A Semianalytical Interface Model of Soil Freeze/Thaw and Permafrost Evolution

A physically based one-dimensional sharp-interface model of active layer evolution and permafrost thaw is presented. This computationally efficient, semianalytical, nonequilibrium solution to soil freeze-thaw problems in partially saturated media is proposed as a component of hydrological models to describe seasonal ground ice, active layer evolution, and changes in permafrost temperature and extent. The model is developed and validated against the analytical Stefan solution and a finite volume coupled heat and mass transfer model of freeze-thaw in unsaturated porous media. Unlike analytic models, the interface model provides a nonequilibrium solution to the heat equation while permitting a wide range of temporally variable boundary conditions and supporting the simulation of multiple interfaces between frozen and unfrozen soils. The model is implemented for use in discontinuous permafrost peatlands where soil properties are highly dependent on soil ice content and infiltration capacity is high. It is demonstrated that the model is suitable for the representation of variably saturated active layer and permafrost evolution in cases both with and without a talik.