Modeling the observed proton aurora and ionospheric convection responses to changes in the IMF clock angle: 1. Persistence of cusp proton aurora

We employ a numerical model of cusp ion precipitation and proton aurora emission to fit variations of the peak Doppler-shifted Lyman-a intensity observed on 26 November 2000 by the SI-12 channel of the FUV instrument on the IMAGE satellite. The major features of this event appeared in response to two brief swings of the interplanetary magnetic field (IMF) toward a southward orientation. We reproduce the observed spatial distributions of this emission on newly opened field lines by combining the proton emission model with a model of the response of ionospheric convection. The simulations are based on the observed variations of the solar wind proton temperature and concentration and the interplanetary magnetic field clock angle. They also allow for the efficiency, sampling rate, integration time and spatial resolution of the FUV instrument. The good match (correlation coefficient 0.91, significant at the 98% level) between observed and modeled variations confirms the time constant (about 4 min) for the rise and decay of the proton emissions predicted by the model for southward IMF conditions. The implications for the detection of pulsed magnetopause reconnection using proton aurora are discussed for a range of interplanetary conditions.

Ionospheric convection response to slow, strong variations in a northward interplanetary magnetic field: A case study for January 14, 1988

We analyze ionospheric convection patterns over the polar regions during the passage of an interplanetary magnetic cloud on January 14, 1988, when the interplanetary magnetic field (IMF) rotated slowly in direction and had a large amplitude. Using the assimilative mapping of ionospheric electrodynamics (AMIE) procedure, we combine simultaneous observations of ionospheric drifts and magnetic perturbations from many different instruments into consistent patterns of high-latitude electrodynamics, focusing on the period of northward IMF. By combining satellite data with ground-based observations, we have generated one of the most comprehensive data sets yet assembled and used it to produce convection maps for both hemispheres. We present evidence that a lobe convection cell was embedded within normal merging convection during a period when the IMF By and Bz components were large and positive. As the IMF became predominantly northward, a strong reversed convection pattern (afternoon-to-morning potential drop of around 100 kV) appeared in the southern (summer) polar cap, while convection in the northern (winter) hemisphere became weak and disordered with a dawn-to-dusk potential drop of the order of 30 kV. These patterns persisted for about 3 hours, until the IMF rotated significantly toward the west. We interpret this behavior in terms of a recently proposed merging model for northward IMF under solstice conditions, for which lobe field lines from the hemisphere tilted toward the Sun (summer hemisphere) drape over the dayside magnetosphere, producing reverse convection in the summer hemisphere and impeding direct contact between the solar wind and field lines connected to the winter polar cap. The positive IMF Bx component present at this time could have contributed to the observed hemispheric asymmetry. Reverse convection in the summer hemisphere broke down rapidly after the ratio |By/Bz| exceeded unity, while convection in the winter hemisphere strengthened. A dominant dawn-to-dusk potential drop was established in both hemispheres when the magnitude of By exceeded that of Bz, with potential drops of the order of 100 kV, even while Bz remained northward. The later transition to southward Bz produced a gradual intensification of the convection, but a greater qualitative change occurred at the transition through |By/Bz| = 1 than at the transition through Bz = 0. The various convection patterns we derive under northward IMF conditions illustrate all possibilities previously discussed in the literature: nearly single-cell and multicell, distorted and symmetric, ordered and unordered, and sunward and antisunward.

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.

The statistical cusp: a flux transfer event model

Traditionally, the cusp has been described in terms of a time-stationary feature of the magnetosphere which allows access of magnetosheath-like plasma to low altitudes. Statistical surveys of data from low-altitude spacecraft have shown the average characteristics and position of the cusp. Recently, however, it has been suggested that the ionospheric footprint of flux transfer events (FTEs) may be identified as variations of the “cusp” on timescales of a few minutes. In this model, the cusp can vary in form between a steady-state feature in one limit and a series of discrete ionospheric FTE signatures in the other limit. If this time-dependent cusp scenario is correct, then the signatures of the transient reconnection events must be able, on average, to reproduce the statistical cusp occurrence previously determined from the satellite observations. In this paper, we predict the precipitation signatures which are associated with transient magnetopause reconnection, following recent observations of the dependence of dayside ionospheric convection on the orientation of the IMF. We then employ a simple model of the longitudinal motion of FTE signatures to show how such events can easily reproduce the local time distribution of cusp occurrence probabilities, as observed by low-altitude satellites. This is true even in the limit where the cusp is a series of discrete events. Furthermore, we investigate the existence of double cusp patches predicted by the simple model and show how these events may be identified in the data.

