Ice records provide a qualitative rather than a quantitative indication of the trend of climate change. Using the bulk aerodynamic method and degree day model, this study quantified ice mass loss attributable to sublimation/evaporation (S/E) and meltwater on the basis of integrated observations (1960-2006) of glacier-related and atmospheric variables in the northeastern Tibetan Plateau. During 1961-2005, the average annual mass loss in the ice core was 95.33 +/- 20.56 mm w.e. (minimum: 78.97 mm w.e. in 1967, maximum: 146.67 mm w.e. in 2001), while the average ratio of the revised annual ice accumulation was 21.2 +/- 7.7% (minimum: 11.0% in 1992, maximum 44.8% in 2000). A quantitative formula expressing the relationship between S/E and air temperature at the monthly scale was established, which could be extended to estimation of S/E changes of other glaciers in other regions. The elevation effect on alpine precipitation determined using revised ice accumulation and instrumental data was found remarkable. This work established a method for quantitative assessment of the temporal variation in ice core mass loss, and advanced the reconstruction of long-term precipitation at high elevations. Importantly, the formula established for reconstruction of S/E from temperature time series data could be used in other regions.

Lakes are commonly accepted as a sensitive indicator of regional climate change, including the Tibetan Plateau (TP). This study took the Ranwu Lake, located in the southeastern TP, as the research object to investigate the relationship between the lake and regional hydroclimatological regimes. The well-known Budyko framework was utilized to explore the relationship and its causes. The results showed air temperature, evapotranspiration and potential evapotranspiration in the Ranwu Lake Basin generally increased, while precipitation, soil moisture, and glacier area decreased. The Budyko space indicated that the basin experienced an obviously drying phase first, and then a slightly wetting phase. An overall increase in lake area appears inconsistent with the drying phase of the basin climate. The inconsistency is attributable to the significant expansion of proglacial lakes due to glacial melting, possibly driven by the Atlantic Multidecadal Oscillation. Our findings should be helpful for understanding the complicated relationships between lakes and climate, and beneficial to water resources management under changing climates, especially in glacier basins.

Numerous endorheic lakes in the Qinghai-Tibet Plateau (QTP) have shown a dramatic increase in total area since 1996. These expanding lakes are mainly located in the interior regions of the QTP, where permafrost is widely distributed. Despite significant permafrost degradation due to global warming, the impact of permafrost thawing on lake evolution in QTP has been underexplored. This study investigated the permafrost degradation and its correlation with lake area increase by selecting four lake basins (Selin Co, Nam Co, Zhari Namco, and Dangqiong Co) in QTP for analysis. Fluid-heat-ice coupled numerical models were conducted on the aquifer cross-sections in these four lake basins, to simulate permafrost thawing driven by rising surface temperatures, and calculate the subsequent changes in groundwater discharge into the lakes. The contribution of these changes to lake storage, which is proportional to lake area, was investigated. Numerical simulation indicates that from 1982 to 2011, permafrost degradation remained consistent across the four basins. During this period, the active layer thickness first increased, then decreased, and partially transformed into talik, with depths reaching up to 25 m. By 2011, groundwater discharge had significantly risen, exceeding 2.9 times the initial discharge in 1988 across all basins. This increased discharge now constitutes up to 17.67 % of the total lake water inflow (Selin Co). The dynamic lake water budget further suggests that groundwater contributed significantly to lake area expansion, particularly since 2000. These findings highlight the importance of considering permafrost thawing as a crucial factor in understanding the dynamics of lake systems in the QTP in the context of climate change.

Snow cover variation significantly impacts alpine vegetation dynamics on the Tibetan Plateau (TP), yet this effect under climate change remains underexplored. This study uses a survival analysis model to assess the influence of snow on vegetation green-up dynamics, while controlling for key temperature and water availability factors. This analysis integrates multi-source data, including satellite-derived vegetation green-up dates (GUDs), snow depth, accumulated growing degree days (AGDD), downward shortwave radiation (SRAD), precipitation, and soil moisture. Our survival analysis model effectively simulated GUD on the TP, achieving an R of 0.62 (p < 0.01), a root mean square error (RMSE) of 11.20 days, and a bias of -1.41 days for 2020 GUD predictions. It outperformed both the model excluding snow depth and a linear regression model. By isolating snow's impact, we found it exerts a stronger influence on vegetation GUD than precipitation in snow-covered areas of the TP. Furthermore, snow depth effects varied seasonally: a 1-cm increase in preseason snow depth reduced green-up rates by 8.48% before 156(th) day but increased them by 4.74% afterward. This indicates that deeper preseason snow cover delays GUD before June, but advances it from June onward, rather than having a uniform effect. These findings highlight the critical role of snow and underscore the need to incorporate its distinct effects into vegetation phenology models in alpine regions.

