A thermomechanical model for frost heave and subglacial frozen fringe

Abstract

Ice-infiltrated sediment, or frozen fringe, is responsible for phenomena such as frost heave, ice lenses and metres of debris-rich ice under glaciers. Understanding frozen fringes is important as frost heave is responsible for damaging infrastructure at high latitudes and sediment freeze-on at the base of glaciers can modulate subglacial friction, influencing the rate of global sea level rise. Here we describe the thermomechanics of liquid water flow and freezing in ice-saturated sediments, focusing on the conditions relevant for subglacial environments. The force balance that governs the frozen fringe thickness depends on the weight of the overlying material, the thermomolecular force between ice and sediments across liquid premelted films and the water pressure required by Darcy flow. We combine this mechanical model with an enthalpy method that conserves energy across phase change interfaces on a fixed computational grid. The force balance and enthalpy model together determine the evolution of the frozen fringe thickness and our simulations predict frost heave rates and ice lens spacing. Our model accounts for premelting at ice–sediment contacts, partial ice saturation of the pore space, water flow through the fringe, the thermodynamics of the ice–water–sediment interface and vertical force balance. We explicitly account for the formation of ice lenses, regions of pure ice that cleave the fringe at the depth where the interparticle force vanishes. Our model results allow us to predict the thickness of a frozen fringe and the spacing of ice lenses in subaerial and subglacial sediments.