Heavy Grazing Altered the Biodiversity–Productivity Relationship of Alpine Grasslands in Lhasa River Valley, Tibet

Grazing is a crucial anthropogenic disturbance on grasslands. However, it is unknown how livestock grazing affects the relationship between biodiversity and productivity of alpine grasslands in Tibet. We carried out a grazing-manipulated experiment from 2016 to 2019 with grazing intensity levels of null (control, grazing exclusion, C.K.), moderate grazing [1.65 standardized sheep unit (SSU) per hectare, M.G.], and heavy grazing (2.47 SSU per hectare, H.G.) on a typical alpine grassland in the Lhasa River Basin, central Tibet. We measured aboveground biomass (AGB), species assembly (alpha and beta diversity indices), and soil nutrients’ availability. The results showed that grazing differently affected plant community in different treatments. Notably, the total dissimilarity value between C.K. and H.G. is 0.334. Grazing decreased the Shannon–Wiener index, increased the Berger–Parker index from 2016 to 2018 significantly, and decreased AGB and total soil nitrogen (STN) significantly. Our results also showed that the grazing affected the relationship between AGB and diversity indices and soil nutrients, including soil organic carbon (SOC) and total soil phosphorus (STP). Specifically, AGB decreased with increasing SOC and STP in all treatments, and heavy grazing changed the positive relationships between AGB, STP, and Shannon–Wiener index to negative correlations significantly compared with grazing exclusion. There was a significant negative correlation between Berger–Parker and Shannon–Wiener indices under each treatment. The general linear models showed that H.G. altered the relationship between diversity and productivity of grassland in central Tibet, and AGB and Shannon–Wiener index positively correlated in C.K. but negatively correlated in H.G. Our study suggests that H.G. caused a negative relationship between plant diversity and productivity. Therefore, sustainable grazing management calls for a need of better understanding the relationship between biodiversity and productivity of alpine grassland in central Tibet.


INTRODUCTION
Livestock grazing is a crucial anthropogenic disturbance on grasslands , reshaping plant community assembly and affecting ecosystem functions (Cingolani et al., 2005). Grazing reduces plant species diversity in dry forests (Schulz et al., 2019) but increases it in grazed areas (Gonzalez-Hernandez et al., 2020). According to the intermediate disturbance hypothesis (Mcnaughton, 1983), the plant community productivity and diversity were the highest at light or moderate grazing intensity. However, over-grazing may reduce species richness (Cardinale et al., 2012). Moreover, grazing can affect plant community structure by changing species composition (Lin et al., 2010), where graminoid and herb species can be inhibited (Rada et al., 2014). Grazing has little influence on nestedness but can increase species turnover (Adler et al., 2001). Furthermore, heavy grazing can result in high species replacement (Sanaei et al., 2018). However, how different grazing intensities affect plant community assembly is still unclear.
Aboveground biomass (AGB) is an important indicator of the quantitative characteristics of ecosystem functions to an anthropogenic disturbance (Ma et al., 2010). Many studies have demonstrated that grazing significantly affected AGB and primary productivity at different scales (Wu, 2013;Su et al., 2017;Liu et al., 2020). Notably, heavy grazing seriously declined grassland productivity and plant coverage in grasslands and intensified land desertification . However, the effect of grazing intensities on alpine grasslands is still under debates. Therefore, a better understanding of the impact of grazing disturbance on AGB could help us to enhance ecosystem stability when facing external disturbances (Rusina et al., 2013). It should examine the balance between ecosystem function and ecosystem service for rational rangeland management (Dias et al., 2017).
The biodiversity-productivity relationship of grassland is a controversial topic in the field of ecosystem function and stability research (Axmanova et al., 2013;Liang et al., 2016;Wu et al., 2019a). The relationship is complex and is affected by many factors that are monotonous linear (positively or negatively) or have no relation (Almufti et al., 1977;Mcnaughton, 1983;Gough et al., 1994;Gillman and Wright, 2006;Hector and Bagchi, 2007;Cardinale et al., 2012;Zhang R. et al., 2018). A positive linear relationship between biodiversity and productivity has been found in the Eurasian Steppe (Bai et al., 2007) and a weak positive relation in the Qinghai-Tibetan Plateau (Fayiah et al., 2019). Zhang et al. (2011) found that the relationship can be changed with special research scales in the temperate steppe of Inner Mongolia. Wang et al. (2013) demonstrated no significant correlation between plant species richness and community productivity in the Tibetan semi-arid regions where available water is limited. How plant diversity regulates vegetation productivity is receiving wide attention to Erfanzadeh et al. (2015) and Grace et al. (2016). However, many studies rarely examine the impact of environmental factors and ways of utilization of the relationship shapes (White et al., 2019). Therefore, it is necessary to examine the mechanisms behind the different biodiversity-productivity relationships (Gross and Cardinale, 2007).
The Qinghai-Tibetan Plateau (QTP) is characterized by high altitude and high fluctuation in temperature, intense radiation, and fragile ecosystems (Gao et al., 2010;Liu et al., 2020). It covers more than 20% of China's grassland (Cao et al., 2004). Alpine grasslands were an essential resource of pastoral production activities (Yi et al., 2019) affected by many biotic and abiotic factors (He et al., 2005;Zhu et al., 2016). At large scales, alpine grasslands on this plateau are thought to be sensitive to climate change (Wu et al., 2019b). At a regional or local scale, they are disturbed by livestock grazing activities (Bai et al., 2007;Milchunas et al., 2012;Erfanzadeh et al., 2015). The relationship between plant diversity and productivity is influenced by environmental variables differently at various spatial scales in Tibetan grasslands (Ma et al., 2010). Xiong et al. (2014) found that grazing exclusion improved forage yield and has negative influences on species diversity in alpine meadows. Wu et al. (2014b) discovered that AGB increased with precipitation in growing season and decreased with growing season temperature in alpine pastures across the northern Tibetan Plateau. They even found a unimodal pattern between AGB and species richness within steppes but no significant relationship neither in meadows nor desert steppes (Wu et al., 2014a). Previous studies focused on the impact of grazing or fencing on the productivity-species diversity relationship in the northern part of the Qinghai-Tibetan Plateau. Besides, rare efforts were done to explore how grazing intensity affects the alpine grasslands in central Tibet. Therefore, it is necessary to consider how different grazing intensities shape the relationship between plant diversity and productivity of grasslands in central Tibet.
Specifically, in the current study, we aim to (1) explore how the diversity-productivity relationship is driven by grazing management, (2) examine whether the relationship between plant diversity and productivity is monotonous linear (positively or negatively) or not related, and finally, and (3) inspect whether the relationship has been altered by heavy grazing. We examined the differences in plant community diversity, AGB, and soil nutrient content between three grazing intensity levels. Finally, we analyzed the relationship between plant diversity and productivity under the different grazing intensities.

