Introduction

Filaments are regions of dense cool plasma floating in the corona that are supported by magnetic fields. They are observed in absorption as dark features on the solar disk in H\(\alpha\) (6562.8 Å) and in emission as prominences above the solar limb. At the end of its life, a filament disappears by slow fading or exhibits a transient eruption. Before it disappears or erupts, small-scale blobs observed in H\(\alpha\) in a filament often show a larger standard deviation of the line-of-sight (LOS) velocity (Seki et al. 2017, 2019b). During eruption phase, the entire body of a filament ascends at a velocity of 100–1000 km \(\text {s}^{-1}\) (Parenti 2014).

Filament eruptions are often associated with coronal mass ejections (CMEs), which are observed by coronagraphs such as the Large Angle and Spectrometric Coronagraph (LASCO) (Brueckner et al. 1995). Some CMEs exhibit a three-part structure Illing and Hundhausen (1985) consisting of a leading edge, faint coronal cavity, and dense core. Others exhibit more complex forms, appearing as narrow jets or global eruptions, which are called halo CMEs (Webb 2015). The core of a CME is believed to originate from the filament mass if the CME is associated with a filament eruption. Because the Sun is concealed by an occulting disk in coronagraph observations of CMEs, it is difficult to observe the early evolution of CMEs. Investigations of the evolution of filament eruptions and their association with CMEs are expected to clarify the early evolution of CMEs.

CMEs often produce severe geomagnetic storms, which expose the Earth to a potential risk of adverse socioeconomic impacts such as a widespread blackout (Boteler 2001). A CME associated with a polar crown filament eruption reportedly caused a severe geomagnetic disturbance (Dst \(\sim\) − 200 nT) 3 days after the eruption (McAllister et al. 1994, 1996). Therefore, to mitigate the socioeconomic impacts of geomagnetic disturbances, it is essential to predict the occurrence of CMEs and their arrival to the Earth. And to reveal the relationship between CMEs and the eruption or disappearance of filaments is important for the better prediction.

However, the reported proportion of filament disappearances or eruptions that are associated with CMEs ranges from \(\sim\) 10 to \(\sim\) 90%. Hori and Culhane (2002) studied 50 prominence eruptions observed at 17 GHz by the Nobeyama Radioheliograph (Nakajima et al. 1994) and found that 92% of them were associated with CMEs. Seki et al. (2019a) investigated 43 filament disappearances in H\(\alpha\) data observed by the Solar Dynamics Doppler Imager (SDDI) (Ichimoto et al. 2017) on the Solar Magnetic Activity Research Telescope (SMART) (UeNo et al. 2004) at Hida Observatory, Kyoto University, and found that 50% of them were associated with CMEs. McCauley et al. (2015) studied 904 filament and prominence eruptions observed in He(II) (304 Å) by the Atmospheric Imaging Assembly (AIA) (Lemen et al. 2011) and found that 73% of them were associated with CMEs. In contrast, Al-Omari et al. (2010) automatically classified 7332 filament/prominence eruptions reported by the National Centers for Environmental InformationFootnote 1 as events associated or not associated with CMEs and found that only 17% of them were associated with CMEs. [For a more detailed summary of previous studies on the filament–CME association, see Table 1 in Al-Omari et al. (2010)]. It is supposed that the discrepancy among these results could depend on how to select the events.

In this study, we aim to reveal what parameters control the association between the filament eruptions and CMEs. We investigate the relationships between physical parameters that characterise filament eruptions, i.e., the length, velocity during eruption, and direction of eruption, and the CME association. Several studies have shown that these parameters are well correlated with the CME association of filament eruptions (Seki et al. 2019a; Gilbert et al. 2000; Morimoto and Kurokawa 2003). Gilbert et al. (2000) studied 54 prominence eruptions observed above the limb in H\(\alpha\). They defined “eruptive prominences” as those in which all or part of the material escaped from the solar gravitational field and “active prominences” as those in which none of the material appeared to escape. They found that eruptive prominences clearly had a larger apparent velocity (the velocity projected on the plane of the sky) above 1.10 solar radii than active prominences did and that eruptive prominences were more strongly associated with CMEs (94%) than active ones (46%). Our previous study (Seki et al. 2019a) found that filament eruptions are more likely to be associated with CMEs if the filament length exceeds 150 Mm, the maximum radial velocity exceeds 140 km \(\text {s}^{-1}\), or their direction is inclined by less than 48° with respect to the solar normal. Thus, in the present study, we focus on these three parameters of filament eruptions and investigate how the association rate varies with respect to them. Note that, in contrast to our previous work that we investigated the tendency of CME association with respect to individual physical parameters, the present study aims to improve the predictability of CME association by combining those parameters. In “Data” section, we provide a description of the data we utilised. In the succeeding section, the results will be provided, followed by summary and discussion.

