Precise time-dependent measurements of the Z = 2 component in the cosmic radiation provide crucial information about the propagation of charged particles through the heliosphere. The PAMELA experiment, with its long flight duration (15th June 2006 - 23rd January 2016) and the low energy threshold (80 MeV/n) is an ideal detector for cosmic ray solar modulation studies. In this paper, the helium nuclei spectra measured by the PAMELA instrument from July 2006 to December 2009 over a Carrington rotation time basis are presented. A state-of-the-art three-dimensional model for cosmic-ray propagation inside the heliosphere was used to interpret the time-dependent measured fluxes. Proton-to-helium flux ratio time profiles at various rigidities are also presented in order to study any features which could result from the different masses and local interstellar spectra shapes.

The local interstellar spectra (LIS's) for galactic cosmic rays (CR's) cannot be directly observed at the Earth below certain energies, because of solar modulation in the heliosphere. With Voyager 1 crossing the heliopause in 2012, in situ experimental LIS data below 100 MeV/nuc can now constrain computed galactic CR spectra. Using galactic propagation models, galactic electron, proton and light nuclei spectra can be computed, now more reliably as very LIS's. Using the Voyager 1 observations made beyond the heliopause, and the observations made by the PAMELA experiment in Earth orbit for the 2009 solar minimum, as experimental constraints, we simultaneously reproduced the CR electron, proton, Helium and Carbon observations by implementing the GALPROP code. Below about 30 GeV/nuc solar modulation has a significant effect and a comprehensive three-dimensional (3D) numerical modulation model is used to compare the computed spectra with the observed PAMELA spectra possible at these energies. Subsequently the computed LIS's can be compared over as wide a range of energies as possible. The simultaneous calculation CR spectra with a single propagation model allows the LIS's for positrons, Boron and Oxygen to also be inferred. This implementation of the most comprehensive galactic propagation model (GALPROP), alongside a sophisticated solar modulation model to compute CR spectra for comparison with both Voyager 1 and PAMELA observations over a wide energy range, allows us to present new self-consistent very LIS's (and expressions) for electrons, positrons, protons, Helium, Carbon, Boron and Oxygen for the energy range of 3 MeV/nuc to 100 GeV/nuc.

Heliospheric modulation conditions were unusually quite during the last solar minimum activity between Solar Cycles 23/24. Fortunately, the PAMELA space-experiment measured six-month averaged Galactic positron spectra for the period July 2006 to December 2009, over an energy range of 80 MeV to 30 GeV, which is important for solar modulation. The highest level of Galactic positrons was observed at Earth during the July-December 2009 period. A well-established, comprehensive three-dimensional (3D) numerical model is applied to study the modulation of the observed positron spectra. This model had been used previously to understand the modulation of Galactic protons and electrons also measured by PAMELA for the same period. First, a new very local interstellar spectrum for positrons is constructed, using the well-known GALPROP code together with the mentioned PAMELA observations. The 3D model is used to distinguish between the dominant mechanisms responsible for the heliospheric modulation of Galactic positrons, and to understand the effect of particle drift during this unusual minimum in particular, which is considered diffusion dominant, even though particle drift still had a significant role in modulating positrons. Lastly, the expected intensity of Galactic positrons during an A>0 polarity minimum, with similar heliospheric conditions than for 2006-2009, is predicted to be higher than what was observed by PAMELA for the 2006-2009 unusual minimum.

With the cosmic ray observations made by the Voyager 1 spacecraft outside the dominant modulating influence of the heliosphere, the comparison of computed galactic spectra with experimental data at lower energies is finally possible. Spectra for specifically protons, Helium and Carbon nuclei, computed by galactic propagation models, can now be compared with observations at low energies from Voyager 1 and at high energies from the PAMELA space detector at Earth. We set out to reproduce the Voyager 1 observations in the energy range of 6 MeV/nuc to 60 MeV/nuc, and the PAMELA spectrum above 50 GeV/nuc, using the GALPROP code, similarly to our previous study for Voyager 1 electrons. By varying the galactic diffusion parameters in the GALPROP plain diffusion model, specifically the rigidity dependence of spatial diffusion, and then including reacceleration, we compute spectra simultaneously for galactic protons, Helium and Carbon.We present new local interstellar spectra, with expressions for the energy range of 3 MeV/nuc to 100 GeV/nuc, which should be of value for solar modulation modeling.

Given that the bulk of the galactic cosmic rays (CRs) originates from transient point-like sources, such as supernova remnants, the flux of the CR primary component measured at Earth depends on the local source history. Whereas the secondary omponent shows little to no variation due to nearby sources, and instead, should depend on the global distribution of CRs. Steady-state, rotational symmetric models (2D) of CR propagation assume smeared-out CR sources in the Galaxy and cannot properly take into account the influence of nearby point sources. In this paper, we infer evidence of nearby sources by illustrating that a 2D propagation model will not describe CR primaries and secondaries equally well. We adapted the 2D version of the GALPROP code to a compute-cluster environment and perform parameter studies comparing CR spectra with primary and secondary CR data separately. Doing a parameter study, one may expect different best fit values looking at the primary and secondary CR components separately, as it is unlikely that the source history mimicked by the 2D models coincides with the real local source history. We find that the primaries and secondaries are fitted with differing best fit parameters, with the results indicating that the 2D model being more suited to modelling secondary CRs. The insufficient modelling of primary CRs can be contributed to the possible presence of local CR sources.