Loss Of Relativistic Electrons Research Articles

Inclusion of Nonresonant Effects Into Quasi‐Linear Diffusion Rates for Electron Scattering by Electromagnetic Ion Cyclotron Waves

AbstractElectromagnetic ion cyclotron (EMIC) waves are a key plasma mode affecting radiation belt dynamics. These waves are important for relativistic electron losses through scattering and precipitation into Earth's ionosphere. Although theoretical models of such resonant scattering predict a low‐energy cut‐off of ∼1 MeV for precipitating electrons, observations from low‐altitude spacecraft often show simultaneous relativistic and sub‐relativistic electron precipitation associated with EMIC waves. Recently, nonresonant electron scattering by EMIC waves has been proposed as a possible solution to the above discrepancy. We employ this model and a large database of EMIC waves to develop a universal treatment of electron interactions with EMIC waves, including nonresonant effects. We use the Green's function approach to generalize EMIC diffusion rates foregoing the need to modify existing codes or recompute empirical wave databases. Comparison with observations from the electron losses and fields investigation mission demonstrates the efficacy of the proposed method for explaining sub‐relativistic electron losses by EMIC waves.

BARREL Observations of Microburst Events With a Slowly‐Varying Component

AbstractElectron microburst precipitation has been shown to have significant potential for depletion of the outer radiation belt. We present observations from the Balloon Array for Radiation belt Relativistic Electron Losses (BARREL) of six (five unique and one dual‐balloon observation) microburst events, each containing minutes to hours of persistent microbursts. We find that each event included a slowly‐varying smooth precipitation component underlying the bursty component. The smooth component has not yet been fully characterized in the literature; we have written a program to identify microburst events and quantify the relative contributions of each component in the BARREL data. In all six events analyzed, the smooth component contributed more to the total X‐ray counts measured, indicating that the smooth component could contribute significantly to radiation belt loss.

Modeling the Global Distribution of Chorus Wave Induced Relativistic Microburst Spatial Characteristics

