Open Access
Issue
A&A
Volume 690, October 2024
Article Number A60
Number of page(s) 12
Section The Sun and the Heliosphere
DOI https://doi.org/10.1051/0004-6361/202450705
Published online 27 September 2024

© The Authors 2024

Licence Creative CommonsOpen 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 energetic particle (SEP) events result from acceleration processes associated with both solar flares and coronal mass ejections (CMEs). SEPs propagate in interplanetary space mostly along interplanetary magnetic field (IMF) lines before they are observed by spacecraft located in the heliosphere (see Desai & Giacalone 2016; Reames 2021, and references there in). A two-class paradigm classifies the SEP events into impulsive or gradual events. The impulsive SEP events are assumed to be associated with solar flares and type III radio bursts. They are limited in duration, reach low peak intensities, and typically have narrow emission cones (see, e.g., Reames 2021). Gradual SEP events are most energetic and are assumed to be associated with CMEs and type II radio bursts. They can last for several days, achieving elevated peak fluxes, and in general have a broad emission cone (see, e.g., Desai & Giacalone 2016). Nonetheless, this paradigm has proven to be a simplified view, and it is further complicated by the association of both strong flares and CMEs with SEPs (Cane et al. 2010; Papaioannou et al. 2016).

The fact that SEPs are driven by CMEs was first discussed by Kahler et al. (1978), who also indicated the close relation of the CME speed and the peak proton flux of SEPs by highlighting the fact that fast CMEs are more likely to drive shocks that are capable of accelerating energetic particles (see also Kahler 1982, 2001). This correlation was widely investigated and verified since then (e.g., Gopalswamy et al. 2002; Cane et al. 2010; Richardson et al. 2014; Papaioannou et al. 2016; Paassilta et al. 2017; Kihara et al. 2020). The interpretation is that shocks driven by fast CMEs are more likely to efficiently accelerate particles since theoretically, the acceleration rate depends on the speed of the shock with respect to the upstream medium (Lee et al. 2012).

Routine CME observations by the Large Angle and Spectrometric Coronagraph (LASCO; Brueckner et al. 1995) on board the Solar and Heliospheric Observatory (SOHO) have been performed since 1997. These observations have revealed that there is no one-to-one correspondence between X-ray flares and CMEs because many more flares than CMEs are observed. However, most CMEs are associated with some level of X-ray flare emission. In particular, Yashiro et al. (2006) (their Fig. 1) demonstrated that the stronger the flare in terms of its peak flux in the 1–8 Å soft X-ray (SXR) wavelength band as measured routinely by the GOES satellites, the more likely it is to be associated with a CME. This results from that fact that solar eruptive events (i.e., flares and CMEs) do not occur in isolation, but in concert, as a consequence of corresponding changes in the coronal magnetic field. It is noteworthy that recently, Li et al. (2021) reported a critical study of the flare–CME relation, showing that it depends on the flare class and the size of the active region (AR) of the source. Moreover, the solar origin of SEPs can almost always (∼94%, Papaioannou et al. 2016) be associated with the occurrence of both SXR flares and CMEs. Some additional SEP events generated at (or even beyond) the west limb of the Sun (Cane et al. 2010) could only be associated with CMEs because flare observations are not possible when they occur on the far side of the Sun. As a result, a wealth of statistical studies indicate that SXR flare peak fluxes, near-Sun CME speeds, the achieved peak proton intensities, and the fluences of the resulting SEP events are related. From recent studies investigating these correlations, it was shown that the most prominent correlation is found between the SEP peak proton flux and the speed of the CME (i.e., cc = 0.57 for E > 10 MeV), with a tendency to decrease at higher particle energies (i.e., 0.40 for E > 100 MeV, see, e.g. Papaioannou et al. 2016).

Scaling relations of the peak proton flux of SEPs (IP) to the speed of the CME (VCME) have been proposed by several authors. In particular, Kahler (2001) showed that the relation of the peak proton flux of SEPs at 2 and 20 MeV depends on the VCME in the form I P V CME 4.36 $ I_{P} \propto V_{\mathrm{CME}}^{4.36} $ and I P V CME 4.83 $ I_{P} \propto V_{\mathrm{CME}}^{4.83} $, respectively. Investigating 130 SEP events at an integral energy of E > 10 MeV and 88 SEP events at an integral energy of E > 100 MeV associated with CMEs originating at western longitudes (i.e., W20–W87°), Belov (2017) found dependences described as I P V n 4.02 ± 0.39 $ I_{P} \propto V_{n}^{4.02\pm0.39} $ and I P V n 3.01 ± 0.50 $ I_{P} \propto V_{n}^{3.01\pm0.50} $, respectively, with V n = V CME 1000 $ V_{n} = \frac{V_{\mathrm{CME}}}{1000} $1. Lario & Karelitz (2014) showed that the SEP peak intensity versus the CME speed in a good approximation follows a triangular distribution. This method directly provided upper limits to the peak proton particle intensity that can be observed in the prompt component of the SEP events. In particular, these authors showed that for three energy ranges spanning from 9–80 MeV, the resulting upper limit dependence scales with I P V CME γ $ I_{P} \propto V_{\mathrm{CME}}^{\gamma} $ with γ ranging from 4.90–5.63 (see their Fig. 5). Takahashi et al. (2016) showed on theoretical grounds that the upper limit for the peak proton flux of E > 10 MeV is proportional to the CME speed as I P V CME 5 $ I_{P} \propto ~V_{\mathrm{CME}}^{5} $.

In our previous study (Papaioannou et al. 2023, – hereafter part I), we presented the dependence of SEPs on flare parameters. In particular, we investigated the relation between the GOES 1–8 Å SXR peak flux (FSXR) of the parent flare versus the SEP event peak proton flux (IP) and fluence (FP). We showed that a direct estimation of the upper limit SEP fluence spectra based on FSXR alone is possible and leads to a quantification of the radiation environment. The present follow-up study deals with the dependence of SEPs on the properties of the associated CMEs. With this purpose, we use the catalog presented in detail in part I (Appendix C) and start with the analysis of the relations between VCME and IP, and we consider the dependence of VCME on FSXR. Based on the findings of Takahashi et al. (2016), showing that the upper limit for the peak proton flux of E > 10 MeV depends on the CME speed as I P V CME 5 $ I_{P} \propto V_{\mathrm{CME}}^{5} $, we derive upper limits and scaling relations among the CME speeds (VCME) and the achieved SEP peak flux (IP) at each integral energy (from E > 10 to E > 100 MeV). We then extend these relations to incorporate the fluence (FP) of SEPs. Additionally, we deduce the upper-limit fluence spectra of SEPs based on VCME, whereas in part I, the SEP fluence spectra were obtained based on FSXR.

