Issue 
J. Space Weather Space Clim.
Volume 8, 2018
Measurement, Specification and Forecasting of the Solar Energetic Particle Environment and GLEs



Article Number  A06  
Number of page(s)  10  
DOI  https://doi.org/10.1051/swsc/2017040  
Published online  06 February 2018 
Research Article
The virtual enhancements − solar proton event radiation (VESPER) model
^{1}
Institute of Accelerating Systems and Applications,
P.O. Box 17214,
GR10024
Athens, Greece
^{2}
Space Applications & Research Consultancy (SPARC),
Kleissovis 8,
GR10677
Athens, Greece
^{3}
Department of Physics, National and Kapodistrian University of Athens,
Athens, Greece
^{4}
European Research and Technology Centre, European Space Agency,
Noordwijk, Netherlands
^{*} sagiamini@sparc.gr
Received:
22
June
2017
Accepted:
4
December
2017
A new probabilistic model introducing a novel paradigm for the modelling of the solar proton environment at 1 AU is presented. The virtual enhancements − solar proton event radiation model (VESPER) uses the European space agency's solar energetic particle environment modelling (SEPEM) Reference Dataset and produces virtual timeseries of proton differential fluxes. In this regard it fundamentally diverges from the approach of existing SPE models that are based on probabilistic descriptions of SPE macroscopic characteristics such as peak flux and cumulative fluence. It is shown that VESPER reproduces well the dataset characteristics it uses, and further comparisons with existing models are made with respect to their results. The production of timeseries as the main output of the model opens a straightforward way for the calculation of solar proton radiation effects in terms of timeseries and the pairing with effects caused by trapped radiation and galactic cosmic rays.
Key words: Radiation environment / SEP / modelling / radiation / missions
© S. AminalragiaGiamini et al., Published by EDP Sciences 2018
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.
1 Introduction
Solar particle events (SPE) are of particular interest for spacecraft design and radiation estimations. Protons with energies higher than ∼1 MeV have the ability to cause typical semiconductor (Si, GaAs, etc.) devices to malfunction through the effects of total ionizing dose (TID) such as errors caused by Single Event Upsets (SEUs) (Dodd & Massengill, 2003) which cause the excitation of electrons triggering false signals. Furthermore lasting damage can be caused by protons, and heavier nuclei, which result in atomic displacements in the crystalline lattices. Such defects can seriously hamper device function or even cause failure. Underestimation of such risks can degrade the performance of systems aboard spacecrafts and cause important loss of mission lifetime. Such considerations become even more important for manned missions, especially outside the protective shielding of Earth's magnetic field especially with regard to manned missions to the Moon or Mars (Cucinotta et al., 2010). On the other hand, overprotection against these risks in the form of excessive and unneeded shielding, results in increased cost, design complexity and severely reduced payloads.
In order to have reliable and costeffective spacecraft design and implement new space technologies, accurate and reliable models for estimating the various radiation risks (Crosby et al., 2008) have been developed over the last years. Specifically for SPEs, evidence has been shown by Xapsos et al. (2006) that they are a SelfOrganised Critical (SOC) phenomenon. This suggests that it may not be possible to predict the size/intensity or the occurrence of SPEs over long periods. It underscores that along with important work in the theoretical understanding (Vainio et al., 2009), SPE modelling efforts will continue to use techniques from the field of probability theory. Statistical, probabilistic models provide predictions of the cumulative and maximum fluence or peak flux values to take place over the course of a mission as a function of confidence level over a range of energies seen at the Earth. The selection of the physical variable of the SPE to be modelled (fluence/peak flux), the data that are being used, and the selection of suitable mathematical functions that provide satisfactory matching with the measured proton flux/fluence distributions have resulted in the development of various probabilistic solar proton models, King (Stassinopoulos & King, 1974), JPL (Feynman et al., 1993), ESP (Xapsos et al., 1998, 1999), Nymmik (Nymmik, 1999). These widely used models are usually based on a Poisson description regarding the event frequency distribution and they usually assume lognormal distributions and/or power function distributions for the modelled variable(s). The lognormal distribution describes the large events well but underestimates the probability of smaller events. On the other hand, the powerlaw functions describe the smaller events well but overestimate the probability of larger events. Jiggens & Gabriel (2009) have questioned that the occurrence of SPEs is really random. The authors chose alternative descriptions (based on Lèvy or on timedependent Poisson distributions) and showed that the resulting probability models do provide better fits for the available event lists. The success on using nonPoisson descriptions suggests that there is an inherent memory in the system and the use of Poissonian description in SPE models may not be the optimum selection. Later on, Jiggens et al. (2012) developed a new modelling methodology − the “Virtual Timelines” method (VTM). The novel idea behind this method was to allow for the incorporation of the duration of SPEs which are nonpointlike in nature but this cannot be explicitly considered when sampling the frequency distribution. Using VTM, it was shown that both waiting times and event duration distributions are fitted better when using nonPoisson distributions. In this work we present a probabilistic approach which is based on the modelling of solar proton flux timeseries. The virtual enhancements − solar proton event radiation (VESPER) model utilizes the European space agency's solar energetic particle environment modelling (SEPEM) (Crosby et al., 2015) Reference Dataset v2.0 (RDSv2) (Heynderickx et al., 2017) and employs a datadriven approach to produce Virtual TimeSeries (VTS) of differential proton fluxes. In turn, the statistical analysis of the VESPER outputs, the timeseries, provides distribution functions used in the probabilistic estimation of quantities such as the Cumulative Fluence and Peak Flux over the duration of a space mission. The fact that the model produces full timeseries to describe solar proton flux enhancements presents unique advantages, since such an output can be used to produce radiation effects timeseries (e.g. ionizing or nonionizing dose) as well as to be combined with similar outputs of radiation belt models, such as the AE9/AP9 (Ginet et al., 2013) models, that also produce timeseries outputs of either fluxes/fluences or radiation effects. Furthermore, by applying a Magnetospheric Shielding model one can derive proton cutoff rigidity values for a userdefined orbit that crosses earth's magnetosphere, and map the output VTS along this orbit. This can provide a full description of the SPE environment in terms of proton flux timeseries. Thus, we are confident that the VESPER model will prove valuable to the scientists and engineers that design future missions. In what follows, we present the basic premises of the VESPER model, the way it is constructed and we evaluate its output and compare it with those of the JPL, ESP and VTM models. It is shown that characteristic statistical properties of VESPER outputs (timeseries) are consistent with those of the RDSv2 database it uses and that there is good agreement with the aforementioned models, especially with the VTM, despite the fundamental differences in their approach.
2 Data
2.1 SEPEM reference dataset
We have used the SEPEM RDS v2.0 dataset (RDSv2) comprised of proton differential flux timeseries. The dataset spans more than 40 years, 19740701–20151231, and contains processed data from the NOAA energetic particles sensor (EPS), part of the space environment monitor (SEM) package onboard GOES and earlier SMS satellites. These data have been crosscalibrated to find the effective energy using data from the goddard medium energy (GME) instrument onboard the IMP8 spacecraft. A description of the crosscalibration is available in Sandberg et al. (2014). Finally, the spectra for each time stamp were rebinned into a set of 11 logarithmically spaced energy bins. In this work we have used the first 10 energy bins shown in Table 1. The rebinned data are continuous with 5 min time intervals without time gaps or missing values with ∼4 × 10^{5} measurements in total. It is worth mentioning that the RDSv2 integral fluxes have been recently validated using STEREOA/HET and STEREOB/HET as well as IMP8/GME integral proton flux data by Rodriguez et al. (2017).
For the purposes of the model construction, we have introduced the new variable LF_{2} which can be defined as the 2nd moment of the flux spectrum along the energy axis in logarithmic space. (1) Here, f(E) [cm^{−2} MeV^{−1} s^{−1} sr^{−1}] denotes the differential proton flux intensity and E [MeV] the proton energy. The minimum and maximum energies used here are those of the first and tenth channel of RDSv2 (see Tab. 1). The calculation of the LF_{2} series for the entirety of the dataset results in the creation of a single timeseries that characterizes RDSv2. The LF_{2} variable is defined to serve the goals of the model while still being firmly based solely on the flux timeseries. As it will be shown, it plays a key role in the construction of VESPER, and it is also used for the event definition process.
The first 10 energy bins and logmean energies of the RDSv2 proton flux dataset.
2.2 SPE definition
In this study, we define as Solar Particle Event the enhancements of the LF_{2} values above their timevarying background level that last for at least 12 h and peaked at a value at least two times larger than the temporary background level. In order to apply this SPE definition, we have employed a novel scheme to determine the time varying background level which is governed by the galactic cosmic ray (GCR) modulation. A “clump” is defined as the group of at least 3 consecutive measurements above a test threshold. For a given timewindow, the background is defined as the threshold level for which the number of “clumps” is maximized. This scheme produces a time varying threshold which follows the GCR modulation and allows an accurate identification of the flux enhancements. Since this extraction of the SPEs is applied on the LF_{2} timeseries, the proton flux enhancements across all energies are taken into account.
In Figure 1, we present the LF_{2} timeseries, the calculated background level and the extracted SPE flux enhancements. Using this event definition, we define 241 SPEs and we calculate their respective:

