Methods: Stealth CMEs represent a particular class of solar eruptions that are clearly distinguished in coronagraph observations, but they don't have a clear source signature. A particular type of stealth CMEs occurs in the trailing current sheet of a previous ejection, therefore, we used the 2.5D MHD package of the code MPI-AMRVAC to numerically simulate consecutive CMEs by imposing shearing motions onto the inner boundary. The initial magnetic configuration consists of a triple arcade structure embedded into a bimodal solar wind, and the sheared polarity inversion line is found in the southern loop system. The mesh was continuously adapted through a refinement method that applies to current carrying structures. We then compared the obtained eruptions with the observed directions of propagation of an initial multiple coronal mass ejection (MCME) event that occurred in September 2009. We further analysed the simulated ejections by tracking the centre of their flux ropes in latitude and their total speed. Radial Poynting flux computation was employed as well to follow the evolution of electromagnetic energy introduced into the system. Results: Changes within 1\% in the shearing speed result in three different scenarios for the second CME, although the preceding eruption seems insusceptible to such small variations. Depending on the applied shearing speed, we thus obtain a failed eruption, a stealth, or a CME driven by the imposed shear, as the second ejection. The dynamics of all eruptions are compared with the observed directions of propagation of an MCME event and a good correlation is achieved. The Poynting flux analysis reveals the temporal variation of the important steps of eruptions. For the first time, a stealth CME is simulated in the aftermath of a first eruption, through changes in the applied shearing speed.

Coronal mass ejections (CMEs) are the main drivers of geomagnetic disturbances, but the effects of their interaction with Earth's 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 simulate several magnetic reconnection processes in the low solar chromosphere/photosphere, the radiation cooling, heat conduction and ambipolar diffusion are all included. Our numerical results indicate that both the high temperature($ \gtrsim 8\times10^4$~K) and low temperature($\sim 10^4$~K) magnetic reconnection events can happen in the low solar atmosphere ($100\sim600$~km above the solar surface). The plasma $\beta$ controlled by plasma density and magnetic fields is one important factor to decide how much the plasma can be heated up. The low temperature event is formed in a high $\beta$ magnetic reconnection process, Joule heating is the main mechanism to heat plasma and the maximum temperature increase is only several thousand Kelvin. The high temperature explosions can be generated in a low $\beta$ magnetic reconnection process, slow and fast-mode shocks attached at the edges of the well developed plasmoids are the main physical mechanisms to heat the plasma from several thousand Kelvin to over $8\times10^4$~K. Gravity in the low chromosphere can strongly hinder the plasmoind instability and the formation of slow-mode shocks in a vertical current sheet. Only small secondary islands are formed; these islands, however, are not well developed as those in the horizontal current sheets. This work can be applied for understanding the heating mechanism in the low solar atmosphere and could possibly be extended to explain the formation of common low temperature EBs ($\sim10^4$~K) and the high tenperature IRIS bombs ($\gtrsim 8\times10^4$) in the future.

In this paper we investigate, by means of two-dimensional magnetohydrodynamic simulations, the impact of temperature-dependent resistivity and thermal conduction on the development of plasmoid instabilities in reconnecting current sheets in the solar corona. We find that the plasma temperature in the current sheet region increases with time and it becomes greater than that in the inflow region. As secondary magnetic islands appear, the highest temperature is not always found at the reconnection $X$-points, but also inside the secondary islands. One of the effects of anisotropic thermal conduction is to decrease the temperature of the reconnecting $X-$points and transfer the heat into the $O-$points, the plasmoids, where it gets trapped. In the cases with temperature-dependent magnetic diffusivity, $\eta \sim T^{-3/2}$, the decrease in plasma temperature at the $X-$points leads to: (i) increase in the magnetic diffusivity until the characteristic time for magnetic diffusion becomes comparable to that of thermal conduction; (ii) increase in the reconnection rate; and, (iii) more efficient conversion of magnetic energy into thermal energy and kinetic energy of bulk motions. These results provide further explanation of the rapid release of magnetic energy into heat and kinetic energy seen during flares and coronal mass ejections. In this work, we demonstrate that the consideration of anisotropic thermal conduction and Spitzer-type, temperature-dependent magnetic diffusivity, as in the real solar corona, are crucially important for explaining the occurrence of fast reconnection during solar eruptions.

