Context. Solar Energetic Particles (SEPs) with energy in the GeV range can propagate to Earth from their acceleration region near the Sun and produce Ground Level Enhancements (GLEs). The traditional approach to interpreting and modelling GLE observations assumes particle propagation only parallel to the magnetic field lines of interplanetary space, i.e. it is spatially 1D. Recent measurements by PAMELA have characterised SEP properties at 1 AU for the ~100 MeV-1 GeV range at high spectral resolution. Aims. We model the transport of GLE-energy solar protons through the Interplanetary Magnetic Field (IMF) using a 3D approach, to assess the effect of the Heliospheric Current Sheet (HCS) and drifts associated to the gradient and curvature of the Parker spiral. The latter are influenced by the IMF polarity. We derive 1 AU observables and compare the simulation results with data from PAMELA. Methods. We use a 3D test particle model including a HCS. Monoenergetic populations are studied first to obtain a qualitative picture of propagation patterns and numbers of crossings of the 1 AU sphere. Simulations for power law injection are used to derive intensity profiles and fluence spectra at 1 AU. A simulation for a specific event, GLE 71, is used to compare with PAMELA data. Results. Spatial patterns of 1 AU crossings and the average number of crossings are strongly influenced by 3D effects, with significant differences between periods of A+ and A- polarities. The decay time constant of 1 AU intensity profiles varies depending on the polarity and position of the observer, and it is not a simple function of the mean free path as in 1D models. Energy dependent leakage from the injection flux tube is particularly important for GLE energy particles, in many cases resulting in a roll-over in the fluence spectrum.

Solar energetic particles (SEPs), accelerated during solar eruptions, propagate in turbulent solar wind before being observed with in situ instruments. In order to interpret their origin through comparison with remote-sensing observations of the solar eruption, we thus must deconvolve the transport effects due to the turbulent magnetic fields from the SEP observations. Recent research suggests that the SEP propagation is guided by the turbulent meandering of the magnetic fieldlines across the mean magnetic field. However, the lengthening of the distance the SEPs travel, due to the fieldline meandering, has so far not been included in SEP event analysis. This omission can cause significant errors in estimation of the release times of SEPs at the Sun. We investigate the distance travelled by the SEPs by considering them to propagate along fieldlines that meander around closed magnetic islands that are inherent in turbulent plasma. We introduce a fieldline randow walk model which takes into account the physical scales associated to the magnetic islands. Our method remedies the problem of the diffusion equation resulting in unrealistically short pathlengths, and the fractal dependence of the pathlength of random walk on the length of the random-walk step. We find that the pathlength from the Sun to 1 au can be below the nominal Parker spiral length for SEP events taking place at solar longitudes 45E to 60W, whereas the western and behind-the-limb particles can experience pathlengths longer than 2 au due to fieldline meandering.

Forecasting Solar Energetic Particle (SEP) fluence, as integrated over an SEP event, is an important element when estimating the effect of solar eruptions on humans and technology in space. Current real-time estimates are based on SEP measurements at a single location in space. However, the interplanetary magnetic field corotates with the Sun approximately 13$^\circ$ each day with respect to Earth, thus in 4 days a near-Earth spacecraft will have changed their connection about 60$^\circ$ from the original SEP source. We estimate the effect of the corotation on particle fluence using a simple particle transport model, and show that ignoring corotation can cause up to an order of magnitude error in fluence estimations, depending on the interplanetary particle transport conditions. We compare the model predictions with STEREO observations of SEP events.

Insights into the processes of Solar Energetic Particle (SEP) propagation are essential for understanding how solar eruptions affect the radiation environment of near-Earth space. SEP propagation is influenced by turbulent magnetic fields in the solar wind, resulting in stochastic transport of SEPs to Earth. Multi-spacecraft observations suggest that the cross-field propagation shapes the SEP fluxes at Earth strongly. However, modelling SEP cross-field transport as spatial diffusion has been shown to be insufficient without use of unrealistically large cross-field diffusion coefficients. Recent work has shown that the early-time propagation of energetic particles across the mean field direction in turbulent fields is not diffusive, as the particles propagating along meandering field lines. This early-time transport mode results in fast access of the particles across the mean field direction, in agreement with the SEP observations. In this work, we demonstrate the significance of turbulence strength on evolution of the SEP radiation environment near Earth. We calculate the transport parameters with a turbulence transport model, parametrised by the SEP parallel scattering mean free path at 1~AU, $\lambda_{\parallel}^{*}$, and show that the parallel and cross-field transport are connected, with conditions resulting in slow parallel transport corresponding to wider events. We find a scaling $\sigma_{\phi,\, \mathrm{max}}\propto (1/\lambda_{\parallel}^{*})^{1/4}$ for the Gaussian fitting of the longitudinal distribution of maximum intensities. The longitudes with highest intensities are shifted towards the west for strong scattering conditions. Our results emphasise the importance of understanding both the SEP transport and the interplanetary turbulence conditions for modelling and predicting the SEP radiation environment at Earth.