duration times;

(sum of the LF_{2} values over the event duration);

waittimes.
We note that we use the convention used by SEPEM where the duration of an event is the same for all energy channels. This way all flux timeseries have common start and end time points. The WaitTimes are defined as the time intervals between the end of an Event and the start of the next one.
Figure 2 shows the differential proton flux series at the ten different energies and the calculated LF_{2} for two SPEs. Here, the way the LF_{2} variable takes into account the flux enhancements across all energy ranges is clear. For Event B, the flux of the 1st energy channel at the beginning of day 2 (Fig. 2b − red arrow) is two orders of magnitude higher than the flux of Event A at its onset (Fig. 2a − red arrow). However, whereas at the onset of Event A contributions from all energy channels are strong, for Event B the higher energy channels barely rise above background level. This results in the calculated LF_{2} values at these time points being comparable precisely due to the spectral profiles of the flux timeseries.
Fig. 1
The LF_{2} timeseries calculated for the entire SEPEM reference dataset v2.0. The background is denoted by the green line and the flux enhancements denoted by red. 
Fig. 2
The timeseries for two defined SPEs. Event A: 16Feb1984 to 24Feb1984, Event B: 11Jul1982 to 14Jul1982. (a) and (b) flux timeseries of the ten energy channels, (c) and (d) timeseries of the LF_{2} variable calculated from the fluxes. 
3 Structure of the VESPER model
The VESPER model implements a probabilistic approach to model the occurrence and characteristics of SPEs with the goal of creating virtual timeseries of differential proton fluxes. This is achieved by creating Virtual Events (VEs) from measured Events extracted from the database, which we dub “seeds”. The profiles/timeseries of these measured Events are rescaled in time and flux which essentially amounts to “stretching” or “compressing” timeseries in their two dimensions, time and intensity, with the goal of creating VEs with new profiles which are however based on existing ones. The rescaling of the SPE flux series is performed in the time and intensity scales in such a manner that it preserves the statistical properties of the SPEs (macroscopic) variables (such as the peak flux and fluence distributions) as well as the inherent spectral properties of the flux series during SPEs. The timeline of the SPE VTS depends on the fitting of the Duration and the WaitTime cumulative distribution functions (CDFs) while the creation of the virtual SPEs is based on the observed dependence between the SLF_{2} and the Duration of the extracted events (see Fig. 4a). The process by which SPE VTS are constructed consists of three main steps:

