Huinan Zheng

Modeling of outer radiation belt electrons by multidimensional diffusion process

[1] In this study, we develop a new code which uses a hybrid finite difference (HFD) method to solve a 2-D bounce-averaged energy-pitch angle diffusion equation in the outer radiation belts. We implement the numerical algorithm by adopting a split operator technique, an implicit scheme for dealing with diagonal diffusion coefficients, and an alternative direction implicit scheme for off-diagonal (or cross) diffusion coefficients. We show that our code HFD can successfully overcome the unstable problem when the large and rapidly varying cross diffusion coefficients are included. Particularly, we examine whether resonant interactions with obliquely propagated whistler mode chorus and hiss waves result in a net acceleration or loss of relativistic electrons in the outer radiation belt L = 4.5. Numerical simulations show that chorus waves can yield significant acceleration of relativistic electrons, while hiss waves are primarily responsible for scattering equatorially mirroring electrons toward the loss cone.

CRRES observation and STEERB simulation of the 9 October 1990 electron radiation belt dropout event

We examine and simulate the electron radiation belt dropout event on 9 October 1990. CRRES observations show that significant depletions of electron fluxes occurred at energies similar to 0.1- 1.0 MeV beyond 6 R-E and at energies >similar to 0.4 MeV within 6 R-E. The three- dimensional kinetic radiation belt model STEERB is used to simulate this dropout event, taking into account the magnetopause shadowing, adiabatic transport, radial diffusion, and plume and chorus wave- particle interactions. Our results show that STEERB code can basically reproduce the observed depletion of similar to 0.1-1.0 MeV electron fluxes throughout the outer radiation belt, suggesting that the competition and combination of all these physical mechanisms can well explain this electron radiation belt dropout event. Citation: Su, Z., F. Xiao, H. Zheng, and S. Wang (2011), CRRES observation and STEERB simulation of the 9 October 1990 electron radiation belt dropout event, Geophys. Res. Lett., 38, L06106, doi: 10.1029/2011GL046873.

Radiation belt electron dynamics driven by adiabatic transport, radial diffusion, and wave-particle interactions

The adiabatic transport process is introduced into our recently developed three-dimensional physics-based electron radiation belt model (STEERB, Storm-Time Evolution of Electron Radiation Belt) via adopting a time-varying Hilmer-Voigt geomagnetic field. The current STEERB model contains more complete physical processes: adiabatic transport, radial diffusion, and various in situ wave-particle interactions. In particular, the influence of adiabatic transport on storm time radiation belt electron dynamics is investigated by some idealized simulations. It is found that the adiabatic transport alone (without plume hiss and electromagnetic ion cyclotron (EMIC) waves) is unable to reproduce the observed main phase loss of energetic outer radiation belt electron fluxes in the presence of a strong chorus-driven acceleration process. However, these adiabatic and nonadiabatic processes for radiation belt electron dynamics are coupled to each other. The adiabatic transport, together with radial diffusion and cyclotron resonant interactions with chorus, plume hiss, and EMIC waves, contributes significantly to the main phase loss and the recovery phase enhancement of energetic electron fluxes. In the absence of adiabatic transport, the energetic outer radiation belt electron fluxes are found to be overestimated by a factor of 5-30 over all the pitch angles during the main phase and to be underestimated by a factor of 2-5 at larger pitch angles (alpha(e) > 50 degrees) during the recovery phase. These numerical results suggest that the adiabatic transport in a time-varying geomagnetic field model should be incorporated into the future radiation belt models for space weather application.

Nonstorm time dynamics of electron radiation belts observed by the Van Allen Probes

Storm time electron radiation belt dynamics have been widely investigated for many years. Here we present a rarely reported nonstorm time event of electron radiation belt evolution observed by the Van Allen Probes during 21-24 February 2013. Within 2 days, a new belt centering around L=5.8 formed and gradually merged with the original outer belt, with the enhancement of relativistic electron fluxes by a factor of up to 50. Strong chorus waves (with power spectral density up to 10(-4)nT(2)/Hz) occurred in the region L>5. Taking into account the local acceleration driven by these chorus waves, the two-dimensional STEERB can approximately reproduce the observed energy spectrums at the center of the new belt. These results clearly illustrate the complexity of electron radiation belt behaviors and the importance of chorus-driven local acceleration even during the nonstorm times. Key Points A rarely reported nonstorm time event of RB reformation observed by RBSP Formation of a new belt near the outer boundary of the original outer belt Importance of chorus-driven local acceleration: observation and simulation

