B. A. Emery

Magnetosphere‐ionosphere‐thermosphere coupling: Effect of neutral winds on energy transfer and field‐aligned current

The assimilative mapping of ionospheric electrodynamics (AMIE) algorithm has been applied to derive the realistic time-dependent large-scale global distributions of the ionospheric convection and particle precipitation during a recent Geospace Environment Modeling (GEM) campaign period: March 28-29, 1992. The AMIE outputs are then used as the inputs of the National Center for Atmospheric Research thermosphere-ionosphere general circulation model to estimate the electrodynamic quantities in the ionosphere and thermosphere. It is found that the magnetospheric electromagnetic energy dissipated in the high-latitude ionosphere is mainly converted into Joule heating, with only a small fraction (6%) going to acceleration of thermospheric neutral winds. Our study also reveals that the thermospheric winds can have significant influence on the ionospheric electrodynamics. On the average for these 2 days, the neutral winds have approximately a 28% negative effect on Joule heating and approximately a 27% negative effect on field-aligned currents. The field-aligned currents driven by the neutral wind flow in the opposite direction to those driven by the plasma convection. On the average, the global electromagnetic energy input is about 4 times larger than the particle energy input.

If the Sun is so quiet, why is the Earth ringing? A comparison of two solar minimum intervals

coronal holes lingered even as the sunspots disappeared. Consequently, for the months surrounding the WHI campaign, strong, long, and recurring high-speed streams in the solar wind intercepted the Earth in contrast to the weaker and more sporadic streams that occurred around the time of last cycle’s WSM campaign. In response, geospace and upper atmospheric parameters continued to ring with the periodicities of the solar wind in a manner that was absent last cycle minimum, and the flux of relativistic electrons in the Earth’s outer radiation belt was elevated to levels more than three times higher in WHI than in WSM. Such behavior could not have been predicted using sunspot numbers alone, indicating the importance of considering variation within and between solar minima in analyzing and predicting space weather responses at the Earth during solar quiet intervals, as well as in interpreting the Sun’s past behavior as preserved in geological and historical records.

Interhemispheric asymmetry of the high-latitude ionospheric convection pattern

The assimilative mapping of ionospheric electrodynamics technique has been used to derive the large-scale high-latitude ionospheric convection patterns simultaneously in both northern and southern hemispheres during the period of January 27-29, 1992. When the interplanetary magnetic field (IMF) Bz component is negative, the convection patterns in the southern hemisphere are basically the mirror images of those in the northern hemisphere. The total cross-polar-cap potential drops in the two hemispheres are similar. When Bz is positive and |By| > Bz, the convection configurations are mainly determined by By and they may appear as normal “two-cell” patterns in both hemispheres much as one would expect under southward IMF conditions. However, there is a significant difference in the cross-polar-cap potential drop between the two hemispheres, with the potential drop in the southern (summer) hemisphere over 50% larger than that in the northern (winter) hemisphere. As the ratio of |By|/Bz decreases (less than one), the convection configuration in the two hemispheres may be significantly different, with reverse convection in the southern hemisphere and weak but disturbed convection in the northern hemisphere. By comparing the convection patterns with the corresponding spectrograms of precipitating particles, we interpret the convection patterns in terms of the concept of merging cells, lobe cells, and viscous cells. Estimates of the “merging cell” potential drops, that is, the potential ascribed to the opening of the dayside field lines, are usually comparable between the two hemispheres, as they should be. The “lobe cell” provides a potential between 8.5 and 26 k V and can differ greatly between hemispheres, as predicted. Lobe cells can be significant even for southward IMF, if |By| > |Bz|. To estimate the potential drop of the “viscous cells,” we assume that the low-latitude boundary layer is on closed field lines. We find that this potential drop varies from case to case, with a typical value of 10 kV. If the source of these cells is truly a viscous interaction at the flank of the magnetopause, the process is likely spatially and temporally varying rather than steady state.

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.

An overview of the early November 1993 geomagnetic storm

This paper describes the development of a major space storm during November 2-11, 1993. We discuss the history of the contributing high-speed stream, the powerful combination of solar wind transients and a corotating interaction region which initiated the storm, the high-speed flow which prolonged the storm and the near-Earth manifestations of the storm. The 8-day storm period was unusually long; the result of a high-speed stream (maximum speed 800 km/s) emanating from a distended coronal hole. Storm onset was accompanied by a compression of the entire dayside magnetopause to within geosynchronous Earth orbit (GEO). For nearly 12 hours the near-Earth environment was in a state of tumult. A super-dense plasma sheet was observed at GEO, and severe spacecraft charging was reported. The effects of electrons precipitating into the atmosphere penetrated into the stratosphere. Subauroral electron content varied by 100% and F layer heights oscillated by 200 km. Equatorial plasma irregularities extended in plumes to heights of 1400 km. Later, energetic particle fluxes at GEO recovered and rose by more than an order of magnitude. A satellite anomaly was reported during the interval of high energetic electron flux. Model results indicate an upper atmospheric temperature increase of 200°K within 24 hours of storm onset. Joule heating for the first 24 hours of the storm was more than 3 times that for typical active geomagnetic conditions. We estimate that total global ionospheric heating for the full storm interval was ∼190 PJ, with 30% of that generated within 24 hours of storm onset.

