Hydrologic-land surface models (H-LSMs) offer a physically-based framework for representing and predicting the present and future states of the extensive high-latitude permafrost areas worldwide. Their primary challenge, however, is that soil temperature data are severely limited, and traditional model validation, based only on streamflow, can show the right fit to these data for the wrong reasons. Here, we address this challenge by (1) collecting existing data in various forms including in-situ borehole data and different large-scale permafrost maps in addition to streamflow data, (2) comprehensively evaluating the performance of an H-LSM with a wide range of possible process parametrizations and initializations, and (3) assessing possible trade-offs in model performance in concurrently representing hydrologic and permafrost dynamics, thereby pointing to the possible model deficiencies that require improvement. As a case study, we focus on the sub-arctic Liard River Basin in Canada, which typifies vast northern sporadic and discontinuous permafrost regions. Our findings reveal that different process parameterizations tend to align with different data sources or variables, which largely exhibit inconsistencies among themselves. We further observe that a model may fail to represent permafrost occurrence yet seemingly fit streamflows adequately. Nonetheless, we demonstrate that accurately representing essential permafrost dynamics, including the active soil layer and insulation effects from snow cover and soil organic matter, is crucial for developing high-fidelity models in these regions. Given the complexity of processes and the incompatibility among different data sources/variables, we conclude that employing an ensemble of carefully designed model parameterizations is essential to provide a reliable picture of the current conditions and future spatio-temporal co-evolution of hydrology and permafrost.

The climate in Northwest China (NWC) has undergone a warming and wetting trend (WWT) since the 1980s, which has attracted considerable attention from the scientific and policy communities. However, the majority of previous studies have focused on overall effects of WWT, and very few have examined how land surface system responds to climate warming or wetting trend, respectively. For this purpose, this study uses the Community Land Model (CLM5) driven by the Chinese Meteorological Forcing Dataset (CMFD) to conduct four modeling experiments: a control experiment (CTRL) and three sensitivity experiments, in which the annual trend of air temperature (NonWarm), precipitation (NonWet), and both (NonWWT) are removed from the CMFD from 1979 to 2018. Compared to CTRL, the land hydrological variables (i.e. soil moisture, runoff and evapotranspiration) show a visible reduction in magnitude, interannual variability, as well as annual trend in NonWet, while they are enhanced in NonWarm. In both NonWarm and NonWet, the magnitude and trend of both net radiation and sensible heat fluxes increase, with a more pronounced change in NonWWT. Further analysis indicates that the land surface processes are more sensitive to wetting trend than to warming trend. Among all land surface hydrological variables and energy variables, runoff and snow cover fraction are the most susceptible to climate change. Overall, the effects of climate change in Ta and Pr on surface hydrological variables are non-linearly offsetting, while the effects on surface energy budgets are non-linearly superimposed. Compared to warming trend, wetting trend plays a larger impact on the variability of land surface processes in NWC.

Tibetan Plateau (TP) is known as the Third Pole of the Earth. Any changes in land surface processes on the TP can have an unneglectable impact on regional and global climate. With the warming and wetting climate, the land surface of the TP saw a darkening trend featured by decreasing surface albedo over the past decades, primarily due to the melting of glaciers, snow, and greening vegetation. Recent studies have investigated the effects of the TP land surface darkening on the field of climate, but these assessments only address one aspect of the feedback loop. How do these darkening-induced climate changes affect the frozen ground and ecosystems on the TP? In this study, we investigated the impact of TP land surface darkening on regional frozen ground and ecosystems using the state-of-the-art land surface model ORCHIDEE-MICT. Our model results show that darkening-induced climate changes on the TP will lead to a reduction in the area of regional frozen ground by 1.1x104 +/- 0.019x104 km2, a deepening of the regional permafrost active layer by 0.06 +/- 0.0004 m, and a decrease in the maximum freezing depth of regional seasonal frozen ground by 0.06 +/- 0.0016 m compared to the scenario without TP land surface darkening. Furthermore, the darkening-induced climate change on the TP will result in an increase in the regional leaf area index and an enhancement in the regional gross primary productivity, ultimately leading to an increase in regional terrestrial carbon stock by 0.81 +/- 0.001 PgC. This study addresses the remaining piece of the puzzle in the feedback loop of TP land surface darkening, and improves our understanding of interactions across multiple spheres on the TP. The exacerbated regional permafrost degradation and increasing regional terrestrial carbon stock induced by TP land surface darkening should be considered in the development of national ecological security barrier.

Arctic land is characterized by a high surface and subsurface heterogeneity on different scales. However, the effects of land surface model resolution on fluxes and soil state variables in the Arctic have never been systematically studied, even though smaller scale heterogeneities are resolved in high-resolution land boundary condition datasets. Here, we compare 210 km and 5 km setups of the land surface model JSBACH3 for an idealized case study in eastern Siberia to investigate the effects of high versus low-resolution land boundary conditions on simulating the interactions of soil physics, hydrology and vegetation. We show for the first time that there are differences in the spatial averages of the simulated fluxes and soil state variables between resolution setups. Most differences are small in the summer mean, but larger within individual months. Heterogeneous soil properties induce large parts of the differences while vegetation characteristics play a minor role. Active layer depth shows a statistically significant increase of +20% in the 5 km setup relative to the 210 km setup for the summer mean and +43% for August. The differences are due to the nonlinear vertical discretization of the soil column amplifying the impact of the heterogeneous distributions of soil organic matter content and supercooled water. Resolution-induced differences in evaporation fluxes amount to +43% in July and are statistically significant. Our results show that spatial resolution significantly affects model outcomes due to nonlinear processes in heterogenous land surfaces. This suggests that resolution needs to be accounted in simulations of land surface models in the Arctic.

