Temperature changes affect the nonlinear consolidation process in soils, and there is limited associated theoretical research. In this study, the governing equations for nonlinear consolidation and thermal conduction are developed, and a mathematical model for one-dimensional nonlinear thermal consolidation in saturated clay under the impeded drainage boundary is established, where the temperature-dependent compressibility and permeability are considered. Meanwhile, the finite-difference solutions for nonlinear consolidation and the analytical solutions for thermal conduction are obtained, respectively. Furthermore, the proposed model's reasonableness is verified by comparison with other theoretical models. Based on this, the impact of several factors on nonlinear thermal consolidation behaviors is investigated. With a rise in temperature increment (Delta T), the dissipation rate of excess pore water pressure (EPWP) accelerates in the later consolidation stage, and the final settlement becomes larger. In addition, the EPWP dissipation rate grows remarkably with an increasing impeded drainage boundary parameter (mu). In particular, the impeded drainage boundary can be degraded into a drainage boundary when the value of mu becomes large (e.g., mu = 100 m-1). Increasing preconsolidation pressure (pcR) results in a reduction in settlement, and the maximum values of EPWP decline with a rising linear loading time (tc). Overall, this study contributes to the accurate prediction of the nonlinear consolidation process taking the thermal effect into account. This study might provide a basis for the analysis of soil consolidation features in geotechnical projects that involve thermal effects, which have been growing in number in recent decades. For instance, an approach that combines thermal treatment with surcharge preloading has started to be employed for soft soil reinforcement. When this method is used in practical engineering projects, the changes in EPWP and settlement can be predicted based on the model developed in this study. In addition, this study reveals that increasing the temperature by 60 degrees C can lead to an increase in the final settlement of saturated clay by approximately 25%, compared with ambient temperature treatment. Besides, compacted clay, which typically serves as a bottom engineering barrier, might be subject to consolidation deformation due to varied temperatures and external loading. The model presented could contribute to properly predicting the porosity change in compacted clay under this scenario.