Artem Smirnov

and 10 more

Radial diffusion is one of the dominant physical mechanisms driving acceleration and loss of radiation belt electrons. A number of parameterizations for radial diffusion coefficients have been developed, each differing in the dataset used. Here, we investigate the performance of different parameterizations by Brautigam and Albert (2000), Brautigam et al (2005), Ozeke et al. (2014), Ali et al. (2015, 2016); Ali (2016), and Liu et al. (2016) on long-term radiation belt modeling using the Versatile Electron Radiation Belt (VERB) code, and compare the results to Van Allen Probes observations. First, 1-D radial diffusion simulations are performed, isolating the contribution of solely radial diffusion. We then take into account effects of local acceleration and loss showing additional 3-D simulations, including diffusion across pitch-angle and energy, as well as mixed diffusion. For the L* range studied, the difference between simulations with Brautigam and Albert (2000), Ozeke et al. (2014), and Liu et al. (2016) parameterizations is shown to be small, with Brautigam and Albert (2000) offering the best agreement with observations. Using Ali et al. (2016)’s parameterization tended to result in a lower flux at 1 MeV than both the observations and the VERB simulations using the other coefficients. We find that the 3-D simulations are less sensitive to the radial diffusion coefficient chosen than the 1-D simulations, suggesting that for 3-D radiation belt models, a similar result is likely to be achieved, regardless of whether Brautigam and Albert (2000), Ozeke et al. (2014), and Liu et al. (2016) parameterizations are used.
In this study we investigate two distinct loss mechanisms responsible for the rapid dropouts of radiation belt electrons by assimilating data from Van Allen Probes A and B and Geostationary Operational Environmental Satellites (GOES) 13 and 15 into a 3-D diffusion model. In particular, we examine the respective contribution of electromagnetic ion cyclotron (EMIC) wave scattering and magnetopause shadowing for values of the first adiabatic invariant μ ranging from 300 to 3000 MeV G. We inspect the innovation vector and perform a statistical analysis to quantitatively assess the effect of both processes as a function of various geomagnetic indices, solar wind parameters, and radial distance from the Earth. Our results are in agreement with previous studies that demonstrated the energy dependence of these two mechanisms. Loss from L* = 4 to L* = 4.8 is dominated by EMIC wave scattering (μ ≥ 900 MeV G) and may amount to between 10%/hr to 30%/hr of the maximum value of phase space density (PSD) over all L shells for fixed first and second adiabatic invariants. Magnetopause shadowing is shown to deplete electrons across all energies, mostly between L* = 5 and L* = 6.6, resulting in loss from 50%/hr to 70%/hr of the maximum PSD. We also identify a boundary located between L* = 3.5 and L* = 5.2 clearly separating the regions where each mechanism dominates. Nevertheless, during times of enhanced geomagnetic activity, both processes can operate beyond such location and encompass the entire outer radiation belt.