Energetic Electron Precipitation Driven by Electromagnetic Ion Cyclotron Waves from ELFIN’s Low Altitude Perspective

We review comprehensive observations of electromagnetic ion cyclotron (EMIC) wave-driven energetic electron precipitation using data collected by the energetic electron detector on the Electron Losses and Fields InvestigatioN (ELFIN) mission, two polar-orbiting low-altitude spinning CubeSats, measuring 50-5000 keV electrons with good pitch-angle and energy resolution. EMIC wave-driven precipitation exhibits a distinct signature in energy-spectrograms of the precipitating-to-trapped flux ratio: peaks at >0.5 MeV which are abrupt (bursty) (lasting ∼17 s, or \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$\Delta L\sim 0.56$\end{document}ΔL∼0.56) with significant substructure (occasionally down to sub-second timescale). We attribute the bursty nature of the precipitation to the spatial extent and structuredness of the wave field at the equator. Multiple ELFIN passes over the same MLT sector allow us to study the spatial and temporal evolution of the EMIC wave - electron interaction region. Case studies employing conjugate ground-based or equatorial observations of the EMIC waves reveal that the energy of moderate and strong precipitation at ELFIN approximately agrees with theoretical expectations for cyclotron resonant interactions in a cold plasma. Using multiple years of ELFIN data uniformly distributed in local time, we assemble a statistical database of ∼50 events of strong EMIC wave-driven precipitation. Most reside at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$L\sim 5-7$\end{document}L∼5−7 at dusk, while a smaller subset exists at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$L\sim 8-12$\end{document}L∼8−12 at post-midnight. The energies of the peak-precipitation ratio and of the half-peak precipitation ratio (our proxy for the minimum resonance energy) exhibit an \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$L$\end{document}L-shell dependence in good agreement with theoretical estimates based on prior statistical observations of EMIC wave power spectra. The precipitation ratio’s spectral shape for the most intense events has an exponential falloff away from the peak (i.e., on either side of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$\sim 1.45$\end{document}∼1.45 MeV). It too agrees well with quasi-linear diffusion theory based on prior statistics of wave spectra. It should be noted though that this diffusive treatment likely includes effects from nonlinear resonant interactions (especially at high energies) and nonresonant effects from sharp wave packet edges (at low energies). Sub-MeV electron precipitation observed concurrently with strong EMIC wave-driven >1 MeV precipitation has a spectral shape that is consistent with efficient pitch-angle scattering down to ∼ 200-300 keV by much less intense higher frequency EMIC waves at dusk (where such waves are most frequent). At ∼100 keV, whistler-mode chorus may be implicated in concurrent precipitation. These results confirm the critical role of EMIC waves in driving relativistic electron losses. Nonlinear effects may abound and require further investigation.


Earth's Radiation Belts
Relativistic electron fluxes in Earth's radiation belts pose a significant hazard to satellites (Horne et al. 2013) and astronauts, especially during magnetic storms (Gonzalez et al. 1994;Baker et al. 2018). These fluxes wax and wane in response to upstream solar wind variations (Baker et al. 1987;Reeves et al. 1998Reeves et al. , 2003, reflecting a delicate competition between acceleration, transport, and loss within the magnetosphere in ways that still defy accurate forecasting . Additionally, energetic electrons precipitate to the highlatitude mesosphere and lower thermosphere, where they can create NO x and H O x , ozonedestroying catalysts (Jackman et al. 1980;Thorne 1980;Randall et al. 2005). Even though NO x has a short lifetime when in sunlight, during polar winter darkness it can last for days to weeks and can be brought by vertical winds to 30 − 50 km altitude, near the stratospheric peak of the ozone layer. Under such conditions it can catalytically convert ozone and contribute significantly to ozone destruction. In particular, relativistic electrons, those at energies of hundreds of keV to several MeV, pose the most significant threat to space assets. These electrons often have both sufficient fluxes and energy to penetrate through spacecraft and space suits, causing deep dielectric charging or high levels of radiation exposure. Additionally, when > 1MeV electrons are scattered by magnetospheric processes and precipitate to the upper stratosphere, they can produce NO x directly in the regions of dominant ozone concentration where they can be most destructive (Baker et al. 1987).
The transient nature of the high energy electron flux and precipitation complicates forecast efforts. The trapped flux is greatly affected by local acceleration (as opposed to transport) via electron resonant interaction with intense electromagnetic whistler-mode waves (Chen et al. 2007; Thorne et al. 2013;Li et al. 2014;Allison and Shprits 2020); by heating during injections or during radial diffusion via intense ultra-low-frequency (ULF) waves (Elkington et al. 2004;Hudson et al. 2012;Mann et al. 2016;Sorathia et al. 2017); by magnetopause shadowing at the dayside leading to rapid electron losses (Shprits et al. 2006;Turner et al. 2012;Mann et al. 2016;Olifer et al. 2018;Sorathia et al. 2018;Pinto et al. 2020); by field-line scattering at the nightside (Sergeev and Tsyganenko 1982;Artemyev et al. 2013); and by wave-driven precipitation into the upper atmosphere Blum et al. 2015b;Shprits et al. 2017). During active times, energetic electrons in the radiation belts can occasionally be sourced adiabatically by inward diffusion from outer L-shells where the electron phase space density is high, due to enhanced electric fields (Dessler and Karplus 1961;Kim and Chan 1997) or ULF waves (Elkington et al. 2003;Mann et al. 2016). They can also be sourced from even lower energy, ∼ 10 keV, perpendicularly anisotropic magnetospheric electrons at larger distances. This can occur by rapid injections (transport and simultaneous acceleration) of the source electrons to the inner magnetosphere, leading to tens to hundreds of keV seed electrons Gabrielse et al. 2017). Waves, in particular whistler-mode chorus, are also excited by these anisotropic, highly unstable source electrons in the magnetosphere. These waves can cause further rapid acceleration of the seed electrons to even higher, relativistic energies . As the energetic electrons are transported into closed drift shells they further interact with waves in the ULF, ELF and VLF range to cause both local acceleration at low L-shells (Elkington et al. 1999;Horne et al. 2005;Ma et al. 2015;Li et al. 2016;Thorne 2010; Thorne et al. 2013) and scattering into the loss cone Millan and Thorne 2007). Drift-shell splitting and other dynamic effects such as solar wind compression pulses can cause further wave excitation and scattering of particles into the loss cone (Jin et al. 2022;Ukhorskiy et al. 2014).
Acceleration and loss of radiation belt electrons can occur simultaneously, over a wide range of temporal scales (from seconds to weeks) and spatial scales (across different local times, L-shells and latitudes). It is evident that precipitating electron fluxes also result from a dynamic competition between acceleration, transport, and loss processes. The two main waves responsible for such losses are whistler-mode chorus and EMIC waves. The efficacy of such wave-driven losses must then be studied in the geophysical context (geomagnetic activity, plasma environment and location) within which these waves occur. Scattering of > 1 MeV electrons by chorus waves in a plasma of realistic (on the order of 1 − 10cm −3 ) density requires cyclotron resonance at high magnetic latitudes along the magnetic field lines (|λ| > 30 • ; see Shprits et al. 2006;Summers et al. 2007;Agapitov et al. 2018)). However, at those latitudes the average whistler-mode wave intensity is observed to be weak (Meredith et al. 2012;Agapitov et al. 2013Agapitov et al. , 2018Wang et al. 2019). Thus, local (i.e., unrelated to magnetopause shadowing) and rapid losses of > 2 MeV electrons are typically attributed to their resonant scattering by EMIC waves (Usanova et al. 2014;Omura and Zhao 2012;Mourenas et al. 2016;Drozdov et al. 2017;Qin et al. 2019). This interaction is also believed to be a key process controlling the dynamics of the relativistic (as low as hundreds of keV) electron fluxes in Earth's radiation belts (Thorne and Kennel 1971). EMIC waves can be very effective in electron scattering both in the quasi-linear diffusion regime (Ni et al. 2015a;Drozdov et al. 2017) and in the nonlinear resonant interaction regime (Albert and Bortnik 2009;Kubota et al. 2015;Grach and Demekhov 2020;Bortnik et al. 2022;Grach et al. 2022). In fact, both event and statistical studies have shown that precipitating fluxes of relativistic electrons at low-altitudes correlate well with equatorial EMIC wave bursts (Blum et al. 2015a;Capannolo et al. 2019;. However, the transient nature of EMIC wave emissions (Blum et al. 2016 and the effects of hot ions on the EMIC wave dispersion relation (Cao et al. 2017a;Chen et al. 2019) complicate the evaluation of the relative contribution of EMIC waves to multi-MeV, MeV, and sub-MeV electron losses in the radiation belts. This is a question that remains open for both observational and theoretical reasons, which we further detail below.

