A Study of a Compound Solar Eruption with Two Consecutive Erupting Magnetic Structures

, , and

Published 2018 June 8 © 2018. The American Astronomical Society. All rights reserved.
, , Citation Suman K. Dhakal et al 2018 ApJ 860 35 DOI 10.3847/1538-4357/aac028

Download Article PDF
DownloadArticle ePub

You need an eReader or compatible software to experience the benefits of the ePub3 file format.

0004-637X/860/1/35

Abstract

We report a study of a compound solar eruption that was associated with two consecutively erupting magnetic structures and correspondingly two distinct peaks, during impulsive phase, of an M-class flare (M8.5). Simultaneous multi-viewpoint observations from SDO, GOES and STEREO-A show that this compound eruption originated from two pre-existing sigmoidal magnetic structures lying along the same polarity inversion line. Observations of the associated pre-existing filaments further show that these magnetic structures are lying one on top of the other, separated by 12 Mm in height, in a so-called "double-decker" configuration. The high-lying magnetic structure became unstable and erupted first, appearing as an expanding hot channel seen at extreme ultraviolet wavelengths. About 12 minutes later, the low-lying structure also started to erupt and moved at an even faster speed compared to the high-lying one. As a result, the two erupting structures interacted and merged with each other, appearing as a single coronal mass ejection in the outer corona. We find that the double-decker configuration is likely caused by the persistent shearing motion and flux cancellation along the source active region's strong-gradient polarity inversion line. The successive destabilization of these two separate but closely spaced magnetic structures, possibly in the form of magnetic flux ropes, led to a compound solar eruption. The study of the compound eruption provides a unique opportunity to reveal the formation process, initiation, and evolution of complex eruptive structures in solar active regions.

Export citation and abstract BibTeX RIS

1. Introduction

Solar eruptions, manifested as flares and coronal mass ejections (CMEs), are the most magnificent events happening on the solar atmosphere. They are the main sources of disturbances in interplanetary space and the driver of space weather near the Earth. Flares are sudden and localized emission of radiation over the entire electromagnetic spectrum. They are routinely identified and classified based on their intensity profile of soft X-ray (SXR) emission (Fletcher et al. 2011). CMEs, which have wider effects, are the eruption of magnetized plasma from the solar atmosphere identified as arc-like brightenings in coronagraph images in the outer corona. Flares and CMEs, often occurring together, are manifested as a single energy release process, during which a tremendous amount of stored magnetic energy is released from the solar atmosphere (Forbes 2000). The general model of eruptive flares (flares with CMEs) describes the solar eruption as the eruption of a magnetic flux rope (MFR). An MFR is a coherent structure of magnetic field twisting around a central axis that carries a large amount of free magnetic energy. A typical solar eruption has three phases in the X-ray profiles: (1) the precursor phase, (2) the impulsive phase (IP), and (3) the gradual phase (see Zhou et al. 2016 for a discussion). This three-phase evolution scenario has been also clearly observed in CMEs and is associated with three distinct acceleration phases of them (e.g., Zhang et al. 2001).

It has been suggested that an MFR can exist before the eruption (van Ballegooijen & Martens 1989; Chen et al. 1997) or it could form during it (Antiochos et al. 1999; Lynch et al. 2008). If an MFR exists before the eruption, the question is how and where it could form and how to detect it. It has been suggested that an MFR could form by different mechanisms prior to its eruption, such as shearing motions (Amari et al. 2000, 2003; Jacobs et al. 2009), flux cancellation at the polarity inversion line (PIL) (van Ballegooijen & Martens 1989; Aulanier et al. 2010; Green et al. 2011), flux emergence (Fan & Gibson 2003; Leake et al. 2013), and confined flaring (Patsourakos et al. 2013; Chintzoglou et al. 2015). The detection of an MFR is difficult as it has low plasma density, and also due to the technological inability to measure magnetic field directly in the corona. However, there is indirect evidence that may suggest the existence of MFRs in the form of a cavity and prominence (associated with filaments when seen on the solar disk), sigmoidal bright loops in X-rays (Rust & Kumar 1996; Green et al. 2007), AIA/SDO hot channel structures (Cheng et al. 2011; Zhang et al. 2012), dips or bald patches in filament channels (Lites 2005; López Ariste et al. 2006), and nonlinear force-free field (NLFFF) extrapolation results obtained for the photospheric boundary conditions (e.g., Chintzoglou et al. 2015).

