Binbin Ni

Resonant scattering of outer zone relativistic electrons by multiband EMIC waves and resultant electron loss time scales

To improve our understanding of the role of electromagnetic ion cyclotron (EMIC) waves in radiation belt electron dynamics, we perform a comprehensive analysis of EMIC wave-induced resonant scattering of outer zone relativistic (>0.5 MeV) electrons and resultant electron loss time scales with respect to EMIC wave band, L shell, and wave normal angle model. The results demonstrate that while H+-band EMIC waves dominate the scattering losses of ~1–4 MeV outer zone relativistic electrons, it is He+-band and O+-band waves that prevail over the pitch angle diffusion of ultrarelativistic electrons at higher energies. Given the wave amplitude, EMIC waves at higher L shells tend to resonantly interact with a larger population of outer zone relativistic electrons and drive their pitch angle scattering more efficiently. Obliquity of EMIC waves can reduce the efficiency of wave-induced relativistic electron pitch angle scattering. Compared to the frequently adopted parallel or quasi-parallel model, use of the latitudinally varying wave normal angle model produces the largest decrease in H+-band EMIC wave scattering rates at pitch angles  ~5 MeV. At a representative nominal amplitude of 1 nT, EMIC wave scattering produces the equilibrium state (i.e., the lowest normal mode under which electrons at the same energy but different pitch angles decay exponentially on the same time scale) of outer belt relativistic electrons within several to tens of minutes and the following exponential decay extending to higher pitch angles on time scales from <1 min to ~1 h. The electron loss cone can be either empty as a result of the weak diffusion or heavily/fully filled due to approaching the strong diffusion limit, while the trapped electron population at high pitch angles close to 90° remains intact because of no resonant scattering. In this manner, EMIC wave scattering has the potential to deepen the anisotropic distribution of outer zone relativistic electrons by reshaping their pitch angle profiles to “top-hat.” Overall, H+-band and He+-band EMIC waves are most efficient in producing the pitch angle scattering loss of relativistic electrons at ~1–2 MeV. In contrast, the presence of O+-band EMIC waves, while at a smaller occurrence rate, can dominate the scattering loss of 5–10 MeV electrons in the entire region of the outer zone, which should be considered in future modeling of the outer zone relativistic electron dynamics.

Modeling inward diffusion and slow decay of energetic electrons in the Earth's outer radiation belt

©2015. American Geophysical Union. All Rights Reserved. A new 3-D diffusion code is used to investigate the inward intrusion and slow decay of energetic radiation belt electrons (>0.5MeV) observed by the Van Allen Probes during a 10day quiet period on March 2013. During the inward transport, the peak differential electron fluxes decreased by approximately an order of magnitude at various energies. Our 3-D radiation belt simulation including radial diffusion and pitch angle and energy diffusion by plasmaspheric hiss and electromagnetic ion cyclotron (EMIC) waves reproduces the essential features of the observed electron flux evolution. The decay time scales and the pitch angle distributions in our simulation are consistent with the Van Allen Probe observations over multiple energy channels. Our study suggests that the quiet time energetic electron dynamics are effectively controlled by inward radial diffusion and pitch angle scattering due to a combination of plasmaspheric hiss and EMIC waves in the Earths radiation belts.

Formation of energetic electron butterfly distributions by magnetosonic waves via Landau resonance

Radiation belt electrons can exhibit different types of pitch angle distributions in response to various magnetospheric processes. Butterfly distributions, characterized by flux minima at pitch angles around 90°, are broadly observed in both the outer and inner belts and the slot region. Butterfly distributions close to the outer magnetospheric boundary have been attributed to drift shell splitting and losses to the magnetopause. However, their occurrence in the inner belt and the slot region has hitherto not been resolved. By analyzing the particle and wave data collected by the Van Allen Probes during a geomagnetic storm, we combine test particle calculations and Fokker-Planck simulations to reveal that scattering by equatorial magnetosonic waves is a significant cause for the formation of energetic electron butterfly distributions in the inner magnetosphere. Another event shows that a large-amplitude magnetosonic wave in the outer belt can create electron butterfly distributions in just a few minutes.

Analytical approximation of transit time scattering due to magnetosonic waves

