Premium
Modeling radiation belt dynamics using a 3‐D layer method code
Author(s) -
Wang C.,
Ma Q.,
Tao X.,
Zhang Y.,
Teng S.,
Albert J. M.,
Chan A. A.,
Li W.,
Ni B.,
Lu Q.,
Wang S.
Publication year - 2017
Publication title -
journal of geophysical research: space physics
Language(s) - English
Resource type - Journals
eISSN - 2169-9402
pISSN - 2169-9380
DOI - 10.1002/2017ja024143
Subject(s) - van allen radiation belt , computational physics , physics , geomagnetic storm , population , radiation , diffusion , code (set theory) , hiss , van allen probes , pitch angle , earth's magnetic field , electron , geophysics , statistical physics , magnetosphere , computer science , optics , magnetic field , nuclear physics , demography , set (abstract data type) , quantum mechanics , sociology , thermodynamics , programming language
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 A L ‐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.