Next Article in Journal
Quantifying the Potential Vegetation Distribution under Climate Change: The Case of Cryptomeria fortunei in Dongting Lake Watershed, China
Next Article in Special Issue
A Review of Eco-Product Value Realization and Eco-Industry with Enlightenment toward the Forest Ecosystem Services in Karst Ecological Restoration
Previous Article in Journal
Identifying Nutrient Export Hotspots Using a Spatially Distributed Model in Boreal-Forested Catchments
Previous Article in Special Issue
Review on Driving Factors of Ecosystem Services: Its Enlightenment for the Improvement of Forest Ecosystem Functions in Karst Desertification Control
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identifying the Landscape Security Pattern in Karst Rocky Desertification Area Based on Ecosystem Services and Ecological Sensitivity: A Case Study of Guanling County, Guizhou Province

School of Karst Science, Guizhou Normal University, Guiyang 550003, China
*
Author to whom correspondence should be addressed.
Submission received: 7 February 2023 / Revised: 8 March 2023 / Accepted: 16 March 2023 / Published: 19 March 2023

Abstract

:
Ecological environmental security in karst areas is an issue of global concern. Identifying the ecological landscape security pattern (ELSP) is key to promoting environmental protection and alleviating the land development and utilization impacts. Ecological sources (ESs) and ecological corridors (ECs) are important bases for constructing an ELSP. We used five influencing factors (land use type, digital elevation model (DEM), rocky desertification degree, normalized difference vegetation index (NDVI) and slope) to obtain the distribution of the importance and sensitivity values of ecosystem services in Guanling County, Guizhou Province. The probability of the connectivity index (PC) was calculated, and the ES was extracted by combining the importance of ecosystem services, ecological sensitivity, and landscape connectivity. According to the topographic and geomorphological characteristics of Guanling County, seven indicators of elevation, slope, landscape type, degree of stone desertification, distance from rivers, distance from settlements, and distance from roads were selected as resistance factors for the outward expansion of the ESs to calculate the comprehensive resistance surface of Guanling County. Based on the gravity model, an interaction matrix between 10 ESs was constructed, and the magnitude of the interaction forces between the source sites was quantitatively evaluated to distinguish the important ECs and general ECs. The study showed that the total length of the ECs in Guanling County was 509.78 km, and the core area of Guanling County was large, accounting for 65.73% of the ecological landscape area. By assessing the importance of ecosystem services, ecological sensitivity, and landscape connectivity, 10 ES and 45 EC were obtained based on ArcGIS10.8, which constituted the landscape security pattern of Guanling County by ESs and ECs. Suggestions were proposed for a planning layout that will benefit the ecological restoration of Guanling County and environmental protection of the karst region according to the study area characteristics.

1. Introduction

With the acceleration of urbanization, frequent human activities have led to the fragmentation of ecosystems. Ecological security issues have increasingly become the focus of attention of experts and researchers at home and abroad. The construction of regional security patterns is an effective way to improve the ecological environment and achieve regional sustainable development. Yu (1995) [1] proposed the theory of landscape security patterns based on the landscape ecological planning method advocated by Forman [2]. This theory satisfied the theoretical requirements of the reasonable regulation of ecological processes in ecological security research and became an effective way to guide the theory and practice of landscape ecology spatial pattern-ecological process coupling [3]. Ecological security patterns focus on biodiversity conservation, landscape restoration, ecosystem service supply, etc. It is of great significance to analyze the interaction between key elements and ecological processes for ecological security such as ecosystem integrity, reasonable structure, and normal function [4]. At present, the theory of landscape security patterns has been widely used in empirical research. For example, Li et al. (2015), using Changzhou as an example, used the minimum cumulative resistance model (MCR) to construct the resistance surface, determine the quantity and pattern optimization of urban ecological land, and solve urban environmental problems [5]. Based on the distance cost analysis method, Su et al. (2016) constructed an urban regional composite ecological security pattern (ESP) to protect the survival and habitat security of important vegetation, wild animals, and human beings [6]. Liang et al. (2018) comprehensively analyzed the priority protection area and the minimum cost EC, established a more representative, connected, and efficient EC network system, and proposed a new framework for a protection area network composed of protection priority and EC [7]. Peng et al. (2018) evaluated the ecological land risk, quantified the value of ecological land, identified the ESs of Shenzhen, and constructed EC to provide a theoretical framework for the study of urban ecological security patterns [8]. A large number of ES and EC studies have also been carried out in karst areas [9,10,11,12]. However, the research on ELSPs has mainly focused on cities, watersheds, and other regions, and few in karst areas. The construction method of ELSP is divided into three steps: the selection of ESs, the construction of comprehensive resistance surfaces, and the extraction of ECs [13]. The traditional method is subjective when selecting the ESs and usually selects areas with the highest ecological values as the ES, which lacks a certain scientific rigor. Therefore, some researchers have introduced more scientific methods such as the importance of ecosystem services [14], morphological spatial pattern analysis (MSPA) [15], and landscape connectivity analysis to identify ES.
Ecosystem services are the direct or indirect contributions of ecosystems to human well-being. They reflect the human need for ecosystems and are the frontiers of ecological, geographic and economic research [16]. At present, ecosystem services have made many contributions to human development. Researchers at home and abroad have studied the function, quantitative evaluation, and ecosystem service flow of ecosystem services [17,18,19]. Ecological sensitivity is mainly used in ecological environmental protection and spatial planning. Researchers have also combined ecological sensitivity evaluation with spatial security patterns, making it an important indicator for the selection of ESs [20,21,22]. Ecosystem services and ecological sensitivity assessments have gradually been applied to the study of regional security patterns, but little attention has been given to ecological protection and ecological security pattern construction in karst areas.
The term karst originates from the karst plateau of the former Yugoslavia, coupled with the atmosphere, hydrosphere, and biosphere to form the natural ecological environment of karst. It is widely distributed, accounting for 10% of the total Earth area [23,24]. The distribution of karst areas in China exceeds 1.24 million km2, which is concentrated in Yunnan, Guizhou, Guangxi, etc., among which Guizhou karst has a continuous distribution and a wide area of carbonate rocks. Due to the lack of water and soil components in the ecological environment, the Guizhou karst area has low ecological capacity, poor environmental stability, and high sensitivity to variation, and is prone to disasters such as soil erosion, rocky desertification, drought, and flood. The special geomorphic conditions result in complex landforms and strong spatial heterogeneity in the karst areas. At present, the research into environmental governance and rocky desertification control in karst regions has achieved much success [25], and is mainly concentrated on the environment monitoring of karst [26,27,28,29], the spatial distribution of rocky desertification [30,31,32], the evolution process and patterns of rocky desertification [33,34,35], driving factor analysis [36,37,38,39,40], and rocky desertification control, etc. [41,42,43]. Guanling County in Guizhou Province is a typical karst landform area in China. With the acceleration of urbanization and the development of a large number of urban construction facilities, the originally fragile karst environment has been continuously affected by human activity and economic development, and the landscape pattern in the area has changed dramatically [44]. The vegetation coverage has been reduced, cave cracks have developed, water leakage has become serious, the terrain has become extremely fragmented, landscape connectivity has declined, and the ecological environment has become fragile [45]. It is of great significance for the sustainable development of the ecological environment in this karst region to construct a security pattern and optimize its spatial layout. This paper identifies the ecological source of Guanling County based on the importance of ecosystem services, ecological sensitivity, and landscape connectivity, and constructs a comprehensive resistance surface of the study area from the aspects of landscape, terrain, human interference, etc. ArcGIS spatial analysis was used to calculate the potential ecological corridor of Guanling County, and identified the important ecological corridor based on the gravity model to construct the landscape security pattern of Guanling County and provide effective scientific support for the planning and construction of the study area.

