We present a new quantitative technique that determines the times and durations of substorm expansion and recovery phases and possible growth phases based on percentiles of the rate of change of auroral electrojet indices. By being able to prescribe different percentile values, we can determine the onset and duration of substorm phases for smaller or larger variations of the auroral index or indeed any auroral zone ground‐based magnetometer data. We apply this technique to the SuperMAG AL (SML) index and compare our expansion phase onset times with previous lists of substorm onsets. We find that more than 50% of events in previous lists occur within 20 min of our identified onsets. We also present a comparison of superposed epoch analyses of SML based on our onsets identified by our technique and existing onset lists and find that the general characteristics of the substorm bay are comparable. By prescribing user‐defined thresholds, this automated, quantitative technique represents an improvement over any visual identification of substorm onsets or indeed any fixed threshold method.

Abstract. We perform two superposed epoch analyses of the auroral evolution during substorms using the FUV instrument on the Imager for Magnetopause-to-Aurora Global Explorer (IMAGE) spacecraft. The larger of the two studies includes nearly 2000 substorms. We subdivide the substorms by onset latitude, a measure of the open magnetic flux in the magnetosphere, and determine average auroral images before and after substorm onset, for both electron and proton aurora. Our results indicate that substorms are more intense in terms of auroral brightness when the open flux content of the magnetosphere is larger, and that magnetic flux closure is more significant. The increase in auroral brightness at onset is larger for electrons than protons. We also show that there is a dawn-dusk offset in the location of the electron and proton aurora that mirrors the relative locations of the region 1 and region 2 current systems. Superposed epoch analyses of the solar wind, interplanetary magnetic field, and geomagnetic indices for the substorms under study indicate that dayside reconnection is expected to occur at a faster rate prior to low latitude onsets, but also that the ring current is enhanced for these events.

We present the first multievent study of the spatial and temporal structuring of the aurora to provide statistical evidence of the near‐Earth plasma instability which causes the substorm onset arc. Using data from ground‐based auroral imagers, we study repeatable signatures of along‐arc auroral beads, which are thought to represent the ionospheric projection of magnetospheric instability in the near‐Earth plasma sheet. We show that the growth and spatial scales of these wave‐like fluctuations are similar across multiple events, indicating that each sudden auroral brightening has a common explanation. We find statistically that growth rates for auroral beads peak at low wave number with the most unstable spatial scales mapping to an azimuthal wavelength λ≈ 1700–2500 km in the equatorial magnetosphere at around 9–12 RE. We compare growth rates and spatial scales with a range of theoretical predictions of magnetotail instabilities, including the Cross‐Field Current Instability and the Shear Flow Ballooning Instability. We conclude that, although the Cross‐Field Current instability can generate similar magnitude of growth rates, the range of unstable wave numbers indicates that the Shear Flow Ballooning Instability is the most likely explanation for our observations.

Using the total radiation belt electron content calculated from Van Allen Probe phase space density, the time‐dependent and global response of the outer radiation belt during storms is statistically studied. Using phase space density reduces the impacts of adiabatic changes in the main phase, allowing a separation of adiabatic and nonadiabatic effects and revealing a clear modality and repeatable sequence of events in storm time radiation belt electron dynamics. This sequence exhibits an important first adiabatic invariant (μ)‐dependent behavior in the seed (150 MeV/G), relativistic (1,000 MeV/G), and ultrarelativistic (4,000 MeV/G) populations. The outer radiation belt statistically shows an initial phase dominated by loss followed by a second phase of rapid acceleration, while the seed population shows little loss and immediate enhancement. The time sequence of the transition to the acceleration is also strongly μ dependent and occurs at low μ first, appearing to be repeatable from storm to storm.