The pre-erupting structures (MFRs) are kept in equilibrium by a downward force due to the pressure and magnetic tension from the overlying magnetic field (Chen 1989). The instability of MFRs could result in three scenarios: (i) failed eruption, (ii) partial eruption, and (iii) full eruption of the MFR. In a failed eruption, the MFR rises up as it is erupting and then stops at a certain maximum height (Song et al. 2014). A failed eruption could be due to many factors such as strong overlying field pressure (Sun et al. 2015), an asymmetric overlying coronal field above the filament (Liu et al. 2009), or not having enough energy for a successful eruption (Shen et al. 2011). In partial eruptions, (a) reconnection within the MFR (Gibson & Fan 2006) separates it into two parts, or (b) the partial emergence and reconnection with the coronal dipole field may result in the split of the MFR into two parts (Leake et al. 2014). In partial eruptions, the upper part of the MFR erupts and its lower part remains in place. Non-erupting filaments and the immediate re-formation of long-lived X-ray sigmoids after eruption were considered to be partial eruptions (Tang 1986; Gilbert et al. 2000; Gibson et al. 2002). The full eruption is the case where the MFR fully escapes from the Sun, leaving behind a post-flare loop arcade in a potential magnetic configuration.

In some instances, multiple solar eruptions occur within a short period of time; these are known as sympathetic eruptions. They may occur at different active regions (ARs) or within the same AR with multiple PILs. It is debatable whether these consecutive eruptions are related to each other (Biesecker & Thompson 2000). Interrelation among sympathetic eruptions has been shown in the simulation work of Török et al. (2011). In their 3D magnetohydrodynamic (MHD) simulation, the eruption of one filament can weaken the overlying magnetic pressure sufficiently, leading to the eruption of the other filaments. There is a high chance of interaction among CMEs when they occur in close proximity and/or within a short period of time. The study by Gopalswamy et al. (2002) found that interaction among CMEs could be an important factor for solar energetic particle (SEP) production. For a complex AR with multiple PILs, MFRs could be formed along the different PILs and they could erupt consecutively (e.g., Chintzoglou et al. 2015). There are also cases where, instead of having a single MFR along the PIL, two MFRs could be formed along the same PIL (e.g., Liu et al. 2012). The two MFRs arrange themselves in a double-decker (DD) configuration: one lying above the other along the same PIL. Unlike a partial eruption, in a DD configuration the lower and upper structures are observed to exist well before the commencement of the eruption. The two branches of the DD could interact with each other (e.g., Liu et al. 2012; Zhu et al. 2015).

Besides sympathetic eruptions, another type of inter-related eruption is known as a homologous event. CMEs originating from the same region, associated with homologous flares, having similar coronographic appearance and similar dimming in extreme ultraviolet imaging are known as homologous CMEs (Zhang & Wang 2002). Their study may also provide information on the physical mechanism of free-energy build-up and the triggering mechanisms of CMEs (Nitta & Hudson 2001; Zhang & Wang 2002; Vemareddy 2017).

In this paper, we analyze solar eruption SOL2012-03-10T16:50:00L303C17, which is associated with a single flare (M8.5) and a single CME. However, the flare contains two peaks in its IP separated by only 12 minutes, and the eruption contains two different erupting structures sharing the same PIL. We name this type of eruption a "compound eruption," i.e., two closely connected magnetic structures erupting consecutively within a short time interval, forming a single flare and a single CME. The compound eruption studied in this paper originated from NOAA AR 11429 on 2012 March 10. Study of the photospheric magnetic structure showed that the shearing and converging motions led to flux cancellation and new connectivity in the corona. These processes resulted in the formation of two MFR candidates, which were arranged in a DD configuration. The paper is structured as follows. The instruments and data are described in Section 2. The compound eruption is described in Section 3. The evolution of the AR prior to the eruption and the pre-eruptive configuration are analyzed in Section 4, and the discussion and conclusion are presented in Section 5.

2. Instruments

The data used in this study are primarily from the Atmospheric Imaging Assembly (AIA; Lemen et al. 2012) and Helioseismic and Magnetic Imager (HMI; Schou et al. 2012), onboard the Solar Dynamics Observatory (SDO; Pesnell et al. 2012). The AIA provides full disk images of the Sun with high cadence (12 s) and high spatial resolution 1farcs2. The 10 passbands of the AIA, covering the temperature range from 60,000 K to 20 MK, observe the Sun from the photosphere, chromosphere, and corona up to 1.3 R. The HMI provides line-of-sight magnetic field observations at a cadence of 45 s and full vector magnetic data at a cadence of 720 s. We used observations of the AR 11429 from 2012 March 9, 04:00 UT to 2012 March 10, 19:00 UT. There were brief moments when the Earth obscured (SDO eclipse) the view of the SDO on each day, therefore data were absent around that time period. We also used imaging from the Extreme Ultraviolet Imager (EUVI; Wuelser et al. 2004) and white light images of the Sun–Earth Connection Coronal and Heliospheric Investigation (SECCHI; Howard et al. 2008) onboard the Solar Terrestrial Relations Observatory (STEREO; Kaiser et al. 2008). The STEREO mission has two satellites; one moves ahead (STEREO-A) and the other behind (STEREO-B) the Earth. During the observational period, AR 11429 was behind the limb for STEREO-B and could be seen through STEREO-A and the SDO. In addition, we used observations from the Large Angle Spectrometric Coronagraph (LASCO; Brueckner et al. 1995) on board the Solar and Heliospheric Observatory (SOHO), and the Solar X-ray Imager (SXI; Hill et al. 2005 and Pizzo et al. 2005) on board the Geostationary Operational Environmental Satellite (GOES).