In an attempt to estimate upper limits of the extreme events that can stem from our host star, Gopalswamy et al. (2010) presented calculations based on a hypothetical AR with the largest reported AR area (i.e., 5000 millionths of a solar hemisphere, msh) and the maximum measured sunspot magnetic field (B = 6100 G). These authors then estimated the potential energy of the AR to be 1036 erg, which could produce an SXR flare of ∼X1000 class (i.e., 10−1 W/m2). Consequently, the estimated maximum CME speed associated with the highest solar flare class, taking into account an upper limit of 26% of the potential energy being converted into CME kinetic energy, was found to be 7200 km s−1. This exceeds the highest measured CME speed by the SOHO/LASCO by a factor of ∼2 (Gopalswamy 2011). By considering the observations used in this work (Sect. 2), the obtained scaling relations (Sect. 3), and the upper limits our Sun can produce (Sect. 4), we estimate the fastest expected CMEs, the worst-case SEP proton fluxes and fluences, and the corresponding SEP spectrum. Implications for the effects of extreme CMEs on the radiation environment and the limits of recent fluence reconstructions on VCME are described and discussed.

2. Data sets

We focus on the relations between SEPs and CMEs. We used a well-defined catalog of 65 well-connected (W20-90°) SEP events that were recorded between 1984 and 2017 and extended from E > 10 to E > 100 MeV. For each event, we first identified the prompt peak intensity (in units of protons cm2 sr1 s1), defined as the maximum intensity observed shortly after the onset of the event in situ. In this way, energetic storm particles (or the ESP component) were excluded. Furthermore, we calculated the omnidirectional fluence (cm2) (integration over time) and tabulated all results (see Appendix C in part I). Because SOHO/LASCO measured linear CME speeds only starting in 1997, the sample used in this work was reduced from 65 to 38 SEP events. The CME speeds and widths were taken from the online CDAW CME catalog2 (Yashiro et al. 2004).

3. Scaling relations

3.1. Soft X-ray flare flux, coronal mass ejection speed, and peak proton fluxes

As a first step, we studied the scaling relations between the CME speed (VCME) and the prompt peak proton flux (IP) for integral energies E > 10 MeV, E > 30 MeV, E > 60 MeV, and E > 100 MeV. Figure 1 shows IP for E > 10 MeV as a function of VCME. The correlation coefficient between IP and VCME for the SEP events that reach an integral energy of E > 10 MeV is cc = 0.58, assuming a linear regression obtained from the reduced major axis (RMA, see discussion in Papaioannou et al. 2023) method, leading to I P V CME γ $ I_{P} \propto V_{\mathrm{CME}}^{\gamma} $, with γ = 5.09 ± 0.78. Additionally, the gray shaded envelope in Fig. 1 provides the estimated error3.

thumbnail Fig. 1.

IP ∝ VCMEγ relation for E > 10 MeV, with γ = 5.09 ± 0.78. The dashed black line corresponds to the upper limit of IP in terms of VCME with γ = 5. The upper point in our sample corresponds to the SEP event on 8 November 2000 (see the relevant appendix in part I) and is used for the scaling. It is labeled with a red dot. In addition, the event on 15 April 2001 is highlighted with a magenta dot (see details in the text).

As discussed in Takahashi et al. (2016), an upper limit for this relation is given by I P V CME 5 $ I_{P} \propto V_{\mathrm{CME}}^{5} $, when passing through the uppermost point of the employed sample. This scaling relation is in general deduced by employing three assumptions: (a) that the CME mass (MCME) is equal to the sum of the gravitational stratified AR corona (see Eq. (1) of Takahashi et al. 2016), (b) that the CME kinetic energy (ECME) is proportional to the total energy released during a flare, which is also a fraction f of the AR magnetic field energy (see e.g. Emslie et al. 2012; Papaioannou et al. 2023) (see Eq. (2) of Takahashi et al. 2016), and (c) that the total kinetic energy of SEPs (Ep) is proportional to the flare energy, and the duration of the proton flux enhancement is determined by the CME propagation timescale tCME ∝ L/VCME (L is the length scale of the flaring AR), which leads to a scaling relation of I P V CME 5 $ I_{P} \propto V_{\mathrm{CME}}^{5} $. Based on our sample, we derived that this upper limit lies at I P = 10 12.05 · V CME 5 $ I_{P} = 10^{-12.05} \cdot V_{\mathrm{CME}}^{5} $ (dashed black line in Fig. 1). The event that marks the upper limit fit is the event on 8 November 2000 SEP (indicated by a red dot in Fig. 1). We further highlighted the ground-level enhancement event of 15 April 2001 (GLE60) in magenta (see details in Appendix B).

A scatter plot between VCME and the flare SXR peak flux (FSXR) is presented in Fig. 2. The RMA regression is quantified as VCME F SXR 0.43 ± 0.08 $ F_{\mathrm{SXR}}^{0.43\pm 0.08} $ with a correlation coefficient of cc = 0.43. Takahashi et al. (2016) suggested that the upper limit of this relation is given by V CME F SXR 1 / 6 $ V_{\mathrm{CME}}\propto F_{\mathrm{SXR}}^{1/6} $. Following the description by Takahashi et al. (2016), based on our sample, we find the relation to be V CME = 1.3 × 10 4 · F SXR 1 / 6 $ V_{\mathrm{CME}} = 1.3 \times 10^{4} \cdot F_{\mathrm{SXR}}^{1/6} $ (dashed black line in Fig. 2), passing the point with the highest VCME = 3387 km/s corresponding to the event on 10 November 2004.

thumbnail Fig. 2.

Relation between CME velocity (VCME) and the flare SXR peak flux (FSXR). The solid black line represents the linear RMA regression fit V CME F SXR α $ V_{\mathrm{CME}}\propto F_{\mathrm{SXR}}^{\alpha} $ with α = 0.43 ± 0.08. The dashed black line corresponds to the upper limit of VCME in terms of FSXR. See text for details.

The results presented to date assumed statistical relations among the SXR peak flux of flares (FSXR), the speed of the CME (VCME), and the peak proton flux (IP) at an integral energy of E > 10 MeV (e.g. Takahashi et al. 2016). Solar flares and CMEs are both drivers of SEP events (e.g. Cane et al. 2010), and a correlation between VCME, FSXR, and IP has therefore often been put forward (e.g. Papaioannou et al. 2016). Nonetheless, previous findings have solely been derived from an E > 10 MeV sample of SEPs. Although it was neglected in most studies, we further investigate solar scaling relations here for the integrated E > 30 Mev, E > 60 Mev, and E > 100 MeV energy channels based on the same method. As a result, the scaling relations I P V CME γ $ I_{P} \propto V_{\mathrm{CME}}^{\gamma} $ are obtained for each of these integral energies (see Appendix A). Our results are presented in Fig. A.1, where we display the IPVCME relations for E > 30 MeV (top panel), E > 60 MeV (middle panel), and E > 100 MeV (bottom panel). Table 1 summarizes the slopes obtained with the RMA regression fits for each case.