Data

We selected events from the SMART/SDDI Filament Disappearance CatalogueFootnote 2 (hereafter, the catalogue) (Seki et al. 2019a). The unique advantage of the SDDI is its wide wavelength coverage, which makes it possible to determine the LOS velocity of erupting filaments up to 400 km \(\text{s}^{-1}\). The catalogue lists 43 filament/prominence disappearances observed by SDDI from 2016 May 1 to 2019 June 18, in which filaments/prominences totally disappeared at the H\(\alpha\) line centre. We selected 28 of these events that had a credibility value of 2 or 3 for CME association in the catalogue (description of “credibility” will be provided later). That is, we used only events whose CME association or non-association is fairly clear. Note that some of the events were excluded from our analysis, even though their credibility was 2 or 3, because (1) terrestrial clouds covered the target filaments, and it was impossible to estimate their precise LOS velocities (Nos. 001, 007, 018, 021, 022, 035, and 043), or (2) the length of the target filament could not be measured due to the lack of observation (No. 029). Most of the selected events (26 of 28) are filament eruptions, and two events (on May 24 and June 20 in 2016) are prominence eruptions observed on the solar limb. Hereafter, we refer to these 28 events simply as filament eruptions.

The credibility value indicates how credible the CME association of an event is. We made a movie containing solar full-disk images of each event in the H\(\alpha\) line centre or in 304 Å and LASCO C2 running difference images. While watching each movie, we examined the directional and temporal association of each filament eruption with the CME and assigned a credibility value based on our judgement. The actual movies used for examination can be accessed at the catalogue webpage (click the credibility column). Our judgement is based on (1) \(\varDelta T\), which is the difference between the time when a CME was first observed in LASCO and the time of total disappearance of the filament in H\(\alpha\) centre (same as FD_end_time in the catalogue), and (2) \(\varDelta \phi\), which is the difference between the position angle of the filament and the central position angle of the CME. Table 1 describes how we determined the credibility 2 or 3 of the CME association on the basis of \(\varDelta T\) and \(\varDelta \phi\). The credibility 1 was labeled to events which were difficult to determine one-to-one correspondence. More specifically, we labeled the credibility 1 (1) if two filaments disappeared within 1 h, the difference between their \(\varDelta \phi\)’s was within 15°, and they were diagnosed as being associated with the same CME, or (2) if there were flares located within 30° from the central position angle of the CME and within a few hours prior to CME occurrence. For example, the events No. 002 and No. 003 in the catalogue were categorised as the credibility of 1 because two filaments disappeared within 1 h, and it was ambiguous which of these events was actually associated with one CME. Another example is No. 012, in which a CME could be attributable to a C-class flare occurred in an active region rather than to the filament eruption of the interest, and thus we concluded the credibility of this event as 1.

Table 1 Criteria for determining the credibility of a CME association on the basis of \(\varDelta T\) and \(\varDelta \phi\)

Table 2 shows the selected 28 events with their CME associations and physical parameters.

Table 2 Filament eruptions used in this study.