random sampling of the analytical fits on the CDFs to produce Virtual WaitTimes and Duration values;

calculation of the Virtual SLF_{2} using the SLF_{2}Duration linear dependence in the loglog space and “seed” event selection;

scaling of the seed Event to the selected Duration and SLF_{2} values to create a VE.
Each step is explained and described in detail in the following paragraphs.
Fig. 4
Cross plots of Duration versus SLF_{2}. The red points correspond to the calculated values based on the analysis of RDSv2. (a) The dashed lines denote the vicinity of the VE which is used to calculate an SLF_{2} value, adjust the Duration value and select a seed for the Virtual Event. (b) Black dots correspond to calculated values for the VEs. 
3.1 Sampling of duration and waittime cumulative distribution functions
The fitting and the random sampling of the WaitTimes and Duration CDFs of the extracted SPEs determines the occurrence of the VEs within the examined timespan of a mission duration, e.g. 10 years. The occurrence frequency of SPEs is strongly dependent with solar cycle phases, i.e. active or quiet, thus for the modelling of WaitTimes between SPEs we have employed separate descriptions for each solar phase and derived two distinct distribution functions. However, since no correlation has been demonstrated between Event duration or intensity and solar phases the SLF_{2} and Duration values are taken from the list of defined SPEs from all years of the RDSv2 regardless of whether the period was a solar active or quiet one. The solar cycles were split into 7 active/maximum and the remaining (∼4) quiet/minimum years, inline with the work of Feynman et al. (1990). For the solar active phases, there were 2.5 years before the peak times and 4.5 years after. The peak times for cycles 20–22 are taken from Xapsos et al. (1998) as 1968.9, 1979.9 and 1989.9, for cycle 23 it was taken as 2000.6 in order to ensure the inclusion of January 2005 in the solar maximum period. For the last cycle the situation cannot be absolutely clear until the active period is complete, however 2009 should still be solar minimum and it's likely that 2011–2017 will be solar maximum. In Figure 3, the experimental and the fitted CDFs are presented. For all variables, we have employed exponential cutoff powerlaw functions which provide a good fit/approximation of the measured data. At 95% confidence bounds, for the Durations fit the RMSE value is equal to 0.081, and for the WaitTimes in solar active and solar quiet phases equal to 0.128 and 0.142, respectively. It is evident, that for a given mission length, e.g. 7 years, the sum of the sampled Durations of the Virtual Events and their respective WaitTimes must be equal to that mission length. As a consequence, the sampled WaitTimes and Durations determine the total number of VEs that will occur within a given mission duration. WaitTimes have quite higher values than Durations and therefore play a dominant role since their underestimation or overestimation can heavily affect the output of the model, especially the resulting Fluence, for a given mission duration.
Fig. 3
CDFs for (a) the Duration and (b) the WaitTime variables and their fits with exponential cutoff powerlaws. 
3.2 SLF_{2} calculation, Duration adjustment and seed Event selection
After the random selection of the Virtual Duration value from the fitted CDF, a Virtual SLF_{2} value is calculated and a seed event is selected. This triplet (Duration, SLF_{2} and seed event) is later on used in the construction of the VE. The SLF_{2} variable has a strong correlation with the SPE Durations which is much higher than those between Durations and SPE Fluences of any energy bin. In this step we employ the singular free parameter in our model, dubbed tolerance, for which we have used the value 0.1. Tolerance (τ) is used to calculate two factors (1 + τ, 1−τ), essentially percentages, applied to the selected log(Duration) to calculate an upper and a lower bound for the subsequent calculation of SLF_{2} and the selection of a seed event. These bounds in turn define the vicinity of a virtual event and dictate how wide or narrow it will be as seen in Figure 4a. For a randomly sampled Virtual Duration (VD) and tolerance value (τ) we determine three points, defined by (cf. Fig. 4a):

