Vegetation greening across the Tibetan Plateau, a critical ecological response to climate warming and land-cover change, affects soil hydrothermal regimes, altering soil moisture (SM) and soil temperature (ST) dynamics. However, its effects on SM-ST coupling remain poorly understood. Using integrated field measurements from a vegetation-soil (V-S) network, reanalysis, and physics-based simulations, we quantify responses of SM, ST, and their coupling to vegetation changes across the Upper Brahmaputra (UB) basin, southern Tibetan Plateau. Results show that strong positive SM-ST correlations occur throughout 0-289 cm soil layers across the basin, consistent with the monsoon-driven co-occurrence of rainy and warm seasons. Spatially, SM-ST coupling strength exhibits pronounced spatial heterogeneity, demonstrating strongest coupling in central basin areas with weaker intensities in eastern and western regions. Overall, vegetation greening consistently induces soil warming and drying: as leaf area index (LAI) increases from 20 % to 180 % of its natural levels, SM (0-160 cm) declines by 15 % to 29 % due to enhanced evapotranspiration and root water uptake. Mean ST simultaneously increases by 1.4 +/- 0.9 degrees C. Crucially, sparsely vegetated regions sustain warming (1.4-2.1 degrees C), while densely vegetated areas transition from initial warming to gradual cooling. These findings advance our understanding of soil hydrothermal dynamics and their broader environmental impacts, improving climate model parameterizations and informing sustainable land management strategies in high-altitude ecosystems.

Thawing permafrost alters climate not only through carbon emissions but also via energy-water feedback and atmospheric teleconnections. This review focuses on the Tibetan Plateau, where strong freeze-thaw cycles, intense radiation, and complex snow-vegetation interactions constitute non-carbon climate responses. We synthesize recent evidence that links freeze-thaw cycles, ground heat flux dynamics, and soil moisture hysteresis to latent heat feedback, monsoon modulation, and planetary wave anomalies. Across these pathways, both observational and simulation studies reveal consistent signals of feedback amplification and nonlinear threshold behavior. However, most Earth system models underrepresent these processes due to simplifications in freezethaw processes, snow-soil-vegetation coupling, and cross-seasonal memory effects. We conclude by identifying priority processes to better simulate multi-scale cryosphere-climate feedback, especially under continued climate warming in high-altitude regions.

In order to investigate the frost-heaving characteristics of wintering foundation pits in the seasonal frozen ground area, an outdoor in-situ test of wintering foundation pits was carried out to study the changing rules of horizontal frost heave forces, vertical frost heave forces, vertical displacement, and horizontal displacement of the tops of the supporting piles under the effect of groundwater and natural winterization. Based on the monitoring condition data of the in-situ test and the data, a coupled numerical model integrating hydrothermal and mechanical interactions of the foundation pit, considering the groundwater level and phase change, was established and verified by numerical simulation. The research results show that in the silty clay-sandy soil strata with water replenishment conditions and the all-silty clay strata without water replenishment conditions, the horizontal frost heave force presents a distribution feature of being larger in the middle and smaller on both sides in the early stage of overwintering. With the extension of freezing time, the horizontal frost heave force distribution of silty clay-sand strata gradually changes from the initial form to the Z shape, while the all-silty clay strata maintain the original distribution characteristics unchanged. Meanwhile, the peak point of the horizontal frost heave force in the all-silty clay stratum will gradually shift downward during the overwintering process. This phenomenon corresponds to the stage when the horizontal displacement of the pile top enters a stable and fluctuating phase. Based on the monitoring conditions of the in-situ test, a numerical model of the hydro-thermo-mechanical coupling in the overwintering foundation pit was established, considering the effects of the groundwater level and ice-water phase change. The accuracy and reliability of the model were verified by comparison with the monitoring data of the in-situ test using FLAC3D finite element analysis software. The evolution of the horizontal frost heaving force of the overwintering foundation pit and the change rule of its distribution pattern under different groundwater level conditions are revealed. This research can provide a reference for the prevention of frost heave damage and safety design of foundation pit engineering in seasonal frozen soil areas.

Accurately modeling soil-fluid coupling under large deformations is critical for understanding and predicting phenomena such as slope failures, embankment collapses, and other geotechnical hazards. This topic has been studied for decades and remains challenging due to the nonlinear responses of geotechnical structures, which typically result from plastic yielding and finite deformation of the soil skeleton. In this work, we comprehensively summarize the theory involved in the soil-fluid coupling problem. Within a finite strain framework, we employ an elasto-plastic constitutive model with linear hardening to represent the solid skeleton and a nearly incompressible model for water. The water content influences the behavior of the solid skeleton by affecting its cohesion. The governing equations are discretized by material point method and two sets of material points are employed to independently represent solid skeleton and fluid, respectively. The proposed method is validated by comparing simulation results with experimental results for the impact of water on dry soil and wet soil. The capability of the method is further demonstrated through two cases: (1) the impact of a rigid body on saturated soil, causing water seepage, and (2) the filling of a ditch, which considers the erosion of the foundation. This work may provide a versatile tool for analyzing the dynamic responses of fluid and solid interactions, considering both mixing and separation phenomena.

