Abstract: Deep groundwater recharge is a complex process involving coupled thermal-hydrodynamic-mechanical (THM) interactions. Focusing on a deep groundwater recharge test site in the Yunhe District, this study developed a THM-coupled numerical model that integrates wellbore-reservoir flow dynamics, thermal effects, and formation deformation. The model was validated and calibrated using measured water level data. We systematically analyzed the spatiotemporal distribution of reservoir pressure and temperature during recharge and predicted the spatiotemporal evolution of soil deformation induced by groundwater recharge under given mechanical parameters. Results indicate that during recharge processes, pressure propagates rapidly through the reservoir in the form of pressure waves. The extent of thermal influence expands with recharge duration, reaching a maximum horizontal distance of 23.8 m by the end of the recharge period. Recharge induces uplift in the overlying strata and compression settlement in the underlying strata. Under the given formation parameters, the maximum uplift within the formation is 0.036 m, while the maximum surface uplift is 0.016 m, and the maximum compression settlement in the underlying strata is 0.015 m. The established model effectively characterizes the multi-field coupled response in deep groundwater recharge, providing a scientific basis for optimizing recharge strategies and controlling land subsidence.