Formation of a Double-decker Magnetic Flux Rope in the Sigmoidal Solar Active Region 11520
In this paper, we address the formation of a magnetic flux rope (MFR) that erupted on 2012 July 12 and caused a strong geomagnetic storm event on July 15. Through analyzing the long-term evolution of the associated active region observed by the Atmospheric Imaging Assembly and the Helioseismic and Magnetic Imager on board the Solar Dynamics Observatory, it is found that the twisted field of an MFR, indicated by a continuous S-shaped sigmoid, is built up from two groups of sheared arcades near the main polarity inversion line half day before the eruption. The temperature within the twisted field and sheared arcades is higher than that of the ambient volume, suggesting that magnetic reconnection most likely works there. The driver behind the reconnection is attributed to shearing and converging motions at magnetic footpoints with velocities in the range of 0.1–0.6 km s. The rotation of the preceding sunspot also contributes to the MFR buildup. Extrapolated three-dimensional non-linear force-free field structures further reveal the locations of the reconnection to be in a bald-patch region and in a hyperbolic flux tube. About two hours before the eruption, indications for a second MFR in the form of an S-shaped hot channel are seen. It lies above the original MFR that continuously exists and includes a filament. The whole structure thus makes up a stable double-decker MFR system for hours prior to the eruption. Eventually, after entering the domain of instability, the high-lying MFR impulsively erupts to generate a fast coronal mass ejection and X-class flare; while the low-lying MFR remains behind and continuously maintains the sigmoidicity of the active region.
Subject headings:Sun: corona — Sun: coronal mass ejections (CMEs) — Sun: magnetic fields — Sun: filaments, prominences
Online-only material: animations, color figures
A magnetic flux rope (MFR) is defined as a current channel with a set of magnetic field lines wrapping more than once around the central axis. Such a structure often erupts from the Sun as a coronal mass ejection (CME; the largest-scale eruption in the solar system), which releases a large quantity of magnetized plasma into the interplanetary space with a velocity of hundreds km s, even up to 3000 km s (Yashiro et al., 2004; Zhang & Dere, 2006; Chen, 2011). The magnetized plasma is still organized frequently as a coherent MFR when arriving at the Earth, as indicated by features such as the magnetic field rotation, the drop of the solar wind density and proton temperature, and the low plasma beta in the in situ observations (Burlaga et al., 1981).
Several lines of evidence imply that the MFR exists in the corona prior to the CME eruption. Sigmoid, a forward or reversed S-shaped emission pattern in soft X-ray (SXR) and extreme ultraviolet (EUV) passbands, often appears in CME-productive active regions (ARs) (Rust & Kumar, 1996; Canfield et al., 1999). The straight section in the middle of the sigmoid is believed to be strong evidence of the MFR existing in the corona (e.g., Sterling et al., 2000; Liu et al., 2007; McKenzie & Canfield, 2008; Liu et al., 2010; Savcheva et al., 2012b). Filaments are another piece of evidence of the existence of the MFR, which includes magnetic dips that are able to collect cool material against gravity (Mackay et al., 2010; Guo et al., 2010; Su et al., 2011; Su & van Ballegooijen, 2012). Filament channels are even thought to be the body of the MFR because rotation motions are often observed at the bottom of dark cavities, the cross sections of filament channels at the solar limb (Low & Hundhausen, 1995; Guo & Wu, 1998; Gibson et al., 2004; Wang & Stenborg, 2010; Li et al., 2012). Moreover, taking advantage of the observed photospheric vector magnetic field at the bottom boundary, the MFR can be reconstructed by extrapolation techniques using the assumption of a non-linear force free field (NLFFF); this has also indicated preexistence of the MFR (e.g., Yan et al., 2001; Canou et al., 2009; Guo et al., 2010; Cheng et al., 2010, 2013b; Su et al., 2011; Jiang et al., 2013, 2014; Inoue et al., 2013).
If the MFR really exists in the corona prior to eruption, the question arises of when and where the MFR is built up. Two possibilities are proposed theoretically. One is that the MFR is generated in the convection zone and partly emerges into the corona by buoyancy (Fan, 2001; Martínez-Sykora et al., 2008). Manchester et al. (2004) noted that when the primary axis approaches the photosphere, the MFR may split into two parts by the reconnection with the surrounding fields, which only allows the upper part of the MFR to ascend to the corona (also see Magara, 2006; Archontis & Török, 2008; Leake et al., 2013). The other possibility is that the MFR is built up directly in the corona. Through imposing shearing and/or vortex motions to the different polarities, an initial potential field is gradually sheared. Converging flows then initiate the reconnection near the polarity inversion line (PIL), which converts sheared fields into twisted ones (van Ballegooijen & Martens, 1989; Amari et al., 2003). Depending on the specific locations of the reconnection, the MFR is created either through flux cancellation in the photosphere prior to the eruption (Aulanier et al., 2010; Amari et al., 2011; Xia et al., 2014) or via tether-cutting and flare reconnection in the corona during the eruption (Moore et al., 2001; Antiochos et al., 1999; Lynch et al., 2008; Karpen et al., 2012).
Observationally, the MFR is also conjectured either to stem from emergence from below the photosphere or to build up in the corona. Lites (2005) studied the properties of the vector magnetograms associated with two filaments and found a concave up geometry underneath the filaments. Okamoto et al. (2008) examined a sequence of vector magnetograms of AR 10953 and found that two abutting opposite-polarity regions with horizontally strong but vertically weak magnetic fields grew laterally and then narrowed. The directions of the horizontal magnetic fields along the PIL gradually reversed from a normal polarity to an inverse one. Both concave up geometry and reversed polarity suggest that the MFR may come from below the photosphere. However, Vargas Domínguez et al. (2012) recently provided an opposite interpretation for the photospheric evolution characteristics of the AR 10953. Through comparing with the numerical results by MacTaggart & Hood (2010), they stated that magnetic cancellation is also able to produce the lateral growing and then narrowing of the opposite polarities, as well as the reversal of the horizontal field direction. Moreover, through analyzing the temperature structure of a sigmoid, Tripathi et al. (2009) discovered that the plasma in the J-shaped arcades can have a higher temperature than that in the S-shaped flux if both are simultaneously visible. They argued that it is most likely that the J-shaped arcades are reconnecting to the S-shaped flux. The reconnection at the same time heats the plasma, which afterwards enters a cooling phase. Green & Kliem (2009) and Green et al. (2011) supported the conjecture that the reconnection is mostly associated with flux cancellation, although only part of the cancelled flux may be injected into the MFR. It is worth noting that the MFR can even be formed during a confined flare and be destabilized in a subsequent major eruption (Patsourakos et al., 2013; Song et al., 2014). This is similar to the conjecture of Guo et al. (2013), who argued that the quasi-separatrix layer reconnection in the interface between the MFR and the surrounding fields, indicated by a series of confined flares, has an important role in injecting self-helicity to the MFR.
Although previous works have displayed elementary results on the formation of the MFR, their objects for study are filaments or sigmoids, not the MFR itself. Recently, using the data from the Atmospheric Imaging Assembly (AIA; Lemen et al., 2012) on board the Solar Dynamics Observatory (SDO), Zhang et al. (2012) and Cheng et al. (2013a) found that the MFR directly exists as an elongated EUV channel structure appearing in the high temperature AIA passbands at 131 and 94 Å. In the impulsive acceleration phase, the MFR is further enhanced by flare reconnection, and the morphology evolves from the sigmoidal shape to the semicircular one. A number of further observations have indicated that this hot channel is actually the MFR that plays an important role in forming and accelerating the CME (Zhang et al., 2012; Cheng et al., 2013a, 2014; Li & Zhang, 2013). In this paper, we pay our attention to the long-term formation process of a hot channel in the sigmoidal AR NOAA 11520. We find that (1) the MFR is most likely formed through coexisting reconnection at a bald patch (BP) and at a hyperbolic flux tube, driven by photospheric shearing and converging flows; (2) a second high-lying MFR (the hot channel) locates vertically above a filament-associated low-lying MFR; the whole structure thus constitutes a double-decker MFR system that stably exists for about two hours prior to the eruption. The instruments are presented in Section 2. Observations of the MFR formation are displayed in Section 3, followed by the causes of the MFR formation in Section 4. In Section 5, we give our summary and discussions.
The data used in this study are mainly from the AIA (Lemen et al., 2012) and the Helioseismic and Magnetic Imager (HMI; Schou et al., 2012), both of which are on board the SDO. The AIA includes ten passbands, six of which image the solar corona almost simultaneously with an unprecedented high cadence (12 seconds) and high spatial resolution (1.2″), covering the temperature range from 0.06 MK to 20 MK (O’Dwyer et al., 2010). The HMI observes the vector magnetic field of the full solar photosphere with approximately the same spatial resolution (1.0″) as the AIA but with a cadence of 12 minutes. The X-Ray Telescope (XRT; Golub et al., 2007) aboard Hinode (Kosugi et al., 2007) images the hot plasma in the corona. The Geostationary Operational Environmental Satellite (GOES) and the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI; Lin et al., 2002) spacecraft register the SXR and hard X-ray (HXR) fluxes of solar flares. In addition, the Large Angle and Spectrometric Coronagraph (LASCO; Brueckner et al., 1995) on board the Solar and Heliospheric Observatory (SOHO) and the Sun-Earth Connection Coronal and Heliospheric Investigation (SECCHI; Howard et al., 2008) on board the Solar Terrestrial Relations Observatory (STEREO-A and STEREO-B) provide the EUV and white-light images of the CME.
|July 7 00:00 UT||The AR 11520 rotated from backside to the east limb.|
|July 11 00:00 UT||A sigmoidal structure appeared, the core field mainly consists of two sets of hot sheared arcades (Figure 1a–1c, Figure 6g, and Figure 8a).|
|July 11 00:00 – July 12 03:00 UT||The sheared arcades transformed to the continuous S-shaped field lines (Figure 1d–1f, Figure 6h, and Figure 8b) driven by the shear and convergence flows near the main neutral line (Figure 7).|
|July 12 14:00 UT||A diffuse and elongated high-lying channel appeared (Figure 2a–2c and Figure 3).|
|July 12 16:10 UT||The hot channel started to expand and rise (Figure 2g–2i), resulting in a CME (Figure 4) and X1.4 class flare, but the filament survived (Figure 5c).|
|July 13 14:00 UT||The low-lying MFR and associated filament could be identified with the disappearance of flare loops (Figure 5f and Figure 8d).|
|July 15 06:00 UT||The high-lying MFR-associated magnetic cloud arrived at the Earth.|
|July 15 19:00 UT||The Dst index peaks at nT .|
3. Observations of the MFR
3.1. Formation of double-decker MFR
On 2012 July 15, a magnetic cloud reached the Earth at 06:00 UT and then gave rise to a strong geomagnetic storm event with the minimal Dst index of –127 nT (19:00 UT). Inspecting the AIA data carefully, we find that the magnetic cloud was related to a strong solar eruption that occurred on July 12 and produced an X1.4 class flare and a fast halo CME. The flare was located at the heliographic coordinates S17W08. The corresponding source region is NOAA active region (AR) 11520. The flare SXR flux started to rise gradually at 14:50 UT, then rapidly increased from 16:10 UT and peaked at 16:49 UT.
In order to investigate the origin of the eruption, we examine the AIA images of the AR in all EUV passbands from July 6 to July 13. The overall evolution of the event is summarized in Table 1. It is found that the AR appeared at the solar limb on July 7 and progressively displayed a sigmoidal structure from July 11 onward. From the XRT Ti-poly and the AIA 94 Å images (Figure 1a–1b), one can see that initially the sigmoid was mainly composed of two groups of J-shaped arcades, whose straight arms were located at the two sides of the PIL with the elbows crossing the PIL at the opposite ends. From 03:00 UT on 2012 July 12, the two groups of J-shaped arcades gradually transformed to continuous S-shaped field lines (Figure 1d–1e), which is a strong indication that an MFR has formed (e.g., Savcheva et al., 2012a). In the AIA 335 Å images, we only see an overall outer envelope of the sigmoid, in which the detailed evolution of the coronal structure is difficult to detect (Figure 1c and 1f). These results show that the formation of the MFR by reconnection of J-shaped arcades mainly occurs in the core field of the sigmoid, where the plasma has been heated to the high temperatures (10 MK); while for the envelope of the sigmoid, the plasma temperature is relatively lower (2.5 MK), thus not indicating any strong reconnection. Moreover, in the AIA 304 Å passband, we notice that a J-shaped filament was visible from June 6 onward, located in the right part of the sigmoid. Its presence also indicates high shear or even twist of the sigmoid.
A closer inspection of the data leads us to the conclusion that a high-lying MFR and a low-lying MFR coexisted above the same PIL of the AR for two hours prior to the eruption. From 14:00 UT on July 12, we can see a diffuse and hot EUV channel structure in the AIA 94 Å passband (elongated S-shaped; outlined by the red dotted line in Figure 2c). The temperature of the channel is higher than 6 MK because it can only be seen in the hot temperature passbands (i.e., AIA 94 and 131 Å and XRT Ti-poly) but not in the other lower temperature ones (Figure 2d). From 16:00 UT, the hot channel showed an accelerating expansion (Figure 2g–2i), with the elbows expanding eastward and westward and the middle moving to the south. Because the channel is diffuse and less bright than the flare signatures, only the animation of the AIA 94 Å images that accompanies Figure 2 permits a full appreciation of the hot channel’s shape and dynamics. While the hot channel erupted, the filament (as seen in Figure 2e) stayed in place; it did not show any displacement larger than the slight changes in position seen during the preceding days. Such partial sigmoid eruptions are not uncommon (Pevtsov, 2002) and have been interpreted as a partial eruption of an MFR, whose top part has become unstable while the bottom part is tied to the photosphere in a BP (Gilbert et al., 2001; Gibson & Fan, 2006, 2008; Green & Kliem, 2009). While this is a plausible scenario in general, here it runs into difficulties because the HMI vector field data presented below show a BP section of the PIL only under the shorter left part of the filament (see, e.g., Figure 9a below). Moreover, at least some motion of the filament was likely to occur if it was part of the same erupting MFR as the hot channel; however, no motion in the southward direction of the eruption was seen. Next we consider the arcade of loops that dominated the middle and left part of the sigmoid in the AIA 94 Å images prior to the eruption of the hot channel (marked by white arrows in Figure 2c). These loops did not show any systematic or significant change in position or shape while the hot channel erupted (see Figure 2g–2i and the accompanying animation). Such a behavior is clearly at variance with the interpretation that a single MFR erupted in part. Finally, we consider the perspective of STEREO-B (Figure 3a) and find that the filament was blocked behind the solar limb, indicating it is low-lying; while the counterpart of the hot channel in the EUV 195 Å passband, probably from the emission of the Fe xxiv line at 192 Å (Milligan & McElroy, 2013), is located at a height of 90 Mm above the limb in the period of 14:00–16:00 UT, which is far above typical heights of active-region filaments (e.g., Tandberg-Hanssen, 1995). Thus, we are led to conclude that the hot channel was independent of the magnetic flux of the filament-associated low-lying MFR and constituted a double-decker MFR system already prior to the onset of the impulsive phase of the eruption. Note that the identification of the structure of the filament with a low-lying MFR is further supported by the NLFFF modeling in Section 4.3.
3.2. Partial Eruption of double-decker MFR
The explosive eruption of the high-lying MFR commenced at 16:10 UT, which resulted in the quick formation and acceleration of a CME and the rapid enhancement of the flare emission at various wavelengths. The CME first appeared as a limb event in the field-of-view (FOV) of the SECCHI coronagraphs COR1-A and COR1-B at 16:25 UT, and as a halo CME in the FOV of LASCO/C2 at 17:00 UT (Figure 4). As for the flare, we note that the most remarkable feature from 15:00 UT to 16:10 UT is the brightening of the right arcade of the sigmoid and of the footpoints of the double-decker MFR, which are closely related to the slow rise and accelerating expansion of the high-lying MFR (Figure 2g–2i and attached high cadence movies). After 16:10 UT, the footpoints brightenings extended along the direction of the PIL. At 16:25 UT, the brightenings formed two ribbons in the chromosphere, which separated from each other afterwards as seen in the movies of the AIA 304, 1600, and 1700 Å passbands (see details in attached high cadence movies of Figure 2). From the AIA 94 Å and XRT Ti-poly images, we can see the arcade of flare loops, the ends of which corresponded well to the two separating ribbons (Figure 5a–5b). Underneath the flare loops, the whole filament stayed quietly with the inferred low-lying MFR regardless of the eruption of the high-lying MFR. Moreover, we find that the CME in the COR1 images (Figure 4) mainly consisted of a bright front and a relatively diffuse cavity. We cannot, however, detect a clear appearance of a bright core that is generally thought to be erupted filament material and located at the bottom part of the cavity (Illing & Hundhausen, 1983; Vourlidas et al., 2013). This seems to also support that the filament did not erupt, or at least did not fully erupt. Mostly, the whole filament structure was not influenced much by the high-lying MFR eruption.
With the cooling and disappearance of the flare loops, the signatures of the low-lying MFR again showed up and redisplayed the sigmoidicity of the AR on July 13. In the high temperature passbands (AIA 94 Å and XRT Ti-poly; Figure 5d and 5e), one can obviously see that some S-shaped field lines became visible in the AR. The survived filament was located in the middle of the sigmoid as seen in the AIA 304 Å passband (Figure 5f).
4. Origin of the MFR
4.1. DEM Properties of MFR
Thanks to the multi-passband multi-temperature imaging ability of AIA, differential emission measure (DEM) structures of the plasma in the AR can be constructed. The observed flux for each AIA passband is given by =d, where denotes the temperature response function of passband . In order to reduce the error of the DEM inversion, we first calibrate six near-simultaneous AIA EUV images to the data level 1.5 using the SolarSoft routine “aia_prep.pro” and then degrade the resolution to 2.4 by the routine “rebin.pro”, thus guaranteeing a good coalignment accuracy of 0.6″ between the images in different passbands (Aschwanden et al., 2013). Using the routine “xrt_dem_iterative2.pro” as proposed by Weber et al. (2004) and Golub et al. (2004), we reconstruct the DEM in each pixel. The validation of this inversion method and uncertainties of the DEM results can be found in Cheng et al. (2012).
With the DEM in each pixel, we calculate the emission measure (EM) at the different temperature intervals () through the formula EM()=d to construct the two-dimensional maps of the plasma EM. Figure 6 shows the EM structure of the AR in three temperature ranges. One can see that the emission in the ambient volume of the AR is dominated by cool plasma (1–2 MK; upper row), whereas the emission of the sigmoidal structure is mostly from warm plasma (3–4 MK; middle row). Hot plasma (8–10 MK) contributes the significant emission in the sigmoid center (bottom row).
The EM maps of the AR provide important clues for understanding the formation of the MFR. From the EM maps in the early phase of the sigmoid, e.g., at 03:00 UT on July 11 (Figure 6d), one can only see an indication of the sigmoidal emission pattern at the warm temperature; while at the hot temperature, two groups of clearly sheared sigmoidal arcades have already appeared. Their maximum EM is 10 cm, being comparable to that of the loops in the warm temperature range. With the time elapsing, a clear sigmoidal emission pattern appeared at 07:00 UT on July 12 in the warm and hot temperature ranges, possibly due to the fact that magnetic flux carrying heated plasma was added to the continuous S-shaped field (Figure 6h). This demonstrates the gradual pre-eruption evolution of the AR, forming an MFR by reconnection. The low-lying MFR is very well captured in the EM maps, and a trace of the high-lying MFR can be identified as well in Figure 6i by the correspondence with the diffuse S-shaped structure marked by the red line in the AIA 94 Å image in Figure 2c. Finally, we note that the peak EM of the low-lying MFR at 10 MK decreased in Figure 6i (to 10 cm); this may correspond to the temporarily reduced rate of flux cancellation (Figure 7a).
4.2. Shearing and Converging Flows and Sunspot Rotation
In this section, we study the photospheric properties during the MFR formation. We first plot the evolution of the line-of-sight magnetic flux of the sigmoidal AR in Figure 7a. During 40 hours before the eruption, the positive flux increased from 3.610 Mx to 4.010 Mx, while the negative flux decreased slightly from 1.710 Mx to 1.510 Mx. Overall, there is no significant flux emergence in the period of the MFR formation (Table 1 and Figure 7a). After the eruption, the positive flux approximately kept a constant value of 4.010 Mx, whereas the negative flux quickly decreased to 1.010 Mx at 00:00 UT on July 14. The decrease is most likely to be attributed to magnetic cancellation along the main PIL as shown in Figure 7f–7i.
In order to investigate the driver of the magnetic cancellation, we calculate the transverse velocity of flux elements in the photosphere through applying the differential affine velocity estimator (DAVE) to a temporal sequence of HMI line-of-sight magnetograms. This technique assumes an affine velocity profile within a small window and then executes variational principles to minimize the deviations in the magnitude of the magnetic induction equation (Schuck, 2005, 2006). The only free parameter is the window size, which is set to 10 pixels in our calculation. According to the tests by Chae & Sakurai (2008), this window size ensures a relative error less than 20%.
We find that the MFR formation in the sigmoid is associated with three main types of photospheric motions: shearing, converging, and rotational flows (see Figure 7f–7i and the accompanying movie). In the center of the sigmoid there are strong flows toward the PIL, mainly from the positive-polarity side. These consist of the moat flow of the preceding sunspot and of the motion of the sunspot as a whole toward the PIL, both of which persist throughout the time range analyzed. These flows lead to strong magnetic cancellation in the middle of the sigmoid, between the straight legs of the initial double-J pair of arcades. Under the left part of the sigmoid, positive flux is intruding along the PIL toward the center of the AR with a velocity of 0.3 km s; while the adjacent negative flux shows motions in the opposite direction at a velocity of 0.2 km s during part of the time. This leads to strong shear flows at a section of the PIL that likely cause significant cancellation, since the polarities are closely butted against each other while possessing a strong relative motion. The shearing motion stretched the loops to form the two sets of J-shaped arcades (Figure 7b and 7f). The converging flows initiated the reconnection between their heads and tails, producing the twisted field lines as suggested by Martens & Zwaan (2001). With the continuous photosphere-driven reconnection, more and more J-shaped arcades are converted to the twisted flux and came into being the MFR (Figure 7c and 7g).
In addition to the shearing and converging motions, the rotation of the sunspot also has a significant role in building up the MFR. From Figure 7g and 7h, one can see that the preceding sunspot, where the right footpoints of both conjectured MFRs were anchored, has an obvious rotation motion (indicated in the red circles). The rotation angular velocity (3 h) results in a maximal flow velocity of 0.6 km s at the edge of the AR where the footpoints of the MFR are located. The clockwise rotation twists the coronal field rooted in this sunspot in the right-handed sense, thus supporting the formation of the MFR. It is however difficult to identify the origin of the rotation, which can be due to either a vortex motion in the photosphere or the emergence of a strongly twisted flux tube, although the probability of the latter is small for a mature AR.
4.3. Non-Linear Force-Free Field Modeling
We use an optimization algorithm proposed by Wheatland et al. (2000) and implemented by Wiegelmann (2004) to extrapolate the three-dimensional (3D) NLFFF structure of AR 11520. Due to the plasma pressure being dominated by the magnetic pressure in the corona (1; Gary, 2001), the coronal magnetic field generally satisfies the force-free criteria, i.e., = and =0, where is a constant along each field line. For the magnetic field above ARs, varies in space (Wiegelmann & Neukirch, 2002). The optimization algorithm minimizes the objective function = through iteration and thus approaches the solution of the NLFFF equations (Wiegelmann, 2004). Before the extrapolation, we apply a preprocessing procedure to the bottom boundary vector data. This removes most of the net force and torque that otherwise generally results in an inconsistency between the forced photospheric magnetic field and the force-free assumption in the NLFFF models (Wiegelmann et al., 2006).
We compute a time sequence of the 3D NLFFF structure of the AR, covering a period of three days with a cadence of 1 hour. The 3D magnetic field structures at four instants are shown in Figure 8a–8h. One can see that at 03:00 UT on July 11, the AR included three sets of field lines: two sets of strongly sheared arcades near the PIL and the overlying constraining field. The right and left arcades correspond very well to the observed double-J sigmoid if seen from above (Figures 8a and 2a–b). By 07:00 UT on July 12, part of the two groups of arcade field lines may have reconnected, as manifested by continuous sigmoidal field lines (Figure 8b), which form a weakly twisted flux rope. At 15:00 UT on July 12, the twist in the rope reached a maximum, since more and more flux was added to the rope by the ongoing reconnection and the rotation of the sunspot (Figure 8c and 8g). At 12:00 UT on July 13, in spite of the significant decrease of the twist as a consequence of the eruption, the remaining twist still preserved the sigmoidal structure (Figure 8d and 8h). Overall, the time sequence of 3D NLFFF structures successfully reproduces the evolution of the sigmoid, including the formation and twisting of an MFR before the eruption, as well as the survival of a weakly twisted MFR after the eruption.
The NLFFF structures provide the significant clues to the question whether the stable filament resides in a magnetic arcade or in a low-lying MFR. The HMI vector data show that the horizontal field in the section of the PIL under the left half of the sigmoid pointed in the inverse direction (from the negative to the positive side of the PIL). This BP signifies the topology of an MFR that must exist at least in this section of the PIL. The magnetic field modeling suggests that this MFR extended along the whole length of the filament. We conjecture that the sigmoid involved two MFRs, at least in the final two hours before the eruption during which the hot channel was seen. One MFR was high-lying and the other was low-lying (holding the filament); both existed simultaneously above the PIL, constituting a stable double-decker MFR system. From the AIA 94 and 131 Å and XRT images, it is obvious that the two branches of the double-decker system had very closely located footpoints at both ends but rather different lengths (Figure 2c).
Furthermore, to study the properties of the reconnection during the MFR formation in detail, we plot three representative field lines DF, BC, and AE in Figure 9a and 9b. DF and BC denote the right and left sheared arcades, respectively. AE shows the S-shaped field of the MFR. To further examine the properties of the reconnection at the different locations, we take three north-south oriented cross sections at =, , and (Figure 9a). The distributions of (the total current density normalized by the total magnetic field) at =, , and are shown in Figure 9c–9e, respectively. Here, the current density is given by , where =4 10 G m A. One can see that at the cross section = is mainly concentrated at the BPs (Figure 9a–9c), where the horizontal photospheric field components show inverse polarity (Figure 9a) and the field lines are concave up with their bottom points touching the photosphere (Figure 9b). At the cross section =, the distribution of displays an X-shape at the height of 4 Mm, suggesting a hyperbolic flux tube (HFT) in 3D; this indicates that the reconnection mainly takes place in the corona (Figure 9d). At the cross section =, the location of high coronal current density ascends to a higher position (15 Mm) and is mostly concentrated inside the MFR; thus this current has likely a role in heating the MFR. Based on the above properties, we argue that different types of reconnection, i.e., BP, HFT, and internal reconnection, probably exist simultaneously during the formation of the MFR, either generating the twist or heating the plasma.
We also compare the distribution of the currents with the emission pattern in the sigmoid. Figure 8i–8l show the distributions of the current density integrated along the line of sight. It can be seen that the concentration of the currents is mostly along the MFR axis with the largest magnitude appearing at the regions corresponding to the BPs and the HFT. This fact infers that it is the currents that heat up the plasma inside and around the sigmoidal MFR, thus making a sigmoidal emission pattern. Moreover, one can notice that the current density near the main PIL increases before the eruption (Figure 8i–j) and decreases afterwards (Figure 8k–l). This can be qualitatively explained by the partial eruption of the configuration which released part of the coronal currents.
In addition, we note that the NLFFF model is difficult to reproduce the high-lying MFR in spite of its significant success in reproducing the long-term evolution of the sigmoidal AR. Before the eruption, both the two sets of sheared arcades and the continuous S-shaped and twisted flux bundle are successfully reproduced by the NLFFF model (Figure 8a and 8b). Even after the eruption, the NLFFF model is still powerful to reconstruct the MFR with a decreased twist that strongly resembles the surviving filament (Figure 8d). However, the NLFFF model does not fully succeed in reconstructing the complex magnetic structure near the eruption. The extrapolated 3D structure at 15:00 UT on July 12 only shows a low-lying twisted structure (Figure 8c), but misses the high-lying MFR that separated from the low-lying twisted field and appeared as the elongated hot channel structure in the AIA 131 and 94 Å passbands. A possible reason for this partial failure can be ascribed to the proximity of the footpoints of the two twisted coronal structures (Figure 2c and Figure 7). For a successful reconstruction of both structures, all four footpoint areas must be well resolved by the vector magnetogram, which obviously was not realized.
5. Summary and Discussion
In this paper, we investigate the origin of an MFR in the sigmoidal AR that erupted at about 16:10 UT on 2012 July 12 and produced an interplanetary magnetic cloud that caused the strong geomagnetic storm event on July 15. We perform a detailed analysis of the AR’s evolution in the three days leading up to the eruption. This includes tracking its morphological evolution, diagnosing the DEM properties, characterizing the motions of magnetic footpoints in the photosphere, and constructing a time sequence of 3D NLFFF models of the coronal structure.
A particularly interesting finding is that two MFRs have been formed above the same PIL of the sigmoidal AR and constituted a stable double-decker MFR system for at least two hours prior to the eruption. The concept of the double-decker MFR was recently proposed by Liu et al. (2012) to explain two vertically separated filaments over the same PIL. These authors payed much attention to the identification of the double-decker filament and to the discussion of the mechanisms of the partial eruption of the system. In the present work, we present a second case that supports their new conjecture for the magnetic structure of some CME source regions and additionally concentrate on the formation process of the double-decker MFR. It is found that during a period of the first 40 hours prior to the eruption, an evolving sigmoid manifested the formation of one MFR, most likely the result of reconnection between two groups of sheared arcades near the main PIL. The driver of the reconnection is attributed to the shearing and converging photospheric flows in the vicinity of the PIL, as derived with the DAVE technique. The distribution of the current layers as indicated by the NLFFF extrapolation suggests that the reconnection happens simultaneously at the BPs, i.e., photospheric flux cancellation, and in the HFT in the corona, i.e., tether-cutting. In the present event both worked at the same time in the process of converting the sheared arcades to the twisted field. Besides the shearing and converging flows, the rotation of the leading sunspot probably also played a role in forming the MFR, similar to the vortex motions used for building up twist in some simulations (Amari et al., 2003; Török & Kliem, 2003; Aulanier et al., 2005). A set of continuous S-shaped hot threads indicates that an MFR structure was formed about half a day before the eruption. This low-lying MFR also hosted a filament and remained stable in the eruption.
From about two hours before the eruption we find evidence for the existence of a second MFR in the form of a hot channel (Zhang et al., 2012; Cheng et al., 2013a), which was located above the first MFR. Only the second MFR erupted in the event studied here. The second MFR could be observed by the AIA high-temperature and XRT passbands. It is also imaged in the 195 Å passband, indicating emission in the Fe xxiv line blend in this channel ( MK; Milligan & McElroy, 2013). The location above the first MFR excludes its formation by emergence from below the photosphere. The high temperature suggests an important role for reconnection in the formation. This and the proximity of the two MFR in the stable phase prior to the eruption suggest that the double-decker MFR system may have formed by a splitting of the initial and low-lying MFR through the internal reconnection. Such a splitting must be considered as a tentative interpretation, since an MFR is a coherent large-scale structure that generally possesses a considerable degree of stability against perturbations. The low-lying MFR is perhaps most evident from the typically long-lasting stability of quiescent filaments. However, indications that MFR can split or even completely disintegrate do exist, both in observations and numerical modeling. Prominences often show two branches that are clearly separated in height (e.g., Liu et al., 2012; Su et al., 2011). The disintegration of a sigmoid by flux dispersal was described in Tripathi et al. (2009). The vertical splitting of an MFR in the evolution to an eruption of only the upper part was found in the numerical modeling of the 1997 May 2 eruption (Kliem et al., 2013). This evolution was driven by photospheric flows converging at the PIL and enforcing flux cancellation, and it involved slow tether-cutting reconnection with the ambient field in an HFT formed between the splitting parts of the rope. The flux added between the two parts gradually stabilizes (destabilizes) the lower (upper) part of the splitting MFR. Kliem et al. (2013) further demonstrated that a double-decker MFR system can be in stable equilibrium if the overlying field is sufficiently strong. They also found that the configuration admits a partial eruption, with only the upper branch erupting and the bottom branch remaining stable, very similar to the eruption of the high-lying hot channel and the stability of the low-lying filament. In the present event, strong perturbations of the previously formed MFR in the sigmoid were given by the intrusion of flux along the PIL under the left part of the sigmoid and by the rotation of the sunspot under the right part of the sigmoid, both of which affected the footpoint regions of either MFR. The conjectured internal reconnection must have occurred high enough in the corona (e.g., the second case of Figure 2 in Gilbert et al., 2001) so that the remaining MFR could still support the filament and continuously maintain the sigmoidal pattern of the AR.
The partial eruption of the double-decker MFR system possesses some similarities to but also differences from a partial eruption of a single MFR (Gilbert et al., 2001; Gibson & Fan, 2006, 2008). In that case, the reconnection happens in the interface between the MFR and the surrounding fields, e.g., near the crossing point of a kinked MFR (Tripathi et al., 2013), which results in the escape of the upper part of the MFR with the lower part remaining at the original place. However, one should note that the driver of this internal reconnection is attributed to the helical kink instability, which writhes the MFR axis and generates a current sheet around the MFR (e.g., Kliem et al., 2010). In contrast, for the case of the double-decker MFR, photospheric shearing and converging motions drive the reconnection. Therefore, the timescale of the separation is also different; it is nearly instantaneous in the partial eruption model but lasts for hours in the present case.
In the early eruption phase, the morphology of the high-lying MFR varied from an S-shape to a loop-shape, which is very similar to the linear feature in erupting sigmoids (e.g., Moore et al., 2001; McKenzie & Canfield, 2008; Liu et al., 2010; Aulanier et al., 2010; Green et al., 2011; Zharkov et al., 2011). The high similarity suggests that the linear feature is most likely the MFR itself rather than the current shell above the MFR as suggested by Aulanier et al. (2010). The difficulty in identifying the linear feature with the MFR in previous studies can be attributed to the unavailability of the high temporal and spatial resolution data.
As the high-lying MFR slowly ascends to a height where the background field declines rapidly enough, the torus instability probably triggers and initiates the impulsive acceleration of the MFR eruption (Kliem & Török, 2006; Török & Kliem, 2005; Fan & Gibson, 2007; Aulanier et al., 2010; Olmedo & Zhang, 2010; Savcheva et al., 2012a; Cheng et al., 2013b, 2014; Dudík et al., 2014). The initial brightening at the footpoints of the MFR in all EUV and UV passbands is consistent with an enhancement of internal reconnection by the commencing eruption. This is followed by a seamless transition to the much more rapid reconnection in the flare current sheet, in a mutual feedback with the unstable MFR, which not only formed the flare loops further constraining the filament but also produced the high-energy particles that stream down along the newly reconnected loops to generate two well-observed flare ribbons.
Finally, we find that the NLFFF model of the coronal field, obtained by extrapolation from a sequence of HMI vector magnetograms, succeeds in simulating important aspects in the long-term quasi-static evolution of the sigmoidal AR. The model reproduces the formation of a twisted, sigmoidal flux rope from the highly sheared arcades in very good agreement with the observed coronal structures, and it also resembles the weaker sigmoidal structure after the partial eruption quite well. On the other hand, the double-decker magnetic configuration suggested by the coronal data could not be found. We conjecture that this results from insufficient resolution of the structure in the magnetogram because the footpoints of the two branches were located in close proximity. Moreover, such a sequence of static models may capture complex evolutions only if a much higher time resolution is realized. These facts suggest that more advanced models or MHD simulations should be developed in the future to deal with complex pre-eruption magnetic structures and their dynamic evolution.
- slugcomment: Accepted for publication in ApJ
- Amari, T., Aly, J.-J., Luciani, J.-F., Mikic, Z., & Linker, J. 2011, ApJ, 742, L27
- Amari, T., Luciani, J. F., Aly, J. J., Mikic, Z., & Linker, J. 2003, ApJ, 585, 1073
- Antiochos, S. K., DeVore, C. R., & Klimchuk, J. A. 1999, ApJ, 510, 485
- Archontis, V., & Török, T. 2008, A&A, 492, L35
- Aschwanden, M. J., Boerner, P., Schrijver, C. J., & Malanushenko, A. 2013, Sol. Phys., 283, 5
- Aulanier, G., Démoulin, P., & Grappin, R. 2005, A&A, 430, 1067
- Aulanier, G., Török, T., Démoulin, P., & DeLuca, E. E. 2010, ApJ, 708, 314
- Brueckner, G. E., Howard, R. A., Koomen, M. J., Korendyke, C. M., et al. 1995, Sol. Phys., 162, 357
- Burlaga, L., Sittler, E., Mariani, F., & Schwenn, R. 1981, J. Geophys. Res., 86, 6673
- Canfield, R. C., Hudson, H. S., & McKenzie, D. E. 1999, Geophys. Res. Lett., 26, 627
- Canou, A., Amari, T., Bommier, V., Schmieder, B., Aulanier, G., & Li, H. 2009, ApJ, 693, L27
- Chae, J., & Sakurai, T. 2008, ApJ, 689, 593
- Chen, P. F. 2011, Living Reviews in Solar Physics, 8, 1
- Cheng, X., Ding, M. D., Guo, Y., Zhang, J., Jing, J., & Wiegelmann, T. 2010, ApJ, 716, L68
- Cheng, X., Zhang, J., Ding, M. D., Liu, Y., & Poomvises, W. 2013a, ApJ, 763, 43
- Cheng, X., Zhang, J., Ding, M. D., Olmedo, O., Sun, X. D., Guo, Y., & Liu, Y. 2013b, ApJ, 769, L25
- Cheng, X., Zhang, J., Saar, S. H., & Ding, M. D. 2012, ApJ, 761, 62
- Cheng, X., et al. 2014, ApJ, 780, 28
- Dudík, J., Janvier, M., Aulanier, G., Del Zanna, G., Karlický, M., Mason, H., & Schmieder, B. 2014, ArXiv e-prints
- Fan, Y. 2001, ApJ, 554, L111
- Fan, Y., & Gibson, S. E. 2007, ApJ, 668, 1232
- Gary, G. A. 2001, Sol. Phys., 203, 71
- Gibson, S. E., & Fan, Y. 2006, ApJ, 637, L65
- —. 2008, Journal of Geophysical Research (Space Physics), 113, 9103
- Gibson, S. E., Fan, Y., Mandrini, C., Fisher, G., & Demoulin, P. 2004, ApJ, 617, 600
- Gilbert, H. R., Holzer, T. E., & Burkepile, J. T. 2001, ApJ, 549, 1221
- Golub, L., Deluca, E., Austin, G., Bookbinder, J., et al. 2007, Sol. Phys., 243, 63
- Golub, L., Deluca, E. E., Sette, A., & Weber, M. 2004, in Astronomical Society of the Pacific Conference Series, Vol. 325, The Solar-B Mission and the Forefront of Solar Physics, ed. T. Sakurai & T. Sekii, 217
- Green, L. M., & Kliem, B. 2009, ApJ, 700, L83
- Green, L. M., Kliem, B., & Wallace, A. J. 2011, A&A, 526, A2
- Guo, W. P., & Wu, S. T. 1998, ApJ, 494, 419
- Guo, Y., Ding, M. D., Cheng, X., Zhao, J., & Pariat, E. 2013, ApJ, 779, 157
- Guo, Y., Schmieder, B., Démoulin, P., Wiegelmann, T., Aulanier, G., Török, T., & Bommier, V. 2010, ApJ, 714, 343
- Howard, R. A., et al. 2008, Space Sci. Rev., 136, 67
- Illing, R. M. E., & Hundhausen, A. J. 1983, J. Geophys. Res., 88, 10210
- Inoue, S., Hayashi, K., Shiota, D., Magara, T., & Choe, G. S. 2013, ApJ, 770, 79
- Jiang, C., Feng, X., Wu, S. T., & Hu, Q. 2013, ApJ, 771, L30
- Jiang, C., Wu, S. T., Feng, X., & Hu, Q. 2014, ApJ, 780, 55
- Karpen, J. T., Antiochos, S. K., & DeVore, C. R. 2012, ApJ, 760, 81
- Kliem, B., Linton, M. G., Török, T., & Karlický, M. 2010, Sol. Phys., 266, 91
- Kliem, B., & Török, T. 2006, Phys. Rev. Lett., 96, 255002
- Kliem, B., Török, T., Titov, V. S., Lionello, R., Linker, J. A., Liu, R., Liu, C., & Wang, H. 2013, ApJ submitted
- Kosugi, T., et al. 2007, Sol. Phys., 243, 3
- Leake, J. E., Linton, M. G., & Török, T. 2013, ApJ, 778, 99
- Lemen, J. R., et al. 2012, Sol. Phys., 275, 17
- Li, L. P., & Zhang, J. 2013, A&A, 552, L11
- Li, X., Morgan, H., Leonard, D., & Jeska, L. 2012, ApJ, 752, L22
- Lin, R. P., Dennis, B. R., Hurford, G. J., Smith, D. M., et al. 2002, Sol. Phys., 210, 3
- Lites, B. W. 2005, ApJ, 622, 1275
- Liu, C., Lee, J., Yurchyshyn, V., Deng, N., Cho, K.-s., Karlický, M., & Wang, H. 2007, ApJ, 669, 1372
- Liu, R., Kliem, B., Török, T., Liu, C., Titov, V. S., Lionello, R., Linker, J. A., & Wang, H. 2012, ApJ, 756, 59
- Liu, R., Liu, C., Wang, S., Deng, N., & Wang, H. 2010, ApJ, 725, L84
- Low, B. C., & Hundhausen, J. R. 1995, ApJ, 443, 818
- Lynch, B. J., Antiochos, S. K., DeVore, C. R., Luhmann, J. G., & Zurbuchen, T. H. 2008, ApJ, 683, 1192
- Mackay, D. H., Karpen, J. T., Ballester, J. L., Schmieder, B., & Aulanier, G. 2010, Space Sci. Rev., 151, 333
- MacTaggart, D., & Hood, A. W. 2010, ApJ, 716, L219
- Magara, T. 2006, ApJ, 653, 1499
- Manchester, IV, W., Gombosi, T., DeZeeuw, D., & Fan, Y. 2004, ApJ, 610, 588
- Martens, P. C., & Zwaan, C. 2001, ApJ, 558, 872
- Martínez-Sykora, J., Hansteen, V., & Carlsson, M. 2008, ApJ, 679, 871
- McKenzie, D. E., & Canfield, R. C. 2008, A&A, 481, L65
- Milligan, R. O., & McElroy, S. A. 2013, ApJ, 777, 12
- Moore, R. L., Sterling, A. C., Hudson, H. S., & Lemen, J. R. 2001, ApJ, 552, 833
- O’Dwyer, B., Del Zanna, G., Mason, H. E., Weber, M. A., & Tripathi, D. 2010, A&A, 521, A21
- Okamoto, T. J., et al. 2008, ApJ, 673, L215
- Olmedo, O., & Zhang, J. 2010, ApJ, 718, 433
- Patsourakos, S., Vourlidas, A., & Stenborg, G. 2013, ApJ, 764, 125
- Pevtsov, A. A. 2002, Sol. Phys., 207, 111
- Rust, D. M., & Kumar, A. 1996, ApJ, 464, L199
- Savcheva, A., Pariat, E., van Ballegooijen, A., Aulanier, G., & DeLuca, E. 2012a, ApJ, 750, 15
- Savcheva, A. S., van Ballegooijen, A. A., & DeLuca, E. E. 2012b, ApJ, 744, 78
- Schou, J., et al. 2012, Sol. Phys., 275, 229
- Schuck, P. W. 2005, ApJ, 632, L53
- —. 2006, ApJ, 646, 1358
- Song, H. Q., Zhang, J., Cheng, X., Chen, Y., Liu, R., Wang, Y. M., & Li, B. 2014, ApJ, 784, 48
- Sterling, A. C., Hudson, H. S., Thompson, B. J., & Zarro, D. M. 2000, ApJ, 532, 628
- Su, Y., Surges, V., van Ballegooijen, A., DeLuca, E., & Golub, L. 2011, ApJ, 734, 53
- Su, Y., & van Ballegooijen, A. 2012, ApJ, 757, 168
- Tandberg-Hanssen, E., ed. 1995, Astrophysics and Space Science Library, Vol. 199, The nature of solar prominences
- Török, T., & Kliem, B. 2003, A&A, 406, 1043
- —. 2005, ApJ, 630, L97
- Tripathi, D., Kliem, B., Mason, H. E., Young, P. R., & Green, L. M. 2009, ApJ, 698, L27
- Tripathi, D., Reeves, K. K., Gibson, S. E., Srivastava, A., & Joshi, N. C. 2013, ApJ, 778, 142
- van Ballegooijen, A. A., & Martens, P. C. H. 1989, ApJ, 343, 971
- Vargas Domínguez, S., MacTaggart, D., Green, L., van Driel-Gesztelyi, L., & Hood, A. W. 2012, Sol. Phys., 278, 33
- Vourlidas, A., Lynch, B. J., Howard, R. A., & Li, Y. 2013, Sol. Phys., 284, 179
- Wang, Y.-M., & Stenborg, G. 2010, ApJ, 719, L181
- Weber, M. A., Deluca, E. E., Golub, L., & Sette, A. L. 2004, in IAU Symposium, Vol. 223, Multi-Wavelength Investigations of Solar Activity, ed. A. V. Stepanov, E. E. Benevolenskaya, & A. G. Kosovichev, 321–328
- Wheatland, M. S., Sturrock, P. A., & Roumeliotis, G. 2000, ApJ, 540, 1150
- Wiegelmann, T. 2004, Sol. Phys., 219, 87
- Wiegelmann, T., Inhester, B., & Sakurai, T. 2006, Sol. Phys., 233, 215
- Wiegelmann, T., & Neukirch, T. 2002, Sol. Phys., 208, 233
- Xia, C., Keppens, R., & Guo, Y. 2014, ApJ, 780, 130
- Yan, Y., Deng, Y., Karlický, M., Fu, Q., Wang, S., & Liu, Y. 2001, ApJ, 551, L115
- Yashiro, S., Gopalswamy, N., Michalek, G., St. Cyr, O. C., Plunkett, S. P., Rich, N. B., & Howard, R. A. 2004, Journal of Geophysical Research (Space Physics), 109, 7105
- Zhang, J., Cheng, X., & Ding, M.-D. 2012, Nature Communications, 3, 747
- Zhang, J., & Dere, K. P. 2006, ApJ, 649, 1100
- Zharkov, S., Green, L. M., Matthews, S. A., & Zharkova, V. V. 2011, ApJ, 741, L35