Ionospheric signatures of pulsed reconnection at the Earth's magnetopause

THE plasma precipitating into the Earth's dayside auroral atmosphere has characteristics which show that it originates from the shocked solar-wind plasma of the magnetosheath1'2. The particles of the magnetosheath plasma precipitate down a funnel-shaped region (cusp) of open field lines resulting from reconnection of the geomagnetic field with the interplanetary magnetic field3. Although the cusp has long been considered a well defined spatial structure maintained by continuous reconnection, it has recently been suggested4–6 that reconnection instead may take place in a series of discontinuous events; this is the ‘pulsating cusp model’. Here we present coordinated radar and satellite observations of a series of discrete, poleward-moving plasma structures that are consistent with the pulsating-cusp model.

Measuring ion temperatures and studying the ion energy balance in the high-latitude ionosphere

Data are presented for a nighttime ion heating event observed by the EISCAT radar on 16 December 1988. In the experiment, the aspect angle between the radar beam and the geomagnetic field was fixed at 54.7°, which avoids any ambiguity in derived ion temperature caused by anisotropy in the ion velocity distribution function. The data were analyzed with an algorithm which takes account of the non-Maxwellian line-of-sight ion velocity distribution. During the heating event, the derived spectral distortion parameter (D∗) indicated that the distribution function was highly distorted from a Maxwellian form when the ion drift increased to 4 km s−1. The true three-dimensional ion temperature was used in the simplified ion balance equation to compute the ion mass during the heating event. The ion composition was found to change from predominantly O4 to mainly molecular ions. A theoretical analysis of the ion composition, using the MSIS86 model and published values of the chemical rate coefficients, accounts for the order-of-magnitude increase in the atomic/molecular ion ratio during the event, but does not successfully explain the very high proportion of molecular ions that was observed.

Transient reconnection: Search for ionospheric signatures

The concept of magnetic reconnection originated with the suggestion by Giovanelli [1946] that particles could be energized during solar flares near nulls in the magnetic field. Hoyle [1949] subsequently proposed that such a process could also act at nulls between the geomagnetic field and the interplanetary magnetic field (IMF) to generate the energized particles responsible for auroral displays. However, the idea of the interconnection of the two magnetic fields, as we know it today, was first presented by Hoyle's student, Dungey [1953, 1961].

Characteristics of the high-latitude trough

The EISCAT radar has provided data for a comprehensive study of the high-latitude trough in electron concentration, which occurs in the auroral zone. In this paper the characteristics of the trough are illustrated, the method of its formation is outlined and important features of the trough are described. A large upward velocity along the geomagnetic field line is shown to play a significant role in the formation of the trough. The large ion-neutral difference velocities which initiate the formation of the trough may also drive the plasma into a non-thermal state which should be taken into account during the analysis of incoherent scatter data.

Analysis of incoherent scatter radar data from non-thermal F-region plasma

A procedure is presented for fitting incoherent scatter radar data from non-thermal F-region ionospheric plasma, using theoretical spectra previously predicted. It is found that values of the shape distortion factor D∗, associated with deviations of the ion velocity distribution from a Maxwellian distribution, and ion temperatures can be deduced (the results being independent of the path of iteration) if the angle between the line-of-sight and the geomagnetic field is larger than about 15–20°. The procedure can be used with one or both of two sets of assumptions. These concern the validity of the adopted model for the line-of-sight ion velocity distribution in the one case or for the full three-dimensional ion velocity distribution function in the other. The distribution function employed was developed to describe the line-of-sight velocity distribution for large aspect angles, but both experimental data and Monte Carlo simulations indicate that the form of the field-perpendicular distribution can also describe the distribution at more general aspect angles. The assumption of this form for the line-of-sight velocity distribution at a general aspect angle enables rigorous derivation of values of the one-dimensional, line-of-sight ion temperature. With some additional assumptions (principally that the field-parallel distribution is always Maxwellian and there is a simple relationship between the ion temperature anisotropy and the distortion of the field-perpendicular distribution from a Maxwellian), fits to data for large aspect angles enable determination of line-of-sight temperatures at all aspect angles and hence, of the average ion temperature and the ion temperature anisotropy. For small aspect angles, the analysis is restricted to the determination of the line-of-sight ion temperature because the theoretical spectrum is insensitive to non-thermal effects when the plasma is viewed along directions almost parallel to the magnetic field. This limitation is expected to apply to any realistic model of the ion velocity distribution function and its consequences are discussed. Fit strategies which allow for mixed ion composition are also considered. Examples of fits to data from various EISCAT observing programmes are presented.

Analysis of incoherent scatter spectra from non-Maxwellian plasma

