Search

During the second Balloon Array for Radiation Belt Relativistic Electron Losses (BARREL) campaign two solar energetic proton (SEP) events were observed. Although BARREL was designed to observe X-rays created during electron precipitation events, it is sensitive to X-rays from other sources. The gamma lines produced when energetic protons hit the upper atmosphere are used in this paper to study SEP events. During the second SEP event starting on 7 January 2014 and lasting ∼3 days, which also had a solar energetic electron (SEE) event occurring simultaneously, BARREL had six payloads afloat spanning all magnetic local time (MLT) sectors and L values. Three payloads were in a tight array (∼2 h in MLT and ∼2 ΔL) inside the inner magnetosphere and at times conjugate in both L and MLT with the Van Allen Probes (approximately once per day). The other three payloads mapped to higher L values with one payload on open field lines for the entire event, while the other two appear to be crossing from open to closed field lines. Using the observations of the SEE and SEP events, we are able to map the open-closed boundary. Halford et al. (2015) demonstrated how BARREL can monitor electron precipitation following an interplanetary shock created by a coronal mass ejection (ICME-shock) arrival at Earth, while in this study we look at the SEP event precursor to the arrival of the ICME-Shock in our cradle-to-grave view: from flare, to SEE and SEP events, to radiation belt electron precipitation.

Radiation belt protons in the kinetic energy range 24 to 76 MeV are being measured by the Relativistic Electron Proton Telescope on each of the two Van Allen Probes. Data have been processed for the purpose of studying variability in the trapped proton intensity during October 2013 to August 2015. For the lower energies (≲32 MeV), equatorial proton intensity near L = 2 showed a steady increase that is consistent with inward diffusion of trapped solar protons, as shown by positive radial gradients in phase space density at fixed values of the first two adiabatic invariants. It is postulated that these protons were trapped with enhanced efficiency during the 7 March 2012 solar proton event. A model that includes radial diffusion, along with known trapped proton source and loss processes, shows that the observed average rate of increase near L = 2 is predicted by the same model diffusion coefficient that is required to form the entire proton radiation belt, down to low L, over an extended (∼103 year) interval. A slower intensity decrease for lower energies near L = 1.5 may also be caused by inward diffusion, though it is faster than predicted by the model. Higher-energy (≳40 MeV) protons near the L = 1.5 intensity maximum are from cosmic ray albedo neutron decay. Their observed intensity is lower than expected by a factor ∼2, but the discrepancy is resolved by adding an unspecified loss process to the model with a mean lifetime ∼120 years.

The loss of protons in the outer part of the inner radiation belt (L = 2 to 3) during the 6 April 2000 solar energetic particles event has been investigated using test particle simulations that follow full Lorentz trajectories with both magnetic and electric fields calculated from an empirical model. The electric fields are calculated as inductive fields generated by the time-changing magnetic field, which is achieved by time stepping analytic magnetic fields. The simulation results are compared with proton measurements from the highly elliptical orbit satellite for three different energy ranges (8.5–35 MeV, 16–40 MeV, and 27–45 MeV) as well as previous modeling work done. In previous work, inner zone radiation belt loss during geomagnetic storms has been modeled by simulating field line curvature scattering in static magnetic field snapshots with no electric field. The inclusion of the inductive electric field causes an increase in loss to lower L shells, improving the agreement with the satellite data.

