Search

Magnetometers are a key component of heliophysics research providing valuable insight into the dynamics of electromagnetic field regimes and their coupling throughout the solar system. On satellites, magnetometers provide detailed observations of the extension of the solar magnetic field into interplanetary space and of planetary environments. At Earth, magnetometers are deployed on the ground in extensive arrays spanning the polar cap, auroral and sub-auroral zone, mid- and low-latitudes and equatorial electrojet with nearly global coverage in azimuth (longitude or magnetic local time—MLT). These multipoint observations are used to diagnose both ionospheric and magnetospheric processes as well as the coupling between the solar wind and these two regimes at a fraction of the cost of in-situ instruments. Despite their utility in research, ground-based magnetometer data can be difficult to use due to a variety of file formats, multiple points of access for the data, and limited software. In this short article we review the Open-Source Python library GMAG which provides rapid access to ground-based magnetometer data from a number of arrays in a Pandas DataFrame, a common data format used throughout scientific research.

We present analytic expressions for ULF wave-derived radiation belt radial diffusion coefficients, as a function of L and Kp, which can easily be incorporated into global radiation belt transport models. The diffusion coefficients are derived from statistical representations of ULF wave power, electric field power mapped from ground magnetometer data, and compressional magnetic field power from in situ measurements. We show that the overall electric and magnetic diffusion coefficients are to a good approximation both independent of energy. We present example 1-D radial diffusion results from simulations driven by CRRES-observed time-dependent energy spectra at the outer boundary, under the action of radial diffusion driven by the new ULF wave radial diffusion coefficients and with empirical chorus wave loss terms (as a function of energy, Kp and L). There is excellent agreement between the differential flux produced by the 1-D, Kp-driven, radial diffusion model and CRRES observations of differential electron flux at 0.976 MeV—even though the model does not include the effects of local internal acceleration sources. Our results highlight not only the importance of correct specification of radial diffusion coefficients for developing accurate models but also show significant promise for belt specification based on relatively simple models driven by solar wind parameters such as solar wind speed or geomagnetic indices such as Kp. Key Points Analytic expressions for the radial diffusion coefficients are presented The coefficients do not dependent on energy or wave m value The electric field diffusion coefficient dominates over the magnetic

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.

Electromagnetic ion cyclotron (EMIC) waves are believed to be an important source of pitch angle scattering driven relativistic electron loss from the radiation belts. To date, investigations of this precipitation have been largely theoretical in nature, limited to calculations of precipitation characteristics based on wave observations and small-scale studies. Large-scale investigation of EMIC wave-driven electron precipitation has been hindered by a lack of combined wave and precipitation measurements. Analysis of electron flux data from the POES (Polar Orbiting Environmental Satellites) spacecraft has been suggested as a means of investigating EMIC wave-driven electron precipitation characteristics, using a precipitation signature particular to EMIC waves. Until now the lack of supporting wave measurements for these POES-detected precipitation events has resulted in uncertainty regarding the driver of the precipitation. In this paper we complete a statistical study comparing POES precipitation measurements with wave data from several ground-based search coil magnetometers; we further present a case study examining the global nature of this precipitation. We show that a significant proportion of the precipitation events correspond with EMIC wave detections on the ground; for precipitation events that occur directly over the magnetometers, this detection rate can be as high as 90%. Our results demonstrate that the precipitation region is often stationary in magnetic local time, narrow in L, and close to the expected plasmapause position. Predominantly, the precipitation is associated with helium band rising tone Pc1 waves on the ground. The success of this study proves the viability of POES precipitation data for investigating EMIC wave-driven electron precipitation.

Very low frequency (VLF) transmissions propagating between the conducting Earth's surface and lower edge of the ionosphere have been used for decades to study the effect of space weather events on the upper atmosphere. The VLF response to these events can only be quantified by comparison of the observed signal to the estimated quiet time or undisturbed signal levels, known as the quiet day curve (QDC). A common QDC calculation approach for periods of investigation of up to several weeks is to use observations made on quiet days close to the days of interest. This approach is invalid when conditions are not quiet around the days of interest. Longer-term QDCs have also been created from specifically identified quiet days within the period and knowledge of propagation characteristics. This approach is time consuming and can be subjective. We present three algorithmic techniques, which are based on either (1) a mean of previous days' observations, (2) principal component analysis, or (3) the fast Fourier transform (FFT), to calculate the QDC for a long-period VLF data set without identification of specific quiet days as a basis. We demonstrate the effectiveness of the techniques at identifying the true QDCs of synthetic data sets created to mimic patterns seen in actual VLF data including responses to space weather events. We find that the most successful technique is to use a smoothing method, developed within the study, on the data set and then use the developed FFT algorithm. This technique is then applied to multiyear data sets of actual VLF observations.