We present a new global model of the solar corona, including the low corona, the transition region and the top of chromosphere. The realistic 3D magnetic field is simulated using the data from the photospheric magnetic field measurements. The distinctive feature of the new model is incorporating the MHD Alfven wave turbulence. We assume this turbulence and its non-linear dissipation to be the only momentum and energy source for heating the coronal plasma and driving the solar wind. The difference between the turbulence dissipation efficiency in coronal holes and that in closed field regions is because the non-linear cascade rate degrades in strongly anisotropic turbulence in coronal holes (no inward propagating wave), thus resulting in colder coronal holes with the bi-modal solar wind originating from them. The detailed presentation of the theoretical model is illustrated with the synthetic images for multi-wavelength EUV emission compared with the observations from SDO AIA and Stereo EUVI instruments for the Carrington rotation 2107.

This study aims to provide a reference to different magnetic field models and reconstruction methods for interplanetary coronal mass ejections (ICMEs). In order to understand the differences in the outputs of those models and codes, we analyze 59 events from the Coordinated Data Analysis Workshop (CDAW) list, using four different magnetic field models and reconstruction techniques; force-free fitting (Goldstein,1983,Burlaga,1988,Lepping et al.,1990), magnetostatic reconstruction using a numerical solution to the Grad-Shafranov equation (Hu and Sonnerup, 2001), fitting to a self-similarly expanding cylindrical configuration (Marubashi and Lepping, 2007) and elliptical, non-force free fitting (Hidalgo,2003). The resulting parameters of the reconstructions for the 59 events are compared statistically, as well as in selected case studies. The ability of a method to fit or reconstruct an event is found to vary greatly: the Grad-Shafranov reconstruction is successful for most magnetic clouds (MCs) but for less than 10% of the non-MC ICMEs; the other three methods provide a successful fit for more than 65% of all events. The differences between the reconstruction and fitting methods are discussed, and suggestions are proposed as to how to reduce them. We find that the magnitude of the axial field is relatively consistent across models but not the orientation of the axis of the ejecta. We also find that there are a few cases for which different signs of the magnetic helicity are found for the same event when we do not fix the boundaries, illustrating that this simplest of parameters is not necessarily always well constrained by fitting and reconstruction models. Finally, we look at three unique cases in depth to provide a comprehensive idea of the different aspects of how the fitting and reconstruction codes work.

In 2010 May 23-24, SDO observed the launch of two successive coronal mass ejections (CMEs), which were subsequently tracked by the SECCHI suite onboard STEREO. Using the COR2 coronagraphs and the heliospheric imagers (HIs), the initial direction of both CMEs is determined to be slightly west of the Sun-Earth line. We derive the CME kinematics, including the evolution of the CME expansion until 0.4 AU. We find that, during the interaction, the second CME decelerates from a speed above 500 km/s to 380 km/s the speed of the leading edge of the first CME. STEREO observes a complex structure composed of two different bright tracks in HI2-A but only one bright track in HI2-B. In situ measurements from Wind show an "isolated" ICME, with the geometry of a flux rope preceded by a shock. Measurements in the sheath are consistent with draping around the transient. By combining remote-sensing and in situ measurements, we determine that this event shows a clear instance of deflection of two CMEs after their collision, and we estimate the deflection of the first CME to be about 10 degrees towards the Sun-Earth line. The arrival time, arrival speed and radius at Earth of the first CME are best predicted from remote-sensing observations taken before the collision of the CMEs. Due to the over-expansion of the CME after the collision, there are few, if any, signs of interaction in in situ measurements. This study illustrates that complex interactions during the Sun-to-Earth propagation may not be revealed by in situ measurements alone.

We present a numerical investigation of the coronal evolution of a coronal mass ejection (CME) on 2005 August 22 using a 3-D thermodynamics magnetohydrodynamic model, the SWMF. The source region of the eruption was anemone active region (AR) 10798, which emerged inside a coronal hole. We validate our modeled corona by producing synthetic extreme ultraviolet (EUV) images, which we compare to EIT images. We initiate the CME with an out-of-equilibrium flux rope with an orientation and chirality chosen in agreement with observations of a H-alpha filament. During the eruption, one footpoint of the flux rope reconnects with streamer magnetic field lines and with open field lines from the adjacent coronal hole. It yields an eruption which has a mix of closed and open twisted field lines due to interchange reconnection and only one footpoint line-tied to the source region. Even with the large-scale reconnection, we find no evidence of strong rotation of the CME as it propagates. We study the CME deflection and find that the effect of the Lorentz force is a deflection of the CME by about 3 deg/Rsun towards the East during the first 30 minutes of the propagation. We also produce coronagraphic and EUV images of the CME, which we compare with real images, identifying a dimming region associated with the reconnection process. We discuss the implication of our results for the arrival at Earth of CMEs originating from the limb and for models to explain the presence of open field lines in magnetic clouds.