Coronal mass ejection (CME)-shock compression of the dayside magnetopause has been observed to cause both prompt enhancement of radiation belt electron flux due to inward radial transport of electrons conserving their first adiabatic invariant and prompt losses which at times entirely eliminate the outer zone. Recent numerical studies suggest that enhanced ultra-low frequency (ULF) wave activity is necessary to explain electron losses deeper inside the magnetosphere than magnetopause incursion following CME-shock arrival. A combination of radial transport and magnetopause shadowing can account for losses observed at radial distances into L = 4.5, well within the computed magnetopause location. We compare ULF wave power from the Electric Field and Waves (EFW) electric field instrument on the Van Allen Probes for the 8 October 2013 storm with ULF wave power simulated using the Lyon–Fedder–Mobarry (LFM) global magnetohydrodynamic (MHD) magnetospheric simulation code coupled to the Rice Convection Model (RCM). Two other storms with strong magnetopause compression, 8–9 October 2012 and 17–18 March 2013, are also examined. We show that the global MHD model captures the azimuthal magnetosonic impulse propagation speed and amplitude observed by the Van Allen Probes which is responsible for prompt acceleration at MeV energies reported for the 8 October 2013 storm. The simulation also captures the ULF wave power in the azimuthal component of the electric field, responsible for acceleration and radial transport of electrons, at frequencies comparable to the electron drift period. This electric field impulse has been shown to explain observations in related studies (Foster et al., 2015) of electron acceleration and drift phase bunching by the Energetic Particle, Composition, and Thermal Plasma Suite (ECT) instrument on the Van Allen Probes.

Access of solar and galactic cosmic rays to the Earth's magnetosphere is quantified in terms of geomagnetic cutoff rigidity. Numerically computed grids of cutoff rigidities are used to model cosmic ray flux in Earth's atmosphere and in low Earth orbit. In recent years, the development of more accurate dynamic geomagnetic field models and an increase in computer power have made a real-time data-driven geomagnetic cutoff computation extending over the inner magnetosphere possible. For computational efficiency, numerically computed cutoffs may be scaled to different altitudes and directions of arrival using the known analytic variation of cutoff in a pure dipole magnetic field. This paper is a presentation of numerical techniques developed to compute effective cutoff rigidities for space weather applications. Numerical tests to determine the error associated with scaling vertical cutoff rigidities with altitude in a realistic geomagnetic field model are included. The tests were performed to guide the development of spatial grids for modeling cosmic ray access to the inner magnetosphere and to gain a better understanding of the accuracy of numerically modeled cutoffs.

Balloon-borne instruments detecting radiation belt precipitation frequently observe oscillations in the millihertz frequency range. Balloons measuring electron precipitation near the poles in the 100 keV to 2.5 MeV energy range, including the MAXIS, MINIS, and most recently the Balloon Array for Relativistic Radiation belt Electron Losses balloon experiments, have observed this modulation at ULF wave frequencies. Although ULF waves in the magnetosphere are seldom directly linked to increases in electron precipitation since their oscillation periods are much larger than the gyroperiod and the bounce period of radiation belt electrons, test particle simulations show that this interaction is possible. Three-dimensional simulations of radiation belt electrons were performed to investigate the effect of ULF waves on precipitation. The simulations track the behavior of energetic electrons near the loss cone, using guiding center techniques, coupled with an MHD simulation of the magnetosphere, using the Lyon-Fedder-Mobarry code, during a coronal mass ejection (CME)-shock event on 17 March 2013. Results indicate that ULF modulation of precipitation occurs even without the presence of electromagnetic ion cyclotron waves, which are not resolved in the MHD simulation. The arrival of a strong CME-shock, such as the one simulated, disrupts the electric and magnetic fields in the magnetosphere and causes significant changes in both components of momentum, pitch angle, and L shell of radiation belt electrons, which may cause them to precipitate into the loss cone.

The outer radiation belt consists of relativistic (>0.5 MeV) electrons trapped on closed trajectories around Earth where the magnetic field is nearly dipolar. During increased geomagnetic activity, electron intensities in the belt can vary by orders of magnitude at different spatial and temporal scales. The main phase of geomagnetic storms often produces deep depletions of electron intensities over broad regions of the outer belt. Previous studies identified three possible processes that can contribute to the main‐phase depletions: adiabatic inflation of electron drift orbits caused by the ring current growth, electron loss into the atmosphere, and electron escape through the magnetopause boundary. In this paper we investigate the relative importance of the adiabatic effect and magnetopause loss to the rapid depletion of the outer belt observed at the Van Allen Probes spacecraft during the main phase of 17 March 2013 storm. The intensities of >1 MeV electrons were depleted by more than an order of magnitude over the entire radial extent of the belt in less than 6 h after the sudden storm commencement. For the analysis we used three‐dimensional test particle simulations of global evolution of the outer belt in the Tsyganenko‐Sitnov (TS07D) magnetic field model with an inductive electric field. Comparison of the simulation results with electron measurements from the Magnetic Electron Ion Spectrometer experiment shows that magnetopause loss accounts for most of the observed depletion at L>5, while at lower L shells the depletion is adiabatic. Both magnetopause loss and the adiabatic effect are controlled by the change in global configuration of the magnetic field due to storm time development of the ring current; a simulation of electron evolution without a ring current produces a much weaker depletion., Key Points Main‐phase depletions are caused by magnetopause lossesLosses are enabled by a diamagnetic effect due to storm time ring currentBoth the third and the second invariants are violated in the loss process