Study Area and Experimental Design
The research station locates in the Linzhou County (29 • 45 ∼30 • 08 N, 90 • 51 ∼91 • 28 E, 3,900 m.a.s.l.), approximately 60 km away from Lhasa, the capital of the Tibet Autonomous Region, China (Figure 1). Table 1 shows community productivity, plant species richness, and soil nutrients. At this station, annual precipitation ranges from 343 to 466 mm, over 70% of which falls in the summer monsoon season, hereafter named the plant-growing season. Mean air temperature during the plant-growing season ranges from 10.6 to 12.6 • C (Chen et al., 2014). The vegetation is alpine steppe dominated by Pennisetum centrasiaticum, Tripogon bromoides, Kobresia  All data are measured in the field and lab, and the values are average of multiple years, which are expressed as mean ± standard error. C.K., grazing exclusion; M.G., moderate grazing; H.G., heavy grazing. Soil nutrients include the contents of soil organic carbon (SOC), total soil nitrogen (STN), and total soil phosphorus (STP).
pygmaea, and Carex atrofusca. However, degradation indicator species are Stipa capillacea, Astragalus strictus, Euphorbia kansui, and Ligularia rumicifolia. The soil is poor with an average depth of about 20-30 cm. The grassland at this station was fenced and excluded from livestock grazing since May 2011. In 2016, we built 15 plots of 10 m × 10 m within the fenced area, grouped randomly into three different treatments. According to the grassland carrying capacity in the Lhasa River Basin , the grazing intensity level was manipulated at null (control, grazing exclusion, C.K.), moderate grazing [1.65 standardized sheep unit (SSU) per hectare, M.G.], and heavy grazing (2.47 SSU per hectare, H.G.). Nine fixed quadrats of 0.5 m × 0.5 m, at least 0.5 m away from the fence lines, were located in each plot. The grazing was carried out once per month, between July 15th and September 17th, from 2016 to 2019. Grazing duration is 8 h/day each time, with two female healthy sheep (native sheep) for M.G. and four female sheep for H.G., respectively.

