Ecological Corridors Analysis Based on MSPA and MCR Model—A Case Study of the Tomur World Natural Heritage Region

: The rapid urbanization process, accompanied by the transformation of high-intensive land development and land use, constantly encroaches on habitat patches, making them increasingly fragmented and isolated, which directly inﬂuences the regional landscape pattern and sustainable development. Taking the Tomur World Natural Heritage region as the study area, the morphological spatial pattern analysis (MSPA) method and landscape index method were used to extract the ecological source areas of great signiﬁcance to the construction of ecological corridors. Then, using the minimum cumulative resistance model ( MCR ), the comprehensive resistance surface was constructed and the potential corridors were generated by the minimum cost path method. Finally, according to the gravity model, the important corridors of the study area were designed. Results showed that the MSPA method and MCR model can be used in combination to identify the potential ecological corridors in the study area and clarify the priority of landscape element protection in the study area, which can provide guidance to construct the ecological network and provide reference for other regions as well.


Introduction
The rapid urbanization process, accompanied by the transformation of high-intensive land development and land use, constantly encroaches on habitat patches, making them increasingly fragmented and isolated, which directly influences the regional landscape pattern and sustainable development. Landscape fragmentation reduces the patch area of internal habitat, hinders the operating and regulating ability of normal landscape ecological processes, and damages ecological corridors [1,2]. Landscape connectivity is reduced, which blocks the migration of species and the spread of the material energy and impairs the health of ecosystem integrity, leading to changes in ecosystem services [3]. Therefore, connecting isolated broken ecological patches and stepping stones through potential corridors within the borders can improve the impact of fragmented landscapes on biodiversity and the connectivity of landscape and promote the exchanges of genetic material and species between patches, which would effectively improve the service functions of natural ecosystems and have an important ecological significance [4,5].
Since the 1990s, many scholars have tried to build and analyze ecological corridors at different spatial scales from the perspective of landscape ecology, and many models and methods have been produced [6][7][8]. Among them, the most common methods include the minimum cumulative resistance Based on the MSPA method, this paper takes the Tomur World Natural Heritage site as the research area and identifies and extracts the landscape types of the core areas with the best ecological functions in the study area. According to the integral index of connectivity (IIC), the probability of connectivity (PC), and the delta of PC (d PC ) in the landscape index, the patches in the core area would be quantitatively evaluated so as to select the ecological source in the study area. Ecological corridors are generated through the MCR model using the minimum path method, and the relative importance of corridors is determined based on the gravity model. Then, the patches with good intermediary function are identified as stepping stones through the intermediary centrality, and the potential corridor in the research area is planned, so as to provide a reference for constructing the regional landscape ecological pattern with a complete structure.

Study Area
The Tomur region includes the counties of Zhaosu, Wensu, and Baicheng in Xinjiang, China. The highest peak of the mountain of Tianshan, located in the three counties, is one of the constituent areas of Xinjiang Tianshan World Natural Heritage. East-west mountain systems such as Halketawu, Khantengri, and Tomur cross the study area and divide it into the south and the north. The geographic coordinates are 79 • 27 10"~83 • 04 20" E, 40 • 48 57"~43 • 22 05" N, with a total area of 40,884.07 km 2 ( Figure 1A). The study area belongs to the temperate continental climate, and the climate is significantly affected by the terrain. The north slope is wet, while the south slope is dry. Precipitation is concentrated in the summer and winter, and precipitation in the north slope is significantly higher than that in the south slope. The difference of natural geographical environment, such as the climate of the south and north slope, forms a different vertical natural landscape belt. There are seven large rivers on the south and north slopes, which are mainly fed by melting water from snow and ice [31]. Based on the MSPA method, this paper takes the Tomur World Natural Heritage site as the research area and identifies and extracts the landscape types of the core areas with the best ecological functions in the study area. According to the integral index of connectivity (IIC), the probability of connectivity (PC), and the delta of PC (dPC) in the landscape index, the patches in the core area would be quantitatively evaluated so as to select the ecological source in the study area. Ecological corridors are generated through the MCR model using the minimum path method, and the relative importance of corridors is determined based on the gravity model. Then, the patches with good intermediary function are identified as stepping stones through the intermediary centrality, and the potential corridor in the research area is planned, so as to provide a reference for constructing the regional landscape ecological pattern with a complete structure.

