Thawing-triggered slope failures and landslides are becoming an increasing concern in cold regions due to the ongoing climate change. Predicting and understanding the behaviour of frozen soils under these changing conditions is therefore critical and has led to a growing interest in the research community. To address this challenge, we present the first mesh-free smoothed particle hydrodynamics (SPH) computational framework designed to handle the multi-phase and multi-physic coupled thermo-hydro-mechanical (THM) process in frozen soils, namely the THM-SPH computational framework. The frozen soil is considered a tri-phase mixture (i.e., soil, water and ice), whose governing equations are then established based on u-p-T formulations. A critical-state elasto-plastic Clay and Sand Model for Frozen soils (CASM-F), formulated in terms of solid-phase stress, is then introduced to describe the transition response and large deformation behaviour of frozen soils due to thawing action for the first time. Several numerical verifications and demonstrations highlight the usefulness of this advanced THM-SPH computational framework in addressing challenging problems involving thawing-induced large deformation and failures of slopes. The results indicate that our proposed single-layer, fully coupled THM-SPH model can predict the entire failure process of thawing-induced landslides, from the initiation to post-failure responses, capturing the complex interaction among multiple coupled phases. This represents a significant advancement in the numerical modelling of frozen soils and their thawing-induced failure mechanisms in cold regions.

Frozen soil, covering most of the Tibetan Plateau (TP), critically influences land surface and climate simulations. Although some studies have made advancements in simulations, further investigation into the distinct mechanisms underlying relevant parameterization schemes remains essential. This study compares two frozen soil permeability schemes in Noah-MP (NY06: high-permeability; Koren99: low-permeability) to elucidate their distinct hydrological mechanisms. Although significant disparities exist in the simulation of soil water and ice content between the two schemes in permafrost regions, the simulated soil water content in the shallow layer exhibits similarity. Their underlying physical processes behind this similarity differ fundamentally: Koren99 relies on cross-seasonal ice melt recharge, whereas NY06 depends more on current-season precipitation and snowmelt. With greater soil depth, soil water differences progressively propagate downward, amplifying variations in hydraulic conductivity, and soil memory effects become increasingly dominant. Meanwhile, the Koren99 scheme more effectively impedes bottom-up melting water transport than top-down effect. However, the aforementioned disparities are not apparent in seasonally frozen soil. Notable disparities also exist in simulated evapotranspiration and surface runoff over permafrost regions, particularly during the summer months. This research investigates the differences in water transport within frozen soil over the TP, elucidates the distinct hydrological mechanisms underlying different frozen soil permeability schemes, and highlights that similar soil hydrothermal simulations are associated with different physical processes, leading to varying degrees of effectiveness in soil memory. Furthermore, this research elucidates the dual role of soil ice (permeability restriction and water storage) in hydrological processes, providing a theoretical basis for improving frozen soil parameterization.

The hysteresis effect of unfrozen water during freeze-thaw cycles greatly influences the hydrothermal properties of soil. To better understand the hysteresis behavior of unfrozen water in the soil, this study utilized frequency domain reflectometry to measure the unfrozen water content variations in silty clay under both stepwise and rapid temperature change modes. The hysteresis effect of unfrozen water in soil was analyzed, also the underlying mechanism was revealed. The results indicate that unfrozen water content variations are consistent across the two temperature change modes, with hysteresis observed in both scenarios. This effect was more noticeable during the rapid temperature change mode, and soil samples with higher initial moisture content froze earlier and thawed more slowly in this mode. The hysteresis phenomena are mainly influenced by the ice crystal metastable nucleation, the blockage effect of pore ice crystallization, and the pore water pressure changes during phase transition. The main cause of unfrozen water hysteresis in soil during the initial freezing phase is the metastable nucleation process. In the later stages of freezing, the hysteresis effect is primarily driven by changes in capillary water curvature, induced by the blockage effect of pore ice crystallization, and shifts in pore water pressure during the ice-water phase transition. Also, a hysteresis model was proposed and validated against experimental data and existing models, demonstrating good performance and accurately predicting unfrozen water content under varying temperature conditions. This research enhances the understanding of the mechanism responsible for the hysteresis effect of unfrozen water content in frozen soil.

