The Origin, Early Evolution and Predictability of Solar Eruptions

Coronal mass ejections (CMEs) were discovered in the early 1970s when space-borne coronagraphs revealed that eruptions of plasma are ejected from the Sun. Today, it is known that the Sun produces eruptive flares, filament eruptions, coronal mass ejections and failed eruptions; all thought to be due to a release of energy stored in the coronal magnetic field during its drastic reconfiguration. This review discusses the observations and physical mechanisms behind this eruptive activity, with a view to making an assessment of the current capability of forecasting these events for space weather risk and impact mitigation. Whilst a wealth of observations exist, and detailed models have been developed, there still exists a need to draw these approaches together. In particular more realistic models are encouraged in order to asses the full range of complexity of the solar atmosphere and the criteria for which an eruption is formed. From the observational side, a more detailed understanding of the role of photospheric flows and reconnection is needed in order to identify the evolutionary path that ultimately means a magnetic structure will erupt.


Introduction
Our Sun is a dynamic star, exhibiting a range of large-scale eruptive activity that, over the last roughly 150 years of study, has led to various activity categories being developed. These include eruptive flares, filament eruptions, coronal mass ejections (CMEs) and failed eruptions. The modern understanding is that this list represents different manifestations of a disruption and reconfiguration of the coronal magnetic field, which ultimately leads to an ejection of magnetised plasma that may, or may not, propagate into the heliosphere. In this review we use the general term CME to encapsulate all eruptive activity types, and focus on the physical processes that are involved.
Coronal mass ejections are the sudden release of 10 15 -10 16 g of chromospheric and coronal plasma into the heliosphere. Their speed ranges from between a few tens to a few thousands km s −1 , with an average value of ∼490 km s −1 (Webb and Howard 2012). Fast CMEs (that have speeds greater than the background solar wind through which they are propagating) can drive shocks ahead of them. CMEs are identified through the use of coronagraph data in which they are seen via the detection of Thomson scattered photons in outward moving plasma structures that have a higher electron density than the background corona ( Fig. 1). Their plane-of-sky configuration varies from event to event and includes amorphous blobs, a "3-part" structure consisting of a bright front, dark cavity and bright central core (Illing and Hundhausen 1985), and helical shapes. The plasma traces out the structure of the magnetic field and the ejected plasma is carried away from the Sun as it is "frozen-in" to the magnetic field. Indeed studying the magnetic field is central to understanding CMEs since the energy to power CMEs, which can be up to the order of 10 32 ergs, can only come from the conversion of free magnetic energy into other forms (Forbes 2000). The free magnetic energy is the excess energy stored in the coronal field above the energy of the current-free potential field configuration. The free magnetic energy is that available to be released and is stored in the form of electric currents in the pre-eruptive non-potential magnetic structures (for a review see, e.g., Low 1996;Priest and Forbes 2002;Vršnak 2006;Aulanier 2014;Schmieder et al. 2015).
When CMEs are formed at the solar limb many can be seen to originate from coronal streamers (Hundhausen 1993). However, for the majority of CMEs, the eruption includes a low-lying sheared core field, which may or may not contain filament plasma, and that cannot be observed with a coronagraph. Therefore instruments that image the solar disk must be used to study the origin of CMEs in addition to coronagraph data. It is these low altitude observations that showed CMEs are closely related to solar flares (e.g. Vršnak 2008;Schmieder et al. 2015). During a solar flare, free magnetic energy is converted into energetic particles, heat and electromagnetic radiation and the close temporal and spatial CME-flare relationship shows that they are the product of the same disruption and reconfiguration of the coronal magnetic field (e.g. Temmer et al. 2010, and references therein). A relationship between CME occurrence and flare size has been noted, with observations revealing that the likelihood of a flare having an associated CME increases with flare magnitude and that all flares above X5 in the GOES classification occur in concert with a CME (Andrews 2003;Fig. 1 A coronal mass ejection imaged by the SOHO/LASCO C2 coronagraph on 4th January 2002. The CME is revealed due to Thomson scattered photospheric light and the ejection can be seen in the upper left of the image to have a bright core in addition to the surmounting structure. The black disk is a result of the occulting disk, which has a radius of 1.7R sun ). However, it should be pointed out that CMEs and flares do not occur on a one-to-one basis.
CMEs are common events and their occurrence rate varies with the phase of the solar cycle. On average one CME is formed per day at cycle minimum with around five per day at cycle maximum (Webb and Howard 2012) when active regions are more abundant on the solar disk. Active regions can be the source of CMEs during their emergence phase, when current carrying magnetic field is being brought into the corona (Wang et al. 1994;Leka et al. 1996), and through the decay phase when fragmentation and flux dispersal occur, redistributing the magnetic field into the quiet sun. However, CME characteristics can be very different at different evolutionary stages of an active region. From compact, fast CMEs that occur during the emergence phase to quiescent filament eruptions that occur when active regions are at the end of their lives, and finally streamer blowouts associated with the largest-scale coronal structures/or quiescent filaments. The eruptive structure that forms the CME can come from within an active region and also between active regions. Two examples of how the CME rate varies with an active region's evolutionary phase were presented in  and Green et al. (2002). These studies followed NOAA 7978 and NOAA AR 8100, respectively, for several months and showed that these active regions were CME productive throughout their lifetime (whereas flare activity disappeared as the decay phase progressed). The findings indicate that CMEs do not seem to be so dependent on the magnetic field conditions that are necessary for flares to occur. Flares, especially major flares, occurred when each active region was young and flux emergence had recently taken place. CMEs on the other hand occurred at all evolutionary stages, including during the decay phase when there was no new major flux emergence and the energy density of the magnetic field had reduced. These studies revealed that CMEs are not as dependent on strong field and magnetic complexity as are flares. Still, active regions that have large and complex magnetic field configurations can be the origin of many CMEs during their emergence phase. For example, NOAA active region 9236 produced 17 CMEs in 5 days  and NOAA active region 8100 produced 15 CMEs in 7 days (Green et al. 2002).
However, these large and dramatic regions should not distract from the fact that CMEs are formed in magnetic structures that span a range of spatial scales. The size-scale of the preeruption structure can be as large as 100s thousands of km long (polar crown filaments) but the structure can be as small as only 10s of thousands of km across (Mandrini et al. 2005).
Some CMEs detected in coronagraph data show no signatures in the lower atmosphere, earning them the name "stealth CMEs". They originate at altitudes of >0.1R , where R is the solar radius, (Robbrecht et al. 2009) in magnetic structures where the field strength is weak and the plasma density low, meaning that observational signatures are less likely to be seen (Howard and Harrison 2013). A study of the Sun at solar minimum suggests that stealth CMEs may make up 30% of all CMEs that occur .  found in their study after solar maximum in 2014 that ten out of the 16 studied events could be classed as stealth CMEs.
Taking into account both the remote sensing and in situ data it has been realised that a CME in the inner corona can be composed of some or all of the following: -a shock that moves ahead of the CME -a frontal structure with high plasma density (n e ≈ 10 14 m −3 ), coronal temperature (∼2 MK) and high magnetic field strength (10 −4 T) -a cavity with low plasma density (n e ≈ 10 13 m −3 ), coronal temperature (∼2 MK) and field strength of few 10 −4 T -a prominence core with high plasma density (n e ≈ 10 17 m −3 ), low temperature (∼80 000 K) and high field strength (few 10 −3 T) -a post-eruption arcade with a plasma temperature of ∼10 MK and a few 10 −4 T field strength Studies of CMEs have found a new relevance in the last three decades through the realisation that they can drive severe space weather at Earth (Gosling 1993;Green and Baker 2015). Parameters that are important for determining the level to which a CME will be geo-effective upon its impact with the magnetosphere include the strength of the southward component of the CME's magnetic field (B z ), the speed and plasma density of the CME (V CME , ρ) and the CME's dynamic pressure (ρ V 2 CME ). This means that fast CMEs with a strong and sustained southward magnetic field component and/or a high plasma density are likely to be the ones that have the most significant effect. Historical case studies of such events include the Carrington event of 1859 when a CME reached the Earth in just 17.5 hours (Cliver and Svalgaard 2004) giving it an average velocity along the Sun-Earth line of ≈2300 km s −1 , and the March 1989 event which was exceptionally bright in coronagraph data, implying it had a high plasma density although it also had a much lower average Sun-Earth line speed of around 770 km s −1 (Feynman and Hundhausen 1994). These CMEs produced geomagnetic storms with a disturbance storm time (Dst) index estimated to be −850 nT ) and −548 nT (Cliver and Svalgaard 2004) respectively. It is worth noting that the Dst value during the 1859 event has been a topic of controversy, with Tsurutani et al. (2003) proposing that −1760 nT may have been reached.
Especially important in the creation of severe geomagnetic storms appear to be CMEs that follow one another in quick succession Dumbović et al. 2015;Temmer and Nitta 2015;Vennerstrom et al. 2016). The Carrington event is one example of such a case. In these events, the preceding CME may "pre-condition" the space between the Sun and the Earth so that the following CME is embedded in a faster solar wind stream and so experiences less drag; allowing it to maintain a high velocity out to 1 AU (Temmer et al. 2017).
The impact of a CME extends beyond its direct interaction with the magnetosphere as shocks driven by fast CMEs (roughly above around 800 km s −1 , Kahler 2001) can accelerate particles, producing an additional space weather effect when these solar energetic particles (SEPs) arrive at the Earth (Reames 2013). Meaning that when CME predictions are being developed, speed and spatial extent of the CME are important quantities to predict too. As are the upstream wind conditions ahead of the shock that can affect the shock geometry.
The very significant role that CMEs play in the creation of space weather has naturally led to a desire to predict their time of arrival at Earth. Once a CME has been launched, it can take from less than a day to more than three days before it arrives. Since the fastest CMEs can be the most geo-effective, the ability to predict their occurrence before they happen would significantly improve forecasting lead times. In this review paper we draw together the wide-ranging observational (Sect. 2.1) and theoretical (Sect. 2.2) work on CMEs to assess the current status of understanding how eruptive structures are formed and how energy is injected into the magnetic field, over timescales of weeks to hours before the eruption. We focus on processes that occur during the energy build-up phase that bring the magnetic field to a point where a "driving" mechanism occurs to drive the rapid expansion of the structure. From this, we look at current and future directions of predicting CMEs based on both observations and modelling (Sect. 3).

Observations
Coronal mass ejections are observed using a wide range of imaging and spectroscopic data, spanning wavelengths from radio to X-ray and altitudes from the photosphere to the outer corona. These observations reveal that in many cases, the onset of a CME will be marked by the eruption of a filament/prominence in the lower solar atmosphere. Indeed, over 70% of CMEs begin this way (Munro et al. 1979). The filament may be seen as a bright core in coronagraph data in CMEs that have a 3-part structure. Lower coronal observations of CME source regions can also show the rise of plasma structures that are emitting at EUV and/or soft X-ray wavelengths. Such rising structures have been seen in Yohkoh/SXT data (e.g. Moore et al. 2001;Green and Kliem 2014), SDO/AIA channels (e.g. Liu et al. 2010;Zharkov et al. 2011) and STEREO/EUVI data (Bein et al. 2011;Patsourakos et al. 2013). In the wake of a CME, the reconfiguration results in a post-eruption arcade, which is also known as a flare arcade. When the CME source region is located at the solar limb, hot plasma structures referred to as plasmoids or flux ropes have been observed (Shibata et al. 1995;Cheng et al. 2011;Reeves and Golub 2011). Once the eruption is well underway, the footpoints of the expanding and erupting magnetic volume can sometimes be identified through EUV or soft X-ray regions which have dimmed due to the reduction in plasma density (Rust and Hildner 1976;Sterling and Hudson 1997;Zarro et al. 1999). Figure 2 provides examples of the variety of lower corona signatures that indicate the occurrence of a CME.
In this section we review observational aspects of the coupled evolution of the photosphere and the corona in the long-term (Sect. 2.1.1) and short-term (Sect. 2.1.2) leading up to the eruption and during the rapid acceleration phase of the eruption itself (Sect. 2.1.3).

