The interstellar medium (ISM) of star-forming galaxies is magnetized and turbulent. Cosmic rays (CRs) propagate through it, and those with energies from ∼ GeV − TeV are likely subject to the streaming instability, whereby the wave damping processes balances excitation of resonant ionic Alfv én waves by the CRs, reaching an equilibrium in which the propagation speed of the CRs is very close to the local ion Alfv én velocity. The transport of streaming CRs is therefore sensitive to ionic Alfv én velocity fluctuations. In this paper we systematically study these fluctuations using a large ensemble of compressible MHD turbulence simulations. We show that for sub-Alfv énic turbulence, as obtains for a strongly magnetized ISM, the ionic Alfv én velocity probability density function (PDF) is determined solely by the density fluctuations from shocked gas forming parallel to the magnetic field, and we develop analytical models for the ionic Alfv én velocity PDF up to second moments. For super-Alfv énic turbulence, magnetic and density fluctuations are correlated in complex ways, and these correlations as well as contributions from the magnetic fluctuations sets the ionic Alfv én velocity PDF. We discuss the implications of these findings for underlying "macroscopic" diffusion mechanisms in CRs undergoing the streaming instability, including modeling the macroscopic diffusion coefficient for the parallel transport in a sub-Alfv énic plasma. We also describe how, for highly-magnetized turbulent gas, the gas density PDF, and hence column density PDF, can be used to access information about ionic Alfv én velocity structure from observations of the magnetized ISM.