The Van Allen Probes spacecraft have provided detailed observations of the energetic particles and fields environment for coronal mass ejection (CME)-shock-driven storms in 2012 to 2013 which have now been modeled with MHD test particle simulations. The Van Allen Probes orbital plane longitude moved from the dawn sector in 2012 to near midnight and prenoon for equinoctial storms of 2013, providing particularly good measurements of the inductive electric field response to magnetopause compression for the 8 October 2013 CME-shock-driven storm. An abrupt decrease in the outer boundary of outer zone electrons coincided with inward motion of the magnetopause for both 17 March and 8 October 2013 storms, as was the case for storms shortly after launch. Modeling magnetopause dropout events in 2013 with electric field diagnostics that were not available for storms immediately following launch have improved our understanding of the complex role that ULF waves play in radial transport during such events.

Measurements of inner radiation belt protons have been made by the Van Allen Probes Relativistic Electron-Proton Telescopes as a function of kinetic energy (24 to 76 MeV), equatorial pitch angle, and magnetic L shell, during late 2013 and early 2014. A probabilistic data analysis method reduces background from contamination by higher-energy protons. Resulting proton intensities are compared to predictions of a theoretical radiation belt model. Then trapped protons originating both from cosmic ray albedo neutron decay (CRAND) and from trapping of solar protons are evident in the measured distributions. An observed double-peaked distribution in L is attributed, based on the model comparison, to a gap in the occurrence of solar proton events during the 2007 to 2011 solar minimum. Equatorial pitch angle distributions show that trapped solar protons are confined near the magnetic equator but that CRAND protons can reach low altitudes. Narrow pitch angle distributions near the outer edge of the inner belt are characteristic of proton trapping limits.

[1]Relativistic electron intensities in Earth's outer radiation belt can vary by multiple orders of magnitude on the time scales ranging from minutes to days. One fundamental process contributing to dynamic variability of radiation belt intensities is the radial transport of relativistic electrons across their drift shells. In this paper we analyze the properties of three-dimensional radial transport in a global magnetic field model driven by variations in the solar wind dynamic pressure. We use a test particle approach which captures anomalous effects such as drift orbit bifurcations. We show that the bifurcations lead to an order of magnitude increase in radial transport rates and enhance the energization at large equatorial pitch angles. Even at quiet time fluctuations in dynamic pressure, radial transport at large pitch angles exhibits strong deviations from the diffusion approximation. The radial transport rates are much lower at small pitch angle values which results in a better agreement with the diffusion approximation.

[1] Observations of geomagnetically trapped 27–45 MeV protons following the November 2003 magnetic storm show a gradual intensity rise that is interpreted as a direct measurement of the cosmic ray albedo neutron decay (CRAND) source strength. The intensity rise is simulated by combining the detector response function with a model CRAND source, obtained by drift-averaging neutron intensity from Monte Carlo simulation of cosmic ray interactions in the atmosphere. The simulation, for 2.4