‘Date and time’ is the start time of a filament eruption and is defined as the first observation of a dark feature in H\(\alpha\) − 0.5 Å. It is equivalent to ‘FD_start_time’ in the catalogue. ‘CME’ indicates whether a filament eruption is associated with a CME. \(V_{\text{r}\_\max}\), \(V_{\text{r}\_\text{fin}}\), \(V_{\text{pos}}\), and \(V_{\max}\) are determined as follows. We manually tracked and measured the position and LOS velocity of a blob at the apex of the filament that was present until its total disappearance in H\(\alpha\). Then, we constructed its three-dimensional velocity as a function of time. \(V_{\text{r}\_\max}\) is the maximum radial (or ascending) velocity during the eruption, whilst \(V_{\max}\) is the maximum magnitude of three-dimensional velocity. \(V_{\text{pos}}\) is equal to \(\sqrt{V_{x}^2 + V_{x}^2}\), where \(V_{x}\) and \(V_{y}\) are the velocities of the filament in the east–west and south–north directions on the plane of the sky, at the time of \(V_{\text{r}\_\max}\), respectively. \(V_{\text{r}\_\text{fin}}\) is the radial velocity at the last observation of a filament in H\(\alpha\). L is the length of a filament measured at the same time as ‘date and time’. The projection effect is corrected according to the location of the filament on the solar disk. \(\varTheta\) is the inclination angle between the direction of the filament velocity at the time of \(V_{\text{r}\_\max}\) and the solar normal (see Figure B on the catalogue webpage). For further details of how these values were determined, see our previous paper (Seki et al. 2019a).

Results

Figure 1 displays the CME associations according to L (vertical axis) and \(V_{\text{r}\_\max}\) (top left), \(V_{\max}\) (top right), or \(V_{\text{pos}}\) (bottom left) on a logarithmic scale.

Fig. 1
figure 1

Plots of filament eruptions according to \(V_{\text{r}\_\max}\) (top left), \(V_{\max}\) (top right), or \(V_{\text{pos}}\) (bottom left) and filament length, L, on a common logarithmic scale. \(V_{0}\) and \(L_{0}\) correspond to the typical velocity (100 km \(\text {s}^{-1}\)) and typical length (100 Mm) of filaments, respectively. Open circles and crosses represent events with and without CMEs, respectively. The solid lines are described in the text

Here, the length and velocities are normalised by \(L_{0} = 100\) Mm and \(V_{0} =\) 100 km \(\text {s}^{-1}\), respectively. We can see the tendency that the longer and faster filaments are more likely to be associated with CMEs. The solid lines in the panels are drawn by the following relationships;

$$\left( \frac{V_{\text{r}\_\max}}{V_{0}}\right) \times \left( \frac{L}{L_{0}}\right) ^{0.96}= 0.80,$$
(1)
$$\left( \frac{V_{\max}}{V_{0}}\right) \times \left( \frac{L}{L_{0}}\right) ^{1.1}= 1.4,$$
(2)
$$\left( \frac{V_{\text{pos}}}{V_{0}}\right) \times \left( \frac{L}{L_{0}}\right) ^{0.72}= 0.85.$$
(3)

They were determined by using the algorithm of Linear Support Vector Classification implemented in LIBLINEAR (Fan et al. 2008) (for further explanation of the algorithm, see Appendix and Fan et al. 2008). In the top left panel (\(V_{\text{r}\_\max}\)), 27 events out of 28 (96%) were correctly classified into the two groups of filament eruptions with (open circles) and without (crosses) CMEs, whilst in the other cases, 21 (for \(V_{\max}\)) and 25 (\(V_{\text{pos}}\)) events were correctly separated. If we make a separation so that the number of correctly classified events can be maximised, 27 (\(V_{\text{r}\_\max}\)), 24 (\(V_{\max}\)), and 26 (\(V_{\text{pos}}\)) events will be correctly classified (the separations not shown in the figure). Thus, a better prediction of the CME association could be obtained by using \(V_{\text{r}\_\max}\) rather than using \(V_{\max}\) or \(V_{\text{pos}}\) at least with our limited number of the events, 28. This result suggests the advantage of measuring the radial velocity of filament eruptions. It also suggests that measuring both the velocity and the length of filaments should contribute to the better prediction of CME occurrence. The three-dimensional velocity observation provides a better capability for predicting the occurrence of CMEs, whilst the H\(\alpha\) imaging observations without Doppler measurements (\(V_{\text{pos}}\)) can still contribute to it.

Figure 2 displays histograms of events with (grey) and without (black) CMEs with respect to the left-hand sides (LHS) of Eqs. (1)–(3). These histograms also demonstrate that the CME association is better identified when \(V_{\text{r}\_\max}\) is used than when \(V_{\max}\) or \(V_{\text{pos}}\) is used. In the top left panel, we see a clear bimodal distribution, which is less clear in the other cases. To confirm the bimodality quantitatively, we introduced a statistic, D-value, which is defined as

