An anomalous warm weather event in the Antarctic McMurdo Dry Valleys on 18 March 2022 created an opportunity to characterize soil biota communities most sensitive to freeze-thaw stress. This event caused unseasonal melt within Taylor Valley, activating stream water and microbial mats around Canada Stream. Liquid water availability in this polar desert is a driver of soil biota distribution and activity. Because climate change impacts hydrological regimes, we aimed to determine the effect on soil communities. We sampled soils identified from this event that experienced thaw, nearby hyper-arid areas, and wetted areas that did not experience thaw to compare soil bacterial and invertebrate communities. Areas that exhibited evidence of freeze-thaw supported the highest live and dead nematode counts and were composed of soil taxa from hyper-arid landscapes and wetted areas. They received water inputs from snowpacks, hyporheic water, or glacial melt, contributing to community differences associated with organic matter and salinity gradients. Inundated soils had higher organic matter and lower conductivity (p < .02) and hosted the most diverse microbial and invertebrate communities on average. Our findings suggest that as liquid water becomes more available under predicted climate change, soil communities adapted to the hyper-arid landscape will shift toward diverse, wetted soil communities.

Surface soil moisture (SSM) is a key limiting factor for vegetation growth in alpine meadow on the Qinghai-Tibetan Plateau (QTP). Patches with various sizes and types may cause the redistribution of SSM by changing soil hydrological processes, and then trigger or accelerate alpine grassland degradation. Therefore, it is vital to understand the effects of patchiness on SSM at multi-scales to provide a reference for alpine grassland restoration. However, there is a lack of direct observational evidence concerning the role of the size and type of patches on SSM, and little is known about the effects of patches pattern on SSM at plot scale. Here, we first measured SSM of typical patches with different sizes and types at patch scale and investigated their patterns and SSM spatial distribution through unmanned aerial vehicle (UAV)-mounted multi-type cameras at plot scale. We then analyzed the role of the size and type of patchiness on SSM at both patch and plot scales. Results showed that: (1) in situ measured SSM of typical patches was significantly different (P < 0.01), original vegetation patch (OV) had the highest SSM, followed by isolate vegetation patch (IV), small bare patch (SP), medium bare patch (MP) and large bare patch (LP); (2) the proposed method based on UAV images was able to estimate SSM (0-40 cm) with a satisfactory accuracy (R-2 = 0.89, P < 0.001); (3) all landscape indices of OV, with the exception of patch density, were positively correlated with SSM at plot scale, while most of the landscape indices of LP and IV showed negative correlations (P < 0.05). Our results indicated that patchiness intensified the spatial heterogeneity of SSM and potentially accelerated the alpine meadow degradation. Preventing the development of OV into IV and the expansion of LP is a critical task for alpine meadow management and restoration.

Soil freeze-thaw state influences multiple terrestrial ecosystem processes, such as soil hydrology and carbon cycling. However, knowledge of historical long-term changes in the timing, duration, and temperature of freeze-thaw processes remains insufficient, and studies exploring the combined or individual contributions of climatic factors-such as air temperature, precipitation, snow depth, and wind speed-are rare, particularly in current thermokarst landscapes induced by abrupt permafrost thawing. Based on ERA5-Land reanalysis, MODIS observations, and integrated thermokarst landform maps, we found that: 1) Hourly soil temperature from the reanalysis effectively captured the temporal variations of in-situ observations, with Pearson' r of 0.66-0.91. 2) Despite an insignificant decrease in daily freeze-thaw cycles in 1981-2022, other indicators in the Qinghai-Tibet Plateau (QTP) changed significantly, including delayed freezing onset (0.113 d yr- 1), advanced thawing onset (-0.22 d yr- 1), reduced frozen days (-0.365 d yr- 1), increased frozen temperature (0.014 degrees C yr- 1), and decreased daily freeze-thaw temperature range (-0.015 degrees C yr- 1). 3) Total contributions indicated air temperature was the dominant climatic driver of these changes, while indicators characterizing daily freeze-thaw cycles were influenced mainly by the combined effects of increased precipitation and air temperature, with remarkable spatial heterogeneity. 4) When regionally averaged, completely thawed days increased faster in the thermokarstaffected areas than in their primarily distributed grasslands-alpine steppe (47.69%) and alpine meadow (22.64%)-likely because of their stronger warming effect of precipitation. Locally, paired comparison within 3 x 3 pixel windows from MODIS data revealed consistent results, which were pronounced when the thermokarst-affected area exceeded about 38% per 1 km2. Conclusively, the warming and wetting climate has significantly altered soil freeze-thaw processes on the QTP, with the frozen soil environment in thermokarstaffected areas, dominated by thermokarst lakes, undergoing more rapid degradation. These insights are crucial for predicting freeze-thaw dynamics and assessing their ecological impacts on alpine grasslands.