Relative Impact of EMIC Waves on Relativistic Electron Precipitation
Trapped fluxes of electrons of relativistic energy, associated with EMIC and chorus wavedriven precipitation, are monotonically decreasing with energy. EMIC waves are expected to be mainly responsible for multi-MeV electron scattering, but exhibit a minimum energy of precipitation potentially extending down to several hundred keV under some specific wave and plasma conditions (Summers and Thorne 2003;Cao et al. 2017a;. All other wave modes, and in particular whistler-mode chorus waves, are more effective at pitchangle scattering tens of keV particles, with progressively reduced efficacy at hundreds of keV (except under very special circumstances, such as nonlinear scattering and ducted wave propagation away from the equator (Horne and Thorne 2003;Miyoshi et al. 2020;Zhang et al. 2022)). Thus, the precipitating-to-trapped electron flux ratio, plotted as a function of energy, should be a good indicator of the scattering mechanism. In particular it would be ideal for separating the precipitation induced by EMIC waves from that by other types of waves. However, due to the scarcity of low-altitude satellite data with good pitch-angle and energy resolution in the relevant (sub-MeV to multi-MeV) energy range, this has been difficult in the past, and the relative contribution of the various wave modes to sub-MeV electron losses remains an open question.
This has not been for the lack of trying, as the contribution of energy and pitch-angle scattering by various waves to the overall flux levels and spectra is critically important for modeling and predicting space weather. Theoretical modelling using diffusion theory of EMIC waves, whistler-mode chorus and hiss does a decent job in predicting the high equatorial-pitch-angle flux decay as a function of activity indices, such as Kp or Dst, over many days to weeks Ma et al. 2015;Drozdov et al. 2017;Mourenas et al. 2017;Pinto et al. 2019). Similarly, an overall agreement has been found between diffusion theory predictions and observations, from the combination of low-altitude POES satellites and equatorial Van Allen Probes for the flux decay due to chorus waves at low equatorial-pitch-angles Reidy et al. 2021;Mourenas et al. 2021). However, the relative contribution of each mode to relativistic electron acceleration, precipitation, and short-term flux-evolution has been more difficult to pin down. This is particularly true for the relative contribution of EMIC waves to relativistic (hundreds of keV to MeV range) electron scattering. This relative contribution has been surmised recently from Versatile Electron Radiation Belt modeling comparisons with Van Allen Probes data. It was shown that EMIC waves are critical for the flux evolution of > 1 MeV electrons, but in order to explain these electrons' equatorial-pitch-angle spectra over a broad range of L-shells, a combination of EMIC waves and whistler-mode chorus and hiss waves was needed (Drozdov et al. 2015(Drozdov et al. , 2020. Indeed it has been theoretically and observationally shown that the efficacy of MeV electron precipitation by EMIC waves is enhanced in the presence of whistler-mode waves even when the two wave types are operating at different local times, because EMIC waves alone cannot cause precipitation of the most abundant, high pitch-angle electrons Mourenas et al. 2016;Zhang et al. 2017;Qin et al. 2019).
Another example of an important unanswered question in the same area is the origin of microburst precipitation. This is a particularly intense, short-lived, electron precipitation phenomenon, lasting on the order of 1s or less. It is thought to contribute significantly to the overall energetic electron losses at 100s of keV to several MeV (Blake et al. 1996;O'Brien et al. 2004;Blum et al. 2015b;Greeley et al. 2019;Hendry et al. 2019). While the sub-MeV energy range of this precipitation could be consistent with whistler-mode resonant interactions with electrons, especially since the broad spatial extent of microbursts at dawn overlaps with the typical location of whistler-mode chorus (Douma et al. 2017;Shumko et al. 2018;Zhang et al. 2022), theoretical studies and observations suggest that intense EMIC waves can also drive relativistic electron microbursts, especially longer duration ones, at dusk (Blum et al. 2015a;Zhang et al. 2016a;Kubota and Omura 2017). As the two wave modes (chorus and EMIC waves) are both able to scatter electrons in the hundreds of keV to ∼ 1 MeV range, overlap in their spatial distribution, and both can exhibit a bursty nature, it remains unclear which wave mode dominates relativistic microburst precipitation.
Thus, the relative contribution of EMIC waves to scattering loss of sub-MeV to few MeV electrons compared to the contribution by other waves is still an open question. This question is important for accurately modeling short-term variations of both radiation belt fluxes and the atmospheric response to relativistic electron precipitation. Part of the difficulty in addressing this question can be attributed to the previous lack of energy-resolved and pitch-angle resolved spectra of precipitating electrons in the 10 s of keV to a few MeV range that would enable a quantitative validation of theoretical models of diffusion. This situation has changed with the recent launch of the ELFIN mission, which provides for the first time such data over a wide range of L-shells and local times. ELFIN's multi-year dataset allows us now to accurately compare precipitating-to-trapped electron flux ratios and, thus, to infer electron diffusion rates as a function of energy (Kennel and Petschek 1966;Li et al. 2013). This is especially useful for comparisons of theoretical expectations of such diffusion rates and measurements of those rates using ELFIN observations, especially when combined with conjugate, equatorial measurements of waves and plasma parameters at THEMIS, Van Allen Probes, Arase and MMS. This represents a significant improvement over the otherwise massive and previously well-utilized POES dataset (Rodger et al. 2010;Yando et al. 2011;Yahnin et al. 2016;Capannolo et al. 2019;Qin et al. 2018Qin et al. , 2020, which can only provide limited pitch-angle and integral (high-) energy spectra. It is also an improvement over the dataset from the 1972-076B mission (Imhof et al. 1977), which had similar pitch-angle and energy resolution as ELFIN but was not accompanied by conjugate equatorial missions.
ELFIN was proposed as a focused investigation to address the question of whether EMIC waves, the primary candidate for pitch-angle scattering of relativistic electrons, can be definitively proven to be responsible for such scattering using the advantages offered by its new dataset. In this paper we aim to achieve that objective and exemplify the salient features that accompany such scattering. We will address this objective using ELFIN together with its numerous fortuitous conjunctions with equatorial spacecraft and ground observations. We first review, below, the properties of EMIC waves and their interaction with relativistic electrons. We next discuss how chorus waves may be also implicated in the scattering and precipitation of such electrons and how to differentiate the effects of these two wave types. Next, we present the first comprehensive ELFIN measurements of EMIC wave-driven electron precipitation. We discuss the observed features of the precipitating electron fluxes that indicate nonlinear resonant interaction of EMIC waves with electrons, compare precipitating electron energy spectra at high energy-resolution with theoretical expectations, and provide the first statistical distributions of EMIC wave-driven precipitation and its properties.

Generation
EMIC waves were first postulated to be excited by a low density, high energy population of hot ions which achieve cyclotron resonance with the ion cyclotron wave of a cold, dense ion background by appropriately Doppler-shifting the wave's frequency in their own frame through streaming along the magnetic field (Cornwall 1965;Cornwall et al. 1970). Such conditions prevail near the plasmapause where drift-shell splitting of ring current ions, or fresh ion injections, or magnetopause compressions of ambient, low density hot plasma may acquire perpendicular anisotropy. Portions of this (hot) ion distribution having a field-aligned streaming velocity that can thus attain cyclotron resonance with the waves can liberate the free energy available in their anisotropy to achieve wave growth (Kennel and Petschek 1966;Cornwall et al. 1970. This resonance condition is: ω − k v i, = n ci . Here, n = +1 is the relevant harmonic resonance number corresponding to first order resonance; ci and v i are the ion cyclotron angular frequency and ion velocity, respectively; ω and k are the wave angular frequency and wave number, respectively; and the parallel symbol denotes components along the ambient magnetic field. Electromagnetic waves of the background (cold, presumed dominant) plasma population propagating opposite to the beam (k v i, < 0) over a range of frequencies near, say, 0.5 ci and with wave vectors satisfying the cold plasma dispersion relation can thus become unstable. The dispersion relation of the cold component for parallel propagation (assuming that ions are protons and that the hot anisotropic component has a sufficiently low density to make a negligible contribution to the plasma dielectric response) is: Alfvén speed, pi the plasma frequency and other symbols having their usual meaning (note that: c/V A ≡ pi / ci ). This is a monotonic function of ω, approaching the Alfvén wave dispersion relation in the low frequency, MHD limit (ω → 0). In the high frequency limit, ω → ci as k → ∞, which means that in this limit the waves are absorbed by the plasma and cannot propagate -this is the ion cyclotron resonance. At intermediate frequencies, though, when the cyclotron growth provided by the hot component exceeds cyclotron damping by the cold component, the waves can grow. At oblique propagation, the dispersion relation near ci is only slightly modified, becoming: Using this dispersion relation, the aforementioned cyclotron resonance condition can be recast as Maximum growth occurs for parallel propagation since at oblique propagation the resonant velocity decreases and the ion cyclotron damping by the cold component prevails quickly, due to that component's high density and low temperature. It is evident from the above resonance condition that EMIC wave generation depends critically on the ratio pi / ci = c/V A (or equivalently on f pe /f ce , the ratio of electron cyclotron and plasma frequencies that we use more commonly below) and on 1 − x, the latter denoting the proximity of the wave frequency to the ion cyclotron frequency. These parameters determine the resonance energy and its proximity to the free energy available in the velocity-distribution's anisotropy. Typical EMIC wave excitation requires that this resonance energy be low enough for the waves to resonate with anisotropic hot magnetospheric ions in the few keV to 10 s of keV range. Hence, the larger the aforementioned frequency ratios are (the smaller V A is), the easier it is for EMIC waves to resonate with the free energy source of hot ions typically available. Because the Alfvén speed increases rapidly away from the equator along a field line, conditions at the geomagnetic equator favor such wave excitation. At high-density equatorial regions that are far enough from Earth so the geomagnetic field is also low, such as near the plasmapause, or within plasmaspheric plumes, pi / ci can increase (V A can decrease) sufficiently for EMIC waves to be excited if anisotropic hot ions are also present.
And indeed, EMIC waves are often excited near the post-noon, dusk and pre-midnight sectors where the cold, dense plasmaspheric bulge and plume (Horne and Thorne 1993) are intersected by the drift-paths of (hot) ring current ions exhibiting velocity space anisotropies. The cold, dense background plasma there is critical for lowering the resonance energy into the range where there is a sufficient number of hot ions with large enough anisotropy. This situation occurs in that sector, especially during storm times, according to case studies (Kozyra et al. 1997;Jordanova et al. 1998). However, statistical studies have also revealed that banded, low-frequency electromagnetic waves exist at other local times as well (Anderson et al. 1992;Erlandson and Ukhorskiy 2001;Fraser et al. 2010;Min et al. 2012;Meredith et al. 2014;Allen et al. 2016;Paulson et al. 2017). These waves have amplitudes 0.1 − 10nT, are typically left-hand polarized and field-aligned near the equator, and can extend from the Alfvén mode at low frequencies upwards to the local ion cyclotron frequency . They too can be identified as EMIC waves. Further supporting this identification is that in the presence of a multi-component plasma, typically H + with a few percent of either H e + , O + , or both, such waves are observed to split into the classical EMIC wave distinct bands (a H + , H e + and O + band), each between their respective gyrofrequency and the one below it, except that the O + band extends continuously below the O + gyrofrequency down to the Alfvén branch (e.g., (Cornwall and Schulz 1971;Young et al. 1981;Horne and Thorne 1993)).
While the highest amplitude waves are most frequently observed at the duskside equator in the H e + band with left-hand circular polarization and nearly field-aligned propagation, lower amplitude waves are also routinely observed at the dawnside equator except in the H + band with linear or elliptical polarization and occasional oblique propagation (Min et al. 2012). They are also seen further away from the equator, where they become oblique, likely due to their propagation, and are eventually (at high enough latitudes) Landau damped. Hence off-equatorial waves are seen with lower amplitudes and occurrence rates. However, such waves can occasionally also be ducted. Then they can propagate nearly-field-aligned and evade damping, thus reaching the ionosphere and the ground (Kim et al. 2010) where they are detected (Engebretson et al. 2008(Engebretson et al. , 2015 as continuous magnetic pulsations of Type 1 (Pc1, 0.2-5 Hz) or Type 2 (Pc2, 5-10 Hz). Other means of evading damping are mode conversion to the R-mode, and tunneling near the bi-ion frequency, just below the respective ion gyrofrequency . Substorm-related, freshly injected, anisotropic ring current ions drifting duskward from midnight and interacting with the plasmaspheric bulge or plume are most often responsible for exciting the EMIC waves seen at the duskside (Cornwall and Schulz 1971;Chen et al. 2010;Morley et al. 2010). However, solar wind compressions can also cause hot ions drifting in the inner magnetosphere with preexisting moderate (marginally stable) anisotropy to attain (through betatron acceleration) an enhanced anisotropy, one that exceeds the threshold for EMIC wave growth. This excitation mechanism is often credited for EMIC wave observations at the dayside, at pre-and post-noon (Anderson and Hamilton 1993;Arnoldy et al. 2005;Usanova et al. 2008Usanova et al. , 2010McCollough et al. 2010). However, prolonged, quiet time EMIC wave activity over a broad range of local times in the dayside (pre-and post-noon), but over a narrow L-shell range, is attributed to the large anisotropy of freshly supplied ions from the nightside plasmasheet by injections (that can persist at low occurrence rates at large distances even during geomagnetically quiet conditions). Such anisotropy develops at the dayside due to differential drifts at different energies and can excite EMIC waves near an expanded plasmasphere (Anderson et al. 1996;Engebretson et al. 2002).

Interaction with Electrons
Relativistic electron pitch-angle scattering due to their resonant interaction with H + EMIC waves was first considered by Thorne and Kennel (1971) and Lyons and Thorne (1972). Horne and Thorne (1998) calculated the minimum resonance energies for a multi-ion plasma (H + , H e + , O + ) inside and outside the plasmapause during storms. Summers et al. (1998) addressed relativistic effects, showing that even with such corrections, electrons in gyroresonance with EMIC waves undergo nearly pure pitch-angle (but not much energy) diffusion. Summers and Thorne (2003) demonstrated that such interactions can rarely result in scattering of electrons at or below 1 MeV. Such conditions arise only for pe / ce ≥ 10 (or equivalently ci / pi ≡ V A /c ≤ 4.3) which occurs near and just inside the dusk plasmapause, most often at storm times (where ce is the unsigned electron cyclotron angular frequency). To understand why, we discuss below the fundamental characteristics of this interaction.
Electrons can resonate with an EMIC wave by overtaking (moving in the same direction, but faster than) the wave if they have a sufficiently high (relativistic) speed to Dopplershift the very low sub-ion-cyclotron wave frequency to the very high, electron cyclotron frequency. The left-hand circularly polarized (in time) EMIC wave electric field vector tip carves a right-handed helical wave-train in space. From the viewpoint (in the frame) of the guiding center of an electron able to overtake the wave's helix crests and valleys, the electric field vector tip rotates now (in time) in a right-handed way, opposite to that in the (ion or plasma) rest frame. This polarization reversal has the potential to put the electron, also gyrating in a right-handed sense, in cyclotron resonance with the EMIC wave. The generic electron cyclotron resonance condition is: ω − k v e, = n ce /γ . Here v e is the velocity of the electron, which in our case has a projection along the magnetic field that is in the same direction as the wave's projection along the field (k v e, > 0), γ is the relativistic correction factor (the Lorentz factor) and n = −1, for first order anomalous cyclotron resonance. (Anomalous, because due to the aforementioned overtaking, the sense of polarization experienced by the electron is opposite to that of the wave in the plasma rest frame.) Since ω ce , the resonance condition becomes simply: k v e, ∼ ce /γ . The electron resonance energy obtained from the aforementioned cold plasma dispersion relation of ion cyclotron waves and from the above electron cyclotron resonance condition (Thorne and Kennel 1971), simplified for a proton-electron plasma is E R /E 0 = γ − 1, where E 0 is the electron rest mass, and (based on the solution of the above two equations) γ is given by: The minimum resonance energy, E R,min , is obtained for zero pitch angle, α = 0, for a given total energy. The most common situation of parallel propagation (ξ = 1) serves as a case-inpoint. Moreover, for a fixed wavenumber, θ = 0 also minimizes the resonance energy. We see that E R,min , corresponding to γ 2 − 1 = ce / pi √ 1 − x/x, is a monotonic function of x = ω/ ci and ce / pi , so the closer ω gets to ci the lower the E R,min . This is seen more clearly if the resonant velocity in the resonance condition above can be simply recast as resonance energy: E R (k) = E 0 1 + ( ce /ck) 2 − 1 . This is minimum for the maximum unstable wavenumber, which (based on the cold plasma dispersion relation, seen earlier) corresponds to the maximum ω/ ci , closest to 1.
For fixed ω/ ci ∼ 0.5, E R (k) falls off with L-shell as a power law in the plasmasphere, due to the magnetic field decreasing faster than the square root of the density (Sheeley et al. 2001;Ozhogin et al. 2012). At the plasmapause, E R,min increases abruptly with L-shell (outward) by about an order of magnitude, to ∼ 10 MeV as the density drops by 1-2 orders of magnitude (Cornwall 1965;Thorne and Kennel 1971). Therefore, E R,min has a local minimum (near ∼ 1 MeV) just at the interior of the plasmapause. This situation remains true for EMIC resonances with heavy ions, when those are included in the dispersion relation (Summers and Thorne 2003).
However, incorporating thermal effects in the cold plasma dispersion relation complicates this picture (Chen et al. 2011). When even a fraction of the low-energy ions has a significant temperature (10s to 100s of eV), as is often observed (Lee and Angelopoulos 2014), the dispersion relation is significantly modified: the waves can propagate through their respective cyclotron frequencies and the stop bands can vanish (Silin et al. 2011;Chen et al. 2011;Lee et al. 2012). While the dispersion relation becomes more complex, and the wave frequency is not limited to just below, or between the ion gyrofrequencies as the case may be, heavy cyclotron damping by the cold species near those frequencies severely limits wave propagation away from the source, even when the waves are, in principle, unstable due to an exceedingly strong anisotropy of the hot ions. These conditions cause excessive wave damping at large wavenumbers, those with kc/ pi higher than ≈ 1. Yet, the resonance condition, expressed above as the resonance energy as a function of wavenumber, E R (k), still applies regardless of the dispersion relation, and shows that there is still a lower limit to the minimum resonance energy, the one for the maximum wavenumber permitted for propagation, even with warm plasma effects accounted for. This realization simplifies the analysis: approximating the maximum wavenumber that can be attained under the presence of thermal effects as kc/ pi ≈ 1, we obtain a similar, monotonic dependence of the resonance energy on ce / pi as for the cold plasma approximation: E R,min ∼ E 0 1 + ( ce / pi ) 2 − 1 . We will compare this relationship with data, later in the paper.
A parametric analysis of the instability for multi-species plasmas including warm plasma effects confirms that the maximum unstable wavenumber rarely results in E R,min below 1 MeV: this only occurs for conditions of large pe / ce (15 to 100), and a large hot species anisotropy A 2 (Chen et al. 2011. In those cases, electrons of energy as low as ∼ 500 keV may be able to resonate with and be scattered by waves of sufficiently high frequency. (Note that diffusion rates still peak at energies near E peak higher than E R,min , corresponding to frequencies near peak-wave-power frequency that are lower than the maximum observed frequency of wave propagation linked to E R,min .) In particular, H e + EMIC waves which are most easily able to resonate with 1MeV electrons in cold plasma theory are strongly suppressed by cyclotron absorption; warm plasma effects cause the H + band to resonate more readily with 0.5 − 1 MeV electrons than the H e + band ). Thus, even though warm plasma effects allow EMIC wave spectra to reach closer to and even cross the cyclotron frequency, consistent with some observations, at least in the context of quasi-linear theory E R,min still remains most often above 1 MeV except in rare cases of high density regions such as plumes at high L-shells  or at low L-shells for compressed plasmaspheric conditions during the storm main phase (Cao et al. 2017a;Qin et al. 2019). Observations of precipitating electrons from POES, albeit with instruments of limited energy and pitch-angle resolution , have shown that > 0.7 MeV electron precipitation can indeed be observed at POES, preferentially when equatorial spacecraft in close conjunction with POES confirm the existence of plasma conditions favorable for 0.7 − 1MeV electron scattering by H + waves. Note, though, that POES does not have differential energy channels to finely resolve the peak in precipitation as a function of energy, hence these results should be considered suggestive, not conclusive evidence for the operation of EMIC waves. Additionally, many counter-examples were also found (when theoretically expected precipitation from waves was not observed, or was unexpected but observed) suggesting that quasi-linear theory alone may not be able to fully explain these observations. Further supporting the latter suggestion is that on occasion, relativistic electron precipitation (∼ 1 MeV) events can occur on timescales of a few seconds or less (Imhof et al. 1992;Lorentzen et al. 2000Lorentzen et al. , 2001O'Brien et al. 2004;Douma et al. 2017), whereas the usual timescales of quasi-linear diffusion are on the order of many minutes to hours (Albert 2003;Li et al. 2007;Ukhorskiy et al. 2010;Ni et al. 2015a). Such short-lived precipitation can often extend down to hundreds of keV. These counter-examples cast doubt on the ability of EMIC waves to fully explain the observations, at least when studied in the quasi-linear regime even when hot plasma effects are incorporated into the theory.
Nonlinear treatments of EMIC wave interaction with relativistic electrons have also resulted in some successes in interpreting observations of rapid sub-MeV electron precipitation. Early work initially showed that nonlinear interaction with moderate amplitude, fixed frequency waves in a dipole field typically leads to scattering towards large pitch angles, away from the loss cone (Albert and Bortnik 2009). However, the observed departures of the EMIC waveforms from a constant frequency and the presence of a magnetic field gradient near the equator can cause phase trapping of resonant electrons and result in very rapid pitch-angle scattering and precipitation of MeV and even sometimes sub-MeV energies (Omura and Zhao 2012;Kubota et al. 2015;Hendry et al. 2017;Nakamura et al. 2019;Grach et al. 2021). This effect can be enhanced by diffusive scattering by large amplitude EMIC waves, which may transport electrons directly into the loss cone from intermediate (∼ 30 • ) pitch-angles . For realistic EMIC waveforms having sufficiently steep edge effects, or equivalently having a few wave periods in a single packet, even sub-MeV nonresonant electrons can be pitch-angle scattered, when their interaction occurs over a small number of gyroperiods (Chen et al. 2016;An et al. 2022). Additionally, bounce resonance of near-equatorially mirroring, hundreds of keV energy electrons with EMIC waves can also result in moderate pitch-angle scattering and contribute to precipitation at those energies (Cao et al. 2017b;Blum et al. 2019). But since hundreds of keV electrons can also interact with chorus waves, which may occur simultaneously with EMIC waves, either at different local times  or even at the same location when driven by ULF pulsations Bashir et al. 2022), an unambiguous determination of the distinct (let alone independent) EMIC wave contribution to the precipitation can be difficult.

Identification in Precipitation Spectra
Previous studies have presented suggestive evidence of telltale signatures of EMIC wavedriven relativistic electron precipitation. This was achieved either using in-situ magnetospheric observations of depletion of near-field-aligned flux (in velocity-space) concurrent with EMIC wave enhancements (Usanova et al. 2014;Zhang et al. 2016c;Bingley et al. 2019;Adair et al. 2022), or by identifying local minima in radial profiles of the phasespace density at L-shells consistent with simultaneous ground observations of EMIC waves , or through observations of simultaneous precipitation of 10s of keV protons and MeV electrons (Imhof et al. 1986;Miyoshi et al. 2008;Hendry et al. 2017;Capannolo et al. 2019). However, chorus waves can also scatter and cause precipitation of electrons of hundreds of keV to MeV Ma et al. 2016b;Miyoshi et al. 2020;Zhang et al. 2022). The relative contribution of chorus and EMIC waves was not addressed in those studies (e.g., see discussions in ).
Noting that typical chorus wave scattering is most effective at tens of keV rather than at hundreds of keV, a monotonically decreasing precipitating-to-trapped flux ratio as a function of energy would favor a chorus wave scattering interpretation over an EMIC wave scattering one. Conversely, that ratio increasing with energy, particularly when peaking at 1 MeV or greater, would favor the EMIC wave scattering interpretation, since EMIC waves are most effective scatterers at > 1 MeV electron energies. However, electron spectra of sufficiently high resolution in energy and pitch-angle to make the above distinction were not available in prior studies, which were mostly based on POES data (Evans and Greer 2004;Yahnin et al. 2017;Capannolo et al. 2018Capannolo et al. , 2019. Such high resolution spectra, obtained at a low-altitude (ionospheric) satellite, especially when combined with equatorial or ground-based measurements of the EMIC waves, are critical for determining if such waves are responsible for relativistic electron scattering, and for addressing the physical mechanism of the scattering process (quasi-linear, nonlinear, resonant or nonresonant, etc). Such measurements are needed not only to identify but also to quantify EMIC wave-driven precipitation and its role in radiation belt dynamics and magnetosphere-atmosphere coupling.
EMIC wave resonant interactions with electrons can be attributed to (and studied as) one of two processes: quasi-linear diffusion toward the loss cone (Kennel and Petschek 1966;Lyons 1974) and fast nonlinear phase trapping transport toward the loss cone (Albert and Bortnik 2009;Kubota et al. 2015;Grach et al. 2021;Bortnik et al. 2022;Grach et al. 2022). The relative importance and occurrence rate of these two regimes of wave-particle interaction for EMIC wave scattering has not been addressed yet, even though there is consensus from observations that EMIC waves Saikin et al. 2015;Zhang et al. 2016c) are often sufficiently intense to resonate with electrons nonlinearly . The strongest losses associated with quasi-linear diffusion, those in the strong diffusion limit, have (by definition) loss-cone fluxes comparable to trapped fluxes, those next to the loss-cone edge (Kennel and Petschek 1966). However, nonlinear electron interaction may exceed the strong diffusion limit and produce loss-cone fluxes higher than trapped fluxes (Grach and Demekhov 2020). Distinguishing these two precipitation regimes requires electron flux measurements at fine pitch-angle resolution near and within the loss cone, which is possible with energetic particle detectors of modest angular resolution observing from low altitudes.
The recently launched, ELFIN CubeSat twins, ELFIN A and ELFIN B, provide a new dataset of precipitating electrons that is very helpful for addressing the above questions related to the process and efficiency of EMIC wave resonant scattering of energetic electrons. Their energetic particle detector for electrons (EPDE) measures the full 180 • pitch-angle distribution of electron fluxes with approximately 22.5 • resolution, over the energy range 50 − 5000 keV sampled at 16 logarithmically-spaced energy channels of width E/E < 40%. Thus, they can resolve perpendicular (locally trapped), precipitating, and backscattered fluxes with good pitch-angle and energy resolution . Due to ELFIN's altitude, 300-450 km, the locally-trapped (perpendicular) flux measured corresponds to particles that are most often in the drift loss cone, i.e., destined to be lost before they complete a full circle around the Earth due to the variation of the geomagnetic field magnitude with geographic longitude. Near the longitude of the south-Atlantic anomaly, in the northern hemisphere the perpendicular fluxes are still inside the bounce loss cone (they will precipitate in the south) but even in that case, intense fluxes generated locally in the same hemisphere above, at the equator, will still provide valuable information on EMIC wave scattering at a rate faster than a quarter-bounce period, and are therefore valuable to retain. However, at most longitudes the measured perpendicular fluxes still correspond to electrons outside the local maximum bounce loss cone, meaning that such electrons have had a chance to drift in longitude for some time. In this paper we simply refer to perpendicular fluxes as trapped, meaning at least quarter-bounce trapped, or locally trapped. In addition to the energetic particle detector for electrons, EPDE, instrument, which forms the basis of the investigations discussed in this paper, the ELFIN satellites carry an energetic particle detector for ions, EPDI, instrument measuring total ions with the same energy and pitch angle range and resolution as EPDE. It was calibrated in orbit late during the mission operations and delivered good data for the last three months of operations on both -A and -B satellites, before their deorbit (due to atmospheric drag) in mid-September 2022. Both satellites also carry a research-grade fluxgate magnetometer, FGM, instrument that operated flawlessly throughout the mission, and has returned low-noise (< 0.3 nT/ √ Hz at 1 Hz) data for several years on both spacecraft, up until deorbit. Its data are presently being calibrated. The EPDI and FGM data complement the analysis from the primary instrument, EPDE, as intended by the ELFIN mission design . They constitute high value datasets to be used in the future, in conjunction with the EPDE data to further study the electron scattering by EMIC waves, its relative contribution to precipitation compared to other mechanisms, and its origin at various magnetospheric plasma regions.
In Sect. 3, below, we present ELFIN examples of EMIC wave-driven electron precipitation. We show the salient features of that precipitation and its difference from whistler-mode precipitation, consistent with the prior discussion in the subsection above. In Sect. 4, we also incorporate in our analysis ancillary observations from other assets, such as conjugate measurements from ground-based stations and equatorial spacecraft. These are providing a regional context for the observed ELFIN precipitation (equatorial density and magnetic field), independent confirmation of the trapped particle fluxes and information on the occurrence and properties of EMIC waves that may be responsible for the observed precipitation. In Sect. 5, we take a statistical approach to the study of ELFIN's observed precipitation spectra attributed to EMIC waves. We show that their spectral properties, such as the peak precipitation energy and the slope of precipitating-to-(locally-)trapped flux ratio as a function of energy, as well as the spatial distribution of the inferred EMIC wave power are all consistent with expectation from theory and equatorial observations of these waves. We also find evidence of nonlinear interactions that can be further explored with the new dataset at hand.

ELFIN Examples of EMIC Wave-Driven Electron Precipitation
Moving along their low-Earth (∼ 400 km altitude), ∼ 90 min period orbits, the ELFIN Cube-Sats can, in principle, record up to four science zones per orbit (each zone including traversals of the near-Earth plasma sheet, outer radiation belt, plasmasphere, and inner belt). However, power and telemetry constraints demand judicious selection of (typically) 4-12 such zones per day per spacecraft. Choice of science zones (planned and scheduled with a weekly cadence) depends on conjunction availability with other missions and ground stations, or uniformity of coverage in time, MLT and hemisphere. For several months after ELFIN's launch in September 2018, there were conjunctions with the near-equatorial, dual satellite, ∼ 9 hour period Van Allen Probes mission (Mauk et al. 2013). During the four years (2018-2022) of ELFIN operations, there have been multiple conjunction periods with the equatorial Time History of Events and Macroscale Interactions during Substorms (THEMIS) mission (three spacecraft, roughly on a string-of-pearls orbital configuration when in the inner magnetosphere) on a ∼ 1 day orbital period, and ∼ 13R E apogee (Angelopoulos 2008)), as well as with the near-equatorial Exploration of energization and Radiation in Geospace (ERG) spacecraft (also known as ARASE; (Miyoshi et al. 2018)), and with the near-equatorial, ∼ 3 day period, four closely-separated satellite Magnetospheric Multiscale (MMS) mission (Burch et al. 2016)). Additionally there have been very useful ELFIN conjunctions with ground-based magnetometer stations providing magnetic field measurements in the EMIC wave frequency range. Such stations often detect equatorial EMIC waves propagating down to the ionosphere and associated with relativistic electron precipitation (Usanova et al. 2014;Yahnin et al. 2017).

EMIC Wave-Driven Versus Whistler-Mode Wave-Driven Electron Precipitation
Two ELFIN A science zone crossings in Fig. 1, one at the nightside/dawn flank (left) and the other at dayside/dusk flank (right), show precipitation signatures representative of the outer radiation belt and auroral zone, under moderately active conditions. (See comprehensive plots of these crossings also in the standard ELFIN overview plots at https://elfin.igpp.ucla. edu under Science → Summary Plots. Navigate to the appropriate time using the drop-down menus and then use mouse-clicks to zoom in and out in time on daily-plots and on science zone overview plots respectively.) To orient the reader, we first describe the magnetospheric regions ELFIN A traversed and discuss the predominant scattering mechanisms responsible for sub-MeV precipitation. We then explain how EMIC wave-driven MeV electron precipitation can be recognized within this otherwise fairly typical context. Figure 1 shows energy-time spectrograms of (locally) trapped fluxes (averaged over nearperpendicular pitch angles, those outside the local bounce loss cone, also referred to as "trap" in plots), precipitating fluxes (averaged over pitch-angles inside the local loss cone, and also referred to as "prec" in plots), as well as precipitating-to-trapped (or "prec-totrap" in plots) flux ratios. ELFIN A travelled on a post-midnight/dawn meridian from high Overview of two consecutive ELFIN A science zone crossings: one at the nightside/dawn sector (left) and the other at the dayside/dusk sector (right, primed panel letters). From top to bottom shown are 3 energyspectrograms (a-c), 2 pitch-angle spectrograms (d-e) and the satellite's L-shell and magnetic latitude (f) computed using the international geophysical reference field (IGRF) model. All spectrograms show products derived from the number-flux of electrons (measured in individual sectors, in units of 1/cm 2 /s/sr/MeV) averaged over the selected pitch-angle and energy range as follows: The energy spectrograms in (a) and (a ) are for locally trapped electrons (only pitch angles outside the loss cone and anti-loss cone, near perpendicular to the local field line direction, were included); those in (b) and (b ) are for precipitating electrons (with pitches in the loss cone); and those in (c) and (c ) are precipitating-to-trapped flux spectra ratios formed from the panels right above. The pitch-angle spectrograms in (d-e) and (d -e ) are average fluxes in two broad energy ranges: a low energy range, 80-270 keV, and a high energy range, 0.4 -3.4 MeV. The horizontal lines demarcate 90deg (vertically centered solid line), the loss cone (the other solid line) and the anti-loss cone (the dashed line). Horizontal color bars above Panel (a) represent magnetospheric regions identified based on the data and discussed in the main text. Arrows in Panels (c ) and (e ) represent spectral features also discussed in the main text to low L-shells as depicted in the bottom panel of that figure and as demarcated in the annotations at the bottom. (Note that L-shells and magnetic latitudes/longitudes throughout this paper are computed using the international geophysical reference model, IGRF.) The satellite was in the plasma sheet prior to ∼11:05:40 UT, and moved to the outer radiation belt and plasmasphere soon thereafter.
The plasma sheet is operationally identified by above-background (locally) trapped fluxes which do not exceed energies ∼ 200 keV, concurrent with a precipitating-to-trapped flux ratio of around one. The latter ratio is expected for energetic electrons that have been isotropized by the small field-line curvature radius and the low equatorial field intensity of the plasma sheet (Birmingham 1984;Büchner and Zelenyi 1989;Delcourt et al. 1994;Lukin et al. 2021). At ∼11:05:55 UT (shortly after, but still near the transition region between the plasma sheet and inner magnetosphere), ELFIN A observed a distinct, dispersive feature of the precipitating-to-trapped flux ratio: the lowest energy approaching the ratio of ∼ 1 increased with proximity to Earth (as L-shell decreased), a trend marked by an inclined dotted arrow in the figure. This energy versus L-shell dispersion is typical of the isotropy (magnetic latitude) boundary (see statistical analysis of such ELFIN observations in Wilkins et al. 2023). For a given energy, this boundary is the ionospheric magnetic latitude poleward of which electrons of that energy become isotropic (due to field-line scattering) and equatorward of which the same energy electrons are anisotropic (due to the strong intensity and large radius of curvature field in the inner magnetosphere that prevents such scattering (Sergeev and Tsyganenko 1982;Dubyagin et al. 2002;Sergeev et al. 2012)). At progressively higher energy this isotropy boundary typically resides at progressively lower latitude (corresponding to a lower L-shell), because electrons of higher gyroradius can still be fieldline scattered at the stronger and less curved equatorial magnetic field at this lower latitude (see Wilkins et al. 2023). This isotropy boundary signature in the precipitating-to-trapped ratio can be used as an additional identifier of the transition from the plasma sheet to the outer radiation belt.
Subsequently, between 11:05:55 UT and 11:08:00 UT ELFIN A was in the outer radiation belt: both the trapped flux magnitude at all energies increased, and the maximum energy of the (statistically significant) trapped fluxes also increased as the L-shell decreased. There, a low-level precipitation (ratio < 0.1) is intermittently punctuated by bursty, intense precipitation (ratio reaching ∼ 1). In this (post-midnight/dawn) science zone crossing, the precipitating-to-trapped flux ratio was highest at the lowest energies. This is expected for electron resonant interactions with whistler-mode (chorus) waves (Kennel and Petschek 1966). Such spectra are typical at ELFIN's outer belt crossings at post-midnight or dawn, especially at times of moderate to high geomagnetic activity.
At ∼11:07:40 UT, ELFIN A started to enter the plasmasphere, as evidenced by the decrease in the trapped electron fluxes and the simultaneous weakening of the precipitating fluxes. Shortly after that time, at ∼11:08:10 UT, the trapped fluxes at ∼ 200 − 300 keV decreased below instrument background level. This is a distinct feature of the plasmasphere, where 100 s of keV electrons are efficiently scattered by plasmaspheric hiss (also whistler-mode) waves Ripoll et al. 2016). This scattering, occurring at L ∼ 3.5 − 4.3, forms a slot region between the outer and inner belts that is nearly devoid of energetic electrons (Lyons andThorne 1972, 1973;Ma et al. 2016aMa et al. , 2017Mourenas et al. 2017). We thus interpret the ELFIN observations as entry into the slot region at ∼11:08:10 UT. The precipitating-to-trapped flux ratio in the plasmasphere (which was seen inside of L ∼ 4.5 in this event) is almost nil because (i) the flux of > 200 keV electrons near the losscone has decreased considerably, and (ii) hiss and lightning-generated whistler-mode wave power at frequencies that can attain cyclotron-resonance with ∼ 50 − 200 keV electrons is weak in the 17-24 and 00-05 MLT sectors inside the plasmasphere (Agapitov et al. , 2018Li et al. 2015), where the much higher plasma density compared to outside the plasmasphere further reduces the scattering efficiency of such whistler-mode waves (Mourenas and Ripoll 2012). Figure 1, right column, shows ELFIN A post-noon/dusk observations acquired about 1.5 hours later. In this science zone crossing, ELFIN A moved from low to high L-shells (see Panel (f )) and traversed the regions discussed in the previous crossing but in reverse order. The plasmasphere was traversed first, and the plasmapause was encountered at ∼12:24:10 UT (L ∼ 4.5) as evidenced by the transition from low to high trapped fluxes and from low to high precipitating-to-trapped flux ratio of ∼ 200 − 300 keV electrons. Between 12:23:40 UT and 12:26:40 UT, ELFIN A traversed the outer radiation belt, as evidenced by the significant fluxes of > 200 keV electrons. In that period it observed whistler-mode (chorus) wavedriven electron scattering, as implied by the bursty nature of the precipitation and by the precipitating-to-trapped flux ratio being largest at the lowest energies. After 12:26:45 UT ELFIN A was likely magnetically conjugate to the equatorial post-noon/dawn plasma sheet (and outside the outer radiation belt) as suggested by the reduction in the flux of trapped electrons of energy > 200 keV. But because in that local time sector the equatorial magnetic field is both stronger than and not as curved as that at the nightside plasma sheet at a similar L-shell, we do not expect plasma sheet field-line scattering to be significant, explaining the dearth of < 200 keV precipitation at ELFIN A at that time.
While in the outer radiation belt but near the plasmapause, ELFIN A also observed a distinctly different type of precipitation: bursts of precipitation at relativistic (E ≥ 1 MeV) energies. These bursts are demarcated by black arrows in the fifth panel (the one depicting the precipitating-to-trapped flux ratio). The most intense bursts occurred at ∼12:25:05UT and ∼12:25:28UT, but less intense ones are also evident around the same time period (all occurred within a ∼1 min interval, or L≈1.5). The MeV range precipitating-to-trapped flux ratio peaks are clearly separated in energy from the lower-energy peaks of chorus-driven precipitation discussed previously. The MeV range bursts are accompanied by increases in both trapped and precipitating fluxes from near-or below-background levels outside the bursts to above such levels within the bursts. We attribute this type of precipitation to EMIC wave scattering. Aside from the fact that EMIC waves resonate with such high energy electrons and are abundant in this region of space (at or near the post-noon and dusk plasmapause), the spectral properties observed by ELFIN A are also consistent with this interpretation: It is known that equatorial fluxes of relativistic electrons are typically very anisotropic and thus have a low flux level near the loss cone (Gannon et al. 2007;Chen et al. 2014;Shi et al. 2016;Zhao et al. 2018). Scattering by EMIC waves transports electrons to smaller pitch-angles in velocity-space over a broad range of pitch-angles, increasing both trapped and precipitating fluxes near (and on either side of) the loss cone. Since low-altitude spacecraft, like ELFIN A, measure locally mirroring fluxes that still map to low equatorial pitch angles, they experience EMIC wave scattering as flux increases at both locally trapped (i.e., locally mirroring or perpendicular) and precipitating electrons. In Sect. 4, we provide more examples of such relativistic electron precipitation, using conjunctions with ground-based and equatorial observatories measuring the relevant EMIC waves directly, further supporting our EMIC wave-driver interpretation of these precipitation signatures seen at ELFIN.

The Latitudinally Localized and Intense (Bursty) Nature of EMIC Wave-Driven Electron Precipitation
The event presented next is a prototypical observation of EMIC wave-driven electron precipitation at ELFIN with accompanying wave measurements on a conjugate platform. Figure 2 shows ELFIN A measurements in the same format as Fig. 1. The event occurred on 2021-04-29 in the outer radiation belt at dawn (where EMIC waves are also observed quite often, see (Zhang et al. 2016c)). The energy-and time-localized precipitation ratio (Panel (c)) at ∼13:21:05 UT marks the relativistic precipitation burst of interest. It lasted anywhere from 7 to 14 s (2.5-5 spins) the range depending on the intensity used for its definition. That ratio peaked in energy at ∼ 1 MeV, while both precipitating and (locally) trapped fluxes increased at the same time. These are all indicative of EMIC wave scattering Ni et al. 2015a;Mourenas et al. 2016). In particular, whistler-mode chorus waves (otherwise also abundant in this region of space) preferentially scatter lower energy electrons (see energy distributions of EMIC and whistler-mode wave scattering rates in Summers et al. 2007;Shprits and Ni 2009)) and cannot be responsible for these observations. Note that Panels (d) and (e), which show pitch-angle spectrograms during the event, demonstrate that the precipitation was evident far from the edge of the loss cone, close to the downgoing direction (whereas upgoing electons are interpreted as reflected electrons from the atmosphere below). The pitch-angles that enter the precipitating flux energy spectro-grams are, by selection, centered at >22.5 degrees from the edge of the loss cone, providing clean separation of precipitating and trapped fluxes. About 35 min prior to these ELFIN A observations at L ∼ 5, the equatorial THEMIS E spacecraft traversed the inner magnetosphere approximately in the radial direction away from Earth at a similar MLT as ELFIN (Fig. 3, Panel (a)) and detected intense wave activity at 0.5 -1 Hz frequency. The observed waves were propagating near-parallel to the background magnetic field (Panel (b)), had left-hand circular polarization (Panel (c)) and were seen to peak just below the H e + gyrofrequency (dashed lines, calculated using the local magnetic field), allowing us to identify them as H e + -band EMIC waves. The waves were seen in a locally measured magnetic field of (−66, 48, 227)nT in GSM coordinates, when THEMIS E was about 2000 km away from the magnetic equator. They were observed near the plasmapause: a region of density gradient (Panel (d)) exhibiting more than an order of magnitude change throughout the event (from 500 to 15 per cm 3 ) and short-scale density variations of a factor of 2 on timescales as short as tens of seconds. The EMIC wave emission was seen from L ∼ 4.65 to L ∼ 4.88, i.e. close to the L-shell where ELFIN detected the ∼ 1 MeV electron precipitation. The narrow L-shell range of EMIC observation at THEMIS E, L ≈0.23 (a 6 min duration x 4 km/s THEMIS E's radial velocity) is traversed by ELFIN A at 7.9 km/s in 9.6 sec, or 3.4 spin periods. This is similar to the duration of ELFIN-A's observation of precipitation (7-14 s, as discussed earlier), and explains the short lifetime (the bursty apparent nature) of the electron precipitation observed from ELFIN's ionospheric vantage-point. The EMIC wave emission was not only localized in L-shell but, in this particular case, was likely also relatively short-lived: THEMIS A and THEMIS D which traversed the same region of space as THEMIS E about 50 and 70 minutes prior, respectively, did not observe the emission. Neither did ELFIN A, which traversed the same region in its prior and subsequent orbits (at 02:40 and 05:40 UT) observe similar type of precipitation. This limits the likely duration of the EMIC wave-driven precipitation to <1.5 hours. As we shall see in the following subsection, in other instances the precipitation can last a lot longer (many hours), confirming our assumption of the spatial interpretation of the EMIC wave power variability at THEMIS E.
It is instructive to examine several features of the EMIC waves in space that influence the precipitation signatures at ELFIN A. We continue with our assumption that at least the gross features of plasmaspheric density gradients and EMIC wave power last longer than a few minutes (their crossing time by THEMIS E) and are organized by the plasmapause in azimuthal sheets. We also note that the measured flow velocity and ExB velocity in the radial direction at THEMIS E are both 4 km/s consistent with the spacecraft velocity through the plasma. (Corotation, about 2.3 km/s in the azimuthal direction, is ignored here as it is both smaller than the satellite velocity and because the structures of interest vary mostly in the radial direction.) Examining the time-series EMIC wave data in one component on the plane of polarization (Panel (e)) we see that the emissions consist of bursts of <30 s duration (<23 wave periods) including transient bursts (such as the one around 03:35:25UT marked by an arrow) as short as 7 s (or 5 wave periods). The crossing time of such 7-30 s structures (if spatial at the equator) by ELFIN A in the ionosphere (scaled from the estimates in the previous paragraph) would be 0.19-0.8 s which is much smaller than ELFIN's nominal spin period (2.8 s), and around 1.-4.5 azimuthal spin sectors (each of 16 spin sectors lasts 0.175 s). Examination of the longer waveforms in further detail (Panel (f)) reveals phase skips every 5-15 s, or 4-11 wave periods (marked by the arrows and their separation times), mapping to 0.13-0.40 s, or 0.75-2.3 spin phase sectors at ELFIN. These EMIC wave structures (jumps in coherence and amplitude) from 5-30 s, corresponding to 20 -120 km at THEMIS, are smaller than the H e + inertial length (d H e+ ≈124-264 km in this case for Observations from THEMIS E at the equator at an MLT and UT near those of the science zone crossing by ELFIN A depicted in Fig. 2. Panels (a-d) show about an hour of data centered around an EMIC wave emission: power spectral density of the magnetic field measured by the fluxgate magnetometer, FGM, instrument (a), wave normal angle (b), ellipticity (c), electron density inferred from the spacecraft potential computed on-board by the electric field instrument, EFI, and processed on the ground using the measured electron temperature by the electrostatic analyzer, ESA, instrument (d). The black dashed line is the H e + gyrofrequency. Panel (e) shows ∼2.5 min of data from a single magnetic field component in a field-alignedcoordinate system (FAC), B X,F AC . It is oriented perpendicular to the average magnetic field direction (hence, near the plane of polarization) and lies on a plane also containing the sunward direction. Panel (f) is an expanded, ∼1 min long, view of the same quantity as above it. Arrows in Panels (e) and (f) are discussed in the main text N e ≈ 50-100/cc and a H e + number density ∼ N e /16) and the thermal proton gyroradius (100 km for 60 keV H + ) and much smaller than the expected (parallel) EMIC wavelength (λ ≈ 4πd H e+ ). They may be due to oblique propagation and interference of waves from different source locations in the presence of plasma density inhomogeneity of a comparable gradient scale (Panel (d)). As bouncing and drifting energetic electrons near cyclotron resonance encounter these EMIC wave structures, both along and across the field, they should experience short-lived coherent interactions with an ensemble of waves of a range of wavelengths, amplitudes and phases that, for small amplitudes, would appear like turbulence. (Even if wave-field temporal variations occur and are partly responsible for the wave observations, above, the energetic electrons stream and drift so fast across them that they are effectively stationary in the Earth frame.) Their resultant interaction with electrons would then be describable by quasi-linear theory. Field-aligned electrons are organized on drift shells with equatorial cross-sections encircling Earth but centered towards dawn, since field lines are stretched out further from Earth at post-midnight. Conversely, waves excited due to ion cyclotron resonance near density gradients (like the plasmapause) are arranged on distorted drift shells with circular equatorial cross-sections having centers near Earth but displaced towards dusk due to the plasmasphere buldging toward dusk. The intersection of resonant electron drift-shells and plasmaspheric density enhancements filled with EMIC waves would be sheet-like structures: longer in azimuth and thin in radial distance. A good fraction of that layer of waves could be also interacting with resonant drifting electrons, resulting in precipitation over several spins. We would expect the EMIC wave-driven electron precipitation to be organized in thin azimuthal layers of thickness roughly consistent with the aforementioned variations in wave power at THEMIS E over a 6min in duration, i.e., over a L ≈0.23 as discussed above. The significant variations in power on spatial scales 7-30 s at THEMIS E could result in abrupt enhancements in near-field-aligned equatorial fluxes on that time scale. Mapping at ELFIN A to 1-4.5 spin sectors, they can be either inside the loss cone or outside it (both are quite close to the equatorial loss cone), depending on ELFIN's spin phase and look direction. Hence the perpendicular-to-trapped flux ratio would be timealiased at ELFIN, exhibiting large increases (even above unity) or decreases due to abrupt changes in equatorial wave power, even when that power is consistent with the quasi-linear regime. Two arrows in Fig. 1(e ), of the event discussed in the previous subsection, show instances when the average flux in the loss-cone dominates the precitating-to-trapped flux ratio (the ratio exceeds one) due to flux enhancements in a single sector -clearly aliased. In fact, because there are more spin-phase sectors in the loss cone than outside (typically 6 versus 4 in each full spin) a random distribution (in time) of 1-5 sector-long flux enhancements would result in more flux ratios exceeding one than below one, statistically. That bias can be normalized away, however, in statistical studies. Case studies can rely on the consistency between consecutive trapped fluxes, twice per spin, to ensure aliasing has been minimized (e.g., Fig. 1 in Zhang et al. (2022)). Figure 4 shows energy-spectra of the number flux of precipitating and trapped electrons, averaged over the 6 s (four ELFIN A half-spins) when the dominant precipitation attributed to EMIC waves occurred, in the event of Fig. 2 that was discussed earlier in this subsection. The measured spectrum of the precipitating flux (dotted thin red line) shows a peak near 1 MeV. This peak is even more pronounced when the average of the precipitation determined from data collected 6 s before and 6 s after the dominant precipitation interval (the trend, dashed thin red line) is removed, revealing the net contribution to the precipitation from just the EMIC waves (solid thick red line). (We interpret the trend as most likely due to low-level hiss waves.) The ratio of the average measured precipitating-to-trapped flux (not shown) also peaks near 1 MeV at a high value, ∼ 0.4, consistent with the color spectrogram in Fig. 2(c), that depicts this ratio for individual half-spins peaking at ∼ 0.7. After Trends are average fluxes from 6 s immediately before and 6 s immediately after the strongest precipitation interval detrending the average of the precipitating flux in Fig. 4, the detrended precipitating flux has an even more clear peak near 1 MeV (see the solid thick red line), corresponding also to a stronger peak near 1 MeV for the detrended precipitating to un-detrended trapped flux ratio appropriate for comparisons with quasi-linear theory (Kennel and Petschek 1966). The detrended precipitating and trapped fluxes have a very similar energy spectrum peaked near 1 MeV (compare solid thick red and blue curves), suggesting that EMIC waves may be responsible for both flux increases compared to the trend, possibly through nonlinear transport from higher pitch-angles . As we will see later, in about 50% of the time in our database of EMIC events the un-detrended flux ratio actually peaks and exceeds one above 1 MeV. However, given the significant spatial variability of the EMIC wave field and its electron interaction region in the magnetosphere, the ratio exceeding unity must be viewed with caution, because of the temporal aliasing effects arising from latitudinally narrow regions of precipitation lasting a few spin sectors, which prevents the precipitating and (locally) trapped fluxes to be measured simultaneously.
Nevertheless, it is evident from the above discussion that the precipitating flux during EMIC events exhibits a strong peak and a high precipitating-to-trapped flux ratio in the 0.5 − 1.5 MeV energy range, making this a hallmark of EMIC wave-driven precipitation in spectra that are well-resolved in energy and pitch-angle.
Assuming that there exist cases when time aliasing does not affect the precipitating-totrapped flux ratio, that ratio exceeding one would signify the presence of nonlinear EMIC wave-relativistic electron interactions. This is because a ratio >1 cannot be explained by quasi-linear diffusion, which has an upper limit of precipitation, the strong diffusion limit (see (Kennel and Petschek 1966)), that necessitates that precipitating and trapped fluxes be equal. However, nonlinear resonant interactions can indeed result in loss cone fluxes that are greater than permitted by quasi-linear theory. This could be due to the phase trapping effect (Kubota et al. 2015;Kubota and Omura 2017;Grach and Demekhov 2020;Grach et al. 2021). In this mechanism, very intense EMIC waves can interact resonantly with electrons initially located well above the loss-cone edge and transport them in phase space directly into the loss cone. Such nonlinear trapping results in a large pitch-angle change, 20 • − 40 • , during a single resonant interaction. Therefore, electrons with large equatorial pitch-angles exhibiting increasingly larger flux at fixed energy (due to the typically strong perpendicular anisotropy of relativistic electrons (Ni et al. 2015b)) can be transported all the way into the loss cone directly and without simultaneously enhancing the trapped flux near the edge of the equatorial loss-cone that corresponds to the only trapped electron population visible at ELFIN's altitude. This has been already demonstrated for a case study of EMIC wave-driven precipitation by Grach et al. (2022) who used simultaneous equatorial observations of the EMIC waves and modeling of the wave-particle interactions to demonstrate the precipitation ratio should exceed unity, as was indeed observed on ELFIN. In that case the ratio exceeded unity for 3 consecutive spins (albeit at different energies) which bolsters the case for nonlinear scattering. Similar case-by-case studies of the details of the ELFIN particle distributions are needed, hand-in-hand with modeling, to verify the presence of nonlinear effects and separate them from temporal aliasing. Statistical studies of the problem can either rely on multiple consecutive spins with similar signatures, or probabilistic analysis of the cluster of strong precipitation events (ratio >1) after removal of trend and biases. Thus, the fine energy and pitch-angle resolution of ELFIN's energetic electron measurements from low altitude allow us to also identify and study the properties of nonlinear interactions of EMIC waves and relativistic electrons, which are likely important at times of intense EMIC waves.
We see that although EMIC wave bursts are often relatively short-lived and spatially limited (Blum et al. 2016;Meredith et al. 2014), their intensity, either inferred from spacecraft wave statistics or surmised from ELFIN's precipitation measurements, implies that their time-and MLT-averaged wave power could be sufficient to lead to strong electron losses at high energy Drozdov et al. 2020;Ross et al. 2022), especially during long-lasting events extending over many hours in MLT during active times (Engebretson et al. 2015).

Evolution of Long-Lasting EMIC Wave-Driven Electron Precipitation
As we saw in the previous section, relativistic electron precipitation driven by EMIC waves is often localized in L-shell due to the localization of the EMIC wave excitation and distribution in the magnetosphere (Blum et al. 2016. From ELFIN's low-altitude vantage point, such a spatial localization is evidenced as temporally localized, transient or bursty precipitation (as in Fig. 2 and in events to be discussed subsequently in Sects. 4.1, 4.2). However, EMIC waves can occasionally persist for hours, as revealed in prior studies from combinations of ground-based and multi-spacecraft observations (Engebretson et al. 2015;Blum et al. 2020). From its low-altitude, ∼ 90 min period orbit, ELFIN can traverse the same MLT and L-shell region repeatedly, and thus it too can identify long-lasting EMIC waves, as well as monitor and study their gradual evolution in intensity and latitude using their clear electron precipitation signatures. Figure 5 shows such an event over three ELFIN B orbits (lasting ∼ 3 hours).
At first (Fig. 5, Panels (a-c)), ELFIN B observed EMIC wave-driven electron precipitation at 02:42:50 UT, at L-shell∼ 5.4 (see Panel (j)), lasting ∼ 8.4 seconds ( L ≈ 0.1). This is evident based on the precipitating-to-trapped flux ratio, showing that the most efficient precipitation lies in the range [0.5, 2] MeV, where that ratio is ∼1 -precipitating fluxes at these energies are comparable to trapped. An orbit later (Fig. 5, Panels (df)), ELFIN B passed over approximately the same MLT (at a distance of MLT ≈ 0.1) and observed again strong precipitation (at 04:15:25-04:15:35 UT) at L-shell∼ 5.2, over a somewhat broader energy range, now [0.4, 2] MeV, and lasting much longer, ∼ 22.4 seconds ( L ≈ 0.25). The third time around (Fig. 5, Panels (g-i)), EMIC wave-driven precipitating fluxes were seen again at ELFIN B (at 05:48:20-05:48:35 UT) at L-shell∼ 5.0. They were comparable to trapped fluxes and extended over an even greater range in energy, [0.3., 1.5] MeV, and L-shell ( L ≈ 0.5). There is a small MLT evolution between the three orbits (from the first to the third crossing the MLT of the precipitation event changed by MLT ≈ 0.35), but this is well within the expected EMIC wave azimuthal extent of many  Figs. 1 and 2, i.e., they are energy spectrograms of trapped fluxes, precipitating fluxes and precipitating-to-trapped flux ratio hours in MLT in the equatorial magnetosphere (Blum et al. 2016. The location of the center-time of the emissions moves closer to Earth in each encounter (from L ≈ 5.4, to 5.2 to 5.0 in IGRF; the absolute values and their variation are, naturally, model-dependent, but nonetheless dominated by the Earth's internal magnetic moments, to zero order).) Given their similarity in L-shell, we conclude that these observations are likely due to continuous EMIC wave activity from roughly the same region in space, where either the EMIC wave free energy source (e.g. drifting ions) or the density gradient that enables resonant interactions with that source was evolving in time (moving closer to Earth and expanding in radial extent, in this case). Therefore, using ELFIN's repeated passes over a long-lasting (∼ 3 hours) relativistic electron precipitation event also allows us to infer characteristics of the EMIC wave spatial location, intensity, extent, and temporal evolution at timescales of an orbit period or, on occasion, even faster (due to the occasional availability of data from two satellites or from two science zones at the same MLT, in the north and south hemispheres, on each orbit). The longevity of some EMIC wave events coupled with their anticipated presence over many hours in MLT, near the plasmaspheric boundary where conditions for their excitation prevail, suggest that EMIC waves may be significant contributors to outer belt electron losses during active times.

Studies of EMIC Wave-Driven Precipitation with ELFIN and Its
Conjunctions with Ancillary Datasets Figure 6 shows an overview of ELFIN A observations on 2 November 2020. There is a clear peak of electron precipitation above 1 MeV at ∼15:19:00 UT. The precipitating-to-trapped flux ratio is ∼ 1 for energies ∼ 1 − 3 MeV and decreases with decreasing energy at energies < 1 MeV. Enhanced precipitation is observed during 6 consecutive ELFIN half-spins (∼ 8.5 seconds), the temporal localization again likely being due to ELFIN A's crossing of flux-tubes mapping to a spatially localized region of EMIC waves near the equator at L-shell∼ 5.5. The typical scale of such wave emissions in the radial direction at the equator is about 0.5R E (Blum et al. 2016, consistent with the observed L-shell extent in this event of L ≈ 0.4. ELFIN A's trajectory projections to the north and south hemisphere, shown in Fig. 7, demonstrate the existence of several nearby ground-based magnetometer stations. These could potentially provide high-resolution magnetic field measurements that can reveal EMIC waves of interest during this pass. (The time of strong precipitation at ELFIN A is denoted by a thick trace superimposed on its projections in that figure.) Thus, although there is no equatorial spacecraft conjugate to ELFIN A at this time, it is still possible to obtain information on the presence and properties of the EMIC waves associated with the observed precipitation using such stations (Usanova et al. 2014), if the EMIC waves managed to propagate to the ground. And indeed, stations PG3, PG4, and SPA measured low frequency, banded emissions consistent with EMIC waves. We elect to work with data from SPA (Fig. 8) which was closest to ELFIN A at the time of interest. Using the T89 magnetic field model (Tsyganenko 1989) we plot the equatorial helium gyrofrequency for the L-shells pertaining to these stations. The observed wave emission, exhibiting an upper limit f max just below the H e + gyrofrequency, f cH e , evidently corresponds to helium band EMIC waves. To further confirm the resonance condition for these waves, we use an empirical plasma density model (Sheeley et al. 2001) to estimate the equatorial plasma frequency f pe ; its ratio to the electron cyclotron frequency f ce (from the aforementioned use of T89) is about 8.5 − 9. This is typical of the plasma trough, where ELFIN A was located during the subject precipitation (between the plasmasphere and the inner edge of the plasma sheet, marked atop the trapped electron energy spectrogram in Fig. 7(a) using the same criteria as explained earlier for the events of Fig. 1). If equatorial EMIC wave measurements are available in individual case studies, even in the absence of direct measurements of ion composition, an estimate of the H e + to H + concentration ratio can be inferred from the stop band between f cH e and H + -band wave power Qin et al. 2019). In the absence of conjugate equatorial wave measurements, we adopt here a statistical treatment, using typical H e + concentrations, ≈ 5% and 10% Ross et al. 2022). We find that the above f pe /f ce ratio and the observed value of f max /f cH e ≈ 0.97 − 0.99 (Fig. 8) can result in a theoretical estimate of the minimum resonance energy of EMIC waves. To demonstrate this, we plot the theoretical minimum resonance energy for a wide range of f max /f cH e (vertical axis) and f pe /f ce (horizontal axis) in Fig. 9. We then denote the range of the latter two parameters expected from observations as the gray areas between two horizontal and two vertical lines. The intersection of those two gray areas defines the region in this two-parameter space that corresponds to the expected range of parameters and solutions for the resonance energy as depicted by the plot's contours. This resonance energy estimate, Fig. 7 ELFIN A projections to the ionosphere in the north and south for the event of Fig. 6. Diamonds and asterisks mark the start and end times of the trajectories; crosses are 1 min tickmarks; thick traces denote times of intense relativistic electron precipitation identified from Fig. 6(c) as a putative EMIC wave-driven precipitation event Fig. 8 Magnetic field spectra from the ground-based station at SPA in the Antarctic. As shown in Fig. 7, SPA is in close conjunction with ELFIN-A during this event. Superimposed in the spectra are H e + and O + equatorial gyrofrequencies (horizontal lines) using the magnetic field at their equatorial projection, inferred from the T89 model. The vertical magenta lines mark the time interval of ELFIN A' science zone crossing during this event between 1.5 and 2.5 MeV, is quite close to the moderate and strong electron precipitation energies observed at ELFIN A.

Confirming H e + -Band EMIC Wave Resonance with Electrons
To make the last point clear, we marked in yellow and red in Fig. 9 the energies where moderate and strong precipitation was observed, respectively. This determination was based on the ratio R, depicted in Panel (a), being 0.5 > R > 0.3 and R > 0.5, for moderate (yellow) and strong (red) precipitation, respectively. Transferring the energy range of these two categories onto the contours of Panels (b) and (c) readily depicts our observational assessment of moderate and strong precipitation within the resonance energy contours. We can thus see that the expected resonance energy based on the most likely values of the two aforementioned parameters, f max /f cH e and f pe /f ce (the intersection of their respective grayed area bounds), overlaps with the observed moderate and strong precipitation (yellow and red highlighted contours), as one would expect from quasi-linear theory. and of the f pe /f ce ratio for a H e + ion concentration of 5%; and same as the panel above, but for a 10% H e + concentration. Red and yellow colors depict the electron energy ranges for which ELFIN A measured strong precipitation and moderate precipitation (R>0.5 and 0.5>R>0.3, respectively). The plot shows that resonance energies exhibiting moderate and strong precipitation at ELFIN are consistent with the range of parameters f max /f cH e , f pe /f ce inferred from in-situ measurements at conjugate platforms (at the intersection of the corresponding grayed areas), for a reasonable range of H e + concentrations Figure 10 shows an overview of ELFIN A observations on 6 Dec 2020 exhibiting a clear peak of electron precipitation at E > 1 MeV around 20:20:40 UT (marked by a left-leaning arrow), a putative EMIC wave scattering event. At that time, the precipitating-to-trapped flux ratio peaked (at 1) near ∼ 1 − 2 MeV, and decreased with decreasing energy at < 1MeV. Strong relativistic electron precipitation at energies ∼ 1 MeV was also seen at four earlier times in the same ELFIN science zone crossing (denoted by down-pointing arrows). Each of the five instances of intense precipitation was transient, lasting for a few (2-8) spins (∼ 5 − 20 seconds), had a peak ratio 1 at ∼ 0.5 − 2 MeV and exhibited a decreasing ratio with decreasing energy below its peak.

Confirming H + -Band EMIC Wave Resonance with Electrons
As in the event we examined previously, ELFIN A's magnetic projections to the north and south hemispheres (Fig. 11) reveal that several ground-based stations were magnetically conjugate to ELFIN A, and might be able to provide magnetic field data that can ex- Fig. 10 Overview of ELFIN A observations during a H + band EMIC wave-driven event, on 6 December 2020, in a format similar to that of Fig. 1 plore the presence and properties of EMIC waves during this interval. And indeed, stations KEV, OUL, and SOD of the Finnish pulsation magnetometer network do show evidence for narrow-banded waves, likely EMIC waves, that started abruptly with a broadband burst, probably due to a nightside injection. In Fig. 12, we elect to show data only from SOD, which was closest to ELFIN A during the 20:20:40 UT burst. Although the wave emission was not well defined spectrally on the ground at the time of ELFIN A's science zone crossing (demarcated by the two vertical magenta lines in that figure), it became clearly defined a few minutes after (at ∼20:25:00 UT): it had peak-power and maximum frequencies (f peak and f max ), that varied (rose) considerably over the next 2 hours, both within the range 0.2-0.7 Hz. T89 magnetic field model-based O + , H e + and H + gyrofrequencies at the equatorial location conjugate to SOD are over-plotted in the above figure. They reveal that the observed wave intensity peaked below the expected equatorial H + gyrofrequency (≈ 0.7H z) and was suppressed as that frequency was approached from below. The lack of a local power minimum (a stop band) near the helium gyrofrequency and just above it leads us to surmise that there is likely negligible helium concentration in the wave source region Chen et al. 2019), and that the observed emission is likely due to an equatorial source of H + EMIC waves. We take note of the delay in appearance of the EMIC wave on the ground relative to the broadband burst that initiated the wave activity and also relative to the overhead passage of ELFIN A. We interpret this delay as due to the fact that it takes some time for the duct that facilitates unimpeded EMIC wave propagation to the ground to establish itself along the entire flux tube. The broadband waves, in the Pi2 and Pi1 range, are (at least part of the way) compressional and Alfvénic and arrive to the ionosphere from  Magnetic field power spectral density (arbitrary units) from the Finland ground-based station at SOD, located as shown in Fig. 11. Superimposed H + , H e + and O + equatorial gyrofrequencies (horizontal lines) using the equatorial magnetic field conjugate to these stations. The magenta vertical lines bracket the time interval of ELFIN A's science zone crossing during this event the equatorial magnetosphere without the need for ducting. Thus, due to the need for duct formation for EMIC waves to propagate far from the equator, ground EMIC wave detection is delayed relative to the equatorial injection that likely initiated the hot anisotropic ions, the equatorial EMIC wave emission, the associated electron precipitation, the initial broadband waves on the ground, and the density perturbation that eventually led to the duct.
Equatorial observations of spacecraft potential from THEMIS-A, -D and -E (Angelopoulos 2008), which crossed the same L-shell as ELFIN A two hours later and a few MLT hours away (at ∼20 MLT), provide estimates of the equatorial plasma density (and plasma frequency). At the times when the locally measured magnetic field at THEMIS was ∼ 45 − 50nT, consistent with the aforementioned H + gyrofrequency, the inferred equatorial density from the spacecraft potential was 2.5-3.0 cm −3 , such that the inferred equatorial f pe /f ce ratio was about 10 − 12. Together with the observed value of f max /f cp ∼ 0.6 − 0.9 from the ground-observatories, this gives a theoretical estimate of the wave minimum resonance energy for electrons of ≈ 0.3 -2 MeV as marked by the cross-section of the hori- , and of the f pe /f ce ratio for 0% helium concentration; and same as the panel above but for a 2.5%H e + concentration. Red and yellow colors depict the electron energy ranges for which ELFIN A measured strong and moderate precipitation (R>0.5 and 0.5>R>0.3, respectively). The plot shows that resonance energies exhibiting moderate and strong precipitation at ELFIN are consistent with the range of parameters f max /f cp and f pe /f ce inferred from in-situ observations at conjugate platforms (f max /f cp ∼0.6-0.9 and f pe /f ce ∼10-12) for a reasonable range of H e + concentrations zontal and vertical highlighted regions in Fig. 13. This estimate is quite close to the energies of moderate and strong electron precipitation observed by ELFIN, depicted in the same figure with red and yellow areas within the wider f max /f cp versus f pe /f ce parameter space. Our ELFIN observations of peak precipitation combined with ground based estimates of the wave frequency and equatorial spacecraft estimates of plasma density are therefore consistent with a H + band EMIC wave-driven interpretation in the context of quasi-linear theory.

EMIC Wave-Driven Precipitation in the Context of TEC Maps
As discussed above, the equatorial electron density is one of the most important parameters affecting EMIC wave generation and resonance with relativistic electrons. When that density is large, as in the case of plasmaspheric plumes (Fraser et al. 2005;Usanova et al. 2013;Halford et al. 2015), the EMIC wave resonant excitation by hot anisotropic ions is favored (Horne and Thorne 1993;Chen et al. 2009Chen et al. , 2010, while the energy of electrons resonant with EMIC waves also decreases down to energies of peak flux in the outer radiation belt, 0.5-2 MeV, leading to intense precipitation of such relativistic electrons (Summers and Thorne 2003;Summers et al. 2007). Thus, magnetospheric density enhancements, collocated with azimuthally drifting energetic ions, have been commonly reported as favorable sites for EMIC wave growth and relativistic electron precipitation. Such sites are expected around boundaries between the nominal ring current region (filled by hot ions injected from the plasma sheet Runov et al. 2015)) and the plasmasphere (Thorne and Kennel 1971;Horne and Thorne 1993;Chen et al. 2010). It is challenging, however, to definitively establish the connection between EMIC wave-driven electron precipitation and density enhancements solely based on fortuitous conjunctions between low-altitude and equatorial spacecraft. This is because density ducts (spatially limited enhancements) can be highly-localized and the equatorial spacecraft in those fortuitous conjunctions are not always in the optimal location to reveal the spatial profile of the pertinent density enhancement responsible for wave excitation. Mapping uncertainties due to imperfect magnetic field models further complicate such conjunction-based studies.
An alternate approach to exploring the density variations at play during ELFIN measurements of relativistic electron precipitation is to use total electron content (TEC) measurements of the ionosphere. These are obtained using the phase delay of radio signals transmitted from Global Navigation Satellite Systems (GNSS) satellites (moving along circular orbits at an altitude of ∼ 20, 000 km) to ground-based receivers. That phase delay allows estimates of the altitude-integrated electron density (Davies 1965;Vo and Foster 2001) along the line-of-sight. Data from multiple propagation rays during a finite time-step (from tens of seconds to minutes) at a wide range of propagation angles are assimilated through a tomographic reconstruction process to produce TEC maps over a wide area on the ground. TEC maps reveal well the spatial and temporal variations of plasmaspheric density (Heise et al. 2002;Belehaki et al. 2004;Lee et al. 2013), including those at adjacent density enhancement structures, like plasmaspheric plumes (Foster et al. 2002;Walsh et al. 2014). Moreover, TEC data also reveal quite well the nightside region of enhanced hot (ring-current energy) ions, which, is known to the ionospheric community as the mid-latitude ionospheric trough (MIT; see discussion on its formation in Aa et al. (2020)). This topside ionosphere, subauroral latitude region is recognized in satellite data by its enhanced ionospheric electron temperature and its prominent density reduction. The latter is also commonly captured in TEC data (see: (Yizengaw and Moldwin 2005;Weygand et al. 2021)). Because ring current ions can provide free energy for EMIC wave generation in the pre-midnight sector, a significant fraction of EMIC wave-driven precipitation is observed in the nightside region (Capannolo et al. 2022(Capannolo et al. , 2023Carson et al. 2021), precisely where the MIT develops. Using TEC data to study the correlation between the TEC-derived density gradients with relativistic electron precipitation can be quite advantageous: The magnetic projection of the low-altitude ELFIN satellites downwards, onto ionospheric TEC maps is highly accurate compared to upwards, to the magnetic equator. Moreover, the large availability and wide spatial coverage (in MLT and L-shell) of TEC observations provides a large dataset with which a correlation between precipitation events and plasma density boundaries can be investigated. In the following, we present first results from such studies and exemplify advances that can be made in the future using a similar approach. Figure 14 shows two relativistic electron precipitation events on successive orbits of ELFIN A (∼ 90 min apart, at the same MLT). (Note that Panel (e), showing MLAT and L-shell, applies to both events, as the event start/stop times have been selected to match each-other in the coordinates depicted on this panel.) The events show clear signatures of EMIC wave-driven scattering at around 11:30:00 UT and 13:02:30 UT, respectively. At those times, ELFIN observed strong, transient increases in trapped and precipitating electron fluxes at an L-shell ∼ 6.5. The precipitating-to-trapped flux ratio peaked at relativistic energies (∼ 1MeV). Both bursts were observed between the inner edge of the plasma sheet (the earthward boundary of isotropic plasma sheet fluxes in the energy range of < 200 keV, at ∼11:29:23 UT and ∼13:02:05 UT in the two crossings, respectively) and the plasmasphere (where trapped fluxes of energy ∼200 keV fall below those at >500 keV energy due to effective scattering of <500 keV electrons by plasmasheric hiss (Ma et al. 2016a;Mourenas et al. 2017), at and after 11:30:55 UT and 13:03:25 UT in the two crossings, respectively). The precipitation was likely spatially localized ( L ∼ 0.4 for both events). The similarity of the EMIC-driven precipitation signatures in two ELFIN orbits suggests here too that EMIC wave generation persists for at least 1.5 hours (as in the cases presented earlier, in Fig. 5, and as was previously reported by Engebretson et al. (2015), Blum et al. (2020)). Figure 15 shows projections of ELFIN A orbits onto TEC maps in the northern hemisphere (provided by MIT Haystack via the Madrigal database (Rideout and Coster 2006;Coster et al. 2013;Vierinen et al. 2015)). The spatial resolution of the TEC map is 1 by 1 degree. The images are geographic projections but the magnetic longitude is denoted by the two overlaid MLT meridians, for MLT = 0 and MLT = 6 respectively. The ELFIN A orbits map magnetically near the midnight meridian, slightly post-midnight (MLT ≈ 0-1). This is to the west of the westernmost boundary of the MIT, which can be identified as the blue regions of TEC depletion. The intense TEC values at the dayside are due to enhanced TEC due to the sunlight. Note that during the second ELFIN A orbit, at ∼13:02:00 UT, the westernmost boundary of the MIT is not as clearly visible, due to statistical uncertainties from lack of high-latitude TEC stations in the Northern Pacific. This TEC quality degradation occurs as the Earth rotates counterclockwise, viewed from the north, by 22.5 degrees between the two successive orbits, and the American sector moves eastward. The ELFIN A orbit and the sunlit-enhanced TEC at the dayside, which are both nearly-fixed in magnetic local time coordinates, appear to rotate clockwise (westward) by the same amount in the geographic coordinate system of this figure.
The MIT region's western edge evolves as the season changes between the summer solstice and the autumnal equinox, moving from dawn to dusk, according to previous statistical studies (Aa et al. 2020). As expected for the season of our event, close to autumnal equinox, the MIT is observed in Fig. 15 preferentially at post-midnight to dawn (MLT ∈ [0, 6]). It has also been previously established through statistical investigations of the MIT in conjunction with near-equatorial spacecraft measurements that at each longitude the minimum of TEC viewed as a function of latitude within the MIT corresponds to the equatorial plasmapause location (Shinbori et al. 2021;Heilig et al. 2022), whereas the westernmost MIT boundary (roughly along a magnetic meridian) is magnetically conjugate to the plasmaspheric plume (Heilig et al. 2022). Thus, both ELFIN A orbits map magnetically near the westernmost MIT boundary. This boundary should be interpreted as the overlap between ring current ions (resulting in the MIT's formation in the first place, (Aa et al. 2020)) and the cold plasma density region associated with the plasmaspheric plume (Heilig et al. 2022). Finally, the poleward TEC boundary of the MIT maps to the inner edge of the plasma sheet, as it is attributed to local, ionospheric electron enhancements brought-about by plasma sheet electron precipitation (Aa et al. 2020).
During its 6-min-long science-zone crossing in each orbit, ELFIN A was moving equatorward and was initially in a region of localized TEC enhancement associated with plasma sheet electron precipitation, as also evident by the electron spectrograms in Fig. 14, Panels (a) and (a ), and denoted by horizontal blue bars above those panels. It traversed the inner edge of the plasma sheet at approximately 11:29:23 UT and 13:02:05 UT, respectively Fig. 14 Overview of ELFIN A observations on 10 September 2020, which is being examined together with concurrent TEC observations. Two consecutive northern hemisphere, equatorward science zone crossings at midnight/pre-midnight (MLT ∼0.75) are depicted in three panels per crossing (a-c; and f-g), arranged in time to have a common L-shell and magnetic latitude, shown in Panel (e). Each crossing's three panels have the same format as the top three panels in Figs. 1 and 2, i.e., they show energy spectrograms of trapped fluxes, precipitating fluxes and precipitating-to-trapped flux ratio. A fourth panel for each crossing (Panels (d) and (h), respectively) show the TEC values at the satellite ionospheric projection (after averaging over 3 degrees each in latitude and longitude, and over 20min in time) (the times when the precipitating-to-trapped flux ratio stopped being isotropic and trapped fluxes exceeding background levels were still below < 200keV those panels). Observations of EMIC wave-driven precipitation, which ensued and are demarcated along its orbit projections in Fig. 15 by thick magenta lines, appear to map into the green or blue TEC region, between the poleward and equatorward MIT boundaries and near the west MIT boundary. The TEC values along the satellite tracks on Fig. 15 are transferred as line plots into Fig. 14 (see Panels (d) and (d )) to facilitate comparison with the ELFIN A precipitation signatures in the panel stacks right above them, and with the region identifications depicted by the horizontal color bars at the top of each stack. At ∼11:31:20 UT and ∼13:03:55 UT, the times of ELFIN A's crossings of the plasmapause, as inferred from in-situ measurements in Fig. 14 and shown by the black color bars, ELFIN A was located on the poleward edge of, but near the MIT's TEC minima in Fig. 14(d) and (d ) and also evident in Fig. 15. These minima are known to map to the equatorial plasmapause location (Shinbori et al. 2021;Heilig et al. 2022)) and therefore are consistent with our approximate identification of the plasmapause from in-situ energetic electron measurements. Therefore, comparison of TEC and ELFIN A-measured electron precipitation demonstrates that regions of EMIC wave-driven electron precipitation indeed correlate well with the TEC minimum and its poleward gradient within the MIT region, i.e., they are located just outside the plasmapause. Similar comparisons with TEC maps can provide two-dimensional ionospheric and plasmaspheric context information for ELFIN's locally trapped and precipitating flux measurements.

Statistical Results
To analyze the spatial and spectral properties of EMIC wave-driven precipitation statistically, we examined all ELFIN A&B observations from 2019-2021. Based on the aforementioned telltale signatures of EMIC wave-driven precipitation, we applied an operational definition for such a precipitation event to be a peak in the precipitating-to-trapped flux ratio at > 0.5 MeV. We thus identified ∼ 50 relativistic electron precipitation events similar to those discussed in the previous sections. We specifically excluded from the event database those which were within the expected location, or were consistent with, the electron isotropy boundary. The latter is distinctly evident by its location at the inner edge of the plasma sheet and the characteristic energy dispersion exhibited by the precipitating-to-trapped flux ratio (Imhof et al. 1977(Imhof et al. , 1979Yahnin et al. 1997); it has been recently studied extensively using ELFIN data in its own right (Wilkins et al. 2023) and in contrast to EMIC waves (Capannolo et al. 2023). The putative EMIC wave-driven precipitation events thus selected correspond to a total of ∼ 310 spins, with an average event duration of ∼ 6 spins or ∼ 17 s. They are most often located at L∼6 and have typical L ∼ 0.56, consistent with prior reports of the radial extent of the EMIC waves at the equator (Blum et al. 2016. This is in accordance with our earlier assertion (in Sect. 3.2) that the bursty nature of the subject precipitation at ELFIN is due to the radial localization of the EMIC waves and of their electron-interaction region in the magnetosphere.

Spatial Distribution of All EMIC Wave-Driven Events
Figure 16(a) shows that most events occurred at around L ∼ 5 − 7 consistent with their expected location within the outer radiation belt. The events occurred predominantly at the midnight and dusk sectors, even though the ELFIN database of science zone crossings has uniform coverage in local time. Their probability distribution in MLT -L space, shown in Fig. 16(b), reveals that while events are indeed most probable at pre-midnight and dusk at L ∼ 5 − 7, a second class of events is present at midnight, post-midnight and dawn at higher L-shells, L ∼ 8 − 12. This bares close resemblance to the distribution of EMIC waves in the magnetosphere (see e.g., (Min et al. 2012)).
All events have a large precipitating-to-trapped flux ratio, near one (Panel (c)). In fact, for almost half of the events that ratio exceeds unity. Full-spin resolution data have been used (thus trapped fluxes represent an average of flux accumulations from two look directions, one before and one after the precipitation measurement). Only points with relative error of the ratio, R, R/R < 50% were included herein. We thus find statistically that the above ratio exceeds unity for a very considerable fraction of all events. This suggests that nonlinear effects may be considerable in our database of EMIC wave-driven electron scattering process. However, because transient precipitation on a sub-spin resolution, lasting only a few spin sectors, may be common, careful consideration and removal of aliasing needs to  (d) and (e) depict theoretical estimates of the respective quantities based on statistical wave spectra (derived as discussed in text). Dashed and dotted blue lines in Panel (e) depict theoretical lower limits of the resonance energy estimated based on wave cyclotron damping at kc/ pi ∼ 1 and kc/ pi ∼ 2, respectively take place in a statistical or event-by-event analysis of this database in a future study to definitively identify events that are due to nonlinear scattering.
It is possible that few, some, or many of our events are due to nonlinear scattering due to spatially-localized or short-coherence scale, high-amplitude wave emissions. In fact, it has been known for a while that individual EMIC wave-packets often reach sufficient amplitudes to allow faster, nonlinear wave-particle interactions (Engebretson et al. 2008;Albert and Bortnik 2009;Pickett et al. 2010;Omura and Zhao 2012;Grach et al. 2021). However, despite the presence of such nonlinearly scattered, bursty precipitation events in our database, it is still quite likely that such events statistically conform to a diffusive treatment. The argument for this is an analogy with nonlinear whistler-mode chorus waves. Chorus wave-packets having large amplitude (thus, in the nonlinear regime) can have short wavepacket lengths or exhibit strong, random phase jumps, or both, resulting in phase decoherence. This allows a diffusive description of the resulting electron phase space density evolution, despite the large wave amplitudes at play Artemyev et al. , 2022Mourenas et al. 2022). And for EMIC waves too, it is likely that a prevalence of short EMIC wave-packets exhibiting strong and random wave phase jumps across or within them (e.g., see various examples of short packets in  and ; some examples have also been seen in Fig. 3 and discussed in Sect. 3.2) can permit a diffusive description of the scattering process despite the nonlinear nature of the scattering. In this approach the average precipitating fluxes and the average wave amplitudes incorporate the combined effects of linear and nonlinear resonant interaction regimes. We therefore proceed in our analysis of the precipitation energy spectra with a diffusive paradigm in mind, and study the minimum resonance energy, peak precipitation energy and the typical wave amplitude variation with frequency that correspond to the energy-spectrum of the precipitation, next. Figure 16(d) shows that the energy of the peak precipitating-to-trapped flux ratio E * lies typically in the range 1 to 3 MeV. We interpret E * as the energy of electrons resonating with the most intense EMIC waves. This is because the typical spectrum of EMIC waves on the ground or in space (as seen for example in Fig. 3(a) or Fig. 8) has a peak intensity at a certain frequency, f peak . Then, electrons of energy E peak resonating with the most intense EMIC waves at that frequency should also exhibit a peak in their precipitating-to-trapped flux ratio. The higher frequency part of the EMIC wave spectrum (at f > f peak ), has a lower wave intensity, but it will still lead to electron precipitation at energies below E peak , down to the minimum resonance energy E R,min corresponding to f max , the maximum wave frequency of appreciable wave power.
To better approximate the minimum resonance energy for the purpose of statistical studies of its spatial distribution, we define E * min as the half-way point in precipitation intensity below its peak. E * min is thus the energy at half-peak of the measured precipitating-to-trapped electron flux ratio at energies lower than that of peak precipitation (E * min < E * ). (For some EMIC spectra without appreciable wave power at frequencies higher than their peak intensity, it would be E * min = E * , but that is a rare occasion.) We note that warm plasma effects may limit scattering to frequencies below f max (energies above E * min ) but also that such effects may also limit propagation and suppress the average wave power to below its f max attainable in each individual event time series. Thus the half-way point appears as a reasonable energy selection in precipitation data to represent the theoretical E R,min expected from a statistical average of f max in available datasets.
In summary, E * is our measured estimate for the theoretical E peak corresponding to the point of maximum precipitation by EMIC wave-driven waves of peak wave power at f peak , and E * min is our measured proxy for the theoretical minimum resonance energy E R,min corresponding to the maximum frequency of appreciable wave power f max in statistical averages of such wave power observations. In other words, E * min is an estimate of the minimum resonance energy E R,min corresponding to significant wave-driven electron scattering toward the loss-cone. Figure 16(e) shows E * min for each spin in our database of events. Evidently, its average lies in the range of 0.5 to 1.5 MeV.
Both E * and E * min tend to decrease, on average, as L increases (despite the large scatter, which is in part due to uncertainties in L-shell determination). This is expected, due to the decrease of the minimum resonance energy for cyclotron resonance with distance from Earth: Since the dipole field intensity fall-off with distance (∼ 1/L 3 ) is faster than that of the square-root of the density (Sheeley et al. 2001), f pe /f ce increases with distance, and the minimum resonance energy for electrons (which is monotonically and inversely dependent on f pe /f ce , as discussed in Sect. 2.2) decreases with distance. We can see this behavior even in case studies with multiple EMIC wave-driven bursts, such as that of Fig. 10: the energy of peak precipitating-to-trapped ratio, i.e., the energy of the most efficient precipitation, decreases at progressively larger L-shells.
To assess the trends revealed in the averages shown, we compare them with estimates for E * and E * min determined independently, from published statistical averages of EMIC wave spectra providing peak-power and maximum frequency (Zhang et al. 2016c) and from a model-based estimation of the f pe /f ce ratio (Summers and Thorne 2003). Our independent estimates resulted in the dashed red lines in Panels (d) and (e). Specifically, our procedure was as follows: Using an empirical model for the plasmaspheric plume density n e ∼ 1300(3/L) 4.83 cm −3 at L > 4 (Sheeley et al. 2001), as appropriate for the dusk sector where most of our events were observed, and assuming that the resonance is mainly with hydrogen band EMIC waves (based on previous statistical observations Mourenas et al. 2016Mourenas et al. , 2017Zhang et al. 2016cZhang et al. , 2017), we get for the resonance energy of precipitating fluxes E[MeV] = [(1 + C/L 1.17 ) 1/2 − 1]/2. The coefficient C is to be determined from EMIC wave and background plasma characteristics. If conjugate equatorial wave measurements were available, an estimate of the H e + concentration could be inferred from the observed stop band between f cH e and H + -band wave power Qin et al. 2019). Here, in the absence of conjugate equatorial wave measurements, we adopt a statistical treatment, using typical H e + concentrations Ross et al. 2022). For typical peak power parameters Zhang et al. 2016b) f peak /f cp 0.4, hydrogen band EMIC waves, and 98% proton +2% helium ions (or f peak /f cp = 0.43 and 92% protons +8% helium ions which result in similar values), we get C 246. The corresponding theoretical estimate of the typical peak resonance energy, E peak is shown as a dashed red curve in Panel (d). Next, we used a factor C/2.8 instead of C, to get our theoretical estimate of the minimum resonance energy E R,min corresponding to significant wave-driven electron scattering toward the loss-cone. This is shown as a dashed red line in Panel (e) along with the observationally determined energies at half-max precipitatingto-trapped electron flux ratio E * min and their average. The factor C/2.8 corresponds to electron cyclotron resonance with hydrogen band EMIC waves of frequency f/f cp ∼ 0.56, well above the peak-power frequency f peak /f cp 0.4, and taken to be a reasonable approximation for f max /f cp in statistical averages of wave power. At that frequency, wave power is still finite albeit an order of magnitude lower than at peak-power frequency (Zhang et al. 2016c. This choice is consistent with the smaller precipitating-to-trapped flux ratio at E R,min (compared to that at E peak ), despite the theoretically anticipated increase in diffusion rate with decreasing energy for otherwise constant wave power Ni et al. 2015a;Summers and Thorne 2003). As we see from the panels under discussion, there is a reasonably good match between the theoretical expectation for E peak and E R,min independently derived from EMIC wave-power statistics alone (the dashed red lines), and the averages of E * and E * min from observations of the precipitating-to-trapped flux ratio at ELFIN (the solid black lines).
In addition, we show in Panel (e) two lower limits of E R,min estimated based on wave cyclotron damping by cold/cool ions near the proton gyrofrequency (see Sect. 2.2) at kc/ pi ∼ 1 (dashed blue line) and at kc/ pi ∼ 2 (dotted blue line). These two estimates fit well the average E * min values (solid black line) and the lower limit of E * min data points at L = 4 − 10, respectively. This suggests a significant effect of wave cyclotron damping in determining the wave spectral shape at high frequencies f > 0.6 f cp .

