“…Under these circumstances, numerical codes such as the Codeon Multiflo for transient and steady-state well simulation (MULTIFLO from Lichtner et al, 2000, with application in Ferguson andWoodbury, 2005;Majorowicz et al, 2006) and VS2DH (Healy and Ronan, 1996, with application in Su et al, 2004;Essaid et al, 2008;Duque et al, 2010) facilitate the simulations of heat transport in the subsurface layer under variable boundary conditions. Previously, there have been no studies that have estimated the climate change impacts on aquifer thermal regimes by considering the dynamic nature of the boundary conditions and unsteady groundwater flow.…”