Search for: All records

Magnetopause shadowing (MPS) effect could drive a concurrent dropout of radiation belt electrons and ring current protons. However, its relative role in the dropout of both plasma populations has not been well quantified. In this work, we study the simultaneous dropout of MeV electrons and 100s keV protons during an intense geomagnetic storm in May 2017. A radial diffusion model with an event‐specific last closed drift shell is used to simulate the MPS loss of both populations. The model well captures the fast shadowing loss of both populations atL* > 4.6, while the loss atL* < 4.6, possibly due to the electromagnetic ion cyclotron wave scattering, is not captured. The observed butterfly pitch angle distributions of electron fluxes in the initial loss phase are well reproduced by the model. The initial proton losses at low pitch angles are underestimated, potentially also contributed by other mechanisms such as field line curvature scattering.

Last closed drift shell (LCDS) has been identified as a crucial parameter for investigating the magnetopause shadowing loss of radiation belt electrons. However, drift orbit bifurcation (DOB) effects have not been physically incorporated into the LCDS calculation. Here we calculate event‐specific LCDS using different approaches to dealing with the DOB effects, that is, tracing field lines ignoring DOB, tracing test particles rejecting field lines with DOB, and tracing particles including field lines with DOB, and then incorporate them into a radial diffusion model to simulate the fast electron dropout observed by Van Allen Probes in May 2017. The model effectively captures the fast dropout at highL*and exhibits the best agreement with data when LCDS is calculated by tracing test particles with DOB more physically included. This study represents the first quantitative modeling of the DOB effects on radiation belt magnetopause shadowing loss via a more physical specification of LCDS.

During the 9 March 2018 event with two consecutive interplanetary shocks compressing the dayside magnetosphere, the azimuthal mode structure and frequency spectrum of ultra low frequency magnetic pulsations are resolved using a cross‐spectral analysis based on high‐fidelity multi‐probe Magnetospheric Multiscale Mission (MMS) magnetometer data. The results based on the MMS 4 and MMS 3 pair of measurements show that shock arrival leads to low mode () magnetic fluctuations in the Pc4‐5 regimes, and smaller spatial scale fluctuations implied by the dominant high mode numbers are observed after both shock signatures hit and passed the magnetosphere. Detailed evolution of the mode structure is also shown for the first shock to reveal the development of high mode structure from a bump‐on‐tail distribution atto a dominant peak atin about 10 min. In addition, an interesting change of sign infrom negative to positive is observed as MMS crosses ∼11 MLT pre‐noon, which is consistent with the picture of wave generation by dayside magnetopause compression and then anti‐sunward propagation. For both shocks, the contribution of higher frequency waves (Pc‐4 range compared with Pc‐5) to the total wave power is found to be negligible before and after the shock impact, but it becomes more significant during the shock impact.

Ultra‐low‐frequency (ULF) waves are known to radially diffuse hundreds‐keV to few‐MeV electrons in the magnetosphere, as the range of drift frequencies of such electrons overlaps with the frequencies of the waves, leading to resonant interactions. The theoretical framework for this process is described by analytic expressions of the resonant interactions between electrons and toroidal and poloidal ULF wave modes in a background magnetic field. However, most expressions estimate the radial diffusion rates based on estimates of the power of ULF waves that are obtained either from spacecraft close to the equatorial plane or from the ground. In this study, using multiyear measurements from the THEMIS and Arase missions, we present a statistical analysis of the distribution of ULF wave power in magnetic latitude and local time and show that the wave power of the radial and azimuthal components of the magnetic field increases away from the magnetic equator. Our result could have significant implications for the radial diffusion rates as currently estimated.

Angular response functions are derived for four electron channels and six proton channels of the SEM‐2 MEPED particle telescopes on the POES and MetOp satellites from Geant4 simulations previously used to derive the energy response. They are combined with model electron distributions in energy and pitch angle to show that the vertical 0° telescope, intended to measure precipitating electrons, instead usually measures trapped or quasi‐trapped electrons, except during times of enhanced pitch angle diffusion. A simplified dynamical model of the radiation belt electron distribution near the loss cone, as a function of longitude, energy, and pitch angle, that accounts for pitch angle diffusion, azimuthal drift, and atmospheric backscatter is fit to sample MEPED electron data atL = 4during times of differing diffusion rates. It is then used to compute precipitating electron flux, as function of energy and longitude, that is lower than would be estimated by assuming that the 0° telescope always measures precipitating electrons.

Radiation belt electrons undergo frequent acceleration, transport, and loss processes under various physical mechanisms. One of the most prevalent mechanisms is radial diffusion, caused by the resonant interactions between energetic electrons and ULF waves in the Pc4‐5 band. An indication of this resonant interaction is believed to be the appearance of periodic flux oscillations. In this study, we report long‐lasting, drift‐periodic flux oscillations of relativistic and ultrarelativistic electrons with energies up to ∼7.7 MeV in the outer radiation belt, observed by the Van Allen Probes mission. During this March 2017 event, multi‐MeV electron flux oscillations at the electron drift frequency appeared coincidently with enhanced Pc5 ULF wave activity and lasted for over 10 h in the center of the outer belt. The amplitude of such flux oscillations is well correlated with the radial gradient of electron phase space density (PSD), with almost no oscillation observed near the PSD peak. The temporal evolution of the PSD radial profile also suggests the dominant role of radial diffusion in multi‐MeV electron dynamics during this event. By combining these observations, we conclude that these multi‐MeV electron flux oscillations are caused by the resonant interactions between electrons and broadband Pc5 ULF waves and are an indicator of the ongoing radial diffusion process during this event. They contain essential information of radial diffusion and have the potential to be further used to quantify the radial diffusion effects and aid in a better understanding of this prevailing mechanism.