Table 1.

Slopes of the relations obtained for the peak proton flux (IP) to the CME speed (VCME), and the corresponding correlation coefficients for each integral energy.

The power-law index γ presents a relatively slight increase with energy and varies between 5.09 ± 0.78 at E > 10 MeV up to 5.52 ± 1.02 at E > 100 MeV, but also shows a consequent slight increase in the uncertainties. At the same time, the correlation coefficients of the IPVCME relation seem to decrease with energy.

3.2. Establishing the relations

According to Takahashi et al. (2016), who followed the argumentation by Emslie et al. (2012) that the CME mass is the sum of the mass within gravitationally stratified AR corona, the kinetic energy of CMEs is proportional to the flare energy (ECME ∝ Eflare) and the energetic proton flux in response to the SXR class of flares can be estimated under the assumption that FSXR is roughly proportional to the total energy released during flares, that is, FSXR ∝ Eflare (consistent with the observational findings in Emslie et al. 2012). VCME then scales with FSXR as

V CME , upper = V 0 · F SXR 1 / 6 ( k m / s ) , $$ \begin{aligned} V_{\rm CME, upper} = V_{0} \cdot F_{\rm SXR}^{1/6} \mathrm {(km/s)}, \end{aligned} $$(1)

with FSXR normalized in units of 1 W/m2 and V0 = 1.3 ⋅ 104 km/s, as derived from our sample. Equation (1) is the dashed black line in Fig. 2.

Moreover, as outlined in Takahashi et al. (2016) and above, when we assume that the total kinetic energy of protons in SEP events (EP) is proportional to the flare energy (Eflare) and that the duration of the proton flux enhancement is determined by the CME propagation timescale tCME, it follows that EP ∝ IP ⋅ tCME ∝ Eflare. Therefore, IP is scaled with VCME as

I P , upper V CME 5 . $$ \begin{aligned} I_{P,\mathrm {upper}} \propto V_{\rm CME}^{5}. \end{aligned} $$(2)

The relations for each integral energy are positioned to run through the strongest SEP events of our sample, so that we can discuss the upper limits of the CME velocity (VCME, upper) and peak proton flux (IP,upper) in each plot of Figs. 1 and A.1.

4. Estimating the coronal mass ejection speed, peak proton flux, and fluence for extreme solar events

4.1. Upper limits

4.1.1. Estimating VCME, upper

Figure 3 is similar to Fig. 2 and shows a scatter plot between VCME and FSXR. The RMA fit (solid black line) is embedded into the gray error band and the upper limit is deduced from Eq. (1) and presented as the dashed black line. We investigated the upper limit of the VCME based on the SXR peak flux of the associated solar flare together with Eq. (1). Considering the FSXR values estimated by Cliver et al. (2022) of X400 ± 200 (i.e., 4 × 10−2 W/m2 ± 2 × 10−2 W/m2) for the AD774/775 SEP event and X28 ± 14 (i.e., 2.8 × 10−3 W/m2 ± 1.4 × 10−3 W/m2) for GLE05, we obtained V CME , upper = 7600 + 534 827 $ V_{\mathrm{CME, upper}} =\,{\sim} 7600 \substack{+534 \\ -827} $ km/s and V CME , upper = 4880 + 342 532 $ V_{\mathrm{CME, upper}} =\,{\sim} 4880\substack{+342 \\ -532} $ km/s for GLE05. (for more details about the FSXR ranges, see Cliver et al. 2022 and Papaioannou et al. 2023).

thumbnail Fig. 3.

Similar to Fig. 2 with the addition of Eq. (5) (blue line) and Eq. (6) (red line).

In the next step, known published scaling relations of VCME were investigated and compared with our findings. Gopalswamy (2018) showed that the highest expected CME speed might be as high as ∼7200 km/s. This estimate was based on extreme solar conditions, which differ from those of the current Sun. Interestingly, Fig. 7(a) of Gopalswamy (2018) provides an estimation of VCME and an upper limit (VCME, upper) based on the magnetic potential energy (MPE [erg]) of active regions (ARs), giving the corresponding empirical relations as

V CME = 748 · l o g ( M P E ) + 636 ( k m / s ) , $$ \begin{aligned} V_{\rm CME} = 748 \cdot \mathrm {log(MPE)} + 636 \mathrm {(km/s)}, \end{aligned} $$(3)

which represents the whole sample of CMEs and ARs considered by Gopalswamy (2018) and

V CME , upper = 1136 · l o g ( M P E ) + 1557 ( k m / s ) , $$ \begin{aligned} V_{\rm CME, upper} = 1136 \cdot \mathrm {log(MPE)} + 1557 \mathrm {(km/s)}, \end{aligned} $$(4)

which is extracted from the highest values of the Goplaswamy sample, leading to the upper limit VCME, upper.

According to Emslie et al. (2012), the bolometric flare energy is related to MPE by FTSI = 0.025 ⋅ MPE, or vice versa, MPE = 40 ⋅ FTSI. Using our 38 SEP events (see Appendix C of part I), we first substituted FSXR in Eq. (1) of Cliver et al. (2020), which translated FSXR to FTSI for each flare. From the above relation, the MPE was estimated. After they were calculated for each case, the expected CME speeds and their upper limit in our sample were derived using Eqs. (3) and (4). The obtained data were subsequently used to obtain fits of VCME (VCME, upper) versus FSXR, which were then added to Fig. 2, and we obtained the representation in Fig. 3. Thereby, the obtained fits are given by

V CME = 6162 · F SXR 0.158 [ 6.2 × 10 3 · F SXR 1 / 6 ] ( k m / s ) , $$ \begin{aligned} V_{\rm CME} = 6162 \cdot F_{\rm SXR}^{0.158} \,\,\,\left[\sim 6.2 \times 10^{3} \cdot F_{\rm SXR}^{1/6}\right]\mathrm {(km/s)}, \end{aligned} $$(5)

with FSXR normalized in units of 1 W/m2 and 6.2 × 103 in km/s. This relation provides the VCME ∝ FSXR relation presented as a blue line in Fig. 3 and

V CME , upper = 8734 · F SXR 0.124 [ 8.7 × 10 3 · F SXR 1 / 8 ] ( k m / s ) , $$ \begin{aligned} V_{\rm CME, upper} = 8734 \cdot F_{\rm SXR}^{0.124} \,\,\,\left[{\sim } 8.7 \times 10^{3} \cdot F_{\rm SXR}^{1/8}\right]\mathrm {(km/s)}, \end{aligned} $$(6)

