1. Introduction
Building an ecological civilization is a “millennium plan” for sustainable development in China, where the ecological environment has always been a key concern [
1]. With rapid social and economic development, the ecological environment is still deteriorating globally, large ecological patches that maintain ecosystem stability are gradually fragmenting the landscape, and patch connectivity is being reduced, thereby greatly hindering species migration and material energy spread, which are serious threats to ecosystem structure and function [
2]. The SRYR is an important water conservation and livestock farming base in China, and includes numerous ecological patches. The President of China, Xi Jinping, once emphasized at a symposium on ecological protection and high-quality development of the Yellow River Basin: “The Yellow River Basin is an important ecological barrier and an important economic zone in China” [
3]. In recent years, soil erosion, water and soil loss, human activities, rodent browsing, and other phenomena have plagued the SRYR. Consequently, the SRYR ecosystem structure has lost its inherent balance, and suffered functional decline and weakened recovery ability. Therefore, it is urgent to construct and optimize the SRYR ecological network, to scientifically and effectively promote patch connectivity. Additionally, ecological network construction and optimization are highly significant for maintaining ecological security [
4], optimizing ecological patterns [
5], and improving ecosystem quality [
6].
Based on landscape ecology, island biogeography, and population theory, ecological networks comprehensively analyze the distribution and connection of ecological patches in space [
7]. Since the 1990s, ecological network research has involved all ecosystem aspects, including energy flow, material cycles, information transmission, and ecological network structure and composition [
8]. For example, Marc et al. [
9] measured local, regional, and inter-sample network diversity (α-, γ-, and β-diversity) to describe how ecological interactions change over space and time. Isadora et al. [
10] developed a spatial model that identifies and prioritizes riparian corridors to improve landscape connectivity. Ecological network construction simplifies ecological patches in a region into ecological nodes to build an ecological corridor and ecological network. Presently, ecological network construction methods include models such as graph theory [
11], landscape suitability [
12], minimum consumption distance [
13], current theory [
14], and thermodynamic law [
15]. Commonly used software include ConeforSensinode, Circulitscape, Guidos, Zzonation, and Marxan. Among them, the commonly used ecological network construction method is the least cumulative resistance model (MCR). The MCR model construction is mainly about source selection and resistance surface construction.
For source selection, considering the impact of habitat quality and human activities, Gao et al. [
16] extracted ecological sources based on ecosystem service function and ecological sensitivity to construct ecological resistance surface. Yu et al. [
17] selected Dengkou County, a typical ecologically fragile area, as an ecological source area and improved ecological network stability by optimizing the spatial layout of ecological nodes. However, most current studies selected scenic parks or large nature reserves with good habitat patches as ecological sources, although this approach is somewhat subjective. In recent years, a morphological spatial pattern analysis (MSPA) method focusing on structural connections has gradually been integrated into ecological network construction and analysis. Based on Ritters research, Vogt et al. [
18] combined the convolution algorithm with the mathematical morphological mapping algorithm proposed by Soille [
19], and proposed a new method for a landscape connectivity analysis based on the principles of expansion, corrosion, and open–close operation, i.e., morphological spatial pattern analysis. This algorithm can divide the binary image into seven non-overlapping categories (namely, core area, bridge area, loop, branch, edge area, pore, and island patch). Then, the landscape types that are important to maintain patch connectivity are identified, which increases the scientific rigor of the selection of ecological sources and ecological corridors. For example, using the methods of morphological spatial pattern analysis (MSPA) and landscape connectivity, Xiao et al. [
20] combined the graphic theory and quantitative analysis to evaluate the spatio-temporal pattern and network connectivity changes of ecological networks in Zhengzhou. Yi et al. [
14], based on a morphological spatial pattern analysis and circuit theory, focused on the importance of human activities in tropical southwest China to the optimization of the Asian elephant ecological network.
Construction of the resistance surface represents the degree of interference encountered when the target species moves between patches, and it will seriously influence the ecological corridor and ecological network research outputs [
21]. Presently, scholars constructed the resistance surface based on various methods, such as expert scoring, entropy weighting, landscape development intensity index, and biological behavior resistance estimation. Based on the TOPSIS model of entropy weight, Li [
22] constructed an evaluation model of the eco-geological environmental carrying capacity. Li et al. [
23] took the Sichuan-Yunnan ecological barrier as a typical national complex ecological barrier area, and proposed to construct a sustainable Sichuan-Yunnan ecological barrier based on the cycle theory and future land growth changes. Some scholars modified the resistance surface according to the actual situation, to scientifically judge and simulate the potential ecological corridor. Yu Gao [
24] proposed a landscape resistance surface construction method based on a habitat quality assessment, and compared it with a resistance surface constructed using the entropy coefficient and expert scoring methods, and found it more suitable for ecological network research in the scattered Changzhou landscape. However, due to differences in land nutrients and environmental elevations, there may also be differences between the same land use types. Currently, most studies are based on professional knowledge and overall rating of some land use types to construct the landscape resistance surface, which leads to heavy dependence of the landscape resistance surface on the grade coefficient. The MCR model can solve this problem well. Moreover, in general, the combination of MSPA and MCR has been applied to ecological networks in urban landscapes in the central and eastern regions of China, but it has rarely been used in the field of natural landscapes and biological protection in the northwest inland areas.
Although MSPA can identify patches that are important for maintaining landscape connectivity, it still requires assistance from the overall connectivity index (IIC), possibility connectivity index (PC), and equivalent connection proposed by Pascual-Hortal et al. [
25]. In addition to patch abundance and spatial arrangement, these indices combine the dispersal specificity of plant habitats. Wu et al. [
26] took the Guangdong-Hong Kong-Macao Greater Bay Area as an example, and found that the overall ecological connectivity of ecological networks at all scales showed a gradual upward trend, and the overall connectivity index (IIC) and the possible connectivity index (PC) gradually increased with the increase in the maximum dispersal distance of species. Javier Babi Almenar [
27] integrated a landscape index analysis, including the overall connectivity index (IIC), probable connectivity index (PC), and equivalent connectivity index (EC) to show that from 1999 to 2007, habitat fragmentation and loss increased ecological connectivity in Luxembourg, Western Europe. Although the MCR model can judge and simulate the potential ecological corridor by constructing the regional cumulative resistance surface, for corridor relative importance, it is necessary to analyze the interaction strength between patches through a gravity model [
28]. This method mainly combines the network structure index and gravity model to obtain the important patch node rank classification and potential corridor suitability analysis through quantitative calculation, to make the research results more consistent with ecological principles. L. Thiault [
29] evaluated an ecological network of marine-protected areas established on Moorea and French Polynesia through a progressive BACIPS method. For the ecological network evaluation index, scholars also measured the ecological service value based on the probability of occurrence of a certain species in ecological patches [
30]. However, this method lacks consideration of the spatial relationship between landscape ecological elements and is therefore unsuitable.
In summary, this paper took the SRYR alpine grassland as the research object, and based on the MSPA method, identified and extracted the core area landscape type with the best ecological function in the research area. According to the overall connectivity (IIC), probable connectivity (PC), and patch importance (dPC) in the landscape index, core area patches were quantitatively evaluated, to select the ecological sources. The least-LCPs method was used to generate the ecological corridor through the MCR model, and patch interaction intensity was determined based on the gravity model. Then, through betweenness centrality, we identified patches with a better mediating effect as stepped stones to identify potential corridors, and to build the SRYR ecological network. Our research results can provide a basis for the construction and planning of the SRYR ecological network, and also have a certain reference value for ecosystem protection in other regions.
5. Conclusions
In this paper, the SRYR is taken as the research area. From the perspective of ecological landscape connectivity, the ecological source region is identified based on the MSPA method, and the land cover ecological network of the SRYR is constructed and optimized by combining the MCR model. The conclusions are as follows.
(1) The core landscape area of the SRYR was 99,560.85 km2, accounting for 80.53% of the total grassland area, which was mainly distributed in the northeast region, with relatively large marginal areas and void patches. Ring roads, branch lines, and connecting bridges were mainly distributed in the western region. The island had the smallest landscape area.
(2) The dPC values of different ecological sources in the SRYR were significantly different. The number of ecological sources in the west was much less than that in the east, and the northern and southern regions lack the distribution of ecological sources. The patches with large dPC values were located in the east, mainly distributed in cultivated land and wetland, while the dPC value was less than 1 in the western ecological source area. The patches were mainly regional, small, fragmented patches, mainly distributed in bare land.
(3) The minimum cumulative resistance of the ecological network in the SRYR decreased from west to east, and the northwestern region showed the highest resistance, with a value of 4.5. The eastern part had the lowest resistance, with a resistance value of 1.16. Meanwhile, 45 potential ecological corridors were identified based on the MCR model, among which the important corridors were mainly distributed in the eastern part.
(4) In total, 190 planned ecological corridors were obtained by combining 10 core area patches with increased betweenness centrality, which optimized the ecological network of the SRYR. The optimized network structure index was much higher than before.
The results show that the ecological network based on the MCR model is poor, and the eastern and western parts lack connectivity. The optimized ecological network effectively improves the connectivity of the whole ecological patches in the SRYR, and promotes the material exchange and energy flow among the core regions. The results of this study provide an important basis for the sustainable development of the SRYR, and provide a reference for the research and protection of fragile ecosystems. However, the influences of the spatial scale, time scale, and resistance factor on the research results still need to be further discussed.