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.

A statistical analysis of the relationship between solar energetic particles (SEPs) and properties of solar flares and coronal mass ejections (CMEs) is presented. SEP events during solar cycle 23 are selected which are associated with solar flares originating on the visible hemisphere of the Sun and at least of magnitude M1. Taking into account all flares and CMEs that occurred during this period, the probability for the occurrence of an SEP event near Earth is determined. A strong rise of this probability is observed for increasing flare intensities, more western locations, larger CME speeds and halo CMEs. The correlations between the proton peak flux and these solar parameters are derived for a low (>10 MeV) and high (>60 MeV) energy range excluding any flux enhancement due to the passage of fast interplanetary shocks. The obtained correlation coefficients are: 0.55+-0.07 (0.63+-0.06) with flare intensity and 0.56+-0.08 (0.40+-0.09) with the CME speed for E>10 MeV (E>60 MeV). For both energy ranges, the correlations with flare longitude and CME width are very small or non-existent. Furthermore, the occurrence probabilities, correlation coefficients and mean peak fluxes are derived in multi-dimensional bins combining the aforementioned solar parameters. The correlation coefficients are also determined in different proton energy channels ranging from 5 to 200 MeV. The results show that the correlation between the proton peak flux and the CME speed decreases with energy, while the correlation with the flare intensity shows the opposite behavior. Furthermore, the correlation with the CME speed is stronger than the correlation with the flare intensity below 15 MeV and becomes weaker above 20 MeV. Excluding the flux enhancements due to interplanetary shocks, only a small but not very significant change is observed in the correlation between the peak flux below 7 MeV and the CME speed.

Shortly after the seminal paper {\sl "Self-Organized Criticality: An explanation of 1/f noise"} by Bak, Tang, and Wiesenfeld (1987), the idea has been applied to solar physics, in {\sl "Avalanches and the Distribution of Solar Flares"} by Lu and Hamilton (1991). In the following years, an inspiring cross-fertilization from complexity theory to solar and astrophysics took place, where the SOC concept was initially applied to solar flares, stellar flares, and magnetospheric substorms, and later extended to the radiation belt, the heliosphere, lunar craters, the asteroid belt, the Saturn ring, pulsar glitches, soft X-ray repeaters, blazars, black-hole objects, cosmic rays, and boson clouds. The application of SOC concepts has been performed by numerical cellular automaton simulations, by analytical calculations of statistical (powerlaw-like) distributions based on physical scaling laws, and by observational tests of theoretically predicted size distributions and waiting time distributions. Attempts have been undertaken to import physical models into the numerical SOC toy models, such as the discretization of magneto-hydrodynamics (MHD) processes. The novel applications stimulated also vigorous debates about the discrimination between SOC models, SOC-like, and non-SOC processes, such as phase transitions, turbulence, random-walk diffusion, percolation, branching processes, network theory, chaos theory, fractality, multi-scale, and other complexity phenomena. We review SOC studies from the last 25 years and highlight new trends, open questions, and future challenges, as discussed during two recent ISSI workshops on this theme.