skip to main content


Title: Application of a Modified Spheromak Model to Simulations of Coronal Mass Ejection in the Inner Heliosphere
Abstract

The magnetic fields of interplanetary coronal mass ejections (ICMEs), which originate close to the Sun in the form of a flux rope, determine their geoeffectiveness. Therefore, robust flux rope‐based models of CMEs are required to perform magnetohydrodynamic (MHD) simulations aimed at space weather predictions. We propose a modified spheromak model and demonstrate its applicability to CME simulations. In this model, such properties of a simulated CME as the poloidal and toroidal magnetic fluxes, and the helicity sign can be controlled with a set of input parameters. We propose a robust technique for introducing CMEs with an appropriate speed into a background, MHD solution describing the solar wind in the inner heliosphere. Through a parametric study, we find that the speed of a CME is much more dependent on its poloidal flux than on the toroidal flux. We also show that the CME speed increases with its total energy, giving us control over its initial speed. We further demonstrate the applicability of this model to simulations of CME‐CME collisions. Finally, we use this model to simulate the 12 July 2012 CME and compare the plasma properties at 1 AU with observations. The predicted CME properties agree reasonably with observational data.

 
more » « less
NSF-PAR ID:
10375434
Author(s) / Creator(s):
 ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Space Weather
Volume:
18
Issue:
5
ISSN:
1542-7390
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    We provide exact analytical solutions for the magnetic field produced by prescribed current distributions located inside a toroidal filament of finite thickness. The solutions are expressed in terms of toroidal functions, which are modifications of the Legendre functions. In application to the MHD equilibrium of a twisted toroidal current loop in the solar corona, the Grad–Shafranov equation is decomposed into an analytic solution describing an equilibrium configuration against the pinch-effect from its own current and an approximate solution for an external strapping field to balance the hoop force. Our solutions can be employed in numerical simulations of coronal mass ejections (CMEs). When superimposed on the background solar coronal magnetic field, the excess magnetic energy of the twisted current loop configuration can be made unstable by applying flux cancellation to reduce the strapping field. Such loss of stability accompanied by the formation of an expanding flux rope is typical for the Titov & Démoulin eruptive event generator. The main new features of the proposed model are as follows: the filament is filled with finiteβplasma with finite mass and energy, the model describes an equilibrium solution that will spontaneously erupt due to magnetic reconnection of the strapping magnetic field arcade, and there are analytic expressions connecting the model parameters to the asymptotic velocity and total mass of the resulting CME, providing a way to connect the simulated CME properties to multipoint coronograph observations.

     
    more » « less
  2. Graduated Cylindrical Shell (GCS) model is a widely used tool to determine direction, kinematic and orientation properties of Coronal Mass Ejections (CME) using multi-viewpoint observations from SOHO and STEREO A&B coronagraphs. In this study, we estimate the subjective uncertainties typically seen while deriving these CME properties by comparing the GCS model results reported in multiple studies and catalogs for 56 CMEs. We find that the GCS estimates of latitude, longitude, and tilt show an average uncertainty of 5.7, 11.2, and 24.7 degrees with standard deviation of 5.5, 12.7, and 19.7 degrees respectively. We found that the uncertainties in estimated latitudes are correlated with uncertainties in estimated longitude, tilt, and speed, showing that some CMEs are inherently difficult to fit than others. We then introduced these uncertainty values in our 3-D magnetohydrodynamic flux rope based modified spheromak CME model to figure out their consequences for space weather prediction. We find that much better CME observations are required to reliably predict magnetic field of CMEs at 1 AU using flux rope based models, since the uncertainties in estimated GCS values can result in large differences in 1 AU signatures, especially for CMEs launched away from the Sun-Earth line. 
    more » « less
  3. 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
  4. 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
  5. Abstract

    The trajectories of coronal mass ejections (CMEs) are often seen to deviate substantially from a purely radial propagation direction. Such deviations occur predominantly in the corona and have been attributed to “channeling” or deflection of the eruptive flux by asymmetric ambient magnetic fields. Here, we investigate an additional mechanism that does not require any asymmetry of the preeruptive ambient field. Using magnetohydrodynamic numerical simulations, we show that the trajectories of CMEs through the solar corona can significantly deviate from the radial direction when propagation takes place in a unipolar radial field. We demonstrate that the deviation is most prominent below ∼15Rand can be attributed to an “effectiveI×Bforce” that arises from the intrusion of a magnetic flux rope with a net axial electric current into a unipolar background field. These results are important for predictions of CME trajectories in the context of space-weather forecasts, as well as for reaching a deeper understanding of the fundamental physics underlying CME interactions with the ambient fields in the extended solar corona.

     
    more » « less