We investigate the dynamics of relativistic electrons in the Earth’s outer radiation belt by analyzing the interplay of several key physical processes: electron losses due to pitch angle scattering from electromagnetic ion cyclotron (EMIC) waves and chorus waves, and electron flux increases from chorus wave-driven acceleration of ~100-300 keV seed electrons injected from the plasma sheet. We examine a weak geomagnetic storm on April 17, 2021, using observations from various spacecraft, including GOES, Van Allen Probes, ERG/ARASE, MMS, ELFIN, and POES. Despite strong EMIC- and chorus wave-driven electron precipitation in the outer radiation belt, trapped 0.1-1.5 MeV electron fluxes actually increased. We use theoretical estimates of electron quasi-linear diffusion rates by chorus and EMIC waves, based on statistics of their wave power distribution, to examine the role of those waves in the observed relativistic electron flux variations. We find that a significant supply of 100-300 keV electrons by plasma sheet injections together with chorus wave-driven acceleration can overcome the rate of chorus and EMIC wave-driven electron losses through pitch angle scattering toward the loss cone, explaining the observed net increase in electron fluxes. Our study emphasizes the importance of simultaneously taking into account resonant wave-particle interactions and modeled local energy gradients of electron phase space density following injections, to accurately forecast the dynamical evolution of trapped electron fluxes.

We have identified for the first time an energy-time dispersion of precipitating electron flux in a pulsating aurora patch, ranging from 6.7 keV to 580 keV, through simultaneous in-situ observations of sub-relativistic electrons of microburst precipitations and lower-energy electrons using the LAMP sounding rocket launched from the Poker Flat Research Range in Alaska. Our observations reveal that precipitating electrons with energies of 180-320 keV were observed first, followed by 250-580 keV electrons 0-30 ms later, and finally, after 500-1000 ms, 6.7-14.6 keV electrons were observed. The identified energy-time dispersion is consistent with the theoretical estimation that the relativistic electron microbursts are a high-energy tail of pulsating aurora electrons, which are caused by chorus waves propagating along the field line.

The relativistic (~500 keV–2 MeV) and ultra-relativistic (~> 2 MeV) electron three-belt structure can originate from the partial depletion of the preexisting outer belt and replenishment of the new belt without direct electron injections. These processes may be related to radial diffusion caused by ULF (Ultra-Low Frequency) waves and local acceleration by VLF (Very-Low Frequency) waves. In this study, we reported a three-belt event with several hundred keV electrons based on Arase observation. The partial depletion of the outer belt and subsequent formation of a new belt due to source and seed electron (~30–500 keV) injection may be the primary mechanisms responsible for the three-belt structure. This discovery showed that the three-belt structure might not be limited to relativistic and ultra-relativistic electrons. Our study provided evidence for the existence of the source and seed electron three-belt structure, which could help improve our understanding of the radiation belt configuration.

A document by Mykhaylo Shumko. Click on the document to view its contents.

Electromagnetic ion cyclotron (EMIC) waves are believed to cause the loss of relativistic electrons from the outer radiation belt into the atmosphere due to pitch angle scattering. However, it is still unclear whether all EMIC waves can scatter relativistic electrons or which conditions are favorable for pitch angle scattering by EMIC waves. In this study, we performed a two-year data analysis of EMIC waves and relativistic electron precipitation (REP) caused by EMIC waves, from 1 November 2016 to 31 October 2018. EMIC waves were observed using a ground-based magnetometer installed at Athabasca (ATH), Canada. REP events were identified from very low-frequency radio waves propagated from the transmitters at the NDK and NLK stations (North Dakota and Seattle, USA, respectively) to the receiver installed at ATH. The magnetic local time dependence of EMIC waves showed higher occurrence rates in the dawn sector. In contrast, EMIC waves accompanied by REP were localized in the dusk sector and were likely to occur during geomagnetic substorms. We found that EMIC waves accompanied by REP were associated with the main phase of geomagnetic storms and occurred inside the plasmapause. These results suggest that the EMIC waves that cause REP occur in the overlap region between the ring current and dense cold plasma during the main phase of geomagnetic storms. This is consistent with previous studies describing that the electron resonant energy with EMIC waves is lower in regions with high plasma density.

We present X-ray observations of the upper atmospheric density disturbance caused by the explosive eruption of the Hunga Tonga-Hunga Ha’apai (HTHH) volcano on 15 January 2022. From 14 January to 16 January, the Chinese X-ray astronomy satellite, Insight-HXMT, was observing the supernova remnant Cassiopeia~A. The X-ray data obtained during Earth’s atmospheric occultations allowed us to measure neutral densities in the altitude range of $\sim$90–150\,km. The density profiles above 110\,km altitude obtained before the major eruption are in reasonable agreement with expectations by both GAIA and NRLMSIS~2.0 models. In contrast, after the HTHH eruption we found a severe density depletion at 500–1,000\,km away from the epicenter, and a relatively weak depletion extending up to $\sim7,000$\,km for over 8\,hr after the eruption. In addition, density profiles showed wavy structures with a typical length scale of either $\sim$20\,km (vertical) or $\sim$1,000\,km (horizontal). This may reflect gravity waves triggered by the volcanic eruption.

We analyzed time-of-flight (TOF) data from the Arase satellite to investigate temporal variations of O2+, NO+, and N2+ at 19.2 keV/q in the inner magnetosphere for 6.5 years from the solar declining to rising phases. Molecular ion counts were estimated by subtracting the background contamination of oxygen counts. While the number of clear molecular events was small, the estimated molecular ion counts exhibited good correlation with the solar wind dynamic pressure and SYM-H index. Long-term variations of molecular ions were different from that of oxygen ions. Additionally, we discuss the importance of the solar wind dynamic pressure in causing the escape of molecular ions into the magnetosphere through an increase in the convection electric field, which causes different evolutions of oxygen ions and molecular ions.

