The Physical Processes of CME/ICME Evolution

As observed in Thomson-scattered white light, coronal mass ejections (CMEs) are manifest as large-scale expulsions of plasma magnetically driven from the corona in the most energetic eruptions from the Sun. It remains a tantalizing mystery as to how these erupting magnetic fields evolve to form the complex structures we observe in the solar wind at Earth. Here, we strive to provide a fresh perspective on the post-eruption and interplanetary evolution of CMEs, focusing on the physical processes that define the many complex interactions of the ejected plasma with its surroundings as it departs the corona and propagates through the heliosphere. We summarize the ways CMEs and their interplanetary CMEs (ICMEs) are rotated, reconfigured, deformed, deflected, decelerated and disguised during their journey through the solar wind. This study then leads to consideration of how structures originating in coronal eruptions can be connected to their far removed interplanetary counterparts. Given that ICMEs are the drivers of most geomagnetic storms (and the sole driver of extreme storms), this work provides a guide to the processes that must be considered in making space weather forecasts from remote observations of the corona.


Introduction
Coronal Mass Ejections were first observed from space with the coronagraph onboard NASA's Seventh Orbiting Solar Observatory (OSO-7) on 14 December 1971, Tousey (1973) as bright transients expelled through the coronagraph field of view over a period of minutes to hours. Subsequent space-borne coronagraphs such as the Coronagraph/Polarimeter (or C/P) on the Solar Maximum Mission (SMM) (MacQueen et al. 1980) and the Large Angle Spectrometric Coronagraph (LASCO) onboard the Solar and Heliosphere Observatory (SOHO) (Brueckner et al. 1995) have observed thousands of CMEs from which their characteristics are documented over more than two solar cycles (e.g., Hundhausen 1993;). An analysis of these data has led to a very detailed understanding of the structure and evolution of CMEs in the corona, which has been summarized in several reviews (e.g., Hundhausen et al. 1984;Kahler 1987;Hundhausen 1987;Kahler 1992;Gosling 1993a). Related reviews also discuss the pre-event conditions leading to CMEs , while the theoretic underpinnings of CME initiation are treated in e.g., Forbes (2000), Forbes et al. (2006). The structure of interplanetary CMEs (ICMEs) is well described in Gopalswamy (2006) and by Kilpua et al. in this issue. For the purpose of this paper, we define CMEs as transients that occur within the field-of-view of classical coronagraphs, which may extend to maximum of 30 R (solar radii) in the case of the LASCO C3 coronagraph. Transients beyond this range, and certainly beyond orbit of Mercury, we classify as ICMEs.
A fresh look at the literature on CMEs and ICMEs is timely given the enormous advances that have occurred in the past decade. We have just passed the ten year anniversary of the launch of the twin Solar Terrestrial Relations Observatory (STEREO) spacecraft (Kaiser et al. 2008), which provide continuous multi-viewpoint white-light observations of CMEs from Sun to Earth . STEREO also provides in situ measurements including both the plasma and the magnetic field made with the In-situ Measurements of Particles and CME Transients (IMPACT) instrument (Luhmann et al. 2008;. With these capabilities, STEREO was ideally designed to connect solar eruptions with their solar wind disturbances. Of similar importance is the introduction of massively-parallel supercomputers in the late 1990's, which have allowed the first realistic three-dimensional magnetohydrodynamic (MHD) simulations of specific CME events (e.g., Odstrcil et al. 2005;Taktakishvili et al. 2009). A review by Webb and Howard (2012) summarizes observational studies of CMEs utilizing data from STEREO among many other spacecraft while a companion paper by Chen (2011) summarizes recent advances in numerical modeling of CME initiation. We complement these works by reviewing the results of both observational studies and numerical simulations to garner a more complete understanding of the physical processes governing the evolution of CMEs in the corona and ICMEs in the solar wind. We do so by following the sequence of events that affect the structure and velocity of CMEs as they depart the closed fields of the corona to pass through the depths of interplanetary space.
As we take inventory of relevant physical processes, we first consider the structure and appearance of CMEs and ICMEs from the low corona to 1 AU in Sect. 2. We then cover those processes that are most pronounced in close proximity to the Sun, where CMEs make their way through the highly structured magnetic field of the low corona. In this regard, we consider deflections and rotations (in Sects. 3 and 4, respectively). In Sect. 5, we examine the kinematic evolution followed by a discussion of the impact of ICMEs on the surrounding heliosphere in Sect. 6. We then cover the processes that affect the magnetic structure of ICMEs in Sect. 7, most significantly reconnection, which causes magnetic ejecta to erode. In Sect. 8, we examine the charge state composition of CMEs with a special emphasis on filament material. We consider the effects of CME-CME interaction and the formation of complex ejecta in Sect. 9, and in Sect. 10, we describe a Sun-to-Earth simulation of the Bastille Day CME event. Finally, in Sect. 11, we summarize the salient points of CME and ICME evolution. As we elucidate these processes, we endeavor to first introduce them as they were discovered, and then describe supporting work in a historical narrative. Understanding the evolution of CMEs and ICMEs in the corona and solar wind is fundamental to explaining and potentially predicting the complex structures observed throughout the heliosphere. This fundamental goal of heliospheric physics has great practical value in space weather forecasting capability, which is addressed by Kilpua et al. in this issue. Similarly the topic of solar energetic particles (SEPs) will be covered in two complementary texts, one by Dalla and Klein and another by Schwadron et al.

Structure of CMEs to ICMEs
Manifest in coronagraph images, CMEs are seen in Thomson-scattered white-light, where the brightness reflects the electron density near the plane of the sky (Billings 1966). While CMEs may show many morphologies, the simplest configuration is a three-part structure: a bright leading loop enclosing a dark low-density cavity, which contains a high-density core (e.g., Hundhausen 1993;). This three-part structure can usually be traced directly to a progenitor of the same form that may exist at small-scale confined to an active region or exist at global-scale contained within a helmet streamer. The entire system, with a mass in the range of 10 15 -10 16 g (Colaninno and Vourlidas 2009) may lift out of the corona gradually over a period of hours with speeds less than 100 km s −1 , or they may be impulsively ejected with speeds approaching 3000 km s −1 (Hundhausen et al. 1994), with kinetic energy approaching 10 33 ergs. The journey continues beyond the view of classical coronagraphs where the corresponding ejecta in the solar wind are identified as ICMEs that may or may not bare a clear connection to a coronal counterpart (e.g., Gopalswamy et al. 1998).
Disturbances in the solar wind have long been associated with various eruptive phenomena, such as flares, eruptive prominences, type II radio bursts, and later with coronal mass ejections. Since the 1970's efforts have been made to predict the arrival time and impact speed of the eruption-related solar wind disturbances (e.g., De Young andHundhausen 1971, 1973;Steinolfson and Dryer 1978;Wu et al. 1979). However, at that time, CMEs were not recognized as a major source of interplanetary disturbances and research was primarily focused on the propagation of flare-driven shocks. With time, the CME and more precisely, the ICME ejecta and sheath came to be seen as the source of the most significant solar wind disturbances and consequently the source of non-recurrent geomagnetic storms (e.g., Gosling 1993b). This transformation in understanding was reflected in the evolution of modeling efforts from flare-driven shocks (also called blast waves) to CMEs (Wu et al. 1981;Wei 1982;Smart and Shea 1985;Smith and Dryer 1990;Wei and Dryer 1991;Gosling and Riley 1996;Vandas et al. 1996, e.g., and references therein; see also reviews by Pizzo 1985 andDryer et al. 1988).
The connection between CME and ICME was made abundantly more clear through STEREO observations (e.g., Möstl et al. 2009;Liu et al. 2010a). An example of a particularly well observed CME/ICME pair is found in the 12-18 December 2008 event, which serves to highlight a range of physical processes governing the evolution of these phenomena (e.g., Byrne et al. 2010;Liu et al. 2010a,b;DeForest et al. 2011DeForest et al. , 2013Howard and DeForest 2012). In this case, STEREO data is utilized to make a convincing connection between the solar eruption and the disturbance observed at Earth. The CME is induced by a prominence eruption in the northern hemisphere (see left panel of Fig. 1), which started between 03-04 UT on 12 December 2008. The prominence material (visible in EUVI at 304 Å) is well aligned with the CME core. The CME slowly rotates and expands toward the ecliptic plane, and seems fully developed in COR2. The basic structure of the CME remains organized out to at least the field of view of HI1. In HI2 of STEREO A, we see a dark cavity bracketed by structures with enhanced densities. The time-elongation maps shown in Fig. 1 are produced by stacking running difference intensities of COR2, HI1 and HI2 within a slit along the ecliptic plane. Two features corresponding to the CME leading and trailing edges can be identified up to 50°elongation for both STEREO A and B. Intermittent ones between the two tracks, probably associated with the CME core, are also seen but later disappear presumably due to the expansion of the ICME. Figure 2 shows the plasma and magnetic signatures of the corresponding ICME, which passed Wind on 17 December 2008. The shaded region identifies the magnetic ejecta of the CME.

ICMEs and Magnetic Clouds
The observed structure and evolution of CMEs/ICMEs, as illustrated in Figs. 1, 2 and 3, can be described by physical models, which can reproduce observed properties and explain and predict their evolution. For CMEs with a common three-part structure, the magnetic field is commonly taken to be of the form of a twisted flux rope contained within an interior plasma cavity (e.g., Gibson et al. 2010). The core of the structure is typically considered to be filament material that was supported by the magnetic field above the system's photospheric polarity inversion line (PIL) prior to the eruption. While long-standing, it is worth noting recent work by Howard et al. (2017) questions the filament-core connection for some events. Flux ropes have often been invoked as a theoretical construction corresponding to the pre-event plasma cavity, which contains the free energy necessary to drive CMEs (e.g., Low 2001;Török and Kliem 2003;Fuller et al. 2008). The pre-event magnetic field supporting the filament can also be well described by highly sheared magnetic arcades crossing the PIL (e.g., Mikić et al. 1988;Steinolfson 1991;Antiochos et al. 1999;Lynch et al. 2008;van der Holst et al. 2009). The left column shows the CME/ICME evolution observed by STEREO A (left) and STEREO B (right) near simultaneously. From top to bottom, the panels display the composite images of EUVI at 304 Å and COR1 showing the nascent CME (indicated by the arrow), combined COR1 and COR2 images of the fully developed CME, and running difference images from HI1 and HI2 when the ICME is far away from the Sun. The crosses mark the locations of the CME leading and trailing edges obtained from the time-elongation map. The positions of the Earth and Venus are labeled as E and V. Right: Time-elongation maps constructed from running difference images of COR2, HI1 and HI2 along the ecliptic plane for STEREO A (upper) and B (lower). The arrows indicate two tracks associated with the CME. The vertical dashed lines show the MC interval observed at Wind, and the horizontal dashed line marks the elongation angle of the Earth. Adapted from Liu et al. (2010a) Upon eruption, these arcades invariably neck off and reconnect to form erupting magnetic flux ropes attached to the Sun at both ends. Regardless of the simulated CME initiation process, the magnetic structure expelled from the corona is almost universally considered to be a twisted structure that can be characterized as a flux rope. It is here that we begin our evaluation of the physical processes that guide the evolution of CMEs from the low corona to interplanetary space. Flux rope models have been shown to self-consistently reproduce many observed properties of CMEs, including the three-part density structure (e.g., Gibson and Low 1998;Manchester et al. 2004b;). In a similar vein, the graduated cylindrical shell (GCS) model (Thernisien et al. 2009;Vourlidas et al. 2011;Colaninno and Vourlidas 2015) provides a geometric representation of the CME cavity that is consistent with an idealized flux rope. Parameters for the model are determined by fitting two or three nearly simultaneous multi-viewpoint observations derived from STEREO Sun Earth Connection Coronal and Heliospheric Investigation (SECCHI), SOHO/LASCO C2 and C3 Fig. 2 The MC observed at Wind corresponding to the 12 December 2008 CME (after Liu et al. 2010a). From top to bottom, the panels show the proton density, bulk speed, proton temperature, and magnetic field strength and components, respectively. The shaded region indicates the MC interval, and the hatched area shows the predicted arrival times (with uncertainties) of the ICME leading and trailing edges. The horizontal lines mark the corresponding predicted velocities at 1 AU. The dotted line denotes the expected proton temperature from the observed speed and Solar Dynamics Observatory (SDO) Atmospheric Imaging Assembly (AIA). The GCS model is capable of describing a wide range of quantities including the bulk velocity mass distribution and three-dimensional trajectory of the CME (e.g., Shi et al. 2015). Figure 3 shows application of the GCS model to the 16 June 2010 CME event, which is a relatively slow CME occurring with the eruption of a quiescent filament. Here, coronagraph images from STEREO/SECCHI and LASCO coronagraphs are shown along with GCS model represented as green circular lines describing the location of a three-dimensional flux rope in the shape of a crescent as seen in the COR2-A field of view of Fig. 3. The model in this case shows evidence for a variety of physical processes we will discuss in following sections, including super-radial expansion and rotation within the first 5 R . In several examples (e.g., Liu et al. 2010b;Vourlidas et al. 2011;Nieves-Chinchilla et al. 2012;Isavnin et al. 2014;Shi et al. 2015;Schmidt et al. 2016), the GCS model has shown how the early 3-D evolution of CMEs can be well described as magnetic flux ropes that are prone to both deflection and rotation.
Flux ropes ejected from the solar corona during CMEs may travel through interplanetary space largely intact, and careful examination suggests that they can be connected to the magnetic structures observed at 1 AU (e.g., Yurchyshyn et al. 2007;Démoulin 2008;Möstl et al. 2008;Liu et al. 2010aHoward and DeForest 2012;Manchester et al. 2014a;. The magnetic fields associated with ICMEs often retain a coherent structure resembling a flux rope. Referred to as magnetic clouds (MCs) (Burlaga 1981(Burlaga , 1988Lepping et al. 1990;Burlaga et al. 1995), these ICMEs are characterized by high magnetic field strength with a smooth rotation (e.g., south to north or east to west) of the field direction, low ion temperature, low plasma beta (typically less than 0.1). The rotation of the field is suggestive of a flux rope geometry (e.g., Lepping et al. 1990;Hu and Sonnerup 2002;Liu et al. 2008a), while the occasional presence of counter-streaming electrons suggest the magnetic field remains attached to the Sun at both ends (e.g., Gosling et al. 2001). The charge state composition of MCs shows elevated ionization states, which are suggestive of flare heated material being ejected with the CME (e.g., Neugebauer and Goldstein 1997;Lepri and Zurbuchen 2004;Zurbuchen and Richardson 2006). Magnetic clouds are also distinguished by their large-scale, and their passage past the Earth that may last 7 to 48 hours, with an average of approximately 21 hours (Lepping et al. 2006). The time and speed indicates that near Earth, the average radial width of MCs is about 0.22 AU and about 1.33 AU at 10 AU (Liu et al. 2006a). MCs are particularly likely to be detected in the near-Earth solar wind when a CME originates within 30°of disk center (Gopalswamy et al. 2001a), indicative of the longitudinal size of these truly global-scale heliospheric disturbances. It should be noted that apart from size, the signatures of ICMEs do not usually occur simultaneously and few ICMEs have all of them.
The relative proportion of ICMEs that appear as MC events has historically shown great variation through the solar cycle. At solar minimum, nearly all ICMEs at Earth can be identified as MCs (Cane and Richardson 2003;Richardson and Cane 2004b), while at solar maximum only ≈ 15% can be. Averaged over the solar cycle, MCs comprised ≈ 30% of ICMEs (Gosling 1990). The cycle dependence reflects several different aspects of CMEs including their place of origin, orientation and mutual interaction. At solar minimum, a greater majority of CMEs originate from streamer blowouts and quiescent filament eruptions. These eruptions are more prone to produce slow CMEs, which are less likely to interact with one another owing to lower eruption rates. Also at solar minimum, CMEs erupt more often at low latitude, providing a greater opportunity for near-central impacts for observing spacecraft, which are more likely to register MC signatures. In contrast, at solar maximum, high latitude eruptions result in off-center ICME in situ measurements which are less likely to register the field line rotations of a flux rope. Also at solar maximum, more CMEs originate from active regions where high eruption rates are prone to produce complex interacting ICMEs, as will be discussed in Sect. 9. However, more recent analysis suggest that nearly all ICMEs have flux rope structure (e.g., Owens et al. 2005;, and that even plasma-dense ejecta may be fit with flux ropes (Marubashi et al. 2015). ICME-related signatures can also continue well beyond the MC boundaries (e.g., Richardson and Cane 2010b;Kilpua et al. 2013). Manchester and Zurbuchen (2006) found that open field lines deflected around the ejected flux rope can have plasma and magnetic characteristics of a MC at latitudes extending beyond the ejected flux rope.

Characteristics and Causes of CME Deflection
CME deflection is the departure from a radial trajectory that commonly occurs with significant in-course changes in direction (e.g., Vandas et al. 1996;Wang et al. 2004;Lugaz et al. 2011;Kay et al. 2013Kay et al. , 2016. Figure 4 shows a clear example of deflection for the 2 November 2008 event where the CME's change in latitude is obvious when comparing the STEREO-B COR1 to COR2 images. A survey by Isavnin et al. (2014) found a maximum total CME deflection (from Sun to Earth) in latitude was 49 • and almost 30 • in longitude. These deflections can be attributed to two primary causes: First, magnetic forces produced by the background corona (e.g., MacQueen et al. 1986;Kilpua et al. 2009;, including the active region of origin . Second, the background solar wind flow pattern can inhibit the latitudinal expansion of the CME in the corona (e.g., ) and the wind also interacts with ICMEs farther out in the heliosphere (e.g., Wang et al. 2004). Magnetic forces control the deflection low in the corona, while the importance of kinematic interactions increases at larger heliospheric distances. In the case of kinematic interactions, the deflecting forces are also ultimately magnetic in nature (see e.g., discussion . This CME deflected quickly to the ecliptic and was observed at STEREO-A as a well-defined magnetic cloud a few days later. The event is studied in detail in Kilpua et al. (2009). Panels (a)-(c) adapted from Kilpua et al. (2009) in Isavnin et al. 2014). This deflection occurs because of ICME interactions with the ambient solar wind that pileup plasma and drape magnetic field at the edges of the CME ejecta, as will be discussed in detail in Sects. 5 and 6.  argue that CMEs tend to deflect toward the region of the lower magnetic energy density by the combined effect of the magnetic pressure and tension forces. The authors presented a theoretical method to account for this effect, which was used statistically by  to confirm the deflection toward the magnetic energy minimum.  demonstrated this deflection property using the Forecasting CMEs Altered Trajectory (ForeCAT) model (see also Kay et al. 2013). This tool propagates CMEs using a drag-based empirical model that takes into account magnetic forces as well as CME expansion. The results of  show a wide range of deflection and also illustrate circumstances when forces are insufficiently strong to fully deflect CMEs toward the energy minimum.
Consider CME deflection from the corona/heliosphere system in the minimum energy state, which has open flux extending from coronal holes separated by a streamer belt with a heliospheric current sheet (HCS) extension. The magnetic field in coronal holes is typically stronger than that found in the surrounding closed flux systems, which provides a magnetic gradient that readily deflects CMEs. The importance of coronal holes for deflection is featured in many studies (e.g., . Gopalswamy et al. (2009a) and Mohamed et al. (2012) estimated the magnitude and direction of the resultant magnetic force exerted by all coronal holes present on the solar disk during the time of the CME eruption (defined as the CHIP parameter in Mohamed et al. 2012) and compared this with CME trajectories and in situ observations. It was found that CMEs tend to move away from the coronal holes and the CMEs that erupt close to disk center but had large CHIP parameters produced generally complex ICMEs or driverless shocks in the near-Earth solar wind. In addition, several recent studies have demonstrated that CMEs are also deflected by strong magnetic fields in the CME source active region (e.g., . Figure 5 shows examples of deflection in both circumstances: CMEs originating from a lowlatitude active region and from a coronal hole boundary depicted in Fig. 5a, b, respectively. These results, calculated with ForeCAT, indicate that deflections in latitude and longitude may reach 30 • to 40 • , respectively, with the magnitude being inversely related to CME speed and mass. The work of Lugaz et al. (2011) on CME deflection presents a slightly different configuration, that of an anemone active region, where the Lorentz force drives the deflection.   shown for CMEs originating from two locations: (a) low-latitude active region and (b) coronal hole boundary. Here, the dots represent the deflected location of the CMEs where the size of the dots are proportional to mass ranging from 10 14 to 10 15 g, and the speed of the CMEs ranging from 300 to 1500 km s −1 is indicated by color. It is shown that deflection in latitude and longitude may reach 30 • to 40 • , respectively, with deflection being inversely related to CME speed and mass The CME deflection in latitude is constrained by the location of the streamer belt/HCS and this deflection occurs predominantly close to the Sun in the neighborhood of the streamers. However, longitudinal deflections are largely controlled by kinematic interactions (e.g., Wang et al. 2004) that occur at larger distances in the corona and heliosphere. The main source of longitudinal deflection is caused by interaction of ICMEs with the Parker spiral structured solar wind that occurs when there is sufficient speed difference with the ICME (Wang et al. 2004). The resulting magnetic forces and the direction of the deflection are different for slow and fast ICME populations: Slow ICMEs are deflected westward when they are pushed from behind by the faster wind, while fast ICMEs deflect eastward as they are decelerated by the slower wind head. We note however that ambient Parker spiral field by itself may not be able to deflect a ICME by more than a few degrees as even for a slow ICME its kinetic energy density is about two orders of magnitude higher than the magnetic energy density of the Parker field. Finally, interactions between multiple CMEs/ICMEs, in particular when they collide, can also cause longitudinal deflections (e.g., Shen et al. 2012;.

Deflection Dependence on the Background Corona
Based on the above-described studies, the rate and amount of CME/ICME deflection is controlled by the strength and distribution of the background magnetic field, and the mass, size, and speed of the CME/ICME relative to the solar wind. Hence, both the global configuration of the Sun's magnetic field and intrinsic CME properties have crucial importance on the degree and direction of deflection, which we now further quantify. For example, Xie et al. (2009) showed that during solar minimum slow CMEs deflected toward the ecliptic and the streamer belt by strong polar magnetic fields, while fast CMEs are deflected less, sometimes also away from the streamer belt away from the streamer belt, confirming earlier results by MacQueen et al. (1986). Similarly Wang et al. (2011) found 62% of CMEs deflected towards equator with an average deflection angle of 22 • , and only 5% of CMEs deflected towards poles, with an average deflection angle of 16 • . Also a case study by Kilpua et al. (2009) found that slower and wider CMEs deflected toward the equator while the faster and narrower CME propagated radially from its source active region. It was suggested that slow and wider CMEs cannot penetrate through the overlying coronal fields, but are channeled toward the streamer belt. These findings are consistent with the ForeCAT model  showing that slow, wide and low-mass CMEs deflect most as shown in Fig. 5.
The rate of CME deflection is clearly fastest close to the Sun where magnetic forces dominate. The analysis of 14 CMEs by Isavnin et al. (2014) showed that about 60% of the total evolution from the Sun to Earth orbit takes place in the corona, i.e., within about first 20-30 R from the Sun, in particular, the rate being highest in the low corona (< 5 R ).  came to similar conclusions using the ForeCAT model; the majority of CME deflection occurs within 10 R from the Sun.  also pointed out that deflections are confined closer to the Sun in the case of stronger background magnetic fields. For strong fields studied in their paper, the deflections occurred primarily below 2 R from the solar surface. This result is consistent with  who found a positive correlation with the rate of the deflection and the strength of the magnetic energy density gradient. The intrinsic magnetic polarity of the CME flux rope relative to the background coronal magnetic field determines whether and where the magnetic reconnection can occur and consequently deflects the CME (e.g., see simulation work by Chané et al. 2005;Zuccarello et al. 2012a;Zhou and Feng 2013). If the CMEs intrinsic field is parallel to the ambient field, the CME deflects equator-ward while when the fields are antiparallel, the CME is likely to deflect poleward (Zhou and Feng 2013).

Characteristics and Examples of CME/ICME Deflection
There are several examples in the literature where deflections of CMEs have changed the expected ICME impact at Earth and the expected geoeffectivity. For example, Zhou et al. (2006) showed that intrinsically high-latitude CMEs can drive strong storms; nearly 30% of the Earth-encountered ICMEs associated with high-latitude polar crown filament disappearances investigated in their study caused at least moderate space weather effects. CMEs have a strong tendency to deflect from high latitudes toward the equator in particular near solar minimum (e.g., Kilpua et al. 2009;Byrne et al. 2010;Isavnin et al. 2014). This behavior is expected, as at this time, the global magnetic field of the Sun is relatively close to a dipole field and two large polar coronal holes dominate the field structure. These coronal holes can effectively guide CMEs toward the ecliptic, consistent also with the suggestion by  as the minimum energy region streamer belt/HCS is relatively flat and confined close to the equator during slow solar activity period. Near solar maximum the configuration of the Sun's global magnetic field and the distribution of magnetic energy density is more complex than near solar minimum and CMEs deflect less and more randomly. Poleward deflection occurs mainly during the times of solar maximum when large low-latitude polar coronal holes are present.
As discussed in the beginning of this section, longitudinal deflection may cause CMEs/ICMEs to deviate from the Sun-Earth line or toward it. According to study by Gopalswamy et al. (2009b), almost ten percent of large geospace storms are caused by CMEs that originate close to the limb of the Sun. In such cases, the ICME sheath is typically the primary driver of the storm (see also Huttunen et al. 2002), but cases have been reported where clear ejecta signatures upstream of the Earth and geomagnetic activity have been associated with limb CMEs (e.g., Schwenn et al. 2005;Cid et al. 2012;. For example, the CME studied by  was initially heading toward the STEREO-B, but it deflected in the heliosphere and arrived at the Earth instead, which was located about 35°away from the STEREO-B at that time. An opposite case is featured e.g., in ,  where the authors studied a CME Fig. 6 Difference between the observed and calculated transit times (O − C) of 55 fast CMEs (v > 500 km s −1 ) presented as a function of central meridian distance (CMD). The calculated transit times are based on the Drag-Based Model (DBM) using the solar wind speed of w = 450 km s −1 and the drag parameter Γ = 0.2 (for details see Vršnak et al. 2013) that originated from an active region near the solar disk center and significant geomagnetic response was expected. However, multi-spacecraft observations and modeling demonstrated that this CME deflected almost 40°in longitude and caused only minimal space weather effects. A similar event illustrating the magnitude of CME/ICME longitudinal deflection as reported by . Gopalswamy et al. (2009a) also demonstrated that some halo CMEs that erupt near disk center appear to be associated with driverless shocks near Earth orbit (i.e., an interplanetary shock that is not followed by the discernible driver or a slow-fast stream interaction region) due to significant deflection in longitude away from the Sun-Earth line. In these cases the ICME shock and sheath is detected in situ at Earth, but the driving ejecta is missed completely.
Finally, we illustrate the eastward-deflection effect on fast CMEs by presenting in Fig. 6 the dependence of the difference between the observed and calculated transit times (O − C) on the central meridian distance (CMD) of their source region. The sample contains 55 fast CMEs, whose transit times are calculated using the Drag-Based Model (DBM) by using the solar wind speed of w = 450 km s −1 and the drag parameter Γ = 0.2 (for details see Vršnak et al. 2013). Note that O − C > 0 means that the ICME arrived later than expected for the chosen DBM parameters (in the presented case the presumed solar wind speed is too high, causing underestimation of transit times and leading to the average O − C of 9.6 h). The scatter-plot includes the quadratic least squares fit, characterized by the correlation coefficient of cc = 0.24 and the F-test statistical significance of P > 99%. The fit shows that O − C values are on average larger for CMEs launched closer to the solar limb, meaning that such CMEs are slower than those launched close to the disc center, perhaps indicating that the flank speed is lower than at the nose of the CME. Furthermore, the effect is larger for the eastern-hemisphere events than for western ones, demonstrating the eastward deflection of fast events.

CME Rotation
Erupting prominences (or filaments) frequently exhibit a rotation about their rise direction as they ascend in the corona, which leads to a deviation from their original orientation on the solar surface. Strongly rotating filaments display a characteristic "inverse γ " shape, which develops when the legs of the filament cross along the line of sight of the observer (Fig. 7a, b). The direction of the rotation is determined by the sign of helicity of the source region, though there seem to be exceptions (Muglach et al. 2009). When viewed from above,  Török and Kliem (2005). Colored field lines outline the core of a kink-unstable flux rope with positive helicity (right-handed twist), surrounded by green ambient field lines. (c): Top view on the simulation. The flux-rope core has rotated clockwise by almost 90°with respect to its initial orientation. (d): Fully erupting prominence on 9 April 2008 observed by STEREO Ahead at 10:55 and 11:25 UT. The colored strands were used for a 3D reconstruction. The right panels show field lines from a simulation of the event . (e): Prominence rotation vs. heliocentric height of the prominence's leading edge. Panels (a)-(c) adapted from Green et al. (2007). Panels (d)-(e) adapted from  clockwise (anti-clockwise) rotation is typically observed for filaments erupting from source regions with positive (negative) helicity (Green et al. 2007;Fig. 7c). This suggests the conversion of twist into writhe in a kink-unstable magnetic flux rope as a possible mechanism of the rotation . Figure 7 shows a confined (or failed) eruption, but significant rotations occur also in cases where the filament fully erupts as part of a CME (Fig. 7d). While filament material outlines only a small fraction of a CME, it is believed to be located close to the axis of the CME flux rope. Thus, the observed rotation of a filament implies that the whole CME flux rope rotates as well, even though outer flux surfaces of the rope may rotate at a smaller rate.
Since the magnetic orientation of an ICME upon arrival at Earth is one of the main parameters that determine its geo-effectiveness (as discussed by Kilpua et al. in this issue), it is important to understand the physical mechanisms that cause the rotation, and to quantify the total rotation that CMEs/ICMEs undergo during their travel to Earth, particularly in light of the desire to develop methods to reliably predict the sign and magnitude of B z at 1 AU. However, obtaining reliable measurements of the total rotation directly from the observations is difficult for several reasons. First, complete observational coverage of the CME and ICME propagation from Sun to Earth is not always available. Second, rotation is hard to recognize in coronagraph images, especially if the CME is observed above the solar limb and no filament is present, though measurements could be obtained for halo CMEs using flux-rope fitting and three-dimensional (3D) reconstruction techniques (e.g., Yurchyshyn et al. 2007;Liu et al. 2010b;Vourlidas et al. 2011). Third, most of the rotation may take place before the CME enters a coronagraphs' field of view, which requires EUV observations of the associated filament eruption, with the filament material exhibiting a coherent shape for a sufficiently long time. Even if the latter condition is met, reliable measurements of the rotation can only be obtained if the filament either erupts toward the observer or if 3D reconstructions of its shape can be made.
A rare example of direct rotation measurements is shown in Fig. 7d and e. Using observations from both STEREO spacecraft,  were able to measure the rotation of the prominence that erupted during the "Cartwheel" CME on 9 April 2008, which showed a strong rotation of about 115°up to a height of about 1.5 R above the surface. Afterward, the rotation direction seemed to reverse slowly. In a similar investigation of an erupting quiescent polar-crown prominence, Thompson (2011) obtained a rotation of at least 90°. Also in this case, the rotation saturated before the prominence became invisible. These results suggest that CME rotation saturates (or even reverses) already in the low corona. In contrast, the studies by, e.g., Yurchyshyn et al. (2007), Liu et al. (2010b), Lynch et al. (2010), and Vourlidas et al. (2011) suggest that CME/ICMEs still exhibit a significant rotation at larger coronal heights and in interplanetary space.
Given the various observational limitations mentioned above, the total rotation of CMEs is typically estimated by comparing the orientation of the pre-eruptive structure on the Sun with the orientation of the axis of the magnetic cloud at 1 AU. Pre-eruptive orientation is inferred from the polarity inversion line of the source region or from filament, sigmoid or coronal loop observations, while the orientation at 1 AU is obtained by fitting 3D fluxrope models to the 1D in situ data, given that the circumstances allow such fitting. Such comparisons of orientation, albeit hampered by some uncertainty, suggest that in many, if not most, cases the total rotation of the CME/ICME remains relatively small. However, as reviewed by Démoulin (2008), cases with total rotations larger than 30°are not uncommon, and very large values of about 120°or more have been reported (e.g., Rust et al. 2005;Dasso et al. 2007;Liu et al. 2008bIsavnin et al. 2014;Vemareddy et al. 2016).
Numerical simulations and theoretical considerations have been employed to understand the physical causes of CME rotation and to explain the wide range of rotations observed. Török and Kliem (2003) showed that erupting flux ropes undergo a significant rotation (of more than 90°) when driven by vortex flows at their foot points (used to model rotating sunspots), and suggested that the strong rotation occurred due to the development of the ideal MHD kink instability. Starting from analytical flux-rope models, Fan and Gibson (2004) and Török and Kliem (2005) reported similarly large rotations (up to 120°) as a result of the same instability. Using a similar flux-rope model, Fan (2016) obtained an even larger rotation of almost 180°, i.e., a full reversal of the magnetic field vector at the front of the flux rope, in a recent simulation of the 13 December 2006 event. Isenberg and Forbes (2007) suggested that the presence of an external shear field surrounding a pre-eruptive flux rope (i.e., of an ambient field component pointing along the axis of the rope) provides a different mechanism for the rotation of a flux rope, once the rope leaves its equilibrium state and rises in the corona. The Lorentz forces that arise from the interaction of the flux-rope current with the external shear field cause a rotation that acts in the same direction as the kink instability.  confirmed the suggestion by Isenberg and Forbes in a series of MHD simulations, in which they varied the twist within a flux rope, the strength of the external shear field, and the decrease of the ambient potential field with height (by changing the distance between the photospheric polarities that generate that field). Changing these parameters, they obtained flux-rope rotations in a wide range of about (35-170)°. Their study revealed a number of interesting results: (1) For small distances between the polarities, the dominant contribution to the rotation comes from the external shear, even for strongly kink-unstable flux ropes; (2) a moderate rotation in erupting flux ropes is always present due to the conversion of flux-rope twist into writhe (Török et al. 2010), even in the absence of a shear field and of the kink instability; (3) the amount of rotation depends on the slope of the ambient field: the slower the field drops with height, the more the rope rotates. This effect is, however, relatively weak as long as the distance between the polarities is smaller or comparable to the distance of the flux-rope foot points (as it is typically the case on the Sun), but becomes significant if the former distance dominates.
In contrast to these simulations, in which the eruption starts from a fully-developed flux rope, Lynch et al. (2009) modeled CME rotation starting from a sheared arcade, using the breakout model (Antiochos et al. 1999;Lynch et al. 2008). Initially, the arcade rose slowly, without a clear sign of rotation. When it reached a height of about 2 (heliocentric) R , flare reconnection set in and a twisted flux rope was formed. From that point on, the rope rotated at a constant rate of about 30°R −1 , and reached a rotation angle of about 50°at a height of 3.5 R . The twist in the rope was low, clearly below the threshold of the kink instability. The authors explained the rotation by the effect of the tension force associated with the sigmoidal shape of the erupting field lines, which acted to straighten out those field lines once the eruption was underway. As the field lines became straight, the rotation still continued due to the angular momentum imparted in the early phase of the eruption. Lynch et al. (2009) did not report how the rotation evolved beyond 3.5 R , so it is not clear at which point a saturation of the rotation may have occurred. Their simulation results agree qualitatively very well with the case observed by Vourlidas et al. (2011), where apparently little or no rotation of the CME occurred below 2 R , and an almost constant rotation (at a rate twice as large as in the simulation) was seen above that height.
Other mechanisms that have been suggested based on observational or numerical studies are the straightening of an initially strong S-shape of a flux rope during its eruption (Török et al. 2010;); reconnection of an erupting flux rope with the ambient magnetic field (Jacobs et al. 2009;Shiota et al. 2010;Cohen et al. 2010;Thompson 2011;Lugaz et al. 2011); and the alignment of the CME flux rope with the heliospheric current sheet (e.g., Yurchyshyn 2008).
Overall, these studies reveal a substantial number of mechanisms that can cause CMEs to rotate about their direction of propagation. Moreover, as discussed in detail in , several mechanisms can contribute simultaneously or successively to the rotation in a complicated, parameter-dependent manner, making the quantitative prediction of CME rotation a very challenging task. Specifically, predictions of the sign of B z at 1 AU merely based on the pre-eruptive orientation of the ejecta have to be taken with care, given the significant number of events that exhibit a large total rotation. Lynch et al. (2009) suggested that the amount of total rotation my be predicted by the degree of "sigmoidality" of the pre-eruptive field, but this approach needs to be tested using observations, and it will likely underestimate the total rotation if other mechanisms such as the kink instability are involved. If most of the rotation typically occurs low in the corona, as suggested by some observations and simulations, then measurements of the rotation of erupting filaments may be utilized to improve (shorter-term) B z predictions, but the practicability of such an approach needs to be tested as well. Attempts to systematically predict the rotation (and deflection) of CMEs in the corona, based on the properties of the background magnetic field and using analytical flux-rope models, have been recently developed and tested with two observed events by Kay et al. ( , 2016.

A Brief Overview of CME Kinematics
Here, we examine the details of when, where, and how CMEs/ICMEs accelerate/decelerate in interplanetary space and quantify the kinematics behavior. To that end, observations by the STEREO spacecraft have provided unprecedented opportunities to observed kinematic behavior of ICMEs with multiple views that enable accurate measurements of ICMEs over a large distance. Liu et al. (2010a) have developed a triangulation technique to determine CME/ICME Sun-to-Earth kinematics based on the wide-angle imaging observations from STEREO. The technique initially assumes a relatively compact CME structure simultaneously seen by the two spacecraft. It has no free parameters and can give CME/ICME kinematics (both propagation direction and radial velocity) as a function of distance from the Sun continuously out to 1 AU. This capability is key to probing CME propagation and interaction with the inner heliosphere. Later,  and Liu et al. (2010b) realize that the same idea can be applied by assuming CME geometry as a spherical front attached to the Sun. In this case, what is seen by a spacecraft is the segment tangent to the line of sight. The triangulation concept has proven to be a useful tool for determining CME Sun-to-Earth kinematics and connecting imaging observations with in situ signatures (e.g., Liu et al. 2010aMishra and Srivastava 2013). Details of when, where, and how CMEs accelerate/decelerate in interplanetary space can be quantified using the kinematics derived with the triangulation method. Figure 8 shows the comparison of Sun-to-Earth propagation profiles characteristic of slow, fast and intermediate-speed CMEs. The 7 March 2012 CME, with a peak speed of more than 2000 km s −1 shows a typical three-phase speed profile developed by : an impulsive acceleration for a CME up to 10-15 R , followed by a rapid deceleration for the ICME out to about 50 R , and thereafter a nearly constant speed (or gradual deceleration). The predicted speed at the Earth overestimates the observed speed by about 250 km s −1 , which may be partly due to the large longitudinal separation between the two STEREO spacecraft (227 • , i.e., both observing the CME from behind the Sun). Note that the observed speed at the Earth is the average solar wind speed in the sheath between the shock and ejecta, which is usually a little smaller than the shock speed at 1 AU. Hence, a better agreement may be achieved if the shock speed (rather than the sheath speed) is used for the comparison. The 12 December 2008 CME with a peak speed of about 700 km s −1 has a similar profile but different cessation distances for the acceleration and deceleration. Compared with the 7 March 2012 CME, its speed first increased with a lower rate up to about 20 R , then decreased out to 80-90 R , and thereafter became roughly constant, similar to the evolution described by  for a CME of similar speed. The cessation distance of ICME rapid deceleration (about 50 and 85 R for the 7 March 2012 and 12 December 2008 CMEs, respectively) is much shorter than the average cessation distance of 0.76 AU (163 R ) inferred indirectly by Gopalswamy et al. (2001b). The Sun-to-Earth propagation profile of the 16 June 2010 CME/ICME exhibited only two phases: an acceleration with an even slower rate up to 25-30 R followed by a nearly invariant speed at about 390 km s −1 . The predicted speeds at the Earth for the latter two cases are well confirmed by the in situ measurements at 1 AU; the longitudinal separation between the two STEREO spacecraft is 86.3 • for the 12 December 2008 CME and 143.6 • for the 16 June 2010 event.
In examining the question of CME/ICME kinematics and momentum transfer between the ejecta and the solar wind, there are many aspects that must be considered: First, the ambient solar wind is often highly variable in time and space (e.g., Kilpua et al. 2012;. Second, that the Lorentz force driving the CME can sometimes significantly contribute to the CME dynamics up to large distances from the Sun (e.g., Vršnak et al. 2004;). Third, CME-CME interactions including the affects of CME preconditioning the solar wind, (e.g., Lugaz et al. 2005b;Maričić et al. 2014;, which will be taken up in Sect. 9. All of these phenomena affect the CME/ICME propagation and significantly contribute to the uncertainties in predicting the arrival time and impact speed of ICMEs at Earth. Generally, the CME/ICME propagation can be divided into several phases. The CME initiation phase is most often characterized by swelling and slow rising motion of the preeruptive structure, which is usually interpreted as an evolution through a series of quasiequilibrium states (see, e.g., Vršnak 2008, and references therein). When the slowly evolving structure reaches an unstable state (so called "loss of equilibrium" mechanism of Forbes and Isenberg 1991), the structure starts to accelerate, driven by some form of ideal MHD instabilities such as the kink or torus modes (e.g., and references therein Kliem 2003, 2005;Kliem and Török 2006;Démoulin and Aulanier 2010;Olmedo and Zhang 2010;Török et al. 2010). Then follows the take-off stage, which is most often characterized by accelerations on the order of 100 m s −2 , and lasts for ∼ 1 h, so the majority of CMEs achieve velocities in the range 100-1000 km s −1 (e.g., Bein et al. 2011;). However, sometimes the take-off phase is characterized by an extremely impulsive acceleration, achieving peak values on the order of 10 km s −2 , but lasting only for several minutes (e.g., Bein et al. 2011). In such events the gradual pre-eruptive evolution is frequently not observed. In fact, in some events a confined-flare type of process leads to fast reconfiguration of the pre-eruptive structure into an unstable configuration and erupts immediately after being formed (Aurass et al. 1999).
On the other side of the "spectrum" of event types are very gradual events, characterized by weak accelerations that result in velocities on the order of 100 km s −1 .
The take-off phase imparts the initial velocity, which can be categorized as slow (speeds below 400 km s −1 ) intermediate (speeds between 400 and 1000 km s −1 ) and fast (speeds above 1000 km s −1 ), and it is this speed that largely determines the kinematical evolution in the ensuing "main acceleration" phase. Most generally, in the upper corona and heliosphere, CMEs that are slower than the ambient solar wind continuously accelerate, whereas those that are faster than the wind decelerate (Lindsay et al. 1999;Gopalswamy et al. 2001b;Jones et al. 2007) such that CMEs tend to approach the speed of the ambient solar wind. This acceleration/deceleration tendency is already observable in the upper corona (Moon et al. 2002;Vršnak et al. 2004) and was found by comparing the distribution of coronagraphic CME velocities with the ICME speeds measured in situ at 1 AU (Gopalswamy et al. , 2001b. In the heliospheric range, the velocity profile was directly confirmed by various types of measurements, e.g., by white-light heliospheric imagers such as those made by Coriolis instrument (e.g., Reiner et al. 2005b,a;Tappin 2006;) onboard the Solar Mass Ejection Imager (SMEI) (Eyles et al. 2003;Jackson et al. 2004) and STEREO-HI (e.g., , by tracking interplanetary radio type II bursts (e.g., Reiner et al. 2005a,b;Liu et al. 2008b), by analyzing radio interplanetary scintillation (e.g., Manoharan et al. 2001;Manoharan 2006), or by comparing coronagraphic and 1 AU CME/ICME speeds (e.g., Manoharan and Mujiber Rahman 2011).

Empirical Models of CME/ICME Kinematics
In parallel with progress of numerical propagation models, a number of empirical methods and analytical physics-based models were developed. Although, quite clearly, the future of space weather forecasting lies in numerical modeling, the performance of empirical and analytical methods is currently comparable to, or even slightly better than that of numerical methods. On the other hand, these methods are easy to handle and could be promptly adjusted to the new data incoming in the course of the CME/ICME propagation, i.e., the predictions could be easily refreshed instantaneously with receiving the newest input data. Furthermore, the empirical and physics-based analytical models educe our comprehension of the CME/ICME dynamics, and thus, could help in advancing the numerical methods. In this section, we review a physical background of this approach to the space weather forecasting, which is generally based on empirical studies of the CME/ICME kinematics and/or theoretical considerations of their dynamics.
Results of various statistical studies of the CME/ICME Sun-Earth transit times, reflecting the described general characteristics of the CME/ICME kinematics, can be used to establish empirical methods for predicting the ICME arrival. The simplest one, based on a statistical study of a sample of geoeffective ICMEs, was proposed by Brueckner et al. (1998), saying: "In many cases, the travel time between the explosion on the Sun and the maximum geomagnetic activity is about 80 hours." To explain this, so-called "Brückner's 80-hour rule", corresponding to the Sun-Earth average speed of 500 km s −1 , one has to bear in mind that the Sun-Earth travel time for a typical slow solar wind of the velocity 400 km s −1 is about 100 h. On the other hand, most of geoeffective ICMEs are relatively fast, and for the speed of 1000 km s −1 , the transit time is 40 h. Thus, it can be concluded that the 80-hour rule is a consequence of deceleration of fast CMEs in the ambient solar wind.
A somewhat more detailed empirical method, based on the correlation of the Sun-Earth transit times, tt, and the coronagraphic speeds of CMEs, V CME is illustrated in Fig. 9a. A similar plot is presented in Vršnak and Žic (2007), and for additional examples see also In Fig. 9a, the arrival time for a CME of a given coronagraphic speed can be predicted by employing the presented power-law fit. In Fig. 9b, the residuals O − C are presented to show the span of the prediction errors as a function of V CME . Inspecting the graph one finds that the errors can be as large as 40 h and tend to be larger for slower CMEs (however the relative errors are similar for the whole range of V CME ). The standard deviation of the O − C distribution is 18 h. The errors can be somewhat reduced by taking into account the sourceregion location, the CME mass, the associated-flare importance, etc., but not significantly. Note that the same procedure can be applied by using, e.g., the coronal shock velocities inferred from the radio type II burst dynamic spectra.
The scatter-plot in Fig. 9a shows that the power-law fit is characterized by tt ∝ V −0.41 , which is distinctively different from tt ∝ V −1 that would be expected in the case of V CME = const. This means that slow/fast CMEs have shorter/longer transit times than they would have if propagating at constant speed. Again, this can be attributed to the adjustment of the CME/ICME propagation to the ambient solar wind.
The statistical relationships between the arrival time and various CME parameters measured from coronagraphic observations can be used to employ more sophisticated prediction methods. For example, Sudar et al. (2016) applied the "neural network" approach to determine the most probable transit time using the coronagraphic CME speed and the central meridian distance of the source-region as the input parameters. The analysis, involving a sample of 153 CME-ICME pairs, showed that the average tt error is 12 h. The tt (V CME ) dependence showed a typical drag-like trend, with acceleration turning to deceleration at V CME 500 km s −1 , consistent with the "Brückner's 80-hour rule" (note that similar holds for the scatter-plot shown in Fig. 9a). Furthermore, the results clearly demonstrate larger tt for larger CMD, as well as the eastward/westward deflection of fast/slow CMEs. Note that an analogous technique can be applied to forecast the geoeffectiveness of CMEs (e.g., Valach et al. 2009;Uwamahoro et al. 2012;. The acceleration/deceleration tendency of slow/fast CMEs was quantified by Lindsay et al. (1999) and Gopalswamy et al. ( , 2001b by defining a linear relationship between the CME acceleration expressed in m s −2 and the CME speed expressed in km s −1 as a = 2.193 − 0.0054V CME , which can be also expressed as a = −0.0054(V CME − 406). The range over which the deceleration occurs was estimated to r < 0.76 AU. Later on, Manoharan et al. (2004), Manoharan (2006), and Gopalswamy (2009) extended the a(V CME ) relationship to 2nd-degree polynomial forms. The described kinematical forecasting technique, employing various forms of the a(V CME ) relationship, was applied also by, e.g., Manoharan et al. (2004), Gopalswamy et al. (2005), Reiner et al. (2005b), Manoharan andMujiber Rahman (2011), Salas-Matamoros andKlein (2015), where also forecasting of the CME-driven shocks was included.

Physics-Based Kinematic Models
The statistical tendency showing that fast CMEs decelerate, whereas slow CMEs accelerate, as demonstrated by Lindsay et al. (1999), and confirmed by Gopalswamy et al. (2001b), Moon et al. (2002), Vršnak et al. (2004), , indicates that CMEs tend to adjust their velocity to the ambient solar wind, which could be interpreted as a consequence of "aerodynamic" drag (see, e.g., Cargill et al. 1996;Vršnak 2001;Vršnak et al. 2004Cargill 2004;Owens and Cargill 2004;Manoharan 2006). In particular, Vršnak et al. (2004) have shown that the a(V CME ) relationship for the CMEs observed in the LASCO field of view has a quadratic form, providing a strong evidence for the aerodynamic drag effect. Furthermore, they showed that the Lorentz force is decaying rapidly with the height, implying that the drag becomes a dominant force in the heliospheric dynamics of CMEs.
Following observational facts and theoretical considerations, the so called Drag-Based Model (DBM; Vršnak and Gopalswamy 2002;Vršnak et al. , 2013 was developed to describe analytically the dynamics/kinematics of the heliospheric propagation of CMEs. DBM is based on the assumption that the aerodynamic drag is a dominant force that governs the CME propagation in the interplanetary space. Thus, CME/ICME dynamics are described by the equation of motion of the form a ≡r ≡v = −γ (v − w)|v − w|, where a(t), v(t), and r(t) are the instantaneous CME acceleration, speed, and position, w is the ambient solar wind speed (generally depending on r, but often used to be constant). The drag parameter γ defines the drag "effectiveness", usually being expressed as γ = c d Aρ w /M, where c d is the dimensionless drag coefficient, A is the CME cross-sectional area, ρ w is the ambient solar wind density, and M is the total CME mass (for details see, e.g., Cargill 2004;Vršnak et al. 2013).
The described equation of motion, as well as its modifications (e.g., Tappin 2006;Shanmugaraju and Vršnak 2014;Shi et al. 2015), was employed to study various aspects of the CME/ICME heliospheric propagation. For example, Vršnak and Žic (2007) and Shanmugaraju and Vršnak (2014) studied the role of the solar wind speed on the CME/ICME transit times, showing that the wind speed is an essential parameter for the CME/ICME dynamics, which results in "Brückner's 80hour rule".  analyzed the effects of variable/structured solar wind and showed that it significantly affects the CME/ICME kinematics, but also, they demonstrated that in some cases the driving Lorentz force might contribute to CME/ICME dynamics beyond a distance of 100 R .  studied the role of the CME mass, and found that in the range of distances covered by LASCO, the drag affects more strongly the propagation of light CMEs, and inferred that the Lorentz force is usually stronger in more massive events.  and  employed DBM to investigate how the CME-CME interaction affects the kinematics, whereas Wang et al. (2016b) used it to infer the degree of deflection of a CME/ICME event.  and  compared the DBM and ENLIL ICME arrival-time predictions, and found that results of the two models to be similar accurate.
After being proposed in its initial form by Vršnak and Gopalswamy (2002), the formulation of the basic form of DBM, suited for space-weather forecasting, was completed by Vršnak et al. (2013). In this basic option of DBM, it is taken, approximately, that the CME/ICME cross section increases as A ∝ r 2 , the solar wind density decreases as ρ w ∝ r −2 (consistent with the w = const. approximation), and the mass and drag coefficient are constant. Under this approximation one finds γ = const., which simplifies the equation of motion, providing explicit solutions r(t) and v(t). Obviously, the main advantage of so-defined DBM option is its practical use in real-time forecasting of the ICME arrival time and impact speed. However, serious intrinsic drawbacks are related to approximations applied to the DBM: (1) c d does not vary with distance (e.g., Cargill 2004), (2) ICME virtual-mass effect is not included (e.g., Bein et al. 2013;Cargill 2004;Tappin 2006), (3) solar wind properties are taken to be constant in time (e.g., , (4) and (5) the geometry of the ICME leading edge and location of the source region are not considered, and finally, (6) and (7) the roles of ICME over-expansion (e.g., Riley and Crooker 2004) and the ICME-driven shock  are omitted. Some of these drawbacks were removed by Žic et al. (2015), who included CME geometry and source-region location factors, and also took into account more general forms for ρ w (r) and w(r). Finally, it should be noted that there are also various forms of empirical forecasting methods that rely on relationships between CME-driven shocks and CMEs, CME-associated flares, or type II radio bursts. Analytical models also exit as well as hybrid approaches that combine components of empirical and analytical models (e.g., Feng and Zhao 2006;Zhao and Dryer 2014;Zhao and Feng 2015;Zhao et al. 2016, and references therein).

Mass Increase in CMEs/ICMEs
The momentum transfer between ICMEs and the solar wind can be described as "aerodynamic" drag, however this process is strongly affected by the expanding nature of CMEs and ICMEs as shown by Siscoe and Odstrcil (2008). ICMEs tend to expand only slightly faster than the spherically diverging flow of the surrounding solar wind. Consequently, the ejecta has a nearly constant (Sun-centered) angular size while the transverse size increases linearly with distance from the Sun. Siscoe and Odstrcil (2008) found the expansion speed of the ICME flux rope can be greater than the non-radial deflection speed of the solar wind, so that impacted solar wind is unable to move completely around the ICME, and accumulates in the sheath. The deflection flow pattern around an ICME is shown in Fig. 10a, which depicts the non-radial velocity and the trajectory of solar wind parcels ahead of a moderately fast simulated ICME. The impacted solar wind is unable to move around the flux rope, and as a result, the plasma is swept up in the ICME sheath, increasing the total mass of the disturbance.
This snow-plow effect has been difficult to discern in near-sun coronagraph observations because the swept up mass increase is difficult to distinguish from the sustained outflow of dense plasma with the CME (Bein et al. 2013). However, far from the Sun, the  Siscoe and Odstrcil (2008). Color contours show the non-radial velocity while white lines show the trajectory of solar wind parcels ahead of the ICME. The azimuthal deflection speed is insufficient to allow the impacted solar wind to move around the flux rope, and as a result the plasma is swept up in the ICME sheath. (b): plots of the mass and velocity of a simulated ICME as a function of distance from the Sun taken from Lugaz et al. (2005a). The solid line shows the mass of plasma above ambient density, the dashed line shows the mass of low density material while the dot-dash line shows speed of the ICME. Mass and speed follow equal and opposite trends indicating momentum conservation as the solar wind is compressed at the shock and then piles up in the sheath region snow-plow affect has recently been measured in STEREO HI observations (DeForest et al. 2013), which confirms prior predictions made with numerical simulations (Manchester et al. 2004b), which established a large mass increase occurring with CME deceleration close to the Sun. Here, the speed of the CME was found to follow an exponential-like decay that asymptotically approaches the speed of the background solar wind, thus verifying the empirical relationship established by Sheeley et al. (1999) (see also Fig. 6). A study by Lugaz et al. (2005a) explicitly calculated this ICME mass from both 3-D structure and from synthetic coronagraph images, the results of which are shown in Fig. 10a. Here, ICME mass and velocity are plotted as functions of distance from the Sun. The CME/ICME mass increased by more than a factor of four from the low corona to 1 AU with the majority (factor of three) of the increase occurring within 30 R . This plasma is compressed and accelerated at the CME-driven shock where a significant amount accumulates in the sheath region between the flux rope and the shock as shown in Fig. 11. As a result of momentum conservation, ICME velocity follows an equal and opposite decrease in magnitude.  found similar mass increases in their simulation of the 28 October 2003 CME event. In this case, the mass is derived from analysis of synthetic LASCO C2 and C3 coronagraph images, which showed the mass more than doubled to over 4 × 10 16 g. Analysis of ICME data from SMEI have found similar mass increases, and in the case of the 28 October 2003 event, the mass reached the extraordinary level of 13.6 × 10 16 g when the ICME reached 1 AU as reported by . Much of the mass increase is related to the extreme nature of this event, which reached Earth in approximately 20 hours.

CME-Driven Shocks
High initial speed in conjunction with CME/ICME expansion drive fast-mode waves into the surrounding medium, which may steepen into shocks. These disturbances can propagate far beyond the expanding magnetic driver and are typically observed as forward shocks proceeding ejecta (e.g., Reiner et al. 1998;Wood et al. 2012;Temmer and Nitta 2015;Liu et al. 2017). CME-driven shocks were originally identified in coronagraph images as faint arcs observed at the outer edges of fast CMEs. An early example is provided by Sime and Hundhausen (1987), who observed a bright loop at the front of a fast CME identified as a shock from its high speed (1070 km s −1 ), absence of deflections preceding the loop, as well as the absence of stationary legs, and continuously expanding arc. Perhaps the most compelling observational evidence for shocks appearing in LASCO images is presented in Raymond et al. (2000), Mancuso et al. (2002), and . In the case of Raymond et al. (2000) and Mancuso et al. (2002), shocks were observed simultaneously in the low corona (r < 3 R ) by LASCO, the Ultra-violet Coronagraph Spectrometer (UVCS), and associated type II radio bursts. UVCS gave clear spectroscopic evidence for the presence of shock fronts, while type II radio bursts indicated the presence of shock-accelerated electrons. A connection between type II radio bursts and CME-driven shocks has been established by many studies (e.g., Reiner et al. 1998;Liu et al. 2009Reames 2013;Nitta et al. 2014; MHD simulations also provide predictions of the appearance of CME-driven shocks in a variety of observing modes. Pagano et al. (2008) modeled the spectral line signatures of CME-driven shocks as they appear in UVCS, and ,  simulated the appearance of coronal shocks in synthetic coronagraph images derived from simulations of specific CME events. In the case of  the simulated shock is directly compared to LASCO observations of the 28 October 2003 event , which shows agreement with both the general morphology and the quantitative brightness of the shock front as shown in Fig. 12. Here, the increase in brightness at the shock, in both the observed event and the model, is only at a level between 1% and 2%, even though in the model the density increase at the shock is very near the theoretical limit of 5 (for γ = 1.5). The small increase in brightness is due to the relatively short distance along the line-of-sight spanned by the compressed plasma behind the shock front. In this particular example, the shock morphology is largely determined by its interaction with the ambient solar wind and may not be sensitive to the initiation process.
Multidimensional simulations have shown how the geometry of the shock front is affected by the structure of the ambient solar wind. For example, Odstrčil et al. (1996), Riley et al. (1997), and  simulated coronal shocks driven by CMEs initiated at low latitude in solar minimum conditions where the shape of the shock front is strongly affected by latitudinal gradients in solar wind speed and density. Shock fronts propagate significantly faster at high latitude than at low, resulting in a shock surface with a saddle-like structure that is outward-concave in latitude while outward-convex in longitude, as shown in Fig. 11a. Examination of a shock in this configuration by  revealed that in the concave region, the fast-mode shock deflection away from the shock normal necessarily drives converging flows that lead to a significant density enhancement behind the shock as seen in the line plots in Fig. 11c In complement, line plots in Fig. 11d show the shock Mach number, compression ratio and shock angle (between B and the shock normal, θ BN ) as functions of distance from Sun center.

Flux Rope Over-Expansion and Forward-Reverse Shock Pairs
Flux ropes erupting from the Sun have high internal magnetic pressures that drive CME expansion, a process that may persist well beyond 1 AU. Departures from force-free equilibrium are common in ICMEs, providing an overpressure that can drive expansion (Liu et al. 2006a). In fact, ICME expansion is often observed in situ as declining speed profiles (e.g., Liu et al. 2006a), which are consistent with over-expansion or may be a residual effect of initiation. At high helio-latitudes, measurements by the Ulysses spacecraft led to the identification of a new class of disturbances produced by CMEs coasting through the high-speed wind (Gosling et al. 1994). Bounded by a pair of forward and reverse shocks, these events were also defined by low internal densities, temperatures, and hence, pressures. It was inferred, and confirmed by numerical calculations (e.g., Riley et al. 1997Gosling et al. 1998;Reisenfeld et al. 2003) that the initial high internal pressure drove magnetosonic waves into the surrounding plasma in all directions, which then steepened into shock waves. The shocks are divided into two distinct modes: forward shocks that travel ahead of the ejecta, and in contrast, reverse shocks that propagate back toward the Sun in the rest frame of the ejecta, while actually being carried away from the Sun by the solar wind. In the case of a reverse shock, fast low-density plasma passes through the shock to produce slow, dense, hot shocked plasma on the far side (away from the Sun) of shock. Similar processes occur at low latitudes, but because the ejecta is usually immersed within slow and denser material, the interaction is more asymmetric (Gosling 2000) and occurs without the companion reverse shock.
CMEs associated with forward-reverse shock pairs have many distinguishing properties. First, all of these CMEs are found at high heliographic latitude embedded in high-speed streams. Furthermore, these CMEs have an expanding velocity profile with a central speed approximately the same as the fast ambient solar wind. These CMEs commonly exhibit roughly symmetric increases in pressure, density, temperature, and magnetic field strength down stream of the shock pairs, while the minima of these plasma quantities occur near the centers of the disturbances. The magnetic structures associated with forward-reverse shock pairs have the features of magnetic clouds with shock-enhanced field strengths that are typically a factor of 2 or 3 greater than the ambient solar wind field. Counter streaming electrons suggest closed field lines in the CME, while the ion charge state of CME plasma is more typical of the ambient solar wind rather than a more highly ionized ejecta. While this is the common and ideal structure of CMEs with forward-reverse shocks, some examples are much more irregular and asymmetric (Gosling et al. 1998). Manchester and Zurbuchen (2006) suggested an alternative explanation for forwardreverse shock pairs in which the reverse shock forms as a result of deflections of the solar wind caused by the passage of the CME. In this model, fast solar wind overtakes slower plasma forming a reverse shock at high-latitude poleward of the ejected flux rope. Apparent signatures of ICMEs, such as enhanced magnetic field strength, low plasma beta, and field direction rotations, are produced by the plasma flows outside of the flux rope. The salient feature of this mechanism is that the reverse shock can only form laterally beyond the ejecta, as shown in Fig. 11a. Further analysis by Manchester and Zurbuchen (2007) shows that the stream-deflection mechanism reproduces many observed characteristics of high-latitude ICMEs bounded by forward-reverse shock pairs, including: magnetic fields with increased strength, rope like rotations, entrained with low pressure, and low density plasma, all contained between nearly symmetric pressure/density enhancements created by the forward and reverse shocks. In this case, stream deflections naturally produce an expanding velocity profile imbedded in the fast wind, and also explain the fast-wind charge state composition found between the shocks, all as a direct consequence of the interaction of solar wind streams deflected by the passage of an ICME.
Apart from shocks, large disturbances in the solar wind can be driven by the passage of the CME ejecta. One example is the draping of the interplanetary magnetic field (IMF) around the ejecta, which forms strong out-of-the-ecliptic magnetic fields (e.g., Gosling and McComas 1987;Odstrčil and Pizzo 1999a). In the case of field line draping, open magnetic field lines are bent around the ICME ejecta and compressed to have enhanced field strengths that may be two or three times as great as found in the ambient solar wind, as was found in Pioneer observations examined by McComas et al. (1988). A CME simulation showing such draping is shown in Fig. 11. Here, open field lines are colored to show field strength, which is particularly enhanced on the Sun-ward side of the flux rope. Another important disturbance is found trailing fast ICMEs, the large rarefactions that form in the trailing wake of ICMEs (e.g., Gosling and McComas 1987;Riley and Gosling 2007;). Here, the ICME outruns the material following it, forming an expansion wave (or rarefaction region) of low density . Sunward of fast CMEs, the magnetic field is also distorted away from the nominal Parker Spiral configuration into a near radial configuration (e.g., Riley and Gosling 2007;.

MC/ICME Distortion
Although suspected, the substantial distortion that ICMEs undergo was not fully realized until the Ulysses mission, when combined with Geotail, provided the first multi-spacecraft in situ measurements of the same CME at significantly different helio-latitudes and distances from the Sun (Hammond et al. 1995). Liu et al. (2006b) similarly analyzed Ulysses in situ observations to find that MCs are highly flattened. At solar minimum, MCs have curved cross sections that are outward concave, while at solar maximum the MCs are outward convex. The results are indicative that the MCs are deformed by their interaction with the ambient solar wind. White-light observations of ICMEs at large elongation made with SMEI found a relatively small number of ICMEs to be in the shape of outward-concave arcs, consistent with deformation solar wind velocity shear (Kahler and Webb 2007). Liu et al. (2008c) found that at solar minimum eight out of eleven MCs were outward concave as determined by the in situ measurement of the elevation angle of the MC normal with respect to the solar equatorial plane. Finally, as reported by Savani et al. (2010), coronagraphs onboard the STEREO-B observed the distortion of the leading edge of the CME/ICME event of 14 November 2007, which began close to the Sun and progressively increased as the CME moved to large elongation.
Subsequent numerical simulations provided ever-increasingly more accurate explanations for how these deformations evolve as the flux rope transits the space from the Sun to the Earth and beyond (e.g., Riley et al. 1997Odstrčil and Pizzo 1999b;Manchester et al. 2004a). These simulations show the deformation occurring from the exchange of momentum between the ejecta and the surrounding medium. Early on, radial expansion of the ejecta naturally provides a flow that evolves the erupting flux rope in a nearly selfsimilar fashion (e.g., Gibson and Low 1998). This simple evolution breaks downs and the ICME flattens as the solar wind inhibits the expansion of the MC whose radial extent becomes more nearly constant, while simultaneously, the spherically diverging flow causes the transverse extent to increase linearly with distance from the Sun. The aspect ratio of the ejecta, (transverse versus radial size) begins to increases with time, which simulations show to be particularly noticeable beyond 40 R .
Simulations also illustrate the deformation of CMEs by bimodal ambient solar wind typical of solar minimum (e.g., Gosling 1999;Odstrčil and Pizzo 1999b;Riley and Crooker 2004;Manchester et al. 2004a). At the largest scales, ICMEs are deformed into a highly flattened, concave-outward structure by the latitudinal velocity gradient (speed increasing with latitude) that serves to carry the high-latitude portion of the flux rope to greater distances. An example of a 3-D numerical simulation of a CME in a bimodal wind is shown in Fig. 13. Here, the ejected flux rope is shown from two perspectives at 4.5 and 60 hours after initiation in the top and bottom rows, respectively. Close to the Sun, the flux rope reflects its original shape greatly expanded. As the flux rope approaches 1 AU, it appears flattened and outward-convex with field lines connected to the IMF. At solar maximum, the solar wind velocity is relatively slow and unstructured so MC are deformed by aerodynamic drag to a convex shape, which is particularly true for fast CMEs (e.g., .

Flux Rope Erosion
Magnetic reconnection is ubiquitous in the solar wind, from large-scale structures such as the heliospheric current sheet (e.g., Gosling et al. 2005Gosling et al. , 2006 to the smallest scale associated with magnetic turbulence. So it should be no surprise that there are clear signatures that reconnection occurs with magnetic ejecta in interplanetary space (Farrugia et al. 2001). Reconnection between the IMF and the ICME field is most prone to occur where morenearly-opposed fields are brought into contact, which for a highly twisted horizontal flux rope would tend to be either at the front or rear of the magnetic flux rope. Dasso et al. (2006) point out that such reconnection will reduce the magnetic flux of a rope, resulting in Fig. 13 Flux rope structure and evolution in a simulated CME. The magnetic field of a CME simulated with a spheromak flux rope is shown at 4.5 hours (top row) and 60 (bottom row) hours after initiation. (a) and (c): simulation results seen viewed from above (+z axis). (b) and (d): side and oblique views of the simulated CME. Close to the Sun, the flux rope evolves in a nearly self-similar fashion. Far from the Sun, the flux rope is deformed to an outward-convex crescent, while the outer field lines have reconnected with the IMF an imbalance of poloidal flux, a process called flux rope erosion. A schematic of this process is shown in Fig. 14 adapted from Ruffenach et al. (2012). Here, Fig. 14a shows the flux distribution of a cylindrical flux rope prior to reconnection with the IMF, where as Fig. 14b shows the outer layers of the flux rope peeled away by reconnection, leaving an imbalance in poloidal flux.
The analyses of MCs in many studies (e.g., Chané et al. 2005;Dasso et al. 2006Dasso et al. , 2007Ruffenach et al. 2012Ruffenach et al. , 2015Lavraud et al. 2013Lavraud et al. , 2014 have revealed the predicted features of flux rope erosion. Dasso et al. (2006) examined the 21 January 2005 ICME and found a large imbalance in B z that they attribute to magnetic reconnection. Examination of the 9 November 2004 9 ICME by Dasso et al. (2007) revealed a long lasting tail (more than a day) of flux that is interpreted to be the result of magnetic reconnection at the front (and not at the rear) of the MC. A similar tail was found with the 18-20 October 1995 ICME, where the remaining magnetic field is interpreted as a partially reconnected flux rope by Dasso et al. (2006). Ruffenach et al. (2012) and Lavraud et al. (2014) closely examined several ICME events and found evidence (including reconnection jets) that magnetic reconnection occurred on the front side of the MC flux ropes, leading to their erosion. Ruffenach et al. (2012) examined the 19-21 November 2007 event with multi-spacecraft including STEREO A, B, ACE, Wind and THEMIS to establish that the axis of the rope was near the  (Lepping and Wu 2007). Ruffenach et al. (2015) analyzed the magnetic flux content of 50 MCs and found erosion at the front or rear in similar proportions with an average of 40% of the total azimuthal magnetic flux removed by reconnection. 30% of the magnetic clouds showed signs of reconnection.
Numerical simulations have also been used to study magnetic erosion in CMEs (e.g., Schmidt and Cargill 2003;Manchester et al. 2014a,b). The earliest examples are two-dimensional numerical experiments conducted by Schmidt and Cargill (2003), where the cylindrical flux rope was made to move out through open flux and through the current sheet. It was found that reconnection occurs on the lead/trailing edges of the flux rope in the case of fast/slow CMEs, and that the reconnection rate increases with the relative speed of the rope. An example from this simulation is shown in Fig. 15a, where the flux rope is initially traveling at 1.5 times the solar wind speed through the current sheet. Reconnection is manifest by the formation of small islands on opposite sides of the flux rope.
The magnetic reconnection rate is greatest close to the Sun where field strengths and Alfvén speeds are highest, with the rate decreasing with distance from the Sun. This evolutionary trend was shown in Manchester et al. (2014b), where the poloidal flux of the ejected rope is integrated and plotted as a function of time. Here, 14% of the total flux is lost during the first hour after initiation followed by additional 6% loss as the ICME propagates to 1 AU. A more extreme example of CME-IMF reconnection is found in the simulation of the 15 May 2005 ICME (Manchester et al. 2014a) where the poloidal field of the ejected rope is in the opposite direction of the surrounding streamer belt. Figure 15b provides a global view of the reconnection between the MC and the IMF occurring primarily at the north and south extremities of the flux rope. As the ICME approaches 1 AU, the MC is largely reconfigured by reconnection, yet the portion passing Earth retains the B z orientations of the original rope and reproduces in situ observations of total field strength and B z . As the result of reconnection, a long duration of B z -south trails behind the MC.  examines simulations of reconnection between multiple CMEs, which will be discussed in Sect. 9.  Schmidt and Cargill (2003). Here, magnetic reconnection with the IMF occurs near the top and bottom front of the flux rope resulting in the formation of a pair of magnetic islands. (b): the results of the 3-D simulation of Manchester et al. (2014a). In this case, reconnection occurs at both the front and back sides of the flux rope A magnetic flux rope can also be distorted and the flux redistributed to produce localized flux imbalances and give the appearance of flux rope erosion. Manchester et al. (2014b) provides a CME simulation exhibiting an extreme case of flux redistribution that mimics flux rope erosion. In this case, as a fast ICME propagates through the solar wind, it undergoes strong deceleration, as previously discussed in Sect. 5. As a result, there is a tendency for denser filament material to move forward through the decelerating flux rope. Manchester et al. (2014b) have found that such movement is accompanied by azimuthal flows that can transport poloidal flux away from the nose of the ICME producing a large imbalance in B z , giving the appearance of flux rope erosion, while, in fact, the majority of the flux remains intact. In this example, reconnection did occur with the IMF, but only removed less than twenty percent of the toroidal flux of the rope at 1 AU. Evidence of this behavior may be found in the 21 January 2005 ICME. In situ observations of this event show very dense cold plasma with a charge state composition consistent with filament material following immediately behind the ICME sheath, and most significantly the MC is almost entirely unbalanced in poloidal flux.

Plasma Evolution
Plasma that escapes the corona is heated to temperatures near 1.5 MK, which is sufficiently high to partially ionize all heavy elements. The freeze-in location of each ion charge state Fig. 16 Solar wind plasma and magnetic field parameters across an ICME (after Zurbuchen and Richardson 2006). The vertical green line marks the arrival time of the ICME-driven shock, and two vertical pink lines indicate the ICME interval. This event can be classified as a magnetic cloud, showing an enhanced, smoothly-rotating magnetic field. Other ICME signatures are also present: (bidirectional electrons) BDEs (with the interval indicated by the horizontal bar), depressed proton temperatures, declining solar wind speed profiles, He/proton abundance enhancements, enhanced oxygen and iron charge states and Mg/O ratio, cosmic ray depressions ( 60 MeV particle count rate), and geomagnetic storms (indicated by the D st index) (typically within 4 R ) depends both on ionization and recombination rates, which are functions of plasma density and temperature. The in situ charge state composition, therefore, complements remote observations and provides a unique and very important measure of the heating and expansion of the solar wind near the Sun (e.g., Fenimore 1980;Galvin 1997;Henke et al. 1998;Lepri and Zurbuchen 2004;Reinard 2005;Zurbuchen and Richardson 2006). Used in combination, ion freeze-in states can be used to map out the heating and expansion rates of CMEs and as well as the ambient solar wind (e.g., Buergi and Geiss 1986;Esser and Edgar 2001;Lynch et al. 2011;Jin et al. 2012) Numerous studies have shown the value of ion charge states as identifiers of ICMEs, finding that MCs are more likely to exhibit a higher composition of elevated charge states than non-MC ICMEs (e.g., Henke et al. 1998). Figure 16 (adapted from Zurbuchen and Richardson 2006) provides an example of ICME charge state signatures shown in conjunction with the MHD plasma quantities, supra-thermal electron populations and magnetospheric DST. Enhanced helium abundance (He 2+ ) has been observed in ICMEs since the early space era (Hirshberg et al. 1972;Borrini et al. 1982), and roughly 50%-70% ICMEs have enhanced iron charge states up to Fe 16+ and average (Q)Fe > 12 (e.g., Bame et al. 1979;Fenimore 1980;Lepri et al. 2001;. Enhanced charge states of oxygen are also found in ICMEs (e.g., Galvin 1997;Henke et al. 2001) with O 7+ /O 6+ > 1 ( Richardson and Cane 2004a). The elevated charge states of oxygen and iron (Lepri and Zurbuchen 2004) signify an intense heating of plasma close to the Sun, while the cause of the helium enhancement is less clear. Neugebauer and Goldstein (1997) ascribe it to a "sludge removal phenomenon" whereby the helium settling at the base of the solar wind flux tube is cleared out by the CME eruption. Bimodal charge state distributions are also found in ICMEs where both elevated charge states associated with flare heated plasma are found along with ambient solar wind charge states . This condition suggest that either the material expands in such a way that two temperature peaks form (e.g., Gruesbeck et al. 2011) or two plasmas mix (e.g., . A small number of ICMEs have also been identified with intervals of highly elevated levels of unusually low ion charge states, such as He + (e.g., Gosling et al. 1980;Illing and Hundhausen 1985;Burlaga et al. 1998). In some cases, the He + /He ++ ratio may reach 0.1, which is nearly ten thousand times greater than ambient levels. The singly charged He particles imply that these cold events originate from a low temperature region, possibly a cool, dense filament in the solar atmosphere. This assertion has found further support from Lepri and Zurbuchen (2010) who found unusually cold material that they attribute to prominence material that escaped the corona without further heating. Such events remain rare compared to the observed fraction of CMEs associated with prominence eruptions. A possible example is the 21 January 2005 ICME event, where both charge state analysis and simulation (Kozyra et al. 2013;Manchester et al. 2014b) strongly suggest that cold dense prominence material associated with a fast CME impacted the Earth.
Electrons, by contrast, often show strange thermodynamic behavior across the ICME plasma. The electron temperature, T e , is often larger than the proton temperature within ICMEs, with a typical ratio T e /T p > 3 (Liu et al. 2005), so the thermal pressure is dominated by the electron component. However, T e tends to be lower than the proton temperature in the sheath regions of ICMEs with preceding shocks (Skoug et al. 2000), analogous to planetary magnetosheaths. A possible explanation is that protons resonate at their gyro-frequency with magnetic fluctuations induced at the preceding shock and thus gain energy from the fluctuations, while electrons cannot do so. Another explanation is that the shocks are supersonic relative to the protons, but subsonic relative to the electrons, and as a consequence, the protons are preferentially shock heated (Kosovichev and Stepanova 1991;Jin et al. 2013). Another feature of electron thermodynamics is that the polytropic index γ e < 1 inside ICMEs (Gosling 1999), the physical cause of which remains a mystery.
The radial evolution of T α /T p and v αp from 0.3 to 5.4 AU is shown in Fig. 17. The differential streaming in the solar wind (lower panel) shows a fairly steady decrease with distance, probably due to the accumulative effect of Coulomb collisions. In contrast to the solar wind, the ICME v αp decreases to less than 10 km s −1 by 0.4 AU, which indicates the importance of Coulomb collisions inside ICME plasma. The temperature ratio T α /T p of ICMEs is larger than that of the solar wind inside 1 AU but roughly the same magnitude beyond 1 AU. The little variation of the ICME temperature ratio with distance implies the existence of a local mechanism responsible for heating the alphas. One may consider that the free energy associated with the differential streaming is converted into the thermal energy of the plasma. It is not likely that the free energy preferentially goes to alphas through Fig. 17 Radial variations of the temperature ratio and differential speed between alphas and protons for ICMEs and solar wind (after Liu et al. 2006a). For ICMEs, the horizontal bars indicate the bounds of the bins while the error bars show the standard deviation of the parameters within the bins. The solar wind levels are represented by the average values over 0.1 AU bins within 1 AU and 0.5 AU bins beyond. All the data are within ±20 • in latitude Coulomb collisions. In addition, the differential speed within ICMEs is very small outside about 0.4 AU, but the temperature ratio is conserved at much larger distances. These results, in particular the collisional nature of the ICME plasma which has been ignored in previous studies, have significant consequences for modeling the realistic thermodynamics of ICMEs. The readers are directed to Liu et al. (2006a) for more detailed discussions.
Recent advancements in techniques for analyzing the solar wind composition measured have provided the most comprehensive ionic composition data to date, which set solid constraints on the heating and acceleration of the solar wind and ICMEs. These observations are even more powerful by combining them with remote observations of coronal plasma, which reveals details about CME expansion and solar wind source regions (e.g., Gruesbeck et al. 2011;Jin et al. 2012). These data are further complemented by recent simulations (Lynch et al. 2011) of ICME and solar wind (Jin et al. 2012) charge state composition.
Filaments (or more specially the filament channel) are at the heart of the CME energy source, where highly sheared magnetic fields run nearly parallel to the PIL in an energized state (Manchester 2001Fang et al. 2010Fang et al. , 2012. Here, filaments form as cold dense concentrations of plasma suspended in the corona at chromospheric temperatures with densities 100 times greater than the surrounding corona. Filament material is typically taken to be the source of the dense core of the canonical three-part structure that is often wholly expelled in CMEs (Hundhausen 1993;Gibson and Low 1998). Gilbert et al. (2000) found 94% of erupting filaments were in association with CMEs and of those, 84% were found to have a portion of the filament escape the corona and contribute to the ejected CME plasma. However, connecting the components of the CME three-part density structure to in situ counterparts is often not possible. While a pileup of plasma is commonly found ahead of the flux rope, signatures of a cavity are not always clear and filament material is most often not detectable. ICME expansion tends to produce densities (and kinetic temperatures) lower than that of the surrounding solar wind (Liu et al. 2006a). However, there are often density fluctuations within the cloud that attain or exceed the ambient value (Burlaga 1981;Burlaga et al. 1982), which may obscure the signatures of the cavity.
A key to distinguishing different CME structures is the plasma charge state composition, which is determined by the ionization environment. Rapid decreases in mass density and collision rates cause the charge state to become fixed, or frozen in, beyond a few solar radii (e.g., Hundhausen et al. 1968;Ko et al. 1999;Kilpua et al. 2013). Given its relatively cool temperatures, filament material in the corona is dominated by heavy ion charge states such as C 2−3+ , O 2−4+ and Fe 4−6+ , which have been directly observed with SOHO/UVCS (e.g., Lee and Raymond 2012). Filament material can be ejected into interplanetary space where it can be distinguished by these ionization states (Schwenn et al. 1980;Burlaga et al. 1998;Skoug et al. 1999;Yao et al. 2010). Gopalswamy et al. (1998) provides a clear example of the identification of low charge state filament material in the 07-11 February 1997 CME/ICME event. Here, filament material was observed in SOHO/LASCO C2 coronagraph images at the base of the CME, while corresponding dense filament material was found in situ in the Sun-ward side of a magnetic cloud containing relatively high abundances of Fe 5+ and Fe 6+ .
While filament eruptions are commonly associated with CMEs, it is surprising that convincing in situ charge state evidence of filament material is rarely found. Lepri and Zurbuchen (2010) reviewed 10 years of data from SWICS/ACE and found only 4% of ICMEs (11 events of the total) with charge state consistent with filament material. The relative paucity of in situ filament material observations may have been resolved by recent observations. Lee and Raymond (2012) found that low-ionization filament material (as observed by UVCS) covers only a few percent of the total area of CMEs, which is consistent with the low rate of in situ detection. This explanation is further corroborated by remarkable observations by Wood et al. (2016), which found that unlike most plasma in the CME, the filament material does not expand to fill space but fragments into high-density parcels of plasma that ballistically propagate through the heliosphere. Furthermore, in the case of fast ICMEs that strongly decelerate, the dense filament material maintains a high velocity and moves radially outward through the body of the CME, just as first predicted by Manchester et al. (2014b). 9 CME-CME Interaction ICME properties may be strongly influenced as they propagate between the Sun and the Earth, in particular when they interact with other ICMEs or with corotating interaction regions (CIRs) and fast wind streams (Prise et al. 2015;. Two of the main sources of successive CMEs are homologous and sympathetic eruptions. Solar observations have revealed that recurrent CMEs occur from the same active region, often associated with homologous flares (Schmieder et al. 1984;Moon et al. 2003;Liu et al. 2014a,b;Wang et al. 2014a). Sympathetic CMEs, when the eruption of one CME triggers another, is another source of successive CMEs in relatively close angular and temporal separation (see Fig. 18a), revealed by means of observations and simulations to occur regularly (Schrijver and Title 2011;Wang et al. 2016a).
During their Sun-to-Earth propagation, the interaction of successive ICMEs may take a variety of forms: (i) multiple CME-driven shock waves may interact without the associated ejecta interacting, (ii) one shock wave may interact with a preceding magnetic ejecta, or (iii) the successive magnetic ejecta may interact and/or reconnect. Early detections of multiple flares and complex ejecta in the interplanetary space revealed in the 1970s and 1980s that CME-CME interaction is a frequent occurrence. For example, the series of seven flares in 72 hours in early August 1972 have been extensively analyzed (Dryer et al. 1976;Intriligator   . (c): Radial size of simulated CMEs before and during their interaction depending on their initial relative orientation, as adapted from . (d): 3-D view of a MHD simulation of two interacting CMEs from Shen et al. (2013) 1976; Ivanov 1982). Another example is the study by Burlaga et al. (1987), in which the authors describe a variety of compound streams resulting from the interaction of a transient with another transient or with a solar wind stream. The authors concluded that "the interaction between two fast flows, is in general a nonlinear process, and hence, a compound stream is more than a linear superposition of its parts". An early numerical investigation is that of Vandas et al. (1997), who simulated the interaction of a shock wave with a magnetic cloud in 2.5-D MHD. The authors noted that the shock propagation results in a radial compression of the magnetic cloud, a change of its aspect ratio, and acceleration as well as heating of the cloud.
Direct observations of CME-CME interaction became possible in the mid-1990s with the larger field-of-view of LASCO/C3 (up to 32 R ∼ 0.15 AU), which yielded the first reported white-light observations of CME-CME interaction (Gopalswamy et al. 2001c). Interest in this topic grew following this initial report (e.g., see Reiner et al. 2003), and from debates about the association of interacting CMEs with large solar energetic particle (SEP) events . In solar cycle 24 (2008-present), observations with high spatial and temporal resolution by SDO have returned the study of sympathetic eruptions to central stage. The development of heliospheric imaging with the Solar Mass Ejection Imager (SMEI, Eyles et al. 2003;Jackson et al. 2004) and the Heliospheric Imagers (HIs, ) onboard STEREO have led to a large increase in the number of published cases of ICME-ICME interactions being remotely observed in the past 10 years (see Fig. 18b). From the first remote detection in 2001 to routine remote observations in the late 2000s, numerical simulations have been used to fill the gap between the upper corona and the near-Earth space Schmidt and Cargill 2004;Lugaz et al. 2005b;Xiong et al. 2006). In the past decade, the combination of remote observations, in situ measurements and numerical simulations has resulted in a deeper understanding of the physical processes occurring during CME-CME interaction.

Compression and Reconnection Between CMEs
One of the essential aspects of CME-CME interaction is the change in CME properties, such as their speed, size, expansion rate, etc. CME-CME interaction directly affects space weather forecasting, by modifying the hit/miss probability and the expected arrival time of a CME (changes in the CME speed and direction) and by modifying the expected geomagnetic responses (changes in the CME internal magnetic field).
Here, we first discuss the changes in the CME internal properties, such as radial extent and magnetic field strength, associated with the compression and reconnection between CMEs. The radial size of a CME can be derived from remote observations (Savani et al. 2009;Nieves-Chinchilla et al. 2012;. These observations and numerical simulations confirmed that the radial extent of the leading CME plateaus during the main phase (i.e. when the speed of both CMEs changes significantly) of interaction (Schmidt and Cargill 2004;Lugaz et al. 2005bXiong et al. 2006) and this is typically associated with a "pancaking" of the leading CME (Vandas et al. 1997). In Lugaz et al. (2005b), the authors discussed how the shock propagation through the leading CME is the main way in which the expansion slows. In situ measurements at 1 AU of shocks propagating inside CMEs (Lugaz et al. 2015) revealed the compression of the leading CME by the overtaking shock and CME (Wang et al. 2003b;Farrugia and Berdichevsky 2004;Liu et al. , 2014c. It is, as yet, unclear whether this compression changes for cases when the overtaking CME drives a shock or not. What clearly changes is the resulting expansion of the leading CME after the end of the main interaction phase (i.e. after the shock exited the ejecta). In Xiong et al. (2006), the authors found that the leading CME over-expands to return to its expected size; as such the compression is only a temporary state. This was confirmed by a statistical study of the radial size of magnetic ejecta at different heliocentric distances (Gulisano et al. 2010). In one case study, similar results were found: remote observations indicated compression but, a day after the interaction ended, when the CME impacted Earth, the in situ measurements indicated a typical CME size . In numerical simulations with two magnetic ejecta, the rate of over-expansion is found to depend on the rate of reconnection between the two ejecta (see Fig. 18c); as such, it depends on the relative orientation of the two magnetic ejecta (Schmidt and Cargill 2004;), but also probably on their density. Associated with ICME compression is the enhanced magnetic field inside the ICMEs as revealed by in situ measurements.  find that an extremely enhanced field strength, larger than 100 nT, can be produced by CME-CME interactions. CME-CME "cannibalism", i.e. the total coalescence of two CMEs, was reported based on LASCO images by Gopalswamy et al. (2001c). However, because coronagraphs and white-light imagers only give direct information about a CME density structure, it is unclear whether this coalescence is associated with a total reconnection of one CME's magnetic flux or not. Known magnetic reconnection rates and CME magnetic flux makes it unlikely that total reconnection is a common occurrence until beyond Earth's orbit. The potential full coalescence of two ejecta into one was discussed in a few studies Schmidt and Cargill 2004;Chatterjee and Fan 2013) but has not been investigated in detail with realistic reconnection rates. In situ measurements of transients resulting from the interaction of two CMEs sometimes contain a period of weaker and more turbulent magnetic field, higher in density and temperature, resulting in a higher plasma β (Wang et al. 2003a). These periods have been explained by ongoing or past reconnection between the two magnetic ejecta; however multi-spacecraft measurements or a successful Walén test (Paschmann and Sonnerup 2008) would be required to assert that reconnection is occurring. No such case has been reported yet.
In , the authors considered the interaction of two CMEs with axial fields perpendicular to one another by means of MHD simulations, and found that the resulting ejecta appeared as one of a single, extended but seemingly isolated, CME with a complex magnetic field rotation. Such events can sometimes be measured at 1 AU (Dasso et al. 2009;Lugaz and Farrugia 2014) and may be related to these instances of CME-CME cannibalism.

Momentum Transfer Between CMEs
Many studies have investigated the momentum transfer between CMEs during their interaction, both through remote observations and numerical simulations (Schmidt and Cargill 2004;Lugaz et al. 2005b;Shen et al. 2012;Mishra et al. 2015a,b;Colaninno and Vourlidas 2015). CMEs are large-scale magnetized plasma structures propagating in solar wind, and therefore, the momentum transfer between CMEs is complex for a variety of factors: (i) depending on the speed of the CMEs, the interaction between two CMEs may involve zero, one or two CME-driven shocks, some of which may dissipate during the interaction, (ii) the interaction takes at least one Alfvén crossing time of a CME, i.e. typically 8-24 hours, (iii) the CME speeds can change significantly, even at large distances from the Sun, due to their interaction with the solar wind (Manchester et al. 2004a;, (iv) CME-CME interaction is inherently a three-dimensional process, and the changes in kinematics may differ greatly depending on the CME part that is considered (Shen et al. 2012;). Using numerical simulations (see Fig. 18d), it is somewhat possible to control for some of these effects, for example, by performing simulations with or without interaction but with identical CME properties, and by knowing the velocity field in the entire 3-D domain (Shen et al. 2013(Shen et al. , 2016. As CME-CME interaction involves a faster CME overtaking a slower CME, the end result is to homogenize the speed, as was noted from in situ measurements in Burlaga et al. (2002), Farrugia and Berdichevsky (2004) and through simulations by Schmidt and Cargill (2004) and Lugaz et al. (2005b), among others. It has been found to occur independently of the relative speeds of the two CMEs. Based on this finding, one of the main issues related to momentum transfer is to understand what determines the final speed of the transient that was formed by the CME-CME interaction. In an early work, Wang et al. (2005) found that, in the absence of CME-driven shock waves, the final speed is determined by that of the slower ejecta, whereas Schmidt and Cargill (2004) and Lugaz et al. (2005b) found that when the CMEs drive shocks, the final speed is primarily determined by that of the faster ejecta, as the shock's propagation through the first magnetic ejecta accelerates it to a speed similar to that of the second ejecta. Work combining remote observations and numerical simulations lead us to conclude that the final speed of the transient depends on the relative masses of the CMEs, as well as their approach speed (Poedts et al. 2003;Shen et al. 2016).
In addition to changes in velocity, CME-CME interaction may result in the deflection of one CME by another (Xiong et al. 2009;Shen et al. 2012). Combining these works, it appears that the deflection can reach up to 15 • when the two CMEs are initially about 15-20 • apart. Such angular separations are quite frequent between successive CMEs, as it corresponds to a delay of about one day for two CME originating from the same active region (due to solar rotation). This change in direction must be taken into account when deriving the changes in velocity, as done in Shen et al. (2012) and .

Changes in the Shock Properties
In addition to changes in the CME properties, any fast forward shock propagating inside the overtaken magnetic ejecta encounters highly-varying and unusual upstream conditions, affecting the shock properties. Most of what is known about the changes in shock properties was learned from numerical simulations; however, there have been many reported detections of shocks propagating inside a magnetic cloud or magnetic ejecta at 1 AU (Wang et al. 2003c;Collier et al. 2007;Richardson and Cane 2010a;Liu et al. , 2014cLugaz et al. 2015. Vandas et al. (1997) noted that a shock propagates faster inside a magnetic cloud due to the enhanced fast magnetosonic speed inside, which may result in shock-shock merging close to the nose of the magnetic cloud but two distinct shocks in the flanks.  noted that associated with this acceleration, the density jump becomes smaller. Lugaz et al. (2005b) performed an in-depth analysis of the changes in the shock properties, dividing the interaction into four main phases: (i) before any physical interaction, when the shock propagates faster than an identical isolated shock due to the smaller density in the solar wind, (ii) during the shock propagation inside the magnetic cloud, when the shock speed in a rest frame increases and its compression ratio decreases, confirming the findings of , (iii) during the shock propagation inside the dense sheath when the shock decelerates, as pointed out by Vandas et al. (1997), and (iv) the shock-shock merging when, as predicted by MHD theory, a stronger shock forms followed by a contact discontinuity. If the shock is weak or slow enough, it may dissipate as it propagates into the region of higher magnetosonic speed inside the magnetic cloud (Xiong et al. 2006;. High spatial resolution is necessary to resolve weak shocks in MHD simulations, and low resolution may affect the prediction of shock dissipation. The merging or dissipation of shocks was noted by Farrugia and Berdichevsky (2004), when Helios measured four shocks at 0.67 AU and ISEE-3 measured only two shocks later on at 1 AU.
For Solar Cycle 23, Lugaz et al. (2015) identified in total 49 shocks that were propagating within magnetic ejecta. Most such shocks occur toward the back of the ejecta, and shocks tend to be slower as they get closer to the CME front. This can be interpreted as an indication that a number of shocks dissipate inside a CME before exiting it. The two main reasons are that CMEs tend to be expanding and have a decreasing speed profile and that the peak Alfvén speed typically occurs close to the center of the magnetic ejecta. The latter reason means that shocks become weaker as they approach the center of the ejecta. The former reason implies that shocks propagate into higher and higher upstream speeds as they move from the back to the front of the CME. Lugaz et al. (2015) reported cases when the speed at the front of the first CME exceeds the speed of the overtaking shock, i.e. because of the CME expansion, the shock cannot overtake the front of the CME.

Modeling CMEs from Sun to Earth: The Bastille Day Event
MHD simulations of CMEs can be roughly divided into two classes. On the one hand, there are models that consider idealized configurations and are mainly used to test suggested ideas and to investigate specific aspects of CMEs, such as their initiation mechanisms. This kind of simulation uses a simplified set of the MHD equations and numerical domains that typically cover only the lower corona (see Green et al. in this Volume). On the other hand, there are simulations that are designed predominantly to reproduce observed eruptions. Two examples are shown in Figs. 12 and 15, respectively. These semi-realistic simulations vary widely in complexity, depending on the physics included, the phases of the CME/ICME evolution modeled, and the extent of observational constraints and real data included. Many models use observed magnetograms as boundary condition for the magnetic field at the lower boundary, and recently observed surface flows have been included as well (Jiang et al. 2016). Some simulations restrict the calculation to the evolution in the corona (e.g., Zuccarello et al. 2012b;Amari et al. 2014;Fan 2016), while others model the propagation of the ejecta to 1 AU but do not include the coronal evolution. Instead, the latter models start the simulation in the inner heliosphere (typically at around 20-30 R ) and use for the initial ICME some idealized model whose parameters are chosen guided by coronagraph observations (e.g., Shiota and Kataoka 2016). Some of these models, such as ENLIL , which is the only MHD model that is currently used for operational space weather predictions (see Green et al. in this Volume), simplify even further by neglecting the magnetic field of the ICME and initiating its propagation by specifying a cone of constant velocity at the inner boundary of the model.
The currently most advanced simulations model CMEs all the way from the Sun to Earth (e.g., . The most recent ones employ a sophisticated treatment of the energy transfer in the corona, often referred to as "thermodynamic MHD" (Lionello et al. 2009;. Such a treatment is required for a realistic description of the plasma properties in the corona, in particular if one wants to reproduce simultaneously the emission in active regions, the quiet sun, and coronal holes (Lionello et al. 2009). Some models also allow for a separate treatment of the electron and ion temperatures (van der Holst et al. 2014). CMEs are then launched in such a thermodynamic MHD environment, often using an analytical magnetic flux-rope model such as the one by Titov and Démoulin (1999), and propagated into the corona (Lugaz et al. 2011;Downs et al. 2012;Jin et al. 2013) or, by coupling different MHD codes , further out to 1 AU and beyond (Lionello et al. 2013;Manchester et al. 2014a;Jin et al. 2017b). A long-term aim is to use these sophisticated simulations for space-weather predictions by modeling observed CMEs, especially the strongest, most geo-effective events, in real-time, i.e., while the eruption is still on its way to Earth (see, e.g., Jin et al. 2017a). In what follows, we briefly describe the properties and limitations of these modeling efforts, using as an example a recent simulation of the well-known "Bastille Day" eruption. Other examples can be found in the articles mentioned above.
The 14 July 2000 ("Bastille Day") event, one of the largest eruptions of solar cycle 23, occurred in NOAA active region (AR) 9077. It was associated with a fast halo CME, an X5.7 flare, and an intense radiation storm that resulted in one of the 16 GLE events of cycle 23. The CME traveled in the corona at a speed of ≈ 1700 km s −1 (Andrews 2001). The  Linker et al. (2016) shock wave driven by the interplanetary CME/ICME and a large magnetic cloud reached Earth in about 28 and 33 hours, respectively, triggering a large geomagnetic storm (Lepping et al. 2001). Strong eruptions such as the Bastille Day event are perfectly suited to test the capabilities of MHD simulations in reproducing the properties of "extreme" events.
A simulation of this event was performed using the 3D-MHD code MAS developed at Predictive Science Inc. (e.g., Lionello et al. , 2013. Here, we focus on the main properties of the simulation; a more detailed description can be found in Linker et al. (2016) and in a forthcoming publication (T. Török et al., in preparation).
To provide the background environment for the CME, first an MHD solution of the global corona was calculated, using a spherical grid that covers 1 to 20 R , and includes the sonic and Alfvénic critical points. The surface radial magnetic field, B r , was specified by combining a line-of-sight (LOS) SOHO/MDI synoptic map (10/7/2000-6/8/2000) with a LOS MDI magnetogram measured about 30 minutes prior to the flare onset. After calculating a global potential magnetic field based on this boundary condition, a thermodynamic MHD model of the corona is produced, using the equations and coronal-heating specification described by Lionello et al. (2009). The system is relaxed for about 64 hours, until a steady state, including the solar wind, is obtained. The inclusion of thermal conduction, radiative cooling, and empirical coronal heating provides realistic density and temperature distributions, which allows one to produce synthetic EUV or soft X-rays satellite images that can be directly compared to the observations. Figure 19a shows, as an example, a synthetic SOHO/EIT image of the corona obtained from the simulation.
The AR was then energized by inserting several instances of the modified Titov-Démoulin (TDm) flux-rope model (Titov et al. 2014) along the PIL. The original TD model of a force-free coronal flux rope (Titov and Démoulin 1999) has been widely used as an initial condition for CME simulations. Typically, a TD rope out of magnetic equilibrium is used to initiate a CME (e.g., Lugaz et al. 2011). This approach is convenient as it reduces the complexity of the computation and reproduces many observed CME properties sufficiently well, but neglects that in reality CMEs always arise from magnetic fields in an initial state of force balance. Thus, starting simulations from a stable magnetic equilibrium is important for understanding the energy storage and release of CMEs. The TDm model facilitates the construction of such configurations by taking into account the background coronal magnetic field, and was therefore, used in the simulation described here. To account for the highly curved and elongated PIL, along which the strength of the background field strongly varies, a chain of seven overlapping TDm ropes was inserted (Fig. 19b). The axial field orientation was chosen to produce a left-handed flux rope, as suggested by the observations (e.g., Yurchyshyn et al. 2001). The model was further modified to preserve the observed B r distribution when inserting the ropes, and the parameters of each TDm rope were chosen such that the merged flux rope is stable with respect to the kink and torus instabilities (e.g. Kliem and Török 2006). The resulting configuration is close to being force-balanced but is not in thermal equilibrium. Flows and plasma condensations appeared along the flux-rope field lines that possibly arose from thermal non-equilibrium (e.g., Mikić et al. 2013a). The system was relaxed for about one hour, during which time the numerical dissipation (enhanced by the thermal flows) leads to some loss of magnetic energy.
After this additional relaxation, the eruption of the flux rope was triggered by imposing sub-Alfvénic flows at the inner boundary that converge toward the PIL (e.g., Mikić et al. 2013b). These flows were imposed for about 50 minutes, at which point the configuration destabilized. The flows were not directly derived from the observations, but adjusted such that the eruption starts at the western end of the AR and proceeds toward the East, in qualitative agreement with the observations. The full eruption released 1.3 × 10 33 ergs in about 4 minutes; about 31% (4 × 10 32 ergs) was converted into kinetic energy. The simulated CME reached a speed of about 1500 km s −1 , close to the observed speed. The striking morphology of the flare arcade and the halo CME are qualitatively well reproduced in the simulation (see Linker et al. 2016). The simulation also self-consistently (i.e., without further boundary-driving) reproduces a second, "sympathetic" eruption, that was observed four hours after the main eruption in the western part of the AR (T. Török et al., in preparation).
Finally, in order to model the propagation of the ICME to Earth, the coronal simulation was coupled to the heliospheric version of MAS (Lionello et al. 2013), which uses a simpler set of equations and, is therefore, computationally more efficient. The coupling is done by setting the outer boundary of the coronal simulation as the inner boundary of the heliospheric simulation, using time-dependent boundary conditions as described in Lionello et al. (2013). The heliospheric calculation uses a 486 × 264 × 583 nonuniform spherical mesh extending from 19 to 230 R .
The observed ICME decelerated during its propagation, and arrived at Earth at about 19 UT on 15/7/2000 with a speed of ≈ 1100 km s −1 and a maximum field strength of ≈ 50 nT (e.g., . The magnetic-field measurements at 1 AU suggest a left-handed flux rope passed the spacecraft below its axis (Yurchyshyn et al. 2001). Figures 19c, d show the simulated ICME, with field lines colored by the magnitude of B θ ; positive B θ (red) is in the direction of negative B z at Earth. It can be seen that a negative B z arises from the azimuthal fields wrapping the core of the flux rope, and that the flux rope connects back to the inner boundary of the simulation domain. Figure 19e shows that the simulated ICME arrives at 1 AU seven hours later than the real event, with a speed of ∼ 900 km s −1 and a maximum field strength of ∼ 27 nT. The simulated ICME is qualitatively consistent with the flux rope inferred from the observations: a full rotation of B z from negative to positive, and B x , B y > 0. However, it passes Earth about 15-20 • northward of the real event. At Earth's location, the correct (negative) sign of B z is reproduced, but the rotation of B z from negative to positive is not seen, and the field strengths are weaker.
This example illustrates that state-of-the-art CME simulations are capable of qualitatively reproducing many important features of extreme events, such as the rapid, strong energy release, the observed emission and white-light morphology, and the magnetic structure of the ICME. The quantitative differences in the case shown here are relatively small, given that the simulation still uses various simplifications, such as a steady-state corona and heliosphere, and does not employ "trial-and-error" runs to match the observations (such as the CME speed) quantitatively. Yet, the discrepancies are still significant in light of the potential future application of such simulations for space-weather predictions. Further improvements, such as data-driven, time-dependent models of the background corona, and more case studies are required to assess to which accuracy state-of-the-art CME simulations can match the observations at 1 AU.

Conclusions
3-D numerical simulations and wide-angle heliospheric imaging of CMEs and ICMEs now routinely performed by STEREO, have resulted in a new era of study of CME/ICME evolution and their interaction in the heliosphere. Upon eruption, CMEs strongly interact with the surrounding environment from the active region magnetic field, the streamer belt, the solar wind and even other CMEs. A multitude of physical processes govern the transfer of energy, mass, momentum and magnetic flux to and from CMEs and ICMEs as they travel from low corona to interplanetary space, most of which can be well described in the context of extended MHD. CME deflection has an early impact on CMEs, where the trajectory departs from the expected radial direction and is most prevalent in the low corona. Here, magnetic forces can drive deflections in longitude and latitude that can reach 30-40 • for low-mass and slow-speed CMEs. Meanwhile, patterns of CME impact on Earth suggest a consistent eastward deflection related to spiral structure of the ambient solar wind. Apart from being deflected, CMEs can also rotate about their rise direction as they propagate in the corona and interplanetary space. Like deflection, rotations are largest in the low corona, but are less predictable given that a number of mechanisms, including reconnection of the ejecta with complex background magnetic fields, can contribute to CME rotation in a complicated manner. At present, it is unclear what determines the total amount of rotation in CMEs, and further efforts are required in order to gain a better understanding of this important problem.
Once free of the low corona where strong acceleration occurs, ICME evolution is foremost dominated by kinematic processes as ICMEs ballistically propagate through the solar wind. As a first approximation, this process is described by aerodynamic drag where slow CMEs are accelerated and fast CMEs decelerated by the interaction with the solar wind. These processes are sufficiently well represented by empirical and numerical models that they can typically predict CME arrival times at Earth to within eight hours or less. Consider next the reaction that ICME propagation has on the heliosphere. CMEs traveling faster than the solar wind drive forward propagating shocks that can extend far beyond the ejecta, and whose geometry can be largely modified by the structure of the ambient solar wind. Simulations show that radial expansion of the ejecta can drive high-latitude reverse shocks, which may also be caused by CME-deflected high-speed streams. The ejecta also expands laterally with the diverging solar wind at a rate that is faster than the flow deflected around the ejecta. As a result, plasma accumulates in a sheath surrounding the flux rope and adds substantially to the overall mass of the disturbance. The interplanetary magnetic field can similarly be compressed as it drapes around the ejecta, leading to significant enhancements in the overall flux and field strength. The IMF can also reconnect with the ejecta, eroding the magnetic rope by removing poloidal flux. Reconnection is common with 30% of flux ropes showing signs of erosion with an average of 40% of the flux removed at 1 AU. In addition to MHD quantities, charge state composition is measured both in situ and spectroscopically, which gives measures of the temperature and thermodynamic behavior of the low corona where the charge state freezes in. These data are particularly important in identifying cold filament material and flare heated plasma.
Many recent investigations have focused on the relation between CME-CME and ICME-ICME interaction and the resulting structures at 1 AU, as well as their exchange of momentum. CME-CME interaction is one of the ways in which large changes in the properties of CMEs may occur during their Sun-to-Earth propagation, in addition to interaction with the solar wind. It often results in a compression and acceleration of the leading CME, although this might only be a temporary effect; it can also result in a partial or total coalescence of the two interacting CMEs (CME-CME "cannibalism"). There remains significant debate about momentum exchange between CMEs and the resulting final speed. In many instances, after their interaction, the two CMEs propagate with a uniform speed profile, similar to that of an isolated CME, but how this speed relates to that of the two CMEs before the interaction is still an area of active research. Other potential effects are a deflection of the two CMEs, and the "disappearance" of shock waves, either through shock-shock merging or as the trailing shock dissipates as it propagates inside the leading magnetic ejecta. About 15% of shocks at 1 AU are found to propagate inside a preceding magnetic ejecta, a clear example of CME-CME interaction in the process of happening.
Numerical simulations have evolved to a state where a semi-realistic modeling of observed eruptions has become possible. State-of-the-art simulations employ measured magnetograms, a realistic treatment of the energy transfer in the corona, and propagate ICMEs to 1 AU or beyond. Using the example of the well-known Bastille Day event, we described the capability of such simulations to reproduce many important observational features of eruptions, which makes them extremely valuable for improving our physical understanding of these complex events. Further improvements in computational power and in the accuracy of the simulations are required, though, before they can be used for forecasting the impacts of CMEs on space weather close to Earth.
Together, these results indicate general pictures and universal features, that regulate the Sun-to-Earth propagation of different types of CMEs. Although exceptions from the norms can occur, depending on the coronal and interplanetary conditions where a CME/ICME propagates, the results are representative of many events. Any theory/model of CME/ICME Sun-to-Earth propagation and space weather forecasting should be guided and constrained by these results.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.