Toroidal flux oscillation as possible cause of geomagnetic excursions and reversals

It is proposed that convection driven dynamos operating in planetary cores could be oscillatory even when the oscillations are not directly noticeable from the outside. Examples of dynamo simulations are pointed out that exhibit oscillations in the structure of the azimuthally averaged toroidal magnetic flux while the mean poloidal field shows only variations in its amplitude. In the case of the geomagnetic field, global excursions may be associated with these oscillations. Long period dynamo simulations indicate that the oscillations may cause reversals once in a while. No special attempt has been made to use most realistic parameter values. Nevertheless some similarities between the simulations and the paleomagnetic record can be pointed out. Crown Copyright © 2008.

Hourly values of X, Y and Z geomagnetic components recorded at Coimbra Station (IAGA code: COI) and local geomagnetic K indices for 2007 to 2014

This study examines the performance of series of two geomagnetic indices and series synthesized from a semi-empirical model of magnetospheric currents, in explaining the geomagnetic activity observed at Northern Hemipshere's mid-latitude ground-based stations. We analyse data, for the 2007 to 2014 period, from four magnetic observatories (Coimbra, Portugal; Panagyurishte, Bulgary; Novosibirsk, Russia and Boulder, USA), at geomagnetic latitudes between 40° and 50° N. The quiet daily (QD) variation is firstly removed from the time series of the geomagnetic horizontal component (H) using natural orthogonal components (NOC) tools. We compare the resulting series with series of storm-time disturbance (Dst) and ring current (RC) indices and with H series synthesized from the Tsyganenko and Sitnov (2005, doi:10.1029/2004JA010798) (TS05) semi-empirical model of storm-time geomagnetic field. In the analysis, we separate days with low and high local K-index values. Our results show that NOC models are as efficient as standard models of QD variation in preparing raw data to be compared with proxies, but with much less complexity. For the two stations in Europe, we obtain indication that NOC models could be able to separate ionospheric and magnetospheric contributions. Dst and RC series explain the four observatory H-series successfully, with values for the mean of significant correlation coefficients, from 0.5 to 0.6 during low geomagnetic activity (K less than 4) and from 0.6 to 0.7 for geomagnetic active days (K greater than or equal to 4). With regard to the performance of TS05, our results show that the four observatories separate into two groups: Coimbra and Panagyurishte, in one group, for which the magnetospheric/ionospheric ratio in QD variation is smaller, a dominantly QD ionospheric contribution can be removed and TS05 simulations are the best proxy; Boulder and Novosibirsk,in the other group, for which the ionospheric and magnetospheric contributions in QD variation can not be differentiated and correlations with TS05 series can not be made to improve. The main contributor to magnetospheric QD signal are Birkeland currents. The relatively good success of TS05 model in explaining ground-based irregular geomagnetic activity at mid-latitudes makes it an effective tool to classify storms according to their main sources. For Coimbra and Panagyurishte in particular, where ionospheric and magnetospheric daily contributions seem easier to separate, we can aspire to use the TS05 model for ensemble generation in space weather (SW) forecasting and interpretation of past SW events.

Time evolution of the South Atlantic Magnetic Anomaly

The South Atlantic Magnetic Anomaly (SAMA) is one of the most outstanding anomalies of the geomagnetic field. The SAMA secular variation was obtained and compared to the evolution of other anomalies using spherical harmonic field models for the 1590-2005 period. An analysis of data from four South American observatories shows how this large scale anomaly affected their measurements. Since SAMA is a low total field anomaly, the field was separated into its nondipolar, quadrupolar and octupolar parts. The time evolution of the non-dipole/total, quadrupolar/total and octupolar/total field ratios yielded increasingly high values for the South Atlantic since 1750. The SAMA evolution is compared to the evolution of other large scale surface geomagnetic features like the North and the South Pole and the Siberia High, and this comparison shows the intensity equilibrium between these anomalies in both hemispheres. The analysis of non-dipole fields in historical period suggests that SAMA is governed by (i) quadrupolar field for drift, and (ii) quadrupolar and octupolar fields for intensity and area of influence. Furthermore, our study reinforces the possibility that SAMA may be related to reverse fluxes in the outer core under the South Atlantic region.

The polar ionosphere at Zhongshan Station on May 11, 1999, the day the solar wind almost disappeared