© 2015. American Geophysical Union. All Rights Reserved. Recent test particle simulations have shown that energetic electrons traveling through fast magnetosonic (MS) wave packets can experience an effect which is specifically associated with the tight equatorial confinement of these waves, known as transit time scattering. However, such test particle simulations can be computationally cumbersome and offer limited insight into the dominant physical processes controlling the wave-particle interactions, that is, in determining the effects of the various wave parameters and equatorial confinement on the particle scattering. In this paper, we show that such nonresonant effects can be effectively captured with a straightforward analytical treatment that is made possible with a set of reasonable, simplifying assumptions. It is shown that the effect of the wave confinement, which is not captured by the standard quasi-linear theory approach, acts in such a way as to broaden the range of particle energies and pitch angles that can effectively resonate with the wave. The resulting diffusion coefficients can be readily incorporated into global diffusion models in order to test the effects of transit time scattering on the dynamical evolution of radiation belt fluxes.

Predominance of ECH wave contribution to diffuse aurora in Earth's outer magnetosphere

©2014. American Geophysical Union. All Rights Reserved. Due to its importance for global energy dissipation in the ionosphere, the diffuse aurora has been intensively studied in the past 40 years. Its origin (precipitation of 0.5-10 keV electrons from the plasma sheet without potential acceleration) has been generally attributed to whistler-mode chorus wave scattering in the inner magnetosphere (R < ∼8 RE), while the scattering mechanism beyond that distance remains unresolved. By modeling the quasi-linear diffusion of electrons with realistic parameters for the magnetic field, loss cone size, and wave intensity (obtained from Time History of Events and Macroscale Interactions during Substorms (THEMIS) observations as a function of magnetospheric location), we estimate the loss cone filling ratio and electron cyclotron harmonic (ECH) wave-induced electron precipitation systematically throughout the entire data set, from 6 RE out to 31 RE (the THEMIS apogee). By comparing the wave-induced precipitation directly with the equatorially mapped energy flux distribution of the diffuse aurora from ionospheric observations (OVATION Prime model) at low altitudes, we quantify the contribution of auroral energy flux precipitated due to ECH wave scattering. Although the wave amplitudes decrease, as expected, with distance from the Earth, due to the smaller loss cone size and stretched magnetic field topology, ECH waves are still capable of causing sufficient scattering of plasma sheet electrons to account for the observed diffuse auroral dissipation. Our results demonstrate that ECH waves are the dominant driver of the diffuse aurora in the outer magnetosphere, beyond ∼8 RE.

Hemispheric asymmetry of the structure of dayside auroral oval

A comprehensive analysis of long-term and multispectral auroral observations made in the Arctic and Antarctica demonstrates that the dayside auroral ovals in two hemispheres are both presented in a two-peak structure, namely, the prenoon 09:00 magnetic local time (MLT) and postnoon 15:00 MLT peaks. The two-peak structures of dayside ovals, however, are asymmetric in the two hemispheres; i.e., the postnoon average auroral intensity is more than the prenoon one in the Northern Hemisphere but less in the Southern Hemisphere. The hemispheric asymmetry cannot be accounted for by the effect of the interplanetary magnetic field By component and the seasonal difference of ionospheric conductivities in the two hemispheres, which were used to interpret satellite-observed real-time auroral intensity asymmetries in the two hemispheres in previous studies. We suggest that the hemispheric asymmetry is the combined effect of the prenoon-postnoon variations of the magnetosheath density and local ionospheric conductivity.

Two types of whistler waves in the hall reconnection region

Whistler waves are believed to play an important role during magnetic reconnection. Here we report the near-simultaneous occurrence of two types of the whistler-mode waves in the magnetotail Hall r ...

Resonant scattering of central plasma sheet protons by multiband EMIC waves and resultant proton loss timescales

