Supplementary Files for: Community heat flow recommendations: suitable basal boundary conditions for Greenland and Antarctica in ISMIP7

Geothermal heat flow (GHF) influences ice sheet thermal conditions, affecting ice flow by sliding and deformation. However, GHF distribution under polar ice sheets remains poorly constrained, with few direct borehole-derived estimates and large discrepancies between glaciological and geophysicalmodels caused by methodological differences and data limitations. As a result, many ice sheet models rely on uniform GHF estimates, ensemble averages or outdated fields that oversimplify reality. The choice of GHF product can lead to significantly different thermal conditions simulated at the ice-bed interface, which affects the projected evolution of ice sheets under climate warming. Therefore, we conducted an expert elicitation survey to identify the most suitable GHF fields for use as basal boundary conditions in ice sheet modelling, particularly for the Ice Sheet Modelling Intercomparison Project for CMIP7 (ISMIP7). GHF fields generally fall into three categories: (1) outdated due to improved data availability, (2) overly simplified parameterisations and (3) current and preferred. For GHF fields that rank highly in the survey, we discuss uncertainty, data dependency and provide guidance on their use in different applications. Finally, we recommend two Antarctic and one Greenlandic GHF field(s) for ISMIP7.