On the Characteristics of Footpoints of Solar Magnetic Flux Ropes during the Eruption
We investigate the footpoints of four erupted magnetic flux ropes (MFRs) that appear as sigmoidal hot channels prior to the eruptions in the Atmospheric Imaging Assembly high temperaure passbands. The simultaneous Helioseismic and Magnetic Imager observations disclose that one footpoint of the MFRs originates in the penumbra or penumbra edge with a stronger magnetic field, while the other in the moss region with a weaker magnetic field. The significant deviation of the axis of the MFRs from the main polarity inversion lines and associated filaments suggests that the MFRs have ascended to a high altitude, thus being distinguishable from the source sigmoidal ARs. The more interesting thing is that, with the eruption of the MFRs, the average inclination angle and direct current at the footpoints with stronger magnetic field tend to decrease, which is suggestive of a straightening and untwisting of the magnetic field in the MFR legs. Moreover, the associated flare ribbons also display an interesting evolution. They initially appear as sporadical brightenings at the two footpoints of and in the regions below the MFRs and then quickly extend to two slender sheared J-shaped ribbons with the two hooks corresponding to the two ends of the MFRs. Finally, the straight parts of the two ribbons separate from each other, evolving into two widened parallel ones. These features mostly conforms to and supports the recently proposed three-dimensional standard CME/flare model, i.e., the twisted MFR eruption stretches and leads to the reconnection of the overlying field that transits from a strong to weak shear with the increasing height.
Subject headings:Sun: corona — Sun: coronal mass ejections (CMEs) — Sun: magnetic fields — Sun: flares
Online-only material: animations, color figures
Coronal mass ejections (CMEs) and flares are among the most spectacular explosive phenomena in our solar system. They are able to release various relativistic particles and eject a large quantity of plasma and magnetic flux with a velocity of hundreds of km s, even up to 3000 km s (Yashiro et al., 2004), into the interplanetary space. When arriving at the Earth, these energetic particles and high-speed magnetized plasma may interact with the magnetosphere and ionosphere, and seriously affect the safety of human high-tech activities in the outer space such as disrupting communications, overloading power grids, presenting a hazard to astronauts, and so on (Webb et al., 1994; Liu et al., 2014).
In the past decades, theoretical solar physicists proposed that a magnetic flux rope (MFR), a coherent magnetic structure with the field lines wrapping around the central axis more than once, is the fundamental structure in the CME/flare dynamical process (e.g., Shibata et al., 1995; Chen, 1996; Titov & Démoulin, 1999; Chen, 2011; Vourlidas et al., 2013). Forbes & Isenberg (1991) studied a cylindrical MFR and found that such helical structure can suddenly loss its equilibrium and rapidly erupt upward. Kliem & Török (2006) considered a freely expanded current torus. With a detailed analysis, they concluded that the MFR can experience a nonlinear expansion if the background magnetic field that constrains the MFR declines fast enough. Such an expansion instability is also named as torus instability. Following the work of Kliem & Török (2006), Olmedo & Zhang (2010) further studied a more realistic situation, in which two footpoints of the MFR are line-tied to the photosphere, well resembling the pre-eruption state of the MFR. They pointed out that the threshold of torus instability not only depends on the decline of the background field but also on the geometrical circularity of the MFR: a ratio of the MFR arc length above the photosphere to its circumference. Besides torus instability, kink instability can also initiate an eruption. When the twist of the MFR exceeds a critical value, the axis of the MFR will rapidly deform, transforming part of the twist into writhe (Török et al., 2004; Fan & Gibson, 2007). Meanwhile, the height of the MFR quickly ascends, triggering the MFR eruption (e.g., Ji et al., 2003; Guo et al., 2010a; Srivastava et al., 2010; Cheng et al., 2014a; Yan et al., 2014).
In addition to ideal MHD instabilities, reconnection-associated mechanisms such as tether-cutting (Moore et al., 2001) and breakout reconnection (Antiochos et al., 1999; Chen & Shibata, 2000; Karpen et al., 2012) are also proposed to interpret the initiation of the MFR. By contrast, these two models do not require that the MFR has to exist prior to the eruption. In the tether-cutting model, the erupted magnetic field initially is comprised of two groups of sheared arcades. The reconnection between the sheared arcades will form a helical magnetic field, i.e., the nascent MFR (e.g., Liu et al., 2013; Chen et al., 2014b). An upward Lorentz force is simultaneously generated to initiate the MFR. On the contrary, the breakout model consists of a quadrupole field configuration with a null point situated at between the central arcade and the overlying large-scale field. Once the reconnection near the null point commences, the overlying field will be opened, leading to the reduction of the downward tension and thus the rise of the central arcade (e.g., Zuccarello et al., 2008; Shen et al., 2012; Lynch & Edmondson, 2013).
No matter which model applies, the flare reconnection below the erupted configuration can be switched on, which rapidly produces the poloidal flux to accelerate the eruption of the MFR. At the same time, the fast reconnection generates energetic electrons, which stream down along newly reconnected field lines and heat the chromosphere, mapping two elongated ribbon-like structures and casuing the flare emission in various electromagnetic passbands (Priest & Forbes, 2002). Different from those ideal MHD mechanisms that suppose a pre-existing MFR, the tether-cutting reconnection forms the nascent MFR in the slow rise phase, which further grows up in the acceleration phase, while the breakout reconnection does not build up the MFR until in the acceleration phase.
Owing to the important role of the MFR in the CME/flare dynamics, the solar community has been searching for the observational evidence of the MFR. In solar active regions (ARs), forward or reversed sigmoidal emission pattern in extreme ultraviolet (EUV) and soft X-ray (SXR) passbands is often seen as a progenitor of the CME (e.g., Rust & Kumar, 1996; Canfield et al., 1999; Sterling et al., 2000; Liu et al., 2007; McKenzie & Canfield, 2008; Liu et al., 2010; Savcheva et al., 2012b), thus being considered as evidence of the MFR preexisting in the corona. On the other hand, filaments correspond to a collection of relatively dense and cool plasma suspended at dips of a helical structure; thus, filaments also serve as an indicator of the MFR (Mackay et al., 2010). If there are no filament materials, the helical MFR structure probably appears as dark filament channels along the polarity inversion line (PIL) as been visible in the EUV passbands (van Ballegooijen et al., 1998; Chen et al., 2014c). When rotating to the solar limb, the cross section of filament channels will manifest as dark cavities. The frequently observed spinning motions (Low & Hundhausen, 1995; Gibson et al., 2004; Wang & Stenborg, 2010; Li et al., 2012), bright ring (Dove et al., 2011), and “lagomorphic” structure of linear polarization (Bak-Stȩślicka et al., 2013) in the cavities, are even thought to be direct evidence of the existence of helical magnetic field. Moreover, the descending motion of filament materials along a helical trajectory is also argued to be strong evidence of a helical structure (e.g., Li & Zhang, 2013b; Cheng et al., 2014c; Yang et al., 2014; Zhang et al., 2015). The intrinsic relationship between filaments and MFRs has further been supported by various case studies, in which the dips of the twisted field lines reproduced by the non-linear force free field (NLFFF) extrapolation basically conform with the filament locations (e.g., Guo et al., 2010b; Savcheva et al., 2012b; Su & van Ballegooijen, 2012; Cheng et al., 2013b; Inoue et al., 2013; Jiang et al., 2014a, b, 2016; Yan et al., 2015).
Recently, Cheng et al. (2011) presented unambiguous evidence of an MFR being formed in the lower corona (also see; Song et al., 2014a, b). It manifests as an EUV blob in the high temperature passbands of the Atmospheric Imaging Assembly (AIA; Lemen et al., 2012). In the other low temperature passbands, however, it appears as a dark cavity. Subsequently, Zhang et al. (2012) and Cheng et al. (2013a) unveiled that, if seen off the solar limb, the MFR will take on a coherent channel-like hot structure that exists prior to the eruption (also see; Patsourakos et al., 2013; Li & Zhang, 2013a, c; Tripathi et al., 2013; Chintzoglou et al., 2015; Joshi et al., 2015; Cheng et al., 2015). Once the hot channel takes off, it gradually evolves to a semicircular loop-like structure in the early rise phase. Afterwards, the loop-like structure quickly drives the formation of the CME and gives rise to the flare emission simultaneously. Furthermore, through a detailed case study, Cheng et al. (2014c) identified that the hot channel can evolve smoothly from the inner into the outer corona by retaining its coherence, whose morphology coincides with the CME cavity in the white-light images. Moreover, Cheng et al. (2014a) and Chen et al. (2014a) documented that the hot channel is initially cospatial with but later on separates from the associated prominence. At the later time, the prominence only occupies the lower part of the hot channel spatially. Nindos et al. (2015) further made a statistical study on the appearance frequency of the hot structures and found that almost half of major eruptive flares contain a hot blob or channel-like configuration. All of these studies support the arguement that the EUV hot blobs and channels are actually the helical MFRs that exist prior to the eruption.
Therefore, it can be said that the EUV hot channels can now be regarded as the most promising evidence of the MFRs. However, previous studies on the MFRs are mostly restricted to their morphology or kinematics, while no attention is paid to the properties of their footpoints. This is mainly due to that it is extremely difficult to find out appropriate hot channels events, in which the footpoints can be identified unambiguously. In the current study, by inspecting tens of solar eruption events, we find out four hot channels events which are well observed and have the discernable footpoints. Furthermore, we detailedly investigate various characteristics at the footpoints of the MFRs, as well as their relation to the flare ribbons. On the other hand, Aulanier et al. (2012) and Janvier et al. (2013) recently extended the two-dimensional (2D) standard CME/flare model (CSHKP; Carmichael, 1964; Sturrock, 1966; Hirayama, 1974; Kopp & Pneuman, 1976) to 3D, which predicts the vertical straightening of the inner legs of the CME and the transfer of the magnetic shear from strong to weak after the eruption. Thus, the revealed observational characteristics can also be used to compare and validate the 3D CME/flare model. In Section 2, we introduce the observing instruments, followed by the event selection in Section 3. In Section 4, we display the results. Finally, we give a summary and discussions in Section 5.
The data sets used in this study are primarily from Solar Dynamics Observatory (SDO; Pesnell et al., 2012). The AIA on board SDO provides images of the solar corona that are acquired almost simultaneously at temperatures ranging from 0.06 MK to 20 MK through 9 passbands (7 EUV and 2 UV passbands) with a temporal cadence of 12 s or 24 s and spatial resolution of 1.2″. The Helioseismic and Magnetic Imager (HMI; Schou et al., 2012) also on board SDO observes the photosphere at 6173 Å and provide light-of-sight magnetograms of the full disk with a temporal cadence of 45 s and spatial resolution of 1.0″. The vector magnetograms are provided by HMI Active Region Patches (HARPs; Bobra et al., 2014) product, in which the 180 ambiguity in the transverse component has been removed. Moreover, we take advantage of the Large Angle and Spectrometric Coronagraph (LASCO; Brueckner et al., 1995) on board Solar and Heliospheric Observatory (SOHO) and GOES to inspect CMEs and monitor the soft X-ray (SXR) 1–8 Å flux of the associated flares, respectively.
3. Event Selection
In order to investigate the footpoints of MFRs, we search for appropriate events from the ISEST Event List
Four events are found to fit the above criteria as shown in Table 1. One can see that the four flares have a long rise time (20 minutes) and a large magnitude (M6.0), and that the four CMEs are full halo ones with velocities of 800 km s. Figure 1 shows the long-term evolution of the positive and negative magnetic fluxes in their source ARs. We find that all eruptions appear in the decay phase of the ARs but at their early stage. Note that, for the SOL2013-04-11T event, the increase of the positive and negative fluxes is due to the emergence of a second bipolar that appears in the northwest of the main bipolar. Moreover, we can see that in a long period, at least one day, before and after the eruption, there is no other big event taking place, which shows that the magnetic environment of the four eruptions is not influenced by other major eruptions.
4.1. Morphology of the MFRs
We find that prior to the eruption, the most apparent feature appearing in the four ARs is a curved cylinder or writhed channel-like structure as seen in the AIA 131 Å and 94 Å images. From Figure 2–9, one can see that all of the channels take on a forward (SOL2012-07-12T and SOL2014-04-18T events) or reversed (SOL2013-04-11T and SOL2014-09-10T events) sigmoidal shape when approaching the eruption, at least in its early phase. Whereas in the AIA 171 Å passband, the ARs mainly display the overlying loops that straddle over and are almost perpendicular to the channels. The channels most likely contain plasmas of high temperatures since they only appear in the AIA high temperature passbands and are almost invisible in other low temperature passbands. Actually, the detailed differential emission measure analyses, which have been applied to the SOL2012-07-12T and SOL2013-04-11T events by Cheng et al. (2014b) and Vemareddy & Zhang (2014), respectively, proved that the channels do have a high temperature of 8 MK. Because of the high temperature property, Zhang et al. (2012) and Cheng et al. (2013a) named the curved cylinder or writhed channel structures as EUV hot channels.
It is worth noting that, although the ARs often take on a sigmoidal shape in the EUV and/or soft X-ray images, it does not mean that the magnetic field lines there must be continuously sigmoidal threads (Titov & Démoulin, 1999; Kliem et al., 2004; Valori et al., 2010; Schmieder et al., 2015; Jiang & Feng, 2016). Alternatively, there could be two groups of sheared arcades, making up a sigmoidal shape as a whole (e.g., Rust & Kumar, 1996; Canfield et al., 1999; McKenzie & Canfield, 2008; Liu et al., 2010). However, the latter explanation is not appropriate for the four hot channels studied here, because they are mostly comprised of continuously sigmoidal threads for the reasons stated in detail below.
As for the SOL2012-07-12T case, Cheng et al. (2014b) studied the formation of the hot channel-like MFR in detail and found that the twisted field lines have been well built up a half day before the eruption. Furthermore, they found that, two hours before the eruption, the whole sigmoid actually is a double-decker MFR structure; the low-lying one is manifested by the filament and the high-lying one is evidenced by the strongly writhed channel-like structure (Figure 2). With the flare beginning, it is clearly seen that the high-lying channel detach from the sigmoidal AR with the middle part moving toward the south. The high-lying hot channel being an MFR is also proved by Dudík et al. (2014) with the MHD simulation, in which they showed that the high-lying sigmoidal threads and their evolution are well reproduced by an unstable MFR and its eruption process.
Figure 3 presents the early evolution of the SOL2013-04-11T sigmoid (also see the attached online movie of Figure 7). One can see that a reversed sigmoid first appears in the core field of the source AR. From 06:45 UT, part of the S-shaped field lines start to detach from the sigmoid and then slowly rise up with the middle part having the fastest speed. After ten minutes, the S-shaped field lines show an M-shape, which quickly evolves into a semicircular shape. This similar evolution behavior to the SOL2011-03-08T event analyzed by Zhang et al. (2012) and Cheng et al. (2013a) strongly suggests that the erupted S-shaped threads are actually the MFR (also see, Vemareddy & Zhang, 2014).
For the SOL2014-04-18T event, the forward sigmoidal channel appears as early as one hour before the eruption. Thanks to the high resolution of the AIA, we can see that the channel structure mainly consists of a group of coherently sigmoidal threads (Figure 4). Through analyzing the IRIS spectral lines formed in the transition region and the chromosphere and the HMI magnetic field in the photosphere, Cheng et al. (2015) concluded that the formation of the continuously sigmoidal threads in the hot channel is mainly via tether-cutting reconnection in the chromosphere, which gradually transforms the sheared arcades into twisted field lines prior to the eruption. Joshi et al. (2015) also studied this event but mostly focused on the flare. They found that the beginning of the flare is triggered by the eruption of the pre-existing sigmoidal hot channel, which is also identified as the MFR (see their Figure 3).
The SOL2014-09-10T hot channel is slightly different from the above three ones. The sigmoidal emission pattern of the source AR is mostly a result of two groups of sheared arcades, which manifest a sigmoid as a whole (Figure 5). However, through further checking the property of the magnetic field at the sigmoid center in the period of five hours before the eruption, Cheng et al. (2015) found that magnetic reconnection has been playing a role in transforming part of sheared arcades into longer threads. Recently, using Grad-Rubin method, Zhao et al. (2016) extrapolated the 3D NLFFF field of the AR two hours before the eruption and also found that some twisted field lines are comparable with the sigmoidal loops with, though, a small twist number (1.5). Actually, in the early phase of the eruption, we do see a set of continuously sigmoidal threads that detach from the source sigmoidal AR and then erupt outward.
By carefully inspecting the orientation of the four hot channels, we also find that it has some deviation from that of the main PILs of the ARs (red dashed lines in Figure 6f–9f). It indicates that the hot channels are most likely located in the corona with a large separation in altitude from the PILs of the source ARs in the photosphere. Owing to the projection effect, the SOL2013-04-11T (SOL2014-04-18T) hot channel in the ARs located at N07E13 (S20W34) deviates from the PIL toward the east (west). Moreover, in the events except the SOL2014-04-18T one, the filaments appear to lie nicely along the main PILs as seen in the AIA 304 Å images. Interestingly, during the CME eruptions, the filaments do not erupt. It shows again that the erupted hot channels are magnetic structures different from the source sigmoids and filaments.
To sum up, we can say that, although the four sigmoids seen in the EUV images may be comprised of two groups of sheared arcades, the hot channels should be the MFRs since they have ascended to a certain height in the corona and thus separated from the sigmoids, at least in the slow rise phase. Taking advantage of the AIA 131 Å and 94 Å images, we identify and depict the boundaries of the four hot channel-like MFRs, as shown (yellow lines) in Figures 6–9.
From Figure 6b, 7b, 8b, 9b, and Table 2, we find that all of hot channels have a similar footpoint separation distance. We introduce a parameter , the ratio of the projected length of the hot channels to the footpoint separation, to roughly estimate the writhe of the axis of the hot channels. The projected lengths of the four channels, which actually are the lower limit of the real lengths, range from 183 to 312 Mm. They are much larger than the separation distances (67–80 Mm) of the two footpoints. It indicates that the hot channels not only detach from the PILs but also have a significant writhe. One can see that the SOL2012-07-12T and SOL2013-04-11T hot channels have a very large , indicating a strong writhe, in addition to the twist in the threads of the hot channels (Figures 6 and 7). By contrast, the SOL2014-04-18T and SOL2014-09-10T hot channels have a relatively small , which is consistent with a weak writhe in the corresponding threads (Figures 8 and 9).
4.2. Overall Properties of Source Regions of the MFRs
We first study the magnetic properties of the ARs that host the hot channels. Here we use the HARP cylindrical equal-area vector field data (Figures 11a, 13a, 15a, and 17a), in which the magnetic field, , has been remapped to a lambert cylindrical equal-area projection, i.e., the Heliocentric Earth Ecliptic Cartesian coordinate system, and includes projection correction and geometric modification (Sun, 2013). Based on the Ampre’s law, the vertical current density can be calculated through /, where =410 G m A. Note that, in order to avoid the effect of noise in the weak fields, the calculation is only done for 50 G. Moreover, assuming that the uncertainty in the transverse field is 100 G, following the formula (Gary & Demoulin, 1995), the error of is estimated to be of the order of 0.02 A m. We also compare the values of calculated from the HARP data with that calculated from the vector data that are reverted with the UNNOFIT code (Janvier et al., 2014), which are shown to be almost the same.
Figures 11b, 13b, 15b, and 17b show the distribution of the vertical current density in the ARs. We find that the strongest current density always appears near the PILs, in particular for the SOL2012-07-12T and SOL2014-09-10T events. This is consistent with previous studies, which stated that magnetic free energy is mainly stored there (e.g., Sun et al., 2012). Moreover, the strong current density at the PILs also indicates that magnetic reconnection probably occurs there, which converts the sheared arcades into the twist field so as to form the MFRs (e.g., Cheng et al., 2014b; Vemareddy & Zhang, 2014). However, it should be noticed that those currents seem not to be related to the erupted hot channels. For example, in the SOL2012-07-12T, SOL2013-04-11T, and SOL2014-04-18T events, the filaments, which usually denote the locations with strong currents along the PILs, do not erupt following the eruption of the hot channels.
Figures 11c, 13c, 15c, and 17c show the distribution of the inclination angle in the ARs. One can see that the inclination angle increases from the sunspot center to outward and reaches its maximum in the penumbra. This conforms with the canopy-type configuration of ARs. Moreover, along the PILs, the inclination angle is also large, which is, however, due to the shear of the magnetic field, as shown in Figures 11a, 13a, 15a, and 17a. Due to the oppositely directed movements of magnetic footpoints (see attached online movies), the original potential field become more and more sheared, thus causing the large inclination angle.
Previous studies found that, in well-isolated ARs, the net current (the sum of direct current and return current ) over the positive polarity is qual to that over the negative one, i.e., the total current over the whole ARs is close to zero (Wheatland et al., 2000; Georgoulis et al., 2012). However, whether the net current over one polarity is also close to zero is debated. Here, we also tentatively calculate the net current over one polarity in the four ARs. Note that, for the forward sigmoidal ARs, which usually have a right-handed twist, the direct (return) current is parallel (antiparallel) with the direction of the magnetic field; while for the reversed sigmoidal ARs that usually have a left-handed twist, the direct (return) current is antiparallel (parallel) with the direction of the magnetic field (Rust & Kumar, 1996; Green et al., 2007; Török et al., 2014). The results are shown in Table 2, from which one can see that for the SOL2012-07-12T and SOL2014-09-10T cases, the net current significantly deviates from zero. While for the SOL2013-04-11T and SOL2014-04-18T cases, the net current almost vanishes.
4.3. Magnetic Properties of Footpoints of the MFRs
We consider the two ends of the hot channels as the footpoints of the MFRs (From Figures 6f, 7f, 8f, and 9f). We find that the two footpoints do not overlap the two ends of the main PILs. The fact is that, for all of hot channels, one footpoint comes from the moss region near the following sunspot, while the other one stems from the penumbra (SOL2014-04-18T) or penumbra edge (SOL2012-07-12T, SOL2013-04-11T, and SOL2014-09-10T) of the preceding sunspot. In the moss region, the magnetic field is weaker and the inclination angle is smaller; while in the penumbra and penumbra edge, the magnetic field is stronger and the inclination angle is larger.
It is noticed that the strong currents appear not only near the PILs, but also at the footpoints of the hot channels with the relatively strong magnetic field (Figures 11b, 13b, 15b, and 17b), the latter of which may be related to the currents of the MFRs. For the SOL2012-07-12T, SOL2014-04-18T, and SOL2014-09-10T (SOL2013-04-11T) events, at the positive (negative) footpoints, there always appears a dominated direct current. Moreover, we find that the distribution of the direct current takes on a curled shape on a large scale, which is surrounded by some sporadic return current, in particular for the SOL2012-07-12T (Figure 11b and 11h), SOL2013-04-11T (Figure 13b and 13h), and SOL2014-09-10T events (Figure 17b and 17h) (also see, Janvier et al., 2014).
In order to calculate the average magnetic field strength, average inclination angle, direct currents, and return currents at the footpoints, we choose two square regions with their centers determined from the ends of the hot channels but slightly shifted so as to cover most of the brightenings at the ends of the hot channels in the early rise phase as seen in the AIA 304 Å and 1600 Å passbands. The size of the square is set to be 30″30″(22 Mm22 Mm), which is big enough to contain all of the footpoints of the hot channels. On the other hand, such a size is also small enough for the sake of avoiding the influence of the surrounding magnetic field. Furthermore, we also compare the selected regions with the synthesized flare ribbons, which are derived through summing up all flare ribbons in the period of the early phase of the flares. We find that for each event, the two selected regions basically cover the hooks of the two J-shaped flare ribbons, in particular, the ribbons with stronger magnetic fields (Figure 10). Finally, it is needed to determine the locations of the footpoints in the coordinate system of the HARPs data. Here, this is done by means of co-aligning the same magnetic features of the footpoints in the light-of-sight magentograms with that in the vertical component of the HARPs vector magnetograms. It is worth noticing that the real size (30″30″) of the two square regions in the light-of-sight magnetograms is slightly larger than the size (22 Mm22 Mm) of the corresponding two square regions in the HARPs vector magnetograms because the former has a projection effect. However, we vary the center and/or size of the region by within 5 Mm as a test and do not find that our results are influenced qualitatively although the values of magnetic parameters are changed slightly. Note that, for the SOL2012-07-12T and SOL2014-09-10 events, there appears a secondary ribbon (Zhang et al., 2014) that extends from the tip of the hook (arrows in Figure 10). With the development of the flare, it further propagates toward a particular direction, forming a long and slender ribbon-like structure (see the attached online movies of Figure 6 and 9).
Table 2 shows the magnetic properties at the two regions where the hot channels root in. We find that at the footpoints with stronger magnetic fields, the average inclination angle tends to be larger with the maximal (minimal) value being 45 (23). The direct current there is calculated to be of the order of 10 A. While at the footpoints with weaker magnetic fields, the average inclination angle tends to be smaller with the maximal (minimal) value being 34 (18). The direct and return current is distributed nearly uniformly without obvious concentrations in some specific regions. Only for the SOL2012-07-12T event, the direct current there is comparable in quantity to that at the footpoints with stronger magnetic fields. One may expect that, the direct current flowing into the MFRs from one footpoint should flow out from the other, or that the direct current at the footpoints with weaker magnetic fields should be roughly equal to that at the footpoints with stronger magnetic fields. However, the real situation is complicated. Some potential reasons that cause the direct current at the footpoints with weaker magnetic fields significantly deviating from that at the other ones probably include: (1) the measurement of the vector field is not precise enough in the weak field region, which leads to a relatively large uncertainty in the current calculation there, and (2) the magnetic field strength is different at two different polarities, which requires that the sizes of regions for integrating the direct currents are also different to ensure the magnetic flux balance. However, observationally, the real sizes of the regions that contain the entire footpoints of MFRs can hardly be determined accurately.
We also study the evolution of the average magnetic field strength, average inclination angle, and direct current at the two footpoints during the eruptions of the hot channels (Figures 12, 14, 16, and 18). The average magnetic field strength has no significant change with the eruptions, as a consequence of no obvious flux emergence or magnetic cancellation appearing. However, we find that the average inclination angle at most of the footpoints tends to decrease. For example, at the two footpoints of the SOL2013-04-11T and SOL2014-04-18T hot channels, the inclination angle significantly decreases. While the inclination angle only slightly decreases at the footpoints of the SOL2014-09-10T hot channel and it even increases at the footpoint of the SOL2012-07-12T hot channel with positive magnetic polarity. We suspect that such an increase is the result of the rapid rotation of the preceding sunspot (see Figure 7 of Cheng et al., 2014b), which compensates for the decrease of the inclination angle caused by the eruption of the hot channel. Moreover, we also find that the direct current at all of the footpoints with stronger magnetic fields experiences an obvious decrease (Figures 12c, 14c, 16c, and 18c). These features indicate that the magnetic field emanating from the footpoints tends to become more straight and less twisted after the eruptions. This is reasonable because the legs of the hot channel-like MFRs are gradually stretched and straightened. At the same time, such a stretching causes an increase of the length of the hot channels. However, due to the conservation of the total twist, the twist per unit length just decreases, leading to the reduction of the direct current.
4.4. Evolution of Footpoints of the MFR Eruptions
The detailed kinematical studies suggest that the MFRs play an important role in forming the CMEs and producing the flare emission (Cheng et al., 2013a, b; Sun et al., 2015). However, during the MFR eruption, the detailed evolution of the magnetic field in the MFR and its relation to the ambient field are still unknown. In this Section, we pay our attention to the chromospheric brightenings caused by the four erupted hot channels, in particular studying their evolution and relation to flare ribbons.
During the eruptions, the overall evolution of the chromospheric brightenings can be divided into three stages. The first stage occurs in the slow rise phase. The brightenings sporadically appear at the footpoints and in the regions below the two elbows of the hot channels (panels a, d, and g of Figure 19–22). It indicates that some pre-heating probably has started, most likely taking place within the threads of the hot channels or between the hot channels and the ambient sheared arcades.
The second stage of chromospheric brightenings starts with the beginning of the impulsive acceleration of the hot channels, also the early rise phase of the flare (panels b, e, and h of Figure 19–22). The brightenings at the two footpoints first extend outward, forming a hook-like shape. At the same time, the brightenings also propagate along the main PILs with opposite orientations at two different polarities. Quickly, two slender ribbon-like brightenings come into being, taking on sheared double-J shape, the straight parts of which are adjacent to each other on either side of the PILs and the two hooks face on each other at the opposite ends. Such characters are especially apparent for the SOL2013-04-11T and SOL2014-04-18T events. For the SOL2012-07-12T and SOL2014-09-10T events, the J ribbon in the negative polarity is much larger than that in the positive polarity, which is primarily due to that the magnetic field is very strong there and thus the brightenings only propagate with a relatively short distance.
The third stage happens in the later rise and decay phase of the flare (panels c, f, and i of Figure 19–22). The primary characteristics are the separation of the two slender ribbons and the gradual disappearance of the two hooks. From the attached online movies of Figures 6–9, one can clearly see that the two straight parts of the J ribbons expand outward with opposite directions but almost perpendicular to the PILs. Based on the CSHKP model, such a separation is mainly caused by the systematic ascending of the reconnection site in the corona. Moreover, from the second to third stage, we find that the two straight parts of the double J structure evolve from greatly sheared to almost parallelized (less sheared). This is most likely the consequence of the flare reconnection that proceeds from the strongly sheared overlying field at a low site to a nearly potential one at a high site (Aulanier et al., 2012; Janvier et al., 2013).
5. Summary and Discussion
In this paper, we investigate the magnetic properties of four erupted hot channel-like MFRs that initially appear as forward or reversed sigmoidal structures prior to the eruption as seen in the AIA 131 and 94 Å passbands. Thanks to the high resolution of the AIA, the continuously sigmoidal threads that comprise the coherent and sigmoidal channels can be seen to detach from the source sigmoidal ARs and then erupt to become the CMEs. The ratio of the projected length of the hot channels to their footpoint separation distance ranges from 2.3 to 4.0, showing that the magnetic fields in the hot channels have a strong writhe.
The axis of the hot channels significantly deviates from the main PILs and associated filaments prior to the eruption. It shows that the channels have ascended to a high altitude in the corona and been likely separated from that of the filaments when approaching the eruptions. Using the optimization algorithm (Wheatland et al., 2000; Wiegelmann, 2004), we also extrapolate the 3D NLFFF structures in the four ARs. Unfortunately, we could not reproduce the field lines that fit spatially with the hot channels. The reasons could be the insufficient resolution of magnetic field measurement, in particular at the footpoints of the hot channels (Cheng et al., 2014b), or the preprocessing over-smoothes the vector field before doing the extrapolation. Moreover, the observed vector magnetograms often contain significant Lorentz and buoyancy forces, which actually do not satisfy the NLFFF assumption (De Rosa et al., 2009), thus challenging the reliability of extrapolated structures. At present, the potential method that is likely able to reproduce the magnetic field resembling the hot channels is the MFR insertion method developed by van Ballegooijen (2004), which inserts an MFR into the potential field and then perform a relaxation to fit the observations (e.g., Su & van Ballegooijen, 2012; Su et al., 2015; Savcheva et al., 2012b, a, 2015). In a follow up work, we plan to utilize this method to model the structures of the four ARs and then compare the results with observations. We will also try the promising Grad-Rubin method, which has been applied to the SOL2014-09-10 event by Zhao et al. (2016) and reproduced twisted field lines that are comparable with the sigmoidal threads of the hot channel.
Through comparing the AIA 131 Å and 94 Å images and HMI light-of-sight magnetograms, we find that one footpoint of the hot channels originates in the penumbra or penumbra edge while the other one is from the moss region. Such an observational fact can be compared with various numerical models. In the simulations of a pre-existing MFR emerging into the potential magnetic arcades (e.g., Fan & Gibson, 2004; Fan, 2010; Leake et al., 2013), the two footpoints of the MFR are exactly located at the centers of two different polarities. The situation is somewhat different in the simulations by Aulanier et al. (2010) and Amari et al. (2010), in which an initially potential and asymmetric bipolar field evolves by means of the shearing motions along the main PIL. The flux cancellation in a bald-patch and/or the reconnection in a quasi-separatrix layer (QSL) transform the sheared arcades into continuously sigmoidal field, finally forming an MFR. The footpoints of the MFRs are located at the outer boundaries of the two polarities while the axis of the MFRs is basically parallel to the main PIL. Recently, Amari et al. (2014) performed a cutting-edge data-driven simulation and completely reproduce that one footpoint of the MFR is located at the edge of the positive polarity and the other one stems from the moss region nearby the negative polarity. Unfortunately, for the event they modeled, there are no AIA observations; thus it is difficult to compare the magnetic field lines with the hot channels.
We find that the net currents in SOL2012-07-12T and SOL2014-09-10T ARs significantly deviate from zero, but in the SOL2013-04-11T and SOL2014-04-18T ARs they almost vanish. Nevertheless, at the footpoints of the four hot channels, we always see a dominated direct current. It indicates that the direct currents in the MFRs are always larger than the return currents irrespective of whether the currents in their source ARs are neutralized or not. In the ARs with significant non-zero net currents, the MFRs may be directly originated from the net currents; while in ARs with almost neutralized net currents, some opposite currents possibly exist elsewhere and can counteract the currents of the MFRs when integrating the net currents over the entire polarity. In fact, in many CME/flare initiation models (e.g., Titov & Démoulin, 1999; van Ballegooijen, 2004; Török & Kliem, 2005; Kliem & Török, 2006; Savcheva et al., 2012b), the existence of the MFR-associated current is required; one can see a current that comes out from one footpoint of the MFR and returns to the other one. However, it is worth noting that the direct currents at the footpoints with weaker magnetic field are almost one order of magnitude smaller than that at the footpoints with stronger magnetic field. This is somewhat different from the result of Janvier et al. (2014), in which both of the footpoints of the erupted MFR originate in the regions with very strong magnetic field and exhibit comparative direct currents. Such a distinction may be due to the fact that, in our cases, the measurement of the vector magnetic field at the footpoints where the magnetic field is weak seems not to be accurate enough, which may lead to a relatively large uncertainty in the current calculation.
We study the evolution of the magnetic field at the footpoints of the hot channels. The average magnetic field strength does not change apparently with the eruptions going on. However, the average inclination angle at most of the footpoints has a decrease. The direct current at all of the footpoints with stronger magnetic field also experiences a decrease. At first glance, these results seem to conflict with the previous results that the photospheric magnetic field near the PIL becomes more horizontal (Wang & Liu, 2010, 2015; Wang et al., 2012; Liu et al., 2012), implying an increase of the inclination angle there, and the direct current tends to increase slightly after the eruptions (Janvier et al., 2014). We argue that the different results may refer to the different physical processes that occur at different locations. With the MFR erupting upward, the magnetic field rooted in the two ends should be stretched to be more and more vertical. Meanwhile, due to the total twist conservation, the twist per unit length decreases. This results in a decrease of both the inclination angle and the direct current. By contrast, the magnetic field near the PIL is first stretched upward and then reconnects. Afterwards, the reconnected field lines shrink toward the photosphere due to magnetic tension, thus leading to an increase in the inclination angle. However, as for the increase of the direct current near the main PIL, it is possibly a result of magnetic reconnection or due to some other unknown mechanisms.
We also investigate the chromospheric brightenings at the footpoints of the hot channel-like MFRs, as well as their evolution and relation to flare ribbons. The overall evolution of the chromospheric brightenings is comprised of three stages. The brightenings in the first stage appear in the slow rise phase and is located at the two footpoints, as well as in the regions below the two elbows of the hot channels. The brightenings at the second stage take on a double J shape with the two hooks at the opposite ends corresponding to the extended footpoints. Once entering the third stage, the two straight parts of double J-ribbons start to separate from each other, evolving from two sheared slender ribbons into parallel widened ones.
The evolution of the chromospheric brightenings provides a perspective on the 3D evolution of the MFR, which is illustrated by a schematic drawing as shown in Figure 23. In the slow rise stage, the reconnection may take place inside the MFR or between the MFR and its ambient strongly sheared field. As a result, the footpoints of the MFR and the strongly sheared field are heated, manifesting themselves as some sporadical brightenings (Figure 23a). At the same time, the reconnection also produces an upward Lorentz force to lift the MFR. Once the MFR ascends to a height where the decline of the background field is fast enough, kink and/or torus instability will be set in (e.g., Török & Kliem, 2005; Kliem & Török, 2006; Olmedo & Zhang, 2010; Savcheva et al., 2012a; Cheng et al., 2013b; Zuccarello et al., 2015; Su et al., 2015). Then, the MFR accelerates upward nonlinearly. In the early acceleration stage (Figure 23b), the overlying field that straddles over the MFR is strongly sheared. When these fields are stretched by the MFR upward motion, their legs approach each other and then reconnect, resulting in two sheared J-shaped brightenings in the chromosphere. This result is consistent with previous findings that hard X-ray footpoint sources and post-flare loops sometimes display a strong shear in the initial phase of the flare (e.g., Su et al., 2007; Ji et al., 2008; Liu et al., 2009; Guo et al., 2012). Actually, the two strongly sheared J-shaped brightenings may represent the footprints of the curved QSLs in the chormosphere as expected by the 3D standard CME/flare model (Aulanier et al., 2012; Janvier et al., 2013, 2014; Dudik et al., 2016). Very recently, through calculating the QSLs in several sigmoidal ARs, Savcheva et al. (2015) and Zhao et al. (2016) confirmed that the straight segments of the two sheared J-ribbons are basically matched by the MFR-related QSLs, thus approving the interpretation of the reconnection between strongly sheared overlying field in the early phase of the eruption.
In the later acceleration stage (Figure 23c), the most apparent evolution is the separation of two flare ribbons. Since the MFR has risen to a relatively high altitude, the overlying field is almost potential. When a row of such overlying field lines is stretched and then reconnects, the heated footpoints will form two parallel ribbon-like brightenings. With the MFR continuously rising, the overlying field being reconnected becomes higher and higher; correspondingly, the separation of their footpoints becomes larger and larger. As a consequence, the increase in the reconnection height is mapped as the separation of two flare ribbons.
Note: () denotes magnetic flux in the positive (negative) polarity of ARs. refers to the imbalance of magnetic fluxes defined as (+)/(–). sign indicates the sign of the magnetic polarity where we calculate direct currents , return current , and net current =+. is the ratio of the net current to total unsigned current and defined as /(–). The quantities and are the projected length and footpoint separation of the hot channels, respectively, and =/. (), ( ), (), (), (), and () represent average magnetic field strength, magnetic flux, direct current, return current, net current, and average inclination angle at the negative (positive) footpoints of the four hot channels, respectively.
- slugcomment: Accepted for publication in ApJ Supplement Series
- Amari, T., Aly, J.-J., Mikic, Z., & Linker, J. 2010, ApJ, 717, L26
- Amari, T., Canou, A., & Aly, J.-J. 2014, Nature, 514, 465
- Antiochos, S. K., DeVore, C. R., & Klimchuk, J. A. 1999, ApJ, 510, 485
- Aulanier, G., Janvier, M., & Schmieder, B. 2012, A&A, 543, A110
- Aulanier, G., Török, T., Démoulin, P., & DeLuca, E. E. 2010, ApJ, 708, 314
- Bak-Stȩślicka, U., Gibson, S. E., Fan, Y., Bethge, C., Forland, B., & Rachmeler, L. A. 2013, ApJ, 770, L28
- Bobra, M. G., Sun, X., Hoeksema, J. T., Turmon, M., Liu, Y., Hayashi, K., Barnes, G., & Leka, K. D. 2014, Sol. Phys., 289, 3549
- Brueckner, G. E., Howard, R. A., Koomen, M. J., Korendyke, C. M., et al. 1995, Sol. Phys., 162, 357
- Canfield, R. C., Hudson, H. S., & McKenzie, D. E. 1999, Geophys. Res. Lett., 26, 627
- Carmichael, H. 1964, NASA Special Publication, 50, 451
- Chen, B., Bastian, T. S., & Gary, D. E. 2014a, ApJ, 794, 149
- Chen, H., Zhang, J., Cheng, X., Ma, S., Yang, S., & Li, T. 2014b, ApJ, 797, L15
- Chen, J. 1996, J. Geophys. Res., 101, 27499
- Chen, P. F. 2011, Living Reviews in Solar Physics, 8, 1
- Chen, P. F., Harra, L. K., & Fang, C. 2014c, ApJ, 784, 50
- Chen, P. F., & Shibata, K. 2000, ApJ, 545, 524
- Cheng, X., Ding, M. D., & Fang, C. 2015, ApJ, 804, 82
- Cheng, X., Ding, M. D., Zhang, J., Srivastava, A. K., Guo, Y., Chen, P. F., & Sun, J. Q. 2014a, ApJ, 789, L35
- Cheng, X., Ding, M. D., Zhang, J., Sun, X. D., Guo, Y., Wang, Y. M., Kliem, B., & Deng, Y. Y. 2014b, ApJ, 789, 93
- Cheng, X., Zhang, J., Ding, M. D., Liu, Y., & Poomvises, W. 2013a, ApJ, 763, 43
- Cheng, X., Zhang, J., Ding, M. D., Olmedo, O., Sun, X. D., Guo, Y., & Liu, Y. 2013b, ApJ, 769, L25
- Cheng, X., Zhang, J., Liu, Y., & Ding, M. D. 2011, ApJ, 732, L25
- Cheng, X., et al. 2014c, ApJ, 780, 28
- Chintzoglou, G., Patsourakos, S., & Vourlidas, A. 2015, ApJ, 809, 34
- De Rosa, M. L., Schrijver, C. J., Barnes, G., Leka, K. D., & others. 2009, ApJ, 696, 1780
- Dove, J. B., Gibson, S. E., Rachmeler, L. A., Tomczyk, S., & Judge, P. 2011, ApJ, 731, L1
- Dudík, J., Janvier, M., Aulanier, G., Del Zanna, G., Karlický, M., Mason, H. E., & Schmieder, B. 2014, ApJ, 784, 144
- Dudik, J., et al. 2016, ArXiv e-prints
- Fan, Y. 2010, ApJ, 719, 728
- Fan, Y., & Gibson, S. E. 2004, ApJ, 609, 1123
- —. 2007, ApJ, 668, 1232
- Forbes, T. G., & Isenberg, P. A. 1991, ApJ, 373, 294
- Gary, G. A., & Demoulin, P. 1995, ApJ, 445, 982
- Georgoulis, M. K., Titov, V. S., & Mikić, Z. 2012, ApJ, 761, 61
- Gibson, S. E., Fan, Y., Mandrini, C., Fisher, G., & Demoulin, P. 2004, ApJ, 617, 600
- Green, L. M., Kliem, B., Török, T., van Driel-Gesztelyi, L., & Attrill, G. D. R. 2007, Sol. Phys., 246, 365
- Guo, Y., Ding, M. D., Schmieder, B., Démoulin, P., & Li, H. 2012, ApJ, 746, 17
- Guo, Y., Ding, M. D., Schmieder, B., Li, H., Török, T., & Wiegelmann, T. 2010a, ApJ, 725, L38
- Guo, Y., Schmieder, B., Démoulin, P., Wiegelmann, T., Aulanier, G., Török, T., & Bommier, V. 2010b, ApJ, 714, 343
- Hirayama, T. 1974, Sol. Phys., 34, 323
- Inoue, S., Hayashi, K., Shiota, D., Magara, T., & Choe, G. S. 2013, ApJ, 770, 79
- Janvier, M., Aulanier, G., Bommier, V., Schmieder, B., Démoulin, P., & Pariat, E. 2014, ApJ, 788, 60
- Janvier, M., Aulanier, G., Pariat, E., & Démoulin, P. 2013, A&A, 555, A77
- Ji, H., Wang, H., Liu, C., & Dennis, B. R. 2008, ApJ, 680, 734
- Ji, H., Wang, H., Schmahl, E. J., Moon, Y.-J., & Jiang, Y. 2003, ApJ, 595, L135
- Jiang, C., Wu, S. T., Feng, X., & Hu, Q. 2014a, ApJ, 780, 55
- —. 2014b, ApJ, 786, L16
- Jiang, C.-W., & Feng, X.-S. 2016, Research in Astronomy and Astrophysics, 16, 015
- Jiang, C.-W., Wu, S.-T., Feng, X.-S., & Hu, Q. 2016, Research in Astronomy and Astrophysics, 16, 018
- Joshi, N. C., Liu, C., Sun, X., Wang, H., Magara, T., & Moon, Y.-J. 2015, ApJ, 812, 50
- Karpen, J. T., Antiochos, S. K., & DeVore, C. R. 2012, ApJ, 760, 81
- Kliem, B., Titov, V. S., & Török, T. 2004, A&A, 413, L23
- Kliem, B., & Török, T. 2006, Phys. Rev. Lett., 96, 255002
- Kopp, R. A., & Pneuman, G. W. 1976, Sol. Phys., 50, 85
- Leake, J. E., Linton, M. G., & Török, T. 2013, ApJ, 778, 99
- Lemen, J. R., et al. 2012, Sol. Phys., 275, 17
- Li, L. P., & Zhang, J. 2013a, A&A, 552, L11
- Li, T., & Zhang, J. 2013b, ApJ, 770, L25
- —. 2013c, ApJ, 778, L29
- Li, X., Morgan, H., Leonard, D., & Jeska, L. 2012, ApJ, 752, L22
- Liu, C., Deng, N., Lee, J., Wiegelmann, T., Moore, R. L., & Wang, H. 2013, ApJ, 778, L36
- Liu, C., Lee, J., Yurchyshyn, V., Deng, N., Cho, K.-s., Karlický, M., & Wang, H. 2007, ApJ, 669, 1372
- Liu, C., et al. 2012, ApJ, 745, L4
- Liu, R., Liu, C., Wang, S., Deng, N., & Wang, H. 2010, ApJ, 725, L84
- Liu, W., Petrosian, V., Dennis, B. R., & Holman, G. D. 2009, ApJ, 693, 847
- Liu, Y. D., et al. 2014, Nature Communications, 5, 3481
- Low, B. C., & Hundhausen, J. R. 1995, ApJ, 443, 818
- Lynch, B. J., & Edmondson, J. K. 2013, ApJ, 764, 87
- Mackay, D. H., Karpen, J. T., Ballester, J. L., Schmieder, B., & Aulanier, G. 2010, Space Sci. Rev., 151, 333
- McKenzie, D. E., & Canfield, R. C. 2008, A&A, 481, L65
- Moore, R. L., Sterling, A. C., Hudson, H. S., & Lemen, J. R. 2001, ApJ, 552, 833
- Nindos, A., Patsourakos, S., Vourlidas, A., & Tagikas, C. 2015, ApJ, 808, 117
- Olmedo, O., & Zhang, J. 2010, ApJ, 718, 433
- Patsourakos, S., Vourlidas, A., & Stenborg, G. 2013, ApJ, 764, 125
- Pesnell, W. D., Thompson, B. J., & Chamberlin, P. C. 2012, Sol. Phys., 275, 3
- Priest, E. R., & Forbes, T. G. 2002, A&A Rev., 10, 313
- Rust, D. M., & Kumar, A. 1996, ApJ, 464, L199
- Savcheva, A., Pariat, E., McKillop, S., McCauley, P., Hanson, E., Su, Y., Werner, E., & DeLuca, E. E. 2015, ApJ, 810, 96
- Savcheva, A., Pariat, E., van Ballegooijen, A., Aulanier, G., & DeLuca, E. 2012a, ApJ, 750, 15
- Savcheva, A. S., van Ballegooijen, A. A., & DeLuca, E. E. 2012b, ApJ, 744, 78
- Schmieder, B., Aulanier, G., & Vršnak, B. 2015, Sol. Phys., 290, 3457
- Schou, J., et al. 2012, Sol. Phys., 275, 229
- Shen, Y., Liu, Y., & Su, J. 2012, ApJ, 750, 12
- Shibata, K., Masuda, S., Shimojo, M., Hara, H., Yokoyama, T., Tsuneta, S., Kosugi, T., & Ogawara, Y. 1995, ApJ, 451, L83
- Song, H. Q., Zhang, J., Chen, Y., & Cheng, X. 2014a, ApJ, 792, L40
- Song, H. Q., Zhang, J., Cheng, X., Chen, Y., Liu, R., Wang, Y. M., & Li, B. 2014b, ApJ, 784, 48
- Srivastava, A. K., Zaqarashvili, T. V., Kumar, P., & Khodachenko, M. L. 2010, ApJ, 715, 292
- Sterling, A. C., Hudson, H. S., Thompson, B. J., & Zarro, D. M. 2000, ApJ, 532, 628
- Sturrock, P. A. 1966, Nature, 211, 695
- Su, Y., Golub, L., & Van Ballegooijen, A. A. 2007, ApJ, 655, 606
- Su, Y., & van Ballegooijen, A. 2012, ApJ, 757, 168
- Su, Y., van Ballegooijen, A., McCauley, P., Ji, H., Reeves, K. K., & DeLuca, E. E. 2015, ApJ, 807, 144
- Sun, J. Q., et al. 2015, Nature Communications, 6, 7598
- Sun, X. 2013, ArXiv e-prints
- Sun, X., Hoeksema, J. T., Liu, Y., Wiegelmann, T., Hayashi, K., Chen, Q., & Thalmann, J. 2012, ApJ, 748, 77
- Titov, V. S., & Démoulin, P. 1999, A&A, 351, 707
- Török, T., & Kliem, B. 2005, ApJ, 630, L97
- Török, T., Kliem, B., & Titov, V. S. 2004, A&A, 413, L27
- Török, T., et al. 2014, ApJ, 782, L10
- Tripathi, D., Reeves, K. K., Gibson, S. E., Srivastava, A., & Joshi, N. C. 2013, ApJ, 778, 142
- Valori, G., Kliem, B., Török, T., & Titov, V. S. 2010, A&A, 519, A44
- van Ballegooijen, A. A. 2004, ApJ, 612, 519
- van Ballegooijen, A. A., Cartledge, N. P., & Priest, E. R. 1998, ApJ, 501, 866
- Vemareddy, P., & Zhang, J. 2014, ApJ, 797, 80
- Vourlidas, A., Lynch, B. J., Howard, R. A., & Li, Y. 2013, Sol. Phys., 284, 179
- Wang, H., & Liu, C. 2010, ApJ, 716, L195
- —. 2015, Research in Astronomy and Astrophysics, 15, 145
- Wang, S., Liu, C., Liu, R., Deng, N., Liu, Y., & Wang, H. 2012, ApJ, 745, L17
- Wang, Y.-M., & Stenborg, G. 2010, ApJ, 719, L181
- Webb, D. F., Forbes, T. G., Aurass, H., Chen, J., Martens, P., Rompolt, B., Rusin, V., & Martin, S. F. 1994, Sol. Phys., 153, 73
- Wheatland, M. S., Sturrock, P. A., & Roumeliotis, G. 2000, ApJ, 540, 1150
- Wiegelmann, T. 2004, Sol. Phys., 219, 87
- Yan, X. L., Xue, Z. K., Liu, J. H., Ma, L., Kong, D. F., Qu, Z. Q., & Li, Z. 2014, ApJ, 782, 67
- Yan, X. L., Xue, Z. K., Pan, G. M., Wang, J. C., Xiang, Y. Y., Kong, D. F., & Yang, L. H. 2015, ApJS, 219, 17
- Yang, S., Zhang, J., Liu, Z., & Xiang, Y. 2014, ApJ, 784, L36
- Yashiro, S., Gopalswamy, N., Michalek, G., St. Cyr, O. C., Plunkett, S. P., Rich, N. B., & Howard, R. A. 2004, Journal of Geophysical Research (Space Physics), 109, 7105
- Zhang, J., Cheng, X., & Ding, M.-D. 2012, Nature Communications, 3, 747
- Zhang, J., Li, T., & Yang, S. 2014, ApJ, 782, L27
- Zhang, J., Yang, S. H., & Li, T. 2015, A&A, 580, A2
- Zhao, J., Gilchrist, S. A., Aulanier, G., Schmieder, B., Pariat, E., & Li, H. 2016, ArXiv e-prints
- Zuccarello, F. P., Aulanier, G., & Gilchrist, S. A. 2015, ApJ, 814, 126
- Zuccarello, F. P., Soenen, A., Poedts, S., Zuccarello, F., & Jacobs, C. 2008, ApJ, 689, L157