skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Decent Estimate of CME Arrival Time From a Data‐Assimilated Ensemble in the Alfvén Wave Solar Atmosphere Model (DECADE‐AWSoM)
Abstract Forecasting the arrival time of Earth‐directed coronal mass ejections (CMEs) via physics‐based simulations is an essential but challenging task in space weather research due to the complexity of the underlying physics and limited remote and in situ observations of these events. Data assimilation techniques can assist in constraining free model parameters and reduce the uncertainty in subsequent model predictions. In this study, we show that CME simulations conducted with the Space Weather Modeling Framework (SWMF) can be assimilated with SOHO LASCO white‐light (WL) observations and solar wind observations at L1 prior to the CME eruption to improve the prediction of CME arrival time. The L1 observations are used to constrain the model of the solar wind background into which the CME is launched. Average speed of CME shock front over propagation angles are extracted from both synthetic WL images from the Alfvén Wave Solar atmosphere Model (AWSoM) and the WL observations. We observe a strong rank correlation between the average WL speed and CME arrival time, with the Spearman's rank correlation coefficients larger than 0.90 for three events occurring during different phases of the solar cycle. This enables us to develop a Bayesian framework to filter ensemble simulations using WL observations, which is found to reduce the mean absolute error of CME arrival time prediction from about 13.4 to 5.1 hr. The results show the potential of assimilating readily available L1 and WL observations within hours of the CME eruption to construct optimal ensembles of Sun‐to‐Earth CME simulations.  more » « less
Award ID(s):
2323303 2027555
PAR ID:
10567951
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Space Weather
Volume:
23
Issue:
1
ISSN:
1542-7390
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract We present observations and modeling results of the propagation and impact at Earth of a high-latitude, extended filament channel eruption that commenced on 2015 July 9. The coronal mass ejection (CME) that resulted from the filament eruption was associated with a moderate disturbance at Earth. This event could be classified as a so-called “problem storm” because it lacked the usual solar signatures that are characteristic of large, energetic, Earth-directed CMEs that often result in significant geoeffective impacts. We use solar observations to constrain the initial parameters and therefore to model the propagation of the 2015 July 9 eruption from the solar corona up to Earth using 3D magnetohydrodynamic heliospheric simulations with three different configurations of the modeled CME. We find the best match between observed and modeled arrival at Earth for the simulation run that features a toroidal flux rope structure of the CME ejecta, but caution that different approaches may be more or less useful depending on the CME–observer geometry when evaluating the space weather impact of eruptions that are extreme in terms of their large size and high degree of asymmetry. We discuss our results in the context of both advancing our understanding of the physics of CME evolution and future improvements to space weather forecasting. 
    more » « less
  2. Abstract Modeling the impact of space weather events such as coronal mass ejections (CMEs) is crucial to protecting critical infrastructure. The Space Weather Modeling Framework is a state‐of‐the‐art framework that offers full Sun‐to‐Earth simulations by computing the background solar wind, CME propagation, and magnetospheric impact. However, reliable long‐term predictions of CME events require uncertainty quantification (UQ) and data assimilation. We take the first steps by performing global sensitivity analysis (GSA) and UQ for background solar wind simulations produced by the Alfvén Wave Solar atmosphere Model (AWSoM) for two Carrington rotations: CR2152 (solar maximum) and CR2208 (solar minimum). We conduct GSA by computing Sobol' indices that quantify contributions from model parameter uncertainty to the variance of solar wind speed and density at 1 au, both crucial quantities for CME propagation and strength. Sobol' indices also allow us to rank and retain only the most important parameters, which aids in the construction of smaller ensembles for the reduced‐dimension parameter space. We present an efficient procedure for computing the Sobol' indices using polynomial chaos expansion surrogates and space‐filling designs. The PCEs further enable inexpensive forward UQ. Overall, we identify three important model parameters: the multiplicative factor applied to the magnetogram, Poynting flux per magnetic field strength constant used at the inner boundary, and the coefficient of the perpendicular correlation length in the turbulent cascade model in AWSoM. 
    more » « less
  3. Abstract EUropean Heliospheric FORecasting Information Asset (EUHFORIA) is a physics‐based data‐driven solar wind and coronal mass ejections (CMEs) propagation model designed for space weather forecasting and event analysis investigations. Although EUHFORIA can predict the solar wind plasma and magnetic field properties at Earth, it is not equipped to quantify the geo‐effectiveness of the solar transients in terms of geomagnetic indices like the disturbance storm time (Dst) index and the auroral indices, that quantify the impact of the magnetized plasma encounters on Earth's magnetosphere. Therefore, we couple EUHFORIA with the Open Geospace General Circulation Model (OpenGGCM), a magnetohydrodynamic model of the response of Earth's magnetosphere, ionosphere, and thermosphere to transient solar wind characteristics. In this coupling, OpenGGCM is driven by the solar wind and interplanetary magnetic field obtained from EUHFORIA simulations to produce the magnetospheric and ionospheric response to the CMEs. This coupling is validated with two observed geo‐effective CME events driven with the spheromak flux‐rope CME model. We compare these simulation results with the indices obtained from OpenGGCM simulations driven by the measured solar wind data from spacecraft. We further employ the dynamic time warping (DTW) technique to assess the model performance in predicting Dst. The main highlight of this study is to use EUHFORIA simulated time series to predict the Dst and auroral indices 1–2 days in advance, as compared to using the observed solar wind data at L1, which only provides predictions 1–2 hr before the actual impact. 
    more » « less
  4. Context.Predicting geomagnetic events starts with an understanding of the Sun-Earth chain phenomena in which (interplanetary) coronal mass ejections (CMEs) play an important role in bringing about intense geomagnetic storms. It is not always straightforward to determine the solar source of an interplanetary coronal mass ejection (ICME) detected at 1 au. Aims.The aim of this study is to test by a magnetohydrodynamic (MHD) simulation the chain of a series of CME events detected from L1 back to the Sun in order to determine the relationship between remote and in situ CMEs. Methods.We analysed both remote-sensing observations and in situ measurements of a well-defined magnetic cloud (MC) detected at L1 occurring on 28 June 2013. The MHD modelling is provided by the 3D MHD European Heliospheric FORecasting Information Asset (EUHFORIA) simulation model. Results.After computing the background solar wind, we tested the trajectories of six CMEs occurring in a time window of five days before a well-defined MC at L1 that may act as the candidate of the MC. We modelled each CME using the cone model. The test involving all the CMEs indicated that the main driver of the well-defined, long-duration MC was a slow CME. For the corresponding MC, we retrieved the arrival time and the observed proton density. Conclusions.EUHFORIA confirms the results obtained in the George Mason data catalogue concerning this chain of events. However, their proposed solar source of the CME is disputable. The slow CME at the origin of the MC could have its solar source in a small, emerging region at the border of a filament channel at latitude and longitude equal to +14 degrees. 
    more » « less
  5. Abstract Flux-rope-based magnetohydrodynamic modeling of coronal mass ejections (CMEs) is a promising tool for prediction of the CME arrival time and magnetic field at Earth. In this work, we introduce a constant-turn flux rope model and use it to simulate the 2012 July 12 16:48 CME in the inner heliosphere. We constrain the initial parameters of this CME using the graduated cylindrical shell (GCS) model and the reconnected flux in post-eruption arcades. We correctly reproduce all the magnetic field components of the CME at Earth, with an arrival time error of approximately 1 hr. We further estimate the average subjective uncertainties in the GCS fittings by comparing the GCS parameters of 56 CMEs reported in multiple studies and catalogs. We determined that the GCS estimates of the CME latitude, longitude, tilt, and speed have average uncertainties of 5.°74, 11.°23, 24.°71, and 11.4%, respectively. Using these, we have created 77 ensemble members for the 2012 July 12 CME. We found that 55% of our ensemble members correctly reproduce the sign of the magnetic field components at Earth. We also determined that the uncertainties in GCS fitting can widen the CME arrival time prediction window to about 12 hr for the 2012 July 12 CME. On investigating the forecast accuracy introduced by the uncertainties in individual GCS parameters, we conclude that the half-angle and aspect ratio have little impact on the predicted magnetic field of the 2012 July 12 CME, whereas the uncertainties in longitude and tilt can introduce relatively large spread in the magnetic field predicted at Earth. 
    more » « less