Noninvasive magnetic resonance temperature maps that are used to monitor thermal ablation of tissue are described. In magnetic resonance images, thermally induced proton nuclear magnetic resonance frequency shifts, and changes in the longitudinal relaxation time produce both phase and magnitude changes in the MR signal. Temperature maps with improved sensitivity are derived from the complex-difference nuclear magnetic resonance signal. Bovine muscle specimens were heated with focused ultrasound to model thermal surgery and create a known thermal distribution to test the method. Resulting MR images acquired in 2 s produce temperature maps with 1 min resolution and 2 degrees C temperature sensitivity. The temperature sensitivity was increased by extending the acquisition to 5 s, by decreasing the receiver bandwidth, and increasing the echo time.