“…For all subglacial sites, heat flow is calculated using a slightly modified version of the temperature gradient approach. We use ice temperature profiles from the Greenland Ice Borehole Temperature Profiles database (Mankoff, 2021; https://github.com/GEUS-PROMICE/greenland_ice_ borehole_temperature_profiles, last access: 8 February 2022). At sites where the ice temperature profile is complete, meaning it spans the full ice thickness from surface to bed (Devon73, PrinceWales05, Meighen67, Agas-siz77/79A/79B/84, Penny96, CampCentury, DYE-3, GISP2, NEEM, HansTausen_Dome, HansTausen_Hare), we fit a second-order temperature-depth function in the bottom 10 % of the borehole (Paterson, 1968;Weertman, 1968;Paterson et al, 1977;Fisher and Koerner, 1984;Gundestrup and Hansen, 1984;Clarke et al, 1987;Fisher et al, 1988;Clausen et al, 2001;Kinnard et al, 2008;Rasmussen et al, 2013;Mac-Gregor et al, 2015).…”