[1] We use current density data from the Active Magnetosphere and Planetary Electrodynamics Response Experiment (AMPERE) to identify the location of maximum region 1 current at all magnetic local times (MLTs). We term this location the R1 oval. Comparing the R1 oval location with particle precipitation boundaries identified in DMSP data, we find that the R1 oval is located on average within 1 ı of particle signatures associated with the open/closed field line boundary (OCB) across dayside and nightside MLTs. We hence conclude that the R1 oval can be used as a proxy for the location of the OCB. Studying the amount of magnetic flux enclosed by the R1 oval during the substorm cycle, we find that the R1 oval flux is well organized by it: during the growth phase the R1 oval location moves equatorward as the amount of magnetic flux increases whereas after substorm expansion phase onset significant flux closure occurs as the R1 current location retreats to higher latitudes. For about 15 min after expansion phase onset, the amount of open magnetic flux continues to increase indicating that dayside reconnection dominates over nightside reconnection. In the current density data, we find evidence of the substorm current wedge and also show that the dayside R1 currents are stronger than their nightside counterpart during the substorm growth phase, whereas after expansion phase onset, the nightside R1 currents dominate. Our observations of the current distribution and OCB movement during the substorm cycle are in excellent agreement with the expanding/contracting polar cap paradigm.

Trans‐ionospheric high frequency (HF: 3–30 MHz) signals experience strong attenuation following a solar flare‐driven sudden ionospheric disturbance (SID). Solar flare‐driven HF absorption, referred to as short‐wave fadeout, is a well‐known impact of SIDs, but the initial Doppler frequency shift phenomena, also known as “Doppler flash” in the traveling radio wave is not well understood. This study seeks to advance our understanding of the initial impacts of solar flare‐driven SID using a physics‐based whole atmosphere model for a specific solar flare event. First, we demonstrate that the Doppler flash phenomenon observed by Super Dual Auroral Radar Network (SuperDARN) radars can be successfully reproduced using first‐principles based modeling. The output from the simulation is validated against SuperDARN line‐of‐sight Doppler velocity measurements. We then examine which region of the ionosphere, D, E, or F, makes the largest contribution to the Doppler flash. We also consider the relative contribution of change in refractive index through the ionospheric layers versus lowered reflection height. We find: (a) the model is able to reproduce radar observations with an root‐median‐squared‐error and a mean percentage error (δ) of 3.72 m/s and 0.67%, respectively; (b) the F‐region is the most significant contributor to the total Doppler flash (∼48%), 30% of which is contributed by the change in F‐region's refractive index, while the other ∼18% is due to change in ray reflection height. Our analysis shows lowering of the F‐region's ray reflection point is a secondary driver compared to the change in refractive index.

Ion drag is known to play an important role in driving neutral thermosphere circulation at auroral latitudes, especially during the main phase of geomagnetic storms. During the recovery phase, the neutrals are known to drive the ions and generate ionospheric electric fields and currents via the disturbance dynamo mechanism. At midlatitudes, the precise interplay between ions and neutrals is less understood largely because of the paucity of measurements that have been available. In this work, we investigate ion‐neutral coupling at middle latitudes using colocated ion drift velocity measurements obtained from Super Dual Auroral Radar Network radars and neutral wind velocity and temperature measurements obtained from the North American Thermosphere Ionosphere Observing Network (NATION) Fabry‐Perot interferometers. We examine one recent storm period on 2–3 October 2013 during both the main phase and late recovery phase. By using ion‐neutral momentum exchange theory and a time‐lagged correlation analysis, we analyze the coupling time scales and dominant driving mechanisms. We observe that during the main phase the neutrals respond to the ion convection on a time scale of ∼84 min which is significantly faster than what would be expected from local ion drag momentum forcing alone. This suggests that other storm time influences are important for driving the neutrals during the main phase, such as Joule heating. During the late recovery phase, the neutrals are observed to drive the ion convection without any significant time delay, consistent with the so‐called “neutral fly wheel effect” or disturbance dynamo persisting well into the late recovery phase.

The existence of Birkeland magnetic field‐aligned current (FAC) system was proposed more than a century ago, and it has been of immense interest for investigating the nature of solar wind‐magnetosphere‐ionosphere coupling ever since. In this paper, we present the first application of deep learning architecture for modeling the Birkeland currents using data from the Active Magnetosphere and Planetary Electrodynamics Response Experiment (AMPERE). The model uses a 1‐hr time history of several different parameters such as interplanetary magnetic field (IMF), solar wind, and geomagnetic and solar indices as inputs to determine the global distribution of Birkeland currents in the Northern Hemisphere. We present a comparison between our model and bin‐averaged statistical patterns under steady IMF conditions and also when the IMF is variable. Our deep learning model shows good agreement with the bin‐averaged patterns, capturing several prominent large‐scale features such as the Regions 1 and 2 FACs, the NBZ current system, and the cusp currents along with their seasonal variations. However, when IMF and solar wind conditions are not stable, our model provides a more accurate view of the time‐dependent evolution of Birkeland currents. The reconfiguration of the FACs following an abrupt change in IMF orientation can be traced in its details. The magnitude of FACs is found to evolve with e‐folding times that vary with season and MLT. When IMF Bz turns southward after a prolonged northward orientation, NBZ currents decay exponentially with an e‐folding time of ∼25 min, whereas Region 1 currents grow with an e‐folding time of 6–20 min depending on the MLT.

Long-lasting Pc5 ultralow frequency (ULF) waves spanning the dayside and extending from L ∼ 5.5 into the polar cap region were observed by conjugate ground magnetometers. Observations from MMS satellites in the magnetosphere and magnetometers on the ground confirmed that the ULF waves on closed field lines were due to fundamental toroidal standing Alfvén waves. Monochromatic waves at lower latitudes tended to maximize their power away from noon in both the morning and afternoon sectors, while more broadband waves at higher latitudes tended to have a wave power maximum near noon. The wave power distribution and MMS satellite observations during the magnetopause crossing indicate surface waves on a Kelvin-Helmholtz (KH) unstable magnetopause coupled with standing Alfvén waves. The more turbulent ion foreshock during an extended period of radial interplanetary magnetic field (IMF) likely plays an important role in providing seed perturbations for the growth of the KH waves. These results indicate that the Pc5 waves observed on closed field lines and on the open field lines of the polar cap were from the same source. Plain Language Summary The Earth's magnetic field lines can oscillate at ultralow frequencies (ULF: 1 mHz to 5 Hz). These natural oscillations of closed magnetic field lines, analogous to vibrations on a stretched string, are also called geomagnetic pulsations or ULF waves. ULF waves play a key role in the transfer of energy from terrestrial space to Earth's upper atmosphere. In this study, we report a long-lasting large spatial scale ULF wave event observed by ground observatories from both hemispheres. Together with satellite measurements in space, we are able to confirm that these waves were driven by upstream turbulent structures due to the interaction between matter and electromagnetic fields emitted from the Sun and the Earth's outer atmosphere and magnetic field.