Long-Term Evolution Prior to Eruption
The energy required to power a CME appears to be built-up in the coronal field over days/weeks prior to the eruption, in an energy storage and release process. This injection of energy is driven by both the emergence of flux (which is non-potential, Leka et al. 1996) and photospheric motions due to convection and differential rotation that can stress the coronal field. In some cases sunspot rotation is observed (e.g., Brown et al. 2003) and/or rotational motions of same-polarity magnetic fragments, one about the other .
Development of non-linear force-free field (NLFFF) modelling techniques have provided methods to reconstruct the coronal magnetic field using photospheric vector magnetograms. In turn the evolution of the energy in the coronal field can be analysed from these extrapolated configurations. The analysis of the magnetically complex NOAA active region 11158 by Sun et al. (2012) followed the energy stored in the coronal field over a period of five days, during which time the region produced multiple flares and CMEs. The study found a rapid increase in the free magnetic energy of the coronal field lasting around four hours as the flux emergence began. This was followed by a more gradual injection of free energy into the coronal configuration during the following ∼two days. The energy was mostly accumulated at low heights (less than 6 Mm) along a polarity inversion line where a low-lying filament formed, and where later eruptions originated. The extrapolated coronal field configuration shows that around the time of an X-class flare approximately 3.4 × 10 31 ergs of energy is released. The NLFFF modelling studies by Gibb et al. (2014), Mackay et al. (2011) followed the evolution of two short-lived small bipolar active regions during their emergence and decay phases using a magnetofrictional method. The simulated coronal field in these studies captured the energy build-up prior to the CMEs that the regions produce. Free energy is found to be built-up during the flux emergence phase and also during the decay phase when the regions start to undergo flux cancellation at their internal polarity inversion lines (flux cancellation is described and discussed at the end of this section).
As well as energetics, the quantity of magnetic helicity has been considered in the evolution of active regions that produce CMEs. Magnetic helicity is the quantity that encapsulates how twisted, sheared and interlinked the magnetic field is. Magnetic helicity is approximately conserved even during non-ideal magnetohydrodynamic (MHD) processes (Berger 1984). For example, magnetic reconnection serves to transport magnetic helicity to larger spatial scales, opposite to energy, which cascades to smaller scales where it can dissipate. Magnetic helicity is injected into the corona via flux emergence, photospheric motions and torsional and shear Alfvén waves and this, along with its conservation property, has led to the suggestion that CMEs may act as a valve to prevent magnetic helicity from endlessly accumulating in the solar atmosphere (Heyvaerts and Priest 1984;Low 1996). Indeed, active regions during their entire lifetime can be the source of many CMEs, which eject a significant amount of magnetic helicity Green et al. 2002). It has even been proposed that CMEs are the unavoidable product of magnetic helicity accumulation in the force-free magnetic field of the corona . The proposal of an intimate link between this topological quantity and a CME spurred a wave of research in this area. New methods to quantify magnetic helicity injection into the solar atmosphere via flux emergence and photospheric flows provided the tools to be able to make these investigations (Démoulin and Berger 2003).
Many studies of the evolution of magnetic helicity in active regions in the time periods around the occurrence of flares and CMEs have been carried out. Nindos and Andrews (2004) studied the coronal magnetic helicity evolution of active regions that were the source of 133 M-class and X-class flares. They found that the magnetic helicity contained in the coronal field of an active region was larger in regions that produced flares that had associated CMEs than for those regions that had CME-less flares. Other studies have looked at magnetic helicity flux, even proposing that a threshold helicity injection rate must be exceeded for a CME to occur (LaBonte et al. 2007). However, in contrast to this, other work has found that major changes in the magnetic helicity content of an active region can occur both before and after a CME (Smyrli et al. 2010). Helicity injection prior to a CME in the Smyrli et al. (2010) study was seen to be due to flux emergence, whereas the changes after the CME were proposed to be due to the torque imbalance between the helicity-depleted corona and the sub-photospheric portions of the flux as suggested in Longcope and Welsch (2000), , Green et al. (2002). Whether or not magnetic helicity is a quantity that can be used to predict CME occurrence will be discussed further in Sect. 3.1.
The role of flux emergence and photospheric shear flows in injecting sufficient helicity into the corona before a CME were studied by Nindos and Zhang (2002) using NOAA active region 9165. The study finds that the horizontal flows did not inject sufficient magnetic helicity into the corona to account for that ejected by CMEs. Instead, flux emergence is invoked as the dominant helicity source. However, studies of smaller and shorter-lived bipolar active regions suggest that during the decay phase of these short-lived regions the photospheric motions caused by convection (super-granular flows) may play an important role in injecting helicity into the corona (Savcheva 2012;Gibb et al. 2014). All the above studies indicate that the magnetic helicity in the corona can be injected via variety of mechanisms, and the mechanism that dominates may depend on the evolutionary stage of the region. Next, we focus on how a ubiquitous process observed in the photosphere known as flux cancellation might prime the corona, ready for an eruption.
Flux cancellation is a common phenomenon that occurs all over the Sun  and is the name given to an observational event consisting of the approach of two opposite polarity fragments that collide and then disappear from the photospheric magnetic field measurement . It's importance for this review paper is that flux cancellation is often observed in the time leading up to a CME , which suggests that it may be important in the creation of an eruptive magnetic field configuration and its destabilisation. The physical interpretation of flux cancellation is that when opposite polarities collide they undergo reconnection in the photosphere or chromosphere (e.g. Litvinenko and Martin 1999). These opposite polarities represent the footpoints of two sheared loops and the low-altitude reconnection reconfigures the magnetic field so that a long loop forms in the corona and also a very short loop that submerges below the photosphere due to its tension force overcoming buoyancy. Shear and free magnetic energy are concentrated along polarity inversion lines where flux cancellation is ongoing for extended periods of time (Welsch 2006), hence it relevance for CME studies.
In some small bipolar active regions, an eruptive configuration is formed via flux cancellation along the internal polarity inversion line during the decay phase of the region, when the flux is fragmenting. This is well observed in bipolar active regions that form S-shaped "sigmoidal" emission structures in soft X-ray or EUV wavebands over a couple of days (e.g. . Both the coronal evolution and the photospheric evolution in these cases match well the scenario proposed by van Ballegooijen and Martens (1989) in which a sheared arcade is transformed into a flux rope. This is also in line with the theoretical findings that sigmoids are created by hot plasma along field lines in the vicinity of a quasiseparatrix layer formed in flux rope topology (Titov and Démoulin 1999).
In summary, the evolution of an active region prior to a CME involves the emergence of magnetic flux into the corona that is then modified by photospheric flows, which eventually disperse the field over an ever wider area. During the lifetime of the active region flux cancellation may be occurring that is able to modify the configuration of the magnetic field and allow free magnetic energy to buildup along polarity inversion lines. CMEs can occur at any point during this evolutionary sequence although the frequency of CMEs is highest at times of new flux emergence (Green et al. 2003).

Short-Term Evolution Prior to Eruption
Collated observations of the evolution of CME source regions in the tens of minutes before an eruption reveal that the magnetic structure evolves through distinct dynamic phases. The majority of CMEs exhibit a slow-rise phase, a phase of rapid acceleration and then a phase where the erupting structure propagates into the heliosphere (Zhang et al. 2001;. The slow rise phase occurs over timescales of minutes to hours and can be observationally identified through the upward motion of a filament or the expansion of coronal loops. The phase of rapid acceleration then occurs and signals that the eruption is in progress. During the propagation phase the CME interacts with the ambient solar wind ). Here we describe the short-term evolution in the hours or minutes leading up to the phase of rapid acceleration.
It has been realised for many years that the magnetic field configuration associated with a sigmoid is highly likely to erupt (Canfield et al. 1999), however an exact observational definition of a sigmoid has not yet been reached. The term is used to describe continuous S-shaped structures and also collections of loops that appear S-shaped overall. In the cases where a continuous S-shaped sigmoid is observed, an eruption follows around within around 14 hours (Green and Kliem 2014). Given the finding that continuous S-shaped sigmoids that form along polarity inversion lines with sustained flux cancellation are indicative of a low altitude flux rope, the data indicate that these flux ropes can form and remain stable on the Sun for several hours. Sigmoids are best observed near disk centre, so that projection effects do not mask their shape. There is a related set of observations at the solar limb that indicate the formation of flux ropes at a higher altitude, formed by coronal reconnection. These flux ropes appear to form from around 7 hours to 10 minutes before their eruption (Patsourakos et al. 2013;Cheng et al. 2011). A recent study identified such a flux rope near disk centre, allowing the photospheric magnetic field evolution to be observed and revealing that the flux rope formation was driven by flux emergence which exhibited rotational motions of the fragments in the leading polarity, one about the other ). These photospheric motions and the corresponding evolution of the coronal field led to magnetic reconnection in the corona and the formation of the flux rope.
As the time of the eruption approaches observations reveal that coronal structures begin to rise, albeit slowly, before the main phase of acceleration. Filaments darken and start their slow ascent from around 3 hours before their successful eruption (Martin 1980). More compact active region structures complete their slow rise phase more rapidly, taking as little as just a few minutes (e.g. Zharkov et al. 2011) before the structure undergoes rapid acceleration.

Evolution of the Eruption in the Lower Corona
Over many decades, observations of the onset and evolution of the eruption of a CME in the lower corona have been drawn together to form the standard model for a CME or an eruptive flare. This model is also known as the CSHKP model after the authors of the early papers (Carmichael 1964;Sturrock 1966;Hirayama 1974;Kopp and Pneuman 1976). For a review of some of the many observational papers that have further supported this model see McKenzie (2002). The standard model of an eruption begins with the presence of a core field, either a sheared arcade or magnetic flux rope, embedded in an arcade field. Such sheared structures and flux ropes are well supported by observations from the photospheric vector magnetic field to EUV and soft X-ray emission structures in the corona. At the start of the eruption, the core field starts to rise, stretching the overlying arcade. A current sheet forms in the arcade below the core field, and reconnection sets in. Observations that support the occurrence of reconnection include the reconfiguration of the coronal field, indications of non-thermal particles from the formation of flare ribbons and hard X-ray emission, and reconnection inflow (McKenzie 2002). These observational signatures indicate that reconnection is closely related in time to the rapid acceleration phase of the CME. Indeed, the reconnection helps facilitate the eruption as it builds poloidal flux around the erupting core field and cuts the tethers of the overlying field. Flare loops form below the reconnection region as the downward product of the reconnection. Some CMEs show a clear rotation of the erupting structure as the eruption begins, indicating the presence of a flux rope and a conversion of twist into writhe (deformation of the axis of the flux rope, Green et al. 2007).
Not all CMEs that are initiated make it to a full eruption (Rust 2003). Instead, their rapid ascension is arrested in the lower corona. This sub-set of events are known as failed eruptions or confined eruptions. The failure of the eruption might be due to the reconnection that occurs below the erupting core field not sufficiently cutting the tethers of the overlying arcade field (Moore et al. 2001) or the field strength in the overlying arcade not dropping off sufficiently rapidly with height . However, failed eruptions still serve a useful purpose in CME studies as they allow an investigation of the triggers versus the drivers of an eruption. In a failed eruption, presumably the trigger occurs, but the physical processes needed to drive the CME into the heliosphere do not (see Sect. 2.2.2).

Theory and Simulations
From a theoretical point of view, there is little doubt that large-scale solar eruptions are magnetically driven, since they originate in the corona where the plasma beta (ratio of gas to magnetic pressure) is small. They are powered by the free magnetic energy that is stored in the corona in non-potential magnetic fields, i.e., in volumetric electric currents. These currents are believed to emerge from the convection zone (e.g., Leka et al. 1996) or be produced by photospheric flows that shear and twist the coronal magnetic field (e.g., Klimchuk and Sturrock 1992). Eruptions occur when the slow, continuous stressing of the coronal field by flux emergence and surface flows reaches some threshold above which magnetic equilibrium cannot be maintained and the field violently erupts (Forbes 2000). While this overall picture is well established, the detailed processes that govern the evolution prior to and during an eruption are not yet fully understood. In this subsection, we thus focus on the following three questions: How is the free energy built up in the corona? (Sect. 2.2.1) How are solar eruptions initiated? (Sect. 2.2.2) What are the physical mechanisms that can propel the magnetic field and plasma into interplanetary space? (Sect. 2.2.3) We will also discuss potential connections between individual eruptions (leading to "sympathetic" events) in Sect. 2.2.4.