3. Observations of the Compound Eruption

NOAA AR 11429 was a mature AR when the eruption under study occurred. It was located near the center of the solar disk (N17E01) and had a very complex magnetic configuration (β/γ/δ). The compound eruption began on 2012 March 10, around 16:54 UT, in the southwest (SW) part of the AR. The X-ray profile of the eruption had three distinct peaks; all three were also shown in the AIA pass-bands (see Figure 1). The flare could be divided into three phases: precursor phase, IP, and gradual phase. However, unlike typical flares (of which time profiles are single peaked), this flare exhibited two peaks during the IP (indicated by two black arrows in Figure 1). AR 11429 had two pre-eruptive hot-channel structures and corresponding filaments on the SW side of the AR. These structures are called hot-channel structures as they are observed in hot passbands of the AIA i.e., 94 Å (6 MK) and 131 Å (10 MK) (e.g., Zhang et al. 2012). These two hot-channel structures (and associated filaments) were lying one above the other along the same PIL in a DD configuration (a detailed study of the pre-eruptive configuration is in Section 4.3). Hereafter, the higher-lying hot-channel structure will be referred to as "hot-channel structure A" (HCSA), and the lower-lying one as "hot-channel structure B" (HCSB). We also refer to the higher-lying filament as "filament-A" and lower-lying filament as "filament-B."

Figure 1.

Figure 1. Three phases of the compound eruption. The black curve shows the SXR time-profile observed by GOES. The colored lines show the intensity plot for different AIA passbands. The vertical-dashed lines indicate the beginning and ending of the impulsive phase. The two peaks in the impulsive phase are indicated by two black arrows.

Standard image High-resolution image

3.1. The Precursor Phase

The precursor phase is defined as the time period during which a small amount of energy is released just prior to the IP (main energy release phase) of the solar eruption (e.g., Zhou et al. 2016). Generally, a small peak is observed in the time profile of the SXR flux in the precursor phase, which marks the initiation or slow rise phase of the associated CME (Zhang & Dere 2006). During the precursor phase of this event, brightenings were observed in between the two filaments, starting at around 16:33 UT. After that moment the whole HCSA became bright in 131 Å (Figure 2(a); an animation accompanies Figure 2). Brightenings could have been due to the interaction between the two filaments. Evidence for such interaction was observed in a previous study of DD structures (Liu et al. 2012). The region near the left elbow of the HCSA was intensely brighter than the rest of the HCSA structure. In AIA 1600 Å there were sporadic brightenings underneath and along the HCSA, which later concentrated more toward the left elbow around 17:08 UT (Figure 2(f)). The study by Cheng & Ding (2016) of erupting sigmoidal structures found that the initial sporadic brightenings of the flare ribbons marked the slow rise phase of the erupting sigmoidal structure. Brightenings in 1600 Å during the precursor phase may indicate the instability and the slow rise of the HCSA in the DD configuration. The slow rise may in turn indicate the initiation -phase of the solar eruption; this is observed as a small peak in the SXR time profile in this event. The intense brightening at the left elbow of the HCSA might be due to the interaction of the rising HCSA with the overlying magnetic fields. Such activities could result in a weakening of the strapping force from the overlying magnetic fields.

Figure 2. Precursor phase of the compound eruption. The left four images show the change of the hot structure in AIA 131 Å during the precursor phase. The right four panels show the flare-ribbon evolution during the pre-flare phase in 1600 Å (blue contours), plotted over HMI/line-of-sight magnetograms (grayscale). The sequence starts on 2012 March 10, 16:50.09 and ends on 2012 March 10, 17:16:09. The video duration is 5 s.

(An animation of this figure is available.)

Video Standard image High-resolution image

3.2. The Impulsive Phase

The IP of a solar eruption is the phase of main energy release. It is characterized by a sharp increase in SXR flux and fast acceleration of the erupting structure (Zhang & Dere 2006). The onset of the IP for the compound eruption of 2012 March 10 started at 17:16 UT with the acceleration of the HCSA (the latter shown with green asterisks in the AIA 131 Å images in Figure 3). The onset of the IP was also observed as an intense and continuous brightening in AIA 1600 Å along the HCSA. As the IP progressed, the HCSA lifted to higher levels, accompanied by the separation of ribbons seen in 1600 Å and simultaneously to the rising and writhing of filament-A (for a better view of the evolution, see the animation showing the filament eruption). This series of events accompanying the acceleration of HCSA occurred during the first IP-peak in the SXR profile (first black arrow on the left in Figure 1).

Figure 3. Evolution of the compound eruption during the impulsive phase. The left four panels show the running difference images for AIA 131 Å and the right four panels show the flare-ribbon evolution, overplotted on the magnetogram, in AIA 1600 Å in blue. The erupting HCSA and HCSB are shown as green and red asterisks respectively. The blue line in panel (a) is the position of the slice used for the study of kinematic evolution of the two erupting structures. The sequence starts on 2012 March 10, 17:16:21 and ends on 2012 March 10, 17:48:21 for the AIA 131 Å left images. For the right images the flare-ribbon evolution starts on 2012 March 10, 17:16:17 and ends on 2012 March 10, 17:48:17. The video duration is 3 s.