The solar wind almost disappeared on May 11,1999: the solar wind plasma density and' dynamic pressure were less than 1 cm(-3) and 0.1 nPa respectively, while the interplanetary magnetic field was northward. The polar ionospheric data observed by the multi-instruments at Zhongshan Station in Antarctica on such special event day was compared with those of the control day (May 14). It was shown that geomagnetic activity was very quiet on May 11 at Zhongshan. The magnetic pulsation, which usually occurred at about magnetic noon, did not appear. The ionosphere was steady and stratified, and the F-2 layer spread very little. The critical frequency of dayside F-2 layer, f(0)F(2), was larger than that of control day, and the peak of f(0)F(2) appeared 2 hours earlier. The ionospheric drift velocity was less than usual. There were intensive auroral E-s appearing at magnetic noon. All this indicates that the polar ionosphere was extremely quiet and geomagnetic field was much more dipolar on May 11. There were some signatures of auroral substorm before midnight, such as the negative deviation of the geomagnetic H component, accompanied with auroral E-s and weak Pc3 pulsation.

Palaeomagnetic study of a subaerial volcanic ridge (Sao Jorge Island, Azores) for the past 1.3 Myr: evidence for the Cobb Mountain Subchron, volcano flank instability and tectonomagmatic implications

We present a palaeomagnetic study on 38 lava flows and 20 dykes encompassing the past 1.3 Myr on S. Jorge Island (Azores ArchipelagoNorth Atlantic Ocean). The sections sampled in the southeastern and central/western parts of the island record reversed and normal polarities, respectively. They indicate a mean palaeomagnetic pole (81.3 degrees N, 160.7 degrees E, K= 33 and A95= 3.4 degrees) with a latitude shallower than that expected from Geocentric Axial Dipole assumption, suggesting an effect of non-dipolar components of the Earth magnetic field. Virtual Geomagnetic Poles of eight flows and two dykes closely follow the contemporaneous records of the Cobb Mountain Subchron (ODP/DSDP programs) and constrain the age transition from reversed to normal polarity at ca. 1.207 +/- 0.017 Ma. Volcano flank instabilities, probably related to dyke emplacement along an NNWSSE direction, led to southwestward tilting of the lava pile towards the sea. Two spatially and temporally distinct dyke systems have been recognized on the island. The eastern is dominated by NNWSSE trending dykes emplaced before the end of the Matuyama Chron, whereas in the central/western parts the eruptive fissures oriented WNWESE controlled the westward growth of the S. Jorge Island during the Brunhes Chron. Both directions are consistent with the present-day regional stress conditions deduced from plate kinematics and tectonomorphology and suggest the emplacement of dykes along pre-existing fractures. The distinct timing and location of each dyke system likely results from a slight shift of the magmatic source.

Исследование тонкой структуры геомагнитного поля в прибрежной зоне курорта уреки

ცოცხალი ორგანიზმი და გარემო წარმოადგენს ურთულეს დინამიურ სისტემებს, რომლებიც მუდმივ და მჭიდრო ურთიერთკავშირში იმყოფებიან. არსებობს მოსაზრება, რომ ძირითად ფიზიკურ ფაქტორად, რომლის საშუალებითაც ხდება გავლენა ცოცხალ ორგანიზმებზე ხორციელდება მზის აქტივობით გამოწვეული დედამიწის მაგნიტური ველის დაძაბულობის ცვლილებით, და ის ფაქტი, რომ ცოცხალ ორგანიზმებზე მოქმედებენ, როგორც ხელოვნური, ასევე ბუნებრივი მაგნიტური ველები.

A three-dimensional magnetometer for motion sensing of a balloon-carried atmospheric measurement package

An instrument is described which carries three orthogonal geomagnetic field sensors on a standard meteorological balloon package, to sense rapid motion and position changes during ascent through the atmosphere. Because of the finite data bandwidth available over the UHF radio link, a burst sampling strategy is adopted. Bursts of 9s of measurements at 3.6Hz are interleaved with periods of slow data telemetry lasting 25s. Calculation of the variability in each channel is used to determine position changes, a method robust to periods of poor radio signals. During three balloon ascents, variability was found repeatedly at similar altitudes, simultaneously in each of three orthogonal sensors carried. This variability is attributed to atmospheric motions. It is found that the vertical sensor is least prone to stray motions, and that the use of two horizontal sensors provides no additional information over a single horizontal sensor