Propagation of energetic particles across the mean field direction in turbulent magnetic fields is often described as spatial diffusion. Recently, it has been suggested that initially the particles propagate systematically along meandering field lines, and only later reach the time-asymptotic diffusive cross-field propagation. In this paper, we analyse cross-field propagation of 1--100 MeV protons in composite 2D-slab turbulence superposed on a constant background magnetic field, using full-orbit particle simulations, to study the non-diffusive phase of particle propagation with a wide range of turbulence parameters. We show that the early-time non-diffusive propagation of the particles is consistent with particle propagation along turbulently meandering field lines. This results in a wide cross-field extent of the particles already at the initial arrival of particles to a given distance along the mean field direction, unlike when using spatial diffusion particle transport models. The cross-field extent of the particle distribution remains constant for up to tens of hours in turbulence environment consistent with the inner heliosphere during solar energetic particle events. Subsequently, the particles escape from their initial meandering field lines, and the particle propagation across the mean field reaches time-asymptotic diffusion. Our analysis shows that in order to understand solar energetic particle event origins, particle transport modelling must include non-diffusive particle propagation along meandering field lines.

Solar and heliospheric cosmic rays provide a unique perspective in cosmic ray research: we can observe not only the particles, but also the properties of the plasmas in which the they are accelerated and propagate, using in situ and high-resolution remote sensing instruments. The heliospheric cosmic ray observations typically require space missions, which face stern competition against planetary and astrophysics missions, and it can take up to decades from the initial concept proposal until the actual observing of the cosmic rays can commence. Therefore it is important to have continuity in the cosmic ray mission timeline. In this overview, we review the current status and the future outlook in the experimental solar and heliospheric research. We find that the current status of the available cosmic ray observations is good, but that many of the spacecraft are near the end of their feasible mission life. We describe the three missions currently being prepared for launch, and discuss the future outlook of the solar and heliospheric cosmic ray missions.

Current particle transport models describe the propagation of charged particles across the mean field direction in turbulent plasmas as diffusion. However, recent studies suggest that at short time-scales, such as soon after solar energetic particle (SEP) injection, particles remain on turbulently meandering field lines, which results in non-diffusive initial propagation across the mean magnetic field. In this work, we use a new technique to investigate how the particles are displaced from their original field lines, and quantify the parameters of the transition from field-aligned particle propagation along meandering field lines to particle diffusion across the mean magnetic field. We show that the initial decoupling of the particles from the field lines is slow, and particles remain within a Larmor radius from their initial meandering field lines for tens to hundreds of Larmor periods, for 0.1-10 MeV protons in turbulence conditions typical of the solar wind at 1~AU. Subsequently, particles decouple from their initial field lines and after hundreds to thousands of Larmor periods reach time-asymptotic diffusive behaviour consistent with particle diffusion across the mean field caused by the meandering of the field lines. We show that the typical duration of the pre-diffusive phase, hours to tens of hours for 10 MeV protons in 1~AU solar wind turbulence conditions, is significant for SEP propagation to 1~AU and must be taken into account when modelling SEP propagation in the interplanetary space.

Context. The intensity profiles of iron and oxygen in Solar Energetic Particle (SEP) events often display differences that result in a decreasing Fe/O ratio over time. The physical mechanisms behind this behaviour are not fully understood, but these observational signatures provide important tests of physical modelling efforts. Aims. In this paper we study the propagation of iron and oxygen SEP ions using a 3D model of propagation which includes the effect of guiding centre drift in a Parker spiral magnetic field. We derive time intensity profiles for a variety of observer locations and study the temporal evolution of the Fe/O ratio. Methods. We use a 3D full orbit test particle model which includes scattering. The configuration of the interplanetary magnetic field is a unipolar Parker spiral. Particles are released instantaneously from a compact region at 2 solar radii and allowed to propagate in 3D. Results. Both Fe and O experience significant transport across the magnetic field due to gradient and curvature drifts. We find that Fe ions drift more than O ions due to their larger mass-to-charge ratio, so that an observer that is not magnetically well connected to the source region will observe Fe arriving before O, for particles within the same range in energy per nucleon. As a result, for the majority of observer locations, the Fe/O ratio displays a decrease in time. Conclusions. We conclude that propagation effects associated with drifts produce a decay over time of the Fe/O ratio, qualitatively reproducing that observed in SEP event profiles.