Vast deserts and sandy lands in the mid-latitudes cover an area of 17.64 x 106 km2, with 6.98 x 106 km2 experiencing seasonal frozen soil (SFG). Freeze-thaw cycles of SFG significantly influence local surface processes in deserts, impacting meteorological disasters such as infrastructure failures and sandstorms. This study investigates the freeze-thaw dynamics of SFG in crescent dunes from three deserts in northern China: the Tengger Desert, Mu Us Sandy Land, and Ulan Buh Desert, over the period from 2019 to 2024.Freezing occurs from November to January, followed by thawing from January to March. The thawing rate (2.72 cm/day) was 1.8 times higher than the freezing rate (1.48 cm/day). The maximum seasonal freezing depth (MSFD) exceeded 0.80 mat all dune slopes, with depths surpassing 1.10 mat the leeward slope and lower slope positions. Soil moisture content, ranging from 1 % to 1.6 %, is critical for freezing, and this threshold varies depending on the dune's mechanical composition. The hardness of frozen desert soil is primarily controlled by moisture, along with temperature and particle size.Temperature initiates freezing, while moisture and particle size control the resulting hardness.These findings shed light on the seasonal freeze-thaw processes in desert soils and have practical implications for agricultural management, engineering design, and environmental hazard mitigation in arid regions.

Alpine wet meadow (AWM), an important wetland type on the Qinghai-Tibet Plateau (QTP), is sensitive to climate change, which alters the soil hydrothermal regime and impacts ecological and hydrological functions in permafrost regions. The mechanisms underlying extreme AWM degradation in the QTP and hydrothermal factors controlling permafrost degradation remain unclear. In this study, soil hydrothermal processes, soil heat migration, and the permafrost state were measured in AWM and extremely degraded AWM (EDAWM). The results showed that the EDAWM exhibited delayed onset of both soil thawing and freezing, shortened thawing period, and extended freezing period at the lower boundary of the active layer. The lower ground temperatures resulted in a 0.2 m shallower active layer thickness in the EDAWM compared with the AWM. Moreover, the EDAWM altered soil thermal dynamics by redistributing energy, modifying soil moisture, preserving soil organic matter, and adjusting soil thermal properties. As for energy budget, a substantial amount of heat in the EDAWM was consumed by turbulent heat fluxes, particularly latent heat flux, which reduced the amount of heat transferred to the ground. Additionally, the higher soil organic matter content in EDAWM decreased the annual mean soil thermal conductivity from 1.42 W m- 1 K-1 in AWM to 1.26 W m- 1 K-1 in EDAWM, slowing down heat transfer within the active layer and consequently mitigating permafrost degradation. However, with continued climate warming, the soil organic matter content in EDAWM will inevitably decline due to microbial decomposition in the absence of new organic inputs. As the soil organic matter content diminishes, soil heat transfer processes will likely accelerate, and the permafrost warming rate may surpass that in undistributed AWM. These findings enhance our understanding of how alpine ecosystem succession influences regional hydrological cycles and greenhouse gas emissions.

Permafrost is undergoing rapid changes due to climate warming, potentially exposing a vast reservoir of carbon to be released to the atmosphere, causing a positive feedback cycle. Despite the importance of this feedback, its specifics remain poorly constrained, because representing permafrost dynamics still poses a significant challenge for Earth System Models (ESMs). This review assesses the current state of permafrost representation in land surface models (LSMs) used in ESMs and offline permafrost models, highlighting both the progress made and the remaining gaps.We identify several key physical processes crucial for permafrost dynamics, including soil thermal regimes, freeze-thaw cycles, and soil hydrology, which are underrepresented in many models. While some LSMs have advanced significantly in incorporating these processes, others lack fundamental elements such as latent heat of freeze-thaw, deep soil columns, and Arctic vegetation dynamics. Offline permafrost models provide valuable insights, offering detailed process testing and aiding the prioritization of improvements in coupled LSMs.Our analysis reveals that while significant progress has been made in incorporating permafrost-related processes into coupled LSMs, many small-scale processes crucial for permafrost dynamics remain underrepresented. This is particularly important for capturing the complex interactions between physical and biogeochemical processes required to model permafrost carbon dynamics. We recommend leveraging advancements from offline permafrost models and progressively integrating them into LSMs, while recognizing the computational and technical challenges that may arise in coupled simulations. We highlight the importance of enhancing the representation of physical processes, including through improvements in model resolution and complexity, as this is a fundamental precursor to accurately incorporate biogeochemical processes and capture the permafrost carbon feedback.

