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.

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.

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.

The distribution and variation of active layer thickness (ALT) are crucial indicators for assessing the stability and environmental conditions of permafrost regions, which significantly impact regional hydrology, ecology, climate change, engineering construction, and disaster risk assessment. Based on the measured ALT data and Stefan equation, this study investigated the spatial distribution characteristics of ALT in the Tuotuo River region and explored the factors influencing its variability. The results showed that the ALT in the Tuotuo River area ranged from 0.15 to 5.18 m, with an average value of 2.65 m. The spatial distribution showed that the ALT was thinner in the southern region, which exhibited strong spatial heterogeneity, while the northeastern region generally had thicker ALT. Additionally, mountain areas tended to have thinner ALT, whereas plains showed thicker ALT. There was a good linear correlation between the simulated and measured ALT values, and the R 2 was up to 80%. The ALT in the Tuotuo River area was mainly controlled by air temperature and surface water thermal conditions. Among all factors, soil water content was identified as the key determinant. Topographic factors influenced ALT distribution and variation mainly through their impact on soil water content.

Increasing Arctic warming rates drive significant environmental change, including permafrost thaw and new groundwater pathway development, thereby increasing groundwater vulnerability to contaminant transport at the thousands of unremediated sites in the circumpolar north. As a first step in assessing hydrogeological controls of Arctic contaminant transport, this study uses numerical modelling to disentangle the impacts of increasing precipitation and air temperature on groundwater flow within the active layer at a high arctic site (63 degrees 30 ' N). The study uses the numerical model SUTRA 4.0 to simulate groundwater flow and energy transport, including dynamic freeze-thaw processes, across an Arctic hillslope under current and future air temperatures due to climate warming. The model domain represents a two-dimensional hillslope terminating in a lake. Two layers implemented in the model represent unconsolidated glacial till and underlying crystalline bedrock. Four simulation cases are examined based on downscaled CMIP5 projections under the high-emissions business as usual scenario: Baseline Conditions (1981-2010), Near-Projections (2011-2040), Mid-Projections (2041-2070), and Far Projections (2071-2100). Climate projections indicate increasing mean annual air temperatures, reducing annual air temperature amplitude, and increasing precipitation. Further, model results show that groundwater flow dynamics are primarily influenced by the coupling of both increased mean annual temperatures and precipitation, with the consequent deepening and prolonged thawing of the active layer allowing for increased groundwater exfiltration to the lake. Sensitivity analysis identifies overburden permeability, overburden residual liquid freezing temperature, and model base temperature as significant parameters that affect model outcomes. Finally, a variable transmissivity assessment provides new insight into active layer groundwater flows.

Ongoing climate change is critically endangering cold regions, with the Arctic warming at nearly four times the global average. This rapid warming is not only accelerating the irreversible thawing of permafrost but is also reshaping the region's topography, vegetation, hydrology, infrastructure integrity, and carbon exchange. The destabilization of the ground through thaw of ice-rich permafrost, known as thermokarst, is increasing to mass-wasting events such as active-layer detachment failures (ALDs), shallow landslides that are becoming increasingly common in the Arctic. In light of these alarming developments, our study employs the Maxent statistical model to analyze ALD distribution, develop a susceptibility map for Alaska and Northwest Territories, Canada, in the current climate, and assess the potential impact to infrastructure. We identified high-susceptibility zones across critical regions, including the Brooks Range, Franklin Mountains, and West Crazy Mountains in Alaska, as well as the Dawson City and Mackenzie River areas in Canada. Particularly concerning is the vulnerability of linear infrastructure: 878 km of roads, 167 km of the Trans-Alaska pipeline, and 140 km of the Norman Wells pipeline are situated in areas of high to very high susceptibility to ALDs. These results highlight the urgent need for proactive strategies and infrastructure planning to deal with the growing threats from permafrost thaw and its wide-ranging effects.

This manuscript presents a comprehensive presentation of ground temperature data collected at 16 nodes of the 121 of the Crater Lake Circumpolar Active Layer Monitoring (CALM) site on Deception Island, Antarctica, from 2008 to early 2022. Each one of the 16 shallow boreholes has been equipped with miniature temperature loggers, providing valuable insights into the thermal regime of the ground at a depth of 50 cm, which corresponds to the mean depth of the top of the permafrost table as observed by annual mechanical probing in the CALM site. Despite a 9-month long gap in data collection during 2017 due to persistent snow cover, the time series remains largely intact, with annual measurements taken every 3 h. The manuscript details the methodologies employed for data collection, including the use of iButton loggers, and outlines the challenges faced in retrieving and processing the data in the harsh Antarctic environment. The cleaned dataset, which consolidates data from various nodes while removing erroneous records, is made freely accessible to the scientific community without any additional processing of the data such as offset corrections or gaps interpolation. This resource is expected to facilitate further research into the thermal dynamics of the active layer and permafrost and its implications for climate change since both are influenced by external factors such as snow cover, air temperature and others. Overall, the presented dataset contributes to the limited body of knowledge regarding Antarctic permafrost and provides a foundation for future investigations into the effects of climate change on frozen ground dynamics. The dataset serves as a vital tool for researchers aiming to model ground thermal behaviour and assess the impacts of environmental changes in polar regions.

Southern boundary areas of high-latitude permafrost regions may represent the future permafrost temperature regimes; therefore, understanding the carbon stocks and their stability in these systems can shed light on the permafrost carbon cycle under a warming climate. In this study, we sampled soils at three sites representing three differing land covers (forest swamp, dry forest, and shrub swamp) located in the southern boundary area of a high-latitude permafrost region and investigated their carbon fractions and the relationships of these fractions with soil physicochemical parameters in the active and permafrost layers. The results show that the proportion of active carbon is higher in permafrost than in the active layer under forest swamp and dry forest, implying that carbon pools in the permafrost are more decomposable. However, in shrub swamp, the active carbon components in the permafrost layer are lower than in the active layer. Soil pH and water content are the most significant factors associated with soil organic carbon concentration both in the active layers and in the permafrost layers. Our results suggest that, although soil organic carbon concentrations largely decrease with depth, the proportion of the forest swamp, dry forest labile carbon is higher in the permafrost layer than in the active layer and that the vertical distribution of labile carbon proportions is related to land covers.

In the context of global warming, understanding the impact of thaw slump on soil hydrothermal processes and its responses to climate is essential for protecting engineering facilities in cold regions. This study aimed to investigate the effect of thaw slump development on active layer soil. We considered the early thaw slump development in the Tibetan Plateau as research object and conducted long-term monitoring of soil hydrothermal activity in the active layer of various parts of the landslide and the regional meteorology. The results showed that thaw slump development shortened the freezing and thawing time of the active layer, increased the freezing and thawing rates of the shallow soil (10-20 cm), and enhanced the heat exchange between the active layer soil and the atmosphere and the heat transfer between the soils. The heat-exchange efficiency of the active layer, from largest to smallest, was headwall > collapsed area > unaffected area (bottom of the slope) > unaffected area (top of the slope). Furthermore, thaw slump development lowered the water storage of the active layer prof ile and weakened the dynamic response of soil water to precipitation. The events of soil water responses and soil water increments were smaller in the landslide area than in the unaffected area. During a co-precipitation event, the overall soil water storage increment (SWSI) of the profile was significantly smaller in the landslide area than in the unaffected area (P < 0.05), with an SWSI of 2.04 mm in the headwall and 1.77 mm in the collapsed area. In addition, thaw slump development altered the mechanism of soil water transport driven by soil temperature changes, which affected soil water redistribution of profile. The study gives ecohydrology-related research in cold climates a scientific foundation, thereby guiding the construction and maintenance of infrastructure projects.

Permafrost is one of the crucial components of the cryosphere, covering about 25% of the global continental area. The active layer thickness (ALT), as the main site for heat and water exchange between permafrost and the external atmosphere, its changes significantly impact the carbon cycle, hydrological processes, ecosystems, and the safety of engineering structures in cold regions. This study constructs a Stefan CatBoost-ET (SCE) model through machine learning and Blending integration, leveraging multi-source remote sensing data, the Stefan equation, and measured ALT data to focus on the ALT in the Qinghai-Tibet Plateau (QTP). Additionally, the SCE model was verified via ten-fold cross-validation (MAE: 20.713 cm, RMSE: 32.680 cm, R2: 0.873, and MAPE: 0.104), and its inversion of QTP's ALT data from 1958 to 2022 revealed 1998 as a key turning point with a slow growth rate of 0.25 cm/a before 1998 and a significantly increased rate of 1.26 cm/a afterward. Finally, based on multiple model input factor analysis methods (SHAP, Pearson correlation, and Random Forest Importance), the study analyzed the ranking of key factors influencing ALT changes. Meanwhile, the importance of Stefan equation results in SCE model is verified. The research results of this paper have positive implications for eco-hydrology in the QTP region, and also provide valuable references for simulating the ALT of permafrost.