Context. Current solar energetic particle (SEP) propagation models describe the effects of interplanetary plasma turbulence on SEPs as diffusion, using a Fokker-Planck (FP) equation. However, FP models cannot explain the observed fast access of SEPs across the average magnetic field to regions that are widely separated in longitude within the heliosphere without using unrealistically strong cross-field diffusion. Aims. We study whether the recently suggested early non-diffusive phase of SEP propagation can explain the wide SEP events with realistic particle transport parameters. Methods. We used a novel model that accounts for the SEP propagation along field lines that meander as a result of plasma turbulence. Such a non-diffusive propagation mode has been shown to dominate the SEP cross-field propagation early in the SEP event history. We compare the new model to the traditional approach, and to SEP observations. Results. Using the new model, we reproduce the observed longitudinal extent of SEP peak fluxes that are characterised by a Gaussian profile with $\sigma=30-50^\circ$, while current diffusion theory can only explain extents of 11$^\circ$ with realistic diffusion coefficients. Our model also reproduces the timing of SEP arrival at distant longitudes, which cannot be explained using the diffusion model. Conclusions. The early onset of SEPs over a wide range of longitudes can be understood as a result of the effects of magnetic field-line random walk in the interplanetary medium and requires an SEP transport model that properly describes the non-diffusive early phase of SEP cross-field propagation.

We investigate the deceleration of Solar Energetic Particles (SEPs) during their propagation from the Sun through interplanetary space, in the presence of weak to strong scattering in a Parker spiral configuration, using relativistic full orbit test particle simulations. The calculations retain all three spatial variables describing particles' trajectories, allowing to model any transport across the magnetic field. Large energy change is shown to occur for protons, due to the combined effect of standard adiabatic deceleration and a significant contribution from particle drift in the direction opposite to that of the solar wind electric field. The latter drift-induced deceleration is found to have a stronger effect for SEP energies than for galactic cosmic rays. The kinetic energy of protons injected at 1 MeV is found to be reduced by between 35 and 90% after four days, and for protons injected at 100 MeV by between 20 and 55%. The overall degree of deceleration is a weak function of the scattering mean free path, showing that, although adiabatic deceleration plays a role, a large contribution is due to particle drift. Current SEP transport models are found to account for drift-induced deceleration in an approximate way and their accuracy will need to be assessed in future work.

The capability to predict the parameters of an SEP event such as its onset, peak flux, and duration is critical to assessing any potential space weather impact. We present a new flexible modelling system simulating the propagation of Solar Energetic Particles (SEPs) from locations near the Sun to any given location in the heliosphere to forecast the SEP flux profiles. SPARX uses an innovative methodology, that allows implementation within an operational framework, to overcome the time constraints of test particle modelling of SEP profiles, allowing the production of near real time SEP nowcasts and forecasts, when paired with appropriate near real time triggers. SPARX has the capability to produce SEP forecasts within minutes of being triggered by observations of a solar eruptive event. The model is based on the test particle approach and is spatially 3D, thus allowing for the possibility of transport in the direction perpendicular to the magnetic field. The model naturally includes the effects of perpendicular propagation due to drifts and drift-induced deceleration. The modelling framework and the way in which parameters of relevance for Space Weather forecasting are obtained are described. The first results from the modelling system are presented. These results demonstrate that corotation and drift of SEP streams play an important role in shaping SEP flux profiles.

To understand the origin of Solar Energetic Particles (SEPs), we must study their injection time relative to other solar eruption manifestations. Traditionally the injection time is determined using the Velocity Dispersion Analysis (VDA) where a linear fit of the observed event onset times at 1 AU to the inverse velocities of SEPs is used to derive the injection time and path length of the first-arriving particles. VDA does not, however, take into account that the particles that produce a statistically observable onset at 1 AU have scattered in the interplanetary space. We use Monte Carlo test particle simulations of energetic protons to study the effect of particle scattering on the observable SEP event onset above pre-event background, and consequently on VDA results. We find that the VDA results are sensitive to the properties of the pre-event and event particle spectra as well as SEP injection and scattering parameters. In particular, a VDA-obtained path length that is close to the nominal Parker spiral length does not imply that the VDA injection time is correct. We study the delay to the observed onset caused by scattering of the particles and derive a simple estimate for the delay time by using the rate of intensity increase at the SEP onset as a parameter. We apply the correction to a magnetically well-connected SEP event of June 10 2000, and show it to improve both the path length and injection time estimates, while also increasing the error limits to better reflect the inherent uncertainties of VDA.