The lower bound pair: LB_{x} = (1 − τ) * ln(VD), LB_{y}.

The middle pair: ln(VD), MV_{y}.

The upper bound pair: UB_{x} = (1 + τ) * ln(VD), UB_{y}.
The lines perpendicular to the linear fit on a loglog scale that cross these three points create the bounds that define the upper and lower parts of the vicinity of the VE as seen in Figure 4a. The distances of all the Events in the vicinity from the fitting line are initially calculated, and their standard deviation is used as the σ (sigma) parameter in a Gaussian distribution sampling. The Gaussian curve is centred at the [ln(VD), MV_{y}] point and is also perpendicular to the loglinear fit. A value is randomly sampled and used as the distance of the VE from the fit. This way the SLF_{2} and the adjusted Duration values for the VE are calculated allowing for the scatter seen in the linear regression on a loglog scale. The resulting pairs for the VEs and their loglog linear fit adequately reproduce the scattering around the loglog linear fit of the actual Events (cf. Fig. 4b). Lastly, a seed Event is randomly chosen from a pair of Events which have been randomly selected one from the upper and one from the lower vicinity of the VE.
3.3 Flux scaling and construction of Virtual Event timeseries
For each VE to be created the flux timeseries of its seed Event are initially interpolated to the selected Virtual Duration creating an “intermediate” timescaled event. We underline here that since the seed events are chosen from the vicinity of the virtual Duration no extreme under or oversampling occurs that could heavily alter the profiles of the time series. Next, the (interpolated) flux timeseries are rescaled to match the virtual SLF_{2} value. The rescaling scheme essentially alters the intensity of each proton spectrum within the event while maintaining its spectral profile. This way the resulting VEs consist of spectra with energy slopes very close to measured ones avoiding the creation of unnaturally abrupt or hard/soft spectra. This process consists of the steps described below.

