The plasma wave instruments on both Voyager spacecraft have observed electron plasma oscillations in the very local interstellar medium (VLISM). The generally accepted explanation of these events is that the electron foreshock of shocks in the VLISM comprise electron beams in the range of 10 to 100 eV that are unstable to Langmuir waves, or electron plasma oscillations. Further, at least some of these events have been tied to solar transients departing the Sun more than a year earlier that evolve as they propagate outward. These disturbances are led by shocks and the impulse of these on the heliospause results in some of the shock impulse continuing into the VLISM. Previously, Voyager 1 had detected the most distant evidence of these transients at about 145 AU. In August 2025 Voyager 2 detected electron plasma oscillations near 140 AU. A simple model of the propagation of this disturbance suggests a transient from the Sun in 2022 as its source, near the beginning of the current solar maximum. New Horizons observed a series of shocks in 2022 – 2023 at heliocentric distances near 55 AU that could be related to the Voyager 2 event. Given these events occur early in solar cycle 25, it is possible additional shocks will be detected by Voyager and enable us to extend the distance over which these disturbances can travel in the VLISM.We further relate some of the transients observed by the Voyager plasma wave instruments to global models of the VLISM density and magnetic field (Fraternale et al., 2026). For example, these models show the increased density and magnetic field associated with the so-called pf2 (pressure front 2) described by Burlaga et al. (2021). We can now show that the 2-3 kHz radio emissions observed by the Voyagers in the early 1980’s, 1990’s, and 2000’s are related to density structures just beyond the heliopause presumed to be associated with global merged interaction regions stemming from very active solar conditions.
The Systematic Study of lower Thermosphere Energetics by a Rocket (SYSTER) is a multi-payload sounding rocket designed to study Joule heating at auroral latitudes. SYSTER provides pathfinder measurements for a future thermospheric ’dipper’ spacecraft, currently being evaluated by the ESA-NASA Lower Thermosphere-Ionosphere Science (ENLOTIS) working group.
Pulsating aurora is a common and long-lasting type of aurora composed of precipitating electrons with energies up to 100s keV to MeV. Since it is more energetic it can deposit energy into the D region and below, and accounts for a large amount of energy transferred from the magnetosphere into the ionosphere. To understand the effects of pulsating aurora on the ionosphere/atmosphere we seek to understand its precipitating energy flux spectra. One method to do this is to use an inversion method outlined by Semeter and Kamalabadi [1], that uses an electron density vs altitude profile from an Incoherent Scatter Radar to find the precipitating energy flux. This method uses both an atmospheric chemistry model, that finds the recombination coefficient, and an ionization rate model. In this study we seek to validate a selection of these models by comparing the resulting energy flux to spacecraft data. We use the ELFIN bounce loss cone electron flux during conjunctions with the Poker Flat Incoherent Scatter Radar (PFISR). Since ELFIN’s electron energy range is from 50 keV to 5 MeV it overlaps with PFISR’s inverted flux range of 1 keV to 500 keV, and as a low earth orbiting spacecraft, its bounce loss cone flux will give the precipitating fluxes. While not a new method, this is the first time this inversion method has been validated with spacecraft data.
Particle acceleration is a commonly observed phenomenon at dipolarization fronts. Many studies have attempted to determine the acceleration mechanism, with betatron acceleration being a major candidate. In previous work, we attempted to match the observed change in electron energy to the change predicted by betatron acceleration, but found that although this worked in some cases, overall betatron acceleration alone could not describe the observed energy spectrum changes. In this work, we attempted to study whether ion acceleration showed similar behavior and whether a quasi-adiabatic correction would be more accurate. On average the betatron acceleration equation overestimated the observed acceleration and the quasi-adiabatic correction did not account for the difference, although there are limitations to this study due to data fidelity. We then turned to study whether our assumptions about the source population having the same phase space density as the cold pre-existing background population in the plasma sheet are valid. We indirectly studied this by comparing the relative abundances of O+ and He++ as proxies for ionospheric and solar wind populations respectively. We found the betatron acceleration equation method performs slightly better when there is a stronger ionospheric component. This suggests that when more plasma containing O+ is present in the dipolarization front, it indicates that the source population is more local and therefore this method of using betatron acceleration is more valid.
This paper outlines the science and basic design choices associated with a mission concept study known as the LaboratOry for the Behavior of the SloT Region (LOBSTR). This mission concept focuses on energetic particles, both electrons and protons, as they impinge upon the slot region in the Van Allen radiation belts around Earth. In particular, it emphasizes the drift dynamics of particles that were not captured by Van Allen Probes. We conceptualize a mission, utilizing state-of-the-art instruments and components, and calculate the mission’s orbit, thrust, and radiation requirements using industry-standard methods. The concept uses two SmallSats in a near-equatorial orbit, with precise orbital timing to capture the desired dynamics. The total radiation dose and the details of the orbital dynamics are examined and found to be within the capabilities of current technology.
We made observations of magnetic field variations in association with pulsating auroras with the magneto-impedance sensor magnetometer (MIM) carried by the Loss through Auroral Microburst Pulsations (LAMP) sounding rocket that was launched at 11:27:30 UT on 5 March 2022 from Poker Flat Research Range, Alaska. At an altitude of 200-250 km, MIM detected clear enhancements of the magnetic field by 15-25 nT in both the northward and westward components. From simultaneous observations with the ground all-sky camera, we found that the footprint of LAMP at the 100 km altitude was located near the center of a pulsating auroral patch. The auroral patch had a dimension of similar to 90 km in latitude and similar to 25 km in longitude, and its major axis was inclined toward northwest. These observations were compared with results of a simple model calculation, in which local electron precipitation into the thin-layer ionosphere causes an elliptical auroral patch. The conductivity within the patch is enhanced in the background electric field and as a result, the magnetic field variations are induced around the auroral patch. The model calculation results can explain the MIM observations if the electric field points toward southeast and one of the model parameters is adjusted. We conclude that the pulsating auroral patch in this event was associated with a one-pair field-aligned current that consists of downward (upward) currents at the poleward (equatorward) edge of the patch. This current structure is maintained even if the auroral patch is latitudinally elongated. Magneto-impedance (MI) sensor was carried by a sounding rocket and first applied to magnetic field measurement in space Magnetic field variations were observed when the footprint of the payload was located near the center of an auroral patch The magnetic field variations are thought to be caused by field-aligned currents flowing at the edges of the pulsating auroral patch
The reprocessing of radiation belt electron flux measurements into phase space density (PSD) as a function of the adiabatic invariants is a widely-used method to address major questions regarding electron energization and loss in the outer radiation belt. In this reprocessing, flux measurements j (α, E) at local pitch angles α, energies E, and optionally magnetometer measurements B, are combined with a global magnetic field model to express the phase space density f (L*) in terms of the third invariant Φ ∝ 1/L* at fixed first and second invariants M and K. While the general framework of the calculation is agreed upon, implementation details vary amongst the literature, and the issue of magnetic field model dependence is rarely addressed. This work reviews the steps of the calculation with lists of commonly used implementation options. For the first time, analysis is presented to display the effect of doing the calculation with different implementation options and with different backing models (including both empirical and MHD-driven models). The results are summarized to inform evaluation of existing results and future efforts calculating and analyzing radiation belt electron phase space density. Three events are analyzed, and while differences are found, the primary structural interpretations of the phase space density analysis exhibit model independence.
Energetic electron precipitation (EEP) associated with pulsating aurora can transfer greater than 30 keV electrons from the outer radiation belt region into the upper atmosphere and can deplete atmospheric ozone via collisions that produce NOx and HOx molecules. Our knowledge of exactly how EEP occurs is incomplete. Previous studies have shown that pitch angle scattering between electrons and lower-band chorus waves can cause pulsating aurora associated with EEP and that substorms play an important role. In this work, we quantify the timescale of chorus wave decay following substorms and compare that to previously determined timescales. We find that the chorus decay e-folding time varies based on magnetic local time (MLT), magnetic latitude, and wave frequency. The fastest decay occurs for lower-band chorus in the 21 to 9 MLT region and compares well to the timescale of Troyer et al. (2022) for energetic pulsating aurora. We are able to further support this connection by modelling our findings in a quasi-linear diffusion simulation. These results provide observations of how chorus waves behave after substorms and add additional statistical evidence linking energetic pulsating aurora to substorm driven lower-band chorus waves.
Pulsating aurora, which consists of diffuse patches blinking on and off, is caused by pitch angle scattering of radiation belt electrons into the loss cone by lower-band chorus waves. Understanding the drivers of pulsating aurora is important as it is a long-lasting and widespread phenomenon, accounting for significant energy transfer from the solar wind into the ionosphere. Substorm injections, which transport electrons from the magnetotail into the inner magnetosphere, are one source of electrons in this region. Injections have been observed simultaneously with pulsating aurora during conjunctions between ground cameras and satellites. In addition, previous work has also shown that substorms can enhance chorus activity (the fundamental process that produces pulsating aurora), providing a mechanism linking substorms to pulsating aurora. To further study this connection, we used the Van Allen Probes and all-sky cameras to look at events where pulsating aurora and substorm injections were observed at different locations in Magnetic Local Time (MLT), rather than focusing only on conjunctions. To make this comparison, we calculated the drift rate of electrons originating from observed injections and projected their motion forward in time until their Magnetic Local Time was the same as the ground camera. When the electrons are located at the same MLT as the ground camera, the pulsating aurora they cause would most likely occur in the field of view of the camera. We compared the time drifting substorm-injected particles arrived at the MLT of the camera to when pulsating aurora was observed. We found several instances where the initiation or intensification of pulsating aurora was accompanied by the arrival of substorm-injected electrons. This observation gives further evidence that pulsating aurora can be enhanced by or occur after substorm injections.
We estimated the altitude of aurora by combining data from all-sky cameras at multiple places which were obtained during the LAMP sounding rocket experiment in Alaska on 5 March 2022. During the launch window of the rocket, three high-speed all-sky cameras were operative at three stations immediately below the trajectory of the rocket: Poker Flat, Venetie and Fort Yukon. The all-sky cameras captured all-sky images with a temporal resolution of 100 Hz (80 Hz for the Fort Yukon case). The method of altitude determination is based on analyses of time-series of the optical intensity obtained from the all-sky cameras in Venetie and Poker Flat covering the downrange area of the rocket trajectory. The estimated altitude of pulsating aurora during the rocket experiment was found to be consistent with that derived from the in-situ observation of precipitating electrons with a model of optical emission, which confirms the feasibility of deriving the emission altitude through correlation analyses using time-series. The estimated altitude of aurora decreased after the expansion onset of the substorm and stayed slightly below 100 km during the interval of pulsating aurora in the recovery phase. In particular, prompt and brief lowering of the auroral emission, well down to around 90 km, was detected during a transition of auroral form from discrete to diffuse which occurred similar to 10 min after the onset. This result implies an existence of a process causing harder electron precipitation operative soon after the start of the expansion phase of auroral substorm.
Turbulent and compressed sheath regions preceding interplanetary coronal mass ejections strongly impact electron dynamics in the outer radiation belt. Changes in electron flux can occur on timescales of tens of minutes, which are unlikely to be captured by a two‐satellite mission. The recently released Global Positioning System (GPS) data set generally has shorter revisit times (at L ∼ 4–8) owing to the large number of satellites in the constellation equipped with energetic particle detectors. Investigating electron fluxes at energies from 140 keV to 4 MeV and sheaths observed in 2012–2018, we show that the flux response to sheaths on a timescale of 6 hr, previously reported from Van Allen Probes (RBSP) data, is reproduced by GPS measurements. Furthermore, GPS data enables derivation of the response on a timescale of 30 min, which further confirms that the energy and L ‐shell dependent changes in electron flux are associated with the impact of the sheath. Sheath‐driven loss is underestimated over longer timescales as the electrons recover during the ejecta. We additionally show the response of electron phase space density (PSD), which is a key quantity in identifying non‐adiabatic loss from the system and electron energization through wave‐particle interactions. The PSD response is calculated from both RBSP and GPS data for the 6 hr timescale, as well as from GPS data for the 30 min timescale. The response is divided based on the geoeffectiveness of the sheaths revealing that electrons are effectively accelerated only during geoeffective sheaths, while loss commonly occurs during all sheaths.
Abstract Two types of filamentary currents (FCs) were observed inside a magnetic flux rope at the magnetopause by the Magnetospheric Multiscale mission. The first FC is identified as an electron vortex, while the other is a reconnecting current sheet. Stochastic electric fields were generated within the FCs, resulting in electron acceleration up to a few keV, similar to recent simulations of electron acceleration inside vortex, which is a second‐order Fermi acceleration. Furthermore, two FCs propagated at different speeds, causing compression in the region between them. Energetic electrons up to 200 keV were detected in the compressed region and displayed a double power‐law spectrum. Observations suggest that the electrons were mainly accelerated by betatron mechanism in the compressed region. The formation, evolution, and interaction of FCs provide a novel mechanism for electron acceleration. These results clearly show the significance of electron‐scale dynamics within flux rope.
Recognition for All: A Way Forward to Enhance Diversity, Equity and Inclusion in Space Physics M. Fraz Bashir,1 Amy M. Keesee,2, 3 Seth G. Claudepierre,4 Michael D. Hartinger,5 Elizabeth A. MacDonald,6 and Allison Jaynes7 Department of Earth, Planetary and Space Sciences, UCLA, Los Angeles, CA, USA∗ Department of Physics and Astronomy, University of New Hampshire, Durham, NH, USA Space Science Center, University of New Hampshire, Durham, NH, USA Department of Atmospheric and Oceanic Sciences, UCLA, Los Angeles, CA, USA Center for Space Plasma Physics, Space Science Institute, Boulder, CO, USA NASA Goddard Space Flight Center, Greenbelt, MD, USA Department of Physics Astronomy, University of Iowa, IA, USA
In this study, we examine particle energization and injections that show energetic electron enhancements at both MMS in the magnetotail and Van Allen Probes in the inner magnetosphere. Observing injections along with a corresponding flow burst allows us to better understand injections overall. Searching for suitable events, we found that only a small number of events at MMS had corresponding injections that penetrated far enough into the inner magnetosphere to observe with Van Allen Probes. With the four suitable events we did find, we compared the energy spectra at the two spacecraft and mapped the boundary of where the injection entered the inner magnetosphere. We found that, among these injections in the inner magnetosphere, the electron flux did not increase above ∼400 keV, similar to previous results, but the corresponding signatures in the tail observed increased fluxes at 600 keV or higher. There does not appear to be a comparable flux increase at Van Allen Probes and MMS for a given event. None of our injections included ion enhancements at Van Allen Probes, but one included an ion injection at geosynchronous orbit in the GOES spacecraft. All of our injections were dispersed at Van Allen Probes, and we were therefore able to map an estimate of the injection boundary. All of the injections occurred in the premidnight sector. Although we found some events where particle energizations in the tail are accompanied by inner magnetospheric injections, we do not find a statistical link between the two.
We present empirical conductance relations that are derived from incoherent scatter radar observations and correlated with all sky imager observations to identify the morphology of the aurora. We use 75,461 events collected using the Poker Flat Incoherent Scatter Radar (PFISR) with associated all sky imagers observations spanning the years 2012–2016. In addition to classifying these events based on auroral morphology, we estimated the Hall and Pedersen conductance and the differential number flux from which the energy flux and the average energy can be calculated. The differential number flux was estimated using the maximum entropy inversion method described in Semeter and Kamalabadi (2005, https://doi.org/10.1029/2004RS003042 ), but now incorporating the Fang et al. (2010, https://doi.org/10.1029/2010GL045406 ) ionization model. The main results of this investigation are the power law equations that describe the median, 90th, and 10th percentile Hall and Pedersen conductance as a function of energy flux and average energy. These power law fits are performed for different auroral morphology including all events, discrete, diffuse, and pulsating auroral events. The median Pedersen conductance is found to be in good agreement with past empirical conductance specifications by Robinson et al. (1987, https://doi.org/10.1029/JA092iA03p02565 ); however, the median Hall conductance from the PFISR observations is found to be larger than the empirical Hall conductance formulas by Robinson et al. (1987, https://doi.org/10.1029/JA092iA03p02565 ). Pulsating aurora is found to be the most frequently occurring auroral morphology. Furthermore, pulsating aurora has an important contribution to Hall conductance since it has higher average energies than discrete aurora. The results from this investigation are applicable to space weather models and may enable better agreement between model‐data comparisons.
In the next decade, there is an opportunity for very high return on investment of relatively small budgets by elevating the priority of smallsat funding in heliophysics. We've learned in the past decade that these missions perform exceptionally well by traditional metrics, e.g., papers/year/\$M (Spence et al. 2022 -- arXiv:2206.02968). It is also well established that there is a "leaky pipeline" resulting in too little diversity in leadership positions (see the National Academies Report at https://www.nationalacademies.org/our-work/increasing-diversity-in-the-leadership-of-competed-space-missions). Prioritizing smallsat funding would significantly increase the number of opportunities for new leaders to learn -- a crucial patch for the pipeline and an essential phase of career development. At present, however, there are far more proposers than the available funding can support, leading to selection ratios that can be as low as 6% -- in the bottom 0.5th percentile of selection ratios across the history of ROSES. Prioritizing SmallSat funding and substantially increasing that selection ratio are the fundamental recommendations being made by this white paper.
Chorus wave emissions are pervasive electromagnetic waves observed throughout the inner magnetosphere, and can intensify significantly during enhanced geomagnetic activity. Chorus waves have been shown to drive a particular type of precipitation that is composed of ring current and radiation belt losses and results in the quasi-periodic optical ionospheric phenomenon called pulsating aurora. Understanding how chorus waves evolve over time following periods of generation and rapid growth (i.e. injection-driven) and how they evolve as they propagate from equatorial sources to higher magnetic latitudes will lead to insights about how they drive energetic electron precipitation.
Electron cyclotron harmonic (ECH) waves along with whistler mode waves are suggested to be responsible for causing the persistent diffuse aurora in Jupiter. In this work, for the first time we systematically analyze the ECH waves in the Jovian inner magnetosphere, which was surveyed by Juno during the later orbits (>25). We find that in the Jovian inner magnetosphere, ECH waves occur in two specific regions—one equatorial and the other off-equatorial, just outside the Io torus. Equatorial ECH waves have higher intensity compared to their off-equatorial counterpart. We also notice an overlap between the region of mid-latitude hot injections and the region of off-equatorial ECH wave occurrence. Finally, we show an event to describe the complex nature of ECH wave growth/damping varying with particle density structures of the injection region at mid-latitude.
The accessibility of the activities and spaces that heliophysics occupies is connected to our ability to maintain a diverse and inclusive group of scientists.Accessibility is often overlooked in the discussions around promoting diversity, equity, and inclusion.There are very few statistics to show how accessible our field is, but larger studies of all STEM students show that the percentage of students who identify as disabled decreases with increasing education level, and that disabled students face a larger financial burden in their training compared to their peers.We recommend assessing the inclusion of those with disabilities in heliophysics in metrics related to the state of the profession.We also recommend promoting policies that make professional gatherings and educational settings more accessible to those with visible and invisible disabilities.