Purpose: Hyperthermia is a cancer treatment technique that could be delivered as a stand-alone modality or in conjunction with chemotherapy or radiation therapy. Noninvasive and real-time temperature monitoring of the heated tissue improves the efficacy and safety of the treatment. A temperature-sensitive acoustic parameter is required for ultrasound-based thermometry. In this paper the amplitude and the energy of the acoustic harmonics of the ultrasound backscattered signal are proposed as suitable parameters for noninvasive ultrasound thermometry. Materials and methods: A commercial high frequency ultrasound imaging system was used to generate and detect acoustic harmonics in tissue-mimicking gel phantoms and ex vivo bovine muscle tissues. The pressure amplitude and the energy content of the backscattered fundamental frequency (p(1) and E-1), the second (p(2) and E-2) and the third (p(3) and E-3) harmonics were detected in pulse-echo mode. Temperature was increased from 26 degrees to 46 degrees C uniformly through both samples. The amplitude and the energy content of the harmonics and their ratio were measured and analysed as a function of temperature. Results: The average p(1), p(2) and p(3) increased by 69%, 100% and 283%, respectively as the temperature was elevated from 26 degrees to 46 degrees C in tissue samples. In the same experiment the average E-1, E-2 and E-3 increased by 163%, 281% and 2257%, respectively. A similar trend was observed in tissue-mimicking gel phantoms. Conclusions: The findings suggest that the harmonics generated due to nonlinear ultrasound beam propagation are highly sensitive to temperature and could potentially be used for noninvasive ultrasound tissue thermometry.