A Comparison of the Galactic Cosmic Ray H, He and C/O Nuclei Spectra Measured Between ~5 and 500 MeV/nuc Beyond 122 AU at Voyager 1 with the Predictions of a Diffusion Model for Propagation in the Galaxy

After the disappearance of lower energy heliospheric particles at Voyager 1 starting on August 25th, 2012, spectra of H, He and C/O nuclei were revealed that resembled those to be expected for galactic cosmic rays. These spectra had intensity peaks in the range of 30-60 MeV, decreasing at both lower energies down to a few MeV and at higher energies up to several hundred MeV. We have modeled the propagation of these particles in the galaxy using an updated Leaky Box Diffusion model which determines the spectra of these components from ~2 MeV to >200 GeV. The key parameters used in the model are a galactic input spectrum ~P^-2.24, the same for all components and independent of rigidity, and a diffusion coefficient that is ~P^0.5 above a lower rigidity and increases ~beta^-1.0 below a lower rigidity ~0.56 GV. These same parameters also fit the high energy H and He data from ~10-200 GeV/nuc from the PAMELA and BESS experiments. The new Voyager spectra for all three nuclei are thus consistent with rigidity spectra ~P^-2.24 from the lowest energies to at least 100 GeV. Deviations from this spectrum can reasonably be attributed to propagation effects. Some deviations between the calculated and newly observed spectra are noted, however, below ~30 MeV/nuc, particularly for C/O nuclei, that could be significant regarding the propagation and sources of these particles.

The Unfolding of the Spectra of Low Energy Galactic Cosmic Ray H and He Nuclei as the Voyager 1 Spacecraft Exits the Region of Heliospheric Modulation

This paper describes the unfolding of the solar modulated galactic cosmic ray H and He nuclei spectra beyond ~105 AU in the heliosheath. Between 2008.0 and 2012.3 when Voyager 1 went from about 105 to 120.5 AU the spectral intensities of these two components between about 30 and 500 MeV/nuc unfolded (increased) in a manner consistent with an average modulation potential decrease ~5 MV per AU as described by a Parker like cosmic ray transport in the heliosphere where the overall modulation is described by a modulation potential in MV. Between 120.5 and 121.7 AU, however, as a result of two sudden intensity increases starting on May 8th and August 25th, 2012, this modulation potential decreased by ~80 MV and spectra resembling possible local interstellar spectra for H and He were revealed. Considering these spectra to be the local interstellar spectra would imply that almost 1/3 of the total modulation potential of about 270 MV required to explain the spectra of these components observed at the Earth must occur in just a 1 AU radial interval in the outer heliosheath. As a result about ~80% of the total modulation potential observed at the Earth at this time occurs in the heliosheath itself. The remaining 20% of the total modulation occurs inside the heliospheric termination shock. The details of these intensity changes and their description by a simple modulation model are discussed.

20-125 mev/nuc cosmic ray carbon nuclei intensities between 2004-2010 in solar cycle #23 as measured near the earth, at voyager 2 and also in the heliosheath at voyager 1 - modulation in a two zone heliospehre

The recovery of cosmic ray Carbon nuclei of energy ~20-125 MeV/nuc in solar cycle #23 from 2004 to 2010 has been followed at three locations, near the Earth using ACE data and at V2 between 74-92 AU and also at V1 beyond the heliospheric termination shock at between 91-113 AU. To describe the observed intensity changes and to predict the absolute intensities measured at all three locations we have used a simple spherically symmetric (no drift) two-zone heliospheric transport model with specific values for the diffusion coefficient in both the inner and outer zones. The diffusion coefficient in the outer zone is determined to be ~5-10 times smaller than that in the inner zone out to 90 AU. For both V1 and V2 the calculated C nuclei intensities agree within an average of pm 10% with the observed intensities. Because of this agreement between V1 and V2 observations and predictions there is no need to invoke an asymmetrical squashed heliosphere or other effects to explain the V2 intensities relative to V1 as is the case for He nuclei. The combination of the diffusion parameters used in this model and the interstellar spectrum give an unusually low overall solar modulation parameter phi = 250 MV to describe the Carbon intensities observed at the Earth in 2009. At all times both the observed and calculated spectra are very closely ~ E1.0 as would be expected in the adiabatic energy loss regime of solar modulation.

Sudden Intensity Increases and Radial Gradient Changes of Cosmic Ray MeV Electrons and Protons Observed at Voyager 1 beyond 111 AU in the Heliosheath

Voyager 1 has entered regions of different propagation conditions for energetic cosmic rays in the outer heliosheath beginning at a distance of about 111 AU from the Sun. This conclusion is based on the fact that the low energy 6-14 MeV galactic electron intensity suddenly increased by ~20% over a time period leg 10 days and the electron radial intensity gradient abruptly decreased from ~19%/AU to ~8%/AU at 2009.7 at a radial distance of 111.2 AU. A sudden radial gradient change was also observed at this time for >200 MeV protons. The gradients were constant during the time period before and after the electron increase. At about 2011.2 at a distance of 116.6 AU a second abrupt intensity increase was observed, this time for both electrons and protons. The increase for electrons was ~25% and occurred over a time period ~15 days or less. For >200 MeV protons the increase at this time was ~5% (unusually large) and occurred over a longer time period ~50 days. Between about 2011.2 and 2011.6, radial intensity gradients ~18%/AU and 3%/AU were observed for electrons and protons, respectively. These gradients were very similar to those observed for these particles before the 1st sudden increase at 2009.7. These large positive gradients observed after 2011.2 indicate that V1, although it has entered a different propagation region, is still within the overall heliospheric modulating region at least up to a time ~2011.6 (118.0 AU). In this paper we will discuss these events in more detail and consider possibilities for their explanation that have recently been suggested.

Recovery of 150-250 MeV/nuc Cosmic Ray Helium Nuclei Intensities Between 2004-2010 Near the Earth, at Voyager 2 and Voyager 1 in the Heliosheath - A Two Zone Helioshpere

The recovery of cosmic ray He nuclei of energy ~150-250 MeV/nuc in solar cycle #23 from 2004 to 2010 has been followed at the Earth using IMP and ACE data and at V2 between 74-92 AU and also at V1 beyond the heliospheric termination shock (91-113 AU). The correlation coefficient between the intensities at the Earth and at V1 during this time period is remarkable (0.921), after allowing for a ~0.9 year delay due to the solar wind propagation time from the Earth to the outer heliosphere. To describe the intensity changes and to predict the absolute intensities measured at all three locations we have used a simple spherically symmetric (no drift) two-zone heliospheric transport model with specific values for the diffusion coefficient in both the inner and outer zones. The diffusion coefficient in the outer zone, assumed to be the heliosheath from about 90 to 120 (130) AU, is determined to be ~5 times smaller than that in the inner zone out to 90 AU. This means the Heliosheath acts much like a diffusing barrier in this model. The absolute magnitude of the intensities and the intensity changes at V1 and the Earth are described to within a few percent by a diffusion coefficient that varies with time by a factor ~4 in the inner zone and only a factor of ~1.5 in the outer zone over the time period from 2004-2010. For V2 the observed intensities follow a curve that is as much as 25% higher than the calculated intensities at the V2 radius and at times the observed V2 intensities are equal to those at V1. At least one-half of the difference between the calculated and observed intensities between V1 and V2 can be explained if the heliosphere is squashed by ~10% in distance (non-spherical) so that the HTS location is closer to the Sun in the direction of V2 compared to V1.