Search

Two of the geomagnetic storms for the Space Weather Prediction Center (SWPC) Geospace Environment Modeling (GEM) challenge [cf. Pulkkinen et al., 2013] occurred after data were first acquired by the Active Magnetosphere and Planetary Electrodynamics Response Experiment (AMPERE). We compare Birkeland currents from AMPERE with predictions from four models for the 4-5 April 2010 and 5-6 August 2011 storms. The four models are: the Weimer [2005b] field-aligned current statistical model; the Lyon-Fedder-Mobarry magnetohydrodynamic (MHD) simulation; the Open Global Geospace Circulation Model MHD simulation; and the Space Weather Modeling Framework MHD simulation. The MHD simulations were run as described in Pulkkinen et al. [2013] and the results obtained from the Community Coordinated Modeling Center (CCMC). The total radial Birkeland current, ITotal, and the distribution of radial current density, Jr, for all models are compared with AMPERE results. While the total currents are well correlated, the quantitative agreement varies considerably. The Jr distributions reveal discrepancies between the models and observations related to the latitude distribution, morphologies, and lack of nightside current systems in the models. The results motivate enhancing the simulations first by increasing the simulation resolution, and then by examining the relative merits of implementing more sophisticated ionospheric conductance models, including ionospheric outflows or other omitted physical processes. Some aspects of the system, including substorm timing and location, may remain challenging to simulate, implying a continuing need for real-time specification.

Recently Gordeev et al. [2015] suggested a method to test global MHD models against statistical empirical data. They showed that four community-available global MHD models supported by the Community Coordinated Modeling Center (CCMC) produce a reasonable agreement with reality for those key parameters (the magnetospheric size, magnetic field and pressure) that are directly related to the large-scale equilibria in the outer magnetosphere. Based on the same set of simulation runs, here we investigate how the models reproduce the global loading-unloading cycle. We found that in terms of global magnetic flux transport, three examined CCMC models display systematically different response to idealized 2h north then 2h south IMFBz variation. The LFM model shows a depressed return convection and high loading rate during the growth phase as well as enhanced return convection and high unloading rate during the expansion phase, with the amount of loaded/unloaded magnetotail flux and the growth phase duration being the closest to their observed empirical values during isolated substorms. Two other models exhibit drastically different behavior. In the BATS-R-US model the plasma sheet convection shows a smooth transition to the steady convection regime after the IMF southward turning. In the Open GGCM a weak plasma sheet convection has comparable intensities during both the growth phase and the following slow unloading phase. We also demonstrate potential technical problem in the publicly-available simulations which is related to post-processing interpolation and could affect the accuracy of magnetic field tracing and of other related procedures.

There has long been interest in knowing the shape and location of the Earth's magnetopause and of the standing fast-mode bow shock upstream of the Earth's magnetosphere. This quest for knowledge spans both the research and operations arenas. Pertinent to the latter, nowcasting and near-term forecasting are important for determining the extent to which the magnetosphere is compressed or expanded due to the influence of the solar wind bulk plasma and fields and the coupling to other magnetosphere-ionosphere processes with possible effects on assets. This article provides an update to a previous article on the same topic published 15 years earlier, with focus on studies that have been conducted, the current status of nowcasting and forecasting of geophysical boundaries, and future endeavors.

To address challenges of assessing space weather modeling capabilities, the CommunityCoordinated Modeling Center is leading a newly establishedInternational Forum for Space WeatherModeling Capabilities Assessment. This paper presents preliminary results of validation of modeled foF2 (F2 layer critical frequency) and TEC (total electron content) during the first selected 2013 March storm event (17 March 2013). In this study, we used eight ionospheric models ranging from empirical to physics-based, coupled ionosphere-thermosphere and data assimilation models. The quantities we considered are TEC and foF2 changes and percentage changes compared to quiet time background, and the maximum and minimum percentage changes. In addition, we considered normalized percentage changes of TEC. We compared the modeled quantities with ground-based observations of vertical Global Navigation SatelliteSystem TEC (provided by Massachusetts Institute of Technology Haystack Observatory) and foF2 data (provided by Global Ionospheric Radio Observatory) at the 12 locations selected in middle latitudes of the American and European-African longitude sectors. To quantitatively evaluate the models’ performance, we calculated skill scores including correlation coefficient, root-mean square error (RMSE), ratio of the modeled to observed maximum percentage changes (yield), and timing error. Our study indicates that average RMSEs of foF2range from about 1 MHz to 1.5 MHz. The average RMSEs of TEC are between ~5 and ~10 TECU (1 TEC Unit= 1016el/m2). dfoF2[%] RMSEs are between 15% and 25%, which is smaller than RMSE of dTEC[%] ranging from30% to 60%. The performance of the models varies with the location and metrics considered.

