skip to main content


The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 11:00 PM ET on Thursday, May 23 until 2:00 AM ET on Friday, May 24 due to maintenance. We apologize for the inconvenience.

Title: Statistics of geomagnetic storms: Global simulations perspective
We present results of 131 geomagnetic storm simulations using the University of Michigan Space Weather Modeling Framework Geospace configuration. We compare the geomagnetic indices derived from the simulation with those observed, and use 2D cuts in the noon-midnight planes to compare the magnetopause locations with empirical models. We identify the location of the current sheet center and look at the plasma parameters to deduce tail dynamics. We show that the simulation produces geomagnetic index distributions similar to those observed, and that their relationship to the solar wind driver is similar to that observed. While the magnitudes of the Dst and polar cap potentials are close to those observed, the simulated AL index is consistently underestimated. Analysis of the magnetopause position reveals that the subsolar position agrees well with an empirical model, but that the tail flaring in the simulation is much smaller than that in the empirical model. The magnetotail and ring currents are closely correlated with the Dst index, and reveal a strong contribution of the tail current beyond 8 R E to the Dst index during the storm main phase.  more » « less
Award ID(s):
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Frontiers in Astronomy and Space Sciences
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    An intriguing aspect of the famous September 2, 1859 geomagnetic disturbance (or “Carrington” event) is the horizontal magnetic (BH) data set measured in Colaba, India (magnetic latitude approximately 20°N). The field exhibits a sharp decrease of over 1,600 nT and a quick recovery of about 1,300 nT, all within a few hours during the daytime. The mechanism behind this has previously been attributed to magnetospheric processes, ionospheric processes or a combination of both. In this study, we outline our efforts to replicate this low‐latitude magnetic field using the Space Weather Modeling Framework. By simulating an extremely high pressure solar wind scenario, we can emulate the low‐latitude surface magnetic signal at Colaba. In our simulation, magnetospheric currents adjacent to the near‐Earth magnetopause and strong Region 1 field‐aligned currents are the main contributors to the large ColabaBH. The rapid recovery ofBHin our simulated scenario is due to the retreat of these magnetospheric currents as the magnetosphere expands, as opposed to ring current dynamics. In addition, we find that the scenario that best emulated the surface magnetic field observations during the Carrington event had a minimum calculated Dst value between −431 and −1,191 nT, indicating that Dst may not be a suitable estimate of storm intensity for this kind of event.

    more » « less
  2. Abstract

    Since the middle of the last decade, UCSD has incorporated magnetic field data in its Institute for Space‐Earth Environmental Research interplanetary scintillation tomographic analysis. These data are extrapolated upward from the solar surface using the Current Sheet Source Surface model (Zhao & Hoeksema, 1995, to provide predictions of the interplanetary field in RTN coordinates. Over the years this technique has become ever more sophisticated, and allows different types of magnetogram data (SOLIS, Global Oscillation Network Group, etc.,) to be incorporated in the field extrapolations. At Earth, these fields can be displayed in a variety of ways, including Geocentric Solar Magnetospheric (GSM) Bx, By, and Bzcoordinates. Displayed daily, the Current Sheet Source Surface model‐derived GSM Bzshows a significant positive correlation with the low‐resolution (few day variation) in situ measurements of the Bzfield. The nano‐Tesla variations of Bzmaximize in spring and fall as Russell and McPherron (1973, have shown. More significantly, we find that the daily variations are correlated with geomagnetic Kp and Dst index variations, and that a decrease from positive to negative Bzhas a high correlation with minor‐to‐moderate geomagnetic storm activity, as defined by NOAA Space Weather Prediction Center planetary Kp values. Here we provide an 11‐year study of the predicted Bzfield, from the extrapolation of the Global Oscillation Network Group‐magnetograms. We provide a skill‐score analysis of the technique's geomagnetic storm prediction capability, which allows forecasts of moderate enhanced geomagnetic storm activity. UCSD and the Korean Space Weather Center currently operate a website that predicts this low‐resolution GSM Bzfield component variation several days in advance.

    more » « less
  3. Abstract

    We report on the solar and interplanetary (IP) causes of the third largest geomagnetic storm (26 August 2018) in solar cycle 24. The underlying coronal mass ejection (CME) originating from a quiescent filament region becomes a 440 km/s magnetic cloud (MC) at 1 au after ∼5 days. The prolonged CME acceleration (for ∼24 hr) coincides with the time profiles of the post‐eruption arcade intensity and reconnected flux. Chen et al. (2019, obtain a lower speed since they assumed that the CME does not accelerate after ∼12 hr. The presence of multiple coronal holes near the filament channel and the high‐speed wind from them seem to have the combined effect of producing complex rotation in the corona and IP medium resulting in a high‐inclination MC. The Dst time profile in the main phase steepens significantly (rapid increase in storm intensity) coincident with the density increase (prominence material) in the second half of the MC. Simulations using the Comprehensive Inner Magnetosphere‐Ionosphere model show that a higher ring current energy results from larger dynamic pressure (density) in MCs. Furthermore, the Dst index is highly correlated with the main‐phase time integral of the ring current injection that includes density, consistent with the simulations. A complex temporal structure develops in the storm main phase if the underlying MC has a complex density structure during intervals of southward IP magnetic field. We conclude that the high intensity of the storm results from the prolonged CME acceleration, complex rotation of the CME flux rope, and the high density in the 1‐au MC.

    more » « less
  4. Abstract

    The occurrence of St. Patrick's Day (17 March) geomagnetic storms during two different years (2013 and 2015) with similar solar flux levels but varying storm intensity provided an opportunity to compare and contrast the responses of the ionosphere‐thermosphere (IT) system to different levels of geomagnetic activity. The evolution of positive ionospheric storms at the southern polar stations Bharati (76.6°S MLAT) and Davis (76.2°S MLAT) and its causative connection to the solar wind driving mechanisms during these storms has been investigated in this paper. During the main phase of both the storms, significant enhancements in TEC and phase scintillation were observed in the magnetic noon/ midnight period at Bharati and Davis. The TEC in the midnight sector on 17 March 2015 was significantly higher compared to that on 17 March 2013, in line with the storm intensity. The TEC enhancements during both the storm events are associated with the formation of the storm‐enhanced densities (SEDs)/tongue of ionization (TOI). The strong and sustained magnetopause erosion led to the prevalence of stronger storm time electric fields (prompt penetration electric field (PPEF)/subauroral polarization streams (SAPS)) for long duration on 17 March 2015. This combined with the action of neutral winds at midlatitudes favored the formation of higher plasma densities in the regions of SED formation on this day. The same was weaker during the 17 March 2013 storm due to the fast fluctuating nature of interplanetary magnetic field (IMF)Bz. This study shows that the duration and extent of magnetopause erosion play an important role in the spatiotemporal evolution of the plasma density distribution in the high‐midlatitude ionosphere.

    more » « less
  5. Abstract

    Previously, Tsurutani and Lakhina (2014, created estimates for a “perfect” interplanetary coronal mass ejection and performed simple calculations for the response of geospace, including. In this study, these estimates are used to drive a coupled magnetohydrodynamic‐ring current‐ionosphere model of geospace to obtain more physically accurate estimates of the geospace response to such an event. The sudden impulse phase is examined and compared to the estimations of Tsurutani and Lakhina (2014, The physics‐based simulation yields similar estimates for Dst rise, magnetopause compression, and equatorialvalues as the previous study. However, results diverge away from the equator.values in excess of 30 nT/s are found as low asmagnetic latitude. Under southward interplanetary magnetic field conditions, magnetopause erosion combines with strong region one Birkeland currents to intensify theresponse. Values obtained here surpass those found in historically recorded events and set the upper threshold of extreme geomagnetically induced current activity at Earth.

    more » « less