The SLF_{2} value of the (intermediate) timescaled event is calculated. The division with the Virtual SLF_{2} value selected from the Gaussian sampling, produces a scaling factor, sf, defined by:

The proton flux spectra are assumed to be of the general form by f (E) = aG (E) where a is the intensity of the spectrum and G(E) the intrinsic energy distribution described by some function, e.g. for powerlaw G(E) = E^{b}_{.} Therefore in the logarithmic space this is expressed as: (3)
Proton spectra can be adequately fitted with exponential or powerlaw functions, however such specific assumptions are not made here, only the very general form is assumed.

Using the general assumption for the spectrum and equation (1), the LF_{2} integral can be analytically expressed as: (4)

The rescaled timeseries is produced by multiplying the timescaled event's LF_{2}(t) values with the scaling factor sf, which can be also expressed as (5)

Using equations (4) and (5) new intensity ln(a ′) values for all spectra are directly derived from the values. This produces the new intensities of the spectra. (6)

Equation (6) provides the factor to directly alter the intensity of the spectrum. Using equation (4) the rescaled fluxes are calculated by adding the factor in logspace to the existing spectrum. (7)
As can be seen in equation (7), only the scaling factor and the LF_{2} value of the spectrum are actually required for the scaling while the use of a specific analytic function that fits the spectrum is not. This property stems from the fact that in the rescaling scheme essentially only the intensity is altered and neither a nor G(E) need be determined explicitly by fitting the spectrum. Using the process described VEs are constructed having been scaled both in time and flux in a consistent way. An example of such a rescaled event is shown in Figure 5.
The output of the process is a set of Virtual Events described by their flux timeseries and their intermediate WaitTimes. For a given mission duration a large number of iterations or “scenarios” are run and from each iteration important quantities of the encountered SEP environment − such as the Cumulative Mission Fluence, the WorstCase SPE Fluence and Peak Flux − are calculated. For a sufficient number of iterations, the values of these quantities are used to construct the resulting probability and cumulative density functions and derive further probabilistic outputs.
Fig. 5
(a) Solar proton event, 27Nov1989 to 05Dec1989 (b) Virtual event rescaled to higher Duration and SLF_{2} values. 
4 Results and discussion
4.1 Model outputs and dataset reproducibility
In what follows, we present outputs of the VESPER model for mission durations of 1 and 7 years during solar active period. From every iteration the sum of all virtual fluxes results in the mission Cumulative Fluence and the maximum value is the Highest Peak Flux. We run 1 50 000 iterations for 1 year mission and 50 000 for 7 year mission to produce the cumulative distributions shown below in Figure 6.
An interesting feature appears in the way particular statistical quantities scale from 1 to 7 year mission. In the Cumulative Fluence CDFs for proton energies of 6 MeV, significant differences − of the order of 20 − appear at the high probability of exceeding = 0.9, while the differences become much smaller − of the order of 3 − at the low probability of exceeding = 0.01. This is reasonable since at low probabilities one would expect few very energetic events to largely determine the cumulative fluence for a mission and this should not differ dramatically whether it is a 1 year or a 7 year period. Conversely, at low probabilities the cumulative fluence is determined largely by many more frequent and less energetic events which will result in very different total values between the two cases of mission duration. This is also reflected in the Peak Flux CDFs, however the differences are less dramatic for high probabilities as only single flux measurements are examined and the effect of the summation of multiple events does not come into play. Also, an interesting feature is that the VESPER outputs result in curves that are not always smooth but show varying curvatures and “knees“ as in the case of Figures 6a, 6c. This is due to the inherent statistic behavior of the RDSv2. It is therefore a feature of the data that is fed into the model and we consider it a strength of the model since it seems to reproduce the statistical behavior of the measurements in the Event list.
We also perform a comparative statistical analysis on the RDSv2 database and the output of the VESPER model. This analysis serves to verify the selfconsistency of the model in the sense that it reproduces basic statistical characteristics of the database it is based upon. Flux intensity and total cumulative fluence over a mission duration are two major considerations for any probabilistic model and thus it is important that the VESPER model reproduces characteristics of its progenitor database. We use a moving threshold approach to evaluate the output in terms of total fluence and median flux.
All measurements of the flux timeseries whose LF_{2} values lie above a threshold are identified and summed to calculate the total fluence and their median flux. The process is repeated for multiple threshold values resulting in the graphs shown in Figure 7. The total fluence is normalized per year to compare 42 years of the RDSv2 to 10 000 years produced by the VESPER. As can be seen in Figure 7, the curves defined by the moving threshold method applied to the virtual timeseries reasonably reproduce those derived by using RDSv2.
Fig. 6
(a) and (b) CDFs of cumulative differential fluence produced from the VESPER model for 1 and 7 year solar active periods, respectively. Each curve presents one of the channels in the 5–200 MeV proton energy range, (c) and (d) same for peak differential flux. 
Fig. 7
Comparison of (a) total fluence and (b) median flux values for RDSv2 and VESPER virtual data derived using a moving threshold on the respective LF_{2} timeseries. Crosses denote RDSv2 and continuous lines virtual data. 
4.2 Comparisons with existing models
We compare the output of our model to three well tested probabilistic models that are widely used for the estimation of probability functions of quantities such as the Cumulative Fluence and the Peak Flux over a mission duration. We have used here the ESP, JPL and VTM models which are available in the online ESA SEPEM application server (http://dev.sepem.oma.be/). Figure 8 shows a comparison of the output of the three models with the resulting Cumulative Fluence and Peak Flux distributions from VESPER using 50 000 scenarios for a mission duration of 7 year in solar active period.
The comparison of the models' outputs demonstrates that VESPER's virtual timeseries result in Fluence and Peak values whose cumulative distributions are of the same broader magnitude with those produced by existing models. VESPER is within an order of magnitude with the ESP model and shows good agreement with the VTM model while JPL results are typically higher as one would expect from this model. Especially for the VTM model we note that for lowprobabilities/high values the two models show a remarkable degree of agreement. For high probabilities VESPER produces typically somewhat lower values, however the differences are not dramatic remaining within a factor of 2.
In the VESPER model we have employed the “virtual timeline” approach introduced by the VTM model, where a duration and a waittime distribution are sampled to determine the occurrence of Virtual Events. However, despite this common element, the two models differ fundamentally in their respective approaches and implementations. The most indicative difference is, of course, that VESPER models flux timeseries of SPEs, while VTM models macroscopic variables such as Fluence and Peak Flux by fitting their distribution functions. Therefore, we consider the good agreement with VTM an important indication of quality of the VESPER model, especially since it has been designed to be solely datadriven and not reproduce the results of any other model.
Fig. 8
Comparison of VESPER CDFs with ESP, JPL and VTM distributions for a 7 year mission in a solar active period. (a), (b) and (c) Cumulative fluence at 6 MeV, 38 MeV and 166 MeV, respectively. (d), (e) and (f) same for Peak flux. 
5 Conclusions
In this work, a new paradigm for the modelling of the solar proton radiation environment is presented. In contrast to the existing SPE models, that are based on probabilistic descriptions of SPE characteristics (peak flux, cumulative fluence, etc.), the VESPER model produces as fundamental output virtual timeseries of proton differential fluxes that can be used afterwards for the derivation of (secondary) products depending on user defined selected energies and flux intensity thresholds.
It is shown that LF_{2} demonstrates a high correlation with Event Duration, in its summed SLF_{2} form, while it also offers a way of scaling individual proton spectra in a consistent way while maintaining their profiles. We show that the model outputs demonstrate a good degree of agreement with the Event list produced using the SEPEM Reference Dataset v2.0. Furthermore the model demonstrates a quite good agreement with the VTM model regarding the cumulative probability functions produced for Cumulative Fluence and Peak Flux.
The production of virtual flux timeseries environment offers new possibilities in SPE modelling and opens new avenues for the creation of a radiation effect system resulting from the combination of different radiation effect timesseries. VESPER's outputs can be readily combined with tools that calculate radiation effects (Seltzer, 1980; Lei et al., 2002) to produce timeseries of the ionizing or nonionizing dose, internal charging and other such effects. Finally, combining these results with the corresponding outputs of radiation belt models, such as AE9 (Ginet et al., 2013), and/or a GCR model (ISO15390, 2004; O'Neill, 2006; Matthiä et al., 2013), we will be able to provide a complete estimate of the radiation effects on a specified orbit, covering all possible sources of radiation and thus vitally contributing to the planning of future missions. Additionally, we plan to integrate magnetospheric shielding effects in VESPER in a timeseries manner. This will allow further indetail modelling of SPE fluxes in magnetospheric regions where magnetospheric shielding effects should not be ignored. This can be used for modelling SPE fluxes and their effects in orbits were these considerations are vital such as lowearth orbits (LEO) and highly elliptical orbits (HEO). Finally, in the future, we will further our studies and possibly update the VESPER model. Among other things we plan to run the model using the RDS database with the background subtracted for each channel and compare its outputs with the current version. In addition, we will investigate in further detail the impact of the tolerance factor and the establishing of the upper and lower limits on the regression.
Acknowledgements
This work has been supported by ESA/ESTEC contract 4000112863/14/NL/HB in the framework of the HERMES project.
The authors and editor thank Michael Xapsos and an anonymous referee for their assistance in evaluating this paper.
References
 Crosby NB, Bothmer V, Facius R, Griessmeier JM., Moussas X, Panasyuk M, Romanova N, Withers P. 2008. Interplanetary space weather and its planetary connection. Space Weather 6: S01003. [Google Scholar]
 Crosby N, et al. 2015. SEPEM: a tool for statistical modeling the solar energetic particle environment. Space Weather 13: 406–426. [CrossRef] [Google Scholar]
 Cucinotta FA, Hu S, Schwadron NA, Kozarev K, Townsend LW, Kim MHY. 2010. Space radiation risk limits and EarthMoonMars environmental models. Space Weather 8: S00E09. [CrossRef] [Google Scholar]
 Dodd PE, Massengill LW. 2003. Basic mechanisms and modeling of singleevent upset in digital microelectronics. IEEE Trans Nucl Sci 50: 583–602. [CrossRef] [Google Scholar]
 Feynman J, Armstrong TP, DaoGibner L, Silverman S. 1990. New interplanetary proton fluence model. J Spacecr Rocket 27: 403. [Google Scholar]
 Feynman J, Spitale G, Wang J, Gabriel SB. 1993. Interplanetary proton fluence model: JPL 1991. J Geophys Res 98: A8. [Google Scholar]
 Ginet GP, O'Brien TP, Huston SL, Johnston WR, Guild TB, et al. 2013. AE9, AP9 and SPM: new models for specifying the trapped energetic particle and space plasma environment. Space Sci Rev 179: 579–615. [CrossRef] [Google Scholar]
 Heynderickx D, Sandberg I, Jiggens P. 2017. SEPEM reference data set (RDS) v2.0. http://test.sepem.eu/help/SEPEM_RDS_v200.zip. [Google Scholar]
 ISO15390. 2004. Space environment (Natural and artificial) − galactic cosmic ray model. Ref no. ISO 15390:2004(E). [Google Scholar]
 Jiggens PTA, Gabriel SB. 2009. Time distributions of solar energetic particle events: are SEPEs really random? J Geophys Res 114: A10105. [CrossRef] [Google Scholar]
 Jiggens PTA, Gabriel SB, Heynderickx D, Crosby N, Glover A, Hilgers A. 2012. ESA SEPEM project: peak flux and fluence model. IEEE Trans Nucl Sci 59: 1066–1077. [CrossRef] [Google Scholar]
 Lei F, Truscott RR, Dyer CS, Quaghebeur B, Heynderickx D, Nieminen R, Evans H, Daly E. 2002. MULASSIS: a Geant4based multilayered shielding simulation tool. IEEE Trans Nucl Sci 49: 2788–2793. [CrossRef] [Google Scholar]
 Matthiä D, Berger T, Mrigakshi AI, Reitz G. 2013. A readytouse galactic cosmic ray model. Adv Space Res 51: 329–338. [CrossRef] [Google Scholar]
 Nymmik RA. 1999. Probabilistic model for fluences and peak fluxes of solar energetic particles. Rad Meas 30: 287–296. [Google Scholar]
 O'Neill PM. 2006. BadhwarO'Neill galactic cosmic ray model update based on advanced composition explorer (ACE) energy spectra from 1997 to present. Adv Space Res 37: 1727–1733. [NASA ADS] [CrossRef] [Google Scholar]
 Rodriguez JV, Sandberg I, Mewaldt RA, Daglis IA, Jiggens P. 2017. Validation of the effect of crosscalibrated GOES solar proton effective energies on derived integral fluxes by comparison with STEREO observations. Space Weather 15: 290–309. [CrossRef] [Google Scholar]
 Sandberg I, Jiggens P, Heynderickx D, Daglis IA. 2014. Cross calibration of NOAA GOES solar proton detectors using corrected NASA IMP8/GME data. Geophys Res Lett 41: 4435–4441. [CrossRef] [Google Scholar]
 Seltzer SM, SHIELDOSE: a computer code for spaceshielding radiation dose calculations, National Bureau of Standards, NBS Technical Note 1116, US Government Printing Office, Washington, 407 DC, 1980 [Google Scholar]
 Stassinopoulos EG, King JH. 1974. Empirical solar proton model for orbiting spacecraft applications. IEEE Trans Aerosp Electron Syst 10: 442–450. [CrossRef] [Google Scholar]
 Vainio R, et al. 2009. Dynamics of the earth's particle radiation environment. Space Sci Rev 147: 187–231. [Google Scholar]
 Xapsos MA, Summers GP, Burke EA. 1998. Probability model for peak fluxes of solar proton events. IEEE Trans Nucl Sci 45: 2948–2953. [Google Scholar]
 Xapsos MA, Summers GP, Barth JL, Stassinopoulos EG, Burke EA. 1999. Probability model for worst case solar proton event fluences. IEEE Trans Nucl Sci 46: 1481–1485. [Google Scholar]
 Xapsos MA, Stauffer C, Barth JL, Burke EA. 2006. Solar particle events and selforganized criticality: are deterministic predictions of events possible? IEEE Trans Nucl Sci 53: 1839–1843. [CrossRef] [Google Scholar]