(An animation of this figure is available.)

Video Standard image High-resolution image

The HCSB began rising at around 17:23 UT. Its fast acceleration followed at 17:29 UT as seen in AIA 131 Å (shown as red asterisks in the Figure 3). With the acceleration and rise motion of the HCSB, new sets of flare ribbons were observed which were accompanied by the rising filament-B (see the animation showing the filament eruption). The second peak in the SXR profile during the IP (second black arrow in Figure 1) is associated with the rise of the HCSB. The eruptions of the two hot-channel structures of the DD configuration were separated in time and thus were observed as two peaks in the IP of the SXR profile.

To quantitatively study the rising motion of the erupting structures, we chose an image slit, shown by the blue line in Figure 3(a), and made a time-stacking plot along this slit. In the stack-plot (see Figure 4(a)) we can clearly identify the eruption of two hot-channel structures of the DD configuration. The position of the outer edge of the HCSA is shown as green asterisks in the plot in Figure 4(a). The velocities obtained from these points indicate that the HCSA has undergone two acceleration phases. The first corresponds to the first peak of the IP, and the HCSA accelerated the second time when the HCSB erupted. The interaction between HCSA and HCSB seems to be responsible for the second acceleration phase of the HCSA. The HCSA had the second acceleration phase around 17:28 UT. The second acceleration phase also coincided with an increase in the normalized intensity of AIA 304 Å (Figure 4(b)). The velocity time profile of the HCSA and the intensity profile of the flux in AIA 304 Å shows that the fast acceleration of the HCSB started at 17:28 UT, 12 minutes after that of the HCSA. Afterwards the two structures appeared to be moving with the same velocity. The eruptions of the two structures were observed as a single CME in the SOHO LASCO/COR2 and the SECCHI STEREO/COR2 coronagraph images (Figure 5).

Figure 4.

Figure 4. Kinematic evolution of erupting structures. (a) The time-stacking plot obtained from the running difference of AIA 131 Å images (the position of the slice is shown in Figure 3(a)). The green (red) asterisks mark the positions of the outer edge of the erupting HCSA (HCSB) in the DD configuration. The two erupting structures are clearly resolved. (b) The black curve is the GOES SXR profile and blue curve shows the normalized AIA 304 Å intensity. The green curve shows the velocity variation of the HCSA during the impulsive phase; it has two acceleration phases. The curve in yellow is the derivative of the GOES flux. The vertical line shows the start time of the second acceleration phase of the HCSA and it corresponds to the second increase in the intensity of AIA 304 Å.

Standard image High-resolution image
Figure 5.

Figure 5. CME resulting from the compound eruption observed from three different points of view. The three panels from left to right show the COR2 image from SECCHI/STEREO-B, the C2 image from LASCO/SOHO, and the COR2 images from SECCHI/STEREO-A, respectively.

Standard image High-resolution image

4. Observations of the Pre-eruption Evolution

The AR 11429 was very active in terms of number and intensity of solar eruptions since it appeared on the eastern limb; it had produced 48 C-class, 15 M-class, and 3 X-class flares during its transit across the solar disk. This super-activity seems to be associated with its very complex magnetic configuration, δ-spot, and anti-Hale configuration (see Elmhamdi et al. 2014 for the evolution of the magnetic configuration). The following subsections describe the pre-eruption evolution and magnetic configuration before the onset of the compound eruption on 2012 March 10.

4.1. Evolution of the Photospheric Magnetic Field

To study the evolution of the magnetic configuration of AR 11429 we made use of preprocessed cutouts of the line-of-sight (LOS) HMI magnetograms. The preprocessing resolves the 180° azimuth ambiguity and remaps helioprojective images into a cylindrical equal area (CEA) projection, where each pixel has same surface area (Hoeksema et al. 2014). We focused on the time range between March 9, 04:00 UT to March 10, 19:00 UT. This large and complex AR could be divided into two sub-regions, NE and SW (shown as boxes in Figure 6(a)). Both sub-regions contained a long and sharp strong-gradient PIL (shown in green lines in Figure 6, obtained from an image gradient operation as detailed in Zhang et al. 2010). The compound eruption on March 10 occurred in the SW sub-region of the AR, whose evolution is shown in Figure 6. The PIL in the SW sub-region was located in between two very strong opposite polarities. We measured the centroid of positive (negative) flux greater (smaller) than 500 G (−500 G) in the SW sub-region. During the evolution of AR 11429, the distance between the centroid of the two opposite polarities (shown by red lines in Figures 6(a) and (g)) in SW sub-region decreased from ∼29.53 Mm to ∼25.86 Mm, indicating the convergence motion of the two poles. This convergence was accompanied by fragmentation and diffusion of the magnetic flux, which resulted in further collision and cancellation of opposite flux along the PIL (see the animation accompanying the Figure 6). Flux cancellation occurred along the entire PIL; however it was most noticeable for the chunk of positive flux on the left side of the SW sub-region (shown with a yellow circle in Figure 6(b)). The cancellation along the PIL suggests an increase of the free magnetic energy and has been associated with the formation of MFRs (van Ballegooijen & Martens 1989; Green et al. 2011). In addition, we observed shearing motions along the PIL (the direction of motion is shown with an orange arrow in Figure 6(c)), which seemed to change the semi-circular shape of the PIL to a linear-like shape. Shearing motions are also thought to increase the free magnetic energy (e.g., Cheng et al. 2014) in ARs. Thus, the flux cancellation and shearing motions suggest an enhancement of the free magnetic energy in the region and create the magnetic structure pending eruption.

