Abstract:Subglacial lakes are isolated, cold-temperature and high-pressure water environments hidden under ice sheets, which might host extreme microorganisms. Here, we use two-dimensional direct numerical simulations in order to investigate the characteristic temperature fluctuations and velocities in freshwater subglacial lakes as functions of the ice overburden pressure, $p_i$, the water depth, $h$, and the geothermal flux, $F$. Geothermal heating is the unique forcing mechanism as we consider a flat ice-water interface. Subglacial lakes are fully convective when $p_i$ is larger than the critical pressure $p_*\approx 2848$ dbar, but self organize into a lower convective bulk and an upper stably-stratified layer when $p_i < p_*$, because of the existence at low pressure of a density maximum at temperature $T_d$ greater than the freezing temperature $T_f$. For both high and low $p_i$, we demonstrate that the Nusselt number $Nu$ and Reynolds number $Re$ satisfy classical scaling laws provided that an effective Rayleigh number $Ra_{eff}$ is considered. We show that the convective and stably-stratified layers at low pressure are dynamically decoupled at leading order because plume penetration is weak and induces limited entrainment of the stable fluid. From the empirical power law equation for $Nu$ with $Ra_{eff}$, we derive two sets of closed-form expressions for the variables of interest, including the unknown bottom temperature, in terms of the problem parameters $p_i$, $h$ and $F$. The two predictions correspond to two limiting regimes obtained when the effective thermal expansion coefficient is either approximately constant or linearly proportional to the temperature difference driving the convection.

An exact solution to the temperature equation in a column of ice and bedrock

Topological rigidity of dimensional reduction to three dimensions

Precise Solution to 3D Coupled Thermohydromechanical Problems of Layered Transversely Isotropic Saturated Porous Media

A Three-Dimensional Fictitious Domain Method for Direct Numerical Simulations of Particle-Laden Flows with Heat Transfer

An Asymptotic Local Solution for Axi-Symmetric Heat Conduction Problems with Temperature Discontinuity on the Boundary

The Cooling Box Problem: Convection with a quadratic equation of state

Turbulent Convection in Subglacial Lakes

A Numerical Model for Water and Heat Transport in Freezing Soils with Nonequilibrium Ice‐water Interfaces

Closed Solutions for Transient Heat Transport in Geological Media: New Development, Comparisons, and Validations

Equilibrium states of the ice-water front in a differentially heated rectangular cell

Application of the heat equation to the study of underground temperature

Modeling the Link Between Air Convection and the Occurrence of Short‐Term Permafrost in a Low‐Altitude Cold Talus Slope

On mass conservation for ice sheets

Analytical solution for surrounding rock temperature of cold-region tunnel considering phase change

Formulation of approximate analytical solutions for a heat transfer model for an ice layer during microwave heating

Turbulent heat exchange between water and ice at an evolving ice-water interface

Freezing of Water in a Slab with Boundary Conditions of the Third Kind

Homogenization of the Stefan problem, with application to maple sap exudation

A Study of the One-Dimensional Heat-Conduction Equation with Radiation

Explicit and exact solutions concerning the Antarctic Circumpolar Current with variable density in spherical coordinates

Absence of superheating for ice Ih with a free surface : a new method of determining the melting point of different water models