Context: Diffusive shock acceleration in the solar corona can accelerate solar energetic particles to very high energies. Acceleration efficiency is increased by entrapment through self-generated waves, which is highly dependent on the amount of accelerated particles. This, in turn, is determined by the efficiency of particle injection into the acceleration process. Aims: We present an analysis of the injection efficiency at coronal shocks of varying obliquity. We assessed injection through reflection and downstream scattering, including the effect of a cross-shock potential. Both quasi-thermal and suprathermal seed populations were analysed. We present results on the effect of cross-field diffusion downstream of the shock on the injection efficiency. Methods: Using analytical methods, we present applicable injection speed thresholds that were compared with both semi-analytical flux integration and Monte Carlo simulations, which do not resort to binary thresholds. Shock-normal angle and shock-normal velocity were varied to assess the injection efficiency with respect to these parameters. Results: We present evidence of a significant bias of thermal seed particle injection at small shock-normal angles. We show that downstream isotropisation methods affect the theta-dependence of this result. We show a non-negligible effect caused by the cross-shock potential, and that the effect of downstream cross-field diffusion is highly dependent on boundary definitions. Conclusions: Our results show that for Monte Carlo simulations of coronal shock acceleration a full distribution function assessment with downstream isotropisation through scatterings is necessary to realistically model particle injection. Based on our results, seed particle injection at quasi-parallel coronal shocks can result in significant acceleration efficiency, especially when combined with varying field-line geometry.

Drifts are known to play a role in galactic cosmic ray transport within the heliosphere and are a standard component of cosmic ray propagation models. However, the current paradigm of Solar Energetic Particle (SEP) propagation holds the effects of drifts to be negligible, and they are not accounted for in most current SEP modelling efforts. We present full-orbit test particle simulations of SEP propagation in a Parker spiral interplanetary magnetic field which demonstrate that high energy particle drifts cause significant asymmetric propagation perpendicular to the interplanetary magnetic field. Thus in many cases the assumption of field aligned propagation of SEPs may not be valid. We show that SEP drifts have dependencies on energy, heliographic latitude, and charge to mass ratio, that are capable of transporting energetic particles perpendicular to the field over significant distances within interplanetary space, e.g. protons of initial energy 100 MeV propagate distances across the field on the order of 1 AU, over timescales typical of a gradual SEP event. Our results demonstrate the need for current models of SEP events to include the effects of particle drift. We show that the drift is considerably stronger for heavy ion SEPs due to their larger mass to charge ratio. This paradigm shift has important consequences for the modelling of SEP events and is crucial to the understanding and interpretation of in-situ observations.

Solar energetic particles (SEPs) have been observed to easily spread across heliographic longitudes, and the mechanisms responsible for this behaviour remain unclear. We use full-orbit simulations of a 10 MeV proton beam in a turbulent magnetic field to study to what extent the spread across the mean field can be described as diffusion early in a particle event. We compare the full-orbit code results to solutions of a Fokker-Planck equation including spatial and pitch angle diffusion, and of one including also propagation of the particles along random-walking magnetic field lines. We find that propagation of the particles along meandering field lines is the key process determining their cross-field spread at 1 AU at the beginning of the simulated event. The mean square displacement of the particles an hour after injection is an order of magnitude larger than that given by the diffusion model, indicating that models employing spatial cross-field diffusion cannot be used to describe early evolution of an SEP event. On the other hand, the diffusion of the particles from their initial field lines is negligible during the first 5 hours, which is consistent with the observations of SEP intensity dropouts. We conclude that modelling SEP events must take into account the particle propagation along meandering field lines for the first 20 hours of the event.

Drifts in the Parker spiral interplanetary magnetic field are known to be an important component in the propagation of galactic cosmic rays, while they are thought to be negligible for Solar Energetic Particles (SEPs). As a result they have so far been ignored in SEP propagation modelling and data analysis. We examine drift velocities in the Parker spiral within single particle first-order adiabatic theory, in a local coordinate system with an axis parallel to the magnetic field. We show that, in the presence of scattering in interplanetary space, protons at the high end of the SEP energy range experience significant gradient and curvature drift. In the scatter-free case, drift due to magnetic field curvature is present. The magnitude of drift velocity increases by more than an order of magnitude at high heliographic latitudes compared to near the ecliptic; it has a strong dependence on radial distance r from the Sun, reaching a maximum at r~1 AU at low heliolatitudes and r~10 AU at high heliolatitudes. Due to the mass over charge dependence of drift velocities, the effect of drift for partially ionised SEP heavy ions is stronger than for protons. Drift is therefore likely to be a considerable source of cross field transport for high energy SEPs.