The ionospheric signatures of flux transfer events and solar wind dynamic pressure changes

The generation of flow and current vortices in the dayside auroral ionosphere has been predicted for two processes ocurring at the dayside magnetopause. The first of these mechanisms is time-dependent magnetic reconnection, in “flux transfer events” (FTEs); the second is the action of solar wind dynamic pressure changes. The ionospheric flow signature of an FTE should be a twin vortex, with the mean flow velocity in the central region of the pattern equal to the velocity of the pattern as a whole. On the other hand, a pulse of enhanced or reduced dynamic pressure is also expected to produce a twin vortex, but with the central plasma flow being generally different in speed from, and almost orthogonal to, the motion of the whole pattern. In this paper, we make use of this distinction to discuss recent observations of vortical flow patterns in the dayside auroral ionosphere in terms of one or other of the proposed mechanisms. We conclude that some of the observations reported are consistent only with the predicted signature of FTEs. We then evaluate the dimensions of the open flux tubes required to explain some recent simultaneous radar and auroral observations and infer that they are typically 300 km in north–south extent but up to 2000 km in longitudinal extent (i.e., roughly 5 hours of MLT). Hence these observations suggest that recent theories of FTEs which invoke time-varying reconnection at an elongated neutral line may be correct. We also present some simultaneous observations of the interplanetary magnetic field (IMF) and solar wind dynamic pressure (observed using the IMP8 satellite) and the ionospheric flow (observed using the EISCAT radar) which are also only consistent with the FTE model. We estimate that for continuously southward IMF ( ≈ 5 nT) these FTEs contribute about 30 kV to the mean total transpolar voltage (∼30%).

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.

Ionospheric ion upwelling in the wake of flux transfer events at the dayside magnetopause

The effects of flux transfer events (FTE) on the dayside auroral ionosphere are studied, using a simple twin-vortex model of induced ionospheric plasma flow. It is shown that the predicted and observed velocities of these flows are sufficient to drive nonthermal plasma in the F region, not only within the newly opened flux tube of the FTE, but also on the closed, or "old" open, field lines around it. In fact, with the expected poleward neutral wind, the plasma is more highly nonthermal on the flanks of, but outside, the open flux tube: EISCAT observations indicate that plasma is indeed driven into nonthermal distributions in these regions. The nonthermal plasma is thereby subject to additional upforce due to the resulting ion temperature anisotropy and transient expansion due to Joule heating and also to ion accelerations associated with the FTE field aligned current system. Any upflows produced on closed field lines in the vicinity of the FTE are effectively bunched-up in the "wake" of the FTE. Observations from the AMPTE-UKS satellite at the magnetopause reveal ion upflows of energy ∼100 eV flowing out from the ionosphere on closed field lines which are only found in the wake of the FTE. Such flows are also only found shortly after two, out of all the FTEs observed by AMPTE-UKS. The outflow from the ionosphere is two orders of magnitude greater than predicted for the "classical" polar wind. It is shown that such ionospheric ion flows are only expected in association with FTEs on the magnetopause which are well removed from the sub-solar point-either towards dusk or, as in the UKS example discussed here, towards dawn. It is suggested that such ionospheric ions will only be observed if the center of the FTE open flux tube passes very close to the satellite. Consequently, we conclude the ion upflows presented here are probably driven by the second of two possible source FTEs and are observed at the satellite with a lag after the FTE which is less than their time-of-flight.

The modelled occurrence of non-thermal plasma in the ionospheric F-region and the possible consequences for ion outflows into the magnetosphere

A global, time-dependent, three-dimensional, coupled ionosphere-thermosphere model is used to predict the spatial distribution of non-thermal plasma in the F-layer. It is shown that, even for steady-state conditions with Kp as low as 3, the difference between the ion and neutral velocities often exceeds the neutral thermal speed by a factor, D', which can be as large as 4. Theoretically, highly non-Maxwellian, and probably toroidal, ion velocity distributions are expected when D' exceeds about 1.5. The lack of response of the neutral winds to sunward ion drifts in the dawn sector of the auroral oval cause this to be the region most likely to contain toroidal distributions. The maximum in D' is found in the throat region of the convection pattern, where the strong neutral winds of the afternoon sector meet the eastward ion flows of the morning sector. These predictions are of interest, not only to radar scientists searching for non-thermal ionospheric plasma, but also as one possible explanation of the initial heating and upward flows of ions in the cleft ion fountain and nightside auroral oval, both of which are a major source of plasma for the magnetosphere.

