-Solar energetic particles (SEPs) constitute an important component of the radiation environment in interplanetary space. Accurate modeling of SEP events is crucial for the mitigation of radiation hazards in spacecraft design. In this study we present two new statistical models of high energy solar proton fluences based on ground level enhancement (GLE) observations during solar cycles 19-24. As the basis of our modeling, we utilize a four parameter double power law function (known as the Band function) fits to integral GLE fluence spectra in rigidity. In the first model, the integral and differential fluences for protons with energies between 10 MeV and 1 GeV are calculated using the fits, and the distributions of the fluences at certain energies are modeled with an exponentially cut-off power law function. In the second model, we use a more advanced methodology: by investigating the distributions and relationships of the spectral fit parameters we find that they can be modeled as two independent and two dependent variables. Therefore, instead of modeling the fluences separately at different energies, we can model the shape of the fluence spectrum. We present examples of modeling results and show that the two methodologies agree well except for a short mission duration (1 year) at low confidence level. We also show that there is a reasonable agreement between our models and three well-known solar proton models (JPL, ESP and SEPEM), despite the differences in both the modeling methodologies and the data used to construct the models.

This paper reviews Vlasov-based numerical methods used to model plasma in space physics and astrophysics. Plasma consists of collectively behaving charged particles that form the major part of baryonic matter in the Universe. Many concepts ranging from our own planetary environment to the Solar system and beyond can be understood in terms of kinetic plasma physics, represented by the Vlasov equation. We introduce the physical basis for the Vlasov system, and then outline the associated numerical methods that are typically used. A particular application of the Vlasov system is Vlasiator, the world’s first global hybrid-Vlasov simulation for the Earth’s magnetic domain, the magnetosphere. We introduce the design strategies for Vlasiator and outline its numerical concepts ranging from solvers to coupling schemes. We review Vlasiator’s parallelisation methods and introduce the used high-performance computing (HPC) techniques. A short review of verification, validation and physical results is included. The purpose of the paper is to present the Vlasov system and introduce an example implementation, and to illustrate that even with massive computational challenges, an accurate description of physics can be rewarding in itself and significantly advance our understanding. Upcoming supercomputing resources are making similar efforts feasible in other fields as well, making our design options relevant for others facing similar challenges.

For decades, monochromatic large‐scale ultralow frequency (ULF) waves with a period of about 30 s have been observed upstream of the quasi‐parallel bow shock. These waves typically propagate obliquely with respect to the interplanetary magnetic field (IMF), while the growth rate for the instability causing the waves is maximized parallel to the magnetic field. It has been suggested that the mechanism for the oblique propagation concerns wave refraction due to the spatial variability of the suprathermal ions, originating from the E × B drift component. We investigate the ULF foreshock under a quasi‐radial IMF with Vlasiator, which is a newly developed global hybrid‐Vlasov simulation solving the Vlasov equation for protons, while electrons are treated as a charge‐neutralizing fluid. We observe the generation of the 30 s ULF waves and compare their properties to previous literature and multipoint Time History of Events and Macroscale Interactions during Substorms (THEMIS) spacecraft observations. We find that Vlasiator reproduces the foreshock ULF waves in all reported observational aspects. We conclude that the variability of the density and velocity of the reflected back streaming ions determines the large‐scale structure of the foreshock, which affects the wave frequency, wavelength, and oblique propagation. We conclude that the wave refraction may also be at work for radial IMF conditions, which has earlier been thought of as an exception to the refraction mechanism due to the small E × B drift component. We suggest that additional refraction may be caused by the large‐scale spatial variability of the density and velocity of the back streaming ions.

We present results from a first study of the local reconnection rate and reconnection site motion in a 2D‐3V global magnetospheric self‐consistent hybrid‐Vlasov simulation with due southward interplanetary magnetic field. We observe magnetic reconnection at multiple locations at the dayside magnetopause and the existence of magnetic islands, which are the 2‐D representations of flux transfer events. The reconnection locations (the X lines) propagate over significant distances along the magnetopause, and reconnection does not reach a steady state. We calculate the reconnection rate at the location of the X lines and find a good correlation with an analytical model of local 2‐D asymmetric reconnection. We find that despite the solar wind conditions being constant, the reconnection rate and location of the X lines are highly variable. These variations are caused by magnetosheath fluctuations, the effects of neighboring X lines, and the motion of passing magnetic islands.