Study Area
The Tomur region includes the counties of Zhaosu, Wensu, and Baicheng in Xinjiang, China. The highest peak of the mountain of Tianshan, located in the three counties, is one of the constituent areas of Xinjiang Tianshan World Natural Heritage. East-west mountain systems such as Halketawu, Khantengri, and Tomur cross the study area and divide it into the south and the north. The geographic coordinates are 79°27 '10 "~ 83°04' 20" E, 40°48 '57 "~ 43°22' 05" N, with a total area of 40,884.07 km 2 ( Figure 1A). The study area belongs to the temperate continental climate, and the climate is significantly affected by the terrain. The north slope is wet, while the south slope is dry. Precipitation is concentrated in the summer and winter, and precipitation in the north slope is significantly higher than that in the south slope. The difference of natural geographical environment, such as the climate of the south and north slope, forms a different vertical natural landscape belt. There are seven large rivers on the south and north slopes, which are mainly fed by melting water from snow and ice [31]. The vegetation distribution in the study area has obvious vertical differentiation characteristics. In the south, a complete set of seven vertical natural belts have developed, with the altitude from low to high successively being the warm temperate desert belt, temperate desert steppe belt, mountain steppe belt, subalpine meadow belt, alpine meadow belt, alpine cushion vegetation belt, and snow and ice belt. In the north, the vegetation coverage is high, and the soil is mainly chestnut and chernozem. Human activities mainly include grazing, tourism, and mining. Due to the lack of reasonable protection and management measures, unreasonable human activities lead to vegetation degradation and soil erosion. In the south, there is a large proportion of saline-alkali land, bare rock, fertile land, and desert, with moisture soil, meadow soil, saline soil, sand soil, and brown desert soil. The vegetation distribution in the study area has obvious vertical differentiation characteristics. In the south, a complete set of seven vertical natural belts have developed, with the altitude from low to high successively being the warm temperate desert belt, temperate desert steppe belt, mountain steppe belt, subalpine meadow belt, alpine meadow belt, alpine cushion vegetation belt, and snow and ice belt. In the north, the vegetation coverage is high, and the soil is mainly chestnut and chernozem. Human activities mainly include grazing, tourism, and mining. Due to the lack of reasonable protection and management measures, unreasonable human activities lead to vegetation degradation and soil erosion. In the south, there is a large proportion of saline-alkali land, bare rock, fertile land, and desert, with moisture soil, meadow soil, saline soil, sand soil, and brown desert soil. The loose soil structure and poor consolidation strength lead to low vegetation coverage and distribution of shallows, and vegetation degradation and desertification are easy to occur, encountering mudslides, heavy rain, and other natural disasters [32].

