PDF (3.058 MB) - Journal of Space Weather and Space Climate

1 downloads 0 Views 3MB Size Report
Received 31 January 2017 / Accepted 16 May 2017. ABSTRACT ..... 30-Nov-2011–05-Jan-2012 b ..... class flare over the eastern solar limb for event 66.
J. Space Weather Space Clim., 7, A14 (2017) DOI: 10.1051/swsc/2017013  M. Paassilta et al., Published by EDP Sciences 2017

OPEN

RESEARCH ARTICLE

ACCESS

Catalogue of 55–80 MeV solar proton events extending through solar cycles 23 and 24 Miikka Paassilta1,*, Osku Raukunen1, Rami Vainio1, Eino Valtonen1, Athanasios Papaioannou2, Robert Siipola1, Esa Riihonen1, Mark Dierckxsens3, Norma Crosby3, Olga Malandraki2, Bernd Heber4, and Karl-Ludwig Klein5 1

Department of Physics and Astronomy, University of Turku, 20014, Finland Corresponding author: [email protected] Institute for Astronomy, Astrophysics, Space Applications and Remote Sensing (IAASARS), National Observatory of Athens, I. Metaxa & Vas. Pavlou St. GR-15236 Penteli, Greece Royal Belgian Institute for Space Aeronomy, Avenue Circulaire 3, 1180 Uccle, Belgium Institut für Experimentelle und Angewandte Physik, Christian-Albrechts-Universität zu Kiel, 24118 Kiel, Germany LESIA, Observatoire de Paris, Section de Meudon, 5 place Jules Janssen, 92195 Meudon Cedex, France

* 2 3 4 5

