Ablation of sloping ice faces into polar seawater

Abstract

The effects of the slope of an ice–seawater interface on the mechanisms and rate of ablation of the ice by natural convection are examined using turbulence-resolving simulations. Solutions are obtained for ice slopes $\unicode[STIX]{x1D703}=2^{\circ }{-}90^{\circ }$, at a fixed ambient salinity and temperature, chosen to represent common Antarctic ocean conditions. For laminar boundary layers the ablation rate decreases with height, whereas in the turbulent regime the ablation rate is found to be height independent. The simulated laminar ablation rates scale with $(\sin \unicode[STIX]{x1D703})^{1/4}$, whereas in the turbulent regime it follows a $(\sin \unicode[STIX]{x1D703})^{2/3}$ scaling, both consistent with the theoretical predictions developed here. The reduction in the ablation rate with shallower slopes arises as a result of the development of stable density stratification beneath the ice face, which reduces turbulent buoyancy fluxes to the ice. The turbulent kinetic energy budget of the flow shows that, for very steep slopes, both buoyancy and shear production are drivers of turbulence, whereas for shallower slopes shear production becomes the dominant mechanism for sustaining turbulence in the convective boundary layer.