key: cord-0786651-es0x1dig authors: Shi, Yan; Gao, Jay; Li, Xilai; Li, Jiexia; Brierley, Gary title: Effects of disturbances on aboveground biomass of alpine meadow in the Yellow River Source Zone, Western China date: 2022-03-08 journal: Ecol Evol DOI: 10.1002/ece3.8640 sha: c0c592c936a670c291cd97a5d8108539642f8f05 doc_id: 786651 cord_uid: es0x1dig A field experiment quantifies the impacts of two external disturbances (mowing‐simulated grazing and number of pika) on aboveground biomass (AGB) in the Yellow River Source Zone from 2018 to 2020. AGB was estimated from drone images for 27 plots subject to three levels of each disturbance (none, moderate, and severe). The three mowing severities bear a close relationship with AGB and its annual change. The effects of pika disturbance on AGB change were overwhelmed by the significantly different AGB at different mowing severities (−.471 < r < −.368), but can still be identified by inspecting each mowing intensity (−.884 < r < −.626). The impact of severe mowing on AGB loss was more profound than that of severe pika disturbance in heavily disturbed plots, and the joint effects of both severe disturbances had the most impacts on AGB loss. However, pika disturbance made little difference to AGB change in the moderate and non‐mowed plots. Mowing intensity weakens the relationship between pika population and AGB change, but pika disturbance hardly affects the relationship between mowing severity and AGB change. The effects of both disturbances on AGB were further complexified by the change in monthly mean temperature. Results indicate that reducing mowing intensity is more effective than controlling pika population in efforts to achieve sustainable grazing of heavily disturbed grassland. decreased grassland aboveground biomass (AGB) (Dong & Sherman, 2015) , which not only governs the quantity of forage production but also reflects livestock carrying capacity (Tao et al., 2005; Zhu et al., 2005) . These factors occur alongside climate changes, collectively influencing the condition of grazing adapted ecosystems and impacting upon the sustainability of prevailing land-use practices in this region (Harris, 2010) . Grasslands on the QTP have been grazed for millennia (Li, 1997) . Although grazing can improve grassland health by promoting the diversity of plant species (Wang et al., 2005; Wilcox et al., 2017) and stimulating grass growth by boosting the root system (Yan et al., 2020) , it adversely influences grassland vegetation after plant biomass is removed and consumed by livestock (Biondini et al., 1998; Hik & Jefferies, 1990) . In recent decades, grasslands in the Yellow River source zone of the QTP have been overgrazed, at an estimated ratio of 67.88% in 2010, or an equivalent of 6.52 million sheep units overloaded (Zhang et al., 2014) . Overgrazing is regarded as the main cause of grassland degradation, reducing plant biodiversity and productivity (Dong et al., 2020; Török et al., 2018) . Conversely, grazing exclusion is effective at restoring degraded alpine grassland (Liu et al., 2020; Wu et al., 2019) . Determining and implementing an appropriate grazing intensity is a key management strategy to support the sustainable use of grasslands (Deak et al., 2017; Li et al., 2017) . Grazing intensity refers to the cumulative effects of livestock on grassland, mainly represented by forage intake, over a given time period in a given area (Holechek et al., 1998) . It directly affects vegetation structure, coverage, and biomass volume (Toth et al., 2018) . Although plateau pikas (Ochotona curzoniae) may not devour as much grassland biomass as livestock, they do more damage to grassland roots and play a more important role in enlarging bare patches within the grassland (Chen et al., 2017) . If controlled under an appropriate population, pikas are ecosystem engineers, exerting beneficial effects on plant richness and soil mixing (Pang & Guo, 2018; Yu et al., 2017; Zhou et al., 2018) . In fact, they have been considered as a keystone animal of grassland ecosystems (Hogan, 2010) . However, at an average density of 429 individuals/km 2 on the QTP, they consume about 136 million tons of fresh grass annually, enough to feed 150 million sheep (Sun et al., 2016) . Apart from direct consumption of fresh forage, pikas also affect AGB via burrowing and gnawing, both of which can damage grassland by piling the loosened soil atop healthy grasses and creating bare patches, thus accelerating grassland degradation which is hard for plant regrowth and recovery. Li, Perry, et al., 2013; Qin et al., 2020; Wu et al., 2015) . All ecosystems are impacted by global climate change (McMichael et al., 2004) , but alpine grasslands on the QTP are especially sensitive to changing temperature and precipitation Zha et al., 2005; Zhong et al., 2019) . In this region, warmer temperatures can prolong the growing season by advancing the start date and delaying the end date of growth, boosting AGB production (Sun, 2020; Zhang et al., 2013; ; Zhang, Kong, et al., 2018) . However, the dry climate is a major limiting control upon grass growth (Dai et al., 2019; Gao et al., 2016) . From 2004 to 2016, an increase in air temperature by 0.44°C while precipitation remained stable in northern Qinghai Province was accompanied by a significant increase in belowground biomass of the alpine grassland, but AGB hardly changed (Dai et al., 2019) . In general terms, grazing decreases AGB but climate warming increases AGB , both of which are further compounded by pika disturbance. Thus, the achievement of long-term sustainable grazing requires consideration of multiple stressors, assessing impacts of both livestock grazing and pika disturbance, and their interactions, in light of climate change (Brotherton & Joyce, 2015; Li et al., 2019 Li et al., , 2020 Pech et al., 2007) . To date, researchers have studied the effects of different levels of plateau pika disturbance on plant biomass (Pang & Guo, 2018) and the effect of grazing on grassland biomass productivity . However, there are limited understandings of how the inter-annual fluctuations in AGB are related to different levels of grazing and pika disturbances. This study quantifies the singular and combined effects of artificial grazing and pika disturbance severities on AGB and its changes in an alpine grassland on the QTP, assessing the relative importance of both disturbances. The specific objectives are to (a) quantify the effects of three pika disturbance severities (controlled by population) and three mowing-simulated grazing disturbance severities on AGB and its change, (b) evaluate the importance of pika and mowing disturbances on AGB change individually and jointly, and (c) to contextualize the influence of the disturbances on AGB against the backdrop of climate change. The study site was located on a flat floodplain of the Yellow River near the Henan County town center on the QTP (101°47′E, 34°44′N) ( Figure 1 ). It has an altitude of 3743 m above sea level with a plateau continental climate. Annual temperature ranges from −1.3 to 1.6°C, and annual precipitation ranges from 597.1 to 615.5 mm (http://data. cma.cn/). This area has been affected by pikas and overgrazed for decades Zhang, Fan, et al., 2018) . Livestock grazing was simulated via mowing in this study (Tälle et al., 2015) . This enables careful and specific regulation by setting the mower blade to different heights (Jerrentrup et al., 2014) . Pika disturbance severity was made proportional to its population. Both pika (P) and mowing (M) disturbances were set at three levels (none (n), medium (m), and high (h)). Of the nine treatments, one served as the reference, two were subject to the sole effect of pika and mowing, and the remaining four plots were subject to joint influences of both disturbances (Figure 2a) . Each treatment was replicated three times randomly in three rows, resulting in a total of 27 plots, 25 m × 30 m in size. All plots were spatially disconnected with a 5 m corridor between them. Each plot was fenced with a 50 cm high iron sheet above the ground to prevent pikas from roaming among the plots ( Figure 2) . A mesh fence of 50 cm was installed underground (mesh diameter of 2.5 cm × 2.5 cm), deeper than pika's burrowing depth of 20-30 cm to prevent pika intrusion and loss (Wei et al., 2007; Xiao et al., 1981) . Plateau pikas are widely scattered in the alpine grassland of the QTP (Hoffmann et al., 2013) . Small mammals used in this study weighted 130-200 g. Based on previous research in an alpine meadow near this study (Han et al., 2008; Sun, 2008) , the highest pika density (P h ) was set at about 200 individuals/ha, equivalent to 14 individuals per plot. The medium severity (P m ) was set at seven individuals per plot, and no pika plots served as the reference (P n ). Before setting the pika population, traps were deployed in each plot for one month to capture and remove existing pikas, and the reference plots were further treated with poisons to eradicate any residual pikas. All the trapped pikas were later released to the predesignated plots at a sex ratio of 1:1 (male vs female) for the P h plots and 3:4 for the P m plots (Dobson et al., 1998) . Pika number and sex ratio were reassessed annually, and any deviation from the experiment design was remedied accordingly. Grazing disturbance was simulated at three levels of high, moderate, and none by machinery mowing at the beginning of July and August of each year. High mowing severity (M h ) was achieved by clipping almost all the AGB by setting the mower blade to the lowest level (about 2 cm above the ground), medium mowing severity (M m ) was achieved by removing half of the plant (about 10 cm above the ground). M n plots were not mowed. Residues of the mowed plants were removed from the plots. Blade height was set based on the one-third principle to reduce stress to the remaining plants, namely, the mowed grass height was less than one-third of the total grass height at one time (Law et al., 2016) . Thus, it would take 2-4 times (2-3 days interval) for the mowed grass to reach the designed height. To model the effect of mowing on AGB quantitatively, mowing severity was quantified at a scale of 0.0 (M n )-1.0 (M h ), and a value of 0.5 was assigned to M m . The 27 The images were not radiometrically corrected, even though this F I G U R E 1 Experiment site location and treatments (1, 2, 3 are replications) could enhance the accuracy of VI-derived AGB (De Carvalho et al., 2013) . This omission is permissible because VI can be reliably estimated from non-calibrated high-resolution images so long as they are obtained under the same radiometric conditions (Bolten et al., 2015) . After eight 1 m × 1 m quadrats in each plot were selected year. Non-vegetated areas (e.g., denudated patches, pika burrows, and mounds) were mapped from the drone images using Support Vector Machine (SVM) in ArcMap at an overall accuracy higher than 94.0%, with a Kappa value over 0.89. They were used to assess the relationship between disturbance severity and AGB change. AGB was estimated from the drone images using the Red Green Blue Vegetation Index (RGBVI) (Equations 1-4) proposed by Shi et al. (2021) . After the RGBVI of the pixels (about 250) corresponding to the biomass-harvested quarter on the ground was extracted from each RGBVI image, it was used to construct the AGB estimation model from the ground truth data. where n is the number of samples, B is the average measured AGB of the reference plots (G N ), and b is the measured quadrats AGB. t is the total number of pixels in a sampled quarter; VI i is the vegetation index value of pixel i; j denotes mowing intensity; (1 − F j ) k is the grazing intensity reversal factor; k is a constant determined via nonlinear regression (k = 4 produced the best result for the study area, but can be slightly different for other areas). The estimated AGB from the constructed model was modified twice, first by the mowing ratio F (Equations 2 and 3) which can reflect the vegetation canopy vertical properties. The second modification is to remove biomass of the bare ground area by the RGBVI of non-vegetated pixels (Equation 4). They were determined from a newly created mask layer. The pixels were assigned a new value of zero in the output metrics and treated as the ground truth dataset. The coefficient of determination (R 2 ) and root mean square error (RMSE), employed to evaluate the model accuracy, were 0.88 and 19.9 g m −2 , respectively. All analyses, including model regression, ANOVA (Analysis of Variance), the least significant difference (LSD), and Pearson correlation coefficient (r), were implemented via Python scripting (version 3.8). ANOVA was introduced to check the potential differences of AGB between the treatments, and the least significant difference (LSD) was used to test the differences. Annual AGB changes of 2018-2019 and 2019-2020 were calculated to explore the individual and joint effects of the two disturbances. To evaluate the effects of disturbance severities on AGB change, the correlation between disturbance severities and AGB change was assessed using the Pearson correlation coefficient (r) (Hauke & Kossowski, 2011) . r was also tested at each severity of both disturbances (pika population and mowing intensity). Furthermore, to quantify the exact effect of the increasing disturbance severity on annual AGB change, relative AGB change (%) between disturbance severities was calculated (Equations 5 and 6). Field experiment design (a is the design of treatments, b is the installation of fence and weather station) where D denotes AGB change caused by pika (P) or mowing (M) disturbance. In total, two disturbance gradients were assessed: none-tomoderate (D n,m ) in Equation 5 and moderate-to-high (D m,h ) in Equation 6. The former represents the relative AGB change of moderate disturbance severity (m) corresponding to non-disturbance severity (n), and the latter expresses the relative AGB change of severe disturbance (h) relative to moderate disturbance (m). Monthly temperature averaged from −15℃ (December 2019) to 15℃ (July 2020) over the study period ( Figure 3a ). Monthly mean precipitation ranged from 70 to 170 mm, concentrated in summer and early autumn (from June to Sep), and it was much lower in other months. Figure 4 shows that all treatments always had a consistent rank of AGB and a similar level of significant difference throughout the three years (LSD, p < .01). In 2018, the reference plot (P n M n ) had the highest AGB of 289.47 g m −2 (equivalent to 217.11 kg dry matter per plot). The other two non-mowed plots (P m M n and P h M n ) had an equivalent AGB of around 280 g m -−2 , while the moderately mowed plots (M m ) achieved a much lower AGB of about 190 g m −2 , and the heavily mowed plots (M h ) yielded the lowest AGB of around 150 g m −2 . Moreover, AGB was not significantly different between pika disturbance severities at the same mowing intensity, but it varied significantly between the three mowing intensities at every pika disturbance level. This finding still holds true in the other 2 years. The AGB rank and the significant difference were the same regardless of the absolute quantity. 3.3.1 | Annual AGB change under sole disturbance The quantitative relationship between disturbance severity (pika population and mowing intensity) and AGB changes is presented in Table 3 . Pika population had a poor relationship with AGB changes in all treatments, with a Pearson coefficient (r) of −.471 and −.368 for two periods, respectively. However, at a given mowing severity, pika population had a close negative relationship with AGB change, achieving a coefficient from −.779 to −.626 in the first period, and from −.884 to −.733 in the second period. Similar results were also obtained if judged by r between bare area and AGB, and between bare area change and AGB change (Table 4) . Such a relationship is due to the strong relationship (r > .7) between pika population and bare area change (Table 5) . Mowing intensity had different correlations with AGB change in the two periods (Table 3) . It was hardly correlated with AGB change in all treatments in the first period (r = .12), but highly correlated in the second period (r = −.849). Such different impacts in two periods can be explained by the abnormal AGB drop in the reference plot in the first period. However, the relationship between mowing intensity and AGB change was less influenced by pika disturbance severity, given that r at each pika population was equivalent to that of all treatments (.090 < r < .328 in the first period, −.939 < r < −.911 in the second period). Overall, mowing intensity weakened the relationship between pika population and AGB change, but pika disturbance hardly affected the relationship between mowing severity and AGB change. Note: Numbers in the bracket represent actual differences of AGB change (g m −2 ). between mowing disturbance severities Zhang, Kong, et al., 2018) , but the storage of rainfall in the soil can promote plant growth in the next year (Guo et al., 2015; Mueller et al., 2016) . In this study, mowing directly affected the absolute quantity of AGB, and indirectly influenced its pace of change over time. Mowing can promote grass biomass by stimulating plant growth. Accordingly, the mowed meadow can produce more accumulated biomass than the un-mowed meadow, but only slightly (Bernhardt-Römermann et al., 2011; Ten Cate & Gusewell, 2003) . Due to the regulated clipping height, the effect of mowing on plant height is universal throughout all the plots. As the mowed plots have a lower AGB with less tissue for photosynthesis, they are more subject to the effect of climate stresses than un-mowed plots, and it takes a longer time for the remaining biomass to reach a higher level (Bernhardt-Römermann et al., 2011) . Furthermore, if the grass cover was spatially uniform, then the clipped biomass from a plot would also be spatially uniform, but less if the grassland contained more bare patches Qin et al., 2020) . This is demonstrated by the close negative linear relationship between bare area change and AGB change at every mowing severity (Table 4) . Generally, disturbed grassland has less biomass than undisturbed grassland (Fu et al., 2012; Qin et al., 2014; Zhang et al., 2016) . However, mowing had an opposite effect on annual AGB change in the two periods ( Figure 5 ). This result can be explained by the decrease of AGB in the reference plot which had low AGB and short vegetation height in 2019. All the mowed plots had the same mowing intensity over 3 years, but the amount of clipped grass varied from year to year because different plants had their own natural heights in response to climate conditions (Pöyry et al., 2006) . Vertically, the allocation of grassland biomass has a pyramid-like distribution, with a higher biomass density closer to the ground (Tackenberg, 2007) , 40%-60% of the total biomass was distributed in 0-10 cm above the ground, about 30% in 20-30 cm, and less than 20% over the 30 cm height in a meadow in Ireland (Ten Cate & Gusewell, 2003) . Thus, the mowed plots contained more biomass than un-mowed plots, resulting in less annual AGB change. In comparison, the un-mowed plots had different vegetation heights every year, and taller plants had more biomass, leading to a negative annual AGB change in the climate-stressed period. Obviously, mowing was unable to increase annual AGB change. On the contrary, AGB change of the mowed plots was sensitive to climate variability which is the key control upon alpine grassland productivity (Benot et al., 2014; Liu & Chen, 2000; Piseddu et al., 2021; Yang et al., 2012) . Inevitably, mowing-simulated grazing differs from livestock grazing in that livestock preference for palatable species results in an uneven height of grass (Isselstein et al., 2007) . In contrast, mowing leads to a uniform height. Besides, AGB grazed by livestock is mostly recycled back into the ecosystem through excretion (Long et al., 2008) . Also, livestock trampling can influence grass growth and expand bare area by enhancing soil bulk density, also decreasing hydraulic conductivity and respiration rate (Chai et al., 2019) . Generally, livestock grazing has more positive influences on grassland biomass than mowing, but the differences between the two are rather small (Talle et al., 2016) . A comparison between livestock grazing and mowing in a mountain grassland reveals that they have little differences in changing plant species diversity, but livestock grazing has twofold higher AGB loss than mowing under drought tress (Deléglise et al., 2015) . Additionally, mowing based on clipping height is different from livestock grazing which is based on a stable intake amount of forage (Xin-wei et al., 2005) . It is recommended that artificial grazing simulated by mowing should use the removed grass height or biomass as the indicator for artificial grazing intensity. Due to pika's complex feeding habit such as gnawing roots, seeds, and caching forage (Kang et al., 2019) , pika-consumed biomass cannot be factored in the AGB estimation. However, biomass loss associated with bare area enlargement caused by pika can be accurately estimated (Shi et al., 2021) , owing to the strong relationship of bare area change with AGB change (Table 4) . Generally, an adult pika devours 77.3 g of fresh grass/day (Fan et al., 1999) , equivalent to 11.37 kg dry matter per annum. In this study, the biomass loss caused by pika (M n of P n,m gradient for example) was at 13.02-17.88 g m −2 (equivalent to 1.40-1.92 kg per pika annually). This figure is much lower than the directly consumed biomass, but the loss is non-recoverable as pika likely destroyed the grass root system via its burrowing activities (Yu et al., 2017; Zhang et al., 2003) . Furthermore, the fenced plots may not encompass sufficient palatable species of grasses, forcing pika to seek more food via digging and gnawing (Guo et al., 2012; Hogan, 2010) . However, the doubling of pika population from 7 to 14 individuals per plot did not double the annual AGB change because multiple pikas can share food and burrow (Dobson et al., 1998; Fan et al., 1999) . On the other hand, pika population still had a positive relationship with bare area change (Table 5 ), but the increase in the bare area slows down at a higher pika population (Sun et al., 2016) . A good correlation existed between pika population (bare area) and AGB change at each mowing severity (Tables 3 and 4 ). However, these relationships disappeared when all the plots were included in the regression. This and AGB that was decreased mainly by livestock, but their work did not explore the relationship to each grazing intensity. Thus, the effect of pika disturbance on alpine grassland biomass was subject to mowing severity that overwhelmed the relationship between pika disturbance and AGB change. Conversely, some authors note that pika is a keystone rodent for increasing grassland biomass (Liu et al., 2017; Smith & Foggin, 1999) , with herbs occupying more than 90% of bare patches Wei et al., 2007) . Despite this increase, however, most of the newly grown biomass was inedible for livestock, and some even poisonous (Wen et al., 2010) . Further study is required to investigate the change of plant species in different pika disturbance severities. A more severe disturbance level can amplify the effect of another disturbance on AGB decrease, and the joint effects of both mowing and pika disturbances can amplify the negative effect of either disturbance on AGB change. In particular, a high pika disturbance level boosted biomass loss as the mowing intensity rose from moderate to high. In greater numbers, pikas face more competition for food and living space, which drives them to create more bare patches (e.g., gnawing at grassroots), leading to more biomass loss . Also, the impact of pika on landscape fragmentation was higher than mowing in this experiment in another study . These effects were amplified in the second period because bare area expanded rapidly once it was initiated (Li, Perry, et al., 2013 Piao et al., 2011; Xu et al., 2013) . On the other hand, plant height can impact pika burrowing activities in that lower plants can lead to a higher burrow density and more bare patches , because pikas will dig more holes to hide and escape from predators rapidly (Fan & Zhang, 1996; Wangdwei et al., 2013) . Furthermore, pikas will cut down plants to create an open view for safety checking and danger warning (Smith & Foggin, 1999) . Additionally, mowing intensity has significant impacts on grassland AGB under different climate conditions. In this study, the relative importance of rainfall was the highest at high mowing intensity and the contribution of temperature was greatest at low mowing intensity (Bernhardt-Römermann et al., 2011) . However, research by Liu et al. (2018) reported that mowing was able to decrease the damaged area by 14% in a grassland in Inner Mongolia, but differences between varied mowing intensities were not considered in this work. In this study, the obtained results demonstrate that severe mowing is conducive to rodent outbreaks. Severe mowing led to more loss of AGB than severe pika disturbance in heavily disturbed plots. Severe pika disturbance caused no obvious differences in AGB change between moderately mowed and intact plots, but had much higher negative effects on the AGB change in severely mowed plots. Therefore, reducing mowing severity is the primary choice to improve AGB in the plots heavily disturbed by mowing and pika, and pikas do not pose a threat if the grassland is mowed at a low intensity. In reality, livestock and pika have coexisted for thousands of years in grazing-adapted alpine grassland (Badingquiying et al., 2018) . However, both grazing and rodent disturbances had negative impacts on grassland biomass (Chen et al., 2007; Pang & Guo, 2018) , and grazing has been reported as an important factor in causing pika outbreaks (Foster et al., 2014; Su et al., 2016) . Thus, based on findings from this study, it is more effective to reduce grazing intensity than control the pika population to restore the severely disturbed grassland. It will be interesting to see if grazing combined with pika disturbance will have the same effects on AGB change in a natural alpine grassland as those found in this controlled experimental study. Inevitably, alpine grassland AGB change is subjective to the complex interplay and associated uncertainties of multiple stressors; besides, livestock and pika disturbances, including interactions between climate and land-use change (Carlyle et al., 2014; Hovenden et al., 2014) , environmental legacy effects (De Boeck et al., 2018; Gong et al., 2020) , and stages of degradation or recovery of grasslands also play important roles in AGB and its change. This controlled study quantified the singular and joint effects of mowing and pika disturbances on alpine grassland AGB changes on the QTP for the first time. Pika-induced AGB change through bare area expansion was much smaller than their daily consumption, but their damage on grassland was hardly recoverable, and this effect can be amplified with time. Moreover, pika's effects on AGB change were always overwhelmed by mowing severities, but it can be identified by inspecting it into each mowing intensity. Mowing intensity weakened the relationship between pika population and AGB change, but pika disturbance hardly affected the relationship between mowing severity and AGB change. The joint effect of both disturbances amplified the negative effect of individual disturbances on AGB. The damage of high mowing severity on AGB loss was larger than high pika severity in the heavily disturbed grassland. Pika disturbances induced less difference in AGB change if mowing did not exceed the moderate level. Thus, reducing mowing intensity is more effective than controlling pika population for boosting AGB production in heavily disturbed plots. Based on the findings of this study, we recommend that reducing livestock grazing intensity rather than controlling pika population should be the primary method to improve AGB in the real livestock grazing grasslands on the QTP. This study was financially supported by the Project of the Science The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper. The biomass and climate data are available in Dryad: https://doi. org/10.5061/dryad.p5hqb zkr2. https://orcid.org/0000-0003-4520-1458 Gary Brierley https://orcid.org/0000-0002-1310-1105 Summer habitat use of plateau pikas (Ochotona curzoniae) in response to winter livestock grazing in the alpine steppe Qinghai-Tibetan Plateau Feasibility study of using non-calibrated UAV-based RGB imagery for grassland monitoring: Case study at the Rengen Long-term Grassland Experiment (RGE) Stronger short-term effects of mowing than extreme summer weather on a subalpine grassland Explaining grassland biomass -the contribution of climate, species and functional diversity depends on fertilization and mowing frequency Grazing intensity and ecosystem processes in a northern mixed-grass prairie Extreme climate events and wet grasslands: Plant traits for ecological resilience Response of grassland biomass production to simulated climate change and clipping along an elevation gradient Effects of yak and Tibetan sheep trampling on soil properties in the northeastern Qinghai-Tibetan Plateau Changes in vegetation photosynthetic activity trends across the Asia-Pacific region over the last three decades The contribution of plateau pika disturbance and erosion on patchy alpine grassland soil on the Qinghai-Tibetan Plateau: Implications for grassland restoration Model analysis of grazing effect on above-ground biomass and above-ground net primary production of a Mongolian grassland ecosystem Value of ecosystem services in China Thirteen-year variation in biomass allocation under climate change in an alpine Kobresia meadow, northern Qinghai-Tibetan Plateau Legacy effects of climate extremes in alpine grassland Radiometric normalization of temporal images combining automatic detection of pseudo-invariant features from the distance and similarity spectral measures, density scatterplot analysis, and robust regression The effects of micro-habitats and grazing intensity on the vegetation of burial mounds in the Kazakh steppes Drought-induced shifts in plants traits, yields and nutritive value under realistic grazing and mowing managements in a mountain grassland Social and ecological influences on dispersal and philopatry in the plateau pika (Ochotona curzoniae) Enhancing sustainability of grassland ecosystems through ecological restoration and grazing management in an era of climate change on Qinghai-Tibetan Plateau Enhancing the resilience of coupled human and natural systems of alpine rangelands on the Qinghai-Tibetan Plateau Distribution and cause of plateau pika (Ochotona curzoniae) burrows in Henan Mongolian Autonomous county, Qinghai Province Foraging behavior of Ochotona curzoniae and Ochotona daurica and their adaptation to habitat Rodent Pest Management in the Qinghai-Tibet Alpine Meadow Ecosystem Effects of large native herbivores on other animals Modeling the responses of grassland vegetation coverage to grazing disturbance in an alpine meadow Climatic change controls productivity variation in global grasslands Legacy effects of precipitation amount and frequency on the aboveground plant biomass of a semi-arid grassland Changes in vegetation parameters and soil nutrients along degradation and recovery successions on alpine grasslands of the Tibetan plateau Contrasting responses of gross primary productivity to precipitation events in a water-limited and a temperature-limited grassland ecosystem Response of alpine meadow communities to burrow density changes of plateau pika (Ochotona curzoniae) in the Qinghai-Tibet Plateau Rodent damage assessment on the plateau pika Rangeland degradation on the Qinghai-Tibetan plateau: A review of the evidence of its magnitude and causes Comparison of values of Pearson's and Spearman's correlation coefficients on the same sets of data Increases in the net above-ground primary production of a salt-marsh forage grass: A test of the predictions of the herbivore-optimization model Mammals of China The plateau pika: A keystone engineer on the Tibetan Plateau Grazing intensity: Critique and approach Seasonal not annual rainfall determines grassland biomass response to carbon dioxide Climate changes in temperature and precipitation extremes in an alpine grassland of Central Asia Effects of livestock breed and grazing intensity on biodiversity and production in grazing systems. 1. Nutritive value of herbage and livestock performance. Grass and Forage Science Grazing intensity affects insect diversity via sward structure and heterogeneity in a long-term experiment Food habits and its seasonal changes of plateau pika (Ochotona curzniae) in Gannan meadow Selecting turfgrasses and mowing practices that reduce mowing requirements The rangeland degradation in North China and its preventive strategy Micro-scale fragmentation of the alpine meadow landscape on the Qinghai-Tibet Plateau under external disturbances. Catena Short-term effects of experimental trampling on alpine grasslands in Shangri-la Patterns and thresholds of grazing-induced changes in community structure and ecosystem functioning: Species-level responses and the critical role of species traits Rangeland degradation on the Qinghai-Tibet plateau: Implications for rehabilitation. Land Degradation & Development Restoration prospects for Heitutan degraded grassland in the Sanjiangyuan Livestock grazing impacts on plateau pika (Ochotona curzoniae) vary by species identity The interaction between poisonous plants and soil quality in response to grassland degradation in the alpine region of the Qinghai-Tibetan Plateau One-year grazing exclusion remarkably restores degraded alpine meadow at Zoige, eastern Tibetan Plateau Climatic warming in the Tibetan Plateau during recent decades Relationships between plateau pika (Ochotona curzoniae) densities and biomass and biodiversity indices of alpine meadow steppe on the Qinghai-Tibet Plateau China Effects of nitrogen addition and mowing on rodent damage in an Inner Mongolian steppe The yak grazing system on the Qinghai-Tibetan plateau and its status Climate change alters interannual variation of grassland aboveground productivity: Evidence from a 22-year measurement series in the Inner Mongolian grassland Global climate change Impacts of warming and elevated CO 2 on a semi-arid grassland are nonadditive, shift with precipitation, and reverse over time Plant production, and carbon and nitrogen source pools, are strongly intensified by experimental warming in alpine ecosystems in the Qinghai-Tibet Plateau Plateau pika disturbances alter plant productivity and soil nutrients in alpine meadows of the Qinghai-Tibetan Plateau Effects of plateau pika disturbance levels on the plant diversity and biomass of an alpine meadow Population dynamics and responses to management of plateau pikas Ochotona curzoniae Changes in satellite-derived vegetation growth trend in temperate and boreal Eurasia from 1982 to Mowing and warming effects on grassland species richness and harvested biomass: Meta-analyses Different responses of plants and herbivore insects to a gradient of vegetation height: An indicator of the vertebrate grazing intensity and successional age Effects of plateau pikas' foraging and burrowing activities on vegetation biomass and soil organic carbon of alpine grasslands Responses of typical grasslands in a semi-arid basin on the Qinghai-Tibetan Plateau to climate change and disturbances Improved estimation of aboveground biomass of disturbed grassland through including bare ground and grazing intensity The plateau pika (Ochotona curzoniae) is a keystone species for biodiversity on the Tibetan plateau Effects of livestock grazing on rodents in grassland ecosystems Effects of burrowing plateau Pika (Ochotona curzoniae) densities on primary productivity and soil resource characteristics in alpine meadow Plateau pika population survey and its control threshold in the alpine meadow ecosystems of the Tibetan plateau Effects of rodent activities on primary productivity and soil physical characteristics in alpine meadow The response of vegetation dynamics of the different alpine grassland types to temperature and precipitation on the Tibetan Plateau Delayed autumn leaf senescence date prolongs the growing season length of herbaceous plants on the Qinghai-Tibetan Plateau A new method for non-destructive measurement of biomass, growth rates, vertical biomass distribution and dry matter content based on digital image analysis Grazing vs. mowing: A meta-analysis of biodiversity benefits for grassland management The conservation benefit of mowing vs grazing for management of speciesrich grasslands: A multi-site, multi-year field experiment Remote sensing of crop production in China by production efficiency models: Models comparisons, estimates and uncertainties Impact of mowing regime on species richness and biomass of a limestone hay meadow in Ireland Vegetation type and grazing intensity jointly shape grazing effects on grassland biodiversity Livestock type is more crucial than grazing intensity: Traditional cattle and sheep grazing in short-grass steppes. Land Degradation & Development Effect of grazing on the plant functional group diversity and community biomass and their relationship along a precipitation gradient in Inner Mongolia Steppe Grass height and transmission ecology of Echinococcus multilocularisin Tibetan communities Demographic responses of plateau pikas to vegetation cover and land use in the Tibet Autonomous Region Soil erosion and vegetation succession in alpine Kobresia steppe meadow caused by plateau pika-A case study of Nagqu County The construction of Grassland Degradation Index for Alpine Meadow in Qinghai-Tibetan Plateau Asynchrony among local communities stabilises ecosystem function of metacommunities Limitation of grassland productivity by low temperature and seasonality of growth Impacts of grazing exclusion on productivity partitioning along regional plant diversity and climatic gradients in Tibetan alpine grasslands Impacts of burrows and mounds formed by plateau rodents on plant species diversity on the Qinghai-Tibetan Plateau The distribution of the plateau pika and its effect on grass Kobresia pygmaea in the Tianjun and Kangyang regions Study on the forage intake and dry-matter digestion coefficient in the different season in HaiNan QingHai Temperature and vegetation seasonality diminishment over northern lands Grazing significantly increases root shoot ratio but decreases soil organic carbon in Qinghai-Tibetan Plateau grasslands: A hierarchical meta-analysis. Land Degradation & Development Effects of fertilization and mowing on community biomass compensation in the eastern alpine meadow of Tibetan plateau Soil disturbance and disturbance intensity: Response of soil nutrient concentrations of alpine meadow to plateau pika bioturbation in the Qinghai-Tibetan Plateau Grassland productivity in an alpine environment in response to climate change Green-up dates in the Tibetan Plateau have continuously advanced from 1982 to 2011 Livestockcarrying capacity and overgrazing status of alpine grassland in the Three-River Headwaters region A method for calculating the suitable monthly carrying capacity of seasonal pasture-Taking heriheng Village Vegetation phenology on the Qinghai-Tibetan Plateau and its response to climate change Grazing induced changes in plant diversity is a critical factor controlling grassland productivity in the Desert Steppe Responses of alpine vegetation and soils to the disturbance of plateau pika (Ochotona curzoniae) at burrow level on the Qinghai-Tibetan Plateau of China Effects of grazing and climate warming on plant diversity, productivity and living state in the alpine rangelands and cultivated grasslands of the Qinghai-Tibetan Plateau Burrowing rodents as ecosystem engineers: The ecology and management of plateau zokors Myospalax fontanierii in alpine meadow ecosystems on the Tibetan Plateau Climate change trends and impacts on vegetation greening over the Tibetan Plateau Impact of plateau pikas (Ochotona curzoniae) on soil properties and nitrous oxide fluxes on the Qinghai-Tibetan Plateau Estimating net primary productivity of terrestrial vegetation based on GIS and RS: As case study in Inner Mongolia Effects of disturbances on aboveground biomass of alpine meadow in the Yellow River Source Zone, Western China