Assessing the Impact of Village Development on the Habitat Quality of Yunnan Snub-Nosed Monkeys Using the INVEST Model

Simple Summary The Yunnan snub-nosed monkey is one of the most endangered species on the IUCN Red List. The study of its population and habitat quality is important in identifying opportunities for balancing socio-economic development against species conservation in the area’s villages. Such balances are important to protecting and improving habitat diversity and biodiversity. Our habitat quality analysis indicates that increases in socio-economic developments in the villages around the habitat area have decreased both the habitat area and the habitat quality over time. This has resulted in a decline in biodiversity persistence, resilience, and breadth. It also has exacerbated the risk of declining species populations, potentially to extinction. Though focused on the Yunnan snub-nosed monkey, our approach toward the assessment of habitat quality based on species habitat suitability introduces a new perspective for assessing village development impacts on the habitat quality for the conservation of other species. Abstract The habitats of the already endangered Yunnan snub-nosed monkey (Rhinopithecus bieti) are degrading as village economies develop in and around these habitat areas, increasing the depopulation and biodiversity risk of the monkey. The paper aims to show the areas of these monkeys’ high-quality habitats that are at highest risk of degradation by continued village development and hence be the focus of conservation efforts. Our analysis leveraged multiple tools, including primary component analysis, the InVEST Habitat-Quality model, and GIS spatial analysis. We enhanced our analysis by looking at habitat quality as it relates to the habitat suitability for the monkey specifically, instead of general habitat quality. We also focused on the impact of the smallest administrative scale in China—the village. These foci produced a clearer picture of the monkeys’ and villages’ situations, allowing for more targeted discussions on win–win solutions for both the monkeys and the village inhabitants. The results show that the northern habitat for the monkey is currently higher quality than the southern habitat, and correspondingly, the village development in the north is lower than in the south. Hence, we recommend conservation efforts be focused on the northern areas, though we also encourage the southern habitats to be protected from further degradation lest they degrade beyond the point of supporting any monkeys. We encourage developing a strategy that balances ecological protection and economic development in the northern region, a long-term plan for the southern region to reduce human disturbance, increase effective habitat restoration, and improve corridor design.

Habitat-Quality model to study the influence of village development on the habitat quality of Yunnan snub-nosed monkey distribution area in 2018, exploring these three challenges: (1) Analyze and categorize the development status of villages; (2) Analyze the spatial distribution of habitat quality in Yunnan snub-nosed monkey distribution area; (3) With the results of those analyses, determine the impact of village development on the habitat quality of the Yunnan snub-nosed monkey population.

Study Area
The study area ( Figure 1) is in the Three Parallel River region of northwest Yunnan Province (between 29.020 N, 98.030 E in the north and 25.053 N, 99.022 E in the south, with the elevation varying from 1200 m to 5500 m). It is one of the most ecologically significant areas of China in terms of biodiversity, covering approximately 17,000 km 2 across seven counties in Yunnan (Deqin, Weixi, Lanping, Shangri-La, Lijiang, Jianchuan, and Yunlong), with a total human population of about 1,182,500. (1) Analyze and categorize the development status of villages; (2) Analyze the spatial distribution of habitat quality in Yunnan snub-nosed monkey distribution area; (3) With the results of those analyses, determine the impact of village development on the habitat quality of the Yunnan snub-nosed monkey population.

Study Area
The study area ( Figure 1) is in the Three Parallel River region of northwest Yunnan Province (between 29.020 N, 98.030 E in the north and 25.053 N, 99.022 E in the south, with the elevation varying from 1200 m to 5500 m). It is one of the most ecologically significant areas of China in terms of biodiversity, covering approximately 17,000 km 2 across seven counties in Yunnan (Deqin, Weixi, Lanping, Shangri-La, Lijiang, Jianchuan, and Yunlong), with a total human population of about 1,182,500.  The Yunnan snub-nosed monkey lives at a very high elevation (mainly above 3000 m) and is an endangered species on the IUCN Red List [32,33]; 15 groups of these monkeys have been identified previously. We obtained information on the monkey population size as of 2016, with the monkey population being approximately 3000 individuals at that time [6,32].

Land Use and Land Cover
Land Use and Land Cover (LULC) data in 2018 were obtained from a supervised classification on SPOT-5 images (Institute of Forest Inventory and Planning, Yunnan, 2012) with ground-truthing by the Conservation Information Centre of the Nature Conservancy's China program. All data were geo-corrected in ERDAS 9.2 with a root-mean-square (RMS) error < 1. LULC types (such as Armand pine and hemlock) were assigned one of five habitat categories, based on the Yunnan vegetation classification system and the monkey's habitat preferences. These five habitat categories, in declining quality, were optimal habitat, suboptimal habitat, suitable habitat, unfavorable habitat, and highly unfavorable habitat [5].