Figure 6. Evolution of the photospheric magnetic field in NOAA AR 11429. Images are taken from the line-of-sight magnetograms of SDO/HMI and are shown in a cylindrical equal-area heliographic projection. Green lines are the strong-gradient polarity inversion line. Panel (a) shows the entire AR, and white and black boxes show the NE and SW sub-regions respectively. Panels (b)–(g) show the evolution of the magnetic configuration in the SW part of the AR. An animation of these panels is available starting on 9 March 2012, 03:59:45 until 2012 March 10, 17:26:45. The video duration is 31 s. The red lines in panel (a) and (g), show the distance between the centroid of two strong bipoles. The yellow circle in panel (b) is the location of prominent flux cancellation. The orange arrow in panel (c) is the direction of motion of negative flux.

(An animation of this figure is available.)

Video Standard image High-resolution image

4.2. Evolution of the Pre-eruption Structure in the Corona

During our observing period of 39 hr between March 9, 04:00 UT and March 10, 19:00 UT, many small and transient brightenings and confined flares occurred in the AR, including three confined C-class flares. This type of activity was mainly confined to the left side of the SW sub-region. There were also sporadic brightenings at the right side of the SW sub-region. Specifically, during that period (see Figure 7) we observed the brightenings of coronal structures lying along the PIL. The intermittent brightenings in the AIA 131 Å along the PIL indicated the presence of low-lying coherent coronal structures, which brightened up during the small confined flares. As the AR evolved, the length of the hot structure grew in size. We transformed the helioprojective images to a Carrington projection and estimated the length of hot-channel structure. The length was ∼45 Mm on 2012 March 9, at 17:21 UT and ∼76 Mm on 2012 March 10, at 02:31 UT (see Figures 7(b) and (e)) as revealed by confined flares. In the corona, the length of the non-potential magnetic field lines could increase due to the tether cutting mechanism during small flare events when new connectivity could form. In the SW sub-region of the AR, converging and shearing motions may have facilitated the tether cutting reconnection and thus increased the length of the hot-channel structure. Chintzoglou et al. (2015) observed a similar kind of flare activity, small transient brightenings and confined flares, during the formation of MFRs.

Figure 7.

Figure 7. Evolution of the coronal hot-structure before the eruption in NOAA AR 11429. AIA 131 Å passband data are shown in red/white, and the HMI data are shown in white (positive flux) and black (negative flux). In panels (b) and (e) the coronal hot structures are also traced with dark lines.

Standard image High-resolution image

4.3. Observations of the Double-decker Configuration

Long before the eruption on March 10, an underlying filament seemed to exist along the SW PIL. It was revealed when the hot post-flare loops, the product of an earlier eruption on March 9, at ∼03:30 UT, faded away. The existence of the filament shortly after the eruption suggests that it either did not erupt or erupted only partially. Around one hour prior to the compound eruption on March 10, we observed two filaments in the SW part of the AR, clearly seen from the SDO point of view (POV) (Figure 8). In the SDO POV, the two filaments seemed to be in very close proximity to each other; one filament (filament-A) was lying slightly off the SW PIL and the other (filament-B) was lying along it. Coronal structures above the solar disk were observed to deviate more and more from the PIL as they rose higher in the corona (Cheng & Ding 2016), suggesting that filament-A was probably lying higher in the corona than filament-B. To determine the geometry of the two filaments, we made use of the STEREO-A observations of them. For a specific point in one image, the scc_measure IDL routine (available in SolarSoftware Package; Freeland & Handy 1998) calculates the 3D LOS as observed by one vantage point (e.g., SDO) and then draws the projection of this line into second image (e.g., STEREO at an angle of 109fdg7). The projected line from the second satellite is known as the epipolar line (see Inhester 2006). Then the user can identify the location of features selected in the first image (SDO) onto the second image (STEREO) along this epipolar line. The 3D positions of the two filaments determined from such multiple viewpoint measurements show that they are indeed in a DD configuration (see Figure 8). The separation along the vertical direction between the two filament was found to be ∼12 Mm on 2012 March 10, at 16:02 UT. The study by Liu et al. (2012) observed a similar (∼13 Mm) separation between the two filaments in a DD configuration.