with FSXR normalized in units of 1 W/m2 and 8.7 × 103 in km/s. This relation provides the VCME, upper ∝ FSXR relation presented as a red line in Fig. 3.

Table 2 shows that VCME, upper for the extreme cases of AD774/775 and GLE05, derived from Eq. (1) (Fig. 3; dashed black line), leads to CME speeds that are higher by ∼1.13–1.32 times than those obtained by Eq. (6) (solid red line in Fig. 3). Nonetheless, Fig. 3 shows that Eq. (6) includes all 38 SEP events of our sample, except for the event on 10 November 2004, and it provides an alternative upper limit for the CME speed compared to that proposed by Takahashi et al. (2016) (i.e., the dashed black line in Fig. 3). Because of the difference of the obtained slopes (i.e., red versus dashed black line), the difference is clearly larger for stronger, that is, extreme flares. In other words, the stronger the flare in terms of FSXR, the larger the difference for the different estimates of VCME, upper. At the same time, Eq. (5) (solid blue line in Fig. 3) has a similar slope as Eq. (1), but with a lower scaling. Thus, the blue line underestimates the CME speed, especially for M- and X-class flares of our sample, but provides a statistically deduced representation of the expected VCME, in agreement with Fig. 7(a) of Gopalswamy (2018).

Table 2.

Upper limit CME speeds (VCME, upper, [km/s]) for the SEP event on AD774/775 and GLE05 derived in this work (i.e., from Eq. (1) and Eq. (6)), while Eq. (5) provides the VCME for each event.

4.1.2. Peak fluxes and fluences of solar energetic particles driven by VCME

We investigated the relation between IP, FP, and VCME further. In particular, the upper limit peak proton flux was calculated using the I P = V energy · V CME 5 $ I_{P} = V_{\mathrm{energy}} \cdot V_{\mathrm{CME}}^{5} $ relations (where Venergy are the coefficients used in the dashed black lines in Fig. 1 and Fig. A.1 for all integral SEP energies employed in this work). Similar to part I, from the established upper-limit relation of the peak proton flux, IP,upper, the corresponding upper limit of the fluence FP,upper can also be retrieved. In this case, as a function of VCME as

F P , upper = F P , energy · ( V energy · V CME 5 ) δ , $$ \begin{aligned} F_{P,\mathrm {upper}} = F_{P, \mathrm {energy}} \cdot (V_{\rm energy} \cdot V_{\rm CME}^{5})^{\delta } , \end{aligned} $$(7)

with FP,energy and δ directly taken from Table A.1 of Papaioannou et al. (2023)4 and Venergy being VE10 = 10−12.05, VE30 = 10−12.55, VE60 = 10−13.14, and VE100 = 10−13.65. Venergy for each integral energy we considered was scaled with (km/s)−5.

The corresponding outputs are presented in Appendix B. Figure B.1 shows the FP versus VCME relation obtained for our sample of 38 SEP events for which CME information was available, using the RMA regression fit (solid black lines) for each integral energy embedded in a gray error-envelope. The dashed black lines correspond to the upper limits from Eq. (7).

In order to investigate the effect of VCME, upper and VCME on the calculation of IP and FP, we used the outputs from Eq. (1), Eq. (6), and Eq. (5) (see Table 2). Using these values (including upper and lower limits), the peak proton flux was derived employing the I P V CME 5 $ I_{P} \propto V_{\mathrm{CME}}^{5} $ relations (Figs. 1 and A.1), and the fluence was calculated by substituting VCME (or VCME, upper) in Eq. (7).

The obtained results are provided for the upper-limit peak proton flux and fluence in Tab. B.1, showing that the upper limit IP, which was calculated with VCME, upper from Eq. (1) and Eq. (6), differs by a factor of ∼2.5. The corresponding FP,upper differs by a factor of ∼3.0. At the same time, the obtained IP and FP are lower than the upper limits by a factor of ∼40 when VCME from Eq. (5) was used as input.

4.2. VCME and its upper limit

The upper limit of VCME obtained from Eq. (1) leads to a maximum value of ∼7600 km/s for an X425 SXR flare associated with the AD774/775 SEP event. This CME speed agrees with the value of ∼7200 km/s reported by Gopalswamy (2011). However, this high CME speed requires extreme conditions (i.e., B = 6100 G over the entire sunspot region) that are difficult to reconcile with the current Sun. Moreover, the RMA fit for the same case (see the solid black line in Fig. A.1) leads to values of ∼12 000 km/s (for the same SXR flare of X425), which is even more difficult to present with the inherent limitations imposed by the energetics of our Sun, and thus is unrealistic. Moreover, it was shown that VCME from Eq. (6) provides maximum values for the SXR flare associated with the AD774/775 event of 5300–6100 km/s, but this decreases to ∼3300–3900 km/s when using Eq. (5) (see Table 2 and Fig. 3).

During the modern era of CME measurements, the Sun produced an SEP event on 4 November 2003, which was associated with an X28 (X43.25) solar flare and a CME with a linear speed of 2657 km/s6. However, to determine the CME speed of this event, only three points were used in the LASCO field of view. Additionally, the CME emerged under distorted conditions, which complicates the speed estimates. Moreover, Gopalswamy et al. (2005) suggested that the maximum VCME of our host star may not be much higher than ∼3000 km/s.

Figure 4 shows our results for the FP - VCME dependence for E > 30 MeV. In each panel, we added the fluence of the extreme and rare SEP events found in the cosmogenic radionuclide records (e.g., Miyake et al. 2012; Mekhaldi et al. 2015; Brehm et al. 2021; Mekhaldi et al. 2021): the AD993 (upper left), AD774/775 (upper right), 660 BCE (lower left), and 7176 BCE (lower right). For these selected events, we used the recent fluence reconstructions reported by Koldobskiy et al. (2023). Thus far, when we used the obtained fluence range for each of these events (Y-axis in each panel), the VCME range leading to this fluence can be directly obtained based on the upper limit scaling law (dashed black line). Similar figures were constructed for E > 60 MeV (Fig. C.1) and E > 100 MeV (Fig. C.2). When the upper limit is used, the range of the CME speed for these extreme SEPs falls within a mid-mean range of 3241–5255 km/s based on all integral energies (i.e., see details in Appendix C and Table C.1).

thumbnail Fig. 4.

FP-VCME relations for E > 30 MeV. The magenta ribbons on the Y-axis correspond to the FP range for AD993, AD774/775, 7176 BCE, and 660 BCE (in the clockwise direction) as published by Koldobskiy et al. (2023). The ribbons on the X-axis show the estimated VCME range for the events based on the RMA fit (solid black line) and the upper limit (worst-case scenario; dashed black lines).

