Recurring Homologous Solar Eruptions in NOAA AR 11429
Abstract
We present the study of three homologous solar eruptions from NOAA active region (AR) 11429 over four days. This large and complex AR divided into two relatively simple sub-regions: northeast (NE) and southwest (SW). Recurrent eruptions occurred from the SW sub-region over different evolutionary phases, which provided a unique opportunity to isolate the physical processes responsible for solar eruptions. Persistent shearing and convergence of opposite magnetic polarities led to continuous flux cancellation along the SW polarity inversion line (PIL). A filament persistently lying along the SW-PIL was observed to survive each eruption, which suggests the partial eruption of the magnetic system. Further, following the first and second eruptions, a sigmoidal magnetic structure of similar morphology was reformed along the SW-PIL. The photospheric motion of magnetic flux continuously injected and stored the negative helicity in the partially erupted magnetic system and built up the magnetic free energy for the successive eruptions. These results suggest that the shearing motion and magnetic flux cancellation of opposite fluxes were: (1) the dominant factor, irrespective of the evolutionary phase, that contributed to the recurrent homologous eruption, and (2) the key processes of forming the erupting structure, likely a magnetic flux rope, and its long-lasting continuation results in reformation of identical erupting structure. The study also finds that similar magnetic topology could result in the magnetic reconnection at the same location, and such flares during the precursor phase would help in the eruption by decreasing the constraint of overlying magnetic field.
1 Introduction
Solar eruptions, manifested as flares and Coronal mass ejections (CMEs), are the most spectacular phenomena happening on the solar corona. Besides being magnificent, these are the main sources of the disturbances in the interplanetary space and also the space weather effects near the Earth. Flares and CMEs (often occur together) are probably the manifestation of a single energy release process, during which a tremendous amount of stored magnetic energy is released from the solar atmosphere (see, Forbes 2000; Zhang et al. 2001). The general model of solar eruptions, also known as CSHKP model (initially developed by Carmichael 1964; Sturrock 1966; Hirayama 1974; Kopp & Pneuman 1976), describes the process as the eruption of magnetic flux rope (MFR) through magnetic reconnection. An MFR is a highly sheared and twisted magnetic field structure and carries a large amount of free magnetic energy. It is kept under equilibrium by the magnetic pressure of the overlying magnetic fields and its instability makes it to rise. While rising it drags the overlying magnetic field lines, making them to come close to each other in an anti-parallel manner below it and form a current sheet. The magnetic reconnection at current sheet results in flares and the new connection removes the overlying constrain over MFR and accelerate it (see, Forbes 2000; Chen 2011 for review). As long as there is an erupting MFR, the general theory describes the evolution of solar eruptions very well. However, it does not address the initial formation of MFRs and onset of solar eruptions.
Many agree with the notion that the MFR is an essential feature for solar eruption. However, there is a debate regarding when, how and where it is formed. It can exist before (e.g., Kliem & Török 2006; Zhang et al. 2012) or form during solar eruption (e.g., Antiochos et al. 1999; Cheng et al. 2011). Its detection before the eruption is a difficult task, due to low plasma density and technological insufficiency to measure magnetic field directly in the corona. Nevertheless, its existence is suggested through 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 (e.g., Chintzoglou et al. 2015). Several mechanisms are associated with its formation, such as flux emergence (e.g. Fan & Gibson 2003; Leake et al. 2013), shearing motion between opposite magnetic fluxes (e.g., Amari et al. 2000, 2003; Jacobs et al. 2009), magnetic flux cancellation along the polarity inversion lines (PILs) (e.g., van Ballegooijen & Martens 1989; Aulanier et al. 2010; Green et al. 2011), and confined flaring (e.g., Patsourakos et al. 2013, Chintzoglou et al. 2015).
A typical solar eruption has three phases: (1) the precursor phase (PP), (2) the impulsive phase (IP), and (3) the gradual phase (GP; see, Zhou et al. 2016 for discussion); these are associated with three distinct acceleration phases of CMEs (Zhang et al. 2001). PP is important to understand the initiation and validating different models of solar eruptions. In the past, different observational signatures have been identified in precursor phase for certain individual events, such as: soft X-rays (SXR) emission (Fárník et al. 1996; Fárník & Savy 1998), hard X-rays burst (Harra et al. 2001), UV/EUV brightenings (e.g., Joshi et al. 2011; Awasthi et al. 2014), type II radio bursts (e.g., Klassen et al. 2003; Liu et al. 2007), prominence oscillation (Chen et al. 2008) and MFR oscillation (Zhou et al. 2016). However, due to weak or entirely absent observational activities during precursor phase its current understanding is very poor compared to IP and GP (Hudson 2011).
Most of the time, single and independent solar eruption happens on the Sun. However, there are cases where multiple eruptions occur in a short time-span. Such eruptions are known as sympathetic eruptions. These could occur in different active regions (ARs) or at different locations within the same complex AR. Although interrelation among these is a debatable issue (see, Biesecker & Thompson 2000), but causal relationship have been identified in many cases (e.g., Schrijver & Title 2011; Török et al. 2011; Liu et al. 2009).
Besides sympathetic eruptions, multiple eruptions could happen from the same local region consecutively. Consecutive 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). Homologous eruptions are possible only if there are mechanisms that continuously store magnetic energy in the same local region. Therefore, their studies are important to understand the long-term magnetic energy build-up in the corona. In the past, they have been studied either during magnetic flux emergence phase and were considered due to flux emergence (e.g. Nitta & Hudson 2001; Chatterjee & Fan 2013) or flux decaying phase and were considered due to shearing motion and magnetic flux cancellation (e.g., Li et al. 2010; Vemareddy 2017a). The persistent photospheric horizontal motion of magnetic structure along the PIL was also considered to produce homologous eruptions (e.g., Romano et al. 2015; Romano et al. 2018). Past studies also suggest that homologous eruptions can be triggered by similar mechanisms, such as moving magnetic features (Zhang & Wang 2002), coronal null point magnetic configuration (e.g., DeVore & Antiochos 2008) and shearing motion and magnetic reconnection (e.g., Vemareddy 2017a)
In this paper, we analyze three homologous eruptions, viz., SOL2012-03-07T01:05, SOL2012-03-09T03:22, and SOL2012-03-10T16:50. These eruptions originated from the same location in the complex AR 11429, but at different magnetic flux evolutionary phase. An MFR of similar morphology was formed, along the same PIL, before each of the three eruptions. Magnetic flux cancellation led by the shearing and converging motion and new connectivity in the corona were responsible for its formation. A confined flare was observed during the PP of each eruption, which helped in the eruption by weakening the constraint of overlying magnetic fields. The paper is structured as follows. Instrument and data are described in Section 2. The precursor phase and each eruption are described in Section 3. Evolution of the AR before each eruption and pre-eruptive coronal magnetic structures are described in Section 4, and the discussion and conclusion are presented in Section 5.
2 Instruments
This study mainly used the data 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 unprecedented temporal cadence (12 s) and spatial resolution (). The different spectral bands of AIA observe the Sun from the photosphere, chromosphere, and corona up to 0.5 above the solar limb. The HMI provides the photospheic line-of-sight and vector magnetic field observations at a cadence of 45 s and 720 s respectively. The SDO has an inclined geosynchronous orbit around the Earth. During two seasons of each year, the view of the SDO gets obscured by the Earth temporarily each day for two-three weeks (SDO eclipse; Pesnell et al. 2012). In this study, we focused on the SDO observation during the solar front disk passage of the NOAA AR 11429 from 2012 March 6, 00:00 UT to 2012 March 10, 19:00 UT. There were some data-gaps during the observational period due to SDO eclipse. In addition, we used observations from the Large Angle Spectrometric Coronagraph (LASCO; Brueckner et al. 1995) on board the Solar and Heliospheric Observatory (SOHO), white light images of the Sun-Earth Connection Coronal and Heliospheric Investigation (SECCHI; Howard et al. 2008) onboard Solar Terrestrial Relations Observatory (STEREO; Kaiser et al. 2008), H image from Big Bear Solar Observatory (BBSO; Denker et al. 1999), X-ray Telescope (XRT; Golub et al. 2007) onboard the Hinode (Kosugi et al. 2007), 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 Homologous Eruptions
The NOAA AR 11429 was initially emerged on the backside and appeared on the eastern limb of the Sun on 2012 March 4. During its front disk transition, the AR acquired very complex magnetic-configuration (), with many small magnetic polarities (see Elmhamdi et al. 2014 for more details about the evolution of magnetic configuration). However, as discussed in Chintzoglou et al. 2015 and Dhakal et al. 2018, it could be divided into two relatively simpler sub-regions namely: northeast (NE) and southwest (SW) sub-regions, more discussed in the section 4.1. Also, it had an anti-hale magnetic configuration i.e., the leading magnetic polarity (positive) was opposite of the normal hemispheric trend (negative). Both and anti-hale magnetic configuration are identifiable with severe flaring activity (e.g., Zirin 1970; Tanaka 1991; Sammis et al. 2000). Indeed, it was super-flare productive. Most of the flares were observed before it crossed the central meridian on March 9. By the end of March 9, it had produced 3 X-class, 11 M-class, and 30 C-class flares. Afterward, there were only 2 M-class flares and 7 C-class flares. The events under study occurred in the SW sub-region in the interval of 4 days. These homologous eruptions are summarized in Table 1.
Events | Initiation Timea | Locationb | Associated Flare (Peak Time)c | LASCO C2 Appearanced | LASCO Speed (km/s)e |
CME1 | 2012/03/07 01:05 | N15 E26 | X1.3 (01:14 UT) | 2012/03/07 01:30 UT | 1825 |
CME2 | 2012/03/09 03:22 | N15 W03 | M6.3 (03:53 UT) | 2012/03/09 04:26 UT | 950 |
CME3 | 2012/03/10 16:50 | N17 W24 | M8.4 (17:44 UT) | 2012/03/10 18:00 UT | 1296 |
aInitiation time of the solar eruption
bflare location
cpeak time of GOES X-ray flux
dfirst LASCO C2 appearance
elinear speed of CME in LASCO C2, obtained from https://cdaw.gsfc.nasa.gov/CME_list/UNIVERSAL/2012_03/univ2012_03.html
The first solar eruption occurred on 2012 March 7 at 01:05 UT when the AR was in the north-east quadrant of the solar disk. On March 7, within one hour, there were two solar eruptions viz., SOL2012-03-07T00:02 (shown in red dotted line in Figure 1(a)) and SOL2012-03-07T01:05 (shown in black line in Figure 1(a)) from the NE and SW sub-regions respectively (see Chintzoglou et al. 2015). Here, we study the eruption from the SW sub-region. It was associated with X1.3 class flare and its location was N15E26. Associated flare ribbons were concentrated mainly on the positive and negative fluxes of the SW sub-region (see Figure 1(d) for observed flare ribbon in AIA 1600 Å) and were observed to move away on either side of the PIL during the MFR eruption. The position and separation of the flare ribbons indicated that the MFR erupted from the SW-PIL. The eruption led to the coronal mass depletion in the region and observed as coronal dimming in different AIA passbands (e.g. Tian et al. 2012; see Figure 1(g) for observed dimming in AIA 193 Å). Later it was observed as halo CME, moving with 1825 km/s (CDAW CME catalog), in LASCO white light observation at 01:30 UT. The morphology of the CME at a single time snapshot is not clear due to contamination from the earlier CME. However, its morphology was clear in STEREO-A FOV as eastward-moving CME and in STEREO-B FOV as westward-moving CME (see upper panels of the Figure 2).