On 31 May 2013 several rising-tone electromagnetic ion-cyclotron (EMIC) waves with intervals of pulsations of diminishing periods (IPDP) were observed in the magnetic local time afternoon and evening sectors during the onset of a moderate/large geomagnetic storm. The waves were sequentially observed in Finland, Antarctica, and western Canada. Co-incident electron precipitation by a network of ground-based Antarctic Arctic Radiation-belt Dynamic Deposition VLF Atmospheric Research Konsortia (AARDDVARK) and riometer instruments, as well as the Polar-orbiting Operational Environmental Satellite (POES) electron telescopes, was also observed. At the same time POES detected 30-80 keV proton precipitation drifting westwards at locations that were consistent with the ground-based observations, indicating substorm injection. Through detailed modelling of the combination of ground and satellite observations the characteristics of the EMIC-induced electron precipitation were identified as: latitudinal width of 2-3° or ΔL=1 Re, longitudinal width ~50° or 3 hours MLT, lower cut off energy 280 keV, typical flux 1×104 el. cm-2 sr-1 s-1 >300 keV. The lower cutoff energy of the most clearly defined EMIC rising tone in this study confirms the identification of a class of EMIC-induced precipitation events with unexpectedly low energy cutoffs of

We present case studies of the evolution of magnetic wave amplitudes and auroral intensity through the late growth phase and the expansion phase of the substorm cycle. We present strong evidence that substorm-related auroral enhancements are clearly and demonstrably linked to ULF wave amplitudes observed at the same location. In most cases, we find that the highest correlations are observed when the magnetometer time series is advanced in time, indicating that the ULF wave amplitudes start to grow before measured auroral intensities, though interestingly this is not always the case. Further we discuss the four possible reasons that may be able to explain both the timing and the high correlations between these two phenomena, including: a simple coincidence, an artifact of instrumental effects, the response of the ionosphere to magnetic waves and auroral particle precipitation, and finally that ULF waves and auroral particle precipitation are physically linked. We discount coincidence and instrumental effects since in the studies presented here they are unlikely or in general will contribute negligible effects, and we find that the ionospheric response to waves and precipitation can explain some, but not all of the results contained within this paper. Specifically, ionospheric response to substorm waves and auroral precipitation cannot explain that the result that previous studies have shown, that onset of ULF wave activity and the onset of auroral particle precipitation occur at the same time and in the same location. This leaves the possibility that ULF waves and auroral particles are physically linked.

Using a discrete wavelet transform with a Meyer wavelet basis, we present a new quantitative algorithm for determining the onset time of Pi1 and Pi2 ULF waves in the nightside ionosphere with ∼20- to 40-s resolution at substorm expansion phase onset. We validate the algorithm by comparing both the ULF wave onset time and location to the optical onset determined by the Imager for Magnetopause-to-Aurora Global Exploration (IMAGE)–Far Ultraviolet Imager (FUV) instrument. In each of the six events analyzed, five substorm onsets and one pseudobreakup, the ULF onset is observed prior to the global optical onset observed by IMAGE at a station closely conjugate to the optical onset. The observed ULF onset times expand both latitudinally and longitudinally away from an epicenter of ULF wave power in the ionosphere. We further discuss the utility of the algorithm for diagnosing pseudobreakups and the relationship of the ULF onset epicenter to the meridians of elements of the substorm current wedge. The importance of the technique for establishing the causal sequence of events at substorm onset, especially in support of the multisatellite Time History of Events and Macroscale Interactions During Substorms (THEMIS) mission, is also described.

Despite the characterization of the auroral substorm more than 40 years ago, controversy still surrounds the processes triggering substorm onset initiation. That stretching of the Earth's magnetotail following the addition of new nightside magnetic flux from dayside reconnection powers the substorm is well understood; the trigger for explosive energy release at substorm expansion phase onset is not. Using ground-based data sets with unprecedented combined spatial and temporal coverage, we report the discovery of new localized and contemporaneous magnetic wave and small azimuthal scale auroral signature of substorm onset. These local auroral arc undulations and magnetic field signatures rapidly evolve on second time scales for several minutes in advance of the release of the auroral surge. We also present evidence from a conjugate geosynchronous satellite of the concurrent magnetic onset in space as the onset of magnetic pulsations in the ionosphere, to within technique error. Throughout this time period, the more poleward arcs that correspond to the auroral oval which maps to the central plasma sheet remain undisturbed. There is good evidence that flows from the midtail crossing the plasma sheet can generate north-south auroral structures, yet no such auroral forms are seen in this event. Our observations present a severe challenge to the standard hypothesis that magnetic reconnection in stretched magnetotail fields triggers onset, indicating substorm expansion phase initiation occurs on field lines that are close to the Earth, as bounded by observations at geosynchronous orbit and in the conjugate ionosphere.

In this paper, we present case studies of the optical and magnetic signatures of the characteristics of the first minute of substorm expansion phase onset observed in the ionosphere. We find that for two isolated substorms, the onset of magnetic pulsations in the 24–96 s period wavelet band are colocated in time and space with the formation and development of small-scale optical undulations along the most equatorward preexisting auroral arc prior to auroral breakup. These undulations undergo an inverse spatial cascade into vortices prior to the release of the westward traveling surge. We also present a case study of a multiple activation substorm, whereby discrete onsets of ULF wave power above a predetermined quiet time threshold are shown to be associated with specific optical intensifications and brightenings. Moreover, in the multiple activation substorm event, we show that neither the formation of the small-scale undulations nor the formation of similar structures along a north–south aligned arc is sufficient to produce auroral breakup associated with expansion phase onset. It is only ∼10 min after these two disparate activation regions initiate that auroral breakup and the subsequent formation of a westward traveling surge occur. We discuss the implications of these results in terms of the triggering mechanisms likely to be occurring during these specific events.