This is a companion study to Liang et al. (2014) which reported a “reversed” energy-latitude dispersion pattern of ion precipitation in that the lower energy ion precipitation extends to lower latitudes than the higher-energy ion precipitation. Electromagnetic ion cyclotron (EMIC) waves in the central plasma sheet (CPS) have been suggested to account for this reversed-type ion precipitation. To further investigate the association, we perform a comprehensive study of pitch angle diffusion rates induced by EMIC wave and the resultant proton loss timescales at L = 8–12 around the midnight. Comparing the proton scattering rates in the Earths dipole field and a more realistic quiet time geomagnetic field constructed from the Tsyganenko 2001 (T01) model, we find that use of a realistic, nondipolar magnetic field model not only decreases the minimum resonant energies of CPS protons but also considerably decreases the limit of strong diffusion and changes the proton pitch angle diffusion rates. Adoption of the T01 model increases EMIC wave diffusion rates at > ~ 60° equatorial pitch angles but decreases them at small equatorial pitch angles. Pitch angle scattering coefficients of 1–10 keV protons due to H+ band EMIC waves can exceed the strong diffusion rate for both geomagnetic field models. While He+ and O+ band EMIC waves can only scatter tens of keV protons efficiently to cause a fully filled loss cone at L > 10, in the T01 magnetic field they can also cause efficient scattering of ~ keV protons in the strong diffusion limit at L > 10. The resultant proton loss timescales by EMIC waves with a nominal amplitude of 0.2 nT vary from a few hours to several days, depending on the wave band and L shell. Overall, the results demonstrate that H+ band EMIC waves, once present, can act as a major contributor to the scattering loss of a few keV protons at lower L shells in the CPS, accounting for the reversed energy-latitude dispersion pattern of proton precipitation at low energies (~ keV) on the nightside. The pitch angle coverage for H+ band EMIC wave resonant scattering strongly depends on proton energy, L shell, and field model. He+ and O+ band EMIC waves tend to cause efficient scattering loss of protons at higher energies, thereby importantly contributing to the isotropic distribution of higher energy (> ~ 10 keV) protons at higher L shells on the nightside where the geomagnetic field line is highly stretched. Our results also suggest that scattering by H+ band EMIC waves may significantly contribute to the formation of the reversed-type CPS proton precipitation on the dawnside where both the wave activity and occurrence probability is statistically high.

Occurrence characteristics of outer zone relativistic electron butterfly distribution: A survey of Van Allen Probes REPT measurements

Using Van Allen Probes Relativistic Electron Proton Telescope (REPT) pitch angle resolved electron flux data from September 2012 to March 2015, we investigate in detail the global occurrence pattern of equatorial (|λ| ≤ 3°) butterfly distribution of outer zone relativistic electrons and its potential correlation with the solar wind dynamic pressure. The statistical results demonstrate that these butterfly distributions occur with the highest occurrence rate ~ 80% at ~ 20–04 magnetic local time (MLT) and L > ~ 5.5 and with the second peak (> ~ 50%) at ~ 11–15 MLT of lower L shells ~ 4.0. They can also extend to L = 3.5 and to other MLT intervals but with the occurrence rates predominantly < ~25%. It is further shown that outer zone relativistic electron butterfly distributions are likely to peak between 58° and 79° for L = 4.0 and 5.0 and between 37° and 58° for L = 6.0, regardless of the level of solar wind dynamic pressure. Relativistic electron butterfly distributions at L = 4.0 also exhibit a pronounced day-night asymmetry in response to the Pdyn variations. Compared to the significant L shell and MLT dependence of the global occurrence pattern, outer zone relativistic electron butterfly distributions show much less but still discernable sensitivity to Pdyn, geomagnetic activity level, and electron energy, the full understanding of which requires future attempts of detailed simulations that combine and differentiate underlying physical mechanisms of the geomagnetic field asymmetry and scattering by various magnetospheric waves.

Solar cycle variations of trapped proton flux in the inner radiation belt

Trapped proton population in the inner radiation belt is highly dense, posing a potential danger to astronauts and man-made space assets traversing through this region. While being significantly stable within timescales up to hundreds of days, inner zone proton fluxes can exhibit considerable solar cycle variations, which has not been investigated comprehensively yet. To analyze the long-term variation of the South Atlantic Anomaly (SAA), we adopt the proton flux data measured by NOAA 15 from 1999 through 2009 and perform statistical analyses on the basis of reasonable Gaussian fits. We report that the variation of the peak proton flux in the SAA is anticorrelated with that of F-10.7 during a solar cycle. There also exists a phase lag of 685 days between the solar F-10.7 flux and the proton flux. Similar features are seen for changes of the SAA distribution area, which in addition shows a rapid decrease during the solar maximum and a slow increase during the solar minimum. We also find that the region where the proton flux peaks drifts westward year by year with larger drift rates during the solar minimum. The peak region shifts southward during the solar maximum but in the opposite direction during the solar minimum with higher shift speed. Enhancements in solar wind dynamic pressure can favor the north-south drift of the SAA.