A set of direct shear tests on the soil-geotextile interface (SGI) were conducted using a temperature-controlled constant normal stiffness (CNS) direct shear apparatus. This was done in order to evaluate the effects of normal stiffness, initial normal stress, soil water content, and temperature on SGI shear behavior and microdeformation patterns. The observations indicate that all shear stress-shear displacement curves demonstrate strain-hardening characteristics, with SGI cohesion and friction angle increasing at higher normal stiffness and lower temperatures. At freezing conditions, water content significantly affects the interface friction angle, while this effect is minimal at positive temperatures. Normal stress increases with higher water content, lower temperatures, and higher normal stiffness. Shear stress initially rises with normal stress before decreases, with a more pronounced rise under sub-zero conditions. Normal stress shrinkage shows a positive correlation with normal stiffness. Micro-deformation analysis of soil particles at the interface indicates significant strain localization within the shear band, which is less pronounced under sub-zero temperatures compared to positive temperatures. These patterns of normal displacement vary across analysis points within the shear band, with the macroscopic normal displacement reflecting a cumulative effect of these microscopic variations.

Freeze-thaw cycles in seasonally frozen soil affect the boundary conditions of aqueducts with pile foundations, consequently impacting their seismic performance. To explore the damage characteristics and seismic behaviour of aqueduct bent frames in such regions, a custom testing apparatus with an integrated cooling system was developed. Two 1/15 scale models of reinforced concrete aqueduct bent frames with pile foundations were constructed and subjected to pseudo-static testing under both unfrozen and frozen soil conditions. The findings revealed that ground soil freezing has minimal impact on the ultimate bearing capacity and energy dissipation of the bent frame-pile-soil system, but significantly enhances its initial stiffness. Additionally, the frozen soil layer exerts a stronger embedding effect on the pile cap, ensuring the stability of the pile foundation during earthquakes. However, under large seismic loads, aqueduct bent frames experience greater damage and residual deformation in frozen soil compared to unfrozen soil conditions. Therefore, the presence of a seasonally frozen soil layer somewhat compromises the seismic performance of aqueduct bent frames. Subsequently, a finite element model considering pile-soil interaction (PSI) and frozen soil hydro-thermal effects was developed for aqueduct bent frames and validated against experimental results. This provides an effective method for predicting their seismic behaviors in seasonally frozen soil regions. Furthermore, based on the seismic damage characteristics of aqueduct bent frame with pile foundations observed in pseudo-static tests, a novel selfadaptive aqueduct bent frame system was designed to mitigate the adverse effects of seasonally frozen soil layer on seismic performance. This system is rooted in the principle of balancing resistance with adaptability, rather than solely depending on resistance. The seismic performance of this innovative system was then discussed, providing valuable insights for future seismic design of reinforced concrete aqueduct bent frames with pile foundations in seasonally frozen soil regions.

In cold regions, the strength and deformation characteristics of frozen soil change over time, displaying different mechanical properties than those of conventional soils. This often results in issues such as ground settlement and deformation. To analyze the rheological characteristics of frozen soil in cold regions, this study conducted triaxial creep tests under various creep deviatoric stresses and established a corresponding Discrete Element Method (DEM) model to examine the micromechanical properties during the creep process of frozen clay. Additionally, the Burgers creep constitutive model was used to theoretically validate the creep deformation test curves. The research findings indicated that frozen clay primarily exhibited attenuated creep behavior. Under low confining pressure and relatively high creep deviatoric stress, non-attenuated creep was more likely to occur. The theoretical model demonstrated good fitting performance, indicating that the Burgers model could effectively describe and predict the creep deformation characteristics of frozen clay. Through discrete element numerical simulations, it was observed that with the increase in axial displacement, particle displacement mainly occurs at both ends of the specimen. Additionally, with the increase in creep deviatoric stress, the specimen exhibits different deformation characteristics, transitioning from volumetric contraction to expansion. At the same time, the vertical contact force chains gradually increase, the trend of particle sliding becomes more pronounced, and internal damage in the specimen progresses from the ends toward the middle.

