No Arabic abstract
Coronal Mass Ejections (CMEs) are major drivers of extreme space weather conditions, this being a matter of serious concern for our modern technologically-dependent society. Development of numerical approaches that would simulate CME generation and propagation through the interplanetary space is an important step towards our capability to predict CME arrival times at Earth and their geo-effectiveness. In this paper, we utilize a data-constrained Gibson--Low (GL) flux rope model to generate CMEs. We derive the geometry of the initial GL flux rope using the Graduated Cylindrical Shell (GCS) method. This method uses multiple viewpoints from STEREO A & B Cor1/Cor2, and SOHO/LASCO C2/C3 coronagraphs to determine the size and orientation of a CME flux rope as it starts to erupt from the Sun. A flux rope generated in this way is inserted into a quasi-steady global magnetohydrodynamics (MHD) background solar wind flow driven by SDO/HMI line-of-sight magnetogram data, and erupts immediately. Numerical results obtained with the Multi-Scale Fluid-Kinetic Simulation Suite (MS-FLUKSS) code are compared with STEREO and SOHO/LASCO coronagraph observations in particular in terms of the CME speed, acceleration, and magnetic field structure.
Coronal mass ejections (CMEs) are the main drivers of geomagnetic disturbances, but the effects of their interaction with Earths magnetic field depend on their magnetic configuration and orientation. Fitting and reconstruction techniques have been developed to determine the important geometrical and physical CME properties. In many instances, there is disagreement between such different methods but also between fitting from in situ measurements and reconstruction based on remote imaging. Here, we compare three methods based on different assumptions for measurements of thirteen CMEs by the Wind spacecraft from 1997 to 2015. These CMEs are selected from the interplanetary coronal mass ejections catalog on https://wind.nasa.gov/ICMEindex.php due to their simplicity in terms of 1) small expansion speed throughout the CME and 2) little asymmetry in the magnetic field profile. This makes these thirteen events ideal candidates to compare codes that do not include expansion nor distortion. We find that, for these simple events, the codes are in relatively good agreement in terms of the CME axis orientation for six out of the 13 events. Using the Grad-Shafranov technique, we can determine the shape of the cross-section, which is assumed to be circular for the other two models, a force-free fitting and a circular-cylindrical non-force-free fitting. Five of the events are found to have a clear circular cross-section, even when this is not a pre-condition of the reconstruction. We make an initial attempt at evaluating the adequacy of the different assumptions for these simple CMEs. The conclusion of this work strongly suggests that attempts at reconciling in situ and remote-sensing views of CMEs must take in consideration the compatibility of the different models with specific CME structures to better reproduce flux ropes.
We present a first-principles-based coronal mass ejection (CME) model suitable for both scientific and operational purposes by combining a global magnetohydrodynamics (MHD) solar wind model with a flux rope-driven CME model. Realistic CME events are simulated self-consistently with high fidelity and forecasting capability by constraining initial flux rope parameters with observational data from GONG, SOHO/LASCO, and STEREO/COR. We automate this process so that minimum manual intervention is required in specifying the CME initial state. With the newly developed data-driven Eruptive Event Generator Gibson-Low (EEGGL), we present a method to derive Gibson-Low (GL) flux rope parameters through a handful of observational quantities so that the modeled CMEs can propagate with the desired CME speeds near the Sun. A test result with CMEs launched with different Carrington rotation magnetograms are shown. Our study shows a promising result for using the first-principles-based MHD global model as a forecasting tool, which is capable of predicting the CME direction of propagation, arrival time, and ICME magnetic field at 1 AU (see companion paper by Jin et al. 2016b).
We report on a numerical investigation of two coronal mass ejections (CMEs) which interact as they propagate in the inner heliosphere. We focus on the effect of the orientation of the CMEs relative to each other by performing four different simulations with the axis of the second CME rotated by 90 degrees from one simulation to the next. Each magneto-hydrodynamic (MHD) simulation is performed in three dimensions (3-D) with the Space Weather Modeling Framework (SWMF) in an idealized setting reminiscent of solar minimum conditions. We extract synthetic satellite measurements during and after the interaction and compare the different cases. We also analyze the kinematics of the two CMEs, including the evolution of their widths and aspect ratios. We find that the first CME contracts radially as a result of the interaction in all cases, but the amount of subsequent radial expansion depends on the relative orientation of the two CMEs. Reconnection between the two ejecta and between the ejecta and the interplanetary magnetic field (IMF) determines the type of structure resulting from the interaction. When a CME with a high inclination with respect to the ecliptic overtakes one with a low inclination, it is possible to create a compound event with a smooth rotation in the magnetic field vector over more than 180 degrees. Due to reconnection, the second CME only appears as an extended tail, and the event may be mistaken for a glancing encounter with an isolated CME. This configuration differs significantly from the one usually studied of a multiple-magnetic cloud event, which we found to be associated with the interaction of two CMEs with the same orientation.
We study interplanetary coronal mass ejections (ICMEs) measured by probes at different heliocentric distances (0.3-1AU) to investigate the propagation of ICMEs in the inner heliosphere and determine how the generic features of ICMEs change with heliospheric distance. Using data from the MESSENGER, Venus Express and ACE spacecraft, we analyze with the superposed epoch technique the profiles of ICME substructures, namely the sheath and the magnetic ejecta. We determine that the median magnetic field magnitude in the sheath correlates well with ICME speeds at 1 AU and we use this proxy to order the ICMEs at all spacecraft. We then investigate the typical ICME profiles for three categories equivalent to slow, intermediate and fast ICMEs. Contrary to fast ICMEs, slow ICMEs have a weaker solar wind field at the front and a more symmetric magnetic field profile. We find the asymmetry to be less pronounced at Earth than at Mercury, indicating a relaxation taking place as ICMEs propagate. We also find that the magnetic field intensities in the wake region of the ICMEs do not go back to the pre-ICME solar wind intensities, suggesting that the effects of ICMEs on the ambient solar wind last longer than the duration of the transient event. Such results provide an indication of physical processes that need to be reproduced by numerical simulations of ICME propagation. The samples studied here will be greatly improved by future missions dedicated to the exploration of the inner heliosphere, such as Parker Solar Probe and Solar Orbiter.
Stealth coronal mass ejections (CMEs) are events in which there are almost no observable signatures of the CME eruption in the low corona but often a well-resolved slow flux rope CME observed in the coronagraph data. We present results from a three-dimensional numerical magnetohydrodynamics (MHD) simulation of the 1--2 June 2008 slow streamer blowout CME that Robbrecht et al. (2009) called the CME from nowhere. We model the global coronal structure using a 1.4 MK isothermal solar wind and a low-order potential field source surface representation of the Carrington Rotation 2070 magnetogram synoptic map. The bipolar streamer belt arcade is energized by simple shearing flows applied in the vicinity of the helmet streamers polarity inversion line. The flows are large scale and impart a shear typical of that expected from the differential rotation. The slow expansion of the energized helmet streamer arcade results in the formation of a radial current sheet. The subsequent onset of expansion-induced flare reconnection initiates the stealth CME while gradually releasing the stored magnetic energy. We present favorable comparisons between our simulation results and the multiviewpoint SOHO-LASCO (Large Angle and Spectrometric Coronagraph) and STEREO-SECCHI (Sun Earth Connection Coronal and Heliospheric Investigation) coronagraph observations of the preeruption streamer structure and the initiation and evolution of the stealth streamer blowout CME.