Modelling the Build-up of Free Magnetic Energy
Non-potential magnetic fields are at the epicenter of solar eruptive behavior. Eruptions may originate from large-scale relatively weak magnetic fields as in the case of streamer blowouts or quiescent prominence eruptions (Hundhausen 1993;Lynch et al. 2010) or they may originate from strong active region magnetic fields. Whatever the circumstances, solar eruptions originate from strongly non-potential fields where the photospheric magnetic field is nearly parallel to the polarity inversion line (PIL), a so-called magnetic shear (e.g. Forbes 2000;Schrijver 2009), which stands in contrast to a potential field that runs perpendicular to the PIL and has no free energy. There is enormous evidence for the existence of highly sheared magnetic fields associated with CMEs and large flares. At the photosphere, magnetic shear in CME-productive active regions is directly measurable with vector magnetographs (e.g. Hagyard et al. 1984;Zirin and Wang 1993;Falconer et al. 2002;Yang et al. 2004;Sun et al. 2012) and the associated shear flows can be observed along the photospheric PIL (e.g. Strous et al. 1996;Yang et al. 2004;Deng et al. 2006;Sun et al. 2012).
With such observations in mind, MHD simulations of solar eruptions have long invoked shear flows prescribed as boundary conditions to energize magnetic fields to produce largescale solar eruptions. Many examples exist including simulations by Steinolfson (1991), Mikic and Linker (1994), Wolfson (1995), Antiochos et al. (1999), Amari et al. (2003a), Lynch et al. (2008), , which are successful in reproducing many characteristics of CMEs, filament eruptions and flares. CME-like eruptions have also been modelled with prescribed rotational flows that twist up the field (e.g., Amari et al. 1996;Török and Kliem 2003;Aulanier et al. 2005, see also Table 1). A significant limitation of these numerical models is that because of the enormous differences in plasma conditions, they all treat the corona as being disconnected from the solar interior. While this simplification greatly reduces computational costs, it also removes any possibility of elucidating the physical processes that by their very nature transport free magnetic energy and magnetic helicity from the solar interior to the corona.
In contrast, simulations of magnetic flux emerging from the convection zone to the corona set the stage to understand how active region magnetic fields are organized and energized. Global-scale simulations of the solar convection zone have shown how magnetic fields are amplified by the dynamo to produce buoyant loops that capture active region characteristics such as Joy's law and sunspot asymmetries (e.g. Spruit and van Ballegooijen 1982;Abbett et al. 2001;Fan and Fang 2014;Nelson et al. 2014). However, constraints in timescale and numerical resolution prevent these simulations from reaching beyond 95% of the solar radius. To incorporate the upper layers of the solar atmosphere, fully compressible MHD models are applied on small scales in Cartesian domains (e.g. Shibata et al. 1989;Matsumoto et al. 1993) that illustrate magnetic flux buoyantly rising through solar photosphere to the corona. Simulations that followed (Manchester 2001;Fan 2001;Magara and Longcope 2003;Török et al. 2014) found that under these circumstances, shear flows develop spontaneously where horizontal flows reverse direction across emerging loop structures. An analytical model by Manchester and Low (2000) first showed that these shear flows are in response to the Lorentz force, which occurs as magnetic fields rise through the stratified atmosphere and are deformed by the intense pressure gradient of the surrounding plasma. The shearing motions are of the form of large-amplitude shear Alfvén waves where magnetic tension drives the flows that draw the magnetic field along the direction parallel to the polarity inversion line and transport magnetic flux and energy from the convection zone into the corona. Manchester (2003Manchester ( , 2007, Manchester et al. (2004a) found that such shear flows are capable of driving eruptions of magnetic arcades and emerging flux ropes and proposed this mechanism as an energy source and initiation mechanism for CMEs.
Early simulations of flux emergence treated the solar interior as a simple polytropic layer, which roughly captured the temperature structure and density but was free of convection. In recent years, realistic solar models with granular convective and coronal atmospheres have been developed with the inclusion of radiative transfer and heat conduction. Examples include Stein and Nordlund (2006), Abbett (2007) Rempel and Cheung (2014).
Following this work, came realistic simulations of flux emergence and active region formation in a fully convective layer (e.g. Fang et al. 2010Fang et al. , 2012bCheung et al. 2010;Rempel and Cheung 2014). In the convection zone, high-plasma beta magnetic flux is severely distorted as distinct parts of the flux system rise and fall with the large-scale convective flows. As flux passes through the photosphere, it is shredded by the granular convection only to later reorganize as like polarities coalesce into kiloGauss pores and sunspots (Fang et al. 2012b;Cheung et al. 2010;Rempel and Cheung 2014). Even with this complex structure the Lorentz force continues to drive both shear and rotational flows that transport magnetic flux, helicity and free magnetic energy from the convection zone just as found in earlier simulations. In the case of Fang et al. (2012a), it was also found that while rotational flows provide free energy in the area of sunspots, shear flows dominate the energy build-up at the PIL. Fang et al. (2012a) further found that converging motions, flux cancellation and tether-cutting reconnection can combine to build-up the magnetic shear and free energy in the corona necessary for eruptive events. In this case, nearly 40% of the magnetic energy is free energy, which is near the amount required to produce eruptions (Amari et al. 2003a;. Recent observations are proving that the ubiquitous shear and rotational flows so strongly associated with CMEs produced during an active region's emergence phase are driven by the Lorentz force. Analysis of high-resolution vector magnetograms (Georgoulis et al. 2012) indicate that photospheric electric currents generate the Lorentz force that drives shear flows along active region PILs during this phase. Similar analysis by Su et al. (2008) has revealed that the Lorentz force is the driver of rotational flows in sunspots. In these cases, the flows are shown to produce a build-up of magnetic energy with magnetic stress passing from the convection zone into the corona to form non-potential fields. These observational results both verify prior numerical simulations and confirm a fundamental mechanism for energizing coronal magnetic fields before eruption.
Simulation results of the development of shear flows, sheared magnetic fields, Poynting flux and an erupting flux rope are shown in Fig. 3. Panel (a) shows the results of a  Manchester et al. (2004a) shows the eruption of the flux rope resulting from shear flows. At the photosphere, the vertical magnetic field strength is shown in color and the horizontal direction is shown with vectors. Here, the current sheet is shown with a gray iso-surface, forming a sigmoid structure that runs nearly parallel to the highly sheared PIL 2.5 D simulation of flux emergence that shows the development of sheared magnetic fields (Manchester 2001). Panel (b) shows the Poynting flux associated with horizontal flows passing into the corona (Fang et al. 2012a) while Panel (c) shows the development of a highly sheared magnetic field resulting from flux emergence (Török et al. 2014). Panel (d) shows the magnetic field of an erupting magnetic flux rope, where the eruption is triggered by shear flows driven by the Lorentz force.

Initiation Mechanisms of Solar Eruptions
As described in the previous sections, solar eruptions are powered by the free magnetic energy that is stored in current-carrying magnetic fields in the corona. Prior to an eruption, Fig. 4 Non-linear force-free field extrapolation of NOAA AR 10930 before it produced an X3.4 flare and CME (from Schrijver et al. 2008a). Red iso-surfaces show electric currents. Coloured (white) field lines are open (closed). Closed field lines arch over the currents, stabilizing the sheared or twisted core magnetic field the magnetic configuration essentially consists of two parts. One part is the current-carrying (sheared and/or twisted) core field, which is always located above a polarity inversion line and which may or may not harbor a filament. This field naturally seeks to expand and must therefore be stabilized by the other part of the configuration, namely by the largely currentfree (potential) ambient field that surrounds the core field (often referred to as the "strapping field"), as illustrated in Fig. 4. The idea underlying many eruption-models is that an eruption can be initiated by any physical mechanism that is capable of bringing this configuration to a state at which stable equilibria cease to exist and the system violently erupts, thereby releasing some of the stored free energy. The initiation mechanism could be one that either increases the current in the core field or decreases the stabilizing magnetic tension of the ambient field, or both (see, e.g., Aulanier 2014).
However, even though eruptions have been observed extensively for many years, the specific build-up and trigger mechanisms, and the mechanisms that drive CMEs remain rather elusive. This is, to a large extent, due to the persistent lack of routine magnetic field measurements in the corona. Analytical models and numerical simulations provide a useful tool for our attempts to close this gap. These approaches typically employ idealized initial configurations and a reduced set of the MHD equations, both of which are strong simplifications of reality. However, such a simplified approach has two significant advantages. On the one hand, it often allows one to isolate in the model the specific mechanism of interest. On the other hand, it keeps the computational costs moderate and thus provides the possibility of performing extensive parametric studies. Such studies are useful because they allow, for example, an assessment of physical quantities that are difficult or impossible to measure in the corona, such as magnetic twist (e.g., Kliem et al. 2012). In the past decades, idealized analytical models and numerical simulations have been widely used to test and study suggested eruption mechanisms, and have thereby greatly increased our understanding of how solar eruptions work. We note that significantly more complex MHD simulations of eruptions exist as well. These simulations include observed data and more of the physics, and they are aiming to reproduce specific eruptions with as much realism as presently possible (see Manchester et al. in this issue). Table 1 shows a compilation of eruption models of the so-called "storage and release" type, together with references and a compilation of review papers on the subject. 1 Storageand-release models assume that the free energy required to power an eruption is slowly (i.e., quasi-statically) accumulated in the corona over hours or days, and then rapidly released during the eruption. "Directly driven" models, such as thermal blasts during flares or flux injection from the solar interior, appear outdated and are not considered here (for discussions of their validity see, e.g., Forbes 2000, Klimchuk 2001, Linker et al. 2003; but see also Chen 2017). Furthermore, the table does not include eruption scenarios that have, to the best of our knowledge, not yet been tested using numerical simulations. Such scenarios are, for instance, the triggering of eruptions by the merging of magnetic flux ropes (e.g., Pevtsov et al. 1996) or by the unloading of prominence mass (e.g., Low 1996;Seaton et al. 2011), as well as the conjecture that CMEs occur in order to avoid an over-accumulation of magnetic helicity in the corona (e.g., Low 1994).
The mechanisms are grouped according to their predominant role in the eruption process, i.e., whether they act as a "trigger" or as a "driver". 2 By driver we mean here any mechanism that can account for the rapid acceleration and huge expansion of plasma and magnetic field observed during eruptions, i.e., a mechanism that can produce a CME. As a trigger, on the other hand, we consider mechanisms that are capable of initiating an eruption, but can by themselves not produce a CME. The role of a trigger is rather to bring the system to a state at which a driver can take over and complete the eruption.
Such a strict distinction between trigger and driver is somewhat artificial, but we use it here for the following reasons. First, rise profiles of eruptions sometimes clearly show two distinct phases that cannot be fitted with one functional form (see Sect. 2.1.2 and e.g., Sterling et al. 2007a,b;Schrijver et al. 2008b), indicating the occurrence of two distinct physical mechanisms. Second, it has been demonstrated in MHD simulations that certain mechanisms (such as tether-cutting and magnetic breakout) do not seem sufficient for producing a full eruption (e.g., Karpen et al. 2012). Finally, some mechanisms do not appear to be capable of producing a CME, just by their nature. The helical kink instability, for example, is expected to lead to a (relatively moderate) deformation of the shape of an unstable flux rope, but not to its huge expansion.
We note that some of the processes listed as a trigger mechanism in Table 1 (flux emergence, sunspot rotation, flux cancellation) typically act on time-scales of days, much longer than the typical duration of the slow rise phase preceding eruptions. Their preliminary role is to accumulate free energy in the corona (see Sect. 2.2.1). However, as demonstrated by the numerical simulations referenced in the table, they can also, on shorter time scales, induce changes to the system that push it towards an eruption (an example would be flux emergence into an existing pre-eruptive configuration). In that sense it appears justified to consider these processes not just as build-up, but also as trigger mechanisms.
At present, it seems that there may exist only two driving mechanisms in solar eruptions, namely the torus instability and the "flare reconnection" (for a somewhat different view, see Aulanier 2014). The latter is the reconnection that takes place across the vertical current sheet below a CME and produces the flare signatures. The torus instability is an ideal MHD instability, closely related to the processes known as "flux-rope catastrophe" and "loss of equilibrium" (see Kliem et al. 2014a for a detailed discussion). In this instability, the driver responsible for the acceleration of the ejecta is the hoop Table 1 Compilation of suggested physical mechanisms of the "storage-and-release" type for the triggering and driving of solar eruptions. "FE" stands for flux emergence. Note that (i) we do not distinguish here between 2D and 3D simulations; (ii) we do not include simulations that model CMEs by starting with an "out-of-equilibrium" flux-rope configuration (e.g., Manchester et al. 2004b;Tóth et al. 2007;Cohen et al. 2009;Pagano et al. 2013), even though they are compatible with the storage-and-release paradigm; and that (iii) each article is referenced only once, even if the model or scenario it describes involves more than one mechanism. The bottom row contains related review articles Mechanism References
As apparent from Table 1, many more trigger mechanism than driving mechanisms seem to exist. Rather than describing all the trigger mechanisms that have been suggested, we only give a few general remarks here, noting that more detailed explanations can be found in the references compiled in Table 1. First of all, trigger mechanisms should not be thought of as acting in isolation or being mutually exclusive. It is likely that several of them work together, simultaneously or successively, in starting an eruption (e.g., Williams et al. 2005).
Unfortunately it is often difficult, if not impossible, to pin down which particular mechanism or combination of mechanisms triggered a specific eruption. Most mechanisms cannot be observed directly, so indirect evidence such as EUV brightening or the slow rise of a filament must be used. Also, the same observational manifestations may be produced by different mechanisms. This lack of detailed knowledge greatly hampers our capabilities to forecast the onset of eruptions (see Sect. 3), as different trigger mechanisms may have very different onset conditions. Moreover, even if one could predict that a certain trigger mechanism will take place, it would still be hard to predict whether an eruption will produce a CME or remains confined (e.g., Moore et al. 2001;Guo et al. 2010) as we know little about the conditions under which trigger mechanisms can succeed in switching on driving mechanisms. In fact, it seems that, at least in some events, several confined eruptions are required before a full eruption can occur (e.g., Panesar et al. 2015;Chintzoglou et al. 2015;. Overall, observations and simulations have helped us to develop a relatively clear qualitative picture of how a typical solar eruption works. Slow magnetic reconnection acting below (potentially also above) a sheared or twisted core field results in a slow rise of this field. During this slow rise, the core field successively detaches from the photosphere, reconnection occurs at increasingly larger heights in the corona, and the field successively transforms from a structure that was initially probably more like an arcade into a flux rope. Once the rope has been lifted to a sufficiently large height, the torus instability sets in and starts to rapidly accelerate the rope upward. More or less simultaneously, the flare reconnection sets in and supports the rapid acceleration. What is still needed, though, is a better understanding of the onset conditions, efficiency, and interplay of the underlying mechanisms. In the next sub-section, we will use some general considerations and an analytical approach to dig deeper into the physics behind driving mechanisms and the acceleration of CMEs.

Eruption Driving Mechanisms and Physical Background
The free magnetic energy contained in the coronal field can be expressed in the form W = LI 2 , where I is the total electric current flowing through the system and L is the self-inductance that includes information on the entire magnetic-field geometry (e.g., Garren and Chen 1994;Chen 1989;Chen and Krall 2003;Jackson et al. 2004;Žic et al. 2007, and references therein). Depending on the desired complexity of the model, an otherwise intricate coronal magnetic configuration is represented in the analytical approach as a simple line current, a flux rope, a flux rope embedded in a magnetic arcade, etc. In the process of an eruption, the free magnetic energy is explosively released and transformed into kinetic and gravitational potential energy. In addition, if the eruption is accompanied by a flare, part of the energy is consumed for various thermal and non-thermal processes (heating, convective flows, non-thermal particle populations, etc.). Finally, as the speed of the eruption increases, progressively more and more energy is spent for the work against the "aerodynamic" drag. Let's now go on to consider the maximum velocity and acceleration that a CME can attain.
Bearing in mind energy conservation, one can relate the kinetic energy density and the magnetic field energy density as ρv 2 /2 < B 2 /2μ, where ρ, v, B, and μ are the plasma massdensity, flow speed, magnetic field strength, and permeability, respectively. This relation can be rewritten as where v A is the Alfvén speed within the erupting structure. Consequently, it is expected that eruptions from active regions should attain higher speeds than those occurring in quiet-sun regions, which is fully consistent with observations (e.g., Vršnak et al. 2005;Bein et al. 2011, 2012, andreferences therein).
The CME acceleration is driven by the Lorentz force, or formulated in a more simplified way, by the gradient of the magnetic pressure ∇(B 2 /2μ). Taking an order of magnitude estimate for the magnetic pressure as ∇ (B 2 where d is the size of the structure and f L is the Lorentz force per unit volume, the acceleration can be expressed as a = f L /ρ = B 2 /2dρμ, i.e., a ≈ v 2 A /2d. This implies that the acceleration of a CME is limited by the intrinsic Alfvén speed, and that the acceleration should be larger in compact eruptions than in extended ones. Indeed, in eruptions of active-region structures, where the Alfvén speed is high and the source region is rather compact, the acceleration is usually much larger than in eruptions from quiet-sun regions (e.g. eruptions of quiescent prominences, stealth CMEs; see, e.g., Vršnak et al. 2005;Bein et al. 2012;Howard and Harrison 2013). Note that the relation a ≈ v 2 A /2d implies that in the case of very compact sources, d < 100 Mm, accelerations can achieve values on the order of 10 km s −2 , such as is observed in the most impulsive events Bein et al. 2011).
Using these estimates for the maximum velocity and acceleration, one can also roughly estimate the time and distance over which an eruption will be accelerated. The acceleration time can be expressed as t a ≈ v/a, implying t a ≈ d/v A = t A , where t A represents the Alfvén travel time over the eruptive structure. The acceleration, a, speed, v, and the acceleration distance, d a , can be related as a ≈ √ 2ad a , implying d a ≈ v 2 /2a = d, i.e., the acceleration distance should be comparable with the size of the eruptive structure. Indeed, both types of relationships t a ∝ d and d a ∝ d, as well as t a ∝ 1/a, are found in observations Bein et al. 2011). Furthermore, it is well known that eruptions originating from active regions (smaller size, higher Alfvén speed) are more impulsive (stronger acceleration over a shorter acceleration time) than those launched from quiet-sun regions (Zhang and Dere 2006;Bein et al. 2011), consistent with t a ≈ d/v A and a ≈ v 2 A /2d. Characteristic values for active-region and quiescent-prominence eruptions, based on the presented order-of-magnitude relations, are displayed in Table 2.
Let us now consider some general aspects related to the evolution of an eruption. The selfinductance L is primarily dependent on the size d of the current-carrying structure, L ∝ d (Jackson et al. 2004). Thus, as an erupting structure expands, the self-inductance increases, and it can be approximately taken that d is proportional to the heliospheric distance r, i.e., d ∝ r. Under this assumption L ∝ r. On the other hand, the self-inductance of the magnetic structure is by definition related to the electric current as Φ = LI, where Φ is the magnetic flux associated with the current I (see, e.g., Batygin and Toptygin 1962;Žic et al. 2007). In the absence of dissipative processes that could lead to magnetic reconnection, the flux Φ stays preserved, implying that the electric current has to decrease as I ∝ 1/r. This means that the free magnetic energy of the system, W = LI 2 , decreases, being transformed into the kinetic energy and the work done against drag and gravity.
In line with the magnetic origin of an eruption, the CME acceleration is driven by the Lorentz force, or formulated in a simplified way, by the Ampère force F = I × B. Consequently, in the absence of reconnection, the driving force of a CME rapidly decreases with distance from the Sun, because both the electric current I and the magnetic field B protruding through the magnetic structure decrease (for the dependence of the magnetic field on height see, e.g., Dulk and McLean 1978;Gary 2001;Vršnak et al. 2002Vršnak et al. , 2004aLiu 2008). Indeed, observations indicate that the driving force decreases approximately as r −α , where the exponent ranges between α ≈ 1 and 2 (Vršnak , 2006Vršnak et al. 2004c). In other words, the inductive effects and the decrease of the magnetic field with height, act to prevent a very strong acceleration happening over a large time and distance range.
The general physical principles that govern the eruption process described here should be incorporated in any quantitative model, no matter what the specific magnetic configuration is. Most often, the eruptive structure is represented by a semi-toroidal magnetic flux rope that is fixed at both ends in the photosphere. It is worth noting that the flux-rope concept has been in use since before the discovery of CMEs, i.e., in early models of solar flares and eruptive prominences (e.g., Gold and Hoyle 1960;Ellison et al. 1961;Gold 1962;Byrne et al. 1964;Hyder 1967;Anzer and Tandberg-Hanssen 1970). Gradually, flux-rope models developed from very simple ones (e.g., Kuperus and Raadu 1974;Mouschovias and Poland 1978;Sakurai 1976;Anzer 1978;van Tend and Kuperus 1978;Anzer and Pneuman 1982;Pneuman 1980;Vršnak 1984;Steele and Priest 1989) to more complex forms, where various relevant effects have been sequentially taken into account (e.g., Chen 1989;Vršnak 1990;Titov and Démoulin 1999;Török and Kliem 2003). In addition, the development of numerical techniques has provided more demanding studies, beyond those provided by the analytical approach. Now we focus on the implications of the basic physical principles that have just been discussed as applied to the flux rope concept of an eruption. To do so, let us consider the forces acting on a semi-toroidal magnetic flux rope, sketched in Fig. 5a. In such a configuration, several effects have to be considered. First, the curvature of the flux rope axis causes the Lorentz self-force (also known as "hoop force"). The contribution of the selfforce to the equation of motion was firstly estimated by applying the "virtual work" principle (e.g., Shafranov 1966;Garren and Chen 1994). Hereafter, we will follow the procedure proposed by Vršnak (1984) and Vršnak (1990), where the self-force is estimated considering directly the Lorentz force components related to the magnetic pressure gradient and tension (sketched in the upper right and left corners of Fig. 5a).
Next, one has to take into account a diamagnetic effect, which comes from the influence of eddy electric currents induced at the solar "surface". In the simplest way, this effect can be roughly estimated by considering the effect of a line-current located at the height h above the conducting surface, as done by Kuperus and Raadu (1974). Furthermore, the flux rope concept includes an overlying coronal field (B c ), and the influence of this field must be taken into account. The effect depends on the orientation of the background field-in the case of the Kuperus-Raadu configuration (Kuperus and Raadu 1974), the force is directed downwards, whereas in the case of the Kippenhahn-Schlüter configuration (Kippenhahn and Schlüter 1957), the force is directed upwards.
Finally, one has to consider the effect of gravity. Taking into account buoyancy, the acceleration reads a g = −g(ρ − ρ w )/ρ, where g is the local gravitational acceleration and ρ and ρ w represent the density of the CME and the ambient plasma, respectively. Note that in the case ρ ≈ ρ w the gravitation can be neglected, and if ρ < ρ w then a g becomes positive (i.e., acts upwards) (Low et al. 1982).
Let us now consider the stability of a flux-rope and the conditions under which it will erupt. Putting all the previously mentioned effects together leads to an expression for the net acceleration (for details of the derivation see Vršnak 2008Vršnak , 2016: whereR t , Λ, and H represent the major radius of the semi-toroidal flux rope, the length of its axis, and the height of the axis summit, respectively, all normalized with respect to the footpoint half-separation (following the notation presented in Fig. 5a, one can writẽ The parameter X represents the ratio of the poloidal and longitudinal magnetic field components, X = B φ /B , at the surface of the flux rope (Fig. 5a). The first three terms on the right-hand side of (1), written within the square brackets, represent the effects of the magnetic pressure gradient, magnetic tension, and diamagnetic effect, respectively. The fourth term is due to the background field, whereas the last term considers gravity. In (1), we took into account the previously inferred relation I ∝ 1/λ, i.e., we express the current as I = I 0 λ 0 /λ, where I 0 is the value of the axial current when the flux-rope axis is at the height h = d (i.e., R t = h = d, λ = R t π ,R t = H = 1, Λ = π ) and we considered that the mass of the flux rope is M = const. The constant C L = μI 2 0 π/4M determines the overall "strength" of the eruption, and it shows that the CME acceleration is proportional to I 2 0 . Taking into account typical values of M = 10 12 -10 13 kg (e.g., Webb and Howard 2012) and I 0 = 10 10 -10 11 A (e.g., Hofmann et al. 1992;Filippov et al. 2015, and references therein), it can be estimated that the value of C L should typically range between 100 and 1000 m s −2 .
Furthermore, we assumed that B c (background coronal magnetic field surrounding the rope) is uniform over the rope. So the flux conservation B c r t λ = const., implies B c = B c0 r 0t λ 0 /r t λ. We note that in the absence of reconnection the eruptive structure should expand in a self-similar way, i.e., r t ∝ R t (Vršnak 2008), which means B c = B c0 R 0t λ 0 /R t λ = B c0 π/R t Λ. Furthermore, we introduced a constant C c = I 0 B c0 π 3 d/M, where we have taken into account I λ = I 0 λ 0 , i.e., I = I 0 π/Λ. Taking typical values of d = 100 Mm and B c0 = 0.1-1 G (i.e., 10 −5 -10 −4 T), together with the previously used values for I 0 and M, one finds C c could be as large as 5 × 10 4 m s −2 . In the case where the effect of a decrease of the background magnetic field with height is considered, B c0 should be replaced by, e.g., −n t , as is often used in treatments of the torus instability (e.g., Kliem and Török 2006;Olmedo and Zhang 2010, and references therein). Under this approximation, the last term in (1) should be replaced by C c /Λ 2R1+n t .
To solve (1), one has to specify the relationship between the geometrical parameters R t , H , and Λ. Using Fig. 5a, where we have included the angle ζ that represents the inclination of a line connecting a flux-rope footpoint and the center of curvature of the rope axis, one finds the parametric relationshipsR t = 1/ cos ζ , H = (1 + sin ζ )/ cos ζ , and Λ = (π + 2ζ )/ cos ζ . Taking into account that the poloidal-to-axial field ratio can be expressed as X ≡ B φ /B = 2 πr t /Λ, (1) constitutes an explicit expression for the acceleration as a function of the angle ζ , and from that, one can establish a functional dependence a(H ).
Various solutions to (1) are given in Fig. 6. Several a(H ) curves are shown, following an increasing C L ( Fig. 6a and b) and X ( Fig. 6c and d) sequence. The increase of C L and/or X can be caused by various evolutionary processes related to the photospheric and subphotospheric motions. For example, twisting motions at the flux-rope footpoints would lead to a transport of poloidal flux into the rope (thus increasing X) and increase the electric current (thus increasing C L ). A similar effect may also be expected in the case of reconnection below the rope, which tends to increase the poloidal flux in the rope (see e.g., Vršnak 2016). The value of C L could also be increased by emerging flux, since the increase of the total magnetic flux encircled by the current circuit associated with the rope would lead to an enhancement of the current. Finally, C L might be increased by a mass leakage through the flux-rope legs, since C L ∝ 1/M.
The intercepts of the a(H ) curves with the x-axis are equilibrium points (a = 0). If the slope of the a(H ) curve, ∂a/∂H ≡ ω 2 , at the equilibrium point is characterized by ω 2 < 0 the equilibrium is stable (points denoted in Fig. 6b as H 1 , H 2 , H 3 , and H 4 ). If the rope is moved from this position, it will oscillate around it, with a period of P = 2π/ω. Note that from the observational point of view, oscillations of coronal structures are especially well visualized in the case of prominence oscillations initiated by disturbances coming from distant eruptions, a phenomenon being reported as "winking filaments" a half of century ago by Ramsey and Smith (1966). Inspecting Figs. 6b and 6d, one finds that as the value of C L or X increases, the equilibrium position shifts to larger heights (see the insets in Figs. 6b and 6d) and the oscillation period increases (the slope ∂a/∂H becomes less steep).
The x-axis intercepts of a(H ) curves characterized by ∂a/∂H > 0 represent unstable equilibria. If the rope is displaced from a stable equilibrium to an unstable one (e.g., by a disturbance from a distant eruption), it will erupt. Thus, ropes that are characterized by a(H ) curves that have both types of x-axis intercepts are in a metastable state. Note that with increasing C L or X, the distance between stable and unstable equilibria decreases, i.e., it becomes easier to push the rope from the stable to the unstable point and trigger the eruption. Finally, it should be noted that the a(H ) curve achieves the metastable form only for a sufficiently large value of X. This is illustrated in Fig. 6d, where the a(H ) curve for X = 2 (black curve) is shown, which is characterized by only one (stable) equilibrium point. The curve calculated for X = 2.5 (green curve) already contains the unstable equilibrium at H ≈ 2. Thus, a value of X at which the magnetic structure transforms from stable to metastable configuration lies between 2 and 2.5; a more detailed analysis shows that for the considered combination of parameters (C L = 144 m s −2 , r t /R t = 0.1, C c = 1.3 × 10 4 m s −2 ) this transition occurs at X = 2.2.
The increase of C L or X leads to another essential physical situation. As the a(H ) curves gradually shift "upwards", at a given critical value of C * L or X * , the stable and unstable equilibrium points merge (see the red curve in Fig. 6), after which there is no equilibrium possible. The local minimum in the a(H ) curve just touches the x-axis (in Fig. 6b this critical height is denoted as H * ), i.e., the whole a(H ) curve is located in the a ≥ 0 region. Consequently, the rope erupts, dynamically trying to find a new equilibrium. In the case when the gravity is neglected (C g = 0) a new equilibrium does not exist, and the rope erupts. Note that the critical value C * L is lower for a higher value of X. This is illustrated in Fig. 6c, where the loss-of-equilibrium curve for X = 5 (full-red line; C * L = 144 m s −2 ) is compared with the analogous curve for X = 6 (dotted-red line; C * L = 143 m s −2 ). An example of a metastable flux rope containing filament plasma, gradually evolving towards the loss-ofequilibrium stage and eruption, is presented in Fig. 5b. Filament F1 was slowly rising until reaching a critical height, after which it erupted (Fig. 5c). Besides causing a flare below it, the eruption destabilized and triggered the eruption of the nearby metastable flux rope that was embedding the filament F2 (for details see Miklenic et al. 2009).
For a flux rope characterized by a given value of the poloidal-to-axial field ratio X, a critical value C * L at which the eruption starts depends on the relative contribution of the background magnetic field and the gravity. In other words, it depends on the values of C c and C g , which is illustrated in Figs. 7a and 7b, respectively.
In Fig. 7a we present four loss-of-equilibrium a(H ) curves calculated for different values of C c > 0 and the corresponding values of C * L , applying fixed values X = 5 and C g = 0. The graph illustrates that the peak acceleration a m is larger for a larger value of C c , i.e., larger C * L , and that a m → 0 for C c → 0. The value of a m increases if larger values of X or lower values of C g are applied. Note that in all cases shown in Fig. 7a the eruption starts at the same height, very close to H = 1. This is consistent with conclusions drawn by (e.g., Chen 1989;Vršnak 1990) and confirmed by observations (e.g., Vršnak et al. 1991;Chen and Krall 2003). This property is a consequence of the hoop-force behaviour, achieving a maximum at H = 1, since at this height the radius of curvature attains minimum,R t = 1.
For comparison, besides the four a(H ) curves based on the Kuperus-Raadu configuration (C c > 0) shown in Fig. 7a, one example of the Kippenhahn-Schlüter case is presented (gray-dotted curve; C c = −700 m s −2 ). It reveals a completely different behaviour than the Kuperus-Raadu configuration. First, the equilibrium point is shifted to much lower heights than in the Kuperus-Raadu case. Furthermore, the metastable form of the a(H ) curve is possible only for C g > 0, and in that case there is always an upper-equilibrium position present (the x-axis intercept of the a(H ) curve located at H ≈ 1 in the presented case). Thus, the Kippenhahn-Schlüter configuration can result only in failed eruptions, and will not be considered hereafter.
To illustrate the effect of gravity, Fig. 7b shows four loss-of-equilibrium a(H ) curves for different values of C g and the corresponding values of C * L , applying fixed values d = 70 Mm, X = 5, and C c = 1.3 km s −2 . The three curves calculated with C g = 0 show an upper equilibrium point, unlike the curve calculated with C g = 0. Thus, in C g = 0 cases, the eruption should stop at an upper equilibrium. As a matter of fact, after loosing equilibrium, the rope should accelerate upwards, reach maximum speed at the upper equilibrium, after which it should be decelerated until being stopped. Then, the restoring force at the upper equilibrium would start accelerating it downwards-in the absence of drag the rope should oscillate around the upper equilibrium position. Taking into account the drag effect, the oscillations would be damped until the rope settles at the upper equilibrium (for observations see, e.g., Mrozek 2008).
In Fig. 7b we also show the case C g = 0.3 m s −2 with C * L = 144 m s −2 (red-dotted curve), to compare it with the curve calculated for the same C * L = 144 m s −2 and C g = 0 (red-full curve). The graph shows that the gravity effect shifts the a(H ) curve downwards, and since the effect is weaker at larger heights due to decreasing g, this leads to the formation of the upper equilibrium position. In this respect, let us also note that the gravity effect weakens with increasing d, i.e., it is smaller for larger ropes. For example, for d two times larger than the applied value of d = 70 Mm, the upper equilibrium disappears in the C g = 0.3 m s −2 case.
The effect of the decreasing background magnetic field is illustrated in Fig. 7c, where we show a(H ) curves calculated using again C * L = 144 m s −2 , C c = 1.3 km s −2 , X = 5, and C g = 0, but now applying several different values of the exponent n in the dependence B c = B 0cR −n t . The red curve, representing the case n = 0, is the same as the red curves in Figs. 7a and 7b. The graph shows that the peak acceleration increases with increasing n, i.e., it is larger for eruptions occurring in coronal regions where the decrease of the magnetic field is steeper. The rate of decrease is dependent on the flux distribution at the photosphere. Furthermore, the acceleration peak systematically shifts to lower heights.
Although the effects presented here are similar to those governing the torus instability, they do not describe the same physical mechanism. If one would like to consider the torusinstability effect, the pressure-gradient term in the hoop force, i.e., the first term in (1), should be replaced by a new term, estimated using a new difference between the magnetic field strength at the inner and outer boundary of the torus. The difference can be estimated from the derivative of B c = B 0cR −n t , and from that one gets the pressure-gradient term of the form 2 nC L /ΛR n+1 t . Note that this approach is different from that usually applied (e.g., Kliem and Török 2006), where the equation of motion is derived employing the virtualwork principle. Consequently, the equation of motion is of a somewhat different form than that used in torus instability studies. For example, the "standard" torus instability equation of motion includes explicitly the self-inductance, whereas it is herein included implicitly. In addition, (1) includes the diamagnetic effect, Ampére force, and gravity, which are not usually considered in torus instability studies.
The torus instability effect is illustrated in Fig. 7d, where a(H ) curves for several values of n are presented. Comparing this graph with the graphs presented in Fig. 7a-c, one finds a considerably different behaviour of the flux rope. First, the loss-of-equilibrium height H * is increasing with decreasing n, and at the same time, the peak acceleration a m decreases.
Second, in the case of n > 1 the a(H ) profile is more sharply peaked than profiles shown in Fig. 7a-c, and the peak accelerations are higher. Thus, eruptions driven by the torus instability should be more impulsive than those driven by the hoop force that represents a form of the kink instability. Finally, note that an eruption is not possible if the value of n is too low, due to the trend a m → 0 and H * → ∞ for decreasing n. For example, for the set of parameters used for Fig. 7d, the critical value of n is ≈0.8.
The analysis presented here shows that under the conditions considered it is difficult to get accelerations higher than 1000 m s −2 . For example, even taking C L as high as 1000 m s −2 in the hoop-force option (at the same time requiring extremely large C c = 100 km s −2 ), one gets a peak acceleration up to 600 m s −2 . This limitation is related to the conservation of magnetic flux associated with the electric current in the eruptive structure, in conjunction with the inductivity effects. As previously demonstrated, these effects cause a decrease of the current as the size of an eruption increases, which in turn leads to a rapid decrease of the driving Lorentz force. Thus, in order to get a really high acceleration, the flux-conservation constraint has to be broken. This can be attributed either to an abundant flux emergence, as proposed by Chen (1989), or by eruption-related reconnection resulting in the flare energy release, as suggested by Vršnak (2008). Since in the former case the emergence rate should be much higher than ever observed, the latter option seems to be more appropriate (Kunkel and Chen 2010). It is also important to note that since the open field configuration holds more energy than the closed state, it has been postulated that magnetic reconnection is necessary for a successful eruption to occur (Aly 1984;Sturrock 1991).
The role of reconnection in the dynamics of a flux rope eruption was analysed by Vršnak (2008) and Vršnak (2016), providing an insight into the physics of the CME/flare relationship. Reconnection results in two effects that can significantly enhance and prolong the CME acceleration. First, reconnection supplies the rope with additional poloidal flux (see Fig. 4 in Vršnak 2008), which enhances the hoop force, providing higher peak accelerations. Furthermore, the poloidal-flux supply reduces the inductive decrease of the axial current in the later phases of the eruption, i.e., the acceleration phase is prolonged. The second important effect of reconnection is the weakening of the magnetic tension of the arcade field overlying the flux rope by detaching this field from the photosphere.
In Vršnak (2008) and Vršnak (2016) the reconnection effects were analysed assuming that reconnection starts simultaneously with the force imbalance, attains maximum rate when the rope is at a certain height, and ends at twice this height. Various combinations of the height range and peak reconnection rate were considered to find out how these parameters affect the eruption kinematics (see Fig. 2 in Vršnak 2016). It was demonstrated that the acceleration time profile is tightly synchronized with the reconnection-rate time profile, whereas the velocity time profile is synchronized with the cumulative reconnected flux. This explains the observed synchronization of the CME acceleration phase and the impulsive phase of the associated flare (Kahler et al. 1988;Neupert et al. 2001;Zhang et al. 2001Zhang et al. , 2004Shanmugaraju et al. 2003;Vršnak et al. 2004b;Zhang and Dere 2006;Temmer et al. 2008Temmer et al. , 2010. It was also shown that the peak acceleration is higher for higher peak reconnection rate. Consequently the maximum speed of the eruption is higher for larger total reconnected flux, consistent with observations (Qiu and Yurchyshyn 2005). It was also demonstrated that the highest observed accelerations, reaching ≈10 km s −2 , could be achieved if the total reconnected flux is a few times larger than the initial poloidal flux of the rope, corresponding to ≈10 22 Mx. Such amounts are frequently observed in flares (e.g., Qiu et al. 2004Qiu et al. , 2007Qiu and Yurchyshyn 2005;Miklenic et al. 2007Miklenic et al. , 2009. The sequence of events presented here can be brought together in a relatively simple eruption scenario. If the coronal structure contains a flux rope and this rope is characterized by weak electric current and/or low degree of twist (i.e., low values of C L and X, respectively), it is stable. If the rope is displaced from the equilibrium position, it will oscillate around the equilibrium position. If C L or X (or both) start to increase due to photospheric twisting, shearing, or flux emergence, the structure gradually evolves through a series of quasi-equilibrium states. The equilibrium height in the a(H ) curve increases, which is reflected in a slow rise of the structure. At a given amount of twist, the structure becomes metastable, and if pushed strongly enough to reach the point of unstable equilibrium, it will erupt.
If the structure is sufficiently sheared or twisted, the gradual rise of the structure caused by increasing C L and X ends up with a loss of equilibrium, and consequently, an eruption. The threshold value of X (the ratio of poloidal to axial field) is lower for structures characterized by larger C L . It is important to note that observations of the flux swept out by flare ribbons, compared to in situ flux measurements of the associated magnetic cloud, indicate that the pre-eruption structure is weakly twisted, consistent with EUV and soft X-ray observations. The poloidal flux is added via the flare reconnection (Qiu et al. 2007;Hu et al. 2014). Under certain circumstances, when the a(H ) curve shows the presence of an upper equilibrium, the eruption might be halted, showing damped oscillations around the upper equilibrium position. Otherwise, the eruption continues into the heliosphere. The acceleration is stronger, and speed higher, for larger values of C L and X. Since C L is related to the current (C L ∝ I 2 ) and X with the self-inductance (L ∝ X), this is equivalent to saying that the acceleration and speed are higher for structures containing larger amount of free energy W = LI 2 . Furthermore, since C L can be also expressed as C L ∝ B 2 /ρ, we find that acceleration and speed are higher for structures characterized by a higher Alfvén speed. This is consistent with the order-of-magnitude relationships anticipated at the beginning of this section.
To conclude this sub-section, the analysis presented demonstrates that events characterized by a max > 1 km s −2 can occur only if the eruption is accompanied by strong magnetic reconnection below the erupting rope. Indeed, the most powerful/impulsive eruptions are associated with powerful/impulsive flares, supporting such a conclusion. Note that the reconnection process can go on for a long time, on the order of a day, as evidenced by growing post-eruption loop systems (e.g., , and references therein) and signatures of current sheets in the wake of the CME (e.g., Webb et al. 2003;Lin et al. 2005;Ciaravella and Raymond 2008;). Finally, the acceleration can be additionally enhanced by a transport of twist into the expanded summit of the rope (Jockers 1978;Browning and Priest 1983) and by writhing of the flux-rope axis that decreases the radius of curvature (e.g., ).

Sympathetic Eruptions
The previous discussions in this section focussed on single, isolated eruptions. However, solar eruptions can occur relatively close to each other in both space and time, and may thus interact low in the corona or later on in interplanetary space. Such interactions are discussed in Manchester et al. in this Volume; see also the recent review by Lugaz et al. (2016). Here we focus on a special type of interaction, namely those in which one eruption appears to initiate another one. Such events were first mentioned more than 80 years ago (Richardson 1936), and are now commonly referred to as sympathetic eruptions.
Sympathetic eruptions occur within a relatively short period of time (tens of minutes to several hours) at different locations on the Sun. The individual eruptions may originate in one, typically large and complex source region (e.g., Liu et al. 2009) or in different source regions that are often, but not always, adjacent. 3 Research on sympathetic eruptions revolves around two main questions. Is the close timing of individual eruptions purely coincidental, or does some causal link exist between them? And if it does, what is the physical nature of such a link?
Prior to the launch of SDO, research on the timing of flare/CME events focused mainly on the flares (Fritzova-Svestkova et al. 1976;Pearce and Harrison 1990;Bumba and Klvana 1993;Bagalá et al. 2000;Biesecker and Thompson 2000;Moon et al. 2002;Wheatland and Craig 2006). Many of these studies considered waiting-times (i.e., the time interval between two successive events) and the location of the event, both of which are easier to establish for flares than for CMEs . These predominantly statistical studies have provided strong evidence that a causal relationship between sympathetic eruptions does indeed exist. This was supported by detailed case studies (e.g., Bumba and Klvana 1993;. Sympathetic CMEs were studied less frequently though (e.g., Simnett and Hudson 1997;Cheng et al. 2005;Jiang et al. 2008), and the evidence for their existence appeared to be weaker than for flares . However, since flares and CMEs are often merely different observational manifestations of the same eruption process, this discrimination is somewhat artificial. The physical connections between sympathetic eruptions were attributed to, e.g., large-scale convective motions that led to simultaneous flux emergence in the flaring regions (Bumba and Klvana 1993) or waves triggered by a previous eruption (e.g., Ramsey andSmith 1966, Khan andHudson 2000; see also Mullan 1976 for sympathetic stellar flares). Later it was suggested that these connections are of a magnetic nature, such as surges or jets produced by an adjacent eruption Jiang et al. 2008) or changes in the background field due to expansion or reconnection (e.g., DeVore and Antiochos 2005; Jiang et al. 2011;Yang et al. 2012).
The launch of SDO spurred a new series of investigations. This was triggered by the complex sympathetic eruptions that took place on 1-2 August 2010 and that were observed in unprecedented detail by AIA. Schrijver and Title (2011) analyzed the coronal magnetic field during the time period of the eruptions (using potential field extrapolations), and found evidence for connections between all source regions involved via structural features such as separators and quasi-separatrix layers (QSLs; Priest and Forbes 1992;Démoulin et al. 1996). They concluded that the structural properties of the large-scale coronal field play an important role in the genesis of sympathetic eruptions (see also Schrijver et al. 2013). Other cases of sympathetic eruptions observed by AIA can be found in Shen et al. (2012Shen et al. ( , 2013, Joshi et al. (2016), and .
The 1-2 August 2010 events also inspired new modelling efforts (for earlier attempts see, e.g., Ding et al. 2006;Peng and Hu 2007). Figure 8 shows an MHD simulation by , which was designed to qualitatively reproduce a sub-set of these events (see also Lynch and Edmondson 2013). The simulation suggests two purely magnetic mechanisms for the links between successive eruptions. Figure 8 panel (a) shows the initial configuration, consisting of two flux ropes (F2, F3) located in the lobes of a pseudostreamer and one rope (F1) next to it. Pseudostreamers or pseudostreamer-like configurations appear to be prone to producing sympathetic eruptions; they often harbor twin-filaments that erupt successively (Panasenco and Velli 2010), which was the case also in the 1-2 August 2010 events. In the simulation, the eruption of F1 triggers the formation of a thin current  . The eruption of F1 triggers the subsequent eruptions of F2 and F3 due to the reconnection events R1 and R2, respectively layer around the separator line above the pseudostreamer lobes, across this current sheet breakout-like reconnection (R1) transfers flux from the right lobe to the left one (panel (b)). This reduction of flux, and thus of stabilizing magnetic tension, triggers the torus instability in F2, which subsequently erupts (panel (c)). Finally, the flare-reconnection (R2) below F2 redistributes flux from the left lobe to the growing flare arcade (panel (d)). This new flux reduction now destabilizes F3, completing the sequence of eruptions. Note that these two similar mechanisms can work independently, and their applicability is not restricted to pseudostreamers. They can work also in a quadrupolar configuration, as hypothesized by DeVore and Antiochos (2005). What is required for their occurrence is merely a flux system that contains adjacent lobes overlaid by a separator line and a sufficiently strong perturbation to set things in motion. An important and yet largely unexplored question is under which conditions such reconnections will be actually strong enough to trigger a subsequent eruption (see ). Very recently, Jin et al. (2016) undertook the first steps in answering this question, by investigating and quantifying the effects CMEs can have on adjacent and remote magnetic fields.
The simulation shown in Fig. 8 constituted a significant step forward in understanding the physical links between sympathetic eruptions. However, it was designed as a "proof of concept" model, using a simplified magnetic configuration. Employing instead a potential field source surface (PFSS) model based on observed magnetograms,  performed a detailed topological analysis of the coronal magnetic field during the time of the 1-2 August 2010 events, and found support for the occurrence of the mechanisms suggested in . In order to further substantiate this scenario, MHD simulations with more realistic magnetic fields need to be employed to test whether sympathetic eruptions can occur also under less idealized and symmetric conditions (Mikic et al. 2011). Such simulations will allow us also to investigate the long-range effects of eruptions (e.g., Schrijver and Higgins 2015), as recently demonstrated by Jin et al. (2016). Furthermore, whether other mechanisms, such as waves, can serve as a possible link between eruptions can be tested too. At present, the observational and numerical studies that have been performed since the launch of SDO strongly support the conjecture that the mechanisms at work in sympathetic eruptions are magnetic in nature. These mechanisms, by means of reconnection, act predominantly by reducing the magnetic tension in the source region of an eruption, allowing the core field of the region to erupt.

Predictability of Solar Eruptions
Significant progress has been made in understanding the physics that underlies CME occurrence, as well as the observational characteristics of eruptions. So, where do we stand in our ability to predict CMEs for use in space weather forecasting? There are two basic goals and approaches: (1) to predict the occurrence and timing of an eruption (before it has happened) and (2) to predict the erupting magnetic field configuration that leaves the Sun, and its arrival time at Earth which is related to the impact of an eruption. In this section, we summarise the current status of our ability to predict the occurrence of a CME based on observations (Sect. 3.1) and based on MHD modelling (Sect. 3.2).

Predictions Based on Observations
The evolution of the coronal and chromospheric plasma structures in CME source regions can sometimes provide an indication of an impending eruption, as summarised in Green and Baker (2015). However, the timescales over which observations reveal that a CME might be imminent can vary. For example, filaments darken, broaden and slowly rise in the minutes prior to their eruption (for a review see Martin 1980) whereas streamers brighten over around a day before they lift off as a CME (Illing and Hundhausen 1986). Ultimately, a variety of observations may need to be utilised to both identify that a magnetic structure is evolving toward an eruptive configuration, and give a warning of when the eruption is imminent. Here we summarise the large body of observational work that is relevant to CME prediction.
The observational approach to predicting CMEs broadly involves two aspects; identifying the existence of a potentially eruptive magnetic field configuration and identifying the path to its destabilisation. Perhaps the longest-known eruptive magnetic field configuration is that which supports a filament. Feynman and Martin (1995) investigated the role that flux emergence may play in the destabilisation of filaments and found that 19/30 of the eruptions studied were observed to have new flux emergence in the vicinity of the filament. The flux emerged in a time period of a few days before the CME. Filaments that did not have nearby flux emergence did not erupt. The study showed that the orientation of the newly emerging flux was important, with bipoles that were orientated in such a way as to be favourable for reconnection being key for CME occurrence.
Sigmoids are another coronal configuration that is associated with eruptive activity. Sigmoids are S-shaped coronal emission structures seen in soft X-ray and EUV observations (see Sect. 2.1.1), which form in a sub-set of CME source regions. These configurations are one of the most successful indicators of an eruption (Canfield et al. 1999) and have become synonymous with CME occurrence. Given the observation that sigmoidal active regions are highly likely to produce a CME, the question can be asked "how much time elapses between the sigmoid formation and the CME?". Green and Kliem (2014) found that around 5 to 14 hours elapsed between formation and eruption in the small sample of sigmoid regions studied. This study used the time of the first observation of a continuous S-shaped emission structure (rather than an overall S shape formed in a collection of loops) to indicate that a sigmoid had formed. Canfield et al. (2007) found 107 active regions with overall sigmoidal shape in a study of all the active regions imaged in partial-frame SXT data. This means that around 6% of active regions observed in the partial frame Yohkoh/SXT data during 1991 to 2001 were sigmoidal. Although Sterling (2000) has a slightly more optimistic finding, stating that "These observations indicate that the pre-eruption sigmoid patterns are more prominent in SXRs than in EUV, and that sigmoid precursors are present in over 50% of CMEs." So, whilst sigmoids are a useful proxy, not all active regions evolve to produce a sigmoidal configuration. In addition, sigmoids are more likely to be seen in active regions when they are near disk centre and viewed from above, when projection effects do not mask the S shape. As discussed in Sects. 2.1.1 and 2.1.2, a sigmoidal emission structure can indicate the presence of a flux rope magnetic field configuration.
The observation of a continuous sigmoid therefore gives a only few hours warning that the region will produce a CME. Sigmoids can form in decaying bipolar active regions though, that appear to show a systematic evolution toward this shape during their lifetime and this can be exploited to give a longer lead time for CME prediction. For sigmoids that form in decaying bipolar active regions, the process of formation takes a couple of days. During this time the coronal structures evolve due to flux cancellation in the photosphere. Flux cancellation can take the coronal field through three observational stages; 1) increase of shear along the polarity inversion line, 2) the formation of two sets of J-shaped loops and 3) the formation of a sigmoid Kliem 2009, 2014). Accurately identifying stage 2 would increase the CME forecast time by up to a day. Sigmoids are relatively easy to identify when present in small and bipolar regions and viewed from above. They could be harder to identify in complex multipolar active regions with several polarity inversion lines and strong EUV and soft X-ray emission from adjacent structures. However, it is still possible to observe sigmoids in these environments (e.g. Zharkov et al. 2011) where they also seem to be associated with flux cancellation.
The interpretation that some sigmoids indicate a flux rope topology has led to studies of flux cancellation that probe the amount of flux that builds into the rope. This can potentially be used to investigate the stability of the structure. The flux cancellation associated with the formation of a sigmoid was studied in detail for the small bipolar NOAA active region 10977 . The amount of flux cancelled was measured from the time that the active region had fully emerged to the time of the first eruption from the region. The results showed that at the time of the eruption the ratio of flux in the rope to that in the overlying field of the active region had an upper value of 1:1.5. Another study, in an active region that wasn't sigmoidal found that at the time of a filament eruption the ratio of the flux contained in the rope to the flux of the overlying arcade field was 1:0.9 (Yardley et al. 2016). However, there were previous eruptions from this region and the study wasn't able to track the flux cancellation during the entire time during which the filament was forming, and so may represent a lower limit to the flux that was actually contained in the rope. These observational results can be compared to a series of studies that use the flux rope insertion method (van Ballegooijen 2004) to investigate the stability of flux ropes embedded in an overlying arcade. Stable flux ropes have been produced that have a ratio of flux in the rope to flux in the overlying arcade ranging from 1.15 to 1.9 ( 2008;Savcheva and van Ballegooijen 2009;Su et al. 2009;Savcheva 2012). The Savcheva (2012) study compared the modelled results to the observational deductions of flux rope flux, which found that around 60% to 100% of the flux involved in the flux cancellation built into the rope.
Flux cancellation at the photosphere may be an indicator of the formation of an eruptive configuration but cancellation is a ubiquitous process that takes place all over the photosphere. The challenge is to recognise when and where this process is important for CME occurrence, which can only be done when studied in concert with the evolution of the coronal field. Another aspect is that a small number of sigmoids appear more transiently in active regions. These sigmoids appear to be formed via reconnection in the corona (rather than via flux cancellation and reconnection in the photosphere or chromosphere). Sigmoids formed in this way have been observed around 2 hours prior to their eruption (Cheng et al. 2014; although they tend to flash into view and fade, only to reappear again around eruption onset making them harder to use for CME forecasting. An alternative and complementary approach is to study the occurrence and predictability of solar flares in order to make a CME forecast. This approach is relevant due to the relationship between flares and CMEs (see Sect. 2.1.3) and has the added advantage that the datasets used to study flares include observations of the photosphere, which are available over a much longer time period than X-ray and EUV imaging, and coronagraph observations are. Sammis et al. (2000) analysed eight years of active region observations in order to establish whether a relation between an active region's magnetic classification and the occurrence of large flares (for a description of the active region classes see Hale et al. 1919). They find that active regions classified as βγ δ produce significantly more large flares than other regions of comparable size. Each active region of magnetic class βγ δ and a size greater than 1000 µh revealed a probability of nearly 40% of producing flares of GOES class X1 or greater. Flares of this size have around a 90% chance of being associated to a CME . More recently, surface magnetic parameters were combined with the amount of stored (free) magnetic energy reconstructed from local non-linear force-free field modelling, and tested against the flare productivity of active regions . Though this study revealed that the free magnetic energy clearly differentiates flaring active regions from non-flaring active regions, the temporal variation of an region's free magnetic energy did not exhibit a clear pre-flare pattern.
Studies that apply the physical knowledge of an active region's morphology and magnetic properties, to parameterize the solar data in order to quantify flare/CME productivity, mostly use the photospheric line-of-sight magnetic field. However, vector field measurements can also be used. For the parameterizations based on the magnetic properties of active regions, one can further distinguish between those that relate directly to the observed photospheric field (such as magnetic flux, length and strength of main polarity inversion line, etc.) and those that are derived after reconstruction of the coronal field distribution from the vector photospheric data (such as total magnetic energy, free magnetic energy, etc.), see Barnes et al. (2016). In general, it seems that quantities that are integrated over active region scales are stronger determinants for upcoming flare/CME activity than are localized variations on small scales (e.g. review by Schrijver 2009). In addition, the flaring history of an active region plays a decisive role ). There exist statistical methods that predict flaring probabilities using solely the flaring history of the active region under study. This type of statistical model makes use of the information that the past history of flare occurrence in an active region is an important indicator of its future flare productivity and that the solar flare size distribution follows a power-law (Wheatland 2004).
One of the key features used to describe the magnetic characteristics of an active region with respect to its flare/CME productivity is the polarity-inversion line. Schrijver (2007) conducted a statistical study of the magnetic properties of active regions associated with almost 300 M-and X-class flares, in comparison with 2,500 randomly selected active regions.
The main finding of the study was that active regions that are the source of large flares all have a pronounced polarity inversion line with a strong gradient across it. Calculating the unsigned flux within 15 Mm of such strong field polarity inversion lines is shown to be a good indicator for the occurrence of a large flare.
A series of papers led by Falconer used the polarity inversion line along with a photospheric proxy for the free magnetic energy stored in the coronal magnetic field to investigate whether there exists a relation to the region's CME productivity. The work used vector magnetic field data to derive various parameters that describe the global non-potentiality of the active region's magnetic field. These parameters are the length of the strongly sheared field across the polarity inversion line, the net electric current and the degree of twist. They showed that all three parameters are significantly correlated with the active region's CME productivity (Falconer et al. 2002). Note that in their study, they focussed on active regions with one well-defined, dominant polarity-inversion line. In a follow up study, Falconer et al. (2006) extended this work to a larger sample of active regions, and also extended the set of parameters derived from the vector field data. Three of the six parameters calculated characterize the total non-potentiality of the active region (length of the strong-shear segment of polarity inversion line, length of the main polarity inversion line with a strong field-gradient across it, and the total net electric current). Two of them describe the overall twist in the active region (the net-current α and best constant α derived from a linear force-free field) using the twist parameter α, which is determined from the net current divided by the magnetic flux. One parameter describes the active region size (the total magnetic flux). They find that the magnetic twist and the size are not correlated, i.e. they are separate quantities describing the active region's characteristics, but both show a strong relation to the active region's CME productivity, suggesting that a combination of total flux and total twist would provide the best metric for forecasting. In addition, they find that the total free magnetic energy in an active region is more relevant to its CME productivity than is the overall twist (or its large-scale helicity) alone.  carried out a study where only a single metric was used to provide a probability forecast of flaring occurrence in an active region. The authors defined a quantity called the "effective connected magnetic field", B eff , that is calculated by connectivity matrices of field-line lengths and fluxes between positive and negative flux elements. They derived B eff for 298 active regions (93 X-and M-flaring, 205 non-flaring) from MDI/SOHO data over a 10 yr period during solar cycle 23, finding that B eff provides a robust measure to distinguish flaring from non-flaring regions. Leka and Barnes (2003a,b) conducted a thorough study into whether an active region will be flare-productive or flare-quiet. They derived various parameters, together with their statistical moments from vector magnetic field data. These parameters include the magnetic flux distribution, horizontal spatial gradients of the field, vertical current, twist parameter α, current helicity, shear angles, and photospheric excess magnetic energy. Leka and Barnes (2007) examined a set of 496 active regions observed in more than 1200 magnetograms. Statistical tests based on linear discriminant analysis are applied to the numerous photospheric magnetic parameters listed above, and compared to epochs of non-flaring. The most accurate predictor for the occurrence of large flares (defined as larger than M1), is found to be the total excess photospheric magnetic energy (in excess of the potential field). Other top performers identified are parameters that measure global properties of the field, i.e. ones that are integrated over the full active region. The analysis of Leka and Barnes (2007) shows that a number of parameters are important to quantify an active region as flare-productive, but many of these quantities are strongly correlated. For instance, large active regions as measured by the total magnetic flux also tend to have large vertical currents, significant excess energy, and significant current helicity. Thus, the choice of which few variables to use is not unique. The discriminant analysis in Leka and Barnes (2007) gives a success rate of 80% for the occurrence of flares of GOES class C1 or larger. This is compared to a success rate of 70% when all regions are uniformly classified as "flare quiet". For flares of class M1 and larger, the success rate obtained is 93%, which is only slightly larger than the rate of 91% that is obtained for the uniform flare-quiet labelling. This is an intrinsic problem related to the low occurrence rate of large flares, as low event occurrences typically lead to large false alarm rates. Based on these outcomes, the authors conclude that the state of the photospheric magnetic field at any given time does not have a large influence on whether an active region will be flare productive or not. This is in line with the conclusion in Schrijver (2009), that so far none of these various measures (either single or combined) has been identified as being particularly well correlated with the flare/CME activity of an active region.
In order to directly compare the relative success of existing flare prediction algorithms, an interagency "All-Clear workshop" was held in Colorado in 2009. The outcomes of the workshop are summarised in Barnes et al. (2016). The workshop brought together people working on a number of existing algorithms. Previously, it had been difficult to compare the outcome of different flare forecasting algorithms, as they use different analysis techniques and data sets. The analysis carried out for the All-Clear workshop was based on a common data set, using line-of-sight magnetic field and continuum intensity maps from MDI/SOHO, and applied standard verification metrics to evaluate the performance of the different methods. In total, 11 different algorithms were applied including the ones that are described above, i.e. the methods of , Falconer et al. (2002Falconer et al. ( , 2008, Barnes (2003a,b, 2007), Schrijver (2007) and Wheatland (2004). In addition, a number of further algorithms were included, which are now described. The method of Qahwaji (2008, 2009) is based on a feature-recognition technique to automatically derive the McIntosh classification from white-light images and subsequently uses machine learning to make a probability forecast based on the sunspot class. The method by Bloomfield et al. (2012) uses historical flaring rates related to the McIntosh class of the source active region to make forecasts using Poisson probabilities. Yuan et al. ( , 2011 apply machine learning techniques to provide flare forecasts from characteristic magnetic parameters, such as the total unsigned magnetic flux, length of the strong-gradient polarity inversion line, and total magnetic energy dissipation. The method by McAteer et al. (2010) is based on the fractal dimension of the magnetic flux concentrations in an active region to determine its flare productivity. The algorithm by Higgins et al. (2011) and  uses automatic characterization of various magnetic parameters of the active region to apply flare forecasts based on cascade correlation network techniques. The main outcome of the comprehensive workshop study detailed in Barnes et al. (2016) is that no method was clearly outperforming the others. This may be due to the fact that there are strong correlations between the different parameters that are used in the different prediction schemes. For the prediction of flares of class M and larger, the set of methods studied tends to give a slightly positive skill score, but none of the methods achieved large skill score values. All skill scores were 0.2. For context, a perfect forecast would have a skill score of 1 and forecasts that have positive skill score values indicate that the method is more successful than the reference method being used.
An interesting case study against which to examine and test ideas around predicting flares/CMEs based on observations is NOAA active region 12192. This region appeared on the solar disk in October 2014 and it hosted the largest sunspot seen by that point in solar cycle 24. From its size and complexity (βγ δ) it was expected to be the source of a significant amount of free magnetic energy and high-intensity flares, therefore catching the attention of space weather forecasters. Indeed the region produced six GOES X-class flares during its disk passage but no successful eruptions occurred (Thalmann et al. 2015;Sun et al. 2015). Although the magnetic field of NOAA region 12192 was complex there was a separation of the main positive and negative magnetic polarities with time. No convergence or collision of polarities was seen during its first disk passage. Small bipoles between the main spots (called "serpentine field") were observed indicating that the flux emergence process was still taking place. Overall, the active region configuration was that of one main bipole, with at least 2 new bipolar emergences on the periphery of the sunspot region or near the centre, later in time. The region was studied by Sun et al. (2015) who compared NOAA active region 12192 to two other flare-productive regions (11429 and 11158). They constructed NLFFF extrapolations and calculated the free energy in the coronal field, the decay index of the coronal field and the squashing factor. They find that the region's core field had a low level of non-potentiality and that the overlying field was relatively strong. The lack of colliding polarities (and therefore opportunities for magnetic reconnection) may be the reason that a strongly non-potential field didn't build up along any polarity inversion line in the active region and the lack of flux cancellation would have resulted in no tether cutting (which removes overlying restraining field).
The suggestion that CMEs may act as a valve for magnetic helicity that has accumulated in the corona appears to lead intuitively to the idea that CMEs might occur once a threshold in this quantity has been reached. Can this be used to forecast the occurrence of a CME? Although this idea is by no means accepted or proven, some work has indicated that magnetic helicity is a necessary but not sufficient condition for an eruption (Amari et al. 2003b) whilst others show that magnetic helicity is not necessary for CMEs ). It is worth noting that CMEs originate from a wide range of coronal magnetic structures, from bright points with a flux (φ) of around 10 20 Mx, to large active regions, where φ = 10 22 Mx, giving a potentially large range in magnetic helicity values (since helicity scales with φ 2 ). Therefore, there will not be a one-size-fits-all value for magnetic helicity in order to create an eruptive configuration. Instead, the helicity of an eruptive configuration will vary from case-to-case. A study by Pariat et al. (2017) takes a different approach into the potential utility of magnetic helicity to forecast and eruption. Pariat et al. (2017) study the evolution of stable and eruptive magnetic configurations using the flux emergence simulations of Leake et al. (2013. The relative magnetic helicity of the computed field is decomposed into that of the current carrying component and that of the potential component. They find that the eruptive configurations are not necessarily those with the highest magnetic helicity value but that an eruptivity indicator exists in the ratio of these two components of the relative magnetic helicity. Opening, for the first time, an avenue for investigating how magnetic helicity may be used for CME forecasting.

Predictions Based on MHD Models
MHD models provide a potentially powerful future tool in the quest for a predictive capability of CMEs. However, at present, the MHD numerical models described in Sect. 2.2 cannot be used for predicting the time of an eruption. This is in part because most of these models do not provide quantitative onset criteria that could then be applied to observations or magneticfield extrapolations of potential CME source regions. Even those that do (i.e. the kink and torus instability models) can at present only provide rough estimates of the onset criteria. This is because the exact thresholds of these instabilities depends on the detailed structure of the magnetic-field configuration, which is not sufficiently well known. Moreover, exact thresholds have so far only been derived for rather idealized model configurations rather than those exhibited by the real Sun. A complementary approach is to create a simulation that is driven by line-of-sight magnetograms. Using real data provides an opportunity to create a realistic simulation so that the coronal field can be monitored for its eruptive potential. This is the approach used by Gibb et al. (2014) in a study of NOAA active region 10977 that employs a magnetofrictional method driven by line-of-sight magnetograms. The study found that a realistic magnetic field line reconstruction was not possible around the time that the active region erupted suggesting that the approach captured well the destabilisation of the coronal field. Since the top boundary of the simulation box is closed, the flux rope that had formed and started to erupt wasn't able to escape the computational volume. That is, the computed field was trying to relax by ejecting the flux rope but the boundary conditions did not allow this to happen. Consequently, this led to the formation of highly twisted structures in the computational volume.
Apart from the eruption onset, one would also like to forecast the arrival time and impact at the Earth of the CME, once the eruption is on its way. CMEs typical take from 0.7 to 4 days to travel from the Sun to the Earth. This is an amount of time that allows an MHD model, based on solar data, to be run. In contrast, predicting the effects of energetic-particle events created by CME-induced shocks is much more challenging, as these particles can reach the Earth in less than an hour. We therefore consider here only to what extent MHD models can be employed for predicting the time and impact of space-weather disturbances caused by CMEs. The most important quantities for the latter are the speed, duration, plasma density, and magnetic field strength and orientation of the CME. Sophisticated models should also include the effects of the shock and the sheath that precede the CME flux-rope in fast events.
As discussed in Jin et al. (2017), present CME forecasting models can be divided into different categories (see also Siscoe and Schwenn 2006;Messerotti et al. 2009). The first categories contain empirical forecasting models (e.g., Gopalswamy et al. 2001;Schwenn et al. 2005), kinematic models such as the cone model (e.g., Zhao et al. 2002) and analytical models such as the drag-model ). These models have been widely used to predict the arrival time, orientation, and speed of CMEs (and of associated shocks), but they cannot provide predictions of the plasma density and the magnetic field (though see, for example,  for attempts to predict the maximum magnetic field strength in an CME from its predicted speed). To this end, the use of models of the third category, namely of MHD simulations, appears inevitable.
However, MHD simulations that employ observed data and continuously model an eruption from its onset in the low corona out to 1 AU (see Manchester et al. in this issue) are limited to using data from the Earth's point-of-view. Therefore, the Earth-directed component of the CME velocity is not accurately determined. Moreover, these simulations cannot yet be performed in a constantly evolving and accurate solar wind. For this, in situ solar wind measurements are needed that provide real data to adjust the models to. This limits the capability of these MHD models, but efforts are being made to develop such models for future operational use (e.g., Jin et al. 2017). Presently, ENLIL (Odstrcil 2003) is the only MHD model that is used for operational forecasting (at NOAA and the UK's Met Office Space Weather Operations Centre for example). ENLIL produces estimates of the solar wind speed, density, and temperature and magnetic field, and provides 1-4 day advance warning of disturbances that will produce geomagnetic storms, such as CMEs. To be applicable in practice, simplifications had to be made to the model. For example, ENLIL simulates the solar wind plasma and magnetic field only beyond 30 R , i.e., the corona is neglected. To feed the model, the solar wind plasma flow and radial magnetic field at 30 R are provided by the semi-empirical WSA model (Wang and Sheeley 1990;Arge and Pizzo 2000). If an eruption is to be simulated, the cone model is used to set up the location, direction, and speed of the CME (e.g., Odstrcil et al. 2005). Also, while the solar wind magnetic field is included in the model, CMEs are initiated without an internal magnetic field. Therefore, the model cannot yet be used to predict the southward component of the magnetic field during CMEs. However, observations of the CME source region can be used to give first order approximation of the magnetic configuration of the erupting structure. For example, through the detection of a flux rope and determination of the orientation of its axis and its chiral properties (e.g. ).

Bringing Observational and Theoretical Predictions Closer Together
There is much to be learnt from the complementary approaches of observational CME studies and theoretical work, including numerical simulations. In terms of the most promising approaches to predicting CME occurrence ahead of time, combining current observational and theoretical knowledge, there appear to be three main avenues: (i) Flux imbalance between a flux rope and overlying arcade (ii) Torus instability due to the rate of decay of the arcade overlying a flux rope (iii) Helicity proxy derived from the current carrying and potential field components For (i), observations of flux cancellation can be used to investigate the flux built into the rope versus that in the overlying field. However, more work is needed to understand how the shear and spatial extent of the cancellation in the active region influence what proportion of cancelled flux is built into the rope, and what the flux imbalance criteria is. For (ii), the gradient of the overlying field can be investigated using a potential field extrapolation, but simulations so far describe rather idealised configurations and the exact stability criteria varies on a case-to-case basis. For (iii), the theoretical knowledge is still developing, but early indications suggest that there might be a critical value for the ratio of helicity in the current carrying field component to the helicity in potential component. Still, none of the above are able to give the desired several days notice of the occurrence of a CME. In addition, approaches (i) and (ii) are relevant for the flux rope concept of CMEs rather than a sheared arcade.

Summary
This review aims to provide a summary of the wide range of observational and theoretical work on coronal mass ejections (CMEs). These eruptions are known to be a driver of space weather, including the most severe space weather events. In response to this, an increasing number of developed and developing nations provide space weather forecasts to help mitigate the risk of economic and societal disruption resulting from the technology affected by space weather. The severe impact that results from the arrival of a CME of high velocity, southward directed magnetic field and enhanced plasma density means that there is a significant interest in forecasting these eruptions ahead of time.
Currently the physical processes at play after CME onset can be brought together in the standard CHSKP model. This captures and explains CME observations well, once the eruption is underway. However, the energy required to power a CME (and any associated flare) cannot be supplied to the corona on the timescale of the dynamic event and instead it must be built up and stored in the coronal magnetic field in the hours or days beforehand. We need to extend the standard model further back in time to understand the physical processes that create an eruptive magnetic field configuration and the associated currents, as well as the processes that trigger and drive the eruption. CME occurrence requires more than just estimating the free magnetic energy stored in the coronal field-the free energy is a necessary but not sufficient condition. The exact magnetic field configuration and its stability need to be better characterised.
Many flare/CME forecasts currently focus on using a snapshot of the configuration of the photosheric magnetic field to provide a probabilistic forecast. Improvements to this technique are limited by the physical origins of a CME being related to the stability of the coronal field, when a trigger switches a driver on, which is unlikely to be captured in the photospheric data. CME triggers themselves can be monitored over many days using both photospheric and coronal data. But this must be supplemented with data to reveal the likelihood of a CME occurring, using either a NLFFF extrapolation or an eruptivity signature(s) determined using atmospheric observations. With no regular direct measurements of the coronal magnetic field, observations are limited to inferences obtained from plasma structures.
To move on from probabilistic forecasting, to the capability of being able to predict CME onset time, we first need a clear understanding of the criteria that encapsulates the eruptive status of the magnetic field. At what point does a CME driver switch on? Parametric modelling studies that can identify the loss of stability criteria are required and this could be determined from more realistic magnetic field configurations and their evolution. Developing the capability to conduct data driven simulations, relevant for active regions from their birth to their dispersal and decay into the background, field will greatly help this aim.