Impacts of human pressure and climate on biodiversity–multifunctionality relationships on the Qinghai–Tibetan Plateau

Many studies have investigated the effects of environmental context on biodiversity or multifunctionality in alpine regions, but it is uncertain how human pressure and climate may affect their relationships. Here, we combined the comparative map profile method with multivariate datasets to assess the spatial pattern of ecosystem multifunctionality and further identify the effects of human pressure and climate on the spatial distribution of biodiversity–multifunctionality relationships in alpine ecosystems of the Qinghai–Tibetan Plateau (QTP). Our results indicate that at least 93% of the areas in the study region show a positive correlation between biodiversity and ecosystem multifunctionality across the QTP. Biodiversity–multifunctionality relationships with increasing human pressure show a decreasing trend in the forest, alpine meadow, and alpine steppe ecosystems, while an opposite pattern was found in the alpine desert steppe ecosystem. More importantly, aridity significantly strengthened the synergistic relationship between biodiversity and ecosystem multifunctionality in forest and alpine meadow ecosystems. Taken together, our results provide insights into the importance of protecting and maintaining biodiversity and ecosystem multifunctionality in response to climate change and human pressure in the alpine region.

1 Introduction

One of the primary directions for ecology is to determine how biodiversity is generated and maintained across broad macroecological patterns and temporal dynamics (Chase, 2010; Guo et al., 2019), which can enhance our overall awareness of the underlying mechanisms shaping biodiversity, the development of ecological theories, and the conservation of biodiversity (Zhou et al., 2008; Morlon et al., 2011; Guo et al., 2019). Most studies have demonstrated that biodiversity and ecosystem functions are closely related (Jing et al., 2015; Gross et al., 2017; Hu et al., 2021). More importantly, their relationships can further influence ecosystem services related to human well-being (Cameron et al., 2018), including the quantity and quality of the food we eat, the water we drink, and the air we breathe and the control of disease (Wall et al., 2015). However, with increasing human pressures and climate change (Yao et al., 2022), the relationships between biodiversity and multifunctionality may be subject to great uncertainty in response to these changes. There is therefore an urgent need to improve comprehensive information on biodiversity–multifunctionality relationships at spatially explicit scales in the context of global change.

At present, field experiments (Delgado-Baquerizo et al., 2020), belt transect surveys (Jing et al., 2015; Hu et al., 2021), and meta-analyses (Zhang R. et al., 2022) have shown that biodiversity plays an important role in driving ecosystem functioning at different trophic levels. The multiple lines of evidence suggest that the relationships between biodiversity and ecosystem functioning can be strongly predicted by the biodiversity of plant and soil microbiota (Maestre et al., 2012; Meyer et al., 2018; Albrecht et al., 2021) and that there is the potential for the synergistic and complementary influence of plant (Fanin et al., 2018; Valencia et al., 2018) and soil microbial diversity (Zhang and Zhang, 2016; Delgado-Baquerizo et al., 2020) on ecosystem functioning. For example, an extensive climatic gradient study on the Qinghai–Tibetan Plateau (QTP) showed that the observed differences in ecosystem multifunctionality between sites were largely (45%) due to the joint influence of above- and belowground biodiversity (Jing et al., 2015). As such, the aridity gradient study in northern China showed that soil multifunctionality has a strong positive association with plant species richness in less arid regions and with microbial diversity in more arid regions, respectively (Hu et al., 2021). Yet, in the context of human-induced global climate change, biodiversity and ecosystem multifunctionality may be subject to potential ecological consequences under uncertain environmental change (Wagg et al., 2014; Jing et al., 2015; Delgado-Baquerizo et al., 2020; Hu et al., 2021). For instance, human-induced land-use changes, such as desertification, urbanization, deforestation, and agriculture, can cause profound destruction of biodiversity and ecosystem functions (Haddad et al., 2015; Tsiafouli et al., 2015; Wall et al., 2015). Similarly, increased soil erosion caused by extreme water (Wall et al., 2015) and wind events (Garrison et al., 2003) can alter the intrinsic structure of dust storms and the distribution of soil organisms and pathogens (Quinton et al., 2010), with significant impacts on the decomposition function of ecosystems. Although these works provide multidimensional perspectives to demonstrate the importance of environmental context in maintaining biodiversity and ecosystem multifunctionality, there is still a lack of macro-ecological perspectives and determinations of the impacts of human pressure and climate on their relationships across specific ecological gradients, such as the QTP.

The QTP region is one of the world’s biodiversity hotspots and contains numerous endemic amphibians, mammals, birds, and plants, especially along the eastern Himalayas and the Hengduan Mountains (Fu and Shen, 2017; Deng et al., 2020; Yu et al., 2021). A previous report indicated that the Hengduan Mountains are home to 8,439 species of seed plants and that a total of 358 species of passerine birds can be found along the gradient of forest ecosystems that appears to be warm and humid in the eastern Himalayas (Päckert et al., 2020; Yu et al., 2020; Mao et al., 2021). These unique ecological conditions combined with the environmental background have made the QTP a “natural experimental site,” triggering a lot of scientific research, mainly on plant phenology (Shen et al., 2022), cryospheric environment transformation (Kang et al., 2010), water supply changes (Immerzeel et al., 2010), alpine forest growth (Silva et al., 2016), alpine ecosystem vegetation shifts (Zhang et al., 2013), ecosystem function (Zhang H. et al., 2022), and ecosystem services (Jiang et al., 2020). In comparison, although there is important research on the effects of climatic variation and biodiversity on ecosystem functions at the plot level (Jing et al., 2015; Wang Y. et al., 2021), there is a lack of studies describing the spatial distribution of ecosystem multifunctionality and biodiversity–multifunctionality relationships and how this relationship responds to human pressure and climate change along different ecological gradients.

Here, we used multivariate data (including total N, total P, total K, alkali-hydrolyzable N, available P, available K, cation exchange capacity, soil organic matter, carbon stock density of above- and belowground living biomass, soil microbial biomass carbon and nitrogen, vascular plant species richness, and hydrothermal conditions) based on the averaging approach, principal component analysis (PCA) approach, threshold approach, and Comparison Map Profile Method to assess the impacts of human pressure and climate on the spatial distribution of biodiversity–multifunctionality relationships across different ecosystems of the QTP. Our main objectives were to (1) assess the spatial distribution of ecosystem multifunctionality based on averaging, PCA, and threshold approaches across the QTP and validate the robustness of the assessment results; (2) identify the spatial distribution of biodiversity–multifunctionality relationships across the QTP; (3) assess the impacts of human pressure and climate on the spatial distribution of biodiversity–multifunctionality relationships across different ecosystems of the QTP. We provide insights that can guide future conservation and maintenance of biodiversity and ecosystem multifunctionality in response to climate change and human pressure in the alpine region.

2 Materials and methods

2.1 Study area

Being the largest area of ice and snow storage apart from the Arctic and Antarctic, the QTP possesses a mean elevation of over 4,000 m and an area of ~2.50 × 10⁶ km² (Zhang et al., 1999; Wu and Zhang, 2010; Kou et al., 2020), with a representative dry and cold climate, mainly influenced by the Indian summer monsoon and East Asian monsoon. Because of the gradual decrease of the East Asian monsoon with the altitudinal gradient from southeast to northwest, the climate pattern correspondingly represents the identical gradient change, with ranges of mean annual precipitation from 1,000 to 100 mm and ranges of mean annual temperature from 7.4°C to -4.1°C (Kou et al., 2020; Sun et al., 2020). These geographic variations of climate and topographic conditions conjointly characterize the particular ecosystem type on the QTP, spanning forest, alpine meadow, alpine steppe, and alpine desert steppe from southeast to northwest (Figure 1). As the primary ecosystem throughout the QTP, the alpine steppe is chiefly dominated by Stipa purpurea and Carex moorroftii, whereas in the alpine meadow, Kobresia pygmaea and Kobresia humilis are the dominant species (Ding et al., 2017). In brief, these exchanges of matter, energy, and organisms across different ecosystems shape the particularity of QTP around the world (Scherer-Lorenzen et al., 2022).

