In-situ observations from temperature gradient measurements are sparse and go along with several uncertainties like climatic changes or hydrothermal circulation (Burton-Johnson et al., 2020).To establish continent-wide heat flow models, one must refer to indirect methods using geophysical or geological data. The results differ immensely, for example, between magnetic and seismological data (e.g., An et al., 2015b;Martos et al., 2017), and the underlying assumptions cannot be easily combined (Lösing et al., 2020). While these different approaches usually take the respective sensitivity ranges into account, simplifications like a definition of laterally constant thermal parameters are not considered in the assessment. Recently, Shen et al. (2020) demonstrated how an updated seismic tomography model completely changed the estimated GHF and its uncertainties compared to earlier studies by Shapiro and Ritzwoller (2004). In ice-covered regions, different geophysical models show no consensus on magnitude and spatial distribution of heat flow (Rezvanbehbahani et al., 2019;Van Liefferinge, 2018).