Villages and Rural Roads
Data on 2572 villages, related rural roads, and their socio-economic status were obtained from the National Geographic Information Resources Catalog Service system (https://www.webmap.cn/main/do?method=index, accessed on 13 December 2016) and Yunnan Digital Rural Network (http://www.ynszxc.gov.cn, accessed on 6 May 2018).

Principal Component Analysis
We selected 2572 villages located in the study area, obtaining 30 indicators for principal component analysis ( Table 1). The 30 indicators were grouped into six categories: the natural resources, population, economy, infrastructure, energy, and educational factors. We used the principal component analysis [34] in R v.2.14.1 to calculate eigenvalues, using eigenvalues greater than 1 to identify the principal components of development factors of the village. This downscaled 30 variables into 8 principal components. We used weightedsum method with these 8 principal components to calculate a comprehensive score for each village, which was used to evaluate the development level of village. Villages were then assigned into one of seven grades (I, II, III, IV, V, VI, and VII) using the Equal Interval Breaks Method ( Table 2). Villages of grade I had the lowest village development, villages of grade VII had the highest.

Plots
The 344 plots were obtained from the 2017 Forest Resources Survey Data of Yunnan Forestry Survey and Design Institute, which used the field survey method. The plots were laid out as squares, each with an area of 4 km 2 .
By regularly updating the survey plot data every 5 years, dynamic changes in forest resource growth and decline were obtained. We used the ecological quality formula of W i X i to calculate a comprehensive score for each plot, where: i is the evaluation index (1, 2, and 3 for types I, II, and III, respectively); X i is the type score value of each evaluation index; W i is the weight of each evaluation index.
We then assigned the type score value of each evaluation index according to the classification criteria ( Table 3). The weights (as a percentage) were determined by via the expert scoring method for each evaluation index against the whole.
According to the comprehensive score of the plots, we used Equal Interval Natural Breaks Method to divide the range to four ecological quality grades (Table 4).

Habitat Quality Evaluation and Spatial Analysis
InVEST-Habitat Quality model [35] was used to evaluate snub monkey's habitat quality based on LULC map from 2018 and biodiversity threat factors including villages, village roads, other non-forestry land, economic forest, cropland, and artificial construction [28,36]. We analyzed the correlation between villages, rural roads, and plots to get maximum distance thresholds of villages and rural roads in ArcGIS 10.6 and R [3,8]. Finally, we scored the parameters related to threat factor weights, suitability, and sensitivity of land use types and obtained threat factor attributes ( Table 5) and sensitivity of land use types to threats (Table 6). We also applied Spatial Autocorrelation (Moran's I) and Hot Spot Analysis (Getis-Ord Gi*) to explore the spatial variation characteristics of the influence of villages on Yunnan snub-nosed monkey habitat quality [37,38].

Impact of Village Development on Habitat Quality
We used the kernel density in ArcGIS to analyze the spatial pattern of the socioeconomic development of the villages [39,40]. Then the impact of villages on the habitat quality of Yunnan snub-nosed monkeys was evaluated based on GIS spatial analysis.

Analysis the Development of Villages
Our village development analysis revealed that more than half the villages (59.84%) were of Grade II (the second-lowest economic development grade). Grade III (just slightly higher economically developed than Grade II) represented 24.49% of the villages. The remaining Grades represented less than 7% of the villages each (Table 7).  Figure 2 shows higher density development of villages in the southern part of the study area than in the north, with the development trend showing a multi-center radiation pattern. This corresponds directly with a lower quality monkey habitat in the south than in the north as shown in Figure 3. pattern. This corresponds directly with a lower quality monkey habitat in the south than in the north as shown in Figure 3.    Figure 4 shows that habitat quality is higher in the northern part of the study ar than in the southern part, that the development of villages in the north is less impactf on the habitat quality of snub-nosed monkeys than in the south, that ecological stabili is poor, and that ecological protection needs to be prioritized.  Figure 4 shows that habitat quality is higher in the northern part of the study area than in the southern part, that the development of villages in the north is less impactful on the habitat quality of snub-nosed monkeys than in the south, that ecological stability is poor, and that ecological protection needs to be prioritized.  The mean value of the habitat quality index for Yunnan snub-nosed monkeys in 2018 was 0.4679, representing a medium habitat quality. Areas of very poor-quality habitat represented over half the total area (7759.39 km 2 , or 51.94%). These very poor-quality habitats were mainly located in the outer and southern parts of the region and had villages densely distributed with socio-economic activities significantly degrading the habitat. The areas of good and excellent habitat quality were less than 30% of the total area and were mainly located in the northern part of the region. Those areas had more suitable habitat areas, less village distribution, and less human interference, resulting in slower habitat degradation (Table 8). Differing socio-economic development in the villages has caused differing degrees of habitat destruction, but overall socio-economic development has resulted in more low-quality habitat area and less high-quality habitat area. The mean value of the habitat quality index for Yunnan snub-nosed monkeys in 2018 was 0.4679, representing a medium habitat quality. Areas of very poor-quality habitat represented over half the total area (7759.39 km 2 , or 51.94%). These very poor-quality habitats were mainly located in the outer and southern parts of the region and had villages densely distributed with socio-economic activities significantly degrading the habitat. The areas of good and excellent habitat quality were less than 30% of the total area and were mainly located in the northern part of the region. Those areas had more suitable habitat areas, less village distribution, and less human interference, resulting in slower habitat degradation (Table 8). Differing socio-economic development in the villages has caused differing degrees of habitat destruction, but overall socio-economic development has resulted in more low-quality habitat area and less high-quality habitat area.