Interplanetary coronal mass ejections (ICMEs) have complex magnetic and density structures, which is the result of their interaction with the structured solar wind and with previous eruptions. ICMEs are revealed by in situ measurements and in the past five years, through remote-sensing observations by heliospheric imagers. However, to understand and analyze these observations often requires the use of numerical modeling. It is because no instruments can yet provide a simple view of ICMEs in two or three dimensions. Numerical simulations can be used to determine the origin of a complex ejecta observed near Earth, or to analyze the origin, speed and extent of density structures observed remotely. Here, we review and discuss recent efforts to use numerical simulations of ICMEs to investigate the magnetic topology, density structure, energetics and kinematics of ICMEs in the interplanetary space. After reviewing existing numerical models of ICMEs, we first focus on numerical modeling in support of the SMEI and STEREO observations. 3-D simulations can help determining the origins of the fronts observed by SECCHI and SMEI, especially for complex events such as the January 24-25, 2007 CMEs. They can also be used to test different methods to derive ICME properties from remote observations, to predict and explain observational effects, and to understand the deceleration and deformation of ICMEs. In the last part, we focus on the numerical investigation of non-magnetic cloud ejecta. We discuss how simulations are crucial to understand the formation of non-twisted ejecta and the formation of complex ejecta due to the interaction of multiple ICMEs.

We discuss how simultaneous observations by multiple heliospheric imagers can provide some important information about the azimuthal properties of Coronal Mass Ejections (CMEs) in the heliosphere. We propose two simple models of CME geometry that can be used to derive information about the azimuthal deflection and the azimuthal expansion of CMEs from SECCHI/HI observations. We apply these two models to four CMEs well-observed by both STEREO spacecraft during the year 2008. We find that in three cases, the joint STEREO-A and B observations are consistent with CMEs moving radially outward. In some cases, we are able to derive the azimuthal cross-section of the CME fronts, and we are able to measure the deviation from self-similar evolution. The results from this analysis show the importance of having multiple satellites dedicated to space weather forecasting, for example in orbits at the Lagrangian L4 and L5 points.

In this work we describe our implementation of a thermodynamic energy equation into the global corona model of the Space Weather Modeling Framework (SWMF), and its development into the new Lower Corona (LC) model. This work includes the integration of the additional energy transport terms of coronal heating, electron heat conduction, and optically thin radiative cooling into the governing magnetohydrodynamic (MHD) energy equation. We examine two different boundary conditions using this model; one set in the upper transition region (the Radiative Energy Balance model), as well as a uniform chromospheric condition where the transition region can be modeled in its entirety. Via observation synthesis from model results and the subsequent comparison to full sun extreme ultraviolet (EUV) and soft X-Ray observations of Carrington Rotation (CR) 1913 centered on Aug 27, 1996, we demonstrate the need for these additional considerations when using global MHD models to describe the unique conditions in the low corona. Through multiple simulations we examine ability of the LC model to asses and discriminate between coronal heating models, and find that a relative simple empirical heating model is adequate in reproducing structures observed in the low corona. We show that the interplay between coronal heating and electron heat conduction provides significant feedback onto the 3D magnetic topology in the low corona as compared to a potential field extrapolation, and that this feedback is largely dependent on the amount of mechanical energy introduced into the corona.