The second solar eruption started on 2012 March 9, at 03:22 UT when the AR was located around the central meridian of the solar disk. It was associated with an M6.3 class flare (see Figure 1(b)) at N15W03. The flare ribbons, similar to the ribbons of the previous eruption, evolved and concentrated on the opposite magnetic flux in the SW sub-region (see Figure1(e)). This again indicated that the MFR had erupted from the SW-PIL. A dimming similar to the previous eruption was observed in EUV observation (see Figure 1(h)) and eventually observed as halo CME in LASCO C2 at 04:26 UT. In the LASCO FOV CME moved with a linear velocity of 950 km/s. STEREO-A data was not available during that period, but like the earlier eruption it appeared as westward-moving CME in STEREO-B FOV (see middle panels of the Figure 2).
The third solar eruption started on 2012 March 10, at 16:50 UT when the AR was at the north-west quadrant. It was associated with an M8.4 class flare (see Figure 1(c)) at N17W24. The evolution and separation of flare ribbons were similar to the previous two eruptions and they were primarily concentrated on the opposite magnetic flux of the SW sub-region (see Figure1(f)). Therefore, the third MFR also erupted from the SW-PIL. Observed EUV dimming was also similar to the previous two eruptions (see Figure 1(i)). It appeared as halo CME in LASCO C2 at 18:00 UT, moving with a linear velocity of 1296 km/s. Like earlier eruptions, it was observed as eastward-moving CME in STEREO-A, and westward-moving CME in STEREO-B (see the lower panels of Figure 2) with a similar CME morphology.
All three solar eruptions were originated from the SW-PIL of the AR 11429. The direction of successive CMEs seemed to move towards north-side in STEREO FOV and this could be due to the change in magnetic flux distribution with evolution. Nevertheless, the overall morphology of the CMEs appeared similar in each recurring eruptions. Also, the X-ray profile, evolution of flare ribbons, and EUV dimming associated with these eruptions were nearly identical to be homologous eruptions. The time difference between the first and the second eruption was 50 hours, and the time difference between the second and the third eruption was 37 hours. The next section discusses the evolution of the AR that led to identical solar eruptions.

