Understanding soil organic carbon (SOC) distribution and its environmental controls in permafrost regions is essential for achieving carbon neutrality and mitigating climate change. This study examines the spatial pattern of SOC and its drivers in the Headwater Area of the Yellow River (HAYR), northeastern Qinghai-Xizang Plateau (QXP), a region highly susceptible to permafrost degradation. Field investigations at topsoils of 86 sites over three summers (2021-2023) provided data on SOC, vegetation structure, and soil properties. Moreover, the spatial distribution of key permafrost parameters was simulated: temperature at the top of permafrost (TTOP), active layer thickness (ALT), and maximum seasonal freezing depth (MSFD) using the TTOP model and Stefan Equation. Results reveal a distinct latitudinal SOC gradient (high south, low north), primarily mediated by vegetation structure, soil properties, and permafrost parameters. Vegetation coverage and above-ground biomass showed positive correlation with SOC, while soil bulk density (SBD) exhibited a negative correlation. Climate warming trends resulted in increased ALT and TTOP. Random Forest analysis identified SBD as the most important predictor of SOC variability, which explains 38.20% of the variance, followed by ALT and vegetation coverage. These findings likely enhance the understanding of carbon storage controls in vulnerable alpine permafrost ecosystems and provide insights to mitigate carbon release under climate change.

Uneven displacement of permafrost has become a major concern in cold regions, particularly under repeated freezing-thawing cycles. This issue poses a significant geohazard, jeopardizing the safety of transportation infrastructure. Statistical analyses of thermal penetration suggest that the problem is likely to intensify as water erosion expands, with increasing occurrences of uneven displacement. To tackle the challenges related to mechanical behavior under cyclic loading, the New Geocell Soil System has been implemented to mitigate hydrothermal effects. Assessment results indicate that the New Geocell Soil System is stable and effective, offering advantages in controlling weak zones on connecting slopes and reducing uneven solar radiation. Consequently, the New Geocell Soil System provides valuable insights into the quality of embankments and ensures operational safety by maintaining displacement at an even level below 1.0 mm. The thermal gradient is positive, with displacement below 6 degrees C/m, serving as a framework for understanding the stability of the subgrade. This system also enhances stress and release the sealing phenomenon.

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.

This study highlights the results of a palaeoecological analysis conducted on five permafrost peatlands in the northern tundra subzone along the Barents Sea coast in the European Arctic zone. The depth of the peat cores that were sampled was approximately 2 m. The analysis combined data on the main physical and chemical soil properties, radiocarbon dating, botanical composition, and mass fraction of polycyclic aromatic hydrocarbons (PAHs). The concentrations of 16 PAHs in peat organic layers ranged from 140 to 254 ng/g, with an average of 182 ng/g. The peatlands studied were dominated by PAHs with a low molecular weight: naphthalene, phenanthrene, fluoranthene, pyrene, chrysene. The vertical distribution patterns of PAHs along the peat profile in the active layer and permafrost were determined. PAHs migrating down the active layer profile encounter the permafrost barrier and accumulate at the boundary between active layer and permafrost layer. The deep permafrost layers accumulate large amounts of PAHs and PAH derivatives, which are products of lignin conversion during the decomposition of grassy and woody vegetation during the Holocene climate optima. The total toxic equivalency concentration (TEQ) was calculated. Peatlands from the Barents Sea coast have low toxicity for carcinogenic PAHs throughout the profile. TEQ ranged from a minimum of 0.1 ng/g to a maximum of 13.5 ng/g in all peatlands investigated. For further potential use in Arctic/sub-Arctic environmental studies, PAH indicator ratios were estimated. In all investigated sections and peatland horizons, the most characteristic ratios indicate the petrogenic (natural) origin of PAHs.

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.

