A. C. Kellerman

Energetic, relativistic, and ultrarelativistic electrons: Comparison of long‐term VERB code simulations with Van Allen Probes measurements

In this study, we compare long-term simulations performed by the Versatile Electron Radiation Belt (VERB) code with observations from the Magnetic Electron Ion Spectrometer and Relativistic Electron-Proton Telescope instruments on the Van Allen Probes satellites. The model takes into account radial, energy, pitch angle and mixed diffusion, losses into the atmosphere, and magnetopause shadowing. We consider the energetic (>100 keV), relativistic (~0.5–1 MeV), and ultrarelativistic (>2 MeV) electrons. One year of relativistic electron measurements (μ = 700 MeV/G) from 1 October 2012 to 1 October 2013 are well reproduced by the simulation during varying levels of geomagnetic activity. However, for ultrarelativistic energies (μ = 3500 MeV/G), the VERB code simulation overestimates electron fluxes and phase space density. These results indicate that an additional loss mechanism is operational and efficient for these high energies. The most likely mechanism for explaining the observed loss at ultrarelativistic energies is scattering by the electromagnetic ion cyclotron waves.

Combined convective and diffusive simulations: VERB-4D comparison with 17 March 2013 Van Allen Probes observations: VERB-4D

This study is focused on understanding the coupling between different electron populations in the inner magnetosphere and the various physical processes that determine evolution of electron fluxes at different energies. Observations during the 17 March 2013 storm and simulations with a newly developed Versatile Electron Radiation Belt-4D (VERB-4D) are presented. Analysis of the drift trajectories of the energetic and relativistic electrons shows that electron trajectories at transitional energies with a first invariant on the scale of ~100 MeV/G may resemble ring current or relativistic electron trajectories depending on the level of geomagnetic activity. Simulations with the VERB-4D code including convection, radial diffusion, and energy diffusion are presented. Sensitivity simulations including various physical processes show how different acceleration mechanisms contribute to the energization of energetic electrons at transitional energies. In particular, the range of energies where inward transport is strongly influenced by both convection and radial diffusion are studied. The results of the 4-D simulations are compared to Van Allen Probes observations at a range of energies including source, seed, and core populations of the energetic and relativistic electrons in the inner magnetosphere.

Wave-induced loss of ultra-relativistic electrons in the Van Allen radiation belts

The dipole configuration of the Earths magnetic field allows for the trapping of highly energetic particles, which form the radiation belts. Although significant advances have been made in understanding the acceleration mechanisms in the radiation belts, the loss processes remain poorly understood. Unique observations on 17 January 2013 provide detailed information throughout the belts on the energy spectrum and pitch angle (angle between the velocity of a particle and the magnetic field) distribution of electrons up to ultra-relativistic energies. Here we show that although relativistic electrons are enhanced, ultra-relativistic electrons become depleted and distributions of particles show very clear telltale signatures of electromagnetic ion cyclotron wave-induced loss. Comparisons between observations and modelling of the evolution of the electron flux and pitch angle show that electromagnetic ion cyclotron waves provide the dominant loss mechanism at ultra-relativistic energies and produce a profound dropout of the ultra-relativistic radiation belt fluxes.

Multi-MeV Electron Loss in the Heart of the Radiation Belts

Significant progress has been made in recent years in understanding acceleration mechanisms in the Earths radiation belts. In particular, a number of studies demonstrated the importance of the local acceleration by analyzing the radial profiles of Phase Space Density (PSD) and observing building up peaks in PSD. In this study, we focus on understanding of the local loss using very similar tools. The profiles of PSD for various values of the first adiabatic invariants during the previously studied January 17, 2013 storm are presented and discussed. The profiles of PSD show clear deepening minimums consistent with the scattering by electromagnetic ion-cyclotron (EMIC) waves. Long-term evolution shows that local minimums in PSD can persist for relatively long times. During considered interval of time the deepening minimums were observed around L* = 4 during January 17, 2013 storm and around L* = 3.5 during March 1, 2013 storm. This study shows a new method that can help identify the location, magnitude and time of the local loss and will help quantify local loss in the future. This study also provides additional clear and definitive evidence that local loss plays a major role for the dynamics of the multi-MeV electrons.

