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.

💡 Research Summary

The paper investigates the recovery of 20‑125 MeV/nucleon carbon nuclei intensities during solar cycle 23 (2004‑2010) as observed at three distinct locations: near Earth (ACE spacecraft), Voyager 2 (V2) between 74‑92 AU, and Voyager 1 (V1) beyond the heliospheric termination shock (91‑113 AU). The authors aim to reproduce the temporal evolution of the carbon fluxes and to predict absolute intensities at all three points using a very simple, spherically symmetric, two‑zone heliospheric transport model that neglects particle drifts.

Model structure
The heliosphere is divided into an inner zone (0‑90 AU) and an outer zone (90‑120 AU). Each zone is assigned its own spatial diffusion coefficient, K₁ for the inner region and K₂ for the outer region. By fitting the observed data, the authors find that K₂ must be roughly 5‑10 times smaller than K₁, implying that particle diffusion is strongly suppressed in the outer heliosphere (the region just inside the heliosheath). The model solves the one‑dimensional Parker transport equation in spherical coordinates, incorporating solar‑wind convection, adiabatic energy loss, and spatial diffusion. The outer boundary condition is the local interstellar spectrum (LIS) for carbon, while the inner boundary is set by the measured Earth intensity.

Key parameter – modulation potential
A single modulation potential, φ, of 250 MV is adopted to connect the LIS to the Earth‑based spectrum in 2009. This value is unusually low for a solar minimum, reflecting the exceptionally weak solar modulation during that period. The low φ, together with the diffusion coefficients, yields an energy spectrum that follows a power law ∝ E⁻¹·⁰ across the entire 20‑125 MeV/nuc range, which is characteristic of a regime dominated by adiabatic energy losses rather than diffusion‑driven spectral shaping.

Results and validation
When the model is run with the above parameters, the calculated carbon intensities at V1 and V2 agree with the measured values to within an average of ±10 %. This level of agreement is achieved without invoking any asymmetry of the heliosphere (e.g., a “squashed” shape) or additional physical effects. In contrast, earlier studies of helium nuclei required such asymmetries to reconcile V1 and V2 observations. The model also reproduces the Earth‑based intensity time profile, including the peak in 2009, confirming that the same diffusion parameters and low φ can simultaneously explain data from the inner heliosphere out to the heliosheath.

Physical interpretation
The suppression of diffusion in the outer zone suggests that the plasma conditions (magnetic turbulence, field strength, and solar‑wind speed) beyond ~90 AU are markedly different from those inside the termination shock, leading to longer residence times for cosmic‑ray particles. Consequently, the recovery of carbon intensities observed at V1 and V2 is delayed relative to the Earth, a behavior captured naturally by the two‑zone diffusion contrast. The E⁻¹·⁰ spectral shape across the whole energy band confirms that adiabatic cooling dominates the modulation process for these low‑energy carbon nuclei.

Implications

1. Simplicity vs. realism – The success of a drift‑free, spherically symmetric two‑zone model demonstrates that, for carbon nuclei in the 20‑125 MeV/nuc range, complex three‑dimensional structures or drift effects are not required to explain the observations. This greatly reduces the number of free parameters and simplifies forecasting efforts.
2. Species dependence – The fact that carbon does not require a squashed heliosphere, whereas helium does, highlights that modulation can be species‑dependent, likely because of differing rigidity spectra and charge‑to‑mass ratios.
3. Solar minimum conditions – The derived φ = 250 MV provides a quantitative benchmark for the unusually weak solar modulation during the 2009 minimum, useful for calibrating other cosmic‑ray transport studies and for space‑weather modeling.
4. Future missions – The model framework can be readily extended to interpret data from New Horizons, Parker Solar Probe, and upcoming interstellar probes, especially when combined with more detailed measurements of solar‑wind turbulence in the outer heliosphere.

Conclusions
The authors conclude that a two‑zone diffusion model with K₂ ≈ K₁/5‑10 and a low modulation potential of 250 MV accurately reproduces carbon nuclei intensities from the Earth out to beyond the termination shock during the recovery phase of solar cycle 23. No additional asymmetries or drift effects are needed, and the observed E⁻¹·⁰ spectra confirm that adiabatic energy loss is the dominant modulation mechanism for these particles. This work provides a parsimonious yet robust description of low‑energy carbon cosmic‑ray transport and sets a baseline for future comparative studies of other species and heliospheric conditions.