Browsing by Author "Chan, A.A."

Now showing 1 - 4 of 4
  • Thumbnail Image
    Item
    An efficient and positivity-preserving layer method for modeling radiation belt diffusion processes
    (Wiley, 2016) Tao, X.; Zhang, L.; Wang, C.; Li, X.; Albert, J.M.; Chan, A.A.
    An efficient and positivity-preserving layer method is introduced to solve the radiation belt diffusion equation and is applied to study the bounce resonance interaction between relativistic electrons and magnetosonic waves. The layer method with linear interpolation, denoted by LM-L (layer method-linear), requires the use of a large number of grid points to ensure accurate solutions. We introduce a monotonicity- and positivity-preserving cubic interpolation method to be used with the Milstein-Tretyakov layer method. The resulting method, called LM-MC (layer method-monotone cubic), can be used to solve the radiation belt diffusion equation with a much smaller number of grid points than LM-L while still being able to preserve the positivity of the solution. We suggest that LM-MC can be used to study long-term dynamics of radiation belts. We then develop a 2-D LM-MC code and use it to investigate the bounce resonance diffusion of radiation belt electrons by magnetosonic waves. Using a previously published magnetosonic wave model, we demonstrate that bounce resonance with magnetosonic waves is as important as gyroresonance; both can cause several orders of magnitude increase of MeV electron fluxes within 1ﾠday. We conclude that bounce resonance with magnetosonic waves should be taken into consideration together with gyroresonance.
  • Thumbnail Image
    Item
    Dependence of EMIC wave parameters during quiet, geomagnetic storm, and geomagnetic storm phase times
    (Wiley, 2016) Halford, A.J.; Fraser, B.J.; Morley, S.K.; Elkington, S.R.; Chan, A.A.
    As electromagnetic ion cyclotron (EMIC) waves may play an important role in radiation belt dynamics, there has been a push to better include them into global simulations. How to best include EMIC wave effects is still an open question. Recently many studies have attempted to parameterize EMIC waves and their characteristics by geomagnetic indices. However, this does not fully take into account important physics related to the phase of a geomagnetic storm. In this paper we first consider how EMIC wave occurrence varies with the phase of a geomagnetic storm and the SYM-H, AE, and Kp indices. We show that the storm phase plays an important role in the occurrence probability of EMIC waves. The occurrence rates for a given value of a geomagnetic index change based on the geomagnetic condition. In this study we also describe the typical plasma and wave parameters observed in L and magnetic local time for quiet, storm, and storm phase. These results are given in a tabular format in the supporting information so that more accurate statistics of EMIC wave parameters can be incorporated into modeling efforts.
  • Thumbnail Image
    Item
    Effects of magnetic drift shell splitting on electron diffusion in the radiation belts
    (Wiley, 2016) Zheng, Liheng; Chan, A.A.; O’Brien, T.P.; Tu, W.; Cunningham, G.S.; Albert, J.M.; Elkington, S.R.
    Drift shell splitting in the presence of pitch angle scattering breaks all three adiabatic invariants of radiation belt electron motion and produces new diffusion terms that fully populate the diffusion tensor in the Fokker-Planck equation. The Radbelt Electron Model (REM) solves such a Fokker-Planck equation and is used to investigate the phase space density sources. Our simulation results and theoretical arguments suggest that drift shell splitting changes the phase space location of the source to smaller L shells, which typically reduces outer zone phase space density enhancements, and this reduction has a limit corresponding to two-dimensional local diffusion on a curved surface in the phase space.
  • Thumbnail Image
    Item
    Modeling radiation belt dynamics using a 3-D layer method code
    (Wiley, 2017) Wang, C.; Ma, Q.; Tao, X.; Zhang, Y.; Teng, S.; Albert, J.M.; Chan, A.A.; Li, W.; Ni, B.; Lu, Q.; Wang, S.
    A new 3-D diffusion code using a recently published layer method has been developed to analyze radiation belt electron dynamics. The code guarantees the positivity of the solution even when mixed diffusion terms are included. Unlike most of the previous codes, our 3-D code is developed directly in equatorial pitch angle (α0), momentum (p), and L shell coordinates; this eliminates the need to transform back and forth between (α0,p) coordinates and adiabatic invariant coordinates. Using (α0,p,L) is also convenient for direct comparison with satellite data. The new code has been validated by various numerical tests, and we apply the 3-D code to model the rapid electron flux enhancement following the geomagnetic storm on 17 March 2013, which is one of the Geospace Environment Modeling Focus Group challenge events. An event-specific global chorus wave model, an AL-dependent statistical plasmaspheric hiss wave model, and a recently published radial diffusion coefficient formula from Time History of Events and Macroscale Interactions during Substorms (THEMIS) statistics are used. The simulation results show good agreement with satellite observations, in general, supporting the scenario that the rapid enhancement of radiation belt electron flux for this event results from an increased level of the seed population by radial diffusion, with subsequent acceleration by chorus waves. Our results prove that the layer method can be readily used to model global radiation belt dynamics in three dimensions.