We combine wavelet analysis and data fusion to investigate geomagnetically induced currents (GICs) on the Mäntsälä pipeline and the associated horizontal geomagnetic field, BH, variations during the late main phase of the 17 March 2013 geomagnetic storm. The wavelet analysis decomposes the GIC and BH signals at increasing ‘scales’ to show distinct multi-minute spectral features around the GIC spikes. Four GIC spikes > 10 A occurred while the pipeline was in the dusk sector – the first sine-wave-like spike at ~16 UT was ‘compound.’ It was followed by three ‘self-similar’ spikes two hours later. The contemporaneous multi-resolution observations from ground-(magnetometer, SuperMAG, SuperDARN), and space-based (AMPERE, TWINS) platforms capture multi-scale activity to reveal two magnetospheric modes causing the spikes. The GIC at ~16 UT occurred in two parts with the negative spike associated with a transient sub-auroral eastward electrojet that closed a developing partial ring current (PRC) loop, whereas the positive spike developed with the arrival of the associated mesoscale flow-channel in the auroral zone. The three spikes between 18-19 UT were due to bursty bulk flows (BBFs). We attribute all spikes to flow-channel injections (substorms) of varying scales. We use previously published MHD simulations of the event to substantiate our conclusions, given the dearth of timely in-situ satellite observations. Our results show that multi-scale magnetosphere-ionosphere activity that drives GICs can be understood using multi-resolution analysis. This new framework of combining wavelet analysis with multi-platform observations opens a research avenue for GIC investigations and other space weather impacts.

A document by Bhagyashree Waghule. Click on the document to view its contents.

Dr. Jennifer Gannon passed away suddenly on 2024 May 2 in Greenbelt, MD. Dr. Gannon had served as an editor for the Space Weather journal since April 2019, and she was the longest-serving editor on the current board, having started under the previous editor-in-chief, Dr. Delores Knipp. Abstract content goes here

We provide high-resolution maps of quasi-static Poynting flux (PF) in each hemisphere based on nine-satellite years of Defense Meteorological Satellite Program (DMSP) data. Conjugate comparisons from ~850 km reveal more quasi-static PF arriving in the northern hemisphere (NH) than the southern hemisphere (SH). This tendency is clear in the dawn-dusk sectors and during intervals when Kp < 3, which accounts for ~80% of the study interval. Summer-to-summer comparisons indicate this asymmetry is partially associated with more NH solar illumination, which supports stronger NH field-aligned currents (FAC). Differing hemispheric FAC configurations may also play a role. Our findings support and broaden earlier reports of similar NH preference for the deposition of Alfvenic PF. Regionally the NH has stronger dusk-region PF, while the SH has stronger mid-morning PF. We find PF deposition in the near-cusp regions that rivals and often exceeds the PF intensity in the auroral zones.

In this study, a new high-latitude empirical model is introduced, named for Auroral energy Spectrum and High-Latitude Electric field variabilitY (ASHLEY). This model aims to improve specifications of soft electron precipitations and electric field variability that are not well represented in existing high-latitude empirical models. ASHLEY consists of three components, ASHLEY-A, ASHLEY-E and ASHLEY-Evar, which are developed based on the electron precipitation and bulk ion drift measurements from the Defense Meteorological Satellite Program (DMSP) satellites during the most recent solar cycle. On the one hand, unlike most existing high-latitude electron precipitation models, which have assumptions about the energy spectrum of incident electrons, the electron precipitation component of ASHLEY, ASHLEY-A, provides the differential energy fluxes in the 19 DMSP energy channels under different geophysical conditions without making any assumptions about the energy spectrum. It has been found that the relaxation of spectral assumptions significantly improves soft electron precipitation specifications with respect to a Maxwellian spectrum (up to several orders of magnitude). On the other hand, ASHLEY provides consistent mean electric field and electric field variability under different geophysical conditions by ASHLEY-E and ASHLEY-Evar components, respectively. This is different from most existing electric field models which only focus on the large-scale mean electric field and ignore the electric field variability. Furthermore, the consistency between the electric field and electron precipitation is better taken into account in ASHLEY.

Poynting flux calculated from LEO spacecraft in-situ ion drift and magnetic field measurements is an important measure of energy exchange between the magnetosphere and ionosphere. Defense Meteorological Satellite Program (DMSP) spacecraft provide an extensive back-catalog of ion drift and magnetic perturbation measurements, from which quasi-steady Poynting flux could be calculated. However, since DMSP are operations-focused spacecraft, data must be carefully preprocessed for research use. We describe an automated approach for calculating earthward Poynting flux focusing on pre-processing and quality control. We produce a Poynting flux dataset using 9 satellite-years of DMSP F15, F16 and F18 observations. To validate our process we inter-compare Poynting flux from different spacecraft using more than 2000 magnetic conjunction events. We find no serious systematic differences. We further describe and apply an equal-area binning technique to obtain average spatial patterns of Poynting flux, magnetic perturbation, electric field and ion drift velocity. We perform our analysis using all components of electric and magnetic field and comment on the adverse consequences of the typical single-electric-field-component DMSP Poynting flux approximation on inter-spacecraft agreement. Including full-field components significantly increases the relative strength of near-cusp Poynting flux and increases the integrated high-latitude Poynting flux by ~25%

For southward interplanetary magnetic field (IMF) during local summer, the hemispherically integrated Poynting flux estimated by FAST-satellite-derived empirical models is significantly larger for the northern hemisphere (NH) than for the southern hemisphere (SH). In order to test whether the difference is statistically significant, the model uncertainties have been estimated by dividing the data sets for each hemisphere into two nonintersecting sets, and separately constructing the model using each of the four sets. The model uncertainty appears to be smaller than the estimated asymmetry. The asymmetry is mostly absent when the IMF is northward, except there is some evidence that it may actually reverse during local winter. The phenomena is coupled with what appears to be a more distinct two-cell convection pattern in the NH, and a possibly greater cusp contribution in the SH. All this suggests an effect of magnetosphere-ionosphere-thermosphere coupling, probably related to asymmetries in Earth’s geomagnetic field.

Atmospheric drag describes the main perturbing force of the atmosphere on the orbital trajectories of near-Earth orbiting satellites. The ability to accurately model atmospheric drag is critical for precise satellite orbit determination and collision avoidance. Assuming we know atmospheric winds and satellite velocity, area and mass, the primary sources of uncertainty in atmospheric drag include mass density of the space environment and the spacecraft drag coefficient, CD. Historically, much of the focus has been on physically or empirically estimating mass density, while CD is treated as a fitting parameter or fixed value. Presently, CD can be physically modeled through energy and momentum exchange processes between the atmospheric gas particles and the satellite surface. However, physical CD models rely on assumptions regarding the scattering and adsorption of atmospheric particles, and these responses are driven by atmospheric composition and temperature. Modifications to these assumptions can cause CD to change by up to ~40%. The nature and magnitude of these changes also depend on the shape of the spacecraft. We can check the consistency of the CD model assumptions by comparing densities derived from satellite drag measurements and computed CD values for satellites of different shapes orbiting in the same space environment. Since all of the satellites should see the same density, offsets in the derived densities should be attributable to inconsistencies in the CD model. Adjusting the CD model scattering assumptions can improve derived density consistency among the different satellites and inform the physics behind CD modeling. In turn, these efforts will help to reduce uncertainty in CD, leading to improved atmospheric drag estimates.