Bounce-averaged advection and diffusion coefficients for monochromatic electromagnetic ion cyclotron wave: Comparison between test-particle and quasi-linear models

The electromagnetic ion cyclotron (EMIC) wave has long been suggested to be responsible for the rapid loss of radiation belt relativistic electrons. The test-particle simulations are performed to calculate the bounce-averaged pitch angle advection and diffusion coefficients for parallel-propagating monochromatic EMIC waves. The comparison between test-particle (TP) and quasi-linear (QL) transport coefficients is further made to quantify the influence of nonlinear processes. For typical EMIC waves, four nonlinear physical processes, i.e., the boundary reflection effect, finite perturbation effect, phase bunching and phase trapping, are found to occur sequentially from small to large equatorial pitch angles. The pitch angle averaged finite perturbation effect yields slight differences between the transport coefficients of TP and QL models. The boundary reflection effect and phase bunching produce an average reduction of >80% in the diffusion coefficients but a small change in the corresponding average advection coefficients, tending to lower the loss rate predicted by QL theory. In contrast, the phase trapping causes continuous negative advection toward the loss cone and a minor change in the corresponding diffusion coefficients, tending to increase the loss rate predicted by QL theory. For small amplitude EMIC waves, the transport coefficients grow linearly with the square of wave amplitude. As the amplitude increases, the boundary reflection effect, phase bunching and phase trapping start to occur. Consequently, the TP advection coefficients deviate from the linear growth with the square of wave amplitude, and the TP diffusion coefficients become saturated with the amplitude approaching 1 nT or above. The current results suggest that these nonlinear processes can cause significant deviation of transport coefficients from the prediction of QL theory, which should be taken into account in the future simulations of radiation belt dynamics driven by the EMIC waves.

Ultra-low-frequency wave-driven diffusion of radiation belt relativistic electrons

Van Allen radiation belts are typically two zones of energetic particles encircling the Earth separated by the slot region. How the outer radiation belt electrons are accelerated to relativistic energies remains an unanswered question. Recent studies have presented compelling evidence for the local acceleration by very-low-frequency (VLF) chorus waves. However, there has been a competing theory to the local acceleration, radial diffusion by ultra-low-frequency (ULF) waves, whose importance has not yet been determined definitively. Here we report a unique radiation belt event with intense ULF waves but no detectable VLF chorus waves. Our results demonstrate that the ULF waves moved the inner edge of the outer radiation belt earthward 0.3 Earth radii and enhanced the relativistic electron fluxes by up to one order of magnitude near the slot region within about 10 h, providing strong evidence for the radial diffusion of radiation belt relativistic electrons.

Intense duskside lower band chorus waves observed by Van Allen Probes: Generation and potential acceleration effect on radiation belt electrons

Local acceleration driven by whistler mode chorus waves largely accounts for the enhancement of radiation belt relativistic electron fluxes, whose favored region is usually considered to be the plasmatrough with magnetic local time approximately from midnight through dawn to noon. On 2 October 2013, the Van Allen Probes recorded a rarely reported event of intense duskside lower band chorus waves (with power spectral density up to 10(-3)nT(2)/Hz) in the low-latitude region outside of L=5. Such chorus waves are found to be generated by the substorm-injected anisotropic suprathermal electrons and have a potentially strong acceleration effect on the radiation belt energetic electrons. This event study demonstrates the possibility of broader spatial regions with effective electron acceleration by chorus waves than previously expected. For such intense duskside chorus waves, the occurrence probability, the preferential excitation conditions, the time duration, and the accurate contribution to the long-term evolution of radiation belt electron fluxes may need further investigations in future.

Magnetohydrodynamic simulation of the interaction between two interplanetary magnetic clouds and its consequent geoeffectiveness: 2. Oblique collision

The numerical studies of the interplanetary coupling between multiple magnetic clouds (MCs) are continued by a 2.5-dimensional ideal magnetohydrodynamic (MHD) model in the heliospheric meridional plane. The interplanetary direct collision (DC) / oblique collision (OC) between both MCs results from their same/different initial propagation orientations. Here the OC is explored in contrast to the results of the DC (Xiong et al., 2007). Both the slow MC1 and fast MC2 are consequently injected from the different heliospheric latitudes to form a compound stream during the interplanetary propagation. The MC1 and MC2 undergo contrary deflections during the process of oblique collision. Their deflection angles of

GENERATION OF TYPE III SOLAR RADIO BURSTS IN THE LOW CORONA BY DIRECT AMPLIFICATION

|\delta \theta_1|

MHD simulation of the formation and propagation of multiple magnetic clouds in the heliosphere

and