$$D \equiv \frac{|\mu _1 - \mu _2|}{\sigma },$$
(4)

where \(\mu _1\) and \(\mu _2\) are the averages of the two normal distributions fitted to the events with and without CMEs, and \(\sigma\) is equivalent to

$$\sigma = \sqrt{\frac{\sigma _1^2 + \sigma _2^2}{2}},$$
(5)

where \(\sigma _1\) and \(\sigma _2\) are their standard deviations. The D-value represents the distance between the means of two normal distributions relative to their standard deviations. These distributions can be regarded as being separated if the D-value is larger than 2 (Ashman et al. 1994; Muratov and Gnedin 2010; Carlson et al. 2018). The means and standard deviations of two normal distributions (with and without CMEs) for each case of Eqs. (1)–(3) are summarised in Table 3 together with the D-values. We obtained D values of 2.3, 1.8, and 1.7 for \(V_{\text{r}\_\max}\), \(V_{\max}\), and \(V_{\text{pos}}\), respectively. Only \(V_{\text{r}\_\max}\) exhibits the D-value larger than 2. The better bimodality when \(V_{\text{r}\_\max}\) is used is confirmed quantitatively.

Fig. 2
figure 2

Histograms of the LHS of Eqs. (1)–(3) on a common logarithmic scale. The dark and light grey bars correspond to the events without and with CMEs, respectively. The bars are stacked. The size of a bin is 0.33

Table 3 Summary statistics of the fitted normal distributions and the corresponding D values

Figure 3 shows the CME association with respect to the radial velocity of the last observation (\(V_{\text{r}\_\text{fin}}\)) and L (left panel) or a common log of the LHS of Eq. (1) (right panel) for each filament eruption. Open circles and crosses denote events with and without CMEs, respectively. Most of the filament eruptions (80%) with negative \(V_{\text{r}\_\text{fin}}\) (grey area), i.e., events in which the filaments fall back to the Sun, were not associated with CMEs. In addition, 77% of the filament eruptions with positive \(V_{\text{r}\_\text{fin}}\) and L larger than 70 Mm were associated with CMEs. Note that the filaments with the smaller (larger) value of the LHS of Eq. (1) similarly tend to have smaller (larger) \(V_{\text{r}\_\text{fin}}\) (see the right panel).

Fig. 3
figure 3

Plot of filament eruptions according to \(V_{\text{r}\_\text{fin}}\) and L (left) or common log of the product of normalised \(V_{\text{r}\_\max}\) and normalised L to the power of 0.92 (right). Open circles and crosses have the same meaning as in Fig. 1. Grey area corresponds to negative \(V_{\text{r}\_\text{fin}}\)

From Fig. 3, we can also recognise that there are exceptional events that were associated with CMEs despite their negative \(V_{\text{r}\_\text{fin}}\)’s (− 6.40 km \(\text {s}^{-1}\) and − 10.6 km \(\text {s}^{-1}\)). We speculate that the blobs which escaped the solar gravity and erupted into the interplanetary space became invisible in H\(\alpha\), and we tracked a part of the filament that fell back to the solar surface.

Figure 4 shows the CME association according to the LHS of Eq. (1) (vertical axis) and \(\varTheta\), the inclination angle (angle from the solar normal) of the velocities (horizontal axis). We found that 82% of the filament eruptions with directions that were inclined by more than 45° from the solar normal were not associated with CMEs, and 71% of those with their \(\varTheta\)’s smaller than 45° were associated with CMEs. Thus, the inclination angle of eruptions will provide a clue for forecasting the CME occurrence. Note that the LHS of Eq. (1) of the filament eruptions with their \(\varTheta\)’s larger than 45° seldom exceeds − 0.097 (= log10(0.80)). Figure 4 also shows that 86% of the events associated with CMEs have their \(\varTheta\)’s smaller than 45°, while only 36% of the non-associated ones have \(\varTheta\) smaller than 45°. These results are consistent with the work of Gopalswamy et al. (2003), in which they defined two types of prominence eruptions, radial and transverse events, according to the eruptive motion of the prominence observed on the solar limb. In their statistical study, they found that 94% of the prominence eruptions associated with CMEs were radial events (86% in our study) and that 76% of the transverse events were not associated with CMEs (82% in our study).