AbstractThe full spatiotemporal distribution of chorus wave‐induced relativistic electron microburst is modeled for chorus waves originated from different L shells and MLTs, based on the newly developed numerical precipitation model (Kang et al., 2022, https://doi.org/10.1029/2022gl100841). The wave‐particle interaction process that induces each microburst is analyzed in detail, and its relation to the chorus wave propagation effects is explained. The global distribution of maximum precipitation fluxes and scale sizes of relativistic microbursts is then obtained by modeling chorus waves at different L‐shells and local times. The characteristics of dawn and midnight sector microbursts have little difference, but the noon sector has much larger maximum flux and much smaller full width at half maximum, which may be due to dayside's low electron flux in the Landau resonance range. This suggests the controlling effect of keV electrons on the MeV electron precipitation intensity and properties and the overall relativistic electron loss in the outer radiation belt.

New Chorus Diffusion Coefficients for Radiation Belt Modeling

AbstractWhistler mode chorus is an important magnetospheric wave emission playing a major role in radiation belt dynamics, where it contributes to both the acceleration and loss of relativistic electrons. In this study we compute bounce and drift averaged chorus diffusion coefficients for 3.0 < L* < 6.0, using the TS04 external magnetic field model, taking into account co‐located near‐equatorial measurements of the wave intensity and fpe/fce, by combining the Van Allen probes measurements with data from a multi‐satellite VLF wave database. The variation of chorus wave normal angle (WNA) with spatial location and fpe/fce is also taken into account. We find that chorus propagating at small WNAs has the dominant contribution to the diffusion rates in most MLT sectors. However, in the region 4 ≤ MLT < 11 high WNAs dominate at intermediate pitch angles. In the region 3 < L* < 4, the bounce and drift averaged pitch angle and energy diffusion rates during active conditions are primarily larger than those in our earlier models by up to a factor of 10 depending on energy and pitch angle. Further out, the results are similar. We find that the bounce and drift averaged energy and pitch angle diffusion rates can be significantly larger than the new model in regions of low , where the differences can be up to a factor of 10 depending on energy and pitch angle.

Kinetic Modeling of Radiation Belt Electrons With Geant4 to Study Energetic Particle Precipitation in Earth's Atmosphere

AbstractWe present a new model designed to simulate the process of energetic particle precipitation, a vital coupling mechanism from Earth's magnetosphere to its atmosphere. The atmospheric response, namely excess ionization in the upper and middle atmosphere, together with bremsstrahlung X‐ray production, is calculated with kinetic particle simulations using the Geant4 Monte Carlo framework. Mono‐energy and mono‐pitch angle electron beams are simulated and combined using a Green's function approach to represent realistic electron spectra and pitch angle distributions. Results from this model include more accurate ionization profiles than previous analytical models, deeper photon penetration into the atmosphere than previous Monte Carlo model predictions, and predictions of backscatter fractions of loss cone electrons up to 40%. The model results are verified by comparison with previous precipitation modeling results, and validated using balloon X‐ray measurements from the Balloon Array for RBSP Relativistic Electron Losses mission and backscattered electron energy and pitch angle measurements from the Electron Loss and Fields Investigation with a Spatio‐Temporal Ambiguity‐Resolving CubeSat mission. The model results and solution techniques are developed into a Python package for public use.

Spatio‐Temporal Characteristics of IPDP‐Type EMIC Waves on April 19, 2017: Implications for Loss of Relativistic Electrons in the Outer Belt

AbstractTo understand the mechanism of the increased frequency of intervals of pulsations of diminishing periods (IPDPs), we analyzed IPDP‐type electromagnetic ion cyclotron (EMIC) waves that occurred on 19 April 2017, using ground and satellite observations. Observations by low‐altitude satellites and ground‐based magnetometers indicate that the increased IPDP frequency is caused by an inward (i.e., Earthward) shift of the EMIC wave source region. The EMIC wave source region moves inward along the mid‐latitude trough, which we used as a proxy for the plasmapause location. A statistical analysis shows that increases in the IPDP frequency showed a positive correlation with polar cap potentials. These results suggest an enhanced convection electric field causes an inward shift of the source region. The inward shift of the source region allows EMIC waves to scatter relativistic electrons over a wide range of radial distances during the IPDP event. This mechanism suggests that IPDP‐type EMIC waves are more likely to scatter relativistic electrons than other EMIC waves. We also show that the decreased phase‐space density of relativistic electrons in the outer radiation belt is consistent with the extent of the source region and the resonant energy of EMIC waves, implying a possible contribution of EMIC waves to outer radiation belt loss during the main phase of geomagnetic storms.

MLT Dependence of Relativistic Electron Scattering Into the Drift Loss Cone: Measurements From ELFIN‐L on Board Lomonosov Spacecraft

AbstractThere have been a number of theories proposed concerning the loss of relativistic electrons from the radiation belts. However, direct observations of loss were not possible on a number of previous missions due to the large field of view of the instruments and often high‐altitude orbits of satellites that did not allow researchers to isolate the precipitating electrons from the stably trapped. We use measurements from the ELFIN‐L suit of instruments flown on Lomonosov spacecraft at LEO orbit, which allows us to distinguish stably trapped from the drift loss cone electrons. The sun‐synchronous orbit of Lomonosov allows us to quantify scattering that occurred into the loss cone on the dawn‐side and the dusk‐side magnetosphere. The loss at MeV energies is observed predominantly on the dawn‐side, consistent with the loss induced by the chorus waves. The companion data publication provides processed measurements.

Parametric analysis of pitch angle scattering and losses of relativistic electrons by oblique EMIC waves

This study analyzes the effects of electromagnetic ion cyclotron (EMIC) waves on relativistic electron scattering and losses in the Earth’s outer radiation belt. EMIC emissions are commonly observed in the inner magnetosphere and are known to reach high amplitudes, causing significant pitch angle changes in primarily >1 MeV electrons via cyclotron resonance interactions. We run test-particle simulations of electrons streaming through helium band waves with different amplitudes and wave normal angles and assess the sensitivity of advective and diffusive scattering behaviors to these two parameters, including the possibility of very oblique propagation. The numerical analysis confirms the importance of harmonic resonances for oblique waves, and the very oblique waves are observed to efficiently scatter both co-streaming and counter-streaming electrons. However, strong finite Larmor radius effects limit the scattering efficiency at high pitch angles. Recently discussed force-bunching effects and associated strong positive advection at low pitch angles are, surprisingly, shown to cause no decrease in the phase space density of precipitating electrons, and it is demonstrated that the transport of electrons into the loss cone balances out the scattering out of the loss cone. In the case of high-amplitude obliquely propagating waves, weak but non-negligible losses are detected well below the minimum resonance energy, and we identify them as the result of non-linear fractional resonances. Simulations and theoretical analysis suggest that these resonances might contribute to subrelativistic electron precipitation but are likely to be overshadowed by non-resonant effects.

Evidence of H+‐Band EMIC Waves in the Inner Radiation Belt Observed by Van Allen Probes During Magnetic Storms

AbstractElectromagnetic ion cyclotron (EMIC) waves play an important role in the losses of ring current ions and outer radiation belt relativistic electrons. While EMIC waves have been reported to occur frequently in the outer belt and slot region, there is little information about the occurrence of EMIC waves in the inner radiation belt. On basis of the magnetometer observations from Van Allen Probes, we analyze two representative events of H+‐band EMIC waves that occurred in the inner zone on December 21 and June 22, 2015. We find that these waves occurred synchronously accompanying with 5–35 keV proton dispersive injections at the dawn and night sectors. During the periods of EMIC waves, the perpendicular temperature of injected protons increased, while the corresponding parallel component decreased. By analyzing the relation between the temperature anisotropy and parallel plasma beta for the two cases, we suggest that the evident enhancements of the temperature anisotropy (T⊥/T∥) of 5–35 keV protons may provide the energy source to excite the H+‐band EMIC waves. A total of 16 left‐hand H+‐band EMIC wave events have been found in the inner radiation belt during the 4‐year (2013–2016) interval and all occur during magnetic storms. The statistical results implicated that left‐hand H+‐band EMIC waves in the inner radiation belt is associated with storm‐time ion injections.

The dependence of radial diffusion coefficients on solar/interplanetary drivers

Radial diffusion driven by Ultra Low Frequency (ULF) waves is very important for magnetospheric dynamics, as it contributes to both acceleration and loss of relativistic electrons in the outer Van Allen radiation belt, as well as to other planetary magnetospheres. The dependence of ULF wave power spectral density and radial diffusion coefficients (DLL) on solar wind parameters has already been investigated but their dependence on the various solar and interplanetary drivers is poorly studied. In this study, we conduct a statistical analysis of radial diffusion coefficients DLL to look into their dependence on different interplanetary drivers, i.e. Interplanetary Coronal Mass Ejections (ICMEs) and Stream Interaction Regions (SIRs), which originate from different solar drivers: CMEs and coronal holes. We use the ”SafeSpace” database (https://synergasia.uoa.gr/modules/document/?course=PHYS120), which includes radial diffusion coefficients DLL and ULF wave power spectral density, created using magnetic and electric field measurements by the THEMIS satellites in the 2011–2019 time period. We study how the properties of these solar wind drivers influence the behavior of radial diffusion coefficients. The results indicate significant differences between drivers of different solar origin at the ratio of the electric over the magnetic DLL component (DLLE over DLLB), especially when the solar wind dynamic pressure is maximized due to shock existence. Furthermore, differences are observed between SIRs with shock and SIRs without shock. This feature introduces a significant energy dependence to the radial diffusion coefficients, which is further depending on the radial distance and the different first adiabatic invariant (μ) values. These results emphasize the importance of the energy dependence of the magnetic component of the DLL, especially during the occurrence of interplanetary shocks, which is neglected by most semi-empirical models used in radiation belt simulations.

Enhancement in energy loss of relativistic electrons in ferroelectric PVDF above phase transition

Enhancement in energy loss of relativistic electrons in ferroelectric PVDF above phase transition

The “SafeSpace” database of ULF power spectral density and radial diffusion coefficients: dependencies and application to simulations

Abstract. Radial diffusion has been established as one of the most important mechanisms contributing to both the acceleration and loss of relativistic electrons in the outer radiation belt, as well as to the supply of particles to the inner radiation belt. In the framework of the “SafeSpace” project, we have used 9 years (2011–2019) of multi-point magnetic and electric field measurements from THEMIS A, D and E satellites to create a database of radial diffusion coefficients (DLL) and ultra-low-frequency (ULF) wave power spectral densities (PSDs) spanning an L∗ range from 3 to 8. In this work we investigate the dependence of the DLL on the various solar wind parameters, geomagnetic indices and coupling functions, as well as the L-shell, during the solar cycle 24. Moreover, we discuss the uncertainties introduced on the estimation of DLL time series by the partial azimuthal coverage provided by in situ measurements. Furthermore, we investigate, via a superposed analysis, the dependence of the DLL on solar wind drivers. We show, for the first time to the best of our knowledge, that the interplanetary coronal mass ejection (ICME)-driven disturbances accompanied by high solar wind pressure values combined with intense magnetospheric compression can produce DLLB values comparable to or even greater than the ones of DLLE. This feature cannot be captured by semi-empirical models and introduces a significant energy dependence on the DLL. Finally, we show the advantages of using DLL time series by means of numerical simulations of relativistic electron fluxes performed with the Salammbô code and significant deviations in the predictions of several semi-empirical models depending on the level of geomagnetic activity and L-shell.

Runaway electron generation during tokamak start-up

Tokamak start-up is characterized by low electron densities and strong electric fields, in order to quickly raise the plasma current and temperature, allowing the plasma to fully ionize and magnetic flux surfaces to form. Such conditions are ideal for the formation of superthermal electrons, which may reduce the efficiency of ohmic heating and prevent the formation of a healthy thermal fusion plasma. This is of particular concern in ITER where engineering limitations put restrictions on the allowable electric fields and limit the prefill densities during start-up. In this study, we present a new 0D burn-through simulation tool called STREAM (STart-up Runaway Electron Analysis Model), which self-consistently evolves the plasma density, temperature and electric field, while accounting for the generation and loss of relativistic runaway electrons. After verifying the burn-through model, we investigate conditions under which runaway electrons can form during tokamak start-up as well as their effects on the plasma initiation. We find that Dreicer generation plays a crucial role in determining whether a discharge becomes runaway-dominated or not, and that a large number of runaway electrons could limit the ohmic heating of the plasma, thus preventing successful burn-through or further ramp-up of the plasma current. The runaway generation can be suppressed by raising the density via gas fuelling, but only if done sufficiently early. Otherwise a large runaway seed may have already been built up, which can avalanche even at relatively low electric fields and high densities.

On the Magnetosphere‐Ionosphere Coupling During the May 2021 Geomagnetic Storm

AbstractOn 12 May 2021 the interplanetary doppelgänger of the 9 May 2021 coronal mass ejection impacted the Earth's magnetosphere, giving rise to a strong geomagnetic storm. This paper discusses the evolution of the various events linking the solar activity to the Earth's ionosphere with special focus on the effects observed in the circumterrestrial environment. We investigate the propagation of the interplanetary coronal mass ejection and its interaction with the magnetosphere—ionosphere system in terms of both magnetospheric current systems and particle redistribution, by jointly analyzing data from interplanetary, magnetospheric, and low Earth orbiting satellites. The principal magnetospheric current system activated during the different phases of the geomagnetic storm was correctly identified through the direct comparison between geosynchronous orbit observations and model predictions. From the particle point of view, we have found that the primary impact of the storm development is a net and rapid loss of relativistic electrons from the entire outer radiation belt. Our analysis shows no evidence for any short‐term recovery to pre‐storm levels during the days following the main phase. Storm effects also included a small Forbush decrease driven by the interplay between the interplanetary shock and subsequent magnetic cloud arrival.

Hot Plasma Effects on Electron Resonant Scattering by Electromagnetic Ion Cyclotron Waves

AbstractResonant scattering by electromagnetic ion cyclotron (EMIC) waves is one of the most effective mechanisms of relativistic electron losses in Earth’s inner magnetosphere. Low‐altitude spacecraft measurements, however, often show that the energy range of precipitating electrons is wider than theoretical predictions based on the cold plasma dispersion of EMIC waves. To explain this discrepancy, we examine the diffusion rates of EMIC waves by including hot plasma effects in their dispersion relation. Using the observed ion distribution functions, we investigate the hot plasma effects on the EMIC wave dispersion for a wide frequency range. We develop analytical equations for hot plasma effects on EMIC dispersion, and apply this model to diffusion rate evaluations. We show that hot ion effects tend to increase the minimum resonant energy for the frequency range around wave intensity maxima, but can decrease the minimum resonant energy for the higher‐frequency part of wave spectra.

Early-Time Non-Equilibrium Pitch Angle Diffusion of Electrons by Whistler-Mode Hiss in a Plasmaspheric Plume Associated with BARREL Precipitation

In August 2015, the Balloon Array for Radiation belt Relativistic Electron Losses (BARREL) observed precipitation of energetic (<200 keV) electrons magnetically conjugate to a region of dense cold plasma as measured by the twin Van Allen Probes spacecraft. The two spacecraft passed through the high density region during multiple orbits, showing that the structure was spatial and relatively stable over many hours. The region, identified as a plasmaspheric plume, was filled with intense hiss-like plasma waves. We use a quasi-linear diffusion model to investigate plume whistler-mode hiss waves as the cause of precipitation observed by BARREL. The model input parameters are based on the observed wave, plasma and energetic particle properties obtained from Van Allen Probes. Diffusion coefficients are found to be largest in the same energy range as the precipitation observed by BARREL, indicating that the plume hiss waves were responsible for the precipitation. The event-driven pitch angle diffusion simulation is also used to investigate the evolution of the electron phase space density (PSD) for different energies and assumed initial pitch angle distributions. The results show a complex temporal evolution of the phase space density, with periods of both growth and loss. The earliest dynamics, within the ∼5 first minutes, can be controlled by a growth of the PSD near the loss cone (by a factor up to ∼2, depending on the conditions, pitch angle, and energy), favored by the absence of a gradient at the loss cone and by the gradients of the initial pitch angle distribution. Global loss by 1-3 orders of magnitude (depending on the energy) occurs within the first ∼100 min of wave-particle interaction. The prevalence of plasmaspheric plumes and detached plasma regions suggests whistler-mode hiss waves could be an important driver of electron loss even at high L-value (L ∼6), outside of the main plasmasphere.

Conjugate Observation of Magnetospheric Chorus Propagating to the Ionosphere by Ducting

AbstractWhistler‐mode chorus waves are critical for driving resonant scattering and loss of radiation belt relativistic electrons into the atmosphere. The resonant energies of electrons scattered by chorus waves increase at increasingly higher magnetic latitudes. Propagation of chorus waves to middle and high latitudes is hampered by wave divergence and Landau damping but is promoted otherwise if ducted by density irregularities. Although ducting theories have been proposed since the 1960s, no conjugate observation of ducted chorus propagation from the equatorial magnetosphere to the ionosphere has been observed so far. Here we provide such an observation, for the first time, using conjugate spacecraft measurements. Ducted chorus waves maintain significant wave power upon reaching the ionosphere, which is confirmed by ray‐tracing simulations. Our results suggest that ducted chorus waves may be an important driver for relativistic electron precipitation.

Role of Ducting in Relativistic Electron Loss by Whistler‐Mode Wave Scattering

AbstractResonant interactions of energetic electrons with electromagnetic whistler‐mode waves (whistlers) contribute significantly to the dynamics of electron fluxes in Earth's outer radiation belt. At low geomagnetic latitudes, these waves are very effective in pitch angle scattering and precipitation into the ionosphere of low equatorial pitch angle, tens of keV electrons and acceleration of high equatorial pitch angle electrons to relativistic energies. Relativistic (hundreds of keV), electrons may also be precipitated by resonant interaction with whistlers, but this requires waves propagating quasi‐parallel without significant intensity decrease to high latitudes where they can resonate with higher energy low equatorial pitch angle electrons than at the equator. Wave propagation away from the equatorial source region in a non‐uniform magnetic field leads to ray divergence from the originally field‐aligned direction and efficient wave damping by Landau resonance with suprathermal electrons, reducing the wave ability to scatter electrons at high latitudes. However, wave propagation can become ducted along field‐aligned density peaks (ducts), preventing ray divergence and wave damping. Such ducting may therefore result in significant relativistic electron precipitation. We present evidence that ducted whistlers efficiently precipitate relativistic electrons. We employ simultaneous near‐equatorial and ground‐based measurements of whistlers and low‐altitude electron precipitation measurements by ELFIN CubeSat. We show that ducted waves (appearing on the ground) efficiently scatter relativistic electrons into the loss cone, contrary to non‐ducted waves (absent on the ground) precipitating only keV electrons. Our results indicate that ducted whistlers may be quite significant for relativistic electron losses; they should be further studied statistically and possibly incorporated in radiation belt models.

On the Similarity and Repeatability of Fast Radiation Belt Loss: Role of the Last Closed Drift Shell

AbstractProperly characterizing fast relativistic electron losses in the terrestrial Van Allen belts remains a significant challenge for accurately simulating their dynamics. In particular, magnetopause shadowing losses can deplete the radiation belt within hours or even minutes, but can have long‐lasting impacts on the subsequent belt dynamics. By statistically analyzing seven years of data from the entire Van Allen Probes mission in the context of the last closed drift shell, here we show how these losses are much more organized and predictable than previously thought. Once magnetic storm electron dynamics are properly analyzed in terms of the location of the last closed drift shell, not only is the loss shown to be repeatable but its energy‐dependent spatio‐temporal evolution is also revealed to follow a very similar pattern from storm to storm. Employing an energy‐dependent ULF wave radial diffusion model, we further show for the first time how the similar and repeatable fractional loss of the pre‐storm electron population in each storm can be reproduced and explained. Empirical characterization of this loss may open a pathway toward improved radiation belt specification and forecast models. This is especially important since underestimates of loss can also create unrealistic sources in models, creating phantom electron radiation and leading to the prediction of an overly harsh radiation environment.

Interaction between runaway electrons and internal kink in a post-disruption plasma

Direct interactions between the n = 1 (n is the toroidal mode number) resistive internal kink instability and relativistic runaway electrons (REs) in a post thermal quench toroidal plasma are numerically investigated. A recently developed hybrid model, where the runaway current is treated as a fluid component, is adopted to study how REs affect the internal kink stability and the associated eigenmode structure. The RE-modified internal kink perturbation is in turn superposed with the equilibrium field, in order to study the effects of 3D fields on the drift orbits of REs and consequently on the RE confinement and loss. Results are compared with that assuming the single fluid model for the runaway beam. It is found that REs destabilize the resistive internal kink mode, leading to a better recovery (compared to the standard single fluid model) of the mode growth rate scaling of γ ∼ S −1/3 at lower values of Lundquist number S. Furthermore, REs significantly modify the internal kink eigenfunction inside or near the q = 1 surface. In particular, a new narrow layer forms within the resistive layer, where a strong peaking of the perturbed parallel current develops. The RE contribution is also found to substantially enhance the coupling among poloidal harmonics, in particular between m = 1 and its nearest sidebands. All these changes to the perturbation structure consequently affect the RE drift orbits in the presence of the internal kink instability. Compared to the fluid model, 3D perturbations computed with the more consistent hybrid model lead to less loss of relativistic electrons in the RE beam, when the perturbation level reaches that of the equilibrium field. At lower (but still large) perturbation amplitude (∼10% of the equilibrium field), the internal kink mode has little effect on the RE loss. REs launched within the q = 1 surface always stay confined within the plasma, independent of the field perturbation amplitude (up to the level of the equilibrium field), the beam model adopted (fluid versus hybrid), or the particle traveling direction. On the other hand, the individual drift orbit trajectory is qualitatively sensitive to these factors.