An ionospheric conductance model based on ground magnetic disturbance data

An attempt is made to construct an improved ionospheric conductance model employing ground magnetic disturbance data as input. For each of the different regions in the auroral electrojets specified by different combinations of horizontal (ΔH) and vertical (ΔZ) magnetic perturbations, as well as by magnetic local time (MLT), an empirical relationship is obtained between the ionospheric conductance deduced from Chatanika radar data and magnetic disturbances from the nearby College magnetic station. The error involved in the empirical formula is generally of the order of 20–50%. However, some sectors are so poorly covered that uncertainty estimates cannot be made. The new formulas are applied to an average magnetic disturbance distribution to deduce the average conductance distribution. This is compared with a conductance model based on electron precipitation data [Hardy et al., 1987], finding good agreement in terms of the magnitude and distribution pattern. Combining our empirical relationships with the empirical formulas proposed by Robinson et al. [1987], the average energy and energy flux of precipitating electrons are also estimated. Notable similarities exist between the global distribution patterns of these and those obtained by Hardy et al. [1985]. It is proposed that the present conductance model can be used to complement more direct measurements in order to obtain the global distribution needed to study the large-scale electrodynamics of the polar ionosphere. Several interesting characteristics about the auroral electrojet system are apparent from the empirical relationship: (1) For a given magnitude of ΔH, the electric field is relatively stronger in the eastward electrojet region than in the westward electrojet region. (2) The electric field plays a greater role in the intensification of electrojet current than the ionospheric conductance does in the poleward half of the westward electrojet, whereas the opposite trend is apparent in the equatorward half. However, no such different roles of the electric field and conductance is noticeable in the eastward electrojet region. (3) The auroral conductance enhancements tend to be largest around midnight, due to more intense particle fluxes there. (4) The mean particle energy depends on MLT but is relatively insensitive to magnetic activity.

Transformation of high‐latitude ionospheric F region patches into blobs during the March 21, 1990, storm

Discrete F region electron density enhancements of a factor of 2 or more have been observed in the high-latitude ionosphere. These enhancements have been termed patches if they occur within the polar cap and blobs if they occur outside of the polar cap. It is important to understand the formation and evolution of these structures because they are associated with large phase and amplitude scintillation in transionospheric radio signals. Blobs are generally thought to result from the breakup of patches as they exit the polar cap; however, this process has not previously been observed. Detailed study of high-latitude ionospheric plasma transport is generally difficult because of the sparseness (spatial and temporal) of electron density and velocity observations. In this paper, we present electron density enhancements measured from the Qaanaaq Digisonde, the Millstone Hill incoherent scatter radar, and the DMSP F8 satellite during a 5-hour interval of the March 21, 1990, storm period and show definitively how a patch is transformed into a blob. We present a new trajectory analysis package that is capable of using ionospheric convection patterns to determine the motion of ionospheric plasma over a period of several hours. The new package uses convection patterns from the Assimilative Mapping of Ionospheric Electrodynamics (AMIE) technique to track the motion of observed patches from one site to another and thus determines where the measured electron density enhancements originated and where they went after being observed. The trajectory analysis also establishes that there is a direct connection between the enhancements observed by the different instruments at different locations. In this case, within ∼4 hours, plasma observed by a Digisonde near the pole is convected through 35° of latitude to the northeastern United States, where it is observed by the Millstone Hill radar, then roughly equal portions are transported westward to Alaska and eastward to Scandinavia where they are observed by the DMSP satellite. This study demonstrates that the changing convection pattern can significantly distort the patch shape and trajectory, and illustrates the high degree of mixing of ionospheric plasma by convection. The changing convection pattern leads to the simultaneous existence of a boundary blob and a subauroral blob which are both observed by the Millstone Hill radar. This work is very relevant to our future ability to specify and forecast ionospheric conditions at high latitudes. It represents a critical step from a merely qualitative ability to model the evolution of patches and blobs to a quantitative ability.

