Van Allen radiation belts consist of relativistic electrons trapped by Earth's magnetic field. Trapped electrons often drift azimuthally around Earth and display a butterfly pitch angle distribution of a minimum at 90° further out than geostationary orbit. This is usually attributed to drift shell splitting resulting from day–night asymmetry in Earth's magnetic field. However, direct observation of a butterfly distribution well inside of geostationary orbit and the origin of this phenomenon have not been provided so far. Here we report high-resolution observation that a unusual butterfly pitch angle distribution of relativistic electrons occurred within 5 Earth radii during the 28 June 2013 geomagnetic storm. Simulation results show that combined acceleration by chorus and magnetosonic waves can successfully explain the electron flux evolution both in the energy and butterfly pitch angle distribution. The current provides a great support for the mechanism of wave-driven butterfly distribution of relativistic electrons.

The recent launching of Van Allen probes provides an unprecedent opportunity to investigate variations of the radiation belt relativistic electrons. During the 17–19 March 2013 storm, the Van Allen probes simultaneously detected strong chorus waves and substantial increases in fluxes of relativistic (2 − 4.5 MeV) electrons around L = 4.5. Chorus waves occurred within the lower band 0.1–0.5fce (the electron equatorial gyrofrequency), with a peak spectral density ∼10−4 nT2/Hz. Correspondingly, relativistic electron fluxes increased by a factor of 102–103 during the recovery phase compared to the main phase levels. By means of a Gaussian fit to the observed chorus spectra, the drift and bounce‐averaged diffusion coefficients are calculated and then used to solve a 2‐D Fokker‐Planck diffusion equation. Numerical simulations demonstrate that the lower‐band chorus waves indeed produce such huge enhancements in relativistic electron fluxes within 15 h, fitting well with the observation.

[1] We report correlated observation of enhanced electromagnetic ion cyclotron (EMIC) waves and dynamic evolution of ring current proton flux collected by Cluster satellite near the location L = 4.5 during March 26-27, 2003, a nonstorm period (D st > À10). Energetic (5-30 keV) proton fluxes are found to drop rapidly (e.g., a half hour) at lower pitch angles, corresponding to intensified EMIC wave activities. By adopting a Gaussian fit to the observed spectra of EMIC waves, we present two-dimensional (2D) numerical simulations which demonstrate that EMIC wave can yield such decrements in proton flux within 30 minutes, consistent with the observational data. The current result provides a further understanding of ring current dynamics driven by wave-particle interaction under different geomagnetic activities.

[1] Simultaneous observations of enhanced fast magnetosonic (MS) waves and distinct proton ring distributions collected by Cluster satellite near the location L = 4-5 on 28 May 2005 are analyzed to study instability of MS waves. A sum of subtracted bi-Maxwellian components is utilized to fit the observed proton (2-10 keV) ring distributions. A ray-tracing simulation is performed to calculate the local growth rate and path-integrated gain of MS waves. Peak growth rates are found to occur at the multiples of proton gyrofrequency mainly in the range 70-120 Hz, and wave gain lies in 40-80 dB, comparable to the observation. Moreover, MS waves primarily locate within a few degrees of the geomagnetic equator and propagate either into or out of the plasmasphere through the plasmapause. The current results provide observational support for instability of fast magnetosonic waves generated by the proton ring distribution.

During the small storm on 14-15 April 2014, Van Allen Probe A measured a continuously distinct proton ring distribution and enhanced magnetosonic (MS) waves along its orbit outside the plasmapause. Inside the plasmasphere, strong MS waves were still present but the distinct proton ring distribution was falling steeply with distance. We adopt a sum of subtracted bi-Maxwellian components to model the observed proton ring distribution and simulate the wave trajectory and growth. MS waves at first propagate toward lower L shells outside the plasmasphere, with rapidly increasing path gains related to the continuous proton ring distribution. The waves then gradually cross the plasmapause into the deep plasmasphere, with almost unchanged path gains due to the falling proton ring distribution and higher ambient density. These results present the first report on how MS waves penetrate into the plasmasphere with the aid of the continuous proton ring distributions during weak geomagnetic activities.