Analysis of the Main Factors Controlling the Formation of Subaerial Taliks, Using a One-Dimensional Mathematical Model. A Case Study for the Shestakovka River Basin (Central Yakutia)

Received March 9, 2023; revised September 7, 2023; accepted October 2, 2023

This study presents a mathematical model of heat transfer in a subaerial talik. The model is based on the concepts presented in classical works on permafrost, as well as on the results of geological and geophysical research carried out in the Shestakovka River basin (Central Yakutia). This model is based on the solution of the classical Stefan problem on the moving of the phase transition boundaries for a multilayer and multiphase medium. The solution was calculated on a unstructured mesh. When the phase boundaries move, thawed or frozen layers of soil are formed or wedged out. The layers include: snow cover, seasonally thawed soil, seasonally frozen and frozen sand deposits, as well as soil-vegetative layer. Published empirical relationships were used to calculate thermophysical coefficients, which are presented in this article. Simple variants of the model were considered to clarify the contribution of various factors to the process of formation and evolution of taliks. It has been established that the presence of snow cover and soil-vegetative layer have the most significant effect on the formation of taliks. Calculations show that taliks are formed in the first years of the modeled period, in the presence of snow and the absence of soil-vegetative layer. The soil-vegetative layer, depending on its composition and moisture content (ice content), can prevent the formation and development of taliks. The authors do not consider cases where shrubs contribute to snow accumulation. The humidity and iciness of the layer of sand sediments located in Central Yakutia have practically no effect on this process.