Thermal safety indices for diagnostic ultrasound beams are calculated under the assumption that the sound propagates under linear conditions. A non-axisymmetric finite difference model is used to solve the KZK equation, and so to model the beam of a diagnostic scanner in pulsed Doppler mode. Beams from both a uniform focused rectangular source and a linear array are considered. Calculations are performed in water, and in attenuating media with tissue-like characteristics. Attenuating media are found to exhibit significant nonlinear effects for finite-amplitude beams. The resulting loss of intensity by the beam is then used as the source term in a model of tissue heating to estimate the maximum temperature rises. These are compared with the thermal indices, derived from the properties of the water-propagated beams. ©2000 American Institute of Physics.