4 Evolution of AR 11429 leading to homologous eruptions
We used the preprocessed cutouts of SHARP data products (hmi.sharp_cea_720s) to study the photospheric evolution of the magnetic flux distribution in the AR. Preprocessing resolves the 180∘ azimuthal uncertainty and remaps helioprojective images into a cylindrical equal-area (CEA) projection, where each pixel has the same surface area (Hoeksema et al. 2014). Although there were several episodes of new magnetic flux emergence in the early evolutionary period, no major emergence was observed after March 8, as is evident from the time-profile of magnetic flux of the AR shown in the Figure 3 (also, see the movie accompanying Figure 4). Here, we used the CEA maps of HMI LOS-magnetogram to analyze the magnetic flux. The apparent 12-hour periodicity in the plots of magnetic flux (in Figure 3) is due to instrumental effect and related to spacecraft orbital velocity, also there is a variation of the number of pixels contributing to the low-to-moderate magnetic field values as AR moves from central meridian to the limb (see, Hoeksema et al. 2014). Solar eruptions under consideration had occurred at the different evolutionary phases of the AR (eruptions are shown by dashed-vertical lines in Figure 3). The first eruption had occurred when the total magnetic flux was increasing i.e. during the magnetic flux emergence phase. The second eruption happened at the end of the flux emergence phase, and the third eruption happened during the flux decay phase of the AR. In the following, we describe the evolution of the AR, focusing on SW sub-region, in the period of 2012 March 6 to March 11.

4.1 Evolution of the Photospheric Magnetic Field
Initial compactness and complexity of the magnetic flux distribution of the AR decreased with its evolution. It was divided into two simpler sub-regions: NE and SW based on the presence of two erupting MFRs on 2012 March 7 (Chintzoglou et al. 2015) and are shown in red and cyan box respectively in Figure 4(a). The distinction between the two sub-regions become more clear in the later evolutionary period. Long and sharp strong-gradient PILs, shown in green asterisks in the left panels of Figure 4 (obtained from an image gradient operation as discussed in Zhang et al. 2010), were present throughout the observational period. The plasma velocity maps show the persistent shearing motion between opposite magnetic fluxes along the SW-PILs (see the middle panels of Figure 4; also see the movie accompanying Figure 4 for the plasma velocity). These maps were derived using DAVE4VM method with an apodization window of 19 pixels (Schuck 2008). Shearing motions along the PILs are well known to convert potential field to sheared field arcade (e.g., van Ballegooijen & Martens 1989). The direction of horizontal magnetic fields, which were aligned along the length of PILs (see the right panels of Figure 4), indicated that there was indeed strong non-potential shear along the PILs. The SW-PIL was lying between two large areas of opposite magnetic polarity i.e., P-SW and N-SW (see Figure 4(d)). The distance between the centroid of opposite magnetic poles was decreasing with time (see the plot in black-asterisks in Figure 5). This showed that the opposite magnetic poles were continuously converging towards each other. The centroid-distance was decreased by 1 Mm between first and second solar eruptions (from March 7 to March 9) and by 2 Mm between second and third eruptions (from March 9 to March 10). Though the convergence was observed in both flux emerging and decaying phase, the rate of convergence was increased after the flux emergence had been stopped. Flux cancellation was not obvious from the plot of the magnetic flux of the entire AR and it was hard to isolate the SW sub-region completely. However, after March 6 it was possible to isolate the negative polarity of the SW sub-region. Therefore we encompassed the negative polarity with an ellipse and ellipse was varying with time to make sure that it fully encompassed the negative polarity all the time, and analyzed it. The plot in red asterisks in Figure 5 shows the time-profile of negative flux in the SW sub-region. The negative flux was decreasing gradually, indicating a continuous flux cancellation in the SW sub-region. Though the rate of convergence was increased after the second eruption on March 9, there was no significant change in the slope of the negative magnetic flux, which indicates that the flux cancellation rate was almost the same during and after the flux emergence phase.
To understand the high eruptivity of the AR 11429, we analyzed certain non-potential parameters: viz, non-neutralized current, length of sheared PIL, helicity (), and energy injection () rate (for computational details of these parameters see, Vemareddy 2015; Vasantharaju et al. 2018; Vemareddy 2019a). The upper panel of Figure 6 shows the non-neutralized currents in positive (negative) polarity in blue (red). A strong non-neutralized current was observed in the AR 11429 throughout the observational period. Also, long sheared PILs were observed throughout the evolution of the AR (shown in grey in the upper panel of Figure 6). The stress in the field lines due to plasma flows generates non-neutralized current in an AR (e.g., Török et al. 2014; Vemareddy et al. 2015; Vemareddy 2017b). A higher level of deviation ( 1.5) from the current neutrality condition and long sheared PILs suggests that the AR was very CME productive (e.g., Liu et al. 2017; Vemareddy 2019b). There was a significant decrease in the length of sheared PIL after March 9, and maybe that was the reason for the decrease in the flare activity after it. Further, continuous shearing and converging motions of opposite magnetic fluxes were injecting and storing magnetic helicity and magnetic energy in the solar corona (see the middle and lower panels of Figure 6). The rate of helicity and energy injection was almost constant after March 7. The initial high helicity and energy rates were due to major magnetic flux emergence during the early evolutionary phase of the AR. Continuous shearing and converging motion of the opposite magnetic fluxes were continuously injecting magnetic helicity and energy in the AR 11429 and were responsible for recurrent solar eruptions. The homologous eruptions can be understood as a cycle of storage and release of magnetic helicity and energy from the AR (Vemareddy, 2017a).