Glaciers provide multiple ecosystem services (ES) to human society. Due to the continued global warming, the valuation of glacier ES is of urgent importance because this knowledge can support the protection of glaciers. However, a systematic valuation of glacier ES is still lacking, particularly from the perspective of ES contributors. In this study, we introduce the concept of emergy to establish a methodological framework for accounting glacier ES values, and take the Tibetan Plateau (TP) as a case study to comprehensively evaluate the spatiotemporal characteristics of glacier ES during the early 21st century. The results show that the total glacier ES values on the TP increased from 2.36E+24 sej/yr in the 2000s to 2.40E+24 sej/yr in the 2010s, with an overall growth rate of 1.6%. The values of the various services in the 2010s are ranked in descending order: climate regulation (1.59E+24 sej/yr, 66.1%), runoff regulation (4.40E+23 sej/yr, 18.4%), hydropower generation (1.88E+23 sej/ yr, 7.8%). Significantly higher glacier ES values were recorded in the marginal TP than in the endorheic area. With the exception of climate regulation and carbon sequestration, all other service values increased during the study period, partially cultural services, which have experienced rapid growth in tandem with social development. The results of this study will help establish the methodological basis for the assessment of regional and global glacier ES, as well as a scientific basis for the regional protection of glacier resources.

As a key component of the cryosphere, permafrost is sensitive to climate change, but mapping permafrost, especially in the Tibetan Plateau, has been challenging due to the heterogeneous mountainous landscape and limited representativeness of ground observations. Using 155 compiled ground observations and more than 20,000 rock glacier records, we developed a machine learning model to map the distribution of permafrost and produce an improved permafrost zonation index (PZI) map. The model was applied by incorporating several control variables, including terrain (elevation and relief), soil (bulk density, clay, coarse fragments, sand, and silt), and temperature (MAAT, FDD, and TDDT) to estimate the PZI at a 1-km resolution in the southern Tibetan Plateau. Excluding glaciers and lakes, the area of permafrost estimated by the new map is approximately 103.5 x 103 km2, accounting for 47.8% of the total area of the region. The result was assessed with various datasets and compared with existing permafrost maps and achieved higher accuracy compared with previous studies. The overall classification accuracy was 96.1% in high plain areas and 84.4% in mountain areas. The results demonstrated the substantial potential for improving mapping permafrost and understanding the periglacial environment with rock glacier inventories and machine learning, especially in complex terrain and climate.

Due to the impact of climate change, significant alterations in snowmelt have already occurred, which have been demonstrated to play a crucial role in photosynthetic carbon sequestration processes in vegetation. However, the effect of changes in snowmelt on light use efficiency (LUE) of grassland remain largely unknown in the permafrost region of Qinghai-Tibetan Plateau (QTP). By utilizing remote sensing data from 2000 to 2017, we conducted an analysis on the spatial and temporal patterns of LUE for various types of permafrost and grassland on the QTP. The LUE of the growing season was 1.1588 g CMJ(-1), displaying variations among different ecosystems: alpine steppe of seasonally frozen ground (ASS) > alpine meadow of seasonally frozen ground (AMS) > alpine meadow of permafrost (AMP) > alpine steppe of permafrost (ASP). Furthermore, our study demonstrated that decreasing snowmelt during the growing season had a negative impact on LUE through meteorological factors, elucidating its influence on LUE for approximately 40.65%, 34.06%, 41.05%, and 32.68% of ecosystems studied. Reduced snowmelt indirectly affects LUE by lowering air temperatures, vapor pressure deficit and solar radiation, while replenishing soil moisture. Additionally, changes in snowmelt can directly affect LUE by reducing the insulating properties of snow cover. Therefore, when estimating gross primary productivity (GPP) using remote sensing data based on LUE, it is essential to consider the impact of snowmelt. This will better represent vegetation phenology's response to climate change.