Solar energetic particle (SEP) cutoffs at geosynchronous orbit are sensitive to moderate geomagnetic activity and undergo daily variations due to the day–night asymmetry of the magnetosphere. At geosynchronous orbit, cutoff rigidity also has a large directional dependence, with the highest cutoff rigidity corresponding to ions arriving from magnetic east and lowest cutoff rigidity corresponding to ions incident from the west. Consequently, during geomagnetically quiet periods, the SEP flux observed by an eastward facing particle detector is significantly lower than observed by a westward facing particle detector. During geomagnetically disturbed periods the cutoff is suppressed allowing SEPs access well inside of geosynchronous, so that the east–west SEP flux ratio approaches unity. Variations in the east–west SEP flux ratio observed by GOES Energetic Particle Sensors (EPS) have recently been reported by Rodriguez et al. (2010) . In NOAA’s operational processing of EPS count rates into differential fluxes, the differential flux is treated as isotropic and flat over the energy width of the channel. To compare modeled SEP flux with GOES EPS observations, the anisotropy of the flux over the EPS energy range and field of view must be taken into account. A technique for making direct comparisons between GOES EPS observations and SEP flux modeled using numerically computed geomagnetic cutoffs is presented. Initial results from a comparison between modeled and observed flux during the 6–11 December 2006 SEP event are also presented. The modeled cutoffs reproduce the observed flux variations well but are in general too high.