Dependence on AE, MLT and Concurrence with Whistler-Mode Waves
Let us now investigate some properties of the highly relativistic, strong electron precipitation EMIC wave-driven events, those with R = j prec /j trap > 1/2 at E * > 1 MeV, of the predominant spatial category of events, the ones with L < 7. These represent the most efficient Fig. 17 Properties of the most efficient EMIC wave-driven precipitation events, those exhibiting highly relativistic (E * >1 MeV), strong (R = jprec/jtrap > 1/2) electron precipitation, observed at L < 7. (a) Fraction of these most efficient events in each AE * bin divided by total number of most efficient events (AE * is the maximum AE in the preceeding 3 hours). (b) Same as (a) but as a function of MLT. (c) Fraction of most efficient events that also have 0 R < 1/3) in the three low-energy range categories listed in the abscissa. The rest (those with R > 1/3), representing concurrent low-energy moderate or strong precipitation, can still be a significant fraction of the most efficient EMIC wave-driven precipitation events (15 -35%, depending on energy range category). (d) Fraction of most efficient EMIC wave-driven precipitation events satisfying two different criteria for R = j prec /j trap at lower energy as defined in annotations. Implications are discussed in the main text subset of the previously discussed strong precipitation (R>0.5) EMIC wave-driven events, those with peak energy >1 MeV at L-shell L < 7, but otherwise very similar. Figure 17(a) shows that when binned as a function of AE * (the maximum AE in the preceding 3 hours), the fraction of these most efficient events over the total number of events within our database strongly increases from 250 nT to ∼ 1500 nT and is quite significant for AE * > 600 nT. It is noteworthy that the increase of this fraction with AE * happens despite the rapidly decreasing probability of intervals with large AE * , for values above ∼ 200 nT. Thus the smaller number of events in the highest AE * bin (1750 − 2500 nT) is likely due to the decrease in occurrence rate of such extremely high AE * periods. We interpret the result of Fig. 17(a) as a consequence of the fact that EMIC wave power (which determines the efficiency of electron scattering) is itself strongly dependent on geomagnetic activity, owing to activitydependent injections of anisotropic hot ions which provide the free energy for EMIC wave generation (Chen et al. 2010(Chen et al. , 2011. Figure 17(b) shows the fraction of the most efficient EMIC wave-driven electron precipitation events (normalized to all events within L < 7) as a function of MLT. It also reflects the relative occurrence probability of such events with MLT, since ELFIN's science zone collections are uniformly distributed in local time. This is much higher in the 18-24 MLT sector than elsewhere. We suggest this is because it is towards this sector that highly anisotropic hot ions drift, after being produced by nightside injections (which are known to peak in occurrence rate around the pre-midnight sector). The probability of efficient precipitation is much weaker at 0-16 MLT, likely because of the reduced anisotropy of the aforementioned hot ion populations (since these ions have had a chance to be modified by EMIC wave scattering at dusk). Thus when the injected ions reach the 12-16 MLT range they can generate only weaker EMIC waves, if any. Conversely, the 0-6 MLT sector can still harbor direct injections from the magnetotail, albeit at a lower probability than at pre-midnight, which could explain the reduced probability of both EMIC wave generation and subsequent efficient electron precipitation in that sector. Figure 17(c) shows the fraction of the aforementioned most efficient EMIC wave-driven events that have a low-energy j prec /j trap ratio that is clearly low, 0 R < 1/3. This condition is applied in three low-energy ranges: E 100 keV, E∈ [100,200] keV, and E∈ [100,300] keV. The figure reveals that 85% of events with strong precipitation ratio at E * > 1 MeV are accompanied by a significantly weaker precipitation ratio, < 1/3, at E 100 keV. However, the unity-complement of the third category (E∈ [100,300] keV) also shows that ∼35% of the most efficient EMIC wave-driven events still have j prec /j trap > 1/3 at some energy between 100 and 300 keV. Therefore, moderate to strong precipitation is still present at low energies (up to a few hundred keV), for quite a significant fraction of these most efficient, highly relativistic precipitation events.
So, we next investigate whether the precipitation at low (100 − 300 keV) energies we observed here during the most efficient EMIC wave-driven events could be due to electron scattering by whistler-mode chorus or hiss waves that may occur simultaneously with EMIC waves. Since the whistler-mode pitch-angle diffusion rate (and, thus, the precipitation efficiency) decreases with energy, E, the precipitating-to-trapped flux ratio R = j prec /j trap should first decrease and reach a minimum somewhere above 100 keV, before increasing again at relativistic energies due to cyclotron resonance with EMIC waves (e.g., see (Mourenas et al. 2021). We note that such intense sub-MeV precipitation has been reported previously (Hendry et al. 2017(Hendry et al. , 2019Capannolo et al. 2021), but in the past it was not possible to separate the effects of whistler-mode wave and EMIC wave scattering by studying at sufficient energy resolution and extent this R − E relationship. Here, we explore this relationship next, in Fig. 17(d). The figure shows, on the right, the fraction of all events with j prec /j trap > 1/2 at E * > 1 MeV for which R(200 keV)< 0.7 × R(100 keV) OR R(300 keV)< 0.7 × R(100 keV); in other words, where the spectral slope of the ratio R is decidedly negative (0.7 being an arbitrary but reasonable choice). Such a negative slope with energy at the energy range E < 300 keV likely corresponds to chorus or hiss wave-driven precipitation simultaneous with EMIC wave precipitation at MeV energies. We see that the depicted fraction is 22%, indicating that about one fifth of the most efficient EMIC wave events exhibit signatures of simultaneous whistler-mode chorus precipitation below 300 keV.
The OR condition has been used above to ensure inclusive accounting of a broad range of low energies and an increased statistical significance of the results. However, it is still possible that sub-MeV precipitation by EMIC waves due to the aforementioned nonresonant (Chen et al. 2016;An et al. 2022), or bounce resonant (Cao et al. 2017b;Blum et al. 2019) interactions, might extend down to the 300 keV range of energies. This could provide a positive R − E slope in the 300 keV range but retain a negative slope in the 200 keV range. Therefore, we examine a more stringent criterion, applying the AND condition to the above quantities. In other words, to examine bona fide concurrent whistler-mode and EMIC wave scattering, we investigate the fraction of putative EMIC wave-precipitation events (those with j prec /j trap > 1/2 at E * > 1 MeV) which exhibit both R(200 keV)< 0.7 × R(100 keV) AND R(300 keV)< 0.7 × R(100 keV). Figure 17(d), left, depicts this fraction. It shows that whistler-mode wave-driven precipitation may be present at E ≤ 300 keV for only ∼ 6.5%, or one sixteenth of all events with j prece /j trap > 1/2 at E * > 1 MeV at L < 7. These results suggest that at least ∼ 78% and possibly up to ∼ 93.5% of all intense EMIC wave-driven precipitation events extending down to 200 keV are likely not due to concurrent whistlermode and EMIC wave scattering, and deserve further scrutiny.

Consistency of Precipitation with Diffusion Theory
The statistically significant and moderately efficient precipitation observed at ∼200 keV in the presence of strong precipitation at highly relativistic energies (E * > 1 MeV) by EMIC waves, does not seem, at first glance, to be similarly due to resonant scattering by EMIC waves: Based on Van Allen Probes statistics , EMIC waves with sufficiently high frequencies (relatively close to the ion, and in particular the proton gyrofrequency) to provide resonant scattering and precipitation at such low energies (Kennel and Petschek 1966;Mourenas et al. 2022) indeed have very low wave power (at least 10 times lower than that at peak wave power (Zhang et al. 2016c). In addition, hot plasma corrections to the wave dispersion relation become necessary in the close vicinity of the ion gyrofrequency, where they can suppress the expected resonant scattering Cao et al. 2017a). Nonresonant scattering by EMIC waves with very sharp edges may be implicated (Chen et al. 2016;An et al. 2022), but the actual presence of such sharp wave-packet edges still has to be verified experimentally and its effects have yet to be compared with theory and simulations. Bounce resonance scattering can, in principle, provide precipitation at sub-MeV energies (Cao et al. 2017b;Blum et al. 2019), but the associated scattering rate is quite small for quasi-parallel EMIC waves and, thus, further statistical investigation of the most effective oblique EMIC waves is needed to evaluate the relative contribution of this mechanism. Nevertheless, some useful insights into the possible origin of the moderate sub-MeV electron precipitation that accompanies strong relativistic electron precipitation by EMIC waves, can be gained from a careful examination of the energy spectrum characteristics of EMIC wave-driven precipitation in our database. We do this next. Figure 18(a) shows the average trapped electron flux (in black) and the average precipitating-to-trapped flux ratio j prec /j trap (in red) as a function of energy for the dominant category of all highly relativistic electron strong precipitation events (at L < 7 with j prec /j trap > 1/2 at E * > 1 MeV). The flux ratio, denoting precipitation efficiency, increases approximately as γ 2 (where γ = 1 + E/mc 2 is the Lorentz factor) as kinetic energy E increases (blue dotted line). It starts from a low, approximately constant average value j prec /j trap ∼ 1/8 at E < 200 keV and rises to j prec /j trap > 1/2 at 1 MeV. It attains a high, approximately constant average value j prec /j trap ≈ 0.85 for E > 1.5 MeV. The low values of j prec /j trap ∼ 1/8 at the low energies could be partly due to the simultaneous presence of chorus wave-driven precipitation during a small fraction ≈ 6.5% of these events (as suggested by results in Fig. 17(d)). Since such chorus events should give a higher j prec /j trap ratio at E < 100 keV , whereas all the other events without chorus waves should give a slightly smaller j prec /j trap in the same energy range, the resulting average j prec /j trap can be nearly constant at 50 − 200 keV. The flattening of the slope of j prec /j trap versus energy at E ≥ 1.5 MeV is consistent with the most common minimum cyclotron resonance energy with intense EMIC waves being E R,min > 1 MeV based on wave statistics from satellites (Summers and Thorne 2003;Kersten et al. 2014;Cao et al. 2017a;Ross et al. 2021). It is therefore likely that, with some energies below and some above E R,min in the range ∼ 1.5 MeV to 4 MeV in our dataset, the averaging of j prec /j trap over all events lowers its value (to well below the peak ratio value of ≈ 1) and flattens its spectrum. Fig. 18 Spectral properties of EMIC wave-driven electron precipitation. (a) Average spectrum of trapped fluxes (black curve) and of precipitating-to-trapped flux ratio (red curve) for events with j prec /j trap > 1/2 at E * > 1 MeV. The dotted blue curve shows a best least-squares fit to the average flux ratio at low energy. (b) Average precipitating-to-trapped flux ratio for the events of Panel (a), plotted as a function of the normalized energy E/E * , where E * is the peak precipitating-to-trapped flux ratio energy for each event (a proxy for the minimum resonance energy E R,min with the most intense waves). Best least-squares fits to the average flux ratio are shown below the peak, j prec /j trap ≈ γ 2 (E)/γ 2 (E * ) ∼ 0.065 γ 2 (E) for E * = E * ∼ 1.45 MeV (dashed blue curve), and above the peak, j prec /j trap ≈ 1.01 · (E/E * ) −1.01 (dashed red curve). Note that if the last point (E/E * ≈ 4.84) with a small number of measurements was discarded, the best fit would become j prec /j trap ≈ 0.92 · (E/E * ) −0.8 (dotted red curve) Normalizing the energy to that of our estimate of maximum precipitation, E*, should help better organize the spectra and eliminate the aforementioned apparent flattening. This is done in Fig. 18(b), which shows the average ratio j prec /j trap as a function of E/E * . For each event, the energy E has been normalized to the event's E * at maximum j prec /j trap , the approximate energy for cyclotron resonance with the most intense EMIC waves. A leastsquares fit to the data at E ≥ E * (red dashed line) shows that j prec /j trap ∼ (E/E * ) −1 . This agrees well with the prediction of quasi-linear diffusion theory, which can be expressed as where D αα is the electron bounce-averaged quasi-linear pitch-angle diffusion rate at the loss-cone angle α LC and, again, γ is the Lorentz factor (Kennel and Petschek 1966;Ni et al. 2015a;Mourenas et al. 2016Mourenas et al. , 2021Mourenas et al. , 2022. This suggests that electron precipitation driven by EMIC waves is described well by a quasi-linear diffusion treatment, which was developed under the assumption of low amplitude waves exhibiting low coherence (random phases) leading to slow phase space evolution compared to the wave period.
The precipitation ratio falloff with normalized energy at E/E * > 1 in Fig. 18(b) is consistent with our assertion that the energy E * likely corresponds to the energy for cyclotron resonance at the frequency of the most intense EMIC waves. This is because, we contend, at lower energies, E/E * < 1, the observed rapid decrease of j prec /j trap with progressively decreasing energy, down to E/E * ∼ 0.1, likely corresponds to cyclotron resonance with less intense EMIC waves at progressively higher frequencies (Zhang et al. 2016c;Ross et al. 2021), or to electron nonresonant scattering by the steep edge of some EMIC wave packets (Chen et al. 2016;An et al. 2022), both of which result in a reduced efficiency of scattering compared to that at peak wave power. Thus, viewed as a function of decreasing energy starting from the highest energies, while that ratio increases exponentially at first (above E * ) due to the increase in the diffusion rate, it starts to decrease just as quickly later (below E * ) due to the fast wave power decrease -the turning point in that behavior being consistent with E * , the resonance energy at peak wave power.
To check whether the decrease of j prec /j trap with decreasing energy below E * is actually consistent with EMIC wave power observations, it is useful to fit that ratio below E * (blue dashed line in Fig. 18(b)). This fit quantifies the ratio's observed energy dependence at E < E * as j prec /j trap ∼ γ 2 (E)/γ 2 (E * ) ∼ 0.065 γ 2 (E), where we have replaced E * with its mean value E * ≈ 1.45 MeV in our dataset (E * varies between ∼ 1 MeV and ∼ 3 MeV). With this replacement, the best fit in Fig. 18(b) is identical to the best fit in Fig. 18(a). The fit value for j prec /j trap ∼ 1 actually occurs at an energy E ∼ E * . For a peak wave power nearly flat over a significant range of frequencies, E * represents the lowest energy at which cyclotron resonance with the highest power waves is achieved, because the diffusion rate D αα (E) that characterizes the behavior on the high-energy side decreases toward higher energy even for a flat wave power spectrum B 2 w (f ) . In that case, E * corresponds to the highest frequency over the nearly flat region of peak wave power, henceforth simply denoted f peak , for which the diffusion rate D αα (E R,min

Consistency of Observed Precipitation with Equatorial EMIC Wave Power
Next, using comparisons with observations, we investigate the simplest explanation for the j prec /j trap dependence on E/E * , at E < E * : cyclotron resonance with progressively less intense EMIC waves at higher frequencies. This is done in two steps: First, assuming that this precipitation is due to quasi-parallel hydrogen band EMIC waves we employ quasilinear theory to infer from the fit to the ELFIN measurements, described earlier, the statistical EMIC wave-power ratio B 2 w (f )/B 2 w (f peak ). And second, we compare this inferred wave power ratio with published statistics of EMIC wave power B 2 w (f ) directly observed by the Van Allen Probes near the equator in 2012-2016 (Zhang et al. 2016c). As a reminder, nonlinear resonant scattering by an ensemble of independent, short-range, and large amplitude wave packets that may partake in the measured average wave spectral shape should produce average j prec /j trap energy spectra similar to that caused by classical diffusion by waves with the average wave power. In other words, the diffusive formalism, even though borrowed from quasi-linear theory, also applies for such nonlinear interactions if they participate statistically in the process. Yet, since we focus now on the precipitation of sub-MeV electrons, at E < E * , which can reach cyclotron resonance only with high frequency waves, at f > f peak , of much lower amplitudes than at peak wave power, the contribution from nonlinear interactions is likely much smaller than in the case of multi-MeV electron precipitation.
Towards the first step, we start from the full quasi-linear expressions for the precipitating to trapped flux ratio at ELFIN j prec (α)/j trap (α trap ) ( (Kennel and Petschek 1966) and ). Pitch angles are referenced to the equator and α trap is such that ln(sin α trap / sin α LC ) ∼ 1/20 (where α LC is the equatorial loss-cone angle). Averaging j prec over equatorial pitch-angles α < α LC , we get j prec j trap where z 0 = 2α LC / √ D αα τ B , and τ B ∼ γ /(γ 2 − 1) 1/2 is the bounce period (Schulz and Lanzerotti 1974). This integral can be approximated as: j prec /j trap = 0.9/z 0 with less than 20% error for z 0 ∈ [0.9, 8], corresponding to j prec /j trap ∈ [0.11, 1]. At low energies, E < E * , and for quasi-parallel left-hand-polarized hydrogen band EMIC waves ) with a monotonically decreasing power B 2 w (f ) toward higher frequencies f > f peak in Van Allen Probe statistics (Zhang et al. 2016c), the most efficient wavedriven pitch-angle diffusion of electrons near the loss-cone through cyclotron resonance should occur for E(f ) ∼ E R,min (f ) at similar latitudes λ R close to the equator for all E in that range . Bounce-averaging ) the full expression of the local pitch-angle diffusion rate in the cold plasma approximation (Summers and Thorne 2003;Su et al. 2012) then gives at an equatorial pitch-angle α α LC .
Combining the theoretical scaling laws for j prec /j trap and D αα and the fit to ELFIN observations j prec /j trap ∼ γ 2 in Fig. 18, gives us the estimate of the wave-power spectrum consistent with these observations: To obtain the mapping between energy and frequency, we utilize the cyclotron resonance condition coupled to the cold plasma dispersion relation. Its expression for E < E * and f > f peak is combined with its expression for E * and f peak , the (highest-)frequency of peak EMIC wave-power (f peak /f cp ∼ 0.37 − 0.41, as seen in statistical wave observations from the Van Allen Probes (Zhang et al. 2016c)). We also assume an ion composition with >94% protons, as appropriate for when hydrogen band waves are present Ross et al. 2022). This yields a second order equation for f/f cp , with solution f/f cp 2/(1 + √ 1 + 4C), where C = (f cp /f peak )(f cp /f peak − 1)E(1 + E)/ (E * (1 + E * )). The resonant frequency f (E)/f cp , expressed as a function of the resonance energy in the above equation, allows us, by substitution in the previous equation, to obtain the EMIC wave-power ratio B 2 w (f )/B 2 w (f peak ) = B 2 w (E)/B 2 w (E * ) inferred from ELFIN, with an error smaller than 40% for j prec /j trap ∈ [0.11, 1]. Note that the above normalizations of j prec /j trap to its level at E * and of EMIC wave power to its level at f peak /f cp allow us to eliminate the dependencies of j prec /j trap and D αα on f pe /f ce and on the non-normalized EMIC wave power at the equator conjugate to ELFIN, which are unknown. Here we only had to assume that the ensemble of j prec /j trap measurements from ELFIN are statistically compliant to an ensemble of EMIC wave power observations on another, equatorial, platform in order to convert the energy spectrum of the precipitation to a frequency spectrum of wave power that should be consistent with that precipitation.
We proceed, now, to the second step in our consistency check between the ELFIN observations of electron precipitation and independently collected EMIC wave-power spectra, in order to assess the validity of quasi-linear diffusion in describing the precipitation down to low energies. Figure 19(a) compares the wave power ratio B 2 w (f )/B 2 w (f peak ) inferred above from ELFIN statistics of j prec /j trap (solid blue curve) with the measured hydrogen band EMIC wave power ratio from Van Allen Probes statistics when f pe /f ce > 15 (Zhang et al. 2016c) in different MLT sectors (black, green, magenta, and red curves). In the high density plasmaspheric plume or at the plasmapause, as implied by the f pe /f ce ratio criterion used in this subset of wave data, the cyclotron resonance condition can indeed be satisfied for the optimum parameters E * ∼ 1.45 MeV in Figs. 18(a,b) and f peak /f cp ∼ 0.37 from the above wave statistics. Figure 19(c) indicates that 200 keV electrons are then in resonance with waves at f/f cp ∼ 0.8 near the equator. We see in Fig. 19(a) that the observed Fig. 19 (a) EMIC wave power ratio B 2 w (f )/B 2 w (f peak ) inferred, using quasi-linear theory, from ELFIN statistics of precipitating-to-trapped electron flux ratio j prec /j trap (solid blue curve), and compared with statistical EMIC wave power ratios obtained from Van Allen Probes 2012-2016 observations in four different MLT sectors when f pe /f ce > 15 (black, green, magenta, and red curves). Specifically, the solid blue curve was derived from the best fit to j prec /j trap in Fig. 18(a)  statistical EMIC wave power ratios in the 12-22 MLT sector (magenta and red curves) agree quite well with the power ratios inferred from ELFIN data (blue curve with dashed lines indicating uncertainties) over two decades in power, in the frequency range f/f cp ∼ 0.37 up to f/f cp ∼ 0.8 − 0.95. Therefore, in high density regions (those with f pe /f ce > 15), the observed EMIC wave power ratio is consistent with that inferred from ELFIN precipitation of j prec /j trap from ∼ 1.45 MeV down to 200 keV, sometimes even down to 100 keV.
In regions of lower plasma density, with 5 < f pe /f ce < 15, EMIC wave statistics from the Van Allen Probes show that f peak /f cp ∼ 0.41 (Zhang et al. 2016c), corresponding to a higher E * ∼ 2.5 MeV for f pe /f ce ∼ 13 − 15. For these parameters, Fig. 19(b) shows that in the 12-22 MLT sector the observed EMIC wave power ratios in the range f/f cp ∼ 0.41 to f/f cp ∼ 0.90 − 0.95 (Zhang et al. 2016c) at Van Allen Probes (magenta and red curves) are consistent with those inferred from ELFIN statistics (blue curves). In the above f pe /f ce range, the frequencies where agreement prevails correspond to j prec /j trap from ∼ 2.5 MeV down to ∼ 100 − 200 keV.
It is evident from Figs. 19(a,b) that the inferred and observed EMIC wave power ratios agree well, and up to high frequencies (low resonance energies), only in the 12-22 MLT sector, i.e., near dusk. This is quite consistent with the fact that most of the EMIC wavedriven events at L < 7 in our database occurred near dusk ( Fig. 17(b)). It is also consistent with the predominance of similar events at dusk in Firebird-II statistics (Capannolo et al. 2021). In the 0-3 MLT sector (black lines in Figs. 19(a,b)) the observed EMIC wave power is consistent with EMIC wave-driven electron precipitation down to ∼ 300 keV only when f pe /f ce > 15. Few events exist in this sector in our database. In the 4-12 MLT sector (green lines) the observed EMIC wave power is too weak to drive significant electron precipitation below ∼ 1 MeV, consistent with the absence of precipitation events in that sector in our database. (Note: Since EMIC wave power spectra depend weakly on L in Van Allen Probes statistics (Zhang et al. 2016c;Ross et al. 2021), the present results should hold for L ∈ [4,7].) Therefore, we see that the observed, moderately efficient electron precipitation at energies as low as ∼ 200 − 300 keV could simply be due to quasi-linear scattering by moderate intensity EMIC waves at high frequencies, up to f/f cp ∼ 0.8 − 0.9, provided that such waves are present during a majority of these events with a similar average intensity as in statistical empirical wave models. Such waves may still be sufficiently below the proton gyrofrequency to evade hot plasma effects, at least to first order (Chen et al. 2011;Ross et al. 2021). Assuming a maximum wavenumber that can be attained in the presence of hot plasma effects given by kc/ pi ∼ 2 (corresponding to a typical cold ion temperature of ∼ 10 eV in a plasmaspheric plume or just within the plasmasphere) for left-hand-polarized hydrogen band waves in a plasma with more than 94% protons, the minimum electron energy for cyclotron resonance can indeed be as low as ∼ 150 keV to 350 keV for f pe /f ce ∼ 15 to 30, when a sufficient transient hot H + temperature anisotropy A > 2.3 generates such high-frequency waves (Chen et al. 2011. A typical precipitation event of this kind has been analyzed in Sect. 4.2. Another important point is that the measured ratios j prec /j trap ∼ 1/8 − 1 from low energy, E E * , to E * actually correspond to a quasi-linear diffusion close to the strong diffusion regime with characteristic time scales for partial or full loss-cone filling on the order of a quarter of a bounce period τ B /4 ∼ 0.25 s (Kennel and Petschek 1966). Such fast diffusive time scales are therefore still consistent with the time scales of the observed electron precipitation at ELFIN. Thus, even in the quasi-linear regime it is possible for temporal or spatial wave power variations encountered by drifting and bouncing electrons to explain the short-lived (sub-spin), bursty nature of the precipitation at ELFIN.
The present results agree with test particle hybrid simulations suggesting that lowamplitude, high-frequency hydrogen band EMIC waves could be the main cause of 0.55 MeV electron precipitation (Denton et al. 2019). Several previous studies have also provided hints of a prevalence of hydrogen band waves in driving sub-MeV relativistic electron precipitation Qin et al. 2018;. Even when such high-frequency waves are not evident at apparently conjugate measurements at the equator, this may be due to mapping uncertainties -the waves could still be present 0.5 Earth radii away where the ionospheric electron precipitation measurements actually map to. Further analysis of EMIC wave data and sub-MeV precipitation with ELFIN's twin spacecraft, providing higher spatio-temporal resolution of these phenomena, will be able to further study these points in greater detail. In addition, one cannot totally rule out a possible role of helium band EMIC waves in driving some sub-MeV electron precipitation, in spite of the expected strong cyclotron damping very close to the helium gyrofrequency Cao et al. 2017a). A similar study as in Fig. 19, but focusing on helium band waves, could be performed to check this point. This would require extensive numerical calculations including hot plasma effects in the dispersion relation .
As a final note, we comment on the global contribution of EMIC waves to sub-MeV electron scattering. Although bursts of EMIC waves can drive relatively intense sub-MeV electron precipitation in a narrow MLT sector, their very weak time-and MLT-averaged power at high frequencies should prevent them from contributing significantly to the global loss rates of sub-MeV electrons up to high equatorial pitch-angles α ∼ 90 • , above and beyond what can already be done by typical chorus waves in the dawnside trough Boynton et al. 2017;Agapitov et al. 2018;Drozdov et al. 2019Drozdov et al. , 2020Miyoshi et al. 2020;Ross et al. 2021;Zhang et al. 2022). This is because at a fixed energy E ∼ E * , cyclotron resonance between electrons of equatorial pitch-angle α and hydrogen band EMIC waves at a frequency f indeed corresponds to a scaling cos 2 α/ cos 2 α LC ∼ (f cp /f − 1)/(f cp /f peak − 1)f peak /f . This implies that electrons of higher α reach resonance with waves of higher frequency f > f peak (e.g. see first standalone equation in Mourenas et al. (2016)). Therefore, the steep decrease of EMIC wave power from f/f cp ∼ 0.4 to ∼ 0.9, by two orders of magnitude in Fig. 19, implies a rapid decrease of the quasi-linear pitch angle diffusion rate D αα (α) away from the loss-cone (Summers and Thorne 2003;Mourenas et al. 2016). This is consistent with pitch angle bite-out signatures (stronger electron loss near the loss-cone than at higher pitch angles) produced by the most intense EMIC waves (Usanova et al. 2014;Adair et al. 2022) at energies of ∼ 1 − 4 MeV. Such signatures demonstrate the inefficacy of EMIC waves in depleting the relativistic electron flux at equatorial pitch angles far from the loss cone, absent intense chorus waves Zhang et al. 2017). The same arguments hold at sub-MeV energies, except that cyclotron resonant EMIC waves have higher frequencies than at ∼ 1 − 4 MeV and, therefore, much lower amplitudes, making them even less efficient than typical intense chorus waves in the dawnside trough in driving the precipitation of electrons near the loss-cone.