The Tibetan Plateau (TP) has experienced accelerated warming in recent decades, especially in winter. However, a comprehensive quantitative study of its long-term warming processes during daytime and nighttime is lacking. This study quantifies the different processes driving the acceleration of winter daytime and nighttime warming over the TP during 1961-2022 using surface energy budget analysis. The results show that the surface warming over the TP is mainly controlled by two processes: (a) a decrease in snow cover leading to a decrease in albedo and an increase in net downward shortwave radiation (snow-albedo feedback), and (b) a warming in tropospheric temperature (850 - 200 hPa) leading to an increase in downward longwave radiation (air warming-longwave radiation effect). The latter has a greater impact on the spatial distribution of warming than the former, and both factors jointly influence the elevation dependent warming pattern. Snow-albedo feedback is the primary factor in daytime warming over the monsoon region, contributing to about 59% of the simulated warming trend. In contrast, nighttime warming over the monsoon region and daytime/nighttime warming in the westerly region are primarily caused by the air warming-longwave radiation effect, contributing up to 67% of the simulated warming trend. The trend in the near-surface temperature mirrors that of the surface temperature, and the same process can explain changes in both. However, there are some differences: an increase in sensible heat flux is driven by a rise in the ground-atmosphere temperature difference. The increase in latent heat flux is associated with enhanced evaporation due to increased soil temperature and is also controlled by soil moisture. Both of these processes regulate the temperature difference between ground and near-surface atmosphere.

Long-term, high-resolution soil moisture (SM) is a vital variable for understanding the water-energy cycle and the impacts of climate change on the Qinghai-Tibet Plateau (QTP). However, most existing satellite SM data are only available at coarse scale (+/- 25 km) and suffer a lot from data gaps due to satellite orbit coverage and snow cover, especially on the QTP. Although substantial efforts have been devoted to downscale SM utilizing multiple soil moisture indices (SMIs) or diverse machine learning (ML) methods, the potentials of different SMIs and ML approaches in SM downscaling on the complex plateau remain unclear, and there is still a necessity to obtain an accurate, long-term, high-resolution and seamless SM data over the QTP. To address this issue, this study generated the long-term, high-accuracy and seamless soil moisture dataset (LHS-SM) over the QTP during 2001-2020 using a two-step downscaling method (first downscaling then merging). Firstly, the daily SM data from the Climate Change Initiative program of the European Space Agency (ESA CCI) was downscaled to 1 km utilizing five ML approaches. Then, a dynamic data merging method that considers spatiotemporal nonstationary error was applied to derive the final LHS-SM data. The performance of fifteen SMIs was also assessed and the optimal indexes for downscaling were identified. Results indicated that the shortwave infrared band-based indices had better performance than the near infrared band-based and energy-based indices. The generated LHS-SM data exhibited satisfying accuracy (mean R = 0.52, ubRMSE = 0.047 m(3)/m(3)) and certain improvement to the ESA CCI SM data both at station and network scales. Compared with existing 1 km SM datasets, the LHS-SM data also showed the best performance (mean R = 0.62, ubRMSE = 0.047 m(3)/m(3)), while existing datasets either failed to fully characterize the spatial details or had some data gaps and unreasonable distributions. Strong spatial heterogeneity was observed in the SM dynamics during 2001-2020 with the southwest and northeast showing a dry gets wetter scheme and the southeast presenting a wet gets drier trend. Overall, the LHS-SM dataset gained its added values by compensating the drawbacks of existing 1 km SM products over the QTP and was much valuable for many regional applications.

Monitoring and modelling surface deformation are crucial components of understanding the freeze-thaw process and preventing disasters in permafrost regions. However, previous methods had limitations that inhibited the interpretation of freeze-thaw deformation, such as a lack of physical meaning, an inability to reflect the physical freeze-thaw process and consideration of only a single external factor's impact on permafrost deformation. This study proposes an improved degree-day model (IDM) for quantitatively isolating surface deformation using interferometric synthetic aperture radar (InSAR) technology over permafrost. We considered the effect of soil moisture variation on permafrost deformation and incorporated interannual variation in the freeze-thaw process due to climate change. By applying small baseline subset (SBAS) technology to Sentinel-1 InSAR measurements over the Wudaoliang permafrost region on the Qinghai-Tibet Plateau from 2018 to 2019, we estimated long-term and seasonal permafrost deformation. The reliability of InSAR results was validated using in situ measurements, with root mean square errors (RMSEs) less than 10 mm. The results showed that the average linear deformation rates in 2018 and 2019 were -3.8 mm a-1 and -11.0 mm a-1, respectively, and the maximum seasonal deformations were 15.7 mm and 13.2 mm, respectively. Compared with the original degree-day model (ODM), the method used in this study produced smaller residual deformations of 6.9 mm and 6.4 mm, highlighting its ability to improve a quantitative description of permafrost deformation.