Relations between basal condition, subglacial hydrological networks and geothermal flux in Antarctica

When modelling the Antarctic ice sheet, the velocity of the ice flow is linked to its temperature. Depending on the thermal rate, the flow rate may vary between deformation and sliding. In this study, we focus on the geothermal flux because it is the least well-known component of the heat equation, and because it constrains the temperature at the bottom of the ice sheet. We used available geological data to build a map of the geothermal flux, which was found to increase from 51 mW/m2 in East Antarctica to 68 mW/m2 in West Antarctica. These values were integrated in the computation of a basal temperature map. The available map of hydrological networks clearly shows more melted areas in West Antarctica than in the earlier results. So we suggest that the model should be forced with higher geothermal flux values, over 85 mW/m2 in this sector. This increase is in good agreement with published results which found a geothermal flux three times higher in West Antarctica. Finally, we computed the bottom melt rate over the ice sheet area which has a mean value of 3.5 mm/yr resulting in a lost of melted ice equal to 1% of the total mass balance.