- 1Planetary Environmental and Astrobiological Research Laboratory (PEARL), School of Atmospheric Sciences, Sun Yat-sen University, Zhuhai, China
- 2Key Laboratory of Tropical Atmosphere-Ocean System, Sun Yat-sen University, Ministry of Education, Zhuhai, China
- 3CAS Center for Excellence in Comparative Planetology, Hefei, China
- 4Institute of Space Science and Applied Technology, Harbin Institute of Technology, Shenzhen, China
- 5CAS Key Laboratory of Geospace Environment, University of Science and Technology of China, Hefei, Anhui, China
Solar eruptions often show the rotation of filaments, which is a manifestation of the rotation of erupting magnetic flux rope (MFR). Such a rotation of MFR can be induced by either the torque exerted by a background shear-field component (which is an external cause) or the relaxation of the magnetic twist of the MFR (an internal cause). For a given chirality of the erupting field, both the external and internal drivers cause the same rotation direction. Therefore, it remains elusive from direct observations which mechanism yields the dominant contribution to the rotation. In this paper, we exploit a full MHD simulation of solar eruption by tether-cutting magnetic reconnection to study the mechanism of MFR rotation. In the simulation, the MFR’s height–rotation profile suggests that the force by the external shear-field component is a dominant contributor to the rotation. Furthermore, the torque analysis confirms that it is also the only factor in driving the counterclockwise rotation. On the contrary, the Lorentz torque inside the MFR makes a negative effect on this counterclockwise rotation.
1 Introduction
Solar eruptions are the most spectacular phenomena that occur in the solar corona, often (but not always) releasing coronal mass ejections (CMEs) into interplanetary space. A magnetic flux rope (MFR), characterized by a twisted and writhed topological structure of magnetic fields, is generally accepted as a basically magnetic configuration underlying the phenomenon of CMEs [1]. A southward orientation of its front-side magnetic field may cause severe geoeffectiveness of the space environment [2,3]. During its eruption and propagation, a rotational motion is also observed, which modulates the orientation of the MFR ([4–6], and references therein). Therefore, unveiling the physics that triggers the rotation of erupting MFR is a key issue in understanding the space weather and its terrestrial effects.
Several theories of how the flux rope rotates are competing in the last decades. For an MFR existing prior to an eruption, when its twist exceeds a threshold, a helical kink instability (KI [7]) may occur. Sigmoids observed in Soft-X-ray and rotating filaments in EUV channels are interpreted to be outlining a MFR undergoing KI [4]. Observationally, the rotation direction of a filament apex has a strong one-to-one relationship with the filament helicity (chirality): For positive/negative helicity, filaments (with sinistral/dextral chirality) rotate clockwise (CW)/counterclockwise (CCW) [5,6]. Numerical simulations [8] find this relationship to be consistent with the conversion of twist into writhe under the ideal MHD constraint of helicity conservation, suggesting the rotation mechanism as the relaxation of tension in the twisted field. Isenberg and Forbes [9] propose an alternative rotation mechanism: a Lorentz force exerted on the MFR legs from the misalignment between the MFR’s toroidal current and the external toroidal field drives the MFR rotating. It cannot be easily distinguished which mechanism is responsible for the rotation, because both of these two mechanisms writhe the flux rope axis in a similar manner from observations. For a given chirality, the rotation driven by the external shear-field component yields the same rotation direction as that driven by the helical kink instability [8,9].
To determine which mechanism dominates during the flux rope rotation, Kliem et al. [8] carried out a comparative study of these two mechanisms: Both the forces by twist and shear are potentially significant contributors to the rotation. Their contribution can be disentangled if the rotation and rise profiles in the observations are simultaneously compared with model calculations. The profile of MFR’s rotation vs. height in the simulation shows that the twist-driven rotation tends to saturate in the low corona (a height range up to several times of the distance between the footpoint of the MFR), while, the shear-driven rotation distributes across a larger height range. For parameters characteristic of filament source regions, the shear field is usually present and dominates the rotation in the corona even if the twist is sufficiently high to trigger the helical kink instability. For a considerable rotation angle (angles of order 90° and higher), this shear-driven component is required.
Recently, Jiang et al. [10] performed a fully three-dimensional (3D) magnetohydrodynamic (MHD) simulation demonstrating a runaway tether-cutting reconnection alone initiates a solar eruption. In the simulation, a MFR with dextral chirality is built up with a reverse-S shape through reconnection of the sheared arcades. While rising, the simulated MFR rotates CCW, which compares favorably with the observation rule of filament chirality and rotation direction [11].
Lacking the regular observations of the coronal magnetic field, we employ the same simulation result as presented in Zhou et al. [11] as the observation substitution to look into the structure evolution and force distribution of the MFR to explore the underlying mechanism of the MFR rotation. In this direction, we also calculate the torque associated with the Lorentz force. The numerical model is described in Section 2. A whole picture of the simulated eruption and a comparative analysis of height–rotation profile with the result of Kliem et al. [8] are given in Section 3, and a systematic survey of the global Lorentz torque is given in Section 4. Discussion and summary are presented in Section 5.
2 MHD simulation
Jiang et al. [10] performed a high-accuracy, fully 3D MHD simulation that examined a fundamental mechanism for the initiation of solar eruptions. They found that a current sheet can be slowly formed within a bipolar magnetic field as driven by quasi-static shearing motion at the photosphere. Then once magnetic reconnection sets in at the current sheet, it triggers and drives the eruption impulsively. In this process, an MFR is built up from the sheared arcade during the eruption, similar to the classic tether-cutting model. In this paper we study the evolution of the erupting MFR by using a similar simulation run like the one in Jiang et al. [10] but with a slightly lower resolution such that the evolution of the MFR can be followed longer. The exactly similar run is used by Zhou et al. [11] to demonstrate the CCW rotational motion of the MFR.
The simulation solves the full set of MHD equations including both solar gravity and plasma pressure. The initial condition consists of a bipolar potential magnetic field and a hydrostatic plasma (ρ = 2.3 × 10−15 g cm−3) stratified by solar gravity with typical coronal temperature (T = 106K). To energize the coronal field, we applied shearing flows along the PIL, which are implemented by rotating the two magnetic polarities at the bottom surface in the same CCW direction. The surface flow is stopped once an eruption is triggered (the time is referred to as t = 0 in the following). During the quasi-static evolution phase driven by the shearing motion, a current sheet is gradually built up. Because we used no explicit resistivity in the MHD model, magnetic reconnection is triggered when the current sheet is sufficiently thin with its width close to the grid resolution, owing to the implicit, grid-dependent numerical resistivity. The whole size of the computational volume extends as (−32, −32, 0)
3 Observation features of simulated eruption
As shown in Figure 1 (and accompanied Movie 1), the simulation demonstrates a typical solar eruption leading to a CME, as seen in many well-observed eruptions in different view angles [6,11]. The core structure of evolving magnetic field is transformed from the pre-eruptive sheared arcades to twisted MFR that subsequently expands, and then transforms into a semi-circular shape stretching outward its overlying field, finally displays a clear inverse γ shape, indicative of a writhing motion of the apex (Figure 1A). From the top view (Figure 1B), as the reconnection kicks in, an inversed S-shaped MFR is formed on-the-fly. After that, the MFR’s axis shows a significant CCW rotation, straightening, and even flipping in shape from reversed S to forward S.
FIGURE 1. Structural evolution of the eruption. (A) 3D perspective view of magnetic field lines colored by the force-free factor. Here the field lines are traced at fixed footpoints on the bottom surface, and they represent the core structure of the MFR and its ambient magnetic field. At t = 20ts, a red horizontal line gives the location of the slice in Figure 3A. (B) Top view of the structure is shown in (A). The unit of time is 105 s. Also, see Supplementary Movie S1 for a high-cadence evolution of the eruption process.
For the first time, this simulation reproduces the predictions for the relationship between the chirality of the erupting field, sigmoid shape, and projected MFR’s apex rotation in the tether-cutting model [5], i.e., for positive (negative) helicity the filament apex rotates CW (CCW) when viewed from above, and both the filament and related sigmoid take on a forward (reverse) S shape. This is consistent with that from observations [6,11].
Through this CCW rotation, the original inverse-S-shaped MFR is straightened to a straight loop, implying a reduction of writhe [12,13]. Furtherly, the initial increase of the flux rope twist by the transformation of writhe helicity is directly confirmed based on this simulation [11], which contradicts the expectation of the kink instability since this instability transforms twist into writhe. So instead we turn to an alternative rotation mechanism that is controlled by an external shear-field component [9]. Both of them yield the same rotation direction for a given chirality, still, we can tell which rotation mechanism dominates by comparing both the observed height–rotation profile with the corresponding curves from models [8].
A powerful tool—the height–rotation profile proposed by Kliem et al. [8] - can be used to diagnose the rotation mechanism of this simulated eruption. Here we use the rotation of MFR’s axis as a proxy of MFR’s rotation. The determination of this modeled MFR’s axis can be referred to Zhou et al. [11]. The rotation angle (Φ) is measured as the difference in angle between the tangent vector at the axis top, and its initial direction at t0, while the apex height (h) of MFR’s axis is normalized by the initial apex height (h0 = 2.75l0) of the axis (Figure 2A).
FIGURE 2. (A) The distribution of twist number (at t = 20ts) on the normal plane that is perpendicular to the tangent vector of MFR’s axis (the red curve) at the apex point (black filled diamond, h = 6.32h0). The bottom plane (z = 0) shows the Bz field in greyscale. (B) Comparison of flux rope rotation as a function of normalized apex height. The black-filled diamond marks out their relation at t = 20ts.
The profile of rotation vs. height then can be obtained in this simulated solar eruption (Figure 2B). At the beginning (t = 12), the MFR is formed with its axis apex at the initial height (h0) and the rotation angle is 0°. A distinct monotonic increase in the rotation angle could be found as the MFR rise up. When it reaches the height of 19h0, the axis rotates 85°. Apparently, there is no saturation of the rotation under the height of 19 h0. The parametric study by Kliem et al. [8] shows that twist-driven rotation alone tends to saturate at lower heights (
4 Force analysis
The height–rotation profile implies that the force by an external shear-field component makes a significant contributor to the rotation. This mechanism can easily be understood in a simplified scenario: a flux rope is settled in an ambient field. The misalignment between the flux rope’s toroidal current and the shear field yields antisymmetric sideways Lorentz forces on opposite side legs, forcing the flux rope to rotate. Gas pressure and gravity have a negligible effect on the magnetic structure in the corona. Thus, in this simulation, the only force to be considered is the Lorentz force F = ∫J ×BdV.
As all the structures are centrosymmetric with respect to the origin on each horizontal plane (z = z0), the MFR and background can be distinguished from the twist distribution (Figure 3A) and the vertical (z) component of the Lorentz torque M = ∫r × (J ×B)dV drives the structures on the horizontal plane to rotate. The positive/negative value in z direction means a CCW/CW rotation. In the simulation, the magnetic field B and current density J consist of those of both MFR and its background field. We assume that the background field B0 is a force-free magnetic arcade before the formation of the flux rope, which is thus taken as the pre-eruption field at t = 0 since the flux rope only begin to form after the eruption starts. Since the force-free field B0 is determined by the magnetic field at the bottom boundary (corresponding to the photospheric field), which remains almost unchanged during the fast eruption, we further assume that the background field B0 does not change with time. Then, to isolate shear-driven Lorentz torque from the total Lorentz torque, the magnetic field and current at time t can be instantly decomposed into B = B0 + Bf and J = J0 + Jf where Bf and Jf represent the MFR’s magnetic field and electric current respectively. J0 is current associated with the background magnetic field B0 at the initial time (t = 0) when the MFR has not yet formed, therefore it is considered to be constant with time. The original data from the simulation has been remapped to a uniform grid, then the total “Lorentz torque density” M integrated on a z-constant plane can be decomposed as:
FIGURE 3. Torque analysis at t =20ts. (A) Distribution of twist number on the horizontal slice at z = 12 l0. The red cross symbol is the center of rotation. (B) Distributions of each z-direction moment component on the z-constant plane (z = 12 l0). (C) The z-direction component of Lorentz torque induced by the external field (M0, grey line), MFR itself (M3, green line), their interaction force (M4, red line), and the total force (M, blue line) integrated on each horizontal slice. The black vertical dashed lines show the location of the left panel.
Of which, the components of “Lorentz torque density” consist of the background (M0), MFR itself (M3), and their interaction (M1 and M2). In particular, M2 refers to the torque induced by the current of the MFR with the background field, while M1 refers to that induced by the background current with the magnetic field of the MFR. To carry out the computation, the Lorentz torque is integrated on a z-constant plane. Due to the original data from the simulation being remapped to a uniform grid, Lorentz torque is the volume integration of the “Lorentz torque density”
A comparison between the contributions of these components is shown in Figure 3C. Of which, the blue profile shows that the total Lorentz torque (M) is always positive along the height, consistent with the observed CCW rotation. The Lorentz torque (M2, the orange dash line in Figure 3C), caused by the misalignment of the MFR-carried current and the sheared ambient field, is the dominant factor for the rotation, not only it has a greater influence (with a larger absolute value than the others), but also it is the only factor (with a positive value) that drives a CCW rotation. Compared to the other components in Figure 3C, the M0 part is almost zero (its mean value
5 Discussion and summary
With the aid of Jiang et al. [10]’s simulation, we investigated the rotation mechanism of erupting MFR formed in the runaway tether-cutting reconnection model. The relationship between chirality, sigmoid shape, and apex rotation in this simulation is in line with that in observations. As the MFR ascends during the eruption, its apex rotates CCW. The rotation angle of simulated MFR monotonically increases to 85° and shows no sign of saturation at a relatively large height
Further on, the force analysis confirms the driving force of the rotation originates solely from external sources. The decomposition of the z direction component of Lorentz torque shows clearly that the part induced by the MFR’s current and the background field, M2 = ∫r × (Jf ×B0)dV, is the only positive component to drive the MFR rotating CCW. This sideways Lorentz force continuously adjusts the loop current’s orientation to align with the external shear field. On the contrary, the z direction component of inner Lorentz torque, M3 = ∫r × (Jf ×Bf)dV, has a negative value, hindering the MFR’s CCW rotation. Therefore, the strength of the external shear field serves as the dominant factor determining the CCW rotation. Previous numerical simulations [8] suggest the KI yields the same rotation direction for the given chirality of the erupting field. Our simulation has indicated the occurrence of the chirality-corresponding rotation does not necessarily reflect a relaxation of tension in the twisted field, and it quite suggests the opposite process: an enhancement of tension in the twisted field. This has also been validated by our recent study [11]: a quantitative measurement in this simulation clearly shows conversion of writhe into twist during rotation.
It is worth noting that our research is based on the MHD simulation of CME initiation by tether-cutting reconnection in a single magnetic arcade rather than a pre-existing unstable MFR. Therefore, although the inference from this model is consistent with the observed sigmoid–chirality relationship, it does not rule out other possibilities, because the initial magnetic topology being that of a sheared arcade is still intensely debated [14]. Moreover, there also exist some other mechanisms, e.g., magnetic reconnection of rope field lines with the ambient field [15–18], asymmetric deflection of the rising flux rope induced by adjacent coronal holes [19] and/or local alignment with the heliospheric current sheet [20], can cause the rotation of a MFR. All these theories are needed to be tested with a measurement of the coronal magnetic field in future works, to understand fully the mechanism of filament rotation during solar eruptions.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author contributions
ZZ leads this work, CJ performs the numerical simulations and all contribute to the study.
Funding
This work is supported by the B-type Strategic Priority Program XDB41000000 funded by the Chinese Academy of Sciences. The authors also acknowledge support from the National Natural Science Foundation of China (NSFC 42004142,42274203), Open Research Program of CAS Key Laboratory of Geospace Environment, Science and Technology Project 202102021019 in Guangzhou, the Fundamental Research Funds for the Central Universities (grant No. HIT.BRETIV.201901), and Shenzhen Technology Project JCYJ20190806142609035.
Acknowledgments
The authors wish to express their special thanks to the referee for suggestions and comments which led to the improvement of the paper.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphy.2023.1119637/full#supplementary-material
References
1. Vourlidas A, Lynch BJ, Howard RA, Li Y. How many CMEs have flux ropes? Deciphering the signatures of shocks, flux ropes, and prominences in coronagraph observations of CMEs. Solar Phys (2013) 284:179–201. doi:10.1007/s11207-012-0084-8
2. Gosling JT. The solar flare myth. J Geophys Res Space Phys (1993) 98:18937–50. doi:10.1029/93JA01896
3. Webb DF, Forbes TG, Aurass H, Chen J, Martens P, Rompolt B, et al. Material ejection. Solar Phys (1994) 153:73–89. doi:10.1007/BF00712493
4. Gibson SE, Fan Y, Török T, Kliem B. The evolving sigmoid: Evidence for magnetic flux ropes in the corona before, during, and after CMES. Space Sci Rev (2006) 124:131–44. doi:10.1007/s11214-006-9101-2
5. Green LM, Kliem B, Török T, van Driel-Gesztelyi L, Attrill GDR. Transient coronal sigmoids and rotating erupting flux ropes. Sol Phys (2007) 246:365–91. doi:10.1007/s11207-007-9061-z
6. Zhou Z, Liu R, Cheng X, Jiang C, Wang Y, Liu L, et al. The relationship between chirality, sense of rotation, and hemispheric preference of solar eruptive filaments. Astrophysical J (2020) 891:180. doi:10.3847/1538-4357/ab7666
7. Török T, Kliem B, Titov VS. Ideal kink instability of a magnetic loop equilibrium. Astron Astrophysics (2004) 413:L27–L30. doi:10.1051/0004-6361:20031691
8. Kliem B, Török T, Thompson WT. A parametric study of erupting flux rope rotation. Modeling the “cartwheel CME” on 9 april 2008. Solar Phys (2012) 281:137–66. doi:10.1007/s11207-012-9990-z
9. Isenberg PA, Forbes TG. A three-dimensional line-tied magnetic field model for solar eruptions. Astrophysical J (2007) 670:1453–66. doi:10.1086/522025
10. Jiang C, Feng X, Liu R, Yan X, Hu Q, Moore RL, et al. A fundamental mechanism of solar eruption initiation. Nat Astron (2021) 5:1126–38. doi:10.1038/s41550-021-01414-z
11. Zhou Z, Jiang C, Liu R, Wang Y, Liu L, Cui J. The rotation of magnetic flux ropes formed during solar eruption. Astrophysical J Lett (2022) 927:L14. doi:10.3847/2041-8213/ac5740
12. Török T, Berger MA, Kliem B. The writhe of helical structures in the solar corona. Astron Astrophysics (2010) 516:A49. doi:10.1051/0004-6361/200913578
13. Liu R, Kliem B, Török T, Liu C, Titov VS, Lionello R, et al. Slow rise and partial eruption of a double-decker filament. I. Observations and interpretation. Astrophysical J (2012) 756:59. doi:10.1088/0004-637X/756/1/59
14. Patsourakos S, Vourlidas A, Török T, Kliem B, Antiochos SK, Archontis V, et al. Decoding the pre-eruptive magnetic field configurations of coronal mass ejections. Space Sci Rev (2020) 216:131. doi:10.1007/s11214-020-00757-9
15. Jacobs C, Roussev II, Lugaz N, Poedts S. The internal structure of coronal mass ejections: Are all regular magnetic clouds flux ropes? Astrophysical J Lett (2009) 695:L171–5. doi:10.1088/0004-637X/695/2/L171
16. Cohen O, Attrill GDR, Schwadron NA, Crooker NU, Owens MJ, Downs C, et al. Numerical simulation of the 12 May 1997 CME Event: The role of magnetic reconnection. J Geophys Res (Space Physics) (2010) 115:A10104. doi:10.1029/2010JA015464
17. Thompson WT. Strong rotation of an erupting quiescent polar crown prominence. J Atmos Solar-Terrestrial Phys (2011) 73:1138–47. doi:10.1016/j.jastp.2010.07.005
18. Vourlidas A, Colaninno R, Nieves-Chinchilla T, Stenborg G. The first observation of a rapidly rotating coronal mass ejection in the middle corona. Astrophysical J Lett (2011) 733:L23. doi:10.1088/2041-8205/733/2/L23
19. Panasenco O, Martin S, Joshi AD, Srivastava N. Rolling motion in erupting prominences observed by STEREO. J Atmos Solar-Terrestrial Phys (2011) 73:1129–37. doi:10.1016/j.jastp.2010.09.010
Keywords: magnetic fields, magnetohydrodynamics (MHD), methods: numerical, Sun: corona, Sun: flares
Citation: Zhou Z, Jiang C, Yu X, Wang Y, Hao Y and Cui J (2023) The mechanism of magnetic flux rope rotation during solar eruption. Front. Phys. 11:1119637. doi: 10.3389/fphy.2023.1119637
Received: 09 December 2022; Accepted: 16 February 2023;
Published: 23 February 2023.
Edited by:
Abhishek Kumar Srivastava, Indian Institute of Technology (BHU), IndiaReviewed by:
David MacTaggart, University of Glasgow, United KingdomSanjay Kumar, Patna University, India
Copyright © 2023 Zhou, Jiang, Yu, Wang, Hao and Cui. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Zhenjun Zhou, emhvdXpoajdAbWFpbC5zeXN1LmVkdS5jbg==; Chaowei Jiang, Y2hhb3dlaUBoaXQuZWR1LmNu