Abstract: Deep groundwater recharge is a complex process involving coupled thermal-hydrodynamic-mechanical (THM) interactions. Aiming at a deep groundwater recharge test site in the Yunhe District test area, this study established a THM coupling numerical model integrating wellbore-reservoir flow dynamics, thermal effects and formation deformation. The model was validated and calibrated using measured water level data. The spatiotemporal distribution characteristics of reservoir pressure and temperature during recharge were systematically analyzed, and the spatiotemporal evolution of soil deformation induced by groundwater recharge was predicted 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.