Figure 8. Two erupting filaments/prominences in a DD configuration. The upper filament (A) is marked with green asterisks and the lower filament (B) with red asterisks. The upper three panels show the filaments at different time in STEREO-A EUVI 171 Å. The lower three panels show the two filaments in SDO AIA 304 Å. The sequence starts on 2012 March 10, 16:31:32 and ends on 2012 March 10, 18:08:32. The video duration is 31 s.

(An animation of this figure is available.)

Video Standard image High-resolution image

Since there was no apparent flux emergence in the sub-region during the period of our observation, one can refute the idea of the emergence of lower magnetic structure below a higher-lying magnetic structure to form a DD configuration. Also, both parts of the DD were visible well before the initiation of the compound eruption, therefore the present case was different from the partial-eruption scenario where a single flux system divides into two parts during the eruption. It is possible that the flux cancellation and shearing motions played a role in the formation of the DD structures. Considering the initially stable configuration of the two structures, it is also possible that the axial currents in the lower and upper magnetic structures (or branches) flow in the same direction with respect to each other and the interaction between the filaments may be attractive in nature due to the ${\boldsymbol{J}}\,\times \,{\boldsymbol{B}}$ force as suggested by Kliem et al. (2014). Liu et al. (2012) observed plasma being transferred from the lower to the upper filament. This so-called filament–filament interaction in the DD was also reported in the study by Zhu et al. (2015). Here, we also observed intermittent brightenings between the two filaments in AIA/SDO, which may suggest a possible interaction between the two flux systems. From the POV of the SDO the separation between the filaments was so small that any plasma mass transfer could not be discerned easily. However, we speculate that there may be mass transfer between the two structures.

Distinguishing the DD structures clearly in the AIA passbands against the strong on-disk background emission was a very challenging task. However, intermittent brightenings occurring at different times show two different hot-channel structures along the SW PIL (see Figure 9). The temperature map, obtained using the differential emission measure (DEM) method (following Cheng et al. 2012; see Figures 9(a) and (b)), showed that the temperature of these structures was >7 MK. Hot-channel structures were also visible in the SXR images obtained by SXI/GOES (see Figures 9(e) and (f)). The foot-points of the two hot-channel structures were in close proximity to each other; however they were clearly different. Comparing the locations of the two hot-channels (AIA 94 Å; Figures 9(c) and (d)) with the two filaments (seen in AIA 304 Å; lower panels in Figure 8), we found that the structure marked with green asterisks (HCSA) corresponds well to filament-A and that with red asterisks (HCSB) corresponds to filament-B.

Figure 9.

Figure 9. Images illustrating the two hot-channel structures of the DD configuration in the SW region of the AR 11429. Panels (a) and (b) show the temperature map obtained using the DEM method. Panels (c) and (d) show the images in AIA 94 Å. Panels (e) and (f) show two hot coronal structures in X-ray images from GOES/SXI. The green and red asterisks show the high-lying and low-lying structures respectively.

Standard image High-resolution image

4.4. Coronal Magnetic Field Inferred from NLFFF Extrapolation

Solar eruptions are magnetic in origin so in order to better understand the associated mechanisms, it is crucial to study the magnetic configuration of the ARs (with the aid of modeling). The location of the AR 11429 on the solar disk made it a suitable candidate to study the coronal magnetic field through magnetic field extrapolation. In the solar corona the magnetic pressure is much higher than the plasma pressure and thus we can assume that the magnetic field lines satisfy the force-free field approximation (see Wiegelmann & Sakurai 2012), i.e., ${\boldsymbol{J}}\parallel {\boldsymbol{B}}$ and ${\rm{\nabla }}\times {\boldsymbol{B}}=\alpha {\boldsymbol{B}}$. Here, α is the force-free parameter and is constant along each magnetic field line; however it varies over different field lines in the case of NLFFFs. To get the NLFFF extrapolation of coronal magnetic field, we use an optimization numerical approach developed by Wheatland et al. (2000) and later extended by Wiegelmann (2004). In this iterative optimization process, a "penalty function" ${\rm{L}}\,={\int }_{V}w(x,y,z)[{{\boldsymbol{B}}}^{-2}| ({\rm{\nabla }}\times {\boldsymbol{B}})\times {\boldsymbol{B}}{| }^{2}$ + $| {\rm{\nabla }}\cdot {\boldsymbol{B}}{| }^{2}]{d}^{3}{\rm{x}}$ is minimized, where w is a weighting function (Wiegelmann 2004). For the bottom boundary of the domain, we use HMI vector magnetogram cutouts. These data were preprocessed to make them suitable for the force-free condition (Wiegelmann et al. 2006). The NLFFF extrapolation was done for the entire area containing AR 11429 at every hour from 2012 March 9, 04:00 UT to 2012 March 10, 17:00 UT.