Three‐dimensional data assimilation and reanalysis of radiation belt electrons: Observations of a four‐zone structure using five spacecraft and the VERB code

Obtaining the global state of radiation belt electrons through reanalysis is an important step toward validating our current understanding of radiation belt dynamics and for identification of new physical processes. In the current study, reanalysis of radiation belt electrons is achieved through data assimilation of five spacecraft with the 3-D Versatile Electron Radiation Belt (VERB) code using a split-operator Kalman filter technique. The spacecraft data are cleaned for noise, saturation effects, and then intercalibrated on an individual energy channel basis, by considering phase space density conjunctions in the T96 field model. Reanalysis during the CRRES era reveals a never-before-reported four-zone structure in the Earths radiation belts during the 24 March 1991 shock-induced injection superstorm: (1) an inner belt, (2) the high-energy shock-injection belt, (3) a remnant outer radiation belt, and (4) a second outer radiation belt. The third belt formed near the same time as the second belt and was later enhanced across keV to MeV energies by a second particle injection observed by CRRES and the Northern Solar Terrestrial Array riometer network. During the recovery phase of the storm, the fourth belt was created near L*=4RE, lasting for several days. Evidence is provided that the fourth belt was likely created by a dominant local heating process. This study outlines the necessity to consider all diffusive processes acting simultaneously and the advantage of supporting ground-based data in quantifying the observed radiation belt dynamics. It is demonstrated that 3-D data assimilation can resolve various nondiffusive processes and provides a comprehensive picture of the electron radiation belts.

Electric field control ofEregion coherent echoes: Evidence from radar observations at the South Pole

Characteristics and formation mechanisms of E region plasma irregularities at high latitudes are investigated using observations with the newly deployed Super Dual Auroral Radar Network (SuperDARN) radar at the South Pole Antarctic station (SPS) near a magnetic latitude (MLAT) of 75°S. It is shown that E region echo occurrence at SPS exhibits a diurnal variation that is significantly different from those at auroral and polar cap latitudes. Moreover, analysis of major spectral populations also showed a distinct and previously unreported diurnal pattern. The plasma drift velocity estimates are derived at E region ranges of SPS, leveraging the SPS radars position well within the MLAT region where SuperDARN convection estimates are well constrained by the data. It is shown that E region irregularity occurrence increases when the convection direction is within the SPS field of view and/or when the plasma drift component is comparable with the nominal ion-acoustic speed Cs of 350 m/s. This is the expected behavior for irregularities generated directly by the modified two-stream plasma instability (MTSI). On the other hand, irregularity velocity dependence on convection velocity showed an unexpected saturation at velocity values smaller than nominal Cs. It is demonstrated that the convection velocity at which irregularity velocity starts to differ from the convection component and to approach a maximum value is dependent on the magnetic aspect angle. Moreover, the maximum velocity value itself also depends on the aspect angle. The observed behavior is discussed in context of recent models that involve evolving aspect angles as a key characteristic of MTSI saturation.

EMIC wave parameterization in the long-term VERB code simulation

Electromagnetic ion cyclotron (EMIC) waves play an important role in the dynamics of ultrarelativistic electron population in the radiation belts. However, as EMIC waves are very sporadic, developing a parameterization of such wave properties is a challenging task. Currently, there are no dynamic, activity-dependent models of EMIC waves that can be used in the long-term (several months) simulations, which makes the quantitative modeling of the radiation belt dynamics incomplete. In this study, we investigate Kp, Dst, AE indices, solar wind speed and dynamic pressure as possible parameters of EMIC wave presence. The EMIC waves are included in the long-term simulations (one year, including different geomagnetic activity) performed with the Versatile Electron Radiation Belt (VERB) code, and we compare results of the simulation with the Van Allen Probes observations. The comparison shows that modeling with EMIC waves, parameterized by solar wind dynamic pressure, provides a better agreement with the observations among considered parameterizations. The simulation with EMIC waves improves the dynamics of ultrarelativistic fluxes and reproduces the formation of the local minimum in the phase space density profiles.