Incoherent scatter data from non-thermal F-region ionospheric plasma are analysed, using theoretical spectra predicted by Raman et al. It is found that values of the semi-empirical drift parameter D∗, associated with deviations of the ion velocity distribution from a Maxwellian, and the plasma temperatures can be rigorously deduced (the results being independent of the path of iteration) if the angle between the line-of-sight and the geomagnetic field is larger than about 15–20 degrees. For small aspect angles, the deduced value of the average (or 3-D) ion temperature remains ambiguous and the analysis is restricted to the determination of the line-of-sight temperature because the theoretical spectrum is insensitive to non-thermal effects when the plasma is viewed along directions almost parallel to the magnetic field. This limitation is expected to apply to any realistic model of the ion velocity distribution, and its consequences are discussed. Fit strategies which allow for mixed ion composition are also considered. Examples of fits to data from various EISCAT observing programmes are presented.

Pressure-driven magnetopause motions and attendant response on the ground

The terrestrial magnetopause suffered considerable sudden changes in its location on 9–10 September 1978. These magnetopause motions were accompanied by disturbances of the geomagnetic field on the ground. We present a study of the magnetopause motions and the ground magnetic signatures using, for the latter, 10 s averaged data from 14 high latitude ground magnetometer stations. Observations in the solar wind (from IMP 8) are employed and the motions of the magnetopause are monitored directly by the spacecraft ISEE 1 and 2. With these coordinated observations we are able to show that it is the sudden changes in the solar wind dynamic pressure that are responsible for the disturbances seen on the ground. At some ground stations we see evidence of a “ringing” of the magnetospheric cavity, while at others only the initial impulse is evident. We note that at some stations field perturbations closely match the hypothesized ground signatures of flux transfer events. In accordance with more recent work in the area (e.g. Potemra et al., 1989, J. geophys. Res., in press), we argue that causes other than impulsive reeonnection may produce the twin ionospheric flow vortex originally proposed as a flux transfer even signature.

The influence of anisotropic F region ion velocity distributions on ionospheric ion outflows into the magnetosphere

The contribution to the field-aligned ionospheric ion momentum equation, due to coupling between pressure anisotropy and the inhomogeneous geomagnetic field, is investigated. We term this contribution the “hydrodynamic mirror force” and investigate its dependence on the ion drift and the resulting deformations of the ion velocity distribution function from an isotropic form. It is shown that this extra upforce increases rapidly with ion drift relative to the neutral gas but is not highly dependent on the ion-neutral collision model employed. An example of a burst of flow observed by EISCAT, thought to be the ionospheric signature of a flux transfer event at the magnetopause, is studied in detail and it is shown that the nonthermal plasma which results is subject to a hydrodynamic mirror force which is roughly 10% of the gravitational downforce. In addition, predictions by the coupled University College London-Sheffield University model of the ionosphere and thermosphere show that the hydrodynamic mirror force in the auroral oval is up to 3% of the gravitational force for Kp of about 3, rising to 10% following a sudden increase in cross-cap potential. The spatial distribution of the upforce shows peaks in the cusp region and in the post-midnight auroral oval, similar to that of observed low-energy heavy ion flows from the ionosphere into the magnetosphere. We suggest the hydrodynamic mirror force may modulate these outflows by controlling the supply of heavy ions to regions of ion acceleration and that future simulations of the effects of Joule heating on ion outflows should make allowance for it.

Vertical ion flow in the polar ionosphere

Outflowing ions from the polar ionosphere fall into two categories: the classical polar wind and the suprathermal ion flows. The flows in both these categories vary a great deal with altitude. The classical polar wind is supersonic at high altitude: at ∼3 RE geocentric, the observed polar wind is H+ dominated and has a Mach number of 2.5–5.1. At 400–600 km, thermal and suprathermal upward O+ ion fluxes frequently occur at the poleward edge of the nightside auroral oval during magnetically active times. Above 500 km, ions are accelerated transverse to the local geomagnetic field. At 1400 km, transversely accelerated ions are frequently observed in winter nights but rarely appear in the summer. In the dayside cleft above ∼2000 km, ions of all species are transversely heated and upwell with significant number and heat fluxes, forming a cleft ion fountain as they convect across the polar cap. Upwelling ions are observed most (least) frequently in the summer (winter). At yet higher altitudes, energetic (>10 eV to several kiloelectron volts) upflowing H+ and O+ ions are frequently observed, their active time occurrence frequency being as high as 0.7 at auroral latitudes and 0.3 in the polar cap. Their composition, intensity, and angular characteristics vary quantitatively with solar activity, being O+ dominant and more intense near solar maximum. Their resulting ion outflow is dominated by ions below 1 keV and reaches 3.5×10^26 O+ and 7×10^25 H+ ions s^{−1} at magnetically active times (Kp≥5) near solar maximum. In comparison, the estimated polar wind ion outflow at times of moderate solar activity is 7×10^25H+ and 4×10^24 He+ ions s^{−1}. The estimated <10-eV cleft ion fountain flow is 3.8×10^25 O+ and 8.6×10^23 H+ ions s^{−1} near solar maximum.

