Quantifying Temperature-sliding Inconsistency in Thermomechanical Coupling: A Comparative Analysis of Geothermal Heat Flux Datasets at Totten Glacier

Abstract. Rapid sliding of ice sheets requires warm basal temperatures and lubricating basal meltwater, whereas slow velocities typically correlate with a frozen bed. However, ice sheet models often infer basal sliding by inverting surface velocity observations with the vertical structure of temperature and hence rheology held constant. If the inversion is allowed to freely vary sliding over the model domain, then inconsistencies between the basal thermal state and ice motion can arise lowering simulation realism. In this study, we propose a new method that quantifies inconsistencies when inferring warm and cold-bedded regions of ice sheets. This method can be used to evaluate the quality of ice sheet simulation results without requiring any englacial or subglacial measurements. We apply the method to evaluate simulation results for Totten Glacier using an isotropic 3D full-Stokes ice sheet model with eight geothermal heat flux (GHF) datasets and compare our evaluation results with inferences on basal thermal state from radar specularity. The rankings of GHF datasets based on inconsistency are closely aligned with those using the independent specularity content data. Examples of the method utility are 1. an inconsistency characterizing overcooling with all GHFs near the western boundary of Totten Glacier between 70° S–72° S, where there is a bedrock canyon and fast surface ice velocities, which suggests that GHF is low in all published datasets; 2. an overheating inconsistency in the eastern Totten Glacier with all GHFs that leads to overestimation of ice temperature due, in this case, to an unrealistically warm surface temperature. Our approach opens a new avenue for assessing the self-consistency and reliability of ice sheet model results and GHF datasets, which may be widely applicable.