Shallow geothermal potential needs to be evaluated when introducing the ground source heat pump (GSHP) system. Accurate evaluation of shallow geothermal potential can be done by estimating underground thermal properties with precision. The data previously established from measuring underground thermal properties, however, have limitations in that the locations where the measurements were made were not identical. Therefore, this study aimed to develop a hybrid model for estimating UGT using the finite element method (FEM) and the kriging method. This study was conducted in three steps: (i) data collection; (ii) estimation of the UGT in 216 non -measurement locations using FEM and the kriging method; and (iii) estimation of the shallow geothermal potential using the geothermal potential (G.POT) method. The average mean absolute percentage error of the hybrid model (3.67%) was lower than those of the other estimation models using FEM and the kriging method (FEM: 4.62%; Kriging method: 6.34%). Thus, the analysis showed that the developed hybrid model was the most suitable method for estimating UGT in all the length ranges (65-250 m) where the ground heat exchanger was to be installed. The hybrid model can help policy makers evaluate shallow geothermal potential at the national level. (C) 2017 Elsevier B.V. All rights reserved.