Data Resource
The main data adopted in this study are the 2015 land use data of Tomur region, the DEM (Digital Elevation Model) data of the Tomur region (http://www.gscloud.cn/), and the road network data of the Tomur region (https://www.openstreetmap.org/). The digital maps were rectified by ground control points to project into WGS_1984_UTM_Zone_44N coordinates. ENVI (The Environment for Visualizing Images, manufactured by American company Exelis Visual Information Solutions) software, which is the ideal software for the visualization, analysis, and presentation of all types of digital imagery (https://envi.software.informer.com/), was used to further enhance remote sensing image data in the study area, so as to improve image data quality and prepare for the next step of data processing. Then, ENVI software was used to monitor and classify the image data in the study area based on visual interpretation, and the Kappa value was 0.84 after an accuracy test of the classification results, which is of very good agreement [33].
The land use classification map of the study area with a grid size of 30 × 30 m was finally obtained ( Figure 1B). According to the actual situation and research purpose of the research area, the land use types in the research area are divided into six types: cultivated land, woodland, grassland, waters, artificial land, and other land. The original resolution of the data is 30 m. All data were used in raster format with the same spatial resolution (i.e., grid cell size corresponding to 1000 m in each raster).

Landscape Pattern Analysis Based on the MSPA Method
The MSPA classification routine starts by identifying core areas, based on user-defined rules for defining connectivity and edge width [34]. Different from the traditional method in which the nature reserves or forest parks were selected directly as the patches or corridors, MSPA is used to identify the important habitat patches and corridors which play an important role in landscape connectivity at the pixel level [35].
After the classification of land use, woodland, grassland, and waters are extracted to be the foreground, and artificial land and other land as the background, a series of image processing methods are used to divide the foreground into seven non-overlapping categories (namely, core, bridge, edge, branch, loop, islet and performation) ( Table 1), and then categories that are important for maintaining connectivity are identified, which increases the scientific nature of the ecological source and ecological corridor selection [36,37].
Then, using the Guidos Toolbox [38], MSPA analysis was conducted on raster-grid data, core patches were obtained through erosion calculation, and the other 6 categories were obtained through a series of mathematical operations. Thus, 7 categories with different functions can be obtained, and the analysis of the results can be conducted. Table 1. Morphological spatial pattern analysis (MSPA) categories and ecological implications.

MSPA Categories Ecological Implication
Core It can be used as the "source" of a variety of ecological processes, most of which are forest parks with large patch areas and large forest farms, etc., which are of great significance for species reproduction and biodiversity protection.

Bridge
The narrow and long areas connecting the patches of different core areas, and which have the characteristics of ecological corridors, which are mostly green belts, which are conducive to the migration of species and the connection of landscape within the territory.

Edge
The transition zone between the marginal zone of the core area and the peripheral nongreen landscape area, which can reduce the impact brought by the external environment and human disturbance, usually the peripheral forest zone of forest parks and large forest farms.

Branch
Only one end is connected to the main patch, mainly an extension of the green space, which is the channel for species diffusion and energy exchange with the peripheral landscape.

Loop
The internal channel of material and energy exchange in the same core area is the shortcut of material and energy exchange in the core area.

Islet
Small patches, which are independent of each other and have low connectivity, are less likely to communicate with other patches in terms of material and energy, and are mostly small green Spaces in urban or rural areas.

Perforation
As a transition region, the edge effect also exists between the core patch and its inner non-green space.

Evaluation of Landscape Connectivity in the Study Area
The level of landscape connectivity in a region can quantitatively characterize whether a certain landscape type is suitable for species exchange and migration, which is of great significance for biodiversity protection and ecosystem balance [39].
At present, in the aspect of landscape connectivity evaluation, the integral index of connectivity (IIC, Equation (1)), the probability of connectivity (PC, Equation (2)), and the delta of PC (d PC , Equation (3)) are commonly used as the important indicators of landscape pattern and function, which can reflect well the degree of connection between core patches in the regional level [40].
where n represents the total number of patches in the region; a i and a j are the areas of patch i and j, respectively; nl ij represents the number of connections between patch i and patch j; P* ij is the maximum product of all path probabilities between patch i and patch j; A L is the total landscape area of the study area; and IIC represents the integral index of connectivity. Please note that 0 ≤ IIC ≤ 1, and if IIC = 0, there is no connection between habitat patches; whereas if IIC = 1, the whole landscape is a habitat patch. PC represents the possible connectivity index of a patch in the landscape of the study area; 0 ≤ PC ≤ 1, the larger the PC value is, the higher the connectivity degree of the patch is; d PC is the change (in %) of the connectivity index after removing one patch. Conefor software was used to evaluate the landscape connectivity of patches in the core areas in the study area through the IIC, PC, and d PC landscape index, and the 10 patches in the core areas with a d PC value greater than 5 were regarded as the source of biological species development and proliferation [41].

Ecological Corridor Construction Based on the MCR Model
The minimum cumulative resistance model (MCR) is characterized by simplicity of construction, extensibility of elements, and universality of application. By calculating the minimum cumulative resistance distance between source and the target to determine the path, it can better reflect the movement possibility and tendency of material energy and biological species in landscape patches [42]. After modification by several experts, the formula of MCR model is as follows: where D ij represents the spatial distance from the source point j to the space unit i and R i represents the resistance coefficient of space unit i. The key to the establishment of the MCR model is the selection of the source and the construction of the resistance surface system. The core process is the competitive diffusion process of the source to the surrounding space units under the resistance surface system [43].
In this paper, combined with the analysis and evaluation results of MSPA and landscape connectivity, 10 sources of ecological processes were selected according to the d PC value of patches in the core area, and resistance factors such as elevation, slope, land type, distance from railway, and distance from other roads were selected to construct the comprehensive resistance surface through the "comprehensive weighted index and method". With reference to relevant research results, the AHP (analytical hierarchy process) method, which is achieved through a pairwise comparison method between the elements for each hierarchical level [30], was used to determine the resistance score and weight of each factor, and a resistance system was constructed, which was divided into 5 resistance scores. The higher the score, the greater the resistance of the diffusion process of biological species [44]. A grid unit size of 30 × 30 m was adopted, and the integrated resistance surface was obtained by the grid calculator as the cost data of the MCR model.
Then, using the ArcGIS software platform and the cost distance module in distance, the minimum cumulative cost distance between each pixel and the nearest unit on the cost surface was calculated using the ecological source and comprehensive resistance surface. Then, the cost path module in distance was used to calculate the minimum cost path from the source to the target patch; thus, 45 potential corridors in the study area were generated. The interaction matrix between patches of 10 ecological source areas was calculated using the gravity model. The interaction intensity between patches of different source areas was evaluated in a quantitative way so that the relative importance of potential corridors in the region could be judged more scientifically.
The gravity model formula is as follows: where G ij is the interaction intensity between patch i and patch j, and N i and N j are the weight coefficients of patch i and patch j, respectively. D ij is the standardized resistance value of potential corridors between patch i and patch j. P i is the overall resistance value of patch i; S i is the area of patch i; and L ij is the cumulative resistance value of potential corridors between patch i and patch j. L max is the maximum resistance value of all corridors in the study area. In summary, the methodology steps are as follows (Figure 2).

Landscape Pattern Analysis of the Tomur Region Based on MSPA
Using the Guidos Toolbox, MSPA analysis was conducted on raster-grid data, core patches were obtained through corrosion calculation, and the remaining six types of landscapes were obtained through a series of mathematical operations such as expansion and reconstruction and skeleton extraction. Thus, seven types of landscapes with different functional types can be obtained ( Figure  3), and their analysis results can be counted (Table 2).

Landscape Pattern Analysis of the Tomur Region Based on MSPA
Using the Guidos Toolbox, MSPA analysis was conducted on raster-grid data, core patches were obtained through corrosion calculation, and the remaining six types of landscapes were obtained through a series of mathematical operations such as expansion and reconstruction and skeleton extraction. Thus, seven types of landscapes with different functional types can be obtained (Figure 3), and their analysis results can be counted (Table 2).

Evaluation of Landscape Connectivity
Conefor software was used to evaluate the landscape connectivity of patches in the core areas in the study area through the IIC, PC, and dPC landscape index, and the 10 patches in the core areas with a dPC value greater than 5 were regarded as the source of biological species development and proliferation (Table 3).

Evaluation of Landscape Connectivity
Conefor software was used to evaluate the landscape connectivity of patches in the core areas in the study area through the IIC, PC, and d PC landscape index, and the 10 patches in the core areas with a d PC value greater than 5 were regarded as the source of biological species development and proliferation (Table 3).

Ecological Corridor Construction Based on the MCR Model
Combined with the analysis and evaluation results of MSPA and landscape connectivity, 10 sources of ecological processes were selected according to the d PC value of patches in the core area, and resistance factors such as elevation, slope, land type, distance from railway, and distance from other roads were selected to construct the comprehensive resistance surface through the "comprehensive weighted index and method" (Figure 4).

Ecological Corridor Construction Based on the MCR Model
Combined with the analysis and evaluation results of MSPA and landscape connectivity, 10 sources of ecological processes were selected according to the dPC value of patches in the core area, and resistance factors such as elevation, slope, land type, distance from railway, and distance from other roads were selected to construct the comprehensive resistance surface through the "comprehensive weighted index and method" (Figure 4).  The higher the score, the greater the resistance of the diffusion process of biological species (Table 4). A grid unit size of 30 × 30 m was adopted, and the integrated resistance surface was obtained by a grid calculator as the cost data of the MCR model.   The higher the score, the greater the resistance of the diffusion process of biological species (Table 4). A grid unit size of 30 × 30 m was adopted, and the integrated resistance surface was obtained by a grid calculator as the cost data of the MCR model.

Resistance Factor Classification Index Evaluation Weight
Then, based on ArcGIS software platform and the cost distance module in distance, the minimum cumulative cost distance between each pixel and the nearest unit on the cost surface was calculated using the ecological source and comprehensive resistance surface. Using the cost path module in distance, 45 potential corridors in the study area were generated ( Figure 5).
Using the gravity model, the interaction matrix between patches of 10 ecological source areas was calculated ( Table 5). The interaction intensity between patches of different source areas was evaluated in a quantitative way so that the relative importance of potential corridors in the region could be judged more scientifically. Table 4. Score and weight of resistance factors.

Resistance Factor Classification Index Evaluation Weight
Elevation Then, based on ArcGIS software platform and the cost distance module in distance, the minimum cumulative cost distance between each pixel and the nearest unit on the cost surface was calculated using the ecological source and comprehensive resistance surface. Using the cost path module in distance, 45 potential corridors in the study area were generated ( Figure 5). Using the gravity model, the interaction matrix between patches of 10 ecological source areas was calculated ( Table 5). The interaction intensity between patches of different source areas was evaluated in a quantitative way so that the relative importance of potential corridors in the region could be judged more scientifically.

Analysis of Landscape Pattern Based on MSPA Method
According to Figure 2 and Table 2, the core area in the study area is about 40721.71 ha, accounting for 84.27% of the total area of woodland. It is mainly concentrated in the northern and southern regions of the study area, while the central region is scarce, leading to poor landscape connectivity in the southern and northern regions of the study area, which is not conductive to the diffusion and communication of organisms and substances between patches in the core area. The area of connection bridge is about 62.81 ha, accounting for 0.13% of the total area of woodland, which is of great significance for the communication of matter and energy between core patches in the study area. As the transition area between the core area and the peripheral nongreen landscape area, the marginal area accounts for about 14.42% of the total area of the forest land and has the edge effect, which can reduce the impact caused by the interference of external factors. Islets were mainly scattered in the central and marginal areas of the study area, accounting for about 0.01% of the total area of woodland. The branch lines with certain connectivity in the core patches and the perforation in the patches accounted for 0.87% and 0.30% of the total forest area, respectively.

Analysis of Landscape Connectivity
It can be seen from Table 3 and Figure 3 that 10 sources of ecological processes were selected according to the d PC value of the core area, showing that the core patches with good connectivity are mainly distributed in the northern part of the study area. Compared with the central and southern areas, the habitat patches in the northern part are more suitable for the migration of biological species and the exchange of material and energy, which can better provide habitats and habitats for species. However, the integral connectivity in the study area is very poor, with a serious fault between the south and the north. On the whole, it is very necessary to construct stepping-stone green patches in the central and southern regions, so as to maintain the sound development of the ecological service function and the balance of the ecosystem in the study area, which needs to be protected.

Analysis of Ecological Corridors Construction
As can be seen from Figure 4 and Table 5, the interaction matrix between patches generated based on the gravity model can measure the strength of the relationship between patches by quantization, and then the importance of corridors between patches to the ecosystem in the region can be inferred. For example, the interaction between patches 4 and 5 was the strongest, and so was the association between these two patches; it is more likely that the biological species can overcome the resistance to migrate, and the communication and transmission of material and energy between the two plaques are simpler; therefore, the protection of the corridor should be strengthened to avoid damage caused by the expansion of construction land.
Due to the farthest distance between patches 6 and 8, communication is least likely between these two patches. Considering the weak connection between the southern and northern regions of the study area, it is necessary to build a bridge for material and energy exchange between the northern and southern regions to ensure the integral connectivity of the study area. Therefore, the corridor between patches 6 and 8 is particularly important, the protection of which is beneficial to the species migration between the north and south regions of the study area.
The results show that the MSPA method and MCR model can be organically integrated, which is a new idea of ecological corridor planning in the world natural heritage area and can provide reference for the construction of an ecological network. When using the MSPA method to analyze the landscape, the research scale is relatively sensitive, and the change of resolution of the grid cell usually influences the analysis results of MSPA. Choosing the appropriate scale is an important content of MSPA research. The area studied in this paper is the Tomur World Natural Heritage area, which has a small area and a large particle size, will lead to the loss of some small patches with good connectivity. Therefore, the particle size is set as 30 m. In addition, the setting of ecological resistance value has a great influence on the construction of an ecological network, but there is no accepted standard at present. Due to the lack of detailed ecological information in the study area, the construction of ecological corridors has not specifically considered the living characteristics of a species. This paper attempts to improve the connectivity of the habitat patches in the study area and the ecological environment quality and avoid the ecological system damage caused by urban expansion, so as to realize the organic unity of protection and growth. Finally, when Conefor software is used to calculate landscape connectivity, the connectivity distance threshold is set as 500 m, and the connectivity probability is 0.5. The setting of the threshold will have a certain impact on the calculation of connectivity index, and more detailed research is needed on the scientific nature of the threshold size.

Conclusions
Taking Tomur as the research area and from the perspective of ecological landscape planning, this paper establishes important ecological corridors based on the MSPA method and MCR model to protect the living environment of ecological patches and species. The MSPA method was adopted to analyze the woodland landscape in the study area. In the selection of an ecological source, not only the value of ecological services was considered, but the structural connection was also emphasized. The landscape types of the core areas with important ecological significance were identified by morphological principles. The importance of the core patches in the study area was quantitatively evaluated by PC and IIC landscape indexes, which changed the subjectivity of the previous artificial ecological source selection to some extent and improved the selection method of the ecological source. Based on the MCR model, multiple resistance factors were comprehensively considered, and the ecological resistance surface was constructed by the comprehensive index weighting method to generate the least resistance corridor between the ecological sources in the study area. According to the strength of interaction between the sources, the Gravity model was used to evaluate the importance of the corridor, and the priority protection order of the corridor was quantitatively analyzed. As an important species location and habitat, the core patch is an important functional node for constructing an ecological network and plays an important role in species migration and diffusion. On the basis of protecting the natural base, reasonable corridor planning suggestions are put forward: (1) The existing corridors, especially the important corridors, should be protected, and the gaps should be repaired to enhance the connectivity of corridors; (2) In areas where there is a lack of connectivity between the relatively important core areas, an appropriate connection bridge corridor or island corridor should be added according to site conditions to promote material exchange and energy flow between the core areas and enhance species diversity.