Most simulations of coronal mass ejections (CMEs) to date either focus on the interplanetary propagation of a giant plasma "blob" without paying too much attention to its origin and to the formation process or they focus on the complex evolution of the coronal magnetic field due to (sub-)photospheric motions which result in an eruption. Here, we present global simulations of CMEs where coronal motions are used to produce a realistic evolution of the coronal magnetic field and cause an eruption. We focus on active region 10069, which produced a number of eruptions in late August 2002, including the August 24, 2002 CME - a fast (~2000 km/s) eruption originating from W81-, as well as a slower eruption on August 22, 2002 (originating from W62). Using a three-dimensional magneto-hydrodynamic (MHD) simulation of these ejections with the Space Weather Modeling Framework (SWMF), we show how a realistic initiation mechanism enables us to study the deflection of the CME in the corona and in the heliosphere. Reconnection of the erupting magnetic field with that of neighboring streamers and active regions modify the solar connectivity of the field lines connecting to Earth and change the expected solar energetic particle fluxes. Comparing the results at 1 AU of our simulations with in situ observations by the ACE spacecraft, we propose an alternate solar origin for the shock wave observed at L1 on August 26.

We present general considerations regarding the derivation of the radial distances of coronal mass ejections (CMEs) from elongation angle measurements such as those provided by SECCHI and SMEI, focusing on measurements in the Heliospheric Imager 2 (HI-2) field of view (i.e. past 0.3 AU). This study is based on a three-dimensional (3-D) magneto-hydrodynamics (MHD) simulation of two CMEs observed by SECCHI on January 24-27, 2007. Having a 3-D simulation with synthetic HI images, we are able to compare the two basic methods used to derive CME positions from elongation angles, the so-called "Point-P" and "Fixed-Phi" approximations. We confirm, following similar works, that both methods, while valid in the most inner heliosphere, yield increasingly large errors in HI-2 field of view for fast and wide CMEs. Using a simple model of a CME as an expanding self-similar sphere, we derive an analytical relationship between elongation angles and radial distances for wide CMEs. This relationship is simply the harmonic mean of the "Point-P" and "Fixed-Phi'' approximations and it is aimed at complementing 3-D fitting of CMEs by cone models or flux rope shapes. It proves better at getting the kinematics of the simulated CME right when we compare the results of our line-of-sights to the MHD simulation. Based on this approximation, we re-analyze the J-maps (time-elongation maps) in January 26-27, 2007 and present the first observational evidence that the merging of CMEs is associated with a momentum exchange from the faster ejection to the slower one due to the propagation of the shock wave associated with the fast eruption through the slow eruption.

We discuss how some coronal mass ejections (CMEs) originating from the western limb of the Sun are associated with space weather effects such as solar energetic particles (SEPs), shock or geo-effective ejecta at Earth. We focus on the August 24, 2002 coronal mass ejection, a fast (~ 2000 km/s) eruption originating from W81. Using a three-dimensional magneto-hydrodynamic simulation of this ejection with the Space Weather Modeling Framework (SWMF), we show how a realistic initiation mechanism enables us to study the deflection of the CME in the corona and the heliosphere. Reconnection of the erupting magnetic field with that of neighboring streamers and active regions modify the solar connectivity of the field lines connecting to Earth and can also partly explain the deflection of the eruption during the first tens of minutes. Comparing the results at 1 AU of our simulation with observations by the ACE spacecraft, we find that the simulated shock does not reach Earth, but has a maximum angular span of about 120$^\circ$, and reaches 35$^\circ$ West of Earth in 58 hours. We find no significant deflection of the CME and its associated shock wave in the heliosphere, and we discuss the consequences for the shock angular span.

We numerically model the coronal mass ejection (CME) event of October 28, 2003 that erupted from active region 10486 and propagated to Earth in less than 20 hours causing severe geomagnetic storms. The magnetohydrodynamic (MHD) model is formulated by first arriving at a steady state corona and solar wind employing synoptic magnetograms. We initiate two CMEs from the same active region, one approximately a day earlier that preconditions the solar wind for the much faster CME on the 28th. This second CME travels through the corona at a rate of over 2500 km s$^{-1}$ driving a strong forward shock. We clearly identify this shock in an image produced by the Large Angle Spectrometric Coronagraph (LASCO) C3, and reproduce the shock and its appearance in synthetic white light images from the simulation. We find excellent agreement with both the general morphology and the quantitative brightness of the model CME with LASCO observations. These results demonstrate that the CME shape is largely determined by its interaction with the ambient solar wind and may not be sensitive to the initiation process. We then show how the CME would appear as observed by wide-angle coronagraphs onboard the Solar Terrestrial Relations Observatory (STEREO) spacecraft. We find complex time evolution of the white-light images as a result of the way in which the density structures pass through the Thomson sphere. The simulation is performed with the Space Weather Modeling Framework (SWMF).