The loss of protons in the outer part of the inner radiation belt (L = 2 to 3) during the 17 March 2015 geomagnetic storm was investigated using test particle simulations that follow full Lorentz trajectories with both magnetic and electric fields calculated from an empirical model. The simulation results presented here are compared with proton pitch angle measurements from the Van Allen Probe satellites Relativistic Electron Proton Telescope (REPT) instrument before and after the coronal mass ejection-shock-driven storm of 17–18 March 2015, with minimum Dst =− 223 nT, the strongest storm of Solar Cycle 24, for four different energy ranges with 30, 38, 50, and 66 MeV mean energies. Two simulations have been run, one with an inductive electric field and one without. All four energy channels show good agreement with the Van Allen Probes REPT measurements for low L (L 2.4. A previous study using the Highly Elliptical Orbiter 3 spacecraft also showed improved agreement when including the inductive electric field but was unable to compare effects on the pitch angle distributions.

As part of the International Heliospheric Year, the Whole Heliosphere Interval, Carrington Rotation 2068, from March 20 to April 16, 2008 was chosen as an internationally coordinated observing and modeling campaign. A pair of solar wind structures identified as Corotating Interaction Regions (CIR), characteristic of the declining phase of the solar cycle and solar minimum, was identified in solar wind plasma measurements from the ACE satellite. Such structures have previously been determined to be geoeffective in producing enhanced outer zone radiation belt electron fluxes, on average greater than at solar maximum. MHD fields from the Coupled Magnetosphere–Ionosphere–Thermosphere (CMIT) model driven by ACE solar wind measurements at L1 have been used to drive both 2D and 3D weighted test particle simulations of electron dynamics for the CIR subset of the month-long CMIT fields. Dropout in electron flux at geosynchronous orbit and enhancement during recovery phase, characteristic of CIR-driven storms, is seen in these moderate (Dstmin=−56, −33 nT) events, while the two CIRs were characterized by increased solar wind velocity in the 650–750 km/s range. The first beginning March 26 produced a greater enhancement in IMF Bz southward and stronger magnetospheric convection, leading to a greater radiation belt electron response at GOES. This study provides the first comparison of 2D and 3D particle dynamics in MHD simulation fields, incorporating the additional diffusive feature of Shebansky orbit trapping of electrons in the magnetic minima on the dayside above and below the equatorial plane. Overall loss occurs during the main phase for 2D and 3D simulations, while incorporation of plasmasheet injection in 2D runs produces a moderate enhancement for the March 26–30 storm, less than observed at GOES, and recovery to initial flux levels as seen for the April 4–7 storm.

Test-particle trajectories are computed in fields from a global MHD magnetospheric model simulation of the 29 October 2003 Storm Commencement to investigate trapping and transport of solar energetic electrons (SEEs) in the magnetosphere during severe storms. SEEs are found to provide a source population for a newly formed belt of > 10 MeV electrons in the Earth's inner zone radiation belts, which was observed following the 29 October 2003 storm. Energy and pitch angle distributions of the new belt are compared with results previously obtained [Kress, B.T., Hudson, M.K., Looper, M.D., Albert, J., Lyon, J.G., Goodrich, C.C., 2007. Global MHD test particle simulations of > 10 MeV radiation belt electrons during storm sudden commencement. Journal of Geophysical Research 112, A09215, doi:10.1029/2006JA012218], where outer belt electrons were used as a source for the new belt.

Discovery of the Van Allen radiation belts by instrumentation flown on Explorer 1 in 1958 was the first major discovery of the Space Age. A view of the belts as distinct inner and outer zones of energetic particles with different sources was modified by observations made during the Cycle 22 maximum in solar activity in 1989–1991, the first approaching the activity level of the International Geophysical Year of 1957–1958. The dynamic variability of outer zone electrons was measured by the NASA–Air Force Combined Radiation Release and Effects Satellite launched in July 1990. This variability is caused by distinct types of heliospheric structure which vary with the solar cycle. The largest fluxes averaged over a solar rotation occur during the declining phase from solar maximum, when high-speed streams and co-rotating interaction regions (CIRs) dominate the inner heliosphere, leading to recurrent storms. Intense episodic events driven by high-speed interplanetary shocks launched by coronal mass ejections (CMEs) prevail around solar maximum when CMEs occur most frequently. Only about half of moderate storms, defined by intensity of the ring current, lead to an overall flux increase, emphasizing the need to quantify loss as well as source processes; both increase when the magnetosphere is strongly driven. Three distinct types of acceleration are described in this review: prompt and diffusive radial transport, which increases energy while conserving the first invariant, and local acceleration by waves, which change the first invariant. The latter also produce pitch angle diffusion and loss, as does outward radial transport, especially when the magnetosphere is compressed. The effect of a dynamic magnetosphere boundary on radiation belt electrons is described in the context of MHD-test particle simulations driven by measured solar wind input.

The prompt trapping of solar energetic particles (SEPs) in the inner magnetosphere around L = 2 –2.5, including protons and heavier ions, has been observed at both the Cycle 22 and 23 solar maxima, in association with high-speed interplanetary shocks and storm sudden commencements (SSCs). Recent observations include the Bastille Day 2000 CME-driven storm as well as two in November 2001, which produced a long-lived new proton belt, as well as trapping of heavy ions up to Fe in all three cases. A survey of such events around the most recent solar maximum, including high altitude measurements from Polar and HEO satellites along with low altitude measurements from `the Solar Anomalous and Magnetospheric Particle Explorer (SAMPEX)', indicates similarities to the well-studied March 24, 1991 SSC event. In this event, electrons and protons in drift resonance with a magnetosonic impulse were transported radially inward. A requirement for such shock-induced acceleration is a high-speed CME-shock at 1 AU, which launches a perturbation with comparable velocity inside the magnetosphere. Secondly, there must be a source population which is drift-resonant with the impulse. The CME-shock itself is a source of solar energetic particles, both protons and heavy ions, with higher fluxes and harder spectra associated with faster moving CMEs. Arrival of the interplanetary shock compresses and changes the magnetosphere topology, leading to a reduction of the geomagnetic cutoff, initially around L = 4 for SEP protons. This effect is modeled using a 3D Lorentz integration of SEP trajectories in electric and magnetic fields taken from the Lyon–Fedder–Mobarry (LFM) global MHD code, with solar wind input parameters taken from spacecraft measurements upstream from the bow shock, carried out for the November 24, 2001 SEP event. The results indicate that an enhancement in solar wind dynamic pressure for this event plays a role in the observed injection of ions to low L-values, to form a new proton belt which has lasted for more than 2 years.

[1] Radiation levels in Earth's outer electron belt (L ≳ 2.5) vary by orders of magnitude on the time scales ranging from minutes to days. Multiple acceleration and loss processes operate across the belt and compete in defining its global variability. One such process is the drift orbit bifurcation effect. Caused by coupling of the drift and bounce motions, it breaks the second adiabatic invariant of radiation belt electrons producing their transport in radius and pitch angle. In this paper we investigate implications of drift orbit bifurcations to the global state and variability of the outer electron belt. For this purpose we use three-dimensional test particle simulations of electron guiding center motion in a realistic magnetic field model. We show that even at most quiet solar wind conditions bifurcations affect a broad range of the belt penetrating inside the geosynchronous orbit. This has an important practical implication for the analysis of experimental particle data: since the third adiabatic invariant is undefined for bifurcating orbit, the electron phase space density cannot be expressed in terms of three adiabatic invariants. We show that long-term transport of electrons due to drift orbit bifurcations is a complex combination of large ballistic jumps and small-amplitude diffusion in the second invariant and radial location. To model long-term transport, we derive an empirical map of the second invariant and radial jumps at bifurcations. The map can also be implemented by other radiation belt models, which cannot directly account for the physics of drift orbit bifurcations. Drift orbit bifurcations can produce electron losses through the magnetopause escape and through pitch angle scattering into the atmospheric loss cone. Most electrons, however, can stay quasi-trapped in the bifurcation regions for very long time periods. The pitch angle and radial transport due to drift orbit bifurcations lead to their meandering back and forth across the region producing mixing and recirculation of particle populations with different initial conditions. We show that this recirculation can greatly amplify electron energization by radial diffusion. Compared to the diffusion alone, the combined action of radial diffusion and drift orbit bifurcations can double electron energization at each recirculation cycle. Our results suggest that drift orbit bifurcations can play an important role in the buildup of increased electron fluxes in the storm recovery phase.

[1] A survey of 27 to 45 MeV proton measurements from the HEO-3 satellite during the years 1998 through 2005 has been taken to describe variability in the outer part of the inner radiation belt and slot region (L = 2 to 3). Rapid (∼1-day) changes are described as injection or loss events, characterized respectively by Gaussian or exponential L dependencies. The radial extent of both event types is correlated to the minimum Dst of associated magnetic storms, while the injection magnitude is correlated to the flux of associated interplanetary solar proton events. Changes in the maximal L of observed trapped protons are consistent with trapping limits estimated from magnetic field line curvature. The inward extent and energy independence of the observed loss events are inconsistent with field line curvature induced scattering in a static magnetic field. However, time-dependent geomagnetic cutoff suppression, observed during magnetic storms, may be the cause of significant losses. Drift resonance with electric field impulses caused by rapid magnetospheric compression is the likely cause of both solar proton injections and radial shifts of preexisting trapped protons.

[1] Prior to 2003, there are two known cases where ultrarelativistic (≳10 MeV) electrons appeared in the Earth's inner zone radiation belts in association with high speed interplanetary shocks: the 24 March 1991 and the less well studied 21 February 1994 storms. During the March 1991 event electrons were injected well into the inner zone on a timescale of minutes, producing a new stably trapped radiation belt population that persisted for ∼10 years. More recently, at the end of solar cycle 23, a number of violent geomagnetic disturbances resulted in large variations in ultrarelativistic electrons in the inner zone, indicating that these events are less rare than previously thought. Here we present results from a numerical study of shock-induced transport and energization of outer zone electrons in the 1–7 MeV range, resulting in a newly formed 10–20 MeV electron belt near L ∼ 3. Test particle trajectories are followed in time-dependent fields from an MHD magnetospheric model simulation of the 29 October 2003 storm sudden commencement (SSC) driven by solar wind parameters measured at ACE. The newly formed belt is predominantly equatorially mirroring. This result is in part due to an SSC electric field pulse that is strongly peaked in the equatorial plane, preferentially accelerating equatorially mirroring particles. The timescale for subsequent pitch angle diffusion of the new belt, calculated using quasi-linear bounce-averaged diffusion coefficients, is in agreement with the observed delay in the appearance of peak fluxes at SAMPEX in low Earth orbit. We also present techniques for modeling radiation belt dynamics using test particle trajectories in MHD fields. Simulations are performed using code developed by the Center for Integrated Space Weather Modeling.