On 21 January 2005, a moderate magnetic storm produced a number of anomalous features, some seen more typically during superstorms. The aim of this study is to establish the differences in the space environment from what we expect (and normally observe) for a storm of this intensity, which make it behave in some ways like a superstorm. The storm was driven by one of the fastest interplanetary coronal mass ejections in solar cycle 23, containing a piece of the dense erupting solar filament material. The momentum of the massive solar filament caused it to push its way through the flux rope as the interplanetary coronal mass ejection decelerated moving toward 1 AU creating the appearance of an eroded flux rope (see companion paper by Manchester et al. (2014)) and, in this case, limiting the intensity of the resulting geomagnetic storm. On impact, the solar filament further disrupted the partial ring current shielding in existence at the time, creating a brief superfountain in the equatorial ionosphere—an unusual occurrence for a moderate storm. Within 1 h after impact, a cold dense plasma sheet (CDPS) formed out of the filament material. As the interplanetary magnetic field (IMF) rotated from obliquely to more purely northward, the magnetotail transformed from an open to a closed configuration and the CDPS evolved from warmer to cooler temperatures. Plasma sheet densities reached tens per cubic centimeter along the flanks—high enough to inflate the magnetotail in the simulation under northward IMF conditions despite the cool temperatures. Observational evidence for this stretching was provided by a corresponding expansion and intensification of both the auroral oval and ring current precipitation zones linked to magnetotail stretching by field line curvature scattering. Strong Joule heating in the cusps, a by-product of the CDPS formation process, contributed to an equatorward neutral wind surge that reached low latitudes within 1–2 h and intensified the equatorial ionization anomaly. Understanding the geospace consequences of extremes in density and pressure is important because some of the largest and most damaging space weather events ever observed contained similar intervals of dense solar material.

Dynamics of the dayside magnetosphere and proton radiation belt was analyzed during unusual magnetic storm on 21 January 2005. We have found that during the storm from 1712 to 2400 UT, the subsolar magnetopause was continuously located inside geosynchronous orbit due to strong compression. The compression was found to be extremely strong from 1846 to 2035 UT when the dense plasma of fast erupting filament produced the solar wind dynamic pressure Pd peaked up to >100 nPa and, in the first time, the upstream solar wind was observed at geosynchronous orbit during almost 2 hours. Under the extreme compression, the outer magnetosphere at L > 5 was pushed inward and the outer radiation belt particles with energies of several tens of keV moved earthward, became adiabatically accelerated and accumulated in the inner magnetosphere at L 20%, which is well appropriate for erupting filaments and which is in agreement with the upper 27% threshold for the He/H ratio obtained from Cluster measurements.

We examine the uncertainties in two plasma parameters from their true values in a simulated asymmetric corona. We use the Corona Heliosphere (CORHEL) and Magnetohydrodynamics Around the Sphere (MAS) models in the Community Coordinated Modeling Center (CCMC) to investigate the differences between an assumed symmetric corona and a more realistic, asymmetric one. We were able to predict the electron temperatures and electron bulk flow speeds to within ± 0.5 MK and ± 100 km s−1, respectively, over coronal heights up to 5.0 R⊙ from Sun center. We believe that this technique could be incorporated in next-generation white-light coronagraphs to determine these electron plasma parameters in the low solar corona. We have conducted experiments in the past during total solar eclipses to measure the thermal electron temperature and the electron bulk flow speed in the radial direction in the low solar corona. These measurements were made at different altitudes and latitudes in the low solar corona by measuring the shape of the K-coronal spectra between 350 nm and 450 nm and two brightness ratios through filters centered at 385.0 nm/410.0 nm and 398.7 nm/423.3 nm with a bandwidth of ≈ 4 nm. Based on symmetric coronal models used for these measurements, the two measured plasma parameters were expected to represent those values at the points where the lines of sight intersected the plane of the solar limb.

