A Study on Spatiotemporal Changes of Ecological Vulnerability in Yunnan Province Based on Interpretation of Remote Sensing Images

: The inherent ecological environment of mountainous regions is highly fragile, and the degree of sustainable development is low. There has not yet been a multi-phase ecological vulnerability evaluation (EVE) study based on remote sensing (RS) and GIS for mountainous provinces, for which there is an urgent need to establish a system that is appropriate, practicable and easily operated and applied. In this study, an integrated “RS and GIS + multi-phase


Introduction
The study of vulnerability has become a focal topic and a meaningful analysis tool in the field of sustainable science and global environmental changes in recent years [1].Sustainable development goals (SDGs) mean common development among countries worldwide today, which includes both the sustainability of the ecological environment and the sustainability of economic and social development.Among them, the sustainability of the ecological environment is the foundation of the sustainability of economic and social development.In terms of the SDGs, vulnerability (especially the vulnerability of the ecological environment) is often an obstacle or in opposition to achieving sustainability (especially with regard to economic and social development).To enhance regional sustainable development capabilities, it is necessary to find ways to reduce vulnerability.The concept of "vulnerability" originated from research on natural calamities and impoverishment [2,3].In 1981, Timmerman expanded this study to the field of geography [4].It has then generally been applied in various disciplines: ecology, economics, sociology, etc.In particular, research on ecological environment vulnerability (EEV) and its quantitative evaluation has become a significant direction for research on the SDGs [5,6].Most of the existing empirical evaluation studies on ecological vulnerability focus on specific ecological areas, such as karst areas [7,8], northern agricultural and pastoral interlaced zones [9,10], coastal zones [11][12][13], river basins [14,15], reservoir areas [16] and high cold areas [17].In general, for mountainous areas, while some scholars have conducted empirical analyses [18], the depth and breadth of research in this area are still significantly insufficient.From a global perspective, mountainous areas account for approximately 30% of the earth; from China's perspective, it is a mountainous country, with mountainous areas accounting for over two-thirds [19].Mountainous areas are a unique natural-human synthesis with specific slopes and altitudes [20], and their extremely obvious feature is the EEV system.The innate fragility of mountainous systems leads to the instability of the ecosystems and the hardships in resource development and determines that once mountainous ecosystems are destroyed, they lack the ability to restore to their original state, thereby restricting the use of mountainous resources and the progress of mountainous economy and society.In terms of general mountainous regions, rural residents are constrained by poor natural resources and ecological, economic and social conditions, which might result in irrational human activities and rough development, further leading to the destruction of the ecological environment.Moreover, the continual occurrence of natural calamities and a decrease in income levels could increase the social-economic-ecological vulnerability.This has led to a vicious cycle of mountainous vulnerability → irrational exploitation → ecological damage → socioeconomic regression → mountainous vulnerability, causing the entire mountainous ecological, economic and social system to exhibit enormous fragility.For example, considering Yunnan province, located on the southwestern border of China, the mountainous area accounts for 94%.The ecological environment is highly vulnerable [21], and the economy remains correspondingly weak.Among the 129 counties (cities, districts) in Yunnan, 88 are classified as nationally impoverished, making Yunnan the representative mountainous province with the highest number of impoverished counties among all the provinces in China.Therefore, in-depth research on the ecological vulnerability of mountainous regions has important practical significance.
The core content of ecological vulnerability research is ecological vulnerability evaluation (EVE), which refers to the quantitative, qualitative or integrated study and identification of the degree of vulnerability in regional ecosystems or ecological environments.The aim of the study is to investigate the causes, mechanisms and change rules underlying ecosystem or regional ecological environment fragility, so that reasonable resource utilization methods and measures for ecological protection and restoration can be proposed toward ensuring the regional coordinated advancement of resources, environment, society and economy.In terms of the system of indicators and data sources of EVE, the progress and literature can be generally classified into four categories.The first is "existing statistical data and survey + quantitative measurement of EEV in the current situation or a certain year".For example, Liu Yi et al. (2010) conducted a natural disaster vulnerability evaluation of 31 provinces in China using DEA [22].Lu Yaling et al. (2010) analyzed the autocorrelation of EVE and spatial aspects for the Bohai Rim Region [23], and Liu Gang (2013) conducted a disaster risk vulnerability evaluation of communities [24].The second category is "the quantitative measurement of eco-environmental vulnerability for several years based on existing surveys and statistical data".For example, Li Yonghua et al. (2015) carried out a study on the spatiotemporal distinction of ecological vulnerability in Chaoyang County, Liaoning, China, from 2003 to 2012 [25].Nguyen (2016) conducted a study on the eco-environment vulnerability in Hue province, Vietnam, for the years 1989, 2003 and 2014 [26], and Liu Yujie et al. ( 2016) evaluated water resource vulnerability [27].The third category is "existing statistical data and survey + quantitative evaluation of vulnerability level of composite systems over several years".For example, Petrosillo et al. (2006) conducted a quantitative modeling study on the vulnerability of composite systems based on tourism activities [28].Li Bo (2012) conducted a spatiotemporal analysis of the vulnerability of the human-sea resource and environment system [29].Chen Jia et al. ( 2016) carried out a study on the vulnerability of the social ecosystem in Yulin City, Shaanxi, China, from 2000 to 2011 [30], and Wen Xiaojin et al. (2016) conducted a socioecological vulnerability evaluation of mountainous cities from 1997 to 2013 [31].The fourth category is "RS image interpretation + quantitative measurement of ecological vulnerability in a certain year or some years".For example, Zhang et al. (2007) studied the eco-environment vulnerability of Fuzhou [32].Xu Qingyong et al. ( 2013) comprehensively evaluated the eco-environment vulnerability of the Pearl River Delta from 2004 to 2008 [33].Hu Baoqing et al. (2014) comprehensively evaluated the eco-environment vulnerability of Guangxi karst areas based on GIS [8].Sahoo et al. (2016) evaluated the vulnerability of ecoenvironment [34].Wang Xuemei et al. (2016) conducted EVE of the Weigan River Basin [35].França et al. (2022) studied the ecological vulnerability of Brazil [36].From the perspective of the EVE index system, many researchers have conducted studies from various focuses and aspects, generally selecting various specific evaluation indicators from four framework systems and models: the Sensitivity-Resilience-Pressure (SRP) conceptual model [19], Natural Ecology-Human Socioeconomy [8], Pressure-State-Response (PSR) [37,38] and Vulnerability Scoping Diagram (VSD) [39,40].The core of the evaluation process is to determine the evaluation method.Currently, the primary approaches used include the comprehensive index method (CIM) [8,[41][42][43], the analytic hierarchy process (AHP) [33,44] and the cluster analysis method [8,33].Recently, EVE based on GIS and RS has generally been applied [8,[32][33][34][35][36]45,46].This is because, through RS image interpretation and GIS technology, information on LULC can be objectively obtained, which results in higher accuracy of the data on eco-environment vulnerability.In the existing research literature, there is currently no comprehensive EVE for certain consecutive years based on RS and GIS for mountainous provinces.Correspondingly, an appropriate, practicable, easily operated and applied regional EVE system for mountainous regions has not been established, which is not beneficial for carrying out the SDGs of mountainous areas.
Therefore, for this study, Yunnan province, the southwest frontier mountainous province, was selected as the empirical study region to conduct comprehensive EVE of mountainous provinces based on RS and GIS for numerous consecutive years.The reason why Yunnan province was chosen as the study area is not only because its mountainous area accounts for 94%, but also because it is located at the headstream or upstream region of six world-famous rivers, namely, the Yangtze River, Pearl River, Lancang River, Red River, Nu River and Irrawaddy River.Its location is extremely significant, and it is an ecological security barrier.Furthermore, Yunnan plays an indispensable and imperative role in ensuring the ecological security of these six major rivers and has significance as an international ecological security barrier.Strengthening ecological environment protection and striving to become the VECC are major strategic tasks of direct concern to both the whole country and Yunnan province.In January 2015, President Xi visited Yunnan and proposed the strategic position of striving to become the VECC for Yunnan.In January 2020, President Xi visited Yunnan again and demanded that Yunnan strive to continuously make new progress in building the VECC.For many years, Yunnan province has firmly established the concepts of "protecting the eco-environment is protecting productivity, improving the eco-environment is developing productivity", placing eco-civilization construction in a prominent position in the overall work and always adhering to the policy of "ecological priority and green development".As a result of work toward realizing the VECC in Yunnan province, there have been many unprecedented achievements [47].So, does Yunnan's VECC work alleviate and reduce the inherent ecological vulnerability and enhance the potential for sustainable development?What exactly is the effect?There is an urgent need to conduct accurate scientific evaluation.
Compared with existing study achievements, the primary characteristic and contribution of this paper consist in an approach, involving a comprehensive combination of "RS and GIS + LUCC with multi phases + practically quantitative theory and methods of EVE" based on the features of mountainous provinces and results of RS images interpretation of LULC from Yunnan province over five periods (i.e., 1980, 1990, 2000, 2010 and 2020).According to the research needs and data availability, this paper establishes a comprehensive evaluation indicator system, a grading system for various levels of DEV, and quantitative classification standards for EVE in the mountainous provinces.Furthermore, it promotes qualitative and quantitative integrated analysis methods for EVE in mountainous regions.At the same time, specific evaluation practices were carried out using county-level regions as evaluation units to quantitatively measure the ecological vulnerability level (EVL) in Yunnan province.The ecological vulnerability levels of 129 counties from 1980 to 2020 were determined, and based on these, the spatiotemporal evaluation of EVL for the whole province and 129 counties over the past 40 years was conducted, aiming to offer a reference for alleviating the DEV and boosting regional sustainable development strategic planning and management in mountainous provinces.In theory, this study not only promotes the progression and advancement of EVE research but also provides fundamental ideas and approaches for boosting regional actions toward achieving the SDGs, which is theoretically helpful for further enrichment and innovative development of the ecological vulnerability system.In actuality, methods are also established for qualitative and quantitative integrated analysis of DEV and its dynamic change trend in mountainous areas.

Overview of the Study Area
Yunnan is located in the southwestern border area of China (21 • 8 32 N-29 • 15 8 N and 97 • 31 39 E-106 • 11 47 E).It is one of the provinces with the most neighboring countries and the longest border in China.So far, it has 16 prefecture-level cities under its jurisdiction, and there are 129 county-level administrative units [48] (Figure 1).
Yunnan is a representative mountainous province mainly composed of mountains amounting to 94% being mountainous area [49].The terrain of the whole province is characterized by being high in the northwest and low in the southeast (Figure 1).The average elevation is about 2000 m, and the difference in elevation between the highest point and the lowest point is 6663.6 m.The basic characteristics of Yunnan's landform are mainly mountainous areas and large terrain slopes.In Yunnan province, more than three-quarters of the land has a ≥15 • slope, and nearly two-thirds of the land has a ≥25 • slope [50], indicative of the steep terrain of Yunnan.The terrain characteristics dominated by mountainous areas cause the "innate" ecological environment to be exceedingly vulnerable, and cultivated land resources remain restricted.The long-term unreasonable exploitation and use of mountainous land resources can easily result in the incompatibility of humans and land in the mountainous regions, causing soil erosion and environment deterioration and further leading to serious disasters.Currently, major eco-environmental issues such as rocky desertification in the province are still severe, with a trend of frequent and prone occurrences of geological disasters, such as drought, floods, wind and hail, low-temperature freezing, landslides and mudslides [51].
occurrences of geological disasters, such as drought, floods, wind and hail, low-temperature freezing, landslides and mudslides [51].The special geographical location and natural environment make Yunnan province a barrier to ecological security.The country is focused on preserving the ecological environment in Yunnan, and there are positive expectations that Yunnan will strive to become the VECC.This is not only a major strategic task undertaken by Yunnan but also represents the foundation for sustainable development in the province.It is necessary to implement an ecologically friendly development approach according to local conditions and take a win-win path toward protection and development.

RS Data Acquisition and Interpretation
Data for the five phases of the RS image (i.e., 1980, 1990, 2000, 2010 and 2020) used in this paper were obtained from the website of the Resource and Environmental Science and Data Center of the Chinese Academy of Sciences (CAS) (URL: "https://www.resdc.cn/"and accessed on 23 August 2022).Images with winter cloud cover of less than 10% were selected for interpretation (Table 1).The entire interpretation work is based on a unified land use classification system and it was combined with RS interpretation symbols in the ArcGIS software environment to interpret the LULC types from 5 phases of RS images through human-machine interaction, thus obtaining a vector database for the 5 phases of LULC (Figure 2).The special geographical location and natural environment make Yunnan province a barrier to ecological security.The country is focused on preserving the ecological environment in Yunnan, and there are positive expectations that Yunnan will strive to become the VECC.This is not only a major strategic task undertaken by Yunnan but also represents the foundation for sustainable development in the province.It is necessary to implement an ecologically friendly development approach according to local conditions and take a win-win path toward protection and development.

RS Data Acquisition and Interpretation
Data for the five phases of the RS image (i.e., 1980, 1990, 2000, 2010 and 2020) used in this paper were obtained from the website of the Resource and Environmental Science and Data Center of the Chinese Academy of Sciences (CAS) (URL: "https://www.resdc.cn/"and accessed on 23 August 2022).Images with winter cloud cover of less than 10% were selected for interpretation (Table 1).The entire interpretation work is based on a unified land use classification system and it was combined with RS interpretation symbols in the ArcGIS software environment to interpret the LULC types from 5 phases of RS images through human-machine interaction, thus obtaining a vector database for the 5 phases of LULC (Figure 2).Referring to Xu Xinliang et al. [52] and Liu Jiyuan et al. [53][54][55][56], and considering the real conditions and research needs of Yunnan, five-phase LULC classification systems in Yunnan were determined to comprise 6 first-level land use types and 12 second-level land use types.Figure 2 provides a detailed description of the entire process of RS image interpretation to acquire the five-phase LULC map of Yunnan.Finally, a vector database of LULC was generated, and ArcGIS was used to compile the five-phase LULC map (Figure 3), and values for the classifications of LULC area were calculated (Table 2).Referring to Xu Xinliang et al. [52] and Liu Jiyuan et al. [53][54][55][56], and considering the real conditions and research needs of Yunnan, five-phase LULC classification systems in Yunnan were determined to comprise 6 first-level land use types and 12 second-level land use types.Figure 2 provides a detailed description of the entire process of RS image interpretation to acquire the five-phase LULC map of Yunnan.Finally, a vector database of LULC was generated, and ArcGIS was used to compile the five-phase LULC map (Figure 3), and values for the classifications of LULC area were calculated (Table 2).This study is based on the awareness that regional ecological vulnerability is the result of "innate" natural ecological vulnerability (INEV) and the land use ecological vulnerability (LUEV) and land cover ecological vulnerability (LCEV) formed by human acquired socioeconomic activities.Consequently, the EVE of research is mostly a comprehensive assessment that takes county-level administrative regions as the assessment unit and, based on RS image interpretation of LULC information, organically integrates the three dimensions of INEV, LUEV and LCEV in the counties, evaluating the comprehensive conditions of ecological vulnerability within the counties.The quantitative multi-indicator comprehensive evaluation method and qualitative comprehensive analysis method are organically combined, and the degree of INEV, LUEV and LCEV in each county over 5 periods is computed.Furthermore, these three indicators are seamlessly combined to compute the level of overall ecological vulnerability (OEV) across all counties over the past 40 years in order to acquire a comprehensive assessment and understanding of 129 counties and comprehensively measure and compare the situations of ecological vulnerability within these counties.
(2) Evaluation Indicator System The evaluation indicator system to be built is primarily aimed at regional EVE.Referring to Zhang Fengrong et al. [57] (2003), the indicator system is classified into three levels (Table 3).They are generally determined based on the demands of assessment work and the accessibility of existing data from the 129 counties.
× 100 Socioeconomic Statistical Yearbook Considering the regional background, this paper takes the G YPUA of the county with the highest grain yield in the province as the relative optimal value.Sowing Area of Grain Crops (S AGC ) Soil Erosion Area (S EA ) I SEA = P SEA /Maximum P SEA × 100

Existing Thematic Surveys
The larger the proportion of soil erosion area, the higher the vulnerability.This paper takes the proportion of soil erosion area in the county with the largest relative area in the province as the relative extreme value.Total Land Area (T LA ) P SEA (Proportion of Soil Erosion Area) = S EA /T LA × 100% Land Area A bio = 511.2642[58][59][60] Calculate According to the Interpretation Results of RS Images It depends on the background of the region.According to the counties in Yunnan province with the best ecological protection, the relative optimal value was determined.

I ESV
Value of Ecological Service per Unit Land Area (V ES ) × 100 Where, V ES is calculated according to Xie Gaodi et al. [61,62] Calculate According to the Interpretation Results of RS Images Considering the regional background, the counties with the best V ES were used to determine the relative optimal value.
The INEV index is a fundamental indicator of the ecological vulnerability of the entire system, directly determining whether the county can sustain its inherent role and function in terms of ecology.Considering the actual conditions in various regions and the accessibility of existing data, the main indicators selected for assessing the natural ecological vulnerability status are the index of mountain area (I MA ), the index of steep-slope area with slope ≥ 25 • (I SSA ), the index of high-altitude area (I HAA ) and the index of average annual rainfall (I AAR ).
The ecological vulnerability of land use mainly reflects the excessive development and use status of land as well as the level of resource transformation.For the former, the main indicators selected are the index of over-reclaimed rate (I OR ) and the index of bare land area (I BLA ).For the latter, the main indicators selected are the index of effective irrigated area of cultivated land (I EI ) and the index of grain yield per unit area (I GY ).
The ecological vulnerability of land cover mainly reflects the level of environmental protection and ecological value in various regions.For this study, the main indicators selected are the index of forest coverage rate (I FC ), the index of soil erosion area (I SEA ), the index of biological richness conversion (I BRC ) and the index of ecosystem services value (I ESV ).Of these, I BRC is a significant index for evaluating eco-environment conditions, indirectly reflecting the richness of biological abundance in the evaluated areas [58][59][60].It is positively correlated with the ecological benefits of resource utilization; the higher the I BRC , the lower the ecological vulnerability in the region.I ESV refers to the various benefits that people directly or indirectly obtain from the ecosystem.Sustainable utilization of resources and optimal ecological value are the fundamental aims of regional development strategies.The higher the I ESV , the lower the DEV.
It might be pointed out that for the purpose of convenient contrast and analysis, a unified requirement has been proposed that all values of indicators could be converted into a value from 0 to 100, where values close to 100 indicate the highest vulnerability but lowest sustainability; values close to 0 indicate the lowest vulnerability but highest sustainability.Therefore, a process for transforming the index is also needed (Table 3).
(3) Assessment Criteria for EVE Indicators The assessment criteria are also known as the threshold.Developing assessment criteria for each indicator is a very complex task; as a result, no uniform assessment criterion for EVE has been established so far, representing an issue that needs to be constantly explored [63].Hence, the data and values of assessment criteria may be obtained through various approaches: First, some indicators (such as I SEA ) that have already been stipulated by the nation can be evaluated using standards formulated by the nation, local or industry.Second, for some restrictive indicators such as I OR , it is necessary to use the bottom line or warning value acquired from the scientific study as the assessment criteria and reference.Third, for some indicators, their average value for the country (or the research area) or the target value of regional planning can often be used as the assessment criterion.In addition, there are also some indicators that need to be comprehensively analyzed according to specific regional conditions, such as I FC and I BRC , etc.In this study, the assessment criteria of the indicators of ecological vulnerability in Yunnan's counties were determined (Table 3).

Comprehensive Method for EVE
To quantitatively evaluate regional ecological vulnerability as a whole, four indicators of QMCE are proposed, namely, D INEV , D EVLU , D EVLC and Degrees of Overall Ecological Vulnerability (D OEV ), which reflect the degree of INEV, LUEV, LCEV and overall ecological vulnerability (OEV), respectively.The "multi-index comprehensive evaluation method" [64] was used to determine the calculation method of these four comprehensive indexes.
In Equation ( 1), w 11 , . .., w 14 are the weight values of the I MA , I SSA , I HAA and I AAR , respectively.The higher the D NEV value, the more pronounced the degree of INEV and the lower the degree of sustainability.
The four indicators that characterize D EVLU are I OR , I BLA , I EI and I GY .On the basis of separately calculating these indicator values, the D EVLU value can be quantitatively computed using the following equation: In Equation ( 2), w 21 , . .., w 24 are the weight values of I OR , I BLA , I EI and I GY , respectively.The higher the D EVLU value, the more pronounced the degree of LUEV and the lower the degree of sustainability.
There are also four indicators that characterize D EVLC , namely, I FC , I SEA , I BRC and I ESV .On the basis of separately calculating these index values, the D EVLC value may be quantitatively computed by using the following equation: In Equation (3), w 31 , . .., w 34 are the weight values of the I FC , I SEA , I BRC and I ESV , respectively.The higher the D EVLC value, the more pronounced the LCEV and the lower the degree of sustainability.
The OEV is the organic combination of values of the aforementioned D INEV , D EVLU and D EVLC .On the basis of separately calculating the values of these three comprehensive indicators, the D OEV value of the study area can be quantitatively calculated as follows: In Equation ( 4), w 1 , w 2 and w 3 are the weight values of D INEV , D EVLU and D EVLC , respectively.The higher the D OEV value, the more pronounced the OEV level and the lower the sustainability level. (

5) Method for determining indicator weights and resulting values
The methods for determination of weight mainly include Principal Component Analysis (PCA), AHP and Delphi Method (or Expert Consultation Method), of which the latter approach is more commonly used to determine the weight coefficient.This study uses the frequently used Delphi method (expert consultation method) to determine the weight coefficients of each item.The specific approach is as follows: In November 2022, we organized 16 experts to assign weights to the various indicators of ecological vulnerability assessment in the above-mentioned regions.After providing feedback on the probability estimation results, the experts scored the weight of each indicator in the second and third rounds, gradually narrowing the scattered weights.Finally, we obtained a more coordinated and consistent weight value for each factor.After corresponding statistical processing, we obtained the weight value of each level of indicator.After the corresponding processing, the weight values of various levels of indicators were obtained (Table 4).After calculating the DEV of the study region, it is also necessary to use this as a basis to classify the regional ecological vulnerability, aiming toward qualitatively assessing the DEV at different levels and combining qualitative and quantitative research results to better provide a scientific reference and guidance for implementing regional measures to meet the SDGs.
For a long time, there has been no specialized exploration or research on the regional ecological vulnerability classification system and standards.Here, based on our awareness and experience gained through serial surveys in various regions, combined with real conditions of the background and with reference to the rating system built by Yang Renyi et al. [21] (2021) and Yang Zisheng et al. [60] (2007), the regional DEV is divided into five levels (Table 5).In addition, the grading standard and the main meaning of each grade are shown (Table 5).
Table 5. Grading standard and main meaning of regional ecological vulnerability.

Grades of Ecological
Vulnerability D OEV Meaning

Very Slightly Vulnerable <35
The DEV of regional development is very low; regional resource development and utilization activities have not caused significant impact or damage to the ecological environment; it can ensure the ecological sustainability of regional development.

Lowly Vulnerable 35~45
The DEV in regional development is not high; regional resource development and utilization activities have caused a certain degree of impact and damage; by adopting general measures, the ecological sustainability of regional development can be ensured.

Moderately Vulnerable 45~55
The DEV in regional development is relatively high; regional resource development and utilization activities have caused significant impact and damage; effective measures need to be taken to ensure the ecological sustainability of regional development.

Highly Vulnerable 55~65
The DEV in regional development is high; regional resource development and utilization activities have caused great impact and damage; strong measures need to be taken to ensure the ecological sustainability of regional development.

Very Highly Vulnerable ≥65
The DEV in regional development is very high, and the degradation and deterioration of the ecological environment are particularly prominent; regional resource development and utilization activities have caused tremendous impact and damage.It is necessary to fundamentally reverse the ways of resource utilization and regional development and take significant measures to ensure the ecological sustainability of regional development.

Results and Analysis
According to the above comprehensive EVE method, D INEV , D EVLU , D EVL and D OEV in 1980, 1990, 2000, 2010 and 2020 were quantitatively calculated.The county-level evaluation map of INEV (Figure 4), LUEV (Figure 5), LCEV (Figure 6) and OEV (Figure 7) in Yunnan province from 1980 to 2020 was compiled.Based on the calculation results, the spatiotemporal changes in INEV, LUEV, LCEV and OEV of Yunnan province can be analyzed.

The Spatiotemporal Evolution Characteristics of INEV and Its Spatial Difference
The D INEV of Yunnan province is 58.21, which corresponds to the "highly vulnerable" level.It means that for Yunnan province as a whole, the D INEV is high.This is the first characterization of the INEV in Yunnan province.This feature is mainly caused by the high and steep mountain terrain characteristics.
The second remarkable characteristic is that the spatial difference in D INEV value is large in the province and the "innate" natural ecological environment of the high mountain and canyon area in northwest Yunnan is the most fragile.In the region, Nujiang prefecture, Diqing prefecture and Lijiang city have D INEV values of 81.02, 84.51 and 69.79, respectively, which correspond to the "very highly vulnerable" level, followed by the northeast mountainous area of Yunnan, where most counties have D INEV values of more than 55, with the highest being over 75 (i.e., Dongchuan, Qiaojia), and most counties are classified as "highly vulnerable" and "very highly vulnerable".In addition, the D INEV values are also high in central and southeastern Yunnan, with most counties reaching over 55, and the region of southern Yunnan and southwestern Yunnan has a relatively low D INEV , for example, with a D INEV of 42.23 in Dehong prefecture.Considering the overall situation of the 129 counties in the province, most have a high D INEV .Table 6 and Figure 4 show that 16.28% of the counties were classified as "very highly vulnerable", 36.43% as "highly vulnerable", 26.36% as "moderately vulnerable" and 13.18% as "low vulnerable", and "very slightly vulnerable" counties accounted for only 7.75%.In other words, 79.07% of counties were classified at the "moderately vulnerable" level or above.The reason why the INEV of Yunnan province is of great spatial difference is mainly due to the significant spatial difference and the superposition effect of the four indexes I MA , I SSA , I HAA and I AAR .classified as "highly vulnerable" and "very highly vulnerable".In addition, the DINEV values are also high in central and southeastern Yunnan, with most counties reaching over 55, and the region of southern Yunnan and southwestern Yunnan has a relatively low DINEV, for example, with a DINEV of 42.23 in Dehong prefecture.Considering the overall situation of the 129 counties in the province, most have a high DINEV.Table 6 and Figure 4 show that 16.28% of the counties were classified as "very highly vulnerable", 36.43% as "highly vulnerable", 26.36% as "moderately vulnerable" and 13.18% as "low vulnerable", and "very slightly vulnerable" counties accounted for only 7.75%.In other words, 79.07% of counties were classified at the "moderately vulnerable" level or above.The reason why the INEV of Yunnan province is of great spatial difference is mainly due to the significant spatial difference and the superposition effect of the four indexes IMA, ISSA, IHAA and IAAR.  Figure 5 shows the spatiotemporal changes in DEVLU, and the main characteristics of its grades for the whole province and 129 counties from 1980 to 2020 will now be dis-

The Spatiotemporal Evolution Characteristics of LUEV
Figure 5 shows the spatiotemporal changes in D EVLU , and the main characteristics of its grades for the whole province and 129 counties from 1980 to 2020 will now be discussed.First, the DEVLU values of Yunnan have gradually dropped from 1980 to 2020, and the corresponding level of LUEV was reduced from the "moderately vulnerable" to the "lowly vulnerable" level.The average value of DEVLU was 53.50 in 1980, and it dropped to 43.94 in 2020, with a net decrease of 9.56.It means that in the past 40 years, significant progress has been achieved as a result of the measures to protect and construct the ecological environment in land use in Yunnan province.In particular, due to the implementation of crucial ecological preservation measures in the past decade, such as comprehensive management of sloping cultivated land, the average annual decrease in DEVLU of Yunnan province from 2010 to 2020 was 0.76%, and the grade of LUEV improved to the level of "lowly vulnerable".
Second, from 1980 to 2020, the LUEV in counties has changed significantly, and the general rule is that the number of counties with high DEVLU significantly decreased while the number of counties with low DEVLU significantly increased.Figure 5 and Table 7 show that over the 40-year period from 1980 to 2020, the number of "very highly vulnerable" counties decreased from 14 to 1, with a net decrease of 92.86%; the number of "highly vulnerable" counties decreased from 41 to 11; the number of "moderately vulnerable" counties decreased from 53 to 41, with a net decrease of 22.64%; the number of "lowly vulnerable" counties increased from 18 counties to 49 counties, with a net increase of 1.72 times; and the number of "very slightly vulnerable" counties increased from 3 to 27, with a net increase of 8 times.This indicates that the many ecological protection measures that were taken in the past 40 years have been effective, with obvious effects observed.First, the D EVLU values of Yunnan have gradually dropped from 1980 to 2020, and the corresponding level of LUEV was reduced from the "moderately vulnerable" to the "lowly vulnerable" level.The average value of D EVLU was 53.50 in 1980, and it dropped to 43.94 in 2020, with a net decrease of 9.56.It means that in the past 40 years, significant progress has been achieved as a result of the measures to protect and construct the ecological environment in land use in Yunnan province.In particular, due to the implementation of crucial ecological preservation measures in the past decade, such as comprehensive management of sloping cultivated land, the average annual decrease in D EVLU of Yunnan province from 2010 to 2020 was 0.76%, and the grade of LUEV improved to the level of "lowly vulnerable".
Second, from 1980 to 2020, the LUEV in counties has changed significantly, and the general rule is that the number of counties with high D EVLU significantly decreased while the number of counties with low D EVLU significantly increased.Figure 5 and Table 7 show that over the 40-year period from 1980 to 2020, the number of "very highly vulnerable" counties decreased from 14 to 1, with a net decrease of 92.86%; the number of "highly vulnerable" counties decreased from 41 to 11; the number of "moderately vulnerable" counties decreased from 53 to 41, with a net decrease of 22.64%; the number of "lowly vulnerable" counties increased from 18 counties to 49 counties, with a net increase of 1.72 times; and the number of "very slightly vulnerable" counties increased from 3 to 27, with a net increase of 8 times.This indicates that the many ecological protection measures that were taken in the past 40 years have been effective, with obvious effects observed.Third, the regional difference in LUEV is large.Although the values of D EVLU have decreased to different degrees from 1980 to 2020, the annual average decreases are quite varied.Overall, the annual average decrease in D EVLU value of central, southern and southwestern Yunnan is relatively large.The annual average decrease in D EVLU value for most counties in these areas is more than 0.60, and the decrease in values reached more than 1.00 in Chenggong, Hongta, Guandu, Xishan, Qilin, Wuhua and Anning, while smaller annual decreases were observed in the northwest, northeast and southeast Yunnan, where the average annual decrease in D EVLU was less than 0.40 in most of these counties, with Yanjin, Daguan, Fugong, Ninglang, Shangri-la, Lanping and Lushui dropping below 0.30.This is mainly related to the condition of land over-exploitation (mainly land overcultivation and bare land) and the level of resource transformation (mainly E IR and G YPUA ).

The Spatiotemporal Evolution Characteristics of LCEV
Figure 6 shows that from 1980 to 2020, the spatiotemporal evaluation of D EVLC and its grades in the whole province and all counties also exhibit the following three remarkable characteristics.
First, the D EVLC values in Yunnan have obviously dropped from 1980 to 2020.The corresponding level of vulnerability dropped from the "highly vulnerable" to the "very slightly vulnerable" level.In 1980, the average D EVLC value of the whole province was 55.79.By 2020, the D EVLC value decreased to 30.01, the net decrease value reached 25.78 and the annual average decrease was 1.16%, which is equivalent to 2.6 times the average annual decrease in D EVLU .This means that over the past 40 years, Yunnan province has achieved remarkable results in improving land cover (especially forest cover rate) and ecological service value while controlling soil and water loss.In particular, 2000 marked the beginning of measures for the conversion of cultivated land to woodland and grassland (CCWG), rocky desertification land consolidation, afforestation and greening of barren mountains and other major ecological construction projects with comprehensive effect.From 2000 to 2020, the average annual decrease in D EVLC of Yunnan province was 1.70%, with the reduction of LCEV to the "very slightly vulnerable" level.Second, the LCEV in counties has changed greatly from 1980 to 2020.The basic rule is that the number of counties with high DEVLC greatly decreased while the number of counties with low DEVLC greatly increased.Figure 6 and Table 8 show that over the 40-year period from 1980 to 2020, the number of "very highly vulnerable" counties decreased from 48 to 2, with a net decrease of 95.83%; the number of "highly vulnerable" counties decreased from 33 to 8, with a net decrease of 75.76%; the number of "moderately vulnerable" counties decreased from 29 to 24, with a net decrease of 17.24%; the number of "lowly vulnerable" counties increased from 15 to 29, with a net increase of 93.33%; and the number of "very slightly vulnerable" counties increased from 4 to 66, with a net increase of 15.50 times.This shows that in the past 40 years, the counties have taken numerous effective measures for ecological construction, such as CCWG, rocky desertification land consolidation, afforestation and greening of barren mountains.Third, the change in the LCEV presents remarkable regional variations in geographical space.In the past 40 years, the values of DEVLC have decreased to differing degrees, while the regional difference in annual average reduction is very significant.In general, Second, the LCEV in counties has changed greatly from 1980 to 2020.The basic rule is that the number of counties with high D EVLC greatly decreased while the number of counties with low D EVLC greatly increased.Figure 6 and Table 8 show that over the 40-year period from 1980 to 2020, the number of "very highly vulnerable" counties decreased from 48 to 2, with a net decrease of 95.83%; the number of "highly vulnerable" counties decreased from 33 to 8, with a net decrease of 75.76%; the number of "moderately vulnerable" counties decreased from 29 to 24, with a net decrease of 17.24%; the number of "lowly vulnerable" counties increased from 15 to 29, with a net increase of 93.33%; and the number of "very slightly vulnerable" counties increased from 4 to 66, with a net increase of 15.50 times.This shows that in the past 40 years, the counties have taken numerous effective measures for ecological construction, such as CCWG, rocky desertification land consolidation, afforestation and greening of barren mountains.Third, the change in the LCEV presents remarkable regional variations in geographical space.In the past 40 years, the values of D EVLC have decreased to differing degrees, while the regional difference in annual average reduction is very significant.In general, the annual average decrease in D EVLC values of northwest, southwest, south and central Yunnan is relatively large, and the annual average decrease in D EVLC values in most counties in these regions is more than 1.20, reaching 1.50 or more in Jiangcheng, Ning'er, Jingdong, Simao, Shangri-la, Zhenyuan, Jinggu, Yingjiang, Dayao, Ninglang, Jinghong, Fugong, Mojiang, Yunlong, Weixi, Mengla, Yongping, Nanhua, Tengchong, Yulong, Lufeng, Lushui, Gongshan, Deqin, Yongren, Zhenkang and Chuxiong.The annual average decrease in D EVLC was less than 0.90 in most counties in northeastern, eastern and southeastern Yunnan, where it fell below 0.60 in Luxi, Mi'le, Shilin, Fuyuan, Luoping, Zhenxiong, Zhaoyang, Kaiyuan, Dongchuan and Qiubei.This is mainly related to the state of forest cover rate, soil erosion, biological abundance and ecological service value in different periods.the annual average decrease in DEVLC values of northwest, southwest, south and central Yunnan is relatively large, and the annual average decrease in DEVLC values in most counties in these regions is more than 1.20, reaching 1.50 or more in Jiangcheng, Ning er, Jingdong, Simao, Shangri-la, Zhenyuan, Jinggu, Yingjiang, Dayao, Ninglang, Jinghong, Fugong, Mojiang, Yunlong, Weixi, Mengla, Yongping, Nanhua, Tengchong, Yulong, Lufeng, Lushui, Gongshan, Deqin, Yongren, Zhenkang and Chuxiong.The annual average decrease in DEVLC was less than 0.90 in most counties in northeastern, eastern and southeastern Yunnan, where it fell below 0.60 in Luxi, Mi le, Shilin, Fuyuan, Luoping, Zhenxiong, Zhaoyang, Kaiyuan, Dongchuan and Qiubei.This is mainly related to the state of forest cover rate, soil erosion, biological abundance and ecological service value in different periods.First, the DOEV in Yunnan obviously decreased, with the whole OEV grade changing from the "highly vulnerable" to the "lowly vulnerable" level.The average DOEV of Yunnan province in 1980 was 55.88 and decreased to 44.47 in 2020, with a net decrease of 11.41.The OEV level of the whole correspondingly decreased from the "highly vulnerable" to the "lowly vulnerable" level.This means that the measures have achieved significant results in the past 40 years, which is also an important achievement following the vigorous implementation of the VECC strategy in Yunnan province.

The Spatiotemporal Evolution Characteristics of OEV
Second, the regional differences in OEV level and DOEV changes are large.From 1980 to 2020, in terms of the changes in the OEV level of 16 prefectures, Diqing and Zhaotong were both downgraded from the "very highly vulnerable" to the "moderately vulnerable" First, the D OEV in Yunnan obviously decreased, with the whole OEV grade changing from the "highly vulnerable" to the "lowly vulnerable" level.The average D OEV of Yunnan province in 1980 was 55.88 and decreased to 44.47 in 2020, with a net decrease of 11.41.The OEV level of the whole correspondingly decreased from the "highly vulnerable" to the "lowly vulnerable" level.This means that the measures have achieved significant results in the past 40 years, which is also an important achievement following the vigorous implementation of the VECC strategy in Yunnan province.
Second, the regional differences in OEV level and D OEV changes are large.From 1980 to 2020, in terms of the changes in the OEV level of 16 prefectures, Diqing and Zhaotong were both downgraded from the "very highly vulnerable" to the "moderately vulnerable" level.Nujiang, Lijiang, Qujing, Chuxiong, Lincang, Honghe and Wenshan were all downgraded from the "highly vulnerable" to the "moderately vulnerable" level.Kunming, Yuxi, Dali, Baoshan and Pu'er were downgraded from the "moderately vulnerable" to the "lowly vulnerable" level.Xishuangbanna was downgraded from the "moderately vulnerable" to the "very slightly vulnerable" level, and Dehong was downgraded from the "lowly vulnerable" to the "very slightly vulnerable" level.
It is worth pointing out that the actual D OEV values in the two counties that maintained at "very highly vulnerable" and "lowly vulnerable" levels without any level changes have shown a slow downward trend.
The results of the above-mentioned dynamic evolution caused the number of counties with different OEV levels to correspondingly change from 1980 to 2020.The basic rule is that the number of counties with high D OEV significantly decreased while the number of counties with low D OEV significantly increased (Table 9).From 1980 to 2020, the number of "very highly vulnerable" counties decreased from 14 to 1, with a net decrease of 92.86%; the number of "highly vulnerable" counties decreased from 61 to 8, with a net decrease of 86.89%; the number of "moderately vulnerable" counties increased from 44 to 50, with a net growth of 13.64% and the number of "lowly vulnerable" counties increased from 10 to 55, with a net increase of 4.50 times.In 1980, there were no "very slightly vulnerable" counties; by 2020, there were 15 "very slightly vulnerable" counties, accounting for 11.63%.

Spatial Differences in OEV of Yunnan Province
On the basis of the calculation results for 2020, the average D OEV , on the whole, has dropped to 44.47, indicating that the OEV level of the province has decreased to a "lowly vulnerable" level.However, there are significant regional differences in the D OEV value within the province.In general, the D OEV value is lower in the southern, southwestern, western and central regions of Yunnan and higher in the northwest high mountain and canyon areas, the northeastern mountainous areas and some karst areas.From the countylevel perspective, the OEV level is still relatively high, with 50 counties (38.76%) being "moderately vulnerable".Deqin, Daguan, Huize, Zhenxiong, Qiaojia, Yongshan, Xichou and Malipo are classified as "highly vulnerable".Dongchuan, known as the "mudslide museum", remains at the "very highly vulnerable" level.

Analysis of the Reasons for the Spatiotemporal Evolution of OEV
The main reason why the D OEV in Yunnan has gradually dropped over the past 40 years is that Yunnan province has paid great attention to ecological construction and environmental protection and established a strategy of "constructing ecological province and prioritizing the environment".The tremendous efforts in protection and governance have obviously boosted the continual upgrading of the ecological environment in the province, evidently increased F CR and gradually enhanced E SV .From a provincial perspective, over the past 40 years, except for a slight decrease in the E IR due to urbanization, industrialization and various infrastructure constructions occupying paddy fields, the other ecological vulnerability indicators of LULC have shown significant improvement (  However, considering all things, the D OEV in Yunnan remains relatively high, which is an overall reflection of the "innate" situations of the natural environment, such as terrain, excessive long-term land development and resource transformation (consolidation) level.In terms of "innate" terrain, "there are many steep slopes and high mountains with lots of rocks, climbing when going out" is a basic reflection of the terrain characteristics of most areas in Yunnan province.About 77% of the land area has a >15 • slope, and nearly 40% is steep slopes with a >25 • slope.Over 56% of the counties have a >30% land ratio of steep slope, and over one-fifth of counties have a >50% land ratio of steep slope [50].The result of long-term excessive development and use of land is particularly reflected in overreclamation (including deforestation and cultivating on steep slopes) and the abundance of mountainous regions maintaining primitive, extensive and backward practices, accounting for the considerable proportion of slope-cultivated land [49].According to the main data bulletin [65], slope cultivated land of >15 • and ≤25 • represents 27.22% and steep slope cultivated land of >25 • (including terraced fields above 25 • ) 18.64%.The area of suitable cultivated land (including suitable cultivated land in existing cultivated land and unused land) was 4.7076 million hectares in 2020 in the whole province, representing a suitable cultivation rate of just 12.25%.However, the actual reclamation rate in 2020 reached 14.04%.Each county has varying degrees of excessive development and use.Overall, 68.22% of counties have an O RR of over 10%, and around 25% of counties have an O RR of over 20%.This is a core reason behind the generally weak ecological environment for Yunnan's socioeconomic development.Furthermore, the measures for resource transformation (consolidation) and ecological environment preservation have not maintained.For instance, the E IR value is not large.The value of E IA in the whole province in 2020 was less than 1.40 million hectares, and the value of E IR was only about 25%.In addition, due to the long-term impact of different intensities and methods of land development and utilization, most counties have various distributions of B LAR .In 2020, the proportion of B LA in the province was 2.10%.Luoping, Guangnan, Xichou, Qiubei and Yanshan have more than 10% of their land classified as bare.
The obvious areal differentiations of the current OEV level within the province are related to the topographic pattern, LULC situation and ecological environment construction status in the various regions.Specifically, there are obvious areal distinctions in the 3 aspects and 12 indicators that reflect the regional OEV level mentioned above.For instance, in terms of landform pattern, the terrain of the high mountain and canyon area in the northwest and the middle high mountain area in the northeast is steep.The ratio of steep land with slopes greater than 25 • in most counties is over 50%; however, in Gongsha, Fugong and Deqin, this ratio reaches 91.18%, 88.31% and 80.88%, respectively.The ratio of steep land with slopes greater than 25 • in most counties of central and eastern Yunnan is less than 30%.From the perspective of F CR in various counties, there are 50 counties with an F CR < 40%, whereas 9 counties have an F CR < 20%.In terms of the distribution of B LAR , the values of B LAR in Wenshan and Diqing are about 8%.From the perspective of county-level units, there are 55 counties in the province with a B LAR of over 1%, among which 6 counties have a B LAR of 5% to 10%, and 5 counties have a B LAR of over 10%.From the perspective of E IR , over 95% of counties have an E IR below 50%, with about two-thirds of counties having an E IR below 30%, and 25 counties having an E IR below 10%.In addition, there are significant differences in I BR and V ES .About two-fifths of counties have I BR < 120, and 13 counties have I BR < 100.About two-fifths of counties have an average V ES of <1.2 million CNY/km 2 and 20 counties have an average V ES value of <1.0 million CNY/km 2 .

Discussion
In this study, based on the five phases of LULC data obtained from RS image interpretation in Yunnan province, the three dimensions of ecological vulnerability of the INEV, and the "acquired" LUEV and LCEV formed by human social and economic activities, were organically integrated, and comprehensive evaluation of the OEV was carried out.Overall, from 1980 to 2020, Yunnan province gradually paid increasing attention to ecological environment preservation.Among some indicators related to LUEV and LCEV, except for a slight decrease in the E IR due to urbanization, industrialization and various infrastructure constructions occupying paddy fields, the other indicators reflect significant improvements.The O RR , the B LAR , F CR , I BR and V ES in the vast majority of counties have all significantly improved, which is the basic reason why the D OEV of Yunnan province and most counties has gradually decreased over the past 40 years.This also fully reflects the significant achievements made in ecological environment preservation with advancements from 1980 to 2020, especially the important achievements in Yunnan in vigorously implementing the strategy of VECC.
However, because of the inherent terrain and other natural environmental conditions, as well as the comprehensive impact of long-term unreasonable development, the OEV level on the whole is still relatively high.The factor of terrain is the most important in terms of "innate" natural circumstances.Yunnan has 77% land area with a >15 • slope, and about two-fifths of the land area consists of steep slopes with more than 25 • .Over 56% of counties have a >30% land ratio of steep slope, and over one-fifth of counties have a >50% land ratio of steep slope.Therefore, the "innate" ecological vulnerability is relatively high.In terms of "acquired" factors, the high DEV of Yunnan province is first and foremost the result of excessive land development and utilization, particularly reflected in over-reclamation (including deforestation and cultivating on steep slopes) [60], and the abundance of mountainous regions maintaining primitive, extensive and backward practices has resulted in a considerable proportion of slope-cultivated land.The area of suitable cultivated land (including suitable cultivated land in existing cultivated land and unused land) was 4.7076 million hectares in 2020 in the whole province, with only 12.25% being suitable for cultivation.Each county has varying degrees of excessive development and use.Overall, 68.22% of counties have an O RR of over 10%, and about 25% of counties have an O RR of over 20%.This is the root reason behind the generally weak ecological environment for Yunnan's social-economic development.In addition, the measures of resource transformation (consolidation) and ecological environment preservation have not maintained, manifested in the low E IR , with about two-thirds of counties having E IR values of less than 30%.There are significant regional differences in F CR , with some counties having F CR values < 20%.In addition, there is a significant difference between the I BR and the average V ES in different regions.Consequently, it is extremely important to carry out great ecological preservation projects according to local conditions, especially in weak areas with high DEV [66].To promote meeting the SDGs in Yunnan, it is urgent to further significantly reduce ecological vulnerability.Therefore, it is urgent to increase efforts in ecological civilization construction throughout the province, effectively implementing the strategy of VECC [67] and taking practical and effective ecological environment protection and construction measures according to local conditions, especially in terms of developing and utilizing various resources reasonably based on the principles of ecological suitability and environmental friendliness, to steadily reduce the D OEV in various regions, ensuring the coordination of the "population, resources, environment, and economic development" system and the sustainability of regional development.For counties evaluated as being "moderately vulnerable", "highly vulnerable" and "very highly vulnerable", it is necessary to essentially change the mode of resource development and utilization and take ecological and environmental protection and construction measures to significantly reduce ecological vulnerability.
Overall, ecological vulnerability assessment is one of the core contents of vulnerability research, and the development of RS and GIS technology has provided convenient conditions for this study.However, based on existing research results, there is currently no comprehensive multi-phase ecological vulnerability assessment research based on RS and GIS for mountainous provinces.Correspondingly, a suitable, feasible and easy-to-use regional ecological vulnerability assessment system for mountainous areas has not been established, including the congenital natural ecological vulnerability evaluation and the "acquired" ecological vulnerability evaluation based on multi-phase LUCC; thus, it is difficult to judge and obtain the dynamic changes in the degree of ecological vulnerability in mountain areas.And this deficiency is the distinctive feature and innovation of this study.The research results of this study can be more conducive to the implementation of SDGs and the selection of countermeasures in mountainous areas.
On the other hand, in terms of LUCC based on RS, the research on the spatiotemporal evolution characteristics of ecological vulnerability in all counties also represents a further refinement and extension of former vulnerability research.LUCC often has a significant influence on ecosystems, biogeochemistry cycles, climate change, biodiversity, etc. [68,69], which in turn affects regional ecological vulnerability.Therefore, since the 1990s, LUCC has been the core of global environmental change research [70,71], and study achievements continue to emerge.Recently, the study on LUCC has advanced to becoming an interdisciplinary and independent field termed land change science (LCS) [72], which has led to the formation of the land system science (LSS) research group organized by the global land program (GLP) project [73][74][75].The ultimate goal of both LCS and LSS research is to reduce ecological vulnerability and enhance the sustainability of land use and regional development.Therefore, the research in this paper can, to a certain extent, promote the further refinement of LUCC research based on RS and the application of LCS and LSS in the field of regional EVE, sustainable land use and regional development.
Conducting a comprehensive evaluation of county ecological vulnerability and dynamic change research for a mountainous province is indeed a challenging task with extreme difficulty.Especially when collecting, organizing and analyzing basic data from various counties, we often encounter many obstacles.For example, some indicator data from a few counties in specific years may be incomplete or have obvious inaccuracies.This is why we are taking multiple measures simultaneously to overcome this obstacle: Firstly, we will make full use of remote sensing image interpretation results to analyze and supplement missing and obviously inaccurate data.The second is to find ways to collect statistical data from various counties over the years and objectively analyze and determine data that are in line with reality.The third is to conduct necessary research on the ground and specifically implement the relevant situation.By combining these approaches, we can effectively overcome the obstacles encountered during the research process.

Conclusions
This study was based on five periods of LULC data (1980, 1990, 2000, 2010 and 2020) obtained from RS image interpretation, with the county-level region as the assessment unit, adopting organic combination methods of QMCE and QCA.Based on the establishment of an evaluation indicator system of three dimensions of INEV and "acquired" LUEV and LCEV in 129 counties of the province, the D INEV , D EVLU and D EVLC were calculated for 129 counties in 1980,1990,2000, 2010 and 2020.The D OEV of 129 counties from 1980 to 2020 was computed by organically combining the three indicators; this was used to analyze the characteristics of the OEV of each county in the province.The results indicated the following: (1) The INEV of the province was high, and the D INEV value reached 58.21, corresponding to the "highly vulnerable" level.This feature is mainly caused by the high and steep mountain terrain characteristics.In terms of 129 counties, 79.07% of counties were classified as "moderately vulnerable" or above.The regional difference in INEV is large.In the whole province, the "innate" natural ecological environment is the most fragile in the high mountains and canyons of northwest Yunnan, and most counties belong to the "very highly vulnerable" level; the second is the mountainous area of northeast, and most counties belong to the levels of "highly vulnerable" or "very highly vulnerable".Most counties in central and southeastern Yunnan are classified as "highly vulnerable", while the D INEV is relatively low in south and southwest Yunnan.
(2) The D EVLU and D EVLC resulting from acquired human socioeconomic activities were higher in 1980 (40 years ago), and the average D EVLU and D EVLC values in Yunnan were 53.50 and 55.79, which correspond to the classification of "moderately vulnerable" and "highly vulnerable", respectively.After 40 years of ecological protection and construction, the average D EVLU and D EVLC values in Yunnan province fell to 43.94 and 30.01, respectively, in 2020, which are "lowly vulnerable" and "very slightly vulnerable", respectively.Over the past 40 years, the average D EVLU and D EVLC values in Yunnan province dropped by 0.45% and 1.16%, respectively.From 1980 to 2020, the values of D EVLU and D EVLC dropped to various degrees in each county.The basic rule of change of vulnerability grades is that the number of counties with high D EVLU and D EVLC greatly decreased while the number of counties with low D EVLU and D EVLC significantly increased.
(3) From 1980 to 2020, the OEV of Yunnan province and counties decreased significantly.The average D OEV of Yunnan province decreased by 0.51% annually, and the OEV level changed from being "highly vulnerable" in 1980 to "lowly vulnerable" in 2020.This indicates the obvious success of ecological environment preservation, with advancements made from 1980 to 2020, which is also an important achievement from vigorously implementing the strategy of VECC in Yunnan province in recent years.From 1980 to 2020, the D OEV of 129 counties dropped to varying degrees.In terms of the change in the OEV levels, 127 counties (98.45%) have seen a decrease in the level of OEV in the last 40 years.Although two counties (1.55%) had no change in OEV level, their actual D OEV values showed a gradual but somewhat small decrease.During the past 40 years, the basic rule is that the number of counties with high D OEV significantly decreased while the number of counties with low D OEV significantly increased.
(4) There is a significant regional difference in the D OEV within the province.In general, the OEV level is lower in the south, southwest, west and center while being higher in the northwest high mountain and valley areas, the northeast mountainous area and the east and southeast karst areas.
(5) Considering all things, the D OEV of Yunnan province remains high.In 2020, 38.76% of the province's counties were "moderately vulnerable".Deqin, Daguan, Huize, Zhenxiong, Qiaojia, Yongshan, Xichou and Malipo are considered "highly vulnerable", and Dongchuan was classified as "very highly vulnerable".It is time to further enhance the establishment of ecological civilization of the whole province and to take effective measures to preserve and build a beautiful environment in accordance with local conditions, so as to steadily reduce the D OEV in all areas.

Shortcomings and Prospects
Due to the complexity of the regional ecological vulnerability assessment system and the difficulty in obtaining some basic data, this study is not perfect either.For example, in constructing a county-level evaluation index system, due to deficiencies in basic data, various natural disaster data for each county were not obtained; therefore, natural disaster situations were not included in the county-level evaluation index system.For example, when conducting a unified and comprehensive assessment of county ecological vulnerability, typical thematic analysis has not yet taken into account the unique characteristics of some counties.In addition, the evaluation indicators constructed in this study do not involve the discussion of the human footprint, which needs to be addressed in future research.According to Williams et al. [76] (2020), changes in terrestrial human footprints will lead to sustained loss of intact ecosystems, and greater efforts are urgently needed to alleviate the pressure on terrestrial ecosystems from humans.
Therefore, in future research, the following studies are planned to be carried out in this field: The first is to further strengthen the collection and analysis of basic information and data (such as natural disaster data) in each county and to establish a more comprehensive and systematic evaluation system for county ecological vulnerability.The second is to carry out characteristic typical thematic research based on the special characteristics of some counties, so as to organically combine the unified county evaluation on the surface with the typical research on the points, and better provide a solid basis for the implementation of SDGs and the selection of countermeasures in mountainous counties.The third is to conduct research on the impact of human footprint on ecological vulnerability and further expand the ecological vulnerability assessment system.In addition, existing ecological vulnerability information should be combined with other environmental indicators based on monitoring, such as species monitoring based on standardized surveys.Tulloch et al. [77] (2023) showed that different monitoring and management strategies have long-term impacts on endangered animals; Bayraktarov et al. [78] (2021) argue that strengthening active management and monitoring of protected areas to track and report long-term trends across species is important for preventing biodiversity decline.Therefore, in future ecological vulnerability research, it can be combined with the monitoring of additional environmental indicators.

Figure 1 .
Figure 1.Geographical location and DEM map of the research region: (a) location of the province, (b) distribution of all counties and (c) DEM map of the province.

Figure 1 .
Figure 1.Geographical location and DEM map of the research region: (a) location of the province, (b) distribution of all counties and (c) DEM map of the province.

( 1 )
D INEV .D INEV in the county is reflected through four indicators: I MA , I SSA , I HAA and I AAR .Consequently, based on separate calculations of these indicator values, the D INEV value can be calculated using the following equation:

Figure 4 .
Figure 4. Classification map of INEV in Yunnan province.

Diversity 2023 , 26 Figure 5 .
Figure 5. EVE of land use in Yunnan province from 1980 to 2020.

Figure 5 .
Figure 5. EVE of land use in Yunnan province from 1980 to 2020.

Figure 6 .
Figure 6.EVE of land cover in Yunnan from 1980 to 2020.

Figure 6 .
Figure 6.EVE of land cover in Yunnan from 1980 to 2020.

3. 3 .
Figure 7 indicates the changes in D OEV and OEV levels in various regions from 1980 to 2020, primarily exhibiting two features.

3. 3 . 1 .
Figure 7 indicates the changes in DOEV and OEV levels in various regions from 1980 to 2020, primarily exhibiting two features.

Table 1 .
Detailed information on the RS images of the five different phases.

Table 1 .
Detailed information on the RS images of the five different phases.

Table 2 .
The classified area of land use/land cover in Yunnan province in 1980, 1990, 2000, 2010 and 2020.

-Level Land Use Types Second-Level Land Use Types Land Use Area of Various Types (Unit: 10,000 Hectares)
2.3.Methods of EVE 2.3.1.Basic Ideas and Indicator System for EVE (1) Basic Ideas for EVE

Table 3 .
EVE: indicator system, computing approaches, data acquisition techniques and evaluation criteria.LA )I MA = (M AR − Minimum M AR )/MinimumM AR × 100 The Second National Land Survey in Yunnan Province Dam Area Special Survey It depends on the background of the region.It takes the minimum M AR in Yunnan province as the relative optimal value.Total Land Area (T LA ) M AR (Mountain Area Rate) = M LA /T LA × 100% 1.2 I SSA ≥25 • Steep Slope Area (S SA ) I SSA = S SAR /Maximum S SAR × 100 Special Survey of Land Area in Different Climatic Zones and Slopes in Yunnan Province Considering the regional background and rural development needs, the S SAR in the province's relatively largest county is taken as the relative minimum value.The closer the S SAR is to 0, the better the I SSA .Total Land Area (T LA ) S SAR (Steep Slope Area Rate) = S SA /T LA × 100% 1.3 I HAA High-Altitude Area (H A ) I HAA = H AR /Maximum H AR × 100 Special Survey of Land Area in Different Climatic Zones and Slopes in Yunnan Province Considering the background and the needs of rural development, the H AR in the province's relatively largest county is the relative minimum.The closer the H AR is to 0, the better the I HAA is.Total Land Area (T LA ) H AR (High-Altitude Area Rate) = H A /T LA × 100% 1.4 I AAR Average Annual Rainfall (A AR ) I AAR =100 − A AR /Relatively Optimal Value of A AR × 100% Yunnan Agricultural Climate Dataset Considering the regional background, after removing extreme values, the average annual rainfall of relatively large counties in the province is taken as the optimal value.Degrees of Ecological Vulnerability of Land Use (D EVLU ) 2.1 I OR Land Suitable Reclamation Rate (L SRR ) I OR = O RR /Maximum O RR × 100 Land Suitability Evaluation Actual reclamation rate ≤ suitable reclamation rate (i.e., O RR = 0).Actual Land Reclamation Rate (A LRR ) O RR (Over-Reclaimed Rate) = (A LRR − L SRR )/L SRR × 100% RS Image Interpretation 2.2 I BLA Bare Land Area (B LA ) I BLA = B LAR /Maximum B LAR × 100 RS Image Interpretation 0 Total Land Area (T LA ) B LAR (Bare Land Area Rate) = B LA /T LA × 100% 2.3 I EI Paddy Field Area (P EA ) I EI = 100 − E IR /Maximum E IR × 100 RS Image Interpretation The higher the E IR , the lower the vulnerability.Cultivated Area (C A ) E IR (Effective Irrigated Rate of Cultivated Land) = P EA /C A × 100% 2.4 I GY Total Yield of Grain Crops (T YGC )

Table 4 .
The weight of indicators in various levels for regional ecological vulnerability evaluation.

Table 6 .
Statistics of counties with INEV in Yunnan province.

Table 6 .
Statistics of counties with INEV in Yunnan province.

Table 7 .
Number of counties with different LUEV grades in Yunnan province from 1980 to 2020.

Table 8 .
Number of counties with different LCEV grades in Yunnan province from 1980 to 2020.

Table 8 .
Number of counties with different LCEV grades in Yunnan province from 1980 to 2020.

Table 9 .
The number of counties with different OEV grades in the province from 1980 to 2020.

Table 10 )
from 1980 to 2020: the over-reclaimed rate (O RR ) decreased from 17.69% to 14.64%, with a decrease of 17.24%; the bare land area rate (B LAR ) decreased from 15.27% to 10.49%, with a decrease of 31.30%; the forest coverage rate (F CR ) increased from 25.02% to 49.05%, with an increase of 96.04%; the index of biological richness (I BR ) increased from 108.36 to 129.06, with an increase of 19.10% and the ecological services value per unit area (V ES ) increased from 96.09 to 129.63 CNY/hectare, with an increase of 34.90%.The O RR , B LAR , F CR , I BR and V ES in the vast majority of counties have shown significant improvement.This is the basic reason why the OEV of Yunnan province and most counties is gradually decreasing.

Table 10 .
Changes in LULC ecological vulnerability indicators in Yunnan province from 1980 to 2020.