Bibliography
|
Notice:
|
Found 11 entries in the Bibliography.
Showing entries from 1 through 11
| 2019 |
|
Simulations of Electron Energization and Injection by BBFs Using High-Resolution LFM MHD Fields We study electron injection and energization by bursty bulk flows (BBFs), by tracing electron trajectories using magnetohydrodynamic (MHD) field output from the Lyon-Fedder-Mobarry (LFM) code. The LFM MHD simulations were performed using idealized solar wind conditions to produce BBFs. We show that BBFs can inject energetic electrons of few to 100 keV from the magnetotatail beyond -24 RE to inward of geosynchronous, while accelerating them in the process. We also show the dependence of energization and injection on the initial relative position of the electrons to the magnetic field structure of the BBF, the initial pitch angle, and the initial energy. In addition, we have shown that the process can be nonadiabatic with violation of the first adiabatic invariant (μ). Further, we discuss the mechanism of energization and injection in order to give generalized insight into the process. Eshetu, W.; Lyon, J.; Hudson, M.; Wiltberger, M.; Published by: Journal of Geophysical Research: Space Physics Published on: 02/2019 YEAR: 2019 DOI: 10.1029/2018JA025789 |
| 2018 |
|
Pitch Angle Scattering of Energetic Electrons by BBFs Field line curvature scattering by the magnetic field structure associated with bursty bulk flows (BBFs) has been studied, using simulated output fields from the Lyon-Fedder-Mobarry global magnetohydrodynamic code for specified solar wind input. There are weak magnetic field strength (B) regions adjacent to BBFs observed in the simulations. We show that these regions can cause strong scattering where the first adiabatic invariant changes by several factors within one equatorial crossing of energetic electrons of a few kiloelectron volts when the BBFs are beyond 10RE geocentric in the tail. Scattering by BBFs decreases as they move toward the Earth or when the electron energy decreases. For radiation belt electrons near or inside geosynchronous orbit we demonstrate that the fields associated with BBFs can cause weak scattering where the fractional change of the first invariant (μ0) within one equatorial crossing is small, but the change due to several crossings can accumulate. For the weak scattering case we developed a method of calculating the pitch angle diffusion coefficient Dαα. Dαα for radiation belt electrons for one particular BBF were calculated as a function of initial energy, equatorial pitch angle, and radial location. These Dαα values were compared to calculated Dαα for a dipole field with no electric field. We further compared Dαα values with that of stretched magnetic fields calculated by Artemyev et al. (2013, https://doi.org/10.5194/angeo-31-1485-2013) at r≈7RE. Results show that scattering by BBFs can be comparable to the most highly stretched magnetic field they studied. Eshetu, W.; Lyon, J.; Hudson, M.; Wiltberger, M.; Published by: Journal of Geophysical Research: Space Physics Published on: 10/2018 YEAR: 2018 DOI: 10.1029/2018JA025788 |
| 2017 |
|
Simulated prompt acceleration of multi-MeV electrons by the 17 March 2015 interplanetary shock Prompt enhancement of relativistic electron flux at L = 3-5 has been reported from Van Allen Probes Relativistic Electron Proton Telescope (REPT) measurements associated with the 17 March 2015 interplanetary shock compression of the dayside magnetosphere. Acceleration by \~ 1 MeV is inferred on less than a drift time scale as seen in prior shock compression events, which launch a magetosonic azimuthal electric field impulse tailward. This impulse propagates from the dayside around the flanks accelerating electrons in drift resonance at the dusk flank. Such longitudinally localized acceleration events produce a drift echo signature which was seen at >1 MeV energy on both Van Allen Probe spacecraft, with sustained observations by Probe B outbound at L = 5 at 2100 MLT at the time of impulse arrival, measured by the Electric Fields and Waves instrument. MHD-test particle simulations are presented which reproduce drift echo features observed in the REPT measurements at Probe B, including the energy and pitch angle dependence of drift echoes observed. While the flux enhancement was short-lived for this event due to subsequent inward motion of the magnetopause, stronger events with larger electric field impulses, as observed in March 1991 and the Halloween 2003 storm, produce enhancements which can be quantified by the inward radial transport and energization determined by the induction electric field resulting from dayside compression. Hudson, Mary; Jaynes, Allison; Kress, Brian; Li, Zhao; Patel, Maulik; Shen, Xiaochen; Thaller, Scott; Wiltberger, Michael; Wygant, John; Published by: Journal of Geophysical Research: Space Physics Published on: 09/2017 YEAR: 2017 DOI: 10.1002/2017JA024445 17 March 2015; MeV electron acceleration; Radiation belt; test-particle simulation; Van Allen Probes |
|
The 17 March 2015 St. Patrick\textquoterights Day Storm is the largest geomagnetic storm to date of Solar Cycle 24, with a Dst of -223 nT. The magnetopause moved inside geosynchronous orbit under high solar wind dynamic pressure and strong southward IMF Bz causing loss, however a subsequent drop in pressure allowed for rapid rebuilding of the radiation belts. The 17 March 2013 storm also shows similar effects on outer zone electrons: first a rapid dropout due to inward motion of the magnetopause followed by rapid increase in flux above the pre-storm level early in the recovery phase and a slow increase over the next 12 days. These phases can be seen in temporal evolution of the electron phase space density measured by the ECT instruments on Van Allen Probes. Using the Lyon-Fedder-Mobarry global MHD model driven by upstream solar wind measurements, we simulated both St. Patrick\textquoterights Day 2013 and 2015 events, analyzing LFM electric and magnetic fields to calculate radial diffusion coefficients. These coefficients have been implemented in a radial diffusion code, using the measured electron phase space density following the local heating as the initial radial profile and outer boundary condition for subsequent temporal evolution over the next 12 days, beginning 18 March. Agreement with electron phase space density at 1000 MeV/G measured by the MagEIS component of the ECT instrument suite on Van Allen Probes was much improved using radial diffusion coefficients from the MHD simulations relative to coefficients parametrized by a global geomagnetic activity index. Li, Zhao; Hudson, Mary; Patel, Maulik; Wiltberger, Michael; Boyd, Alex; Turner, Drew; Published by: Journal of Geophysical Research: Space Physics Published on: 06/2017 YEAR: 2017 DOI: 10.1002/2016JA023846 March 2013; March 2015; radial diffusion; Radiation belt; Van Allen Probes |
| 2016 |
|
The 17\textendash18 March 2015 storm is the largest geomagnetic storm in the Van Allen Probes era to date. The Lyon-Fedder-Mobarry global MHD model has been run for this event using ARTEMIS data as solar wind input. The ULF wave power spectral density of the azimuthal electric field and compressional magnetic field is analyzed in the 0.5\textendash8.3 mHz range. The lowest three azimuthal modes account for 70\% of the total power during quiet times. However, during high activity, they are not exclusively dominant. The calculation of the radial diffusion coefficient is presented. We conclude that the electric field radial diffusion coefficient is dominant over the magnetic field coefficient by one to two orders of magnitude. This result contrasts with the dominant magnetic field diffusion coefficient used in most 3-D diffusion models. Li, Zhao; Hudson, Mary; Paral, Jan; Wiltberger, Michael; Turner, Drew; Published by: Journal of Geophysical Research: Space Physics Published on: 07/2016 YEAR: 2016 DOI: 10.1002/2016JA022508 March 2015; radial diffusion; radial diffusion coefficient; Radiation belt; ULF waves; Van Allen Probes |
| 2015 |
|
Magnetohydrodynamic modeling of three Van Allen Probes storms in 2012 and 2013 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\textendashFedder\textendashMobarry (LFM) global magnetohydrodynamic (MHD) magnetospheric simulation code coupled to the Rice Convection Model (RCM). Two other storms with strong magnetopause compression, 8\textendash9 October 2012 and 17\textendash18 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. Paral, J.; Hudson, M.; Kress, B.; Wiltberger, M.; Wygant, J.; Singer, H.; Published by: Annales Geophysicae Published on: 08/2015 YEAR: 2015 DOI: 10.5194/angeo-33-1037-2015 |
|
Modeling CME-shock driven storms in 2012 - 2013: MHD-test particle simulations The Van Allen Probes spacecraft have provided detailed observations of the energetic particles and fields environment for 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 pre-noon 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 (Hudson et al., 2014). Modeling magnetopause dropout events in 2013 with electric field diagnostics that were not available for storms immediately following launch has improved our understanding of the complex role that ULF waves play in radial transport during such events. Hudson, M.; Paral, J.; Kress, B.; Wiltberger, M.; Baker, D.; Foster, J.; Turner, D.; Wygant, J.; Published by: Journal of Geophysical Research: Space Physics Published on: 01/2015 YEAR: 2015 DOI: 10.1002/2014JA020833 |
| 2014 |
|
Simulated magnetopause losses and Van Allen Probe flux dropouts Three radiation belt flux dropout events seen by the Relativistic Electron Proton Telescope soon after launch of the Van Allen Probes in 2012 (Baker et al., 2013a) have been simulated using the Lyon-Fedder-Mobarry MHD code coupled to the Rice Convection Model, driven by measured upstream solar wind parameters. MHD results show inward motion of the magnetopause for each event, along with enhanced ULF wave power affecting radial transport. Test particle simulations of electron response on 8 October, prior to the strong flux enhancement on 9 October, provide evidence for loss due to magnetopause shadowing, both in energy and pitch angle dependence. Severe plasmapause erosion occurred during ~ 14 h of strongly southward interplanetary magnetic field Bz beginning 8 October coincident with the inner boundary of outer zone depletion. Hudson, M.; Baker, D.; Goldstein, J.; Kress, B.; Paral, J.; Toffoletto, F.; Wiltberger, M.; Published by: Geophysical Research Letters Published on: 02/2014 YEAR: 2014 DOI: 10.1002/2014GL059222 |
| 2012 |
|
Radiation belt 2D and 3D simulations for CIR-driven storms during Carrington Rotation 2068 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\textendashIonosphere\textendashThermosphere (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\textendash750 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\textendash30 storm, less than observed at GOES, and recovery to initial flux levels as seen for the April 4\textendash7 storm. Hudson, M.; Brito, Thiago; Elkington, Scot; Kress, Brian; Li, Zhao; Wiltberger, Mike; Published by: Journal of Atmospheric and Solar-Terrestrial Physics Published on: 07/2012 YEAR: 2012 DOI: 10.1016/j.jastp.2012.03.017 |
| 2006 |
|
In an MHD particle simulation of the September 1998 magnetic storm the evolution of the radiation belt electron radial flux profile appears to be diffusive, and diffusion caused by ULF waves has been invoked as the probable mechanism. In order to separate adiabatic and nonadiabatic effects and to investigate the radial diffusion mechanism during this storm, in this work we solve a radial diffusion equation with ULF wave diffusion coefficients and a time-dependent outer boundary condition, and the results are compared with the phase space density of the MHD particle simulation. The diffusion coefficients include contributions from both symmetric resonance modes (ω ≈ mωd, where ω is the wave frequency, m is the azimuthal wave number, and ωd is the bounce-averaged drift frequency) and asymmetric resonance modes (ω ≈ (m \textpm 1)ωd). ULF wave power spectral densities are obtained from a Fourier analysis of the electric and magnetic fields of the MHD simulation and are used in calculating the radial diffusion coefficients. The asymmetric diffusion coefficients are proportional to the magnetic field asymmetry, which is also calculated from the MHD field. The resulting diffusion coefficients vary with the radial coordinate L (the Roederer L-value) and with time during different phases of the storm. The last closed drift shell defines the location of the outer boundary. Both the location of the outer boundary and the value of the phase space density at the outer boundary are time-varying. The diffusion calculation simulates a 42-hour period during the 24\textendash26 September 1998 magnetic storm, starting just before the storm sudden commencement and ending in the late recovery phase. The differential flux calculated in the MHD particle simulation is converted to phase space density. Phase space densities in both simulations (diffusion and MHD particle) are functions of Roederer L-value for fixed first and second adiabatic invariants. The Roederer L-value is calculated using drift shell tracing in the MHD magnetic field, and particles have zero second invariant. The radial diffusion calculation reproduces the main features of the MHD particle simulation quite well. The symmetric resonance modes dominate the radial diffusion, especially in the inner and middle L region, while the asymmetric resonances are more important in the outer region. Using both symmetric and asymmetric terms gives a better result than using only one or the other and is better than using a simple power law diffusion coefficient. We find that it is important to specify the value of the phase space density on the outer boundary dynamically in order to get better agreement between the radial diffusion simulation and the MHD particle simulation. Fei, Yue; Chan, Anthony; Elkington, Scot; Wiltberger, Michael; Published by: Journal of Geophysical Research Published on: 12/2006 YEAR: 2006 DOI: 10.1029/2005JA011211 |
| 2002 |
|
MHD/particle simulations of radiation belt dynamics Particle fluxes in the outer radiation belts can show substantial variation in time, over scales ranging from a few minutes, such as during the sudden commencement phase of geomagnetic storms, to the years-long variations associated with the progression of the solar cycle. As the energetic particles comprising these belts can pose a hazard to human activity in space, considerable effort has gone into understanding both the source of these particles and the physics governing their dynamical behavior. Computationally tracking individual test particles in a model magnetosphere represents a very direct, physically-based approach to modeling storm-time radiation belt dynamics. Using global magnetohydrodynamic models of the Earth\textendashSun system coupled with test particle simulations of the radiation belts, we show through two examples that such simulations are capable of capturing the outer zone radiation belt configuration at a variety of time scales and through all phases of a geomagnetic storm. Such simulations provide a physically-based method of investigating the dynamics of the outer radiation zone, and hold promise as a viable method of providing global nowcasts of the radiation environment during geomagnetically active periods. ELKINGTON, S; Hudson, M.; Wiltberger, M.J; Lyon, J.; Published by: Journal of Atmospheric and Solar-Terrestrial Physics Published on: 04/2002 YEAR: 2002 DOI: 10.1016/S1364-6826(02)00018-4 Shock-Induced Transport. Slot Refilling and Formation of New Belts. |
1