The thermal coupling between the atmosphere and the subsurface on the Qinghai-Tibetan Plateau (QTP) governs permafrost stability, surface energy balance, and ecosystem processes, yet its spatiotemporal dynamics under accelerated warming are poorly understood. This study quantifies soil-atmosphere thermal coupling ((3) at the critical 0.1 m root-zone depth using in-situ data from 99 sites (1980-2020) and a machine learning framework. Results show significantly weaker coupling in permafrost (PF) zones (mean (3 = 0.42) than in seasonal frost (SF) zones (mean (3 = 0.50), confirming the powerful thermal buffering of permafrost. Critically, we find a widespread trend of weakening coupling (decreasing (3) at 66.7 % of sites, a phenomenon most pronounced in SF zones. Our driver analysis reveals that the spatial patterns of (3 are primarily controlled by surface insulation from summer rainfall and soil moisture. The temporal trends, however, are driven by a complex and counter-intuitive interplay. Paradoxically, rapid atmospheric warming is the strongest driver of a strengthening of coupling, likely due to the loss of insulative snow cover, while trends toward wetter conditions drive a weakening of coupling by enhancing surface insulation. Spatially explicit maps derived from our models pinpoint hotspots of accelerated decoupling in the eastern and southern QTP, while also identifying high-elevation PF regions where coupling is strengthening, signaling a loss of protective insulation and increased vulnerability to degradation. These findings highlight a dynamic and non-uniform response of land-atmosphere interactions to climate change, with profound implications for the QTP's cryosphere, hydrology, and ecosystems.

Under the action of freeze-thaw cycles, the internal temperature and water distribution of slope soils in cold regions change significantly, which directly affects the stability of slopes. In order to study the differences in hydrothermal reactions at different depths and their impacts on the stability of slopes. This study establishes both a freeze-thaw model and a hydrothermal coupling model, combining field measurements with numerical simulations to examine the dynamic changes in hydrothermal characteristics within the slope. The results indicate that the variation in slope temperature with depth can be divided into three stages: initial freezing, stable freezing, and thawing. In the freezing stage, the negative temperature gradient drives water to migrate towards the freezing front, forming segregated ice and inducing frost heave. In the thawing stage, the latent heat released by the phase change in segregated ice promotes water to move towards the slope toe, increasing the water content there and indirectly exacerbating the risk of slope instability. The heat and moisture transfer in frozen soil slopes shows non-linear and dynamic characteristics. The unique process of one-way freezing and two-way thawing makes the thawing rate 1.35 times that of the freezing rate, and this asymmetric characteristic is the key to understanding the mechanism of slope instability.

Study region: The source area of the Yangtze River, a typical catchment in the cryosphere on the Tibet Plateau, was used to develop and validate a distributed hydrothermal coupling model. Study focus: Climate change has caused significant changes in hydrological processes in the cryosphere, and related research has become hot topic. The source area of the Yangtze River (SAYR) is a key catchment for studies of hydrological processes in the cryosphere, which contains widespread glacier, snow, and permafrost. However, the current hydrological modeling of the SAYR rarely depicts the process of glacier/snow and permafrost runoff from the perspective of coupled water and heat transfer, resulting in distortion of simulations of hydrological processes. Therefore, we developed a distributed hydrothermal coupling model, namely WEP-SAYR, based on the WEP-L (Water and energy transfer process in large river basins) model by introducing modules for glacier and snow melt and permafrost freezing and thawing. New hydrological insights for the region: In the WEP-SAYR model, the soil hydrothermal transfer equations were improved, and a freezing point equation for permafrost was introduced. In addition, the glacier and snow meltwater processes were described using the temperature index model. Compared to previously applied models, the WEP-SAYR portrays in more detail glacier/ snow melting, dynamic changes in permafrost water and heat coupling, and runoff dynamics, with physically meaningful and easily accessible model parameters. The model can describe the soil temperature and moisture changes in soil layers at different depths from 0 to 140 cm. Moreover, the model has a good accuracy in simulating the daily/monthly runoff and evaporation. The Nash-Sutcliffe efficiency exceeded 0.75, and the relative error was controlled within +/- 20 %. The results showed that the WEP-SAYR model balances the efficiency of hydrological simulation in large scale catchments and the accurate portrayal of the cryosphere elements, which provides a reference for hydrological analysis of other catchments in the cryosphere.

Accurately quantifying the impact of permafrost degradation and soil freeze-thaw cycles on hydrological processes while minimizing the reliance on observational data are challenging issues in hydrological modeling in cold regions. In this study, we developed a modular distributed hydro-thermal coupled hydrological model for cold regions (DHTC) that features a flexible structure. The DHTC model couples heat-water transport processes by employing the conduction-advection heat transport equation and Richard equation considering ice-water phase change. Additionally, the DHTC model integrates the influence of organic matter into the hydrothermal parameterization scheme and includes a subpermafrost module based on the flow duration curve analysis to estimate cold-season streamflow sustained by subpermafrost groundwater. Moreover, we incorporated energy consumption due to ice phase changes to the available energy, enhancing the accuracy of evaporation estimation in cold regions. A comprehensive evaluation of the DHTC model was conducted. At the point scale, the DHTC model accurately replicates daily soil temperature and moisture dynamics at various depths, achieving average R-2 of 0.98 and 0.87, and average RMSE of 0.61degree celsius and 0.03 m(3)m(-3), respectively. At the basin scale, DHTC outperformed (Daily: R-2 = 0.66, RMSE = 0.75 mm; Monthly: R-2 = 0.90, RMSE = 15.7 mm) the GLDAS/FLDAS Noah, GLDAS/VIC, and PML-V2 models in evapotranspiration simulation. The DHTC model also demonstrated reasonable performance in simulating daily (NSE = 0.70, KGE = 0.84), monthly (NSE = 0.86, KGE = 0.90), and multi-year monthly (NSE = 0.97, KGE = 0.93) streamflow in the Source Regions of Yangtze River. DHTC also successfully reproduced the snow depth in basin-averaged time series and spatial distributions (RMSE = 0.86 cm). The DHTC model provides a robust tool for exploring the interactions between permafrost and hydrological processes, and their responses to climate change.

Considering the impact of the subgrade water level and freeze-thaw cycles, experiments were conducted on ballast track subgrade mud pumping. The study analyzed the migration of water and fine particles, as well as the characteristics of mud formation during the mud pumping process of the ballast track subgrade under cyclic loading. The research findings indicate that, during the initial loading stage at ambient temperature, moisture migrates upwards from the bottom. As dynamic loading is continuously applied, the internal pore water pressure in the subgrade soil gradually dissipates, resulting in a decrease in the pore water pressure gradient and a stabilization of the moisture content in each soil layer. When the water level is positioned in the middle of the subgrade, the upper soil is in an unsaturated state with a relatively low volumetric water content of approximately 26%. Fine particle migration does not occur, and the effective stress at the subgrade surface is much greater than zero, thus preventing mud pumping. When the water level is at the top of the subgrade, particle migration is more pronounced. The effective stress at the subgrade surface rapidly decreases to below 0 under the action of the load, resulting in mud pumping phenomena. Compared to unidirectional freezing, freeze-thaw loading results in a slower descent rate of the freezing front and a greater amount of moisture migration. Under thawing conditions, the upper soil layer of the subgrade melts before the lower soil layer, forming a frozen soil interlayer. Due to the isolation effect of the frozen soil interlayer, the upper soil layer retains a higher moisture content. Under the action of the load, the effective stress at the subgrade surface rapidly develops into negative values, making it more susceptible to mud pumping.

Introduction: Permafrost and seasonally frozen soil are widely distributed on the Qinghai-Tibetan Plateau, and the freezing-thawing cycle can lead to frequent phase changes in soil water, which can have important impacts on ecosystems.Methods: To understand the process of soil freezing-thawing and to lay the foundation for grassland ecosystems to cope with complex climate change, this study analyzed and investigated the hydrothermal data of Xainza Station on the Northern Tibet from November 2019 to October 2021.Results and Discussion: The results showed that the fluctuation of soil temperature showed a cyclical variation similar to a sine (cosine) curve; the deep soil temperature change was not as drastic as that of the shallow soil, and the shallow soil had the largest monthly mean temperature in September and the smallest monthly mean temperature in January. The soil water content curve was U-shaped; with increased soil depth, the maximum and minimum values of soil water content had a certain lag compared to that of the shallow soil. The daily freezing-thawing of the soil lasted 179 and 198 days and the freezing-thawing process can be roughly divided into the initial freezing period (November), the stable freezing period (December-early February), the early ablation period (mid-February to March), and the later ablation period (March-end of April), except for the latter period when the average temperature of the soil increased with the increase in depth. The trend of water content change with depth at all stages of freezing-thawing was consistent, and negative soil temperature was one of the key factors affecting soil moisture. This study is important for further understanding of hydrothermal coupling and the mechanism of the soil freezing-thawing process.

Hydro -thermal coupling is the essence of the freeze -thaw process, and theoretical studies of this coupled process have been hot topics in the field of frozen soil. Darcy's law of unsaturated soil water flow, heat conduction theory, and relative saturation and solid -liquid ratio are based on this paper. According to the principle that the cumulative curve of particle gradation of canal foundation soil is similar to soil -water properties. A soil -water characteristic curve is derived using the cumulative particle gradation curve. VG model is then used to fit soil -water characteristic curves to obtain the canal foundation soil's hydraulic characteristic parameters, and the established hydro -thermal coupling model is modified to reflect canal foundation soil hydro -thermal evolution more objectively. A closed system one-way freezing test method is used to verify the feasibility of the proposed method in this part. The results show that the optimal parameters of the VG model of the subsoil are a = 0.06, n = 1.2, and m = 0.17, and the temperature and water fields obtained from the simulation are in good agreement with the measured data, showing the utility of the hydro -thermal coupling model in predicting hydraulic parameters. Analysis of the multi -field interaction mechanism and dynamic coupling process of the canal foundation soil during freezing and thawing. This has great importance for preventing freezing damage in canals and protecting agricultural safety.

Landslide damage to soil graben slopes in seasonal freezing zones is a crucial concern for highway slope safety, particularly in the northeast region of China where permafrost thawing is significant during the spring. The region has abundant seasonal permafrost and mostly comprises powdery clay soil that is susceptible to landslides due to persistent frost and thaw cycles. The collapse of a slope due to thawing and sliding not only disrupts highway operations but also generates lasting implications for environmental stability, economic resilience, and social well-being. By understanding and addressing the underlying mechanisms causing such events, we can directly contribute to the sustainable development of the region. Based on the Suihua-Beian highway graben slope landslide-management project, this paper establishes a three-dimensional finite element model of a seasonal permafrost slope using COMSOL Multiphysics 6.1 finite element numerical analysis software. Additionally, the PDE mathematical module of the software is redeveloped to perform hydrothermal-coupling calculations of seasonal permafrost slopes. The simulation results yielded the dynamic distribution characteristics of the temperature and seepage field on the slope during the F-T process. The mechanism behind the slope thawing and sliding was also unveiled. The findings provide crucial technical support for the rational analysis of slope stability, prevention of sliding, and effective control measures, establishing a direct linkage to the promotion of sustainable infrastructure development in the context of transportation and roadway engineering.

To provide reliable input information for the load design and extraction of lunar soil water ice samples, it is necessary to study the water content distribution and water migration of simulated lunar soil water ice samples. On this basis, the temperature field model and the hydrothermal coupling relationship are proposed. The temperature field model was constructed by combining energy conservation and Fourier's heat transfer law. The coupling relationship was established, and the hydrothermal coupling model was obtained by testing the unfrozen water content using the nuclear magnetic resonance method. Finite element software was used to solve the model numerically, and the water migration rule of the soil water ice samples at different ambient temperatures were analyzed. Thin-wall drilling tests were carried out on the simulated lunar soil water ice samples to obtain water content data for different locations, and the simulation results were verified. Due to the migration effect of the cold end of the water, the closer we tested to the edge of the sample, the higher the water content was. The higher the ambient temperature was, the more pronounced the water migration phenomenon of the whole sample was. These research results provide a basis for sampling scheme design.