Phase estimation for thermal Gaussian states