The second-order statistics of the temperature and thermal stresses are evaluated in an axisymmetrically heated functionally graded annular disc of variable thickness with spatially random heat transfer coefficients (HTCs) on the heat dissipation surfaces. This annular disc is assumed to have arbitrary variations in the HTCs and material composition only along the radial direction. The randomness in the HTCs is considered to be a random field. The stochastic temperature field is analysed by considering the annular disc to be multilayered one with stepwise thickness variation, which has constant material properties and spatially constant but random HTCs in each layer. In order to evaluate the statistics, the Monte Carlo simulation method is employed with analytical solutions for the deterministic temperature and thermal stresses. The analytical solution for the thermal stresses is obtained thorough the piecewise power function approximation for Young's modulus. Numerical results demonstrate the effects of the material composition profile and thickness variation on the statistics of the temperature and thermal stresses.