Soil organic carbon (SOC) in the active layer (0-2 m) of the Tibetan Plateau (TP) permafrost region is sensitive to climate change, with significant implications for the global carbon cycle. Environmental factors-including parent material, climate, vegetation, topography, soil, and human activities-inevitably drive SOC variations. However, vegetation and climate are likely the two most influential factors impacting SOC variations. To test this hypothesis, we conducted experiments using 31 environmental variables combined with the recursive feature elimination (RFE) algorithm. These experiments showed that RFE retained all vegetation variables [Land cover types (LCT), normalized difference vegetation index (NDVI), leaf area index (LAI), and gross primary productivity (GPP)] as well as two climate variables [Moisture index (MI) and drought index (DI)], supporting our hypothesis. We then analyzed the relationship between SOC and the retained vegetation and climate variables using random forest (RF), Shapley additive explanations (SHAP), and GeoDetector models to quantify the independent and interactive drivers of SOC distribution and to identify the optimal conditions for SOC accumulation. The RF model explained 68% and 42% of the spatial variability in SOC at depths of 0-1 m and 1-2 m, respectively, with SOC stocks higher in the southeast and lower in the northwest. Additionally, SOC stock at 0-1 m was significantly higher (p 0.05). Spearman correlation coefficients results indicated that NDVI, LAI, GPP, and MI had highly significant positive correlations with SOC (p < 0.01), whereas DI had a highly significant negative correlation with SOC (p < 0.01). SHAP analysis revealed environmental thresholds for SOC variations, with notable shifts at NDVI (0.40), LAI (7), GPP (250 g C m(-)(2) year(-)(1)), MI (0.40), and DI (0.50). The spatial distribution of these thresholds aligns with the 400 mm equivalent precipitation line. Additionally, GeoDetector results emphasized that interactions between climate and vegetation factors enhance the explanatory power of individual variables on SOC variations. The swamp meadow type, with an NDVI range of 0.73-0.84, LAI range of 11.06-15.94, and MI range of 0.46-0.56, was identified as the most favorable environment for SOC accumulation. These findings are essential for balancing vegetation and climate conditions to sustain SOC levels and mitigate climate change-driven carbon release.

The fine-scale controls of active layer dynamics remain poorly understood, particularly at the southern boundary of continuous permafrost. We examined how environmental conditions associated with upland tundra heath, open graminoid fen, and palsa/peat plateau landforms affected active layer thermal regime (timing, magnitude, and rate of thaw) in a subarctic peatland in the Hudson Bay Lowlands, Canada. A significant increase in active layer thaw depth was evident between 2012 and 2024. Within-season thaw patterns differed among landforms, with tundra heath exhibiting the highest thaw rates and soil temperatures, succeeded by fen and palsa. Air temperature mediated by soil properties, topography, and vegetation affected thaw patterns. The increased thermal conductivity of gravel/sandy tundra heath soils exerted a more pronounced influence on thaw patterns relative to fens and palsas, both of which had a thicker organic layer. Near-surface soil moisture was the lowest in tundra, followed by palsas, and fens. Increased soil moisture impeded active layer thaw, likely due to a combination of soil surface evaporation and meltwater percolation. These findings elucidate the relationship between the biophysical properties of landform features and climate, revealing their role in influencing active layer thaw patterns in a subarctic ecosystem.

Soil thermal conductivity (STC) plays a crucial role in regulating the energy distribution of both the surface and underground soil layers. It is widely applied in various fields, including engineering design, geothermal resource development and climate change research. A rapid and accurate estimation of STC remains a key focus in the study of soil thermodynamic parameters. However, the methods for estimating STC and their distinct characteristics have yet to be systematically reviewed. In this study, we used bibliometrics to comprehensively and systematically review the literature on STC, focusing on knowledge graph characteristics to analyze the development trend of calculation schemes. The main conclusions drawn from the study are as follows: (1) In recent years, most studies have been focused on soil thermal characteristics and their main contributing factors, the soil hydrothermal process in the Qinghai-Tibet Plateau, geothermal equipment and numerical simulations, and the exploration of geothermal resources. (2) A systematic review of various schemes indicates that no single scheme is universally applicable to all soil types. Moreover, a single parameterization scheme fails to meet the practical requirements of land surface process models. We evaluated the advantages and disadvantages of the traditional heat conduction schemes, parameterization schemes, and machine learning-based schemes and the findings suggest that a comprehensive scheme that integrates these three different schemes for STC simulations should be urgently developed.

Reanalysis is a valuable potential data source for permafrost studies. The latest-generation reanalysis of the Japanese Reanalysis for three quarters of a century (JRA-3Q) benefits from improved snow and soil schemes and demonstrates encouraging performance for soil temperature in permafrost regions compared to its predecessor, JRA-55, and other state-of-the-art reanalyses. We find JRA-3Q to have an overall mean annual air temperature bias of-0.17 degrees C, with-0.55 degrees C in permafrost regions. The snow depth was underestimated by-5.5 cm. In permafrost regions, the mean annual ground temperature bias was about-0.09 degrees C. The estimated permafrost area from JRA-3Q is between 10.8 and 15.8 x 106 km2. The active layer thickness is substantially overestimated by about 0.65 m. The JRA-3Q soil temperature exhibits a pronounced warm bias in Alaska, which is very likely due to the overestimated snow insulation and simplified soil organic content. The decoupled energy conservation parameterization (DECP) method employed in the JRA-3Q soil scheme restricts its suitability for the interpretation of detailed permafrost phenomena, such as zero-curtain effects. This DECP method is used in many stateof-the-art land surface models; our results demonstrate the need for additional contributions to improve the representation of permafrost-specific processes.