The Arctic has been warming much faster than the global average, known as Arctic amplification. The active layer is seasonally frozen in winter and thaws in summer. In the 2017 Arctic Boreal Vulnerability Experiment (ABoVE) airborne campaign, airborne L- and P- band synthetic aperture radar (SAR) was used to acquire a dataset of active layer thickness (ALT) and vertical soil moisture profile, at 30 m resolution for 51 swaths across the ABoVE domain. Using a thawing degree day (TDD) model, ALT=K root TDD, we estimated ALT along the ABoVE swaths employing the 2-m air temperature from ERA5. The coefficient (K) calibrated has an R2=0.9783. We also obtained an excellent fit between ALT and K root(TDD/theta) where theta is the soil moisture from ERA5 (R2=0.9719). Output based on shared-social economic pathway (SSP) climate scenarios SSP 1-2.6, SSP 2-4.5, and SSP 5-8.5 from seven global climate models (GCMs), statistically downscaled to 25-km resolution, was used to project the impacts of climate warming on ALT. Assuming ALT=K root TDD, the projections of UKESM1-0-LL GCM resulted in the largest projected ALT, up to about 0.7 m in 2080s under SSP5-8.5. Given that the mean observed ALT of the study sites is about 0.482 m, this implies that ALT will increase by 0.074 to 0.217 m (15% and 45%) in 2080s. This will have substantial impacts on Arctic infrastructure. The projected settlement Iset (cm) of 1 to 7 cm will also impact the infrastructure, especially by differential settlement due to the high spatial variability of ALT and soil moisture, given at local scale the actual thawing will partly depend on thaw sensitivity of the material and potential thaw strain, which could vary widely from location to location.

Climate change is transforming the ice-free areas of Antarctica, leading to rapid changes in terrestrial ecosystems. These areas represent <0.5% of the continent and coincide with the most anthropogenically pressured sites, where the human footprint is a source of contamination. Simultaneously, these are the locations where permafrost can be found, not being clear what might be the consequences following its degradation regarding trace element remobilisation. This raises the need for a better understanding of the natural geochemical values of Antarctic soils as well as the extent of human impact in the surroundings of scientific research stations. Permafrost thaw in the Western Antarctic Peninsula region and in the McMurdo Dry Valleys is the most likely to contribute to the remobilisation of toxic trace elements, whether as the result of anthropogenic contamination or due to the degradation of massive buried ice and ice-cemented permafrost. Site-specific locations across Antarctica, with abandoned infrastructure, also deserve attention by continuing to be a source of trace elements that later can be released, posing a threat to the environment. This comprehensive summary of trace element concentrations across the continent's soils enables the geographical systematisation of published results for a better comparison of the literature data. This review also includes the used analytical techniques and methods for trace element dissolution, important factors when reporting low concentrations. A new perspective in environmental monitoring is needed to investigate if trace element remobilisation upon permafrost thaw might be a tangible consequence of climate change.

This computational study focuses on the thermo-hydro-mechanical simulations of the behaviors of freezing soils used for artificial ground freezing (AGF) in a metro project. Leveraging the experimental and field data available in the literature, we simulate the sequential freezing and excavation of a twin tunneling that occurred in months during the actual construction of the tunnel. A thermo-hydro-mechanical model is developed to capture the multi-physical rate-dependent behaviors triggered by phase transitions, as well as the creeping and secondary consolidation of the soil skeleton and the ice crystals. We then calibrate the material models and establish the THM finite element model coupled with the rate-dependent multi-physical models, which may accurately predict the surface heave induced by ground freezing throughout the project. To showcase the potential of using simulations to guide the AGF, we simulate the scenario where a simultaneous freezing scheme is employed as an alternative to the actual sequential scheme design. We then compared the simulated performance with the recorded results obtained from the sequential scheme. Finally, parametric studies on the effect of ground temperature, the porosity of the frozen soil, and the intrinsic elastic modulus of the solid skeleton are conducted. The maximum surface heave is inferred from finite element simulations to quantify the sensitivity and the impact on the safety of AGF operations.

