Sensitivity of the results of modeling of seasonal ground freezing to selection of parameterization of the snow cover thermal conductivity

The relationship between the results of calculations of the dynamics of the temperature regime of the in freezing and thawing soil profile with the heating effect of the snow cover is considered. To analyze this connection, two coupled models are used: the model of formation and degradation of snow cover in winter and the model of heat transfer and soil moisture transport in underlying vadoze zone profile. Parametrization of the influence of the snow cover, which at each calculated moment of time has the current average density and depth, on the dynamics of the temperatures of the soil profile is due to the use of its specific thermal resistance, which depends on its current depth and the thermal conductivity coefficient. The coefficient of thermal conductivity of the snow cover is related with its density using six different published empirical relationships. Modeling of heat transfer in freezing and thawing soil is carried out on the example of the field site for monitoring the thermal regime located on the territory of the Zvenigorod Biological Station of Moscow State University. It is shown that the well-known relationships give similar curves for the dynamics of the depth of seasonal freezing, including the degradation of the seasonal freezing layer in the spring period, with the same dynamics of the snow cover. However, the maximum penetration depth of the zero isotherm differs significantly for different snow conductivity-snow density relationships. The tested six relationships were divided into three groups. Minimal freezing is provided by the Sturm model and the effective medium model. The average and rather poorly differentiating freezing from each other is given by the Pavlov, Osokin et al. and Jordan relationships. The greatest value of the freezing depth is obtained with using Pavlov’s relationship with a temperature correction.