4.3. Spectrum based on VCME

The top panel of Fig. 5 provides the obtained integral fluence spectra of the AD774/775 event for the three different estimates of the CME speed based on the FSXR (see Table 2). In particular, the upper limit of VCME was obtained by Eq. (1) (dashed black line in Fig. 3), Eq. (5) (solid blue line in Fig. 3), and Eq. (6) (solid red line in Fig. 3). Thereafter, the FP − VCME relations (i.e., RMA fits in Fig. B.1 and Eq. (7) for the upper limits) were employed to identify the expected fluence at the respective integral energies. The solid black, blue, and red lines were obtained via the inverse power-law fit of the estimated fluences per integral energy of the uppermost points included in Table B.1. These points were calculated using Eq. (7). For each of the three fits, the shaded area provides the 1σ error. Estimates of the energy-dependent fluences of the AD774/775 event, independently obtained by Usoskin et al. (2021) and Koldobskiy et al. (2023), are included as blue and black squares, respectively.

thumbnail Fig. 5.

Fluence spectra of SEPs obtained from the scaling relation FP-VCME for AD774/775. Each solid line is the inverse power-law fit to the obtained integral fluence values acquired when using VCME from Table 2 (top panel) and the range of CME speeds obtained from Table C.1 (bottom panel). The derived AD774/775 fluences by Usoskin et al. (2021) and Koldobskiy et al. (2023) are shown as filled blue and black squares, respectively.

The obtained integral fluence spectra displayed in the top panel of Fig. 5 are driven by the associated VCME. The reconstructions by Usoskin et al. (2021) and Koldobskiy et al. (2023) seem to fall above the estimated fluence spectra obtained from an X425±175 flare converted into VCME using Eq. (5) (i.e., solid blue line) and below the upper limit ( V CME 5 $ V_{\mathrm{CME}}^{5} $) spectra presented in this work (i.e., solid black and red lines) when converting the same flare magnitude using Eq. (1) and Eq. (6), respectively. Therefore, the obtained upper limit (i.e., the worst case) fluence spectra in these two later cases (solid black and red lines, respectively) seem to restrict the actual fluence-reconstructed values in the case of AD774/775. Nonetheless, for higher energies (i.e., E > 430 MeV and E > 1000 MeV), the difference between the actual estimates of the fluence (black squares) and the obtained fluence spectra (solid black line) exceeds one order of magnitude, which suggests that a CME speed of 7602 + 532 829 $ \substack{+532 \\ -829} $ km/s overestimates the fluence at these energies.

The bottom panel of Fig. 5 is similar to the top panel. However, the obtained fluence spectrum here assumes a VCME range of 3421 km/s (lower limit) to 5482 km/s (upper limit) as obtained in Appendix C (magenta band). It shows that the independently obtained fluence estimates of Usoskin et al. (2021) and Koldobskiy et al. (2023) fall well within the estimated fluence range for energies between E > 30 MeV and E > 1 GeV. When we assume that the AD774/775 event is one of the strongest events ever recorded (based on our current knowledge), the fluence spectra depicted in the bottom panel of Fig. 5 offer an upper limit (worst case) fluence, including the independently obtained fluence values, based on a CME speed of 5482 km/s (see Appendix C for more details). Thus, in turn, a CME speed of ≤5500 km/s could be restrictive for our current Sun.

5. Conclusions

We investigated the dependence of SEP events on VCME. In particular, the scaling relations that describe IP, VCME, and FP seem to be consistent with statistical relations obtained by observations. These were extended from E > 10 MeV up to E > 100 MeV. Based on our analysis, we derived the maximum expected CME speed associated with the largest estimated FSXR flare that was unleashed by the Sun (i.e., X600) to range between ∼3950 and ∼8134 km/s, depending on the underlying relation (see Table 2). However, a CME speed as high as ∼8200 km/s could only result from exceptional and most likely unrealistic conditions (see Gopalswamy 2018).

Moreover, the use of Eq. (5) showed that a CME speed of ∼3950 km/s agrees relatively well with CME observations of the extreme case of 4 November 2001 (VCME = 2660 km/s). This VCME provides an upper-limit fluence spectrum that is consistent with observations for E > 30 MeV (blue line in Fig. 5, top panel). However, at higher energies from E > 60 MeV to E > 430 MeV, this CME speed leads to an underestimation of the fluence. At the same time, the use of Eq. (6) showed that a CME speed of ∼5380 km/s provides an upper-limit spectrum comparable with higher energies (i.e., E > 100 MeV, E > 200 MeV & E > 430 MeV), but it overestimates the respective fluence. Nonetheless, at E > 30 MeV, an upper-limit estimation of the fluence for AD774/775 based on VCME ≤ 4000 km/s seems to represent the obtained fluence better (Usoskin et al. 2021).

The scaling relations presented in this study provide a direct estimate of the upper-limit peak flux (IP) and fluence (FP) based on the VCME5 relation driven by the associated CME speed. The correlation coefficients of FSXR and VCME are roughly similar in strength (see Papaioannou et al. 2023, and this work). However, FSXR possibly provides more robust estimates of IP and FP than the CME speed because (a) the CME speed measurement is more uncertain and may be strongly variable during the early phase close to Sun (Zhang et al. 2001; Vršnak 2008), with peak accelerations as high as 10 km/s2 (Bein et al. 2011; Veronig et al. 2018), whereas the flare peak SXR flux (FSXR) is well measured7, and (b) the conversion of the FSXR into the VCME (Fig. 3) leads to three CME speed estimates for the same flare, depending on which relation is considered. Based on this and the scaling relations presented in part I, (Papaioannou et al. 2023), FSXR scales almost linearly (γ = 5/6) with IP and FP, whereas VCME scales with a γ = 5. As a result, the errors in determining the CME speed are strongly enhanced pertaining to the dependent quantity upon calculations. This spread of values obtained for fluences driven by the same CME speed (see Fig. 4) could also be explained by the apparent lack of VCME values that extend to ≥4000 km/s, which causes a critical gap that is not covered by (any of) the sample(s) used (in any such study). Additionally, one of the theoretical assumptions is that the relation tCME ∝ L/VCME holds true, where L is the AR size. However, for completeness, we note that when a constant length L0 is assumed instead, the scaling relation changes to FP ∝ VCME7. Although plausible for the acceleration of particles whose length scales are independent of the AR size, this scaling law would in fact mean that even slower CMEs would cause extreme fluences (orders of magnitude higher) at the respective energies employed in this work. Based on the measurements of CMEs of the last 27 years, this seems highly unlikely.

