Quantifying Degradation Classiﬁcations on Alpine Grassland in the Lhasa River Basin, Qinghai-Tibetan Plateau

: The Qinghai-Tibetan Plateau (QTP) has the world’s largest alpine grassland ecosystem. The QTP ecosystem is extremely fragile and su ﬀ ers continuous degradation. An accurate determination of the status of alpine grassland is the ﬁrst crucial step in monitoring its degradation. A novel method combining ﬁeld survey with remote sensing information based on ecological indicators is proposed. The degradation classiﬁcation of alpine grassland was identiﬁed by multivariate hierarchical analysis based on 270 ﬁeld plots. The spatial pattern of alpine grassland degradation was mapped by determining remote sensing variables that corresponded to ﬁeld indicators of the degradation classiﬁcation system. The results showed that clustering analysis divided the degradation classiﬁcation of alpine grassland into ﬁve classes: Non-Degraded (ND), Slightly Degraded (SLD), Moderately Degraded (MD), Severely Degraded (SD), and Extremely Degraded (ED). The most signiﬁcant factors for alpine grassland degradation included the dominance of Cyperaceae plants, soil total nitrogen content, soil organic carbon content, soil total carbon content, soil bulk density, soil pH, dominance of miscellaneous plants, and elevation among all 17 variables. The assessment and mapping of alpine grassland degradation provide an important basis for alpine grassland protection and management, particularly at a large scale. Alpine grasslands with different classes of degradation require different corresponding restoration and management measures [43], and the quantitative grading of degradation contributes to a more scientific alpine grassland management. The ecological thresholds


Introduction
Grassland is an important part of terrestrial ecosystems, with natural vegetation types including steppe, meadow, marsh, tundra, savanna, desert, and woodland [1]. Grassland plays important ecological and socio-economic roles, such as affecting biodiversity by evolving grass species that shape grassland environments [2] and storing carbon as carbon sinks and impacting livestock by producing forage [3]. Grassland degradation manifests in the reduction of grassland productivity, deterioration of the environment, and loss of biodiversity, suffering from natural processes and human activities [4]. Grassland degradation is a global problem as demonstrated by the fact that nearly half of the world's grasslands are degraded [5]. Additionally, 90% of grassland in China has been degraded to some extent [6]. The alpine grassland of the Qinghai-Tibetan Plateau (QTP) has suffered from increasingly severe degradation in recent decades due to livestock over-grazing, rodent damage, and climate changes [3,7]. Herbivore (mainly yaks and sheep) grazing affects the productivity and the function of

Study Area
The Lhasa River basin is located in the southern QTP (90 • 05 -93 • 20 E, 29 • 20 -31 • 15 N) and has an area of 3.26 × 10 4 km 2 ( Figure 1). The terrain declines in altitude from north to south, and the elevation ranges from 3598 m to 7074 m, with geographical characteristics of mountains and valley plains. The study area is within the plateau temperate semi-arid monsoon climate zone. The mean annual temperature ranges from −7.1 • C to 9.2 • C, and the multiyear average precipitation ranges from 340 mm to 700 mm. The vegetation types mainly include alpine meadow, alpine steppe, mountain shrub grassland, shrub, and sparse mountain forest [22]. The soil types are mainly leptosols, cambisols, gleysols, and phaeozems. This area has suffered from severe grassland degradation due to the long-term lack of scientific management and the under-developed utilization pattern of grassland resources, which includes over-grazing and rodent damage. The alpine grassland ecosystem has been destroyed and sustainable development of alpine grassland livestock husbandry has been restricted.

Field Data
We selected 64 field sites of representative alpine grassland that covered all types of alpine grassland in the study area. Five subsites representing "Non-Degraded" (ND), "Slightly Degraded" (SLD), "Moderately Degraded" (MD), "Severely Degraded" (SD), and "Extremely Degraded" (ED) alpine grassland were designated at each site in reference to the alpine grassland coverage described by Wang et al. [23]. One plot (1m × 1 m) was set randomly in each subsite, and the plots were separated by at least 10 m buffer. Twenty-four sites did not have five subsites (22 of them had three subsites, and two had two subsites), yielding a total of 270 plots across the region.
All field data were measured during the growing season between July and August in 2017. The latitude, longitude, and elevation of the subsites were recorded by a Global Positioning System (GARMIN GPSMAP ® 631sc). The slope and aspect were measured using an inclinometer and a compass. In each plot, vegetation parameters were measured, including total vegetation coverage, bare land coverage, species richness, coverage of each species, height of each species, and dominance of plant functional groups. The total vegetation coverage and bare land coverage were calculated by fisheye lens images. Species richness was simply the total number of species per plot. The coverage of each species was calculated from the number of individuals of each one. The height of each species was calculated as the average from five random individuals. The plant groups in plots were apportioned into four functional groups containing Cyperaceae plants, Poaceae plants, miscellaneous plants, and inedible plants [24]. The dominance index is based on Importance Values (IVs) of species, and the proportion of the total IVs was calculated as follows [25]: where pi is the important value in the plot, S is the total number of species, Ci is the coverage of species i, Hi is the height of species i, and Fi is the frequency of species i. Soil samples (0-20 cm depth) were collected in each plot and used for measurements of the soil moisture content, soil bulk density, soil pH, soil total carbon, soil organic carbon, soil total nitrogen,