Scattered power from non-thermal, F-region plasma observed by EISCAT—evidence for coherent echoes?

Three rapid, poleward bursts of plasma flow, observed by the U.K.-POLAR EISCAT experiment, are studied in detail. In all three cases the large ion velocities (> 1 kms−1) are shown to drive the ion velocity distribution into a non-Maxwellian form, identified by the characteristic shape of the observed spectra and the fact that analysis of the spectra with the assumption of a Maxwellian distribution leads to excessive rises in apparent ion temperature, and an anticorrelation of apparent electron and ion temperatures. For all three periods the total scattered power is shown to rise with apparent ion temperature by up to 6 dB more than is expected for an isotropic Maxwellian plasma of constant density and by an even larger factor than that expected for non-thermal plasma. The anomalous increases in power are only observed at the lower altitudes (< 300 km). At greater altitudes the rise in power is roughly consistent with that simulated numerically for homogeneous, anisotropic, non-Maxwellian plasma of constant density, viewed using the U.K.-POLAR aspect angle. The spectra at times of anomalously high power are found to be asymmetric, showing an enhancement near the downward Doppler-shifted ion-acoustic frequency. Although it is not possible to eliminate completely rapid plasma density fluctuations as a cause of these power increases, such effects cannot explain the observed spectra and the correlation of power and apparent ion temperature without an unlikely set of coincidences. The observations are made along a beam direction which is as much as 16.5° from orthogonality with the geomagnetic field. Nevertheless, some form of coherent-like echo contamination of the incoherent scatter spectrum is the most satisfactory explanation of these data.

Non-Maxwellian ion velocity distributions observed using EISCAT

Recent observations from the EISCAT incoherent scatter radar have revealed bursts of poleward ion flow in the dayside auroral ionosphere which are consistent with the ionospheric signature of flux transfer events at the magnetopause. These bursts frequently contain ion drifts which exceed the neutral thermal speed and, because the neutral thermospheric wind is incapable of responding sufficiently rapidly, toroidal, non-Maxwellian ion velocity distributions are expected. The EISCAT observations are made with high time resolution (15 seconds) and at a large angle to the geomagnetic field (73.5°), allowing the non-Maxwellian nature of the distribution to be observed remotely for the first time. The observed features are also strongly suggestive of a toroidal distribution: characteristic spectral shape, increased scattered power (both consistent with reduced Landau damping and enhanced electric field fluctuations) and excessively high line-of-sight ion temperatures deduced if a Maxwellian distribution is assumed. These remote sensing observations allow the evolution of the distributions to be observed. They are found to be non-Maxwellian whenever the ion drift exceeds the neutral thermal speed, indicating that such distributions can exist over the time scale of the flow burst events (several minutes).

Paleomagnetism and rock magnetism of the Neoproterozoic Itajai Basin of the Rio de la Plata craton (Brazil): Cambrian to Cretaceous widespread remagnetizations of South America

A detailed rock magnetic and paleomagnetic study was performed on samples from the Neoproterozoic Itajai Basin in the state of Santa Catarina, Brazil, in order to better constrain the paleogeographic evolution of the Rio de la Plata craton between 600 and 550 Ma. However, rock magnetic properties typical of remagnetized rocks and negative response in the fold test indicated that these rocks carried a secondary chemical remanent magnetization. After detailed AF and thermal cleaning, almost all samples showed a normal polarity characteristic remanent magnetization component close to the present geomagnetic field. The main magnetic carriers are magnetite and hematite, probably of authigenic origin. The mean paleomagnetic pole of the ltajai Basin is located at Plat= -84 degrees, Plong = 97.5 degrees (A95 = 2 degrees) and overlaps the lower Cretaceous segment of the apparent polar wander path of South America, suggesting a cause and effect with the opening of the South Atlantic Ocean. A compilation of remagnetized paleomagnetic poles from South America is presented that highlights the superposition of several large-scale remagnetization events between the Cambrian and the Cretaceous. It is suggested that some paleomagnetic poles used to calibrate the APWP of Gondwana at Precambrian times need to be revised; the indication of remagnetized areas in southern South America may offer some help in the selection of sites for future paleomagnetic investigations in Precambrian rocks. (C) 2011 International Association for Gondwana Research. Published by Elsevier B.V. All rights reserved.

Upper limit on the cosmic-ray photon flux above 10(19) eV using the surface detector of the Pierre Auger Observatory