Field‐aligned beams of suprathermal electrons, known as “strahl,” are a frequently observed constituent of solar wind plasma. However, the formation and interplanetary evolution of the strahl electron populations has yet to be fully understood. As strahl electrons travel away from the Sun, they move into regions of decreasing magnetic field strength and thus are subject to adiabatic focusing. However, the widths of strahl pitch angle distributions observed at 1 AU are significantly broader than expected. Previous investigations have found that the average observed strahl pitch angle width actually increases with heliocentric radial distance. This implies that strahl electrons must be subjected to some form of pitch angle scattering process or processes, details of which as of yet remain elusive. In this paper, we use Cassini electron measurements to examine strahl beams across a distance range of approximately 8 AU, from its Earth Flyby in 1999 until its insertion into orbit around Saturn in 2004. We find that, in general, there is a relatively constant rate of broadening of strahl pitch angle distributions with distance between ∼1 and 5.5 AU. Our results from beyond this distance indicate that the strahl population is likely to be completely scattered, presumably to form part of the halo. We find multiple energy dependences at different radial distances implying that there are multiple strahl scattering mechanisms in operation.

We present results from the first systematic survey of proton and electron pitch angle distributions in the magnetotail, based on Cluster CIS and PEACE data binned by proton plasma β (βp). The proton distributions conform to the canonical picture of magnetotail ions ‐ a boundary layer made up of Earthward streaming and bidirectional field‐aligned particles, consistent with recent observations of time‐varying beamlets, which gives way to a broadly isotropic central plasma sheet when βp ∼ 3. The electron distributions are significantly different from the canonical picture. A “boundary layer” made up of bidirectional field‐aligned electrons is observed to values of βp as high as 17. This boundary quickly gives way to perpendicular‐dominated electrons close to the neutral sheet. Hence, our results suggest that, on average, there is no extended, isotropic electron plasma sheet and that the proton plasma sheet is not routinely encountered until higher βp than commonly assumed.

Abstract. Dawn-dusk asymmetries are ubiquitous features of the coupled solar-wind-magnetosphere-ionosphere system. During the last decades, increasing availability of satellite and ground-based measurements has made it possible to study these phenomena in more detail. Numerous publications have documented the existence of persistent asymmetries in processes, properties and topology of plasma structures in various regions of geospace. In this paper, we present a review of our present knowledge of some of the most pronounced dawn-dusk asymmetries. We focus on four key aspects: (1) the role of external influences such as the solar wind and its interaction with the Earth's magnetosphere; (2) properties of the magnetosphere itself; (3) the role of the ionosphere and (4) feedback and coupling between regions. We have also identified potential inconsistencies and gaps in our understanding of dawn-dusk asymmetries in the Earth's magnetosphere and ionosphere.

We investigate how statistical properties of the rate of change R of the surface horizontal magnetic field in the United Kingdom differ during substorm expansion and recovery phases compared with other times. R is calculated from 1‐min magnetic field data from three INTERMAGNET observatories—Lerwick, Eskdalemuir, and Hartland and between 1996 and 2014—nearly two solar cycles. Substorm expansion and recovery phases are identified from the SuperMAG Lower index using the Substorm Onsets and Phases from Indices of the Electrojet method. The probability distribution of R is decomposed into categories of whether during substorm expansion and recovery phases, in enhanced convection intervals, or at other times. From this, we find that 54–56% of all extreme R values (defined as above the 99.97th percentile) occur during substorm expansion or recovery phases. By similarly decomposing the magnetic local time variation of the occurrence of large R values (>99th percentile), we deduce that 21–25% of large R during substorm expansion and recovery phases are attributable to the Disturbance Polar (DP)1 magnetic perturbation caused by the substorm current wedge. This corresponds to 10–14% of all large R in the entire data set. These results, together with asymptotic trends in occurrence probabilities, may indicate the two‐cell DP2 magnetic perturbation caused by magnetospheric convection as the dominant source of hazardous R > 600 nT/min that is potentially damaging to the U.K. National Grid. Thus, further research is needed to understand and model DP2, its mesoscale turbulent structure, and substorm feedbacks in order that GIC impact on the National Grid may be better understood and predicted.