A Perspective on Substorm Dynamics Using 10 Years of Auroral Kilometric Radiation Observations From Wind

We study 10 years (1995-2004 inclusive) of auroral kilometric radiation (AKR) radio emission data from the Wind spacecraft to examine the link between AKR and terrestrial substorms. We use substorm lists based on parameters including ground magnetometer signatures and geosynchronous particle injections as a basis for superposed epoch analyses of the AKR data. The results for each list show a similar, clear response of the AKR power around substorm onset. For nearly all event lists, the average response shows that the AKR power begins to increase around 20 min prior to expansion phase onset, as defined by the respective lists. The analysis of the spectral parameters of AKR bursts show that this increase in power is due to an extension of the source region to higher altitudes, which also precedes expansion phase onset by 20 min. Our observations show that the minimum frequency channel that observes AKR at this time, on average, is 60 kHz. AKR visibility is highly sensitive to observing spacecraft location, and the biggest radio response to substorm onset is seen in the 21:00-03:00 hr local time sector.

New Insights Into the Substorm Initiation Sequence From the Spatio-Temporal Development of Auroral Electrojets

In the present study we examine three substorm events, Events 1-3, focusing on the spatio-temporal development of auroral electrojets (AEJs) before auroral breakup. In Events 1 and 2, auroral breakup was preceded by the equatorward motion of an auroral form, and the ground magnetic field changed northward and southward in the west and east of the expected equatorward flow, respectively. Provided that these magnetic disturbances were caused by local ionospheric Hall currents, this feature suggests that the equatorward flow turned both eastward and westward as it reached the equatorward part of the auroral oval. The auroral breakup took place at the eastward-turning and westward-turning branches in Events 1 and 2, respectively, and after the auroral breakup, the westward AEJ enhanced only on the same side of the flow demarcation meridian. The zonal flow divergence is considered as an ionospheric manifestation of the braking of an earthward flow burst in the near-Earth plasma sheet and subsequent dawnward and duskward turning. Therefore, in Events 1 and 2, the auroral breakup presumably mapped to the dawnward and duskward flow branches, respectively. Moreover, for Event 3, we do not find any pre-onset auroral or magnetic features that can be associated with an equatorward flow. These findings suggest that the braking of a pre-onset earthward flow burst itself is not the direct cause of substorm onset, and therefore, the wedge current system that forms at substorm onset is distinct from the one that is considered to form as a consequence of the flow braking.

Katz Fractal Dimension of Geoelectric Field during Severe Geomagnetic Storms

We are concerned with the time series resulting from the computed local horizontal geoelectric field, obtained with the aid of a 1-D layered Earth model based on local geomagnetic field measurements, for the full solar magnetic cycle of 1996–2019, covering the two consecutive solar activity cycles 23 and 24. To our best knowledge, for the first time, the roughness of severe geomagnetic storms is considered by using a monofractal time series analysis of the Earth electric field. We show that during severe geomagnetic storms the Katz fractal dimension of the geoelectric field grows rapidly.

Time-scale dependence of solar wind-based regression models of ionospheric electrodynamics

The solar wind influence on geospace can be described as the sum of a directly driven component, or dayside reconnection, and an unloading component, associated with the release of magnetic energy via nightside reconnection. The two processes are poorly correlated on short time scales, but exactly equal when averaged over long time windows. Because of this peculiar property, regression models of ionospheric electrodynamics that are based on solar wind data are time scale specific: Models derived from 1 min resolution data will be different from models derived from hourly, daily, or monthly data. We explain and quantify this effect on simple linear regression models of various geomagnetic indices. We also derive a time scale-dependent correction factor that can be used with the Average Magnetic field and Polar current System model. Finally, we show how absolute estimates of the nightside reconnection rate can be calculated from solar wind measurements and geomagnetic indices.

Simultaneous observation of auroral substorm onset in Polar satellite global images and ground-based all-sky images

Substorm onset has originally been defined as a longitudinally extended sudden auroral brightening (Akasofu initial brightening: AIB) followed a few minutes later by an auroral poleward expansion in ground-based all-sky images (ASIs). In contrast, such clearly marked two-stage development has not been evident in satellite-based global images (GIs). Instead, substorm onsets have been identified as localized sudden brightenings that expand immediately poleward. To resolve these differences, optical substorm onset signatures in GIs and ASIs are compared in this study for a substorm that occurred on December 7, 1999. For this substorm, the Polar satellite ultraviolet global imager was operated with a fixed-filter (170 nm) mode, enabling a higher time resolution (37 s) than usual to resolve the possible two-stage development. These data were compared with 20-s resolution green-line (557.7 nm) ASIs at Muonio in Finland. The ASIs revealed the AIB at 2124:50 UT and the subsequent poleward expansion at 2127:50 UT, whereas the GIs revealed only an onset brightening that started at 2127:49 UT. Thus, the onset in the GIs was delayed relative to the AIB and in fact agreed with the poleward expansion in the ASIs. The fact that the AIB was not evident in the GIs may be attributed to the limited spatial resolution of GIs for thin auroral arc brightenings. The implications of these results for the definition of substorm onset are discussed herein.

Using ultra-low frequency waves and their characteristics to diagnose key physics of substorm onset

