ABSTRACT:Inhomogeneous B 0 -magnetic fields generate distortion in magnetic resonance images, particularly those produced using echo planar imaging, and are responsible for signal reduction due to intravoxel dephasing in gradient echo experiments. Such effects increase in magnitude in proportionality with the static field strength, and with the growing use of high-field (3 T and above) systems in medical imaging, it is increasingly important to be able to quantify field inhomogeneities. Here, we describe the implementation and use of a method for rapidly calculating frequency shifts due to spatially varying magnetic susceptibility that is based on an approach previously used to calculate long-range dipolar field effects. The method relies on a simple expression that relates the three-dimensional Fourier transforms of the magnetization distribution and the field, and can naturally include the effect of the sphere of Lorentz. It has been used to evaluate field inhomogeneity in the head due to the variation of magnetic susceptibility with tissue type and to calculate the change in field inhomogeneity that occurs due to small rotations of the head. In addition, this approach has been used to simulate the effect of lung volume changes in generating respiration induced resonant offsets in the brain.