Analytical modelling techniques are frequently used for prediction of the impacts of borehole heat exchangers (BHEs) on the thermal conditions in the subsurface. They are appealing because they are easy to use, compact and computationally efficient. However, analytical formulations also cut down the complexity of the conditions in the field by limited resolution of variations in space and time, and by often simplified specifications of boundary conditions. In this presentation, we focus on the

more »

... efinition of the land surface boundary condition in line source solutions. As BHEs are thin elongated forms that can be approximated by a linear shape function, main attention is on the radial heat transfer, but less on the axial effects. Given the long operation time, however, the sharp physical boundary at the top may exert a substantial influence on the temperature evolution in the shallow ground and thus influence the performance of the ground source heat pump. We examine different formulations with Dirichlet and Cauchy type boundaries included in the line source equation, with and without horizontal groundwater flow. Both homogeneous as well as inhomogeneous implementations are presented, which offer new flexibilities for simulation of nonuniform land surface impacts such as associated with varying land use types. The presented approach also facilitates to account for the accelerated ground heat flux due to global warming and due to urbanization such as observed in so-called subsurface urban heat islands in many cities.