In permafrost regions, the thaw depth strongly controls shallow subsurface hydrologic processes that in turn dominate catchment runoff. In seasonally freezing soils, the maximum expected frost depth is an important geotechnical engineering design parameter. Thus, accurately calculating the depth of soil freezing or thawing is an important challenge in cold regions engineering and hydrology. The Stefan equation is a common approach for predicting the frost or thaw depth, but this equation assumes negligible soil heat capacity and thus exaggerates the rate of freezing or thawing. The Neumann equation, which accommodates sensible heat, is an alternative implicit equation for calculating freeze-thaw penetration. This study details the development of correction factors to improve the Stefan equation by accounting for the influence of the soil heat capacity and non-zero initial temperatures. The correction factors are first derived analytically via comparison to the Neumann solution, but the resultant equations are complex and implicit. Explicit equations are obtained by fitting polynomial functions to the analytical results. These simple correction factors are shown to significantly improve the performance of the Stefan equation for several hypothetical soil freezing and thawing scenarios. Copyright (c) 2015 John Wiley & Sons, Ltd.

Soil temperature regimes were studied in three ecosystems of the north of Western Siberia in the zone of isolated permafrost: the forest ecosystem with gleyic loamy sandy podzol (Stagnic Albic Podzol), the flat-topped peat mound ecosystem with humus-impregnated loamy sandy to light loamy peat cryozem (Histic Oxyaquic Turbic Cryosol (Arenic)), and the peat mound (palsa) ecosystem with oligotrophic destructive permafrost-affected peat soil (Cryic Histosol). Annual temperature measurements in the soil profiles demonstrated that these soils function under different temperature regimes: very cold permafrost regime and cold nonpermafrost regime. The following annual temperature characteristics proved to be informative for the studied soils: sums of above-zero temperatures at the depths of 10 and 20 cm, the maximum depth of penetration of temperatures above 10A degrees C, and the number of days with daily soil temperatures above (or below) 0A degrees C at the depth of 20 cm. On the studied territory, the insulating effect of the snow cover in winter was at least two times more pronounced than the insulating effect of the vegetation cover in summer. Cryogenic soils of the studied region are characterized by the high buffering towards changing climatic parameters. This is explained by the presence of the litter and peat horizons with a very low thermal diffusivity and by the presence of permafrost at a relatively shallow depth with temperature gradients preventing penetration of heat to the permafrost table.