The top panel of Fig. 5 shows that the highest expected fluence for a given VCME cannot (or can only slightly) overcome values demarcated by the red line. In addition, the bottom panel of Fig. 5 highlights that applying the scaling relation F P V CME 5 $ F_{P} \propto V_{\mathrm{CME}}^{5} $ and assuming CME speeds between 3420 km/s and 5480 km/s leads to FP values for one of the largest SEP event observed to date (i.e., AD774/775) that agree well with the values calculated by Usoskin et al. (2021) and Koldobskiy et al. (2023).

We underline that these relations do not necessarily assume a pure solar flare or CME acceleration of SEPs. As noted in the pioneering work by Emslie et al. (2012), there is an interplay between X-ray and SEP emission in complex solar events, including CME generation, so that the energy release is distributed between radiation and accelerated particles. Moreover, scaling relations are inherent of caveats and limitations, as outlined in part I (Papaioannou et al. 2023), but provide valuable content for the estimation of the worst-case radiation environment based on the VCME alone.

Our results may apply to other Sun-like stars. Nonetheless, for our host star, Li et al. (2021) estimated that no more than 50% of solar flares with a magnitude ∼X100 would generate CMEs caused by the strong magnetic confinement exerted on flares/eruptions from the largest ARs in terms of their overall magnetic flux. This finding clearly has consequences for the rate of SEP production. Nevertheless, as highlighted in part I and Sect. 1, scaling relations offer valuable context for a worst-case (upper limit) estimate of the radiation environment because these relations inherently assume flares to be associated with CMEs, and consequently, with SEP events. However, this one-to-one association scenario is unrealistic and not observed on the Sun because tens of thousands of flares result in only a few hundred SEP events (Papaioannou et al. 2016). This finding was recently corroborated by Kahler & Ling (2023), who concluded that scaling relations cannot be used directly without taking this imbalance into account.

We further underlined the difficulties in identifying the upper CME speed and how they impact the resulting IP and FP values. In addition, caution is needed in the case of stellar CMEs because of the observational gap in the stellar regime: While numerous stellar flares have been observed in Sun-like stars, stellar CMEs might be rare and cannot be directly observed so far (see details in Moschou et al. 2019; Leitzinger & Odert 2022). Nonetheless, recent studies demonstrated a different approach for identifying stellar CMEs on cool stars. This approach is based on the sudden dimmings in the extreme-ultraviolet and X-ray emission caused by the CME mass loss (Veronig et al. 2021; Loyd et al. 2022; Notsu et al. 2024), and it shows potential for future research efforts.


1

The normalization employed in that paper does not affect the proportionality.

3

See Papaioannou et al. (2023) for details of the error estimation.

5

See the recent recalibration in Hudson et al. (2024).

7

The main reason is the projection effects, which depend on the location of the source region on the Sun (Paouris et al. 2021).

Acknowledgments

The authors would like to thank the anonymous referee for a critical and constructive reading of the manuscript and for valuable comments that improved the contents of the paper. AP, KH, and DL acknowledge the International Space Science Institute and the supported International Team 441: High EneRgy sOlar partICle Events Analysis (HEROIC). KH acknowledges the support of the DFG priority program SPP 1992 “Exploring the Diversity of Extrasolar Planets (HE 8392/1-1)”. AP and KH also acknowledge the supported International Team 464: The Role Of Solar And Stellar Energetic Particles On (Exo)Planetary Habitability (ETERNAL). AP and DL acknowledge support from NASA/LWS project NNH19ZDA001N-LWS. DL also acknowledges support from project NNH17ZDH001N-LWS. AMV acknowledges the Austrian Science Fund (FWF): project no. 10.55776/I4555. The authors would like to thank Ed Cliver for fruitful discussions that shaped the paper at hand.