With global warming and intensified rainfall, the heat and moisture transfer processes within frozen subgrades beneath asphalt pavements have become increasingly complex, posing risks to highway stability in cold regions. This study developed a multi-physics coupled indoor simulation system based on a typical asphalt highway structure on the Qinghai-Tibet Plateau to examine subgrade responses under solar radiation, wind, and rainfall. Results showed that rainfall shifted the dominant depth of moisture migration from 7 cm to 12 cm, with moisture at 2 cm and 7 cm increasing rapidly by 2.67 % and 1.58 %, respectively. A nonlinear decrease-increase pattern was observed at 12 cm due to the capillary barrier effect. Evaporative latent heat significantly suppressed surface warming, reducing the temperature rise at 2 cm by 59.2 %, and delayed heat transfer to deeper layers (reductions of 54.5 %-64.7 %). A cumulative heat flux prediction model, incorporating solar radiation, evaporation, convection, and surface wetting, showed high accuracy (R-2 = 0.981 and 0.952; relative errors: 4.1 % and 9.6 %). Sensitivity analysis identified the surface wetting rate coefficient (beta) and evaporation attenuation coefficient (gamma) as dominant factors (F > 6.0). These findings improve understanding of rainfall-induced thermal effects and offer guidance for climate-resilient road infrastructure in permafrost regions.

In urban subway construction, shield tunneling near pile groups is common, where additional loads may threaten existing structures. This study establishes multiple 3D nonlinear FDM models with fluid-solid coupling to investigate how tunnel-pile clearances (Hc) affect the mechanical response of low-cap pile groups (2 x2) during side-by-side twin tunneling in composite strata. The advanced CYSoil model, incorporating nonlinearity, strain path dependency, and small strain behavior, is employed to simulate soil response. Results show that tunneling induces up to a similar to 66.7 % reduction in pore water pressure, forming a funnel-shaped seepage pattern. As Hc increases from 0.8D to 2.6D, the low-pressure zone shifts from sidewalls to vault and invert, while maximum displacements reduce by up to 14.04 mm (lateral), 5.28 mm (transverse), and 19.68 mm (vertical). Axial force evolution in piles follows a three-stage decline, i.e., rapid, slow, and moderate, with peak shaft resistance concentrated near the tunnel axis. These findings aid in optimizing tunnel-pile configurations and mitigating geotechnical risks.

A novel thermo-hydro-mechanical-chemical (THMC) coupling model grounded in thermodynamic dissipation theory was established to unravel the intricate behavior of unsaturated sulfate-saline soils during cooling crystallization. The model quantifies energy transfer and dissipation during crystallization and introduces a method to calculate the amount of sulfate crystallization. It intricately captures the interdependencies between crystallization, pore water pressure, crystallization pressure and volumetric expansion, while also accounting for the dynamic feedback of latent heat from phase transitions on heat conduction. The reliability of the model was validated through experimental data. Numerical simulations explored the effects of cooling paths, thermal conductivity, initial salt content and initial porosity on the crystallization behavior and mechanical properties. The model provides theoretical support for optimizing the engineering design and facility maintenance of sulfatesaline soils.

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.

Three-dimensional numerical models are developed to investigate the anti-liquefaction of ordinary (OSCs) and geosynthetic-encased (GESCs) stone columns in sandy soil under sinusoidal loading using the fluid-solid coupling method. The validated models capture and compare the vertical and radial deformation, excess pore water pressure (EPWP), and vertical effective stress of OSC, GESC, and sandy soil. Furthermore, ten essential factors are selected to conduct the parametric study. Numerical results reveal that GESC is more suitable for improving sandy soil and resisting dynamic load considering the deformation and EPWP. The bulging deformation is no longer the primary reason for failure. The partial encasement (e.g., 1-2D, D = column diameter) and short floating and end-bearing GESCs (e.g., 1-2.5D) are not recommended for reinforcing the sandy soil. GESC is more sensitive to low-frequency and high-amplitude loads, with shear and bending, whereas displays a block movement under higher frequency and lower amplitude loading. The change in loading amplitude is more disadvantageous to GESC than loading frequency. GESC with a large diameter cannot effectively resist the dynamic loads.

For many solids, irreversible deformation is often accompanied by changes in the internal structure, impacting the reversible responses, a phenomenon termed elasto-plastic coupling. This coupling has been observed experimentally in various geomaterials, including clayey and sandy soils, as well as hard and soft rocks. Fabric anisotropy, which characterizes the internal structure, is a distinct feature of soils and significantly influences both reversible and irreversible behaviors. In this study, we adopted a coupling formulation based on the framework of anisotropic critical state theory (ACST) to describe the anisotropic elasto-plastic coupling response of soils. The formulation incorporates a deviatoric fabric tensor F, which consistently quantifies the internal structure of soils in both reversible and irreversible range, into a hyperelastic formulation and a plastic model, respectively. A novel evolution rule of F, defined based on the current stress ratio and plastic strain, is proposed, where the direction gradually aligns with the loading direction and the norm achieves different asymptotic values depending on the applied loading paths. This allows for the representation of evolved anisotropy effects on elasticity, dilatancy and strength simultaneously, providing a natural description of elasto-plastic coupling. Within this coupling framework, any anisotropic model within ACST can serve as the plastic platform for developing the elasto-plastic coupling models with anisotropic hyperelasticity. Herein, a bounding surface plastic model is utilized for illustration. The proposed model's performance is demonstrated by especially comparing simulated results to test data on evolving elastic stiffness ratios and overall elastoplastic responses under varying monotonic and cyclic loading conditions.