Relativistic microburst storm characteristics: Combined satellite and ground-based observations

We report a comparison of Solar Anomalous Magnetospheric Particle Explorer detected relativistic electron microbursts and short-lived subionospheric VLF perturbations termed FAST events, observed at Sodankyl Geophysical Observatory, Finland, during 2005. We show that only strong geomagnetic disturbances can produce FAST events, which is consistent with the strong link between storms and relativistic microbursts. Further, the observed FAST event perturbation decay times were consistent with ionospheric recovery from bursts of relativistic electron precipitation. However, the one-to-one correlation in time between microbursts and FAST events was found to be very low (similar to 1%). We interpret this as confirmation that microbursts have small ionospheric footprints and estimate the individual precipitation events to be <4 km radius. In contrast, our study strongly suggests that the region over which microbursts occur during storm event periods can be at least similar to 90 degrees in longitude (similar to 6 h in magnetic local time). This confirms earlier estimates of microburst storm size, suggesting that microbursts could be a significant loss mechanism for radiation belt relativistic electrons during geomagnetic storms. Although microbursts are observed at a much higher rate than FAST events, the ground-based FAST event data can provide additional insight into the conditions required for microburst generation and the time variation of relativistic precipitation

Initial measurements of O-ion and He-ion decay rates observed from the Van Allen probes RBSPICE instrument.

H-ion (∼45 keV to ∼600 keV), He-ion (∼65 keV to ∼520 keV), and O-ion (∼140 keV to ∼1130 keV) integral flux measurements, from the Radiation Belt Storm Probe Ion Composition Experiment (RBSPICE) instrument aboard the Van Allan Probes spacecraft B, are reported. These abundance data form a cohesive picture of ring current ions during the first 9 months of measurements. Furthermore, the data presented herein are used to show injection characteristics via the He-ion/H-ion abundance ratio and the O-ion/H-ion abundance ratio. Of unique interest to ring current dynamics are the spatial-temporal decay characteristics of the two injected populations. We observe that He-ions decay more quickly at lower L shells, on the order of ∼0.8 day at L shells of 3-4, and decay more slowly with higher L shell, on the order of ∼1.7 days at L shells of 5-6. Conversely, O-ions decay very rapidly (∼1.5 h) across all L shells. The He-ion decay time are consistent with previously measured and calculated lifetimes associated with charge exchange. The O-ion decay time is much faster than predicted and is attributed to the inclusion of higher-energy (&gt; 500 keV) O-ions in our decay rate estimation. We note that these measurements demonstrate a compelling need for calculation of high-energy O-ion loss rates, which have not been adequately studied in the literature to date.Key pointsWe report initial observations of ring current ionsWe show that He-ion decay rates are consistent with theoryWe show that O-ions with energies greater than 500 keV decay very rapidly

Machine Learning Interpretability of Outer Radiation Belt Enhancement \& Depletion Events

We investigate the response of outer radiation belt electron fluxes to different solar wind and geomagnetic indices using an interpretable machine learning method. We reconstruct the electron flux variation during 19 enhancement and 7 depletion events and demonstrate a feature attribution analysis on the superposed epoch results for the first time. We find that the intensity and duration of the substorm sequence following an initial dropout determine the overall enhancement or depletion of electron fluxes, while the solar wind pressure drives the initial dropout in both types of events. Further statistical results from a dataset with 71 events confirm this and show a significant correlation between the resulting flux levels and the average AL index, indicating that the observed "depletion" event can be more accurately described as a "non-enhancement" event. Our novel SHAP-Enhanced Superposed Epoch Analysis (SHESEA) method can be used as an insight discovery tool in various physical systems

Formation of Foreshock Transients and Associated Secondary Shocks