The entrance of permafrost tunnels in cold regions is particularly vulnerable to frost damage caused by complex thermal-hydro-mechanical (THM) interactions in unsaturated frozen soils. The effects of temperaturedependent volumetric strain variations across different stratum materials on heat and moisture transport are often neglected in existing THM coupling models. In this study, a novel THM coupled model for unsaturated frozen soil integrating volumetric strain correction is proposed, which addresses bidirectional interactions between thermal-hydraulic processes and mechanical responses. The model was validated through laboratory experiments and subsequently applied to the analysis of the Yuximolegai Tunnel. The results indicate that distinct layered ice-water distribution patterns are formed in shallow permafrost under freeze-thaw cycles, driven by bidirectional freezing and water migration. Critical mechanical responses were observed, including a shift in maximum principal stress from the invert (1.40 MPa, frozen state) to the crown (5.76 MPa, thawed state), and periodic lining displacements (crown > invert > sidewalls). Frost damage risks are further quantified by the spatial-temporal zoning of ice-water content-sensitive regions. These findings advance unsaturated frozen soil modeling and provide theoretical guidance for frost-resistant tunnel design in cold regions.

Frozen soils exhibit unique mechanical behavior due to the coexistence of ice and unfrozen water, making experimental studies essential for engineering applications in cold regions. This review comprehensively examines laboratory investigations on frozen soils under static and dynamic loadings, including uniaxial and triaxial compression, creep, direct shear, and freeze-thaw (F-T) cycle tests. Key findings on stress-strain characteristics, failure mechanisms, and the effects of temperature and time are synthesized. Advancements in microstructural analysis techniques, such as computed tomography (CT), scanning electron microscopy (SEM), nuclear magnetic resonance (NMR), and mercury intrusion porosimetry (MIP), are also summarized to elucidate the internal structural evolution of frozen soils. While significant progress has been made, further efforts are needed to better replicate complex environmental and loading conditions and to fully understand the interactions between multiple influencing factors. Future research should focus on developing novel experimental techniques, establishing standardized testing protocols, and creating a comprehensive database to enhance data accessibility and advance frozen soil research. This review provides critical insights into frozen soil mechanics and supports validating constitutive models and numerical simulations, aiding infrastructure design and construction in cold regions.

The service performance of frozen soil is one of the important factors that needs to be considered in designing and assessing the safety of artificial ground freezing projects. We conducted shear tests on ice-containing frozen soil and assessed soil performance and damage characteristics of the ice-frozen soil interface. On the basis of experimental results, we further investigated the damage of ice-containing frozen soil numerically using the finite-discrete element method. Experimental and numerical results show that temperature, the normal load, and moisture content are the primary factors influencing the mechanical properties of the ice-frozen soil interface. The effects of these parameters on shear strength, shear modulus, cohesion, and angle of internal friction were analyzed and discussed. There was a transition from ductile to brittle behavior at the ice-frozen soil interface with decreasing temperature. Transition occurred at higher temperatures in soils with higher moisture content. Because ice and sand differ in terms of stiffness, fractures appeared first at the ice-frozen sand interface. Under continued loading, the specific form of damage and maximum load-bearing capacity varied as a function of the location of the maximum shear stress zone and the ice in the soil. Our research findings provide valuable theoretical insights for the design and evaluation of the safety of artificial ground freezing engineering projects.

The cumulative plastic deformation and damage evolution of frozen soil-rock mixtures under cyclic loading was studied by a dynamic triaxial instrument with real-time resistivity measurement function. A series of low- temperature cyclic triaxial tests were conducted under varying confining pressures (200 kPa, 500 kPa, 800 kPa), block proportions (0, 30 %, 40 %, 50 %), and dynamic stress ratios (0.4, 0.6, 0.8). The results reveal that the cumulative plastic deformation process can be divided into three stages, such as microcrack closure as the initial stage, crack steady growth as the middle stage, and rapid crack propagation until it fails as the final stage. Under the same number of cycles, the greater the dynamic stress is, the greater the cumulative plastic deformation is. Furthermore, a strong correlation is identified between the resistivity and the cumulative plastic deformation. With the increase of the number of cycles, the cumulative plastic deformation leads to the accumulation of internal damage, and the resistivity gradually increases. Thus, a damage evolution model based on resistivity damage variables is proposed. The model demonstrated an average fitting accuracy of 97.36 % with the experimental data.