Noise statistics identification for Kalman filtering of the electron radiation belt observations I: Model errors

In this study we present a first attempt to identify errors of the 1-D radial diffusion model for relativistic electron phase space density (PSD). In practice, the model error and characteristics of satellite observations are poorly known, which may cause failure of a Kalman filter algorithm. Correct specification of model errors statistics is necessary for the development of the next generation of radiation belt specification models providing the effective PSD reconstruction and hence the prediction and mitigation of space weather effects in the hazardous space environment. The proposed approach to the identification of errors statistics is based on estimating the unknown bias and the covariance matrix of model errors from the sparse CRRES observations over a period of 441 days, from 28 July 1990 to 11 October 1991. With our technique we demonstrate that model errors are biased. Neglecting the bias when applying a data assimilation algorithm to radiation belt electrons can cause significant errors of the PSD estimate during data gaps. Both the identified bias and the covariance matrix of model errors increase with increase of L shell. Sensitivity of the PSD reconstruction to model errors statistics and advances of the improved physical-based model based on the model errors identification are illustrated by a number of representative examples of the PSD reanalysis. Identification of satellite observations characteristics, and filtration and smoothing algorithms are discussed in the companion paper.

Contamination in electron observations of the silicon detector on board Cluster/RAPID/IES instrument in Earth's radiation belts and ring current

Since more than 15 years, the Cluster mission passes through Earths radiation belts at least once every 2 days for several hours, measuring the electron intensity at energies from 30 to 400 keV. These data have previously been considered not usable due to contamination caused by penetrating energetic particles (protons at >100 keV and electrons at >400 keV). In this study, we assess the level of distortion of energetic electron spectra from the Research with Adaptive Particle Imaging Detector (RAPID)/Imaging Electron Spectrometer (IES) detector, determining the efficiency of its shielding. We base our assessment on the analysis of experimental data and a radiation transport code (Geant4). In simulations, we use the incident particle energy distribution of the AE9/AP9 radiation belt models. We identify the Roederer L values, L*, and energy channels that should be used with caution: at 3 600 MeV);at L* similar or equal to 1 and 4-6, the energy channels at 95-400 keV are contaminated by high-energy electrons (>400 keV). Comparison of the data with electron and proton observations from RBSP/MagEIS indicates that the subtraction of proton fluxes at energies similar or equal to 230-630 keV from the IES electron data adequately removes the proton contamination. We demonstrate the usefulness of the corrected data for scientific applications.

Dependence of radiation belt simulations to assumed radial diffusion rates tested for two empirical models of radial transport

Radial diffusion is one of the dominant physical mechanisms that drives acceleration and loss of the radiation belt electrons, which makes it very important for nowcasting and forecasting space weather models. We investigate the sensitivity of the two parameterizations of the radial diffusion of Brautigam and Albert [2000] and Ozeke et al. [2014] on long-term radiation belt modeling using the Versatile Electron Radiation Belt (VERB). Following Brautigam and Albert [2000] and Ozeke et al. [2014], we first perform 1-D radial diffusion simulations. Comparison of the simulation results with observations shows that the difference between simulations with either radial diffusion parameterization is small. To take into account effects of local acceleration and loss, we perform 3-D simulations, including pitch-angle, energy and mixed diffusion. We found that the results of 3-D simulations are even less sensitive to the choice of parameterization of radial diffusion rates than the results of 1-D simulations at various energies (from 0.59 to 1.80 MeV). This result demonstrates that the inclusion of local acceleration and pitch-angle diffusion can provide a negative feedback effect, such that the result is largely indistinguishable simulations conducted with different radial diffusion parameterizations. We also perform a number of sensitivity tests by multiplying radial diffusion rates by constant factors and show that such an approach leads to unrealistic predictions of radiation belt dynamics.