Quantifying drought and wetness fluctuations is of great significance to the regional ecological environment and water resource security, especially in the fragile Qinghai-Tibetan Plateau (QTP). In this paper, the standardized precipitation evapotranspiration index (SPEI) was calculated based on the observed data and China Meteorological Forcing Dataset (CMFD) in the QTP for the period of 1979–2015, and the drought and wetness evolution based on the SPEI series and respective contribution of temperature and precipitation were also analyzed. Results indicated that meteorological stations are mainly concentrated in the eastern part of the plateau, which cannot reflect the drought and wetness trend of the whole QTP. The linear trend and Mann–Kendall test revealed that SPEI series calculated based on CMFD data at 1-, 3-, 6-, 9-, 12-, and 24-month time scales all showed significant upward trend , indicating that the QTP as a whole tended to be wetter. Spatially, the regions with significant drying and increased drought probability were mainly concentrated in the Qaidam Basin and the southern part of the QTP, and the mean contribution rates of temperature and precipitation variability to SPEI trend in these regions were 60% and −11%, respectively. The regions with significant wetting and decreased drought probability were mainly concentrated in the northeast, central, and western parts of the plateau, and the mean contribution rates of temperature and precipitation variability to SPEI trend were −9% and 61% in these regions. From the statistics in different climatic regions, most of the arid and humid regions in the QTP tended to be drier, while the semiarid regions tended to be wetter.

Understanding the tradeoffs/synergies between ecosystems is crucial to effective watershed ecosystem management and sustainable development. In this study, the Integrated Valuation of Ecosystem Services and Tradeoffs (InVEST) and Carnegie Ames Stanford Approach (CASA) models were utilized to estimate four ecosystem services (ESs), including water conservation capacity, soil retention, habitat quality, and carbon storage services, in Qinghai Lake Basin between 2000 and 2018. Local Indicators of Spatial Association (LISA) and tradeoffs/synergies criterion (TSC) were used to reveal the relationships between the ESs. The results show that the water conservation capacity, soil retention, habitat quality, and carbon storage service in Qinghai Lake Basin all increased between 2000 and 2018. TSC and LISA revealed that carbon storage and habitat quality, habitat quality and water conservation capacity, and carbon storage and soil retention had the same relationship, as did habitat quality and soil retention. In addition, LISA showed that the relationships between ESs are mainly based on high high clusters that concentrate in the middle of the basin. The analysis also revealed obvious spatial heterogeneity. This study aims to compensate the research deficiencies that affected previous studies of the Qinghai Lake Basin and provide a point of reference for the sustainable development of the basin.

Accurate simulation of snow cover process is of great significance to the study of climate change and the water cycle. In our study, the China Meteorological Forcing Dataset (CMFD) and ERA-Interim were used as driving data to simulate the dynamic changes in snow depth and snow water equivalent (SWE) in the Irtysh River Basin from 2000 to 2018 using the Noah-MP land surface model, and the simulation results were compared with the gridded dataset of snow depth at Chinese meteorological stations (GDSD), the long-term series of daily snow depth dataset in China (LSD), and China’s daily snow depth and snow water equivalent products (CSS). Before the simulation, we compared the combinations of four parameterizations schemes of Noah-MP model at the Kuwei site. The results show that the rainfall and snowfall (SNF) scheme mainly affects the snow accumulation process, while the surface layer drag coefficient (SFC), snow/soil temperature time (STC), and snow surface albedo (ALB) schemes mainly affect the melting process. The effect of STC on the simulation results was much higher than the other three schemes; when STC uses a fully implicit scheme, the error of simulated snow depth and snow water equivalent is much greater than that of a semi-implicit scheme. At the basin scale, the accuracy of snow depth modeled by using CMFD and ERA-Interim is higher than LSD and CSS snow depth based on microwave remote sensing. In years with high snow cover, LSD and CSS snow depth data are seriously underestimated. According to the results of model simulation, it is concluded that the snow depth and snow water equivalent in the north of the basin are higher than those in the south. The average snow depth, snow water equivalent, snow days, and the start time of snow accumulation (STSA) in the basin did not change significantly during the study period, but the end time of snow melting was significantly advanced.