2. Materials and Methods

2.1. Overview of the Study Area

Guanling County (Figure 1) is located in central Guizhou Province, the eastern ridge slope of the Yunnan-Guizhou Plateau to the south of the Guangxi hilly slope, within the city of Anshun, 105°15′~105°49′ E, 25°19′~26°05′ N, which is adjacent to Zhenning County, Zhenfeng County, Qinglong County, etc. The total area of the county is 1464 km2. The terrain is high in the northwest and low in the southeast. The elevation is between 370 and 1850 m. The mountainous area accounts for 89% of the total area of the county, and the surface is rugged. The climate in the region is mainly a subtropical monsoon humid climate with sufficient heat, concentrated summer rainfall, and serious soil erosion. The mountains in the territory are part of the Wumeng Mountains. The landform types are complex and diverse. Carbonate rocks are widely distributed, and karst development is strong. Guanling County is a typical karst plateau mountainous area and the area of karst landforms encompasses 83.83% of the county. With socioeconomic development and the influence of human activities, the regional terrain is more fragmented, and the ecological environment has been seriously damaged.

2.2. Data Sources

The basic data used in this study include land use data, digital elevation model (DEM) data, rocky desertification data, and normalized difference vegetation index (NDVI) data.
(1)
Land use data: Derived from the results of the third national land survey (the third national land survey) in Guanling County, according to the land use classification system of the third national land survey (land use status classification GB/T21010-2017) [46]; the data are divided into seven categories.
(2)
DEM data: Obtained from the Chinese Academy of Sciences Data Sharing Center (https://www.resdc.cn/), accessed on 17 May 2022, with a resolution of 30 m, and using ArcGIS10.8 surface analysis to extract the slope data.
(3)
NDVI data: Obtained from the National Ecological Science Data Center (http://www.nesdc.org.cn/), accessed on 19 February 2022, with a resolution of 30 m.
(4)
Rocky Desertification data: Rocky desertification data from Xiong et al. [47], and corrected by the visual interpretation of remote sensing and field survey observations.

2.3. Research Methods

2.3.1. Importance and Sensitivity Analysis of Ecosystem Services

The ecosystem service value is the benefit to humankind from ecosystems [48,49,50]. The ecosystem service value per unit area of land use type was obtained according to the ecosystem service value evaluation research of Xie et al. (2008). (Table 1).
Ecological sensitivity is a comprehensive index used to evaluate the regional ecological environmental quality, land use rationality, and economic development. It is the basis of regional ecological environment planning and management [51]. This paper, combined with the principle of data availability and objectivity, chose five indicators (land use data, DEM, rocky desertification data, NDVI and slope) according to the current research [52] to evaluate the ecological sensitivity of the study area (Table 2).

2.3.2. Landscape Connectivity Analysis

Landscape connectivity refers to the degree to which the landscape promotes or hinders movement between patches [53], which can be reflected by the size of the patches, the distance between similar patches, and the presence of corridors [54], and promotes the communication and exchange of species and biodiversity between patches. MSPA is an image processing method based on mathematical morphology principles such as corrosion, expansion, open operations, and closed operations to measure, identify, and segment raster images [55]. This method can suitably identify the important ecological patches in the study area and obtain seven kinds of landscape elements (core, edge, bridge, islet, perforation, loop and branch) for a clear spatial topological relationship between the target pixel set and the structural elements [56]. Based on the Guidos software, this paper examined the landscape structure of Guanling County by MSPA. The forestland, grassland, water area, cultivated land, and garden land in the seven types of land use in the study area were used as the foreground data, and the other types were used as the background data to extract the core areas. All values were binarized and imported into the Guidos Toolbox software for the MSPA analysis of the images using the eight-neighborhood analysis method.
Next, we used Conefor 2.6 software (available at http://www.conefor.org/ (accessed on 2 March 2023), as proposed by Saura and Torne (2009), to calculate the dPC of the core area. Finally, according to their dPC value from high to low, the patches were divided into four levels, and the importance pattern of the Guanling County landscape connectivity was finally obtained. The calculation formula of dPC is:
d P C = P C P C i r e m o v e P C
where PC is the possible connectivity index of the whole landscape when all patches exist in the landscape and  P C i r e m o v e  is the possible connectivity index value of the remaining patch composition landscape after removing patch i. The higher the dPC value, the higher the importance of the patch in landscape connectivity, and the more obvious its core position in the landscape. The calculation formula of PC is:
P C = i = 1 n j = 1 n a i × a j × a i j A L 2
where n is the total number of habitat nodes in the landscape; ai and aj are the areas of plaque i and plaque j, respectively; A is the total area of the study area; and aij is the maximum final connectivity of all paths between plaque i and plaque j.

2.3.3. Construction of Ecological Resistance Surface

The flow of matter and energy needs to overcome a certain resistance, and the resistance encountered in the outward diffusion of ESs constitutes the comprehensive resistance surface [57]. Based on previous studies [58,59], seven indicators were selected from landscape, topography, and human interference including the elevation, slope, landscape type, degree of stone desertification, distance from rivers, distance from settlements and distance from roads as the resistance factors for the outward expansion of ESs, and these resistance factors were assigned resistance values (from 1 to 9, where the higher the value, the higher the degree of resistance) in a hierarchy (Table 3). Using ArGIS10.8, these seven types of factors were superimposed by finding the mean value to obtain the comprehensive resistance surface of Guanling County.

2.3.4. Construction of ELSP

The core steps of ecological security pattern construction include ecological source identification, ecological resistance surface construction, and key ecological corridor extraction. Among these, the ecological source is the starting point for the outward diffusion of species, which can promote the development of ecological processes and is the key area to ensure regional ecological security [60]. The identification of ESs in Guanling County is mainly based on the analysis of the importance of ecosystem services, ecological sensitivity, and landscape connectivity. ECs are corridors that can provide protection for biodiversity and prevent soil erosion and the loss of ecosystem services [61]. They are the pathways of material and energy flow in the region and key factors in maintaining ecological stability, ecological processes, and ecological function in the region [62]. According to the constructed comprehensive resistance surface, the ECs between the sources were extracted. Using the center point of the ESs as the ecological node, the minimum cost distance from each source to other ESs was calculated by the cost distance using the distance analysis tool in ArcGIS, and then the cost path was used to calculate the minimum cost path between the sources. Finally, the key EC in the study area, which were the channels with the smallest resistance value between the sources, were obtained. Based on the gravity model, we quantitatively evaluated the interaction between each ecological source and determined the importance of the potential corridor (the greater the interaction, the higher the importance of the corridor). The calculation formula of gravity model is as follows:
G i j =   N i N j D i j 2 = ln s i P i ln s j P j L i j L m a x 2 = L m a x 2 ln s i ln s j L i j 2 P i P j
where  G i j  is the interaction force between source patches i and j N i   and   N j is the weight value of source patches,  D i j  is the potential corridor resistance value between source patches i and j P i   and   P j  is the resistance value of source patches i and j L i j  is the potential corridor resistance value from source patches i to j L m a x  is the maximum resistance value of all potential corridors in the study area.

3. Results and Analysis

3.1. Land Use Analysis of Guanling County

According to the land use distribution in Guanling County (Figure 2), the largest area was that of cultivated land at 404.34 km2, accounting for 27.94% of the total area; the forestland area was 345.48 km2, accounting for 23.87% of the total area; the grassland area was 246.31 km2, accounting for 17.02% of the total area; the garden land area was 16.38 km2, accounting for 1.13% of the total area; the water body area was the smallest, at only 15.73 km2, accounting for 1.09% of the total area; and the construction land area was 42.31 km2, accounting for 2.92% of the total area. The area of bare land and other land was 376.50 km2, accounting for 26.02% of the total area. There was more cultivated land, bare land, and other land in Guanling County (Table 4), and the ecological environment was poor. Forest protection and the Grain for Green Project are recommended. The water area was relatively small, and thus water resource protection should be improved.

3.2. Analysis of Importance and Ecological Sensitivity of Ecosystem Services

The importance of ecosystem services can be calculated for different ecosystems to analyze regional differentiation rules and clarify the important areas of ecosystem services. According to research on the evaluation of ecosystem service value by Xie et al. (2008), water areas, forestlands and grasslands, garden lands, cultivated lands, and construction lands were assigned an importance value of 5, 4, 3, 2, and 1, respectively (where the higher the value, the higher the value of the ecosystem services). Then, a distribution map of the importance of ecosystem services in Guanling County (Figure 3) was obtained. The area of extremely important and important areas of ecosystem service value was 607.38 km2, accounting for 41.98% of the total area of the county.
Ecological environmental sensitivity is an important factor restricting urban development. The more sensitive the ecological environment, the more likely environmental problems are to occur in the area. Through the superposition of ecological sensitivity evaluation factors, the ecological sensitivity evaluation results were obtained. The natural breakpoint method was used to divide the ecological sensitivity into five grades: insensitive, mildly sensitive, moderately sensitive, highly sensitive, and extremely sensitive (Figure 4). The highly sensitive and extremely sensitive area of Guanling County was 526.5 km2, accounting for 36.37% of the total area of the county, mainly in the northwest, central, east, and south; the insensitive area was 53.5 km2, accounting for 3.70% of the total area of the county, mainly distributed in the western and southern fringes of the study area. The ecological environment of Guanling County is relatively fragile.

3.3. Analysis of Landscape Connectivity and Landscape Types in Guanling County

There are many indicators used to measure landscape connectivity. The integral index of connectivity (IIC) and probability of connectivity (PC) proposed by Pascual-Hortal and Saura (2006) can evaluate landscape connectivity more effectively than other indicators [63]. Both the IIC and PC are landscape connectivity evaluation indices based on graph theory, but the PC has more advantages and is more rational than the IIC index [64]. Guidos software was used to analyze the landscape structure of Guanling County by MSPA, and the core patches of Guanling County were extracted. Combined with Conefor 2.6 software, the dPC index value of each patch was calculated. According to the level of dPC, the landscape connectivity of Guanling County was divided into four levels: extremely high, high, medium, and low. In the landscape type of Guanling County (Figure 5), the core area was an important species habitat, with an area of 608.61 km2, accounting for 65.73% of the ecological landscape area, and was widely distributed. The loop area had a certain buffer effect on species migration, covering an area of 70.59 km2, accounting for 7.62% of the ecological landscape area. The islet area was the smallest and was the supplementary part of the core area. This area was only 11.71 km2, accounting for 1.30% of the ecological landscape area. The perforation distribution in the study area was relatively uniform, with an area of 45.15 km2, accounting for 4.88% of the ecological landscape area. The edge area was the external boundary of the core area, with an area of 120.06 km2, accounting for 12.97% of the ecological landscape area. The connecting bridge, with an area of 31.30 km2, accounting for 3.39% of the ecological landscape area, had a certain effect on the communication between species. The area of the branch line was 38.53 km2, accounting for 4.16% of the ecological landscape (Table 5).
In the analysis of landscape connectivity in Guanling County (Figure 6), the area of very high and highly connected patches was 518.26 km2, accounting for 35.82% of the total area, mainly in the central, western, and eastern parts of the study area. The area of medium and low connectivity patches was 929.04 km2, accounting for 64.20% of the total area, mainly distributed in the northern and southern regions of the study area, with a relatively large area and wide distribution. The landscape patches in Guanling County are relatively fragmented, and the connectivity between the landscapes is poor, which is mainly related to the special topography of the karst region.

3.4. Comprehensive Resistance Surface and ES

We used the reclassification tool in ArcGIS10.8 to derive the resistance surface of seven single factors, and then superimposed these seven types of factors to find the mean value to obtain the comprehensive resistance surface distribution map of Guanling County (Figure 7). The results show that there is a high resistance area mainly located in the north, central, and south of Guanling County, and the patches in this area are fragmented and the landscape connectivity is poor. The low resistance areas are mainly located in the western and southern fringes, where the landscape connectivity is relatively good and mainly consists of woodlands and grasslands. The selection of ES sites was based on the analysis of ecosystem service importance, ecological sensitivity, and landscape connectivity in Guanling County, and the ecosystem service importance, ecological sensitivity, and landscape connectivity were divided into four levels: very high, high, medium and low, respectively. Patches with two levels, very high and high after overlaid analysis with ArcGIS, were selected as ecological source sites. Finally, 10 ESs were obtained, with an area of 387.48 km2, accounting for 26.78% of the total area. As far as the landscape components are concerned, the source sites were mainly woodlands, waters, and grasslands, and the ecological values of construction land and bare land were relatively low, mostly non-source sites. In terms of distribution area, the southwest region of Guanling County had less distribution, while the central, western, and eastern regions had more distribution. The distribution of ESs was fragmented, and there were many fine patches in the area, with poor inter-patch connectivity and high landscape fragmentation, which is very unfavorable for species dispersal.

3.5. EC Identification and ELSP Construction

EC are the links between each ecological source and the pathways of material flow, allowing species to avoid disturbance during migration [65]. According to the results of the ecological source distribution, using the distance analysis tool in ArcGIS10.8, combined with the comprehensive resistance surface of the study area, the minimum consumption distance between the ecological source points in the study area was calculated, and a total of 45 EC were obtained. The total length of the EC was 509.78 km, the maximum resistance value was 2.06 and the minimum resistance value was 0.23. Twenty of these ECs had resistance values above 1. Based on the gravity model, the interaction matrix between the 10 ecological source sites was constructed (Table 6), which could quantitatively evaluate the magnitude of the interaction between the source sites and discern the importance of potential ECs, and the corridors with interaction forces greater than 10 were identified as important ECs, while the rest were general ECs. The interaction between source 6 and source 8 was the largest with 168.1243, which indicates that the material exchange and transportation between the two sources were more convenient and less costly. The interaction between source 3 and 9 was the smallest at 2.2473, which indicates that the exchange between the two sources requires a higher cost distance and a higher difficulty factor for material exchange and transportation. The ELSP is mainly composed of ecological nodes and EC. Under certain social, economic development and ecological protection conditions, this pattern plays an important role in the regional ecosystem. The ELSP of Guanling County is composed of ESs and ECs (Figure 8).

3.6. Discussion

This study combined ecosystem service values and landscape connectivity to determine the ESs, which is more scientific than previous methods and avoids the subjective selection of ESs. According to the particular topography and geomorphology in the study area and considering the influence of various factors on regional development, the five evaluation factors of rocky desertification degree, land use, NDVI, DEM, and slope were selected to analyze the sensitivity of the study area. The selection of indicators takes into account the geomorphological characteristics of the karst region, which provides a more scientific basis for the construction of the comprehensive resistance surface and the extraction of EC in the study area.
According to the results from this research, we suggest that the protection of the core source area should first be strengthened in the process of ecological land use planning. The ES is a key part of the region. The ecological source area in Guanling County is small, with a low proportion and scattered distribution, which is very unfavorable for the exchange of materials in the region. Therefore, it should be protected in the landscape planning process and development, and construction should be prohibited. Second, ECs should be established, native species should be used as much as possible, and corresponding widths should be set according to the source distribution and the topographic and geomorphic characteristics of the study area to promote species migration and minimize artificial facilities. An organic combination of the ECs identified in this study and the original corridor to form an EC system can improve landscape connectivity between patches. In corridor planning, the potential corridors identified in this study can be organically combined with the original corridors; important ECs can be built in combination with the current corridors to enhance the connectivity between the corridors and the source sites; the construction of general ECs can be combined with some scattered fragmented patches using the existing spatial pattern to increase the connectivity between the source sites. Third, we should focus on the restoration of landscape patches. The general landscape of Guanling County is relatively fragmented. In the planning process, we should focus on the reconstruction and restoration of landscape patches with poor patch connectivity and high ecologically sensitive areas to minimize human interference.

4. Conclusions

This research selected five factors of land use type, DEM, slope, NDVI, and rocky desertification to analyze the ecological sensitivity of Guanling County in Guizhou Province with ArcGIS10.8 software. Through the superposition of various factors, the distribution status of the ecological sensitivity of Guanling County was obtained. Based on the natural breakpoint method, the calculation results were divided into five grades: insensitive, mildly sensitive, moderately sensitive, highly sensitive, and extremely sensitive, and the distributions of seven landscape types (core, islet, perforation, branch, bridge, edge, loop) in the study area were identified. The core area was the main landscape in the study area and an important habitat for species. The total core area was 608.61 km2, accounting for 65.73% of the ecological landscape area, with a wide distribution range. Using Conefor 2.6 software, the dPC index value of each patch was calculated, and the distribution of landscape connectivity in Guanling County was obtained and divided into four grades: extremely high, high, medium, and low. The results showed that the very high and high connectivity patches in Guanling County had an area of 518.26 km2, accounting for 35.82% of the total area, and were mainly distributed in the central, western, and eastern parts of the study area; the medium and low connectivity patches had an area of 929.04 km2, accounting for 64.20% of the total area, and were mainly distributed in the northern and southern parts of the study area, with a relatively large area share and wide distribution. The landscape fragmentation of Guanling County was serious, and the connectivity was poor. According to the overlay of ecosystem service importance, ecological sensitivity, and landscape connectivity in Guanling County, 10 ESs were obtained with an area of 387.48 km2, accounting for 26.78% of the total area of the county. Using the distance analysis tool in ArcGIS10.8, combined with the comprehensive resistance surface of the study area, the minimum consumption distance between each ecological source point was calculated, and a total of 45 ECs were obtained. The maximum resistance value was 2.06 and the minimum resistance value was 0.23, among which 20 ECs had resistance values over 1. Based on the gravity model, the interaction matrix between 10 ESs was constructed to quantitatively evaluate the magnitude of the interaction forces between the source sites and identify the important EC and general EC in Guanling County. Together, ESs and ECs constitute the ELSP of Guanling County. In the southwest part of the study area, the connectivity of landscape patches was relatively poor, the distribution of ES was limited, and the sensitivity was high. In future planning and construction efforts, we should focus on the ecological protection of this area and reduce human interference. However, there were still shortcomings in this study: the study of EC in the study area was mainly based on the minimum cost distance, the width of the corridor was determined by many factors together, influenced by topography, climate change, etc. In this study, the article did not study the width of the corridor in detail because we were limited by the data acquisition and, at the same time, there was not enough time to conduct detailed experiments. The width of the corridor has an important impact on the ecological function of the landscape [66]. Future research will address these inadequacies to improve the scientific accuracy of the study.

Author Contributions

Conceptualization, B.Y. and T.L.; methodology, B.Y. and T.L.; software, T.L. and L.K.; validation, B.Y., T.L. and L.K.; formal analysis, B.Y. and T.L.; investigation, T L., S.L. and R.S.; resources, B.Y. and K.X.; data curation, T.L. and L.K.; writing—original draft preparation, B.Y. and T.L.; writing—review and editing, B.Y. and F.Z.; project administration, B.Y.; funding acquisition, B.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (42261056, 42061028); the Key Science and Technology Program of Guizhou Provence: Poverty alleviation model and technology demonstration for eco-industries derivated from the karst desertification control (No. 5411 2017 Qiankehe Pingtai Rencai); the Oversea Expertise Introduction Program for Discipline Innovation of China: Overseas expertise introduction center for South China Karst eco-environment discipline innovation (No. D17016).

Data Availability Statement

Data supporting this study can be obtained from the Chinese Academy of Sciences Data Sharing Center (https://www.resdc.cn/, accessed on 17 May 2022) and the National Ecological Science Data Center (http://www.nesdc.org.cn/, accessed on 19 February 2022).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Yu, K. Security Patterns in Landscape Planning with a Case in South China. Ph.D. Thesis, Harvard University, Cambridge, MA, USA, 1995. [Google Scholar]
  2. Formam, R.T.T. Land Mosaics: The Ecology of Landscape and Regions; Cambridge University Press: Cambridge, UK, 1995; pp. 32–76. [Google Scholar]
  3. LI, H.; Yi, N.; Yao, W.J.; Wang, S.Q.; Li, Z.Y.; Yang, S.H. Shangri-La county ecological land use planning based on landscape security pattern. Acta Ecol. Sin. 2011, 37, 587–596. (In Chinese) [Google Scholar]
  4. Du, Y.Y.; Hu, Y.N.; Yang, Y.; Peng, J. Building ecological security patterns in southwestern mountainous areas based on ecological importance and ecological sensitivity: A case study of Dali Bai Autonomous Prefecture, Yunnan Province. Acta Ecol. Sin. 2017, 37, 8241–8253. (In Chinese) [Google Scholar]
  5. Li, F.; Ye, Y.P.; Song, B.W.; Wang, R.S. Evaluation of urban suitable ecological land based on the minimum cumulative resistance model: A case study from Changzhou, China. Ecol. Model. 2015, 318, 194–203. [Google Scholar] [CrossRef]
  6. Su, Y.; Chen, X.; Liao, J.; Zhang, H.; Wang, C.; Ye, Y.; Wang, Y. Modeling the optimal ecological security pattern for guiding the urban constructed land expansions. Urban For. Urban Green. 2016, 19, 35–46. [Google Scholar] [CrossRef]
  7. Liang, J.; He, X.; Zeng, G.; Zhong, M.; Gao, X.; Li, X.; Li, X.; Wu, H.; Feng, C.; Xing, W.; et al. Integrating priority areas and ecological corridors into national network for conservation planning in China. Sci. Total Environ. 2018, 626, 22–29. [Google Scholar] [CrossRef]
  8. Peng, J.; Yang, Y.; Liu, Y.X.; Hu, Y.N.; Du, Y.Y.; Meersmans, J.; Qiu, S.J. Linking ecosystem services and circuit theory to identify ecological security patterns. Sci. Total Environ. 2018, 644, 781–790. [Google Scholar] [CrossRef] [Green Version]
  9. Dai, L.; Wang, Z.J. Construction and optimization strategy of ecological security pattern based on ecosystem services and landscape connectivity: A case study of Guizhou Province, China. Environ. Sci. Pollut. Res. 2023, 1–17. [Google Scholar] [CrossRef]
  10. Shuai, N.; Hu, Y.; Gao, M.; Guo, Z.; Bai, Y. Construction and optimization of ecological networks in karst regions based on multi-scale nesting: A case study in Guangxi Hechi, China. Ecol. Inform. 2023, 74, 101963. [Google Scholar] [CrossRef]
  11. Wu, Z.; Xiong, K.; Zhu, D.; Xiao, J. Revelation of coupled ecosystem quality and landscape patterns for agroforestry ecosystem services sustainability improvement in the karst desertification control. Agriculture 2023, 13, 43. [Google Scholar] [CrossRef]
  12. Jiang, Y.; Gao, J.; Wu, S.; Jiao, K. Mediation effect as the component to ecosystem? Establishing the chain effect framework of ecosystem services across typical karst basin in China. CATENA 2023, 221, 106761. [Google Scholar] [CrossRef]
  13. Beier, P.; Majkad, R.; Spencer, W.D. Forks in the road: Choices in procedures for designing wildland linkages. Conserv. Biol. 2008, 22, 836–851. [Google Scholar] [CrossRef]
  14. Wang, H.; Ma, X.; Du, Y. Constructing ecological security patterns based on ecological service importance and ecological sensitivity in Guangdong Province. Acta Ecol. Sin. 2021, 41, 1705–1715. (In Chinese) [Google Scholar] [CrossRef]
  15. Yao, C.Y.; An, R.; Dou, C.; Liu, Y.L. Construction and evaluation of forest ecological network in the Three Gorges Reservoir area based on MSPA and MCR models. J. Resour. Environ. Yangtze River Basin 2022, 31, 1953–1962. (In Chinese) [Google Scholar]
  16. Costanza, R.; Kubiszewski, l. The authorship structure of “ecosystem services” as a transdisciplinary field of scholarship. Ecosyst. Serv. 2012, 1, 16–25. [Google Scholar] [CrossRef] [Green Version]
  17. Li, D.; Wu, S.; Liu, L.; Liang, Z.; Li, S. Evaluating regional water security through a freshwater ecosystem service flow model: A case study in Beijing-Tianjin-Hebei region, China. Ecol. Indic. 2017, 81, 159–170. [Google Scholar] [CrossRef]
  18. Darvill, R.; Lindo, Z. The inclusion of stakeholders and cultural ecosystem services in land management trade-off decisions using an ecosystem services approach. Landsc. Ecol. 2016, 31, 533–545. [Google Scholar] [CrossRef]
  19. Shi, Y.; Shi, D.; Zhou, L.; Fang, R. Identification of ecosystem services supply and demand areas and simulation of ecosystem service flows in Shanghai. Ecol. Indic. 2020, 115, 106418. [Google Scholar] [CrossRef]
  20. Zhang, Q.; Zhang, T. Land consolidation design based on an evaluation of ecological sensitivity. Sustainability 2018, 10, 3736. [Google Scholar] [CrossRef] [Green Version]
  21. Chi, Y.; Zhang, Z.; Gao, J.; Xie, Z.; Zhao, M.; Wang, E. Evaluating landscape ecological sensitivity of an estuarine island based on landscape pattern across temporal and spatial scales. Ecol. Indic. 2019, 101, 221–237. [Google Scholar] [CrossRef]
  22. Ersayin, K.; Tagil, S. Ecological sensitivity and risk assessment in the Kizilirmak Delta. Fresenius Environ. Bull. 2018, 26, 6508–6516. [Google Scholar]
  23. Gao, G.L.; Deng, Z.M.; Xiong, K.N.; Su, X.L. The Call and Hope of Karst; Guizhou Science and Technology Press: Guiyang, China, 2003. (In Chinese) [Google Scholar]
  24. Yuan, D.X. Global view on Karst rock desertification and integrating control measures and experiences of China. J. Grassl. Sci. 2008, 9, 19–25. (In Chinese) [Google Scholar]
  25. Lu, Y.R. Karst Development Mechanism and Research Directions of Developing Engineering Construction Effect. J. Earth Sci. 2016, 37, 419–432. (In Chinese) [Google Scholar] [CrossRef]
  26. Xiong, K.N.; LI, J.; Long, M.Z. Features of Soil and Water Loss and Key Issues in Demonstration Areas for Combating Karst Rocky Desertification. J. Geogr. 2012, 67, 878–888. (In Chinese) [Google Scholar]
  27. Yue, Y.; Liu, B.; Wang, K.; Li, R.; Zhang, B.; Zhang, C.; Chen, H. RETRACTED ARTICLE: Using remote sensing to quantify the fractional cover of vegetation and exposed bedrock within a complex landscape: Applications for karst rocky desertification monitoring. Environ. Monit. Assess. 2013, 185, 1025. [Google Scholar] [CrossRef] [Green Version]
  28. Tong, X.; Brandt, M.; Yue, Y.; Horion, S.; Wang, K.; De Keersmaecker, W.; Tian, F.; Schurgers, G.; Xiao, X.; Luo, Y.; et al. Increased vegetation growth and carbon stock in China karst via ecological engineering. Nat. Sustain. 2018, 1, 44–50. [Google Scholar] [CrossRef]
  29. Qi, X.; Zhang, C.; Wang, K. Comparing remote sensing methods for monitoring Karst rocky desertification at sub-pixel scales in a highly heterogeneous Karst Region. Sci. Rep. 2019, 9, 13368. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Jiang, Y.; Li, L.; Groves, C.; Yuan, D.; Kambesis, P. Relationships between rocky desertification and spatial pattern of land use in typical karst area, Southwest China. Environ. Earth Sci. 2009, 59, 881–890. [Google Scholar] [CrossRef]
  31. Bai, X.-Y.; Wang, S.-J.; Xiong, K.-N. Assessing spatial-temporal evolution processes of karst rocky desertification land: Indications for restoration strategies. Land Degrad. Dev. 2013, 24, 47–56. [Google Scholar] [CrossRef]
  32. Yang, Q.; Jiang, Z.; Yuan, D.; Ma, Z.; Xie, Y. Temporal and spatial changes of karst rocky desertification in ecological reconstruction region of Southwest China. Environ. Earth Sci. 2014, 72, 4483–4489. [Google Scholar] [CrossRef]
  33. Yue, Y.; Wang, K.; Zhang, B.; Jiao, Q.; Liu, B.; Zhang, M. Remote sensing of fractional cover of vegetation and exposed bedrock for karst rocky desertification assessment. Procedia Environ. Sci. 2012, 13, 847–853. [Google Scholar] [CrossRef] [Green Version]
  34. Ying, B.; Xiao, S.Z.; Xiong, K.N.; Cheng, Q.W.; Luo, J.S. Comparative studies of the distribution characteristics of rocky desertification and land use/land cover classes in typical areas of Guizhou province, China. Environ. Earth Sci. 2014, 71, 631–645. [Google Scholar] [CrossRef]
  35. Gao, J.; Wang, H. Temporal analysis on quantitative attribution of karst soil erosion: A case study of a peak-cluster depression basin in Southwest China. CATENA 2019, 172, 369–377. [Google Scholar] [CrossRef]
  36. Jiang, Z.; Lian, Y.; Qin, X. Rocky desertification in Southwest China: Impacts, causes, and restoration. Earth Sci. Rev. 2014, 132, 1–12. [Google Scholar] [CrossRef]
  37. Yan, X.; Cai, Y.L. Multi-scale anthropogenic driving forces of karst rocky desertification in southwest China. Land Degrade. Dev. 2015, 26, 193–200. [Google Scholar] [CrossRef]
  38. Zhang, S.; Zhang, Y.; Xiong, K.; Yu, Y.; Min, X. Changes of leaf functional traits in karst rocky desertification ecological environment and the driving factors. Glob. Ecol. Conserv. 2020, 24, e01381. [Google Scholar] [CrossRef]
  39. Peng, X.D.; Dai, Q.H. Drivers of soil erosion and subsurface loss by soil leakage during karst rocky desertification in SW China. Int. Soil Water Conserv. Res. 2022, 10, 217–227. [Google Scholar] [CrossRef]
  40. Guo, B.; Wei, C.; Yu, Y.; Liu, Y.; Li, J.; Meng, C.; Cai, Y. The dominant influencing factors of desertification changes in the source region of Yellow River: Climate change or human activity? Sci. Total Environ. 2022, 813, 152512. [Google Scholar] [CrossRef] [PubMed]
  41. Yan, D.; Zhong, C.J. Characteristic of rocky desertification and comprehensive improving model in karst peak-cluster depression in Guohua, Guangxi, China. Procedia Environ. Sci. 2011, 10, 2449–2452. [Google Scholar] [CrossRef] [Green Version]
  42. Wu, Z.; Zhu, D.; Xiong, K.; Wang, X. Dynamics of landscape ecological quality based on benefit evaluation coupled with the rocky desertification control in South China Karst. Ecol. Indic. 2022, 138, 108870. [Google Scholar] [CrossRef]
  43. Yao, B.; Yue, X.J.; Huang, P.; Li, Y.-H. The Qing-Long model: China provides a solution to the karst rocky desertification challenge. Acta Ecol. Sin. 2022. [Google Scholar] [CrossRef]
  44. Li, Y.X.; Li, J.; Chen, H.; Wang, Z.J. Landscape connectivity evaluation and temporal-spatial characteristics of Guiyang City from 2008 to 2017 based on MSPA and MCR models. J. Ecol. 2022, 41, 1240–1248. (In Chinese) [Google Scholar] [CrossRef]
  45. Wang, S.J.; Li, Y.B.; Li, R.L. Background, evolution and management of karst rocky desertification. J. Quat. Study 2003, 1, 657–666. (In Chinese) [Google Scholar]
  46. General Administration of Quality Supervision; Inspection and Quarantine of the People’s Republic of China, National Standardization Administration of China. GB/T21010-2017; Land Use Status Classification. Standards Publishing House: Beijing, China, 2017. Available online: https://openstd.samr.gov.cn/bzgk/gb/newGbInfo?hcno=224BF9DA69F053DA22AC758AAAADEEAA (accessed on 22 September 2022). (In Chinese)
  47. Xiong, K.N.; Yuan, J.Y.; Fang, Y.; Lan, A.J.; Chen, Q.W.; Zeng, F.Q.; Ying, B.; Zhang, W.W.; Hu, J.; Zhou, Z.F.; et al. Atlas of Comprehensive Prevention and Control of Karst Rocky Desertification in Guizhou (2006–2050); Guizhou People’s Publishing House: Guiyang, China, 2007. (In Chinese) [Google Scholar]
  48. Ahern, J.; Cilliers, S.; Niemelä, J. The concept of ecosystem services in adaptive urban planning and design: A framework for supporting innovation. Landsc. Urban Plan. 2014, 125, 254–259. [Google Scholar] [CrossRef] [Green Version]
  49. Warsame, A.A.; Abdi, A.H. Towards sustainable crop production in Somalia: Examining the role of environmental pollution and degradation. Cogent Food Agric. 2023, 9, 2161776. [Google Scholar] [CrossRef]
  50. Raihan, A.; Pavel, M.I.; Muhtasim, D.A.; Farhana, S.; Faruk, O.; Paul, A. The role of renewable energy use, technological innovation, and forest cover toward green development: Evidence from Indonesia. Innov. Green Dev. 2023, 2, 100035. [Google Scholar] [CrossRef]
  51. Xu, W.; Wang, J.; Zhang, M.; Li, S. Construction of landscape ecological network based on landscape ecological risk assessment in a large-scale opencast coal mine area. J. Clean. Prod. 2021, 286, 125523. [Google Scholar] [CrossRef]
  52. Yang, T.R.; Kuang, W.H.; Liu, W.D.; Liu, A.L.; Pan, T. Optimization of ecological spatial structure of Guanzhong urban agglomeration based on ecological security pattern. J. Geogr. Res. 2017, 36, 441–452. (In Chinese) [Google Scholar]
  53. Taylor, P. Connectivity is a vital element of landscape structure. Oikos 1993, 68, 571–573. [Google Scholar] [CrossRef] [Green Version]
  54. Chen, L.D.; Fu, B.J. The ecological significance and application of landscape connectivity. J. Ecol. 1996, 15, 37–42. [Google Scholar]
  55. Gandini, M.L.; Lara, B.D.; Moreno, L.B.; Cañibano, M.A.; Gandini, P.A. Landscape dynamics of Paspalum quadrifarium grasslands analyzed by Morphological Spatial Pattern Analysis (MSPA). PeerJ 2018. [Google Scholar] [CrossRef] [Green Version]
  56. Xu, F.; Yin, H.W.; Kong, F.H.; Xu, J.G. Developing ecological networks based on mspa and the least-cost path method: A case study in bazhong western new district. J. Ecol. 2015, 35, 6425–6434. (In Chinese) [Google Scholar]
  57. Gao, X.W.; Feng, Z.J.; Ge, J.F. Analysis of landscape security pattern in Western Mountains of Shijiazhuang. In Proceedings of the Geoinformatics 2008 and Joint Conference on GIS and Built Environment: Geo-Simulation and Virtual GIS Environments, Guangzhou, China, 28–29 June 2008; Volume 7145, p. 71451. [Google Scholar] [CrossRef]
  58. Xiao, S.; Wu, W.; Guo, J.; Ou, M.; Pueppke, S.G.; Ou, W.; Tao, Y. An evaluation framework for designing ecological security patterns and prioritizing ecological corridors: Application in Jiangsu Province, China. Landsc. Ecol. 2020, 35, 2517–2534. [Google Scholar] [CrossRef]
  59. Su, X.; Zhou, Y.; Li, Q. Designing ecological security patterns based on the framework of ecological quality and ecological sensitivity: A case study of Jianghan Plain, China. Int. J. Environ. Res. Public Health 2021, 18, 8383. [Google Scholar] [CrossRef]
  60. Wang, Z.; Shi, P.; Zhang, X.; Tong, H.; Zhang, W.; Liu, Y. Research on landscape pattern construction and ecological restoration of Jiuquan city based on ecological security evaluation. Sustainability 2021, 13, 5732. [Google Scholar] [CrossRef]
  61. Liu, Z.; Gan, X.; Dai, W.; Huang, Y. Construction of an ecological security pattern and the evaluation of corridor priority based on ESV and the “importance–connectivity” index: A case study of Sichuan province, China. Sustainability 2022, 14, 3985. [Google Scholar] [CrossRef]
  62. Wu, M.Q.; Hu, M.M.; Wang, T.; Fan, C.; Xia, B.C. Recognition of urban ecological source area based on ecological security pattern and multi—Scale landscape connectivity. Acta Ecol. Sin. 2019, 39, 4720–4731. (In Chinese) [Google Scholar] [CrossRef]
  63. Xiong, C.N.; Wei, H. Lan Mingjuan. Landscape connectivity of urban green space in Chongqing. Chin. J. Ecol. 2008, 5, 2237–2244. (In Chinese) [Google Scholar]
  64. Garcia-Lozano, C.; Varga, D.; Pintó, J.; Roig-Munar, F. Landscape connectivity and suitable habitat analysis for wolves (Canis lupus L.) in the Eastern Pyrenees. Sustainability 2020, 12, 5762. [Google Scholar] [CrossRef]
  65. Yu, K. Security patterns and surface model in landscape ecological planning. Landsc. Urban Plan. 1996, 36, 1–17. [Google Scholar] [CrossRef]
  66. Chen, X.; Peng, J.; Liu, Y.X.; Yang, Y.; Li, G. Construction of ecological security pattern of Yunfu City based on ‘importance-sensitivity-connectivity’ framework. Geogr. Res. 2017, 36, 471–484. (In Chinese) [Google Scholar] [CrossRef]
Figure 1. Geographical location map of Guanling County.
Figure 1. Geographical location map of Guanling County.
Forests 14 00613 g001
Figure 2. Guanling County land use distribution map.
Figure 2. Guanling County land use distribution map.
Forests 14 00613 g002
Figure 3. Importance distribution map of ecosystem services in Guanling County.
Figure 3. Importance distribution map of ecosystem services in Guanling County.
Forests 14 00613 g003
Figure 4. Ecological sensitivity distribution map of Guanling County.
Figure 4. Ecological sensitivity distribution map of Guanling County.
Forests 14 00613 g004
Figure 5. MSPA landscape type map of Guanling County.
Figure 5. MSPA landscape type map of Guanling County.
Forests 14 00613 g005
Figure 6. Guanling County landscape connectivity distribution map.
Figure 6. Guanling County landscape connectivity distribution map.
Forests 14 00613 g006
Figure 7. Guanling County comprehensive resistance surface distribution map.
Figure 7. Guanling County comprehensive resistance surface distribution map.
Forests 14 00613 g007
Figure 8. Landscape security pattern layout of Guanling County.
Figure 8. Landscape security pattern layout of Guanling County.
Forests 14 00613 g008
Table 1. Ecosystem service value of land use types.
Table 1. Ecosystem service value of land use types.
Value of Ecosystem ServicesLand Use TypeImportance of Ecological Services
5Water bodyExtremely important
4Forestland and grasslandHighly important
3Garden plotModerately Important
2Cultivated landSlightly important
1Construction land and other landUnimportant
Table 2. Ecological environment sensitivity the evaluation factor classifications and weights.
Table 2. Ecological environment sensitivity the evaluation factor classifications and weights.
Sensitivity AssignmentNDVIDEM/mSlope/°Land Use TypeRocky Desertification
1≤0.35≤500≤5Construction LandExtremely strong rocky desertification
3(0.35, 0.50](500, 800](5, 15]Other landIntense rocky desertification
5(0.50, 0.65](800, 1100](15, 25]Cultivated landModerate rocky desertification
7(0.65, 0.75](1100, 1500](25, 35]Grassland and Garden landMild rocky desertification
9>0.75>1500>35Forestland and Water bodyNo rocky desertification
Table 3. Ecological resistance surface evaluation system.
Table 3. Ecological resistance surface evaluation system.
Resistance ValueDistance from River/mDistance from the Settlement/mDistance from Road/mDegree of Stone DesertificationLandscape TypeElevation/mSlope/°
1≤500>1000≤1000No rocky desertificationForestland and Water body≤500≤5
3(500, 1000](1000, 800](500, 1000]Mild rocky desertificationGrassland and Garden land(500, 800](5, 10]
5(1000, 1500](800, 500](200, 500]Moderate rocky desertificationCultivated land(800, 1100](10, 15]
7(1500, 2000](200, 500](100, 200]Intense rocky desertificationOther land(1100, 1500](15, 25]
9>2000≤200≤100Extremely strong rocky desertificationCultivated land>1500>25
Table 4. Areas of various land use types in Guanling County.
Table 4. Areas of various land use types in Guanling County.
Land Use TypesArea (km2)Proportion (%)
Forestland345.4823.87
Grassland246.3117.02
Garden land16.381.13
Water body15.731.09
Cultivated land404.3427.94
Construction land42.312.92
Other land376.5026.02
Table 5. Landscape types in Guanling County.
Table 5. Landscape types in Guanling County.
Landscape TypeArea (km2)Proportion of Total Area of Ecological Landscape (%)
Core608.6165.73
Islet11.711.26
Perforation45.154.88
Edge120.0612.97
Loop70.597.62
Bridge31.303.38
Branch38.534.16
Table 6. Interaction matrix between sources based on gravity model calculations.
Table 6. Interaction matrix between sources based on gravity model calculations.
Number of Source12345678910
1023.5564.15394.76674.7658.95567.9549.546743.410619.5628
2 07.01238.55734.946624.623824.646221.94936.489223.4465
3 022.334217.325416.247932.035210.47922.24734.4374
4 019.636812.92732.96078.16832.28633.911
5 010.36216.85257.28892.49143.4039
6 078.5538168.12433.923717.3001
7 030.9580 3.50739.6917
8 05.344239.4447
9 09.4984
10
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ying, B.; Liu, T.; Ke, L.; Xiong, K.; Li, S.; Sun, R.; Zhu, F. Identifying the Landscape Security Pattern in Karst Rocky Desertification Area Based on Ecosystem Services and Ecological Sensitivity: A Case Study of Guanling County, Guizhou Province. Forests 2023, 14, 613. https://0-doi-org.brum.beds.ac.uk/10.3390/f14030613

AMA Style

Ying B, Liu T, Ke L, Xiong K, Li S, Sun R, Zhu F. Identifying the Landscape Security Pattern in Karst Rocky Desertification Area Based on Ecosystem Services and Ecological Sensitivity: A Case Study of Guanling County, Guizhou Province. Forests. 2023; 14(3):613. https://0-doi-org.brum.beds.ac.uk/10.3390/f14030613

Chicago/Turabian Style

Ying, Bin, Ting Liu, Li Ke, Kangning Xiong, Sensen Li, Ruonan Sun, and Feihu Zhu. 2023. "Identifying the Landscape Security Pattern in Karst Rocky Desertification Area Based on Ecosystem Services and Ecological Sensitivity: A Case Study of Guanling County, Guizhou Province" Forests 14, no. 3: 613. https://0-doi-org.brum.beds.ac.uk/10.3390/f14030613

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop