Refine
Document Type
- Article (15)
- Doctoral Thesis (1)
- Postprint (1)
Language
- English (17)
Is part of the Bibliography
- yes (17)
Keywords
- electrons (3)
- magnetosphere (3)
- radiation belts (3)
- ring current electrons (3)
- wave-particle interactions (3)
- L-Asterisk (2)
- inner magnetosphere (2)
- model (2)
- Datenassimilation (1)
- EMIC (1)
Institute
Simulations of the inner magnetospheric energetic electrons using the IMPTAM-VERB coupled model
(2019)
In this study, we present initial results of the coupling between the Inner Magnetospheric Particle Transport and Acceleration Model (IMPTAM) and the Versatile Electron Radiation Belt (VERB-3D) code. IMPTAM traces electrons of 10-100 keV energies from the plasma sheet (L = 9 Re) to inner L-shell regions. The flux evolution modeled by IMPTAM is used at the low energy and outer L* computational boundaries of the VERB code (assuming a dipole approximation) to perform radiation belt simulations of energetic electrons. The model was tested on the March 17th, 2013 storm, for a six-day period. Four different simulations were performed and their results compared to satellites observations from Van Allen probes and GOES. The coupled IMPTAM-VERB model reproduces evolution and storm-time features of electron fluxes throughout the studied storm in agreement with the satellite data (within similar to 0.5 orders of magnitude). Including dynamics of the low energy population at L* = 6.6 increases fluxes closer to the heart of the belt and has a strong impact in the VERB simulations at all energies. However, inclusion of magnetopause losses leads to drastic flux decreases even below L* = 3. The dynamics of low energy electrons (max. 10s of keV) do not affect electron fluxes at energies >= 900 keV. Since the IMPTAM-VERB coupled model is only driven by solar wind parameters and the Dst and Kp indexes, it is suitable as a forecasting tool. In this study, we demonstrate that the estimation of electron dynamics with satellite-data-independent models is possible and very accurate.
Plain Language Summary Radiation belts of the Earth, which are the zones of charged energetic particles trapped by the geomagnetic field, comprise enormous and dynamic systems. While the inner radiation belt, composed mainly of high-energy protons, is relatively stable, the outer belt, filled with energetic electrons, is highly variable and depends substantially on solar activity. Hence, extended reliable observations and the improved models of the electron intensities in the outer belt depending on solar wind parameters are necessary for prediction of their dynamics. The Cluster mission has been measuring electron flux intensities in the radiation belts since its launch in 2000, thus providing a huge dataset that can be used for radiation belts analysis. Using 16 years of electron measurements by the Cluster mission corrected for background contamination, we derived a uniform linear-logarithmic dependence of electron fluxes in the outer belt on the solar wind dynamic pressure.
The Cluster mission has produced a large data set of electron flux measurements in the Earth's magnetosphere since its launch in late 2000. Electron fluxes are measured using Research with Adaptive Particle Imaging Detector (RAPID)/Imaging Electron Spectrometer (IES) detector as a function of energy, pitch angle, spacecraft position, and time. However, no adiabatic invariants have been calculated for Cluster so far. In this paper we present a step-by-step guide to calculations of adiabatic invariants and conversion of the electron flux to phase space density (PSD) in these coordinates. The electron flux is measured in two RAPID/IES energy channels providing pitch angle distribution at energies 39.2-50.5 and 68.1-94.5 keV in nominal mode since 2004. A fitting method allows to expand the conversion of the differential fluxes to the range from 40 to 150 keV. Best data coverage for phase space density in adiabatic invariant coordinates can be obtained for values of second adiabatic invariant, K, similar to 10(2), and values of the first adiabatic invariant mu in the range approximate to 5-20 MeV/G. Furthermore, we describe the production of a new data product "LSTAR," equivalent to the third adiabatic invariant, available through the Cluster Science Archive for years 2001-2018 with 1-min resolution. The produced data set adds to the availability of observations in Earth's radiation belts region and can be used for long-term statistical purposes.
Significant progress has been made in recent years in understanding acceleration mechanisms in the Earth's 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 17 January 2013 storm are presented and discussed. The profiles of PSD show clear deepening minimums consistent with the scattering by electromagnetic ion cyclotron 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 17 January 2013 storm and around L* = 3.5 during 1 March 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.
Models of ring current electron dynamics unavoidably contain uncertainties in boundary conditions, electric and magnetic fields, electron scattering rates, and plasmapause location. Model errors can accumulate with time and result in significant deviations of model predictions from observations. Data assimilation offers useful tools which can combine physics-based models and measurements to improve model predictions. In this study, we systematically analyze performance of the Kalman filter applied to a log-transformed convection model of ring current electrons and Van Allen Probe data. We consider long-term dynamics of mu = 2.3 MeV/G and K = 0.3 G(1/2) R-E electrons from 1 February 2013 to 16 June 2013. By using synthetic data, we show that the Kalman filter is capable of correcting errors in model predictions associated with uncertainties in electron lifetimes, boundary conditions, and convection electric fields. We demonstrate that reanalysis retains features which cannot be fully reproduced by the convection model such as storm-time earthward propagation of the electrons down to 2.5 R-E. The Kalman filter can adjust model predictions to satellite measurements even in regions where data are not available. We show that the Kalman filter can adjust model predictions in accordance with observations for mu = 0.1, 2.3, and 9.9 MeV/G and constant K = 0.3 G(1/2) R-E electrons. The results of this study demonstrate that data assimilation can improve performance of ring current models, better quantify model uncertainties, and help deeper understand the physics of the ring current particles.
Storm Time Depletions of Multi-MeV Radiation Belt Electrons Observed at Different Pitch Angles
(2019)
During geomagnetic storms, the rapid depletion of the high-energy (several MeV) outer radiation belt electrons is the result of loss to the interplanetary medium through the magnetopause, outward radial diffusion, and loss to the atmosphere due to wave-particle interactions. We have performed a statistical study of 110 storms using pitch angle resolved electron flux measurements from the Van Allen Probes mission and found that inside of the radiation belt (L* = 3 - 5) the number of storms that result in depletion of electrons with equatorial pitch angle alpha(eq) = 30 degrees is higher than number of storms that result in depletion of electrons with equatorial pitch angle alpha(eq) = 75 degrees. We conclude that this result is consistent with electron scattering by whistler and electromagnetic ion cyclotron waves. At the outer edge of the radiation belt (L* >= 5.2) the number of storms that result in depletion is also large (similar to 40-50%), emphasizing the significance of the magnetopause shadowing effect and outward radial transport.
Ring current electrons (1–100 keV) have received significant attention in recent decades, but many questions regarding their major transport and loss mechanisms remain open. In this study, we use the four‐dimensional Versatile Electron Radiation Belt code to model the enhancement of phase space density that occurred during the 17 March 2013 storm. Our model includes global convection, radial diffusion, and scattering into the Earth's atmosphere driven by whistler‐mode hiss and chorus waves. We study the sensitivity of the model to the boundary conditions, global electric field, the electric field associated with subauroral polarization streams, electron loss rates, and radial diffusion coefficients. The results of the code are almost insensitive to the model parameters above 4.5 RERE, which indicates that the general dynamics of the electrons between 4.5 RE and the geostationary orbit can be explained by global convection. We found that the major discrepancies between the model and data can stem from the inaccurate electric field model and uncertainties in lifetimes. We show that additional mechanisms that are responsible for radial transport are required to explain the dynamics of ≥40‐keV electrons, and the inclusion of the radial diffusion rates that are typically assumed in radiation belt studies leads to a better agreement with the data. The overall effect of subauroral polarization streams on the electron phase space density profiles seems to be smaller than the uncertainties in other input parameters. This study is an initial step toward understanding the dynamics of these particles inside the geostationary orbit.
Chorus waves play an important role in the dynamic evolution of energetic electrons in the Earth's radiation belts and ring current. Using more than 5 years of Van Allen Probe data, we developed a new analytical model for upper‐band chorus (UBC; 0.5fce < f < fce) and lower‐band chorus (LBC; 0.05fce < f < 0.5fce) waves, where fce is the equatorial electron gyrofrequency. By applying polynomial fits to chorus wave root mean square amplitudes, we developed regression models for LBC and UBC as a function of geomagnetic activity (Kp), L, magnetic latitude (λ), and magnetic local time (MLT). Dependence on Kp is separated from the dependence on λ, L, and MLT as Kp‐scaling law to simplify the calculation of diffusion coefficients and inclusion into particle tracing codes. Frequency models for UBC and LBC are also developed, which depends on MLT and magnetic latitude. This empirical model is valid in all MLTs, magnetic latitude up to 20°, Kp ≤ 6, L‐shell range from 3.5 to 6 for LBC and from 4 to 6 for UBC. The dependence of root mean square amplitudes on L are different for different bands, which implies different energy sources for different wave bands. This analytical chorus wave model is convenient for inclusion in quasi‐linear diffusion calculations of electron scattering rates and particle simulations in the inner magnetosphere, especially for the newly developed four‐dimensional codes, which require significantly improved wave parameterizations.
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 mu ranging from 300 to 3,000 MeV G(-1). 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. We show that EMIC wave scattering tends to dominate loss at lower L shells, and it may amount to between 10%/hr and 30%/hr of the maximum value of phase space density (PSD) over all L shells for fixed first and second adiabatic invariants. On the other hand, magnetopause shadowing is found to deplete electrons across all energies, mostly at higher L shells, resulting in loss from 50%/hr to 70%/hr of the maximum PSD. Nevertheless, during times of enhanced geomagnetic activity, both processes can operate beyond such location and encompass the entire outer radiation belt.
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, and 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 (1year, including different geomagnetic activity) performed with the Versatile Electron Radiation Belt 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.