- 1Solar Physics and Space Plasma Research Center, School of Mathematics and Statistics, The University of Sheffield, Sheffield, United Kingdom
- 2CAS Key Laboratory of Geospace Environment, Department of Geophysics and Planetary Sciences, University of Science and Technology of China, Hefei, China
- 3CAS Center for the Excellence in Comparative Planetology, Hefei, China
- 4Department of Astronomy, Eötvös Loránd University, Budapest, Hungary
Highly twisted magnetic flux ropes, with finite length, are subject to kink instabilities, and could lead to a number of eruptive phenomena in the solar atmosphere, including flares, coronal mass ejections (CMEs) and coronal jets. The kink instability threshold, which is the maximum twist a kink-stable magnetic flux rope could contain, has been widely studied in analytical models and numerical simulations, but still needs to be examined by observations. In this article, we will study twists released by 30 off-limb rotational solar coronal jets, and compare the observational findings with theoretical kink instability thresholds. We have found that: (1) the number of events with more twist release becomes less; (2) each of the studied jets has released a twist number of at least 1.3 turns (a twist angle of 2.6π); and (3) the size of a jet is highly related to its twist pitch instead of twist number. Our results suggest that the kink instability threshold in the solar atmosphere should not be a constant. The found lower limit of twist number of 1.3 turns should be merely a necessary but not a sufficient condition for a finite solar magnetic flux rope to become kink unstable.
1. Introduction
Eruption of solar magnetic flux ropes (see reviews in e.g., Raouafi, 2009; Schrijver, 2009; Chen, 2011; Filippov et al., 2015a; Karpen, 2015) has been considered as one of the main drivers of the so-called “space weather”. According to magnetohydrodynamics (MHD) theories, highly twisted magnetic flux ropes, with finite length, are subject to the kink instability, which will develop and finally lead to a release of energy when the stored twist exceeds a certain threshold. Various theoretical studies have given similar but different estimations of the kink-unstable threshold. The Kruskal-Shafranov limit (Kruskal and Kulsrud, 1958; Shafranov, 1963) suggests a kink-unstable threshold 2π of the total twist angle in axisymmetric toroidal magnetized plasma columns. Further study on line-tying force-free coronal loops with uniform twist by Hood and Priest (1981) suggested a maximum twist angle of 2.5π a kink-stable, cylindrical flux tube might contain. 3D MHD numerical simulations (e.g., Pariat et al., 2009) gave a slightly higher limit of the twist angle, 2.6π, injected into the system for the onset of kink instability and the eruption of a solar coronal jet. Dungey and Loughhead (1954) suggested a kink-unstable threshold of 2l/R, where l and R are the length and radius of the flux rope, respectively. All the above thresholds for kink instabilities are, however, theoretical and therefore somehow idealized. The realistic threshold(s) for flux ropes to become unstable in the solar atmosphere need to be further studied, and the theoretical thresholds to be confirmed or refuted observationally.
Kink-unstable magnetic flux ropes in the solar atmosphere could account for a wide range of observational phenomena. For instances, Hood and Priest (1979) suggested, using theoretical considerations, that the kink instability could be a main cause of solar flares. This hypothesis has then been supported by a number of observational studies (e.g., Pevtsov et al., 1996; Srivastava et al., 2010; Liu et al., 2016d). It has also been suggested that kink instability could be associated with small-scale (nano) flares (e.g., Browning et al., 2008). Meanwhile, plenty of literature is available to present abundant evidence in theories, numerical simulations and observations on how the eruptions of filaments and CMEs are related to kink-unstable magnetic flux ropes (e.g., Rust and Kumar, 1996; Kliem et al., 2004; Török and Kliem, 2005; Williams et al., 2005; Guo et al., 2010; Kumar et al., 2012; Liu et al., 2016c; Cheng et al., 2017; Vemareddy et al., 2017). Recently, Wang et al. (2016) has further identified the magnetic twist inside post-eruption flux ropes in the heliosphere from analyzing 115 magnetic clouds observed at 1 AU. They found the kink-unstable thresholds vary from case to case.
Besides flares and CMEs, rotational solar coronal jets (see reviews, e.g., Shibata et al., 1996; Raouafi et al., 2016) have also been suggested to be linked to kink instabilities. The relationship between rotational jets and kink instability has been further investigated in the context of magnetized astronomical jets, which are in scales of light years (e.g., Giannios and Spruit, 2006; Barniol Duran et al., 2017). Most theories and observations suggest that the rotational motion of solar coronal jets should be a process involving “untwisting” (e.g., Shibata and Uchida, 1986; Jibben and Canfield, 2004; Moreno-Insertis et al., 2008; Liu et al., 2009; Liu et al., 2016b; Pariat et al., 2009; Shen et al., 2012; Fang et al., 2014; Filippov et al., 2015b; Lee et al., 2015). The physical scenario of “untwisting” jets is usually described as follows: a newly emerging (e.g., observations in Liu et al., 2016b; Zheng et al., 2018) or a pre-existing closed flux system (disturbed by footpoint motions, e.g., observations in Chen et al., 2017) reconnects with the ambient open magnetic field, during which twists contained in the closed flux system could be passed into the open fields and are then released during the rotational motion of the associated jet. Following this idea, a natural question may be raised: will all the twists stored in the pre-reconnection flux rope be released during the coronal jet eruption? In other words, can we infer the twist stored in the pre-reconnection flux rope from the number of turns a jet rotates after it emerges from the magnetic reconnection?
We shall also note that, not all solar coronal jets show clear rotational motion during their lifetime. Different studies have slightly different (but similar in principle) explanations of why some solar coronal jets rotate and others do not (e.g., Shibata et al., 1996; Moore et al., 2010; Pariat et al., 2015). For examples, models (e.g., Shibata and Uchida, 1986; Canfield et al., 1996) summarized in Shibata et al. (1996) suggested that magnetic reconnection with a sheared/twisted flux system involved could result in a jet with obvious rotational motion. Such modeling approach was further confirmed by Moore et al. (2010) with observations of a number of X-ray jets with rotational motion (named as “blowout jets”) and without rotational motion (named as “standard jets”). On the other hand, numerical simulations in e.g., Pariat et al. (2015) have shown that in certain circumstances, a non-rotational jet (names as “straight jet”) may precede a rotational jet (named as “helical jet”) and influences the behavior of the rotational jet. Most importantly, in the simulations, the pre-eruption twisted magnetic flux rope is not directly involved in the eruption of the non-rotational jet.
We have found it difficult to directly compare the twist released by a rotational jet and stored in its pre-reconnection flux rope from observations, because: (1) for an off-limb jet, we do not have accurate vector magnetic field observations to investigate the underlying magnetic flux rope; and (2) for an on-disk jet, even though we can study the twist stored in the underlying magnetic flux rope using magnetic field extrapolations (which are not always accurate), it is hard to investigate the rotational motion of the jet using imaging observations and spectral observations with sufficiently enough spectral resolution at the needed temperatures. Fortunately, we could find a way to go ahead for an answer from realistic numerical simulations. Pariat et al. (2016) performed a series of 3D MHD numerical simulations of solar jets in conditions with different plasma-β, where the plasma-β is the ratio of the plasma pressure to the magnetic pressure. It has been found that, under chromospheric and coronal conditions where plasma-β is less than unity, the number of turns a jet rotates is almost the same with the twist injected into the system before eruption. Though the possible scenario of partial eruption was not included in their simulations, we have found, through a detailed observational and numerical study of solar coronal twin jets (Liu et al., 2016a), that the residual twist remaining after the jet eruption is very small when compared to the total twist stored in the pre-eruption magnetic flux rope. Therefore, we may safely conclude that the twist released by a solar coronal jet should be the lower limit of, and most likely be similar to, the twist stored in the pre-reconnection flux rope.
In this article, we study twists released by 30 rotational solar coronal jets observed off-limb from 2010 to 2016, and compare them to kink instability thresholds proposed by theoreticians. The paper is organized as follows: data collection and event selection are presented in section 2; detailed examples of the analysis of two typical coronal jets are shown in sections 3, 4 is devoted to statistical results; conclusions and discussions are given in section 5.
2. Data Collection
The Solar Dynamics Observatory (SDO, Pesnell et al., 2012), lunched in 2010 to a geosynchronous orbit, carries three different scientific instruments, among which one is the Atmospheric Imaging Assembly (AIA, Lemen et al., 2012). Data used in this research was obtained from one of the AIA broadband images at He II 304 Å targeting at plasmas with a characteristic temperature of 0.05 MK. All the images were taken at a cadence of 12 seconds with a pixel size of 0.6″(Lemen et al., 2012). To find usable coronal jet events for the purpose of this research, we performed the following steps to explore and collect data:
• First of all, we used the built-in Heliophysics Event Knowledgebase (HEK, Hurlburt et al., 2012, http://www.lmsal.com/hek/index.html) module in SunPy (SunPy Community et al., 2015) to find all events labeled as “coronal jet” (with abbreviation of “CJ”) from the year 2010 to 2017. These events were identified automatically or manually by various research groups from different institutes. To exclude all on-disk events, we then removed all entries with the central location of the event <1.02 solar radii from the disk center. 173 events have been found during this initial search.
• In the second step, we downloaded the movies associated with all the 173 events using links provided in the HEK searching results. Movies of events without given links in their HEK entries were then generated locally from automatically downloaded SDO/AIA 304 Å image sequences.
• Next, all movies were carefully examined one-by-one. On-disk events, which were not eliminated by the first step, were further removed. Limb events which were not clear in AIA 304 Å images were also abandoned. We have also discarded events, when it was not sure whether they were coronal jets or filament eruptions. After applying all the above selection procedures, only 44 events were kept.
• All SDO/AIA 304 Å data were then downloaded with their original cadence (12 s). All the events would be studied into details one-by-one, which will be demonstrated with examples in section 3. We note 14 events which were either too close to other bright structures (e.g., complex loop systems, prominences, etc.) or too faint to allow us to obtain firm parameters, will not be included in the final statistics (section 4). The first two columns in Table 1 show the times and locations of all 30 coronal jet events studied in this research.
Figure 1 depicts the locations of the apparent source regions on the solar limb of all 30 coronal jets, with colors denoting dates of eruptions. Animations of the SDO/AIA 304 Å observations of all 30 events are available at https://github.com/PyDL/jet-stat-movie.
Figure 1. Location of the apparent source regions on the solar limb of all the studied 30 rotational coronal jets from 2010 to 2017. Colors denote dates of eruptions.
3. Examples of Events
In this section, we will show, using two typical examples, our analysis of the temporal evolution, axial and rotational motion of all jets in details. The main difference between these two examples is the different behaviors in their rotational motions: the rotational motion of the first example jet manifests recurrent quasi-parallel stripes in the running-difference time-distance diagram of the slit perpendicular to its axis, while the second example jet shows sinusoidal-like features.
3.1. Coronal Jet on 27 June 2010
A coronal jet, together with a flaring event at its source region, started to erupt from the north-west limb with its root latitude of ~31°, at around 02:16 UT on the 27th June 2010. After rising up to more than 100 Mm above the solar surface, the jet began to fall back from around 02:38 UT and finally arrived at the solar surface at around 02:57 UT. The visualization of whole evolution of this event is available as the online Movie M1, which, again, was generated from a sequence of base-difference AIA 304 Å images. Apparent “whip-like” motion, which has also been observed in many other jets (e.g., Shibata et al., 1996), could be observed during the very early stage of its eruption. After that, clear signatures of the rotational motion became visible during the ascending phase of this jet. Its rotational motion stopped before the jet reached its maximum projected length (the projection of its real length in the plane of the sky). Figure 2A depicts a snapshot of this jet at 02:38 UT observed by SDO/AIA 304 Å. The green dashed line is a 50-pixel (~22 Mm) wide slit along the jet axis, and the blue dotted line is a 30-pixel (~13 Mm) wide slit perpendicular to the jet axis.
Figure 2. (A) is an off-limb coronal jet observed by SDO/AIA 304 Å passband on 27th June 2010 at 02:38 UT. Green dashed and blue dotted lines are slits along and perpendicular to the jet axis, respectively. (B) is a running difference time-distance plot along the jet axis. Three yellow vertical dashed lines represent the starting, peak and ending time of the jet, respectively. The red horizontal dashed line denotes the maximum projected length of the jet. Green dashed curves indicate trajectories of sample sub-jets used to estimate the average projected axial speed of the jet during the period of its rotational motion. (C) is a running difference time-distance plot perpendicular to the jet axis. Two blue vertical dashed lines denote the estimated starting and ending time of the rotational motion, respectively. Purple solid lines are linear fittings of the inclined stripes. The red arrow marks a stripe as part of the evidence of the rotational motion. See the text for reasons why these stripes are manifestations of rotational motion, instead of “whip-like” motion.
Figure 2B shows the time-distance diagram of the green slit in panel (a) based on running-difference images of SDO/AIA 304 Å observations. We can identify a number of fine structures, as parts of the whole jet, appearing as alternating black and white curves in the time-distance diagram. These fine structures, known as “sub-jets”, erupt successively and are common in many solar coronal jet events (e.g., Liu et al., 2014). Parabolic fittings to sample sub-jets (indicated by green dashed curves in Figure 2B) reveal an average axial speed of 148.0 ± 20.6 km s−1 of the jet during the period of its rotational motion (see next paragraph). The error of the average axial speed is the standard deviation of the linear speeds of the sample sub-jets (green dashed curves). The three vertical dashed lines (yellow) represent the starting, peak and ending time of the jet, respectively. Meanwhile, the horizontal dashed line (red) represents the maximum projected length (~133 Mm) of the jet, which was directly determined from the time-distance diagram.
It is worth noticing that, from Figure 2B, these successive sub-jets reached their maximum projected distances at different times. Via a careful examination on these sub-jets, we have found a trend that: (1) when the maximum projected distance is above about 100 Mm, a sub-jet with a larger maximum projected distance reached its maximum projected distance slightly earlier; and (2) when the maximum projected distance is below about 100 Mm, a sub-jet with a smaller maximum projected distance reached its maximum projected distance earlier. A similar behavior of sub-jets in a rotational coronal jet was also be found in our earlier study (Liu et al., 2014). We conjecture that the above behavior was caused by the following reason: (1) at the beginning of the eruption, a sub-jet erupted earlier only had a slightly higher initial speed than its successive sub-jet and thus reached its maximum projected distance somewhat earlier; and (2) toward the end of the magnetic reconnection, a sub-jet erupted later had a significantly lower initial speed than its previous sub-jet, and thus reached its maximum projected distance earlier. A similar relationship between the initial speeds of sub-jets can be found from column 1 of Table 1 in Liu et al. (2014).
Figure 2C shows the time-distance diagram of the blue slit in Figure 2A, again, based on running-difference images of SDO/AIA 304 Å observations. Lower values in distance correspond to lower latitudes in the blue slit. Several inclined quasi-parallel stripes could be identified, indicating plasma material moving toward higher latitudes. One might wonder whether these movements are the manifestation of either a “whip-like” motion or a rotational motion. We consider this as a representation of a rotational motion of the jet, mainly based on the following: (1) visual check on the temporal evolution of the jet in SDO/AIA 304 Å images suggests a significant rotational motion of the jet (online Movie M1); (2) a “whip-like” motion would lead to the jet moving as a whole. However, what we see here from Figure 2C, is that the location of the jet body stays almost unchanged, while plasma moves from one side to another across the width of the jet; and (3) there are signs of stripes “changing direction” in the time-distance diagram (indicated by the red arrow) showing plasmas travel from the front (back) to the back (front). Time-distance diagrams along slits perpendicular to most jets in this study reveal quasi-parallel trajectories as what we show here in Figure 2C. We suggest, this could have been caused by the combined effect of that: (1) most jets are intensive, and (2) the SDO/AIA 304 Å passband is optically thick, meaning that it would be hard to see movements of material behind the “back” of jets.
The purple solid lines in Figure 2C denote some typical inclined quasi-parallel stripes. The average rotational speed of the jet is then estimated as the average slope of linear fittings of these stripes. We note that, these example quasi-parallel stripes do not indicate the different rotating phases of the same material. In other words, stripes perpendicular to these marked example stripes, which mark the rotational motion of material behind the “back” of the jet, are not seen clearly between them. Again, we suggest, this was caused by the combined effect of that: (1) most jets are intensive, and (2) the SDO/AIA 304 Å passband is optically thick. Further, the average rotational period of the jet is defined as twice of the average duration of these stripes. The average rotational speed and period of this jet are therefore estimated as 176.6±49.8 km s−1 and 5.8±2.3 min, respectively. Taking into account the starting time (02:22:14 UT, blue vertical dashed line on the left) and ending time (02:35:50 UT, blue vertical dashed line on the right), we can now estimate that the jet has rotated 2.3±0.9 turns. We shall note, the ending time of the rotational motion was determined by investigating the time-distance plot (Figure 2C) and the original observations of the event (online movie M1). Clear evidence of rotational motion has vanished in both of the time-distance plot and original observations after the determined ending time. We demonstrate that the error of the determined ending time is <5 min, which introduces an uncertainty of <0.9 of the total number of turns the jet has rotated.
3.2. Another Coronal Jet on 27 June 2010
Unlike the previous example in section 3.1, this second jet erupted without strong flaring signatures at its source region. The jet began to erupt from the north-east limb of the Sun with a root latitude of ~19° at around 05:50 UT on the same day as the first example. There are some other differences between this jet and the first example, namely, this jet: (1) did not show apparent “whip-like” motion during its early stage of eruption; and (2) reached a height which is beyond the FOV of SDO/AIA. However, this jet also showed some signatures of rotational motion during its ascending phase and stopped rotating before it reached its maximum projected length. The jet fell back to the solar surface at around 06:56 UT, having a lifetime of more than an hour. The visualization of whole evolution of this event is available as the online Movie M2, which was generated from a sequence of base-difference AIA 304 Å images. Figure 3A depicts a snapshot of this jet at 06:18 UT observed at the SDO/AIA 304 Å passband. Again, the green dashed line is a 50-pixel (~22 Mm) wide slit along the jet axis, and the blue dotted line is a 30-pixel (~13 Mm) wide slit perpendicular to the jet axis.
Figure 3. (A) similar to Figure 2A, but for another off-limb coronal jet event on 27th June 2010 at 06:18 UT. (B) is a running difference time-distance plot along the jet axis. Again, three yellow vertical dashed lines represent the starting, peak and ending time of the jet, respectively. The red dashed curve denotes the parabolic fitting of the outermost sub-jet. Green dashed curves (together with the red dashed curve) indicate trajectories of sample sub-jets used to estimate the average projected axial speed of the jet during the period of its rotational motion. (C) is a running difference time-distance plot perpendicular to the jet axis. Two blue vertical dashed lines denote the estimated starting and ending time of the rotational motion, respectively. Purple solid curves are sinusoidal fittings of corresponding sinusoidal-like trajectories. These sinusoidal features are evidence of the rotational motion instead of the kink motion (see main text).
From the running-difference time-distance plot in Figure 3B, taken along the green dashed line in Figure 3A during the eruption, we can again find that, there are many “sub-jets,” which could be evidence of successive magnetic reconnections (Liu et al., 2014). The relationship between the times when these sub-jets reached their maximum projected distances is similar to the previous jet. Parabolic fittings to sample sub-jets (indicated by red and green dashed curves in Figure 3B) reveal an average axial speed of 230.6± 29.2 km s−1 of the jet during the period of its rotational motion. Because the jet finally reached a height beyond the FOV of SDO/AIA, we used a parabolic fitting to the trajectory of the outermost sub-jet (red dashed curve), to estimate its maximum projected length. This jet has then been found to reach a maximum projected length of about 320 Mm at around 06:18 UT (middle yellow dashed line in Figure 3B).
Similar to Figures 2C, 3C shows the time-distance diagram taken at the location of the blue slit in Figure 3A, based on the running-difference images of SDO/AIA 304 Å observations. Lower values in distance correspond to higher latitudes in this example. Instead of inclined quasi-parallel stripes, we can find several sinusoidal-like features (e.g., indicated by purple solid curves) in the time-distance diagram in Figure 3C. These sinusoidal-like features clearly manifest the rotational motion, instead of the kink motion, of the jet. In the case of the jet undergoing a kink motion, the whole jet body would move forwards and backwards periodically. This, however, is not found in Figure 3C. To evaluate the rotational period, we have performed sinusoidal fittings to these trajectories:
where, A, ω, x0 and y0 are the amplitude, frequency, phase shift and vertical shift, respectively. We shall note, due to the complexity of the sinusoidal fitting, one needs usually to make a good initial guess of the above parameters as inputs of the fitting, to avoid being trapped in a local minimum of the χ2. For a given trajectory, we have used its average y-value as the initial guess of the vertical shift y0, and Fast Fourier Transform (FFT) on the series of the y-value of the trajectory to determine the initial guesses of the other three parameters. The average rotational period of this jet is 9.8 ± 5.7 min. Considering that the rotational motion started at around 05:51:50 UT (blue vertical dashed line on the left in Figure 3C) and ended at around 06:06:14 UT (blue vertical dashed line on the right in Figure 3C), the total number of turns the jet rotated is then estimated as 1.5 ± 0.8 turns. Similarly to the previous example, the ending time of the rotational motion was determined by investigating the time-distance plot (Figure 3C) and the original observations of the event (online Movie M2). The error of the determined ending time is, again, <5 min, resulting in <0.5 of the total number of turns the jet has rotated. The average rotational speed of this jet is estimated to be 108.1 ± 69.7 km s−1.
4. Statistical Results
Table 1 lists parameters that we have obtained from all 30 rotational coronal jets studied. The first column represents the starting time of each jet, the second column (θ) is the latitude with positive (negative) values for the northern (southern) hemisphere, the third column (LT) is the lifetime and the forth column (Lmax) is the maximum projected length obtained from time-distance plots along their axes similar to Figures 2B, 3B. Jets of which we could not find clear evidence of when they fell back to the solar surface are denoted with “-” with their lifetimes. The fifth column is the average width () of jet, obtained by averaging the distances between two edges of the jet when it reached its maximum projected length. The five remaining columns are the axial speed (va), rotational speed (vr), rotational period (Pr), duration of the rotational motion (Dr), and total number of turns of the rotation (Tr), respectively.
Figure 4 depicts distributions of the lifetime, width, rotational speed and rotational period of the jets studied. Pink bars are the distributions of corresponding parameters, while blue dashed curves are the Gaussian fittings to these distributions. μ and σ are the arithmetic mean and standard deviations, respectively. The lifetime (Figure 4A) of all jets (excluding those that we are not sure when or whether they fell back to the solar surface) ranges from 20 to 80 min, with 73% being within the 1-σ range (25–58 min). The arithmetic mean of the lifetime is found to be about 42 min. Similarly, the average width (Figure 4B) of all jets ranges from 10 to 50 Mm, with an arithmetic mean of about 21 Mm. The distribution of the width matches well with a Gaussian distribution. The widths of about 67% jets lie within the 1-σ range (15–24 Mm). Note that for a perfect Gaussian distribution, the total probability within the 1-σ range is about 68%.
Figure 4. Statistics of the lifetime (A), average width (B), rotational speed (C) and rotational period (D) of the studied rotational coronal jets. Blue dashed curves are Gaussian fitting results. μ and σ are the arithmetic mean and standard deviation, respectively.
Figures 4C,D are the distributions of the rotational speed and rotational period, respectively. The rotational speed of jets ranges from 30 to 210 km s−1, with an arithmetic mean of about 88 km s−1 and about 73% of them lying within the 1-σ range (49–127 km s−1). Out of 30 jets, 29 have a rotational period ranging from 3 to 15 min. The arithmetic mean is about 7 min, with about 77% lying within the 1-σ range (4–10 min). All the studied jets have an axial speed ranging from 80 to 330 km s−1 (Figure 5A). The arithmetic mean is about 145 km s−1, with about 67% of the jets lying within the 1-σ range (85–205 km s−1).
Figure 5. Statistics of the average axial speed during the period of the rotational motion (A), the projected length (A), the duration of rotational motion (B) and the total number of turns (C) of all studied solar coronal jets. The blue dashed curve in (A) is the Gaussian fitting result. μ and σ are the arithmetic mean and standard deviation, respectively. Blue dashed curves in (B–D) are exponential fitting results. r is the exponential index. The integrated probabilities on the left of the black vertical dashed lines are 0.8.
Different from the distributions of lifetime, width, rotational speed and rotational period, the frequency of events decreases with increased projected length (Figure 5B). The distribution of the maximum projected length could be fitted with an exponential function f∝eγh, where the index γ is found to be about -0.01. The projected length of jets ranges from 70 to 320 Mm, with 80% <220 Mm (black dashed line in Figure 5B). The frequency of events also decreases with increased duration of the rotational motion (Figure 5C). Eighty percent of the studied events rotated for <25 min.
Most importantly, the frequency of events also decreases with increased total number of turns of the rotational motion, and could be fitted well with an exponential function with an index γ of –0.85 (Figure 5D). This indicates that, the number of events with more twist release becomes less. All jets are associated with twist angles (Φj = 2πTr) between 2.6 and 9.4π. Among all the 30 studied events, there is not a single event revealing a rotational motion with <1.3 turns (twist angle of 2.6π). Moreover, 80% of the events released twist angles <5.6π (twist numbers <2.8 turns).
Figure 6 shows the dependencies of the total number of turns of the rotational motion on the lifetime (Figure 6A), the maximum projected length (Figure 6B), the average width (Figure 6C), the rotational speed (Figure 6D), the rotational period (Figure 6E), and the duration of the rotational motion (Figure 6F) of the jets investigated. From both the scatters and the correlation coefficients (CCs), we can find that there are neither positive nor negative correlations between the total number of turns of the rotational motion and the lifetime (the projected length, the average width or the rotational speed). Even though, the total number of turns was derived from the rotational period and duration of the rotational motion, the total number of turns has very low (moderate) correlation with the rotational period (duration of the rotational motion).
Figure 6. Correlation between various jet parameters (A) lifetime, (B) maximum projected length, (C) width, (D) rotational speed, (E) rotational period, and (F) duration of the rotational motion and the total number of turns these jets rotate. The black vertical lines are indicating the associated errors of the total number of turns. The orange horizontal lines show the errors of the corresponding parameters.
However, the product of the duration and the speed of the rotational motion has a strong positive correlation (with a CC of 0.72) with the total number of turns (Figure 7A). Considering that the twist a jet may release would be the lower limit of, and very likely similar to, the twist its pre-eruption flux rope contains, this strong correlation may indicate that a flux rope with a higher twist number would result in a jet with either a longer or a faster rotational motion. The linear fitting (blue dashed line in Figure 6A) using the “fitxy.pro” in the SSW package, which also accounts for the error bars, suggests that there might be a lower cut-off value around 1.0 of the total number of turns.
Figure 7. (A) shows the correlation between the product of the duration and speed of the rotational motion and the total number of turns. (B) demonstrates correlation between the twist pitch and volume of jets (see Equation 3). Black vertical lines are errors of the total number of turns. Orange horizontal lines are propagated errors of corresponding parameters. Blue lines are the linear fitting results of the scattered points with error bars taken into account.
The twist angle and twist number of a magnetic flux rope are defined as:
Here, l, r, Bϕ, and Bz are the length, width, azimuthal and axial magnetic field strength of the pre-eruption magnetic flux rope, respectively. Then, the twist pitch lp = 2πl/Φ, represents the length traveled along the axis when the magnetic field rotates for a full turn. Figure 7B shows the relationship between the twist pitch of jets and the volume of jets. The twist pitch (lpj) and volume (V) of a jet are defined as:
respectively. All variables in the above equation have the same meanings as defined in the caption of Table 1. It is shown in Figure 7B, that the volume of a jet is correlated positively very well to its twist pitch lpj, with a CC of 0.77. Besides, the CCs between the twist pitch of jets and their length Lmax and average width are 0.39 and 0.52, respectively. If we consider that the twist pitch is conservative during magnetic reconnections (e.g., Birn and Priest, 2007), the above results may suggest that the size of a jet should not be determined by the total twist stored in its pre-reconnection magnetic flux rope, instead, by the twist pitch of its pre-reconnection magnetic flux rope.
5. Conclusions and Discussions
In this research, using high spatial- and temporal-resolution observations obtained at the SDO/AIA 304 Å passband, we have studied the detailed temporal and spatial evolution, especially the rotational motion, of 30 off-limb rotational solar coronal jets that had erupted between 2010 and 2017. These jets were obtained from the HEK database, and were identified either automatically or manually by different groups from different colleagues, to minimize any possible selection bias.
One of the major findings of this study is that all the rotating jets have rotated at least 1.3 turns during their lifetime. The number (occurrence/frequency) of jets decreases almost exponentially with increased total number of turns they have rotated. Most (80%) of them have rotated <2.8 turns. Note again that when plasma-β is <1 the twist released by a rotational jet is the lower limit of the twist stored in the pre-eruption magnetic flux rope (Pariat et al., 2016). From the above results, we conclude that flux ropes that finally erupted as coronal jet events contain twists of at least 1.3 turns (Φ = 2.6π). This value is highly consistent with the suggested kink instability threshold given by various theories and numerical simulations (e.g., Hood and Priest, 1981; Pariat et al., 2009). However, twists released by the studied jets are different from each other, indicating that the kink instability threshold in the solar atmosphere should not be seen as a constant. Further, the exponential decrease of the number of events with increased total number of turns these jets have rotated, suggests that, the more twist the pre-eruption magnetic flux rope contains, the rarer the event is. Most magnetic flux ropes associated with coronal jets would become unstable before their stored twist number (twist angle) is accumulated to 2.8 turns (5.6π). All the results we report here, suggest that containing a twist number of 1.3 turns should be a necessary but not a sufficient condition for a finite solar magnetic flux rope to become kink unstable.
We have found no clear correlation between the released twist by jets and their measured characteristic parameters including lifetime, maximum projected length, width, rotational speed and rotational period. However, there is a strong positive correlation (with a CC of 0.72) between the released twist and the product of the duration and the speed of the rotational motion, suggesting that pre-eruption magnetic flux ropes with higher twists tend to generate jets that rotate either for longer durations or with faster rotational speeds. On the other hand, we have found very strong positive correlation (with a CC of 0.77) between the jet twist pitch and volume of jets, indicating that a pre-eruption magnetic flux rope with a higher twist pitch would most likely result in a larger jet.
All the jets studied in this research have projected lengths of at least 70 Mm and lifetimes of at least 20 min. Therefore, one might wonder whether the lower limit of the twist (1.3 turns) released by rotating solar coronal jets would be different if we study more coronal jets with shorter length and lifetime? We shall note that: (1) Liu et al. (2018) analyzed four homologous recurrent jets, which had minimum length shorter than 20 Mm and lifetime <10 min. Detailed spectral analyses using high resolution IRIS data of one of these jets have suggested that, it released a twist of 1.3 turns (twist angle of 2.6π). (2) We have found no clear correlation between the twist released by jets and their length or lifetime. And, (3) the linear fitting in Figure 7 suggests a cut-off value of around 1.0 turns of the twist released by jets. Although more observations would be needed to confirm this conjecture, we suggest that the lower limit of the twist released by rotating solar coronal jets would not be significantly different.
The maximum twist of the jets analyzed here have released 4.7 turns (corresponding to a twist angle of 9.4π). The jet which triggered a coronal mass ejection event studied in Liu et al. (2015) was found to release a twist of at least 3.3 turns (twist angle of 6.6π). Liu et al. (2014) investigated in details the rotational motion and kinetic energy sources of a coronal jet erupted in July 2012. That jet studied in Liu et al. (2014) was still rotating even at the end of its descending phase and, finally, released a twist of at least 5.1 turns (corresponding to a twist angle of 10.2π). Throughout the study of a pair of solar coronal twin jets and their preceding jet, Liu et al. (2016a) found the preceding jet rotated for at least 8.9 turns (twist angle of 17.8π) during its lifetime. Considering all these previous studies and our findings in this research, we cannot conclude about an upper limit for the twist released by rotating solar coronal jets or stored by pre-eruption magnetic flux ropes. Instead, we have found that magnetic flux ropes with very high twist numbers (> 2.8 turns or 5.6π), are much less in number.
Dungey and Loughhead (1954) suggested a kink-unstable threshold (in units of radians) of ωl/R, where l and R are the length and radius of the flux rope, respectively. They found the constant ω as 2. If we assume that: (1) the magnetic twist pitch is conserved during the magnetic reconnections which triggered the observed jets (e.g., Birn and Priest, 2007); (2) the twists released by the observed coronal jets are similar to those stored in the pre-eruption magnetic flux ropes (e.g., Liu et al., 2016a; Pariat et al., 2016); and (3) the average width of each observed jet is similar to the diameter of its associated pre-eruption magnetic flux rope, we then have:
Here, all the variables have the same meaning as defined before. After removing a significant outlier (event 18), the correlation between Φj and l/R (defined in Equation 4) is found to be as high as 0.73. Linear fitting between Φj and l/R reveals an average value of about 0.6 of ω, with the maximum value of ω <2. Wang et al. (2016) studied the twist angle and other properties of 115 magnetic clouds observed at 1 AU. They also found an average ω value of 0.6 and a maximum ω value of <2. These highly consistent results between solar coronal jets and inter-planetary magnetic clouds, suggest the high possibility of a universal mechanism of different types of eruptions of magnetic flux ropes in the solar upper atmosphere. Nonetheless, we note that, we have made a number of assumptions to obtain the above results. Future work will be needed to focus on examining the above results without making such assumptions.
We would also stress that all the above findings are related to “large-scale” rotational solar coronal jets. Solar coronal jets with no rotational motion are also common, but are suggested to be resulted from different mechanisms and not directly being related to magnetic flux ropes (e.g., Moore et al., 2010; Pariat et al., 2015; Sterling et al., 2015). On the other hand, we would not extend our findings to the ubiquitous type-II spicules (“small-scale” jets with typical length <10 Mm and lifetime around 10 s, de Pontieu et al., 2007), which are believed to be resulted from magnetic reconnections in the upper photosphere or lower chromosphere, with the plasma-β being higher than unity (e.g., Shibata et al., 2007).
Author Contributions
JL conducted the analysis of the data and drafted the manuscript. YW and RE led the interpretation of the result and commented on the manuscript. All authors reviewed the manuscript.
Funding
RE and JL are grateful to Science and Technology Facilities Council (STFC grant nr ST/M000826/1) UK and to The Royal Society for the support received. YW acknowledges the support received from the National Natural Science Foundation of China (NSFC grant no 41774178 and 41761134088).
Conflict of Interest Statement
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.
Acknowledgments
The Solar Dynamics Observatory is the first mission for NASA's Living With a Star (LWS) Program. Sunpy is an open-source and free community-developed solar data analysis package written in Python.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fspas.2019.00044/full#supplementary-material
Video 1. Online movie M1, generated from a sequence of base-difference AIA 304 A images for the example jet in Section 3.1.
Video 2. Online movie M2, generated from a sequence of base-difference AIA 304 A images for the example jet in Section 3.2.
References
Barniol Duran, R., Tchekhovskoy, A., and Giannios, D. (2017). Simulations of AGN jets: magnetic kink instability versus conical shocks. Month. Notices R. Astron. Soc. 469, 4957–4978. doi: 10.1093/mnras/stx1165
Birn, J., and Priest, E. R. (2007). Reconnection of Magnetic Fields : Magnetohydrodynamics and Collisionless Theory and Observations. New York, NY: Cambridge University Press. doi: 10.1017/CBO9780511536151
Browning, P. K., Gerrard, C., Hood, A. W., Kevis, R., and van der Linden, R. A. M. (2008). Heating the corona by nanoflares: simulations of energy release triggered by a kink instability. Astron. Astrophys. 485, 837–848. doi: 10.1051/0004-6361:20079192
Canfield, R. C., Reardon, K. P., Leka, K. D., Shibata, K., Yokoyama, T., and Shimojo, M. (1996). H alpha Surges and X-Ray Jets in AR 7260. Astrophys. J. 464:1016. doi: 10.1086/177389
Chen, J., Su, J., Deng, Y., and Priest, E. R. (2017). A complex solar coronal jet with two phases. Astrophys. J. 840:54. doi: 10.3847/1538-4357/aa6c59
Chen, P. F. (2011). Coronal mass ejections: models and their observational basis. Living Rev. Sol. Phys. 8:1. doi: 10.12942/lrsp-2011-1
Cheng, X., Guo, Y., and Ding, M. (2017). Origin and structures of solar eruptions I: magnetic flux rope. Sci. China Earth Sci. 60, 1383–1407. doi: 10.1007/s11430-017-9074-6
de Pontieu, B., McIntosh, S., Hansteen, V. H., Carlsson, M., Schrijver, C. J., Tarbell, T. D., et al. (2007). A tale of two spicules: the impact of spicules on the magnetic chromosphere. Publ. Astron. Soc. Jpn. 59, S655–S662. doi: 10.1093/pasj/59.sp3.S655
Dungey, J. W., and Loughhead, R. E. (1954). Twisted magnetic fields in conducting fluids. Aust. J. Phys. 7:5. doi: 10.1071/PH540005
Fang, F., Fan, Y., and McIntosh, S. W. (2014). Rotating solar jets in simulations of flux emergence with thermal conduction. Astrophys. J. 789:L19. doi: 10.1088/2041-8205/789/1/L19
Filippov, B., Martsenyuk, O., Srivastava, A. K., and Uddin, W. (2015a). Solar magnetic flux ropes. J. Astrophys. Astron. 36, 157–184. doi: 10.1007/s12036-015-9321-5
Filippov, B., Srivastava, A. K., Dwivedi, B. N., Masson, S., Aulanier, G., Joshi, N. C., et al. (2015b). Formation of a rotating jet during the filament eruption on 2013 April 10-11. Month. Notices R. Astron. Soc. 451, 1117–1129. doi: 10.1093/mnras/stv1039
Giannios, D., and Spruit, H. C. (2006). The role of kink instability in Poynting-flux dominated jets. Astron. Astrophys. 450, 887–898. doi: 10.1051/0004-6361:20054107
Guo, Y., Ding, M. D., Schmieder, B., Li, H., Török, T., and Wiegelmann, T. (2010). Driving mechanism and onset condition of a confined eruption. Astrophys. J. Lett. 725, L38–L42. doi: 10.1088/2041-8205/725/1/L38
Hood, A. W., and Priest, E. R. (1979). Kink instability of solar coronal loops as the cause of solar flares. Solar Phys. 64, 303–321. doi: 10.1007/BF00151441
Hood, A. W., and Priest, E. R. (1981). Critical conditions for magnetic instabilities in force-free coronal loops. Geophys. Astrophys. Fluid Dyn. 17, 297–318. doi: 10.1080/03091928108243687
Hurlburt, N., Cheung, M., Schrijver, C., Chang, L., Freeland, S., Green, S., et al. (2012). Heliophysics event knowledgebase for the solar dynamics observatory (SDO) and beyond. Solar Phys. 275, 67–78. doi: 10.1007/s11207-010-9624-2
Jibben, P., and Canfield, R. C. (2004). Twist Propagation in Hα Surges. Astrophys. J. 610, 1129–1135. doi: 10.1086/421727
Kliem, B., Titov, V. S., and Török, T. (2004). Formation of current sheets and sigmoidal structure by the kink instability of a magnetic loop. Astron. Astrophys. 413, L23–L26. doi: 10.1051/0004-6361:20031690
Kruskal, M. D., and Kulsrud, R. M. (1958). Equilibrium of a magnetically confined plasma in a toroid. Phys. Fluids 1, 265–274. doi: 10.1063/1.1705884
Kumar, P., Cho, K.-S., Bong, S.-C., Park, S.-H., and Kim, Y. H. (2012). Initiation of coronal mass ejection and associated flare caused by helical kink instability observed by SDO/AIA. Astrophys. J. 746:67. doi: 10.1088/0004-637X/746/1/67
Lee, E. J., Archontis, V., and Hood, A. W. (2015). Helical blowout jets in the sun: untwisting and propagation of waves. Astrophys. J. 798:L10. doi: 10.1088/2041-8205/798/1/L10
Lemen, J. R., Title, A. M., Akin, D. J., Boerner, P. F., Chou, C., Drake, J. F., et al. (2012). The atmospheric imaging assembly (AIA) on the solar dynamics observatory (SDO). Solar Phys. 275, 17–40. doi: 10.1007/s11207-011-9776-8
Liu, J., Erdélyi, R., Wang, Y., and Liu, R. (2018). Untwisting jets related to magnetic flux cancellation. Astrophys. J. 852:10. doi: 10.3847/1538-4357/aa992d
Liu, J., Fang, F., Wang, Y., McIntosh, S. W., Fan, Y., and Zhang, Q. (2016a). On the observation and simulation of solar coronal twin jets. Astrophys. J. 817:126. doi: 10.3847/0004-637X/817/2/126
Liu, J., Wang, Y., Erdélyi, R., Liu, R., McIntosh, S. W., Gou, T., et al. (2016b). On the magnetic and energy characteristics of recurrent homologous jets from an emerging flux. Astrophys. J. 833:150. doi: 10.3847/1538-4357/833/2/150
Liu, J., Wang, Y., Liu, R., Zhang, Q., Liu, K., Shen, C., et al. (2014). When and how does a prominence-like jet gain kinetic energy? Astrophys. J. 782:94. doi: 10.1088/0004-637X/782/2/94
Liu, J., Wang, Y., Shen, C., Liu, K., Pan, Z., and Wang, S. (2015). A solar coronal jet event triggers a coronal mass ejection. Astrophys. J. 813:115. doi: 10.1088/0004-637X/813/2/115
Liu, L., Wang, Y., Wang, J., Shen, C., Ye, P., Liu, R., et al. (2016c). Why is a flare-rich active region CME-poor? Astrophys. J. 826:119. doi: 10.3847/0004-637X/826/2/119
Liu, R., Kliem, B., Titov, V. S., Chen, J., Wang, Y., Wang, H., et al. (2016d). Structure, stability, and evolution of magnetic flux ropes from the perspective of magnetic twist. Astrophys. J. 818:148. doi: 10.3847/0004-637X/818/2/148
Liu, W., Berger, T. E., Title, A. M., and Tarbell, T. D. (2009). AN INTRIGUING CHROMOSPHERIC JET OBSERVED BY HINODE : FINE STRUCTURE KINEMATICS AND EVIDENCE OF UNWINDING TWISTS. Astrophys. J. Lett. 707, 37–41. doi: 10.1088/0004-637X/707/1/L37
Moore, R. L., Cirtain, J. W., Sterling, A. C., and Falconer, D. A. (2010). Dichotomy of solar coronal jets: standard jets and blowout jets. Astrophys. J. 720, 757–770. doi: 10.1088/0004-637X/720/1/757
Moreno-Insertis, F., Galsgaard, K., and Ugarte-Urra, I. (2008). Jets in coronal holes: hinode observations and three-dimensional computer modeling. Astrophys. J. 673, L211–L214. doi: 10.1086/527560
Pariat, E., Antiochos, S. K., and DeVore, C. R. (2009). A MODEL FOR SOLAR POLAR JETS. Astrophys. J. 691, 61–74. doi: 10.1088/0004-637X/691/1/61
Pariat, E., Dalmasse, K., DeVore, C. R., Antiochos, S. K., and Karpen, J. T. (2015). Model for straight and helical solar jets. Astron. Astrophys. 573:A130. doi: 10.1051/0004-6361/201424209
Pariat, E., Dalmasse, K., DeVore, C. R., Antiochos, S. K., and Karpen, J. T. (2016). A model for straight and helical solar jets. II. Parametric study of the plasma beta. Astron. Astrophys. 596:A36. doi: 10.1051/0004-6361/201629109
Pesnell, W. D., Thompson, B. J., and Chamberlin, P. C. (2012). The solar dynamics observatory (SDO). Solar Phys. 275, 3–15. doi: 10.1007/s11207-011-9841-3
Pevtsov, A. A., Canfield, R. C., and Zirin, H. (1996). Reconnection and helicity in a solar flare. Astrophys. J. 473:533. doi: 10.1086/178164
Raouafi, N. E. (2009). Observational evidence for coronal twisted flux rope. Astrophys. J. 691, L128–L132. doi: 10.1088/0004-637X/691/2/L128
Raouafi, N. E., Patsourakos, S., Pariat, E., Young, P. R., Sterling, A. C., Savcheva, A., et al. (2016). Solar coronal jets: observations, theory, and modeling. Space Sci. Rev. 201, 1–53. doi: 10.1007/s11214-016-0260-5
Rust, D. M. and Kumar, A. (1996). Evidence for helically kinked magnetic flux ropes in solar eruptions. Astrophys. J. Lett. 464:L199. doi: 10.1086/310118
Schrijver, C. J. (2009). Driving major solar flares and eruptions : a review. Adv. Sp. Res. 43, 739–755. doi: 10.1016/j.asr.2008.11.004
Shafranov, V. D. (1963). Equilibrium of a toroidal plasma in a magnetic field. J. Nuclear Energy 5:251. doi: 10.1088/0368-3281/5/4/307
Shen, Y., Liu, Y., Su, J., and Deng, Y. (2012). ON A CORONAL BLOWOUT JET: THE FIRST OBSERVATION OF A SIMULTANEOUSLY PRODUCED BUBBLE-LIKE CME AND A JET-LIKE CME IN A SOLAR EVENT. Astrophys. J. 745:164. doi: 10.1088/0004-637X/745/2/164
Shibata, K., Nakamura, T., Matsumoto, T., Otsuji, K., Okamoto, T. J., Nishizuka, N., et al. (2007). Chromospheric anemone jets as evidence of ubiquitous reconnection. Science 318:1591. doi: 10.1126/science.1146708
Shibata, K., Shimojo, M., Yokoyama, T., and Ohyama, M. (1996). Theory and Observations of X-ray Jets. ASP Conf. Ser. 111:29.
Shibata, K., and Uchida, Y. (1986). Sweeping-magnetic-twist mechanism for the acceleration of jets in the solar atmosphere. Solar Phys. 103, 299–310. doi: 10.1007/BF00147831
Srivastava, A. K., Zaqarashvili, T. V., Kumar, P., and Khodachenko, M. L. (2010). Observation of kink instability during small B5.0 solar flare on 2007 June 4. Astrophys. J. 715, 292–299. doi: 10.1088/0004-637X/715/1/292
Sterling, A. C., Moore, R. L., Falconer, D. A., and Adams, M. (2015). Small-scale filament eruptions as the driver of X-ray jets in solar coronal holes. Nature 523:437. doi: 10.1038/nature14556
SunPy Community, Mumford, S. J., Christe, S., Pérez-Suárez, D., Ireland, J., Shih, A. Y., et al. (2015). SunPy-Python for solar physics. Comput. Sci. Dis. 8:014009. doi: 10.1088/1749-4699/8/1/014009
Török, T., and Kliem, B. (2005). Confined and Ejective Eruptions of Kink-unstable Flux Ropes. Astrophys. J. Lett. 630, L97–L100. doi: 10.1086/462412
Vemareddy, P., Gopalswamy, N., and Ravindra, B. (2017). Prominence eruption initiated by helical kink instability of an embedded flux rope. Astrophys. J. 850:38. doi: 10.3847/1538-4357/aa9020
Wang, Y., Zhuang, B., Hu, Q., Liu, R., Shen, C., and Chi, Y. (2016). On the twists of interplanetary magnetic flux ropes observed at 1 AU. J. Geophys. Res. Sp. Phys. 121, 9316–9339. doi: 10.1002/2016JA023075
Williams, D. R., Török, T., Démoulin, P., van Driel-Gesztelyi, L., and Kliem, B. (2005). Eruption of a kink-unstable filament in NOAA active region 10696. Astrophys. J. Lett. 628, L163–L166. doi: 10.1086/432910
Keywords: solar eruptions, solar coronal jets, MHD instabilities, magnetic flux ropes, magnetic twists
Citation: Liu J, Wang Y and Erdélyi R (2019) How Many Twists Do Solar Coronal Jets Release? Front. Astron. Space Sci. 6:44. doi: 10.3389/fspas.2019.00044
Received: 16 January 2019; Accepted: 22 May 2019;
Published: 11 June 2019.
Edited by:
Dipankar Banerjee, Indian Institute of Astrophysics, IndiaReviewed by:
P. Vemareddy, Indian Institute of Astrophysics, IndiaAbhishek Kumar Srivastava, Indian Institute of Technology (BHU), India
Copyright © 2019 Liu, Wang and Erdélyi. 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: Jiajia Liu, jj.liu@sheffield.ac.uk