Mathematical modeling of temporal changes in snow-firn properties in the cold season

Abstract

The paper describes a non-linear dynamic thermomechanical model of the snow–firn development in the high mountain environment during the cold season. The model allows the estimation of the thermal and mechanical development of snow in an inhomogeneous stratified snowpack and in the upper part of the firn layer, the development of the snow–firn structure and transformation of snow into firn. The non-linear mathematical model consists of a heat problem, water-vapor diffusion problem, densification of the snow–firn layer, structural equations and strength conditions of the metamorphosed snow–firn depending on its temperature, structural parameters, density and rate of sublimation. The model includes conditions of snow transformation into firn and its structural changes. These equations are non-linear and coupled. The iterative finite-difference numerical method was used for the calculations and involved a full mathematical model. Numerical experiments were implemented using the complete model as well as real. The numerical computations were made for various meteorological data, temperature gradients and in the wide range of structural parameters. The calculation results were found to be in good agreement with the results of observations.