Permafrost in High Mountain Asia (HMA) is becoming increasingly vulnerable to thaw due to climate change. However, the lack of either in situ ground surface or borehole temperature data beyond the Tibetan Plateau prevents comprehensive assessments of its impact on the regional hydrologic cycle and local cascading hazards. Although past studies have generated estimates of permafrost extent in Central Asia, many are limited to the Tibetan Plateau, excluding the more remote reaches of the Tien Shan, Pamirs, and Himalayas. By leveraging surface temperatures from both the Moderate Resolution Imaging Spectroradiometer (MODIS) and Atmospheric Infra-Red Sounder (AIRS), this study advances further understanding of remotely sensed permafrost occurrence at high altitudes, which are prone to error due to frequent cloud cover. We demonstrate that the fusion of MODIS and AIRS products can accurately estimate long-term thermal regimes of the subsurface, with reported correlation coefficients of 0.773 and 0.560, RMSEs of 0.890 degree celsius and 0.680 degree celsius, and biases of 0.003 degree celsius and 0.462 degree celsius, respectively, for the ground surface and the depth of zero annual amplitude, during a reference period of 2003-2016. Furthermore, we provide a range of possible permafrost extents based on established equations for calculating the temperature at the top of the permafrost to demonstrate temperature sensitivity to soil moisture and snow cover. The MODIS-AIRS product is recommended to be a robust source of ground temperature estimates, which may be sufficient for inferring mountain permafrost presence in HMA. Incorporating the influence of soil moisture and snow depth, although limited by biased estimates, also produces estimates of permafrost regional areas comparable to previously reported permafrost indices. A total permafrost area of 1.69 (+/- 0.32) million km(2) is estimated for the entire HMA, across 15 mountain subregions.

The Tibetan Plateau (TP) is distributed with large areas of permafrost, which have received increasing attention as the climate warms. Accurately modeling the extent of permafrost and permafrost changes is now an important challenge for climate change research and climate modeling in this region. Uncertainty in land use and land cover (LULC), which is important information characterizing surface conditions, directly affects the accuracy of the simulation of permafrost changes in land surface models. In order to investigate the effect of LULC uncertainty on permafrost simulation, we conducted simulation experiments on the TP using the Community Land Model, version 5 (CLM5) with five high-resolution LULC products in this study. Firstly, we evaluated the simulation results using shallow soil temperature data and deep borehole data at several sites. The results show that the model performs well in simulating shallow soil temperatures and deep soil temperature profiles. The effect of different land use products on the shallow soil temperature and deep soil temperature contours is not obvious due to the small differences in land use products at these sites. Although there is little difference in the simulating results of different land use products when compared to the permafrost distribution map, the differences are noticeable for the simulation of the active layer. Land cover had a greater impact on soil temperature simulations in regions with greater land use inconsistency, such as at the junction of bare soil and grassland in the northwestern part of the TP, as well as in the southeast region with complex topography. The main way in which this effect occurs is that land cover affects the net surface radiation, which in turn causes differences in soil temperature simulations. In addition, we discuss other factors affecting permafrost simulation results and point out that increasing the model plant function types as well as carefully selecting LULC products is one of the most important ways to improve the simulation performance of land-surface models in permafrost regions.

Polar amplification appears in response to greenhouse gas forcing, which has become a focus of climate change research. However, polar amplification has not been systematically investigated over the Earth's three poles (the Arctic, Antarctica, and the Third Pole). An index of polar amplification is employed, and the annual and seasonal variations of land surface temperature over the Earth's three poles are examined using MODIS (Moderate Resolution Imaging Spectroradiometer) observations for the period 2001-2018. As expected, the warming of the Arctic is most conspicuous, followed by the Third Pole, and is weakest in Antarctica. Compared to the temperature changes for the global land region, positive polar amplification appears in the Arctic and the Third Pole on an annual scale, whereas Antarctic amplification disappears, with a negative amplification index of -0.72. The polar amplification for the Earth's three poles shows seasonal differences. Strong Arctic amplification appears in boreal spring and winter, with a surface warming rate of more than 3.40 times the global mean for land regions. In contrast, the amplification of the Third Pole is most conspicuous in boreal summer. The two poles located in the Northern Hemisphere have the weakest amplification in boreal autumn. Differently from the positive amplification for the Arctic and the Third Pole in all seasons, the faster variations in Antarctic temperature compared to the globe only appear in austral autumn and winter, and the amplification signal is negative in these seasons, with an amplification index of -1.68 and -2.73, respectively. In the austral winter, the strong negative amplification concentrates on West Antarctica and the coast of East Antarctica, with an absolute value of amplification index higher than 5 in general. Generally, the polar amplification is strongest in the Arctic except from June to August, and Antarctic amplification is the weakest among the Earth's three poles. The Earth's three poles are experiencing drastic changes, and the potential influence of climate change should receive attention.