Field Data
We selected 64 field sites of representative alpine grassland that covered all types of alpine grassland in the study area. Five subsites representing "Non-Degraded" (ND), "Slightly Degraded" (SLD), "Moderately Degraded" (MD), "Severely Degraded" (SD), and "Extremely Degraded" (ED) alpine grassland were designated at each site in reference to the alpine grassland coverage described by Wang et al. [23]. One plot (1 m × 1 m) was set randomly in each subsite, and the plots were separated by at least 10 m buffer. Twenty-four sites did not have five subsites (22 of them had three subsites, and two had two subsites), yielding a total of 270 plots across the region.
All field data were measured during the growing season between July and August in 2017. The latitude, longitude, and elevation of the subsites were recorded by a Global Positioning System (GARMIN GPSMAP ® 631sc). The slope and aspect were measured using an inclinometer and a compass. In each plot, vegetation parameters were measured, including total vegetation coverage, bare land coverage, species richness, coverage of each species, height of each species, and dominance of plant functional groups. The total vegetation coverage and bare land coverage were calculated by fisheye lens images. Species richness was simply the total number of species per plot. The coverage of each species was calculated from the number of individuals of each one. The height of each species was calculated as the average from five random individuals. The plant groups in plots were apportioned into four functional groups containing Cyperaceae plants, Poaceae plants, miscellaneous plants, and inedible plants [24]. The dominance index is based on Importance Values (IVs) of species, and the proportion of the total IVs was calculated as follows [25]: where p i is the important value in the plot, S is the total number of species, C i is the coverage of species i, H i is the height of species i, and F i is the frequency of species i. Soil samples (0-20 cm depth) were collected in each plot and used for measurements of the soil moisture content, soil bulk density, soil pH, soil total carbon, soil organic carbon, soil total nitrogen, soil available nitrogen, soil available phosphorus, and soil available potassium. The soil moisture content and soil bulk density were determined by the oven drying method, according to Thomasson [26]. The soil pH was measured on a 1:1 water/soil suspension using a pH meter [27]. Soil total carbon and soil total nitrogen were measured using the flash dynamic combustion method, gas chromatographic separation, and thermal conductivity detection system (vario EL cube Elementar, Germany), respectively [28]. Soil organic carbon was determined by the dichromate oxidation method [29]. Soil available nitrogen was measured using the alkali-hydrolysis reduction diffusion method [30,31]. Soil available phosphorus was extracted with the ammonium bicarbonate method, and soil available potassium was extracted using the acetamide extraction method [32].