Data Collection and Calculation
Plant composition was surveyed before the grazing experiment each year. Plant species height was measured and averaged as the community height. Species coverage was visually estimated and summed for the community coverage. The numbers of species within each quadrat and plot was used as plant species richness (n) at the quadrat and plot levels, respectively. A quadrate near each fixed quadrat was used to evaluate AGB. All plants in each quadrat were harvested, kept in an envelope, and weighed after being oven-dried at 65 • C for 48 h. After harvesting AGB, three soil samples were collected from each plot. The soil samples were air-dried in a dark room. Finally, the contents of soil organic carbon (SOC), total soil nitrogen (STN), and total soil phosphorus (STP) were measured for further analysis. A given species' frequency (F ri ) was relative to the sum of all species that appeared in nine quadrats at each plot. Similarly, we calculated relative species coverage (C ri ) and height (H ri ) (Wu et al., 2019b). Thus, the Shannon-Wiener diversity index (H ) was calculated based on the relative species dominance (P i ), which is the mean of relative species coverage (C ri ), frequency (F ri ), and height (H ri ) (see Eqs. 1-3). We also calculated the Berger-Parker index (B ), which is the ratio of the max number of dominant species (n max ) to all species numbers (N) (Morris et al., 2014), Eq. 4.
The beta (β) diversity refers to the dissimilarity of species composition among different communities or species replacement rate along the environmental gradients. In this study, the β diversity index was calculated at plot level among treatments according to the approach of Baselga (2010). The binary presence or absence matrices for plant species across 15 plots were used for a measured dissimilarity index. Species turnover was expressed as Simpson dissimilarity (β sim , Eq. 5). Species nestedness was expressed as a nestedness-resultant fraction of Sørensen dissimilarity (β sne , Eq. 6). Species overall beta diversity was expressed as Sørensen dissimilarity (β sor , Eqs. 7-8), where a is the number of species shared in two plots, b is the number of species that only occurred at the species-poorest plot, and c is the number of species only appearing on the species-richest plot (Leprieur et al., 2012;Tonn et al., 2019). (1)

Statistical Analyses
The data analysis was conducted with R statistic software version 4.0.3 for Windows. First, beta diversity indices (β sim , β sne , and β sor ) were calculated by using the "betapart" package. The differences in plant diversity (Shannon-Wiener index and Berger-Parker index), soil nutrients, and AGB between treatments were examined with the Kruskal-Wallis and tested by rank, which is a non-parametric alternative to one-way analysis of variance (ANOVA) when the assumptions of homogeneity of variance and normality are not met (Wu et al., 2019b). The correlation coefficient matrix was calculated and tested the significance by "corrplot" package within each of the three treatments. Finally, we established general linear models (GLM) to analyze the relationships between soil nutrients, plant diversity indices, and community AGB at different grazing intensities.

Comparisons of Community Assembly
Community dissimilarities was dominated by species turnover, with a relatively small contribution of the nestedness component (Figure 2). In 2016, the dissimilarity value of species turnover of C.K. to H.G. is 0.164, which is lower than C.K. to M.G.

Comparisons of AGB and Soil Nutrients
With the increase of grazing intensity, AGB and STN decreased significantly (Figures 4A,C, P < 0.05), whereas STP under H.G. were lower than C.K. significantly ( Figure 4B, P < 0.05), and SOC under M.G. were higher than C.K. and H.G. significantly ( Figure 4D, P < 0.05).
Aboveground biomass and soil nitrogen decreased dramatically with increasing grazing intensity (Figures 4A,C, P < 0.05). STP under H.G. was lower than that under C.K. and M.G. (Figure 4B, P < 0.05). However, there was no significant difference in STP between M.G. and C.K. (Figure 4B, P > 0.05). SOC under M.G. was higher than that under C.K. and H.G. (Figure 4D, P < 0.05). There was no significant difference in SOC between H.G. and C.K. (Figure 4B, P > 0.05).

The Relationship Among Plant Diversity, AGB, and Soil Nutrients
There is a significantly negative correlation between Berger-Parker and Shannon-Wiener indices in all treatments (Figure 5). The relationships between AGB, STP, and Shannon-Wiener index changed from positive to negative correlations significantly under C.K. and H.G. There is no significant relationship between STN, plant diversity, and aboveground biomass among the three treatments.
Bivariate regressions showed similar patterns of AGB under the three treatments, which decreased with the increasing SOC and STP (Figures 6A,C). Grazing altered the slopes of AGB varying along with SOC and STP compared to C.K., namely, the slopes of AGB along with SOC and STP were steeper with increasing grazing intensities. The proportions of the total variance AGB explained by SOC and STP were 36 and 39% under C.K., 16 and 35% under M.G., 17 and 18% under H.G., respectively.

