In this paper,an analytical method for predicting the time-dependent one-dimension temperature field in a multilayered pavement system is proposed on the basis of the Laplace transform as well as of the inverse Laplace transform being resolved numerically by the Gaussian quadrature formula.In the method,first,a time-dependent one-dimension mathematical temperature model of pavement in natural environment is constructed by the heat con- duction equation.Then,the air temperature and the solar radiation intensity in a user-defined time interval are fit- ted by means of the interpolatory trigonometric polynomials on the basis of the discrete least-square approximation,on which a surface boundary condition is determined.Finally,with the aid of the Gaussian quadrature formula,the inverse Laplace transform is solved,so that the analytical solution of the time-dependent one-dimension temperature field can be easily derived.The results are compared with the measured temperatures in the asphalt overlay on the existing cement concrete pavement system in both summer and winter,it is found that the maximum error of the pre- dicted temperature and the measured one at different depth locations is less than 3℃,meaning that the proposed method achieves a high precision in predicting the pavement temperature field.