Fig. 4
figure 4

Plot of filament eruptions according to \(\varTheta\) (inclination angle in steradians and degrees) and a common log of the LHS of Eq. (1). The symbols have the same meaning as in Fig. 1

Figure 5 shows the linear speed of the CMEs in the SOHO/LASCO CME Catalog (Yashiro et al. 2004) against \(V_{\text{r}\_\max}\). The CME speed is expected to be larger than the velocity of the associated filament eruption (Gopalswamy et al. 2003). The black line indicates equal speeds. All the data points but one (\(V_{\text{r}\_\max}\) \(=\) 359 km \(\text {s}^{-1}\)) are above or close to the black line, as expected. However, since there is a large scatter in the ratio of CME linear speed to \(V_{\text{r}\_\max}\), it is difficult to predict CME linear speeds from the maximum radial velocity.

Fig. 5
figure 5

Plot of filament eruptions according to \(V_{\text{r}\_\max}\) and CME linear speed. The black line corresponds to equal speeds

Summary and discussion

In this study, we investigated the relationships between the physical parameters of filament eruptions (three-dimensional velocity, filament length, and direction of eruption) and their CME associations using 28 events observed by SDDI at Hida Observatory. We found that the filament eruptions are well separated into two groups of events, one with and the other without CMEs, according to the product of the normalised maximum ascending velocity (\(V_{\text{r}\_\max} / V_{0}\)) and the normalised filament length (\(L / L_{0}\)) to the power of 0.96, and that among the filament eruptions with \(\left( \frac{V_{\text{r}\_\max}}{V_{0}}\right) \times \left( \frac{L}{L_{0}}\right) ^{0.96} > 0.80\), 93% are associated with CMEs, and 100% of filament eruptions with the product \(< 0.80\) are not associated with CMEs. The apparent velocity and the length of filaments measured in H\(\alpha\) observation could also provide a good measure for predicting the occurrence of CMEs, though the accuracy of the prediction using the apparent velocity is worse than that using the radial velocity. Our results suggest that the three-dimensional velocity, or more specifically the radial velocity derived from it, and the length of the erupting filament are the notable parameters for improving the predictability of CME association. And thus, we suggests the importance of observations of the three-dimensional velocity of filament eruptions for the prediction of CMEs. It should be noted, however, that improvement of statistics, i.e., studies with a larger number of examples, are strongly required to confirm these results.

Here, we propose a possible physical interpretation for the solid line in the top left panel of Fig. 1. This line, which is represented by Eq. (1), successfully separates events into those with and without CMEs. We assume that (1) the cross section of filaments, A, follows the relationship of

$$\left( \frac{A}{A_{0}}\right) = \left( \frac{L}{L_{0}}\right) ,$$
(6)

where \(A_{0}\) is the typical cross section of filaments (100 \(\text {Mm}^{2}\)), and that (2) the average hydrogen density is common among filaments, i.e., 10\(^{11}\) \(\text {cm}^{-3}\), which is a typical value for quiescent prominences (Heinzel et al. 2008). Then, if we regard Eq. (1) as

$$\left( \frac{V_{\text{r}\_\max}}{V_{0}}\right) \times \left( \frac{L}{L_{0}}\right) \sim 0.80,$$
(7)

or \(V_{\text{r}\_\max} \times L \sim 8.0 \times 10^{6}\ \text {km}^{2}\ \text {s}^{-1}\), then its square represents the kinetic energy of an erupting filament, i.e., \(\frac{1}{2}\times\) proton mass \(\times\) density \(\times\) volume \(\times\) \({V_{\text{r}\_\max}}^{2}\) = 5.4 \(\times\) 10\(^{28}\) erg. This relationship could be regarded as the kinetic-energy threshold above which filament eruptions are associated with CMEs. Note that if the length of a filament is 100 Mm, the deduced mass gets 1.7 \(\times\) 10\(^{15}\) g. (Gilbert et al. 2006) reported the masses of 18 prominences, which ranged from (\(1.08 \pm 0.52\))\(\times 10^{14}\) to (2.09 ± 0.80)\(\times 10^{15}\) g. Our assumed “typical” mass is consistent with the reported values.