High‐latitude ionospheric electrodynamics as determined by the assimilative mapping of ionospheric electrodynamics procedure for the conjunctive SUNDIAL/ATLAS 1/GEM period of March 28–29, 1992

During the conjunctive SUNDIAL/ATLAS 1/GEM campaign period of March 28–29, 1992, a set of comprehensive data has been collected both from space and from ground. The assimilative mapping of ionospheric electrodynamics (AMIE) procedure is used to derive the large-scale high-latitude ionospheric conductivity, convection, and other related quantities, by combining the various data sets. The period was characterized by several moderate substorm activities. Variations of different ionospheric electrodynamic fields are examined for one substorm interval. The cross-polar-cap potential drop, Joule heating, and field-aligned current are all enhanced during the expansion phase of substorms. The most dramatic changes of these fields are found to be associated with the development of the substorm electrojet in the post midnight region. Variations of global electrodynamic quantities for this 2-day period have revealed a good correlation with the auroral electrojet (AE) index. In this study we have calculated the AE index from ground magnetic perturbations observed by 63 stations located between 55° and 76° magnetic latitudes north and south, which is larger than the standard AE index by about 28% on the average over these 2 days. Different energy dissipation channels have also been estimated. On the average over the 2 days, the total globally integrated Joule heating rate is about 102 GW and the total globally integrated auroral energy precipitation rate is about 52 GW. Using an empirical formula, the ring current energy injection rate is estimated to be 125 GW for a decay time of 3.5 hours, and 85 GW for a decay time of 20 hours. We also find an energy-coupling efficiency of 3% between the solar wind and the magnetosphere for a southward interplanetary magnetic field (IMF) condition.

A magnetosphere-thermosphere-ionosphere electrodynamics general circulation model

A new simulation model of the thermosphere, the ionosphere, and the magnetosphere is presented. This model, which we call the magnetosphere-thermosphere-ionosphere electrodynamics general circulation model (MTIE-GCM), calculates the three-dimensional structure of the thermosphere and the ionosphere, the two-dimensional magnetospheric plasma convection in the equatorial plane of the magnetosphere, and the couplings among the thermosphere, the ionosphere, and the magnetosphere. The electrodynamic couplings induced by the auroral precipitation, the region 2 field-aligned currents, and the neutral winds are self-consistently computed. The major drawback of the model is the use of a dipole magnetic field in the magnetosphere such that we do not expect to reproduce observations exactly but rather to predict only correct orders of magnitude. The MTIE-GCM is run for solar maximum conditions and moderate magnetic activity in a sample case, where the electric potential is imposed in the polar cap but calculated everywhere else, to study the electrodynamic couplings occurring at high latitudes. The neutral winds increase the north–south component and decrease the east–west component of the ionospheric electric field, corresponding to an enhancement of the shielding effect by 10%. The region 2 field-aligned currents and the distribution of the magnetospheric pressure present smaller modifications, suggesting that the magnetosphere acts partly as a current generator. However, the results could possibly be different if the polar cap electric potential were allowed to change due to the neutral winds.

Assimilative mapping of ionospheric electrodynamics in the thermosphere‐ionosphere general circulation model comparisons with global ionospheric and thermospheric observations during the GEM/SUNDIAL period of March 28–29, 1992

Satellite and ground-based observations from March 28 to 29, 1992, were combined in the assimilative mapping of ionospheric electrodynamics (AMIE) procedure to derive realistic global distributions of the auroral precipitation and ionospheric convection which were used as inputs to the National Center for Atmospheric Research (NCAR) thermosphere-ionosphere general circulation model (TIGCM). Comparisons of neutral model winds were made with Fabry-Perot measurements and meridional winds derived from ionosondes. The peak equatorward winds occurred 1–2 hours later in the model. Gravity waves launched from high-latitude Joule heating sources reached the equator in about 2 hours and agreed with observed variations in the height of the maximum electron density (hmF2) and in the meridional winds. Joule heating events produced minima in the O/N2 ratio that moved equatorward and usually westward in longitudinal strips which lasted about a day. Changes in the O/N2 ratio and in the peak electron density (NmF2) were strongly correlated so the observed daytime NmF2 values for stations near 50° magnetic latitude were generally reproduced by AMIE-TIGCM on the second day of the simulation. The AMIE-TIGCM underestimated the electron density after midnight by up to a factor of 2 in midlatitudes, while the modeled F2 layer was about 35 km lower than the observations at midnight. Shifting the model winds 2 hours earlier at night could double the NmF2 at 0400 LT and increase hmF2 by 20 km. NmF2 could also be increased at night by realistically increasing the TIGCM nighttime downward fluxes of O+ at the upper boundary.