The dipole configuration of the Earth’s magnetic field allows for the trapping of highly energetic particles, which form the radiation belts. Although significant advances have been made in understanding the acceleration mechanisms in the radiation belts, the loss processes remain poorly understood. Unique observations on 17 January 2013 provide detailed information throughout the belts on the energy spectrum and pitch angle (angle between the velocity of a particle and the magnetic field) distribution of electrons up to ultra-relativistic energies. Here we show that although relativistic electrons are enhanced, ultra-relativistic electrons become depleted and distributions of particles show very clear telltale signatures of electromagnetic ion cyclotron wave-induced loss. Comparisons between observations and modelling of the evolution of the electron flux and pitch angle show that electromagnetic ion cyclotron waves provide the dominant loss mechanism at ultra-relativistic energies and produce a profound dropout of the ultra-relativistic radiation belt fluxes.
The dynamic evolution of the Earth’s electron radiation belts is the result of the competition between acceleration and loss processes. Recent observations and modelling provided significant advances in our understanding of acceleration mechanisms operating in the Earth radiation belts1,2,3, whereas the loss mechanisms remain more controversial. Understanding of acceleration provides only ‘half’ of the global picture required for understanding the evolution of electrons, as the dynamics of the radiation belts are not determined by just the acceleration mechanisms, but rather the dynamic battle between acceleration and loss.
The earliest studies of radiation belt dynamics noted that electron fluxes tend to decrease during the main phase of geomagnetic storms4 and these flux dropouts were attributed to reversible processes5. That is, when slow changes occur in the Earth’s magnetic field configuration during geomagnetic storms, electron fluxes are redistributed in radial distance, energy and pitch angle (the angle between the electron velocity vector and the magnetic field). When the magnetic field relaxes back to the pre-disturbance configuration, the electron distribution also reverses to the original state. Therefore, no net loss of electrons occurs. Recent studies6,7,8 showed that non-reversible electron dropouts occur during the main phase of many storms. One plausible theory to explain the non-reversible dropout is scattering by electromagnetic ion cyclotron (EMIC) waves6,9,10, which are often strongly enhanced during geomagnetic storms. After scattering by the waves, the electrons are lost as a result of precipitation into the upper atmosphere. Another mechanism to explain the dropouts has been proposed11. In this mechanism, compression of the outer boundary of the Earth’s magnetosphere during geomagnetic storms causes electrons at large radial distances to be lost to the interplanetary medium. Such loss creates steep radial gradients in electron-phase space density (PSD) that drive outward radial diffusion. The outward radial diffusion transports particles away from the Earth into the region of lower magnetic field and decelerates particles, which in turn propagates the electron loss to lower radial distance.
To differentiate between the two non-reversible processes, it is necessary to examine changes in the electron energy spectrum and pitch angle distribution. Numerical modelling has shown that dropouts observed down to low energy (100’s of kiloelectron volts) cannot be produced by EMIC waves, but are consistent with expectations from loss produced by outward radial diffusion11. Additional strong evidence for the outward diffusion mechanism came from observations of dropouts at energies of about 500 eV that were not accompanied by precipitating electrons, eliminating EMIC wave scattering as a driver of the observed loss12. However, it remained unclear whether higher energy electrons may still be predominantly scattered by EMIC waves.
A recent study has shown that EMIC waves are required to reproduce an unusually narrow remnant belt at energies >4 MeV13 that was observed on the Van Allen Probes and referred to as the ‘storage ring’14. Observations of electron pitch angle distributions also showed telltale signatures of interactions of ultra-relativistic electrons with EMIC waves15 but did not provide observational evidence that EMIC waves can initialize the net loss at all pitch angles and not just change the shape of the pitch angle distribution.
Identifying and separating the effects of different acceleration and loss mechanisms is often a challenging task. Multiple competing acceleration and loss mechanisms usually occur simultaneously. Another complication is that acceleration and loss mechanisms also depend on the seed population of electrons at somewhat lower energy than relativistic or ultra-relativistic energies. At times when seed population fluxes are increased, all acceleration and loss mechanisms will become intensified and distinguishing between them is challenging.
Careful selection of particular conditions during specific geomagnetic storms can help isolate different acceleration and loss mechanisms, allowing us to determine the dominant processes and ultimately provide insight into the physics driving radiation belt acceleration and loss.
Up until recently, there were no accurate and reliable measurements of ultra-relativistic electron populations with full pitch angle resolution. The relativistic electron–proton telescope (REPT) instrument on the Van Allen Probes spacecraft mission16 has provided very detailed information on the pitch angle distributions and the orbit of the spacecraft allows for the measurement of electrons trapped near the geomagnetic equator. However, for a number of storms, the flux of >3 MeV electrons was below the noise level of the high-energy channels of the REPT instrument, and energy and pitch angle distributions could not be inferred for these energies.
Evidence for EMIC wave scattering can come from observations of the electron energy spectrum, which only produces loss of electrons above a certain minimum energy as only very energetic electrons can be in resonance with these waves. Another typical signature of EMIC scattering is the narrowing of the normalized electron pitch angle distribution toward 90° (perpendicular to the local magnetic field), as only electrons with a large parallel energy and consequently small pitch angles can be in resonance with the waves.
Here we report observations of electron flux at different energies and observations of electron pitch angle distributions during 17 January 2013 storm, which are consistent with EMIC wave-induced loss for ultra-relativistic electrons, whereas electrons at relativistic energies are not affected by EMIC wave scattering. Comparison of modelling results to the observations provides further support of our conclusions. The presented observations and modelling show that the ultra-relativistic electron population experiences different loss mechanisms than the relativistic population.
Observations during the 17 January 2013 storm
The environment close to the Earth (radial distances <3–6 Earth radii (RE)) is usually occupied by the cold plasma bubble co-rotating with the Earth, which is referred to as the plasmasphere. Inside the plasmasphere, relativistic electrons usually decay on the scale of 1–10 days due to scattering by whistler-mode plasma waves known as plasmaspheric hiss. Ultra-relativistic electrons can persist for a very long time and scattering by plasmaspheric hiss will be much slower at these energies13.
The October 2012 geomagnetic storm2,3, which occurred 101 days before our considered event, produced an abundance of electrons at energies of Mega-electron Volt (MeV) and multi-MeV. Elevated fluxes of these most energetic electrons in the heart of the outer belt at levels above REPT instruments’ noise floor allowed for fully resolved measurements of pitch angle distributions as a function of energy. After the October storm, the ultra-relativistic (above ∼4 MeV) remnant belt slowly diffused inward to L between ∼3 and 4 (L is the distance from the centre of the Earth to a given magnetic field line in the equatorial plane measured in RE), where the ultra-relativistic electrons cannot be significantly affected by the loss to the magnetopause and outward radial diffusion, and scattering by hiss inside the plasmasphere is weak.
During a moderate storm on 17 January (Fig. 1a,b shows the index of geomagnetic activity Kp inferred from the fluctuations of magnetic field on the ground), flux of relativistic electrons at 1.02 MeV increased (Fig. 1c). The combination of local acceleration and inward radial diffusion moved the inner boundary of the outer belt to lower radial distances (Fig. 1c). The short-lived dropout of fluxes is associated with the reversible changes described in the introduction. As clearly seen in the observations after the adiabatic dropout (see Supplementary Fig. 1 and Supplementary Note 1), fluxes return to pre-storm values.
Ultra-relativistic electrons at 4.2 MeV (Fig. 1d) show a very different evolution, which is a net decrease in flux. This loss of ultra-relativistic electrons is not produced by the loss to the magnetopause, as the multi-MeV electron belt is located deep inside the outer zone, below 4.5 RE, whereas the magnetopause for this event was compressed down to 7.1 RE according to an empirical model17(see also Supplementary Note 2 and Supplementary Figs 2 and 3) and the variation of the global magnetic field in the inner magnetosphere was not significantly large (see Supplementary Note 3 and Supplementary Fig. 4). Such difference in the evolution of fluxes between 1.02 and 4.2 MeV can be explained by the presence of scattering by EMIC waves that affects only electrons above certain minimum threshold energy, which is the minimum energy for which electrons can resonantly interact with EMIC waves.
Another clear piece of evidence for the EMIC wave scattering is provided by observation of the electron pitch angle distributions (Fig. 1e,f). During resonant interactions with EMIC waves, only small pitch angle electrons are scattered into the atmosphere and the distribution should have characteristic signatures with bite-outs near 0° and 180°, corresponding to particles with very high velocity parallel to the magnetic field that were scattered by EMIC waves. Although the pitch angle distribution of the 1.02 MeV electrons is rather wide (Fig. 1e), the ultra-relativistic electron distribution dramatically changes after the flux drop out and becomes very narrow, with sharp gradients around 50° and 130°.
EMIC waves are extremely strong plasma wave emissions that are often observed only in a very narrow localized region of magnetic local time (MLT). Although satellite measurements can only provide information at a particular radial distance and MLT, ground-based arrays of observations of EMIC waves provide a global view of waves.
Observations from an array of stations in Finland show strong EMIC emissions on 17 January, between ∼15 and 18 universal time (see Fig. 2), with a peak intensity in the heart of the outer zone, exactly at the same time when the Van Allen Probes observed a dropout in ultra-relativistic fluxes (Fig. 1d) and prompt narrowing of the pitch angle distributions (Fig. 1f). The clear signature of narrowing of the distributions during the days when EMIC wave activity was seen on the ground clearly show that the pitch angle distributions are produced by resonant scattering with EMIC waves.
Comparison of modelling with observations
To model the evolution of fluxes, we use numerical simulations that include diffusion in radial distance, energy and pitch angle. The numerical code is described in detail in the Methods section. Simulations of relativistic electron radiation belt dynamics including radial transport, and local acceleration and loss due to whistler-mode waves (Fig. 3a) can well reproduce the dynamics of the 1.02 MeV electrons (Fig. 3e). Simulation of 4.2 MeV electrons with the same numerical code also predicts increase in fluxes (Fig. 3b), contrary to observations that clearly show the dropout as discussed above and shown for comparison in Fig. 3f. The introduction of EMIC waves in the simulation during 15–18 universal time on 17 January does not significantly change the dynamics of the 1.02 MeV electrons (Fig. 3c) but dramatically changes the evolution of ultra-relativistic electrons. Inclusion of EMIC waves at the time when they are observed on the ground produces a dropout in fluxes at 4.2 MeV (Fig. 3d), consistent with the evolution of radial profiles of fluxes observed on the Van Allen Probes (Fig. 3f). Results of the simulations for energy channels of 0.46 and 3.4 MeV electrons also agree well with observations and are presented in Supplementary Fig. 5.
At 1.02 MeV, evolution of the electron pitch angle distribution not including (Fig. 4a) and including (Fig. 4c) the effects of EMIC waves are similar, and both simulations produce a wide pitch angle distribution similar to the pitch angle distribution observed on Van Allen Probes, with the exception of the adiabatic changes (see also Supplementary Notes 1 and 2, and Supplementary Figs 2 and 3). The modelled distributions are rather wide. For the 4.2 MeV electrons, similar wide pitch angle distributions are seen when EMIC waves are not taken into account (Fig. 4b). However, the inclusion of EMIC waves (Fig. 4d) produces bite-outs at small pitch angles (near 0° and 180°) and particles are transported from high pitch angles (near 90°) towards the loss cone by whistler-mode waves. The scale of the dropout at L=3.9 RE and the narrow shape of the pitch angle distribution observed on Van Allen Probes (Fig. 4f) is well reproduced by the simulations with EMIC waves (Fig. 4d). The modelled pitch angle distributions at energies of 0.46 and 3.4 MeV also agree with observations and are presented in Supplementary Fig. 6.
The 17 January 2013 geomagnetic storm provided unique conditions that allowed us to show that EMIC-induced wave scattering can produce a very fast loss into the atmosphere at ultra-relativistic energy. Clear differences were observed in the dynamics at different energies, that is, the 1.02 MeV electrons are accelerated, while EMIC wave-induced precipitation depletes the 4.2 MeV electron population. The narrowing of the pitch angle distribution at 4.2 MeV is consistent with EMIC wave scattering, providing additional confirmation of the EMIC wave-induced loss.
Simulations including EMIC waves can reproduce both an increase in flux near 1 MeV energy and a reduction at 4.2 MeV. During this event, an abundance of ultra-relativistic electrons near L=3.5 allowed us to observe the progression of the shape of the pitch angle distribution. Pitch angle distributions provide additional evidence that the loss processes in the radiation belts are different for relativistic and ultra-relativistic energies. Ultra-relativistic electrons at moderate pitch angles have sufficient parallel velocity to change the sense of rotation of the left-hand polarized EMIC waves. In the rest reference frame EMIC waves propagating along the field line have left-hand polarization, whereas in the reference frame of the ultra-relativistic particles the wave becomes right-hand polarized, which is the same as the rotation sense of the electrons. Electrons gyrating around the field line in the same rotation sense and with the same frequency as waves in the electron reference frame can be in resonance with waves and can be effectively scattered.
The presence of EMIC waves provides very efficient loss mechanism for the most energetic electrons in the energy spectrum and may impose an upper limit on the energy of the trapped electron population. Similar processes may occur on the outer planets and exoplanets. The absence of EMIC waves and/or a relatively low, cold plasma density (due to lack of a dense planetary atmosphere) would imply that electrons on exoplanets may be potentially accelerated to extremely high energies, much higher than that of in the Earth’s magnetosphere, without being scattered into the loss cone.
The Versatile Electron Radiation Belt18,19 code was used to perform the simulations described in this manuscript. The Fokker–Plank equation can be written on the single grid using adiabatic invariant coordinates20:
where f is the electron PSD, L* is a form of the third invariant; , where J is the second adiabatic invariant, m0 is electron rest mass, , where μ is the first adiabatic invariant; and are bounce-averaged diffusion coefficients; ; is the Jacobean of the transformation from an adiabatic invariant system (μ, J, Φ) to (L*, V, K); Φ is the third adiabatic invariant, B0 is the magnetic field on the surface of the Earth, RE is the Earth’s radius; and τ is the electron’s lifetime inside the atmospheric loss cone, which is taken to be a quarter of the bounce period. To compute radial transport, we used activity-dependent21 radial diffusion rates. The flux of particles observed on satellites can be inferred from PSD using equation J=f·p2,where j is the differential directional flux and p is the relativistic momentum.
Wave parameters of hiss and chorus waves are based on statistical studies and were used in the previously published long-term simulations22,23. Wave parameters were not specifically adjusted for this event. The simulation includes whistler-mode chorus waves on the day and night side following the Kp-based parametrizations presented in ref. 23. The chorus waves were included above the plasmapause location. Plasmaspheric hiss waves were included inside the plasmasphere based on a recent empirical model24.
The plasmapause location was constructed based on the plasmapause test particle (PTP) simulation model24. In the Versatile Electron Radiation Belt code simulation, the location of the inner boundary of EMIC waves in the plasmaspheric plume and hiss in plume was calculated as the hourly average minimum of the plasmapause location calculated with the PTP model.
We assumed the presence of the plasmaspheric plume base on the PTP model25. The model shows clear presence of the plume during the interval from 17 January 2013, 15:00 h, to 19 January 2013, 08:00 h. The orbits of the Van Allen Probes did not cross the plume at the times when waves are observed at the ground stations in Fig. 2. However, the Van Allen Probes provide measurements of particles that drift around the Earth and are affected by waves in all MLT sectors. As discussed above, particle measurements provide clear evidence of loss by EMIC waves at ultra-relativistic energies.
The spatial distribution of the EMIC waves was assumed to be very localized at only 5% of the electron orbit. It was assumed that EMIC waves can effectively scatter electrons only in the high-density region of the plasmaspheric plume and were included in the model only above the PTP-calculated minimum plasmapause location.
The spectral parameters of the EMIC waves were obtained from the GOES 15 spacecraft measurements. Figure 5 shows the observation of EMIC waves on 17 January 2013 and the observed spectrum. We used the superposition of three Gaussian shapes to fit the observed spectrum for the further calculations of the EMIC wave diffusion coefficients (see Supplementary Table 1). The parameters that are assumed for the calculation of the EMIC wave-induced diffusion coefficients are based on GOES measurements and are presented and discussed in the Supplementary Note 5.
For this period, plasmaspheric hiss waves in the region of the plasmaspheric plume were also included above the calculated location of the plasmapause, assuming that they affect only 5% of drift time. Corresponding diffusion coefficients were scaled by Kp index.
The Fokker–Plank equation was solved using a fully implicit scheme on a 18 × 101 × 100 point grid for the coordinates L*, V and K, respectively, with a time resolution of 1 h. The 7 days of simulation from 14 January 2013 until 21 January 2013 was initialed with PSD profiles that were inferred from observation. Detailed information on how the data were processed is given in Supplementary Note 6.
The grid is uniform in L* and K, but logarithmic in V to resolve both lower energy and ultra-relativistic electrons. The boundaries for L* were set up at 2.1 and 5.5, with a step of 0.2 RE. The V and K coordinates were calculated in the range of energy and pitch angle defined at the upper radial boundary. The K coordinate was calculated on a linear grid of pitch angles from 0.7° to 89.3°. At the outer radial boundary of the code, the V coordinate was calculated on the logarithmic grid using the range of energies from 10 keV to 10 MeV.
The upper boundary conditions for L* were reconstructed from the Van Allen Probes observations using MagEIS and REPT measurements. The lower radial boundary, the upper boundary condition for V (high-energy population) and the upper boundary condition for K (electrons inside the loss cone) were set to zero, to simulate the absence of the electrons on the edge of the grid (see refs 21, 23 for more details). Practically, identical results were obtained if a zero gradient boundary condition is assumed at the higher K boundary representing a zero gradient in the pitch angle distribution at 0°. Constant PSD was prescribed at the lower boundary condition for V representing a balance between convective sources and losses. Although these conditions may not adequately represent the dynamics of the lower-energy electrons, the exact values of fluxes at this boundary does not substantially influence the results. Detailed sensitivity simulations to the assumed lower-energy boundary have been previously presented26.
All the Van Allen Probes data are publicly available at http://www.rbsp-ect.lanl.gov/data_pub/ by the REPT and MagEIS instrument. The GOES data are publicly available at http://satdat.ngdc.noaa.gov/sem/goes/data/. The OMNI data are obtained from http://omniweb.gsfc.nasa.gov/form/dx1.html. The modelling results of the plasmasphere are available on RBSP-ECT Data Portal (http://www.rbsp-ect.lanl.gov/). Kp index of geomagnetic activity was obtained from the GSFC/SPDF OMNIWeb interface at http://omniweb.gsfc.nasa.gov and produced by GFZ, Potsdam.The results of the model simulations and any other relevant data or information that may be necessary to reproduce the presented results is available from the authors upon request.
How to cite this article: Shprits, Y. Y. et al. Wave-induced loss of ultra-relativistic electrons in the Van Allen radiation belts. Nat. Commun. 7:12883 doi: 10.1038/ncomms12883 (2016).
Horne, R. B. et al. Wave acceleration of electrons in the Van Allen radiation belts. Nature 437, 227–230 (2005).
Reeves, G. D. et al. Electron acceleration in the heart of the Van Allen radiation belts. Science 341, 991–994 (2013).
Thorne, R. M. et al. Rapid local acceleration of relativistic radiation-belt electrons by magnetospheric chorus. Nature 504, 411–414 (2013).
Dessler, A. J. & Karplus, R. Some properties of the Van Allen radiation. Phys. Rev. Lett. 4, 271–274 (1960).
Kim, H.-J. & Chan, A. A. Fully adiabatic changes in storm time relativistic electron fluxes. J. Geophys. Res. 102, 22107 (1997).
Green, J. C. & Kivelson, M. G. Relativistic electrons in the outer radiation belt: differentiating between acceleration mechanisms. J. Geophys. Res. 109, A03213 (2004).
Onsager, T. G. et al. Radiation belt electron flux dropouts: local time, radial, and particle-energy dependence. J. Geophys. Res. [Space Phys.] 107, SMP 21–S1 21–11 (2002).
Shprits, Y., Daae, M. & Ni, B. Statistical analysis of phase space density buildups and dropouts. J. Geophys. Res. 117, A01219 (2012).
Thorne, R. M. & Kennel, C. F. Relativistic electron precipitation during magnetic storm main phase. J. Geophys. Res. 76, 4446–4453 (1971).
Sandanger, M., Søraas, F., Aarsnes, K., Oksavik, K. & Evans, D. S. Loss of relativistic electrons: evidence for pitch angle scattering by electromagnetic ion cyclotron waves excited by unstable ring current protons. J. Geophys. Res. 112, A12213 (2007).
Shprits, Y. Y. et al. Outward radial diffusion driven by losses at magnetopause. J. Geophys. Res. [Space Phys.] 111, A11214 (2006).
Turner, D. L., Shprits, Y., Hartinger, M. & Angelopoulos, V. Explaining sudden losses of outer radiation belt electrons during geomagnetic storms. Nat. Phys. 8, 208–212 (2012).
Shprits, Y. Y. et al. Unusual stable trapping of the ultrarelativistic electrons in the Van Allen radiation belts. Nat. Phys. 9, 699–703 (2013).
Baker, D. N. et al. A long-lived relativistic electron storage ring embedded in Earth’s outer Van Allen belt. Science 340, 186–190 (2013).
Usanova, M. E. et al. Effect of EMIC waves on relativistic and ultrarelativistic electron populations: ground-based and Van Allen Probes observations. Geophys. Res. Lett. 41, 1375–1381 (2014).
Baker, D. N. et al. The relativistic electron-proton telescope (REPT) instrument on board the radiation belt storm probes (RBSP) spacecraft: characterization of Earth’s radiation belt high-energy particle populations. Space Sci. Rev. 179, 337–381 (2013).
Shue, J.-H. et al. A new functional form to study the solar wind control of the magnetopause size and shape. J. Geophys. Res. 102, 9497–9511 (1997).
Shprits, Y. Y., Subbotin, D. A., Meredith, N. P. & Elkington, S. R. Review of modeling of losses and sources of relativistic electrons in the outer radiation belt II: local acceleration and loss. J. Atmos. Sol. Terr. Phys. 70, 1694–1713 (2008).
Subbotin, D. A. & Shprits, Y. Y. Three-dimensional modeling of the radiation belts using the Versatile Electron Radiation Belt (VERB) code. Space Weather 7, S10001 (2009).
Subbotin, D. A. & Shprits, Y. Y. Three-dimensional radiation belt simulations in terms of adiabatic invariants using a single numerical grid. J. Geophys. Res. 117, A05205 (2012).
Brautigam, D. H. & Albert, J. M. Radial diffusion analysis of outer radiation belt electrons during the October 9, 1990, magnetic storm. J. Geophys. Res. 105, 291–309 (2000).
Drozdov, A. Y. et al. Energetic, relativistic, and ultrarelativistic electrons: comparison of long-term VERB code simulations with Van Allen probes measurements. J. Geophys. Res. [Space Phys.] 120, 3574–3587 (2015).
Subbotin, D. A., Shprits, Y. Y. & Ni, B. Long-term radiation belt simulation with the VERB 3-D code: comparison with CRRES observations. J. Geophys. Res. [Space Phys.] 116, A12210 (2011).
Orlova, K., Shprits, Y. & Spasojevic, M. New global loss model of energetic and relativistic electrons based on Van Allen Probes measurements. J. Geophys. Res. [Space Phys.] 121, 1308–1314 (2016).
Goldstein, J. et al. Simulation of Van Allen Probes plasmapause encounters. J. Geophys. Res. [Space Phys.] 119, 7464–7484 (2014).
Subbotin, D. A. et al. Simulation of the acceleration of relativistic electrons in the inner magnetosphere using RCM-VERB coupled codes. J. Geophys. Res. 116, A08211 (2011).
Sakaguchi, et al. Simultaneous ground and satellite observations of an isolated proton arc at subauroral latitudes. J. Geophys. Res. 112, A04202 (2007).
This research was supported by the Helmholtz Association Recruiting Imitative programme, NASA grants NNX16AF91G,NNX13AE34G, NNX10AK99G, NNX15AI94G and NNX16AG78G, RBSP-ECT funding through JHU/APL contract 967399 under NASA contract NAS5-01072, JHU/APL contract 922613 (RBSP-EFW), NASA Grant NNX16AF85G, NSF GEM AGS-1203747 and UC Lab Fee grant 12-LR-235337, project PROGRESS funded by EU Horizon 2020 No 637302 and International Space Science Institute (Bern).
The authors declare no competing financial interests.
About this article
Cite this article
Shprits, Y., Drozdov, A., Spasojevic, M. et al. Wave-induced loss of ultra-relativistic electrons in the Van Allen radiation belts. Nat Commun 7, 12883 (2016). https://doi.org/10.1038/ncomms12883
Journal of Geophysical Research: Space Physics (2021)
Journal of Geophysical Research: Space Physics (2021)
Space Weather (2021)
Energetic Electron Precipitation Observed by FIREBIRD‐II Potentially Driven by EMIC Waves: Location, Extent, and Energy Range From a Multievent Analysis
Geophysical Research Letters (2021)
Catalogue of events with cumulative fluxes below as well as comparable and exceeding those for RBSP discovered third radiation belt
Astrophysics and Space Science (2021)