The high-resolution permafrost distribution maps have a closer relationship with engineering applications in cold regions because they are more relative to the real situation compared with the traditional permafrost zoning mapping. A particle swarm optimization algorithm was used to obtain the index eta with 30 m resolution and to characterize the distribution probability of permafrost at the field scale. The index consists of five environmental variables: slope position, slope, deviation from mean elevation, topographic diversity, and soil bulk density. The downscaling process of the surface frost number from a resolution of 1000 m to 30 m is achieved by using the spatial weight decomposition method and index eta. We established the regression statistical relationship between the surface frost number after downscaling and the temperature at the freezing layer that is below the permafrost active layer base. We simulated permafrost temperature distribution maps with 30 m resolution in the four periods of 2003-2007, 2008-2012, 2013-2017, and 2018-2021, and the permafrost area is, respectively, 28.35 x 10(4) km(2), 35.14 x 10(4) km(2), 28.96 x 10(4) km(2), and 25.21 x 10(4) km(2). The proportion of extremely stable permafrost (< -5.0 degrees C), stable permafrost (-3.0 similar to -5.0 degrees C), sub-stable permafrost (-1.5 similar to -3.0 degrees C), transitional permafrost (-0.5 similar to -1.5 degrees C), and unstable permafrost (0 similar to -0.5 degrees C) is 0.50-1.27%, 6.77-12.45%, 29.08-33.94%, 34.52-39.50%, and 19.87-26.79%, respectively, with sub-stable, transitional, and unstable permafrost mainly distributed. Direct and indirect verification shows that the permafrost temperature distribution maps after downscaling still have high reliability, with 83.2% of the residual controlled within the range of +/- 1 degrees C and the consistency ranges from 83.17% to 96.47%, with the identification of permafrost sections in the highway engineering geological investigation reports of six highway projects. The maps are of fundamental importance for engineering planning and design, ecosystem management, and evaluation of the permafrost change in the future in Northeast China.

Downward solar radiation (DSR) and air temperature (Ta) have significant influences on the thermal state of frozen ground. These parameters are also important forcing terms for physically based land surface models (LSMs). However, the quantitative influences of inaccuracies in DSR and Ta products on simulated frozen ground temperatures remain unclear. In this study, three DSR products (CMFD-SR, Tang-SR, and GLDAS-SR) and two Ta products (CMFD-Ta and GLDAS-Ta) were used to force an LSM model in an alpine watershed in Northwest China, to investigate the sensitivity of simulated ground temperatures to different DSR and Ta products. Compared to a control model (CTRL) forced by in situ observed DSR, ground temperatures simulated by the experimental model forced by GLDAS-SR are obviously decreased because GLDAS-SR is much lower than in situ observations. Instead, simulation results in models forced by CMFD-SR and Tang-SR are much closer to those of CTRL. Ta products led to significant errors in simulated ground temperatures. In conclusion, both CMFD-SR and Tang-SR could be used as good alternatives to in situ observed DSR for forcing a model, with acceptable errors in simulation results. However, more care need to be paid for models forced by Ta products instead of Ta observations, and conclusions should be carefully drawn.

Freezing/thawing indices are important indicators of the dynamics of frozen ground on the Qinghai-Tibet Plateau (QTP), especially in areas with limited observations. Based on the numerical outputs of Community Land Surface Model version 4.5 (CLM4.5) from 1961 to 2010, this study compared the spatial and temporal variations between air freezing/thawing indices (2 m above the ground) and ground surface freezing/thawing indices in permafrost and seasonally frozen ground (SFG) across the QTP after presenting changes in frozen ground distribution in each decade in the context of warming and wetting. The results indicate that an area of 0.60 x 10(6) km(2) of permafrost in the QTP degraded to SFG in the 1960s-2000s, and the primary shrinkage period occurred in the 2000s. The air freezing index (AFI) and ground freezing index (GFI) decreased dramatically at rates of 71.00 & DEG;C & BULL;d/decade and 34.33 & DEG;C & BULL;d/decade from 1961 to 2010, respectively. In contrast, the air thawing index (ATI) and ground thawing index (GTI) increased strikingly, with values of 48.13 & DEG;C & BULL;d/decade and 40.37 & DEG;C & BULL;d/decade in the past five decades, respectively. Permafrost showed more pronounced changes in freezing/thawing indices since the 1990s compared to SFG. The changes in thermal regimes in frozen ground showed close relations to air warming until the late 1990s, especially in 1998, when the QTP underwent the most progressive warming. However, a sharp increase in the annual precipitation from 1998 began to play a more controlling role in thermal degradation in frozen ground than the air warming in the 2000s. Meanwhile, the following vegetation expansion hiatus further promotes the thermal instability of frozen ground in this highly wet period.