Remote Sensing Data
Landsat 8 OLI satellite images were acquired on 10 July 2017, corresponding to the field survey time as much as possible under the condition of no clouds. Eight Landsat 8 images (path/row number: 138/038, 138/039, 138/040, 137/038, 137/039, 137/040, 136/038, 136/039) were downloaded from the United States Geological Survey (USGS) website (http://earthexplorer.usgs.gov) at a spatial resolution of 30 m. These images were processed by radiance calibration and atmospheric correction in ENVI 5.1. In this study, we used only two bands: Band 4 (636-673 nm) and Band 5 (851-879 nm). The vegetation indices considered in this case, including the normalized difference vegetation index (NDVI), soil adjusted vegetation index (SAVI), and renormalized difference vegetation index (RDVI) (according to Davidson et al.) were calculated as follows: where NIR is the radiance in the near infrared region band (NIR = Band5), R is the radiance in red band (R = Band 4), and L is the soil correction factor. In this study, the L value equals 0.5 [33]. The elevation data were derived from ASTER Global DEM (GDEM) with a spatial resolution of 30 m (http://srtm.csi.cgiar.org). The soil type was derived from spatial distribution data of soil types in China at a scale of 1/1000000 (http://www.resdc.cn). The coverage map of bare land was derived from the MOD44B Version 6 Vegetation Continuous Fields (VCF) product with a spatial resolution of 250 m (https://lpdaac.usgs.gov/dataset_discovery/modis/modis_products_table/mod44b_v006).

Statistical Analysis
All statistical analyses were performed in R version 3.5.0 with the Vegan and the Rpart libraries [34]. Euclidean distance matrices were used for community dissimilarity analysis, and hierarchical clustering using Ward's method was calculated on this basis [35]. Nonmetric multidimensional scaling (NMDS) analysis was used to test the differences in plant community composition among different degraded levels of alpine grassland. The "envfit" function software package was used to determine which biophysical variables correlated well with ordination space. Seventeen field-derived biophysical variables were selected and normalized to range from zero (0) to one (1) ( Table 1). The main factors affecting the ordination were chosen by excluding collinear factors. Classification and regression tree (CART) analysis was used to characterize the main biophysical factors of subsites classified as ND, SLD, MD, SD, and ED. The tree determined the thresholds of variables by encoding a set of decision rules in the form of if-then statements [20]. The CART model was pruned by reducing its size to the minimum of the cross-validation error [36].

Spatial Pattern of Alpine Grassland Degradation
Based on the main factors of alpine grassland degradation determined by multivariate statistical analysis, the corresponding GIS layers were selected ( Table 2). Among them, the NDVI, SAVI, and RDVI were chosen as vegetation index layers because of their suitability for alpine grassland degradation according to Zha [4]. The indicator thresholds of alpine grassland degradation obtained by CART analysis consummated the degradation classifications at the field scale. The sites with different degradation classes were placed on the GIS layers, and the pixel values of the sites in each layer were extracted. Thus, the ranges of the pixel values of each layer in the different degradation classes were obtained. The grading layers were in accordance with the locations of the plots in the different classes in layers. We divided the degradation class of alpine grassland according to the values ranging from low (1) to high (5). The most frequent score was the overall assessment score of degradation at each pixel scale ( Figure 2).

Spatial Pattern of Alpine Grassland Degradation
Based on the main factors of alpine grassland degradation determined by multivariate statistical analysis, the corresponding GIS layers were selected ( Table 2). Among them, the NDVI, SAVI, and RDVI were chosen as vegetation index layers because of their suitability for alpine grassland degradation according to Zha [4]. The indicator thresholds of alpine grassland degradation obtained by CART analysis consummated the degradation classifications at the field scale. The sites with different degradation classes were placed on the GIS layers, and the pixel values of the sites in each layer were extracted. Thus, the ranges of the pixel values of each layer in the different degradation classes were obtained. The grading layers were in accordance with the locations of the plots in the different classes in layers. We divided the degradation class of alpine grassland according to the values ranging from low (1) to high (5). The most frequent score was the overall assessment score of degradation at each pixel scale ( Figure 2).

Classification of Subsites to Identify Degradation Classes
The cluster analysis showed that the 270 plots could be classified into five groups ( Figure 3). They include Cluster 2, Cluster 5, Cluster 4, Cluster 3, and Cluster 1, represented the ND, SLD, MD, SD, and ED classes of alpine grassland degradation, respectively. Significant differences in the biophysical variables were found among the different degradation classes ( Table 3). The total vegetation coverage of SLD and MD decreased by 9.31% and 13.35%, respectively, compared with that of ND. SD and ED showed significantly decreased total vegetation coverage by 36.84% and 70.69%, respectively, relative to ND. Compared with ND, SLD and MD exhibited significantly increased dominance of Cyperaceae plants by 137.73% and 101.64%. However, the dominance of Cyperaceae plants in SD and ED was lower by 26.89% and 61.21% compared with that in ND, and significantly lower by 69.25% and 83.68% compared with that in SLD. ND and SLD showed no significant difference in the dominance of Poaceae plants, but MD, SD, and ED increased by 14.88%, 116.38%, and 369.74% compared with SLD, while ED exhibited a significant increase. The dominance of inedible plants exhibited a similar trend to the dominance of Poaceae plants in different degradation classes. The dominance of miscellaneous plants in SLD and MD was significantly lower than that of ND by 42.06% and 41.02%, respectively. Bare land cover showed an increasing trend with increasing degradation level as a whole, and the bare land cover was significantly increased in ED compared with ND, SLD, MD, and SD. The soil water content, soil total carbon content, soil organic carbon content, soil total nitrogen content, soil available nitrogen content, and soil available potassium content decreased with increasing degradation level, while the soil total carbon content exhibited a significant decrease. The soil bulk density and soil pH increased with the aggravation of alpine grassland degradation. The soil available phosphorus content was not significantly different in the different degradation classes. The changes in elevation and slope were not highly significant. Consequently, the subsites can be classified into different clusters through a statistical analysis method based on the biophysical variables.

Main Biophysical Factors of Classified Subsites
The NMDS analysis showed the ordination space of all subsites clustered based on the degree of alpine grassland degradation ( Figure 4). The goodness-of-fit of the ordination in linear regression and non-linear regression was high with R 2 values of 0.958 and 0.871, respectively. All the variables were correlated with the ordination pattern of the subsites, with the exception of the soil available phosphorus content and soil available potassium content. Different variables did not have the same influence on the subsite ordination. The subsites plotted toward the ND class were correlated with elevation. SLD and MD were correlated with the dominance of Cyperaceae plants, total vegetation coverage, soil water content, soil total carbon content, soil total nitrogen content, soil organic carbon content, and soil available nitrogen content. These two classes showed the highest correlation with the dominance of Cyperaceae plants. SD and ED were correlated with the dominance of inedible plants, dominance of Poaceae plants, dominance of miscellaneous plants, soil bulk density, soil pH, bare land cover, and slope. They were highly correlated with the dominance of inedible plants, dominance of Poaceae plants, soil bulk density, and soil pH. Collinear variables were found among all of these variables. Bare land cover had significant collinearity with slope. Soil available phosphorus content and soil available potassium content were strongly collinear. Considering collinearity and weak correlations, the main biophysical factors of the classified subsites included all variables except slope, soil available phosphorus content, and soil available potassium content.

Thresholds of Variables Affecting Degradation
A pruned classification tree can encode a set of decision rules to predict the classification of five degraded alpine grassland types, which can be used to determine the main predictor variables and their thresholds ( Figure 5). For example, the right most decision rule can be translated as follows: "when the dominance of Cyperaceae plants is greater than 27.85% and the soil organic carbon content is less than 0.88%, seven subsites were classified as extremely degraded". Figure 5 shows that the main predictor variables included the dominance of Cyperaceae plants, soil total nitrogen content, soil organic carbon content, soil total carbon content, soil bulk density, soil pH, dominance of

Thresholds of Variables Affecting Degradation
A pruned classification tree can encode a set of decision rules to predict the classification of five degraded alpine grassland types, which can be used to determine the main predictor variables and their thresholds ( Figure 5). For example, the right most decision rule can be translated as follows: "when the dominance of Cyperaceae plants is greater than 27.85% and the soil organic carbon content is less than 0.88%, seven subsites were classified as extremely degraded". Figure 5 shows that the main predictor variables included the dominance of Cyperaceae plants, soil total nitrogen content, soil organic carbon content, soil total carbon content, soil bulk density, soil pH, dominance of miscellaneous plants, and elevation, as noncollinear variables with relatively strong effects in the NMDS analysis. The thresholds determined by the classification tree partitioned the degradation classes, obtaining quantitative grading criteria of alpine grassland degradation.  Figure 6 shows the spatial pattern of alpine grassland degradation obtained from quantitative classification based on both fieldwork and remote sensing data. The grassland with an area of 2.26×10 4 km 2 was divided into five levels. The criteria for alpine grassland degradation are presented in Table 2. Of the alpine grassland area, 4.34% was classified as ND, 2.52% as SLD, 49.30% as MD, 27.79% as SD, and 16.06% as ED. The ND was mainly located in the valley of Maizhokunggar County, which is in the southeastern region of the Lhasa River basin, comprising 14.20% of the area. The SD of alpine grassland comprised 5.00% and was found in the Medica wetlands in northeastern Lhari County. The MD was widely distributed throughout the entire basin, occupying 62.36% of Datse, 58.08% of Qushui, 55.64% of Chengguan District of Lhasa city, 52.64% of Damxung, 52.53% of Duilongdeiqin, 50.36% of Nagchu, 47.56% of Sangri, 46.19% of Maizhokunggar, and 42.17% of Lhari. In the county of Lingdrub in this basin, 35.99% of the area was classified as SD, followed by Damxung County, where the severe class occupied 34.08%. The ED of alpine grassland degradation was scattered in the study area, in which Lhari and Maizhokunggar counties had larger proportions of 29.90% and 17.39%, respectively.  Figure 6 shows the spatial pattern of alpine grassland degradation obtained from quantitative classification based on both fieldwork and remote sensing data. The grassland with an area of 2.26 × 10 4 km 2 was divided into five levels. The criteria for alpine grassland degradation are presented in Table 2. Of the alpine grassland area, 4.34% was classified as ND, 2.52% as SLD, 49.30% as MD, 27.79% as SD, and 16.06% as ED. The ND was mainly located in the valley of Maizhokunggar County, which is in the southeastern region of the Lhasa River basin, comprising 14.20% of the area. The SD of alpine grassland comprised 5.00% and was found in the Medica wetlands in northeastern Lhari County. The MD was widely distributed throughout the entire basin, occupying 62.36% of Datse, 58.08% of Qushui, 55.64% of Chengguan District of Lhasa city, 52.64% of Damxung, 52.53% of Duilongdeiqin, 50.36% of Nagchu, 47.56% of Sangri, 46.19% of Maizhokunggar, and 42.17% of Lhari. In the county of Lingdrub in this basin, 35.99% of the area was classified as SD, followed by Damxung County, where the severe class occupied 34.08%. The ED of alpine grassland degradation was scattered in the study area, in which Lhari and Maizhokunggar counties had larger proportions of 29.90% and 17.39%, respectively.

Validity of the Classification Method for Alpine Grassland Degradation
The NMDS is an unconstrained ordination analysis based on any type of distance matrix. Bray-Curtis dissimilarity matrices are commonly used in NMDS ordination [37]. For this multivariate statistical method based on distance matrices, the selection of an appropriate distance matrix directly affects the rationality of the results [38]. In this study, the cophenetic correlation coefficients of the clustering model were compared, and the results showed that the Euclidean distance had a value of 0.72 and a Bray-Curtis distance of 0.58. Thus, the Euclidean distance has a better clustering effect than the Bray-Curtis distance. The Euclidean distance was used in NMDS ordination [9] when it was more robust than Bray-Curtis.
The vegetation index (VI) is the result of an arithmetic operation between the pixel values of two or more spectral bands, which are correlated with plant parameters [39,40], and reveals the status of grassland degradation. The appropriate VIs of the QTP included NDVI, SAVI, and RDVI, which were selected by univariate linear regression analysis, with the larger R 2 values including the percent grass cover and the proportion of unpalatable grass as dependent variables [4]. NDVI is widely used for monitoring natural biological communities and agricultural ecosystems due to its easy accessibility and commonality. However, for sparse vegetation, the soil background contributes greatly to the total spectral reflectance [14], and it is not appropriate to evaluate grassland degradation using only NDVI. SAVI was proposed to compensate for the corresponding soil effect, reducing the significant contribution of soil reflectivity [41]. In high vegetation cover areas, NDVI is easily saturated, and RDVI can avoid the problems applied in the cases of high and low vegetation coverage, thus better adapting to areas with large disparities in coverage [42].

Quantitative Classifications of Alpine Grassland Degradation
The assessment of grassland degradation states is an important prerequisite for restoring and treating degraded grassland. Quantitative classifications of degradation can represent the degradation states more accurately. Alpine grasslands with different classes of degradation require different corresponding restoration and management measures [43], and the quantitative grading of degradation contributes to a more scientific alpine grassland management. The ecological thresholds

Validity of the Classification Method for Alpine Grassland Degradation
The NMDS is an unconstrained ordination analysis based on any type of distance matrix. Bray-Curtis dissimilarity matrices are commonly used in NMDS ordination [37]. For this multivariate statistical method based on distance matrices, the selection of an appropriate distance matrix directly affects the rationality of the results [38]. In this study, the cophenetic correlation coefficients of the clustering model were compared, and the results showed that the Euclidean distance had a value of 0.72 and a Bray-Curtis distance of 0.58. Thus, the Euclidean distance has a better clustering effect than the Bray-Curtis distance. The Euclidean distance was used in NMDS ordination [9] when it was more robust than Bray-Curtis.
The vegetation index (VI) is the result of an arithmetic operation between the pixel values of two or more spectral bands, which are correlated with plant parameters [39,40], and reveals the status of grassland degradation. The appropriate VIs of the QTP included NDVI, SAVI, and RDVI, which were selected by univariate linear regression analysis, with the larger R 2 values including the percent grass cover and the proportion of unpalatable grass as dependent variables [4]. NDVI is widely used for monitoring natural biological communities and agricultural ecosystems due to its easy accessibility and commonality. However, for sparse vegetation, the soil background contributes greatly to the total spectral reflectance [14], and it is not appropriate to evaluate grassland degradation using only NDVI. SAVI was proposed to compensate for the corresponding soil effect, reducing the significant contribution of soil reflectivity [41]. In high vegetation cover areas, NDVI is easily saturated, and RDVI can avoid the problems applied in the cases of high and low vegetation coverage, thus better adapting to areas with large disparities in coverage [42].

Quantitative Classifications of Alpine Grassland Degradation
The assessment of grassland degradation states is an important prerequisite for restoring and treating degraded grassland. Quantitative classifications of degradation can represent the degradation states more accurately. Alpine grasslands with different classes of degradation require different corresponding restoration and management measures [43], and the quantitative grading of degradation contributes to a more scientific alpine grassland management. The ecological thresholds of indicators in degradation assessment provide guidance for management thresholds and have theoretical and practical value for the scientific management of alpine grassland ecosystems [9]. Compared with field surveys, remote sensing is widely used in large-scale grassland degradation because of its higher efficiency [4,44,45], but quantitative classifications at a larger scale are rare in assessment studies.
In this study, the quantitative degradation classifications reflected the spatial pattern of grassland degradation on the basis of field data processed by multivariate statistical analysis.

Analysis of Different Classes of Grassland Degradation Combined with Management Measures
The distribution patterns of grassland degradation have spatial heterogeneity affected by different natural and human factors. The ND areas were mainly distributed in western Maizhokunggar County because this area is characterized by montane shrub grassland with a relatively low level of human activities such as grazing and farming. Most of the SLD areas were located in eastern Lhari County of the study area, and the vegetation type is composed mainly of wetland with a swamp soil. The MD areas were widely distributed in the study area, concentrated in Damxung and Lhari County. Grazing livestock is a heavy attribute for high forage quality with the vegetation types of alpine meadow, temperate grassland and alpine grassland and at higher altitudes. For SD areas, 27% was located in Damxung and 19% in Lingdrub. Sparse grassland was the main vegetation type, and the elevation was low in this area. With the addition of frequent cultivation activities in the surrounding areas, this area suffered from severe human disturbance. The ED areas were distributed sporadically in the study area, and Lhari County had the largest area. Vegetation coverage was low and bare land area was large, with an NDVI value less than 0.08. The soil types were swamp soil and frozen soil.
For different degradation states, specific measures should be implemented and appropriate restoration measures can improve the effectiveness of rangeland restoration [46]. The ranges of indicator thresholds provide a reference for the specific management threshold under different degradation levels. Once the dominance of Cyperaceae plants exceeds 28%, rangeland can be considered slightly degraded or moderately degraded. Enclosure and fencing are needed to maintain the state of rangeland to prevent it from being more seriously degraded. Rangeland can be regarded as extremely degraded when the dominance of Cyperaceae plants is under 28% and the soil total carbon content is under 0.88%. Certain management practices for rangeland restoration are required, including fencing protection, fertilization, and compensatory seed planting. Once the dominance of Cyperaceae plants is under 28%, soil total carbon content exceeds 0.88% and soil bulk density exceeds 0.86, rangeland can be considered severely degraded. Returning reclaimed land to grasslands should be implemented to reduce the human disturbance caused by frequent cultivated land activities [47].

Conclusions
In this study, a quantitative classification of alpine grassland degradation was obtained by a multivariate hierarchical analysis based on field data. Combined with remote sensing information, the scaling classification of degradation from the sample plots scale to the region scale and the spatial pattern of alpine grassland degradation was acquired for the study area. The indicator thresholds of degradation were determined, which contribute to the proceeding of alpine grassland restoration and management measures. The moderately degraded of alpine grassland is the main type (approximately 50%) among the five degradation classes of alpine grassland, and moderately degraded is widely distributed in the study area. When the dominance of Cyperaceae plants is under 28% and the soil total carbon content is under 0.88%, alpine grassland can be considered extremely degraded, and extremely degraded grasslands were found to represent 16% of all alpine grassland. This study not only quantitatively evaluated the degradation situation in this region, but also confirmed the main degradation factors and their corresponding thresholds, which has important significance for pasture production and grazing management.