Delores J. Knipp

Effects of a high‐density plasma sheet on ring current development during the November 2–6, 1993, magnetic storm

The growth and recovery of the November 2–6, 1993 magnetic storm was simulated using a drift-loss ring current model that was driven by dynamic fluxes at geosynchronous orbit as an outer boundary condition. During the storm main phase, a high-density plasma sheet was observed by the Los Alamos National Laboratory geosynchronous satellites to move into and flow around the inner magnetosphere over a period of ∼12 hours [Borovsky et al., 1997; this issue] during the storm main phase. Densities at the leading edge of this structure reached 3 cm−3 as compared with more typical values <1 cm−3. The factor of 3 change in the plasma sheet density from quiet to active times produced a factor of 3 enhancement in the strength of the simulated ring current. In addition, a short-timescale recovery in the Dst index at 1600 UT on November 4 was driven by changes in the outer boundary condition and appeared even in the absence of collisional losses. An overshoot in the minimum Dst* occurred in the simulated ring current compared with observed values at ∼0200 UT on November 4 and is taken as evidence of a loss process not included in the ring current-atmosphere interaction model (RAM). The storm onset was associated with a compression of the entire dayside magnetopause to within geostationary orbit starting at 2307 UT and continuing for a half hour. It is suggested that a possible additional loss may have resulted as ions drifted to the compressed dayside magnetopause. In fact such losses were found in another simulation of the inner magnetosphere for the same storm by Freeman et al. [1996]. The energy supplied to the inner magnetosphere, relative to the total energy input during this magnetic storm, was examined by comparing two widely used energy input functions, the e parameter [Akasofu, 1981] and the F parameter [Burton et al., 1975] against energy input to the ring current model based on geosynchronous plasma observations at the outer boundary. It is found that the e parameter [Akasofu, 1981] overestimates the ring current energy input compared to the drift-loss model by almost an order of magnitude during the main phase. However, the integrated energy input from e, over the 4 day interval of the storm, is in very good agreement with the total energy input inferred from observations. On the other hand, F more closely approximates the magnitude of the ring current energy input alone as calculated in the drift-loss model. An energy budget is constructed for the storm that shows energy inputs from the solar wind and energy dissipation due to ring current buildup and decay, auroral electron precipitation, Joule heating, ion precipitation, and energy storage in the magnetotail in reasonable balance. The ring current energy input accounts for only 15% of the total dissipated energy in this storm interval. A more complete energy budget that extends to November 11, 1993, was compiled by Knipp et al. [this issue].

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.

Magnetospheric dynamics and mass flow during the November 1993 storm

The National Space Weather Program (NSWP) Storm that occurred in November 1993 is examined with the use of plasma and energetic-particle measurements on three satellites in geosynchronous orbit. Geosynchronous orbit affords a powerful perspective on magnetospheric dynamics since both tail and dipole processes can be regularly seen, as well as nightside and dayside processes. The major magnetospheric regions analyzed before, during, and after this storm are the outer plasmasphere, the ion plasma sheet, the electron plasma sheet, and the outer electron radiation belt. Ionospheric outflows into the magnetosphere are also observed, and during the storm the magnetosheath and the low-latitude boundary layer are both seen briefly. The geosynchronous observations indicate that prior to the storm the magnetosphere was very quiet and the outer plasmasphere was filled out to beyond geosynchronous orbit. Extremely large anisotropies were seen in the ion plasma sheet during a compression phase just prior to storm onset. During the storms main phase the drainage of the outer plasmasphere to the dayside magnetopause was observed, a super dense ion plasma sheet was tracked moving around the dipole, and a superdense electron plasma sheet was seen. The anomalously large plasma pressure on the nightside led to a β > 1 situation at geosynchronous orbit. The β > 1 region spread around the dipole with the super dense ion plasma sheet. The magnetic-field tilt angle at geosynchronous orbit indicated that strong cross-tail currents were present very near the Earth. These currents appear to be associated with plasma diamagnetism. Geosynchronous observations indicate that magnetospheric convection was extremely strong. In the electron plasma sheet, severe spacecraft charging occurred. The density of relativistic electrons was observed to peak very early in the storm, whereas the flux of these relativistic electrons peaked much later in the aftermath of the storm.