An overview of Earth's global electric circuit and atmospheric conductivity

The Earth’s global atmospheric electric circuit depends on the upper and lower atmospheric boundaries formed by the ionosphere and the planetary surface. Thunderstorms and electrified rain clouds drive a DC current (∼1 kA) around the circuit, with the current carried by molecular cluster ions; lightning phenomena drive the AC global circuit. The Earth’s near-surface conductivity ranges from 10−7 S m−1 (for poorly conducting rocks) to 10−2 S m−1 (for clay or wet limestone), with a mean value of 3.2 S m−1 for the ocean. Air conductivity inside a thundercloud, and in fair weather regions, depends on location (especially geomagnetic latitude), aerosol pollution and height, and varies from ∼10−14 S m−1 just above the surface to 10−7 S m−1 in the ionosphere at ∼80 km altitude. Ionospheric conductivity is a tensor quantity due to the geomagnetic field, and is determined by parameters such as electron density and electron–neutral particle collision frequency. In the current source regions, point discharge (coronal) currents play an important role below electrified clouds; the solar wind-magnetosphere dynamo and the unipolar dynamo due to the terrestrial rotating dipole moment also apply atmospheric potential differences. Detailed measurements made near the Earth’s surface show that Ohm’s law relates the vertical electric field and current density to air conductivity. Stratospheric balloon measurements launched from Antarctica confirm that the downward current density is ∼1 pA m−2 under fair weather conditions. Fortuitously, a Solar Energetic Particle (SEP) event arrived at Earth during one such balloon flight, changing the observed atmospheric conductivity and electric fields markedly. Recent modelling considers lightning discharge effects on the ionosphere’s electric potential (∼+250 kV with respect to the Earth’s surface) and hence on the fair weather potential gradient (typically ∼130 V m−1 close to the Earth’s surface. We conclude that cloud-to-ground (CG) lightning discharges make only a small contribution to the ionospheric potential, and that sprites (namely, upward lightning above energetic thunderstorms) only affect the global circuit in a miniscule way. We also investigate the effects of mesoscale convective systems on the global circuit.

An overview of Earth’s global electric circuit and atmospheric conductivity

The Earth’s global atmospheric electric circuit depends on the upper and lower atmospheric boundaries formed by the ionosphere and the planetary surface. Thunderstorms and electrified rain clouds drive a DC current (∼1 kA) around the circuit, with the current carried by molecular cluster ions; lightning phenomena drive the AC global circuit. The Earth’s near-surface conductivity ranges from 10−7 S m−1 (for poorly conducting rocks) to 10−2 S m−1 (for clay or wet limestone), with a mean value of 3.2 S m−1 for the ocean. Air conductivity inside a thundercloud, and in fair weather regions, depends on location (especially geomagnetic latitude), aerosol pollution and height, and varies from ∼10−14 S m−1 just above the surface to 10−7 S m−1 in the ionosphere at ∼80 km altitude. Ionospheric conductivity is a tensor quantity due to the geomagnetic field, and is determined by parameters such as electron density and electron–neutral particle collision frequency. In the current source regions, point discharge (coronal) currents play an important role below electrified clouds; the solar wind-magnetosphere dynamo and the unipolar dynamo due to the terrestrial rotating dipole moment also apply atmospheric potential differences. Detailed measurements made near the Earth’s surface show that Ohm’s law relates the vertical electric field and current density to air conductivity. Stratospheric balloon measurements launched from Antarctica confirm that the downward current density is ∼1 pA m−2 under fair weather conditions. Fortuitously, a Solar Energetic Particle (SEP) event arrived at Earth during one such balloon flight, changing the observed atmospheric conductivity and electric fields markedly. Recent modelling considers lightning discharge effects on the ionosphere’s electric potential (∼+250 kV with respect to the Earth’s surface) and hence on the fair weather potential gradient (typically ∼130 V m−1 close to the Earth’s surface. We conclude that cloud-to-ground (CG) lightning discharges make only a small contribution to the ionospheric potential, and that sprites (namely, upward lightning above energetic thunderstorms) only affect the global circuit in a miniscule way. We also investigate the effects of mesoscale convective systems on the global circuit.

Optical characterization of the growth and spatial structure of a substorm onset arc

We present a detailed case study of the characteristics of auroral forms that constitute the first ionospheric signatures of substorm expansion phase onset. Analysis of the optical frequency and along-arc (azimuthal) wave number spectra provides the strongest constraint to date on the potential mechanisms and instabilities in the near-Earth magnetosphere that accompany auroral onset and which precede poleward arc expansion and auroral breakup. We evaluate the frequency and growth rates of the auroral forms as a function of azimuthal wave number to determine whether these wave characteristics are consistent with current models of the substorm onset mechanism. We find that the frequency, spatial scales, and growth rates of the auroral forms are most consistent with the cross-field current instability or a ballooning instability, most likely triggered close to the inner edge of the ion plasma sheet. This result is supportive of a near-Earth plasma sheet initiation of the substorm expansion phase. We also present evidence that the frequency and phase characteristics of the auroral undulations may be generated via resonant processes operating along the geomagnetic field. Our observations provide the most powerful constraint to date on the ionospheric manifestation of the physical processes operating during the first few minutes around auroral substorm onset.

Electron trapping in shear Alfvén waves that power the aurora

Results from 1D Vlasov drift-kinetic plasma simulations reveal how and where auroral electrons are accelerated along Earth’s geomagnetic field. In the warm plasma sheet, electrons become trapped in shear Alfven waves, preventing immediate wave damping. As waves move to regions with larger vTe=vA, their parallel electric field decreases, and the trapped electrons escape their influence. The resulting electron distribution functions compare favorably with in situ observations, demonstrating for the first time a self-consistent link between Alfven waves and electrons that form aurora.

Towards radiocarbon calibration beyond 28ka using speleothems from the Bahamas

We present a new speleothem record of atmospheric Δ14C between 28 and 44 ka that offers considerable promise for resolving some of the uncertainty associated with existing radiocarbon calibration curves for this time period. The record is based on a comprehensive suite of AMS 14C ages, using new low-blank protocols, and U–Th ages using high precision MC-ICPMS procedures. Atmospheric Δ14C was calculated by correcting 14C ages with a constant dead carbon fraction (DCF) of 22.7 ± 5.9%, based on a comparison of stalagmite 14C ages with the IntCal04 (Reimer et al., 2004) calibration curve between 15 and 11 ka. The new Δ14C speleothem record shows similar structure and amplitude to that derived from Cariaco Basin foraminifera (Hughen et al., 2004, 2006), and the match is further improved if the latter is tied to the most recent Greenland ice core chronology (Svensson et al., 2008). These data are however in conflict with a previously published 14C data set for a stalagmite record from the Bahamas — GB-89-24-1 (Beck et al., 2001), which likely suffered from 14C analytical blank subtraction issues in the older part of the record. The new Bahamas speleothem ∆14C data do not show the extreme shifts between 44 and 40 ka reported in the previous study (Beck et al., 2001). Causes for the observed structure in derived atmospheric Δ14C variation based on the new speleothem data are investigated with a suite of simulations using an earth system model of intermediate complexity. Data-model comparison indicates that major fluctuations in atmospheric ∆14C during marine isotope stage 3 is primarily a function of changes in geomagnetic field intensity, although ocean–atmosphere system reorganisation also played a supporting role.

What do cosmogenic isotopes tell us about past solar forcing of climate?

In paleoclimate studies, cosmogenic isotopes are frequently used as proxy indicators of past variations in solar irradiance on centennial and millennial timescales. These isotopes are spallation products of galactic cosmic rays (GCRs) impacting Earth’s atmosphere, which are deposited and stored in terrestrial reservoirs such as ice sheets, ocean sediments and tree trunks. On timescales shorter than the variations in the geomagnetic field, they are modulated by the heliosphere and thus they are, strictly speaking, an index of heliospheric variability rather than one of solar variability. Strong evidence of climate variations associated with the production (as opposed to the deposition) of these isotopes is emerging. This raises a vital question: do cosmic rays have a direct influence on climate or are they a good proxy indicator for another factor that does (such as the total or spectral solar irradiance)? The former possibility raises further questions about the possible growth of air ions generated by cosmic rays into cloud condensation nuclei and/or the modulation of the global thunderstorm electric circuit. The latter possibility requires new understanding about the required relationship between the heliospheric magnetic fields that scatter cosmic rays and the photospheric magnetic fields which modulate solar irradiance.

Temporal and spatial variability of auroral forms in the 10–14 MLT sector: Relationship to plasma convection and solar wind-magnetosphere coupling

Ground-based observations of dayside auroral forms and magnetic perturbations in the arctic sectors of Svalbard and Greenland, in combination with the high-resolution measurements of ionospheric ion drift and temperature by the EISCAT radar, are used to study temporal/spatial structures of cusp-type auroral forms in relation to convection. Large-scale patterns of equivalent convection in the dayside polar ionosphere are derived from the magnetic observations in Greenland and Svalbard. This information is used to estimate the ionospheric convection pattern in the vicinity of the cusp/cleft aurora. The reported observations, covering the period 0700-1130 UT, on January 11, 1993, are separated into four intervals according to the observed characteristics of the aurora and ionospheric convection. The morphology and intensity of the aurora are very different in quiet and disturbed intervals. A latitudinally narrow zone of intense and dynamical 630.0 nm emission equatorward of 75 degrees MLAT, was observed during periods of enhanced antisunward convection in the cusp region. This (type 1 cusp aurora) is considered to be the signature of plasma entry via magnetopause reconnection at low magnetopause latitudes, i.e. the low-latitude boundary layer (LLB I,). Another zone of weak 630.0 nm emission (type 2 cusp aurora) was observed to extend up to high latitudes (similar to 79 degrees MLAT) during relatively quiet magnetic conditions, when indications of reverse (sunward) convection was observed in the dayside polar cap. This is postulated to be a signature of merging between a northward directed IMF (B-z > 0) and the geomagnetic field poleward of the cusp. The coexistence of type 1 and 2 auroras was observed under intermediate circumstances. The optical observations from Svalbard and Greenland were also used to determine the temporal and spatial evolution of type 1 auroral forms, i.e. poleward-moving auroral events occurring in the vicinity of a rotational convection reversal in the early post-noon sector. Each event appeared as a local brightening at the equatorward boundary of the pre-existing type 1 cusp aurora, followed by poleward and eastward expansions of luminosity. The auroral events were associated with poleward-moving surges of enhanced ionospheric convection and F-layer ion temperature as observed by the EISCAT radar in Tromso. The EISCAT ion flow data in combination with the auroral observations show strong evidence for plasma flow across the open/closed field line boundary.

EISCAT observations of ion composition and temperature anisotropy in the high-latitude F-region

The papers by Winser et al. [(1990) J. atmos. terr. Phys.52, 501] and Häggström and Collis [(1990) J. atmos. terr. Phys.52, 519] used plasma flows and ion temperatures, as measured by the EISCAT tristatic incoherent scatter radar, to investigate changes in the ion composition of the ionospheric F-layer at high latitudes, in response to increases in the speed of plasma convection. These studies reported that the ion composition rapidly changed from mainly O+ to almost completely (>90%) molecular ions, following rapid increases in ion drift speed by >1 km s−1. These changes appeared inconsisent with theoretical considerations of the ion chemistry, which could not account for the large fractions of molecular ions inferred from the obsevations. In this paper, we discuss two causes of this discrepancy. First, we reevaluate the theoretical calculations for chemical equilibrium and show that, if we correct the derived temperatures for the effect of the molecular ions, and if we employ more realistic dependences of the reaction rates on the ion temperature, the composition changes derived for the faster convection speeds can be explained. For the Winser et al. observations with the radar beam at an aspect angle of ϕ = 54.7° to the geomagnetic field, we now compute a change to 89% molecular ions in < 2 min, in response to the 3 km s−1 drift. This is broadly consistent with the observations. But for the two cases considered by Häggström and Collis, looking along the field line (ϕ = 0°), we compute the proportion of molecular ions to be only 4 and 16% for the observed plasma drifts of 1.2 and 1.6 km s−1, respectively. These computed proportions are much smaller than those derived experimentally (70 and 90%). We attribute the differences to the effects of non-Maxwellian, anisotropic ion velocity distribution functions. We also discuss the effect of ion composition changes on the various radar observations that report anisotropies of ion temperature.