References

  1. Bein, B. M., Berkebile-Stoiser, S., Veronig, A. M., et al. 2011, ApJ, 738, 191 [NASA ADS] [CrossRef] [Google Scholar]
  2. Belov, A. V. 2017, Geomagn. Aeron., 57, 727 [NASA ADS] [CrossRef] [Google Scholar]
  3. Brehm, N., Christl, M., Adolphi, F., et al. 2021, Nat. Geosci., 14, 10 [NASA ADS] [CrossRef] [Google Scholar]
  4. Brueckner, G. E., Howard, R. A., Koomen, M. J., et al. 1995, Sol. Phys., 162, 357 [NASA ADS] [CrossRef] [Google Scholar]
  5. Cane, H. V., Richardson, I. G., & von Rosenvinge, T. T. 2010, J. Geophys. Res. (Space Phys.), 115, A08101 [Google Scholar]
  6. Cliver, E. W., Hayakawa, H., Love, J. J., & Neidig, D. F. 2020, ApJ, 903, 41 [Google Scholar]
  7. Cliver, E. W., Schrijver, C. J., Shibata, K., & Usoskin, I. G. 2022, Liv. Rev. Sol. Phys., 19, 2 [NASA ADS] [CrossRef] [Google Scholar]
  8. Desai, M., & Giacalone, J. 2016, Liv. Rev. Sol. Phys., 13, 1 [CrossRef] [Google Scholar]
  9. Emslie, A. G., Dennis, B. R., Shih, A. Y., et al. 2012, ApJ, 759, 71 [Google Scholar]
  10. Gopalswamy, N. 2011, Astron. Soc. India Conf. Ser., 2, 241 [Google Scholar]
  11. Gopalswamy, N. 2018, Extreme Events in Geospace (Elsevier), 37 [CrossRef] [Google Scholar]
  12. Gopalswamy, N., Yashiro, S., Michałek, G., et al. 2002, ApJ, 572, L103 [NASA ADS] [CrossRef] [Google Scholar]
  13. Gopalswamy, N., Yashiro, S., Liu, Y., et al. 2005, J. Geophys. Res.: Space Phys., 110, A09S15 [Google Scholar]
  14. Gopalswamy, N., Akiyama, S., Yashiro, S., & Mäkelä, P. 2010, in Magnetic Coupling between the Interior and Atmosphere of the Sun, eds. S. Hasan, & R. J. Rutten (Berlin, Heidelberg: Springer, Berlin Heidelberg), 289 [CrossRef] [Google Scholar]
  15. Hudson, H., Cliver, E., White, S., et al. 2024, Sol. Phys., 299, 39 [NASA ADS] [CrossRef] [Google Scholar]
  16. Kahler, S. W. 1982, J. Geophys. Res., 87, 3439 [NASA ADS] [CrossRef] [Google Scholar]
  17. Kahler, S. W. 2001, J. Geophys. Res., 106, 20947 [NASA ADS] [CrossRef] [Google Scholar]
  18. Kahler, S. W., & Ling, A. G. 2023, ApJ, 956, 24 [NASA ADS] [CrossRef] [Google Scholar]
  19. Kahler, S. W., Hildner, E., & Van Hollebeke, M. A. I. 1978, Sol. Phys., 57, 429 [NASA ADS] [CrossRef] [Google Scholar]
  20. Kihara, K., Huang, Y., Nishimura, N., et al. 2020, ApJ, 900, 75 [NASA ADS] [CrossRef] [Google Scholar]
  21. Koldobskiy, S., Mekhaldi, F., Kovaltsov, G., & Usoskin, I. 2023, J. Geophys. Res. (Space Phys.), 128, e2022JA031186 [NASA ADS] [CrossRef] [Google Scholar]
  22. Lario, D., & Karelitz, A. 2014, J. Geophys. Res. (Space Phys.), 119, 4185 [NASA ADS] [CrossRef] [Google Scholar]
  23. Lee, M. A., Mewaldt, R. A., & Giacalone, J. 2012, Space Sci. Rev., 173, 247 [CrossRef] [Google Scholar]
  24. Leitzinger, M., & Odert, P. 2022, Serb. Astron. J., 205, 1 [NASA ADS] [CrossRef] [Google Scholar]
  25. Li, T., Chen, A., Hou, Y., et al. 2021, ApJ, 917, L29 [NASA ADS] [CrossRef] [Google Scholar]
  26. Loyd, R. O. P., Mason, J. P., Jin, M., et al. 2022, ApJ, 936, 170 [NASA ADS] [CrossRef] [Google Scholar]
  27. Mekhaldi, F., Muscheler, R., Adolphi, F., et al. 2015, Nat. Commun., 6, 1 [CrossRef] [Google Scholar]
  28. Mekhaldi, F., Adolphi, F., Herbst, K., & Muscheler, R. 2021, J. Geophys. Res.: Space Phys., 126, e2021JA029351 [CrossRef] [Google Scholar]
  29. Miyake, F., Nagaya, K., Masuda, K., & Nakamura, T. 2012, Nature, 486, 240 [NASA ADS] [CrossRef] [Google Scholar]
  30. Moschou, S.-P., Drake, J. J., Cohen, O., et al. 2019, ApJ, 877, 105 [NASA ADS] [CrossRef] [Google Scholar]
  31. Notsu, Y., Kowalski, A. F., Maehara, H., et al. 2024, ApJ, 961, 189 [NASA ADS] [CrossRef] [Google Scholar]
  32. Paassilta, M., Raukunen, O., Vainio, R., et al. 2017, J. Space Weather Space Clim., 7, A14 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  33. Paouris, E., Vourlidas, A., Papaioannou, A., & Anastasiadis, A. 2021, Space Weather, 19, e2020SW002617 [CrossRef] [Google Scholar]
  34. Papaioannou, A., Sandberg, I., Anastasiadis, A., et al. 2016, J. Space Weather Space Clim., 6, A42 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  35. Papaioannou, A., Herbst, K., Ramm, T., et al. 2023, A&A, 671, A66 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  36. Reames, D. V. 2021, Solar Energetic Particles. A Modern Primer on Understanding Sources, Acceleration and Propagation (Cham: Springer International Publishing), 978 [Google Scholar]
  37. Richardson, I. G., von Rosenvinge, T. T., Cane, H. V., et al. 2014, Sol. Phys., 289, 3059 [NASA ADS] [CrossRef] [Google Scholar]
  38. Takahashi, T., Mizuno, Y., & Shibata, K. 2016, ApJ, 833, L8 [NASA ADS] [CrossRef] [Google Scholar]
  39. Usoskin, I., Koldobskiy, S., Kovaltsov, G., et al. 2021, in Proceedings of 37th International Cosmic Ray Conference - PoS(ICRC2021), 395, 1319 [Google Scholar]
  40. Veronig, A. M., Podladchikova, T., Dissauer, K., et al. 2018, ApJ, 868, 107 [NASA ADS] [CrossRef] [Google Scholar]
  41. Veronig, A. M., Odert, P., Leitzinger, M., et al. 2021, Nat. Astron., 5, 697 [NASA ADS] [CrossRef] [Google Scholar]
  42. Vršnak, B. 2008, Ann. Geophys., 26, 3089 [CrossRef] [Google Scholar]
  43. Yashiro, S., Gopalswamy, N., Michalek, G., et al. 2004, J. Geophys. Res. (Space Phys.), 109, A07105 [NASA ADS] [CrossRef] [Google Scholar]
  44. Yashiro, S., Akiyama, S., Gopalswamy, N., & Howard, R. A. 2006, ApJ, 650, L143 [NASA ADS] [CrossRef] [Google Scholar]
  45. Zhang, J., Dere, K. P., Howard, R. A., Kundu, M. R., & White, S. M. 2001, ApJ, 559, 452 [Google Scholar]

Appendix A: The IP versus VCME relations for different integral energies

Here, we present the same scatter plots as for Fig. 1 but for IP at energies E> 30 MeV, E> 60 MeV, and E> 100 MeV. The red and magenta dots indicate the events on 8 November 2000 and 15 April 2001, respectively. Our scaling is based on the 8 November 2000 event. This is the event that achieves the highest peak proton flux across all energies considered (i.e. E> 10 - E> 100 MeV). From the interplay of the associated VCME with the achieved peak proton flux, the 8 November 2000 SEP event is clearly distinguished at E> 10 MeV (Fig. 1) and E> 30 MeV (Fig. A.1, top panel). For E> 60 MeV the 15 April 2001 SEP event falls on the upper-limit scaling deduced by the 8 November 2000 SEP event (Fig. A.1, dashed black line in the middle panel). However, for E> 100 MeV, although the achieved peak proton flux on 15 April 2001 is lower than that of the 8 November 2000 SEP event (Fig. A.1, bottom panel; y-axis) this SEP event is associated to a CME with a speed of 1199 km/s - which is lower than the CME speed of 1738 km/s associated with the 8 November 2000 SEP event. As a result, the magenta point seems to differentiate from the upper-limit scaling (dashed black line) by a factor of ∼2.6. Nonetheless, for consistency and in line with part I (see Appendix B of that work) we keep the upper-limit scaling bound to the 8 November 2000 SEP event (red dot) across all energies and propagate the relative error imposed by the differentiation at E> 100 MeV to our calculations.

thumbnail Fig. A.1.

Similar to Fig. 1. From top to bottom, these panels present results for E> 30-; E> 60-, and E> 100 MeV, respectively. In all panels, the red dot corresponds to the 8 November 2000 outstanding large SEP event (see text for further details).

Appendix B: The FP versus VCME relations

Figure B.1 depicts the FP-VCME relation obtained for our sample of 38 SEP events for which CME information was available. The RMA regression fit is presented as a solid black line at each integral energy (i.e. panel) embedded in a gray error-envelope. The dashed black lines represent the upper-limits obtained by using Eq. (7). In each case, the obtained dashed black line yields an upper-limit to the observed fluences.