DISCUSSION
Grasslands in central Tibet are intensively affected by anthropogenic disturbances (grazing, mowing, and construction) and climate change (such as precipitation, temperature, wind  soil, and nutrient) (Zhu et al., 2016). Previous studies were concentrated on the effects of exclusion grazing on grasslands in north Tibet (Spooner et al., 2002;Xiong et al., 2014;Tai et al., 2021). However, it is rarely reported how fencing vs. grazing influences the grassland ecosystem in central Tibet. In the current study, we first demonstrated that grazing alerted alpha and beta diversity compared to grazing exclusion. Next, we detected that grazing activity decreased community AGB, and heavy grazing decreased soil total nitrogen and soil organic carbon. Then, we found that heavy grazing changed the positive relationships between AGB, STP, and Shannon-Wiener index to negative ones compared with grazing exclusion. A negative correlation was also found between Berger-Parker and Shannon-Wiener indices under each treatment. Finally, the GLM showed that heavy grazing altered diversity and productivity in central Tibet grassland. All potential mechanisms and causes are discussed below.

Effects of Grazing Intensity on Plant Community Assembly
Community dissimilarities were dependent on species composition and turnover (Figure 2A), strongly correlated with precipitation under C.K. Herbaceous communities exhibited high rates of temporal species turnover (Bar-Massada and Hadar, 2017), which can be explained by the modern co-existence theory (Chesson, 2000). On a given land, species turnover may be correlated with the different land management. Our results are similar to other scientists (Sanaei et al., 2018) who found that grazing intensity affected vegetation composition by enhancing species replacement. Peper et al. (2011) suggested that grazing should be considered in species composition shifts along ecological gradients. In this study, livestock frequently foraged for family Gramineae and Cyperaceae plants. Sheep grazing behavior can limit the growth of dominant species and increase habitat heterogeneity and species niches (Barot and Gignoux, 2004;Olden and Halme, 2016;Fischer et al., 2019). However, some other species cannot exist and may be lost due to sheep trampling (Su et al., 2005).
Grazing significantly influenced the diversity indices of alpine steppes in central Tibet. Specifically, grazing decreased the Shannon-Wiener index and increased the Berger-Parker index significantly compared to C.K. in 2016 to 2018 (Figure 3). The diversity and dominance of the community have changed greatly under H.G. It indicated that grazing altered the stability of the plant community. These findings were in line with findings that the Shannon-Wiener index under ungrazed sites was greater, but the Simpson dominance index was smaller under over-grazing sites in alpine Kobresia meadows in Tibet (Wei et al., 2008). However, with increasing grazing intensity, vegetation coverage decreased dramatically but had no effects on species richness, Shannon-Wiener index, and Pielou evenness index in Gangcha County, Qinghai Province (Chu et al., 2019). In the southern rangelands of Kenya, heavy grazing decreased species diversity, richness, and evenness (Jawuoro et al., 2017). Dingaan et al. (2016) found no dramatic difference in species richness between L.G. and H.G. sites in Kimberley. Therefore, due to different study areas and research objects, grazing effects on diversity are significantly different.

