Issue |
A&A
Volume 686, June 2024
|
|
---|---|---|
Article Number | A148 | |
Number of page(s) | 12 | |
Section | The Sun and the Heliosphere | |
DOI | https://doi.org/10.1051/0004-6361/202348734 | |
Published online | 05 June 2024 |
Magnetic helicity evolution during active region emergence and subsequent flare productivity
1
School of Earth and Space Sciences, Peking University, Beijing 100871, PR China
e-mail: liting@nao.cas.cn
2
National Astronomical Observatories, Chinese Academy of Sciences, Beijing 100101, PR China
3
School of Astronomy and Space Science, University of Chinese Academy of Sciences, Beijing 100049, PR China
4
State Key Laboratory of Solar Activity and Space Weather, Beijing 100190, PR China
5
Max Planck Institute for Solar System Research, Göttingen 37077, Germany
Received:
25
November
2023
Accepted:
20
March
2024
Aims. Solar active regions (ARs), which are formed by flux emergence, serve as the primary sources of solar eruptions. However, the specific physical mechanism that governs the emergence process and its relationship with flare productivity remains to be thoroughly understood.
Methods. We examined 136 emerging ARs, focusing on the evolution of their magnetic helicity and magnetic energy during the emergence phase. Based on the relation between helicity accumulation and magnetic flux evolution, we categorized the samples and investigated their flare productivity.
Results. The emerging ARs we studied can be categorized into three types, Type-I, Type-II, and Type-III, and they account for 52.2%, 25%, and 22.8% of the total number in our sample, respectively. Type-I ARs exhibit a synchronous increase in both the magnetic flux and magnetic helicity, while the magnetic helicity in Type-II ARs displays a lag in increasing behind the magnetic flux. Type-III ARs show obvious helicity injections of opposite signs. Significantly, 90% of the flare-productive ARs (flare index ≥ 6) were identified as Type-I ARs, suggesting that this type of AR has a higher potential to become flare productive. In contrast, Type-II and Type-III ARs exhibited a low and moderate likelihood of becoming active, respectively. Our statistical analysis also revealed that Type-I ARs accumulate more magnetic helicity and energy, far beyond what is found in Type-II and Type-III ARs. Moreover, we observed that flare-productive ARs consistently accumulate a significant amount of helicity and energy during their emergence phase.
Conclusions. These findings provide valuable insight into the flux emergence phenomena, offering promising possibilities for early-stage predictions of solar eruptions.
Key words: Sun: flares / sunspots
© The Authors 2024
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.
1. Introduction
Solar active regions (ARs) encompass all the observable occurrences that result from the expansion of the Sun’s magnetic field, stretching from the photosphere to the corona, which become evident through emissions across a broad spectrum of wavelengths (van Driel-Gesztelyi & Green 2015; Toriumi & Wang 2019). Due to their strong magnetic field and complicated magnetic configuration, ARs are the principal source of a broad range of solar activities, ranging from small-scale brightenings and jets (Archontis et al. 2010; Chen et al. 2019; Duan et al. 2022; Sun et al. 2023) to large-scale flares and coronal mass ejections (CMEs) (Li & Zhang 2015; Janvier et al. 2016; Sun et al. 2022). Active regions are formed through the process of flux emergence from the convection zone in the Sun (Cheung et al. 2010; Zhang & Low 2001; Fan 2021). The initial stage of creating a simple bipolar AR involves the emergence of an Ω-shaped flux tube. The interaction between the axial magnetic field of the tube and the photosphere leads to the emergence of two magnetic polarities within the AR. As the flux tube continues to emerge, the two main polarities gradually separate, giving rise to the appearance of small magnetic polarities between them (Cheung & Isobe 2014). Despite numerous simulations having been conducted to study the process of AR emergence (Fan & Fang 2014; Fournier et al. 2017), the specific mechanisms of the emergence process and the correlation with subsequent flare and CME productivity remain subjects of ongoing discussions.
There are various parameters that can be utilized to characterize the properties of ARs (Zhang 2001; Chintzoglou et al. 2019; Li et al. 2022). Welsch et al. (2009) classified these parameters into categories of extensive type and intensive type. Extensive-type parameters generally scale with the size of the AR, such as the magnetic flux Φ and magnetic energy E. On the other hand, intensive-type parameters are normalized quantities, for example, the reconnection flux divided by the total flux Φr/Φt and free energy divided by the total energy Ef/Et. The extensive-type parameters usually correlate with flares, while intensive-type parameters are associated with CMEs (Sun et al. 2015; Kazachenko et al. 2022). Several studies have demonstrated that the distribution of some parameters (e.g., magnetic shear, net current) near the polarity inversion line (PIL) can provide valuable insights into the productivity of solar eruptions (Dalmasse et al. 2015; Kontogiannis et al. 2019). However, during the emergence of ARs, the PIL structure may not be unambiguously identified (Stein & Nordlund 2012; Toriumi et al. 2014). Consequently, to investigate the properties of emerging ARs, it becomes necessary to consider parameters calculated over the whole region of ARs. Theoretical studies simulating the magnetic flux tube emergence have suggested that a higher initial twist of the magnetic tube may result in a higher flux emergence rate d(Φ)/dt (Magara 2001; Murray et al. 2006; Fu & Welsch 2016). Abramenko et al. (2017) measured the flux emergence rate for 36 ARs and found that the relationship between d(Φ)/dt and the peak unsigned flux Φp could be described by a power-law index of 0.69 ± 0.10 for the entire sample of ARs. Kutsenko et al. (2021) performed a statistical analysis and found a significant correlation between flare productivity and flux emergence rate, which offers valuable information for predicting flare productivity during the early phase of ARs.
In our study, we investigate the properties of emerging ARs by employing two key parameters: magnetic helicity and magnetic energy. Magnetic helicity serves as a measure of the net linking and the twisting of magnetic field lines within a given volume (Berger 1999; Pevtsov et al. 2014; Valori et al. 2016). It is important to note that magnetic helicity is only gauge invariant under the condition that the normal component Bn of the magnetic field vanishes on the bounding surface (Woltjer 1958; Berger 1999). However, in the solar atmosphere, this condition is not satisfied due to the presence of a strong vertical magnetic component perpendicular to the solar surface. To address this issue, the concept of relative magnetic helicity is introduced (Berger & Field 1984; Yang et al. 2013, 2018). Relative magnetic helicity is defined as the difference between the total magnetic helicity and the magnetic helicity of a reference field. This quantity maintains gauge invariance even when the boundary condition does not match, making it an effective characterization of helicity in the solar atmosphere. For the sake of simplicity, we refer to relative magnetic helicity as “magnetic helicity” throughout this paper. Regarding magnetic energy, it serves as a fundamental parameter that quantifies the amount of energy stored within magnetic fields (Démoulin & Berger 2003; Trujillo Bueno et al. 2004).
Various studies have demonstrated the relationship between magnetic helicity and solar eruptions (Zhang & Flyer 2008; Jing et al. 2015; Park et al. 2021; Wang et al. 2023). LaBonte et al. (2007) conducted a statistical analysis of 48 X-flaring regions and 345 non-X-flaring regions, revealing that a peak helicity flux threshold of 6 × 1036 Mx2 s−1 is required to produce an X-class flare. Park et al. (2010) performed a statistical study of 378 ARs and found a strong correlation between helicity accumulation and flare productivity. They also observed a distinct difference in the helicity injection rate between flaring and non-flaring ARs. Korsós et al. (2020) and Soós et al. (2022) observed differences in the identified periodicities of magnetic helicity between flaring ARs associated with large M- and X-class flares compared to ARs associated with smaller B- and C-class flares. When characterizing CMEs, the use of intensive-type parameters is generally considered more appropriate (Sun et al. 2015; Toriumi & Takasao 2017). Pariat et al. (2017) proposed the ratio of magnetic helicity of the current-carrying magnetic field to the total relative helicity Hj/Hr as a useful indicator for CME eruptions. Numerous subsequent observational cases have supported this proposition (Moraitis et al. 2019; Thalmann et al. 2019; Gupta et al. 2021). Among them, Gupta et al. (2021) examined ten ARs and found that Hj/Hr effectively distinguishes between CME-productive and CME-poor ARs. However, Duan et al. (2023) conducted a statistical analysis of 45 M- and X-class flares from 30 different ARs and found no systematic differences in the value of Hj/Hr between confined and eruptive flares. The ratio of free magnetic energy to potential energy Ef/Ep and normalized current-carrying helicity Hj/Φ2 instead exhibited the highest discriminative ability between confined and eruptive flares.
Numerous simulations and observations have investigated the transport of helicity from the solar interior to the atmosphere during flux emergence (Pevtsov et al. 2003; Yang et al. 2009; Liu et al. 2014; Prior & MacTaggart 2019; Wang et al. 2021). However, there is a lack of large-sample statistical analysis focusing on the helicity evolution during AR emergence. Moreover, there are only a few statistical investigations examining the relationship between emerging ARs and subsequent eruptions (Kutsenko et al. 2021; Liokati et al. 2022). Our work aims to examine the evolution of helicity in emerging ARs and subsequent AR flare productivity. We analyzed 136 samples that we divided into three categories based on the evolution of magnetic helicity and magnetic flux and found the differences in flare productivity among them. This paper is organized as follows. Sections 2 and 3 describe the dataset and analysis methods used in our study, respectively. In Sect. 4, we divide the samples into three distinct categories according to the evolution of magnetic helicity and magnetic flux. Section 5 presents our statistical results in detail. Finally, we discuss and summarize the major results in Sect. 6.
2. Dataset
The data used in our study were obtained from the Helioseismic and Magnetic Imager (HMI; Scherrer et al. 2012) on board the Solar Dynamics Observatory (SDO), which provides magnetograms of the entire solar disk with a spatial resolution of 0.5″ pixel−1. For our analysis, we mainly used the Spaceweather HMI Active Region Patch (SHARP) data (Bobra et al. 2014), which offers the advantage of automated identification and tracking of AR data. The SHARP data also comprises Lambert cylindrical equal-area (CEA) projections of the magnetic field vector, enabling the conversion of the components from [Bx, By, Bz] to the spherical heliographic components [Bϕ, Bθ, Br] with a cadence of 720 s (Sun 2013), which can be used to calculate the magnetic helicity and magnetic energy.
In total, we selected 136 emerging ARs according to the dataset provided in Kutsenko et al. (2021). In the study conducted by Kutsenko et al. (2021), a total of 243 emerging ARs between 2010 May and 2017 December 31 were examined. The starting and ending times of the emergence phase were automatically determined using a two-segment piecewise continuous linear fitting method applied to the unsigned flux curves (Kutsenko et al. 2019). The selection criteria for these samples are as follows: (i) clear identification of the starting and ending times of the emergence phase and (ii) ensuring that the ARs remain within the central meridian distance of less than 60° throughout the entire process. When we used SHARP to access these samples, in some cases the initial stages of emergence for some ARs were not detected by SHARP. We evaluated the ratio between the data-covering interval provided by SHARP and the overall emergence interval and considered that samples with a ratio above approximately 90% can be regarded as possessing nearly complete temporal coverage. Given that there are seven samples with a ratio between 87% and 90%, we manually set the threshold as 87% in order to include seven more samples and develop a more robust analysis. Consequently, based on the additional criterion we introduced, namely (iii) SHARP data has over 87% time range of the emergence process. In this way, we selected 136 samples in total: 61 flaring ARs and 75 non-flaring ARs (shown in Table A.1).
3. Methods
We computed unsigned magnetic flux, magnetic helicity, magnetic energy and flare index for the selected 136 samples within their respective emergence time intervals.
For unsigned magnetic flux Φ, we performed a summation of the absolute magnetic flux density in pixels on the radial magnetic field component Br provided by SHARP. Only pixels with magnetic flux density values exceeding 18 Mx cm−2 were included in the calculation, corresponding to three times the noise level of the HMI magnetogram (Liu et al. 2012). We used the peak value (normally at the final time step) as the peak unsigned magnetic flux Φp.
We employed a helicity flux method to calculate the magnetic helicity H (Berger 1999). The helicity flux equation, Eq. (1), is shown below:
In the equation, Ap represents the magnetic vector potential, while Bn and Bt denote the normal and tangential magnetic field, respectively. Similarly, Vn and Vt respectively indicate the normal and tangential velocity. Typically, the integral surface S for the helicity flux calculation should be a closed hexahedron situated in the corona, with the photosphere serving as the bottom boundary. However, since the contribution of helicity is primarily provided by the photosphere, we considered the integral surface S to be the photosphere itself (Pevtsov et al. 2014). The first term on the right-hand side of the equation is associated with Vn and is referred to as the “emergence term”. This term reflects the helicity contributed by the emergence of magnetic structures from the convection zone. The second term on the right-hand side is associated with Vt and is known as the “shear term”. This term reflects the helicity contributed by the shearing motion of magnetic structures in the photosphere.
To calculate the velocity vector on the photosphere, we used the differential affine velocity estimator for vector magnetograms (DAVE4VM; Schuck 2008) method. We used an apodization window size of 19 pixels following the suggestion by Schuck (2008). In this method, the achieved velocity undergoes an additional adjustment to account for the exclusion of any irrelevant field-aligned plasma flow (Liu & Schuck 2012; Wang et al. 2022). This correction is achieved using the following equation:
Here, V⊥ represents the velocity perpendicular to the magnetic field lines, while V corresponds to the velocity obtained through the application of the DAVE4VM. The corrected velocity V⊥ serves as the crucial component for calculating the helicity flux. Liu & Schuck (2012) and Liu et al. (2014) also used DAVE4VM to calculate the helicity flux in emerging ARs. They used Monte Carlo methods to estimate the numerical errors when calculating the helicity flux, and the results are considered to be within the margin of errors. Thus, in this paper we regard the calculation as reliable.0]Could you please check the link of footnote 1 or provide another one?
Since we could calculate dH/dt for each time step, we integrated dH/dt over the entire duration of emergence. This integration allowed us to obtain the accumulated change in helicity, denoted as ΔH:
To evaluate the eruptive potential of ARs based solely on their structural complexity rather than their area, it is valuable to examine the helicity normalized by the square of the magnetic flux (Toriumi & Park 2022), denoted as ΔH/.
We employed the energy flux method, Eq. (4),
which is similar to the helicity flux method, in order to calculate magnetic energy E. Once we obtained the energy flux dE/rmdt for each time step, we integrated these values to obtain the accumulated energy, denoted as ΔE:
Flare index (FI) was introduced by Abramenko (2005) to assess the flare productivity of an AR, which provides a quantitative measure of the AR’s ability for flaring. The equation to calculate the FI is presented below:
Here, represents the summation of the GOES peak intensities
for a specific flare class (X, M, C, B), where Nj denotes the number of flares within that class. In the above equation, τ signifies the total duration of the AR observation in days, indicating that this parameter is normalized by time. It should be noted that the parameters Φ, H, and E are calculated during emergence, while FI includes the entire period of an AR until we could no longer observe it.
4. Classification of active regions
Once we had calculated the parameters outlined in Sect. 3, namely, the unsigned magnetic flux, helicity, and energy, we could construct the corresponding curves during the emergence phase. From the 136 samples, we found that almost all the unsigned magnetic flux Φ curves exhibit a monotonically increasing trend. This outcome is expected since the magnetic flux naturally increases as the AR emerges. In contrast, the accumulated magnetic helicity ΔH curves manifested different characteristics. Based on the ΔH and Φ curves, we categorized the samples into three distinct groups. Figure 1 illustrates three representative examples for each category. It is worth mentioning that the intensity-normalized helicity represented by the orange curves refers to the normalization of helicity values into the range of (0, 1) rather than the previously mentioned normalized helicity ΔH/.
![]() |
Fig. 1. Unsigned magnetic flux Φ (blacks curve) and accumulated magnetic helicity ΔH (orange curves) evolutions for nine examples. The top, middle, and bottom rows represent the Type-I, Type-II, and Type-III ARs, respectively. The decreasing helicity indicates a negative helicity injection. The correlation coefficients between the two curves are marked in the upper-left corner of each panel. |
The Type-I ARs (the first row in Fig. 1) exhibited a simultaneous growth in the unsigned magnetic flux Φ and accumulated magnetic helicity ΔH. The absolute value of the Pearson correlation coefficient between them exceeds 0.9. This high correlation indicates that the helicity steadily increases with unsigned flux, and it can be inferred that the helicity injection is mainly dominated by a single sign. For Type-II ARs, the helicity curves exhibit a lag behind the Φ curve during the initial phase of emergence. As depicted in the middle row in Fig. 1, the helicity curves demonstrate a slow initial growth followed by a rapid increase. As a result, the absolute value of the correlation coefficient between them is not more than 0.9. It is important to note that the second example of Type-II ARs displays a consistent helicity decrease. This is because this AR is dominated by a negative helicity. We also calculated the ratio between the lag period and the entire period for each Type-II AR, and the results showed that they exhibit a 50% lag phase on average. In this way, the missing data (up to 13% time period) of some samples would not affect the classification for Type-II ARs. Type-III ARs display a significant injection of opposite helicity, resulting in a clear reversal of the helicity curve. Type-I and Type-II ARs may also exhibit slight opposite helicity injections in addition to their dominant helicity sign. However, they are relatively slight in contrast to the total helicity accumulation. To classify an AR as a Type-III, positive and negative signs of helicity must both be considered to constitute at least 30% of the total unsigned helicity. Based on the criterion mentioned above, there are 71 Type-I ARs, 34 Type-II ARs, and 31 Type-III ARs, accounting for 52.2%, 25%, and 22.8% of the total sample, respectively.
Figure 2 presents an example of Br magnetograms of AR 11620 corresponding to Type-I ARs, which provides a visual depiction of the entire process of AR emergence. During the emergence of the AR, a typical scenario involves the presence of a pair of positive and negative magnetic poles. As the AR emerges, the magnetic flux associated with both poles gradually increases, and the poles start to separate from each other. In the area where the poles meet, the PIL structures develop and are characterized by a significant magnetic field gradient and strong shearing motions. The clockwise rotation of the magnetic footpoints at the poles is also observed. This rotation contributes to the shearing term of the helicity injection within the AR. As the emergence process continues, the magnetic field becomes progressively more complicated. Apart from the two main magnetic poles, numerous smaller magnetic poles are also formed. The interaction of opposite polarities leads to ubiquitous magnetic cancellation. Figure 3 depicts the temporal evolution of the calculated parameters corresponding to Fig. 2. In the top row, panel a, of Fig. 3, we observed that the unsigned flux Φ and accumulated energy ΔE exhibited a monotonic increase over time. Panel b displays the time evolution of helicity flux dH/dt, while panel c presents the accumulated helicity ΔH obtained by integrating dH/dt. Notably, we can observe the strong correlation between the ΔH and Φ curves, a characteristic associated with Type-I ARs. Moreover, we observed that the shear term contributes to the majority of the helicity.
![]() |
Fig. 2. Selected HMI images of the normal component, Bn, of the photospheric field of the Type-I AR 11620 showing the emergence process. The length of the horizontal white line corresponds to 100 arcsec. |
![]() |
Fig. 3. Time profiles of the magnetic flux, helicity, and energy for the Type-I AR 11620. Panel (a): time evolution of the unsigned magnetic flux Φ (black) and accumulated energy ΔE (blue). Panel (b): time profile of the helicity flux, d(H)/dt (black). The red and green curves represent the shear and emergence terms, respectively. Panel (c): time profiles of accumulated helicity ΔH. The red and green curves indicate the helicity of the emergence term and shear term, respectively. |
Figure 4 presents another example of AR 11768, classified as a Type-II AR, and illustrates its emergence from the quiet region. Through these magnetograms, we could observe a similar emergence process with Type-I ARs: a bipole magnetic configuration with a PIL structure in the middle. We could also observe the main pole rotation and the emerging smaller magnetic poles. As a result, it is hard to distinguish Type-I and Type-II ARs solely based on the magnetogram. Figure 5 displays the temporal evolution of the calculated parameters for this AR. We observed a monotonic increase in the Φ, while the ΔH shows a lag behind Φ by about 30 h. Moreover, helicity does not begin to increase rapidly until around 60 h after the start of emergence. We found that the when the helicity is not accumulated prominently, the energy accumulates slowly, indicating a small amount of energy injection during the initial phase of emergence for Type-II ARs. When examining the magnetograms from the early phase of emergence (i.e., helicity accumulation is small), we could see a clear increase in magnetic flux (see Figs. 4 and 5a). Consequently, the minimal helicity injection can be attributed to the relatively diminished velocity field during the initial phase of emergence.
Figure 6 shows another example of AR 12521, and its temporal evolution of the calculated parameters is shown in Fig. 7. This example is a Type-III AR characterized by a distinct reversal in the helicity curve. This implies that the helicity injection is primarily driven by positive helicity in the initial stage followed by a transition to negative helicity dominance. From Fig. 6, we could not get conclusive evidence substantiating helicity reversal. The emergence process exhibits similarities with that of the other two types of ARs, rendering their differentiation based solely on magnetograms a challenge. However, the spatial area and the time interval during emergence are relatively smaller than the other two types. Importantly, discerning the source position of helicity reversal proved to be intricate due to the gauge-dependent nature of helicity flux density, which lacks a distinctly defined physical interpretation (Pariat et al. 2005; Liu et al. 2023). We observed that the energy curve also displays alternating positive and negative phenomena, though the trend is different from the helicity. This phenomenon may be attributed to magnetic cancellations occurring between opposite polarities. Notably, among all the cases, the shear terms (represented by the red curves) contribute significantly to the helicity, as opposed to the emerging terms (represented by the green curves). This indicates that the helicity injection is mainly caused by shearing motions rather than emerging twisted structures, which aligns with previous observations (Liu et al. 2014; Vemareddy 2015; Liokati et al. 2022).
5. Statistical results
Figure 8 displays the scatter plots on various parameters for the three types of ARs. Panel a shows the distribution of accumulated helicity ΔH. It is evident that samples with markedly larger ΔH values (above 3 × 1042 Mx2) mainly belong to Type-I ARs. Furthermore, the ΔH values for Type-I ARs are significantly larger than those of the other two types. The mean ΔH for Type-I ARs is 1.81 × 1042 Mx2, approximately three times larger than that of Type-II ARs and about six times larger than that of Type-III ARs. Panel b presents the distribution of accumulated energy ΔE. Notably, samples with significantly larger ΔE values (above 5 × 1032 erg) are predominantly associated with Type-I ARs. Additionally, the ΔE values for Type-I ARs are also larger than those of the other two types. The mean ΔE for Type-I ARs is 1.51 × 1032 erg, which is approximately twice as large as that of Type-II ARs and Type-III ARs. Panel c shows the distribution of normalized helicity. This quantity allowed us to assess the structural complexity of the ARs rather than focusing solely on their area (Démoulin & Pariat 2009). It is evident that Type-III ARs exhibit the lowest normalized helicity compared to the other two types of ARs. In panel d the distribution of emergence time for the three types of ARs is presented. This parameter indicates the interval from the beginning of emergence to the end of the increase of the magnetic flux calculated by Kutsenko et al. (2019), as mentioned in Sect. 2. The analysis reveals that Type-II ARs exhibit a slightly longer mean emergence time compared to Type-I ARs. However, the ARs with a long emergence time, above 100 h, are almost Type-I ARs. Notably, Type-III ARs have the shortest mean emergence time. This observation could be attributed to the fact that Type-III ARs typically have the smallest spatial range, as is evident from panels a and b. Panels e and f depict the distribution of helicity injection rate and the energy injection rate, respectively. Both panels highlight that Type-I ARs exhibit a higher injection rate, which is similar to the result in panels a and b. Based on the above statistical analyses, we could infer substantial differences among the three types of ARs. Moreover, ΔH and ΔE emerge as good indicators for Type-I ARs, while ΔH/ can be used to distinguish Type-III ARs.
![]() |
Fig. 8. Scatter plots of various parameters for the three types of ARs. The horizontal axes represent the peak unsigned magnetic flux Φp, while the vertical axes display the (a) accumulated helicity ΔH, (b) accumulated energy ΔE, (c) normalized helicity ΔH/ |
Figure 9 presents scatter plots related to the flare index and other parameters. Notably, panel a of this figure reveals that the vast majority (90%) of ARs with a flare index greater than six belong to Type-I ARs. In contrast, only one sample (10%) falls into the category of Type-III. In other words, Type-I ARs are the most likely candidates for becoming highly flare-productive ARs. Moreover, it is noteworthy that these highly flare-productive ARs consistently exhibit substantial helicity accumulation, with values exceeding 2.2 × 1042 Mx2. Furthermore, there are several Type-III ARs with a flare index above four, which indicates that Type-III ARs possess a moderate likelihood of becoming flare-productive ARs. In contrast, Type-II ARs demonstrate a notably limited propensity to develop into flare-productive regions. In panel b, a flare index trend similar to that in panel a can be seen. However, we observed that the normalized helicity ΔH/ does not exhibit a strong correlation with the flare index. It appears that the flare index may not be closely associated with normalized helicity. Panel c illustrates the ΔE − ΔH diagram. A quasi-linear relationship with a Pearson coefficient of 0.8 can be observed between helicity and energy accumulation, which is consistent with previous observations (Tziotziou et al. 2012; Liokati et al. 2022). Using the least-squares best logarithmic fits with a Kolmogorov–Smirnov significance level of about 0.8, the relationship between energy and helicity satisfies the following equation:
![]() |
Fig. 9. Scatter plots of the flare index versus the (a) accumulated helicity ΔH, (b) normalized helicity ΔH/ |
Notably, samples with a high flare index (red and orange scatter points) tend to concentrate in the upper-right corner of the diagram, indicating that flare-productive ARs exhibit a larger helicity and energy accumulation. However, non-flaring ARs (indicated by deep blue scatter points with FI = 0) can also be found in the upper-right corner, suggesting that even with significant helicity and energy accumulation, certain ARs may remain inactive. Moderately active ARs (marked as green and light blue scatter points) can be observed in the lower-left corner, indicating that ARs with a relatively smaller helicity and energy injection can still exhibit moderate flare productivity. In contrast to the ΔE − ΔH diagram with a random time duration in Tziotziou et al. (2012), we strictly defined the starting and ending time of the emergence. In this way, the helicity and energy could reflect the exact value accumulated during the emergence of each AR. Moreover, we used data from SDO/HMI, which are widely used nowadays, to develop the helicity and energy calculation, thus providing a supplement to and validation of their work.
6. Summary and discussion
In this work, we selected 136 samples of emerging ARs and calculated their parameters of magnetic flux, helicity, and energy accumulation during emergence. Through analysis of the distribution of the parameters and their relationship with the flare index, we obtained the following important results:
(i) According to the time evolution of magnetic helicity and unsigned magnetic flux, we classified the emerging ARs into three categories. Type-I ARs exhibit a continuous increase in helicity during emergence closely synchronized with the changes in magnetic flux, with a correlation coefficient above 0.9. In Type-II ARs, the helicity shows minimal growth during the initial phase of emergence followed by a rapid increase after a certain period, with a correlation coefficient below 0.9 compared to the magnetic flux. Type-III ARs display evident injection of opposite helicity, where the accumulated value of helicity with a specific sign constitutes at least 30% of the total accumulated unsigned helicity.
(ii) The three types of ARs exhibit notable differences in terms of flare index, magnetic flux, helicity, normalized helicity, and energy. The most prominent distinction lies in the flare index. We observed that 90% of highly flare-productive ARs (FI ≥ 6) belong to Type-I ARs, suggesting that Type-I ARs are the most likely candidates for achieving a high level of flare productivity. In contrast, Type-II ARs show a very low likelihood of becoming flare productive, while Type-III ARs still maintain a moderate probability of becoming flare productive. Moreover, the ΔH and ΔE values for Type-I ARs are predominantly larger than those of the other two types, and significantly Type-III ARs have the lowest ΔH/ value.
(iii) A highly flare-productive AR is certain to accumulate a significant amount of helicity and energy during its emergence phase. However, despite accumulating substantial energy during emergence, some ARs may remain as non-flaring. We also found that there is a quasi-linear relationship between helicity and energy accumulation, which is in accordance with previous studies (Tziotziou et al. 2012; Liokati et al. 2022).
Our findings have two significant implications. Firstly, it provides valuable insights into the physical mechanisms of flux emergence from the perspectives of magnetic helicity and magnetic energy evolution. It also enhances our understanding of the transport of magnetic helicity from the solar interior to the solar atmosphere. Secondly, it helps establish a potential relationship between magnetic parameters during emergence and subsequent flare activities. Being able to identify whether an AR is flare productive or not when it emerges offers the possibility of early-stage predictions for solar eruptions.
The most distinct characteristic of Type-II ARs is their lag between the helicity and the unsigned flux. They constitute 25% of the total number of emerging ARs, indicating that the phenomenon of delayed helicity injection is widespread. This raises the question of why such a delay in helicity accumulation leads to a flare-poor nature. There are two possible explanations. The first is that the coherency of the emerging flux tube is weak in Type-II ARs. At the initial stage of the emergence, because of the weak coherency of the flux tube, the magnetic helicity can be dispersed due to the convection flows at the photosphere, and thus there is only a small increase in the magnetic helicity. However, for Type-I ARs, which may have a strong coherency, the initial magnetic helicity is not dispersed due to the convection flows at the photosphere, and thus it accumulates from the beginning. In conclusion, the difference in the relations of magnetic helicity with magnetic flux for Type-I and Type-II ARs may reflect the coherency of the emerging magnetic flux tubes. The flux tubes with a strong coherency in Type-I ARs tend to produce more flares compared to those with a weak coherency in Type-II ARs. The second explanation is the lack or minimal injection of helicity during the initial phase of emergence. Previous studies as well as ours have shown that helicity injection is primarily driven by shear motions in the photosphere and sub-photosphere (Liu et al. 2014; Liokati et al. 2022). The absence of helicity injection in the early stages of emergence suggests the absence of significant shear motions at the onset of emergence. Toriumi et al. (2023) performed numerical simulations of the emergence of three types of magnetic flux tubes: non-twisted, weakly twisted, and strongly twisted. They placed the three magnetic flux tubes in the convection zone and observed their emerging process and helicity injection upward, respectively. Their results revealed that in the case of non-twisted flux tubes, the onset of helicity injection is delayed compared to the other two cases (see their Fig. 2). The total helicity accumulation for the non-twisted case is also significantly less than the other two cases. Additionally, their model predicts that ARs formed by non-twisted flux tubes are hard to produce above M-class flares. The scenario of the non-twisted case appears plausible when applied to Type-II ARs in our work. During the initial phase of emergence, the lack of prominent shear motions results in minimal helicity injection. A short time later, the footpoint motions induced by convective turbulence leads to the injection of helicity. These ARs formed by a non-twisted flux tube will have a relatively small helicity accumulation, and it will be difficult for them to be flare productive.
Type-III ARs demonstrate a distinctive characteristic of successive injections of opposite helicity signs. They constitute 22.8% (31 out of 136) of the samples, but only one Type-III case can develop into an AR with a flare index greater than six. Zhang et al. (2006) proposed the concept of an upper bound on the total magnetic helicity that a force-free field can contain. Exceeding this upper bound would initiate a non-equilibrium situation, such as CMEs. In this scenario, the continuous injection of opposite helicity can effectively remove a portion of the accumulated helicity, thus inhibiting CMEs. To support this idea, Vemareddy & Démoulin (2017) and Vemareddy (2021, 2022) conducted detailed case studies of several emerging ARs. They found that the helicity flux changed sign during the evolution of these cases and that the magnetic flux ropes and CMEs are absent. Though they are CME-poor, these cases can be associated with numerous C-class and sub-C-class flaring activities. They concluded that the injection of opposite helicity causes the relaxation of the sheared field through magnetic reconnection, resulting in small-scale flares, but it is difficult to form a large-scale flux rope and a large flare. Their investigations might explain why Type-III ARs have a moderate probability of becoming flare productive. Moreover, the injection of opposite magnetic helicity results in a smaller total helicity accumulation (see Fig. 8d), which is another contributing factor to the inhibition of large flares. There are indeed some studies showing that opposite helicity injection can cause a large flare (Park et al. 2012; Jing et al. 2015); however, the detailed triggering mechanism between opposite helicity injection and flares is debated. These results may explain why four out of our 31 Type-III ARs can have a flare index exceeding four.
Through our study, we conducted a large-sample analysis of the relationship between parameters of emerging active regions and subsequent flare eruptions. However, several limitations must be acknowledged in our approach. Firstly, we did not consider the surrounding conditions of ARs during their emergence. In scenarios where an AR emerges in the presence of a pre-existing neighboring AR, the magnetic field structure can become more complicated, potentially leading to a higher probability of flare production (Toriumi et al. 2017). Thus, we focused solely on calculating the parameters of the AR within its emergence region. Secondly, we did not conduct a detailed analysis of how the energy and helicity accumulations trigger and contribute to flare eruptions. Our investigation only focused on a general examination of the relationship between the calculated parameters and the flare index. In future, we would like to investigate the long-term helicity evolution of ARs and their following solar eruptions in order to establish the relationship between AR helicity evolution and flare/CME productivity. Moreover, we will likely conduct some simulations focusing on the eruption potential of emerging ARs, which can provide possibilities for early-stage predictions of solar eruptions.
Acknowledgments
We thank Dr. Jingxiu Wang and Dr. Shin Toriumi for fruitful suggestions. This work is supported by the Strategic Priority Research Program of the Chinese Academy of Sciences (Grant No. XDB0560000 and XDB41000000), the National Key R&D Program of China (2019YFA0405000, 2021YFA1600500, 2022YFF0503800) and the National Natural Science Foundations of China (11973056, 12222306, and 12273060).
References
- Abramenko, V. I. 2005, ApJ, 629, 1141 [NASA ADS] [CrossRef] [Google Scholar]
- Abramenko, V. I., Kutsenko, A. S., Tikhonova, O. I., & Yurchyshyn, V. B. 2017, Sol. Phys., 292, 48 [NASA ADS] [CrossRef] [Google Scholar]
- Archontis, V., Tsinganos, K., & Gontikakis, C. 2010, A&A, 512, L2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Berger, M. A. 1999, Geophysical Monograph Ser. (Washington, D.C.), 111, 1 [NASA ADS] [Google Scholar]
- Berger, M. A., & Field, G. B. 1984, J. Fluid Mech., 147, 133 [NASA ADS] [CrossRef] [Google Scholar]
- Bobra, M. G., Sun, X., Hoeksema, J. T., et al. 2014, Sol. Phys., 289, 3549 [Google Scholar]
- Chen, Y., Tian, H., Zhu, X., et al. 2019, Sci. China E: Technol. Sci., 62, 1555 [CrossRef] [Google Scholar]
- Cheung, M. C. M., & Isobe, H. 2014, Liv. Rev. Sol. Phys., 11, 3 [Google Scholar]
- Cheung, M. C. M., Rempel, M., Title, A. M., & Schüssler, M. 2010, ApJ, 720, 233 [Google Scholar]
- Chintzoglou, G., Zhang, J., Cheung, M. C. M., & Kazachenko, M. 2019, ApJ, 871, 67 [NASA ADS] [CrossRef] [Google Scholar]
- Dalmasse, K., Aulanier, G., Démoulin, P., et al. 2015, ApJ, 810, 17 [NASA ADS] [CrossRef] [Google Scholar]
- Démoulin, P., & Berger, M. A. 2003, Sol. Phys., 215, 203 [CrossRef] [Google Scholar]
- Démoulin, P., & Pariat, E. 2009, AdSpR, 43, 1013 [Google Scholar]
- Duan, Y., Shen, Y., Zhou, X., et al. 2022, ApJ, 926, L39 [NASA ADS] [CrossRef] [Google Scholar]
- Duan, A., Jiang, C., & Feng, X. 2023, ApJ, 945, 102 [NASA ADS] [CrossRef] [Google Scholar]
- Fan, Y. 2021, Liv. Rev. Sol. Phys., 18, 5 [NASA ADS] [CrossRef] [Google Scholar]
- Fan, Y., & Fang, F. 2014, ApJ, 789, 35 [NASA ADS] [CrossRef] [Google Scholar]
- Fournier, Y., Arlt, R., Ziegler, U., & Strassmeier, K. G. 2017, A&A, 607, A1 [Google Scholar]
- Fu, Y., & Welsch, B. T. 2016, Sol. Phys., 291, 383 [NASA ADS] [CrossRef] [Google Scholar]
- Gupta, M., Thalmann, J. K., & Veronig, A. M. 2021, A&A, 653, A69 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Janvier, M., Savcheva, A., Pariat, E., et al. 2016, A&A, 591, A141 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Jing, J., Xu, Y., Lee, J., et al. 2015, RAA, 15, 1537 [NASA ADS] [Google Scholar]
- Kazachenko, M. D., Lynch, B. J., Savcheva, A., Sun, X., & Welsch, B. T. 2022, ApJ, 926, 56 [CrossRef] [Google Scholar]
- Kontogiannis, I., Georgoulis, M. K., Guerra, J. A., Park, S.-H., & Bloomfield, D. S. 2019, Sol. Phys., 294, 130 [NASA ADS] [CrossRef] [Google Scholar]
- Korsós, M. B., Romano, P., Morgan, H., et al. 2020, ApJ, 897, L23 [Google Scholar]
- Kutsenko, A. S., Abramenko, V. I., & Pevtsov, A. A. 2019, MNRAS, 484, 4393 [NASA ADS] [CrossRef] [Google Scholar]
- Kutsenko, A. S., Abramenko, V. I., & Kutsenko, O. K. 2021, MNRAS, 501, 6076 [NASA ADS] [CrossRef] [Google Scholar]
- LaBonte, B. J., Georgoulis, M. K., & Rust, D. M. 2007, ApJ, 671, 955 [CrossRef] [Google Scholar]
- Li, T., & Zhang, J. 2015, ApJ, 804, L8 [NASA ADS] [CrossRef] [Google Scholar]
- Li, T., Sun, X., Hou, Y., et al. 2022, ApJ, 926, L14 [NASA ADS] [CrossRef] [Google Scholar]
- Liokati, E., Nindos, A., & Liu, Y. 2022, A&A, 662, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Liu, Y., & Schuck, P. W. 2012, ApJ, 761, 105 [NASA ADS] [CrossRef] [Google Scholar]
- Liu, Y., Hoeksema, J. T., Scherrer, P. H., et al. 2012, Sol. Phys., 279, 295 [Google Scholar]
- Liu, Y., Hoeksema, J. T., Bobra, M., et al. 2014, ApJ, 785, 13 [NASA ADS] [CrossRef] [Google Scholar]
- Liu, Y., Welsch, B. T., Valori, G., et al. 2023, ApJ, 942, 27 [NASA ADS] [CrossRef] [Google Scholar]
- Magara, T. 2001, ApJ, 549, 608 [CrossRef] [Google Scholar]
- Moraitis, K., Sun, X., Pariat, É., & Linan, L. 2019, A&A, 628, A50 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Murray, M. J., Hood, A. W., Moreno-Insertis, F., Galsgaard, K., & Archontis, V. 2006, A&A, 460, 909 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Pariat, E., Démoulin, P., & Berger, M. A. 2005, A&A, 439, 1191 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Pariat, E., Leake, J. E., Valori, G., et al. 2017, A&A, 601, A125 [CrossRef] [EDP Sciences] [Google Scholar]
- Park, S.-H., Chae, J., & Wang, H. 2010, ApJ, 718, 43 [NASA ADS] [CrossRef] [Google Scholar]
- Park, S.-H., Cho, K.-S., Bong, S.-C., et al. 2012, ApJ, 750, 48 [NASA ADS] [CrossRef] [Google Scholar]
- Park, S.-H., Leka, K. D., & Kusano, K. 2021, ApJ, 911, 79 [NASA ADS] [CrossRef] [Google Scholar]
- Pevtsov, A. A., Maleev, V. M., & Longcope, D. W. 2003, ApJ, 593, 1217 [NASA ADS] [CrossRef] [Google Scholar]
- Pevtsov, A. A., Berger, M. A., Nindos, A., Norton, A. A., & van Driel-Gesztelyi, L. 2014, Space Sci. Rev., 186, 285 [NASA ADS] [CrossRef] [Google Scholar]
- Prior, C., & MacTaggart, D. 2019, J. Plasma Phys., 85, 775850201 [NASA ADS] [CrossRef] [Google Scholar]
- Scherrer, P. H., Schou, J., Bush, R. I., et al. 2012, Sol. Phys., 275, 207 [Google Scholar]
- Schuck, P. W. 2008, ApJ, 683, 1134 [NASA ADS] [CrossRef] [Google Scholar]
- Soós, S., Korsós, M. B., Morgan, H., & Erdélyi, R. 2022, ApJ, 925, 129 [CrossRef] [Google Scholar]
- Stein, R. F., & Nordlund, Å. 2012, ApJ, 753, L13 [NASA ADS] [CrossRef] [Google Scholar]
- Sun, X. 2013, ArXiv e-prints [arXiv:1309.2392] [Google Scholar]
- Sun, X., Bobra, M. G., Hoeksema, J. T., et al. 2015, ApJ, 804, L28 [Google Scholar]
- Sun, Z., Tian, H., Chen, P. F., et al. 2022, ApJ, 939, L18 [NASA ADS] [CrossRef] [Google Scholar]
- Sun, Z., Li, T., Tian, H., et al. 2023, ApJ, 953, 148 [CrossRef] [Google Scholar]
- Thalmann, J. K., Moraitis, K., Linan, L., et al. 2019, ApJ, 887, 64 [Google Scholar]
- Toriumi, S., & Park, S. H. 2022, ArXiv e-prints [arXiv:2204.06010] [Google Scholar]
- Toriumi, S., & Takasao, S. 2017, ApJ, 850, 39 [NASA ADS] [CrossRef] [Google Scholar]
- Toriumi, S., & Wang, H. 2019, Liv. Rev. Sol. Phys., 16, 3 [Google Scholar]
- Toriumi, S., Iida, Y., Kusano, K., Bamba, Y., & Imada, S. 2014, Sol. Phys., 289, 3351 [CrossRef] [Google Scholar]
- Toriumi, S., Schrijver, C. J., Harra, L. K., Hudson, H., & Nagashima, K. 2017, ApJ, 834, 56 [Google Scholar]
- Toriumi, S., Hotta, H., & Kusano, K. 2023, Sci. Rep., 13, 8994 [NASA ADS] [CrossRef] [Google Scholar]
- Trujillo Bueno, J., Shchukina, N., & Asensio Ramos, A. 2004, Nature, 430, 326 [Google Scholar]
- Tziotziou, K., Georgoulis, M. K., & Raouafi, N.-E. 2012, ApJ, 759, L4 [NASA ADS] [CrossRef] [Google Scholar]
- Valori, G., Pariat, E., Anfinogentov, S., et al. 2016, Space Sci. Rev., 201, 147 [Google Scholar]
- van Driel-Gesztelyi, L., & Green, L. M. 2015, Liv. Rev. Sol. Phys., 12, 1 [Google Scholar]
- Vemareddy, P. 2015, ApJ, 806, 245 [NASA ADS] [CrossRef] [Google Scholar]
- Vemareddy, P. 2021, MNRAS, 507, 6037 [NASA ADS] [CrossRef] [Google Scholar]
- Vemareddy, P. 2022, MNRAS, 516, 158 [Google Scholar]
- Vemareddy, P., & Démoulin, P. 2017, A&A, 597, A104 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Wang, J., Yan, X., Kong, D., et al. 2021, A&A, 652, A55 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Wang, Q., Yang, S., Zhang, M., & Yang, X. 2022, ApJ, 929, 122 [NASA ADS] [CrossRef] [Google Scholar]
- Wang, Q., Zhang, M., Yang, S., Yang, X., & Zhu, X. 2023, RAA, 23, 095025 [NASA ADS] [Google Scholar]
- Welsch, B. T., Li, Y., Schuck, P. W., & Fisher, G. H. 2009, ApJ, 705, 821 [NASA ADS] [CrossRef] [Google Scholar]
- Woltjer, L. 1958, Proc. Natl. Acad. Sci., 44, 489 [NASA ADS] [CrossRef] [Google Scholar]
- Yang, S., Zhang, H., & Büchner, J. 2009, A&A, 502, 333 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Yang, S., Büchner, J., Santos, J. C., & Zhang, H. 2013, Sol. Phys., 283, 369 [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]
- Zhang, H. 2001, ApJ, 557, L71 [NASA ADS] [CrossRef] [Google Scholar]
- Zhang, M., & Flyer, N. 2008, ApJ, 683, 1160 [NASA ADS] [CrossRef] [Google Scholar]
- Zhang, M., & Low, B. C. 2001, ApJ, 561, 406 [NASA ADS] [CrossRef] [Google Scholar]
- Zhang, M., Flyer, N., & Low, B. C. 2006, ApJ, 644, 575 [NASA ADS] [CrossRef] [Google Scholar]
Appendix A: Calculated parameters of the dataset
Emerging AR samples and their corresponding parameters.
All Tables
All Figures
![]() |
Fig. 1. Unsigned magnetic flux Φ (blacks curve) and accumulated magnetic helicity ΔH (orange curves) evolutions for nine examples. The top, middle, and bottom rows represent the Type-I, Type-II, and Type-III ARs, respectively. The decreasing helicity indicates a negative helicity injection. The correlation coefficients between the two curves are marked in the upper-left corner of each panel. |
In the text |
![]() |
Fig. 2. Selected HMI images of the normal component, Bn, of the photospheric field of the Type-I AR 11620 showing the emergence process. The length of the horizontal white line corresponds to 100 arcsec. |
In the text |
![]() |
Fig. 3. Time profiles of the magnetic flux, helicity, and energy for the Type-I AR 11620. Panel (a): time evolution of the unsigned magnetic flux Φ (black) and accumulated energy ΔE (blue). Panel (b): time profile of the helicity flux, d(H)/dt (black). The red and green curves represent the shear and emergence terms, respectively. Panel (c): time profiles of accumulated helicity ΔH. The red and green curves indicate the helicity of the emergence term and shear term, respectively. |
In the text |
![]() |
Fig. 4. Same as Fig. 2 but for the Type-II AR 11768. |
In the text |
![]() |
Fig. 5. Same as Fig. 3 but for the Type-II AR 11768. |
In the text |
![]() |
Fig. 6. Same as Fig. 2 but for the Type-III AR 12521. |
In the text |
![]() |
Fig. 7. Same as Fig. 3 but for the Type-III AR 12521. |
In the text |
![]() |
Fig. 8. Scatter plots of various parameters for the three types of ARs. The horizontal axes represent the peak unsigned magnetic flux Φp, while the vertical axes display the (a) accumulated helicity ΔH, (b) accumulated energy ΔE, (c) normalized helicity ΔH/ |
In the text |
![]() |
Fig. 9. Scatter plots of the flare index versus the (a) accumulated helicity ΔH, (b) normalized helicity ΔH/ |
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.