The Huangshui National Wetland Park (HNWP) is a unique national wetland park in a city on the Qinghai-Tibetan Plateau, containing three zones: Haihu, Beichuan, and Ninghu. In this study, a total of 54 soil samples (18 sampling points with depths of 0–10 cm, 10–20 cm, and 20–30 cm) were collected in these three zones, and the contents of heavy metals (Cr, Cd, Cu, Hg, Ni, Pb, Zn, and As) of each sample were determined. The ecological risk of eight kinds of heavy metals was evaluated by using the geo–accumulation index ($I_{geo}$), and the ecological risk–controlling effect of the Xining urban wetlands on heavy metals was explored by comparative analysis, and the possible sources of heavy metals in the soil were analyzed via correlation analysis and principal component analysis (PCA). The results revealed that the total heavy metal concentration order was Haihu > Beichuan > Ninghu zone. As and Cu presented vertical accumulation characteristics in the surface and lower horizon, respectively. Cr, Cd, Hg, Ni, Pb, and Zn accumulated downwards along the depth. On the spatial scale, the enrichments of Cd and Hg brought non-negligible ecological risks in plateau urban wetlands. The results of PCA indicated that soil heavy metals mainly came from compound sources of domestic and atmospheric influences, traffic pollution sources, and industrial pollution sources. The study has revealed that human activities have inevitable negative impacts on wetland ecosystems, while the HNWP provides a significant weakening effect on heavy metal pollution.

High spatial–temporal resolution surface net radiation (RN) data are of great significance to the study of climate, ecology, hydrology and cryosphere changes on the Tibetan Plateau (TP), but the verification of the surface net radiation products on the plateau is not sufficient. In this study, the China Meteorological Administration Global Land Surface Reanalysis Products (CRA/Land) and ECMWF Land Surface Reanalysis version 5 (ERA5-Land) RN data were validated using ground measurements at daily and monthly time scales, and the spatiotemporal patterns were also analyzed. The results indicate the following: (1) CRA/Land overestimated while ERA5-Land underestimated RN, but CRA/Land RN outperformed ERA5-Land in observations at the daily and monthly scale. (2) The CRA/Land RN data had a larger error in the central part and a smaller error in the northeast of the TP, while ERA5-Land showed the opposite. (3) The spatial patterns of RN revealed by CRA/Land and ERA5-Land data showed differences in most regions. The CRA/Land data showed that the RN of the TP had a downward trend during 2000 and 2020 with a slope of −0.112 W·m⁻²/a, while the ERA5-Land data indicated an upward trend with a change rate of 0.016 W·m⁻²/a. (4) Downwelling shortwave radiation (DSR), upwelling shortwave radiation (USR), downwelling longwave radiation (DLR) and upwelling longwave radiation (ULR) are the four components of RN, and the evaluation results indicate that the DSR, DLR and ULR recorded via CRA/Land and ERA5-Land are consistent with the observed data, but the consistency between the USR recorded via CRA/Land and ERA5-Land and the observed data is poor. (5) The inconsistency of the USR data is the main reason for the large differences in the spatiotemporal distribution of CRA/Land and ERA5-Land RN data across the TP.

Alpine wetlands sequester large amounts of soil carbon, so it is vital to gain a full understanding of their land-atmospheric CO₂ exchanges and how they contribute to regional carbon neutrality; such an understanding is currently lacking for the Qinghai—Tibet Plateau (QTP), which is undergoing unprecedented climate warming. We analyzed two-year (2018–2019) continuous CO₂ flux data, measured by eddy covariance techniques, to quantify the carbon budgets of two alpine wetlands (Luanhaizi peatland (LHZ) and Xiaobohu swamp (XBH)) on the northeastern QTP. At an 8-day scale, boosted regression tree model-based analysis showed that variations in growing season CO₂ fluxes were predominantly determined by atmospheric water vapor, having a relative contribution of more than 65%. Variations in nongrowing season CO₂ fluxes were mainly controlled by site (categorical variable) and topsoil temperature (T_s), with cumulative relative contributions of 81.8%. At a monthly scale, structural equation models revealed that net ecosystem CO₂ exchange (NEE) at both sites was regulated more by gross primary productivity (GPP), than by ecosystem respiration (RES), which were both in turn directly controlled by atmospheric water vapor. The general linear model showed that variations in nongrowing season CO₂ fluxes were significantly (p < 0.001) driven by the main effect of site and T_s. Annually, LHZ acted as a net carbon source, and NEE, GPP, and RES were 41.5 ± 17.8, 631.5 ± 19.4, and 673.0 ± 37.2 g C/(m² year), respectively. XBH behaved as a net carbon sink, and NEE, GPP, and RES were –40.9 ± 7.5, 595.1 ± 15.4, and 554.2 ± 7.9 g C/(m² year), respectively. These distinctly different carbon budgets were primarily caused by the nongrowing season RES being approximately twice as large at LHZ (p < 0.001), rather than by other equivalent growing season CO₂ fluxes (p > 0.10). Overall, variations in growing season CO₂ fluxes were mainly controlled by atmospheric water vapor, while those of the nongrowing season were jointly determined by site attributes and soil temperatures. Our results highlight the different carbon functions of alpine peatland and alpine swampland, and show that nongrowing season CO₂ emissions should be taken into full consideration when upscaling regional carbon budgets. Current and predicted marked winter warming will directly stimulate increased CO₂ emissions from alpine wetlands, which will positively feedback to climate change.