Abstract: This paper aims to develop an accurate and efficient numerical model for three-dimensional transient thermal groundwater flow problems. The modified Richards equation and heat transport equation are considered to govern ...