Origins of Rolling, Twisting and Non-Radial Propagation of Eruptive Solar Events
We demonstrate that major asymmetries in erupting filaments and CMEs, namely major twists and non-radial motions are typically related to the larger-scale ambient environment around eruptive events. Our analysis of prominence eruptions observed by the STEREO, SDO and SOHO spacecraft shows that prominence spines retain, during the initial phases, the thin ribbon-like topology they had prior to the eruption. This topology allows bending, rolling, and twisting during the early phase of the eruption, but not before. The combined ascent and initial bending of the filament ribbon is non-radial in the same general direction as for the enveloping CME. However, the non-radial motion of the filament is greater than that of the CME. In considering the global magnetic environment around CMEs, as approximated by the Potential Field Source Surface (PFSS) model, we find that the non-radial propagation of both erupting filaments and associated CMEs is correlated with the presence of nearby coronal holes, which deflect the erupting plasma and embedded fields. In addition, CME and filament motions respectively are guided towards weaker field regions, namely null points existing at different heights in the overlying configuration. Due to the presence of the coronal hole, the large-scale forces acting on the CME may be asymmetric. We find that the CME propagates usually non-radially in the direction of least resistance, which is always away from the coronal hole. We demonstrate these results using both low and high latitude examples.
Eruptive solar events which include a coronal mass ejection (CME), an erupting filament with surrounding cavity and a flare are distinguished from simple confined flares which lack most of these erupting signatures. The term CME is increasingly used in reference to the whole eruptive sequence and not just the ejection of coronal mass. Regardless of definition, there has been widespread recognition, for over three decades, of the close association among CMEs, erupting filaments, coronal cavities and flares.
What do CMEs and erupting filaments have in common? All of these events originate above and around regions where the magnetic field changes polarity, seen in photospheric magnetograms and variously named a neutral line or polarity inversion line or a polarity reversal boundary. We know that erupting filaments, in particular, originate from specific polarity reversal boundaries that are also inside filament channels (Martin, 1990; Gaizauskas, 1998). We also already know they only come from filament channels that have reached maximum development; i.e. have reached maximum magnetic shear along the polarity boundary (Martin, 1998). Therefore, not every polarity reversal boundary is a filament channel.
A filament channel is a volume of space around (and encompassing) a polarity reversal boundary with maximum magnetic shear. A filament is not necessarily present but the channel is a necessary condition for the formation of a filament. The formation of filament channels in the chromosphere has been described by Smith (1968), Foukal (1971), Martin (1990), Gaizauskas et al. (1997), Wang and Muglach (2007) and Martin, Lin, and Engvold (2008) and Martin et al. (2008). In and around active regions or decaying active regions observed in H a filament channel may be recognized by the presence of fibrils aligned along a polarity reversal boundary and the local magnetic field. Because no fibrils cross the polarity reversal boundary in a fully-developed filament channel and fibrils are field-aligned, one may conclude that no magnetic field lines from active region or network magnetic fields cross this boundary in the chromosphere, above or within filaments. This is also true for the low corona. The initial lack of coronal connectivity across polarity reversal boundary inside the filament channel volume, under, across and along the filament is confirmed by Wood and Martens (2003). We find that the coronal cells observed and modeled by Sheeley and Warren (2012), do not cross the polarity reversal boundary within a filament channel at the heights below the filament spine top (Figure 1). Coronal cells originate from the network field concentrations and show the same patterns as chromospheric fibrils because they follow the same filament channel magnetic topology, with a (presumably) strong horizontal component of the field.
As deduced from chromospheric fibrils and coronal cells, the direction of the network magnetic fields associated with chromospheric fibrils and coronal cells changes from upward on one side of a polarity reversal boundary to horizontal along the polarity reversal boundary and to vertically downward on the other side of the polarity boundary. A vital precondition for the formation of a filament channel is a magnetic field with a strong horizontal component (Gaizauskas, 1998) and the chromospheric fibrils in a filament channel trace such a horizontal pattern along its polarity reversal boundary. So far no observational evidence of the existence of an ascending twisted flux rope at the chromospheric level within a filament channel above a polarity reversal boundary has been found; this would require fibrils to systematically change their orientation and cross the polarity reversal boundary. Conversely, were the reverse to be true, emerging, twisted flux tubes playing a major role in the physics of filaments and filament channels, one would expect there to be ubiquitous evidence of polarity crossing fibrils, which would necessarily appear earlier, if there was a flux rope axis submerged in the photosphere (Wood and Martens, 2003).
This is a fundamental observation which allows one to deduce at least one significant moment in the formation of CMEs: the most likely time of initiation of twist in CMEs. There is evidence that once an eruption is underway, a twisted flux rope does indeed form during magnetic reconnection of the stretched, sheared overlying coronal loop system and concurrently with the initiation of flare ribbons in the chromosphere (Wang, Muglach and Kliem, 2009). In other words, the observational evidence suggests that the formation of the flux rope of CMEs occurs after the filament and coronal loop system are already ascending.
Similar to the overlying coronal loops within a coronal loop system, before eruption, filaments rise slowly. Unlike coronal loops and if a filament erupts non-radially, as frequently happens, the top of its spine first bends to one side and evolves into a sideways rolling motion, known as the roll effect. As shown by 304Å observations from SOHO and STEREO and earlier H Doppler observations, the rolling motion unambiguously propagates down the legs of erupting filaments resulting in the large scale twists which are commonly observed in them. The rolling of the spine initiates structural twist with opposite chirality in the two legs (Figure 2). Concurrently, rotational mass motions with the opposite signs of chirality (and helicity) in the two legs are observed to flow downward while the body of the filament ascends. However, in asymmetric eruptions, both legs are not always observed (Martin, 2003; Panasenco and Martin, 2008).
In addition to the observed absence of twist in the pre-eruptive state, further evidence that the energy creating the twist comes from above was found in Doppler shifts (Martin 2003); the rotational motions in the legs of erupting filaments are not only opposite in sign to each other but the twists in both legs are opposite in sign to that required if the observed sense of twist were generated at the feet or in the legs of the erupting filament (Martin 2003, Panasenco and Martin, 2008).
Here we will examine events showing that the combined ascent and initial bending of the filament is non-radial in the same general direction as for the surrounding CME. Also, the non-radial motion of the filament is greater than that of the CME. By considering the global magnetic environment around CMEs, as can be approximated by the Potential Field Source Surface (PFSS) model, we will see that both erupting filaments and their surrounding CMEs are non-radial moving away from nearby coronal holes and toward null points. Due to the presence of the coronal hole, both the local forces on an erupting filament and the global forces on the CME are asymmetric. That the CME propagates non-radially always away from coronal holes was demonstrated in many studies (e.g., Plunkett et al., 2001; Gopalswamy et al., 2003; Cremades and Bothmer, 2004; Cremades, Bothmer and Tripathi, 2006; Gopalswamy et al. 2009; Kilpua et al., 2009; Zuccarello et al., 2012). Prominence deflection and rolling motion during eruptions received less attention, though there has been some work in the same general direction as advocated here (Filippov, Gopalswamy and Lozhechkin, 2001; Martin 2003; Panasenco and Martin, 2008; Bemporad 2009; Panasenco et al., 2011; Pevtsov, Panasenco and Martin, 2012; Liewer, Panasenco and Hall, 2012). The formation of the CME flux rope appears to occur in the early phase of filament eruption, when the rolling motion of the filament is already in progress (Liewer, Panasenco and Hall, 2012).
The roll effect appears very early during the beginning phase of prominence eruptions. We will show that it is caused by an asymmetry in magnetic structure and magnetic flux density on the two sides of the filament. The non-radial direction of the eruption of the whole filament system is caused directly either by the open coronal hole magnetic field near the filament channel (cases studied in this paper), or by other strong magnetic field configurations which might be in the neighborhood of the eruption (strong in reference to the overall energy of the eruption).
The paper is organized as follows: we begin with a general introduction to the coronal structures which lead to eruptive solar events. We then discuss a series of recent eruptions within such structures, which may be referred to generically as either streamers or pseudostreamers, exhibiting non-radial expansion and propagation. We then correlate such behavior with the properties of the ambient corona. We also discuss the three-dimensional structure of such events, focusing on one particular case where the two well separated STEREO viewpoints show that what appears to be a highly kinked, flux rope-like erupting prominence and cavity is a projection effect.
2 Magnetic topologies overlying filament channels
Filament channels and filaments develop along and above polarity reversal boundaries in the photosphere and lie under coronal loop systems also known as coronal arcades. In the instances most commonly studied, these arcades correspond to a helmet-streamer type structures in the outer corona. The outer, weaker fields do not always close over the arcade but may remain open, oppositely directed adjacent fields (what we will refer to as a globally dipolar structure). Recently, filament channels and filaments lying below arcades and belonging to more complex magnetic structures have been shown to be very important for multiple and global solar eruptive events, such as that of 01 August 2010 (Schrijver and Title, 2011; Török et al., 2011). The prominences in that event had a magnetic structure corresponding to a hierarchy of higher order multipolar structures resulting in a pseudostreamer (Wang, Sheeley and Rich, 2007), which was found to contain twin filaments at its base (Panasenco and Velli, 2010). Such twin filaments are topologically connected, sharing a neutral point and a separatrix dome. This was a case in which two polarity reversal boundaries contain between them fields with a polarity opposite to that of the global unipolar configuration surrounding them (tripolar pseudostreamer).
The magnetic field, based on the potential field source surface (PFSS) reconstruction, shows that field lines of the arcades above filaments are often asymmetrically placed with respect to the polarity reversal boundaries, because the overarching magnetic field is stronger at one side. In the PFSS this is shown by the greater density of field lines on one side of the arcade with respect to the other. As we shall see, this is a tell-tale sign for strong rolling and nonradiality of the eruption. The question naturally arises as to the fidelity of the PFSS model in the neighborhood of filaments and filament channels, where the strong shears naturally lead to the formations of currents in the corona, currents which are also required to get to the instability involved in the filament eruption and CME formation. In reality, the PFSS does a good job because at least a component of the axial field along the filament spine is potential in nature. This may be understood as having to do with the finite longitudinal extent of the channel Ð so that the polarities on opposite side of the Êinversion line are not translationally invariant along it. Indeed, this component is recovered also by photospheric shear and flux transport models which correctly describe filament channel formation statistically based on photospheric transport (Yeates, Mackay and van Ballegooijen, 2007). In the following sections we will make use of the PFSS model of Schrijver and DeRosa (2003) to study the coronal magnetic field and its role in determining the trajectory of filaments and CMEs for a whole series of events including even higher order multipoles below streamers and pseudostreamers.
3 Non-radial eruptions with the roll effect
We consider here a set of prominence and CME eruptions which display both the roll effect and significant non-radiality. The different examples paint a picture of the rich variability in motions which the prominence and CME may display, depending on the detailed magnetic field configuration surrounding the eruption and the global coronal context. We use images from the Sun-Earth Connection Coronal and Heliospheric Investigation (SECCHI: Howard et al., 2008) aboard the Solar-Terrestrial Relations Observatory (STEREO: Kaiser et al., 2008) and from the Atmospheric Imaging Assembly (AIA: Lemen et al., 2012) aboard the Solar dynamics Observatory (SDO: Pesnell et al., 2012). We also use information from the Large Angle Spectrometric Coronagraph (LASCO; Brueckner et al 1995) on the Solar and Heliospheric Observatory (SOHO: Domingo et al., 1995).
3.1 2 November 2008 filament eruption
This strongly non-radial eruption (Figure 3) was described in detail by Kilpua et al. (2009). The right panel in Figure 3 shows that the filament channel was located under the northernmost arcade of a double arcade system underneath a pseudostreamer, a large-scale coronal structure which differs from the classical helmet streamer because the magnetic field is unipolar above the dome like closed field regions. The pseudostreamer itself, adjacent to a strong northern coronal hole, expanded outward with considerable equatorward deflection. At the early stage of the eruption, the motion of the filament was directed towards the local null point, O, below the pseudostreamer separatrix (Figure 3).
3.2 12 December 2008 filament eruption
This event, illustrated in Figure 4, is described in detail in the paper by Panasenco et al., 2011. It was observed by both STEREO-A, at the north east limb (middle panel), and STEREO-B, at the north west limb (left panel). Though the spacecraft were separated at the time by , both imagers show the rolling motion of the erupting prominence. The erupting prominence and CME were deviated from a radial trajectories by different degrees. To illustrate this, in Figure 5 we superimpose outlines of the outer boundaries of the CME and the bright core of the prominence plasma for different consecutive moments during the eruption (as observed by the Extreme Ultraviolet Imager (EUVI-A) 304 Å and COR1-A, respectively). The deviation of the CME from radial propagation is about , for the prominence Ð about and the difference between the prominence and CME central direction of the eruption is about for this STEREO-A view. Since the prominence erupts approximately toward STEREO-B, and the separation angle between A and B is close to , we may safely assume that the projection on the sky in the STEREO-A view shows a separation of trajectories which is smaller than the real separation between the propagation directions for prominence and CME, shown in Figure 5.
3.3 26 September 2009 filament eruption
This eruption was observed by both STEREO, at the north east limb by A, and across the disk by B as shown in Figure 7 . The separation angle between the spacecraft was . EUVI-B 304 Å images show the prominence crossing the disk during the eruption. STEREO-A observed this eruption above the limb. The position of the filament relative to the coronal arcade and open field of the coronal hole before the eruption is shown in Figure 6. The asymmetric location of the filament is very pronounced. Figure 7 shows the eruption on 26 September 2009 at 22:45 UT as observed simultaneously by STEREO-A/COR1 and EUVI instruments. The non-radial motion of the erupting filament is much greater than that of the CME. The CME envelope expands, containing the prominence within it, but the two bodies are not rigidly linked to one another. Plasma of the prominence moves approximately along the line which can be described as a tangent line to the limb from where the prominence erupted. Apparently the imbalance in magnetic pressure under the coronal arcade and around the filament during the early stage of the eruption was such that the stronger field on the coronal hole side of the arcade caused such an extensive lateral motion of the filament. The corresponding CME also deviates from radial propagation. Analyzing the projection on the sky as viewed by STEREO-A we found that the difference between the prominence and CME deviations amounted to about . Figure 8 shows an intermediate projection of the eruption on the plane of the sky as observed by the SOHO/EIT and LASCO C2 instruments. This projection viewpoint is midway between those of STEREO-A and B, therefore the apparent difference in the deviation angle for the prominence and CME is smaller, about .
3.4 30 April 2010 filament eruption.
This non-radial filament eruption with the roll effect was observed against the solar disk by the Mauna Loa Solar Observatory (MLSO) on 30 April 2010. The movie of the first 40 min of the eruption (23:25 - 00:05 UT) in the He I spectral line shows that the filament crossed the disk and the neighboring active region in the south-eastern direction. This direction of the eruption was away from the coronal hole north-west of the filament channel (movie available in the supplementary material). The distance which the filament covered over these 40 min was about 160000 km, corresponding to a projected speed of 65 km s. The first four frames in Figure 9 correspond to this period but observed from the different point of view of STEREO-A . The top of the prominence is bending away and to the left of the viewing direction of STEREO-A (Figure 9) creating an apparent kink and twist in the prominence ribbon. However, the apparent kink crossing may be an illusion from viewing the warped ribbon of erupting filament mass seen in a 2-dimensional image. MLSO and STEREO-B observations show the filament spine as straight, without any twisting, crossing the disk. The apparent kink disappears as seen in the final six frames in Figure 9. The images are suggestive of a kink structure until 00:26 UT. With the next 304 Å image (at 00:36 UT) the apparent kink disappears as the eruption proceeds along the southeast extension of the filament channel (final six frames in Figure 9).
Török, Berger and Kliem (2010) have pointed out examples of kink unstable flux ropes where there is no crossing as seen from angles. While this may be the case, additional considerations may come from the properties of the line-tied kink instability. The kink instability requires, for a flux rope, a certain minimum number of turns that the magnetic field must make from end to end, not on the axis of the flux rope configuration, but in its immediate surroundings. Oftentimes this is described in terms of the ratio of the average poloidal flux (i.e. the flux circulating around the axis of the flux rope) to the axial flux (the integrated field along the flux rope axis or spine). The required twist for kink instability is at least 2.5 and depending on the configurations considerably more. This is more than one complete turn of the field in the domain (Hood and Priest, 1979; Velli, Hood and Einaudi, 1990) and is a prerequisite, independent of observed kinking of the spine. If the spine is bent, and there is no tracer of plasma making one complete turn around the spine, the presumed axis of the prominence, the evidence for kinking becomes weak. The introduction of toroidal geometry (and the hoop force) does not change these constraints on the kink instability, i.e the flux ropes of Török, Berger and Kliem (2010) also have field lines wrapping around the axis for more than one complete turn. Because individual field lines are the only real relevant concept once line tying is taken into account (because the concept of a flux surface is lost if lines end in different regions in the photosphere), rather than flux surfaces, plasma tracers along such turning field lines should never disappear no matter what the point of view.
The position of the filament under its coronal arcade before its eruption was very asymmetric. Figure 10 shows the PFSS extrapolation of the magnetic field lines above the filament on 2010 Apr 30 18:04 UT. This asymmetry is thought to lead to asymmetric forces on the filament which cause the rolling motion of the filament at the early stages of the eruption, and create a condition for the non-radial propagation of the corresponding CME (Figure 11). The right image in Figure 11 shows the CME front and position of the prominence at the limb from STEREO-A; the CME non-radial propagation relatively to its point of origin is clearly visible.
4 Radial eruptions with the roll effect
Here we consider cases where there is still significant rolling of the prominence spine before the CME is formed, but subsequently the filament and CME propagate outwards radially. Such cases are usually observed when the pre-eruptive filament is lying at the base of a pseudostreamer magnetic configuration where the arcades below the pseudostreamer have a much greater lateral extension relative to the pseudostreamer stalk which is essentially radial (an example of a non-radial pseudostreamer is shown in Figure 3). The asymmetry in the filament position relative to the overlying coronal arcade might lead to additional force on the filament that could enhance the rolling motion of the erupting filament during the early stage of the eruption. This lateral motion of the filament is directed towards the null point in the magnetic configuration. The radial geometry of the pseudostreamer branches (defined as bundles of open magnetic field lines on opposite sides of the spine-fan projection on the sky) in this case favors the radial propagation for the CME flux rope formed in the later stages of the filament eruption.
4.1 1 August 2010 twin filaments sympathetic eruptions.
Pseudostreamers appear in unipolar regions above multiple polarity reversal boundaries. Some of these polarity reversal boundaries can be filament channels, and when this is the case they can occur as twin filament channels often containing twin filaments. Figure 12 shows twin filaments and the coronal magnetic field configuration above them on 2010 August 01. This event has been described in detail by Panasenco and Velli, 2010; Schrijver and Title, 2011; Török et al., 2011. The slight rolling motion of the erupting filaments toward the null point in the pseudostreamer configuration has been observed and reported in Török et al. (2011). Overall both CMEs, that originated under the pseudostreamer, propagated radially along the axis between the pseudostreamer branches.
4.2 16 May 2007 south polar crown filament eruption.
Figure 13 shows a PFSS extrapolation of the open field lines of the pseudostreamer and coronal loops overlying a filament under the southern lobe of the pseudostreamer on 15 May 2007. On the right column are anaglyph images from EUVI showing an erupting filament on 16 May 2007; the separation angle between STEREO-A and B is 8.14 degrees. To assist the reader, the upper edge of the filament is outlined by a red line and the lower edge by a light blue line. When one watches the 3D movies of this event created from STEREO/EUVI data, one can identify a rolling motion beginning to the east of the middle of the filament. A consequence of this rolling motion is that the bottom edge of the filament crosses gradually over the top edge. The roll appears to begin between the white Xs superposed on the blue line. It propagates in both directions away from the initial site (white Xs) and towards the apparent ends of the prominence as seen against the chromosphere. We also see that this erupting filament is erupting non-radially and rolling northward in the direction of the pseudostreamer null point (left panel in Figure 13).
5 Radial eruptions without the roll effect
The final examples are cases where there is no visible rolling of the filament, and both the filament and CME propagate radially outwards through the corona.
5.1 3D Reconstruction of the erupting prominence on 28 February 2010
A polar crown prominence around the northern coronal hole (CH) erupted over the course of a day on 28 February 2010. The eruption was located on the far-side of the solar disk but it was captured by the EUVI imagers onboard both STEREO satellites. The two EUV viewpoints revealed a very different behavior of the erupting prominence during its ascent which is the reason we chose to analyse its 3D morphology in detail. In particular, EUVI-A shows a long piece of the polar crown prominence lifting off parallel to the surface. In contrast, the EUVI-B images show a very strongly kinked prominence (Figure 14).
It appears, therefore, that the apparent kinked morphology from one viewpoint could be just a projection effect. Since we have observations from two viewpoints, we can quantify this effect because we can localize the prominence in 3D space. Tiepointing is the standard technique for this purpose and has been described before (Liewer et al., 2009; Liewer et al., 2011; Thompson, Kliem, and Török, 2012). The technique relies of the correct identification of the same structure in the two images. Filaments are ideally suited for this because of their highly detailed morphology. Because we are interested in the 3D shape of the large scale structure of the prominence, we concentrate only on its ’backbone’ in the EUVI 304 Å images. We use relatively few points to determine its orientation. We select six points for tiepointing which are shown by the yellow crosses in Figure 14. To aid the reader, we connect some of those points with the dashed lines. We use the IDL routine sunloop which is part of the SECCHI software distribution in Solarsoft. The routine allows the user to select one feature in one viewpoint and then click on the same point (according to the use) on the other image. The program returns the location of the feature in 3D space defined by its longitude, latitude and radial distance from the solar surface. Our results are shown in Figure 15 along with the calculated projections of the prominence backbone for the EUVI-A, and B views and on the surface.
The measurements verify that the appearance of a kink crossing is simply a projection effect. It is probably caused by the narrow shape of the prominence. The prominence is quite thin with a 70 deg longitudinal extent but only 8 deg of latitudinal extent. We cannot make a definitive measurement of the size of the prominence because its footpoints in the lower atmosphere are not visible from both EUVI telescopes. But we can easily trace the highest point along the backbone during the eruption. The resulting 3D height-time plots is shown in Figure 16 (left) where we can see that the prominence is erupting quite radially; the latitude changes by only 4 deg while the point remains within a 10 deg longitudinal range. The height-time profile (Figure 16, (right)) is quite typical for such slow eruptions. We measure a very slow rising during the first seven hours followed by a rapid increase after 10 UT. We fit the height-time plot only for the period of the fast rise using a 2nd degree polynomial. The resulting velocity profile is marked by the straight line of star symbols. We find that the prominence attains a speed of only 50 km s at 14 UT when it leaves the EUVI field of view. The prominence and the corresponding CME propagate radially outwards through the corona (Figure 17).
In summary, the three-dimensional analysis of the 28 February 2010 prominence results in two important findings for this paper. First, we find that the kink morphology can be a result of projection effect which has important ramifications for past identification of kinked prominences. Second, we determine that the eruption is radial without any hint of equatorial deflection even though it is a slow prominent eruption (Figure 17). This is contrary to the prevailing consensus in the field which holds that prominence eruptions, especially during solar minimum, tend to be deflected towards the magnetic equator. Although this is certainly true for some events, it should not be taken for granted for any event.
5.2 Radial eruption on 6 July 2010.
A polar crown prominence around the northern CH erupted on 6 July 2010. The eruption was captured by the EUVI imagers onboard both STEREO satellites and SDO. The erupting prominence appeared to rise very radially and showed very little bending of the prominence spine during this eruption (Figure 18). In Figure 19 we also examined the position of the filament before its eruption on 6 July (SDO/AIA) beneath the coronal fields generated by the PFSS model. The filament is centered symmetrically beneath the arcade of overlying coronal loops. Therefore, as the filament rises, the magnetic pressure of the magnetized filament cavity and overlying arcade should be approximately equal on the two sides of the filament. The absence of the lateral rolling motion of the filament spine is consistent with the deduced presence of the magnetic force balance between the fields on the two sides of the prominence.
The path followed by the CME, and by the erupting filament immersed within it, depends on the macroscopic forces to which the erupting plasma is subject. Both the roll effect and the variable CME trajectories imply that the energy in the background corona is non-negligible compared to that liberated in the ejection. In the opposite case of negligible background energy, one would expect an essentially radial expulsion with a self-similar expanding spherical shock wave à la Sedov (Sedov, 1981). Describing the propagation of a part of a plasma cloud involves considerable difficulties, from the proper definition of the plasma cloud volume, to the separation of plasma cloud from the background field, to understanding the photospheric inertial effects of the parts of the plasma volume which remains connected to the Sun. Indeed, some or all of these problems are also intrinsic to numerical simulations. Consider then a volume of plasma as it is ejected and propagates outward into the solar corona. Parker (1957) focused on the possible shapes the plasma should take depending on the relative importance of kinetic and magnetic pressures, and then included the effects of the outward motion of the center of mass in a field with gradients. Assuming a general shape in the form of a deformable prolate spheroid, and using the virial theorem in its tensor form, he showed that if the internal pressure was dominated by the gas, rather than the magnetic field, the plasmoid would elongate into an infinitely thin shape aligned with the external magnetic field. Pneuman and Cargill (1985) showed that if the plasma (ratio of kinetic to magnetic pressures) inside the body was very small, and it was initially in equilibrium with its surroundings, then it would move outward self-similarly, maintaining its shape. If we assume , the dominant force in determining the motion of the plasma center of mass will be the one arising from the magnetic field stresses, which, in the simplest case of a magnetically isolated plasma body, may be written as
where is the surface bounding the volume of the ejected material and is the external magnetic field, as perturbed by the presence of the ejecting plasma itself. More generally the global force has a contribution coming from the non-diagonal part of the magnetic stress (magnetic tension) whose role is to eliminate the net force parallel to the magnetic field connecting the ejecta to the (arbitrarily defined at this point) outer magnetic field,
Parker (1957), Pneuman and Cargill (1985). For the cases of simple radial magnetic field background, these expressions lead to the so-called melon seed expulsion of diamagnetic plasmoids (see also Rappazzo et al. 2005 for an application to blobs formed at the tips of helmet streamers).
These expressions are important because they show that regions in the corona separating strong and weak magnetic field regions, for example the boundaries of coronal holes, will act as walls repelling the motion of plasmoids from the strong field regions, whereas the neighborhoods of coronal null points (for example, between spine and fan for a pseudostreamer configuration), will be regions towards which an ejecting plasma will be accelerated, as they are regions of vanishingly small magnetic pressure. The dynamics of course is richer, and nonlinearly so, because as a consequence of such motions regions neighboring null points may develop strong current sheets, leading to attraction along a privileged axis (locally orthogonal to the tearing sheet) and repulsive in the plane of the sheet itself (a representative example is given by the sympathetic eruptions simulated in Török et al. 2011, where the breakout of one filament at first inhibits that of the other, but once sufficient flux is removed, the breakout of the second filament also follows). A detailed analysis of the acceleration of the filament and global CME requires separating the components of an eruption. Consider then the global CME as a large-scale plasmoid bubble, dominated by the cavity magnetic field, elongated along the filament axis. In this case one may apply the simplifications of a magnetically isolated system, compared to the rest of the coronal field, which in the case of a CME with small dimensions compared to the global coronal scales would lead to an expression for the force as something like
with the volume of the CME and a constant . This can explain departures from radial expansion due, for example, to the presence of a coronal hole, which occurs on large scales. For the filament itself, local deflections occuring in the neighborhood of the relatively low-lying null points in pseudo-streamer or more general multipolar configurations, and complexities due to field lines not isolating the plasma entirely need to be taken into account. Differences in the acceleration of the filament (or filament spine) and the overall CME could therefore arise quite naturally in the initial phases of the eruption. More generally the expansion of the filament leads not only to macroscopic forces but also moments tending to rotate the prominence spine (responsible for the roll effect). Here the natural response of the filament should lead to deflection and alignment following the neighboring cavity field. As shown in the PFSS extrapolations, the asymmetries in the overlying arcade field correspond to both roll and deflection (Figure 20).
When multiple arcades are present below a pseudo-streamer or a streamer configuration, the various polarity regions form a hierarchy of null points (and possible lines/fan structures). From the PFSS modeling we found that the lateral rolling motion of the filament is directed first towards the closest, local null point above the (asymmetric) coronal arcade before this arcade reconnects creating the CME flux rope. Once reconnection has occurred, the just-born CME flux rope then follows moving non-radially towards the higher, more distant global null points/lines separating larger scale structures (the fan/spine null point of the large-scale pseudo-streamer, or the tip of the helmet streamer for a globally bipolar configuration).
We have considered the global environment around CMEs, as can be approximated by the Potential Field Source Surface (PFSS) model. From our analyses, we have demonstrated the following relationships for eruptive solar events consisting of an erupting filament surrounded by a cavity and an enveloping CME:
Both erupting filaments and their surrounding CMEs are non-radial only in the direction away from a nearby coronal hole and toward local and global null points. Due to the presence of the coronal hole, the global forces on the CME are asymmetric. The CME propagates non-radially in the direction of least resistance as we demonstrate by comparing low latitude and high latitude examples.
The presence of lateral rolling motion in erupting filaments can account for their observed bending and twisting during eruption.
The roll effect can depend, at least in part, on the asymmetries in the magnetic flux densities on the two sides of a polarity reversal boundary in a filament channel under a rising coronal loop system.
The magnitude of the filament deflection from the radial propagation during its eruption is typically greater than the deflection of the corresponding CME. This difference can be explained by the different origins and character of filaments and CMEs: (a) The spine of filaments exist as thin ribbons before their eruption and these ribbons roll, bend, and twist in characteristic ways during eruption. (b) Before their eruption, CME magnetic fields are tunnel-like systems of many fine coronal loops, a topology that allows an ascending loop system to reconnect and form a flux rope via magnetic reconnection below the associated and ascending filament ribbon. Such flux ropes therefore encapsulate the eruptive filaments, which still have some freedom to keep moving along the original direction of eruption, as defined prior to the flux rope formation.
We confirm previous findings that the CME phase of rapid acceleration is initiated concurrently with the magnetic reconnection of a rising coronal loop system.
In events with the roll effect and non-radial filament eruptions, but with radial CMEs, rolling motion and lateral deflection for the erupting filament is observed only during the early stages of the eruption.
CME flux ropes commonly propagate radially for eruptive events associated with pseudostreamers.
Non-radial and radial filament eruptions often show kinked shapes, which may, at least in part, be a consequence of projection effects from the rising ribbons of the filaments together with large scale rolling motions.
Through modeling and comparison with observed events, we anticipate that major twists and non-radial motions in erupting prominences and CMEs will become predictable to the extent that their environments are well-defined and measurable.
We are indebted to the SOHO, STEREO/SECCHI and SDO teams. O. P. and S. M. are supported in this research under the NASA grant NNX09AG27G. The work of M.V. was conducted at the Jet Propulsion Laboratory, California Institute of Technology under a contract from the National Aeronautics and Space Administration. A.V. is supported by NASA contract S-136361-Y to the Naval Research Laboratory. SOHO is a mission of international cooperation between ESA and NASA. The SECCHI data are produced by an international consortium of the NRL, LMSAL and NASA GSFC (USA), RAL and Univ. Bham (UK), MPS (Germany), CSL (Belgium), IOTA and IAS (France).The AIA data used here are courtesy of SDO (NASA) and the AIA consortium.
- Bemporad, A.: 2009, Astrophys. J. 701, 298.
- Brueckner, G. E., Howard, R. A., Koomen, M. J., Korendyke, C. M., Michels, D. J., Moses, J. D. et al.: 1995, Solar Phys. 162, 357.
- Cremades, H., Bothmer, V.: 2004, Astron. Astrophys. 422, 307.
- Cremades, H., Bothmer, V., Tripathi, D.: 2006, Adv. Space Res. 38, 461.
- Foukal, P. : 1971, Solar Phys. 19, 59.
- Filippov, B. P., Gopalswamy, N., Lozhechkin, A. V.: 2001, Solar Phys. 203, 119.
- Gaizauskas, V.: 1998, In: New Perspective on Solar Prominences, IAU Colloquium 167, ASP Conf. Series 150, 257.
- Gaizauskas, V., Zirker, J.B., Sweetland, C., Covács, A.: 1997, Astrophys. J. 479, 448.
- Gopalswamy, N., Shimojo, M., Lu, W., Yashiro, S., Shibasaki, K., Howard, R. A.: 2003, Astrophys. J. 586, 562.
- Gopalswamy, N., Mäkelä, P., Xie, H., Akiyama, S., Yashiro, S.: 2009, J.Geophys. Res. A 114, A22.
- Hood, A. W., Priest, E. R.: 1979, Solar Phys. 64, 303.
- Howard, R. A., Moses, J. D., Vourlidas, A., Newmark, J. S., Socker, D. G., Plunkett, S. P., et al.: 2008, Space Sci. Rev. 136, 67.
- Kaiser, M. L., Kucera, T. A., Davila, J. M., St. Cyr, O. C., Guhathakurta, M., Christian, E.: 2008, Space Sci. Rev. 136, 5.
- Kilpua, E. K. J., Pomoell, J., Vourlidas, A., Vainio, R., Luhmann, J., Li, Y., et al.: 2009, Ann. Geophys. 27, 4491.
- Lemen, J. R., Title, A. M., Akin, D. J., Boerner, P. F., Chou, C., Drake, J. F., et al.: 2012, Solar Phys. 275, 17.
- Liewer, P., Panasenco, O., Hall, J.R.: 2012, Solar Phys., DOI: 10.1007/s11207-012-0145-z.
- Liewer, P.C., de Jong, E.M., Hall, J.R., Howard, R.A., Thompson, W.T., Culhane, J.L., Bone, L., van Driel-Gesztelyi, L.: 2009, Solar Phys. 256, 57.
- Liewer, P.C., Hall, J.R., Howard, R.A., de Jong, E.M., Thompson, W.T., Thernisien, A.: 2011, J. Atm. Solar-Terr. Phys. 73, 1173.
- Martin, S.F.: 1990, In: Rudzjak, V., Tandberg-Hanssen, E. (eds.), Dynamics of Quiescent Prominences, Proc. of the 117th Colloquium of IAU, Hvar, Yugoslavia, Springer-Verlag, 1.
- Martin, S.F.: 1998, Solar Phys. 182, 107.
- Martin, S. F.: 2003, Adv. Space Res. 32, 1883.
- Martin, S.F., Lin, Y., Engvold, O.: 2008, Solar Phys. 250, 31.
- Martin, S.F., Panasenco, O., Engvold, O., Lin, Y.: 2008, Ann. Geophys. 26, 3061.
- Panasenco, O., Martin, S.: 2008, In: Howe, R., Komm, R.W., Balasubramaniam, K.S, Petrie, G.J.D. (eds.), Subsurface and Atmospheric Influences on Solar Activity, ASP Conf. Ser. 383, 243.
- Panasenco, O., Velli, M.M.: 2010, AGU Fall Meeting Abstracts, A1663.
- Panasenco, O., Martin, S., Joshi, A. D., Srivastava, N.: 2011, J. Atmos. Solar-Terr. Phys. 73, 1129.
- Parker, E. N.: 1957, Astrophys. J. Suppl. 3, 51.
- Pesnell, W. D., Thompson, B. J., Chamberlin, P. C.: 2012, Solar Phys. 275, 3.
- Pevtsov, A. A., Panasenco, O., Martin, S. F.: 2012, Solar Phys. 277, 185.
- Plunkett, S. P., Thompson, B. J., St. Cyr, O. C., & Howard, R. A. 2001, J. Atmos. Solar-Terr. Phys. 63, 389.
- Pneuman, G. W., Cargill, P. J.: 1985, Astrophys. J. 288, 653.
- Rappazzo, A. F., Velli, M., Einaudi, G., Dahlburg, R. B.: 2005, Astrophys. J. 633, 474.
- Schrijver, C. J., De Rosa, M. L.: 2003, Solar Phys. 212, 165.
- Schrijver, C. J., Title, A. M.: 2011, J. Geophys. Res. (Space Phys.) 116, 4108.
- Sedov, L. I.: 1981, Similarity and Dimensionality Methods in Mechanics, 9th ed. (English translation MIR, Moscow 1985).
- Sheeley, N. R., Jr., Warren, H. P.: 2012, Astrophys. J. 749, 40.
- Smith, S. F.: 1968, In: Kiepenheuer, K. O. (ed.), Structure and Development of Solar Active Regions, Symposium no. 35, D. Reidel Pub. Co., 267.
- Thompson, W.T., Kliem, B., Török, T.: 2012, Solar Phys. 276, 241.
- Török, T., Berger, M. A., Kliem, B.: 2010, Astron. Astrophys. 516, A49.
- Török, T., Panasenco, O., Titov, V. S., et al.: 2011, Astrophys. J. Lett. 739, L63.
- Velli, M., Hood, A. W., Einaudi, G.: 1990, Astrophys. J. 350, 428.
- Wang, Y.-M., Muglach, K.: 2007, Astrophys. J. 666, 1284.
- Wang, Y.-M., Muglach, K., Kliem, B.: 2009, Astrophys. J. 699, 133.
- Wang, Y.-M., Sheeley, N. R., Jr., Rich, N. B.: 2007, Astrophys. J. 658, 1340.
- Wood, P., Martens, P.: 2003, Solar Phys. 218, 123.
- Yeates, A. R., Mackay, D. H., van Ballegooijen, A. A.: 2007, Solar Phys. 245, 87.
- Zuccarello, F. P., Bemporad, A., Jacobs, C., Mierla, M., Poedts, S., Zuccarello, F.: 2012, Astrophys. J. 744, 66.