The cleft ion fountain: a two-dimensional kinetic model

The transport of ionospheric ions from a source in the polar cleft ionosphere through the polar magnetosphere is investigated using a two-dimensional, kinetic, trajectory-based code. The transport model includes the effects of gravitation, longitudinal magnetic gradient force, convection electric fields, and parallel electric fields. Individual ion trajectories as well as distribution functions and resulting bulk parameters of density, parallel average energy, and parallel flux for a presumed cleft ionosphere source distribution are presented for various conditions to illustrate parametrically the dependences on source energies, convection electric field strengths, ion masses, and parallel electric field strengths. The essential features of the model are consistent with the concept of a cleft-based ion fountain supplying ionospheric ions to the polar magnetosphere, and the resulting plasma distributions and parameters are in general agreement with recent low-energy ion measurements from the DE 1 satellite.

Simple M-factor algorithm for improved estimation of the basic maximum usable frequency of radio waves reflected from the ionospheric F-region

The equations of Milsom are evaluated, giving the ground range and group delay of radio waves propagated via the horizontally stratified model ionosphere proposed by Bradley and Dudeney. Expressions for the ground range which allow for the effects of the underlying E- and F1-regions are used to evaluate the basic maximum usable frequency or M-factors for single F-layer hops. An algorithm for the rapid calculation of the M-factor at a given range is developed, and shown to be accurate to within 5%. The results reveal that the M(3000)F2-factor scaled from vertical-incidence ionograms using the standard URSI procedure can be up to 7.5% in error. A simple addition to the algorithm effects a correction to ionogram values to make these accurate to 0.5%.

The bottomside mid-latitude ionospheric trough

The measured power losses and Doppler shifts of h.f. radio waves propagated over a long, west-east, sub-auroral path are found to exhibit features which cannot be explained by simple predictions and models. Both theory and the limited available data indicate that a bottomside F2-layer depletion should be present below the topside mid-latitude trough. Introducing this into the models (using the mean statistical positions of the trough deduced from Alouette I and II soundings) is shown to explain many of these features. From the Doppler shifts and a simple ray-tracing model the height of the depleted F2-peak inside the trough is deduced to be greater than its value outside the trough by an amount of the order of only 30–80 km.

The wavelet method as an alternative for reducing ionospheric effects from single-frequency GPS receivers

The ionospheric effect is one of the major errors in GPS data processing over long baselines. As a dispersive medium, it is possible to compute its influence on the GPS signal with the ionosphere-free linear combination of L1 and L2 observables, requiring dual-frequency receivers. In the case of single-frequency receivers, ionospheric effects are either neglected or reduced by using a model. In this paper, an alternative for single-frequency users is proposed. It involves multiresolution analysis (MRA) using a wavelet analysis of the double-difference observations to remove the short- and medium-scale ionosphere variations and disturbances, as well as some minor tropospheric effects. Experiments were carried out over three baseline lengths from 50 to 450 km, and the results provided by the proposed method were better than those from dual-frequency receivers. The horizontal root mean square was of about 0.28 m (1 sigma).

Stochastic modelling considering ionospheric scintillation effects on GNSS relative and point positioning

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Second and Third Order Ionospheric Effects on GNSS Positioning: A Case Study in Brazil

The Global Positioning System (GPS) transmits signals in two frequencies. It allows the correction of the first order ionospheric effect by using the ionosphere free combination. However, the second and third order ionospheric effects, which combined may cause errors of the order of centimeters in the GPS measurements, still remain. In this paper the second and third order ionospheric effects, which were taken into account in the GPS data processing in the Brazilian region, were investigated. The corrected and not corrected GPS data from these effects were processed in the relative and precise point positioning (PPP) approaches, respectively, using Bernese V5.0 software and the PPP software (GPSPPP) from NRCAN (Natural Resources Canada). The second and third order corrections were applied in the GPS data using an in-house software that is capable of reading a RINEX file and applying the corrections to the GPS observables, creating a corrected RINEX file. For the relative processing case, a Brazilian network with long baselines was processed in a daily solution considering a period of approximately one year. For the PPP case, the processing was accomplished using data collected by the IGS FORT station considering the period from 2001 to 2006 and a seasonal analysis was carried out, showing a semi-annual and an annual variation in the vertical component. In addition, a geographical variation analysis in the PPP for the Brazilian region has confirmed that the equatorial regions are more affected by the second and third order ionospheric effects than other regions.

High precision GPS network with precise ephemerides and earth body tide model

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)