Substorm onset is marked in the ionosphere by the sudden brightening of an existing auroral arc or the creation of a new auroral arc. Also present is the formation of auroral beads, proposed to play a key role in the detonation of the substorm, as well as the development of the large-scale substorm current wedge (SCW), invoked to carry the current diversion. Both these phenomena, auroral beads and the SCW, have been intimately related to ultra-low frequency (ULF) waves of specific frequencies as observed by ground-based magnetometers. We present a case study of the absolute and relative timing of Pi1 and Pi2 ULF wave bands with regard to a small substorm expansion phase onset. We find that there is both a location and frequency dependence for the onset of ULF waves. A clear epicentre is observed in specific wave frequencies concurrent with the brightening of the substorm onset arc and the presence of "auroral beads". At higher and lower wave frequencies, different epicentre patterns are revealed, which we conclude demonstrate different characteristics of the onset process; at higher frequencies, this epicentre may demonstrate phase mixing, and at intermediate and lower frequencies these epicentres are characteristic of auroral beads and cold plasma approximation of the "Tamao travel time" from near-earth neutral line reconnection and formation of the SCW.

What effect do substorms have on the content of the radiation belts?

Substorms are fundamental and dynamic processes in the magnetosphere, converting captured solar wind magnetic energy into plasma energy. These substorms have been suggested to be a key driver of energetic electron enhancements in the outer radiation belts. Substorms inject a keV "seed" population into the inner magnetosphere which is subsequently energized through wave-particle interactions up to relativistic energies; however, the extent to which substorms enhance the radiation belts, either directly or indirectly, has never before been quantified. In this study, we examine increases and decreases in the total radiation belt electron content (TRBEC) following substorms and geomagnetically quiet intervals. Our results show that the radiation belts are inherently lossy, shown by a negative median change in TRBEC at all intervals following substorms and quiet intervals. However, there are up to 3 times as many increases in TRBEC following substorm intervals. There is a lag of 1-3 days between the substorm or quiet intervals and their greatest effect on radiation belt content, shown in the difference between the occurrence of increases and losses in TRBEC following substorms and quiet intervals, the mean change in TRBEC following substorms or quiet intervals, and the cross correlation between SuperMAG AL (SML) and TRBEC. However, there is a statistically significant effect on the occurrence of increases and decreases in TRBEC up to a lag of 6 days. Increases in radiation belt content show a significant correlation with SML and SYM-H, but decreases in the radiation belt show no apparent link with magnetospheric activity levels.

Shape-Constrained Sparse and Low-Rank Decomposition for Auroral Substorm Detection

An auroral substorm is an important geophysical phenomenon that reflects the interaction between the solar wind and the Earth's magnetosphere. Detecting substorms is of practical significance in order to prevent disruption to communication and global positioning systems. However, existing detection methods can be inaccurate or require time-consuming manual analysis and are therefore impractical for large-scale data sets. In this paper, we propose an automatic auroral substorm detection method based on a shape-constrained sparse and low-rank decomposition (SCSLD) framework. Our method automatically detects real substorm onsets in large-scale aurora sequences, which overcomes the limitations of manual detection. To reduce noise interference inherent in current SLD methods, we introduce a shape constraint to force the noise to be assigned to the low-rank part (stationary background), thus ensuring the accuracy of the sparse part (moving object) and improving the performance. Experiments conducted on aurora sequences in solar cycle 23 (1996-2008) show that the proposed SCSLD method achieves good performance for motion analysis of aurora sequences. Moreover, the obtained results are highly consistent with manual analysis, suggesting that the proposed automatic method is useful and effective in practice.

Quantitative maps of geomagnetic perturbation vectors during substorm onset and recovery

We have produced the first series of spherical harmonic, numerical maps of the time-dependent surface perturbations in the Earth's magnetic field following the onset of substorms. Data from 124 ground magnetometer stations in the Northern Hemisphere at geomagnetic latitudes above 33° were used. Ground station data averaged over 5 min intervals covering 8 years (1998-2005) were used to construct pseudo auroral upper, auroral lower, and auroral electrojet (AU*, AL*, and AE*) indices. These indices were used to generate a list of substorms that extended from 1998 to 2005, through a combination of automated processing and visual checks. Events were sorted by interplanetary magnetic field (IMF) orientation (at the Advanced Composition Explorer (ACE) satellite), dipole tilt angle, and substorm magnitude. Within each category, the events were aligned on substorm onset. A spherical cap harmonic analysis was used to obtain a least error fit of the substorm disturbance patterns at 5 min intervals up to 90 min after onset. The fits obtained at onset time were subtracted from all subsequent fits, for each group of substorm events. Maps of the three vector components of the averaged magnetic perturbations were constructed to show the effects of substorm currents. These maps are produced for several specific ranges of values for the peak |AL*| index, IMF orientation, and dipole tilt angle. We demonstrate an influence of the dipole tilt angle on the response to substorms. Our results indicate that there are downward currents poleward and upward currents just equatorward of the peak in the substorms' westward electrojet.

KEY POINTS

Show quantitative maps of ground geomagnetic perturbations due to substorms Three vector components mapped as function of time during onset and recovery Compare/contrast results for different tilt angle and sign of IMF Y-component.