Effects of Grazing Intensity on Soil Nutrients and AGB
Grazing leads to a reduction in productivity. Our results showed that AGB decreased significantly with the increase of grazing intensity ( Figure 4A). The livestock can directly reduce the AGB by eating a large amount of edible herbage (such as Gramineae, etc.). Removing the plant AGB will directly limit the plant growth and reproduction, and this process will redistribute productivities of photosynthesis (Hardman-Mountford et al., 2013). However, these results were inconsistent with other scientists. Olofsson et al. (2001) thought heavy grazing enhanced primary productivity of summer ranges, while moderate grazing decreased primary productivity in northern Norway. These results conceived that grazing activities are the main factors causing changes in grassland productivity, but others believed that grazing disturbance had a greater indirect effect on AGB via plant diversity (Zhang R. et al., 2018). Therefore, it is necessary to detect whether grazing directly affects plant productivity in a future study.
Grazing may affect soil nutrient parameters (i.e., organic matters, soil total and available nitrogen, and soil total and available phosphorus). Our results showed that grazing had dramatically adverse effects on soil nutrients. Specifically, heavy grazing significantly influenced STN and STP comparing to grazing exclusion, but moderate grazing significantly increased both nutrients. Our results agree with Frank et al. (1995), Steffens et al. (2008), Talore et al. (2016), and Dai et al. (2021). Moreover, Wu et al. (2010) had reported that grazing could influence soil carbon pool and nitrogen sequestration and storage in the Qinghai-Tibetan Plateau. Dong et al. (2012) also found that high grazing intensity decreased soil organic carbon and total nitrogen than low grazing intensity did. A meta-analysis concluded that grazing significantly decreased soil organic carbon, total soil nitrogen, and total soil phosphorus .
Furthermore, overgrazing significantly and directly increased the exposure of soil surface and water erosion storage in grassland ecosystems (Wang et al., 2011). This may be attributed to trampling by heavy grazing that is much greater than other treatments on alpine grassland. Trampling can increase the exposure of topsoil and wind erosion and reduce SOM (Wang et al., 2011). Additionally, topsoil exposure made soil temperatures down, which can accelerate the rate of soil organic matter decay and mineralization, thus declining soil organic matter levels. Soil parameter changes are closely associated with vegetation deterioration or over-grazing (Dong et al., 2012). Soil is the physical support of vegetation growth and the main source of plant nutrients, so they often independently interplay with each other (Li et al., 2014). Su et al. (2005) found that continuous grazing dramatically decreased vegetation coverage, plant litters, soil organic carbon, and soil nitrogen contents. Thus, the vegetation and soil influences restrict each other and coordinate development as a unified interaction network. Therefore, grazing by livestock directly affects vegetation and soil characteristics (Gao et al., 2009), but the soil characteristic change will react on the vegetation (Olofsson et al., 2001), thus forming a complicated vegetation-soil network.

