Comprehensive Ecological Risk Changes and Their Relationship with Ecosystem Services of Alpine Grassland in Gannan Prefecture from 2000–2020

Alpine grassland is one of the most fragile and sensitive ecosystems, and it serves as a crucial ecological security barrier on the Tibetan Plateau. Due to the combined influence of climate change and human activities, the degradation of the alpine grassland in Gannan Prefecture has been increasing recent years, causing increases in ecological risk (ER) and leading to the grassland ecosystem facing unprecedented challenges. In this context, it is particularly crucial to construct a potential grassland damage index (PGDI) and assessment framework that can be used to effectively characterize the damage and risk to the alpine grassland ecosystem. This study comprehensively uses multi-source data to construct a PGDI based on the grassland resilience index, landscape ER index, and grass&ndash:livestock balance index. Thereafter, we proposed a feasible framework for assessing the comprehensive ER of alpine grassland and analyzed the responsive relationship between the comprehensive ER and comprehensive ecosystem services (ESs) of the grassland. There are four findings. The first is that the comprehensive ER of the alpine grassland in Gannan Prefecture from 2000&ndash:2020 had a low distribution in the southeast and a high distribution trend in the northwest, with medium risk (29.27%) and lower risk (27.62%) dominating. The high-risk area accounted for 4.58% and was mainly in Lintan County, the border between Diebu and Zhuoni Counties, the eastern part of Xiahe County, and the southwest part of Hezuo. Second, the comprehensive ESs showed a pattern of low distribution in the northwest and high distribution in the southeast. The low and lower services accounted for only 9.30% of the studied area and were mainly distributed in the west of Maqu County and central Lintan County. Third, the Moran&rsquo:s index values for comprehensive ESs and ER for 2000, 2005, 2010, 2015, and 2020 were &minus:0.246, &minus:0.429, &minus:0.348, &minus:0.320, and &minus:0.285, respectively, thereby indicating significant negative spatial autocorrelation for all aspects. Fourth, ER was caused by the combined action of multiple factors. There are significant differences in the driving factors that affect ER. Landscape index is the first dominant factor affecting ER, with q values greater than 0.25, followed by DEM and NDVI. In addition, the interaction between diversity index and NDVI had the greatest impact on ER. Overall, this study offers a new methodological framework for the quantification of comprehensive ER in alpine grasslands.