Building further upon the high spatial resolution offered by ultrasonic imaging and the high optical contrast yielded by laser excitation of photoacoustic imaging, and exploiting the temperature dependence of photoacoustic signal amplitudes, this paper addresses the question whether the rich information given by multispectral optoacoustic tomography (MSOT) allows to obtain 3D temperature images. Numerical simulations and experimental results are reported on agarose phantoms containing gold nanoparticles and the effects of shadowing, reconstruction flaws, etc. on the accuracy are determined.