A high‐resolution global magnetohydrodynamic simulation is conducted with the Lyon‐Fedder‐Mobarry (LFM) model for idealized solar wind conditions. Within the simulation results high‐speed flows are ...seen throughout the magnetotail when the interplanetary magnetic field (IMF) is southward. Case study analysis of these flows shows that they have an enhancement in BZ and a decrease in density preceding a peak in the flow velocity. A careful examination of the structure within the magnetotail shows that these features are driven by bursts of magnetic reconnection. In addition to the case study, a superposed epoch analysis of flows occurring during a 90 min interval of southward IMF yields statistical properties that are in qualitative agreement with observational analysis of bursty bulk flows (BBFs). For the idealized simulation conditions, the most significant differences with the observational results are a broader velocity profile in time, which becomes narrower away from the center of the current sheet, and a larger density drop after flow passage. The peak BZ amplitude is larger than in observations and precedes the peak in the flow velocity. We conclude that the LFM simulations are reproducing the statistical features of BBFs and that they are driven by spatially and temporally localized reconnection events within the simulation domain.
Key Points
High‐resolution LFM simulations contain BBFs
BBFs have statistical properties similar to observations
BBFs are generated by reconnection within simulation
Explosive magnetotail activity has long been understood in the context of its auroral manifestations. While global models have been used to interpret and understand many magnetospheric processes, the ...temporal and spatial scales of some auroral forms have been inaccessible to global modeling creating a gulf between observational and theoretical studies of these phenomena. We present here an important step toward bridging this gulf using a newly developed global magnetosphere‐ionosphere model with resolution capturing
≲ 30 km azimuthal scales in the auroral zone. In a global magnetohydrodynamic (MHD) simulation of the growth phase of a synthetic substorm, we find the self‐consistent formation and destabilization of localized magnetic field minima in the near‐Earth magnetotail. We demonstrate that this destabilization is due to ballooning‐interchange instability which drives earthward entropy bubbles with embedded magnetic fronts. Finally, we show that these bubbles create localized field‐aligned current structures that manifest in the ionosphere with properties matching observed auroral beads.
Plain Language Summary
The aurora has long been used as a window onto the magnetosphere. However, auroral observations are inherently limited in trying to reconstruct global magnetospheric dynamics from the “magnetic shadow” they cast on Earth. For this reason modeling has been used in tandem with observations to better contextualize and understand the data. Substorms, the violent reconfiguration of the magnetotail and one of the most dynamic magnetospheric phenomena, have been known to be preceded by the formation of bead‐like structures in the aurora. The processes responsible for auroral beading and their causal versus correlative role with substorm onset have remained an enduring mystery. The vast disparity between the spatial scales of auroral beads and those of the global magnetosphere has greatly complicated the use of modeling in unraveling this mystery. We show here for the first time a demonstration of the self‐consistent formation of a magnetospheric configuration that becomes unstable during the period preceding the substorm onset and that this instability manifests in the ionosphere with similar morphology to auroral beads. The global context of the model shows that the magnetospheric processes responsible for beading are not necessarily causal to onset but a consequence of the slow magnetotail reconfiguration that precedes onset.
Key Points
We present the first global magnetosphere simulation to reveal ballooning‐interchange instability of a narrow Bz minimum in the near‐Earth magnetotail
The instability is prominent during the substorm growth phase and generates earthward entropy bubbles with embedded magnetic fronts
The bubbles drive mesoscale ionospheric field‐aligned currents and auroral structures (beads) with properties matching to those observed
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.
Key Points
Simulations of electron flux 8 Oct 2012 produce dropout into L~5.8Magnetopause compression bounding the September storage ring formation modeledModel shows effects of magnetopause inward motion and ULF waves on electrons
We present results from global, three‐dimensional magnetohydrodynamic (MHD) simulations of the solar wind/magnetosphere interaction. These MHD simulations are used to study ultra low frequency (ULF) ...pulsations in the Earth's magnetosphere driven by shear instabilities at the flanks of the magnetopause. We drive the simulations with idealized, constant solar wind input parameters, ensuring that any discrete ULF pulsations generated in the simulation magnetosphere are not due to fluctuations in the solar wind. The simulations presented in this study are driven by purely southward interplanetary magnetic field (IMF) conditions, changing only the solar wind driving velocity while holding all of the other solar wind input parameters constant. We find surface waves near the dawn and dusk flank magnetopause and show that these waves are generated by the Kelvin‐Helmholtz (KH) instability. We also find that two KH modes are generated near the magnetopause boundary. One mode, the magnetopause KH mode, propagates tailward along the magnetopause boundary. The other mode, the inner KH mode, propagates tailward along the inner edge of the boundary layer (IEBL). We find large vortical structures associated with the inner KH mode that are centered on the IEBL. The phase velocities, wavelengths, and frequencies of the two KH modes are computed. The KH waves are found to be fairly monochromatic with well‐defined wavelengths. In addition, the inner and magnetopause KH modes are coupled and lead to a coupled oscillation of the low‐latitude boundary layer. The boundary layer thickness, d, is computed and we find maximum wave growth for kd = 0.5–1.0, where k is the wave number, consistent with the linear theory of the KH instability. We comment briefly on the effectiveness of these KH waves in the energization and transport of radiation belt electrons.
General methods for improving the specification of electron precipitation in global simulations are described and implemented in the Lyon‐Fedder‐Mobarry (LFM) global simulation model, and the quality ...of its predictions for precipitation is assessed. LFM's existing diffuse and monoenergetic electron precipitation models are improved, and new models are developed for lower energy, broadband, and direct‐entry cusp precipitation. The LFM simulation results for combined diffuse plus monoenergetic electron precipitation exhibit a quadratic increase in the hemispheric precipitation power as the intensity of solar wind driving increases, in contrast with the prediction from the OVATION Prime (OP) 2010 empirical precipitation model which increases linearly with driving intensity. Broadband precipitation power increases approximately linearly with driving intensity in both models. Comparisons of LFM and OP predictions with estimates of precipitating power derived from inversions of Polar satellite UVI images during a double substorm event (28–29 March 1998) show that the LFM peak precipitating power is >4× larger when using the improved precipitation model and most closely tracks the larger of three different inversion estimates. The OP prediction most closely tracks the double peaks in the intermediate inversion estimate, but it overestimates the precipitating power between the two substorms by a factor >2 relative to all other estimates. LFMs polar pattern of precipitating energy flux tracks that of OP for broadband precipitation exhibits good correlation with duskside region 1 currents for monoenergetic energy flux that OP misses and fails to produce sufficient diffuse precipitation power in the prenoon quadrant that is present in OP. The prenoon deficiency is most likely due to the absence of drift kinetic physics in the LFM simulation.
Key Points
Electron precipitation models are developed for global magnetosphere simulations
Monoenergetic and diffuse precipitation exhibit nonlinear relations with SW driving
Modeled precipitation power is consistent with estimations from UVI images
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.
Key Points
The mechanism of electron energization and injection is sensitive to energy, position, and pitch angle
Changes in μ due to chaotic motion can have a significant effect on the expected energy gain
The pitch angle distribution of injected electrons at the equator can be peaked at high or low pitch angles
Pitch Angle Scattering of Energetic Electrons by BBFs Eshetu, W. W.; Lyon, J. G.; Hudson, M. K. ...
Journal of geophysical research. Space physics,
November 2018, 2018-11-00, 20181101, Letnik:
123, Številka:
11
Journal Article
Recenzirano
Odprti dostop
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.
Key Points
Pitch angle scattering of energetic electrons by bursty bulk flow fields in MHD simulations with idealized solar wind input is demonstrated
A method of calculating the equatorial pitch angle diffusion coefficient for general electromagnetic fields has been developed
Pitch angle diffusion coefficients are calculated for radiation belt electrons in realistic BBF fields for the first time
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.
Key Points
MHD test particle simulations reproduce magnetopause loss seen by VA Probes
E‐field along satellite path agrees with measured shock impulse for 8 October 2013
Electrons are accelerated by prompt inward radial transport
We investigate the plasmaspheric influence on the resonant mode coupling of magnetospheric ultralow frequency (ULF) waves using the Lyon‐Fedder‐Mobarry (LFM) global magnetohydrodynamic (MHD) model. ...We present results from two different versions of the model, both driven by the same solar wind conditions: one version that contains a plasmasphere (the LFM coupled to the Rice Convection Model, where the Gallagher plasmasphere model is also included) and another that does not (the stand‐alone LFM). We find that the inclusion of a cold, dense plasmasphere has a significant impact on the nature of the simulated ULF waves. For example, the inclusion of a plasmasphere leads to a deeper (more earthward) penetration of the compressional (azimuthal) electric field fluctuations, due to a shift in the location of the wave turning points. Consequently, the locations where the compressional electric field oscillations resonantly couple their energy into local toroidal mode field line resonances also shift earthward. We also find, in both simulations, that higher‐frequency compressional (azimuthal) electric field oscillations penetrate deeper than lower frequency oscillations. In addition, the compressional wave mode structure in the simulations is consistent with a radial standing wave oscillation pattern, characteristic of a resonant waveguide. The incorporation of a plasmasphere into the LFM global MHD model represents an advance in the state of the art in regard to ULF wave modeling with such simulations. We offer a brief discussion of the implications for radiation belt modeling techniques that use the electric and magnetic field outputs from global MHD simulations to drive particle dynamics.
Key Points
Magnetosphere responds as a resonant waveguide to ULF fluctuations in solar wind dynamic pressure
Inclusion of a plasmasphere has a substantial impact on the nature of the simulated ULF waves
Inclusion of a plasmasphere leads to a deeper penetration of azimuthal electric field oscillations
PDF