4.2 Pre-eruptive Structures in the Corona
While strongly influenced by magnetic evolution in the photosphere, solar eruptions occur in the corona. Therefore, coronal magnetic structures and their evolution are one of the most important aspects of solar eruptions. One way to analyze them is through the observation of plasma-emission at different electromagnetic wavelengths. Another way to know about the coronal magnetic field is through the extrapolation of photospheric magnetic field data. However it is not possible to get the exact 3-D magnetic field in the corona through 2-D magnetic flux observations. Therefore observational validation is necessary for the analysis of such results. Here, we used both techniques to complement each other and to have a proper understanding of pre-eruptive coronal structures and environment. To reconstruct the coronal magnetic structure, we carried out NLFFF extrapolation (Wiegelmann & Inhester 2010) at every hour from 2012 March 6, 00:10 UT to 2012 March 10, 19:34 UT. Cutouts of HMI vector magnetograms were taken as the bottom boundary of the computational domain. These were inserted in an extended field of view to weaken the effects of lateral boundaries and were pre-processed to make them suitable for the force-free conditions (Wiegelmann et al. 2006). Computations were performed on a uniformly spaced grid of pixels corresponding to a physical volume of Mm3 (more details in Dhakal et al. 2018).
A filament was observed to lie along the SW-PIL (see Figure 7), indicating sheared and/or twisted field lines (e.g., Priest et al. 1989). Following the eruptions on March 7 and 9, the filament was observed to still exist along the SW-PIL under the fading post-flare arcade in AIA cool passbands (e.g., 304 Å). The appearance of the filament, immediately after the eruption, suggests that either only part of it erupted or it did not erupt at all, thus indicating the partial eruption of the magnetic structure (e.g., Tang 1986; Gilbert et al. 2000; Gibson et al. 2002).
On March 6 at 23:35 UT, one hour prior to the first solar eruption, we observed a coherent hot-channel structure (HCS) to lie along the SW-PIL (see Figure 8(a)). Hot channels are so named 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). The temperature map, obtained through differential emission measure (DEM) method (e.g., Cheng et al. 2012), indicated that it was hotter than 7 MK. It was also observed as sigmoid in X-ray images (see Figure 8(c)). These are called sigmoid because of their forward-S or reverse-S shape (e.g., Rust & Kumar 1996). Both the HCS and X-ray sigmoids are considered as the observational proxy of MFR (e.g., Green et al. 2007; Zhang et al. 2012). NLFFF extrapolation results showed that long, sheared and twisted magnetic field lines were lying along the SW-PIL (see the top panels of Figure 9). Both observational and extrapolation results suggest the existence of an MFR in the SW sub-region of the AR. As discussed in section 3, the HCS was observed to erupt on March 7 around 01:05 UT. A coronal structure was seen to reform by March 9 01:30 UT, in the SW sub-region over the period of two days. Its morphology in EUV and X-ray images was similar to the coronal structure observed on March 6 (see the middle panels of Figure 8). Ensuing the eruption on March 9 03:00 UT, another coronal structure of identical morphology (see the lower panels of Figure 8) was observed to reform in the period of one and half days. The extrapolation results showed that the coronal structure (on both days) was sheared and twisted like the observed structure on March 6 (see Figure 9). This suggests that MFR of identical shape was reformed before each homologous eruption.
A persistent convergence and flux cancellation were observed in the SW sub-region. Also, there was a continuous shearing motion between the opposite magnetic fluxes. Shearing motion and flux cancellation are well known to generate long and twisted field lines along the PILs (e.g., van Ballegooijen & Martens 1989). This suggests that the shearing motion and flux cancellation were responsible for the recurrent formation of MFR along the SW-PIL.