The horizontal magnetic field in the PIL was aligned along its length, indicating the presence of a strong non-potential shear along the PILs (Figure 10(a)). Using the extrapolated 3D cubes and the PARAVIEW 3D visualization tool, we analyzed the evolution of the vector magnetic field in the corona. We found sheared and weakly twisted magnetic field lines above the PILs in the AR. In the NE sub-region, the sheared coronal magnetic field lines remained more or less unchanged throughout the period under study (short sheared lines shown in brown in Figure 10(b)). In the SW sub-region, however, there were different groups of sheared and twisted magnetic field lines. With time, a system of long, sheared, and twisted magnetic field lines above the PIL appeared in the modeled corona data, before the compound eruption on 2012 March 10 (shown in red in Figure 10(b)). The evolution of magnetic field lines was also evident from the evolution of the current density of the AR (more in the following paragraph).

Figure 10.

Figure 10. Magnetic field configuration in the AR 11429 at 15:58 UT on 2012 March 10, immediately before the eruption. The background shows the Bz component of the vector magnetogram with positive (negative) flux in white (black). Panel (a) also shows the direction of horizontal magnetic field at the photosphere in red and green. Panel (b) shows selected magnetic field lines in the AR above the photosphere obtained from the NLFFF extrapolation. Sheared magnetic field lines along the NE PIL are shown in brown. The sheared and twisted magnetic field lines, along the SW PIL, are shown in red. The white lines show the overlying magnetic field in the corona.

Standard image High-resolution image

We calculated the current density in the AR as ${\boldsymbol{J}}\,={\rm{\nabla }}\times B$/μo for the entire extrapolated domain. Along the PILs we observed very high current density. The presence of current density signifies non-potential magnetic fields. Due to Ohmic heating, regions of higher current density are hotter than the neighboring regions. These regions are observed as diffuse and sigmoidal structures in the AIA hot-channel passbands. In the NE sub-region, the high current-density region was more diffuse and remained unchanged during the evolution of the AR. In the SW sub-region, on the other hand, the shape of the high current-density region changed significantly over time (see the upper panels of Figure 11, corresponding to the change in the SW PIL in Figure 6) and also became more concentrated (see the lower panels of Figure 11). The change of the high current-density region from a more diffuse to more concentrated current distribution indicated the formation of long, sheared, and twisted magnetic field lines above the photosphere as displayed in the extrapolation results. This is also seen in Chintzoglou et al. (2015) where the 3D current density before the formation of the MFRs was initially a system of fragmented current channels and over time became a single current channel structure.

Figure 11.

Figure 11. Evolution of coronal current density (obtained from the NLFFF extrapolation) in the AR 11429. Panels (a)–(c) show the normalized-integrated current density along the line of sight in red/white and the Bz component in the background. Panels (d)–(f) show the vertical slice of the current density along the blue line (top left panel). During the evolution, the current density changed from more diffuse to more concentric around the location of the hot channel.

Standard image High-resolution image

The shearing motion and flux cancellation changed the magnetic configuration of the AR 11429. The extrapolation results showed that these changes were favorable for the formation of a long, sheared, and twisted magnetic structure at the SW PIL of the AR. There was one group of sheared and twisted magnetic field lines along the SW PIL, shown in red in Figure 10(b), corresponding to a single concentrated current density in the lower-right panel of Figure 11. The DD configuration was not reproduced in the NLFFF extrapolated result. The lack of a DD configuration in the extrapolation might be due to the fact that its lower and the upper branches formed during the temporal evolution. It is also important to note that the NLFFF extrapolation is producing static equilibria from the snapshot image (coronal magnetic field for the photospheric boundary condition) of the magnetic field, thus the dynamics are not fully captured. Nevertheless, NLFFF extrapolation helps in understanding the magnetic distribution above the photosphere and slow quasi-static evolution to first order.

5. Discussion and Conclusion

Our primary aim for the present study is to understand and explain the two peaks in the IP of the solar eruption from AR 11429 on 2012 March 10. It is a compound eruption, since it involves two closely spaced magnetic flux rope candidates which successively erupted within a short time interval (12 minutes in this event). The compound eruption was associated with an M-class flare. We carried out a detailed study of the compound eruption and the evolution of the AR for 39 hr leading to the eruption. Around an hour prior to the compound eruption, we observed two filaments lying along the same PIL on the SW sub-region of the AR.

The main result of the study is that the two peaks in the IP are caused by the eruption of two hot-channel structures (as seen in 131 Å corresponding to 10 MK plasma emission), which were likely in the form of a flux rope during or even before the eruption. Our analysis suggests that each erupting hot-channel structure was preceded by an independent coherent magnetic structure that destabilized and erupted. Here, we call such an eruption-capable pre-eruption magnetic structure a flux bundle, which is so named to emphasize its becoming coherent and independent long before (i.e., many hours or more than one day) the eruption, but without specifying whether such a structure is a flux rope or not. Thus, both sheared structures and flux ropes can be considered as flux bundles, as long as they become an entity of eruption at a later time. These two flux bundles were lying along the same PIL one above the other in a DD configuration. The instability and fast acceleration started first for the high-lying structure, then was followed by the instability and fast acceleration of the low-lying structure. The velocity of the low-lying structure was greater than that of the high-lying structure; as a result, the second structure reached the first and both interacted with each other to form a compound eruption. The eruption of the two flux bundles and the magnetic reconnection in the current sheet underneath them was responsible for the double peaks in the IP of the associated flare.