One of the major types of solar eruption, coronal mass ejections (CMEs) not only impact space weather, but also can have significant societal consequences. CMEs cause intense geomagnetic storms and drive fast mode shocks that accelerate charged particles, potentially resulting in enhanced radiation levels both in ions and electrons. Human and technological assets in space can be endangered as a result. CMEs are also the major contributor to generating large amplitude Geomagnetically Induced Currents (GICs), which are a source of concern for power grid safety. Due to their space weather significance, forecasting the evolution and impacts of CMEs has become a much desired capability for space weather operations worldwide. Based on our operational experience at Space Weather Research Center at NASA Goddard Space Flight Center (http://swrc.gsfc.nasa.gov), we present here some of the insights gained about accurately predicting CME impacts, particularly in relation to space weather operations. These include: 1. The need to maximize information to get an accurate handle of three-dimensional (3-D) CME kinetic parameters and therefore improve CME forecast; 2. The potential use of CME simulation results for qualitative prediction of regions of space where solar energetic particles (SEPs) may be found; 3. The need to include all CMEs occurring within a ~24 h period for a better representation of the CME interactions; 4. Various other important parameters in forecasting CME evolution in interplanetary space, with special emphasis on the CME propagation direction. It is noted that a future direction for our CME forecasting is to employ the ensemble modeling approach.

We employ Magnetohydrodynamic (MHD) simulations to examine the creation and evolution of plasma vortices within the Earth's magnetosphere for steady solar wind plasma conditions. Very few vortices form during intervals of such solar wind conditions. Those that do remain in fixed positions for long periods (often hours) and exhibit rotation axes that point primarily in the x or y direction, parallel (or antiparallel) to the local magnetospheric magnetic field direction. Occasionally, the orientation of the axes rotates from the x direction to another direction. We compare our results with simulations previously done for unsteady solar wind conditions. By contrast, these vortices that form during intervals of varying solar wind conditions exhibit durations ranging from seconds (in the case of those with axes in the x or y direction) to minutes (in the case of those with axes in the z direction) and convect antisunward. The local-time dependent sense of rotation seen in these previously reported vortices suggests an interpretation in terms of the Kelvin–Helmholtz instability. For steady conditions, the biggest vortices developed on the dayside (about 6 RE in diameter), had their rotation axes aligned with the y direction and had the longest periods of duration. We attribute these vortices to the flows set up by reconnection on the high-latitude magnetopause during intervals of northward Interplanetary Magnetic Field (IMF) orientation. This is the first time that vortices due to high-latitude reconnection have been visualized. The model also successfully predicts the principal characteristics of previously reported plasma vortices within the magnetosphere, namely their dimension, flow velocities, and durations.

[1] We use the global magnetohydrodynamic (MHD) code BATS-R-US to model multipoint observations of Flux Transfer Event (FTE) signatures. Simulations with high spatial and temporal resolution predict that cavities of weak magnetic field strength protruding into the magnetosphere trail FTEs. These predictions are consistent with recently reported multi-point Cluster observations of traveling magnetopause erosion regions (TMERs).

[1] One of the ways to address the science needs of the research community and to enable science progress is to provide community access to modern space science models. The Community Coordinated Modeling Center (CCMC) is a multiagency partnership based at the Goddard Space Flight Center that hosts a set of state-of-the-art space science models ranging from the solar atmosphere to the Earth's upper atmosphere. The CCMC provides a Web-based, no-cost, Runs on Request system, by which the interested scientist can readily request simulations for time intervals of interest. CCMC also provides a tailored Web-based visualization interface for the model output, including near-real-time results from select models. Model outputs have been specifically tailored for easy comparison with observational data to facilitate data analysis and model validation. This paper provides an overview of CCMC activities, with an emphasis on the ionosphere-thermosphere models residing there.

[1] The effect of magnetospheric substorms on the ring current is not completely understood. Using a combination of the University of Michigan's BAT-S-RUS Model and Fok Ring Current Model, we modeled the effects of multiple substorms on the ring current by modeling multiple dipolarizations in the tail. Increasing the number of substorms corresponds to increases in the number of injections into the ring current. The ionospheric potential increases during periods of southward IMF. Energy increases are more dependent on the duration of large ionospheric potential than the number of substorm dipolarizations.

[1] The physics of steady driven magnetic reconnection at Earth's subsolar magnetopause is addressed. Three-dimensional, global magnetohydrodynamics (MHD) simulations of the magnetopause are compared with analytical solutions of the resistive MHD equations [Sonnerup and Priest, 1975] corresponding to magnetic field annihilation driven by an incompressible stagnation point flow. The simulations demonstrate that under steady southward interplanetary magnetic field conditions and when the plasma resistivity is spatially uniform, subsolar magnetopause reconnection occurs in long, thin Sweet-Parker current sheets via a flux pileup mechanism [Sonnerup and Priest, 1975; Priest and Forbes, 1986] (rather than in Petschek slow shock configurations). Magnetic energy piles up upstream of the magnetopause current sheet to accommodate the sub-Alfvenic solar wind inflow. The scaling of the pileup with Lundquist number, S, is consistent (approximately ∝ S1/4) with that predicted by the analytical, incompressible stagnation point flow solutions (though there are small corrections due to plasma compressibility in the simulations). Since there is a finite energy in the magnetosheath available to drive the magnetic pileup (and associated rapid magnetic reconnection), we expect the pileup to saturate and the reconnection rate to drop as the upstream plasma pressure drops to accommodate the pileup. Thus we expect the reconnection to stall, the rate vanishing in the limit S → ∞. We discuss the role of Hall electric fields in allowing the magnetic pileup to saturate before the reconnection begins to stall, permitting Alfvenic reconnection to occur in thin current sheets in the limit S → ∞.

Researching and forecasting the ever changing space environment (often referred to as space weather) and its influence on humans and their activities are model-intensive disciplines. This is true because the physical processes involved are complex, but, in contrast to terrestrial weather, the supporting observations are typically sparse. Models play a vital role in establishing a physically meaningful context for interpreting limited observations, testing theory, and producing both nowcasts and forecasts. For example, with accurate forecasting of hazardous space weather conditions, spacecraft operators can place sensitive systems in safe modes, and power utilities can protect critical network components from damage caused by large currents induced in transmission lines by geomagnetic storms.