4.3 Initiation of the eruption
A confined flare was observed on the left side of the SW sub-region during the PP of each eruption. It was observed as a small bump in the SXR intensity profile, a brightening of HCS, and the appearance of flare ribbons (see Figure 10 for the flare ribbons). The SXR peak of the PP for the first eruption was not observed due to the contamination from another nearby event (as discussed in 3; see Figure 1), nevertheless other signatures were observed. Identical PP could suggest that the triggering mechanisms were the same for the three solar eruptions (e.g., Vemareddy 2017a). The confined flare was observed near the small chunk of positive magnetic flux (positive-chunk), which was almost surrounded by the negative magnetic flux (shown by a red arrow in Figure 11 (a)). Using NLFFF extrapolation results and codes developed by R. Liu and J. Chen (See Liu et al. 2016), we calculated the squashing factor Q in the AR 11429. The contours of high-Q are considered as the proxy of quasi-separatrix layers (QSLs; Titov et al. 2002). Flare ribbons appeared around the positive-chunk, therefore we analyzed the magnetic topology around it. The vertical-slice of Q-map in the middle of it revealed the presence of semicircular QSLs (see Figure 11), extending from one end to the other end of it. The previous study by Polito et al. (2017) of the AR 11429, on March 9, had found 3D dome-shaped QSLs surrounding the positive-chunk. Here, the vertical slice of dome-shaped QSLs is seen as semicircular QSLs. The semicircular shape of QSLs was maintained before each eruption, though its size changed over the period. QSLs are the probable sites of magnetic reconnection (Demoulin et al. 1996, 1997). The whole positive-chunk was very dynamic throughout the evolution of the AR. Also, new flux emergence was observed there before March 9. Therefore, the shearing motion or flux emergence could have triggered the magnetic reconnection there.

