Temperature is an integral part of soil quality in terms of moisture content; coupling between water and heat can render a soil fertile, and plays a role in water conservation. Although it is widely recognized that both water and heat transport are fundamental factors in the quantification of soil mass and energy balance, their computation is still limited in most models or practical applications in the root zone under non-isothermal conditions. This research was conducted to: (a) implement a fully coupled mathematical model that contains the full coupled process of soil water and heat transport with plants focused on the influence of temperature gradient on soil water redistribution and on the influence of change in soil water movement on soil heat flux transport; (b) verify the mathematical model with detailed field monitoring data; and (c) analyze the accuracy of the model. Results show the high accuracy of the model in predicting the actual changes in soil water content and temperature as a function of time and soil depth. Moreover, the model can accurately reflect changes in soil moisture and heat transfer in different periods. With only a few empirical parameters, the proposed model will serve as guide in the field of surface irrigation.