Summary and Discussion
Relativistic electrons in the near-Earth environment are an important contributor to space weather and may play a significant role in charged particle energy input to the atmosphere from space. The question of how such electrons, after being accelerated in near-Earth space by waves or transported into it by, say, radial diffusion, are lost remains open. This question is particularly vexing for the outer radiation belt, at (L<7), which has the highest fluxes of such electrons and is thus most important for space weather. Aside from magnetopause shadowing (the result of magnetospheric compressions changing electron drift paths from trapped to open), and field-line curvature scattering precipitation at the isotropy boundary, both affecting L-shells beyond the outer radiation belt, and whistler-mode wave-driven precipitation which mostly affects sub-MeV electrons, most attention on relativistic electron loss has been placed on these electrons' interaction with EMIC waves. This is because EMIC waves can resonate with electrons of highly relativistic energies > 1 MeV under high-density, low magnetic field plasma conditions that can be realized in the outer edge of the plasmasphere, well inside the outer radiation belt at active times. As they can act on trapped, outer radiation belt electrons, such waves can change significantly the outer radiation belt fluxes and therefore are an important process to include in space weather models. The recent launch of the ELFIN mission, with a goal to determine whether EMIC waves are predominantly responsible for this loss or whether other wave processes are implicated, provided a unique dataset of 50 -5000 keV electrons obtained on a polar, low-altitude orbit, with which we can address this question comprehensively and for the first time. In this paper, we primarily focused on the question of whether EMIC waves can be definitively shown to be responsible for the observed scattering of relativistic electrons, and whether whistler-mode chorus waves (the other candidate for pitch-angle scattering of sub-MeV and up to MeV electrons into the loss cone) simultaneously might be implicated.
After a short review of EMIC wave generation and its interaction with energetic electrons in Sect. 2, we presented in Sect. 3 the first comprehensive examples of typical EMIC wave-driven precipitation of such electrons from ELFIN. The high energy and pitch-angle resolution of ELFIN allowed us to identify and quantitatively investigate the energy range of such precipitation and study its properties using the precipitating-to-trapped flux ratio, R, measured by a single detector on a spinning platform, thus avoiding problems with multiple detector inter-calibration. This ratio increasing with energy and peaking in the range >0.5 MeV is a tell-tale signature of EMIC wave-driven precipitation. In Sect. 3.1 we showed how this spectral signature is differentiated quite well from that of whistler-mode chorus driven precipitation which has a decreasing ratio as a function of energy from as low as 50 keV (the lower energy limit of the ELFIN detector) to >0.5 MeV. (Both spectral types can be observed on occasion simultaneously, as expected, since whistler-mode chorus and EMIC waves can co-exist -in that case, the two spectral shapes can still be well separated.) Next, in Sect. 3.2 we used a case study of typical EMIC wave-driven emissions, accompanied by equatorial measurements of EMIC waves on THEMIS to demonstrate that the bursty nature of the precipitation at ELFIN (lasting only 7-14 s, or 2.5-5 ELFIN spin periods) is consistent with the spatial extent of the EMIC wave region at the equator (lasting 6min). In this case, the equatorial extent of the EMIC region was L ≈0.23 which is consistent with published statistical averages of L ≈0.5 (Blum et al. 2016. Thus, the burstiness of EMIC wave-driven precipitation at polar orbiting, low-altitude spacecraft like ELFIN is likely typical. We have also shown that the burstiness can extend to sub-spin timescales, 1-5 spin sectors, at least in part due to the spatial variability of the wave-power at the equator. This can cause the precipitating-to-trapped flux ratio to be, on occasion, aliased and exceed unity. While there are reported cases when this flux ratio on ELFIN exceeds unity due to nonlinear scattering, an assertion supported by fortuitous measurements of the EMIC wave amplitudes at conjugate observatories , conclusions on the preponderance and significance of nonlinear effects cannot be drawn from the ratio alone without further statistical or detailed case study analyses, both left for the future. However, we have argued that nonlinear effects (from presumably mainly short wave packets -e.g., see various examples of such short packets in  and ) are statistically incorporated into the spectral shapes of precipitation at ELFIN, just as these high-amplitude EMIC wave packets are incorporated into the statistical averages of equatorial wave spectra, so that a direct comparison using a diffusive formalism should be possible.
In Sect. 3.3 we showed that consecutive ELFIN passes over the same MLT region can reveal the EMIC wave-driven precipitation's spatio-temporal evolution in location, extent (in L and L), and intensity. This allows us to study the EMIC wave-generation region in a geophysical context (e.g., during substorms and storms) as well as the potential of the waves for reducing the outer radiation belt flux.
To evaluate the consistency of the observed energy of moderate and strong precipitation at ELFIN (ratios R ∈[0.3,0.5] and R >0.5, respectively) with the resonance energy expected for EMIC wave-driven precipitation, we utilized in Sect. 4 conjunctions with ground based or equatorial assets. Using such data-informed estimates of f pe /f ce and of f max /f ci , with only occasional support from models, we examined the resonance energy that would be consistent with these estimates. We found it to be in agreement with the observed energies of moderate and strong precipitation seen at ELFIN. This demonstrates that the observed precipitation as identified by the ratio R is indeed consistent with theoretical expectation from resonant interactions.
We proceeded, in Sect. 5, to study statistically an ensemble of EMIC wave-driven events observed by ELFIN. These were identified in two years of ELFIN data (2019-2021) as enhancements in the precipitating-to-trapped flux ratio R, peaking at > 0.5MeV. The average event duration on ELFIN, 17 s (6 spins), or L ∼ 0.56, is consistent with published reports of the typical radial extent of EMIC waves in the magnetosphere, L ∼ 0.5 (Blum et al. 2016, validating our assertion that the bursty nature of the precipitation is due to the spatial localization of the EMIC wave interaction region in the magnetosphere. The most populous category of events is those occurring at pre-midnight and dusk, at L ∼5-7. A second class of events was found at midnight, post-midnight and dawn at L ∼8-12. These two categories are roughly co-located with the two main populations of EMIC waves in the magnetosphere (Min et al. 2012), further solidifying our assertion that our observations of relativistic electron precipitation are caused by EMIC waves.
The peak precipitating-to-trapped flux ratio exceeds one for many of the events studied. Considering that such a ratio, if scrutinized and confirmed, could signify nonlinear interactions, the only means of exceeding the quasi-linear strong diffusion limit (e.g., Kubota et al. (2015), Grach and Demekhov (2020)), we examined further how it may arise. We found that on many occasions extreme burstiness of the precipitation, lasting 1-5 spin phase sectors was evidenced in the data. (Each ELFIN spin period has 16 sectors and the loss cone most often contains 6 such sectors.) Examination of EMIC waves at THEMIS E on one event studied revealed that this burstiness is, in fact, to be expected at ELFIN due to the spatial localization of the EMIC wave coherence and amplitude arising from interference of wave packets in the presence of extreme density gradients. This burstiness in precipitation can result in temporal aliasing of ELFIN's precipitating-to-trapped flux ratio, R, resulting in values that can exceed 1 or be well below 1. Careful analysis of individual events together with modeling of measured wave spectra from ancillary datasets can, at times, confirm the nonlinear nature of the scattering (Grach et al. 2021). Either future extensive case-by-case analysis or careful statistical analysis to identify and remove events subject to aliasing is required to establish the preponderance of nonlinear effects. However, we argued that by analogy with similar effects for whistler-mode chorus waves, it is still quite likely that precipitation due to nonlinear EMIC wave scattering can conform to a statistical analysis using a diffusive treatment, because precipitation bursts contribute to the average flux (and flux ratio) at a level commensurate to the contribution of nonlinear wave power bursts in statistical averages of equatorial wave power. We proceeded with such a statistical analysis even though we recognize that in small or large part the precipitation and the associated wave amplitudes considered may incorporate nonlinear effects.
Using this database we have shown (Sect. 5.1) that: • The typical energy of the peak precipitating-to-trapped flux ratio, E * , a measured proxy for the resonance energy at the frequency of peak wave power, E peak , is in the range ∼ 1 − 3 MeV. This is consistent with expectation for EMIC wave-electron resonant interactions, based on the cold plasma dispersion relation and prior statistics of wave-power near peak wave-power, f peak , at the equator. • The above measured energy E * decreases with L-shell as ∼ L −0.5±0.1 . This dependence is also in good agreement with theoretical expectations. • The typical energy of the half-peak of the precipitating-to-trapped flux ratio, E * min , a proxy for the minimum resonance energy E R,min corresponding to significant wave-driven electron scattering toward the loss-cone, is in the range ∼ 0.5 − 1.5 MeV and falls off with Lshell at the same rate as E * . This too is consistent with theoretical expectations based on equatorial wave-power estimates at the maximum frequency below the ion gyrofrequency where wave-power remains significant. This shows that sub-MeV electron precipitation accompanying multi-MeV electron precipitation, is also likely driven by EMIC waves.
We next examined (Sect. 5.2) the properties of the most intense (R >1/2), highly relativistic (having E * > 1 MeV), EMIC wave-driven electron precipitation events, of the most populous category in our database (L<7). We consider these to be the most efficient EMIC wave-driven events in our database. Such events are highly correlated with geomagnetic activity. And, as expected from the statistics of all events, they are predominantly seen near dusk (MLT ∼ 18) with their occurrence rates dropping precipitously (though still finite) at noon and post-midnight. We found that: • About 35% of the time, the most efficient EMIC wave-driven events are still associated with moderate or strong (R >1/3) precipitation of electrons with energies as low as 100-300 keV. Only 6.5% of the time do such events have a definitively negative slope of R versus E, at that low energy range. For most events (at least 78% and potentially as many as 93.5%) that slope is positive. This suggests that simultaneous energetic electron scattering by whistler-mode chorus and EMIC waves (as was presented in Fig. 1, right column, and discussed in Sect. 3.1) are a minority, occurring for only one in every 16 of the highly relativistic, strong EMIC wave-driven events. • The average R versus E spectrum of the most efficient EMIC wave-driven events exhibits an exponential increase with energy (more precisely with the square of the Lorentz factor, R ∼ γ 2 ) up to its peak value R ∼ 1 at E ∼ E * (whose average value is E * ∼ 1.45 MeV). Examined versus energy normalized to the peak-precipitation energy, E/E * , R exhibits an exponential decay with normalized energy away from (on both sides of) unity: (R ∼ (E/E * ) −1 above and R ∼ (E/E * ) +1 below E/E * = 1). This decay is in agreement with quasi-linear diffusion theory for resonant scattering by a typical equatorial EMIC wave power spectrum that falls off away from its peak wave power, which is consistent with observations. • Based on both points above, the majority of the weaker precipitation below 1 MeV is thus probably due to electron resonant interaction with EMIC waves at frequencies f > f peak , i.e., above the peak power frequency of hydrogen band EMIC waves (as was also inferred by the E * min decrease as a function of L-shell in Sect. 5.1 and summarized above). Even though this wave power is much smaller than at f peak , it is still finite (Denton et al. 2019;Zhang et al. 2016c. But because of the significantly higher trapped electron flux at lower energies, even waves of such small power may still produce finite precipitating fluxes at ∼ 200 − 500 keV, as observed. • Very low energy (∼ 100 keV) electron precipitation at times of strong, highly relativistic precipitation, could still be predominantly due to interactions with very low-amplitude high-frequency hydrogen band EMIC waves, excited very near the equatorial gyrofrequency, f cp , by 10 − 100 eV anisotropic ion populations (Teng et al. 2019;. Such waves can more easily experience strong cyclotron damping by low energy ions, likely causing them to be transient and rendering their equatorial observations sparse. Since ∼ 100 keV electrons cannot reach cyclotron resonance with the most frequent hydrogen band EMIC waves (Zhang et al. 2016c), other precipitation mechanisms should also be examined and quantified. Those include nonresonant scattering by EMIC wave-packets with sharp edges (Chen et al. 2016;An et al. 2022), which can scatter electrons well below the minimum resonance energy while still exhibiting a rising R versus E spectrum. However, simultaneous whistler-mode wave-driven precipitation which is very efficient at such low energies, is clearly contributing to this precipitation (as discussed above), even though it is likely a minority, and should also be further considered Ma et al. 2016b).
Our results confirm the crucial role played by EMIC waves in relativistic electron losses in the Earth's outer radiation belt (at L-shells between 5 and 7, where ELFIN detected most EMIC wave-driven precipitation events). We find that the main predictions of the classical resonant scattering models of EMIC wave-driven precipitation (e.g., typical energy of precipitation peaks, the energy spectra and the L-shell dependence of minimum and peak-power resonance energy) are consistent with ELFIN case studies and statistical results presented herein. However, ELFIN's observed fine structure of the precipitation in energy, pitch-angle and time reveal interesting details pertaining to the nature of low-energy precipitation concurrent with strong, highly relativistic electron scattering, and the origin of precipitation exceeding the strong diffusion limit, that deserve further attention. Towards that end, significant new information on EMIC wave-driven precipitation can be obtained by further combining ELFIN electron measurements and theoretical modelling of wave-particle interactions with several ancillary datasets: First, numerous ELFIN conjunctions with ground-based and equatorial spacecraft such as Van Allen Probes, THEMIS, ERG and MMS exist; the equatorial spacecraft can provide direct information on the EMIC waves implicated in the scattering. They can also measure the near-equatorial density and magnetic field, and thus constrain the f pe /f ce ratio. The large dataset of ELFIN (spanning four years of operation, from 2018 to 2022) provides ample opportunities for studies using such conjugate observations. These also present good opportunities to investigate the relationship between the fine structure of EMIC waves at the equator (wave packet coherence (Blum et al. 2016, frequency ranges Zhang et al. 2016c), and hot plasma effects (Cao et al. 2017a;Chen et al. 2019)) and the energy distribution of precipitating electrons at ELFIN. Second, ELFIN's energetic particle detector for ions (EPDI), which measures total ions over the same energy range as electrons (50 -5000 keV), has been calibrated and has collected data for approximately three months prior to the end of the mission on both satellites. And finally, ELFIN's fluxgate magnetometer (FGM) instrument, which provides field-aligned current information and has the capability of measuring waves up to a Nyquist frequency of 5 Hz, has been operating well and has provided data on both satellites for several years (note that sub-spin resolution calibration to reveal EMIC waves is still ongoing). These datasets can provide information on energetic ion scattering and precipitation by the same EMIC waves that scatter energetic electrons in the inner magnetosphere (when the ion resonance energy is sufficiently high at appropriate f pe /f ce ratios), on the location of the EMIC wave excitation relative to large scale magnetospheric current system sources, and on EMIC waves potentially seen at ELFIN at the same time as electron precipitation. Critical open questions are (i) What is the overall contribution of nonlinear interactions and what can we learn about the physics of such interactions from ELFIN's precipitation energy and pitch-angle spectrum? (ii) What effects are responsible for precipitation at sub-MeV energies? (iii) What is the relative contribution of EMIC wave-driven precipitation to the total loss of energetic electrons in the inner magnetosphere and to the total atmospheric energy deposition by energetic particles? There is already good theoretical background for further model development of nonlinear resonant interactions (Albert and Bortnik 2009;Kubota et al. 2015;Kubota and Omura 2017;Grach and Demekhov 2020;Grach et al. 2021) and nonresonant interactions (Chen et al. 2016;An et al. 2022). Such development can proceed on solid grounds only if guided by statistically-derived properties of the observed electron precipitation and ancillary measurements afforded by the above datasets.