Analysis of the Habitat Quality of the Yunnan Snub-Nosed Monkey
The analysis in Table 9 shows 12 monkey groups with a mean habitat quality value of 0.7408, which is relatively high for the overall habitat. The highest habitat quality was found in monkey group 11 (C11), with a value of 0.9047, while groups C3, C6, and C14 had the lowest habitat quality (Table 9).

Impacts of Village Development on Habitat Quality
There were significant spatial differences in habitat quality distribution patterns. The spatial autocorrelation of habitat quality revealed characteristics of weak agglomeration distribution in the space. The statistical parameters of spatial analysis indicated that the global Moran's I index of habitat quality was 0.0907, and the possibility of agglomeration distribution was less than 1%. Figure 5 reveals that the hot spots and sub-hot-spots of habitat quality were in the northern and southern region as opposed to the central region, and the habitat quality changed greatly, as formerly forested land was converted into non-forested land such as crop land and villages or converted into ecological restoration projects. Conversely, the cold spots and sub-cold-spots were mainly distributed in the central region, which is a relatively suitable habitat for Yunnan snub-nosed monkey with great conversation and little change in habitat quality. Comparing Figures 4 and 5, we found that there was significant overlap between the Very Poor habitat quality area from Figure 4 and the Not Significant area from Figure 5.
If habitat or land use changes are representative of genetic, species, or ecosystem changes, then a low habitat quality will mean a decline in the biodiversity in the habitat and will mean habitat change unfavorable for species survival [16]. As the population of a village increases dramatically, people expand its boundaries, increase the use of cropland, and develop unused land to carry the increased population and provide food [7,41]. In addition, to develop the economy, people plant large amounts of economic forests, leading to habitat fragmentation [42]. Therefore, we need to carry out ecological restoration projects in these areas through cropland reforestation. We need to promote the restoration of the Huashan pine hemlock arrow bamboo forests, spruce fir forests, and mixed coniferous and broad-leaved forests to improve habitat quality [4]. Bamboo and lichen are important food sources for Yunnan snub-nosed monkeys [10]. Forest rangers should restore bamboo plants and lichen, which will provide adequate food conditions for the Yunnan snubnosed monkeys. If habitat or land use changes are representative of genetic, species, or ecosystem changes, then a low habitat quality will mean a decline in the biodiversity in the habitat and will mean habitat change unfavorable for species survival [16]. As the population of a village increases dramatically, people expand its boundaries, increase the use of cropland, and develop unused land to carry the increased population and provide food [7,41]. In addition, to develop the economy, people plant large amounts of economic forests, leading to habitat fragmentation [42]. Therefore, we need to carry out ecological restoration projects in these areas through cropland reforestation. We need to promote the restoration of the Huashan pine hemlock arrow bamboo forests, spruce fir forests, and mixed coniferous and broad-leaved forests to improve habitat quality [4]. Bamboo and lichen are important food sources for Yunnan snub-nosed monkeys [10]. Forest rangers

Impacts of Village Development on Yunnan Snub-Nosed Monkeys
Coordinating the relationship between the socio-economic development of the village and the protection of the ecological environment can promote a harmonious symbiosis between humans and Yunnan snub-nosed monkeys. In places where the socio-economic development of villages is lagging or in decline, local villagers currently rely on massive timber cutting and occupying forest land to develop agriculture and grazing to sustain their livelihoods, threatening the survival of monkeys. There were 10 villages (Zanzhuaro, Amu Shiguang, Muguang Aji, Sebu Feiha, Zaopo, Guilong, Cuding, Zengda, Yidoushe, and Lilinong) directly competing with the C3, C6, and C14 monkey groups ( Figure 6). The government should implement poverty alleviation policies for these villages, removing the need of deforestation to drive economic development, instead developing monkey habitat-friendly eco-tourism and minority cultural industries.