Upstream of shocks, the foreshock is filled with hot ions. When these ions are concentrated and thermalized around a discontinuity, a diamagnetic cavity bounded by compressional boundaries, referred to as a foreshock transient, forms. Sometimes, the upstream compressional boundary can further steepen into a secondary shock, which has been observed to accelerate particles and contribute to the primary shock acceleration. However, secondary shock formation conditions and processes are not fully understood. Using particle-in-cell simulations, we reveal how secondary shocks are formed. From 1D simulations, we show that electric fields play a critical role in shaping the shock's magnetic field structure, as well as in coupling the energy of hot ions to that of the shock. We demonstrate that larger thermal speed and concentration ratio of hot ions favor the formation of a secondary shock. From a more realistic 2D simulation, we examine how a discontinuity interacts with foreshock ions leading to the formation of a foreshock transient and a secondary shock. Our results imply that secondary shocks are more likely to occur at primary shocks with higher Mach number. With the secondary shock's previously proven ability to accelerate particles in cooperation with a planetary bow shock, it is even more appealing to consider them in particle acceleration of high Mach number astrophysical shocks.Peer reviewe

Nonresonant scattering of energetic electrons by electromagnetic ion cyclotron waves: spacecraft observations and theoretical framework

Electromagnetic ion cyclotron (EMIC) waves lead to rapid scattering of relativistic electrons in Earth's radiation belts, due to their large amplitudes relative to other waves that interact with electrons of this energy range. A central feature of electron precipitation driven by EMIC waves is deeply elusive: moderate precipitating fluxes at energies below the minimum resonance energy of EMIC waves occur concurrently with strong precipitating fluxes at resonance energies in low-altitude spacecraft observations. Here we expand on a previously reported solution to this problem: nonresonant scattering due to wave packets of finite size. We first generalize the quasi-linear diffusion model to incorporate nonresonant scattering by a generic wave shape. The diffusion rate decays exponentially away from the resonance, where shorter packets lower decay rates and thus widen the energy range of significant scattering. Using realistic EMIC wave packets from $\delta f$ particle-in-cell simulations, we then perform test particle simulations, and demonstrate that intense, short packets extend the energy of significant scattering well below the minimum resonance energy, consistent with our theoretical prediction. Finally, we compare the calculated precipitating-to-trapped flux ratio of relativistic electrons to ELFIN observations, and infer the wave power spectra that are consistent with the measured flux ratio. We demonstrate that even with a narrow wave spectrum, short EMIC wave packets can provide moderately intense precipitating fluxes well below the minimum resonance energy.Comment: 27 pages, 7 figure

Nonlinear Landau resonant interaction between whistler waves and electrons: Excitation of electron acoustic waves

Electron acoustic waves (EAWs), as well as electron-acoustic solitary structures, play a crucial role in thermalization and acceleration of electron populations in Earth's magnetosphere. These waves are often observed in association with whistler-mode waves, but the detailed mechanism of EAW and whistler wave coupling is not yet revealed. We investigate the excitation mechanism of EAWs and their potential relation to whistler waves using particle-in-cell simulations. Whistler waves are first excited by electrons with a temperature anisotropy perpendicular to the background magnetic field. Electrons trapped by these whistler waves through nonlinear Landau resonance form localized field-aligned beams, which subsequently excite EAWs. By comparing the growth rate of EAWs and the phase mixing rate of trapped electron beams, we obtain the critical condition for EAW excitation, which is consistent with our simulation results across a wide region in parameter space. These results are expected to be useful in the interpretation of concurrent observations of whistler-mode waves and nonlinear solitary structures, and may also have important implications for investigation of cross-scale energy transfer in the near-Earth space environment

Opening the Black Box of the Radiation Belt Machine Learning Model

Many Machine Learning (ML) systems, especially neural networks, are fundamentally regarded as black boxes since it is difficult to grasp how they function once they have been trained. Here, we tackle the issue of the interpretability of a high-accuracy ML model created to model the flux of Earth's radiation belt electrons. The Outer RadIation belt Electron Neural net model (ORIENT) uses only solar wind conditions and geomagnetic indices as input. Using the Deep SHAPley additive explanations (DeepSHAP) method, we show that the `black box' ORIENT model can be successfully explained. Two significant electron flux enhancement events observed by Van Allen Probes during the storm interval of 17 to 18 March 2013 and non storm interval of 19 to 20 September 2013 are investigated using the DeepSHAP method. The results show that the feature importances calculated from the purely data driven ORIENT model identify physically meaningful behavior consistent with current physical understanding.Comment: Under revie