Using ultra-low frequency waves and their characteristics to diagnose key physics of substorm onset
© The Author(s) 2017
Received: 3 January 2017
Accepted: 11 October 2017
Published: 24 October 2017
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.
Substorm expansion phase onset is traditionally identified in ionospheric optical observations by the brightening of a pre-existing auroral arc or the formation of a new arc (Akasofu 1964, 1977). This brightening is accompanied by ultra-low frequency (ULF) waves in the Pi2 band (40–150 s period; Jacobs et al. 1964) that develop when the cross-tail current is diverted towards the ionosphere during the dipolarisation of the nightside magnetospheric field and formation of the substorm current wedge (SCW; e.g. Atkinson 1967; McPherron et al. 1973). Following substorm onset, these Pi2s have also been shown to couple to the plasmasphere forming a plasmaspheric cavity mode (Yeoman and Orr 1989) and have been linked to plasma sheet flows in the tail (Kepko and Kivelson 1999; Murphy et al. 2011b), presumably corresponding to transient magnetotail reconnection (Keiling et al. 2006). At substorm onset, the general characteristics of Pi2s are used to determine the large-scale characteristics of the SCW, including the locations of the upward and downward field-aligned current (FAC) elements and the central meridian (e.g. Lester et al. 1983, 1984). However, the precise relationship between Pi2 waves and the first observable changes of auroral arc brightening is limited by a number of observational constraints. Pi2 waves have periods of ~ 40–150 s, which are of the order of the “2-min problem” (Ohtani 2004) surrounding onset. Additionally, the integrated ionospheric currents that magnetometers measure are inherently noisy and so it is challenging to identify the “first” evidence of ULF wave activity when the low-amplitude signals are growing out of a noisy or elevated background. Linking Pi2 to substorm onset is further complicated by the plethora of physical phenomena that have been suggested to produce ULF pulsations in the ionosphere. Shorter period ULF wave bands such as the Pi1B (~ 1–10 s) period band, part of the overall Pi1 band (1–40 s), offer a means to reduce the onset timing uncertainty (e.g. Bösinger 1989; Lessard et al. 2006; Posch et al. 2007) in a localised region close to onset (e.g. Bösinger and Yahnin 1987; Arnoldy et al. 1987) and more importantly are not directly linked to other magnetospheric phenomena.
The long-period part of Pi1/short-period part of Pi2 ULF wave band (hereafter referred to as the long-period Pi1 band or Pi1-2 for brevity) was first investigated by Milling et al. (2008). In this paper, a discrete wavelet transform (DWT) based upon the Pi2 algorithm outlined by Nose et al. (1998) was used to investigate the entire ULF wave spectrum during substorm expansion phase onset. The Meyer wavelet was used for the analysis as it closely resembles the impulsive nature of nightside ULF waves and has excellent timing resolution, and hence is ideal for defining any “onset”. Milling et al. (2008) found that the first ULF wave band to rise above a pre-determined noise threshold was in the long-period Pi1 wavelet band, and when the analysis was extended across all available magnetometers, the long-period Pi1 waves had a clear and coherent onset that spread out from an epicentre in the auroral zone.
Murphy et al. (2009a) presented an automated algorithm based on the Meyer wavelet (christened the automated wavelet estimation of substorm onset and magnetic events—AWESOME) and used it to compare the onset timing and location of long-period Pi1 waves with global auroral intensification times detailed in the Frey substorm listings (Frey et al. 2004; Frey and Mende 2006). Murphy et al. found that the epicentre of ULF wave activity was co-located with the initial location of global auroral intensification observed at substorm expansion phase onset, as determined by the global auroral imaging from the IMAGE-FUV instrument. However, the onset of long-period Pi1 ULF wave activity generally occurred ~ minutes prior to the global auroral intensification during substorm expansion phase onset. Liou and Zhang (2009) argued that the temporal resolution of the IMAGE-FUV instrument meant that the ULF wave onset and substorm onset could occur simultaneously. However, Murphy et al. (2009b) demonstrated that even with a more in-depth consideration of the temporal cadence of each instrument in the study, in general ULF wave onset always precedes global auroral breakup.
Rae et al. (2009a) studied the relationship between smaller scale auroral features observed minutes prior to auroral breakup and the onset of Pi1-2 ULF waves observed by both Milling et al. (2008) and Murphy et al. (2009a, b). After identifying the onset arc, Akasofu (1964) and Rae et al. (2009a, b) showed that the onset of Pi1-2 ULF waves occurred at the same time and in the same location as “auroral beads” (Henderson 1994) which developed azimuthally along the onset arc. Rae et al. (2009b) and Walsh et al. (2010) went on to demonstrate the rapid connectivity of the equatorial magnetosphere to the ionosphere during substorm onset, showing that the Pi1-2 ULF onset observed on the ground and at the conjugate point in the equatorial magnetosphere occurred within the timing uncertainty defined by the AWESOME algorithm. Further studies of the onset of the ULF waves shows that the exponential growth in ULF wave amplitudes occurs immediately prior to the onset of an auroral substorm (Rae et al. 2011). Finally, it has been shown that both the ULF wave amplitudes and auroral brightness grow exponentially at the location of ionospheric substorm onset (Rae et al. 2012; Kalmoni et al. 2015, 2017). Although the relationship between the enhancement of ULF waves and auroral brightness is no longer in doubt (see review by Rae and Watt 2016), the relative timing and localisation of Pi1, Pi1-2 and Pi2 waves across the auroral zone at substorm onset have not yet been investigated.
The arrival time of Pi2 ULF waves has been used by Chi et al. (2009) to infer a magnetospheric source of substorm onset in the magnetotail by investigating the latitudinal dependence of the first peak in Pi2 amplitude observed by ground magnetometers along a meridian. In their paper, Chi et al. (2009) found that there is a strong latitudinal dependence of this first Pi2 peak with stations in the auroral zone observing this peak ~ 1 min prior to both higher and lower latitudinal stations. By considering the travel path preserving the most wave energy using the Tamao travel time (e.g. Tamao 1964), these authors concluded that the distances and timing that were most consistent with the latitudinal profile of the first Pi2 peak were ~ 1–3 min prior at 15–25 RE distances in the magnetotail. This technique is often referred to as “magnetoseismology”, and for a comprehensive review of this, we direct the reader to Menk and Waters (2013). Here, we study not the first notable peak in wave amplitude, but the onset of wave activity as defined by the AWESOME algorithm (Murphy et al. 2009a), and we extend the analysis from Pi1 waves, through the long-period Pi1 band, to Pi2 waves to study all frequencies (see Murphy et al. 2011a).
Hence, as discussed above, there have been two specific drivers of two specific wave bands discussed in the literature; Pi1/2 waves are in general associated with the substorm onset arc that maps in the near-Earth region between stretched and dipolar fields where current disruption plays a key role (e.g. Roux et al. 1991), whereas Pi2 waves are related to near-Earth neutral line physics driven by magnetotail reconnection (e.g. Hones 1976). Determining differences between these two ULF wave bands therefore potentially represents a novel method to distinguish between these two substorm onset mechanisms.
To understand all aspects of ULF wave studies, one must study the timing, latitudinal, longitudinal and frequency-dependent onset characteristics of ULF waves and, to date, no such measurement of all four quantities has been made. In this paper, we study the frequency-dependent onset of ULF waves as a function of time and location. Specifically, we present magnetic observations of a small substorm on 5 March 2008, where we utilise AWESOME to define the 2D ionospheric projection of the first rapid increase in magnetic oscillation amplitude in the Pi1, Pi1-2 and Pi2 ULF wave bands close to the auroral onset. Using this technique, we directly compare both the absolute and relative timing and location characteristics of all three essential ULF wave bands during substorm expansion phase onset. We find that the Pi1-2 ULF wave band occurs first at auroral latitudes, but Pi2 ULF waves may be observed concurrently at lower latitudes, to within the timing uncertainty of the technique. We postulate that Pi1-2 ULF waves are critically important for defining the timing and location of the auroral substorm expansion phase, but that the low-latitude Pi2 onset may provide valuable information on the location of the plasmapause. We discuss these results in terms of previous findings.
Instrumentation and large-scale morphology
AWESOME magnetic timing
The discrete wavelet transform (DWT) utilised by AWESOME uses a Meyer basis wavelet to decompose a non-stationary signal into discrete wavelet coefficients that are localised in time and band limited in frequency in order to objectively define the first time ULF wave amplitudes rise above a noise threshold, which we define as the onset of ULF waves, during magnetic substorms (see Murphy et al. 2009a for details). For each individual magnetometer station, the quietest interval preceding the event during the UT day is selected (i.e. an interval that does not include significant ULF wave variations) and the power in each wavelet period band is calculated during that interval. A threshold is then defined independently for each wavelet band as the mean plus two standard deviations of the power. Onset in a specific ULF wave band is then defined as the initial time at which the wavelet power coefficient rises above the quite time threshold. Our definition ensures a 98% confidence level that the signal which rises above this threshold represents a statistically significant signal above quiet time fluctuations. A discrete onset time is defined as the centre of the interval when the wavelet coefficient first exceeds this threshold, and the uncertainty is defined as plus or minus the temporal width of the coefficient band. Statistically, the ULF power spectra around onset are characteristic of a power law; neither the Pi1, Pi1-2 or Pi2 frequency bands are preferred at substorm onset nor is there any clear break in the ULF wave spectra (Murphy et al. 2011b). The AWESOME algorithm is designed so that the entire ULF wave spectrum is analysed without identifying an a priori ULF wave band for study during any particular event. In this paper, the AWESOME algorithm is utilised to detail the absolute and relative timing in the 12–48 s, 24–96 s and 48–192 s wavelet bands, timing uncertainties for each of the three wavelet bands being ± 8 s, ± 16 s and ± 32 s, respectively.
A case study in relative ULF wave timing
Figure 4a shows the minimum curvature fit to the Pi2 arrival times, relative to the first Pi2 arrival time at RABB at 06:05:08 UT ± 32 s, whereas Fig. 4b shows timed relative to the Pi1-2 arrival time at GILL-RABB at 06:04:20 ± 16 s, and Fig. 4c shows timed relative to the Pi1 arrival time at GILL at 06:06:20 ± 8 s. There are notable similarities and differences between the onset of ULF waves in each of the three Pi1 and Pi2 wave bands, but for each wave band at least one location of the epicentre is well aligned with the onset arc as shown in Fig. 1. In the case of the Pi1/2 wave band, the epicentre is indeed well aligned with the east–west onset arc. We note here that the closed contours at RABB and GILL are part of the same epicentre, as each magnetometer observes the onset signature at the same time. Hence, these closed contours are simply a result of the minimum curvature fit to the data, rather than any physical process. This alignment of the Pi1/2 epicentre is consistent with the results of Murphy et al. (2009a, b), Rae et al. (2009a, b, 2010, 2011) and Walsh et al. (2010), amongst others, that the long-period Pi1/short-period Pi2 waves are initially seen in the ionosphere at the location of substorm onset. Interestingly, however, the times at which these epicentres are initiated are significantly different. Using AWESOME, we can account for the uncertainty in timing, and we demonstrate that in this case study the Pi1-2 wave band is observed ~ 1 min prior to the Pi2 wave band and ~ 2 min prior to the Pi1 wave band. Rae et al. (2009a, b) showed that the onset of Pi1-2 wave powers was commensurate with the development of auroral beads and postulated that the onset of Pi2 waves was associated with auroral breakup. Again, this is consistent with the results presented in this paper.
However, the Pi1, Pi1-2 and Pi2 arrival patterns are quite dissimilar across a wide range of latitudes and longitudes. Indeed, there is a second Pi2 epicentre observed at low latitudes centred around OSAK, where Pi2 onset is observed at the same time as the auroral Pi2 onset location at RABB. The onsets of Pi2 waves at stations between OSAK and RABB (THRF and BRDN) are later; there are sufficient data in this part of the plot that we may have confidence in the contours as drawn. We discuss the implications of this Pi2 onset distribution in the “Discussion” section.
The patterns of onset of the Pi1 waves shown in Fig. 4c differ from both the Pi1-2 and Pi2 ULF wave onsets. Firstly, the AWESOME-determined Pi1 onset does not extend to the lower latitudes that the Pi1-2 and Pi2 onsets did in this case study; the lowest latitude station that observed Pi1 onset was PINE at ~ 52° latitude. Secondly, although the high-latitude Pi1 onset exhibits largely the same spatial onset pattern as the Pi1-2 wave onset, both occurring first at GILL and expanding clearly and coherently away from this epicentre region, there is a delay of ~ 2 min after the Pi1-2 onset.
In summary, it can be observed from Fig. 4 that Pi1-2 waves are initially confined along the substorm onset arc that is azimuthally aligned with the GILL-RABB-FSMI magnetometers, and the onset signature propagates relatively coherently across the rest of the magnetometer array, which is in keeping with the results shown previously by Murphy et al. (2009a), Rae et al. (2009a, b) and Walsh et al. (2010). Pi2 waves occur later in the substorm and their onset signals propagate significantly faster but still coherently across the array from two epicentres, one centred on the location of the substorm onset arc and another at lower latitudes. Finally, Pi1 waves occur even later in the substorm, are confined near onset and propagate relatively slowly compared to other ULF wave signals.
Importance of relative timing of wave onsets in substorm physics
In “A case study in relative ULF wave timing” section, we demonstrate that each ULF wave band has significantly different onset times as a function of radial distance away from the auroral onset region. As in numerous previous studies, the epicentre of Pi1-2 ULF wave onset in this case study occurs contemporaneously with the growing auroral bead signature that designates auroral substorm onset, in this case study near the GILL station. These auroral beads, which mark the onset of the substorm in the ionosphere (e.g. Rae et al. 2010; Murphy et al. 2015; Kalmoni et al. 2015, 2017), have been shown to grow exponentially in brightness at the same time as the onset of Pi1-2 ULF waves (~ 0.05 s−1; Kalmoni et al. 2015, demonstrating that auroral onset and ULF wave activity are inextricably linked through whichever process(es) cause substorm onset, e.g. Murphy et al. (2009a, b, 2015) and Rae et al. (2009a, b, 2010, 2011, 2012).
This Pi1-2 onset occurs first, but is closely followed by the onset of Pi2 ULF wave activity at two separated epicentres close to the RABB-GILL station and at mid-latitudes at OSAK. The Pi2 waves are observed around 1 min later at the onset station at GILL. Finally, the shorter period Pi1 waves are observed in a region around the GILL-RABB stations.
The primary aim of developing the AWESOME technique was to analyse the entire ULF wave spectra and provide an estimate of uncertainty in the onset time, in addition to setting an objective and quantitative onset criteria. We note here that there is no overlap in uncertainty between the onsets of the different ULF wave bands, demonstrating that the onset of Pi1-2 waves occurs prior to Pi2 waves, and both occur prior to the higher frequency Pi1 wave band consistent with previous studies at single stations (Murphy et al. 2009a, b; Rae et al. 2009a, b; Walsh et al. 2010). We conclude that the Pi1-2 waves are linked to the formation and evolution of auroral beads and are therefore onset driven (e.g. Rae et al. 2009a, b), whereas the longer period Pi2 waves may be related to the evolution of the large-scale current system that develops following onset, most likely the substorm current wedge (e.g. McPherron et al. 1973; Rae et al. 2007; Murphy et al. 2013). The timing difference between auroral beads and auroral breakup during which the SCW begins to form is similar to timing difference between the onset of the Pi1-2 and Pi2 wave bands; beads and Pi1/2 waves occurring contemporaneously and auroral breakup and Pi2 waves occurring 2–3 min later. This further supports our conclusion that Pi1-2 waves are linked to the auroral beads and Pi2 waves are linked to the SCW. Presumably in this scenario, Pi1 waves correspond to the later structuring and filamentation of the current systems that develop into auroral breakup (e.g. Arnoldy et al. 1987; Bösinger 1989).
We further analyse the onset times of these three ULF wave bands in order to highlight the precise relative timing between the onset times as a function of latitude and longitude and reveal additional information on the formation and evolution of the substorm.
A number of features stand out in Fig. 5a, most notably that at the auroral and ULF wave onset station Pi1-2 ULF waves clearly lead Pi2 waves by ~ 1 min and across the majority of latitudes and longitudes the Pi1-2 onset leads the Pi2 onset by ~ 1–2 min. However, there appears to be a channel in the ionosphere where Pi2 ULF waves are observed prior to Pi1-2 waves, ranging from TPAS to OSAK and to the west of substorm onset, marked by the blue timing contours in Fig. 5a between 320–345° longitudes.
Figure 5b shows the relative timing between the Pi1-2 and Pi1 ULF wave bands. Figure 5b clearly demonstrates that the onset of Pi1-2 ULF waves occur ~ 1–4 min prior to the higher frequency Pi1 ULF waves, other than at one station (BRDN) which lies within the timing uncertainty contour. Note that the relative timing cannot be calculated below 55° latitude since the Pi1 onsets were only observed over a limited spatial region. The onset of Pi1 ULF waves is delayed longest with respect to the onset of Pi1-2 waves at these lowest observable latitudes.
We present a case study of the relative timing of ULF wave onset across the Pi1 and Pi2 ULF wave bands. These waves are known to play a pivotal role in substorm studies as a diagnostic (Rae et al. 2009a, b, 2010) and a timing mechanism for substorm onset (Milling et al. 2008; Murphy et al. 2009a, b). We find that there are clear epicentres observed in each wave frequency band, and we discuss the ramifications on the generation mechanisms of each of these specific ULF wave bands.
Historically, Pi2 waves have been invoked to carry the current of the substorm current wedge (e.g. McPherron et al. 1973; Rae et al. 2007), related to the onset of NENL reconnection (e.g. Liou et al. 2000; Uozumi et al. 2000, 2004, 2009; Chi et al. 2009). Using magnetoseismology or “Tamao travel times” in a cold plasma, Chi et al. (2009) determined that, by studying a visually determined Pi2 onset time in the ionosphere, constraints on the location of NENL reconnection could be determined. Interestingly, Chi et al. (2009) determined that there were two paths which minimised the time-of-flight of a Pi2: one along the PSBL and one at lower latitudes outside the plasmapause where Alfven speeds are higher than the surrounding plasma sheet. This would lead to the onset of Pi2 waves at two discrete epicentres, in precisely the same configuration as shown in Fig. 4(left). Pi2 waves also appear to be global, being observed throughout the nightside ionosphere, which is consistent with the historical literature suggesting that Pi2 waves are a “global mode” (see Keiling and Takahashi 2011 for a comprehensive review). In this scenario, Pi2 waves at higher latitudes would correspond to the open-closed separatrix shown in Fig. 2b to be situated around 70° latitude (e.g. Blanchard et al. 1995). As shown in Fig. 4(left), this location is not populated with a ground magnetometer, but the closest high-latitude magnetometer station (RABB at 67°) is the first magnetometer to observe the first Pi2 signals, which supports this hypothesis. Moreover, depending upon how stretched the magnetic field is and where NENL reconnection is initiated, the onset of Pi2 waves would potentially be simultaneous or later for the magnetometer station closest to the plasmapause (Chi and Russell 2005; Chi et al. 2009). As travel time magnetoseismology is a strong function of plasma sheet temperature, field line stretching and assumed NENL reconnection location, modelling wave propagation times is beyond the scope of this paper. However, using the results from the literature described above, we suggest that the two epicentres of Pi2 wave activity discovered in this case study, coupled with a meridional channel of early Pi2 arrival times, support the hypothesis of NENL reconnection occurring outside of 15RE in the magnetotail (Fig. 2 of Chi et al. 2009). In this scenario, Pi2 waves are observed either in two epicentres described above, or relatively simultaneously in a single meridian (see Figs. 4a, 5a).
On the other hand, Pi1-2 waves have been exclusively linked to the formation and structuring of the substorm onset arc (e.g. Rae et al. 2009a, b). Specifically, Pi1-2 waves have been linked to the formation of what is termed “auroral beads” after Henderson (1994). These auroral fluctuations have the same temporal periodicity as Pi1-2 ULF waves and initially form azimuthally along the onset arc (Rae et al. 2012). Evident from Fig. 4b is the fact that the Pi1-2 onset occurs initially at GILL, RABB and FSMI stations conjugate to the onset arc. This supports the hypothesis that Pi1-2 waves are the magnetic counterparts of auroral beads. Detailed analysis of these auroral beads (Rae et al. 2010; Kalmoni et al. 2015, 2017) and auroral and plasma sheet dynamics at substorm onset (Rae et al. 2009a; Murphy et al. 2015) has provided strong evidence that auroral beads are the ionospheric manifestation of a plasma instability driven by substorm onset, whether it be in the near-Earth magnetotail (e.g. Lui 1991) or closer to the ionosphere (e.g. Motoba et al. 2012). Since both Pi1-2 ULF waves (e.g. Murphy et al. 2011b; Rae et al. 2011) and auroral beads (Kalmoni et al. 2017) are ubiquitous to substorm onset, the Pi1-2 waves appear to be exclusively related to the substorm onset instability at the inner edge of the plasma sheet (e.g. Samson et al. 1992).
Finally, Pi1 waves occur in minutes following substorm onset during the substorm expansion phase and auroral breakup as the aurora expands poleward and toward the dawn and dusk of the onset region. In a statistical study of 11 years of DMSP electron observations, Wing et al. (2013) showed that following substorm onset there is a sharp increase in wave-driven aurora characterised by broadband electron precipitation (c.f., Newell et al. 2010). Our observations together with those from Wing et al. (2013) indicate that Pi1 waves are related to the filamentary currents generated following substorm onset and localised poleward of the onset arc region (e.g. Arnoldy et al. 1987; Bösinger 1989).
That two separate mechanisms are invoked during substorm onset is not a new concept (e.g. Lui 1991; Pu et al. 1999). Studies often appeal to either one mechanism occurring before the other, such as reconnection-driven current disruption (e.g. Angelopoulos 2008) or via current disruption leading to reconnection subsequently (e.g. Lui 2009). Indeed, these could be independent regions being activated independently (e.g. Murphy et al. 2015). Our results are consistent with the literature that Pi2 waves are driven by reconnection and Pi1-2 waves are driven by current disruption, which would be observational evidence of the theory that current disruption occurs first and mid-tail reconnection occurs subsequently (e.g. Lui 1991; Pu et al. 1999). This is also consistent with the results of Rae et al. (2011) who performed a statistical study of the onset times of Pi1-2 and Pi2 waves at the location of substorm auroral brightening and found that Pi1-2 waves were observed consistently before the signatures in Pi2 waves.
In this paper, we probe the relationship between three distinct ULF wave period bands and the initiation of a small substorm expansion phase onset using the AWESOME algorithm (Murphy et al. 2009a). We determine the relative timing between sudden increases in amplitude, or onsets, of different ULF wave bands at multiple ground-based magnetometer stations, and discuss how two separate theories of substorm onset dynamics can explain the wave characteristics observed.
We find that the onset of long-period Pi1/short-period Pi2 (Pi1-2) ULF waves are the first ULF wave signature observed at substorm onset followed by Pi2s and subsequently by Pi1s. The onset of the Pi1-2 waves is co-located in time and space with the formation of small-scale auroral beads. This is in contrast to the later arrival of Pi2 waves that occur in two distinct epicentres, one at auroral latitudes and one at lower latitudes. The onset of Pi1 waves occurs following both Pi1-2 and Pi2 onsets and is confined to latitudes near substorm onset. To date, there are no realistic travel times calculated for substorm onset, but the most realistic of any travel times is by Chi et al. (2009) using cold plasma theory using the onset time of Pi2 waves in the ionosphere as the marker for “Tamao travel times”. Indeed, Ferdouis and Raeder (2016) stated that the real MHD travel paths are more complicated than the Tamao travel path, but that there is good qualitative agreement between their studies and the predicted Tamao travel times.
We conclude that the differing onset times and spatial expansion of the Pi1, Pi1-2 and Pi2 waves are the result of the physical process which drives each wave in the ionosphere. Pi1-2 waves, being observed at onset and collocated with auroral beads, are the result of the initiation of a plasma instability in the inner magnetosphere which subsequently leads to the development of the substorm expansion phase. Pi2 waves, observed following Pi1-2 waves during the expansion phase, are the result of the development of the SCW and dipolarisation of the tail. Finally, Pi1 waves, occurring following both Pi1-2 and Pi2 waves during the expansion phase and localised to higher latitudes, are the result of filamentary currents and wave aurora developing during expansion phase as the aurora expands poleward and toward the dawn and dusk. These results demonstrate how careful analysis of ULF waves surrounding substorm onset can provide vital information on the physical processes occurring and time history of these processes through substorm onset.
Work on this topic has spanned several years, both in the development of the AWESOME wavelet algorithm and the scientific data exploitation. IJR, KRM, CEJW, IRM and DKM all provided substantive contributions in both areas. NMEK, ZY and CF provided detailed scientific discussions on the resultant analysis. All authors have contributed to the writing and preparation of this manuscript. All authors read and approved the final manuscript.
We thank all optical imagery and ground magnetometer providers. The authors thank NRCan for the provision of data from the Canadian magnetic observatory system (CANMOS) Network and CARISMA is operated by the University of Alberta, funded by the CSA. We acknowledge NASA contract NAS5-02099 and V. Angelopoulos for use of data from the THEMIS Mission and NSF support of GIMNAST through grant AGS-1004736. Specifically, we thank S. Mende and E. Donovan for use of the THEMIS ASI data. Deployment and data retrieval of the THEMIS ASIs were partly supported by CSA contract 9 F007-046101 to the University of Calgary. We thank S. Mende and C.T. Russell for the use of the THEMIS ground-based magnetometer data and NSF for support through grant AGS-1004814. Operational support for NORSTAR is provided by the CSA and University of Calgary, and we thank Brian Jackel and the NORSTAR team for providing the MSP data used in this study. Data reported in this study are available via http://themis.ssl.berkeley.edu and http://www.carisma.ca.
The authors declare that they have no competing interests.
Availability of data and materials
Timing information is available on request from the author. CARISMA data can be obtained from http://www.carisma.ca.
Consent for publication
Ethics approval and consent to participate
IJR is funded in part by Science and Technology Facilities (STFC) Grant ST/N000749/1 and Natural Environmental Research Council (NERC) Grant NE/L007495/1. CEJW is funded in part by STFC Grant ST/M000885/1. ZY is a Marie-Curie COFUND postdoctoral fellow at the University of Liege, co-funded by the European Union. CF is supported by a NERC Independent Research Fellowship NE/N014480/1.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
- Akasofu SI (1964) The development of the auroral substorm. Planet Space Sci 12(4):273–282. doi:10.1016/0032-0633(64)90151-5 View ArticleGoogle Scholar
- Akasofu SI (1977) Magnetospheric substorms. Q J R Astron Soc 18(2):170–187Google Scholar
- Angelopoulos V (2008) The THEMIS mission. Space Sci Rev 141:5–34. doi:10.1007/s11214-008-9336-1 View ArticleGoogle Scholar
- Arnoldy RL et al (1987) Simultaneous measurement of aurora-related, irregular magnetic pulsations at northern and southern high latitudes. J Geophys Res 92(A11):12221–12232. doi:10.1029/JA092iA11p12221 View ArticleGoogle Scholar
- Atkinson G (1967) Current system of geomagnetic bays. J Geophys Res 72(23):6063–6067. doi:10.1029/JZ072i023p06063 View ArticleGoogle Scholar
- Blanchard GT, Lyons LR, Samson JC, Rich FJ (1995) Locating the polar-cap boundary from observations of 6300 angstrom auroral emission. J Geophys Res 100:7855–7862View ArticleGoogle Scholar
- Bösinger T (1989) On the spectral index of the Pi1B power spectrum. Ann Geophys 7(4):375–386Google Scholar
- Bösinger T, Yahnin AG (1987) PI1B type magnetic pulsation as a high time resolution monitor of substorm development. Ann Geophys Ser A 5(4):231–237Google Scholar
- Chi PJ, Russell CT (2005) Travel-time magnetoseismology: magnetospheric sounding by timing the tremors in space. Geophys Res Lett 32:L18108. doi:10.1029/2005GL023441 View ArticleGoogle Scholar
- Chi PJ, Russell CT, Ohtani S (2009) Substorm onset timing via traveltime magnetoseismology. Geophys Res Lett 36:L08107. doi:10.1029/2008GL036574 View ArticleGoogle Scholar
- Donovan EF, Trondsen TS, Cogger LL, Jackel BJ (2003) All-sky imaging within the Canadian CANOPUS and NORSTAR projects. Sodank Geophys Obs Publ 92:109–112Google Scholar
- Ferdousi B, Raeder J (2016) Signal propagation time from the magnetotail to the ionosphere: OpenGGCM simulation. J Geophys Res Space Phys 121:6549–6561. doi:10.1002/2016JA022445 View ArticleGoogle Scholar
- Frey HU, Mende SB (2006) Substorm onsets as observed by IMAGE-FUV. In: Syrjäsuo M, Donovan EF (eds) Proceedings of ICS-8, Univ. of Calgary, Canada, pp 71–75Google Scholar
- Frey HU, Mende SB, Angelopoulos V, Donovan EF (2004) Substorm onset observations by IMAGE-FUV. J Geophys Res 109:A10304. doi:10.1029/2004JA010607 View ArticleGoogle Scholar
- Henderson MG (1994) Implications of viking imager results for substorm models. Univ. of Calgary, CalgaryGoogle Scholar
- Hones EW Jr (1976) The magnetotail: its generation and dissipation. In: Williams DJ (ed) Physics of solar planetary environments: proceedings of the international symposium on solar-terrestrial physics, vol 2. AGU, Washington, D. C., pp 558–571Google Scholar
- Jacobs JA, Kato Y, Matsushita S, Troitskaya VA (1964) Classification of geomagnetic micropulsations. J Geophys Res 69:180–181View ArticleGoogle Scholar
- Kalmoni NME, Rae IJ, Watt CEJ, Murphy KR, Forsyth C, Owen CJ (2015) Statistical characterization of the growth and spatial scales of the substorm onset arc. J Geophys Res Space Phys 120:8503–8516. doi:10.1002/2015JA021470 View ArticleGoogle Scholar
- Kalmoni NME, Rae IJ, Murphy KR, Forsyth C, Watt CEJ, Owen CJ (2017) Statistical azimuthal structuring of the substorm onset arc: implications for the onset mechanism. Geophys Res Lett 44:2078–2087. doi:10.1002/2016GL071826
- Keiling A, Takahashi K (2011) Review of Pi2 models. Space Sci Rev 161(1–4):63–148. doi:10.1007/s11214-011-9818-4 View ArticleGoogle Scholar
- Keiling A, Fujimoto M, Hasegawa H, Honary F, Sergeev V, Semenov VS, Frey HU, Amm O, Reme H, Dandouras I, Lucek E (2006) Association of Pi2 pulsations and pulsed re- connection: ground and cluster observations in the tail lobe at 16RE. Ann Geophys 24:3433–3449. doi:10.5194/angeo-24-3433-2006 View ArticleGoogle Scholar
- Kepko L, Kivelson MG (1999) Generation of Pi2 pulsations by bursty bulk flows. J Geophys Res 104:25021–25034View ArticleGoogle Scholar
- Lessard MR et al (2006) Nature of Pi1B pulsations as inferred from ground and satellite observations. Geophys Res Lett 33:L14108. doi:10.1029/2006GL026411 View ArticleGoogle Scholar
- Lester M et al (1983) Polarization patterns of Pi 2 magnetic pulsations and the substorm current wedge. J Geophys Res 88(A10):7958–7966. doi:10.1029/JA088iA10p07958 View ArticleGoogle Scholar
- Lester M, Hughes WJ, Singer HJ (1984) Longitudinal structure in Pi 2 pulsations and the substorm current wedge. J Geophys Res 89(A7):5489–5494. doi:10.1029/JA089iA07p05489 View ArticleGoogle Scholar
- Liou K, Zhang Y-L (2009) Comment on wavelet-based ULF wave diagnosis of substorm expansion phase onset by K. Murphy et al. J Geophys Res. doi:10.1029/2009JA014207 Google Scholar
- Liou K, Meng C-I, Newell PT, Takahashi K, Ohtani S-I, Lui ATY, Brittnacher M, Parks G (2000) Evaluation of low-latitude Pi2 pulsations as indicators of substorm onset using polar ultraviolet imagery. J Geophys Res 105(A2):2495–2505. doi:10.1029/1999JA900416 View ArticleGoogle Scholar
- Lui ATY (1991) A synthesis of magnetospheric substorm models. J Geophys Res 96(A2):1849–1856. doi:10.1029/90JA02430 View ArticleGoogle Scholar
- Lui ATY (2009) Comment on “Tail reconnection triggering substorm onset”. Science 324(5933):1391. doi:10.1126/science.1167726
- Mann IR et al (2008) The upgraded CARISMA magnetometer array in the THEMIS era. Space Sci Rev 141:413–451. doi:10.1007/s11214-008-9457-6 View ArticleGoogle Scholar
- McPherron RL et al (1973) Satellite studies of magnetospheric substorms on August 15, 1968. 9. Phenomenological model for substorms. Eos Trans AGU 54(4):206Google Scholar
- Mende SB, Harris SE, Frey HU, Angelopoulos V, Russell CT, Donovan E, Jackel B, Greffen M, Peticolas LM (2008) The THEMIS array of ground-based observatories for the study of auroral substorms. Space Sci Rev 141:357–387. doi:10.1007/s11214-008-9380-x View ArticleGoogle Scholar
- Menk FW, Waters CL (2013) Introduction, in magnetoseismology: ground-based remote sensing of earth’s magnetosphere. Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim. doi:10.1002/9783527652051.ch1 View ArticleGoogle Scholar
- Milling DK et al (2008) Ionospheric localization and expansion of long-period Pi1 pulsations at substorm onset. Geophys Res Lett. doi:10.1029/2008GL033672 Google Scholar
- Motoba T, Hosokawa K, Kadokura A, Sato N (2012) Magnetic conjugacy of northern and southern auroral beads. Geophys Res Lett 39:L08108. doi:10.1029/2012GL051599 View ArticleGoogle Scholar
- Murphy KR, Rae IJ, Mann IR, Milling DK, Watt CEJ, Ozeke L, Frey HU, Angelopoulos V, Russell CT (2009a) Wavelet-based ULF wave diagnosis of substorm expansion phase onset. J Geophys Res. doi:10.1029/2008JA013548 Google Scholar
- Murphy KR, Rae IJ, Mann IR, Milling DK, Watt CE, Ozeke L, Frey HU, Angelopoulos V, Russell CT (2009b) Wavelet-based ULF wave diagnosis of substorm expansion phase onset. J Geophys Res 114:A10Google Scholar
- Murphy KR, Rae IJ, Mann IR, Milling DK (2011a) On the nature of ULF wave power during nightside auroral activations and substorms: 1. Spatial distribution. J Geophys Res. doi:10.1029/2010JA015757 Google Scholar
- Murphy KR, Rae IJ, Mann IR, Walsh AP, Milling DK, Kale A (2011b) The dependence of Pi2 waveforms on periodic velocity enhancements within bursty bulk flows. Ann Geophys. doi:10.5194/angeo-29-493-2011 Google Scholar
- Murphy KR, Mann IR, Rae IJ, Waters CL, Frey HU, Kale A, Singer HJ, Anderson BJ, Korth H (2013) The detailed spatial structure of field-aligned currents comprising the substorm current wedge. J Geophys Res Space Phys 118:7714–7727. doi:10.1002/2013JA018979 View ArticleGoogle Scholar
- Murphy KR, Mann IR, Rae IJ, Walsh AP, Frey HU (2015) Inner magnetospheric onset preceding reconnection and tail dynamics during substorms: can substorms initiate in two different regions? J Geophys Res Space Phys 119:9684–9701. doi:10.1002/2014JA019795 View ArticleGoogle Scholar
- Newell PT, Lee AR, Liou K, Ohtani S-I, Sotirelis T, Wing S (2010) Substorm cycle dependence of various types of aurora. J Geophys Res 115:A09226. doi:10.1029/2010JA015331 Google Scholar
- Nishimura Y, Lyons L, Zou S, Angelopoulos V, Mende S (2010) Substorm triggering by new plasma intrusion: THEMIS all-sky imager observations. J Geophys Res 115:A07222. doi:10.1029/2009JA015166 Google Scholar
- Nose M et al (1998) Automated detection of Pi 2 pulsations using wavelet analysis: 1 method and an application for substorm monitoring. Earth Planets Space 50:773–783View ArticleGoogle Scholar
- Ohtani S (2004) Flowbursts in the plasma sheet and auroral substorm onset: observational constraints on connection between midtail and near-earth substorm processes. Space Sci Rev 113:77–96. doi:10.1023/B:SPAC.0000042940.59358.2f.493-509 View ArticleGoogle Scholar
- Peticolas LM et al (2008) The time history of events and macroscale interactions during substorms (THEMIS) education and outreach (E/PO) program. Space Sci Rev 141:557–583. doi:10.1007/S11214-008-9458-5
- Posch JL et al (2007) Statistical observations of spatial characteristics of Pi1B pulsations. J Atmos Sol Terr Phys 69(15):1775–1796. doi:10.1016/j.jastp.2007.07.015 View ArticleGoogle Scholar
- Pu ZY et al (1999) Ballooning instability in the presence of a plasma flow: a synthesis of tail reconnection and current disruption models for the initiation of substorms. J Geophys Res 104(A5):10235–10248. doi:10.1029/1998JA900104 View ArticleGoogle Scholar
- Rae IJ et al. (2007) Pi2 pulsations: field line resonances or a driven response? In: Syrjäsuo M, Donovan EF (eds) International Conference on Substorms (ICS-8), University of Calgary, CanadaGoogle Scholar
- Rae IJ, Watt CEJ (2016) ULF waves above the nightside auroral oval during substorm onset. In: Keiling A, Lee DH, Nakariakov V (eds) Low frequency waves in space plasmas. Wiley, Hoboken. doi:10.1002/9781119055006.ch7 Google Scholar
- Rae IJ et al (2009a) Timing and localization of ionospheric signatures associated with substorm expansion phase onset. J Geophys Res. doi:10.1029/2008JA013559 Google Scholar
- Rae IJ et al (2009b) Near-earth initiation of a terrestrial substorm. J Geophys Res 114:A07220. doi:10.1029/2008JA013771 Google Scholar
- Rae IJ, Watt CEJ, Mann IR, Murphy KR, Samson JC, Kabin K, Angelopoulos V (2010) Optical characterization of the growth and spatial structure of a substorm onset arc. J Geophys Res 115:A10222. doi:10.1029/2010JA015376 Google Scholar
- Rae IJ, Murphy KR, Watt CEJ, Mann IR (2011) On the nature of ULF wave power during nightside auroral activations and substorms: 2. Temporal evolution. J Geophys Res. doi:10.1029/2010JA015762.3771 Google Scholar
- Rae IJ, Watt CEJ, Murphy KR, Frey HU, Ozeke LG, Milling DK, Mann IR (2012) The correlation of ULF waves and auroral intensity before, during and after substorm expansion phase onset. J Geophys Res 117:A08213. doi:10.1029/2012JA017534 Google Scholar
- Roux A, Perraut S, Robert P, Morane A, Pedersen A, Korth A, Kremser G, Aparicio B, Rodgers D, Pellinen R (1991) Plasma sheet instability related to the westward traveling surge. J Geophys Res 96(A10):17697–17714. doi:10.1029/91JA01106
- Russell CT et al (2008) THEMIS ground-based magnetometers. Space Sci Rev 141:389–412View ArticleGoogle Scholar
- Samson JC, Lyons LR, Newell PT, Creutzberg F, Xu B (1992) Proton aurora and substorm intensifications. Geophys Res Lett 19:2167–2170View ArticleGoogle Scholar
- Sibeck DG, Angelopoulos V (2008) THEMIS science objectives and mission phases. Space Sci Rev 141:35–59. doi:10.1007/s11214-008-9393-5 View ArticleGoogle Scholar
- Tamao T (1964) The structure of three-dimensional hydromagnetic waves in a uniform cold plasma. J Geomagn Geoelectr 18:89–114View ArticleGoogle Scholar
- Uozumi T, Yumoto K, Kawano H, Yoshikawa A, Olson JV, Solovyev SI, Vershinin EF (2000) Characteristics of energy transfer of Pi 2 magnetic pulsations: latitudinal dependence. Geophys Res Lett 27:1619View ArticleGoogle Scholar
- Uozumi T et al (2004) Propagation characteristics of Pi 2 magnetic pulsations observed at ground high-latitudes. J Geophys Res 109:A08203. doi:10.1029/2003JA009898 View ArticleGoogle Scholar
- Uozumi T et al (2009) Propagation characteristics of Pi 2 pulsations observed at high- and low-latitude MAGDAS/CPMN stations: a statistical study. J Geophys Res 114:A11207. doi:10.1029/2009JA014163 View ArticleGoogle Scholar
- Walsh AP et al (2010) Comprehensive ground-based and in situ observations of substorm expansion phase onset. J Geophys Res. doi:10.1029/2010JA015748 Google Scholar
- Wing S, Gkioulidou M, Johnson JR, Newell PT, Wang C-P (2013) Auroral particle precipitation characterized by the substorm cycle. J Geophys Res Space Phys 118:1022–1039. doi:10.1002/jgra.50160 View ArticleGoogle Scholar
- Yeoman TK, Orr D (1989) Phase and spectral power of mid- latitude Pi2 pulsations: evidence for a plasmaspheric cavity resonance. Planet Space Sci 37:1367–1383View ArticleGoogle Scholar