Impacts of Village Development on Yunnan Snub-Nosed Monkeys
Coordinating the relationship between the socio-economic development of the village and the protection of the ecological environment can promote a harmonious symbiosis between humans and Yunnan snub-nosed monkeys. In places where the socio-economic development of villages is lagging or in decline, local villagers currently rely on massive timber cutting and occupying forest land to develop agriculture and grazing to sustain their livelihoods, threatening the survival of monkeys. There were 10 villages (Zanzhuaro, Amu Shiguang, Muguang Aji, Sebu Feiha, Zaopo, Guilong, Cuding, Zengda, Yidoushe, and Lilinong) directly competing with the C3, C6, and C14 monkey groups ( Figure 6). The government should implement poverty alleviation policies for these villages, removing the need of deforestation to drive economic development, instead developing monkey habitat-friendly eco-tourism and minority cultural industries. Habitat quality, habitat patch size, and habitat connectivity mutually trade off against each other [43]. The habitat quality of the C3, C6, and C14 monkey groups is relatively low due to the development on villages. This indicates that village development Habitat quality, habitat patch size, and habitat connectivity mutually trade off against each other [43]. The habitat quality of the C3, C6, and C14 monkey groups is relatively low due to the development on villages. This indicates that village development had a strong influence on habitat quality, with the higher development level of a village translating to a lower habitat quality of the monkey group. We should enhance habitat connectivity and build up ecological corridors to promote gene exchange and conserve genetic diversity in these areas where connectivity with other monkey groups is impeded by agriculture and grazing land, roads, and villages. This is especially true for the isolated monkey groups (C3, C6, and C14) [4][5][6]44].

Implications for Conservation
In addition to protecting Yunnan snub-nosed monkeys through traditional measures such as ecological restoration, enhancing habitat connectivity, building corridors, relocation, and poverty alleviation, we need to take additional, creative measures. Illegal hunting is still an important factor threatening species conservation [45]. It is necessary to combat the illegal wildlife trade, curbing wildlife consumption through interventions and/or shifting consumption to more sustainable alternatives [46,47].
The government must provide ecological compensation to villages with high levels of ecological protection and restoration and establish ecological compensation mechanisms. These are a form of compensation for the damage cause to ecological functions and quality during development. The purpose of these forms of compensation is to improve the environmental quality of damaged areas or to be used to create new areas with similar ecological functions and environmental quality that can further protect, and sustainably use, ecosystem services [48,49]. We also need to strengthen science and education on biodiversity and wildlife conservation so that villagers can participate in species conservation and reap the benefits from it. In general, through a variety of measures, such as education, social influence, legal regulation, and behavioral facilitation, emphasis must be on multidisciplinary applications, multi-subject cooperation, and multi-scale considerations for species conservation [46,47].
The InVEST models is relatively mature and superior to transitional methods in spatial expression and dynamic research. We preliminarily discussed the improvement of adding quantitative village factors to the conventional model and have conducted analysis on the development on villages to include the habitat quality for Yunnan snub-nosed monkeys.
Due to the protection of some data and limited access to other data, the study area only has the distribution area of the Yunnan snub-nosed monkey in northwestern Yunnan. For example, data for some areas in Tibet province is lacking. The research is also lack dynamicity, integrity, and comprehensiveness. Data for some areas of Tibet will be supplemented in future work, and the influence of villages on the habitat quality of monkeys in the whole territory will be studied quantitatively. We also plan to improve the InVEST model parameters and strengthen parameter verification to obtain more effective variable factors, such as climate change factors. We also plan to further improve the variables that reflect the socio-economic development of villages, such as roads at all levels, population spatial distribution km grid data, nighttime lighting data, and GDP spatial distribution km grid data sets.

Conclusions
By using the smallest administrative scale in China-the village-we were able to accurately study the impact of village development on habitat quality, with focus on the Yunnan snub-nosed monkey. The modeling of habitat quality using the InVEST-Habitat Quality Model provided information about the relationship between the spatio-temporal distribution of habitat quality and potential biodiversity. The results of the habitat quality analysis indicate that decreases in habitat area and quality over time meant a decline in biodiversity persistence, resilience, and breadth and exacerbated the risk of declining species populations. In addition, the assessment of habitat quality based on species habitat suitability provides a new perspective to assess the impact of village developments on the habitat quality for the conservation of specific species. While we used the Yunnan snub-nosed monkey in our paper, these approaches can be applied for the conservation of other species.