A method is developed to search for air showers initiated by photons using data recorded by the surface detector of the Auger Observatory. The approach is based on observables sensitive to the longitudinal shower development, the signal risetime and the curvature of the shower front. Applying this method to the data, tipper limits on the flux of photons of 3.8 x 10(-3), 2.5 x 10(-3), and 2.2 x 10(-3) km(-2) sr(-1) yr(-1) above 10(19) eV, 2 x 10(19) eV, and 4 x 10(19) eV are derived, with corresponding limits on the fraction of photons being 2.0%, 5.1%, and 31% (all limits at 95% c.l.). These photon limits disfavor certain exotic models of sources of cosmic rays. The results also show that the approach adopted by the Auger Observatory to calibrate the shower energy is not strongly biased by a contamination from photons. (C) 2008 Elsevier B.V. All rights reserved.

RINEX_HO: second- and third-order ionospheric corrections for RINEX observation files

When GNSS receivers capable of collecting dual-frequency data are available, it is possible to eliminate the first-order ionospheric effect in the data processing through the ionosphere-free linear combination. However, the second- and third-order ionospheric effects still remain. The first-, second- and third-order ionospheric effects are directly proportional to the total electron content (TEC), although the second- and third-order effects are influenced, respectively, by the geomagnetic field and the maximum electron density. In recent years, the international scientific community has given more attention to these kinds of effects and some works have shown that for high precision GNSS positioning these effects have to be taken into consideration. We present a software tool called RINEX_HO that was developed to correct GPS observables for second- and third-order ionosphere effects. RINEX_HO requires as input a RINEX observation file, then computes the second- and third-order ionospheric effects, and applies the corrections to the original GPS observables, creating a corrected RINEX file. The mathematical models implemented to compute these effects are presented, as well as the transformations involving the earth's magnetic field. The use of TEC from global ionospheric maps and TEC calculated from raw pseudorange measurements or pseudoranges smoothed by phase is also investigated.

Higher order ionospheric effects in GNSS positioning in the European region

After removal of the Selective Availability in 2000, the ionosphere became the dominant error source for Global Navigation Satellite Systems (GNSS), especially for the high-accuracy (cm-mm) demanding applications like the Precise Point Positioning (PPP) and Real Time Kinematic (RTK) positioning.The common practice of eliminating the ionospheric error, e. g. by the ionosphere free (IF) observable, which is a linear combination of observables on two frequencies such as GPS L1 and L2, accounts for about 99% of the total ionospheric effect, known as the first order ionospheric effect (Ion1). The remaining 1% residual range errors (RREs) in the IF observable are due to the higher - second and third, order ionospheric effects, Ion2 and Ion3, respectively. Both terms are related with the electron content along the signal path; moreover Ion2 term is associated with the influence of the geomagnetic field on the ionospheric refractive index and Ion3 with the ray bending effect of the ionosphere, which can cause significant deviation in the ray trajectory (due to strong electron density gradients in the ionosphere) such that the error contribution of Ion3 can exceed that of Ion2 (Kim and Tinin, 2007).The higher order error terms do not cancel out in the (first order) ionospherically corrected observable and as such, when not accounted for, they can degrade the accuracy of GNSS positioning, depending on the level of the solar activity and geomagnetic and ionospheric conditions (Hoque and Jakowski, 2007). Simulation results from early 1990s show that Ion2 and Ion3 would contribute to the ionospheric error budget by less than 1% of the Ion1 term at GPS frequencies (Datta-Barua et al., 2008). Although the IF observable may provide sufficient accuracy for most GNSS applications, Ion2 and Ion3 need to be considered for higher accuracy demanding applications especially at times of higher solar activity.This paper investigates the higher order ionospheric effects (Ion2 and Ion3, however excluding the ray bending effects associated with Ion3) in the European region in the GNSS positioning considering the precise point positioning (PPP) method. For this purpose observations from four European stations were considered. These observations were taken in four time intervals corresponding to various geophysical conditions: the active and quiet periods of the solar cycle, 2001 and 2006, respectively, excluding the effects of disturbances in the geomagnetic field (i.e. geomagnetic storms), as well as the years of 2001 and 2003, this time including the impact of geomagnetic disturbances. The program RINEX_HO (Marques et al., 2011) was used to calculate the magnitudes of Ion2 and Ion3 on the range measurements as well as the total electron content (TEC) observed on each receiver-satellite link. The program also corrects the GPS observation files for Ion2 and Ion3; thereafter it is possible to perform PPP with both the original and corrected GPS observation files to analyze the impact of the higher order ionospheric error terms excluding the ray bending effect which may become significant especially at low elevation angles (Ioannides and Strangeways, 2002) on the estimated station coordinates.