Effects of Grazing Intensity on the Biodiversity-Productivity Relationship
In grazing ecosystems, soil management profoundly impacts ecosystem stability (Tscharntke et al., 2005). Soil properties are important for formations of vegetation productivity, species richness, and community assembly (Dingaan et al., 2016). Particularly, soil organic carbon, nitrogen, and phosphorus are highly essential for plant growth and productivity, which are affected by physical, chemical, and biological factors in multiple ways (Yuan and Jiang, 2021). In our study, the SOC and STP significantly negatively correlate with AGB in all treatments (Figure 5). It maybe that forage growth and compensatory growth absorb and utilize more soil nutrients (i.e., SOC and STP), which caused more AGB and less soil nutrient content (Yuan et al., 2015), as Fay et al. (2015) also proved that plant diversity and productivity are closely related to soil nutrients. Another reason is that livestock ingested plant aboveground biomass (such as plant litter and green leaves), causing decreasing SOC content (Yan et al., 2020) and increasing AGB . Furthermore, our GLM results also showed that AGB decreased significantly with the increase of the values of SOC and STP. More interesting, grazing dramatically altered the slopes of AGB along with SOC and STP compared with grazing exclusion (Figures 6A,C). The possible reasons are that the C.K. had higher vegetation coverage (62.67 ± 1.17%, Table 1), more litter, and moister topsoil; hence, litter decomposition rate and content of SOC were higher in C.K. than M.G. and H.G. (Dimitrakopoulos, 2010). Meanwhile, litter may hinder plant growth and decrease aboveground biomass (He et al., 2012), as intensive grazing has affected and changed soil nutrients and further ignited to changing of ecosystem functions . Previous studies manifested that grazing caused decreasing of AGB and plant coverage, increasing of soil water evaporation and erosion, and decreasing SOC and STN consequently (Rickart et al., 2013;Lu et al., 2017;Su et al., 2017;. Therefore, plant AGB and soil nutrient availability were directly or indirectly affected by livestock. Grazing exhibited no significant influence on the Shannon-Wiener index and SOC and STN (Figure 5). The finding seems to be similar to Wang et al. (2017), who demonstrated that grazing intensity had no effects on plant composition and soil properties across Tibet. Moreover, the pattern of the Shannon-Wiener index declined with SOC dramatically in grazing exclusion, while the pattern of the Shannon-Wiener index increased with SOC without significance in heavy grazing ( Figure 6D). However, the pattern of the Berger-Parker index increased with SOC significantly in grazing exclusion (Figure 6G), which was opposite with the pattern of the Shannon-Wiener index with SOC in grazing exclusion. The relation varies in the study may be due to the highly significant negative correlation of the Shannon-Wiener index with the Berger-Parker index at all treatments (Figures 5, 6L). Our results correspond to previous studies (Prober and Wiehl, 2012;Wang et al., 2017). One possible reason is that grazing had changed the species competitively. Specifically, heavy grazing reduced the competitive exclusion among plants. Another possible reason is that heavy grazing had relatively higher urine and dung contents and more trampling (Mipam et al., 2021).
Heavy grazing altered the relationship between diversity and productivity. Specifically, the Shannon-Wiener index-AGB relationships tended to be positive on grazing exclusion, but a negative relationship was detected on heavy grazing ( Figure 6J). These outcomes seem to be similar to other studies, namely, the patterns of diversity and biomass display a positive or negative linear trend across different grassland ecosystems of the Qinghai-Tibetan Plateau (Zhang et al., 2011;Milchunas et al., 2012). These variations in the study may be due to unavailable water, temperature, precipitation, and soil nutrient limitation on a different spatial scale for natural grassland (Zobel and Paertel, 2008;Li et al., 2019) and maybe driven by sample effects (Wolda, 1981), density-dependent factors, plant biomass, livestock grazing intensity, or other potential factors on a local scale (He et al., 2005;Ma et al., 2010). Grazing can significantly influence plant diversity and productivity patterns of plant communities (Guo and Berry, 1998). The major differences between the positive linear relationship at grazing exclusion and negative linear relationship at heavy grazing we obtained in this study could be connected with livestock grazing activities (Rillig et al., 2019). The results are similar to other researchers, who found that the patterns of biodiversity-productivity were overall linear models in alpine grassland but changed with temporal-spatial scales (Du et al., 2003). Li et al. (2019) found that a plant community with relatively few individuals of each species and concurrently with more species could be accompanied by bigger plant biomass, so the species richness and Shannon-Wiener index have a positive linear pattern with both above-and belowground biomass. Comparing to our study, with lack of grazing livestock activities, each species in a plant community could fully use natural resources (e.g., light, available water, and soil nutrients) for growth (Zuo et al., 2012). However, when grazing at grassland, the livestock could alter the stability of plant community and plant species' competitiveness or even lead to plant species' niche changing (Pavlu et al., 2021). One possible reason is that the response rate of AGB and species diversity on grazing is different in plant community. In other words, heavy grazing decreased AGB remarkably ( Figure 4A) but gradually changed plant diversity (Figures 2, 3). It indicated that the community stability of diversity and productivity exhibits hysteresis concerning grazing (Cipriotti et al., 2019). Another possible reason is that grazing altered the plant species' competitive exclusion within communities (Huston, 1979;DeMalach et al., 2016) and plant species' synchrony. Plant competition and synchrony may be different in communities (Sanaei et al., 2018). Grazing activities released the plant competition for light and niche breadth (Borer et al., 2014). In our study, heavy grazing decreased AGB significantly (Figure 4A), resulting in decreasing competition among plant species and increasing remnant niche breadth of other species; thus, the poorly competitive species may colonize into the communities. Therefore, heavy grazing can influence plant diversity-productivity relationship at small scales through plant competition and colonization (Fayiah et al., 2019). So, it is critical to reveal the effects of species losses on ecosystem functioning (Mulder et al., 2004).

CONCLUSION
This study suggested that heavy grazing lead to high species replacement and greatly changed the species diversity in plant community. Intensified grazing decreased plant productivity and SOC. Furthermore, a positive relationship between the Shannon-Wiener diversity and AGB was found under the condition of grazing exclusion, but a negative relationship was detected on heavy grazing. In conclusion, heavy grazing altered plant community features, soil nutrients, and the patterns of biodiversity-productivity relationship. However, this study focused on alpine grassland in central Tibet. Therefore, it is critically necessary to set a reasonable grazing intensity based on alpine grassland carrying capacity at different areas in Tibet.

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

AUTHOR CONTRIBUTIONS
MX: concepts, ideas, experimental design, data collection, and analysis. JXW: experimental design, interpretation of the results, and editing. JJW: data collection and analysis. YG: experimental design and data collection. DL: ideas, experimental design, interpretation of the results, and editing. YP: editing and part of data collection. XZ: interpretation of the results. All authors have reviewed and approved the final manuscript.