High-Resolution Spectral and Anisotropy Characteristics of Solar Protons During the GLE N∘73 on 28 October 2021 Derived with Neutron-Monitor Data Analysis

The first ground-level enhancement of the current Solar Cycle 25 occurred on 28 October 2021. It was observed by several space-borne and ground-based instruments, specifically neutron monitors. A moderate count-rate increase over the background was observed by high-altitude polar stations on the South Pole and Dome C stations at the Antarctic plateau. Most of the neutron monitors registered only marginal count-rate increases. Using detrended records and employing a method verified by direct space-borne measurements, we derive the rigidity spectra and angular distributions of the incoming solar protons in the vicinity of Earth. For the analysis, we employed a newly computed and parameterized neutron-monitor yield function. The rigidity spectra and anisotropy of solar protons were obtained in their time evolution throughout the event. A comparison with the Solar and Heliospheric Observatory/Energetic and Relativistic Nuclei and Electron (SOHO/ENRE) experiment data is also performed. We briefly discuss the results derived from our analysis.


Introduction
Eruptive processes on the Sun, e.g., flares and/or coronal mass ejections, result in solar energetic particle (SEP) events, that is, distinct enhancements of fluxes of protons, heavy ions and electrons (e.g. Aschwanden, 2012;Desai and Giacalone, 2016;Klein and Dalla, 2017, and references therein). The most energetic SEP events can last from hours to up to several days (Moraal and McCracken, 2012;Vainio et al., 2009Vainio et al., , 2013Gopalswamy et al., 2014;Raukunen et al., 2018;Usoskin et al., 2020). The energy of SEPs in most cases is in the MeV range, yet in some cases, accelerated protons can gain energy up to the GeV energy range.
A study of SEP properties provides a unique basis to reveal important processes of both SEP acceleration on the Sun and their propagation in the interplanetary space (e.g. Debrunner et al., 1988;Reames, 1999;Gopalswamy et al., 2014;Kocharov et al., 2021). A specific interest is paid to relatively rare events, with SEP energy reaching about GeV/nucleon or even greater values. Solar protons with energy above 300 MeV/nucleon can produce secondary particles in the Earth's atmosphere reaching the ground by generating an atmospheric cascade, whose byproducts are eventually registered by ground-based detectors, e.g., neutron monitors (NMs) (e.g. Dorman, 2004;Kühl et al., 2017;Mishev and Poluianov, 2021, and references therein). This class of events is called ground-level enhancements (GLEs) (for details see Shea and Smart, 1982;Poluianov et al., 2017). The current paradigm of ground level enhancements (GLEs) is that most likely they represent the high-energy tail of SEPs and are related to both solar flares and CMEs (e.g. Cliver, 2016;Desai and Giacalone, 2016;Miroshnichenko, 2018;Anastasiadis et al., 2019;Kocharov et al., 2020, and references therein). GLEs can be conveniently studied with ground-based and space-borne instruments (e.g. Bieber and Evenson, 1995;Simpson, 2000;Bruno et al., 2018), yet the latter usually orbit most of the time in high-rigidity cut-off regions, therefore being not full-time sensitive to the SEPs, whilst the former, specifically high-altitude polar instruments are sensitive to solar protons with energy of about 300 MeV/nucleon (e.g. Mishev and Poluianov, 2021).
On the other hand, GLEs can be registered and studied using the worldwide NM network (Simpson, Fonger, and Treiman, 1953;Hatton, 1971;Stoker, Dorman, and Clem, 2000;Mavromichalaki et al., 2011;Papaioannou et al., 2014), that is, stations located at different geographic regions are sensitive to a different part of the SEP spectra and arrival direction so that using the geomagnetosphere as a spectrometer, one can reveal their characteristics (e.g. Bieber and Evenson, 1995). At present, 73 GLEs have been registered, the first to the fifth by ionization chambers. Starting from GLE N • 5, their registration is carried out by NMs, and the records are stored in the International GLE Database (IGLED) (https://gle.oulu.fi, for details see Usoskin et al., 2020).
GLEs occur sporadically and differ from each other in the shape of the spectra, particle flux, angular distribution, duration, as well as time evolution of their characteristics (e.g. Moraal and McCracken, 2012;Raukunen et al., 2018;Koldobskiy et al., 2021). Therefore, GLEs are studied case-by-case, yet some common features have been recently reported (Kocharov et al., 2015. Here, we investigate the most recent GLE, namely, GLE N • 73 registered on 28 October 2021 (e.g. Papaioannou et al., 2022;Velinov, 2022), and present the rigidity spectra and angular distribution of SEPs, including their dynamical evolution throughout the event, employing NM data analysis.

GLE N • 73 on 28 October 2021
The first GLE event of the current Solar Cycle 25 was observed on 28 October 2021 as a weak increase (below 20% with respect to the galactic cosmic-ray (GCR) background) by several NMs, specifically those located in the low-rigidity cut-off region and several space-borne instruments (for details see Papaioannou et al., 2022). The 28 October 2021 event is associated with a gradual, class X1.0 flare located at S28W01, with peak of soft Xray emission at 15:35 UT, and an asymmetric halo coronal mass ejection (CME), brightest over the southern hemisphere of the Sun according to the Large Angle and Spectrometric Coronagraph (LASCO) C2 observations (SOHO/LASCO C2; Brueckner et al., 1995).
In Figure 1 we present the detrended NM count-rate increases registered by several selected stations during the event. We note that detrended data account for the GCR baseline temporal variability, specifically the contribution of short-time variations of GCRs, most likely due to transients and local anisotropy. Therefore, the detrended data provide smooth and free from transients, including diurnal variations, records of NM count-rate increases due to SEPs (for details see Usoskin et al., 2020).
The peak count-rate increase was registered by low-rigidity cut-off, high-altitude polar NMs, namely those located at the South Pole and French-Italian Dome C (Concordia) research stations, South Pole (SOPO) (5.4%) and South Pole Bare (SOPB) (5.7%), DOMC (7.3%) and DOMB (14%), standard and bare monitors, respectively (see Table 1). The highaltitude polar NMs have greater sensitivity in energy to SEPs compared to the sea-level ones because of the reduced atmospheric attenuation, namely about 300 MeV/nucleon for the former and about 430 MeV/nucleon for the latter (for details see Mishev and Poluianov, 2021, and the discussion therein). In addition, bare NMs are more sensitive to the low-energy part of the GLE-producing SEPs (e.g. Clem and Dorman, 2000;Vashenyuk, Balabin, and Stoker, 2007;Nuntiyakul et al., 2018Nuntiyakul et al., , 2020. Several NM stations observed the event onset at about 15:50 UT (e.g. SOPO and Fort Smith (FSMT)). The event lasted for about 4.5 h, however, the statistically significant signal registered by a sufficient number of NMs, allowing reliable data analysis (see the discussion in , was observed from about 16:00 UT to about 20:00 UT. The event revealed a slow, gradual increase of NM count rates and moderate anisotropy, that is the angular distribution width σ 2 ≈ π , similarly to the event(s) described by Bombardieri et al. (2006) as opposed to very anisotropic events described by Bütikofer et al. (2009). The event occurred on the background of a strong diurnal wave caused by the local anisotropy of GCRs, accordingly, the background was detrended and averaged over two hours before the event onset employing the procedure described in Usoskin et al. (2020). For the data analysis, we consider five-minute integrated detrended records of the NM count-rate increases available at the IGLED. Here, we emphasize that the one-minute records of the NM countrate increases depicted considerably larger fluctuations, not allowing strictly precise data analysis, yet some estimations are presented below.

Analysis of the Neutron-Monitor Data
In this section, we employ a method and model for an analysis of the NM data, based on the algorithm initially developed by Shea and Smart (1982), Cramp et al. (1997), , Bombardieri et al. (2007). The details and applications are given elsewhere (Mishev and Usoskin, 2016;Mishev, Poluianov, and Usoskin, 2017;Mishev et al., 2018Mishev et al., , 2021a. The data analysis of the NM records allowing us to derive SEP spectra, angular distribution, and apparent source position, is based on modeling the response of the global NM network and optimization of the modeled over the experimental data, assuming an appropriate initial guess (Cramp, Humble, and Duldig, 1995).
The method initially developed by our team (e.g. Mishev, Kocharov, and Usoskin, 2014) was recently improved and verified by direct space-borne measurements (for details see Koldobskiy et al., 2019a;Mishev et al., 2021b;Koldobskiy et al., 2021). In this article, we perform the modeling using a newly computed and verified, altitude-dependent NM yield function , which depicted reasonable latitude and altitude surveys, and space-borne records by PAMELA (Payload for Antimatter Matter Exploration and Lightnuclei Astrophysics, Adriani et al., 2017) and AMS-02 (Alpha Magnetic Spectrometer, Aguilar et al., 2021) (for details see the discussions in Lara, Borgazzi, and Caballero-Lopez, 2016;Nuntiyakul et al., 2018;Koldobskiy et al., 2019b). We benefit from the employment of stations located at different altitudes as well as the pairs of standard and bare NMs (South Pole and Dome C) (for details, see e.g., Ruffolo et al., 2006;Bieber et al., 2013). In addition, the optimization by Aleksandrov (1971), Golub and Van Loan (1980), Golub, Hansen, andO'Leary (1999), Mishev, Mavrodiev, andStamenov (2005) to the Levenberg-Marquardt (Levenberg, 1944;Marquardt, 1963) algorithm used in this work, allowed us to reduce unfolding uncertainties and obtain a reliable solution even in the case of ill-posed problem(s), arising from the weak and fluctuated NM responses (e.g., see the discussions in Tikhonov et al., 1995;Mavrodiev, Mishev, and Stamenov, 2004;Aster, Borchers, and Thurber, 2005;Mishev et al., 2021b).
The merit function (D in Equation 1) that is the main criterion for the quality of the fit, represents the residual (Himmelblau, 1972;Dennis and Schnabel, 1996): simultaneously including additional criteria as discussed in Mishev et al. (2021a,b). In Equation 1, N i , corresponds to the difference between the modeled and recorded count-rate increases of the ith NM, and N i to the nominal count-rate increase of the ith NM. Whilst D ≤ 5% for strong events (e.g., see , for weak events it is normally ≈ 10 -15%, and in some cases about 20% (e.g. Mishev et al., 2018). In our model, we can imply a spectral approximation with a modified power-law or an exponential rigidity spectrum as in Cramp et al. (1997), Vashenyuk et al. (2008), Mishev et al. (2021b.
The analytical expression of the rigidity spectrum of SEPs with rigidity P > 1 GV given by a modified power law is: where the flux of particles with rigidity P in [GV] is along the axis of symmetry identified by geographic latitude and longitude and the power-law exponent is γ with the steepening of δγ , J 0 is the particle flux at 1 GV in [m −2 s −1 sr −1 GV −1 ]. At P ≤ 1 GV, the rigidity spectrum of SEPs is approximated with: The exponential shape is described by the expression: where P 0 is a characteristic proton rigidity in [GV]. The pitch-angle distribution (PAD) was assumed to be a Gaussian-like distribution: where α is the pitch angle, σ accounts for the width of the distribution. The propagation of SEPs in the geomagnetosphere, necessary for the computation of the rigidity cut-offs and asymptotic directions of the NMs selected for the analysis (Cooke et al., 1991), is performed implying a superposition of the International Geomagnetic Reference Field (IGRF) geomagnetic model (epoch 2020) as the internal field model (Alken et al., 2021) and the Tsyganenko-89 model as the external field (Tsyganenko, 1989), which assures straightforward and reasonably accurate modeling of SEPs propagation in the Earth's magnetosphere (Kudela and Usoskin, 2004;Kudela, Bučik, and Bobik, 2008;Nevalainen, Usoskin, and Mishev, 2013). The NM stations used for our analysis with their standard acronyms, rigidity cut-offs (P c ), geographic coordinates, altitudes above the sea level, and the type of the instrument are given in Table 1.

Results of the Analysis
Before the unfolding of the NM records, we computed the rigidity cut-offs and asymptotic directions for all NMs from Table 1. An example of NM asymptotic directions in the rigidity range 1 -5 GV, corresponding to the rigidity cut-off of the station and encompassing maximal NM response (DOMC and SOPO are plotted in the range 0.7 -5 GV), is presented in Figure 2, while in the analysis we considered the 1 -20 GV rigidity range.
Subsequently, we examine all possibilities in our model as spectral and PAD functional shapes, that is, by modified power law or exponential rigidity spectra of SEPs (Equations 2-4), as well as by Ellison and Ramaty (1985) spectral form.
The best fit is obtained using a modified power-law rigidity spectrum of SEPs and single-Gaussian PAD, depicted for various stages of the event in Figure 3, the details are presented Table 1 Neutron monitors with standard acronyms, geographic coordinates (latitude and longitude), rigidity cut-off, altitude above the sea level, and type of instrument used for the analysis of GLE N • 73. in Table 2. The angular distribution of SEPs derived in this case is simpler than that of the two previous GLEs (Mishev, Kocharov, and Usoskin, 2014;Adriani et al., 2017;Mishev et al., 2021b). An analysis of 1-min records, available in the NM database (Mavromichalaki et al., 2011), was also performed during the event onset. We note that, in this case, considerably larger fluctuations in NM count rates resulted in greater residuals, therefore not all fit-quality criteria were fulfilled (for details see the discussion in Mishev et al., 2021a). However, we present the results of this analysis in Table 3, which are important in relation to particleacceleration studies (e.g. Klein and Trottet, 2001, and references therein).
The derived SEP spectra are moderately hard, slightly softer compared to the previous two weak events: GLE N • 71 and GLE N • 72 (Plainaki et al., 2014;Adriani et al., 2017;Mishev et al., 2018;Bruno et al., 2018Bruno et al., , 2019Mishev et al., 2021a). The SEP spectra revealed an important steepening, with considerable roll-off, significantly greater than that of GLE N • 71 and GLE N • 72 (Mishev et al., , 2021a. The spectra slowly and steadily softened throughout the event, accordingly, the steepening decreased and vanished in the Table 2 Derived spectral and angular characteristics during GLE N • 73 on 28 October 2021 fitted with a modified power-law rigidity spectrum. The columns correspond to the integration interval (1), particle flux (2), spectrum slope (3), steepening of the spectrum (4), width of the angular distribution (5), anisotropy axis position (6,7), merit function D (8), and normalized to degrees of freedom χ 2 r (9). late phase, which was after 19:00 UT. The PAD was wider compared to beam-like events such as GLE N • 69 and GLE N • 70 (for details see Bütikofer et al., 2009) and broadened out during the event. In addition, it was not complicated, in contrast to GLE N • 71 or GLE N • 72, where Sun-anti-Sun SEP flux was observed (e.g. Adriani et al., 2015;Mishev et al., 2021a). The intensity of the SEP flux gradually increased during the initial and main phase of the event (roughly corresponding to the prompt component (Vashenyuk et al., 2008)), reaching its peak at 18:15 UT (DOMB maximal count-rate increase of about 14%), and decreased afterwards.

Integration interval UT
We assess the quality of the fit by comparison between the modeled and the experimentally measured NM count-rate increases. Selected stations are presented in Figure 4. We note that the quality of the fit is similar for the other NM stations.
Similarly to Mishev et al. (2021b), we present in Figure 5 the contour plot of the sum of variances for the best-fit solutions vs. geographic coordinates. Here, the forward modeling is performed over all the possible apparent source positions. We emphasize that an erroneous apparent source-position determination would distort the derived PADs and spectra (see the discussion in Mishev et al., 2021a). One can see that the results of the forward modeling are satisfactory, i.e. the derived solutions are of reasonable quality, implying that the model well described the experimental NM records. Note that the derived apparent source position does not coincide exactly with the minimum of D because we employed additional fit-quality criteria as discussed in Mishev et al. (2021a,b).
The possibility to fit the SEP spectra with an exponential rigidity spectrum, assuming the same single-Gaussian PAD as derived above, was also examined, details are given in Table 4, and accordingly with the 1-min resolution in Table 5. In this case, the residual, i.e. D was greater compared to the previous case, specifically after the event onset, that is during the main and late phase, yet during the initial phase of the event, the goodness of fit was of the same order. Therefore, during the event initial phase, the SEP spectra could be fitted with either a modified power law with considerable roll-off or an exponential rigidity spectrum. We note that after 17:00 UT, the only possibility to describe the SEP spectra was with a modified power-law rigidity spectrum, since assuming an exponent rigidity spectrum led to D ≈ 40 -50%.

Figure 5
Contour plot of D for the best-fit solutions vs. geographic latitude and longitude during the initial phase of GLE N • 73, that was at 16:30 UT. The small white circle depicts the derived apparent source position, the cross depicts the IMF direction measured by ACE satellite.
The particle fluence (the time and angle integrated intensity of protons) of GLE N • 73 is presented in Figure 6. It is compared with the fluence reconstruction employing the "fast" method (details given in Koldobskiy, Kovaltsov, and Usoskin, 2018;Koldobskiy et al., 2019a;Usoskin et al., 2020). This method allows us to consider each NM as an integrating detector such that its SEP-related count rate is directly proportional to the SEP fluence at some rigidity P (assumed to be constant for a given NM). This method allows us to perform  Table 5 Derived spectral and angular characteristics using the exponential rigidity spectrum with 1-minute time resolution during the event onset of GLE N • 73. a robust reconstruction of the fluence during GLEs, when studying temporal and spatial features is impossible, this is so because the approach does not consider the anisotropy nor the evolution of GLE characteristics throughout the event. One can see that the fluence obtained with the "fast" analysis method agrees well with the full reconstruction.

Moments of the Fitted Proton Distribution
The available set of counting-rate profiles of the NM network is fitted with the intensity of arriving solar protons according to Equations 2 and 5: where the pitch angle α is measured from the axis of symmetry of the intensity distribution. The axis direction is determined in geographic coordinates, which is the instrument coordinate system of the NM network (Table 2), while for convenience of data interpretation, here we depict the axis direction in GSE coordinates (Figures 7(b) and (c)). The zero-order moment of the intensity distribution j (P , α, t) is the omnidirectional intensity: It is proportional to the proton number density N(P , t) = J (P , t)/v(P ), where v(P ) is proton speed. The number density of 1 GV protons is shown in Figure 7e. It reaches a maximum value in about 2.5 hours after the event onset. The next moment is the proton net flux: which is relevant to the proton-production profile at the Sun and is also affected by the interplanetary transport. The time profile of the proton net flux is plotted in Figure 7a. The time evolution of the spectrum slope δγ , which is responsible for the spectrum steepening with increase of rigidity, is shown with a color scale. In the late phase of the event, the spectrum is a pure power law, that is δγ = 0, typical for the delayed component of the GLE in contrast to the prompt component typically dominating the initial phase of the events (e.g. . In panels b and c of Figure 7, we present the GSE latitude ψ and longitude λ, of the direction of the proton flux arrival. We also plot the IMF direction, as measured by the ACE spacecraft (Smith et al., 1998) but shifted in time by the solar-wind transit time from the probe to the Earth's orbit. In Figure 7b, the R L bar additionally illustrates the passage time of a solar-wind structure of the Larmor radius scale of 1 GV protons, T L = R L /U sw , with U sw = 300 km s −1 and B = 4.5 nT. During the first ≈ 2.5 hours of the event, the magneticfield direction exhibited strong variations, including a switchback at t ≈ 160 min, while the direction of the proton flux, S, turned smoothly to the south of the ecliptic plane, and then returned back. In Figure 7d we plot the mean angle between the proton flux axis and the magnetic field, α SB . The mean value was obtained via temporal averages of α SB over a sliding time window extending from the current time t to the time t + 2R L /U sw (the bar over α SB illustrates the averaging time scale). The S-B (proton flux-magnetic field) angle increased from ≈ 50 • at the beginning of the event to 90 • at the flux-maximum time and then returned back. Thus, the GLE-producing solar protons arrived at the vicinity of Earth largely via the crossfield transport, perhaps from southern solar locations.
The average pitch-angle cosine of the GLE-producing solar protons, cos(θ ) = S(P , t)/J (P , t), characterizes the proton streaming, relevant to actual interplanetary transport conditions. This is shown in Figure 7e. The observed anisotropy was moderate, even at the beginning of the event, as discussed above.
The last moment shown here is Q, which is the cumulative net flux of energy E up to current time, i.e. fluence of energy of > 1 MV solar protons passing the Earth's orbit: It can be used as a measure of the GLE power as well as the relative power of its phases. For the time interval shown in Figure 7, 90% of the total energy fluence is accounted for by the prompt component emission.

Comparison with SOHO/ERNE Data
An event associated with GLE N • 73 SEP was observed by the ERNE instrument on board the SOHO spacecraft, which is a platform stabilized in such a way that the particle instrument axis points in the ecliptic plain at the GSE longitude either λ = 45 • , either λ = −45 • . On 28 October 2021 the orientation was at λ = −45 • (45 • W). The High Energy Detector (HED) of SOHO/ERNE has a 120 • -wide field of view, hence covering the plane of the ecliptic in the sector −105 • ≤ λ ≤ 15 • , therefore provides angular measurements of the particles outside the geomagnetosphere (Torsti et al., 1995).
In Figure 8, we compare the particle data of ERNE on board SOHO and NMs on Earth, on the one hand, and the ACE magnetic-field measurements on the other hand. In the solarwind flow, the SOHO is orbiting upstream of the ACE, while the Earth is downstream. For comparison sake of Figure 8, we shifted the ERNE data to a later time by the SOHO-ACE transit time of the solar wind, while the NM profiles are shifted to an earlier time by the solar-wind transit time from the ACE to the Earth's orbit.
Based on magnetic-field data, we selected several time intervals, i.e. t 1 , t 2 , and t 3 (Figure 8). The SEP rise was observed in the time interval t 2 -t 3 , when the interplanetary magnetic field was inside the ERNE/HED field of view (panel c). The rise phase of a deka-MeV proton event exhibited a clear velocity dispersion corresponding to a traveled distance of ≈ 1.2 − 1.5 AU (panel a). Then, during the time interval from 215 min to ≈ 245 min, a decrease was observed in all energy channels, implying the spacecraft entered a distinct magnetic flux tube filled with a different particle population. Shortly after the appearance of the distinct proton population, at around time t 3 , the magnetic-field magnitude and direction  (Figure 8b and c). Those changes may indicate an alteration of magnetic connection to the SEP source apparently situated at the western flank of the CME shock.
The GLE rise and maximum phase were observed in a very different magnetic environment (Figure 8b-e). At around the maximum intensity of the GLE-producing solar protons, from t 1 to t 2 , the interplanetary magnetic field dramatically deviated from the standard, spiral field direction, while the GV proton flux did not follow those fast changes, apparently because of the large Larmor radius of the GV protons.
For further comparison of the NM data with ERNE data, we created a virtual GV proton channel of ERNE/HED by sampling the GLE protons in the HED viewing cone, HED , that is a circular cone of 60 • half-width with axis in the ecliptic plane pointing at λ = −45 • : In Figure 9, such a channel profile is compared with the HED-observed profiles of deka-MeV protons. All profiles of Figure 9 are shifted back in time to a near-Sun particle source by subtracting the time it takes a proton to travel the distance of 1.4 AU, and then adding 8 minutes for possible comparison with data of solar electromagnetic emissions. It can be seen that the time-intensity profile of GeV protons is very different from the common profile of the deka-MeV protons. Note that the time-shifting technique is a kind of velocity-dispersion analysis and cannot completely deconvolve the interplanetary transport effect, especially after the event rise phase (e.g. Kocharov et al., 2015).
HED of SOHO/ERNE comprises both silicon detectors and a scintillator, allowing proton-flux anisotropy measurements in the deka -MeV energy range. For the present analysis of the proton-flux anisotropy, we divide the field of view of HED into the five sectors shown in the inset of Figure 10. Time-intensity profiles of protons in those sectors are plotted in panel a. In panel b, intensities observed in opposite sectors are compared using two anisotropy indices, A NORTH−SOUTH and A SUN−WEST , defined in the figure.
In the 28 October 2021 SEP/GLE event, the anisotropy of relativistic protons was surprisingly low compared to the anisotropy of deka -MeV protons. The anisotropy direction in the deka -MeV range and the anisotropy direction in the GeV range were also different. At lower energies, more particles arrived from the north and west, consistently with the observed direction of the interplanetary magnetic field (Figure 8). In contrast, the relativistic protons arrived preferentially from the south, similarly to the eruption center location on the solar disk.

Discussion and Conclusion
In the analysis presented here, we derived the spectral and angular shapes of high-energy SEPs during the first GLE event of Solar Cycle 25, namely GLE N • 73 that occurred on 28  October 2021. The detailed modeling of the global NM network response and the application of a method verified by direct space-borne measurements, allowed us to obtain precise information about the spectra of the SEPs with the best possible time resolution (1 minute during the event onset) available from the existing data sets, thus revealing with high resolution the characteristics of the prompt component and with lower resolution the delayed component. Several possible shapes of the spectra and PAD were studied, and the best fit was achieved with a modified power law in the spectra and a single Gaussian as the angular distribution, yet satisfactory results were obtained with an exponential rigidity spectrum for the prompt component. We revealed that the rigidity spectra were moderately hard (γ ≈ 4.5) during the event initial phase with a significant steeping δγ ≈ 1.1, which constantly decreased throughout the later stages of the event. The derived SEP spectra gradually softened during the event. Accordingly, the derived PAD was relatively wide in contrast to the two previous GLEs, namely with a distribution width σ 2 of about π . The revealed angular distribution may be a result of a perpendicular transport of the SEPs (e.g. Kocharov et al., 2005;Ruffolo et al., 2008), however, a more detailed modeling turns out to be necessary, and is planned as forthcoming work.
The uncertainties of the derived characteristics were explicitly assessed, and were larger than those for the previous weak GLEs N • 72, namely of the order of 30% for the spectra and PAD and about 15 degrees for the apparent source position, most likely due to larger NM count-rate fluctuations. The assessed confidence limits of the spectra and PAD were slightly greater than the systematic errors of the employed NM yield function, i.e. those related to the atmospheric cascade evolution (e.g. Alves Batista et al., 2019, and references therein).
The velocity-dispersion magnitude observed in the deka -MeV proton channels of ERNE/HED corresponded to the traveled distance of ≈ 1.4 AU. The average distance traveled by particles from the source to the detector depended on the magnetic line length and interplanetary scattering conditions. A traveled distance as small as 1.4 AU allowed only a weak scattering of deka -MeV protons to exist between the Sun and the Earth, corresponding to the proton mean free path > 2.5 AU (Figure 1 of Kocharov et al., 2015). How- ever, no strong anisotropy was seen in NM data: the average pitch-angle cosine of arriving high-energy protons cos(θ ) ≤ 0.35 (Figure 7e). This may be explained by the GV-proton transport across magnetic-field lines with corresponding mean free path ⊥ < . The latter implied a larger value of traveled distance and correspondingly, a somewhat earlier injection time of the 1 GV protons compared to the timing shown in Figure 9. By the onset time of the ERNE-observed deka -MeV proton emission, ≈ 16:10 UT, the CME had already expanded to heliocentric distances > 4R , so the deka -MeV protons could be accelerated by the CME shock in the solar wind. In contrast, emission of the NMobserved high-energy protons started before 15:50 UT. At that time, the CME was still low in the corona, below 1.5R above the chromosphere. This supports the idea of a coronal origin of the prompt component of GLE event(s).
The study and related discussion presented here give a reliable basis for understanding the nature of high-energy SEPs as well as for quantifying the corresponding space-weather effects (e.g. Jiggens et al., 2019).