Cite this article as: AminalragiaGiamini S, Sandberg I, Papadimitriou C, Daglis IA, Jiggens P. 2018. The virtual enhancements − solar proton event radiation (VESPER) model. J. Space Weather Space Clim. 8: A06
All Tables
The first 10 energy bins and logmean energies of the RDSv2 proton flux dataset.
All Figures
Fig. 1
The LF_{2} timeseries calculated for the entire SEPEM reference dataset v2.0. The background is denoted by the green line and the flux enhancements denoted by red. 

In the text 
Fig. 2
The timeseries for two defined SPEs. Event A: 16Feb1984 to 24Feb1984, Event B: 11Jul1982 to 14Jul1982. (a) and (b) flux timeseries of the ten energy channels, (c) and (d) timeseries of the LF_{2} variable calculated from the fluxes. 

In the text 
Fig. 4
Cross plots of Duration versus SLF_{2}. The red points correspond to the calculated values based on the analysis of RDSv2. (a) The dashed lines denote the vicinity of the VE which is used to calculate an SLF_{2} value, adjust the Duration value and select a seed for the Virtual Event. (b) Black dots correspond to calculated values for the VEs. 

In the text 
Fig. 3
CDFs for (a) the Duration and (b) the WaitTime variables and their fits with exponential cutoff powerlaws. 

In the text 
Fig. 5
(a) Solar proton event, 27Nov1989 to 05Dec1989 (b) Virtual event rescaled to higher Duration and SLF_{2} values. 

In the text 
Fig. 6
(a) and (b) CDFs of cumulative differential fluence produced from the VESPER model for 1 and 7 year solar active periods, respectively. Each curve presents one of the channels in the 5–200 MeV proton energy range, (c) and (d) same for peak differential flux. 

In the text 
Fig. 7
Comparison of (a) total fluence and (b) median flux values for RDSv2 and VESPER virtual data derived using a moving threshold on the respective LF_{2} timeseries. Crosses denote RDSv2 and continuous lines virtual data. 

In the text 
Fig. 8
Comparison of VESPER CDFs with ESP, JPL and VTM distributions for a 7 year mission in a solar active period. (a), (b) and (c) Cumulative fluence at 6 MeV, 38 MeV and 166 MeV, respectively. (d), (e) and (f) same for Peak flux. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.