FIGURE 1

Figure 1 The distribution of ecosystem types used in this study. The illustration at the bottom left denotes the location of the Qinghai–Tibetan Plateau around the world. The points of different proportions represent the elevation changes of meteorological sites.

2.2 Climate and human pressure datasets

We employed the interpolation of precipitation and temperature to evaluate the effect of climate on biodiversity–multifunctionality relationships across the QTP. The interpolation of climate data (mainly including monthly mean air temperature and monthly precipitation sum) was conducted on the Anusplin 4.2 software (Centre for Resource and Environmental Studies, Australian National University, Canberra) by using the 107 meteorological sites (Figure 1) that were obtained from the China Meteorological Administration (http://cdc.cma.gov.cn) with a 5-km spatial resolution for the period 1982–2010. To evaluate the effect of human pressure on biodiversity–multifunctionality relationships across the geographical gradient of the QTP, we employed the index of global human modification to reflect it, which was obtained from the Global Human Modification of Terrestrial Systems (https://sedac.ciesin.columbia.edu/data/set/lulc-human-modification-terrestrial-systems) with a 1-km resolution and contained specific classification criterion: none (value = 0), low (0< value ≤ 0.1), moderate (0.1< value ≤ 0.4), high (0.4< value ≤ 0.7), very high (0.7< value ≤ 1), and uncategorized (value = NA).

2.3 Ecosystem function and biodiversity datasets

To assess the ecosystem functions provided by different sources, we obtained total N, total P, total K, alkali-hydrolyzable N, available P, available K, cation exchange capacity, and soil organic matter products from the Land-Atmosphere Interaction Research Group at Sun Yat-Sen University (http://globalchange.bnu.edu.cn/research/soil2), with a 30-arcsec resolution in a 0–0.045-m-depth soil. Meanwhile, we also acquired the above- and belowground living biomass carbon stock density from ORNL DAAC (https://daac.ornl.gov/VEGETATION/guides/Global_Maps_C_Density_2010.html) at a spatial resolution of 300 m. The soil microbial biomass carbon and nitrogen are freely available from the Compilation of Global Soil Microbial Biomass Carbon, Nitrogen, and Phosphorus Data (http://daac.ornl.gov/cgi-bin/dsviewer.pl?ds_id=1264) with a spatial resolution of 0.5°. In this study, all of the ecosystem functions are employed to represent how the ecosystem is functioning in relation to relevant ecosystem services, chiefly providing a proxy of ecosystem function linked to soil fertility, nutrient cycling, potassium cycling, phosphorus cycling, carbon cycling, the composition of carbon pools, the dynamics of organic decomposition, and nitrogen mineralization (Xu et al., 2013; Soliveres et al., 2014). Furthermore, to comprehend the potential level of plant biodiversity across the QTP, we obtained a vascular plant species richness product (Kreft and Jetz, 2007) from the Laboratory for Anthropogenic Landscape Ecology (https://ecotope.org/anthromes/biodiversity/plants/maps/) at a spatial resolution of 1°. All grid products described above were resampled to 5 km spatial resolution to facilitate comparison.

2.4 Assessment of ecosystem multifunctionality based on different approaches

To detect the integrated ability of different forms of ecosystem functions across the QTP, we employed averaging, PCA, and threshold approaches to assess the ecosystem’s multifunctionality. The averaging approach collects the average values of each standardized ecosystem function in each grid; here, we calculated the standardized values using the Z-scores as described in Gross et al. (2017). For the PCA approach, we calculated the PCA axis at a grid-scale using the Principal Component Toolbox of ArcMap 10.7 (Environmental Systems Research Institute, Inc., Redlands, CA, USA), based on the standardized indicators of ecosystem functions. We then summed the PCA axis scores weighted by the contribution rate of each axis in groups into a single index to reflect the integrated ability of ecosystem functions (Meyer et al., 2018). In addition to the two approaches described above, we also conducted the threshold approach to generate the spatial maps of ecosystem multifunctionality. For the local-scale studies, the multiple-threshold approach can be readily realized by the most significant relationship between biodiversity and multifunctionality (that is, the capacity of biodiversity to sustain multiple ecosystem functions at high levels). In contrast, given that the ecosystem’s multifunctionality on a widespread scale, there is unprecedented heterogeneity of threshold in each grid spanning the QTP, which will lead to the calculation of multifunctionality inapplicably using the multiple-threshold approach based on the ecosystem functions without time series; thus, we utilized the recommendation as described by Gamfeldt et al. (2008). They employed a single threshold to quantify ecosystem multifunctionality and selected 50% of the maximum for each ecosystem function as the threshold, following the next two reasons: (1) If an ecosystem can provide 50% of the maximum species for a function, thereby considering that the ecosystem can maintain this function by referencing the median lethal effect concentration in ecotoxicology; and (2) For any single function, on average one species was needed to maintain a specific process at ~60% of the maximum rate. Thus, selecting 50% as the threshold can maximize the probability, i.e., one function at least needs a species to sustain. It is important to note that the maximum represents the range at the top 5% of each ecosystem function, and that is to prevent outliers. After determining the threshold and maximum, we counted the total number of ecosystem functions that passed the prescribed threshold and thereby produced the spatial map of ecosystem multifunctionality with the threshold approach. Although we used three basic approaches to assess ecosystem multifunctionality at the grid scale, each approach has its strengths and weaknesses. To address this issue, we average the three indexes after normalization to synthetically represent ecosystem multifunctionality.

2.5 Characterizing aridity level

We employed the aridity index to introduce the aridity extent across the QTP, which was the ratio between precipitation and potential evapotranspiration (i.e., aridity index = precipitation/potential evapotranspiration). To quantize the potential evapotranspiration at grid scale across the QTP, we utilized the Penman–Monteith method based on the interpolation data of monthly mean air temperature and monthly precipitation sum. The detailed calculation process is well explained in a previous study (Vicente-Serrano et al., 2010). Finally, we adopted the value of 1 − aridity index to reflect the aridity level, where higher index values represent drier conditions, and vice versa (Berdugo et al., 2017).

2.6 Assessment of biodiversity–multifunctionality relationships at grid scale

In this study, we used the Comparison Map Profile Method to examine biodiversity-multifunctionality relationships at each grid. The Comparison Map Profile Method employs the averaged values of cross-correlation (CC) coefficient from scales 1 (9 × 9 pixel moving window) to 10 (99 × 99-pixel moving window), increasing by 10 increments to represent biodiversity–multifunctionality correlations; their computational formula is as follows (Yao et al., 2018; Ding et al., 2019):

$\begin{array}{ll}\text{CC}=\frac{1}{N^2}{\displaystyle \sum _{i=1}^N{\displaystyle \sum _{j=1}^N\frac{(x_{ij}-\overline{x})(y_{ij}-\overline{y})}{{\sigma }_x{\sigma }_y}}}& (1)\end{array}$

with

$\begin{array}{ll}{\sigma }_x^2=\frac{1}{N^2-1}{\displaystyle \sum _{i=1}^N{\displaystyle \sum _{j=1}^N{(x_{ij}-\overline{x})}^2}}& (2)\end{array}$

where $\overline{x}$ and $\overline{y}$ denote the averaged values of two moving windows for the two comparative images of biodiversity and ecosystem multifunctionality. The x_ij and y_ij represent the pixel value of two moving windows in two compared images at row i and column j, respectively. σ_x and σ_y denote the standard deviation corresponding to the two moving windows of two comparative images. The N is the pixel number of moving windows. Usually, high CC values illustrate that the two comparative images possess better synergic relationships, and vice versa.

2.7 Analyses

We first assessed the spatial distribution of ecosystem multifunctionality across the QTP and used the biodiversity dataset to verify the robustness of the results. To detect the relationship between plant species richness and ecosystem multifunctionality at grid scale across the QTP, we applied the Comparison Map Profile Method to quantify the correlation between native species richness (NSR) and multifunctionality index (MI). After obtaining the biodiversity–multifunctionality relationships at the grid scale, we then adopted the “GAM” of smooth function in the R 4.1.1 software to examine the variations of biodiversity–multifunctionality relationships with human pressure in different ecosystems. We also analyzed the data distribution characteristics of biodiversity–multifunctionality relationships, from none to the high extent of human pressure in different ecosystems. Given that very high human pressure is principally distributed in urban land, we did not count the distribution of biodiversity–multifunctionality relationships in this type. Finally, we conducted the bivariate correlations between biodiversity–multifunctionality relationships and the multiyear average of precipitation, temperature, or aridity in different ecosystems by using the “lm” function in the R 4.1.1 software.

3 Result

3.1 Spatial distribution of ecosystem multifunctionality and biodiversity–multifunctionality relationships

The spatial pattern of ecosystem multifunctionality shows that three indexes experienced relatively consistent geographic variations, i.e., with the high value primarily occurring in the southeast (≥ 0.60; Figures 2A–C) and the low value chiefly distributed in the northwest region (≤ 0.20; Figures 2A–C), whereas the ecosystem multifunctionality based on averaging approach owned fewer high-value areas when compared with PCA and threshold approaches in the southeast of QTP (Figure 2A). Additionally, the three ecosystem multifunctionality indexes displayed excellent robustness and experienced significant bivariate correlations (R² = 0.49, 0.49, and 0.48 for averaging, PCA, and threshold approaches, respectively; p < 0.05; Figure 2D) with plant species richness.

FIGURE 2

Figure 2 The spatial distribution of ecosystem multifunctionality based on different approaches and evaluation of robustness. (A–C) The spatial patterns of ecosystem multifunctionality calculated based on the averaging (AVE), principal component analysis (PCA), and threshold (THR) approaches on the Qinghai–Tibetan Plateau, respectively. (D) Bivariate correlations between ecosystem multifunctionality and native species richness of plants. Shading denotes the 95% confidence interval of the regression.

For the biodiversity–multifunctionality correlations at the grid scale across the QTP, we found that more than 90% of the region (93.44%; Figure 3A) possessed a synergic relationship (CC > 0). Meanwhile, there were considerable heterogeneous biodiversity–multifunctionality correlations that possessed high values (CC exceeding or equal to 0.3) chiefly concentrating in the central region of QTP (Figure 3A) and low values (CC ≤ 0.05) primarily distributed in the eastern region of QTP (Figure 3A). For the variations of biodiversity–multifunctionality relationships with elevation, we found a decrease in biodiversity–multifunctionality correlation with elevation ascent in the forest (Figure 3B), yet the alpine meadow experienced the opposite tendency (Figure 3C). Additionally, the alpine steppe and alpine desert steppe underwent no clear change for biodiversity–multifunctionality correlation across altitudinal gradient (Figures 3D, E).

FIGURE 3

Figure 3 The spatial distribution of biodiversity–multifunctionality relationships and its variation with elevation. (A) The spatial distribution of biodiversity–multifunctionality relationships across the Qinghai–Tibetan Plateau. (B–E) Variations of biodiversity–multifunctionality relationships with elevation in the forest, alpine meadow, alpine steppe, and alpine desert steppe ecosystems, respectively.

3.2 The effect of human pressure on biodiversity–multifunctionality relationships

Figure 4 shows the variations of biodiversity–multifunctionality relationships with human pressure in four representative ecosystems across the QTP. There were no obvious linear relationships between biodiversity–multifunctionality correlations and human pressure (Figure 4). The biodiversity–multifunctionality correlations in forest and alpine meadow ecosystems underwent a slight downward trend with the enhancement of human pressure (Figures 4A, C), whereas the alpine steppe ecosystem experienced a first downward and then upward trend change, but the magnitude of both trends was also small (Figure 4E). Compared to the forest and alpine meadow, biodiversity–multifunctionality correlations in alpine desert steppe exhibit an opposite trend change with increased human pressure; similarly (Figure 4G), the magnitude of change in this trend is also minor. For the distribution characteristics of biodiversity–multifunctionality correlations in none, low, moderate, and high degrees of human pressure, we found that this relationship in the forest, alpine meadow, and alpine desert steppe ecosystems displayed obvious gradient transformation from none to high levels of human pressure (Figures 4B, D, H), with weakening transformation in the forest (average values of correlations in corresponding degrees of human pressure were 0.14, 0.11, 0.09, and 0.08, respectively) and alpine meadow (average values of correlations in corresponding degrees of human pressure were 0.18, 0.16, 0.12, and 0.11, respectively) and enhanced transformation in the alpine desert steppe (average values of correlations in corresponding degrees of human pressure were 0.09, 0.12, 0.14, and 0.15, respectively). For the alpine steppe, although low human pressure possessed the strongest biodiversity–multifunctionality correlations, overall transformation trends decreased from none to high levels of human pressure (the average values of correlations in corresponding degrees of human pressure were 0.17, 0.18, 0.15, and 0.14, respectively; Figure 4F).

FIGURE 4

Figure 4 The human pressure in relation to biodiversity–multifunctionality correlations in different ecosystems. (A, C, E, G) Variations of biodiversity–multifunctionality relationships with increased human pressure in the forest, alpine meadow, alpine steppe, and alpine desert steppe ecosystems, respectively. (B, D, F, H) The statistics of biodiversity–multifunctionality relationships in different extents of human pressure for the forest, alpine meadow, alpine steppe, and alpine desert steppe ecosystems, respectively. Shading denotes the 95% confidence interval of the regression. The boxes denote the 10th and 90th percentiles, and the transverse lines in the boxes denote the average.

3.3 The effects of climate on the biodiversity–multifunctionality relationships

Expanding the analysis for variations of biodiversity–multifunctionality relationships with aridity, precipitation, and temperature found that biodiversity–multifunctionality correlations significantly increased with enhancement of aridity in the forest (average value of 1-AI was −0.67; Slope > 0, R²  = 0.22, p < 0.05; Figures 5A, M) and alpine meadow (average value of 1-AI was −0.54; Slope > 0, R² = 0.13, p < 0.05; Figures 5D, M). Compared to the first two ecosystems, there were no obvious changes in biodiversity–multifunctionality correlations with increased aridity in the alpine steppe (average value of 1-AI was 0.25; Slope < 0, R² < 0.05, p < 0.05; Figures 5G, M) and alpine desert steppe (average value of 1-AI was 0.64; Slope < 0, R² < 0.05, p < 0.05; Figures 5J, M). For precipitation, we found that the biodiversity–multifunctionality correlations significantly decreased with the increased precipitation in the forest (Slope < 0, R² = 0.21, p < 0.05; Figure 5B) and alpine meadow ecosystems (Slope < 0, R² = 0.19, p < 0.05; Figure 5E). By contrast, alpine steppe (Slope > 0, R² < 0.05, p < 0.05; Figure 5H) and alpine desert steppe ecosystems (Slope > 0, R² < 0.05, p < 0.05; Figure 5K) have weaker relationships between biodiversity–multifunctionality correlations and precipitation. For the variations of biodiversity–multifunctionality correlations with temperature, forest ecosystem exhibited a significant positive relationship between them (Slope > 0, R² = 0.10, p < 0.05; Figure 5C), whereas alpine meadow experienced the opposite trend (Slope < 0, R² = 0.09, p < 0.05; Figure 5F). Same as precipitation, biodiversity–multifunctionality correlations had no obvious variations with temperature increase in the alpine steppe (Slope < 0, R² < 0.05, p < 0.05; Figure 5I) and alpine desert steppe ecosystems (Slope > 0, R² < 0.05, p < 0.05; Figure 5L).

FIGURE 5

Figure 5 The aridity, precipitation, and temperature in relation to biodiversity–multifunctionality correlations in different ecosystems. (A, D, G, J) The bivariate relationships between aridity and biodiversity–multifunctionality correlations in the forest, alpine meadow, alpine steppe, and alpine desert steppe ecosystems, respectively. (B, E, H, K) The bivariate relationships between the multiyear average of precipitation and biodiversity–multifunctionality correlations in corresponding ecosystems. (C, F, I, L) The bivariate relationships between the multiyear average of temperature and biodiversity–multifunctionality correlations in corresponding ecosystems. Shading denotes the 95% confidence interval of the regression. (M) The distribution characteristics of aridity in corresponding ecosystems. Dark red vertical lines denote the average value of aridity (1-AI).

4 Discussion

4.1 Spatial patterns of biodiversity–multifunctionality relationships and ecosystem multifunctionality across the QTP

Multiple pieces of evidence have demonstrated that biodiversity is positively correlated with ecosystem multifunctionality, even though the potential mechanisms by which biodiversity is related to ecosystem function remain unclear (Ma et al., 2010; Adler et al., 2011; Jing et al., 2015). Similar results were identified across the QTP; our study confirmed that more than 90% of the region possessed positive biodiversity–multifunctionality relationships (Figure 3A), extending previous findings to the grid scale. More importantly, we found that the central region of QTP (mainly involving the Sanjiangyuan region) maintained the highest synergic relationship between biodiversity and ecosystem multifunctionality (Figure 3A). This reflects the fact that long-term ecological conservation has a positive effect on improving biodiversity–multifunctionality relationships (Li et al., 2013; Li et al., 2020). Compared with the central region, the eastern region of QTP (mainly involving northwestern Sichuan) has undergone negative biodiversity–multifunctionality relationships (Figure 3A). This might be due to the expansion of local industries and land use (such as land degradation as a result of overgrazing in this region) triggered by past economic growth models (Zhou and Peng, 2021), disrupting the relationship between biodiversity and ecosystem multifunctionality. In addition to the relationships between biodiversity and ecosystem multifunctionality across different regional scales (Delgado-Baquerizo et al., 2020; Hu et al., 2021), the geographic variation of environmental context influencing the spatial pattern of ecosystem multifunctionality also needs further attention. Our results indicate that the spatial distribution of ecosystem multifunctionality has undergone a distinct geographical gradient from forest to alpine desert steppe across the QTP (Figures 2A–C), as theoretically expected. This finding provides evidence linking the underlying impact of an extensive climatic gradient (with a decreasing trend from forest to alpine desert steppe) across the QTP on the spatial pattern of ecosystem multifunctionality. Since previous works suggest that the geographic variation of climate is one of the main drivers of growth and species distributions for plants along the altitudinal gradient (Sun et al., 2014; Chen J. et al., 2020; Yu et al., 2021), the extent of plant biodiversity will influence the functioning of the ecosystem (Pugnaire et al., 2019); thus, ecosystem multifunctionality is spatially patterned in the same way as climate. Likewise, the climate will also control the matter turnover rates. For example, there are many snow-capped mountains on the alpine steppe and alpine desert steppe of QTP, leading to the turnover efficiency of soil carbon being limited by the persistent low-temperature environment (Ding et al., 2019), thereby determining the low ecosystem multifunctionality principally concentrated in this region. Unlike the alpine steppe and alpine desert steppe, QTP’s forest and alpine meadow were able to maintain a suitable temperature and an adequate water input that will facilitate the normal functioning of ecosystem function (Ban et al., 2022; Li et al., 2022; Lu et al., 2022), ultimately contributing to the maintenance of ecosystem multifunctionality.

4.2 The influence of human pressure and climate on biodiversity–multifunctionality relationships

According to the biodiversity-ecosystem functioning concept, species richness and other facets of biodiversity have a significant influence on ecosystem processes (Scherer-Lorenzen et al., 2022). Although numerous studies have been conducted on the interaction between biodiversity and ecosystem functioning (Soliveres et al., 2014; Jing et al., 2015; Meyer et al., 2018; Delgado-Baquerizo et al., 2020; Hu et al., 2021), these studies have largely focused on the identification of the biodiversity–ecosystem relationship, ignoring the probable disturbance of human pressure on biodiversity–multifunctionality relationships in the alpine region. Fortunately, our study showed that the influence of human pressures on biodiversity–multifunctionality relationships was nonsignificant in all ecosystems of QTP (Figure 4). This might be due to policy implementation to control ecological degradation and enhance human well-being (Chen et al., 2014; Pan et al., 2016; Pan et al., 2017; Wang and Wei, 2022), such as the Regional Ecological Construction and Environmental Protection Plan on the QTP, Tibet Ecological Protection Barrier Protection and Construction Plan, and Overall Plan for Ecological Protection and Construction of Qinghai Sanjiangyuan Nature Reserve (Wang and Wei, 2022), leading to human activities being mainly concentrated in urban and scenic areas and the discontinuous variation of human pressures across the QTP (Wei et al., 2022), thereby offsetting the negative influence of human pressures on biodiversity–multifunctionality relationships on the QTP (Figures 4B, D, F) and ultimately generating the ambiguous pattern between human pressure and biodiversity–multifunctionality relationships (Figure 4). The variation of biodiversity–multifunctionality relationships with different degrees of human pressure in alpine desert steppe ecosystem further supplements the uncertain pattern. We found the synchronicity between biodiversity–multifunctionality correlations and enhancement of human pressure (Figures 4G, H), indicating that positive human activities (such as strengthening of management to combat land degradation and desertification) may improve the background conditions in desert areas and enhance synergic relationships between biodiversity and multifunctionality (Liu et al., 2021; You et al., 2021).

In contrast to human pressure, aridity significantly strengthened the synergic relationship between biodiversity and ecosystem multifunctionality in the forest and alpine meadow ecosystems (Figures 5A, D). In the forest, the biodiversity–multifunctionality relationship was significantly negatively and positively related to precipitation (Figure 5B) and temperature (Figure 5C), respectively. It is possible that there is a steeper decline in ecosystem function, such as soil microorganism activity weakening with increasing soil moisture in forests that receive ample precipitation all year round (Schafer and Kotanen, 2003; Zhang and Zhang, 2016; Chen Y. et al., 2020), possibly leading to a large amount of organic matter not being decomposed by soil microorganisms due to respiration limitations (Medrano et al., 2015; Trugman et al., 2019), thereby exhibiting this consequence of negative correlation between biodiversity–multifunctionality relationship and precipitation (Figure 5B). This can be explained by the fact that biodiversity–multifunctionality correlations gradually increased with temperature (Figure 5C), because warming can both evaporate excess water and enhance the decomposition capacity of soil microorganisms to produce more compounds that can be more easily absorbed by vegetation (Chou and Lan, 2012; Levy et al., 2019; Sun et al., 2020). The positive association of aridity with geographic variation of biodiversity–multifunctionality relationship (Figure 5A) also further contributes to answering the significant question that excess water (Figure 5M) is adverse to ecosystem function performance (such as mineralization and nitrification of nitrogen) and biodiversity maintenance in the forest (Zhao and Li, 2017; Pugnaire et al., 2019; Sun et al., 2020; Ye et al., 2020; Sun et al., 2022). Indeed, the strong positive correlation between the biodiversity–multifunctionality relationship and aridity (Figure 5D) in the alpine meadow is also possibly attributed to the obstacle of excess water (Figures 5E, M). For example, there are some alpine meadow areas that are always under waterlogged conditions, which could lower the soil’s oxygen content and increase the mortality of soil microorganisms (Wang X. et al., 2021; Sun et al., 2022), decreasing the potential for ecosystem function complementarity and, as a result, weakening the biodiversity–multifunctionality relationship in the more humid regions (Figures 5D, M).

4.3 Uncertainties and future directions

Our study underscored the obvious fact that aridity significantly affects the geographic variation of biodiversity–multifunctionality relationships across the forest and alpine meadow ecosystems of QTP. In contrast to links between human pressure or climate and biodiversity–multifunctionality relationships in this study, actual interactions between matter, energy, and organisms are dependent on the temporal and spatial dynamics of the change process in species composition and anthropogenic and natural disturbances across coupled ecosystems (Ratcliffe et al., 2017; Gounand et al., 2018; Ehbrecht et al., 2021; Scherer-Lorenzen et al., 2022), which are not reflected in our spatial analysis of functional links between them. Acknowledging this, future studies should consider potential mechanisms of exchange between matter, energy, and organisms at a plot level because our comprehensive understanding of how different ecosystems respond to changes in resource availability can be greatly improved by reflecting the ecological process in this way (Kulmatiski et al., 2008; Pugnaire et al., 2019). While our findings support much of the mechanistic and experimental research related to the spatial link between ecosystem multifunctionality and biodiversity (Berdugo et al., 2017; Ehbrecht et al., 2021; Ma et al., 2021; Moi et al., 2022), we must emphasize that cross-correlations are not necessarily a sign of causality. In addition, the currently ongoing global pattern of vascular plant diversity provides a preliminary assessment of the potential level of plant diversity. Future work is therefore required to trigger the development of new methodologies and metrics to assess these proxies of plant and soil biodiversity on our planet. Finally, more meteorological sites are not available; therefore, future study should take into account relating the actual to high-precision satellite products.

5 Conclusions

By reflecting extensive changes in the geographic gradient of ecosystem multifunctionality with climatic variation across different ecosystems based on multivariate datasets, we report results that show the spatial distribution of biodiversity–multifunctionality relationships at the alpine scale can inform our comprehensive understanding of species distribution and dynamics. Our results suggest that biodiversity and multifunctionality are synergistically related in most areas of the QTP. We have also shown that human pressure has a negative impact on biodiversity–multifunctionality relationships in the forest, alpine meadow, and alpine steppe ecosystems, but with large contrasts in the alpine desert steppe ecosystem, which contributes to the important implications for desertification control that need to be further implemented in the QTP. An important finding of our study is that there is a significant positive association between aridity and biodiversity–multifunctionality relationships in forest and alpine meadow ecosystems. These results demonstrate that climate has a more significant impact on biodiversity–multifunctionality relationships relative to human pressure and the urgency of targeted management actions to monitor the influence of climate change on biodiversity and ecosystem functions across the QTP, especially in forest and alpine meadow ecosystems.

Data availability statement

The original contributions presented in the study are included in the article/supplementary material. Further inquiries can be directed to the corresponding author.

Author contributions

CY: Conceptualization (lead); formal analysis (lead); methodology (lead); visualization (lead); writing – original draft (lead); writing – review and editing (lead). SW: Funding acquisition (lead); Methodology (supporting); visualization (supporting); writing–review and editing (supporting). YW: Methodology (supporting); validation (supporting); writing–review and editing (supporting). TZ: Methodology (supporting); validation (supporting); writing–review and editing (supporting). RL: Methodology (supporting). All authors contributed to the article and approved the submitted version.

Funding

For financial support, we thank the Second Tibetan Plateau Scientific Expedition and Research Program (Grant No. 2019QZKK0405), the Strategic Priority Research Program of the Chinese Academy of Sciences (Grant No. XDA20020402), and the Fundamental Research Funds for the Central Universities of China.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The reviewer GF declared a shared affiliation with the author TZ to the handling editor at the time of review.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

Adler, P. B., Seabloom, E. W., Borer, E. T., Hillebrand, H., Hautier, Y., Hector, A., et al. (2011). Productivity is a poor predictor of plant species richness. Science 333, 1750–1753. doi: 10.1126/science.1204498

PubMed Abstract | CrossRef Full Text | Google Scholar

Albrecht, J., Peters, M. K., Becker, J. N., Behler, C., Classen, A., Ensslin, A., et al. (2021). Species richness is more important for ecosystem functioning than species turnover along an elevational gradient. Nat. Ecol. Evol. 5, 1582–1593. doi: 10.1038/s41559-021-01550-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Ban, C., Xu, Z., Zuo, D., Zhang, R., Chen, H., Ye, C., et al. (2022). Impact of variability in the hydrological cycle components on vegetation growth in an alpine basin of the southeastern Tibet plateau, China. Hydrol. Res. 53, 124–140. doi: 10.2166/nh.2021.086

CrossRef Full Text | Google Scholar

Berdugo, M., Kefi, S., Soliveres, S., Maestre, F. T. (2017). Plant spatial patterns identify alternative ecosystem multifunctionality states in global drylands. Nat. Ecol. Evol. 1, 1–10. doi: 10.1038/s41559-016-0003

PubMed Abstract | CrossRef Full Text | Google Scholar

Cameron, E. K., Martins, I. S., Lavelle, P., Mathieu, J., Tedersoo, L., Gottschall, F., et al. (2018). Global gaps in soil biodiversity data. Nat. Ecol. Evol. 2, 1042–1043. doi: 10.1038/s41559-018-0573-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Chase, J. M. (2010). Stochastic community assembly causes higher biodiversity in more productive environments. Science 328, 1388–1391. doi: 10.1126/science.1187820

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Y., Ma, S., Jiang, H., Hu, Y., Lu, X. (2020). Influences of litter diversity and soil moisture on soil microbial communities in decomposing mixed litter of alpine steppe species. Geoderma 377, 114577. doi: 10.1016/j.geoderma.2020.114577

CrossRef Full Text | Google Scholar

Chen, J., Yan, F., Lu, Q. (2020). Spatiotemporal variation of vegetation on the qinghai–Tibet plateau and the influence of climatic factors and human activities on vegetation trend, (2000–2019). Remote Sens. 12, 3150. doi: 10.3390/rs12193150

CrossRef Full Text | Google Scholar

Chen, B., Zhang, X., Tao, J., Wu, J., Wang, J., Shi, P., et al. (2014). The impact of climate change and anthropogenic activities on alpine grassland over the qinghai-Tibet plateau. Agric. For. Meteorol. 189, 11–18. doi: 10.1016/j.agrformet.2014.01.002

CrossRef Full Text | Google Scholar

Chou, C., Lan, C.-W. (2012). Changes in the annual range of precipitation under global warming. J. Climate 25, 222–235. doi: 10.1175/JCLI-D-11-00097.1

CrossRef Full Text | Google Scholar

Delgado-Baquerizo, M., Reich, P. B., Trivedi, C., Eldridge, D. J., Abades, S., Alfaro, F. D., et al. (2020). Multiple elements of soil biodiversity drive ecosystem functions across biomes. Nat. Ecol. Evol. 4, 210–220. doi: 10.1038/s41559-019-1084-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Deng, T., Wu, F., Zhou, Z., Su, T. (2020). Tibetan Plateau: an evolutionary junction for the history of modern biodiversity. Sci. China Earth Sci. 63, 172–187. doi: 10.1007/s11430-019-9507-5

CrossRef Full Text | Google Scholar

Ding, J., Chen, L., Ji, C., Hugelius, G., Li, Y., Liu, L., et al. (2017). Decadal soil carbon accumulation across Tibetan permafrost regions. Nat. Geosci. 10, 420–424. doi: 10.1038/ngeo2945

CrossRef Full Text | Google Scholar

Ding, J., Wang, T., Piao, S., Smith, P., Zhang, G., Yan, Z., et al. (2019). The paleoclimatic footprint in the soil carbon stock of the Tibetan permafrost region. Nat. Commun. 10, 4195. doi: 10.1038/s41467-019-12214-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Ehbrecht, M., Seidel, D., Annighofer, P., Kreft, H., Kohler, M., Zemp, D. C., et al. (2021). Global patterns and climatic controls of forest structural complexity. Nat. Commun. 12, 519. doi: 10.1038/s41467-020-20767-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Fanin, N., Gundale, M. J., Farrell, M., Ciobanu, M., Baldock, J. A., Nilsson, M.-C., et al. (2018). Consistent effects of biodiversity loss on multifunctionality across contrasting ecosystems. Nat. Ecol. Evol. 2, 269–278. doi: 10.1038/s41559-017-0415-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Fu, G., Shen, Z. X. (2017). Clipping has stronger effects on plant production than does warming in three alpine meadow sites on the northern Tibetan plateau. Sci. Rep. 7, 16330. doi: 10.1038/s41598-017-16645-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Gamfeldt, L., Hillebrand, H., Jonsson, P. R. (2008). Multiple functions increase the importance of biodiversity for overall ecosystem functioning. Ecology 89, 1223–1231. doi: 10.1890/06-2091.1

PubMed Abstract | CrossRef Full Text | Google Scholar

Garrison, V. H., Shinn, E. A., Foreman, W. T., Griffin, D. W., Holmes, C. W., Kellogg, C. A., et al. (2003). African And Asian dust: from desert soils to coral reefs. BioScience 53, 469–480. doi: 10.1641/0006-3568(2003)053[0469:AAADFD]2.0.CO;2

CrossRef Full Text | Google Scholar

Gounand, I., Little, C. J., Harvey, E., Altermatt, F. (2018). Cross-ecosystem carbon flows connecting ecosystems worldwide. Nat. Commun. 9, 1–8. doi: 10.1038/s41467-018-07238-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Gross, N., Bagousse-Pinguet, Y. L., Liancourt, P., Berdugo, M., Gotelli, N. J., Maestre, F. T. (2017). Functional trait diversity maximizes ecosystem multifunctionality. Nat. Ecol. Evol. 1, 1–9. doi: 10.1038/s41559-017-0132

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, X., Zhou, X., Hale, L., Yuan, M., Ning, D., Feng, J., et al. (2019). Climate warming accelerates temporal scaling of grassland soil microbial biodiversity. Nat. Ecol. Evol. 3, 612–619. doi: 10.1038/s41559-019-0848-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Haddad, N. M., Brudvig, L. A., Clobert, J., Davies, K. F., Gonzalez, A., Holt, R. D., et al. (2015). Habitat fragmentation and its lasting impact on earth’s ecosystems. Sci. Adv. 1, e1500052. doi: 10.1126/sciadv.1500052

PubMed Abstract | CrossRef Full Text | Google Scholar

Hu, W., Ran, J., Dong, L., Du, Q., Ji, M., Yao, S., et al. (2021). Aridity-driven shift in biodiversity–soil multifunctionality relationships. Nat. Commun. 12, 5350. doi: 10.1038/s41467-021-25641-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Immerzeel, W. W., Van Beek, L. P., Bierkens, M. F. (2010). Climate change will affect the Asian water towers. Science 328, 1382–1385. doi: 10.1126/science.1183188

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, W., Lü, Y., Liu, Y., Gao, W. (2020). Ecosystem service value of the qinghai-Tibet plateau significantly increased during 25 years. Ecosystem Serv. 44, 101146. doi: 10.1016/j.ecoser.2020.101146

CrossRef Full Text | Google Scholar

Jing, X., Sanders, N. J., Shi, Y., Chu, H., Classen, A. T., Zhao, K., et al. (2015). The links between ecosystem multifunctionality and above-and belowground biodiversity are mediated by climate. Nat. Commun. 6, 8159. doi: 10.1038/ncomms9159

PubMed Abstract | CrossRef Full Text | Google Scholar

Kang, S., Xu, Y., You, Q., Flügel, W.-A., Pepin, N., Yao, T. (2010). Review of climate and cryospheric change in the Tibetan plateau. Environ. Res. Lett. 5, 015101. doi: 10.1088/1748-9326/5/1/015101

CrossRef Full Text | Google Scholar

Kou, D., Yang, G., Li, F., Feng, X., Zhang, D., Mao, C., et al. (2020). Progressive nitrogen limitation across the Tibetan alpine permafrost region. Nat. Commun. 11, 1–9. doi: 10.1038/s41467-020-17169-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Kreft, H., Jetz, W. (2007). Global patterns and determinants of vascular plant diversity. Proc. Natl. Acad. Sci. 104, 5925–5930. doi: 10.1073/pnas.0608361104

CrossRef Full Text | Google Scholar

Kulmatiski, A., Beard, K. H., Stevens, J. R., Cobbold, S. M. (2008). Plant-soil feedbacks: a meta-analytical review. Ecol. Lett. 11, 980–992. doi: 10.1111/j.1461-0248.2008.01209.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Levy, O., Dayan, T., Porter, W. P., Kronfeld-Schor, N. (2019). Time and ecological resilience: can diurnal animals compensate for climate change by shifting to nocturnal activity? Ecol. Monogr. 89, e01334. doi: 10.1002/ecm.1334

CrossRef Full Text | Google Scholar

Li, C., Sun, H., Liu, L., Dou, T., Zhou, M., Li, W., et al. (2022). The importance of permafrost in the steady and fast increase in net primary production of the grassland on the qinghai–Tibet plateau. CATENA 211, 105964. doi: 10.1016/j.catena.2021.105964

CrossRef Full Text | Google Scholar

Li, J., Yin, H., Wang, D., Jiagong, Z., Lu, Z. (2013). Human-snow leopard conflicts in the sanjiangyuan region of the Tibetan plateau. Biol. Conserv. 166, 118–123. doi: 10.1016/j.biocon.2013.06.024

CrossRef Full Text | Google Scholar

Li, S., Zhang, H., Zhou, X., Yu, H., Li, W. (2020). Enhancing protected areas for biodiversity and ecosystem services in the qinghai–Tibet plateau. Ecosystem Serv. 43, 101090. doi: 10.1016/j.ecoser.2020.101090

CrossRef Full Text | Google Scholar

Liu, Y., Ren, H., Zheng, C., Zhou, R., Hu, T., Yang, P., et al. (2021). Untangling the effects of management measures, climate and land use cover change on grassland dynamics in the qinghai–Tibet plateau, China. Land Degradation Dev. 32, 4974–4987. doi: 10.1002/ldr.4084

CrossRef Full Text | Google Scholar

Lu, Y., Gao, Y., Dungait, J. A., Jia, J., Sun, K., Wang, S., et al. (2022). Understanding how inland lake system environmental gradients on the qinghai-Tibet plateau impact the geographical patterns of carbon and water sources or sink. J. Hydrol. 604, 127219. doi: 10.1016/j.jhydrol.2021.127219

CrossRef Full Text | Google Scholar

Ma, W., He, J. S., Yang, Y., Wang, X., Liang, C., Anwar, M., et al. (2010). Environmental factors covary with plant diversity–productivity relationships among Chinese grassland sites. Global Ecol. Biogeography 19, 233–243. doi: 10.1111/j.1466-8238.2009.00508.x

CrossRef Full Text | Google Scholar

Ma, H., Mo, L., Crowther, T. W., Maynard, D. S., Van Den Hoogen, J., Stocker, B. D., et al. (2021). The global distribution and environmental drivers of aboveground versus belowground plant biomass. Nat. Ecol. Evol. 5, 1110–1122. doi: 10.1038/s41559-021-01485-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Maestre, F. T., Quero, J. L., Gotelli, N. J., Escudero, A., Ochoa, V., Delgado-Baquerizo, M., et al. (2012). Plant species richness and ecosystem multifunctionality in global drylands. Science 335, 214–218. doi: 10.1126/science.1215442

PubMed Abstract | CrossRef Full Text | Google Scholar

Mao, K. S., Wang, Y., Liu, J. Q. (2021). Evolutionary origin of species diversity on the qinghai–Tibet plateau. J. Systematics Evol. 59, 1142–1158. doi: 10.1111/jse.12809

CrossRef Full Text | Google Scholar

Medrano, H., Tomás, M., Martorell, S., Flexas, J., Hernández, E., Rosselló, J., et al. (2015). From leaf to whole-plant water use efficiency (WUE) in complex canopies: limitations of leaf WUE as a selection target. Crop J. 3, 220–228. doi: 10.1016/j.cj.2015.04.002

CrossRef Full Text | Google Scholar

Meyer, S. T., Ptacnik, R., Hillebrand, H., Bessler, H., Buchmann, N., Ebeling, A., et al. (2018). Biodiversity-multifunctionality relationships depend on identity and number of measured functions. Nat. Ecol. Evol. 2, 44–49. doi: 10.1038/s41559-017-0391-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Moi, D. A., Lansac-Tôha, F. M., Romero, G. Q., Sobral-Souza, T., Cardinale, B. J., Kratina, P., et al. (2022). Human pressure drives biodiversity–multifunctionality relationships in large Neotropical wetlands. Nat. Ecol. Evol. 6 (9), 1279–1289. doi: 10.1038/s41559-022-01827-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Morlon, H., Schwilk, D. W., Bryant, J. A., Marquet, P. A., Rebelo, A. G., Tauss, C., et al. (2011). Spatial patterns of phylogenetic diversity. Ecol. Lett. 14, 141–149. doi: 10.1111/j.1461-0248.2010.01563.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Päckert, M., Favre, A., Schnitzler, J., Martens, J., Sun, Y. H., Tietze, D. T., et al. (2020). "Into and out of" the qinghai-Tibet plateau and the Himalayas: centers of origin and diversification across five clades of Eurasian montane and alpine passerine birds. Ecol. Evol. 10, 9283–9300. doi: 10.1002/ece3.6615

PubMed Abstract | CrossRef Full Text | Google Scholar

Pan, Y., Yu, C., Zhang, X., Chen, B., Wu, J., Tu, Y., et al. (2016). A modified framework for the regional assessment of climate and human impacts on net primary productivity. Ecol. Indic. 60, 184–191. doi: 10.1016/j.ecolind.2015.06.017

CrossRef Full Text | Google Scholar

Pan, T., Zou, X., Liu, Y., Wu, S., He, G. (2017). Contributions of climatic and non-climatic drivers to grassland variations on the Tibetan plateau. Ecol. Eng. 108, 307–317. doi: 10.1016/j.ecoleng.2017.07.039

CrossRef Full Text | Google Scholar

Pugnaire, F. I., Morillo, J. A., Peñuelas, J., Reich, P. B., Bardgett, R. D., Gaxiola, A., et al. (2019). Climate change effects on plant-soil feedbacks and consequences for biodiversity and functioning of terrestrial ecosystems. Sci. Adv. 5, eaaz1834. doi: 10.1126/sciadv.aaz1834

PubMed Abstract | CrossRef Full Text | Google Scholar

Quinton, J. N., Govers, G., Van Oost, K., Bardgett, R. D. (2010). The impact of agricultural soil erosion on biogeochemical cycling. Nat. Geosci. 3, 311–314. doi: 10.1038/ngeo838

CrossRef Full Text | Google Scholar

Ratcliffe, S., Wirth, C., Jucker, T., van der Plas, F., Scherer-Lorenzen, M., Verheyen, K., et al. (2017). Biodiversity and ecosystem functioning relations in European forests depend on environmental context. Ecol. Lett. 20, 1414–1426. doi: 10.1111/ele.12849

PubMed Abstract | CrossRef Full Text | Google Scholar

Schafer, M., Kotanen, P. M. (2003). The influence of soil moisture on losses of buried seeds to fungi. Acta Oecol. 24, 255–263. doi: 10.1016/j.actao.2003.09.001

CrossRef Full Text | Google Scholar

Scherer-Lorenzen, M., Gessner, M. O., Beisner, B. E., Messier, C., Paquette, A., Petermann, J. S., et al (2022). Pathways for cross-boundary effects of biodiversity on ecosystem functioning. Trends Ecol. Evol. 37 (5), 454–467. doi: 10.1016/j.tree.2021.12.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Shen, M., Wang, S., Jiang, N., Sun, J., Cao, R., Ling, X., et al. (2022). Plant phenology changes and drivers on the qinghai–Tibetan plateau. Nat. Rev. Earth Environ. 3, 633–651. doi: 10.1038/s43017-022-00317-5

CrossRef Full Text | Google Scholar

Silva, L. C., Sun, G., Zhu-Barker, X., Liang, Q., Wu, N., Horwath, W. R. (2016). Tree growth acceleration and expansion of alpine forests: the synergistic effect of atmospheric and edaphic change. Sci. Adv. 2, e1501302. doi: 10.1126/sciadv.1501302

PubMed Abstract | CrossRef Full Text | Google Scholar

Soliveres, S., Maestre, F. T., Eldridge, D. J., Delgado-Baquerizo, M., Quero, J. L., Bowker, M. A., et al. (2014). Plant diversity and ecosystem multifunctionality peak at intermediate levels of woody cover in global drylands. Glob. Ecol. Biogeogr. 23, 1408–1416. doi: 10.1111/geb.12215

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, H., Niu, Y., Chen, Y. S., Song, B., Liu, C. Q., Peng, D. L., et al. (2014). Survival and reproduction of plant species in the qinghai–Tibet plateau. J. Systematics Evol. 52, 378–396. doi: 10.1111/jse.12092

CrossRef Full Text | Google Scholar

Sun, J., Ye, C., Liu, M., Wang, Y., Chen, J., Wang, S., et al. (2022). Response of net reduction rate in vegetation carbon uptake to climate change across a unique gradient zone on the Tibetan plateau. Environ. Res. 203, 111894. doi: 10.1016/j.envres.2021.111894

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, J., Zhou, T. C., Liu, M., Chen, Y. C., Liu, G. H., Xu, M., et al. (2020). Water and heat availability are drivers of the aboveground plant carbon accumulation rate in alpine grasslands on the Tibetan plateau. Global Ecol. Biogeogr. 29, 50–64. doi: 10.1111/geb.13006

CrossRef Full Text | Google Scholar

Trugman, A. T., Anderegg, L. D., Sperry, J. S., Wang, Y., Venturas, M., Anderegg, W. R. (2019). Leveraging plant hydraulics to yield predictive and dynamic plant leaf allocation in vegetation models with climate change. Global Change Biol. 25, 4008–4021. doi: 10.1111/gcb.14814

CrossRef Full Text | Google Scholar

Tsiafouli, M. A., Thébault, E., Sgardelis, S. P., De Ruiter, P. C., van der Putten, W. H., Birkhofer, K., et al. (2015). Intensive agriculture reduces soil biodiversity across Europe. Global Change Biol. 21, 973–985. doi: 10.1111/gcb.12752

CrossRef Full Text | Google Scholar

Valencia, E., Gross, N., Quero, J. L., Carmona, C. P., Ochoa, V., Gozalo, B., et al. (2018). Cascading effects from plants to soil microorganisms explain how plant species richness and simulated climate change affect soil multifunctionality. Global Change Biol. 24, 5642–5654. doi: 10.1111/gcb.14440

CrossRef Full Text | Google Scholar

Vicente-Serrano, S. M., Beguería, S., López-Moreno, J. I. (2010). A multiscalar drought index sensitive to global warming: the standardized precipitation evapotranspiration index. J. Climate 23, 1696–1718. doi: 10.1175/2009JCLI2909.1

CrossRef Full Text | Google Scholar

Wagg, C., Bender, S. F., Widmer, F., van der Heijden, M. G. (2014). Soil biodiversity and soil community composition determine ecosystem multifunctionality. Proc. Natl. Acad. Sci. 111, 5266–5270. doi: 10.1073/pnas.1320054111

CrossRef Full Text | Google Scholar

Wall, D. H., Nielsen, U. N., Six, J. (2015). Soil biodiversity and human health. Nature 528, 69–76. doi: 10.1038/nature15744

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X., He, Y., Zhang, C., Tian, Y.-A., Lei, X., Li, D., et al. (2021). Physiological and transcriptional responses of phalaris arundinacea under waterlogging conditions. J. Plant Physiol. 261, 153428. doi: 10.1016/j.jplph.2021.153428

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., Liu, M., Chen, Y., Zeng, T., Lu, X., Yang, B., et al. (2021). Plants and microbes mediate the shift in ecosystem multifunctionality from low to high patterns across alpine grasslands on the Tibetan plateau. Front. Plant Sci. 12, 1–11. doi: 10.3389/fpls.2021.760599

CrossRef Full Text | Google Scholar

Wang, S., Wei, Y. (2022). Qinghai-Tibetan plateau greening and human well-being improving: the role of ecological policies. Sustainability 14, 1652. doi: 10.3390/su14031652

CrossRef Full Text | Google Scholar

Wei, Y., Lu, H., Wang, J., Wang, X., Sun, J. (2022). Dual influence of climate change and anthropogenic activities on the spatiotemporal vegetation dynamics over the Qinghai-Tibetan plateau from 1981 to 2015. Earth's Future 10 (5), e2021EF002566. doi: 10.1029/2021EF002566

CrossRef Full Text | Google Scholar

Wu, Q., Zhang, T. (2010). Changes in active layer thickness over the qinghai-Tibetan plateau from 1995 to 2007. J. Geophys. Res.: Atmospheres 115, D09107. doi: 10.1029/2009JD012974

CrossRef Full Text | Google Scholar

Xu, X., Thornton, P. E., Post, W. M. (2013). A global analysis of soil microbial biomass carbon, nitrogen and phosphorus in terrestrial ecosystems. Global Ecol. Biogeography 22, 737–749. doi: 10.1111/geb.12029

CrossRef Full Text | Google Scholar

Yao, T., Bolch, T., Chen, D., Gao, J., Immerzeel, W., Piao, S., et al. (2022). The imbalance of the Asian water tower. Nat. Rev. Earth Environ. 3 (10), 618–632. doi: 10.1038/s43017-022-00299-4

CrossRef Full Text | Google Scholar

Yao, Y., Wang, X., Li, Y., Wang, T., Shen, M., Du, M., et al. (2018). Spatiotemporal pattern of gross primary productivity and its covariation with climate in China over the last thirty years. Global Change Biol. 24, 184–196. doi: 10.1111/gcb.13830

CrossRef Full Text | Google Scholar

Ye, C., Sun, J., Liu, M., Xiong, J., Zong, N., Hu, J., et al. (2020). Concurrent and lagged effects of extreme drought induce net reduction in vegetation carbon uptake on Tibetan plateau. Remote Sens. 12, 2347. doi: 10.3390/rs12152347

CrossRef Full Text | Google Scholar

You, Y., Zhou, N., Wang, Y. (2021). Comparative study of desertification control policies and regulations in representative countries of the belt and road initiative. Global Ecol. Conserv. 27, e01577. doi: 10.1016/j.gecco.2021.e01577

CrossRef Full Text | Google Scholar

Yu, H., Deane, D. C., Zhang, Y., Li, S., Miao, S., Xie, G., et al. (2021). Integrating multiple indices of geobiodiversity reveals a series of regional species-rich areas worthy of conservation in the region of the qinghai-Tibet plateau. Biol. Conserv. 261, 109238. doi: 10.1016/j.biocon.2021.109238

CrossRef Full Text | Google Scholar

Yu, H., Miao, S., Xie, G., Guo, X., Chen, Z., Favre, A. (2020). Contrasting floristic diversity of the hengduan mountains, the Himalayas and the qinghai-Tibet plateau sensu stricto in China. Front. Ecol. Evol. 8, 136. doi: 10.3389/fevo.2020.00136

CrossRef Full Text | Google Scholar

Zhang, T., Barry, R. G., Knowles, K., Heginbottom, J., Brown, J. (1999). Statistics and characteristics of permafrost and ground-ice distribution in the northern hemisphere. Polar Geogr. 23, 132–154. doi: 10.1080/10889379909377670

CrossRef Full Text | Google Scholar

Zhang, R., Tian, D., Chen, H. Y., Seabloom, E. W., Han, G., Wang, S., et al. (2022). Biodiversity alleviates the decrease of grassland multifunctionality under grazing disturbance: a global meta-analysis. Global Ecol. Biogeography 31, 155–167. doi: 10.1111/geb.13408

CrossRef Full Text | Google Scholar

Zhang, H., Ye, C., Wang, S. (2022). Elevation dependence of climate effects on ecosystem multifunctionality states over the qinghai-Tibet plateau. Global Ecol. Conserv. 36, e02066. doi: 10.1016/j.gecco.2022.e02066

CrossRef Full Text | Google Scholar

Zhang, F.-G., Zhang, Q.-G. (2016). Microbial diversity limits soil heterotrophic respiration and mitigates the respiration response to moisture increase. Soil Biol. Biochem. 98, 180–185. doi: 10.1016/j.soilbio.2016.04.017

CrossRef Full Text | Google Scholar

Zhang, G., Zhang, Y., Dong, J., Xiao, X. (2013). Green-up dates in the Tibetan plateau have continuously advanced from 1982 to 2011. Proc. Natl. Acad. Sci. 110, 4309–4314. doi: 10.1073/pnas.1210423110

CrossRef Full Text | Google Scholar

Zhao, N., Li, X. G. (2017). Effects of aspect–vegetation complex on soil nitrogen mineralization and microbial activity on the Tibetan plateau. CATENA 155, 1–9. doi: 10.1016/j.catena.2017.02.025

CrossRef Full Text | Google Scholar

Zhou, J., Kang, S., Schadt, C. W., Garten, C. T. (2008). Spatial scaling of functional gene diversity across various microbial taxa. Proc. Natl. Acad. Sci. 105, 7768–7773. doi: 10.1073/pnas.0709016105

CrossRef Full Text | Google Scholar

Zhou, S., Peng, L. (2021). Applying Bayesian belief networks to assess alpine grassland degradation risks: a case study in northwest sichuan, China. Front. Plant Sci. 12, 773759. doi: 10.3389/fpls.2021.773759

PubMed Abstract | CrossRef Full Text | Google Scholar