The MFR, along SW-PIL, was situated between two big magnetic poles and thus a strong overlying magnetic strapping force can be inferred above it. Strong overlying fields can inhibit successful eruption from ARs (e.g., Sun et al. 2015). The rate at which the strength of overlying magnetic fields decreases with height can be analyzed using the decay-index (n). An MFR becomes torus-unstable if n is greater or equal to the critical value ( = 1.5; Bateman 1978; Kliem & Török 2006). Using NLFFF results we calculated the decay-index at different height as,
(1) |
Here, is the horizontal component of the magnetic field, and z is the radial distance from the solar surface. Figure 12 shows the distribution of super-critical decay-index (n 1.5) at different heights before each eruptions. The closed contours represent the area where horizontal magnetic fields are very weak. An erupting MFR would preferably pass through such contours (e.g., Chintzoglou et al. 2015). In the SW sub-region, the contours changed and moved more toward the positive polarity with evolution. This could be the reason for the observed change in the direction of successive CMEs, as discussed in section 3. Also, it is clear from the decay-index maps that it would have been hard for low-lying MFR to become torus-unstable and erupt from the SW-PIL. Therefore, it is most probable that the PP flare facilitated the eruption of underlying MFR by weakening the constraint of overlying magnetic field.
5 Discussion and Conclusion
Our primary goal of the present study is to understand the recurrent homologous solar eruptions from the AR 11429. There were three homologous eruptions from the southwest (SW) sub-region within 87 hours. The first eruption was associated with X-class flare and the other two eruptions were associated with M-class flares. Nevertheless, the shape of X-ray profiles of the three flares was very similar. The same photospheric location, similar evolution and location of flare ribbons, EUV dimmings, and identical coronographic morphology of the CME invoked them as homologous eruptions.
Our main result of the study is that the continuous flux cancellation due to shearing motions and the convergence of opposite magnetic bipoles along the PILs was the primary factor for the homologous solar eruptions. Flux emergence was observed only in the initial evolutionary period of the AR, flux cancellation was observed throughout the evolution that span all three events, including prior to the first event. While flux emergence was important and relevant in AR 11429, our study suggests that long-term flux cancellation and shearing motions were the most probable mechanisms that produced homologous eruptions. The photospheric motions of magnetic flux injected magnetic helicity and magnetic energy in the corona gradually. Stored helicity and energy released intermittently during solar eruptions. We believe that due to continuous shearing and converging motions of the magnetic flux, there was a cycle of storage and release of magnetic helicity and magnetic energy from the AR. In general, opposite magnetic poles of an emerging bipoles (conjugate pair) move away from each other, the separation distance depends on the flux content of the emerged bipoles (Wang & Sheeley 1989). In the present case, we did not observe the separation between opposite poles in SW sub-region, instead they were converging towards each other throughout the observational period. Therefore, we believe that these were non-conjugate pair. It was accompanied by the persistent shearing motion and magnetic flux cancellation. Chintzoglou et al. (2019) suggests that the convergence of opposite magnetic poles of non-conjugated pair leads to the magnetic flux cancellation and shearing motions, which could produce solar eruptions.
A filament was lying along the SW-PIL throughout the observational period and was observed soon after the eruptions. A filament survives and would be visible after a partial eruption. Partial eruption could happen either by the splitting of a single flux rope during the eruption(e.g., Gibson & Fan 2006; Zhang et al. 2015) or by the eruption of upper part of double-decker system (e.g., Liu et al. 2012; Kliem et al. 2014). Though a double-decker system was observed along the SW-PIL before the eruption on March 10 (see Dhakal et al. 2018), it is hard to say whether a single flux rope split or a double-decker system was present before the eruption on March 7 and March 9. Nevertheless, we believe that there was a partial eruption of the magnetic system from the SW-PIL.
A magnetic structure was formed along the SW-PIL before each eruption. It was observed as a coherent hot-channel structure in AIA EUV and sigmoid in X-ray images. Reforming HCS/sigmoid was almost identical in appearance. HCS and sigmoid are considered as an observational proxy of magnetic flux rope (MFR). Long, sheared and twisted magnetic field lines along the SW-PIL were found in our extrapolation results. Therefore, we believe that an MFR was formed and existed before each eruption. Chatterjee & Fan 2013 studied the emergence of twisted flux rope in magnetohydrodynamic simulation. They suggest that the partial eruption and reformation of MFR can produce homologous eruptions. In their simulation, reformed MFR was appeared as X-ray sigmoid. In the present case, flux emergence was going on before the first eruption. It is possible that the erupting MFR, before the first eruption, could have emerged from the sub-photosphere. However, since there were continuous shearing motion and flux cancellation in the SW sub-region, spanning the three eruptions. Also, after partial eruption an erupting magnetic structure of similar morphology was reformed. Therefore we believe that, most probably the shearing motion and flux cancellation were the primary mechanisms that formed the erupting MFR along the SW-PIL.
At last, we want to discuss about the identical precursor phase (PP) activity in the AR. It was identified as a confined flare surrounding an island of positive flux at the left side of the SW sub-region. A complex magnetic configuration, semicircular quasi-separatrix layers (QSLs), was inferred to exist in the region. The magnetic configuration was almost the same before each eruption, we believe this was due to the reformation of identical magnetic structures in the region. It is arguable that the identical PP was due to reconnection at identically shaped QSLs, either due to flux emergence or shifting of the island itself. As a strong overlying magnetic strapping force, due to overlying magnetic field, was inferred above the MFR lying along the SW-PIL. Our study suggests that the PP flare facilitated the eruption by removing the overlying magnetic strapping force above the MFR.
References
- Amari et al. (2003) Amari, T., Luciani, J. F., Aly, J. J., Mikic, Z., & Linker, J. 2003, ApJ, 585, 1073, doi: 10.1086/345501
- Amari et al. (2000) Amari, T., Luciani, J. F., Mikic, Z., & Linker, J. 2000, ApJ, 529, L49, doi: 10.1086/312444
- Antiochos et al. (1999) Antiochos, S. K., DeVore, C. R., & Klimchuk, J. A. 1999, ApJ, 510, 485, doi: 10.1086/306563
- Aulanier et al. (2010) Aulanier, G., Török, T., Démoulin, P., & DeLuca, E. E. 2010, ApJ, 708, 314, doi: 10.1088/0004-637X/708/1/314
- Awasthi et al. (2014) Awasthi, A. K., Jain, R., Gadhiya, P. D., et al. 2014, MNRAS, 437, 2249, doi: 10.1093/mnras/stt2032
- Bateman (1978) Bateman, G. 1978, MHD instabilities
- Biesecker & Thompson (2000) Biesecker, D. A., & Thompson, B. J. 2000, Journal of Atmospheric and Solar-Terrestrial Physics, 62, 1449, doi: 10.1016/S1364-6826(00)00085-7
- Brueckner et al. (1995) Brueckner, G. E., Howard, R. A., Koomen, M. J., et al. 1995, Sol. Phys., 162, 357, doi: 10.1007/BF00733434
- Carmichael (1964) Carmichael, H. 1964, NASA Special Publication, 50, 451
- Chatterjee & Fan (2013) Chatterjee, P., & Fan, Y. 2013, ApJ, 778, L8, doi: 10.1088/2041-8205/778/1/L8
- Chen (2011) Chen, P. F. 2011, Living Reviews in Solar Physics, 8, 1, doi: 10.12942/lrsp-2011-1
- Chen et al. (2008) Chen, P. F., Innes, D. E., & Solanki, S. K. 2008, A&A, 484, 487, doi: 10.1051/0004-6361:200809544
- Cheng et al. (2011) Cheng, X., Zhang, J., Liu, Y., & Ding, M. D. 2011, ApJ, 732, L25, doi: 10.1088/2041-8205/732/2/L25
- Cheng et al. (2012) Cheng, X., Zhang, J., Saar, S. H., & Ding, M. D. 2012, ApJ, 761, 62, doi: 10.1088/0004-637X/761/1/62
- Chintzoglou et al. (2015) Chintzoglou, G., Patsourakos, S., & Vourlidas, A. 2015, ApJ, 809, 34, doi: 10.1088/0004-637X/809/1/34
- Chintzoglou et al. (2019) Chintzoglou, G., Zhang, J., Cheung, M. C. M., & Kazachenko, M. 2019, ApJ, 871, 67, doi: 10.3847/1538-4357/aaef30
- Demoulin et al. (1997) Demoulin, P., Bagala, L. G., Mandrini, C. H., Henoux, J. C., & Rovira, M. G. 1997, A&A, 325, 305
- Demoulin et al. (1996) Demoulin, P., Henoux, J. C., Priest, E. R., & Mandrini, C. H. 1996, A&A, 308, 643
- Denker et al. (1999) Denker, C., Johannesson, A., Marquette, W., et al. 1999, Sol. Phys., 184, 87, doi: 10.1023/A:1005047906097
- DeVore & Antiochos (2008) DeVore, C. R., & Antiochos, S. K. 2008, ApJ, 680, 740, doi: 10.1086/588011
- Dhakal et al. (2018) Dhakal, S. K., Chintzoglou, G., & Zhang, J. 2018, ApJ, 860, 35, doi: 10.3847/1538-4357/aac028
- Elmhamdi et al. (2014) Elmhamdi, A., Romano, P., Kordi, A. S., & Al-trabulsy, H. A. 2014, Sol. Phys., 289, 2957, doi: 10.1007/s11207-014-0507-9
- Fan & Gibson (2003) Fan, Y., & Gibson, S. E. 2003, ApJ, 589, L105, doi: 10.1086/375834
- Fárník et al. (1996) Fárník, F., Hudson, H., & Watanabe, T. 1996, Sol. Phys., 165, 169, doi: 10.1007/BF00149096
- Fárník & Savy (1998) Fárník, F., & Savy, S. K. 1998, Sol. Phys., 183, 339, doi: 10.1023/A:1005092927592
- Forbes (2000) Forbes, T. G. 2000, J. Geophys. Res., 105, 23153, doi: 10.1029/2000JA000005
- Gibson & Fan (2006) Gibson, S. E., & Fan, Y. 2006, ApJ, 637, L65, doi: 10.1086/500452
- Gibson et al. (2002) Gibson, S. E., Fletcher, L., Del Zanna, G., et al. 2002, ApJ, 574, 1021, doi: 10.1086/341090
- Gilbert et al. (2000) Gilbert, H. R., Holzer, T. E., Burkepile, J. T., & Hundhausen, A. J. 2000, ApJ, 537, 503, doi: 10.1086/309030
- Golub et al. (2007) Golub, L., Deluca, E., Austin, G., et al. 2007, Sol. Phys., 243, 63, doi: 10.1007/s11207-007-0182-1
- Green et al. (2007) Green, L. M., Kliem, B., Török, T., van Driel-Gesztelyi, L., & Attrill, G. D. R. 2007, Sol. Phys., 246, 365, doi: 10.1007/s11207-007-9061-z
- Green et al. (2011) Green, L. M., Kliem, B., & Wallace, A. J. 2011, A&A, 526, A2, doi: 10.1051/0004-6361/201015146
- Harra et al. (2001) Harra, L. K., Matthews, S. A., & Culhane, J. L. 2001, ApJ, 549, L245, doi: 10.1086/319163
- Hill et al. (2005) Hill, S. M., Pizzo, V. J., Balch, C. C., et al. 2005, Sol. Phys., 226, 255, doi: 10.1007/s11207-005-7416-x
- Hirayama (1974) Hirayama, T. 1974, Sol. Phys., 34, 323, doi: 10.1007/BF00153671
- Hoeksema et al. (2014) Hoeksema, J. T., Liu, Y., Hayashi, K., et al. 2014, Sol. Phys., 289, 3483, doi: 10.1007/s11207-014-0516-8
- Howard et al. (2008) Howard, R. A., Moses, J. D., Vourlidas, A., et al. 2008, Space Sci. Rev., 136, 67, doi: 10.1007/s11214-008-9341-4
- Hudson (2011) Hudson, H. S. 2011, Space Sci. Rev., 158, 5, doi: 10.1007/s11214-010-9721-4
- Jacobs et al. (2009) Jacobs, C., Roussev, I. I., Lugaz, N., & Poedts, S. 2009, ApJ, 695, L171, doi: 10.1088/0004-637X/695/2/L171
- Joshi et al. (2011) Joshi, B., Veronig, A. M., Lee, J., et al. 2011, ApJ, 743, 195, doi: 10.1088/0004-637X/743/2/195
- Kaiser et al. (2008) Kaiser, M. L., Kucera, T. A., Davila, J. M., et al. 2008, Space Sci. Rev., 136, 5, doi: 10.1007/s11214-007-9277-0
- Klassen et al. (2003) Klassen, A., Pohjolainen, S., & Klein, K.-L. 2003, Sol. Phys., 218, 197, doi: 10.1023/B:SOLA.0000013034.61996.c4
- Kliem & Török (2006) Kliem, B., & Török, T. 2006, Physical Review Letters, 96, 255002, doi: 10.1103/PhysRevLett.96.255002
- Kliem et al. (2014) Kliem, B., Török, T., Titov, V. S., et al. 2014, ApJ, 792, 107, doi: 10.1088/0004-637X/792/2/107
- Kopp & Pneuman (1976) Kopp, R. A., & Pneuman, G. W. 1976, Sol. Phys., 50, 85, doi: 10.1007/BF00206193
- Kosugi et al. (2007) Kosugi, T., Matsuzaki, K., Sakao, T., et al. 2007, Sol. Phys., 243, 3, doi: 10.1007/s11207-007-9014-6
- Leake et al. (2013) Leake, J. E., Linton, M. G., & Török, T. 2013, ApJ, 778, 99, doi: 10.1088/0004-637X/778/2/99
- Lemen et al. (2012) Lemen, J. R., Title, A. M., Akin, D. J., et al. 2012, Sol. Phys., 275, 17, doi: 10.1007/s11207-011-9776-8
- Li et al. (2010) Li, Y., Lynch, B. J., Welsch, B. T., et al. 2010, Sol. Phys., 264, 149, doi: 10.1007/s11207-010-9547-y
- Lites (2005) Lites, B. W. 2005, ApJ, 622, 1275, doi: 10.1086/428080
- Liu et al. (2009) Liu, C., Lee, J., Karlický, M., et al. 2009, ApJ, 703, 757, doi: 10.1088/0004-637X/703/1/757
- Liu et al. (2007) Liu, C., Lee, J., Yurchyshyn, V., et al. 2007, ApJ, 669, 1372, doi: 10.1086/521644
- Liu et al. (2012) Liu, R., Kliem, B., Török, T., et al. 2012, ApJ, 756, 59, doi: 10.1088/0004-637X/756/1/59
- Liu et al. (2016) Liu, R., Kliem, B., Titov, V. S., et al. 2016, ApJ, 818, 148, doi: 10.3847/0004-637X/818/2/148
- Liu et al. (2017) Liu, Y., Sun, X., Török, T., Titov, V. S., & Leake, J. E. 2017, ApJ, 846, L6, doi: 10.3847/2041-8213/aa861e
- López Ariste et al. (2006) López Ariste, A., Aulanier, G., Schmieder, B., & Sainz Dalda, A. 2006, A&A, 456, 725, doi: 10.1051/0004-6361:20064923
- Nitta & Hudson (2001) Nitta, N. V., & Hudson, H. S. 2001, Geophys. Res. Lett., 28, 3801, doi: 10.1029/2001GL013261
- Patsourakos et al. (2013) Patsourakos, S., Vourlidas, A., & Stenborg, G. 2013, ApJ, 764, 125, doi: 10.1088/0004-637X/764/2/125
- Pesnell et al. (2012) Pesnell, W. D., Thompson, B. J., & Chamberlin, P. C. 2012, Sol. Phys., 275, 3, doi: 10.1007/s11207-011-9841-3
- Pizzo et al. (2005) Pizzo, V. J., Hill, S. M., Balch, C. C., et al. 2005, Sol. Phys., 226, 283, doi: 10.1007/s11207-005-7417-9
- Polito et al. (2017) Polito, V., Del Zanna, G., Valori, G., et al. 2017, A&A, 601, A39, doi: 10.1051/0004-6361/201629703
- Priest et al. (1989) Priest, E. R., Hood, A. W., & Anzer, U. 1989, ApJ, 344, 1010, doi: 10.1086/167868
- Romano et al. (2018) Romano, P., Elmhamdi, A., Falco, M., et al. 2018, ApJ, 852, L10, doi: 10.3847/2041-8213/aaa1df
- Romano et al. (2015) Romano, P., Zuccarello, F., Guglielmino, S. L., et al. 2015, A&A, 582, A55, doi: 10.1051/0004-6361/201525887
- Rust & Kumar (1996) Rust, D. M., & Kumar, A. 1996, ApJ, 464, L199, doi: 10.1086/310118
- Sammis et al. (2000) Sammis, I., Tang, F., & Zirin, H. 2000, ApJ, 540, 583, doi: 10.1086/309303
- Schou et al. (2012) Schou, J., Scherrer, P. H., Bush, R. I., et al. 2012, Sol. Phys., 275, 229, doi: 10.1007/s11207-011-9842-2
- Schrijver & Title (2011) Schrijver, C. J., & Title, A. M. 2011, Journal of Geophysical Research (Space Physics), 116, A04108, doi: 10.1029/2010JA016224
- Schuck (2008) Schuck, P. W. 2008, ApJ, 683, 1134, doi: 10.1086/589434
- Sturrock (1966) Sturrock, P. A. 1966, Nature, 211, 695, doi: 10.1038/211695a0
- Sun et al. (2015) Sun, X., Bobra, M. G., Hoeksema, J. T., et al. 2015, ApJ, 804, L28, doi: 10.1088/2041-8205/804/2/L28
- Tanaka (1991) Tanaka, K. 1991, Sol. Phys., 136, 133, doi: 10.1007/BF00151700
- Tang (1986) Tang, F. 1986, Sol. Phys., 105, 399, doi: 10.1007/BF00172056
- Tian et al. (2012) Tian, H., McIntosh, S. W., Xia, L., He, J., & Wang, X. 2012, ApJ, 748, 106, doi: 10.1088/0004-637X/748/2/106
- Titov et al. (2002) Titov, V. S., Hornig, G., & Démoulin, P. 2002, Journal of Geophysical Research (Space Physics), 107, 1164, doi: 10.1029/2001JA000278
- Török et al. (2011) Török, T., Panasenco, O., Titov, V. S., et al. 2011, ApJ, 739, L63, doi: 10.1088/2041-8205/739/2/L63
- Török et al. (2014) Török, T., Leake, J. E., Titov, V. S., et al. 2014, ApJ, 782, L10, doi: 10.1088/2041-8205/782/1/L10
- van Ballegooijen & Martens (1989) van Ballegooijen, A. A., & Martens, P. C. H. 1989, ApJ, 343, 971, doi: 10.1086/167766
- Vasantharaju et al. (2018) Vasantharaju, N., Vemareddy, P., Ravindra, B., & Doddamani, V. H. 2018, ApJ, 860, 58, doi: 10.3847/1538-4357/aac272
- Vemareddy (2015) Vemareddy, P. 2015, ApJ, 806, 245, doi: 10.1088/0004-637X/806/2/245
- Vemareddy (2017a) —. 2017a, ApJ, 845, 59, doi: 10.3847/1538-4357/aa7ff4
- Vemareddy (2017b) —. 2017b, ApJ, 851, 3, doi: 10.3847/1538-4357/aa9706
- Vemareddy (2019a) —. 2019a, ApJ, 872, 182, doi: 10.3847/1538-4357/ab0200
- Vemareddy (2019b) —. 2019b, MNRAS, 486, 4936, doi: 10.1093/mnras/stz1020
- Vemareddy et al. (2015) Vemareddy, P., Venkatakrishnan, P., & Karthikreddy, S. 2015, Research in Astronomy and Astrophysics, 15, 1547, doi: 10.1088/1674-4527/15/9/011
- Wang & Sheeley (1989) Wang, Y. M., & Sheeley, N. R., J. 1989, Sol. Phys., 124, 81, doi: 10.1007/BF00146521
- Wiegelmann & Inhester (2010) Wiegelmann, T., & Inhester, B. 2010, A&A, 516, A107, doi: 10.1051/0004-6361/201014391
- Wiegelmann et al. (2006) Wiegelmann, T., Inhester, B., & Sakurai, T. 2006, Sol. Phys., 233, 215, doi: 10.1007/s11207-006-2092-z
- Zhang et al. (2012) Zhang, J., Cheng, X., & Ding, M.-D. 2012, Nature Communications, 3, 747, doi: 10.1038/ncomms1753
- Zhang et al. (2001) Zhang, J., Dere, K. P., Howard, R. A., Kundu, M. R., & White, S. M. 2001, ApJ, 559, 452, doi: 10.1086/322405
- Zhang & Wang (2002) Zhang, J., & Wang, J. 2002, ApJ, 566, L117, doi: 10.1086/339660
- Zhang et al. (2010) Zhang, J., Wang, Y., & Liu, Y. 2010, ApJ, 723, 1006, doi: 10.1088/0004-637X/723/2/1006
- Zhang et al. (2015) Zhang, Q. M., Ning, Z. J., Guo, Y., et al. 2015, ApJ, 805, 4, doi: 10.1088/0004-637X/805/1/4
- Zhou et al. (2016) Zhou, G. P., Zhang, J., & Wang, J. X. 2016, ApJ, 823, L19, doi: 10.3847/2041-8205/823/1/L19
- Zirin (1970) Zirin, H. 1970, Sol. Phys., 14, 328, doi: 10.1007/BF00221318