Simulation of Soil Freezing and Thawing for Different Groundwater Table Depths

During freeze–thaw periods, the transformation between phreatic water and soil water will change the soil hydrothermal properties and affect the soil freezing and thawing in shallow groundwater areas. The purpose of this study was to determine the effect of four different groundwater table depths (GTDs) and two soil textures on the process of soil freezing and thawing during two successive freeze–thaw periods using the Simultaneous Heat and Water (SHAW) model. The results show that the frost depth was the maximum when the GTD was 1.0 m, and the maximum frost depths of sandy loam and fine sand were 97.6 and 98.9 cm, respectively. When the GTD was larger than 1.5 m, the maximum frost depth decreased with an increase in GTD, and the maximum frost depth of the soil profile was more sensitive to changes in the air temperature. The frost depth of the soil profile was linear with the square root of the accumulated negative soil surface temperature (ANST) under different GTDs. The ANST was influenced by the phreatic evaporation, and the soil freezing rate increased with an increase in GTD under the same ANST. This research is significant for the rational development of soil water and heat resources and the study of soil water–heat transfer in shallow groundwater areas.