Superfast precipitation of energetic electrons in the radiation belts of the Earth

The study investigates how energetic electron fluxes in Earth’s outer radiation belt precipitate into the atmosphere, challenging the prevailing quasi-linear diffusion framework. Traditionally, wave–particle interactions with whistler-mode waves are treated as diffusive, implying loss-cone fluxes cannot exceed trapped fluxes (strong diffusion limit). The authors pose the question: can energetic electron precipitation exceed this limit, and if so, by what mechanism? They highlight the importance of understanding precipitation for satellite risk, radiation belt dynamics, and atmospheric heating/chemistry. They propose that coherent, nonlinear interactions with intense, very oblique whistler-mode waves may drive faster-than-diffusive transport, potentially overfilling the loss cone and increasing precipitation beyond classical limits.

Prior work established the role of whistler-mode waves in accelerating and scattering radiation-belt electrons and framed precipitation using quasi-linear diffusion, with weak and strong diffusion limits for pitch-angle scattering. However, the random-phase assumption can break for intense waves, enabling nonlinear processes such as trapping and advection, which can transport electrons much faster in energy and pitch angle. Very oblique whistler-mode waves, with strong field-aligned electric fields and elliptical polarization, have been observed by missions such as THEMIS, Cluster, and Van Allen Probes, and are known to enable Landau resonance. Although their impact on precipitation has been discussed, they are often excluded from radiation belt models due to low average magnetic field amplitudes and presumed damping. The literature also contrasts Landau resonance with cyclotron resonance; the latter generally increases pitch angle for sub-relativistic energies and cannot readily explain loss-cone overfilling, whereas Landau resonance naturally couples pitch-angle decrease with energy gain along resonance curves, enabling efficient transport toward the loss cone.

The study combines multi-mission observations with numerical simulations and event/statistical analyses.

• Conjugate observations: Low-altitude ELFIN CubeSats (A/B) measured precipitating and trapped energetic electron fluxes (pitch-angle and energy resolved) at ~400–450 km altitude across latitudes, while THEMIS (probe E) at the equator provided magnetic field, plasma parameters, and whistler-mode wave characteristics. A conjunction on 23 Nov 2020 (~03:01–03:04 UT) at L≈7.5–10 and near-noon MLT enabled linking ELFIN precipitation to equatorial waves observed continuously by THEMIS from 02:00–04:00 UT (frequencies 0.1–0.3 fce). Conjunction criteria accounted for typical whistler source region correlation lengths (~1.4 h MLT, ~1.5 RE radially) and timing (~35 min).
• Wave characterization: THEMIS electric field spectra were converted to magnetic field spectra using cold plasma dispersion due to lack of recent search-coil calibration. Waveforms and hodograms showed strong field-aligned electric fields and elliptical polarization indicative of very oblique propagation (near resonance cone). Electric field amplitudes reached ~10 mV/m.
• Supporting near-equatorial measurements: ERG/Arase MEP-e provided jprec/jtrap in the equatorial region encompassing the loss cone, averaged over ~20 min to boost counting statistics, with trapped flux taken from 5–15° pitch-angle range.
• Numerical simulations: A probabilistic phase-space mapping technique modeled nonlinear Landau resonance with oblique whistler waves. Inputs included observed equatorial electron distributions (THEMIS ESA and SST, interpolated across energy gap), wave spectral power distributions P(B,ω) from THEMIS, and observed fpe/fce. Wave-normal angles θ were set uniformly within [θr−10°, θr−5°] near the resonance cone, consistent with waveform inferences and statistical studies. The model uses analytical expressions for trapping probability Π(E,α), trapping/scattering energy changes ΔEtrap(E,α), ΔEscat(E,α), and computes pitch-angle changes from magnetic moment conservation. An ensemble of 10^6 orbits was propagated for up to 10 minutes with observed bursty wave intensity profiles to map equatorial distributions to ELFIN altitudes.
• Statistics: Five months of ELFIN dayside observations (L<12, MLT 9–15) were surveyed to determine the occurrence of loss-cone overfilling events (jprec/jtrap>1) and their L-shell dependence; selected conjunctions with THEMIS verified association with very oblique whistler-mode waves.

