Abstract. Geothermal heat flow is an important boundary condition for ice sheets, affecting, for example, basal melt rates, but for ice-covered regions, we only have sparse heat flow observations with partly high uncertainty of up to 30 m W m−2. In this study, we first investigate the agreement between such pointwise heat flow observations and solid Earth models, applying a 1D steady-state approach to perform a statistical analysis for the entire Arctic region. We find that most of the continental heat flow observations have a high reliability and agreement to solid Earth models, except a few data points, such as, for example, the NGRIP (North Greenland Ice Core Project) point in central Greenland. For further testing, we perform a conditional simulation with focus on Greenland in which the local characteristics of heat flow structures can be considered. Simple kriging shows that including or excluding the less reliable NGRIP point has a large influence on the surrounding heat flow. The geostatistical analysis with the conditional simulation supports the assumption that NGRIP might not only be problematic for representing a regional feature but likely is an outlier. Basal melt estimates show that such a local spot of high heat flow results in local high basal melt rates but leads to less variation than existing geophysical models.