We report a citizen science-motivated study on the cause of an unusually bright red aurora as witnessed from Hokkaido, Japan during a magnetic storm on December 1, 2023. Such an intense red aurora event has occurred in the Halloween 2003 super storm, but the Dst index peak of this December 2023 storm was only -107 nT. In spite of the moderate storm amplitude, the extremely high solar wind density of >50 /cc caused the aurora oval extension to 53 MLAT (L=2.8). We discuss that the drift loss of the ring current particles across the small-size magnetopause is important, and Hokkaido was at the right position to see the direct effect of the large particle injection of the storm-time substorm.

Near-equatorial measurements of energetic electron fluxes, in combination with numerical simulation, are widely used for monitoring of the radiation belt dynamics. However, the long orbital periods of near-equatorial spacecraft constrain the cadence of observations to once per several hours or greater, i.e., much longer than the mesoscale injections and rapid local acceleration and losses of energetic electrons of interest. An alternative approach for radiation belt monitoring is to use measurements of low-altitude spacecraft, which cover, once per hour or faster, the latitudinal range of the entire radiation belt within a few minutes. Such an approach requires, however, a procedure for mapping the flux from low equatorial pitch angles (near the loss cone) as measured at low altitude, to high equatorial pitch angles (far from the loss cone), as necessitated by equatorial flux models. Here we do this using the high energy resolution ELFIN measurements of energetic electrons. Combining those with GPS measurements we develop a model for the electron anisotropy coefficient, n, that describes electron flux jtrap dependence on equatorial pitch-angle, αeq, jtrap ∼ sinnαeq. We then validate this model by comparing its equatorial predictions from ELFIN with in-situ near-equatorial measurements from Arase (ERG) in the outer radiation belt.

We performed a statistical study of electromagnetic ion cyclotron (EMIC) wave distributions and their coupling with energetic protons in the inner magnetosphere using the Arase satellite data from May 2017 to December 2020. We investigated the energetic proton pitch-angle distributions and partial thermal pressures associated with EMIC waves using inter-calibrated proton data in the energy range of 30 eV/q-187 keV/q. With a cold plasma approximation, we computed the proton minimum resonance energies using the observed EMIC wave frequency and plasma density values. We found that the EMIC waves had left-handed polarization near the magnetic equator close to the threshold of proton cyclotron instability, and propagated to higher latitudes along the field line with polarization reversal. H-EMIC waves showed two peak occurrence regions in the morning and noon sectors at L=7.5-9 outside the plasmasphere. The flux enhancements associated with morning side H-EMIC waves appeared at E<1 keV/q among all pitch angles, while H-EMIC waves in the noon sector exhibited flux enhancement in field-aligned directions at E=1-100 keV/q. He-EMIC waves showed a broad occurrence region from 12 to 20 magnetic local time at L=5.5-8.5 inside the plasmasphere with strong flux enhancements at all pitch-angle ranges at E=1-100 keV/q. The proton minimum resonance energy using the obtained central frequency was consistent with the observed flux enhancements at different peak occurrence regions. We conclude that the free energy sources of EMIC waves in different geomagnetic environments drive the two different types of EMIC waves, and they interact with energetic protons at different energy ranges.

Two-dimensional hybrid particle-in-cell (PIC) simulations are carried out on a constant L-shell (or drift shell) surface of the dipole magnetic field to investigate the generation process of near-equatorial fast magnetosonic waves (a.k.a equatorial noise; MSWs hereafter) in the inner magnetosphere. The simulation domain on a constant L-shell surface adopted here allows wave propagation and growth in the azimuthal direction (as well as along the field line) and is motivated by the observations that MSWs propagate preferentially in the azimuthal direction in the source region. Furthermore, the equatorial ring-like proton distribution used to drive MSWs in the present study is (realistically) weakly anisotropic. Consequently, the ring-like velocity distribution projected along the field line by Liouville’s theorem extends to rather high latitude, and linear instability analysis using the local plasma conditions predicts substantial MSW growth up to +- 27deg latitude. In the simulations, however, the MSW intensity maximizes near the equator and decreases quasi-exponentially with latitude. Further analysis reveals that the stronger equatorward refraction at higher latitude due to the larger gradient of the dipole magnetic field strength prevents off-equatorial MSWs from growing continuously, whereas MSWs of equatorial origin experience little refraction and can fully grow. Furthermore, the simulated MSWs exhibit a rather complex wave field structure varying with latitude, and the scattering of energetic ring-like protons in response to MSW excitation occurs faster than the bounce period of those protons so that they do not necessarily follow Liouville’s theorem during MSW excitation.

In this study, by simulating the wave-particle interactions, we show that sub-relativistic/relativistic electron microbursts form the high-energy tail of pulsating aurora (PsA). Whistler-mode chorus waves that propagate along the magnetic field lines at high latitudes cause precipitation bursts of electrons with a wide energy range from a few keV (PsA) to several MeV (relativistic microbursts). The rising tone elements of chorus waves cause individual microbursts of sub-relativistic/relativistic electrons and the internal modulation of PsA with a frequency of a few Hz. The chorus bursts for a few seconds cause the microburst trains of sub-relativistic/relativistic electrons and the main pulsations of PsA. Our simulation studies demonstrate that both PsA and relativistic electron microbursts originate simultaneously from pitch angle scattering by chorus wave-particle interactions.