The present case is different from the classical eruptive model where a single flux-rope-like structure arises during the solar eruption and a single peak is observed in the SXR profile during the IP. A few previous studies consider the eruption of DD structures. The study by Liu et al. (2012) observed two filaments lying along the same PIL. They observed the instability and eruption of only the higher-lying filament. Zhu et al. (2015) also studied the eruption of the DD configuration. In their study the lower filament rose and merged with the upper filament and then the merged filaments erupted together. Unlike the previous studies, in this work, we observed the two components of the DD configuration that destabilized and erupted at different times and then interacted with each other to form a compound eruption.

Based on the studies mentioned above, we can summarize that there exist three different scenarios for eruptions in a DD configuration (shown in Figure 12). They are: (1) Instability and eruption of only the high-lying magnetic structure (this is the partial eruption of the DD configuration); (2) instability and rise of the low-lying magnetic structure, merging with the high-lying magnetic structure. Later on the merged structure erupts (this is the case of full eruption of the DD configuration); (3) instability and acceleration first of the high-lying magnetic structure followed by the low-lying magnetic structure. The two erupting structures interact with each other forming a compound structure (this is also the case of full eruption). Cases (1) and (2) would be similar to the eruption of a single flux rope with a single peak in the IP. Only case (3) would result in two peaks in the IP.

Figure 12.

Figure 12. Three erupting scenarios of the DD configuration. (1) The upper three panels depict the case where only the upper-branch of the DD erupts. (2) The middle three panels show the merging of the two branches of the DD, which erupt together. (3) The lower three panels show the eruption of the two branches of the DD separately and their interaction after the eruption.

Standard image High-resolution image

Kliem et al. (2014) studied two flux ropes in the DD configuration through MHD simulations. Analytically, they found that the DD configuration remained in stable equilibrium if the toroidal component of external sheared field lines, lying above two flux ropes, had a strength above a certain threshold value. The decrease in the strength of the toroidal component of the overlying sheared arcade would result in the instability of both of the flux ropes. Usually in such conditions the lower flux rope becomes unstable first and results in two cases: (a) merging of the two flux ropes and full eruption of the DD structure, and (b) only the upper flux rope erupts and the lower flux rope gets destroyed due to reconnection with the ambient flux. For the case when the toroidal component of the overlying sheared arcade is above the threshold value, an eruption of the upper flux rope is possible if transfer of magnetic flux and current occurs from the lower to the upper flux rope. This case also results in the eruption of the upper flux rope only (partial transfer of flux from bottom flux rope). However, in this case the lower one remains in place, possibly because it loses free energy by transferring it to the upper flux rope. In our case the interaction between the lower and the upper filament was observed as brightenings between them. During such interactions it is possible that plasma, magnetic flux and currents become transferred from the lower to the upper flux rope. Normally, such a transfer would have resulted in an eruption of only the upper flux rope, leaving the lower one undisturbed. However, the interaction of the erupting upper flux rope with the overlying sheared arcade might have reduced the toroidal component of the latter. It is possible that this interaction resulted in the instability of the lower flux rope, leading to a sequential eruption. This is one of the possible explanations of the observed eruption of the DD in the present case, but other explanations are also possible and should be addressed in a future study.

Finally, we would like to comment on the formation of the DD configuration. The sub-flaring events during the evolution of the AR possibly increased the length of coronal hot structures and formed flux bundles through magnetic reconnection. Wang & Zhang (2007) found that there is a low probability of having open eruption from the magnetic center of an AR due to strong overlying magnetic pressure. The SW PIL was lying between strong magnetic polarities, which implies that there was strong overlying magnetic pressure over the coronal structure above it. Due to that strong overlying magnetic pressure it was difficult for the magnetic structures to erupt. Cheng et al. (2014) observed the formation of a new flux rope below the existing one due to shearing, flux cancellation, and rotation of the leading magnetic polarity. Continuous shearing and flux cancellation in the AR may have resulted in the formation of a second flux rope that was lying low in the corona. Thus, strong overlying magnetic pressure and continuous shearing and flux cancellation below a higher lying flux bundle might have been responsible for the formation of a second lower-lying flux bundle, forming the DD configuration along the same PIL.

The authors are grateful to the anonymous referee for the careful review and valuable comments. HMI and AIA are instruments on board SDO, a mission for NASA's Living With a Star Program. S.D. is supported by GMU Presidential Scholarship Award. G.C. acknowledges support by NASA contract NNG04EA00C (SDO/AIA). J.Z. is supported by US NSF AGS-1249270 and NSF AGS-1156120.

Please wait… references are loading.
10.3847/1538-4357/aac028