Water migration behavior is the main cause of engineering disasters in cold regions, making it essential to understand its mechanisms and the resulting mechanical characteristics for engineering protection. This study examined the water migration process during soil freezing through both experimental and numerical simulations, focusing on the key mechanical outcomes such as deformation and pore water pressure. Initially, a series of controlled unidirectional freezing experiments were performed on artificial kaolin soil under various freezing conditions to observe the water migration process. Subsequently, a numerical model of water migration was formulated by integrating the partial differential equations of heat and mass transfer. The model's boundary conditions and relevant parameters were derived from both the experimental processes and existing literature. The findings indicate that at lower clay water content, the experimental results align closely with those of the model. Conversely, at higher water content, the modeled results of frost heaving were less pronounced than the experimental outcomes, and the freezing front advanced more slowly. This discrepancy is attributed to the inability of unfrozen water to penetrate once ice lenses form, causing migrating water to accumulate and freeze at the warmest ice lens front. This results in a higher ice content in the freezing zone than predicted by the model, leading to more significant freezing expansion. Additionally, the experimental observations of pore water pressure under freeze-thaw conditions corresponded well with the trends and peaks projected by the simulation results.

Alpine permafrost environments are highly vulnerable and sensitive to changes in regional and global climate trends. Thawing and degradation of permafrost has numerous adverse environmental, economic, and societal impacts. Mathematical modeling and numerical simulations provide powerful tools for predicting the degree of degradation and evolution of subsurface permafrost as a result of global warming. A particularly significant characteristic of alpine environments is the high variability in their surface geometry which drives large lateral thermal and fluid fluxes along topographic gradients. The combination of these topography-driven fluxes and unsaturated ground makes alpine systems markedly different from Arctic permafrost environments and general geotechnical ground freezing applications, and therefore, alpine permafrost demands its own specialized modeling approaches. In this work, we present a multi-physics permafrost model tailored to subsurface processes of alpine regions. In particular, we resolve the ice-water phase transitions, unsaturated conditions, and capillary actions, and account for the impact of the evolving pore space through freezing and thawing processes. Moreover, the approach is multi-dimensional, and therefore, inherently resolves the topography-driven horizontal fluxes. Through numerical case studies based on the elevation profiles of the Zugspitze (DE) and the Matterhorn (CH), we show the strong influence of lateral fluxes in 2D on active layer dynamics and the distribution of permafrost.