SEPServer is a three-year collaborative project funded by the seventh framework programme (FP7-SPACE) of the European Union. The objective of the project is to provide access to state-of-the-art observations and analysis tools for the scientific community on solar energetic particle (SEP) events and related electromagnetic (EM) emissions. The project will eventually lead to better understanding of the particle acceleration and transport processes at the Sun and in the inner heliosphere. These processes lead to SEP events that form one of the key elements of space weather. In this paper we present the first results from the systematic analysis work performed on the following datasets: SOHO/ERNE, SOHO/EPHIN, ACE/EPAM, Wind/WAVES and GOES X-rays. A catalogue of SEP events at 1 AU, with complete coverage over solar cycle 23, based on high-energy (~68-MeV) protons from SOHO/ERNE and electron recordings of the events by SOHO/EPHIN and ACE/EPAM are presented. A total of 115 energetic particle events have been identified and analysed using velocity dispersion analysis (VDA) for protons and time-shifting analysis (TSA) for electrons and protons in order to infer the SEP release times at the Sun. EM observations during the times of the SEP event onset have been gathered and compared to the release time estimates of particles. Data from those events that occurred during the European day-time, i.e., those that also have observations from ground-based observatories included in SEPServer, are listed and a preliminary analysis of their associations is presented. We find that VDA results for protons can be a useful tool for the analysis of proton release times, but if the derived proton path length is out of a range of 1 AU < s [ 3 AU, the result of the analysis may be compromised, as indicated by the anti-correlation of the derived path length and release time delay from the associated X-ray flare. The average path length derived from VDA is about 1.9 times the nominal length of the spiral magnetic field line. This implies that the path length of first-arriving MeV to deka-MeV protons is affected by interplanetary scattering. TSA of near-relativistic electrons results in a release time that shows significant scatter with respect to the EM emissions but with a trend of being delayed more with increasing distance between the flare and the nominal footpoint of the Earth-connected field line.

Master athletes maintain high physical activity levels and have better health than age-matched non-athletes. World records show accelerated declines after age 70 in swimming, long-distance running and sprint performance. However, less is known about age-related performance declines in the general master athlete population and whether decline rates differ between disciplines and genders. We interrogated a dataset including all track and field athletes of North Rhine from 2001 to 2014 to assess age-related changes in performance. 27,088 results of athletes between 11 and 89 years of age in 12 disciplines were analyzed by regression statistics. The analyses showed an accelerated decline beyond the age of 70 in sprint, middle- and long-distance running, while in throwing and jumping disciplines the performance continued a linear decline. Patterns of decline differed between men and women. The steepest declines were observed in javelin throw and 400 m (women), and in pole vault and 800 m (men). In conclusion, performance declines in aging depend more on the specific profile of requirements than previously assumed.

We investigate mirror mode structures in the Earth's magnetosheath using our global kinetic Vlasiator simulation, which models ion behavior through their distribution function and treats electrons as a charge‐neutralizing fluid. We follow the evolution of waves as they advect along velocity streamlines through the magnetosheath. We find that mirror mode waves are observed preferentially in the quasi‐perpendicular magnetosheath along velocity streamlines that enter the sheath in the vicinity of the foreshock ULF wave boundary where there are enough initial perturbations in the plasma for the mirror modes to grow, and the plasma properties fulfill the mirror instability condition better than elsewhere in the magnetosheath. We test selection criteria defined by previous studies and show that the spatial extent over which mirror modes occur ranges from much of the magnetosheath on the quasi‐perpendicular side of the subsolar point to very small isolated regions depending on the criteria.

Abstract. The key dynamics of the magnetotail have been researched for decades and have been associated with either three-dimensional (3-D) plasma instabilities and/or magnetic reconnection. We apply a global hybrid-Vlasov code, Vlasiator, to simulate reconnection self-consistently in the ion kinetic scales in the noon-midnight meridional plane, including both dayside and nightside reconnection regions within the same simulation box. Our simulation represents a numerical experiment, which turns off the 3-D instabilities but models ion-scale reconnection physically accurately in 2-D. We demonstrate that many known tail dynamics are present in the simulation without a full description of 3-D instabilities or without the detailed description of the electrons. While multiple reconnection sites can coexist in the plasma sheet, one reconnection point can start a global reconfiguration process, in which magnetic field lines become detached and a plasmoid is released. As the simulation run features temporally steady solar wind input, this global reconfiguration is not associated with sudden changes in the solar wind. Further, we show that lobe density variations originating from dayside reconnection may play an important role in stabilising tail reconnection.