Small modular reactors (SMRs) are an alternative for clean energy solutions in Canada's remote northern communities, owing to their safety, flexibility, and reduced capital requirements. Currently, these communities are heavily reliant on fossil fuels, and the transition to cleaner energy sources, such as SMRs, becomes imperative for Canada to achieve its ambitious net-zero emissions target by 2050. However, applying SMR technology in permafrost regions affected by climate change presents unique challenges. The degradation of permafrost can lead to significant deformations and settlements, which can result in increased maintenance expenses and reduced structural resilience of SMR infrastructure. In this paper, we studied the combined effect of climate nonstationarity in terms of ground surface temperature and heat dissipation from SMR reactor cores for the first time in two distinct locations in Canada's North: Salluit in Quebec and Inuvik in the Northwest Territories. It was shown that these combined effects can make significant changes to the ground thermal conditions within a radius of 15-20 m around the reactor core. The change in the ground thermal conditions poses a threat to the integrity of the permafrost table. The implementation of mitigation strategies is imperative to maintain the structural integrity of the nuclear infrastructure in permafrost regions. The thermal modeling presented in this study paves the way for the development of advanced coupled thermo-hydromechanical models to examine the impact of SMRs and climate nonstationarity on permafrost degradation.

This study examines permafrost thermal regimes and hydrological responses to climate change in the Navarro Valley, Chile's Dry Central Andes, using a decade of ground temperature data (2013-2022) from two rock glaciers-RG1 (3805 m) and RG2 (4047 m)-alongside short-term meltwater conductivity records, meteorological data, and long-term streamflow records. We assess permafrost stability and climatic sensitivity by analyzing thermal offset data (2017-2022) and ground temperature trends. Both sites show sustained warming, but RG1 exhibits accelerated warming (+ 2.84 degrees C/decade), frequent freeze-thaw cycles, and extended thaw periods, indicating a transitional regime. In contrast, RG2 shows fewer freeze-thaw cycles and greater thermal buffering, consistent with cold permafrost. The statistical model overestimated thaw dynamics at RG2, highlighting the importance of field-based data for accurate classification. Hydrological signals at RG1-including cold, mineralized meltwater and rapid ground surface temperature stream coupling-are attributed to thawing and deeper flowpaths. Conductivity data (2014-2015) reveal solute pulses consistent with early melt events and debris interaction. Meanwhile, long-term streamflow trends indicate declining discharge. These findings suggest feedback between permafrost loss and water availability. This study underscores the divergent evolution of adjacent rock glaciers under warming by integrating thermal, hydrological, and climatic data. RG1 shows signs of degradation, while RG2 may act as a temporary refuge. Continued monitoring is essential for managing water security in vulnerable mountain regions like the Dry Andes.Graphical AbstractThis graphical abstract visually summarizes a ten-year monitoring effort of mountain permafrost and glacier hydrology in the Navarro Valley, Dry Andes (32 degrees S), with implications for water security under climate change. The left panel situates the study area within the upper Aconcagua Basin, identifying two instrumented sites within the Tres Gemelos rock glacier complex-RG1 (3805 m) and RG2 (4047 m)-and an automatic weather station. These sites were selected for continuous monitoring of ground temperature and streamflow to assess permafrost behavior in a water-stressed mountain catchment. Moving to the center, the image presents an integrated monitoring framework that links temperature-depth profiles, surface-subsurface thermal dynamics, and discharge records. Key indicators such as freeze-thaw cycle counts and thawed-day metrics are used to classify thermal regimes and detect warming trends. The upper-right panel features a conceptual model that connects permafrost degradation to hydrological responses: RG1, characterized as transitional, shows signs of enhanced shallow flow and seasonal meltwater pulses, while RG2 retains cold, thermally buffered conditions that support greater storage stability. These contrasts are further illustrated by temperature trend graphs, which reveal faster warming at RG1 (+ 2.84 degrees C/decade) compared to RG2 (+ 1.92 degrees C/decade), as well as increased thaw metrics. Below, a long-term streamflow graph (1970-2023) documents declining discharge, visually supported by a field photo of a dry riverbed. The bottom panel summarizes the study's key finding: RG1 and RG2 are evolving along divergent thermal and hydrological trajectories, underscoring the need for high-resolution monitoring to guide water resource planning in an era of warming and drought.