A large-scale traveling ionospheric disturbance during the magnetic storm of 15 September 1999

enhancement of GPS total electron content (� 1.0 � 10 16 m � 2 ). Multipoint and imaging observations of these parameters show that the LSTID moved equatorward over Japan with a velocity of � 400–450 m/s. From a comparison with the Sheffield University Plasmasphere-Ionosphere Model (SUPIM) we conclude that an enhancement (250–300 m/s) of poleward neutral wind (that is propagating equatorward) caused these observational features of the LSTID at midlatitudes. To investigate generation of the LSTID by auroral energy input, we have used auroral images obtained by the Polar UVI instrument, magnetic field variations obtained at multipoint ground stations, and the empirical Joule heating rate calculated by the assimilative mapping of ionospheric electrodynamics (AMIE) technique. Intense auroral energy input was observed at 0800–1100 UT (4–6 hours before the LSTID), probably causing equatorward neutral wind at lower latitudes. It is likely that the poleward wind pulse that caused the observed LSTID was generated associated with the cessation of this equatorward wind. The effect of Lorentz force is also discussed. INDEX TERMS: 0310 Atmospheric Composition and Structure: Airglow and aurora; 2427 Ionosphere: Ionosphere/atmosphere interactions (0335); 2435 Ionosphere: Ionospheric disturbances; 2437 Ionosphere: Ionospheric dynamics; 2788 Magnetospheric Physics: Storms and substorms; KEYWORDS: large-scale traveling ionospheric disturbance, thermosphere–ionosphere coupling, magnetic storm, airglow imaging, GPS network, ionosonde

Polar cap index as a proxy for hemispheric Joule heating

The polar cap (PC) index measures the level of geomagnetic activity in the polar cap based on magnetic perturbations from overhead ionospheric currents and distant field-aligned currents on the poleward edge of the nightside auroral oval. Because PC essentially measures the main sources of energy input into the polar cap, we propose to use PC as a proxy for the hemispheric Joule heat production rate (JH). In this study, JH is estimated from the Assimilative Mapping of Ionospheric Electrodynamics (AMIE) procedure. We fit hourly PC values to hourly averages of JH. Using a data base approximately three times larger than studies, we find a quadratic relationship between JH and PC, differentiated by season. A comparison during the November 1993 storm interval with earlier reported methods using the AE index and the cross polar cap potential, shows that the PC-based Joule heating estimate is as equally accurate. Thus the single station PC index appears to provide a quick estimate of, and is an appropriate proxy for, the hemispheric Joule heating rate.

Ionospheric convection response to changing IMF direction

By combining ground-based and satellite-based measurements of ionospheric electric fields, conductivities and magnetic perturbations, we are able to examine the characteristics of instantaneous, ionospheric convection patterns associated with changing directions of the Interplanetary Magnetic Field (IMF). In response to a rapid southward-to-northward turning of the IMF on 23 July 1983, the ionospheric convection reconfigured over a period of 40 minutes. The configuration changed from a conventional two-cell pattern to a contracted four-cell pattern, with reversed convection cells in the high-latitude dayside, associated with a strong potential drop of about 75 kV. Later, in response to a gradual rotation of the IMF from the +Z through the −Y. toward the −Z direction, the nightside cells disappeared and the dawn cell in the reversed pair wrapped around and displaced the dusk cell until a conventional two-cell pattern was reestablished, largely in accord with the qualitative model of Crooker [1988]. Our results suggest that multiple cells can arise as a result of strong southward to northward transitions in the IMF. They appear to persist for sometime thereafter.

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.