Received 31 January 2017 / Accepted 16 May 2017 ABSTRACT We present a new catalogue of solar energetic particle events near the Earth, covering solar cycle 23 and the majority of solar cycle 24 (1996–2016), based on the 55–80 MeV proton intensity data gathered by the Solar and Heliospheric Observatory/the Energetic and Relativistic Nuclei and Electron experiment (SOHO/ERNE). In addition to ERNE proton and heavy ion observations, data from the Advanced Composition Explorer/Electron, Proton and Alpha Monitor (ACE/EPAM) (near-relativistic electrons), SOHO/EPHIN (Electron Proton Helium Instrument) (relativistic electrons), SOHO/LASCO (Large Angle and Spectrometric Coronagraph) (coronal mass ejections, CMEs) and Geostationary Operational Environmental Satellite (GOES) soft X-ray experiments are also considered and the associations between the particle and CME/X-ray events deduced to obtain a better understanding of each event. A total of 176 solar energetic particle (SEP) events have been identified as having occurred during the time period of interest; their onset and solar release times have been estimated using both velocity dispersion analysis (VDA) and time-shifting analysis (TSA) for protons, as well as TSA for near-relativistic electrons. Additionally, a brief statistical analysis was performed on the VDA and TSA results, as well as the X-rays and CMEs associated with the proton/electron events, both to test the viability of the VDA and to investigate possible differences between the two solar cycles. We find, in confirmation of a number of previous studies, that VDA results for protons that yield an apparent path length of 1 AU < s [ 3 AU seem to be useful, but those outside this range are probably unreliable, as evidenced by the anticorrelation between apparent path length and release time estimated from the X-ray activity. It also appears that even the first-arriving energetic protons apparently undergo significant pitch angle scattering in the interplanetary medium, with the resulting apparent path length being on average about twice the length of the spiral magnetic field. The analysis indicates an increase in high-energy SEP events originating from the far-eastern solar hemisphere; for instance, such an event with a well-established associated GOES flare has so far occurred three times during cycle 24 but possibly not at all during cycle 23. The generally lower level of solar activity during cycle 24, as opposed to cycle 23, has probably caused a significant decrease in total ambient pressure in the interplanetary space, leading to a larger proportion of SEP-associated halo-type CMEs. Taken together, these observations point to a qualitative difference between the two solar cycles. Key words. Solar energetic particles – Space weather – Solar cycles – Solar flares – Coronal mass ejections (CMEs)

1. Introduction Solar energetic particle (SEP) events, large injections of particles into interplanetary space from the Sun (see, e.g., Reames 1999, 2013), constitute an important component of space weather. They are closely associated with solar flare and coronal mass ejection (CME) activity, but the exact relationship between these and the particle acceleration processes remains an object of study. Aside from the valuable role that they have served in developing our understanding of the conditions and phenomena in the Sun and its atmosphere, as well as of particle transport and scattering in interplanetary space, SEPs also pose a threat in the form of tremendous increases in radiation dose rates. These may cause considerable harm to both manned and unmanned space missions, even rendering satellites inoperable. The crew and passengers of aircraft flying at high altitude or

high geographic latitude may also be exposed to a nonnegligible radiation hazard during the most energetic events (e.g. Vainio et al. 2009; Reames 2013; Mishev 2014, and references therein). So far, the Sun has been remarkably quiescent during the current solar cycle, number 24 (e.g. Richardson 2013; Gopalswamy et al. 2015a), which is here considered to have commenced in December 2008. This is in contrast to cycle 23 (1996–2008), which exhibited a great deal more of overall solar activity. A comparison between the two cycles in terms of SEP events and associated phenomena, such as X-ray flares and CMEs, is therefore called for, as it may add to our current knowledge regarding particle release and acceleration processes, as well as the interplanetary transport conditions. The purpose of this paper is twofold. Firstly, it presents a catalogue of energetic (55–80 MeV) solar proton events for

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

J. Space Weather Space Clim., 7, A14 (2017)

the years 1996–2016, supplemented with associated electron, electromagnetic and solar CME observations; secondly, these events are subjected to a statistical analysis, with focus on the methods for obtaining a solar release time for the energetic particles and on a rudimentary comparison of the two solar cycles. In both these regards, this article forms a follow-up study for Vainio et al. (2013), and it is also intended to complement similar existing work done on comparing solar cycles 23 and 24 (e.g. Chandra et al. 2013; Richardson et al. 2016). We have extended the period of interest to cover the years 2011–2016 and also revisited the proton velocity dispersion analysis (VDA), electron event onset times, as well as electromagnetic (soft X-ray) and CME observations for the events listed in Vainio et al. (2013). Based on data from the Energetic and Relativistic Nuclei and Electron (ERNE) experiment (Torsti et al. 1995), the Large Angle and Spectrometric Coronagraph (LASCO; Brueckner et al. 1995), and the Electron Proton Helium Instrument (EPHIN; Müller-Mellin et al. 1995) aboard the Solar and Heliospheric Observatory (SOHO) spacecraft, as well as the Electron, Proton and Alpha Monitor (EPAM; Gold et al. 1998) aboard the Advanced Composition Explorer (ACE) and the X-ray flux measuring instrumentation on the Geostationary Operational Environmental Satellites (GOES) in active operation during the period of interest, our catalogue comprises a total of 176 energetic solar particle events between 1996 and 2016. This time span covers solar cycle 23 in its entirety and more than the first half of solar cycle 24, which allows us to gain some insight into the differences of these two cycles. Particle release time analysis is performed for proton events using both VDA and time-shifting analysis (TSA) so as to obtain estimates for the particle release times near the Sun. We also investigate solar flare and CME activity for associations with the energetic particle releases and present some statistical comparisons and results between the previous and the ongoing solar cycle, along with discussion. Our goal has been to compile a comprehensive listing of solar proton events in the energy range of 55–80 MeV, which occurred during the aforementioned time period. The selection of this energy range was motivated by three principal reasons: (1) the significance of high-energy particles for space weather, (2) the relatively fast post-event decrease of intensities, enabling small, closely spaced events to be distinguished more easily than at low energies and finally (3) the desire to maintain continuity with a number of previous ERNE-related studies. It is important to note, however, that due to a number of gaps in the data coverage, as well as certain inherent limitations of the methodology applied to identifying and determining the events and their onset times, the event catalogue as presented herein is not entirely comprehensive. We nevertheless hope that it will both prove useful for further research and also stimulate it. In addition, we have been motivated in our work by the possibility that the apparent qualitative and quantitative changes from the previous solar cycle to the present one may lead to new discoveries and better understanding of the processes involved in SEP events. The European Space Agency (ESA) is establishing a Space Weather (SWE) Service Network1 in the frame of its Space Situational Awareness (SSA) programme (http://swe.ssa. esa.int/). The goal of this incentive is to support end-users in 1

Session 2 – SSA Space Weather Service Network, 13th European Space Weather Week meeting, 14–18 November 2016, Oostende, Belgium.

a wide range of affected sectors to mitigate the effects of space weather on their systems, reduce costs and improve reliability. The network consists of five Expert Service Centres (ESCs): Solar Weather, Heliospheric Weather, Space Radiation, Ionospheric Weather and Geomagnetic Conditions. The domain of the Space Radiation ESC (R-ESC; http://swe.ssa. esa.int/space-radiation) covers the monitoring, modelling and forecasting of space particle radiation and micro-size particulates in the near-Earth space environment, as well as their effects on technological and biological systems. During the second period of the SWE segment of the SSA programme, new products are being provided to improve space weather services. The SEP catalogue described in this paper was recently integrated as a University of Turku (UTU) federated product under the umbrella of the R-ESC. The structure of this paper is as follows. In Section 2, we introduce the data used in this work and present the event catalogue. We describe our statistical analysis and its results in Section 3, and outline our conclusions and outlook in Section 4. 2. The proton event catalogue for the years 1996–2016 2.1. SOHO/ERNE proton data and event selection

The ERNE experiment (Torsti et al. 1995) aboard the SOHO spacecraft served as the primary instrument in this study. ERNE consists of two particle telescopes, the Low-Energy Detector (LED) and the High-Energy Detector (HED), which together are designed to cover the nominal energy range of some 1 MeV/nucleon to a few hundred MeV/nucleon for ions. For protons and helium ions, the upper limit of the energy range is about 140 MeV/nucleon. For the reasons explained in the previous chapter, the proton energy channel of 54.8–80.3 MeV (average energy 67.7 MeV) was chosen as the principal channel to be investigated. To identify the SEP events that had occurred during the time of interest, we scanned visually through SOHO/ERNE intensity data collected during the period between May 1996 and December 2016. Our event selection criterion was such that the 1-min average intensity in the 54.8–80.3 MeV ERNE proton channel was required to surpass the quiet-time background of the relevant phase of the solar cycle by a factor of about 3. The quiet-time background intensities – which are mostly the product of galactic cosmic rays – were estimated to range from ~ 5 · 104 cm2 s1 sr1 MeV1 in 1996–1997 near the solar minimum to ~ 3 · 104 cm2 s1 sr1 MeV1 in 2001–2003 near the solar maximum during solar cycle 23. For solar cycle 24, the corresponding values in 2009–2010 were ~ 7 · 104 cm2 s1 sr1 MeV1 and ~ 3.5 · 104 cm2 s1 sr1 MeV1 at the end of 2013. The end times for the proton events were estimated, as well, for all energy channels. In our study, an SEP event is defined as having ended when the intensity of protons in the 12.6–13.8 MeV (average 13.3 MeV) energy channel falls below twice that of the pre-event background2 of this energy channel for the first time. Observing the onset and end in 2 This was derived by calculating the mean intensity for the 2 h immediately preceding the event onset. In cases where the data were unavailable or indicated that the previous event had not yet ended, estimated quiet-time background intensity, calculated separately for each semiannual period, was substituted for the pre-event background.

A14-p2

M. Paassilta et al.: Catalogue of 55–80 MeV solar proton events extending through solar cycles 23 and 24 Table 1. Major (longer than ~18 h) continuous ERNE data gaps. All times are UT. 28–31-May-1996 03-Jun-1996 15–16-Feb-1997 25-Feb–04-Mar-1997 25–28-Oct-1997 19–21-Nov-1997 27-Nov-1997 11–12-Dec-1997 04–05-Feb-1998 25–28-Feb-1998 24-Jun–09-Oct-1998 14–21-Nov-1998 01–03-Dec-1998 21-Dec-1998–08-Feb-1999 14–18-Feb-1999 21–26-May-1999 03–04-Jun-1999 19–20-Aug-1999 15–20-Nov-1999 28–29-Nov-1999 01–02-Dec-1999 07–08-Jan-2000 22–25-Feb-2000 13–31-Mar-2000 18–19-Apr-2000 22–24-May-2000 14–15-Jan-2001 02–03-Jul-2001

28–29-Jul-2001 10–17-Aug-2001 01–03-Jan-2002 24–25-Jan-2002 05–12-Feb-2002 24-May-2002 01–02-Jul-2003 14–15-Jan-2003 28-Feb–11-Mar-2003 05-Jul-2003 09–10-Jul-2003 22–24-Oct-2003 19–22-Jan-2004 26-Mar-2004 02-Apr-2004 22–29-Apr-2004 06–12-Aug-2004 08–09-Dec-2004 22–23-Dec-2004 31-Jul–03-Aug-2005 08–09-Mar-2006 13–14-Aug-2006 31-Aug-2006 27–29-Jan-2007 11-Nov-2008 20-Feb-2009 07-Sep-2006 12–15-Sep-2006a

26–27-Sep-2010 10–15-Dec-2010 09–13-Jan-2011 13–14-Jan-2011 07-Apr-2011 17-Jul-2011 21-Jul-2011 30-Nov-2011–05-Jan-2012 b 22–24-Jan-2012 25–27-Jan-2012 28-Jan–10-Feb-2012c 15–16-Feb-2012 04-May-2012 07-May-2012 12–14-May-2012 22-May-2012 07–08-Jul-2012 02–07-Nov-2012 06–07-Dec-2012 09-Dec-2012–31-Jan-2013d 01–05-Feb-2013 06–07-May-2013 10–13-May-2013 08–10-Jun-2013 17–22-Jul-2013 13–14-Aug-2013 07–09-Sep-2013 29–30-Sep-2013

14–15-Oct-2013 29-Oct–04-Dec-2013e 29–31-Jan-2014 21-Mar-2014 21-Sep-2014 02–03-Nov-2014 17–19-Dec-2014 09–15-Jan-2015 24–26-Jan-2015 12–14-Feb-2015a 14-Feb-2015 23–24-Mar-2015 05-Jun-2015 04, 05-Sep-2015 10-Sep-2015 02-Oct-2015 06–09-Nov-2015 14–16-Nov-2015 30-Nov–03-Dec-2015a 07-Dec-2015 10–11-Dec-2015 05–08-Jan-2016 30–31-Mar-2016 22–25-Jul-2016 06–07-Oct-2016 18-Oct-2016 21–22-Dec-2016

a

Some data are avaiable during this period. Data available for less than 20 min on 07-Dec-2011. c Data available for a few minutes on 30-Jan-2012. d Some 12 min of data available on 11-Dec-2012. e Data available for a few minutes on 30-Oct-2013. b

different energy channels is justified by the fact that in a typical event, after the most energetic particles have already arrived at and passed the orbit of the Earth, considerable numbers of slower, less energetic particles are still on the way. To consider only the greater energy channel would mean to exclude these slower protons which nevertheless contribute significantly to the total ion fluence and energy release of the event and may potentially cause harm to spacecraft and exposed human beings, such as astronauts. If another high-energy proton event occurs before the previous one has ended according to the criterion explained above, the onset of the next event is taken as the end of the preceding one. It is important to note that the data collection or transmission from the SOHO spacecraft to the Earth was interrupted on several occasions, most notably between the late June and early October 1998, as well as late December 1998 and early part of February 1999. Any SEP events during such periods will have gone unrecorded, as well as some very minor events in the immediate aftermath of extremely large ones that have caused HED to saturate. However, as the saturation can only occur at the height of large events, its overall effect on event selection and onset timing is negligible. A comprehensive listing of ERNE data gaps is given in Table 1. 2.2. Onset, solar release and event end time determination

The proton event onset times as observed by SOHO/ERNE were determined using the so-called Poisson-CUSUM (Cumulative sum) method. This essentially measures the

statistical quality of a process to decide whether or not the process is under control and determines the moment of time of a failure occurring. Here, the pre-event background level proton flux is analogous to a controlled process, while a flux increase at event onset corresponds to a departure from control. The onset determination algorithm and the criteria used have been previously described in detail in Huttunen-Heikinmaa et al. (2005). When the observed event onset time is known, the velocity dispersion analysis for a given particle species and kinetic energy E can be performed. The observed onset time at 1 AU can be written as ð1Þ tonset ðEÞ ¼ t0 þ 8:33 ½min =AU s b1 ðEÞ; where t0 is the particle release time (in minutes) from the acceleration site, s is the apparent path length (in AU) travelled by the particles and b1(E) is their reciprocal speed in units of c1. Linear fitting of the observed onset times as a function of the inverse speed thus yields an estimate for both t0 and s. The basic assumptions of VDA include a simultaneous release of particles of all energies and the same apparent path length for all particles. The latter assumption requires that the first-arriving particles either undergo no scattering, or that the scattering affects these particles independently of their kinetic energy. The VDA was performed using 20 proton energy channels, spanning the range 1.58 MeV–131 MeV (see Table 2). Of this, LED covers the energies between 1.58 MeV and 12.7 MeV, while HED covers those between 13.8 MeV and 131 MeV,

A14-p3

J. Space Weather Space Clim., 7, A14 (2017) Table 2. ERNE energy channels used for proton VDA. Channel

LED

HED

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20

Energy range [MeV] 1.58–1.78 1.78–2.16 2.16–2.66 2.66–3.29 3.29–4.10 4.10–5.12 5.12–6.42 6.42–8.06 8.06–10.1 10.1–12.7 13.8–16.9 16.9–22.4 20.8–28.0 25.9–32.2 32.2–40.5 40.5–53.5 50.8–67.3 63.8–80.2 80.2–101 101–131

Average energy [MeV] 1.68 1.97 2.41 2.98 3.70 4.71 5.72 7.15 9.09 11.4 15.4 18.9 23.3 29.1 36.4 45.6 57.4 72.0 90.5 108

Reciprocal speed [c1] 16.7 15.5 14.0 12.6 11.3 10.0 9.10 8.15 7.24 6.47 5.59 5.06 4.57 4.11 3.69 3.32 2.99 2.70 2.44 2.26

both providing 10 individual channels. The time resolution was 1 min for all channels. In the majority of cases, one or more data points were discarded from the VDA fitting. This was mainly due to an elevated background from a previous event, rendering no additional enhancement detectable at low energies, or the fact that no noticeable intensity enhancement at the highest energies was observed. A reasonable velocity dispersion relation could not be derived for all events, and in some cases, a statistically acceptable fit exists but the resulting apparent path length is unphysical, namely less or very much more than 1 AU. However, as the continuing evaluation of the performance of VDA is one of the objectives of this work, these cases have not been a priori excluded from the listing or the statistical analysis. The time-shifting analysis (TSA) for protons and electrons was performed by backshifting the previously determined onset (near the Earth) to the vicinity of the Sun. The release time for a particle species with kinetic energy E can be expressed as follows: trel ðEÞ ¼ tonset ðEÞ  8:33 ½min =AUL b1 ðEÞ:

ð2Þ

Here, b1 is again the reciprocal speed of the particles, and L is the length of the magnetic field line connecting the source and the observer, computed from the detected speed of the solar wind uSW during the event as LðuSW Þ ¼ zðrSC Þ  zðr Þ;

ð3Þ

where z(r) is the distance along the Archimedean spiral from the centre of the Sun, rSC is the radial distance of the observing spacecraft from the Sun and r is the solar radius (taken to be 6.957 · 105 km). rSC was approximated by calculating the Earth-Sun (centre-to-centre) distance and subtracting from this the Sun-directional distance of the spacecraft from the centre of the Earth. The solar wind speed, relevant to the calculation of z(r), was determined using data from the Solar Wind Electron,

Proton and Alpha Monitor (SWEPAM) aboard ACE (McComas et al. 1998; http://www.srl.caltech.edu/ACE/ASC/ level2/lvl2DATA_SWEPAM.html), with Wind/SWE data (Ogilvie et al. 1995; http://omniweb.gsfc.nasa.gov/ftpbrowser/ wind_swe_2m.html) substituted for the events for which the former were not available. In each case the data were averaged over a 12-h period, centred on tonset for 54.8–80.3 MeV protons observed by ERNE; in cases where the precise tonset for protons was unknown but known for electrons, the UT calendar day (24 h) during which the onset had most likely occurred was used to obtain L for electron TSA. The Earth-Sun and Earthspacecraft distances were averaged over the period of interest in the same manner. In addition to proton data, we examined electron intensity data recorded by both ACE/EPAM (near-relativistic electrons) and SOHO/EPHIN (relativistic electrons). These were analysed for onset times, and ACE/EPAM data additionally for maximum intensities during the events as well as particle release times, so as to form a more complete picture of each event than would be possible by relying on proton data alone. By default, the ACE/EPAM data (Gold et al. 1998) used in our analysis were sectored 1-min intensities observed by the LEFS603 (Low-Energy Foil Spectrometer) telescope. These data were checked for ion contamination by comparing them visually with proton intensities recorded by the LEMS120 (Low-Energy Magnetic Spectrometer) telescope, as well as with electron intensities recorded by DE30 (Deflected Electrons). If any indications of possible ion contamination were detected in the LEFS60 intensity data, DE30 data were used instead. While ion contamination occasionally presents an issue for electron data recorded by LEFS60, it has the advantage of having a larger geometric factor than DE30 (~0.40 cm2 sr for electrons for the former, compared to ~0.14 cm2 sr for the latter) and was for this reason preferred in our study. A similar electron intensity profile comparison to determine possible ion contamination is described in Malandraki et al. (2000). The authors performed their comparison on data recorded during an SEP event by the Ulysses/ HI-SCALE instrument which is virtually identical to ACE/ EPAM (the latter being the flight spare of Ulysses/HI-SCALE). The electron event onset times for ACE/EPAM were determined by first estimating the yearly background intensity (Ibg), calculating the sliding time-average intensity Iav(t) over 60 min, and then comparing Iav (t) and the corresponding standard deviation to data points following the time window by another 60 min. When five consecutive data points all showed an intensity equalling or exceeding Iav (t) + nr, where n was a user-defined constant (with the default value of 2.0), the time stamp of the first point was defined as the preliminary onset of the event. Gathering the intensity data points covering the range [61,60 + d] minutes before the preliminary onset – where d was a user-defined parameter, here set to either 180 or 240 – and subtracting Ibg from these, an exponential curve of the form I0 (t) = exp(A – Bt) was fitted with the least squares method to the intensity data. I0 (t) was then subtracted from the intensity data so as to remove trends, yielding detrended intensity Idetrend (t) = I(t)  Ibg  I0 (t). Finally, the onset search was performed using Idetrend (t) in a similar manner as that described above, with the exception that the comparison time 3

The numbers in the telescope names denote the angles (in degrees) of the boresight directions of the telescopes measured from the spin axis of the spacecraft.

A14-p4

M. Paassilta et al.: Catalogue of 55–80 MeV solar proton events extending through solar cycles 23 and 24

window was kept fixed at 61–120 min before the preliminary onset. The definitive onset was taken to be the time stamp of the first of five consecutive data points equal to or exceeding the 60-min pre-onset quantity Idetrend (t) + nr. The principal ACE/EPAM energy channel used for determining the onsets was 0.18–0.31 MeV. However, in five cases (events 1, 10, 41, 118 and 141) it was deemed necessary to resort to the 0.10–0.18 MeV energy channel so as to obtain a clear result, due to an elevated pre-event background or weak enhancement of the intensity in the higher energy channel. The energy ranges correspond to mean electron energies of 0.23 MeV and 0.13 MeV and thus mean speeds of 0.73 c and 0.61 c, respectively. These speeds are so high, and the number of available energy channels (four) so small, that attempting to perform VDA for ACE/EPAM electron data was not expected to produce reliable results in most cases. Instead, the onset times were determined for each event and, after discarding cases where no reasonable result for tonset could be obtained, only TSA was performed. Near-relativistic electron event information derived from ACE/EPAM observations is complemented by event onset times determined from SOHO/EPHIN electron data (MüllerMellin et al. 1995). The energy channel of interest was 0.7–3.0 MeV. The onset was defined in a similar, but slightly simpler, manner as for ACE/EPAM data: the average intensity over a specified time window Iav and the corresponding standard deviation r multiplied by a user-defined constant n were compared to immediately following data points, and the time stamp of the first data point to exceed the intensity value of Iav + nr (where in this case n = 4.0) was taken as the time of onset. No detrending or release time analysis was performed for these electron observations, but they are included in this article to continue the analysis of the data sets used in Vainio et al. (2013). The durations of the SEP events were estimated according to the basic criteria explained in Section 2.1 in the following manner. The preliminary event end times were derived by locating the time point where the proton intensity in the 12.6–13.8 MeV energy channel first fell below the pre-defined threshold value (2.0 · the pre-event average background intensity or the quiet-time background intensity, if the former was substantially elevated) after the 55–80 MeV proton event onset. The period between the onset and the preliminary event end was then scanned visually for data gaps and new lowenergy SEP events not visible in high energies. If any were found, the event ending time was shifted to the nearest noon or midnight UT before the next low-energy event or to the beginning of an extensive – i.e. longer than about an hour – data gap, if it was considered possible that significant new proton activity might have occurred during the gap. In cases where another 55–80 MeV proton event occurred before any of these conditions was met, the previous event was regarded as ending at the onset of the next. 2.3. Ion fluences and iron-to-oxygen ratios

Fluences (i.e., differential particle intensities integrated over a given range of energy and time, expressed here in cm2 sr1) of energetic protons and oxygen ions, as well as iron/oxygen ratios, were calculated for each event using ERNE data. For protons, the energy range of interest was from ~10 MeV to 140 MeV (prior to 19 April 2000) or 131 MeV (after 19 April 2000); the two different upper limits follow from a modification to the ERNE onboard software with a data format

change on the date given above. As the bulk of the ion fluence is expected to occur at low energies, however, this difference is unlikely to have a large impact on the results. The estimated quiet-time background was first subtracted from the intensities of each energy channel. Short time intervals with known issues with the ERNE instrument were treated as data gaps, and these, together with actual (non-critical) data gaps, were compensated by logarithmic interpolation using neighbouring ‘‘good’’ data points. Because of the limitations of this method, it was not applied to events during which substantial loss of coverage at the event onset or intensity maximum had occurred. Due to the large uncertainties involved, no fluence estimate is given for these events. When applicable, the proton fluence results were checked against those given in Papaioannou et al. (2016) (based on GOES proton data), and when they differed by more than a factor of about 2, the event in question was subjected to further analysis to determine whether the obtained result was indeed reliable. Any results thus deemed unreliable, usually because of ERNE saturation during the event, are marked to that effect in the catalogue. The energy range of interest for oxygen and iron ions was 5–15 MeV/nucleon. Since the intensities of heavy elements are typically several orders of magnitude smaller than those of protons, the interpolation method used with proton intensities could not be employed here as such. Instead, we first calculated 4 He fluences in the energy range of 5–17 MeV/nucleon (this being the closest equivalent available to the heavy ion energy range mentioned above) in the same manner as proton fluences and then used the ratio between the measured and corrected helium fluences as a correction factor for the measured oxygen fluences. Helium was used as a proxy because it has a similar charge-to-mass ratio to oxygen ion (or the same, assuming both are fully ionized). Finally, the iron-to-oxygen ratios were calculated from the measured event-averaged intensities for these particle species. In cases where no counts for oxygen were measured for the duration of an entire event, the value given in the catalogue for oxgyen fluence is the oxygen onecount upper limit, and the Fe/O ratio is omitted. Conversely, when oxygen counts were measured during an event but iron counts were not, we calculated an upper limit for the Fe/O ratio using one-count upper limit for iron intensity.4 2.4. Soft X-ray and CME data

The tabulated information regarding X-ray flare events is based on GOES satellite data made available by the United States National Oceanic and Atmospheric Administration (NOAA; http://www.ngdc.noaa.gov/stp/space-weather/solardata/solar-features/solar-flares), with the flare/SEP event association deduced by us. The electromagnetic (EM) wavelength band of interest was 0.1–0.8 nm, i.e. soft X-rays. However, only flares with NOAA/GOES magnitude classification of C1.0 or greater (peak flux at least 106 W/m2) were taken into consideration. Flare onset, classification and location were in the majority of cases provided directly by NOAA listings, and the time derivative maximum of the soft X-ray flux was calculated based on the flux data smoothed with 5-min sliding average for use in the particle injection time estimation. In cases where the flare classification and onset information were available but the coordinates unavailable in the NOAA 4 Even though neither oxygen fluences nor the Fe/O ratios are discussed or analysed further in this work, they are of interest here since they are prominently featured in the online version of our SEP event catalogue and serve as the main topic of an upcoming study.

A14-p5

J. Space Weather Space Clim., 7, A14 (2017)

X-ray or H-alpha listings5, the derived position provided in the SolarSoft Latest Events Archive (http://www.lmsal.com/solarsoft/latest_events_archive.html) was given as the location of the flare. If the flare that was most likely associated with an SEP event preceded the observed proton onset by more than 6 h, the event was marked with the note ‘‘tp-tx > 6 h’’ in the tables and the flare omitted from any detailed statistical considerations. These results, while derived independently, were checked against the event information provided in Cane et al. (2010) and Richardson et al. (2014). The SEP event-related flares that occurred at the western solar limb or near enough behind it (western longitude between 90 and ~120) to have a NOAA classification and onset time available are marked as ‘‘long.  90’’ in the catalogue. These cases, including their location, were also confirmed using the listings in Cane et al. (2010) and Richardson et al. (2014). Lastly, SEP events with no clearly identifiable GOES flare carry the note ‘‘No GOES flare’’; note, however, that in most (if not all) cases there is farside flare activity involved, according to the articles mentioned above. The time derivative of the X-ray intensity provides a tool for studying electron acceleration in the solar atmosphere. In solar observations, there is a well-known correlation between hard X-ray and microwave intensity and the time derivative of the soft X-ray intensity during the impulsive phase of an X-ray event. Pointed out by Neupert (1968) and since expounded by several authors, this correlation is commonly called the Neupert effect and is thought to arise when non-thermal accelerated electrons lose energy via bremsstrahlung upon encountering dense surrounding plasma, rapidly heating it (see, e.g., Veronig et al. 2002 for details). Thus, we may reasonably assume that the peak of acceleration of near-relativistic electrons at the Sun occurs approximately when the time derivative of the soft X-ray intensity reaches its local maximum. Provided that open magnetic field lines are present at the acceleration site, this corresponds to the time of release of electrons into interplanetary space. The SOHO LASCO CME Catalog (http://cdaw.gsfc.nasa. gov/CME_list/) was used for identifying and listing the features (first appearance, estimated linear speed, width and central position angle) of the CMEs likely associated with each SEP event. After an estimate of the approximate time of particle injection had been derived using VDA and TSA results and refined using flare and radio frequency data, the CME Catalog was searched for relevant CME activity. If the listing included a CME with an estimated solar surface departure time differing from the radio and flare activity by less than ~1 h, it was analysed further and then selected for our catalogue, if its association with the SEP event appeared reasonably certain in light of other evidence (see Sect. 2.5 for an example). We list the first detection of the CME as observed on the SOHO/LASCO C2 coronagraph, with the time figure rounded to the nearest minute. In a few cases, the CME was not observed on C2; for these events, the first observation is recorded for the C3 (see notes for Table 4). For a detailed introduction to the CME Catalog and discussion, see Gopalswamy et al. (2009). While radio frequency observations are not included in the tables presented in this article, spectral plots of Wind/WAVES data6 were scanned visually during periods of SEP, CME and 5 Events 115, 116, 117, 127, 129, 139, 142, 151, 155, 156, 160, 168, 172, 174 and 175. 6 These are available at http://cdaweb.gsfc.nasa.gov; on the main page, select ‘‘Wind’’ and ‘‘Radio and Plasma Waves (space)’’, then ‘‘Wind Radio/Plasma Wave, (WAVES) Hi-Res Parameters’’.

flare X-ray activity to find out whether or not type III microwave bursts were detected close in time to both the onset of the candidate flare and the surface lift-off of the candidate CME. Such radio bursts are in fact a general characteristic of J 20 MeV proton events (Cane et al. 2002). The presence or absence of type III burst activity was used to determine whether or not there was a plausible association between the candidate flare, candidate CME and the SEP event itself. 2.5. Event analysis case study: event 3

In the following, we briefly describe the analysis performed on a typical SEP event and the derivation of the tabulated quantities. The event concerned is event 3, which occurred on 6 November 1997. All times mentioned below and in the rest of this work are UT, and all particle release time estimates include the approximate light travel time, 500 s. A visual scan of the 1-min ERNE 55–80 MeV proton intensity data identified an enhancement of about two orders of magnitude near noon on 6 November 1997 as a candidate event. After it was subjected to closer study, the PoissonCUSUM method yielded the onset time of 12:37 in this energy channel, as well as the corresponding onset time values for the ERNE VDA channels, which allowed VDA to be performed. The resulting apparent path length is 2.61 ± 0.24 AU and the estimated particle release time of 12:21 ± 00:17. On the other hand, TSA yields a magnetic field line length of 1.20 AU and estimated particle release time of 12:18. For electrons, the onset and release time results are remarkably similar. The onset of ACE/EPAM 0.18–0.31 MeV electrons occurs at 12:24, and that of SOHO/EPHIN 0.7– 3.0 MeV electrons at 12:23. For the former, TSA indicates a particle release time of 12:18. Thus, the available information for both protons and electrons suggests a particle release at or very near 12:20. Turning to the SOHO LASCO CME Catalog, we next searched the time period near noon, 6 November 1997, for CME activity. An immediately obvious candidate for the SEP event-associated CME is a fast halo event first detected by LASCO at 12:11; in fact, no other CME is listed as having been detected for several hours either before or after this one. A linear fit to observation data points presented in the SOHO LASCO CME Catalog gives the speed of the CME as 1556 km/s and the departure time from the solar surface as 11:30. This is corroborated by a visual scan of Wind/WAVES data, which indicates considerable type III radio burst activity beginning at circa 11:52. NOAA listings show that during a few hours within noon of 6 November 1997, two soft X-ray flares with NOAA classification greater than C1.0 were detected by GOES: a C4.7 class flare, beginning at 11:31 and ending at 11:44; and an X9.4 class flare, beginning at 11:49 and ending at 12:01. The location of the X9.4 class flare is given in the listing as 18S 63W, and our analysis indicates the maximum of the time derivative of the soft X-ray flux as having occurred at 11:53. The larger and later of the two flares coincides with the estimated departure from the solar surface of the CME, as well as the derived particle release times for both 55–80 MeV protons and 0.18–0.31 MeV electrons. Therefore, combining the results mentioned above, it seems rather safe to conclude that the SEP event, the CME, the radio activity and the X class flare are very likely associated with one another. In addition, it is noted that the SEP event information given in Cane et al. (2010) is in agreement with this interpretation.

A14-p6

M. Paassilta et al.: Catalogue of 55–80 MeV solar proton events extending through solar cycles 23 and 24

Fig. 1. The particle and X-ray observations related to event 3 (6 November 1997). Shown here are proton and electron intensities (A, the former multiplied by 103), the time derivative of the soft X-ray intensity (B) and the soft X-ray intensity (C). The moment of the maximum of dI/dt for X-ray intensity is marked with a vertical red line. The particle intensities have been backshifted in time by (L/v – 500) seconds, where L = 1.20 AU and v = the mean speed of the particle species of interest.

Finally, the end of this proton event is estimated to have occurred at 22:26 on 13 November 1997. While ERNE proton intensity data are available for only 87.3% of the time between event onset and end, none of the data gaps are fortunately located at the critical rise and peak stages of the event, nor are there long continuous periods of missing data. Thus, estimates for the proton and oxygen fluences and the iron/oxygen ratio could be calculated with satisfactory certainty in this case. An overview of event 3 is presented in Figure 1. It shows backshifted proton and electron intensities juxtaposed with X-ray observations; the time derivative maximum of the soft X-ray intensity is marked with a red line. 2.6. The event catalogue

The proton event catalogue is presented in Tables 3 and 4, provided as Supplementary material, available at http://www. swsc-journal.org/10.1051/swsc/2017013/olm. Table 3 lists the proton events, organized by date, and provides the reference number (ID), date and time of detected 55–80 MeV proton onset, approximate maximum intensity Imax (in units of pfu MeV1, where pfu = cm2 s1 sr1), and event end date and time; the spiral field line length L (in AU) and the TSA-derived particle release time trel; the apparent path length s (in AU), release time t0 together with its standard error, and the square of the sample correlation coefficient R2 derived from VDA for 55–80 MeV protons; ERNE proton and oxygen fluences for > 10 MeV and 5–15 MeV/nucleon (in cm2 sr1), respectively; and iron-to-oxygen ratios for 5–15 MeV/nucleon are also given here. Table 4 lists the corresponding electron events together with flare and CME information. It shows the event reference number (ID), electron event onset date and

time, estimated maximum electron intensity Imax and TSAderived release time t0 for 0.18–0.31 MeV electrons; the onset time for 0.7–3.0 MeV electrons; the time of maximum slope of the soft X-ray intensity (Max. dI/dt), NOAA magnitude classification and solar location of the event-related GOES flare; and the first observation time, speed v (in km/s), and width and position angle (in degrees) of the event-associated CME. All particle release times given in the tables include an additional 500 s to account for the light travel time to the observer. For six proton events, the accompanying electron event could not be identified with satisfactory certainty, mainly due to the overall weakness and slow rise of the electron intensity enhancement. In three cases (events 119, 132 and 144), the onset criterion for ACE/EPAM electron intensity was triggered some 7 h or more before the proton event onset, suggesting that the electrons probably originated from a different injection; electron onset times for these entries were omitted from Table 4. In one case (event 149, 30 September 2013), the onset of the proton event was unknown owing to a data gap, rendering further analysis of the electron and X-ray/CME observations (aside from a rough estimation of the maximum intensity of electrons) impossible. Also note that for two events (numbers 71 and 91, 17 April 2002 and 20 November 2003, respectively), the first near-relativistic electrons to arrive at ACE may be from injections that precede the main proton event. A similar situation, concerning relativistic electrons, may have occurred for SOHO/EPHIN in events 41, 65, 73, 91, 94, 127 and 136. In each of these cases, there appear two or more successive electron intensity rises, but the intensity remains elevated to such an extent for several hours that the onset determination algorithms for the respective instruments are unable to distinguish between the apparent events.

A14-p7

J. Space Weather Space Clim., 7, A14 (2017)

Thus, for these catalogue entries, the electron onset time given refers to the earliest possible moment of actual event onset and is prefixed with ‘‘’’ in Table 4. Altogether there were 34 cases in which no visible disk or over-the-limb GOES flare could be associated with the SEP event with confidence; in another eight cases the likely candidate flare had occurred more than 6 h prior to the proton event onset and was for this reason not listed. In addition, four cases could not be investigated due to missing data, and two (events 56 and 66, discussed further below) remained uncertain in this respect due to very long delays between particle release and event onset near the Earth. This leaves 128 events for which at least some information of the associated GOES flare is available and the flare identification is relatively secure. Thirteen of these flares certainly or very likely occurred at or immediately behind the solar limb; thus, full flare information (time of the maximum of the intensity time derivative, classification, location) is listed for 115 events (72 during 1996–2008, 43 events during 2009–2016). Independent longitude estimation was not attempted for the farside flares. The observation that SEP events and CMEs are associated and that the peak intensities of the events are correlated with CME speed is well established (see, e.g., Kahler et al. 1984; Kahler 2001). It is therefore unsurprising that in our data set, virtually every SEP event that extended to proton energies above 55 MeV was accompanied by an identifiable CME. In events 56 (10 August 2001) and 66 (10 January 2002), the proton intensity increase was extremely slow and prolonged in high energies, and the particle injection that caused these events had possibly occurred on the previous calendar day, if not before. As we could not confirm the timing and indeed the identity of the relevant particle injection to our satisfaction, no CME or flare is listed for these events. We note, however, that Cane et al. (2010) suggest a slow partial halo CME (speed 479 km/s) first observed at 10:30 on 9 August 2001 and a possible farside flare for event 55, and a CME (speed 1794 km/s) first observed at 17:54 on 8 January 2002 together with an M1 class flare over the eastern solar limb for event 66. In contrast, Papaioannou et al. (2016) list event 66 as being associated with a C7 class GOES flare occurring at 9:27 on 10 January 2002 but do not offer a suggestion for the accompanying CME. In addition to these two outstanding cases, SOHO/LASCO data are missing for five events, and there are otherwise insufficient data to make a positive determination as to the particle event-associated CME in three cases. All the other 166 events exhibit unambiguous CME activity. Of these, 136 (91 in 1996– 2008, 45 in 2009–2016) were also known to be associated with a GOES flare on the visible solar disk or over the western solar limb (the cases of missing X-ray intensity data are excluded from these figures). Estimates for both the speed and apparent angular size are available for all CMEs listed in our catalogue, aside from seven in 1996–2008, for which the SOHO LASCO Catalog only reports a lower limit for the size, in addition to the estimated speed. 3. Statistical results and discussion 3.1. VDA results

We investigated 114 events of solar cycle 23 and 62 events of solar cycle 24 (until the end of the year 2016), a total of 176 events. A VDA result was found for all but three events of cycle 23 and 12 events of cycle 24, leaving 161 events

(111 for cycle 23, 50 for cycle 24) available for analysis. The statistical considerations presented in this subsection are similar in nature and given in approximately the same order as in Vainio et al. (2013) so as to facilitate comparison of results and show progress made since the publication of that paper. The algorithm used for onset time determination may completely fail in some circumstances. Such circumstances are, e.g., high pre-event background, very slow rise of the intensities, large sudden fluctuations in the intensities just before the onset and low overall fluxes (insufficient sensitivity of the instrument). Unphysical results may be caused by one or a few onset times given by the algorithm, which are clearly erroneous. Removing these data points from the VDA may give physically more significant results, but of course also makes these particular results more subjective in the sense that one has to decide which data points to remove from the analysis. The cases with no VDA result were omitted from consideration. Additionally, there were nine cases of very long (>5 AU) and seven cases of very short (7 AU are omitted from the subfigures but are included in the line fits.

A14-p10

M. Paassilta et al.: Catalogue of 55–80 MeV solar proton events extending through solar cycles 23 and 24

Fig. 4. The distribution of the ratio of the spiral field line length to the apparent path length obtained from VDA for the 1996–2016 events (161 in total). This quantity represents an estimate of the effective value of l of the first-arriving particles. The left-hand bar in each speed bin represents the solar cycle 23 (1996–2008) events, the middle bar the solar cycle 24 (2009–2016) events and the right-hand bar all events during 1996–2016. The error bars denote the statistical error, only, and the numbers on the abscissa denote the lower limit of the l bin. The last bin includes all the events for which L/s  1.

was performed in Vainio et al. (2013), where the fit line slope was also found to be near unity, possibly indicating a meaningful result. As the statistical significance of the fit was rather low, however, the authors suspected that this might be a fortuitous coincidence. When separate data sets for the solar cycles are considered, this indeed seems to be the case. This evidence suggests that the apparent path length is more dependent on other factors, such as the difficulties in onset determination caused by background intensity and the geometric factors of the ERNE particle telescopes not being constant with respect to incident particle energy, than on the spiral field line length. As the final part of this subsection, we examine the distribution of the ratio of the spiral field line length and the apparent path length, L/s. This quantity can be used to arrive at an estimate for the average value of the pitch angle for the earliest-arriving particles. If we write ds ¼ vdt ¼

vdl dl ¼ ; vl l

ð4Þ

where dl is the length element of the line connecting the particle source and observer and l is the pitch angle cosine of the particles, and then integrate both sides of the equation, we see that l ð5Þ s¼ ;  l where l is the actual distance along the field line between the  is the average value of l for source and the observer and l the earliest detected particles on their path from the Sun to the detecting spacecraft near the Earth. Taking l = L, we  ¼ L=s. can express the effective value for l as leff ¼ l Figure 4 represents the distribution of L/s. Omitting the cases for which L/s > 1 as unphysical, the mean and standard deviation, respectively, are for the data sets as follows: 0.54 and 0.20 (all events), 0.55 and 0.20 (the 1996–2008 events), and 0.52 and 0.20 (the 2009–2016 events). This result is in good agreement with the results given in Vainio et al. (2013). It also underscores the fact, earlier found in simulation studies

and pointed out by e.g. Lintunen & Vainio (2004) and Sáiz et al. (2005), that even the earliest-arriving particles cannot be regarded as moving directly along the magnetic field line, challenging one of the basic assumptions of VDA. Overall, the mean and standard deviation values as such do not reveal any remarkable changes between the solar cycles. Despite the poor statistics of the 2009–2016 period, it would still appear that compared to solar cycle 23 as a whole, there have been very few events in 2009–2016 for which L/s > 0.7; in other words, holding to the interpretation that the ratio L/s indeed corresponds to the average pitch angle cosine, scattering tends to be more prominent in the cycle 24 events, during which time the Sun was much less active than in cycle 23. However, the difference between the distributions fails to meet the quantitative criteria for statistical significance7, so no definite conclusions as to the differences between the two solar cycles can be made from this result alone. 3.2. Flares and solar release times of particles

The second main area of our statistical study was the investigation of particle release times derived with VDA for protons and TSA for electrons, alongside the X-ray signature of flares. With the time derivative maxima of the soft X-ray intensity known, it was possible to approximate the solar release times of relativistic electrons with the times of these maxima and compare them to TSA release times. We also made this comparison for proton VDA release times, and additionally considered the release time difference as a function of both flare longitude and the longitudinal distance of the flare from the Parker spiral footpoint. 7

The statistical tests performed on the cycle 23 and cycle 24 data sets were Mann-Whitney U test (Mann & Whitney 1947) and twosample Kolmogorov-Smirnov test (in the latter, the cases in which L/s > 1 were not considered). The results of these tests are not shown here in detail.

A14-p11

J. Space Weather Space Clim., 7, A14 (2017)

Fig. 5. Difference between apparent electron release time (TSA) and the time derivative maximum of X-ray intensity as a function of flare longitude for the 1996–2016 events. Solar cycle 23 (1996–2008) events are marked in red, solar cycle 24 (2009–2016) events in blue. Negative abscissa values refer to eastern and positive values to western longitude. Note the outlying data point near the lower right-hand corner.

Of the 128 listed flares, 45 (33 during 1996–2008, 12 during 2009–2016) belonged to the X class, 69 (45 during 1996–2008, 24 during 2009–2016) to the M class and 14 to the C class (seven during both 1996–2008 and 2009–2016). Flares that occurred over the western limb – i.e. those listed as ‘‘long.  90’’ – include 1 X class, 10 M class and 2 C class flares; they are included in the comparison of proton VDA release times and flare X-ray flux time derivative maxima (see below in this subsection) but are not considered in detail in the context of electron release times. The largest SEP-related flare in our catalogue was classified X20.0; it occurred on 2 April 2001 and has been a subject of study and discussion in its own right.8 The flares associated with the 2009–2016 events are, on average, somewhat weaker than those associated with the 1996–2008 events: they constitute 27% of all X class flares, 35% of M class flares, and 50% of C class flares in Table 4, compared to the total share of some 34%. Since the TSA calculation of electron release times, as employed in this work, does not take into account the magnetic connection between the particle acceleration site and the observer but instead assumes a direct path for the electrons along the Parker spiral, comparing the TSA release times with the maxima of the time derivative of soft X-ray intensity offers a simple means of studying approximately how well a given event is connected with the Earth. Figure 5 demonstrates this by showing the difference between apparent electron release time and the maximum of the X-ray intensity time derivative as a function of flare longitude. Aside from a handful of outlying data points (these will be discussed further below), the flares are concentrated in the western longitudes, and for the great majority, the delay between the maximum of the X-ray intensity time derivative and the TSA-derived injection time is between zero and ~60 min. Figure 6 shows the 1996–2016 events for which the time delay is between 0 and 80 min and longitude is greater than 30E, with western over-the-limb events included, arranged 8 For instance, see Krucker et al. (2011) for a report on electromagnetic observations of this flare.

Fig. 6. The distribution of the 1996–2016 events with the location of origin more westerly than 30E and the difference between apparent electron release time (TSA) and the time derivative maximum of X-ray intensity within the range of [0, 80] minutes. The error bars denote the statistical error, and the values on the abscissa (eastern longitudes expressed as negative, western longitudes as positive) denote the lower limit of each longitude bin. The last bin contains the western over-the-limb events.

into 20-degree longitude-wise and 20-min time-wise bins. Considering only cases where the delay is between zero and 20 min, the entire western solar hemisphere from 10W to near the limb is almost evenly populated: four successive 20-degree bins, starting from 10 to 30 west, contain 10 (6), 8(5), 10 (7) and 9 (6) events for the entire data set (for solar cycle 23). If events with a delay between zero and 40 min are selected, the longitude range of 50–70 west is somewhat favoured over the others in both solar cycle 23 (15 cases) and the entire data set (20 cases). The flare locations in this longitude range are also expected to have the best magnetic connection with the Earth. The other two most populated bins are 10–30 west with 13 (eight) events and 30–50 west with 15 (8) events for the entire data set (for solar cycle 23). In contrast, events for which the delay is between 40 and 60 min are

A14-p12

M. Paassilta et al.: Catalogue of 55–80 MeV solar proton events extending through solar cycles 23 and 24

Fig. 7. The particle and X-ray observations related to the SEP event of 17 July 2012 (event 138). Shown here are proton and electron intensities (A, the former multiplied by 103), the time derivative of the soft X-ray intensity (B) and the soft X-ray intensity (C). The moment of the maximum of dI/dt for X-ray intensity is marked with a vertical red line. The particle intensities have been backshifted in time by (L/v – 500) seconds, where L = 1.13 AU and v = the mean speed of the particle species of interest.

proportionally more common in all the other bins than in the 50–70-degree one, which additionally does not contain any events with delay greater than 60 min. The statistics for the solar cycle 24 events, when considered separately, are too weak to allow any firm conclusions to be made, but the longitude distribution of the event-associated flares seems broadly similar in that set, as well. However, when events originating from the eastern solar hemisphere are considered, a difference between the solar cycles does appear. For solar cycle 23, our listing includes only one case where the associated flare occurred at an eastern longitude greater than 20, namely event 58 on 24 September 2001 (flare longitude 23E). In contrast, such events have been clearly more common and exhibited a generally greater angular separation from the nominal Earth-connected Parker spiral footpoint in solar cycle 24: event 127 on 22 September 2011 (89E), event 133 on 7 March 2012 (27E), event 151 on 25 October 2013 (73E), event 153 on 28 October 2013 (28E) and event 158 on 25 February 2014 (82E). If the cases where the flare onset precedes the proton event onset by more than 6 h are also considered, event 41 on 21 January 2001 and event 81 on 6 September 2002 (probably 41E and 28E,

respectively) were recorded for solar cycle 23, which are balanced by event 132 on 5 March 2012 and event 147 on 15 May 2013 (probably 54E and 64E, respectively) during solar cycle 24. This observation is made more notable by the fact that the total number of 55–80 MeV proton events recorded by ERNE in cycle 24 is about 2/3 of that in cycle 23. (However, we note that lower-energy SEP events with a far-eastern location of origin are not very uncommon; see, e.g., Cane et al. 2010.) Events with probable sources behind the eastern solar limb are quite rare in our listing, but they seem to follow a broadly similar pattern in that they have been at the very least equally frequent, if not more so, in solar cycle 24 than in solar cycle 23. While the timing of the particle injection causing event 66 (10 January 2002) is not certain (see Sect. 2.6), the suggested source longitude of 120E in Cane et al. (2010) may be possible; however, event 128 (3 November 2011) certainly originated at 152E, according to Richardson et al. (2014). Event 163 (1 September 2014), with which no GOES flare could be associated with confidence and which will be considered by us in more detail in an upcoming study, may be another case originating over or behind the eastern limb.

A14-p13

J. Space Weather Space Clim., 7, A14 (2017)

Fig. 8. Difference between apparent electron release time (TSA) and the time derivative maximum of X-ray intensity as a function of the connection angle (the longitudinal distance of the flare from the Earth-connected Parker spiral footpoint) for the catalogued events, with outlying data points (see text) omitted. The red symbols represent solar cycle 23 (1996–2008) events, and the blue symbols represent solar cycle 24 (2009–2016) events; the fit lines coincide nearly perfectly in this graphic. The lower equations refer to the linear fit for solar cycle 23 events, the upper ones to the linear fit for solar cycle 24 events.

A number of solar cycle 24 events in our catalogue apparently also originated from relatively far behind the western solar limb. Correlating our event listing with that of Richardson et al. (2014), it appears that at least six events in 2009–2016 have probably been associated with sources at longitudes of 130, one (that of event 146 on 24 April 2013) being as distant as 175W. Unfortunately, similar coverage of solar farside flare activity is not available for events prior to the beginning of the Solar Terrestrial Relations Observatory (STEREO) mission in late 2006, so a thorough comparison between the two cycles in this respect is at best difficult. So, while the overall relatively poor statistics compel a cautious approach to making definite conclusions, the evidence nonetheless points to solar cycle 24 being a more favourable period than solar cycle 23 for high-energy events to be detected when the particles involved originate from solar locations that are usually not well magnetically connected with the nearEarth observer. This, in turn, might possibly imply that very wide SEP events have been more characteristic of cycle 24 than cycle 23, but our data set alone is not sufficient to establish such a result. In two events, the X-ray intensity exhibits such a slow rise that its time derivative maximum occurs after the arrival of the first protons at the Earth. The most striking example is the SEP event of 17 July 2012: an M1.7 class flare begins at 12:03, ACE/EPAM detects the onset of an electron event at 14:33, followed by the detection of the onset of a 67.7-MeV proton event by SOHO/ERNE at 15:05. However, the time derivative of the X-ray intensity does not reach its maximum until 16:36, more than four and a half hours after the onset of the flare. This event can be seen as an outlier near the bottom right-hand corner in Figure 5. It is obvious that the electron release must have begun some time between the flare and the electron event onsets, probably close to the TSA-derived (actual) release time of 14:20. Figure 7 shows the event in question as observed in proton (SOHO/ERNE), electron (ACE/EPAM), and X-ray (GOES) intensity, with a common time axis. The extremely slow rise of the latter is quite pronounced here. However,

it should be noted that there is very low-intensity X-ray activity showing peaks in the time derivative of X-ray intensity also earlier than the identified peak, related to the global X-ray maximum of the flare. This could signal the acceleration and release of the first-observed particles in the event. Somewhat similar but much less conspicuous cases with probable particle injections and accompanying considerable X-ray activity preceding the main event-associated flare, consequently resulting in a ‘‘too early’’ particle onset, are event 69 (20 February 2002) for both protons and electrons, as well as events 38 (8 November 2000) and 76 (14 August 2002) for electrons. On the other hand, the outlying points for which the difference in time between the estimated electron release and the flare X-ray intensity reaching its maximum slope is positive and large, aside from the far-eastern events discussed earlier, very likely result from insufficient sensitivity of the electron event onset determination or a high pre-event electron background. These tend to lead to delayed detection of the onset. Approximating the solar rotation as rigid, we can estimate the longitudinal distance of the flare and the footpoint of the Parker spiral that connects the solar surface to the Earth at the SEP event onset, often called the connection angle, as follows: rSC : ð6Þ / ¼ /flare  x u SW

The time difference between the TSA-estimated electron release and the time derivative maximum of the X-ray intensity is plotted as a function of D/ in Figure 8. Note that the outlying data points, for which the connection angle is more than 100 or the absolute value of the time difference is more than 100 min, are omitted from this analysis. Again, the trend is that the better connected a flare/SEP release region is, the shorter is the delay between the electron release time as estimated with TSA and the X-ray intensity time derivative maximum. To investigate quantitatively the level of correlation between the longitudinal distance and the time difference,

A14-p14

M. Paassilta et al.: Catalogue of 55–80 MeV solar proton events extending through solar cycles 23 and 24

Fig. 9. Difference of apparent proton release time (as obtained from VDA) and the time derivative maximum of the soft X-ray intensity as a function of VDA apparent path length s for the catalogued events. The red symbols represent solar cycle 23 (1996–2008) events, and the blue symbols represent solar cycle 24 (2009–2016) events. The lower equations refer to the linear fit for solar cycle 23 events, the upper ones to the linear fit for solar cycle 24 events. One event, for which s > 13 AU, is not shown here.

a regression line was fitted to both solar cycle 23 and 24 data sets. When only the core population, shown in Figure 8, is considered in the fit, the results are similar for both cycles. While the statistical significance of the fit is not particularly great in either of these instances, and a lot of scatter is present in the data, the slopes and the intercepts of the lines are nevertheless sensible in that both are greater than zero and the latter lie in the range of 15–20 min. This indicates that in a very well-connected event, the maximum energy release in the flare/CME, as estimated from the X-ray intensity time derivative maximum9, would tend to occur some 20 min or less before the electron release, as estimated from TSA. Such a result appears reasonable, and it also suggests that using the time derivative maximum of the flare X-ray intensity is a basically valid method for approximating the time of global particle release. Figures 9 and 10 show the difference between apparent VDA proton release time and the time derivative maximum of soft X-ray intensity as a function of apparent path length; in the latter, the events for which the path length is greater than 3.0 AU or the absolute value of the time difference is greater than 200 min are omitted from the analysis. Linear fits to the plotted data points are also shown in these figures. Again, the points are strongly scattered, and the linear fits are not very significant statistically. It is seen that the events, which most likely have a reasonable path length (s falls between 1 AU and some 3 AU), also have fairly small release time differences only weakly dependent on the path length. In all considered data sets, this time difference tends to become increasingly negative as s increases. A plausible explanation is that VDA produces overestimated path lengths particularly for events with large s and consequently proton release times that are too early. 9 The relationship between CME kinematics and X-ray intensity profiles is discussed in Zhang & Dere (2006).

3.3. SEP event-associated CMEs

As the final part of our statistical analysis, we briefly examined the CMEs associated with the 55–80 MeV proton/flare events. The width and speed of the CMEs, as indicators of their kinetic energy and size, were compared between the solar cycles. The relationship between fast CMEs and SEP events has long been recognized. For instance, Reames (1999 and references therein) pointed out the essential nature of CME speed for particle acceleration in gradual SEP events, and Kahler (2001) derived a correlation between CME speed and SEP event peak intensity at the Earth. More recently, Dierckxsens et al. (2015) have estimated the mean speed of CMEs associated with SEP events (ranging in proton energy from a few MeV upwards) to be about 1000 km/s, based on the SOHO LASCO CME Catalog data. While noting that for halo CMEs, the likelihood of an SEP event occurring is more dependent on the magnitude of the CME-associated flare than CME speed, they report that SEP event occurrence probability increases for all CME widths as CME speed increases and that halo CMEs are significantly more likely to be associated with SEP events than non-halo CMEs, as defined in the SOHO LASCO CME Catalog. We first considered the angular sizes and speeds of the SEP-related CMEs listed in our catalogue. Omitting the seven cases with an uncertain size estimate, halo CMEs (extending a full 360 around the disk of the Sun as seen from the Earth) constitute 72 of 101 cases (72%) for the 1996–2008 events, 48 of 59 cases (82%) for the 2009–2016 events and 120 of 160 for all events with a well-documented CME (75%); this result again underscores the fact that halo CMEs are typical for SEP events. The average speed for the 1996–2008 CMEs was 1361 km/s, 1276 km/s for the 2009–2016 events and 1331 km/s for all events. These results are in good agreement with those reported by Papaioannou et al. (2016), where an average speed of some 1390 km/s is given for SEP-related

A14-p15

J. Space Weather Space Clim., 7, A14 (2017)

Fig. 10. Difference of apparent proton release time (as obtained from VDA) and the time derivative maximum of the soft X-ray intensity as a function of VDA apparent path length s for the catalogued events, with all data points that have s > 3 AU or the absolute value of the time difference in excess of 200 min. The red symbols represent solar cycle 23 (1996–2008) events, and the blue symbols represent solar cycle 24 (2009–2016) events. The lower equations refer to the linear fit for solar cycle 23 events, the upper ones to the linear fit for solar cycle 24 events.

Fig. 11. The angular sizes of SEP-associated CMEs for the catalogued events. The left-hand bar in each speed bin represents the solar cycle 23 (1996–2008) events, the middle bar the solar cycle 24 (2009–2016) events and the right-hand bar all events during 1996–2016. The error bars denote the statistical error. Seven CMEs, for which only an upper limit estimate of the angular size is listed, are omitted.

CMEs. The authors point out that most CMEs associated with SEP events are faster than 1000 km/s but that there are some slow (500 km/s) SEP-related CMEs, as well; again, our work confirms these observations. As the formation of a strong shock, facilitated by high CME speed, is thought to be essential for particle acceleration, slow CMEs associated with SEP events might warrant further study. Figure 11 represents the distribution of the angular size of the CMEs, and Figure 12 the distribution of their speed. The dominance of halo CMEs in both periods of interest is obvious. There appears to be a mild preference for what might be

termed ‘‘medium-speed’’ SEP event-associated CMEs (750– 1000 km/s) and an expected clear under-representation of slow (10 MeV GOES proton channel (14 events with maximum intensity greater than 1 pfu for both time periods). 4. Conclusions and outlook Studying SOHO/ERNE data, we identified 176 55–80 MeV proton events that occurred during 1996–2016. We determined their onset times at 1 AU, their maximum intensities, proton and oxygen fluences, and Fe/O ratio, and performed VDA and TSA for the events. The associated energetic electron events were also identified from ACE/EPAM data and their onsets determined. TSA was then performed for the electron events. Lastly, we attempted to identify X-ray flares and CMEs corresponding to the particle events, obtaining estimates for particle solar release times from the maxima of the time derivative of the soft X-ray intensity. This accumulated information is presented in the form of an event catalogue, in two parts (Tables 3 and 4). In all but a few cases, the identification and subsequent analysis of the electron event accompanying the proton event was successful. An associated CME could be identified for practically all of the SEP events considered by us, and the majority were also associated with a flare of NOAA classification C1.0 or greater, occurring on the visible solar disk and less than 6 h before the onset of the proton event near the Earth. We performed a statistical analysis on the VDA and TSA path length results, both for the entire event set as a whole and in terms of comparison between the solar cycles 23 and 24. In addition, the SEP-associated X-ray flares and CMEs were briefly investigated. Comparison of the two solar cycles was made somewhat difficult by the poor statistics of especially the latter cycle, but some conclusions can nevertheless

A14-p17

J. Space Weather Space Clim., 7, A14 (2017)

be made. For the most part, our results are consistent with those presented in Vainio et al. (2013), a previous study involving many of the same topics of interest. We found that proton VDA produces reasonable results – for the apparent path length s and the release time t0 – for most events, demonstrating that its underlying assumptions also appear to be reasonable, even if not absolutely correct. However, the results must be regarded as dubious in cases where the apparent path length is less than 1 AU or considerably more than about 3 AU. Provided that the apparent path length is within these realistic limits, there is no significant correlation between the apparent path length and release time t0, or between the spiral field line length L and the apparent path length s, as derived from VDA for protons. It appears likely, therefore, that s may be more affected by such conditions as the background intensities and the energy dependence of the geometric factors of the ERNE particle telescopes than the spiral field line length. As established in previous studies, these protons cannot be regarded as propagating without substantial scattering in the interplanetary medium. In fact, the average apparent path length for the first-arriving particles exceeds the spiral field line length by a factor of about 2. Furthermore, scattering of protons tended to be slightly more pronounced during solar cycle 24 (average L/s  0.52) than during cycle 23 (average L/s  0.55), but this apparent increase does not meet the criteria for statistical significance and thus cannot be considered a meaningful indication of difference between the solar cycles at present. A greater number of data points – obtainable by increasing the energy range to cover, for instance, all >25 MeV proton events – is needed to decide whether the two cycles might indeed differ in this respect. Other essential conclusions of our study can be summarized as follows: (i) Allowing for inaccuracies in determining the electron event onset, the maximum of time derivative of soft X-ray intensity appears to be a valid indicator for the global solar release time of energetic electrons, as evidenced by the time difference between it and TSAderived moment of electron solar release. As expected, the least values of the time difference tend to occur for the events that, judging by the solar longitude of the associated flare, have the best magnetic connection with the observer. However, a considerable amount of scatter is present in the data. (ii) SEP events with a considerable longitudinal distance between the associated flare and the footpoint of the Sun-Earth magnetic line were noticeably more frequent in cycle 24 than in cycle 23; for three well-documented events in our catalogue, the solar longitude of the flare was more than some 60 east, while no directly comparable case occurred in cycle 23. High-energy events with a far-eastern origin in cycle 24 also generally seem to outnumber similar events in cycle 23. This raises the possibility of an increase in angular width for these events, which could be a feature of the current cycle. (iii) The majority of CMEs associated with SEP events were of the halo type in both cycles, but their proportion was greater during cycle 24. This can be attributed to the lower total pressure in the heliosphere, allowing more CMEs to become halo events. This observation is in concurrence with a result previously reported by Gopalswamy et al. (2015b). As the lower total pressure

also facilitates the forming of a shock11, the CMEs of cycle 24 may have been more efficient in accelerating SEPs towards the Earth from poorly connected sites than those of cycle 23, accounting together with (ii) for the SEP events originating from locations with a poor expected magnetic connection to the observer. However, events with very high particle energies have been generally rare in cycle 24. It would be highly desirable to extend the study to cover the rest of solar cycle 24 as soon as it comes to an end. A comprehensive data coverage of several more years could improve the statistics and would also allow direct comparison between the respective rising and falling phases of the two cycles. Another potentially fruitful line of investigation might be to include energetic particle data from other spacecraft, such as the STEREO probes, so as to obtain a more comprehensive spatial coverage of the SEP events. However, a study of this sort would be practically limited to cover only cycle 24. The results detailed in our work may be freely utilized in future studies. An online version of the catalogue presented in this paper is available at http://swe.ssa.esa.int/web/guest/ utu-srl-federated; in addition to the tabulated information given here, it offers plots of various observable quantities over the events. Acknowledgements. The research described in this paper was supported by ESA Contract 4000113187/15/D/MRP. We would like to acknowledge and express our gratitude to the organizations and teams responsible for maintaining the data sources used in this article (SEPServer, SOHO LASCO CME Catalog, SolarSoft Latest Events Archive, NOAA/Solar-Terrestrial Physics at the National Centers for Environmental Information, Coordinated Data Analysis Web, the ACE Science Center and the GLE database). SOHO LASCO CME Catalog: this CME catalog is generated and maintained at the CDAW Data Center by NASA and The Catholic University of America in cooperation with the Naval Research Laboratory. SOHO is a project of international cooperation between ESA and NASA. O.R. wishes to thank the Vilho, Yrjö and Kalle Väisälä foundation for financial support. A.P. would like to acknowledge gratefully the invitation by R.V. and the support of the University of Turku (UTU) that made his working visit to UTU possible. The editor thanks Ian Richardson and an anonymous referee for their assistance in evaluating this paper.

References Brueckner, G.E., R.A. Howard, M.J. Koomen, C.M. Korendyke, D.J. Michels, et al. The large angle spectroscopic coronagraph (LASCO). Sol. Phys., 162, 357–402, 1995, DOI: 10.1007/BF00733434. Cane, H.V., W.C. Erickson, and N.P. Prestage. Solar flares, type III radio bursts, coronal mass ejections, and energetic particles. J. Geophys. Res. [Space Phys.], 107(A10), SSH 14-1–SSH 14-19, 2002, DOI: 10.1029/2001JA000320. Cane, H.V., I.G. Richardson, and T.T. von Rosenvinge. A study of solar energetic particle events of 1997-2006: their composition and associations. J. Geophys. Res. [Space Phys.], 115, A08101, 2010, DOI: 10.1029/2009JA014848. Chandra, R., N. Gopalswamy, P. Mäkelä, H. Xie, S. Yashiro, et al. Solar energetic particle events during the rise phases of solar cycles 23 and 24. Adv. Space Res., 52, 2102–2111, 2013, DOI: 10.1016/j.asr.2013.09.006. 11

This interpretation assumes similar mass densities in the heliosphere during both solar cycles in question.

A14-p18

M. Paassilta et al.: Catalogue of 55–80 MeV solar proton events extending through solar cycles 23 and 24 Dierckxsens, M., K. Tziotziou, S. Dalla, I. Patsou, M.S. Marsh, N.B. Crosby, O. Malandraki, and G. Tsiropoula. Relationship between solar energetic particles and properties of flares and CMEs: statistical analysis of solar cycle 23 events, Sol. Phys., 290(3), 841–874, 2015, DOI: 10.1007/s11207-014-0641-4. Gold, R.E., S.M. Krimigis, S.E. Hawkins III, D.K. Haggerty, D.A. Lohr, E. Fiore, T.P. Armstrong, G. Holland, and L.J. Lanzerotti, Electron, proton, and alpha monitor on the advanced composition explorer spacecraft. Space Sci. Rev., 86, 541–562, 1998, DOI: 10.1023/A:1005088115759. Gopalswamy, N., S. Akiyama, S. Yashiro, H. Xie, P. Makela, and G. Michalek, The mild space weather in solar cycle 24, 2015a, arXiv preprint arXiv:1508.01603. Gopalswamy, N., H. Xie, S. Akiyama, P. Mäkelä, S. Yashiro, and G. Michalek. The peculiar behavior of halo coronal mass ejections in solar cycle 24. AsJ, 804, L23, 2015b, DOI: 10.1088/2041-8205/804/1/L23. Gopalswamy, N., S. Yashiro, G. Michalek, G. Stenborg, A. Vourlidas, S. Freeland, and R. Howard. The SOHO/LASCO CME catalog. Earth Moon Planet, 104, 295–313, 2009, DOI: 10.1007/s11038-008-9282-7. Huttunen-Heikinmaa, K., E. Valtonen, and T. Laitinen. Proton and helium release times in SEP events observed with SOHO/ERNE. A&A, 442, 673–685, 2005, DOI: 10.1051/0004-6361:20042620. Kahler, S.W. The correlation between solar energetic particle peak intensities and speeds of coronal mass ejections: effects of ambient particle intensities and energy spectra. J. Geophys. Res., 106, 20947–20956, 2001, DOI: 10.1029/2000JA002231. Kahler, S.W., N.R. Sheeley, Jr., R.A. Howard, D.J. Michels, M.J. Koomen, R.E. McGuire, T.T. von Rosenvinge, and D.V. Reames. Associations between coronal mass ejections and solar energetic proton events. J. Geophys. Res., 89, 9683–9693, 1984, DOI: 10.1029/JA089iA11p09683. Krucker, S., C.L. Raftery, and H.S. Hudson, Temperature and density estimates of extreme-ultraviolet flare ribbons derived from TRACE diffraction patterns, Astrophys. J., 734(1), 34, 2011, URL http://stacks.iop.org/0004-637X/734/i=1/a=34. Lintunen, J., and R. Vainio. Solar energetic particle event onset as analyzed from simulated data, A&A, 420, 343–350, 2004, DOI: 10.1051/0004-6361:20034247. Malandraki, O., E.T. Sarris, and P. Trochoutsos. Probing the magnetic topology of coronal mass ejections by means of Ulysses/HI-SCALE energetic particle observations. Ann. Geophys., 18, 129–140, 2000, DOI: 10.1007/s00585-000-0129-4. Mann, H.B., and D.R. Whitney. On a test of whether one of two random variables is stochastically larger than the other. Ann. Math. Stat., 18(1), 50–60, 1947, URL http://www.jstor.org/stable/ 2236101. McComas, D.J., S.J. Bame, P. Barker, W.C. Feldman, J.L. Phillips, P. Riley, and J.W. Griffee. Solar wind electron proton alpha monitor (SWEPAM) for the advanced composition explorer. Space Sci. Rev., 86, 563–612, 1998, DOI: 10.1023/A:1005040232597. Mewaldt, R.A., C.M.S. Cohen, G.M. Mason, T.T. von Rosenvinge, G.C.W. Li Smith, and A. Vourlidas, Investigating the causes of solar-cycle variations in solar energetic particle fluences and composition. In: Proceedings of the 34th International Cosmic Ray Conference, The Hague, The Netherlands, 2015, URL https:// pos.sissa.it/archive/conferences/236/030/ICRC2015_030.pdf. Mishev, A.L. Computation of radiation environment during ground level enhancements 65, 69 and 70 at equatorial region and flight

altitudes, Adv. Space Res., 54, 528–535, 2014, DOI: 10.1016/j.asr.2013.10.010. Müller-Mellin, R., H. Kunow, V. Fleißner, E. Pehlke, E. Rode, et al. COSTEP – comprehensive suprathermal and energetic particle analyser. Sol. Phys., 162, 483–504, 1995, DOI: 10.1007/BF00733437. Neupert, W.M. Comparison of solar X-ray line emission with microwave emission during flares, ApJ, 153, L59, 1968, DOI: 10.1086/180220. Ogilvie, K.W., D.J. Chornay, R.J. Fritzenreiter, F. Hunsaker, J. Keller, et al. SWE, a comprehensive plasma instrument for the wind spacecraft, Space Sci. Rev., 71, 55–77, 1995, DOI: 10.1007/BF00751326. Papaioannou, A., I. Sandberg, A. Anastasiadis, A. Kouloumvakos, M.K. Georgoulis, K. Tziotziou, G. Tsiropoula, P. Jiggens, and A. Hilgers. Solar flares, coronal mass ejections and solar energetic particle event characteristics. J. Space Weather Space Clim., 6, A42, 2016, DOI: 10.1051/swsc/2016035. Reames, D.V. Particle acceleration at the sun and in the heliosphere, Space Sci. Rev., 90(3), 413–491, 1999, DOI: 10.1023/A:1005105831781. Reames, D.V. The two sources of solar energetic particles, Space Sci. Rev., 175, 53–92, 2013, DOI: 10.1007/s11214-013-9958-9. Richardson, I.G. Geomagnetic activity during the rising phase of solar cycle 24. J. Space Weather Space Clim., 3(27), A08, 2013, DOI: 10.1051/swsc/2013031. Richardson, I.G., T.T. von Rosenvinge, and H.V. Cane, 25 MeV solar proton events in cycle 24 and previous cycles, Adv. Space Res., 2016, DOI: 10.1016/j.asr.2016.07.035, in press, URL http:// www.sciencedirect.com/science/article/pii/S0273117716304331. Richardson, I.G., T.T. von Rosenvinge, H.V. Cane, E.R. Christian, C.M.S. Cohen, A.W. Labrador, R.A. Leske, R.A. Mewaldt, M.E. Wiedenbeck, and E.C. Stone. 25 MeV proton events observed by the high energy telescopes on the STEREO A and B spacecraft and/or at earth during the first ~ seven years of the STEREO mission, Sol. Phys., 289, 3059–3107, 2014, DOI: 10.1007/s11207-014-0524-8. Sáiz, A., P. Evenson, D. Ruffolo, and J.W. Bieber, On the estimation of solar energetic particle injection timing from onset times near earth, ApJ, 626, 1131–1137, 2005, DOI: 10.1086/430293. Torsti, J., E. Valtonen, M. Lumme, P. Peltonen, T. Eronen, et al. Energetic particle experiment ERNE. Sol. Phys., 162, 505–531, 1995, DOI: 10.1007/BF00733438. Vainio, R., L. Desorgher, D. Heynderickx, M. Storini, E. Flückiger, et al. Dynamics of the earth’s particle radiation environment. Space Sci. Rev., 147, 187–231, 2009, DOI: 10.1007/s11214-009-9496-7. Vainio, R., E. Valtonen, B. Heber, O.E. Malandraki, A. Papaioannou, et al. The first SEP server event catalogue ~68-MeV solar proton events observed at 1 AU in 1996–2010. J. Space Weather Space Clim., 3(27), A12, 2013, DOI: 10.1051/swsc/2013030. Veronig, A., B. Vršnak, B.R. Dennis, M. Temmer, A. Hanslmeier, and J. Magdalenic´, Investigation of the Neupert effect in solar flares. I. Statistical properties and the evaporation model, A&A, 392, 699–712, 2002, DOI: 10.1051/0004-6361:20020947. Zhang, J., and K.P. Dere. A statistical study of main and residual accelerations of coronal mass ejections, ApJ, 649, 1100–1109, 2006, DOI: 10.1086/506903.

Cite this article as: Paassilta M, Raukunen O, Vainio R, Valtonen E, Papaioannou A, et al. Catalogue of 55–80 MeV solar proton events extending through solar cycles 23 and 24. J. Space Weather Space Clim., 7, A14, 2017, DOI: 10.1051/swsc/2017013.

A14-p19