Abstract
The helium abundance, defined as AHe = nHe/nH × 100, is ∼8.5 in the photosphere and seldom exceeds 5 in fast solar wind. Previous statistics have demonstrated that AHe in slow solar wind correlates tightly with sunspot number. However, less attention is paid to the solar cycle dependence of AHe within interplanetary coronal mass ejections (ICMEs) and comparing the AHe characteristics of ICMEs and solar wind. In this paper we conduct a statistical comparison of helium abundance between ICMEs and solar wind near 1 au with observations of the Advanced Composition Explorer from 1998 to 2019 and find that the ICME AHe also exhibits the obvious solar cycle dependence. Meanwhile, we find that the AHe is obviously higher within ICMEs compared to solar wind, and the means within 37% and 12% of ICMEs exceed 5 and 8.5, respectively. It is interesting to answer where and how the high helium abundance originates. Our statistics demonstrate that 21% (3%) of ICME (slow wind) AHe data points exceed 8.5 around solar maximum, which decreases dramatically near minimum, while no such high AHe values appear in the fast wind throughout the whole solar cycle. This indicates that the high AHe (e.g., >8.5) emanates from active regions as more ICMEs and slow wind originate from active regions around maximum, and it supports that both active regions and quiet-Sun regions are the sources of slow wind. We suggest that the high AHe from active regions could be explained by means of the magnetic loop confinement model and/or photoionization effect.
Export citation and abstract BibTeX RIS
![](https://content.cld.iop.org/journals/0004-637X/925/2/137/revision1/apjac3bbflicense.gif)
Original content from this work may be used under the terms of the Creative Commons Attribution 4.0 licence. Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI.
1. Introduction
Coronal mass ejections (CMEs) are one of the most energetic eruptions in the solar system and can cause disastrous space weather effects (e.g., Gosling et al. 1991; Xu et al. 2019). It is widely accepted that CMEs result from eruption of magnetic flux ropes (e.g., Chen 2011), which can form prior to (e.g., Cheng et al. 2011; Patsourakos et al. 2013) and during (e.g., Song et al. 2014; Jiang et al. 2021) eruptions. The composition (including both the charge states and elemental abundances of heavy ions) provides an important avenue to analyze their eruption process (Lepri & Zurbuchen 2004; Lynch et al. 2011; Gruesbeck et al. 2011, 2012; Song et al. 2015, 2016, 2017; Gu et al. 2020; Song & Yao 2020; Song et al. 2021a) and supports that more CMEs come from active regions around solar maximum, whose heavy ions possess higher charge states and enriched relative abundances (Song et al. 2021b).
Solar wind is a stream of charged particles released from the corona and can be divided into fast and slow streams taking speed as a criterion. It is generally accepted that fast wind originates from coronal holes, and slow wind can emanate from active regions and/or quiet-Sun regions (e.g., Zhao et al. 2017; Fu et al. 2018). Two classes of models have been proposed to explain the origin of solar wind, including the wave-turbulence-driven models, which suggest that the solar wind is released along open magnetic field lines directly (Hollweg 1986; Cranmer et al. 2007; Verdini et al. 2009), and the reconnection loop opening (RLO) models, which propose that the stream escapes through magnetic reconnections between open and closed magnetic field lines (Fisk et al. 1999; Fisk 2003; Woo et al. 2004) and work for solar wind originating from the closed-field region. These models can be examined through analyzing the solar wind composition.
The helium abundance, defined as AHe = nHe/nH × 100, can be employed to analyze the source regions and release processes of both CMEs and solar wind. As hydrogen and helium are the first and second most abundant elements within both interplanetary CMEs (ICMEs) and solar wind (Aellig et al. 2001; Kasper et al. 2007, 2012; McIntosh et al. 2011; Alterman & Kasper 2019; Alterman et al. 2021), respectively, they can be measured in situ with relatively higher cadence and accuracy compared to the other heavier ions.
AHe is ∼8.5 in the photosphere (Asplund et al. 2009) and seldom exceeds 5 in fast solar wind (Kasper et al. 2007, 2012) partly because the first ionization potential (FIP) of helium (24.6 eV) is higher than that of hydrogen (13.6 eV). This makes it harder for the neutral helium atoms to be ionized through thermal ionization and flow out than the hydrogen in the lower solar atmosphere, resulting in the AHe depletion in the corona and solar wind (Laming 2004, 2015). Earlier statistics also demonstrated that the AHe in slow wind correlates tightly with sunspot number (SSN; Aellig et al. 2001; Kasper et al. 2007, 2012; McIntosh et al. 2011; Alterman & Kasper 2019) and can be used to herald the solar cycle onset (Alterman et al. 2021). Compared to the solar wind, the statistical studies reveal that the AHe is obviously higher within ICMEs (Owens 2018; Huang et al. 2020).
It is interesting to explore where and how the high AHe of ICMEs originates. Fu et al. (2020) suggested that the chromospheric evaporation during solar eruptions can supply plasma with high AHe to CMEs through a case study, while this scenario cannot consistently explain that ICMEs usually exhibit the obvious FIP effect, i.e., possessing higher Fe/O, Mg/O, and Si/O than the photosphere and chromosphere (Zurbuchen et al. 2016; Owens 2018; Huang et al. 2020). Therefore, it might not work well for most events. Previous studies showed that the AHe is a linear function of solar wind speed around solar minimum, and the correlation between AHe and SSN is strongest at low speed, while the positive correlation decreases quickly with increasing speed (Kasper et al. 2007, 2012), which can be explained qualitatively by means of the coronal loop confinement model (Fisk 2003; Woo et al. 2004; Endeve et al. 2005; Kasper et al. 2007). This model connects the confinement time of plasma in a loop to its physical size, and these factors in turn influence the temperature profile, gravitational settling, and plasma composition (Kasper et al. 2007). Besides, Song et al. (2021b) reported that Ne/O exhibits opposite solar cycle dependence within ICMEs and solar wind, and they suggested that the photoionization might be responsible for the higher Ne/O ratios within ICMEs around solar maximum.
We anticipate that a statistical comparison of AHe between ICMEs and solar wind could shed more light on whether the loop confinement model and photoionization can be employed to answer the origin of high AHe within ICMEs. This is the major motivation for us to conduct this study. The paper is organized as follows. In Section 2, we introduce the data and ICME catalogs. The statistical results and comparisons are displayed in Section 3. Section 4 provides a discussion and summary.
2. Data and ICME Catalogs
The helium abundance (AHe) data in this study are provided by the Solar Wind Electron Proton Alpha Monitor (SWEPAM; McComas et al. 1998) experiment aboard the Advanced Composition Explorer (ACE), which has orbited around the L1 point since its launch in 1997 (Stone et al. 1998). The SWEPAM observations are made with independent electron and ion instruments through electrostatic analyzers. With the fan-shaped fields of view, SWEPAM can sweep out all pertinent look directions as ACE spins and provide 64 s resolution of the full electron and ion distribution functions (McComas et al. 1998). All of these data can be downloaded at the ACE science center. 5 The yearly SSNs are obtained at the Solar Influence Data Center of the Royal Observatory of Belgium. 6
ICMEs can be identified with different physical features, such as enhanced magnetic field intensity, smoothly changing field direction, low proton temperature, low plasma β, and so on (e.g., Wu & Lepping 2011). Researchers have identified ICMEs through in situ data near 1 au and offered several complete and reliable catalogs using measurements of the ACE (Richardson & Cane 2010; RC catalog), WIND (Chi et al. 2016; Nieves-Chinchilla et al. 2018), and the Solar Terrestrial Relations Observatory (Jian et al. 2018). Every catalog lists the ejecta boundaries of each ICME, which are adopted to analyze the helium abundance within ICMEs. In this paper, the RC catalog 7 is used, as the AHe is provided by the SWEPAM aboard the ACE.
3. Statistical Results and Comparisons
The RC catalog lists a total of 487 ICMEs from 1996 to 2019, and the yearly ICME numbers are displayed with white histograms in Figure 1. To demonstrate the solar cycle dependence of ICME numbers, the yearly SSNs are also plotted here with the black line. The Spearman rank cross-correlation coefficient (ρ) between yearly ICME numbers and SSNs is calculated with the IDL code r_correlate.pro, and ρ > 0.6 implies the meaningful cross-correlation (e.g., Alterman & Kasper 2019). The ρ between yearly numbers of ICMEs and sunspots is 0.85, which is also shown in Figure 1 and illustrates a strong correlation between them.
Figure 1. Yearly numbers of ICMEs at L1 point and sunspots from 1996 to 2019. The white histograms represent the ICME numbers, with light gray depicting those with available AHe. The dark-gray and black regions indicate the ICME numbers with AHe exceeding 5 and 8.5, respectively. The black line displays the yearly SSNs. The Spearman rank cross-correlation coefficient (ρ) between yearly numbers of ICMEs and sunspots is also presented.
Download figure:
Standard image High-resolution imageThe SWEPAM began to provide AHe from 1998 February onward, and data gaps appear occasionally during the normal observations. Thus, the AHe is available just for 420 ICMEs that are depicted with light gray in Figure 1, including 155 (37%) and 50 (12%) ICMEs with their average AHe exceeding 5 and 8.5, which are displayed with the dark-gray and black histograms, respectively. Figure 1 shows that ICMEs with high AHe values mainly appear during solar maximum and descending phases, similar to the relative abundances of the other heavier ions (Song et al. 2021b).
To analyze the solar cycle dependence of ICME AHe, we first calculate the average value of AHe within each ICME with the 64 s resolution SWEPAM data. The average duration of ICMEs is ∼20 hr (e.g., Wu & Lepping 2011), which means that one ICME contains ∼1125 data points on average. Then, the yearly mean, median, and standard deviation of all the ICMEs in each year are calculated. In the meantime, we also calculate the yearly values for slow and fast solar wind. The interaction between fast streams, emanating from coronal holes, and slow streams can produce the corotating interaction regions (CIRs) in interplanetary space (Jian et al. 2006), where the solar wind speeds can be changed obviously. For example, the original slow (<400 km s−1) and fast (>600 km s−1) wind is accelerated and decelerated to ∼500 km s−1 (see Figure 2 in Jian et al. 2006 for an example). If we simply set 400 or 500 km s−1 as the critical value to differentiate fast and slow wind, the original slow/fast wind might be counted as fast/slow wind given the existence of CIRs. Therefore, we set our slow and fast wind threshold at <400 km s−1 and >600 km s−1, i.e., excluding the wind between 400 and 600 km s−1, which could guarantee more that pure slow or fast streams are investigated without contamination from each other. Table 1 lists all of the calculated results for ICMEs, slow wind, and fast wind.
Table 1. The Yearly Sunspot Number, Average Value, Median, and Standard Deviation of AHe within ICMEs, Slow Solar Wind (SSW), and Fast Solar Wind (FSW) in Each Year from 1998 to 2019
Year | 1998 | 1999 | 2000 | 2001 | 2002 | 2003 | 2004 | 2005 | 2006 | 2007 | 2008 | 2009 | 2010 | 2011 | 2012 | 2013 | 2014 | 2015 | 2016 | 2017 | 2018 | 2019 | |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
SSN | 88.3 | 136.3 | 173.9 | 170.4 | 163.6 | 99.3 | 65.3 | 45.8 | 24.7 | 12.6 | 4.2 | 4.8 | 24.9 | 80.8 | 84.5 | 94.0 | 113.3 | 69.8 | 39.8 | 21.7 | 7.0 | 3.6 | |
ICME | mean | 3.918 | 3.968 | 5.682 | 5.294 | 6.073 | 6.194 | 5.016 | 5.848 | 2.902 | 2.856 | 1.622 | 1.266 | 2.986 | 3.532 | 4.831 | 5.082 | 5.540 | 4.577 | 3.084 | 3.477 | 1.919 | 2.045 |
median | 3.476 | 3.660 | 4.860 | 4.919 | 5.479 | 5.474 | 4.252 | 5.875 | 1.786 | 2.856 | 1.622 | 0.731 | 2.056 | 2.832 | 4.023 | 4.176 | 4.585 | 3.686 | 2.266 | 3.165 | 1.655 | 0.931 | |
stddev | 2.313 | 2.458 | 3.130 | 2.608 | 3.487 | 3.399 | 2.952 | 3.605 | 2.039 | 0.463 | 1.206 | 1.293 | 2.179 | 2.645 | 3.475 | 3.686 | 3.269 | 2.434 | 2.088 | 1.742 | 1.183 | 1.788 | |
SSW | mean | 2.651 | 3.252 | 4.092 | 3.705 | 4.576 | 4.776 | 3.486 | 2.884 | 2.171 | 1.950 | 1.345 | 1.382 | 2.110 | 3.107 | 3.358 | 3.007 | 3.844 | 3.550 | 2.383 | 2.382 | 1.759 | 1.528 |
median | 2.280 | 2.910 | 3.840 | 3.340 | 4.310 | 4.230 | 2.990 | 2.610 | 1.920 | 1.620 | 0.960 | 1.090 | 1.540 | 2.950 | 3.100 | 2.810 | 3.770 | 3.490 | 2.150 | 2.190 | 1.510 | 1.070 | |
stddev | 1.713 | 1.991 | 2.022 | 2.183 | 2.404 | 2.763 | 2.207 | 1.662 | 1.357 | 1.413 | 1.227 | 1.080 | 1.952 | 1.768 | 1.922 | 1.673 | 1.592 | 1.608 | 1.326 | 1.300 | 1.208 | 1.270 | |
FSW | mean | 2.284 | 2.425 | 2.885 | 2.098 | 2.362 | 2.171 | 2.567 | 2.514 | 2.917 | 3.118 | 2.372 | 3.101 | 4.036 | 4.099 | 4.084 | 3.227 | 3.201 | 3.510 | 2.489 | 3.219 | 2.836 | 2.717 |
median | 2.240 | 2.310 | 2.610 | 1.450 | 2.290 | 1.940 | 2.390 | 2.330 | 2.870 | 3.030 | 2.330 | 3.070 | 4.000 | 4.010 | 4.000 | 3.090 | 3.060 | 3.610 | 2.360 | 3.150 | 2.760 | 2.635 | |
stddev | 1.509 | 1.292 | 1.956 | 2.105 | 1.029 | 1.130 | 1.573 | 1.114 | 1.341 | 1.116 | 1.063 | 0.576 | 1.766 | 1.505 | 0.850 | 1.540 | 1.184 | 1.351 | 1.019 | 1.009 | 1.298 | 0.904 |
Download table as: ASCIITypeset image
Figure 2 displays the yearly means of AHe within ICMEs (pink), slow wind (blue), and fast wind (red), along with the SSNs (black, also listed in Table 1), and the corresponding Spearman rank coefficients of correlation between AHe and SSN are 0.85, 0.90, and −0.22, which demonstrate quantitatively that the helium abundances within ICMEs also possess the solar cycle dependence, while fast wind does not. As only a few ICMEs were detected in each year around solar minimum, we do not divide the yearly ICMEs into magnetic clouds and noncloud events, or fast and slow ones. Owens (2018) reported that magnetic clouds have higher AHe than noncloud ICMEs. In the meantime, the median and standard deviation are not plotted in Figure 2 for clarity, which could be found in Table 1.
Figure 2. The solar cycle dependence of AHe within ICMEs (pink), slow solar wind (SSW; blue), and fast solar wind (FSW; red). The Spearman rank cross-correlation coefficients (ρ) between yearly AHe and SSNs (depicted with the black line) are also presented.
Download figure:
Standard image High-resolution imageHere the strong correlation between AHe and SSN in the slow wind is consistent with previous statistical results with WIND measurements (Kasper et al. 2007, 2012; Alterman & Kasper 2019; Alterman et al. 2021). One different point is that AHe of solar wind from WIND has comparable values during maxima of solar cycles 23 and 24, although the SSN amplitude of cycle 24 is less compared to the previous cycle (see Figure 1 in Alterman & Kasper 2019). While the ACE measurements demonstrate that the AHe values of both slow wind and ICMEs around cycle 24 maximum (2012–2015) decreased slightly compared to cycle 23 (2000–2003), consistent with the expectation of SSN variation (Alterman et al. 2021).
The quantitative comparison between the two solar maxima shows that the average AHe of slow wind (ICMEs) decreases by 16.3% (11.7%), i.e., from 4.16 (6.16) to 3.48 (5.44), and the t-test (see e.g., Song et al. 2021b) results illustrate that the variation is significant at the 99% confidence level. The discrepancy of AHe variation trends between WIND (Alterman & Kasper 2019) and ACE should mainly result from different selections of velocity intervals. Alterman & Kasper (2019) presented the AHe of solar wind with velocity between 312 and 574 km s−1, while our velocity range is less than 400 km s−1 as mentioned. We can get the same result that the helium abundances of solar wind are comparable during the last two solar maxima if choosing the same velocity interval with Alterman & Kasper (2019) using ACE data (not shown).
Figure 3 presents a series of normalized distributions of AHe data points with 64 s resolution for ICMEs (top), slow wind (middle), and fast wind (bottom). Figure 2 has shown that the AHe of ICMEs and slow wind can change significantly from solar maximum to minimum, while we find that no obvious differences exist between the overall distributions of AHe for solar cycles 23 and 24. Therefore, the two solar cycles (1998–2019) are displayed together in the left panels, and the two solar maximum (2000–2003 and 2012–2015) and minimum (2006–2009 and 2016–2019) phases are presented in the middle and right panels, respectively. The red vertical lines denote the photospheric value of AHe (i.e., 8.5). The quantitative analysis shows that 21% of data points within ICMEs possess a high AHe value (>8.5) around solar maximum, which decreases dramatically to 2.3% around minimum. For the solar wind, only slow wind around maximum has 2.7% of data points that exceed 8.5, and in other cases the high AHe (>8.5) data points could be negligible.
Figure 3. The normalized distributions of AHe data points with 64 s resolution within ICMEs (top), slow solar wind (SSW, middle), and fast solar wind (FSW, bottom). The two solar cycles (1998–2019) are displayed together in the left panels, and the two solar maximum (2000–2003 and 2012–2015) and minimum (2006–2009 and 2016–2019) phases are presented in the middle and right panels, respectively. The corresponding mean and median are given in each panel, and the vertical red lines depict the AHe value in the photosphere.
Download figure:
Standard image High-resolution imageThe mean and median of the AHe distributions are also displayed in each panel of Figure 3. It shows that ICMEs possess the highest AHe compared to slow and fast wind, consistent with Figure 2. Focusing on the solar wind alone, the AHe is more enriched in fast wind compared to slow wind around solar minimum, while slow wind possesses higher AHe than fast wind around maximum. The AHe of fast wind does not exhibit the solar cycle dependence, while slow wind possesses obvious higher AHe around maximum compared to minimum. This implies that slow winds originating from the active region and quiet-Sun region have different properties, as slow winds mostly come from the quiet-Sun region around solar minimum, and more fractions of slow wind come from the active region around maximum.
4. Discussion and Summary
Our study demonstrates that AHe in both ICMEs and slow wind exhibits the positive correlation with SSNs, which indicates that the high AHe (e.g., >8.5) emanates from active regions as more ICMEs and slow wind originate from active regions around solar maximum. In the meantime, no high AHe (>8.5) data points existing in fast wind throughout a solar cycle implies that coronal holes do not emanate plasmas with enriched helium. This allows us to infer how the high AHe originates through comparing the characteristics of active regions and coronal holes.
It is generally accepted that the dominated magnetic field lines in active regions and coronal holes are closed and open, respectively. Some studies displayed that the high AHe could be explained through the coronal loop confinement model (Fisk 2003; Woo et al. 2004; Endeve et al. 2005; Kasper et al. 2007). This model suggested that the properties of coronal loop plasmas can vary over their confinement duration before the plasmas are released through magnetic reconnection between loop and open field lines (Fisk 2003). As mentioned, the confinement time of plasmas in the coronal loop is correlated with its physical size, and these factors in turn influence the temperature profile, gravitational settling, and composition of the plasma. For example, simulations (Killie et al. 2005) displayed that the AHe in closed-field regions can increase obviously in several days owing to the thermal force, which is caused by the energy dependence of the Coulomb cross section. The thermal force can be large in the transition region and seeks to push helium ions toward the corona, causing the buildup of helium in the coronal loop (Killie et al. 2005). These high AHe plasma confined in the coronal loop can become part of CME materials through reconnection occurring in the current sheet beneath CMEs during solar eruptions (Lin & Forbes 2000; Lin et al. 2004; Song et al. 2016), corresponding to the enriched helium within ICMEs. Also, they can be released into slow wind by reconnection between closed loop and open field lines.
The other potential mechanism for high AHe within ICMEs is related to photoionization, similar to explaining high Ne/O ratios within flare regions and ICMEs (Feldman et al. 2005; Zurbuchen et al. 2016). According to the FIP effect (Laming 2004, 2015), the Ne/O should decrease in the corona and ICMEs compared to the photosphere, as the FIP of Ne (21.6 eV) is higher than that of O (13.6 eV). However, measurements showed the opposite case within flare regions and ICMEs (Feldman et al. 2005; Song et al. 2021b). Shemi (1991) has suggested that pre-flare soft X-rays can penetrate through the chromosphere and create a slab-like region with photoionization ratios at the chromosphere base. As the photoionization cross section ratio of Ne and O is 9:4 (Yeh & Lindau 1985), the ionization ratio of Ne is higher than that of O. This leads to plasma being transported into the corona with higher Ne/O, corresponding to the enriched Ne/O in flare regions (Schmelz 1993) and ICMEs (Song et al. 2021b). Meanwhile, more CMEs are accompanied by energetic flares during solar maximum, which indicates that the photoionization might play a more important role around maximum compared to minimum, further enhancing the solar cycle dependence of ICME Ne/O. Similarly, the photoionization cross section of He is larger than that of H (Yeh & Lindau 1985); thus, it is straightforward to speculate that the photoionization leads to the enhanced AHe within ICMEs and flare regions (Feldman et al. 2005). Simulations are needed to examine this conjecture.
In this paper we conducted a statistical comparison of helium abundance between ICMEs, slow solar wind (<400 km s−1), and fast solar wind (>600 km s−1) near 1 au with ACE observations from 1998 to 2019, covering solar cycles 23 and 24. The statistics demonstrated that the AHe of ICMEs and slow wind exhibits the obvious solar cycle dependence, and the AHe within ICMEs is highest compared to slow and fast winds. The averages within 37% and 12% of ICMEs exceed 5 and 8.5, respectively. About 21% (3%) of AHe data points of ICMEs (slow wind) exceed 8.5 around solar maximum, which decreases dramatically around minimum, while no such high AHe values appear in the fast wind throughout the solar cycle. Focusing on the solar wind alone, the AHe is more enriched in fast wind compared to slow wind around solar minimum, while slow wind possesses higher AHe than fast wind around maximum. This indicates that high AHe originates from active regions. Two possible mechanisms are discussed to explain the high AHe within both slow wind and ICMEs around solar maximum.
We thank the anonymous referee for the comments and suggestions that helped to improve the original manuscript. We acknowledge the use of the ICME catalog provided by Richardson & Cane. All the SWICS data are downloaded from the ACE science center. H.S. thanks Drs. Bo Li, Liang Zhao, and Liang Guo (a statistician at Shandong University) for their useful discussions. This work is supported by the NSFC grants U2031109, 11790303 (11790300), and 12073042. H.S. is also supported by the CAS grants XDA-17040507 and the open research program of the CAS Key Laboratory of Solar Activity KLSA202107.
Footnotes
- 5
- 6
- 7