TIME DEPENDENT THERMOSPHERIC NEUTRAL RESPONSE TO THE 2-11 NOVEMBER 1993 STORM PERIOD

Abstract Many satellite and ground-based observations from 2–11 November 1993 werecombined in the Assimilative Mapping of Ionospheric Electrodynamics (AMIE) procedure toderive realistic time dependent global distributions of the auroral precipitation and ionosphericconvection. These were then used as inputs to the Thermosphere–Ionosphere–ElectrodynamicsGeneral Circulation Model (TIEGCM) to simulate the thermospheric and ionospheric responseduring the storm period. The November 1993 storm was an unusually strong storm associatedwith a recurring high speed stream of solar plasma velocity in the declining phase of the solarcycle. Significant gravity waves with phase speeds of about 700 m/s caused by Joule heating werepresent in the upper thermosphere as perturbations to the neutral temperature and wind fields,especially on 4 November. The observed gravity waves in the meridional wind and in the height ofthe electron density peak at several southern hemisphere stations were generally reproduced bythe model using the AMIE high latitude inputs. Both model and observed equatorward windswere enhanced during the peak of the storm at Millstone Hill and at Australian ionosondestations. The observed neutral temperature at Millstone Hill increased about 400 K during thenight on 4 November, returning to normal on 9 November, while the model increased 300 K thefirst night at that location but was still elevated on 11 November. Enhanced westward windsduring the storm were evident in the UARS WIND Imaging Interferometer (WINDII) data. Theenhanced westward winds in the model were largest around 40–45° magnetic latitude at night,and also tended to be largest in the longitudes containing the magnetic poles. The peak westwardwind enhancements at 0 LT reached about 250 m/s at 300 km, and about 100 m/s at 125 km thefirst day of the storm at 40° magnetic latitude. At 20° magnetic latitude, the maximum westwardwind enhancements at 125 km at 0 LT appeared 2–4 days after the major part of the storm,indicating very long time constants in the lower thermosphere. The model showed global averageneutral temperature enhancements of 188 K after the peak of the storm that decayed with time,and which correlated with variations 8 h earlier in the Dst index and in the electric potential dropinput from AMIE. The global average temperature enhancement of 188 K corresponded to apotential drop increase of only about 105 kV. The results showed that the TIEGCM usingrealistic AMIE auroral forcings were able to reproduce many of the observed time dependentfeatures of this long-lived geomagnetic storm. The overall global average exospheric temperaturevariation correlated well with the time variation of the cross-tail potential drop and the Dst indexduring the storm period. However, the enhanced westward winds at mid-latitudes were stronglyrelated to the corrected Joule heating defined by the time dependent AMIE inputs.

Energetics of Magnetic Storms Driven by Corotating Interaction Regions: A Study of Geoeffectiveness

We investigate the energetics of magnetic storms associated with corotating interaction regions (CIRs). We analyze 24 storms driven by CIRs and compare to 18 driven by ejecta-related events to determine how they differ in overall properties and in particular in their distribution of energy. To compare these different types of events, we look at events with comparable input parameters such as the epsilon parameter and note the properties of the resulting storms. We estimate the energy output by looking at the ring current energy along with ionospheric Joule heating derived from the PC and Dst indices. We also include the energy of auroral precipitation, estimated from NOAA/TIROS and DMSP observations. In general, ejecta-driven storms produce more intense events, as parameterized by Dst*, but they are usually not as long lasting, and in most cases deposit less energy. This is observed even for events that have similar input quantities, such as epsilon. This may be related to the high speed of the solar wind, in that an increased magnetosonic Mach number may influence the reconnection rate and therefore the coupling. Additionally, we find the efficiency of the coupling varies greatly from CIR-driven to ejecta-driven storms, with the CIR-driven storms coupling substantially more efficiently, particularly in the recovery phase. The efficiency of coupling (output energy divided by input energy) for CIR-driven storms in recovery phase was double that of ejecta-driven storms.