The continuous increase of livestock production in Inner Mongolia has caused severe degradation of the grassland ecosystems in recent years. Previous grazing experiments have shown a wide range of vegetation responses between the biome types on a global scale, but there is still a lack of sufficient studies to discern the relative responses of a given biome type. We conducted a meta-analysis of vegetation coverage (VC), plant density (PD), total biomass (TB), above-ground biomass (AGB), under-ground biomass (UGB) and Shannon–Weaver Index (SI) in different grassland types in Inner Mongolia obtained under conditions of different grazing intensities and durations. The results showed that grazing decreased VC, TB, AGB, UGB, and PD significantly. Compared to the global and national average values, the negative effects of grazing to steppe biomass in Inner Mongolia were higher than that on the global scale, while less pronounced than that in China. TB of the meadow steppe in Inner Mongolia increased by 40% under moderate grazing intensity and duration because of compensatory growth. SI of the desert and meadow steppe showed negative linear relationships with the grazing intensity in Inner Mongolia. The percentage changes in AGB, PD, and SI to grazing showed quadratic relationships with the mean annual temperature of the experimental year. With increasing mean annual precipitation, the negative effects of grazing on UGB and SI first decreased and then increased, with that of VC and grazing showing a cubic relationship.
Grassland is one of the most widely distributed vegetation types, covering approximately 40% of the global land surface (Simons and Weisser, 2017) and accounting for 20% of terrestrial production (Scurlock and Hall, 1998). Grassland provides numerous ecological services, including species diversity maintenance and soil erosion control, and acts as a carbon sink (Wang et al., 2016; Zhang et al., 2019). In addition, most of the global demands for milk and meat are provided by grassland ecosystem, which is a major objective of grassland husbandry (Boval and Dixon, 2012). Over-grazing damages the vegetation community characteristics of grassland ecosystems and accelerates grassland degradation (Gillson and Hoffman, 2007; Gan et al., 2012), which also threatens grassland ecological health (Fedrigo et al., 2017).
The effects of grazing on grasslands have been widely investigated (Milchunas et al., 2011), and diverse results have been obtained. Some studies reported that grazing induced increases in biomass, as proposed by the grazing optimization hypothesis (McNaughton, 1976; de Mazancourt et al., 1998; Patton et al., 2007). Other studies suggested that grazing may be adverse to the vegetation community characteristics in most grassland types (Hanke et al., 2014; Jemma et al., 2017). Wang et al. (2002) found that grazing decreased species diversity and increased the spatial heterogeneity of the communities. Dowling et al. (2005) showed that grazing has no apparent benefit to the grassland. The majority of the studies on the effects of grazing on grasslands were performed using controlled grazing experiments (James et al., 2001; McIntyre et al., 2019; Silva et al., 2019) because the grazing experiment was considered as the most effective approach for examining the effects of grazing up to now (Dong et al., 2007). However, various independent experiments on the grazing effects on grasslands have produced divergent results because of (i) different grass types (Yan et al., 2013), (ii) different climate conditions (Ren et al., 2012), and (iii) differences in the duration of grazing treatments (Su et al., 2017). These differences, combined with the complex response of grasslands to disturbances, raise challenges in quantifying the effects of grazing. Therefore, it is necessary to integrate the divergent experiments to provide asynchronous results for a given biome. The meta-analysis method can combine data from various independent studies mathematically, and has emerged as a standard approach for synthesizing research findings (Yang, 2015).
The grassland of Inner Mongolia is one of the major parts of the Eurasian steppes, comprising meadow, typical, and desert steppes distributed along a gradient of decreasing precipitation from northeast to southwest. Meanwhile, it is one of the grasslands with the most frequent human activities in the world since it has served as a key resource for supporting the subsistence of local people for thousands of years. The introduction and expansion of livestock production pose major challenges to the grasslands in this region. Currently, the area of degraded grassland is approximately 50% of the total land area of Inner Mongolia, of which 20% is seriously degraded. Numerous individual grazing experiments have been undertaken to ascertain the effects of grazing on grassland productivity and plant community diversity in Inner Mongolia (Silletti and Knapp, 2002; Focht and Pillar, 2003; Mike et al., 2012). Based on the data from these individual grazing experiments, several meta-analyses have examined the responses of grassland productivity and plant community diversity to grazing at the biome scale. For example, Ren et al. (2012) analyzed the effect of grazing intensity on species diversity of typical steppe in Inner Mongolia, and demonstrated that species diversity had a neutral response to grazing intensity. Su (2017) quantified the AGB response to grazing in Inner Mongolia grassland. Yan et al. (2013) revealed the response of grassland production to grazing in China based on a transect survey and published data. However, these studies only examined the effects of grazing on individual indexes, but neglected the comprehensive responses of the grassland ecosystem. Furthermore, the limited data on grazing experiments have also prevented the examination of the heterogeneity of the effect size. One example is the study of Yan, which used transect survey data to compensate for the insufficient grazing experiment data. Therefore, high uncertainties remain in the effects of grazing on vegetation community characteristics at the biome scale. A comprehensive analysis based on more extensive experimental data is needed to examine the effects of grazing on grassland vegetation community characteristics in a given biome type.
In this study, we compiled data from 433 observations in 61 studies to quantify the effect of grazing on the vegetation coverage (VC), plant density (PD), total biomass (TB), above-ground biomass (AGB), under-ground biomass (UGB), and Shannon-Weaver Index (SI) in different grassland types in Inner Mongolia for conditions of different grazing intensities and durations. Our objectives were to: 1) examine the patterns of the vegetation community characteristics of the grassland ecosystem in response to grazing; and 2) investigate the regulatory effects of environmental factors on productivity and species diversity in response to grazing.
2 Materials and methods
We compiled peer-reviewed studies published before May 2018 using the resources of the Web of Science (WOS), Google Scholar, and China National Knowledge Infrastructure (CNKI). The extensive searches were conducted using combinations of the keywords: grazing intensity, grazing experiment, vegetation coverage, plant density, Shannon-Weaver Index, and Inner Mongolia. To find additional relevant studies, we consulted the reference lists of the papers found in the databases.
Figure 1 shows a flow diagram detailing the procedure for the selection of the studies used in the meta-analysis. We included only those studies meeting the following criteria in the meta-analysis: 1) The studies were carried out in the grassland ecosystem of Inner Mongolia, and data for at least one of the variables VC, PD, AGB, UGB, TB, and SI was reported; 2) The data for replicates, means, standard errors, standard deviations and the sample sizes of both control and grazing groups were available; and 3) Grazing intensity, grassland type, and grazing duration were provided in the study. A list of the data sources used in the study is provided as “Supplementary Data” ①.
Since publication bias is inevitable in any meta-analysis, the standard method used to identify publication bias currently is a funnel plot (Mavridis and Salanti, 2014), which makes a corresponding scatter plot using the estimated effect value of a single study as the X-axis and the sample size of each study as the Y-axis. Because the accuracy of effect estimation increases with an increase in the sample size, the effect values of small sample studies in a funnel plot would be distributed at the bottom of the graph, while the effect values of large sample studies in a funnel plot would be distributed in the middle or top of the graph.
We used GetData Graph Digitizer ( http://getdata-graphdigitizer.com) to extract data from papers only providing standard deviations and sample numbers in graphical figures. Standard deviation (SD) was calculated as follows:
where SE is the standard error; and n is the sample number. Following the exclusion of duplicates in the various studies, 435 sets of vegetation community characteristic data and environmental data from 61 studies were included in our analysis.
Meta-analysis was conducted using METAWIN 2.1, STATA 12.0, and Origin 8.0 software. The natural log of the response ratio (ln R) represented the effect size (Hedges et al., 1999; Wang et al., 2011). The random-effect model was adopted for those cases. The effect-size ln R was calculated as follows (Cohen, 1988; Hayley and Jones, 2010).
with a variance of:
where ln R is the effect size; and are the means of the treatment and control groups, respectively; nt is the sample size of the treatment group; nc is the sample size of the control group; St is the standard deviation of the treatment group; and Sc is the standard deviation of the control group.
We calculated the heterogeneity of effect size by subgroup analysis, and the between-group (QB) heterogeneity was calculated as follows:
with degree of freedom df = m–1, m is the number of groups, ki is the number of comparisons in the ith group, wi,j is the weighting factor of each group, ln Ri+ is the individual weighted log effect size, ln R+ is the weighted log effect size, If QB is larger than a critical value, then the independent variable has a significant influence on the response ratio (Hedges, 1984). Statistical significance was tested at the P < 0.05 level.
Grazing effects were estimated as a PC (percentage change) relative to the control, and calculated using equation (5).
2.3 Classification of grazing indices
The study region extends from the Greater Hinggan Mountains in the east to the Juyan Sea in the west, covering semi-humid, semi-arid and arid regions with different climatic and soil hydro-meteorological conditions. With an average elevation of 1000–1200 m, the terrain shows a gradual downward trend from the southwest to the northeast. The average annual precipitation in the study area is about 100–400 mm, while the average annual temperature ranges from –5.0 to 7.0 °C. By 2010, the grassland area of Inner Mongolia was about 86.67 million ha, accounting for 44.08% of the total grassland area in China, making it an important animal husbandry production base (Dai, 2016). Due to the variation in hydrothermal conditions and topography, the natural grassland vegetation in the study area consists of meadow steppes, typical steppes, and desert steppes.
In grazing experiments, grazing intensity refers to the number of grazing livestock per unit area of grassland, which directly affects the stability and sustainable development of the grassland ecosystem (GB/T 34754–2017). In this paper, we divided the grazing intensity into light grazing (0.15–0.93 sheep per ha per half year), moderate grazing (0.94–1.82 sheep per ha per half year) and heavy grazing (1.83–2.89 sheep per ha per half year). Grazing durations ranged widely from a half month up to several years, and the experiments with more years of grazing are inevitably affected by the cumulative grazing influence of previous years. To examine the effects of grazing on vegetation community characteristics not affected by the cumulative grazing years, we included an assessment of the grazing experiments with 1-year grazing duration or less, which were divided into short-term grazing (0–3 months), moderate-term grazing (3–6 months), and long-term grazing (6–12 months). The grassland types in the Inner Mongolian study sites were classified as meadow, typical, and desert steppes.
3.1 Effects of grazing on grassland vegetation community characteristics
Figure 2 shows the percentage changes of VC, TB, AGB, UGB, PD, and SI in response to grazing in Inner Mongolia. AGB was affected by grazing with a reduction of 35.04% (95% CI: –47% to –27%), while VC was reduced by 27% (95% CI: –39% to –6%), TB was reduced by 45% (95% CI: –56% to –34%), UGB was reduced by 14% (95% CI: –23.5% to –4%), and PD was reduced by 15% (95% CI: –2.5% to –30%). Meanwhile, no significant effect was found on the SI. Between-group heterogeneity tests (Table 1) showed that the effect on VC was not significantly different among grazing intensities or grazing durations, but a significant difference was found among grassland types (QB = 25.3576, P = 0.009). The AGB showed significant differences among grazing intensities (QB = 68.0871, P = 0.001). The UGB exhibited significant differences among grassland types (QB = 18.8179, P = 0.005). The SI demonstrated significant differences among grazing durations (QB = 72.7619, P = 0.001).
Figure 3 shows the percentage changes of the vegetation community characteristics in response to differing grazing intensity, grazing duration, and grassland type. The responses of vegetation communities in different types of grassland were obviously different under various grazing intensities and durations. Differences were evident in the percentage changes in VC, AGB, and UGB among the corresponding groups. The TB of grasslands decreased by 20% (95% CI: –50% to 10%), while SI increased by 27% (95% CI: 5% to 35%) with short grazing duration.
Between-group heterogeneity (QB) and probability (P) for grazing effects on VC, PD, TB, AGB, UGB, and SI among different grazing intensities, grassland types, and grazing durations.
3.2 Effect of grazing on VC
Figure 4 shows the percentage changes of VC in response to different grazing intensities and grazing durations. The effect on VC varied among the desert, typical, and meadow steppes (QB = 25.3576, P = 0.009). Light grazing decreased VC by 40% (95% CI: –59% to –17%) in the desert steppes, while it had no significant effect on other types of grassland; VC decreased by 32% (95% CI: –61% to –11%) in the desert steppes under moderate grazing intensity, while the moderate grazing had no significant effect on the VC of other types of grassland. The VC of grassland in Inner Mongolia was significantly decreased when grazing was heavy, with the greatest reduction occurring in the desert steppe (a reduction of 40%; 95% CI: –58% to –18%), followed by the typical and meadow steppes. VC was more sensitive to grazing intensity in the desert steppe.
VC of the desert steppe was reduced by 50% (95% CI: –111% to 5%) after short-term grazing. However, the 95% CI crossed the invalidation line due to the limited data collected for this group. VC of desert steppe and meadow steppe decreased by 65% (95% CI: –130% to –5%) and 15% (95% CI: –30% to 2%), respectively, after the moderate-term grazing, while VC of typical grassland had no obvious effect from the moderate-term grazing. After long-term grazing, the VC of desert steppe decreased by 23% (95% CI: –48% to 35%), while the meadow steppe showed no significant change. VC of meadow steppe declined after short-term grazing, but that was reversed to an increase after long-term grazing.
3.3 Effect of grazing on TB
Figure 5 shows the percentage changes of TB in response to different grazing intensities and grazing durations. Light grazing reduced TB of the grasslands in Inner Mongolia gradually. TB of the desert and typical steppe decreased by 48% (95% CI: –50% to –49%) and 31% (95% CI: –35% to –32%), respectively, under moderate grazing intensities, while the TB of meadow steppe increased by 40% (95% CI: 25% to 50%). TB of the desert and typical steppe declined with the increase of grazing intensity, but TB of meadow steppe was not affected by the heavy grazing intensity. The biggest drop in TB was always in the desert steppe, followed by the typical and meadow steppe.
TB of desert steppe decreased by 44% (95% CI: –45% to –43%), while it increased by 40% (95% CI: 8% to 71%) after the moderate-term grazing in the meadow steppe. TB of the three types of grasslands decreased significantly after moderate-term grazing, and the ranking of the decreasing amplitude was: meadow steppe > desert steppe > typical steppe. With the increasing of the grazing duration, TB of the grasslands in Inner Mongolia further declined. Therefore, it can be seen that the desert steppe is more sensitive to the grazing duration.
3.4 Effect of grazing on AGB
Figure 6 shows the percentage changes of AGB in response to differing grazing durations and grassland types. AGB decreased with increasing grazing intensity (QB = 68.0871, P = 0.001). Heavy grazing caused a significant negative effect on AGB, with lesser effects being found for moderate and light grazing intensities. However, the decreases in AGB caused by light and moderate grazing intensities were similar to short-term grazing (approximately 20%). AGB for the desert steppe under light and moderate grazing intensities declined by 15% (95% CI: –33% to –4%) and 50% (95% CI: –80% to –24%), respectively.
3.5 Effect of grazing on SI
The SI is a diversity index that indicates the adaptation of species to the environment following competition. Figure 7 shows the percentage changes in the SI caused by differing grazing durations (a) and intensities (b). The SI was more profoundly affected by grazing intensity (QB = 5.9352, P = 0.413) than by the grazing duration (QB = 72.7619, P = 0.001). Compared with typical and meadow steppe, the SI for the desert steppe was more sensitive to grazing intensity, decreasing by 8% (95 CI: –17% to –1%), 11% (95 CI: –20% to –3%), and 19% (95 CI: –24% to –9%) under light, moderate, and heavy grazing, respectively. The SI for the meadow steppe was not affected by either light or heavy grazing, but showed a 10% (95 CI: –21% to –4%) reduction under moderate grazing.
After short-term grazing, the SI of Inner Mongolia did not show a significant change under the light grazing intensity, but it was reduced by 7% (95% CI: –10% to –3%) and 9% (95% CI: –15% to –1%) under the moderate and heavy grazing intensities, respectively. The SI of grassland decreased with the increasing of the grazing intensity after moderate-term grazing, as it declined by 21% (95% CI: –31% to –15%), 30% (95% CI: –47% to –19%), and 42% (95% CI: –73% to –25%) under light, moderate, and heavy grazing intensities, respectively.
4.1 Effects of grazing on grassland productivity
Our study revealed that grazing can reduce the grassland productivity of Inner Mongolia. We compared the effects of grazing on productivity at global scales (Milchunas and Lauenroth, 1993) and national scales (Yan et al., 2013) to Inner Mongolia (Table 2). AGB was reduced by 43% in Inner Mongolia, approximately double the reduction on a global scale. Grazing resulted in UGB reduction in Inner Mongolia, but an incremental increase on the global scale. TB was reduced by 38% in Inner Mongolia, while there is no obvious effect of grazing when examined on the global scale (Table 2). This difference can be attributed to the more palatable species to livestock which occur in Inner Mongolia, as this leads to heavy grazing as a common form of land-use. Furthermore, the rapid expansion of the grazing area largely reduced the soil nutrients and vegetation degradation (Bai et al., 2015; Du et al., 2020).
Percentage changes in grazing effects globally, in China, and in Inner Mongolia
Compared to results at the national scales, the decrement of the grassland biomass was only slight in Inner Mongolia. The effect of grazing on the grassland biomass is mainly regulated by the variations in the community composition in each grassland type. The meadow and typical steppe, which account for 79% of the total steppe area in Inner Mongolia (Dong et al., 2020), involved plentiful biodiversity including species with higher grazing tolerance, such as Stipa baicalensis, Chinensis, Stipagrandis, Cleistogenes polyphylla, Carex piriformis, and Spodinpogon sibiricus, et al. Moreover, owing to the broader study area in Yan et al (2013), the smaller number of experiments in the Inner Mongolian steppe was covered in the meta-analysis, which also largely contributed to the difference.
The grasslands in Inner Mongolia are distributed in arid and semi-arid areas, where the competition for water dominates the inter-species relationships. The precipitation in Inner Mongolia is lower than in other temperate steppes (Mónica et al., 1996; Ni, 2004; Su et al., 2017). The low and unevenly distributed precipitation during the summer can have negative effects on the AGB, even if the precipitation increases in the subsequent period (Fang et al., 2005). The grazing by livestock reduces the individual species and AGB in the low precipitation condition, increasing surface exposure. As a result, soil evaporation is enhanced, and then the plant species are inhibited from growing due to water stress. The AGB of meadow steppe with perennial hybrid grass is positively correlated with the growing season precipitation, while the AGB of the desert steppe is significantly negatively correlated with the amount of precipitation (Sala et al., 1988).
Figure 8 shows the relationships between observed factors and the mean annual temperature of the experimental year. The response ratios of UGB (Fig. 8c) and TB (Fig. 8b) showed linear relationships with the mean annual temperature of the experimental year. The ln R of UGB and TB decreased with the increasing temperature. However, no significant relationship was found between AGB (Fig. 8d), TB (Fig. 8b), and precipitation. Furthermore, grassland productivity is widely limited by various nutrients, such as N, P, K, and mi-cronutrients (Jobbágy et al., 2002; Dangal et al., 2016).
It is worth noting that the TB of the meadow steppe increased by 40% under moderate grazing intensity and moderate-term grazing, because the vegetation biomass exhibited compensatory growth (Burke et al., 1997; Ferraro and Martín, 2002; Bai et al., 2012; Fay et al., 2015). Previous studies have reported that compensatory growth occurs only under a suitable grazing intensity (Xue et al., 2010) or in the croplands where strong cultivation management measures are applied (Münkemüller and Johst, 2006). Furthermore, the seasonal distribution of temperature determines the beginning and length of the vegetation growing season. Increased precipitation and temperature may exert positive effects on the biomass production of meadow steppes (Xiao, 1996).
4.2 Effects of grazing on SI
The SI declined obviously with an increasing grazing intensity and duration, except under the light and short-term grazing. Moderate grazing reduced SI of the meadow and desert steppe. This is different from the intermediate disturbance hypothesis (IDH) (McManus et al., 1972), which assumes that species diversity reaches peak levels at intermediate disturbance (Wang et al., 2003), because grazing jeopardizes the regenerative ability of edible herbage (Milchunas et al., 1988). Moreover, grazing reduces the soil organic matter, which leads to the depletion of soil nutrients, and thus reduces the competitiveness of plants (Mazumder, 1998).
The water-related factors are the main driving variables of species diversity in the thermal and temperate zones (Connell, 1979; Watt and Gibson, 1988; Wang et al., 1998). With increasing precipitation, the negative effects of grazing on SI first decreased and then increased (Fig. 9f). Previous studies verified IDH through a large number of experiments and found that the IDH appeared in different grazing stages of different grassland types. For example, it was generally confirmed at the stage of light to moderate grazing intensity in humid meadows, and the stage of moderate to heavy grazing intensity in the dry grassland (Nautiyal et al., 2004; Kreft and Jetz, 2007).
4.3 The relationship between productivity and species diversity
The relationship between productivity and species diversity is crucial for explaining the attributes of vegetation community characteristics (Field et al., 2005), and it is generally believed that higher biodiversity can maintain higher productivity. Our results confirmed that the grassland productivity decreased with the reduction of species diversity because the reduction in the number of species decreased the supply sources of productivity (Fig. 10) (Oesterheld et al., 1999; Bai et al., 2001; Guang et al., 2002). The analysis of the structural equation model showed that the grassland vegetation community characteristics were affected not only by grazing intensity and duration, but also the observed community characteristics were closely related to each other. The SI was negatively linearly correlated with VC and PD. The heavy grazing intensity and duration had direct negative effects on VC, PD, AGB, and SI, while they had indirect negative effects on TB (Fig. 11).
Several studies have analyzed the relationships between productivity and species diversity from different spatial scales and ecosystem types since 1990, and six relationship types, including positive correlation (Zhu et al., 2020), humped-back correlation (Shea et al., 2004), negative correlation (Wang et al., 2001), unitary correlation (Guo and Berry, 1998), no correlation (Al-mufti et al., 1977) and “U” relationships (Guo and Berry, 1998), were found between them (Scheiner and Willig, 2005; Jain et al., 2007). The relationship between species diversity and productivity is known to be driven by multiple factors, such as the climate (Kondoh, 2001), soil (Gillman and Wright, 2006), and density (Waide et al., 1999), and each factor changes with the spatial scale (Mittelbach et al., 2001). On the other hand, due to differences in environmental conditions, species composition, and habitats, differences in grassland productivity are often out of synchronization with changes in the factors that affect species diversity. Therefore, productivity may be a poor predictor of plant diversity (Peter et al., 2011). It is difficult to explain plant function without taking into account of all these elements (Hadgu et al., 2009). Furthermore, some details such as sampling area, location of sampling sites, sampling number, and methods of measuring productivity would also affect the relationship between productivity and species diversity.
4.4 The need for appropriate grazing intensity criteria
The standard for the divisions of the grazing intensity affects the accuracy of the meta-analysis by integrating the groups. In addition, an inappropriate division of grazing intensity may lead to grassland degradation under the guidance of a correct grazing policy (Wolfe et al., 2004). Reconsidering the grazing intensity criteria is necessary for defining the relationship between grazing and the grassland degradation degree. Grazing modifies the vegetation community characteristics by changing species composition (Šimova et al., 2013) and alters grassland productivity by changing VC, vegetation height, and biomass (Adler et al., 2011). Proposed indicators of grassland deterioration include vegetation composition, VC, litter, biomass, SI, and soil nutrients (Guang et al., 2002). Some studies inferred the grazing intensity from the distance to drinking water (Wang et al., 2020) or herder huts (Marcos and Willig, 2004). Other studies classified the grazing intensity according to the feces density (Rodríguez et al., 2003), species composition (Xia et al., 2005), or the proportion of herbivore mammals (Peter, 2012). A uniform standard for the grazing intensity classification has not been established yet. Currently, quantifying the number of herbivore mammals per unit area and unit time is the generally used criterion for classifying grazing intensity, but it cannot reflect the corresponding effect. The moderate grazing intensity in one ecosystem may be defined as heavy grazing intensity in another ecosystem. Therefore, the use of grazing intensity should be redefined based on different grassland types according to plant composition, herbage yield, soil coarse fraction, amount of soil litter, and soil erosion rate (Haynes et al., 2013).
4.5 The influences of climate change on the grassland vegetation community characteristics under different grazing intensities and durations
According to the Fifth IPCC Climate Assessment Report, the rate of the national average annual temperature rise is about 0.21–0.25 °C (10 yr)–1 in the recent 50–60 years, which is higher than the global average over the past 100 years (1909–2011) (IPCC, 2013). Sarah et al. (2015) suggested that vegetation community characteristics and ecosystem processes are strongly affected by the climate change. The increasing temperature will bring forward the spring phenological period and prolong the growing period of grasslands, resulting in a significant increase in AGB, thus improving the grassland productivity (Silva, 2020). However, more frequent meteorological disasters caused by extreme climate, such as extremely high temperatures and droughts in summer, could decrease the grassland biomass and species diversity. The grassland will deteriorate severely assuming a long-term exposure of grasslands to extreme climate (Ali, 2020). The trend of climate warming and drying across Inner Mongolia, as well as the occupation of surface water by grazing, have reduced surface water and hinder the growth of vegetation (Zhao, 2019).
To determine the effects of grazing on vegetation community characteristics not affected by the cumulative grazing years, we analyzed the effects just within grazing durations of one year or less. But the lack of experimental data for more than one year's grazing duration in our study results in certain limitations over the long term. In addition, since many unpublished and grey literature studies with statistically non-significant results are inaccessible and could not be compiled into the database, some limitations of this study remain to some extent. We evaluated the publication bias by funnel plots (Gonnet et al., 2003) (Fig. 12). Most of the samples in this study were distributed at the top of the funnel plot and clustered toward the center, indicating there were a few issues with publication biases in this study. Besides, the different variables that represent the vegetation growth in various studies, such as vegetation height, the Simpson Index, annual net primary production, and other ecosystem trait indexes (Li et al., 2008), need to be considered fully in further studies. Therefore, the meta-analysis related to the effects of grazing on grassland requires further study and exploration.
The Inner Mongolia grassland ecosystem exhibits a pronounced response to grazing. Unlike the intermediate disturbance hypothesis (IDH), SI of the desert and meadow steppe showed a negative linear relationship with the grazing intensity in Inner Mongolia. The meadow steppe in Inner Mongolia had a typical compensatory growth phenomenon with suitable grazing intensity and duration. The negative effects of grazing on steppe biomass in Inner Mongolia were higher than that on the global scale, while less pronounced than that in China.
 ① Please log in the journal's website: www.jorae.cn, find this article and download the supplementary data attachment.