thumbnail Fig. B.1.

Relation between SEP fluence (FP) and CME speed (VCME) for the four integral energy bands of the SEP events, i.e., E> 10-; E> 30-; E> 60-; and E> 100 MeV, respectively. The log-log relations are obtained with RMA regression fitting. The estimated upper limits of FP in terms of VCME are depicted as dashed black lines in each panel.

Table B.1.

Upper limit peak proton fluxes (IP, [pfu]) and fluences (FP, [cm−2]) for the SEP event on AD774/775 and GLE05 derived in this work, for each integral proton energy.

Appendix C: FP-VCME relations for E > 60 MeV and E > 100 MeV

F E > 30 MeV = 10 5.46 · ( 10 12.55 · V CME 5 ) 1.12 = 2.53513 · 10 9 · V CME 5.6 V CME ( E > 30 MeV ) = 34.2768 · F E > 30 MeV 0.178571 F E > 60 MeV = 10 5.51 · ( 10 13.14 · V CME 5 ) 1.08 = 2.08353 · 10 9 · V CME 5.4 V CME ( E > 60 MeV ) = 40.5163 · F E > 60 MeV 0.185185 F E > 100 MeV = 10 5.54 · ( 10 13.65 · v 5 ) 1.03 = 3.0234 · 10 9 · v 5.15 V CME ( E > 100 MeV ) = 45.11 · F E > 100 MeV 0.194175 $$ \begin{aligned} F_{\mathrm{E>30~MeV} } =&10^{5.46} \cdot \left(10^{-12.55}\cdot V_{\mathrm{CME} }^5\right)^{1.12}\\ =&2.53513\cdot 10^{-9}\cdot V_{\mathrm{CME} }^{5.6}\nonumber \\ V_{\mathrm{CME~(E>30~MeV)} } =&34.2768 \cdot F_{\mathrm{E>30~MeV} }^{0.178571} \\&\nonumber \\ \hline \nonumber \\&\nonumber \\ F_{\mathrm{E>60~MeV} } =&10^{5.51} \cdot \left(10^{-13.14}\cdot V_{\mathrm{CME} }^5\right)^{1.08}\\ =&2.08353\cdot 10^{-9}\cdot V_{\mathrm{CME} }^{5.4}\nonumber \\ V_{\mathrm{CME~(E>60~MeV)} } =&40.5163 \cdot F_{\mathrm{E>60~MeV} }^{0.185185}\\&\nonumber \\ \hline \nonumber \\&\nonumber \\ F_{\mathrm{E>100~MeV} } =&10^{5.54} \cdot \left(10^{-13.65}\cdot v^5\right)^{1.03}\\ =&3.0234\cdot 10^{-9}\cdot v^{5.15}\nonumber \\ V_{\mathrm{CME~(E>100~MeV)} } =&45.11 \cdot F_{\mathrm{E>100~MeV} }^{0.194175} \end{aligned} $$(C.1)

thumbnail Fig. C.1.

Same as Fig. 4 but for E> 60 MeV.

thumbnail Fig. C.2.

Same as Figs. 4 and C.1 but for E> 100 MeV.

Table C.1.

The integral fluence FP of the extreme SEP events and the resulting VCME per integral energy.

All Tables

Table 1.

Slopes of the relations obtained for the peak proton flux (IP) to the CME speed (VCME), and the corresponding correlation coefficients for each integral energy.

Table 2.

Upper limit CME speeds (VCME, upper, [km/s]) for the SEP event on AD774/775 and GLE05 derived in this work (i.e., from Eq. (1) and Eq. (6)), while Eq. (5) provides the VCME for each event.

Table B.1.

Upper limit peak proton fluxes (IP, [pfu]) and fluences (FP, [cm−2]) for the SEP event on AD774/775 and GLE05 derived in this work, for each integral proton energy.

Table C.1.

The integral fluence FP of the extreme SEP events and the resulting VCME per integral energy.

All Figures

thumbnail Fig. 1.

IP ∝ VCMEγ relation for E > 10 MeV, with γ = 5.09 ± 0.78. The dashed black line corresponds to the upper limit of IP in terms of VCME with γ = 5. The upper point in our sample corresponds to the SEP event on 8 November 2000 (see the relevant appendix in part I) and is used for the scaling. It is labeled with a red dot. In addition, the event on 15 April 2001 is highlighted with a magenta dot (see details in the text).

In the text
thumbnail Fig. 2.

Relation between CME velocity (VCME) and the flare SXR peak flux (FSXR). The solid black line represents the linear RMA regression fit V CME F SXR α $ V_{\mathrm{CME}}\propto F_{\mathrm{SXR}}^{\alpha} $ with α = 0.43 ± 0.08. The dashed black line corresponds to the upper limit of VCME in terms of FSXR. See text for details.

In the text
thumbnail Fig. 3.

Similar to Fig. 2 with the addition of Eq. (5) (blue line) and Eq. (6) (red line).

In the text
thumbnail Fig. 4.

FP-VCME relations for E > 30 MeV. The magenta ribbons on the Y-axis correspond to the FP range for AD993, AD774/775, 7176 BCE, and 660 BCE (in the clockwise direction) as published by Koldobskiy et al. (2023). The ribbons on the X-axis show the estimated VCME range for the events based on the RMA fit (solid black line) and the upper limit (worst-case scenario; dashed black lines).

In the text
thumbnail Fig. 5.

Fluence spectra of SEPs obtained from the scaling relation FP-VCME for AD774/775. Each solid line is the inverse power-law fit to the obtained integral fluence values acquired when using VCME from Table 2 (top panel) and the range of CME speeds obtained from Table C.1 (bottom panel). The derived AD774/775 fluences by Usoskin et al. (2021) and Koldobskiy et al. (2023) are shown as filled blue and black squares, respectively.

In the text
thumbnail Fig. A.1.

Similar to Fig. 1. From top to bottom, these panels present results for E> 30-; E> 60-, and E> 100 MeV, respectively. In all panels, the red dot corresponds to the 8 November 2000 outstanding large SEP event (see text for further details).

In the text
thumbnail Fig. B.1.

Relation between SEP fluence (FP) and CME speed (VCME) for the four integral energy bands of the SEP events, i.e., E> 10-; E> 30-; E> 60-; and E> 100 MeV, respectively. The log-log relations are obtained with RMA regression fitting. The estimated upper limits of FP in terms of VCME are depicted as dashed black lines in each panel.

In the text
thumbnail Fig. C.1.

Same as Fig. 4 but for E> 60 MeV.

In the text
thumbnail Fig. C.2.

Same as Figs. 4 and C.1 but for E> 100 MeV.

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.