• Direct observation of loss-cone overfilling: ELFIN measured precipitating electron fluxes exceeding trapped fluxes (jprec/jtrap>1) for energies below ~200 keV during 03:01–03:04 UT on 23 Nov 2020 at L≈7.5–10. Pitch-angle spectra showed higher fluxes within the local loss cone than just outside it for ~63, 98, and 139 keV.
• Causative waves: Conjugate THEMIS observations showed continuous strong whistler-mode waves with f≈0.1–0.3 fce, very oblique polarization, and strong field-aligned electric fields; wave packet amplitudes reached ~10 mV/m during 02:00–04:00 UT in the same L and MLT sector.
• Mechanism: Nonlinear Landau trapping by very oblique whistler-mode waves accelerates source electrons (~10–30 keV) to ~60–150 keV while rapidly decreasing pitch angle into the loss cone, producing superfast precipitation that exceeds the strong diffusion limit. Energy gain and pitch-angle decrease follow magnetic moment conservation (E sin^2αeq ≈ const), enabling rapid transport into the loss cone within a single resonant interaction.
• Simulation–observation agreement: Numerical mapping using observed wave and particle inputs reproduced the ELFIN pitch-angle distributions, including loss-cone overfilling at ~60–150 keV and its absence above ~300 keV. The model indicates trapping and release up to mid-latitudes (~30–40°) are needed to reach 60–140 keV at high latitudes.
• Equatorial confirmation: ERG/Arase measured jprec≥jtrap at ~10–20 keV within the loss cone near the equator during 02:30–03:00 UT, consistent with low-energy overfilling expected near the equator.
• Occurrence rates: From five months of ELFIN dayside data, loss-cone overfilling (jprec/jtrap>1) occurred in ~10% of precipitation events (those with jprec/jtrap>0.05) at L∈[10,12], ~5% at L∈[6,9], decreasing to near-zero at L<4; 166 orbits (of 465) and 943 spins showed overfilling.
• Atmospheric impact: ~100 keV precipitating electrons can reach ~75 km altitude, implying significant ionospheric/atmospheric effects; strong, sporadic losses can deplete 10–30 keV source populations, influencing radiation belt dynamics and seed availability for relativistic acceleration.

The findings demonstrate that the classical strong diffusion limit is not an upper bound for precipitating fluxes when coherent, intense, very oblique whistler-mode waves drive nonlinear Landau trapping. This mechanism rapidly moves electrons from moderate pitch angles directly into the loss cone while increasing their energy, explaining observed loss-cone overfilling at ~60–150 keV. The agreement between ELFIN observations, conjugate THEMIS wave measurements, ERG equatorial loss-cone observations, and nonlinear simulations supports this interpretation and rules out cyclotron resonance as the primary cause for the observed feature. The relatively frequent occurrence (5–10% regionally) indicates superfast precipitation can contribute substantially to outer belt electron losses, modulate seed populations, and enhance magnetosphere–atmosphere coupling, with potential implications for atmospheric chemistry and climate. The results argue for incorporating nonlinear interactions with oblique whistlers (and potentially time domain structures on the nightside) into radiation belt models and space weather nowcasts/forecasts.

This work provides direct evidence that energetic electron precipitation from the outer radiation belt can exceed the strong diffusion limit due to nonlinear Landau trapping by intense, very oblique whistler-mode waves. Conjugate spacecraft observations and dedicated nonlinear simulations reproduce loss-cone overfilling at ~60–150 keV and identify the source as accelerated ~10–30 keV electrons. Statistical analysis shows that such events are common enough to be consequential for radiation belt dynamics and atmospheric impacts. Future research should incorporate nonlinear, oblique-wave effects into global radiation belt models, extend occurrence statistics across local times (including nightside) and geomagnetic conditions, improve in-situ wave vector and amplitude characterization (including calibrated magnetic measurements), and pursue multi-point, loss-cone-resolving observations to quantify precipitation energetics and atmospheric consequences.

• Conjunction assumptions: Attribution of ELFIN precipitation to THEMIS-observed waves relies on spatial–temporal conjunction criteria and assumes similar wave properties within ~35 minutes and within ~1.4 h MLT and ~1.5 RE in L.
• Instrumental/data constraints: THEMIS magnetic spectra were inferred from electric field measurements via cold plasma dispersion because search-coil calibration was not available for the interval; wave-normal angles were estimated from electric field waveforms. ERG equatorial loss-cone measurements required time averaging, with poor counting statistics at 12–14 keV.
• Observational coverage: ELFIN statistical analysis focuses on dayside, L<12, MLT 9–15; occurrence rates may differ at other local times (e.g., nightside) and geomagnetic conditions. Equatorial loss-cone measurements are generally challenging due to small loss-cone angles, limiting direct equatorial confirmation to certain energy ranges (<~20 keV).
• Modeling assumptions: The mapping technique employs probabilistic ensembles of observed wave characteristics and assumes phase randomization between resonances. Background density and wave-normal angle profiles are modeled and may introduce uncertainties. Nonlinear interaction characteristics cannot be simply averaged, requiring distributional assumptions that may affect quantitative predictions.