M. Wiltberger

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.

Magnetosphere Sawtooth Oscillations Induced by Ionospheric Outflow

Numerical simulations show that a magnetospheric disturbance is caused by an influx of O+ ions from the ionosphere. The sawtooth mode of convection of Earth’s magnetosphere is a 2- to 4-hour planetary-scale oscillation powered by the solar wind–magnetosphere–ionosphere (SW-M-I) interaction. Using global simulations of geospace, we have shown that ionospheric O+ outflows can generate sawtooth oscillations. As the outflowing ions fill the inner magnetosphere, their pressure distends the nightside magnetic field. When the outflow fluence exceeds a threshold, magnetic field tension cannot confine the accumulating fluid; an O+-rich plasmoid is ejected, and the field dipolarizes. Below the threshold, the magnetosphere undergoes quasi-steady convection. Repetition and the sawtooth period are controlled by the strength of the SW-M-I interaction, which regulates the outflow fluence.

Effects of causally driven cusp O + outflow on the storm time magnetosphere-ionosphere system using a multifluid global simulation

densities in the inner magnetosphere can increase the strength of the ring current, reducing Dst and inflating the magnetosphere. This effect is mostly found for the less energetic outflow specification. O + outflow is found to reduce the access of solar wind ions to the inner magnetosphere, which, through the MI coupling in LFM reduces the precipitating electron power, conductance and field‐aligned currents. The effect outflow has on the cross polar cap potential (CPCP) depends upon two competing factors. The reduction in Region I currents when outflow is present appears to increase the CPCP whilst the inflation of the magnetosphere due to an enhanced ring current decreases the CPCP. Citation: Brambles, O. J., W. Lotko, P. A. Damiano, B. Zhang, M. Wiltberger, and J. Lyon (2010), Effects of causally driven cusp O + outflow on the storm time magnetosphere‐ionosphere system using a multifluid global simulation, J. Geophys. Res., 115,

MHD/particle simulations of radiation belt dynamics

Abstract 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–Sun 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.

Geospace Environment Modeling 2008–2009 Challenge: Ground magnetic field perturbations

helps the users of the modeling products to better understand the capabilities of the models and to choose the approach that best suits their specific needs. Further, metrics!based analyses are important for addressing the differences between various modeling approaches and for measuring and guiding the progress in the field. In this paper, the metrics!based results of the ground magnetic field perturbation part of the Geospace Environment Modeling 2008‐2009 Challenge are reported. Predictions made by 14 different models, including an ensemble model, are compared to geomagnetic observatory recordings from 12 different northern hemispheric locations. Five different metrics are used to quantify the model performances for four storm events. It is shown that the ranking of the models is strongly dependent on the type of metric used to evaluate the model performance. None of the models rank near or at the top systematically for all used metrics. Consequently, one cannot pick the absolute“winner”: the choice for the best model depends on the characteristics of the signal one is interested in. Model performances vary also from event to event. This is particularly clear for root!mean!square difference and utility metric!based analyses. Further, analyses indicate that for some of the models, increasing the global magnetohydrodynamic model spatial resolution and the inclusion of the ring current dynamics improve the models’capability to generate more realistic ground magnetic field fluctuations.

Modeling studies of the impact of high-speed streams and co- rotating interaction regions on the thermosphere-ionosphere

Received 28 November 2011; revised 14 June 2012; accepted 14 June 2012; published 1 August 2012. [1] Changes in the thermosphere-ionosphere system caused by high-speed streams in the solar wind, and the co-rotating interaction regions they engender, are studied using a combination of model simulations and data analysis. The magnetospheric responses to these structures and consequent ionospheric drivers are simulated using the numerical Coupled Magnetosphere-Ionosphere-Thermosphere model and the empirical Weimer 2005 model, finding that the interplanetary magnetic field (IMF) is more important than solar wind speed and density per se in controlling magnetosphere-ionosphere coupling. The NCAR Thermosphere-Ionosphere-Electrodynamics General Circulation Model is then employed to calculate neutral density, nitric oxide cooling, and electron density, for comparison to space-based measurements from the STAR instrument on the CHAMP satellite, the SABER instrument on the TIMED satellite, and GPS occultations from the COSMIC mission, respectively. The recurrent, periodic changes observed under solar minimum conditions during 2008, and particularly during the Whole Heliospheric Interval (March–April of 2008), are simulated by the model and compared to these measurements. Numerical experiments were conducted to elucidate the mechanisms of solar wind and IMF forcing, setting the solar wind speed and density to nominal values, smoothing the IMF, and also setting it to zero. The results confirm the importance of IMF variations, particularly its north-south component (Bz), but also show that when the average Bz values are negative (southward), the interaction with increased solar wind speed amplifies the magnetosphere-ionosphere-thermosphere response. Conversely, during events when Bz is on average positive (northward), even large increases in solar wind speed have small effects on the system.

High-resolution global magnetohydrodynamic simulation of bursty bulk flows

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.

The dependence of the coupled magnetosphere-ionosphere- thermosphere system on the Earth's magnetic dipole moment

[1] The strength of the Earth’s magnetic field changes over time. We use simulations with the Coupled Magnetosphere-Ionosphere-Thermosphere model to investigate how the magnetosphere, upper atmosphere, and solar quiet (Sq) geomagnetic variation respond as the geomagnetic dipole moment M varies between 2⋅10 22 and 10⋅10 22 Am 2 . We find that the magnetopause stand-off distance and the cross-polar cap potential increase, while the polar cap size decreases, with increasing M. Their dependence on M is stronger than predicted by previous studies. We also show for the first time that the shape of the magnetosphere starts to change for M ≤ 4⋅10 22 Am 2 . This may be due to enhanced magnetopause erosion and/or to strong changes in the ionospheric conductance, which affect the field-aligned currents and the magnetic fields they create in the magnetosphere,

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 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.

The effects of seasonal and diurnal variations in the Earth's magnetic dipole orientation on solar wind–magnetosphere‐ionosphere coupling

Received 13 April 2012; revised 19 September 2012; accepted 10 October 2012; published 16 November 2012. [1] The angle m between the geomagnetic dipole axis and the geocentric solar magnetospheric (GSM) z axis, sometimes called the “dipole tilt,” varies as a function of UT and season. Observations have shown that the cross-polar cap potential tends to maximize near the equinoxes, when on average m = 0, with smaller values observed near the solstices. This is similar to the well-known semiannual variation in geomagnetic activity. We use numerical model simulations to investigate the role of two possible mechanisms that may be responsible for the influence of m on the magnetosphere-ionosphere system: variations in the coupling efficiency between the solar wind and the magnetosphere and variations in the ionospheric conductance over the polar caps. Under southward interplanetary magnetic field (IMF) conditions, variations in ionospheric conductance at high magnetic latitudes are responsible for 10–30% of the variations in the cross-polar cap potential associated with m, but variations in solar wind–magnetosphere coupling are more important and responsible for 70–90%. Variations in viscous processes contribute slightly to this, but variations in the reconnection rate with m are the dominant cause. The variation in the reconnection rate is primarily the result of a variation in the length of the section of the separator line along which relatively strong reconnection occurs. Changes in solar wind–magnetosphere coupling also affect the field-aligned currents, but these are influenced as well by variations in the conductance associated with variations in m, more so than the cross-polar cap potential. This may be the case for geomagnetic activity too.