Issue |
A&A
Volume 647, March 2021
|
|
---|---|---|
Article Number | A171 | |
Number of page(s) | 9 | |
Section | The Sun and the Heliosphere | |
DOI | https://doi.org/10.1051/0004-6361/202039944 | |
Published online | 30 March 2021 |
How flux feeding causes eruptions of solar magnetic flux ropes with the hyperbolic flux tube configuration
1
CAS Key Laboratory of Geospace Environment, School of Earth and Space Sciences, University of Science and Technology of China, Hefei 230026, PR China
e-mail: zhangqh@ustc.edu.cn
2
CAS Center for Excellence in Comparative Planetology, University of Science and Technology of China, Hefei 230026, PR China
3
Mengcheng National Geophysical Observatory, School of Earth and Space Sciences, University of Science and Technology of China, Hefei 230026, PR China
4
Collaborative Innovation Center of Astronautical Science and Technology, Hefei, Anhui 230026, PR China
5
School of Atmospheric Sciences, Sun Yat-sen University, Zhuhai, Guangdong 519000, PR China
6
Institute for the Study of Earth, Ocean, and Space, University of New Hampshire, Durham, NH 03824, USA
Received:
20
November
2020
Accepted:
26
January
2021
Coronal magnetic flux ropes are generally considered to be the core structure of large-scale solar eruptions. Recent observations found that solar eruptions could be initiated by a sequence of flux feeding, during which chromospheric fibrils rise upward from below, and merge with a pre-existing prominence. Further theoretical study has confirmed that the flux feeding mechanism is efficient in causing the eruption of flux ropes that are wrapped by bald patch separatrix surfaces. But it is unclear how flux feeding influences coronal flux ropes that are wrapped by hyperbolic flux tubes (HFTs), and whether it is able to cause the flux-rope eruption. In this paper, we use a 2.5-dimensional magnetohydrodynamic model to simulate the flux feeding processes in HFT configurations. It is found that flux feeding injects axial magnetic flux into the flux rope, whereas the poloidal flux of the rope is reduced after flux feeding. Flux feeding is able to cause the flux rope to erupt, provided that the injected axial flux is large enough so that the critical axial flux of the rope is reached. Otherwise, the flux-rope system evolves to a stable equilibrium state after flux feeding, which might be even farther away from the onset of the eruption, indicating that flux feeding could stabilize the rope system with the HFT configuration in this circumstance.
Key words: Sun: filaments / prominences / Sun: flares / Sun: coronal mass ejections (CMEs) / Sun: magnetic fields / Sun: activity
© ESO 2021
1. Introduction
Large-scale solar eruptions include prominence and filament eruptions, flares, and coronal mass ejections (CMEs; Benz 2008; Chen 2011; Parenti 2014; Liu 2020). They are capable of inflicting huge impacts on the solar-terrestrial system (Švestka 2001; Cheng et al. 2014; Shen et al. 2014; Lugaz et al. 2017; Gopalswamy et al. 2018). It is widely accepted that different kinds of large-scale solar eruptions are closely related to each other; they are essentially different manifestations of the same eruptive process of a coronal magnetic flux-rope system (Zhang et al. 2001; Vršnak et al. 2005; van Driel-Gesztelyi & Green 2015; Jiang et al. 2018; Liu et al. 2018; Yan et al. 2020). Therefore, it is of great significance to investigate how the eruption of coronal magnetic flux ropes is initiated. According to the magnetic topology, coronal flux ropes are classified into two types of configurations. First, if the flux rope sticks to the photosphere with a bald patch separatrix surface (BPSS, Titov et al. 1993; Titov & Démoulin 1999; Gibson & Fan 2006) wrapping the flux rope, this is usually called the BPS configuration (Filippov 2013). Second, for the flux-rope system in which the rope is suspended in the corona and wrapped around by a hyperbolic flux tube (HFT), it is called the HFT configuration (Titov et al. 2003; Aulanier et al. 2005; Chintzoglou et al. 2017).
Many theoretical analyses have been carried out to investigate the eruptive mechanism of coronal magnetic flux ropes. Apart from the well-known magnetohydrodynamic (MHD) instabilities (Romano et al. 2003; Török & Kliem 2003; Kliem & Török 2006; Fan & Gibson 2007; Aulanier et al. 2010; Guo et al. 2010; Savcheva et al. 2012), it was also suggested by many previous studies that catastrophes could be responsible for solar eruptions (e.g., Forbes & Isenberg 1991; Isenberg et al. 1993; Lin et al. 2001; Chen et al. 2007; Démoulin & Aulanier 2010; Longcope & Forbes 2014; Kliem et al. 2014). For example, if either the axial flux or the poloidal flux of a flux rope exceeds the corresponding critical value, a catastrophic loss of equilibrium will occur, resulting in the eruption of the rope (Li & Hu 2002; Bobra et al. 2008; Su et al. 2011; Zhang et al. 2016; Zhuang et al. 2018). Alternatively, magnetic reconnection may play a dominant role in flux-rope eruptions, such as the break-out model (Antiochos et al. 1999; Sterling & Moore 2004), tether-cutting model (Moore et al. 2001; Inoue et al. 2015), and flux emergence model (Chen & Shibata 2000; Archontis & Hood 2008). Recently, Zhang et al. (2014) observed that chromospheric fibrils rise from below a quiescent prominence and merge with this prominence, and called this phenomenon the “flux feeding” process. They observed that flux feeding processes occur three times, followed by the eruption of the quiescent prominence. This implies that coronal eruptions could be initiated by flux feeding.
Since accurate measurements of the local conditions in the corona are lacking, the physical essence of flux feeding mechanism is unclear based on observational results alone. To solve this, Zhang et al. (2020; hereafter Paper I) carried out numerical simulations to investigate the scenario of the flux feeding mechanism in BPS configurations. In their simulations, the rising chromospheric fibril is represented by a small flux rope emerging from below a pre-existing coronal magnetic flux rope with a BPS configuration. It was found that flux feeding processes only inject axial magnetic flux into the pre-existing flux rope, whose poloidal flux, however, remains unchanged after flux feeding. If the amount of injected axial flux is large enough so that the total axial flux of the flux rope exceeds its critical value, the eruption of the rope is initiated by the upward catastrophe associated with the critical axial flux; otherwise, the flux rope remains in the stable BPS configuration after flux feeding. Different from the BPS configurations, HFT configurations are usually considered to be the pre-eruptive states of coronal eruptions (Galsgaard et al. 2003; Aulanier et al. 2010; Filippov 2013), indicating that HFT configurations could be metastable or even unstable. Several questions arise: How does flux feeding influence the flux-rope systems with an HFT configuration? Is flux feeding still able to cause eruptions? If the flux feeding mechanism is still efficient, what are the similarities and differences between the BPS and the HFT cases? To resolve these questions, theoretical analyses about flux feeding in HFT configurations are needed, which will shed further light on the physical nature of the flux feeding mechanism.
In this paper we carry out 2.5-dimensional numerical simulations to assess the influence of flux feeding on a flux-rope system with an HFT configuration, especially to investigate how flux feeding causes the rope’s eruption. The paper is arranged as follows. The simulating procedures are introduced in Sect. 2. The simulation results of a typical flux feeding process and the flux-rope eruption caused by this process are presented in Sect. 3. The influence of flux feeding on the rope system and the physical scenario of flux feeding mechanism in HFT configurations are analyzed in Sect. 4. Finally, the discussion and conclusion are given in Sect. 5.
2. Simulating procedures
For the 2.5-dimensional cases in our simulations, all the magnitudes satisfy ∂/∂z = 0, thus the magnetic field can be denoted as
Here Bz is the component of the magnetic field in z-direction; ψ is the magnetic flux function, and the isolines of ψ correspond to the magnetic field lines projected in the x − y plane. In this paper the multi-step implicit scheme (Hu 1989) is used to simulate the evolution of coronal flux-rope systems. The initial state in our simulation is obtained by numerical procedures: first we insert a flux rope into a potential background field from the lower base, and then by adjusting the properties of this rope, a flux-rope system with an HFT configuration is eventually obtained. The basic equations and the detailed simulating procedures are introduced in Appendix A. The magnetic configuration of the initial state is shown in Fig. 1a: a flux rope is suspended above a bundle of closed arcades, resulting in an X-type magnetic structure below the rope. The green curves in Fig. 1a show the boundary of the flux rope and that of the arcades below the rope. The magnetic flux of these arcades per unit length along the z-direction is Φa = 0.069 × 1010 Mx cm−1. The background field is a partially open bipolar field, with a negative and a positive surface magnetic charge located at the lower photosphere within −25 Mm < x < −15 Mm and 15 Mm < x < 25 Mm, respectively. The value of Bz is zero in the background field.
![]() |
Fig. 1. Evolution of the magnetic configuration during a flux feeding process with CE = 2.23. Panel a: initial flux-rope system with an HFT configuration; the green curves indicate the boundary of the major flux rope and that of the arcades below the rope. The times indicated in panels b–i are in units of τA. Different colors depict the magnetic field lines corresponding to different values of ψ. |
As in Paper I, the emerging fibril in the scenario of flux feeding is represented by a small flux rope, and the pre-existing large flux rope is called the major rope for simplicity in the rest of this paper. The emergence of the small rope is achieved by simulating procedures similar to those in Paper I: the small rope emerges from the central region right below the major rope at a constant speed; the emergence begins at t = 0 and ends at t = τE = 60τA, where τA = 17.4 s is the characteristic Alfvén transit time. Thus, the emerged part of the small rope at the base is located within −xE ≤ x ≤ xE, where , hE = a(2t/τE − 1), 0 ≤ t ≤ τE, and a is the radius of the small rope. Based on this, by adjusting ψ, Bz, the velocities vx, y, z, the temperature T, and the density ρ at the base (y = 0, − xE ≤ x ≤ xE), the emergence of the small flux rope is achieved as
Here vE = 2a/τE is the emerging velocity; ψi is the magnetic flux function of the initial state, determined by numerical means (see Eq. (A.10)). It is noteworthy that ψ at the lower base is fixed at ψi except during the emergence of the small rope, so that the normal component of the photospheric magnetic field is unchanged, implying that the lower base corresponds to the photosphere (Lin et al. 2003). The parameter CE determines the magnetic field strength of the small rope; its dimensionless values given in this paper are in units of 0.373 × 1010 Mx cm−1. In the small and the major ropes Bz is positive and the component of the magnetic field in the x − y plane is counterclockwise. An anomalous resistivity is used in our simulations to restrict magnetic reconnection within the regions of current sheets:
Here ηm = 9.95 × 10−2, L0 = 107 m, and v0 = 128.57 km s−1; μ0 is the vacuum magnetic permeability; the critical current density is jc = 4.45 × 10−4 A m−2.
3. Simulation results
3.1. Flux feeding process
With the procedures introduced above we are able to simulate flux feeding processes in HFT configurations. The simulation results of a typical flux feeding process with CE = 2.23 are illustrated in Figs. 1b–i; different colors are used to depict the magnetic fields lines corresponding to different values of the flux function ψ so as to clearly demonstrate the interaction between the major rope and emerging small rope in detail.
At the beginning of the flux feeding process the small rope starts to emerge from the lower boundary, as shown by the brown curve in Fig. 1b. The arcades below the major rope in the initial state are then pushed upward by the emerging small rope, so that the top boundary of these arcades reaches the lower boundary of the major rope (as shown by the green curves in Fig. 1c), and then these arcades reconnect with the magnetic fields of the major rope (Fig. 1d). A comparison between Figs. 1b and d shows that the outer section of the major rope is peeled off by the reconnection. After all the arcades have reconnected with the major rope, the small rope itself also reaches the major rope. During the first half (0 ∼ τE/2 = 30τA) of the flux feeding process, what emerges from the lower base is the top half of the small rope, in which the direction of Bx is opposite to that in the bottom half of the major rope. Thus, a horizontal current sheet forms at the interface between the two ropes. This current sheet can be clearly recognized in Fig. 2a, where the distribution of the current density is shown in blue. As a result of the reconnection within this current sheet, the two ropes gradually merge together, as shown by the magnetic field lines plotted in red and pink in Figs. 1e and f. The magnetic cancelation occurring in this current sheet should decrease the local magnetic pressure below the major rope, so that the major rope gradually descends during the early period of the flux feeding process (see Figs. 1f and g), also see the height profile of the major rope axis plotted in Fig. 2i). During the second half of the flux feeding process (30τA ∼ 60τA), however, it is the bottom half of the small rope that emerges from the lower base. The direction of the magnetic field lines in the small rope is now the same as that in the bottom half of the major rope; there should be no reconnection any more. Therefore, after t = 30τA, magnetic flux is injected from the lower base so that the local magnetic pressure accumulates below the major rope, which causes the major rope to gradually rise, as shown in Figs. 1g and h.
![]() |
Fig. 2. Eruptive process of the resultant flux rope after the flux feeding process with CE = 2.23. Panels a–h: illustrate the magnetic configurations at different times, and the blue color depicts the distribution of the current density. Panel i: evolution of the height of the rope axis, and the times of panels a–h are marked by the vertical dotted lines in panel i. The green curves in panels a–h mark the outer boundaries of the flux ropes and the plasmoids. |
The flux feeding process ends at t = 60τA. The configuration of the resultant flux rope after flux feeding is illustrated in Fig. 1i: the rope sticks to the photosphere, without any arcade below the rope.
3.2. Eruption caused by flux feeding
Further evolution of the resultant major rope after flux feeding indicates that the major rope is caused to erupt by this flux feeding process with CE = 2.23. As shown in Fig. 2, the major rope keeps rising after flux feeding. Due to the line-tying effect from the photosphere the rope does not separate from the lower base immediately; its lower boundary sticks to the photosphere for a certain period (Fig. 2b). As the rope rises, the lower section of the flux rope and the ambient background field are stretched. Eventually, as shown in Figs. 2c and d, the rope is fully detached from the lower base, and magnetic reconnection occurs within the vertical current sheet formed below the rope. The reconnection results in closed arcades below the rope, which can be clearly seen in the lower section of Fig. 2d. Moreover, as shown in Fig. 2e, a plasmoid appears within the current sheet below the rope, and then rises and merges with the major rope (Fig. 2f), which is similar to the observations and simulations in Gou et al. (2019). This process is followed by another similar one, as shown in Figs. 2g and h. These plasmoids should result from the tearing mode instability within the current sheet, which have been investigated in detail in many previous studies (e.g., Bárta et al. 2008; Shen et al. 2011; Ni et al. 2012; Huang et al. 2017).
The height of the major rope axis versus time is plotted in Fig. 2i. As discussed in Sect. 3.1, due to the magnetic cancellation during the early period of the flux feeding process, the major flux rope gradually descends, which is accompanied by the contraction of the background arcades above the major rope (see Figs. 1f and g). This kind of downward movement before the onset of the eruption is often considered as a signature of “coronal implosion” (Hudson 2000), which has been reported in many previous observational studies (e.g., Ji et al. 2004; Veronig et al. 2006; Joshi et al. 2009; Liu et al. 2009; Li et al. 2017). For example, Liu et al. (2009) observed the large-scale contraction of the overlying coronal loops before the onset of the flare, which was attributed to the prolonged pre-heating phase dominated by coronal thermal emissions. Our simulation results suggest another scenario of coronal implosions: internal magnetic cancellation during the early period of the flux feeding process in HFT configurations decreases the local magnetic pressure below the major flux rope, resulting in the downward movement of the rope before the onset of its eruption. For comparison, the height profile of the flux-rope eruption caused by flux feeding in BPS configurations is monotonic, without any downward movement before the eruption (see Paper I).
It has been suggested in many previous studies that solar eruptions could not only be initiated by the gradual variations in photospheric magnetic fields, but in turn might also cause rapid magnetic field changes in the photosphere (Wang et al. 1994; Liu et al. 2005; Wang & Liu 2015; Toriumi & Wang 2019). The issue about this kind of photospheric response to the coronal eruptions is graphically called the “tail wags the dog” problem (Aulanier 2016). This phenomenon has been observed in many previous studies (e.g., Wang et al. 2012; Liu et al. 2012; Petrie 2012; Sun et al. 2017; Castellanos Durán et al. 2018). For example, by measuring the various magnetic parameters of a compact region at the central polarity inversion line (PIL) of an active region, Sun et al. (2017) found that the horizontal component of the photospheric magnetic field in this compact region obviously increase after the flare. In our simulation results, Fig. 3 shows the photospheric distribution of the horizontal magnetic field, Bx, at the central region right below the major rope (we note that Bz is zero in the background field). In the initial state, as shown in Fig. 3e, Bx at the photosphere is negative, consistent with the arcades below the rope in the initial state. During the eruptive process of the resultant major rope after flux feeding, closed arcades form right below the rising rope, which also results in the negative photospheric Bx, as shown in Fig. 3f. Moreover, as the rope rises, more and more field lines are closed by reconnection and pile up above the PIL. As a result of the reconnection-driven contraction of the flare loops, the horizontal photospheric magnetic field Bx increases accordingly (see Figs. 3f–h), which is consistent with the simulation results in Barczynski et al. (2019). Comparing the photospheric Bx at t = 200τA with that in the initial state, the average strength of the photospheric Bx after the eruption increases from 6.27 G to 8.50 G.
![]() |
Fig. 3. Variation of the photospheric magnetic field after the eruption. Panels a–d: illustrate the magnetic configurations at the central region above the PIL, and panels e–h: corresponding horizontal component of the photospheric magnetic fields, Bx. The times shown in panels b–d are in units of τA. |
Moreover, as shown in Fig. 3, the magnetic field strength is on the order of 10 G, indicating that coronal eruptions originating in quiescent regions could also result in the tail wags the dog phenomenon; in other words, coronal eruptions originating in quiescent regions could also cause variations of the photospheric magnetic fields.
4. Initiation of eruptions
In order to investigate how the eruption is initiated, we first investigate the influence of flux feeding on the major flux rope. The magnetic properties of a flux rope could be characterized by its axial magnetic flux, Φz, and its poloidal magnetic flux per unit length along the z-direction, Φp. For the initial major rope shown in Fig. 1a, the axial flux Φz0 = 4.366 × 1019 Mx, and the poloidal flux Φp0 = 1.186 × 1010 Mx cm−1. Here we select the flux rope at t = 80τA (Fig. 1i) as the resultant rope after flux feeding, whose axial flux is Φz = 7.115 × 1019 Mx, and poloidal flux is Φp = 1.116 × 1010 Mx cm−1. The increased axial flux is injected by the flux feeding process, which is similar to the BPS cases investigated in Paper I. Different from the result in the BPS cases, however, the poloidal flux is reduced by 0.070 × 1010 Mx cm−1 after flux feeding.
Since the distribution of the coronal and the photospheric magnetic fields play a dominant role in triggering coronal flux-rope eruptions (Yeates 2014; Yang et al. 2018; Thalmann et al. 2019; Xing et al. 2020), the influence of flux feeding processes on the major rope system should be sensitive to the scale of the magnetic field strength in the small rope. The magnetic parameters of the resultant flux ropes after the flux feeding processes with different CE are listed in Table 1. Larger CE implies stronger magnetic field strength in the small emerging rope, such that more axial flux is injected into the major rope. The deduced poloidal fluxes in the different cases, however, are very close to each other, ΔΦp = 0.069 ∼ 0.070 × 1010 Mx cm−1, which is almost the same as the magnetic flux of the arcades below the rope in the initial state, Φa. This indicates that the reduced poloidal flux should be caused by the interaction between the major rope and the arcades below the rope with the HFT configuration (introduced in Sect. 3.1): as pushed by the rising small flux rope, these arcades reconnect with the major rope so that the outer section of the major rope is peeled off.
Parameters of the resultant flux ropes after the flux feeding processes with different CE.
In the cases with differrent CE the further evolution of the resultant rope also varies, as shown in the last column in Table 1. The eruptive and the non-eruptive cases are well separated: the resultant rope erupts if CE is no smaller than 2.23, and the axial fluxes of the resultant ropes in all the eruptive cases are larger than those in the non-eruptive cases. This indicates that a critical axial magnetic flux of the resultant rope should exist, and the onset of its eruption should be initiated by the upward catastrophe associated with this critical axial flux (e.g., Su et al. 2011; Zhang et al. 2016). The value of the critical axial flux is of the order of 7.1 × 1019 Mx. The simulation results of the critical non-eruptive case with CE = 2.22 are illustrated in Fig. 4: the resultant rope does not keep rising after flux feeding (Figs. 4b and c), but falls back to the photosphere (Figs. 4d and e), resulting in an equilibrium state where the rope sticks to the photosphere (Fig. 4f). Therefore, if the injected axial flux is not large enough so that the critical axial flux of the major rope is not reached, not only is the flux feeding process unable to cause the major rope to erupt, but also the major rope system is transformed from an HFT configuration to a stable BPS configuration. This result might also suggest a possible scenario for confined flares in HFT configurations.
![]() |
Fig. 4. Simulation results of the critical non-eruptive case with CE = 2.22. The meanings of the symbols are the same as those in Fig. 2. |
As demonstrated in Fig. 2, obvious magnetic reconnection occurs in the vertical current sheet below the rising major flux rope. Here we simulate a specific case to distinguish between the role that reconnection plays and that upward catastrophe plays in initiating the flux-rope eruption: CE is also 2.23, the same as the eruptive case shown in Fig. 2, but the magnetic reconnection in the current sheet below the rope is prohibited. To achieve this, simulating procedures similar to those used in, for example, Hu (2001) and Zhang et al. (2016), are used here: after t = 80τA, first, the resistivity is adjusted to zero, and second, the flux function ψ is reassigned along the entire vertical current sheet (if it exists) below the major rope with the initial value ψc = ψi(x = 0, y = 0) at each time step, so that ψ along the current sheet is invariant. With these procedures, both numerical and physical magnetic reconnections below the major rope are prohibited.
The simulation results for this specific case are illustrated in Fig. 5. After flux feeding, the upward catastrophe occurs and the resultant rope rises for a while (Figs. 5a–c), but then stops rising, and eventually, the rope does not erupt but is suspended in the corona (Figs. 5d–f). Without magnetic reconnection, the closed background arcades above the major rope are always rooted at the photosphere. These closed arcades constrain the major rope, and prevent the rope from rising further. For comparison, in the case with reconnection demonstrated in Fig. 2, background arcades above the rope reconnect at the vertical current sheet below the rope, so that the “tether” constraining the rising major rope is removed. As suggested by Green et al. (2018), different kinds of eruptive mechanisms can be classified as “trigger” or “driver” according to the predominant role in the eruption process. In our simulation results, the flux feeding process causes the upward catastrophe, which triggers the flux-rope eruption. Magnetic reconnection and the further evolution of the upward catastrophe act as the drivers of the flux-rope eruption.
![]() |
Fig. 5. Simulation results of the critical eruptive case with CE = 2.23, in which the magnetic reconnection below the major rope is prohibited. The meanings of the symbols are the same as in Fig. 2. The red dotted line in panel g is the corresponding eruptive profile with reconnection, which is the same as that in Fig. 2i. |
5. Discussion and conclusion
In this paper we investigated the influence of flux feeding on coronal magnetic flux ropes with an HFT configuration, especially how flux feeding causes the flux rope to erupt. During the flux feeding process, the small emerging flux rope and the arcades below the major rope with the HFT configuration reconnect with the pre-existing major flux rope, so that axial magnetic flux is injected into the major rope, whereas the poloidal magnetic flux of the major rope is reduced. Flux feeding is able to cause the major rope to erupt, provided that the amount of axial flux injected by flux feeding is large enough for the total axial flux of the major rope to reach its critical value, which is similar to that in the BPS cases investigated in Paper I. During the early period of the flux feeding process, the major rope gradually descends due to the magnetic cancellation, demonstrating a downward movement of the flux rope right before its eruption, which suggests an alternative theoretical scenario for coronal implosions. Moreover, our simulation reproduces the photospheric magnetic response to the eruption of the coronal flux rope: the horizontal component of the photospheric magnetic field increases around the PIL after the eruption, as is often observed in the wake of major solar eruptions.
Although the onset of the eruption in both the BPS and the HFT cases is dominated by the upward catastrophe associated with the critical axial flux, the influence of flux feeding on the flux ropes with an HFT configuration is different from those with a BPS configuration. With BPS the major flux rope sticks to the photosphere, so that the small rope directly interacts with the major rope from the very beginning of the flux feeding process, and the major rope only reconnects with the small emerging rope during the flux feeding process; the non-eruptive resultant rope after flux feeding is still with a BPS configuration. For HFT cases, there is an additional interaction between the major rope and the arcades below it in the initial state. Consequently, not only is the poloidal flux of the major rope reduced due to the reconnection with these arcades, but the HFT configuration also collapses after flux feeding. The flux feeding in HFT configurations results in a dramatic change of topology, with the flux rope either erupting or falling back to a BPS configuration. Moreover, previous theoretical studies found that a flux rope could erupt if either its axial or poloidal magnetic flux exceeds the corresponding critical value. For the BPS cases, each episode of flux feeding injects axial flux into the flux rope while maintaining its poloidal flux, therefore pushing the rope system one step closer toward the eruption, no matter how small the amount of the injected axial flux is. For the HFT cases, however, if the injected axial flux is too little, the reduction in poloidal flux could push the rope system even farther away from the onset of the eruption; in addition, the rope system may also evolve from a metastable HFT configuration to a stable BPS configuration after flux feeding. Therefore, the flux feeding process in HFT configurations does not always favor the eruption of the major rope; it could even stabilize the rope system.
As introduced in Sect. 4, flux feeding cases with different CE are simulated in this paper. The emerging velocity of the small rope, however, is unchanged in different cases. Since the emerging velocity might also have some influence on the evolution of the resultant rope, we will try to simulate flux feeding cases with various emerging velocities in our future work, so as to investigate the role that the kinetic properties of flux feeding processes plays in initiating coronal flux-rope eruptions.
Acknowledgments
We appreciate prof. Jie Zhang for his insightful suggestions and comments. We also thank Helenka Kinnan for her valuable editorial comments. This research is supported by the National Natural Science Foundation of China (NSFC 41804161, 41774178, 41761134088, 41774150, 41842037, and 41574165), the Strategic Priority Program of CAS (XDB41000000 and XDA15017300), and the fundamental research funds for the central universities. We acknowledge for the support from National Space Science Data Center, National Science and Technology Infrastructure of China (www.nssdc.ac.cn).
References
- Antiochos, S. K., DeVore, C. R., & Klimchuk, J. A. 1999, ApJ, 510, 485 [Google Scholar]
- Archontis, V., & Hood, A. W. 2008, ApJ, 674, L113 [NASA ADS] [CrossRef] [Google Scholar]
- Aulanier, G. 2016, Nat. Phys., 12, 998 [Google Scholar]
- Aulanier, G., Pariat, E., & Démoulin, P. 2005, A&A, 444, 961 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Aulanier, G., Török, T., Démoulin, P., & DeLuca, E. E. 2010, ApJ, 708, 314 [Google Scholar]
- Barczynski, K., Aulanier, G., Masson, S., & Wheatland, M. S. 2019, ApJ, 877, 67 [NASA ADS] [CrossRef] [Google Scholar]
- Bárta, M., Vršnak, B., & Karlický, M. 2008, A&A, 477, 649 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Benz, A. O. 2008, Liv. Rev. Sol. Phys., 5, 1 [Google Scholar]
- Bobra, M. G., van Ballegooijen, A. A., & DeLuca, E. E. 2008, ApJ, 672, 1209 [NASA ADS] [CrossRef] [Google Scholar]
- Castellanos Durán, J. S., Kleint, L., & Calvo-Mozo, B. 2018, ApJ, 852, 25 [NASA ADS] [CrossRef] [Google Scholar]
- Chen, P. F. 2011, Liv. Rev. Sol. Phys., 8, 1 [Google Scholar]
- Chen, P. F., & Shibata, K. 2000, ApJ, 545, 524 [Google Scholar]
- Chen, Y., Hu, Y. Q., & Sun, S. J. 2007, ApJ, 665, 1421 [Google Scholar]
- Cheng, X., Ding, M. D., Zhang, J., et al. 2014, ApJ, 789, 93 [NASA ADS] [CrossRef] [Google Scholar]
- Chintzoglou, G., Vourlidas, A., Savcheva, A., et al. 2017, ApJ, 843, 93 [NASA ADS] [CrossRef] [Google Scholar]
- Démoulin, P., & Aulanier, G. 2010, ApJ, 718, 1388 [NASA ADS] [CrossRef] [Google Scholar]
- Fan, Y., & Gibson, S. E. 2007, ApJ, 668, 1232 [Google Scholar]
- Filippov, B. 2013, Sol. Phys., 283, 401 [Google Scholar]
- Forbes, T. G., & Isenberg, P. A. 1991, ApJ, 373, 294 [Google Scholar]
- Galsgaard, K., Titov, V. S., & Neukirch, T. 2003, ApJ, 595, 506 [NASA ADS] [CrossRef] [Google Scholar]
- Gibson, S. E., & Fan, Y. 2006, ApJ, 637, L65 [NASA ADS] [CrossRef] [Google Scholar]
- Gopalswamy, N., Akiyama, S., Yashiro, S., & Xie, H. 2018, J. Atm. Solar-Terr. Phys., 180, 35 [Google Scholar]
- Gou, T., Liu, R., Kliem, B., Wang, Y., & Veronig, A. M. 2019, Sci. Adv., 5, 7004 [Google Scholar]
- Green, L. M., Török, T., Vršnak, B., Manchester, W., & Veronig, A. 2018, Space Sci. Rev., 214, 46 [Google Scholar]
- Guo, Y., Ding, M. D., Schmieder, B., et al. 2010, ApJ, 725, L38 [NASA ADS] [CrossRef] [Google Scholar]
- Hu, Y. Q. 1989, J. Comput. Phys., 84, 441 [NASA ADS] [CrossRef] [Google Scholar]
- Hu, Y. Q. 2001, Sol. Phys., 200, 115 [Google Scholar]
- Hu, Y. Q., Li, X., & Ai, G. X. 1995, ApJ, 451, 843 [Google Scholar]
- Huang, Y.-M., Comisso, L., & Bhattacharjee, A. 2017, ApJ, 849, 75 [Google Scholar]
- Hudson, H. S. 2000, ApJ, 531, L75 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
- Inoue, S., Hayashi, K., Magara, T., Choe, G. S., & Park, Y. D. 2015, ApJ, 803, 73 [NASA ADS] [CrossRef] [Google Scholar]
- Isenberg, P. A., Forbes, T. G., & Demoulin, P. 1993, ApJ, 417, 368 [Google Scholar]
- Ji, H., Wang, H., Goode, P. R., Jiang, Y., & Yurchyshyn, V. 2004, ApJ, 607, L55 [NASA ADS] [CrossRef] [Google Scholar]
- Jiang, C., Feng, X., & Hu, Q. 2018, ApJ, 866, 96 [Google Scholar]
- Joshi, B., Veronig, A., Cho, K. S., et al. 2009, ApJ, 706, 1438 [NASA ADS] [CrossRef] [Google Scholar]
- Kliem, B., & Török, T. 2006, Phys. Rev. Lett., 96 [CrossRef] [Google Scholar]
- Kliem, B., Lin, J., Forbes, T. G., Priest, E. R., & Török, T. 2014, ApJ, 789, 46 [NASA ADS] [CrossRef] [Google Scholar]
- Li, G., & Hu, Y. 2002, Sci. China A: Math., 45, 65 [Google Scholar]
- Li, Y., Sun, X., Ding, M. D., Qiu, J., & Priest, E. R. 2017, ApJ, 835, 190 [NASA ADS] [CrossRef] [Google Scholar]
- Lin, J., Forbes, T. G., & Isenberg, P. A. 2001, J. Geophys. Res., 106, 25053 [NASA ADS] [CrossRef] [Google Scholar]
- Lin, J., Soon, W., & Baliunas, S. L. 2003, New Astron. Rev., 47, 53 [NASA ADS] [CrossRef] [Google Scholar]
- Liu, R. 2020, Res. Astron. Astrophys., 20, 165 [Google Scholar]
- Liu, C., Deng, N., Liu, Y., et al. 2005, ApJ, 622, 722 [NASA ADS] [CrossRef] [Google Scholar]
- Liu, R., Wang, H., & Alexander, D. 2009, ApJ, 696, 121 [NASA ADS] [CrossRef] [Google Scholar]
- Liu, C., Deng, N., Liu, R., et al. 2012, ApJ, 745, L4 [NASA ADS] [CrossRef] [Google Scholar]
- Liu, L., Cheng, X., Wang, Y., et al. 2018, ApJ, 867, L5 [NASA ADS] [CrossRef] [Google Scholar]
- Longcope, D. W., & Forbes, T. G. 2014, Sol. Phys., 289, 2091 [Google Scholar]
- Lugaz, N., Farrugia, C. J., Winslow, R. M., et al. 2017, ApJ, 848, 75 [Google Scholar]
- Moore, R. L., Sterling, A. C., Hudson, H. S., & Lemen, J. R. 2001, ApJ, 552, 833 [Google Scholar]
- Ni, L., Roussev, I. I., Lin, J., & Ziegler, U. 2012, ApJ, 758, 20 [NASA ADS] [CrossRef] [Google Scholar]
- Parenti, S. 2014, Liv. Rev. Sol. Phys., 11, 1 [Google Scholar]
- Petrie, G. J. D. 2012, ApJ, 759, 50 [NASA ADS] [CrossRef] [Google Scholar]
- Romano, P., Contarino, L., & Zuccarello, F. 2003, Sol. Phys., 214, 313 [NASA ADS] [CrossRef] [Google Scholar]
- Savcheva, A. S., van Ballegooijen, A. A., & DeLuca, E. E. 2012, ApJ, 744, 78 [NASA ADS] [CrossRef] [Google Scholar]
- Shen, C., Lin, J., & Murphy, N. A. 2011, ApJ, 737, 14 [NASA ADS] [CrossRef] [Google Scholar]
- Shen, F., Shen, C., Zhang, J., et al. 2014, J. Geophys. Res.: Space Phys., 119, 7128 [Google Scholar]
- Sterling, A. C., & Moore, R. L. 2004, ApJ, 602, 1024 [Google Scholar]
- Su, Y., Surges, V., van Ballegooijen, A., DeLuca, E., & Golub, L. 2011, ApJ, 734, 53 [Google Scholar]
- Sun, X., Hoeksema, J. T., Liu, Y., Kazachenko, M., & Chen, R. 2017, ApJ, 839, 67 [Google Scholar]
- Thalmann, J. K., Linan, L., Pariat, E., & Valori, G. 2019, ApJ, 880, L6 [NASA ADS] [CrossRef] [Google Scholar]
- Titov, V. S., & Démoulin, P. 1999, A&A, 351, 707 [NASA ADS] [Google Scholar]
- Titov, V. S., Priest, E. R., & Demoulin, P. 1993, A&A, 276, 564 [NASA ADS] [Google Scholar]
- Titov, V. S., Galsgaard, K., & Neukirch, T. 2003, ApJ, 582, 1172 [Google Scholar]
- Toriumi, S., & Wang, H. 2019, Liv. Rev. Sol. Phys., 16, 3 [Google Scholar]
- Török, T., & Kliem, B. 2003, A&A, 406, 1043 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Švestka, Z. 2001, Space Sci. Rev., 95, 135 [NASA ADS] [CrossRef] [Google Scholar]
- van Driel-Gesztelyi, L., & Green, L. M. 2015, Liv. Rev. Sol. Phys., 12, 1 [NASA ADS] [CrossRef] [Google Scholar]
- Veronig, A. M., Karlický, M., Vršnak, B., et al. 2006, A&A, 446, 675 [Google Scholar]
- Vršnak, B., Sudar, D., & Ruždjak, D. 2005, A&A, 435, 1149 [Google Scholar]
- Wang, H., & Liu, C. 2015, Res. Astron. Astrophys., 15, 145 [NASA ADS] [CrossRef] [Google Scholar]
- Wang, H., Ewell, M. W., Jr., Zirin, H., & Ai, G. 1994, ApJ, 424, 436 [NASA ADS] [CrossRef] [Google Scholar]
- Wang, S., Liu, C., & Wang, H. 2012, ApJ, 757, L5 [NASA ADS] [CrossRef] [Google Scholar]
- Xing, C., Cheng, X., Qiu, J., et al. 2020, ApJ, 889, 125 [Google Scholar]
- Yan, X., Xue, Z., Cheng, X., et al. 2020, ApJ, 889, 106 [NASA ADS] [CrossRef] [Google Scholar]
- Yang, S., Büchner, J., Skála, J., & Zhang, H. 2018, A&A, 613, A27 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Yeates, A. R. 2014, Sol. Phys., 289, 631 [NASA ADS] [CrossRef] [Google Scholar]
- Zhang, J., Dere, K. P., Howard, R. A., Kundu, M. R., & White, S. M. 2001, ApJ, 559, 452 [Google Scholar]
- Zhang, Q., Liu, R., Wang, Y., et al. 2014, ApJ, 789, 133 [Google Scholar]
- Zhang, Q., Wang, Y., Hu, Y., & Liu, R. 2016, ApJ, 825, 109 [Google Scholar]
- Zhang, Q., Wang, Y., Hu, Y., Liu, R., & Liu, J. 2017, ApJ, 835, 211 [Google Scholar]
- Zhang, Q., Wang, Y., Liu, R., et al. 2020, ApJ, 898, L12 [Google Scholar]
- Zhuang, B., Hu, Y., Wang, Y., et al. 2018, J. Geophys. Res.: Space Phys., 123, 2513 [Google Scholar]
Appendix A: Basic equations and initial preparations
Combining the 2.5-dimensional form of the magnetic field (Eq. (1)), we write the MHD equations in the non-dimensional form as
where
Here ρ and T denote the density and the temperature; vx, vy, vz represent the x-, y-, and z-component of the velocity, respectively; γ = 5/3 is the polytropic index in our simulation; g is the normalized gravity; η is the resistivity. Here is the characteristic ratio of the gas pressure to the magnetic pressure, where ρ0 = 3.34 × 10−13 kg m−3, T0 = 106 K, L0 = 107 m, and ψ0 = 3.73 × 103 Wb m−1 are the characteristic values of density, temperature, length, and magnetic flux function, respectively, which are also the calculating units in the simulation. For the other quantities the corresponding characteristic values are v0 = 128.57 km s−1, t0 = 77.8 s, B0 = 3.37 × 10−4 T, and g0 = 1.65 × 103 m s−2. The numerical domain in our simulation is 0 < x < 200 Mm, 0 < y < 300 Mm, which is discretized into 400 × 600 uniform meshes. At the left side of the domain (x = 0) the symmetric boundary condition is used. The radiation and the heat conduction in the energy equation are neglected.
In order to simulate flux feeding processes in the HFT cases, we must first construct a typical coronal flux-rope system with an HFT configuration. Here we select a partially open bipolar field as the background field, in which a negative and a positive surface magnetic charge located at the photosphere within −b < x < −a and a < x < b, respectively. The background field can then be obtained by the complex variable method (e.g., Hu et al. 1995; Zhang et al. 2017). The background magnetic field can be cast in the complex variable form
where ω = x + iy, and
Here a = 15 Mm, b = 25 Mm, and (y = yN = 34.5 Mm, x = 0) is the position of the neutral point of the partially open bipolar field. There is a neutral current sheet located at (x = 0, y ≥ yN). The magnetic flux function could then be calculated as
and the flux function at the lower base as
where ψc = πψ0; the flux function at the neutral point y = yN can be calculated as
Setting Bz equal to 0 in the background field, the configuration of the background field could then be obtained, as shown in Fig. A.1a. The initial corona is isothermal and static:
![]() |
Fig. A.1. Panel a: magnetic configuration of the background field. Panel b: interim state to construct the initial state from the background field. Panel c: obtained flux-rope system with an HFT configuration, which is the initial state of our simulation. |
Except during the emergence of the small rope, the lower boundary is fixed: the flux function ψ is fixed at ψi given by Eq. (A.10); Bz and the velocity is always 0; the density and the temperature are fixed at their initial values. Increment equivalent extrapolation is used at the right and top boundaries.
With the initial condition and the background configuration, Eqs. (A.1)–(A.5) are solved by the multi-step implicit scheme. First, we let a flux rope emerge from the lower base of the background field, as shown in Fig. A.1b. Multiplying Bz within the rope by a factor larger than 1, we increase the axial magnetic flux of the flux rope. The rope then rises, and magnetic reconnection occurs within the vertical current sheet formed below the rope, resulting in closed arcades below the rope. After that, we adjust the axial flux of the rope again, so that the rope stops rising. Eventually, we let the rope system relax to an equilibrium state with an HFT configuration, as shown in Fig. A.1c. This is the initial state of our simulation. It is noteworthy that the radius of the flux rope in our simulation is finite, so that the thin-rope approximation is not satisfied. Under this circumstance, the initial state could only be obtained by numerical procedures.
All Tables
Parameters of the resultant flux ropes after the flux feeding processes with different CE.
All Figures
![]() |
Fig. 1. Evolution of the magnetic configuration during a flux feeding process with CE = 2.23. Panel a: initial flux-rope system with an HFT configuration; the green curves indicate the boundary of the major flux rope and that of the arcades below the rope. The times indicated in panels b–i are in units of τA. Different colors depict the magnetic field lines corresponding to different values of ψ. |
In the text |
![]() |
Fig. 2. Eruptive process of the resultant flux rope after the flux feeding process with CE = 2.23. Panels a–h: illustrate the magnetic configurations at different times, and the blue color depicts the distribution of the current density. Panel i: evolution of the height of the rope axis, and the times of panels a–h are marked by the vertical dotted lines in panel i. The green curves in panels a–h mark the outer boundaries of the flux ropes and the plasmoids. |
In the text |
![]() |
Fig. 3. Variation of the photospheric magnetic field after the eruption. Panels a–d: illustrate the magnetic configurations at the central region above the PIL, and panels e–h: corresponding horizontal component of the photospheric magnetic fields, Bx. The times shown in panels b–d are in units of τA. |
In the text |
![]() |
Fig. 4. Simulation results of the critical non-eruptive case with CE = 2.22. The meanings of the symbols are the same as those in Fig. 2. |
In the text |
![]() |
Fig. 5. Simulation results of the critical eruptive case with CE = 2.23, in which the magnetic reconnection below the major rope is prohibited. The meanings of the symbols are the same as in Fig. 2. The red dotted line in panel g is the corresponding eruptive profile with reconnection, which is the same as that in Fig. 2i. |
In the text |
![]() |
Fig. A.1. Panel a: magnetic configuration of the background field. Panel b: interim state to construct the initial state from the background field. Panel c: obtained flux-rope system with an HFT configuration, which is the initial state of our simulation. |
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.