As mentioned in “Introduction” section, the CME association rates of filament eruptions reported to date range from \(\sim\) 10 to \(\sim\) 90%. Here, we provide a possible interpretation of this wide range based on our results. We showed that the product of the normalised radial velocity of eruptions and the normalised filament length makes a key contribution to the CME association. The high association rates of 80–90% in the past studies might be attributable to the criteria they used to select the events, under which the prominences have a predominantly large radial velocity and a large size. Gilbert et al. (2000) reported that 94% of eruptive prominences (for the definition, see “Introduction” section) were associated with CMEs. Gopalswamy et al. (2003) also reported that 83% of radial prominence eruptions were associated with CMEs. Their selected prominence eruptions should have had a predominant radial velocity. In addition, Gopalswamy et al. (2003) and Hori and Culhane (2002) detected prominences with the NoRH that has its spatial resolution of 10 arcsec (Nakajima et al. 1994), which is worse than the spatial sampling of the SDDI (1.23 arcsec pixel\(^{-1}\)). Therefore, the selected prominences in these studies seem to have a larger size (e.g., larger than 70 Mm, because 75% of the filaments smaller than 70 Mm were not associated with CMEs according to our result).

The association rate could also depend on whether studies include disk events (filament disappearances) in the records. In contrast to the high association rates (80 to 90%) reported in the studies taking into account only limb events (prominence disappearances) (Gilbert et al. 2000; Gopalswamy et al. 2003; Hori and Culhane 2002), some studies (Pojoga and Huang 2003; Jing et al. 2004; Seki et al. 2019a) in which both disk and limb events were considered manifested the association rate of approximately 40–50%. Pojoga and Huang (2003) reported that 39% of filament and prominence eruptions observed in H\(\alpha\) were associated with CMEs. Jing et al. (2004) reported that 56% of filament eruptions were associated with CMEs by automatically detecting filament disappearances in H\(\alpha\). In our study, considering only credible events, we found that 50% of filament eruptions in H\(\alpha\) were associated with CMEs.

Additionally, the observational wavelengths at which filaments or prominences are detected could also affect the association rate. In H\(\alpha\), as mentioned in the previous paragraph, approximately 40 to 50% of disappearance events were associated with CMEs. By contrast, McCauley et al. (2015) used full-disk solar images in the 171, 193, and 304 Å AIA passbands and reported an association rate of 72%.

The low association rate (17%, Al-Omari et al. 2010) might be attributable to the fact that the authors include ejecta such as surges in addition to filament eruptions in their sample. Among their 7332 events, they introduced 15 “filament types”, including coronal rain, sprays, and surges. In our study, we did not refer to these ejecta as filaments, and we excluded them from our list. Thus, the definition of filaments in that study was different from ours. Moreover, most of their events (\(\sim\) 80%) were smaller than \(\sim\) 70 Mm [see Figure 8 in Al-Omari et al. (2010)]. According to our result (see Fig. 1 or 3), 75% of the eruptions of filaments with lengths smaller than 70 Mm were not associated with CMEs. Assuming that this relation holds for coronal rain, surges, and sprays, \(\sim\) 60% (80% \(\times\) 75%) of all their selected events may not be associated with CMEs in our criteria. Therefore, the low association rate in their study can also be attributed to the event selection criteria; i.e., a significant portion of their events is thought to be located below the threshold line proposed in this paper.

The results of this study can be used to develop a methodology to predict the occurrence of CMEs by measuring the three-dimensional velocities of filament eruptions. Moreover, our previous works suggest that the occurrence of filament eruptions can be predicted prior to their initiation by \(1.3 \pm 0.47\) hour for intermediate filaments, on the basis of the mean and standard deviation of the LOS velocity distribution in filaments (Seki et al. 2017, 2019b). Hence, by using SDDI data and measuring the LOS velocity of filaments, we could predict the occurrence of filament eruptions \(\sim\) 1 h in advance and also, during eruptions, estimate the possibility of CME association before coronagraph observations.