The “Quantitative Assessment of Radiation Belt Modeling” focus group was in place at Geospace Environment Modeling from 2014 to 2018. The overarching goals of this focus group were to bring together the current state‐of‐the‐art models for the acceleration, transport, and loss processes in Earth's radiation belts; develop event‐specific and global inputs of wave, plasma, and magnetic field to drive these models; and combine all these components to achieve a quantitative assessment of radiation belt modeling by validating against contemporary radiation belt measurements. This article briefly reviews the current understanding of radiation belt dynamics and related modeling efforts, summarizes the activities and accomplishments of the focus group, and discusses future directions.

During geomagnetic storms, the magnetic field in the outer belt is known to be significantly distorted by the solar wind, such that the outer belt dynamics can be greatly impacted by the magnetic field topology. In this context, we develop a reduced Fokker‐Planck code that includes the effects related to realistic magnetic field models (by using the existing dedicated LANLGeoMag library) through the steps of preprocessing, postprocessing, and, for the first time, during the computation of the reduced Fokker‐Planck diffusion equation itself. We perform the solutions of the reduced Fokker‐Planck equation in the framework of the geomagnetic storm that occurred from 9 October to 15 October 1990. With the use of Combined Release and Radiation Effects Satellite observations, the magnetic field model is shown to strongly affect the way of conciliating theory with observations (processing steps). More specifically, we explain analytically and numerically why the use of a dipole field can lead to misleading interpretations on the local enhancements (attributed to local acceleration) displayed by the electron distribution function, resulting in inaccurate simulations results at large L‐shells. The consideration of a realistic field does not produce any artificial peaks. With such corrected data sets, a great part of the dynamics can be described by radial transport and is thus better reproduced by the simulations. This crucial importance of the field geometry is further emphasized with the calculation of unidirectional, omnidirectional, and integral electron fluxes, and their accuracy is quantified thanks to dedicated metrics.

Characterizing the azimuthal mode number,m, of ultralow‐frequency (ULF) waves is necessary for calculating radial diffusion of radiation belt electrons. A cross‐spectral technique is applied to the compressional Pc5 ULF waves observed by multiple pairs of GOES satellites to estimate the azimuthal mode structure during the 28‐31 May 2010 storm. We find that allowing for both positive and negativemis important to achieve a more realistic distribution of mode numbers and to resolve wave propagation direction. During the storm commencement when the solar wind dynamic pressure is high, ULF wave power is found to dominate at low‐mode numbers. An interesting change of sign inmoccurred around noon, which is consistent with the driving of ULF waves by solar wind buffeting around noon, creating antisunward wave propagation. The low‐mode ULF waves are also found to have a less global coverage in magnetic local time than previously assumed. In contrast, during the storm main phase and early recovery phase when the solar wind dynamic pressure is low and the auroral electrojet index is high, wave power is shown to be distributed over all modes from low to high. The high‐mode waves are found to cover a wider range of magnetic local time than what was previously assumed. Furthermore, to reduce the 2nπambiguity in resolvingm, a cross‐pair analysis is performed on satellite field measurements for the first time, which is demonstrated to be effective in generating more reliable mode structure of ULF waves during high auroral electrojet periods.

Using the global Lagrangian version of the piecewise parabolic method‐magnetohydrodynamic (PPMLR‐MHD) model, we simulate two consecutive storms in December 2015, a moderate storm on 14–15 December and a strong storm on 19–22 December, and calculate the radial diffusion coefficients (D_LL) from the simulated ultralow frequency waves. We find that even though the strong storm leads to more enhancedB_zandE_φpower than the moderate storm, the two storms share in common a lot of features on the azimuthal mode structure and power spectrum of ultralow frequency waves. For both storms, the totalB_zandE_φpower is better correlated with the solar wind dynamic pressure in the storm initial phase and more correlated withAEindex in the recovery phase.B_zwave power is shown to be mostly distributed in low mode numbers, whileE_φpower spreads over a wider range of modes. Furthermore, theB_zandE_φpower spectral densities are found to be higher at higherLregions, with a strongerLdependence in theB_zspectra. The estimatedD_LLbased on MHD fields shows that inside the magnetopause, the contribution from electric fields is larger than or comparable to that from magnetic fields, and our event‐specific MHD‐basedD_LLcan be smaller than some previous empiricalD_LLestimations by more than an order of magnitude. At last, by validating against in situ observations from Magnetospheric Multiscale spacecraft, our MHD results are found to generally well reproduce the totalB_zfields and wave power for both storms, while theE_φpower is underestimated in the MHD simulations.