This paper reports measurements made by the cosmic ray subsystem onboard Voyager 1 and 2 in the Jovian magnetosphere. Energy spectra of oxygen ions in the energy range 1--20 MeV/nuc between 5 and 20 RJ are presented and used to calculate phase space densities. There is a steep positive radial gradient in the phase space density of the energetic oxygen ions in this region, indicating an inward diffusive flow. Solutions of the diffusion equation assuming a diffusion coefficient D and loss lifetime &tgr; of the forms D = D0Ln and &tgr; = &tgr;0Lm, where D0, &tgr;0, n, and m are constants, and L is the McIlwain parameter, are fit to the radial phase space density profile of oxygen ions with magnetic moments of 680 MeV/nuc-G. The best fits are found to have n+m≈6 and 3400 MeV/nuc-G diffuse inward across 10 RJ is 5¿1021¿1 ions s-1. The observations suggest that oxygen and sulfur ions in the Io plasma torus diffuse radially outward, are nonadiabatically accelerated in some region outside 17 RJ and then diffuse inward and outward from the acceleration region.