Search

Space weather is a multidisciplinary research area connecting scientists from across heliophysics domains seeking a coherent understanding of our space environment that can also serve modern life and society's needs. COSPAR's ISWAT (International Space Weather Action Teams) 'clusters' focus attention on different areas of space weather study while ensuring the coupled system is broadly addressed via regular communications and interactions. The ISWAT cluster "H3: Radiation Environment in the Heliosphere" (https://www.iswat-cospar.org/h3) has been working to provide a scientific platform to understand, characterize and predict the energetic particle radiation in the heliosphere with the practical goal of mitigating radiation risks associated with areospace activities, satellite industry and human space explorations. In particular, present approaches help us understand the physical phenomena at large, optimizing the output of multi-viewpoint observations and pushing current models to their limits. In this paper, we review the scientific aspects of the radiation environment in the heliosphere covering four different radiation types: Solar Energetic Particles (SEPs), Ground Level Enhancement (GLE, a type of SEP events with energies high enough to trigger the enhancement of ground-level detectors), Galactic Cosmic Rays (GCRs) and Anomalous Cosmic Rays (ACRs). We focus on related advances in the research community in the past 10-20 years and what we still lack in terms of understanding and predictive capabilities. Finally we also consider some recommendations related to the improvement of both observational and modeling capabilities in the field of space radiation environment.

Large solar eruptions are often associated with long-duration gamma-ray emission extending well above 100 MeV. While this phenomenon is known to be caused by high-energy ions interacting with the solar atmosphere, the underlying dominant acceleration process remains under debate. Potential mechanisms include continuous acceleration of particles trapped within large coronal loops or acceleration at coronal mass ejection (CME)-driven shocks, with subsequent back-propagation towards the Sun. As a test of the latter scenario, previous studies have explored the relationship between the inferred particle population producing the high-energy gamma-rays, and the population of solar energetic particles (SEPs) measured in situ. However, given the significant limitations on available observations, these estimates unavoidably rely on a number of assumptions. In an effort to better constrain theories of the gamma-ray emission origin, we re-examine the calculation uncertainties and how they influence the comparison of these two proton populations. We show that, even accounting for conservative assumptions related to gamma-ray flare, SEP event and interplanetary scattering modeling, their statistical relationship is only poorly/moderately significant. However, though the level of correlation is of interest, it does not provide conclusive evidence for or against a causal connection. The main result of this investigation is that the fraction of the shock-accelerated protons required to account for the gamma-ray observations is >20-40% for six of the fourteen eruptions analyzed. Such high values argue against current CME-shock origin models, predicting a <2% back-precipitation, hence the computed numbers of high-energy SEPs appear to be greatly insufficient to sustain the measured gamma-ray emission., Comment: Accepted for publication in the Astrophysical Journal

The forecasting of solar energetic particles (SEPs) is a prominent area of space weather research. Numerous forecasting models exist that predict SEP event properties at proton energies <100MeV. One of these models is the SPARX system, a physics-based forecasting tool that calculates >10MeV and >60MeV flux profiles within minutes of a flare being detected. This work describes SPARX-H, the extension of SPARX to forecast SEP events above 300MeV . SPARX-H predicts fluxes in three high energy channels up to several hundred MeV. Correlations between SEP peak flux and peak intensity of the associated solar flare are seen to be weak at high energies, but improved when events are grouped based on the field polarity during the event. Initial results from this new high energy forecasting tool are presented here and the applications of high energy forecasts are discussed. Additionally, the new high energy version of SPARX is tested on a set of historic SEP events. We see that SPARX-H performs best when predicting peak fluxes from events with source locations in well-connected regions, where many large SEP events tend to originate., Comment: Accepted for publication in Space Weather 13 pages, 5 figures

Large solar eruptive events, including solar flares and coronal mass ejections (CMEs), can lead to solar energetic particle (SEP) events. During these events, protons are accelerated up to several GeV and pose numerous space weather risks. These risks include, but are not limited to, radiation hazards to astronauts and disruption to satellites and electronics. The highest energy SEPs are capable of reaching Earth on timescales of minutes and can be detected in ground level enhancements (GLEs). Understanding and analyzing these events is critical to future forecasting models. However, the availability of high energy SEP data sets is limited, especially that which covers multiple solar cycles. The majority of analysis of SEP events considers data at energies <100 MeV. In this work, we use a newly calibrated data set using data from Geostationary Operational Environmental Satellite-high energy proton and alpha detector between 1984 and 2017. Analysis of the SEP events in this time period over three high energy channels is performed, and SEP properties are compared to flare and CME parameters. In addition, neutron monitor (NM) observations are examined for the relevant GLE events. We find that correlations between SEP peak intensity and the CME speed are much weaker than for lower SEP energies. Correlations with flare intensity are broadly similar or weaker. Strong correlations are seen between >300 MeV data and GLE properties from NM data. The results of our work can be utilized in future forecasting models for both high energy SEP and GLE events., Comment: 17 pages, 11 figures. Published in Space Weather

Corotation of particle-filled magnetic flux tubes is generally thought to have a minor influence on the time-intensity profiles of gradual Solar Energetic Particle (SEP) events. For this reason many models solve the focussed transport equation within the corotating frame, thus neglecting corotation effects. We study the effects of corotation on gradual SEP intensity profiles at a range of observer longitudinal positions relative to the solar source. We study how corotation affects the duration and decay time constant of SEP events and the variation of peak intensity with observer position. We use a 3D full-orbit test particle code with time-extended SEP injection via a shock-like source. Unlike with focussed transport models, the test particle approach enables us to switch corotation on and off easily. While shock acceleration is not modelled directly, our methodology allows us to study how corotation and the time-varying observer-shock magnetic connection influence intensity profiles detected at six observers. We find that corotation strongly affects SEP profiles, for a monoenergetic population of 5 MeV protons, being a dominant influence during the decay phase. Simulations including corotation display dramatically shortened durations for western events, compared to those which do not include it. When corotation is taken into account, for both eastern and western events the decay time constant is reduced and its dependence on the scattering mean free path becomes negligible. Corotation reduces the peak intensity for western events and enhances it for eastern ones, thus making the east-west asymmetry in peak intensity stronger, compared to the no-corotation case. Modelling SEP intensity profiles without carefully accounting for corotation leads to artificially extended decay phases during western events and profiles with a similar shape regardless of observer longitudinal position., Comment: A&A Letters in press

Solar Energetic Particle (SEP) acceleration and injection into interplanetary space during gradual SEP events is thought to take place at Coronal Mass Ejection (CME)-driven shocks. Features of measured intensity profiles at 1 au have been attributed to properties of the radial and longitudinal/latitudinal injections at the shock. Focussed transport models are typically used to model acceleration at a CME-shock and subsequent propagation. Test particle simulations are an alternative approach but so far they have been carried out only with instantaneous injection near the Sun. We develop the first temporally extended shock-like injection for our 3D test particle code and investigate how the spatial features of injection affect SEP intensity and anisotropy profiles for observers at 0.3 and 1.0 au. We conduct simulations of a monoenergetic population of 5 MeV protons considering three different radial injection functions and two longitudinal/latitudinal injection functions. We consider a range of scattering conditions with scattering mean free path values ranging from 0.1-1.0 au, and determine intensity and anisotropy profiles at six observers at different longitudinal locations. We find that the radial, longitudinal and latitudinal injection functions play a relatively minor role in shaping the SEP intensity profiles. The dependence of profiles on the value of the scattering mean free path is also weak, unlike what is found from 1D focussed transport models. Spatial factors, such as the time of observer-shock-connection/disconnection and time of shock passage have a much stronger influence on SEP intensities and anisotropies. Persistent anisotropies until shock passage are seen in our simulations. Comparing instantaneous and shock-like injections, we find that the link between duration of injection and of the SEP event is very weak, unlike what is commonly assumed., Comment: A&A in press

There are many difficulties associated with forecasting high-energy solar particle events at Earth. One issue is understanding why some large solar eruptive events trigger ground level enhancement (GLE) events and others do not. In this work we perform 3D test particle simulations of a set of historic GLEs to understand more about what causes these powerful events. Particular focus is given to studying how the heliospheric current sheet (HCS) affects high-energy proton transport through the heliosphere following an event. Analysis of $\geq$M7.0 flares between 1976$-$2020 shows that active regions located closer to the HCS ($<$10$^{\circ}$) are more likely to be associated with a GLE event. We found that modelled GLE events where the source region was close to the HCS also led to increased heliospheric transport in longitude and higher count rates (when the Earth was located in the drift direction). In a model that does not include perpendicular diffusion associated with turbulence, the HCS is the dominant mechanism affecting heliospheric particle transport for GLE 42 and 69, and varying other parameters (e.g. a narrow, 10$^{\circ}$, or wider, 60$^{\circ}$, injection width) causes little change. Overall in our model, the HCS is relevant in 71$\%$ of our analysed GLEs and including it more accurately reproduces observed intensities near Earth. Our simulations enable us to produce model profiles at Earth that can be compared to existing observations by the GOES satellites and neutron monitors, as well as for use in developing future forecasting models., Comment: ApJ accepted

Gamma-ray emission during long-duration gamma-ray flare (LDGRF) events is thought to be caused mainly by $>$300 MeV protons interacting with the ambient plasma at or near the photosphere. Prolonged periods of the gamma-ray emission have prompted the suggestion that the source of the energetic protons is acceleration at a coronal mass ejection (CME)-driven shock, followed by particle back-precipitation onto the solar atmosphere over extended times. We study the latter hypothesis using test particle simulations, which allow us to investigate whether scattering associated with turbulence aids particles in overcoming the effect of magnetic mirroring, which impedes back-precipitation by reflecting particles as they travel sunwards. The instantaneous precipitation fraction, $P$, the proportion of protons that successfully precipitate for injection at a fixed height, $r_i$, is studied as a function of scattering mean free path, $\lambda$ and $r_i$. Upper limits to the total precipitation fraction, $\overline{P}$, were calculated for eight LDGRF events for moderate scattering conditions ($\lambda$=0.1 au). We find that the presence of scattering helps back-precipitation compared to the scatter-free case, although at very low $\lambda$ values outward convection with the solar wind ultimately dominates. For eight LDGRF events, due to strong mirroring, $\overline{P}$ is very small, between 0.56 and 0.93% even in the presence of scattering. Time-extended acceleration and large total precipitation fractions, as seen in the observations, cannot be reconciled for a moving shock source according to our simulations. Therefore, it is not possible to obtain both long duration $\gamma$ ray emission and efficient precipitation within this scenario. These results challenge the CME shock source scenario as the main mechanism for $\gamma$ ray production in LDGRFs., Comment: Accepted for publication in A&A

The Geostationary Orbital Environmental Satellites (GOES) Soft X-ray (SXR) sensors have provided data relating to, inter alia, the time, intensity and duration of solar flares since the 1970s. The GOES SXR Flare List has become the standard reference catalogue for solar flares and is widely used in solar physics research and space weather. We report here that in the current version of the list there are significant differences between the mean duration of flares which occurred before May 1997 and the mean duration of flares thereafter. Our analysis shows that the reported flare timings for the pre-May 1997 data were not based on the same criteria as is currently the case. This finding has serious implications for all those who used flare duration (or fluence, which depends on the chosen start and end times) as part of their analysis of pre-May 1997 solar events, or statistical analyses of large samples of flares, e.g. as part of the assessment of a Solar Energetic Particle forecasting algorithm., Comment: In press in Space Weather (2018)

Understanding the transport of Solar Energetic Particles (SEPs) from acceleration sites at the Sun into interplanetary space and to the Earth is an important question for forecasting space weather. The Interplanetary Magnetic Field (IMF), with two distinct polarities and a complex structure, governs energetic particle transport and drifts. We analyse for the first time the effect of a wavy Heliospheric Current Sheet (HCS) on the propagation of SEPs. We inject protons close to the Sun and propagate them by integrating fully 3D trajectories within the inner heliosphere in the presence of weak scattering. We model the HCS position using fits based on neutral lines of magnetic field source surface maps (SSMs). We map 1 au proton crossings, which show efficient transport in longitude via HCS, depending on the location of the injection region with respect to the HCS. For HCS tilt angles around $30^\circ-40^\circ$, we find significant qualitative differences between A+ and A$-$ configurations of the IMF, with stronger fluences along the HCS in the former case but with a distribution of particles across a wider range of longitudes and latitudes in the latter. We show how a wavy current sheet leads to longitudinally periodic enhancements in particle fluence. We show that for an A+ IMF configuration, a wavy HCS allows for more proton deceleration than a flat HCS. We find that A$-$ IMF configurations result in larger average fluences than A+ IMF configurations, due to a radial drift component at the current sheet., Comment: 16 pages, 8 figures and 1 table, plus appendix of 2 pages, 1 figure and 1 table

Solar energetic particle (SEP) events are known to occur following solar flares and coronal mass ejections (CMEs). However some high-energy solar events do not result in SEPs being detected at Earth, and it is these types of event which may be termed "false alarms". We define two simple SEP forecasting algorithms based upon the occurrence of a magnetically well-connected CME with a speed in excess of 1500 km/s ("a fast CME") or a well-connected X-class flare and analyse them with respect to historical data sets. We compare the parameters of those solar events which produced an enhancement of >40 MeV protons at Earth ("an SEP event") and the false alarms. We find that an SEP forecasting algorithm based solely upon the occurrence of a well-connected fast CME produces fewer false alarms (28.8%) than one based solely upon a well-connected X-class flare (50.6%). Both algorithms fail to forecast a relatively high percentage of SEP events (53.2% and 50.6% respectively). Our analysis of the historical data sets shows that false alarm X-class flares were either not associated with any CME, or were associated with a CME slower than 500 km/s; false alarm fast CMEs tended to be associated with flares of class less than M3. A better approach to forecasting would be an algorithm which takes as its base the occurrence of both CMEs and flares. We define a new forecasting algorithm which uses a combination of CME and flare parameters and show that the false alarm ratio is similar to that for the algorithm based upon fast CMEs (29.6%), but the percentage of SEP events not forecast is reduced to 32.4%. Lists of the solar events which gave rise to >40 MeV protons and the false alarms have been derived and are made available to aid further study., Comment: Solar Physics, in press

Solar Energetic Particles (SEPs) are an important component of Space Weather, including radiation hazard to humans and electronic equipment, and the ionisation of the Earth's atmosphere. We review the key observations of SEPs, our current understanding of their acceleration and transport, and discuss how this knowledge is incorporated within Space Weather forecasting tools. Mechanisms for acceleration during solar flares and at shocks driven by Coronal Mass Ejections are discussed, as well as the timing relationships between signatures of solar eruptive events and the detection of SEPs in interplanetary space. Evidence on how the parameters of SEP events are related to those of the parent solar activity is reviewed and transport effects influencing SEP propagation to near-Earth locations are examined. Finally, the approaches to forecasting Space Weather SEP effects are discussed. We conclude that both flare and CME shock acceleration contribute to Space Weather relevant SEP populations and need to be considered within forecasting tools., Comment: accepted for publication

Solar Energetic Particles (SEPs), a major component of space weather, propagate through the interplanetary medium strongly guided by the Interplanetary Magnetic Field (IMF). In this work, we analyse the implications a flat Heliospheric Current Sheet (HCS) has on proton propagation from SEP release sites to the Earth. We simulate proton propagation by integrating fully 3-D trajectories near an analytically defined flat current sheet, collecting comprehensive statistics into histograms, fluence maps and virtual observer time profiles within an energy range of 1--800 MeV. We show that protons experience significant current sheet drift to distant longitudes, causing time profiles to exhibit multiple components, which are a potential source of confusing interpretation of observations. We find that variation of current sheet thickness within a realistic parameter range has little effect on particle propagation. We show that IMF configuration strongly affects deceleration of protons. We show that in our model, the presence of a flat equatorial HCS in the inner heliosphere limits the crossing of protons into the opposite hemisphere., Comment: 16 pages, 15 figures, accepted for publication in The Astrophysical Journal

In the full-orbit particle simulations of energetic particle transport in plasmas, the plasma turbulence is typically described as a homogeneous superposition of linear Fourier modes. The turbulence evolution is, however, typically a nonlinear process, and, particularly in the heliospheric context, the solar wind plasma is inhomogeneous due to the transient structures, as observed by remote and in-situ measurements. In this work, we study the effects of the inhomogeneities on energetic particle transport by using spatially distributed, superposed turbulence envelopes. We find that the cross-field transport is significantly reduced, when compared to the results obtained with homogeneous turbulence. The reduction can reach an order of magnitude when the enveloping breaks the wave phase coherence along the mean magnetic field direction., Comment: 7 pages, 6 figures. Accepted for publication in ApJ

The trajectories of Solar Energetic Particles (SEPs) in an Interplanetary Magnetic Field (IMF) exhibiting large-scale fluctuations due to footpoint motions originating in the photosphere, are simulated using a full-orbit test-particle code. The cross-field transport experienced by the particles in three propagation conditions (scatter-free, with scattering mean free path lambda=0.3 AU and lambda=2 AU) is characterized in the Parker spiral geometry. The role of expansion of the magnetic field with radial distance from the Sun is taken into consideration in the calculation of particle displacements and diffusion coefficients from the output of the simulations. It is found that transport across the magnetic field is enhanced in the lambda=0.3 AU and lambda=2 AU cases, compared to the scatter-free case. Values of the ratios of perpendicular to parallel diffusion coefficients vary between 0.01 and 0.08. The ratio of latitudinal to longitudinal diffusion coefficient perpendicular to the magnetic field is typically 0.2, suggesting that transport in latitude may be less efficient., Comment: 17 pages, 11 figures

Context: The strong electric fields associated with magnetic reconnection in solar flares are a plausible mechanism to accelerate populations of high energy, non-thermal particles. One such reconnection scenario occurs at a 3D magnetic null point, where global plasma flows give rise to strong currents in the spine axis or fan plane. Aims: To understand the mechanism of charged particle energy gain in both the external drift region and the diffusion region associated with 3D magnetic reconnection. In doing so we evaluate the efficiency of resistive spine and fan models for particle acceleration, and find possible observables for each. Method: We use a full orbit test particle approach to study proton trajectories within electromagnetic fields that are exact solutions to the steady and incompressible magnetohydrodynamic equations. We study single particle trajectories and find energy spectra from many particle simulations. The scaling properties of the accelerated particles with respect to field and plasma parameters is investigated. Results: For fan reconnection, strong non-uniform electric drift streamlines can accelerate the bulk of the test particles. The highest energy gain is for particles that enter the current sheet, where an increasing "guide field" stabilises particles against ejection. The energy is only limited by the total electric potential energy difference across the fan current sheet. The spine model has both slow external electric drift speed and weak energy gain for particles reaching the current sheet. Conclusions: The electromagnetic fields of fan reconnection can accelerate protons to the high energies observed in solar flares, gaining up to 0.1 GeV for anomalous values of resistivity. However, the spine model, which gave a harder energy spectrum in the ideal case, is not an efficient accelerator after pressure constraints in the resistive model are included., Comment: 15 pages, 14 figures. Submitted to Astronomy and Astrophysics

A new algorithm for automatic detection of prominences on the solar limb in 304 A EUV images is presented, and results of its application to SOHO/EIT data discussed. The detection is based on the method of moments combined with a classifier analysis aimed at discriminating between limb prominences, active regions, and the quiet corona. This classifier analysis is based on a Support Vector Machine (SVM). Using a set of 12 moments of the radial intensity profiles, the algorithm performs well in discriminating between the above three categories of limb structures, with a misclassification rate of 7%. Pixels detected as belonging to a prominence are then used as starting point to reconstruct the whole prominence by morphological image processing techniques. It is planned that a catalogue of limb prominences identified in SOHO and STEREO data using this method will be made publicly available to the scientific community., Comment: 13 pages, 6 figures (4 coloured)

The distributions of sunspot longitude at first appearance and at disappearance display an east-west asymmetry that results from a reduction in visibility as one moves from disk centre to the limb. To first order, this is explicable in terms of simple geometrical foreshortening. However, the centre-to-limb visibility variation is much larger than that predicted by foreshortening. Sunspot visibility is also known to be affected by the Wilson effect: the apparent dish shape of the sunspot photosphere caused by the temperature-dependent variation of the geometrical position of the tau=1 layer. In this article we investigate the role of the Wilson effect on the sunspot appearance distributions, deducing a mean depth for the umbral tau=1 layer of 500 to 1500 km. This is based on the comparison of observations of sunspot longitude distribution and Monte Carlo simulations of sunspot appearance using different models for spot growth rate, growth time and depth of Wilson depression., Comment: 18 pages, 10 figures, in press (Solar Physics)