A two-dimensional, time-dependent numerical model is used to calculate the modulation of cosmic rays in the heliosphere. Computations are compared to spacecraft observations in the inner and outer heliosphere. It is shown that the model produces cosmic ray proton intensities compatible to different spacecraft observations on a global scale, at Earth and along both Voyager spacecraft trajectories. The study reveals that when the same modulation parameters, which resulted in compatible intensities along Voyager 1, were assumed along the Voyager 2 trajectory, the model failed to reproduce the observations. The study also found that any change in diffusion parameters alone could not reproduce the cosmic ray observations along Voyager 2 so that changes to the heliospheric geometry were necessary i.e the computed intensities along both Voyager trajectories suggest that the heliosphere is asymmetric. Furthermore, $E >$ 70 MeV and 133-242 MeV proton intensities along Voyager 1 and 2 trajectories are predicted from end of 2012 onwards. It is shown that the computed intensities along Voyager 1 increase with an almost constant rate up to the heliopause. However, the model shows that Voyager 2 is still under the influence of temporal solar activity changes because of its relatively large distance to the heliopause. Along the Voyager 2 trajectory, the intensities remained generally constant for some time and should soon start to increase steadily.

The time-dependent modulation of galactic cosmic rays in the heliosphere is studied by computing intensities using a time-dependent modulation model. By introducing recent theoretical advances in the transport coefficients in the model, computed intensities are compared with Voyager 1, International Monitoring Platform (IMP) 8, and Ulysses proton observations in search of compatibility. The effect of different modulation parameters on computed intensities is also illustrated. It is shown that this approach produces, on a global scale, realistic cosmic-ray proton intensities along the Voyager 1 spacecraft trajectory and at Earth upto ~2004, whereafter the computed intensities recovers much slower towards solar minimum than observed in the inner heliosphere. A modified time dependence in the diffusion coefficients is proposed to improve compatibility with the observations at Earth after ~2004. This modified time dependence led to an improved compatibility between computed intensities and the observations along the Voyager 1 trajectory and at Earth even after ~2004. An interesting result is that the cosmic-ray modulation during the present polarity cycle is not determined only by changes in the drift coefficient and tilt angle of the wavy current sheet, but is also largely dependent on changes in the diffusion coefficients.