Observations and modelling suggest that the fluctuations in magnetised plasmas exhibit scale-dependent anisotropy, with more energy in the fluctuations perpendicular to the mean magnetic field than in the parallel fluctuations and the anisotropy increasing at smaller scales. The scale-dependence of the anisotropy has not been studied in full-orbit simulations of particle transport in turbulent plasmas so far. In this paper, we construct a model of critically balanced turbulence, as suggested by \cite{GoSr1995}, and calculate energetic particle spatial diffusion coefficients using full-orbit simulations. The model uses an enveloped turbulence approach, where each 2-dimensional wave mode with wavenumber $k_\perp$ is packed into envelopes of length $L$ following the critical balance condition, $L\propto k_\perp^{-2/3}$, with the wave mode parameters changing between envelopes. Using full-orbit particle simulations, we find that both the parallel and perpendicular diffusion coefficients increase by a factor 2, compared to previous models with scale-independent anisotropy.

High-energy particles stream during coronal mass ejections or flares through the plasma of the solar wind. This causes instabilities, which lead to wave growth at specific resonant wave numbers, especially within shock regions. These amplified wave modes influence the turbulent scattering process significantly. In this paper, results of particle transport and scattering in turbulent plasmas with excited wave modes are presented. The method used is a hybrid simulation code, which treats the heliospheric turbulence by an incompressible magnetohydrodynamic approach separately from a kinetic particle description. Furthermore, a semi-analytical model using quasilinear theory (QLT) is compared to the numerical results. This paper aims at a more fundamental understanding and interpretation of the pitch-angle scattering coefficients. Our calculations show a good agreement of particle simulations and the QLT for broad-band turbulent spectra; for higher turbulence levels and particle beam driven plasmas, the QLT approximation gets worse. Especially the resonance gap at $\mu=0$ poses a well-known problem for QLT for steep turbulence spectra, whereas test-particle computations show no problems for the particles to scatter across this region. The reason is that the sharp resonant wave--particle interactions in QLT are an oversimplification of the broader resonances in test-particle calculations, which result from nonlinear effects not included in the QLT. We emphasise the importance of these results for both numerical simulations and analytical particle transport approaches, especially the validity of the QLT.

We have developed a simulation model of particle acceleration in coronal shock waves. The model is based on a Monte Carlo method, where particles are traced in prescribed large-scale electromagnetic fields utilizing the guiding center approximation. The particles are scattered in the turbulence according to quasilinear theory, with the scattering amplitude directly proportional to the intensity of Alfv\'en waves at gyro-resonant wavenumbers. The Alfv\'en waves are traced simultaneously with the particles, so that the wave field is propagated outwards from the Sun using WKB propagation supplemented with a phenomenological wavenumber diffusion term and a growth rate computed from the net flux of the accelerated particles. We consider initial wave amplitudes small enough to allow rapid escape of particles from the shock to the ambient medium. Thus, in our model the Alfv\'en waves responsible for the diffusive acceleration of particles are generated by the accelerated particles themselves. In this work, we study the effects of non-constant shock velocity and non-monotonic Alfv\'en velocity on particle acceleration scenarios. We report in particular how the deceleration of a shock affects particle intensity and turbulence power evolution in the vicinity of the shock.

Context: Acceleration in coronal mass ejection driven shocks is currently considered the primary source of large solar energetic particle events. Aims: The solar wind, which feeds shock-accelerated particles, includes numerous ion populations, which offer much insight into acceleration processes. We present first simulations of shock-accelerated minor ions, in order to explore trapping dynamics and acceleration timescales in detail. Methods: We have simulated diffusive shock acceleration of minor ions (3He2+, 4He2+, 16O6+ and 56Fe14+) and protons using a Monte Carlo method, where self-generated Alfv\'enic turbulence allows for repeated shock crossings and acceleration to high energies. Results: We present the effect of minor ions on wave generation, especially at low wavenumbers, and show that it is significant. We find that maximum ion energy is determined by the competing effects of particle escape due to focusing in an expanding flux tube and trapping due to the amplified turbulence. We show the dependence of cut-off energy on the particle charge to mass ratio to be approximately (Q/A)^1.5. Conclusions: We suggest that understanding the acceleration of minor ions at coronal shocks requires simulations which allow us to explore trapping dynamics and acceleration timescales in detail, including evolution of the turbulent trapping boundary. We conclude that steady-state models do not adequately describe the acceleration of heavy ions in coronal shocks.

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.

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.