The inherent sensitivity of multiple gradient-echo sequences to the chemical shift is exploited to rapidly map muscle water frequency shifts caused by ultrasonic heating. The use of multiple echoes is shown to offer several advantages over single gradient-echo approaches previously proposed for temperature measurement. An increase in the effective bandwidth significantly reduces aliasing problems observed with single gradient-echo methods in high temperature applications. Of greater significance is the improved immunity to intrascan motion found for multi-echo versus single echo gradient methods, making the former more attractive for clinical applications. Finally, a sensitivity to the presence of multiple spectral components unavailable with single gradient-echo methods is obtained.