Next Article in Journal
Tensile Properties of Natural Fibre-Reinforced FDM Filaments: A Short Review
Next Article in Special Issue
Evolution and Optimization of an Ecological Network in an Arid Region Based on MSPA-MCR: A Case Study of the Hexi Corridor
Previous Article in Journal
Confidence in Greenhouse Gas Emission Estimation: A Case Study of Formaldehyde Manufacturing
Previous Article in Special Issue
Influence of Seasonal Water Level Fluctuations on Food Web Structure of a Large Floodplain Lake in China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Construction of Green Ecological Network in Qingdao (Shandong, China) Based on the Combination of Morphological Spatial Pattern Analysis and Biodiversity Conservation Function Assessment

1
College of Landscape Architecture and Forestry, Qingdao Agricultural University, Qingdao 266000, China
2
Laoshan Branch of Qingdao Ecological Environmental Bureau, Qingdao 266100, China
*
Author to whom correspondence should be addressed.
Sustainability 2023, 15(24), 16579; https://doi.org/10.3390/su152416579
Submission received: 20 October 2023 / Revised: 29 November 2023 / Accepted: 4 December 2023 / Published: 6 December 2023

Abstract

:
Building urban green ecological network systems and increasing urban and rural landscape connectivity are effective ways to improve urban biodiversity and landscape sustainability. The ecological sources in the main urban area of Qingdao City (Shandong, China) were identified based on morphological spatial pattern analysis (MSPA) combined with a biodiversity conservation function assessment, with the ecological corridors established and the ecological network structure optimized. The results showed that (1) the study area lacked high-quality patches with strong landscape connectivity; (2) the potential green ecological network of the study area was composed of 38 ecological sources, 703 ecological corridors, and 284 ecological nodes, effectively connecting urban and suburban green spaces; (3) after optimization, the green ecological network contained a total of 223 important corridors and 61 key nodes, with significantly increased network connectivity; (4) the optimal ecological corridor width in Qingdao was determined to be 30 m. Our study provided important guidance for the construction of ecological security patterns and scientific evidence to support urban green space planning and sustainable development in Qingdao.

1. Introduction

When planning urban green space systems, it is fundamentally important to improve the urban biodiversity and landscape sustainability of urban green spaces if the spaces are limited in number. The accurate identification of effective environmental strategies used to introduce the abundant biological resources of the suburbs into the city has become a popular area of research [1]. In recent decades, a series of environmental impairments have been caused by the continuous expansion of urban construction land, i.e., the ecosystem service value of landscape patches has decreased; the fragmentation of urban green patches has increased; the landscape connectivity between green spaces has decreased; ecological processes, such as energy flow, matter cycle, and information flow, have weakened; and the communication between humans and nature has gradually faded away [2]. Some researchers have considered biological rhythms to be a crucial factor affecting urban biodiversity [3]. However, core patches and ecological corridors are considered to have the strongest positive effects on biodiversity in high-density urban areas [4]. Therefore, it is crucial to construct a reasonable urban green space network system based on both morphological spatial pattern analysis (MSPA) and a biodiversity conservation function (BCF) assessment to improve biological flow efficiency and enhance the biodiversity and sustainability of urban green spaces.
The construction of green ecological networks will effectively connect fragmented sources through ecological corridors to promote the ecological process among patches, increase urban biodiversity, and improve the sustainability of regional ecosystems [5,6]. Numerous studies have investigated the construction and optimization of urban ecological networks, and the fundamental paradigm “identification of ecological sources → construction of resistance surface → extraction of ecological corridor” has been well established [7]. In particular, during the identification of ecological sources, most studies have directly designated patches with high habitat quality, e.g., nature reserves and scenic spots, as ecological sources, and alternatively identified ecological patches based on indicators, e.g., the importance of ecosystem services and ecological sensitivity [8]. However, most of these studies have largely ignored the consideration of the integrity of patches in a landscape’s space. This problem could be adequately solved by applying MSPA, which has been widely used in environmental investigations of urban planning and development. MSPA can not only identify and classify ecological patches based on spatial morphology, but also reveal the structural connections among these patches [9]. Furthermore, the MSPA approach has shown the advantages of small data requirements and the convenient visual outcomes of this analysis. Therefore, as an important method of connecting ecological sources and to ensure species diffusion, ecological corridors are commonly simulated via various methods, e.g., graph theory, the minimum cumulative resistance (MCR) model, and circuit theory. In particular, the MCR model is used to simulate the process of species diffusion overcoming resistance to quantify the potential movement pattern of organisms and to improve the scientific significance of ecological network construction, and is thus widely applied to determine ecological security patterns in territorial spatial planning [10].
The MSPA method can identify habitat patches that play a crucial role in landscape connectivity in the study area based on the high pixels of satellite images. Indeed, most of the current studies on ecological network construction based on MSPA have identified the ecological sources in combination with landscape connectivity assessments. For example, critical areas such as nature reserves and forests were selected as prospects for MSPA and then large core patches were selected as ecological sources [11]; the MSPA method was used to identify the ecological sources based on patch size and landscape connectivity [12], and the ecological sources were identified based on the importance and ecological sensitivity of ecosystem services combined with the results of MSPA [7]. It is noted that although patches with high habitat quality and connectivity were generated, these studies ignored the importance of biodiversity conservation functions in ecological sources for the survival and migration of organisms among ecological networks. Meanwhile, the viability of urban wildlife populations is closely related to the quality of urban open spaces [13]. Therefore, in our study, we combined the BCF assessment with MSPA to identify the critical ecological sources for promoting biological migration and building a green ecological network to ultimately improve urban biodiversity in Qingdao City, Shandong, China.
Qingdao is one of the most livable cities in eastern China. The natural environment in the city is well known for its unique picturesque scenes throughout the entire urban area, harboring a variety of rare animals and plants. However, as a high-density area, numerous human disturbances, such as continuous expansion of urbanization and industrial production activities, have led to dramatic changes in urban green spaces [14], increased habitat fragmentation, and serious threats to biodiversity in urban parks [15]. In recent years, Qingdao has gradually recognized the importance of the construction of an ecological civilization, and the area of urban green space has been significantly increased. According to the statistical yearbook of Qingdao [16,17], the green coverage rate of built-up areas has been increased from 39.44% in 2016 to 43% in 2022. However, an effective urban green space system has not been established in the city and the lack of ecological corridors is not conducive to animal migration and diffusion. Therefore, the aims of our study were to improve the biodiversity and landscape sustainability of urban green spaces and to build an efficient urban green space network system, ultimately providing scientific evidence to support the construction of urban biodiversity protection and sustainable developmental strategies during the planning of urban green spaces.

2. Materials and Methods

2.1. Study Area

The study area was located in the urban areas of Qingdao, covering a total of 5 districts, i.e., Shinan, Shibei, Licang, Laoshan, and Chengyang, mainly between 35°35′–37°09′ N and 119°30′–121°00′ E (Figure 1). The total study area was about 827.70 km2 with the highest altitude of 1080 m, showing evident maritime characteristics, with an annual precipitation of 600–800 mm [18]. Although the green space coverage rate of the study area was high (about 36%), the green spaces were extremely fragmented in urban areas, with low biodiversity and poor ecological functions. In particular, the Laoshan District, located in the east of the study area, was rich in geological and geomorphic variations with high biodiversity, representing a remarkable ecological environment and a significantly important biological source in Qingdao [19].

2.2. Data Collection and Preprocessing

The satellite Landsat 8 OLI_TIRS image data (spatial resolution 30 m) and elevation data presented in the digital elevation model (DEM) were obtained from Geospatial Data Cloud (https://www.gscloud.cn/, accessed on 7 June 2022). Data of net primary productivity (NPP) were retrieved from the National Earth System Science Data Center (http://www.geodata.cn/, accessed on 9 November 2022). Data of rainfall and air temperature were collected from the Resources and Environmental Science Data Registration and Publication System (https://www.resdc.cn/, accessed on 25 August 2022).
ENVI 5.3 was used to divide the satellite image data of the study area into five categories: farm land, green land, water body, construction land, and unused land. Based on the highly accurate images taken during the same period, a land use and land cover (LULC) classification map (Figure 2) was obtained with a grid size of 30 m × 30 m, with the overall accuracy of the final interpretation being 96% and with a Kappa coefficient of 0.9.

2.3. Identification of Ecological Sources

First, ENVI 5.3 software was used to splice and cut the satellite image data. All green spaces in the study area were obtained via supervised classification and visual interpretation. Then, the MSPA was performed to divide all green patches into seven types of landscape, including core, island, loop, bridge, perforation, edge, and branch, based on spatial arrangement characteristics [20]. The ecological significance of each landscape type is shown in Table 1. The core patches were generally large habitat patches and were considered to be ecological source patches [21].
The patch connectivity index (PC) was used to identify the most important habitat patches and evaluate the fragmentation of the landscape and the connectivity of patches [22]. The patch importance value (dI) was calculated to evaluate the contribution of a single patch to the landscape connectivity index [23]. In this study, both PC and dI were calculated to evaluate the landscape connectivity of core patches using conefor 2.6 software, with the distance threshold set to 2500 m and the connectivity probability set to 0.5, as previously reported [24,25]. The PC and dI were calculated using the following equations, respectively:
P C = i = 1 n j = 1 n a i a j p i j * A L 2
d I = P C P C r e m o v e P C × 100
In Equation (1), n is the number of patches; ai and aj are the areas of patches i and j, respectively; AL is the total landscape area; and pij* is the maximum probability of species diffusion in patches i and j. The PC value is [0, 1], with the greater value representing a higher connectivity of the landscape. In Equation (2), dI is the connectivity index value of a certain landscape and PCremove is the connectivity index after the patch was removed. The greater value of dI indicates the greater importance of landscape elements.
The values of NPP, precipitation, temperature, and altitude indicate the capability of habitat patches to maintain biodiversity [26]. Therefore, the BCF index was calculated to select high-quality ecological patches using the following equations:
N P P m e a n = N P P 2010 S 10
P m e a n = P 2010 s 10
T m e a n = T 2010 s 10
B C F = N P P m e a n × P m e a n × T m e a n × 1 A
where NPPmean is the average annual vegetation NPP, Pmean is the average annual precipitation, Tmean is the average annual temperature, A is the altitude, and 2010s represents the sum of the annual data between 2010 and 2020.
The results of landscape connectivity evaluation and BCF evaluation of the core patches were normalized and overlayed, respectively, to obtain the total patch comprehensive evaluation value. The patches with high comprehensive evaluation values were selected as the ecological sources for the subsequent analyses.

2.4. Construction of Ecological Corridors

Species migration between patches is often hindered by different obstacles. A smaller landscape resistance indicates that fewer obstacles are encountered in the process of species migration and the higher efficiency of ecological flow transmission. Therefore, determination of the landscape resistance value is a pivotal step in extracting ecological corridors. The landscape resistance value is generally determined by the type of LULC, elevation, slope, distance to water body, road, and residential site [27,28]. In this study, ArcGIS was used to overlay the landscape resistance factors to obtain the landscape resistance value in the study area (Table 2), with each resistance factor weighted using the expert scoring method [29].
The MCR model was used to compare the cumulative resistance of grid pixels between patches to obtain low-resistance paths for biological migration and diffusion [30]. Based on the landscape resistance values obtained, the Cost Path tool of ArcGIS was used to calculate the raster pixel resistance between two patches and to obtain the minimum cumulative resistance path, i.e., an optimal ecological corridor. The MCR model was established as the following equation:
M C R = f m i n j = n i = m D i j × R i ,
where MCR represents the minimum cumulative cost of species from one ecological source to another patch in space, Dij represents the diffusion distance of species from patches i to j, and Ri is the landscape resistance value of patch i.

2.5. Ecological Network Evaluation and Optimization

Green ecological networks generally contain different structures due to various factors such as the length, location, and connection mode of ecological corridors. Network closure index (α), network connectivity index (β), and network connectivity rate index (γ) are usually calculated to evaluate the complexity of green ecological network structures [31] using the following equations:
α = l υ + 1 2 υ 5
β = l υ
γ = l 3 υ 2
where l is the number of ecological corridors and v is the number of ecological nodes. In this study, the intersections of the maximum- and minimum-cost paths were identified as ecological nodes [32]. The value of α is used to describe the loop level of an ecological network, i.e., a larger α value suggests more loops of the ecological network and higher mobility of matter and energy. The index β is defined as the average number of connections of each node and was used to measure the access to the ecological network, and a larger β value indicated a more complex ecological network structure. The index γ represents the interconnection of each node in the ecological network, i.e., a larger γ value implicates higher interconnection.
Based on the ecological network evaluation, the relationships between patches were quantitatively analyzed according to the interaction matrix between patches generated by the gravity model, which was calculated using the following equation:
G a b = L m a x 2 × l n S a × l n S b L a b 2 × P a × P b
where Gab (gravity value) is the relative significance between patches a and b; Sa and Sb are the areas of patches a and b, respectively; Pa and Pb are the average landscape resistance values of patch a and patch b; Lab is the cumulative landscape resistance value of the corridor between patches a and b; and Lmax is the maximum cumulative landscape resistance value of all corridors.
An ecological corridor generally contains several ecological nodes, and the connectivity of an ecological corridor can be judged by the density of its ecological nodes. Therefore, in this study, the most important ecological corridors were selected based on the interaction matrix, and kernel density (N) analysis was used to analyze the density of the most important ecological corridors. In this study, the Ns of ecological corridors were divided into five groups, i.e., high kernel density (N ≥ 723), relatively high density (468 ≤ N < 723), medium density (275 ≤ N < 468), relatively low density (109 ≤ N < 275), and low density (N < 109). All the ecological nodes in the region with high density were identified as key nodes. Finally, the green ecological network composed of important ecological corridors and key ecological nodes was identified as the optimized ecological network.

2.6. Determination of Ecological Corridor Width

When the width of an ecological corridor reaches the fundamental requirements, organisms can efficiently pass through the corridor and populations can rapidly increase [33]. In general, a corridor width between 30 and 60 m can meet the requirements of most animal migrations, while a width of 60–100 m is suitable for bird migration [24]. The existence of high-quality patches can improve the efficiency of ecological processes [34]. However, due to the expansion of urban construction, it may be a challenge to increase the width of the ecological corridors [35]. Therefore, 30, 60, 120, 300, and 600 m widths were used to identify the optimal ecological corridor width based on the proportions of different types of LULC.

3. Results

3.1. Identification of Ecological Sources

Using the MSPA method, the distribution of different landscape types was obtained in the study area (Figure 3). The results showed that the total area of green land in the study area was 432.82 km2, with a core area of about 329.73 km2, accounting for 76.18% of the total area. The largest core patch, i.e., the most important biological source, was the forest in Laoshan Mountain in the east of Qingdao, which contained abundant vegetation, with favorable growth conditions and rich biodiversity. The distribution of the Laoshan patch was concentrated with other source patches of built-up areas, which could meet the needs of species migration and energy exchange. These results were consistent with those previously reported [36]. There were few core patches observed in the western urban area, showing relatively scattered distribution, poor landscape connectivity, low efficiency of ecological processes between patches, and hindered ecological flow transmission.
The results of MSPA showed that the average core patch area was 0.14 km² (Figure 3). A total of 50 core patches with an area greater than 0.14 km2 were selected for further extraction. Then, based on the values of dI and BCF (Table 3), a total of 38 patches were identified as the ecological sources, each with a normalized total value greater than −0.5 (Figure 4). These ecological sources were mainly distributed in the east and south of the study area, including Laoshan Mountain, Fushan Mountain, Shiwan Forest Park, and other important ecological protection areas. Patch 34 Laoshan was revealed to have the highest dI and BCF, indicating that this patch contained the most important biological sources.
A large range of variations and generally low values were revealed in the ecological patch dI, ranging from 99.57 (patch 34) to 0.01 (patch 38) with an average of 6.57 (Table 3). In addition to patch 34 Laoshan, high dI values (>11) were revealed in patches 21, 10, and 20, which were mainly located in the central and southern parts of the study area, indicating the lack of high-landscape-connectivity patches in the west of the study area.
The results of the BCF assessment revealed that the areas of higher quality were mainly distributed in the eastern and southern regions, while the densely populated and built-up areas within the city, mainly including Fushan Ecological Park and Zhongshan Park, generally lacked high-quality green patches (Figure 5).

3.2. Construction of Ecological Corridors

Based on the six types of resistance factors (Figure 6) and the MCR model, the landscape resistance was obtained (Figure 7). The results showed that the overall landscape resistance values in the study area were high, with evident variations observed between the east and west regions of the study area. The eastern Laoshan District was revealed to have the smallest resistance value, while the landscape resistance values in the central urban area with a high density of construction land were larger, which was not conducive to material exchange and energy flow. Based on the MCR model, a total of 703 ecological corridors with a total length of 935.29 km were initially generated among 38 ecological sources, and these ecological corridors constituted a potential green ecological network (Figure 8). The distribution of ecological corridors in the central and southern parts was relatively dense, indicating that it is beneficial to build ecological networks with more green patches inside the central urban area.

3.3. Evaluation and Optimization of the Potential Green Ecological Network

Based on the potential green ecological network (Figure 8), a total of 284 potential ecological nodes were identified (Figure 9), with the values of the network connectivity indices α, β, and γ being 0.75, 2.48, and 0.83, respectively. These results showed that there were numerous closed loops in the potential green ecological network, the material cycle and energy flow were flexible, and the anti-interference ability was strong, which is conducive to further optimization.
As the most critical and vulnerable locations of a potential green ecological network, most of these ecological nodes were distributed in green land, while a small number of nodes were distributed in farm land and construction land, which enhanced the node connection function, ensuring the integrity of the entire potential green ecological network system.
Then, based on the gravity model, the ecological corridors with high gravity values were extracted as important corridors (Figure 10). The results showed that the gravity values of the potential ecological corridors were generally small with large variations, ranging from 59.49 to 887,457.89 with an average of 4913.66 and median of 715.24. The gravity values between each of the three patches 16, 29, and 30 and other patches were relatively high. In particular, patch 16 was surrounded by multiple patches, and the distances between patches were largely the same, which could facilitate effective connections among the ecological sources of various regions. Similarly, both patches 29 and 30 were shown to be important for maintaining the overall landscape connectivity in the northwest area.
There were 223 ecological corridors with gravity values greater than 1500, accounting for about 32% of the total number of potential ecological corridors. Therefore, the ecological corridors with gravity values greater than 1500 were identified as important corridors, and the rest were regarded as general corridors. The most important corridors were concentrated in the central urban and suburban transitional areas, indicating that the interactions between urban and suburban patches were important for maintaining the overall function of a green ecological network system.
Kernel density analysis was performed based on the important corridors to identify a total of 61 key nodes (Figure 11 and Figure 12). Among them, 26 key nodes were located at the edges of the ecological source patches, indicating that the ecological source patches were readily affected by the surrounding environment. These key nodes were concentrated in the northwest and south of the study area, which was similar to the areas of patches 16, 29, and 30, again indicating the importance of these ecological source patches in the green ecological network. Another 24 key nodes were located around the roads in a north–south linear distribution, indicating that the green corridors along the roads could provide appropriate living environmental resources in the process of biological migration. The other 11 key nodes were located along the rivers. Adequate water resources around these key nodes could not only provide appropriate habitats and nutritional supplies, but also introduce certain obstacles to biological migration. Therefore, biological corridors should be built at key nodes to reduce the resistance of creatures crossing rivers.
Based on the 223 important corridors and 61 key nodes, the values of network connectivity indices α, β, and γ of the optimized green ecological network were increased by 87%, 48%, and 52% to 1.39, 3.66, and 1.26, respectively, compared with those of a potential green ecological network without optimization. These results indicated that the structure of the optimized green ecological network was more stable with a stronger ability to resist external interference, and the efficiency of ecological flow was higher, which could significantly promote the migration of creatures between ecological sources.

3.4. Determination of Ecological Corridor Width

In our study, 30, 60, 120, 300, and 600 m widths were used to identify the appropriate width of ecological corridors based on the proportions of different types of LULC (Figure 13). The results showed that as the corridor width increased, the proportion of construction land increased, whereas the proportion of green land gradually declined. At 30 m, forest land accounted for the largest portion, and a largely optimal overall land type was obtained. Meanwhile, human interference was minimal, which could generally ensure the survival and migration of various organisms in the corridors. Therefore, it was determined that 30 m was the appropriate width of ecological corridors in our study area.
At the width of 30 m, analysis of landscape types showed that the core area accounted for 58.61% of the total area of the ecological corridors (Table 4). In the construction of ecological corridors, some of the landscape elements, such as unused land, could be transformed into green land, which could increase the size of the core area and accelerate the ecological process. The edge is the transition area between core and background, which accounted for about 22.08% of ecological corridors, indicating that the ecological corridors were susceptible to external interference. Therefore, it was recommended that the boundary of the edge should be strictly protected in order to reduce human interference.

4. Discussion

4.1. Identification of Ecological Sources

Although numerous studies have applied the MSPA method to construct green ecological networks [37,38], there is a lack of ecological source identification methods based on biodiversity [39]. Previous studies have graded ecological protection or drawn a red line based on the importance of biodiversity, but have ignored the dynamic process of biological flow [40]. Biodiversity plays a vital role in maintaining the proper functions and sustainability of ecosystems [41], and the BCF assessment can indirectly reveal the regional biodiversity status by evaluating the patch habitat quality [42]. Given these concerns and the very dense environment in Qingdao, our study further combined the MSPA method and BCF assessment to optimize the ecological source identification method and then construct a green ecological network. The results of ecological source patches were largely consistent with the current ecological resources in the study area [36], indicating the appropriate use of the comprehensive identification method applied in our study.
A total of 38 ecological sources were selected and largely coexisted with the national parks listed in the Qingdao “13th Five-Year” Urban Construction Plan [43], which provided the theoretical basis for the subsequent urban green space system planning. Compared with the 26 ecological patches identified in the green ecological network in Xiamen [44], more ecological sources were selected in our study, showing higher dI values of ecological sources, probably due to the high fragmentation of green land, small patch area, and uneven distribution in Qingdao. Similarly, 10 ecological sources were selected based on patch area and habitat quality in Beijing [45], showing that although the number of ecological sources was smaller, the spatial distribution of ecological sources was similar between Beijing and Qingdao, wherein they were mainly detected in mountainous areas.
Furthermore, the source patches overlapped with the areas of high-BCF spaces, which were mainly distributed in the eastern nature reserve and the southern green park of the study area. This was largely due to these areas being located near the coastline, which mainly had the land-use type of green land. Other factors, such as higher elevation, greater relief of terrain, mild climate, and more precipitation than the interior part of the peninsula, made these patches with high BCF suitable for the survival and reproduction of creatures, although strong BCF was observed only in a few large parks along the coast. However, the dI results revealed the weak connectivity between these high-quality patches and other patches. Therefore, it is important to increase the areas of green land in order to improve the quality and connectivity among ecological patches [46].

4.2. Distribution of Landscape Resistance

Qingdao is a typical peninsula city with a high-density urban area. Our results showed that the areas with high landscape resistance values were concentrated in roads and construction land, while the old city at the top of the peninsula was the most dense with the largest landscape resistance (Figure 6 and Figure 7). The results were consistent with those previously reported [47,48], showing the spatial distribution of landscape resistance in highly urbanized areas [49]. The large areas of parks such as Zhongshan Park and Fushan Park showed high BCF, but the landscape resistance between these patches was high, making it difficult for the animals in the Laoshan area to migrate to the inner green land of the city. The ecological corridors constructed in this study formed a dense ecological network, which was conducive to introducing the rich ecological resources of Laoshan Mountain into the existing high-quality green space in the city and alleviating the lack of ecological resources in the urban area of Qingdao.

4.3. Characteristics of Ecological Corridors

The MCR model was used to construct a potential green ecological network in the urban areas of Qingdao, which was composed of 703 ecological corridors and 284 ecological nodes. These results were consistent with the Qingdao Park Urban Construction Plan (2021–2035) [50], which identified Laoshan Mountain as the ecological core and proposed the spatial pattern of a “two belts and multiple corridors into a network” model. Comparing the potential green ecological network constructed in our study with the other high-density urban areas, it was found that the number of ecological corridors and ecological nodes in each region were quite different. For example, although the area of the Chengdu–Chongqing region is about 200 times bigger than that of our study area, the ecological network of this region contains only 17 ecological sources, 33 ecological corridors, and 29 ecological nodes [51]. This was probably because the distance between the sources is relatively long in the Chengdu–Chongqing region. Meanwhile, the ecological network of Guangzhou is similar to our study area, containing 35 ecological sources, 59 ecological corridors, and 39 ecological nodes, with a total length of 817.4 km [52]. Compared with the potential green ecological network in our study, the number of ecological sources and the total length of ecological corridors was similar to that in Guangzhou, but the numbers of ecological corridors and ecological nodes were relatively small. This was likely because the ecological corridors and ecological nodes in Guangzhou were extracted based on circuit theory, with the high-flow corridors and important nodes in the ecological process selected, whereas the MCR model was applied in our study to simulate the possible paths of biological migration, which could effectively connect the fragmented ecological patches.

4.4. Optimization of Potential Green Ecological Network

By optimizing the potential green ecological network, the priority areas for protection and construction could be accurately identified. The results showed that the values of the α, β, and γ indices (1.39, 3.66, and 1.26) of the optimized green ecological network in our study were higher than those (0.27, 2.03, and 0.72) of the green ecological network in Fuzhou [53], indicating that the green ecological network optimized by the gravity model and kernel density was more stable with higher ecological process efficiency. However, higher values of the α, β, and γ indices (3.206, 6.412, and 2.422) were detected in the optimized ecological network in the Dianchi Lake region [54], which treated the ecological sources as ecological nodes in the calculation of the network structure.
Previous studies identified an ecological corridor width of 100 m simulated in the Yanqing District of Beijing [55]. In our study, a shorter optimal ecological corridor width of 30 m was determined, probably due to the large proportion of construction land and the small number of forest patches in Qingdao. Therefore, it was recommended that the ecological corridors should be widened to the highest level possible, and the high-quality forest patches should be included in the boundary of the ecological corridors in order to improve the efficiency of biological migration [56,57]. In ecological corridors, the edge area can be easily affected by human activities. Compared with the internal species that require a single environmental resource condition, ecological corridors are more conducive to adapting to the migration of marginal species in a changing environment [58].

4.5. Suggestions for Ecological Construction of Qingdao

In Qingdao, it was recommended that the new urban areas should make full use of unused land and other patches with low ecological benefits by transforming them into green land patches, which will be conducive to the sustainable development of the city. However, the characteristics of land were well solidified in the old town area, making it largely impossible to add a large area of green space. Therefore, it is necessary to improve the quality of existing green space patches and to optimize the community structure. Furthermore, food-producing plants could provide sufficient nutrients for organisms [59]. To summarize, during ecological construction, native plants should be included on the basis of existing vegetation to provide adequate food and habitats for organisms in order to enhance the potential attractiveness of urban green patches to suburban organisms [60].
Furthermore, the vegetation in Qingdao is mainly characterized by deciduous broad-leaved forest, and organisms lack food and habitats in both autumn and winter [61]. Therefore, it is recommended that during the planning of ecological corridors, the proportion of evergreen tree species in the community should be increased in order to reduce the seasonal impact on the quality of ecological corridors [62], to improve the attraction of ecological corridors to organisms, and to increase the possibility and success rate of biological migration in the corridors.

5. Conclusions

In this study, ecological sources were identified based on the combination of MSPA and BCF assessment and the green ecological network in the urban area of Qingdao was constructed based on the MCR model to explore the spatial distribution of habitat quality and ecological resistance. The ecological network was further optimized based on the gravity model and kernel density analysis to clarify the spatial distribution and conservation importance of the green ecological network. Our study indicated that the construction of the green ecological network played an important role in promoting the landscape connectivity of ecological patches between urban and suburban areas and maintaining the urban biodiversity in Qingdao.

Author Contributions

Methodology, L.T.; Software, L.T.; Formal analysis, L.T.; Investigation, Y.C.; Resources, Y.C.; Writing—original draft, L.T.; Writing—review & editing, F.C. and H.L.; Supervision, F.C. and H.L. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the Qingdao Science and Technology Foundation for Public Wellbeing. Grant Number: 23-2-8-cspz-10-nsh.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Knapp, S.; Aronson, M.F.J.; Carpenter, E.; Herrera-Montes, A.; Jung, K.; Kotze, D.J.; La Sorte, F.A.; Lepczyk, C.A.; MacGregor-Fors, I.; MacIvor, J.S.; et al. A Research Agenda for Urban Biodiversity in the Global Extinction Crisis. Bioscience 2021, 71, 268–279. [Google Scholar] [CrossRef]
  2. Nor, A.N.M.; Corstanje, R.; Harris, J.A.; Grafius, D.R.; Siriwardena, G.M. Ecological connectivity networks in rapidly expanding cities. Heliyon 2017, 3, e325. [Google Scholar] [CrossRef]
  3. Huang, Z.; Zhou, L.; Chi, Y. Spring phenology rather than climate dominates the trends in peak of growing season in the Northern Hemisphere. Glob. Chang. Biol. 2023, 29, 4543–4555. [Google Scholar] [CrossRef] [PubMed]
  4. Beninde, J.; Veith, M.; Hochkirch, A. Biodiversity in cities needs space: A meta-analysis of factors determining intra-urban biodiversity variation. Ecol. Lett. 2015, 18, 581–592. [Google Scholar] [CrossRef] [PubMed]
  5. Aronson, M.F.; Lepczyk, C.A.; Evans, K.L.; Goddard, M.A.; Lerman, S.B.; MacIvor, J.S.; Nilon, C.H.; Vargo, T. Biodiversity in the city: Key challenges for urban green space management. Front. Ecol. Environ. 2017, 15, 189–196. [Google Scholar] [CrossRef]
  6. Belaire, J.A.; Higgins, C.; Zoll, D.; Lieberknecht, K.; Bixler, R.P.; Neff, J.L.; Keitt, T.H.; Jha, S. Fine-scale monitoring and mapping of biodiversity and ecosystem services reveals multiple synergies and few tradeoffs in urban green space management. Sci. Total Environ. 2022, 849, 157801. [Google Scholar] [CrossRef]
  7. Chen, X.; Peng, J.; Liu, Y.; Yang, Y.; Li, G. Construction of ecological security pattern of Yunfu City based on the framework of “importance-sensitivity-connectivity”. Geogr. Res. 2017, 36, 471–484. [Google Scholar]
  8. Yeo, O.T.S.; Mohd Yusof, M.J.; Maruthaveeran, S.; Mohd Shafri, H.Z.; Saito, K.; Yeo, L.B. ABC of green infrastructure analysis and planning: The basic ideas and methodological guidance based on landscape ecological principle. Urban For. Urban Green. 2022, 73, 127600. [Google Scholar] [CrossRef]
  9. Modica, G.; Praticò, S.; Laudari, L.; Ledda, A.; Di Fazio, S.; De Montis, A. Implementation of multispecies ecological networks at the regional scale: Analysis and multi-temporal assessment. J. Environ. Manag. 2021, 289, 112494. [Google Scholar] [CrossRef]
  10. Wang, J.; Li, Y.; Wang, S.; Li, Q.; Li, L.; Liu, X. Assessment of Multiple Ecosystem Services and Ecological Security Pattern in Shanxi Province, China. Int. J. Environ. Res. Public Health 2023, 20, 4819. [Google Scholar] [CrossRef]
  11. Zhang, Q.; Li, F.; Wang, D.; Li, M.; Chen, D. Analysis on changes of ecological spatial connectivity in Jiangsu Province based on ecological network. Acta Ecol. Sin. 2021, 41, 3007–3020. [Google Scholar]
  12. Wang, Y.; Shen, C.; Jin, X.; Bao, G.; Liu, J.; Zhou, Y. Developing and optimizing ecological networks based on MSPA and MCR model. Ecol. Sci. 2019, 38, 138–145. [Google Scholar]
  13. Penteado, H.M. Effects of open space configurations and development patterns on future urban wildlife habitats and populations. City Environ. Interact. 2023, 19, 100106. [Google Scholar] [CrossRef]
  14. Yu, Q.; Xie, L.; Dong, L.; Cheng, N.; Liu, Y. A conceptual framework for riparian green space planning in mountainous cities. Energy Rep. 2023, 9, 499–513. [Google Scholar] [CrossRef]
  15. Yang, J.; Li, X.; Li, S.; Liang, H.; Lu, H. The woody plant diversity and landscape pattern of fine-resolution urban forest along a distance gradient from points of interest in Qingdao. Ecol. Indic. 2021, 122, 107326. [Google Scholar] [CrossRef]
  16. National Bureau of Statistics of China. Qingdao Statistical Yearbook of 2016; Qingdao Municipal Bureau of Statistics NBS Survey Office: Qingdao, China, 2016. [Google Scholar]
  17. National Bureau of Statistics of China. Qingdao Statistical Yearbook of 2022; Qingdao Municipal Bureau of Statistics NBS Survey Office: Qingdao, China, 2022. [Google Scholar]
  18. Guo, X.; Hu, Y.; Ma, J.; Wang, H.; Wang, K.; Wang, T.; Jiang, S.; Jiao, J.; Sun, Y.; Jiang, X.; et al. Nitrogen Deposition Effects on Invasive and Native Plant Competition: Implications for Future Invasions. Ecotoxicol. Environ. Saf. 2023, 259, 115029. [Google Scholar] [CrossRef]
  19. Li, J.; Jia, K.; Zhang, N.; Wei, X.; Wang, B. Evaluation of ecological protection importance in Qingdao City based on remote sensing and ecological service model. Remote Sens. Technol. Appl. 2021, 36, 1329–1338. [Google Scholar]
  20. Soille, P.; Vogt, P. Morphological segmentation of binary patterns. Pattern Recognit. Lett. 2009, 30, 456–459. [Google Scholar] [CrossRef]
  21. Urban, D.; Keitt, T. Landscape Connectivity: A Graph-Theoretic Perspective. Ecology 2001, 82, 1205–1218. [Google Scholar] [CrossRef]
  22. Velázquez, J.; Gutiérrez, J.; García-Abril, A.; Hernando, A.; Aparicio, M.; Sánchez, B. Structural connectivity as an indicator of species richness and landscape diversity in Castilla y León (Spain). For. Ecol. Manag. 2019, 432, 286–297. [Google Scholar] [CrossRef]
  23. Gao, Y.; Mu, H.; Zhang, Y.; Tian, Y.; David, T.; Li, X. Research on construction path optimization of urban-scale green network system based on MSPA analysis method: Taking Zhaoyuan City as an example. Acta Ecol. Sin. 2019, 39, 7547–7556. [Google Scholar]
  24. Yang, Z.; Jiang, Z.; Guo, C.; Yang, X.; Xu, X.; Li, X.; Hu, Z.; Zhou, H. Construction of ecological network using morphological spatial pattern analysis and minimal cumulative resistance models in Guangzhou City, China. Chin. J. Appl. Ecol. 2018, 29, 3367–3376. [Google Scholar]
  25. Rincón, V.; Velázquez, J.; Pascual, Á.; Herráez, F.; Gómez, I.; Gutiérrez, J.; Sánchez, B.; Hernando, A.; Santamaría, T.; Sánchez-Mata, D. Connectivity of Natura 2000 potential natural riparian habitats under climate change in the Northwest Iberian Peninsula: Implications for their conservation. Biodivers. Conserv. 2022, 31, 585–612. [Google Scholar] [CrossRef]
  26. Yang, W.; Liu, J.; Wang, W.; Liu, X.; Hao, X. Biodiversity monitoring and assessment based on multispectral remote sensing data. Spectrosc. Spect. Anal. 2023, 43, 1282–1290. [Google Scholar]
  27. Milanesi, P.; Holderegger, R.; Bollmann, K.; Gugerli, F.; Zellweger, F. Three-dimensional habitat structure and landscape genetics: A step forward in estimating functional connectivity. Ecology 2017, 98, 393–402. [Google Scholar] [CrossRef]
  28. Betbeder, J.; Laslier, M.; Hubert-Moy, L.; Burel, F.; Baudry, J. Synthetic Aperture Radar (SAR) images improve habitat suitability models. Landsc. Ecol. 2017, 32, 1867–1879. [Google Scholar] [CrossRef]
  29. Huang, Y.; Huang, J.; Liao, T.; Liang, X.; Tian, H. Simulating urban expansion and its impact on functional connectivity in the Three Gorges Reservoir Area. Sci. Total Environ. 2018, 643, 1553–1561. [Google Scholar] [CrossRef]
  30. Lopes, I.J.C.; Biondi, D.; Corte, A.P.D.; Reis, A.R.N.; Oliveira, T.G.S. A methodological framework to create an urban greenway network promoting avian connectivity: A case study of Curitiba City. Urban For. Urban Green. 2023, 87, 128050. [Google Scholar] [CrossRef]
  31. Urban, D.L.; Minor, E.S.; Treml, E.A.; Schick, R.S. Graph models of habitat mosaics. Ecol. Lett. 2009, 12, 260–273. [Google Scholar] [CrossRef]
  32. Liu, X.; Yang, C.; Gao, Y.; Sun, Q.; Wang, S.; Feng, Z.; Liu, X. Identification of ecological security pattern and hierarchical management of ecological control area in Xiamen. Acta Ecol. Sin. 2023, 43, 5357–5369. [Google Scholar]
  33. Zhu, Q.; Yu, K.; Li, D. The width of ecological corridor in landscape planning. Acta Ecol. Sin. 2005, 25, 2406–2412. [Google Scholar]
  34. Díaz-Forestier, J.; Abades, S.; Pohl, N.; Barbosa, O.; Godoy, K.; Svensson, G.L.; Undurraga, M.I.; Bravo, C.; García, C.; Root-Bernstein, M.; et al. Assessing Ecological Indicators for Remnant Vegetation Strips as Functional Biological Corridors in Chilean Vineyards. Diversity 2021, 13, 447. [Google Scholar] [CrossRef]
  35. Chi, Y.; Shi, H.; Feng, A. Typical island landscape ecological network establishmen—A case study of Chongming Island. Mar. Eviromental Sci. 2015, 34, 433–440. [Google Scholar]
  36. Qian, W.; Zhao, Y.; Li, X. Construction of ecological security pattern in coastal urban areas: A case study in Qingdao, China. Ecol. Indic. 2023, 154, 110754. [Google Scholar] [CrossRef]
  37. Staccione, A.; Candiago, S.; Mysiak, J. Mapping a Green Infrastructure Network: A framework for spatial connectivity applied in Northern Italy. Environ. Sci. Policy 2022, 131, 57–67. [Google Scholar] [CrossRef]
  38. Lattanzio, M.S.; Miles, D.B. Ecological divergence among colour morphs mediated by changes in spatial network structure associated with disturbance. J. Anim. Ecol. 2014, 83, 1490–1500. [Google Scholar] [CrossRef]
  39. Van Teeffelen, A.J.A.; Vos, C.C.; Opdam, P. Species in a dynamic world: Consequences of habitat network dynamics on conservation planning. Biol. Conserv. 2012, 153, 239–253. [Google Scholar] [CrossRef]
  40. Karimi, A.; Yazdandad, H.; Reside, A.E. Spatial conservation prioritization for locating protected area gaps in Iran. Biol. Conserv. 2023, 279, 109902. [Google Scholar] [CrossRef]
  41. Hector, A.; Joshi, J.; Lawler, S.; Spehn, E.M.; Wilby, A. Conservation implications of the link between biodiversity and ecosystem functioning. Oecologia 2001, 129, 624–628. [Google Scholar] [CrossRef]
  42. Wang, R.; Gamon, J.A. Remote sensing of terrestrial plant biodiversity. Remote Sens. Environ. 2019, 231, 111218. [Google Scholar] [CrossRef]
  43. Qingdao “13th Five-Year” Urban Construction Plan; Qingdao Urban and rural Construction Committee: Qingdao, China, 2016.
  44. Liu, X.; Zeng, J.; Zeng, P. Construction and Optimization of the Green Space. Chin. Landsc. Archit. 2019, 36, 76–81. [Google Scholar]
  45. Dong, J.; Peng, J.; Liu, Y.; Qiu, S.; Han, Y. Integrating spatial continuous wavelet transform and kernel density estimation to identify ecological corridors in megacities. Landsc. Urban Plan. 2020, 199, 103815. [Google Scholar] [CrossRef]
  46. Beita, C.M.; Murillo, L.F.S.; Alvarado, L.D.A. Ecological corridors in Costa Rica: An evaluation applying landscape structure, fragmentation-connectivity process, and climate adaptation. Conserv. Sci. Pract. 2021, 3, e475. [Google Scholar] [CrossRef]
  47. Liu, Q.; Sun, Y.; Mei, Y.; Jian, Z.; Pan, F.; Zhang, L. Construction and Analysis of Ecological SecurityPattern of Qingdao Based on MSPAand MCR Models. Pol. J. Environ. Stud. 2023, 32, 155–169. [Google Scholar] [CrossRef] [PubMed]
  48. Kang, J.; Qing, Y.; Lu, W. Construction and optimization of the Saihanba ecological network. Ecol. Indic. 2023, 153, 110401. [Google Scholar] [CrossRef]
  49. Wang, Y.; Qu, Z.; Zhong, Q.; Zhang, Q.; Zhang, L.; Zhang, R.; Yi, Y.; Zhang, G.; Li, X.; Liu, J. Delimitation of ecological corridors in a highly urbanizing region based on circuit theory and MSPA. Ecol. Indic. 2022, 142, 109258. [Google Scholar] [CrossRef]
  50. Qingdao Park Urban Construction Plan (2021–2035) (Public notice); Qingdao Municipal Landscape Gardens and Forestry Bureau: Qingdao, China, 2023.
  51. Jiang, J.; Zhang, H.; Huang, Q.; Liu, F.; Li, L.; Qiu, H.; Zhou, S. Diagnosis of Key Ecological Restoration Areas in Territorial Space under the Guidance of Resilience: A Case Study of the Chengdu–Chongqing Region. Land 2023, 12, 973. [Google Scholar] [CrossRef]
  52. Qian, M.; Huang, Y.; Gao, Y.; Wu, J.; Xiong, Y. Ecological network construction and optimization in Guangzhou from the perspective of biodiversity conservation. J. Environ. Manag. 2023, 336, 117692. [Google Scholar] [CrossRef]
  53. Lin, J.; Yang, W.; Yu, K.; Geng, J.; Liu, J. Identification and Construction of Ecological Nodes in the Fuzhou Ecological Corridors. Forests 2022, 13, 1837. [Google Scholar] [CrossRef]
  54. Zhou, S.; Song, Y.; Li, Y.; Wang, J.; Zhang, L. Construction of Ecological Security Pattern for Plateau Lake Based on MSPA–MCR Model: A Case Study of Dianchi Lake Area. Sustainability 2022, 14, 14532. [Google Scholar] [CrossRef]
  55. Kong, Y.; Wang, S. Construction of urban and rural ecological network in Yanqing District of Beijing based on MSPA model. J. Beijing For. Univ. 2020, 42, 113–121. [Google Scholar]
  56. Tian, M.; Chen, X.; Gao, J.; Tian, Y. Identifying ecological corridors for the Chinese ecological conservation redline. PLoS ONE 2022, 17, e271076. [Google Scholar] [CrossRef]
  57. van Schalkwyk, J.; Pryke, J.S.; Samways, M.J.; Gaigher, R. Corridor width determines strength of edge influence on arthropods in conservation corridors. Landsc. Ecol. 2020, 35, 1175–1185. [Google Scholar] [CrossRef]
  58. Fried, J.H.; Levey, D.J.; Hogsette, J.A. Habitat corridors function as both drift fences and movement conduits for dispersing flies. Oecologia 2005, 143, 645–651. [Google Scholar] [CrossRef] [PubMed]
  59. Feng, M.; Dai, A.; Guo, S. Research on the Design Strategies of Wildlife Protection Facilities in Nanjing Zijin Mountain Forest Park. Landc. Plan. Des. 2020, 81–86. [Google Scholar]
  60. Li, X.; Zhang, S.; Huang, R.; Feng, L.; Xu, S.; Liu, B. Diversity and distribution variation of urban spontaneous vegetation with distinct frequencies along river corridors in a fast-growing city. J. Environ. Manag. 2023, 333, 117446. [Google Scholar] [CrossRef] [PubMed]
  61. Peng, J.; Li, Y.; Tian, L.; Liu, Y.; Wang, Y. Vegetation Dynamics and Associated Driving Forces in Eastern China during 1999–2008. Remote Sens. 2015, 7, 13641–13663. [Google Scholar] [CrossRef]
  62. Zhao, X.; Xu, X. Research on Landscape Ecological Security Pattern in a Eucalyptus Introduced Region Based on Biodiversity Conservation. Russ. J. Ecol. 2015, 46, 59–70. [Google Scholar] [CrossRef]
Figure 1. Study area.
Figure 1. Study area.
Sustainability 15 16579 g001
Figure 2. Classification and distribution of 5 types of land use and land cover (LULC) in the study area.
Figure 2. Classification and distribution of 5 types of land use and land cover (LULC) in the study area.
Sustainability 15 16579 g002
Figure 3. Distribution of core patches in the study area based on morphological spatial pattern analysis (MSPA).
Figure 3. Distribution of core patches in the study area based on morphological spatial pattern analysis (MSPA).
Sustainability 15 16579 g003
Figure 4. Distribution of ecological sources based on a total of 38 patches listed in Table 3, and the numbers in the figure represent patch number in Table 3.
Figure 4. Distribution of ecological sources based on a total of 38 patches listed in Table 3, and the numbers in the figure represent patch number in Table 3.
Sustainability 15 16579 g004
Figure 5. Biodiversity conservation function (BCF) assessment of the study area.
Figure 5. Biodiversity conservation function (BCF) assessment of the study area.
Sustainability 15 16579 g005
Figure 6. Distribution of 6 types of resistance factors in the study area.
Figure 6. Distribution of 6 types of resistance factors in the study area.
Sustainability 15 16579 g006
Figure 7. Landscape resistance based on 6 types of resistance factors.
Figure 7. Landscape resistance based on 6 types of resistance factors.
Sustainability 15 16579 g007
Figure 8. Potential green ecological network based on the minimum cumulative resistance (MCR) model.
Figure 8. Potential green ecological network based on the minimum cumulative resistance (MCR) model.
Sustainability 15 16579 g008
Figure 9. Potential ecological nodes based on the potential green ecological network.
Figure 9. Potential ecological nodes based on the potential green ecological network.
Sustainability 15 16579 g009
Figure 10. Optimized ecological corridors based on gravity model. Yellow stars represent Patches 16, 29, and 30.
Figure 10. Optimized ecological corridors based on gravity model. Yellow stars represent Patches 16, 29, and 30.
Sustainability 15 16579 g010
Figure 11. Optimized green ecological network distribution showing the 61 key nodes based on kernel density analysis.
Figure 11. Optimized green ecological network distribution showing the 61 key nodes based on kernel density analysis.
Sustainability 15 16579 g011
Figure 12. Distribution of the 61 key nodes based on kernel density analysis.
Figure 12. Distribution of the 61 key nodes based on kernel density analysis.
Sustainability 15 16579 g012
Figure 13. Proportion of 5 types of land use and land cover (LULC) under different ecological corridor widths.
Figure 13. Proportion of 5 types of land use and land cover (LULC) under different ecological corridor widths.
Sustainability 15 16579 g013
Table 1. MSPA landscape types and their ecological significance.
Table 1. MSPA landscape types and their ecological significance.
Type of LandscapeEcological Significance
CoreHabitat patches with larger areas, which provide larger habitats for species, are generally used as ecological source areas in the ecological network.
IslandLittle isolated patches that are not connected to each other; the possibility of internal material and energy exchange and transfer is lower.
LoopThe ecological corridor connecting the same core areas serves as a channel for energy exchange within the core areas.
BridgeThe ribbon area connecting the core areas plays a corridor role in the ecological network and is an important landscape type for enhancing landscape connectivity.
PerforationMarginal zones of inner core patches with edge effects.
EdgeThe external boundary of the core areas protects the internal core areas from external interference and ensures normal energy flow and material exchange in the core areas.
BranchA linear landscape with only one end connected to bridge, loop, edge, or perforation.
Table 2. Characteristics of landscape resistance factors in this study.
Table 2. Characteristics of landscape resistance factors in this study.
Landscape Resistance FactorContentClassification of Landscape ResistanceValueWeight
Land use typeLand use and land coverGreen land10.2
Water body2
Farm land3
Unused land4
Construction land5
Natural topographic factorAltitude (m)0–20010.15
200–4002
400–6003
600–8004
>8005
Slope (°)0–510.15
5–102
10–153
15–204
>205
Distance to water body (m)0–20010.1
200–4002
400–6003
600–8004
>8005
Human activity factorDistance to road (m)0–10010.2
100–2002
200–3003
300–4004
>4005
Distance to residential area (m)0–20010.2
200–4002
400–6003
600–8004
>8005
Table 3. A total of 38 core patches identified in the study area ranked by the values of dI and BCF.
Table 3. A total of 38 core patches identified in the study area ranked by the values of dI and BCF.
Rank of PatchPatch NumberBCFdINormalized Value of Both BCF and dI
13429,241.5999.5713.74
22116.1412.970.41
3105.3511.900.34
42017.6011.240.29
541.8210.370.23
6128.599.830.19
7155.708.950.13
8115.048.900.12
952.748.870.12
1010.748.800.11
1136819.634.390.00
12139.786.78−0.03
13291.255.92−0.09
1437668.272.20−0.19
1566.344.15−0.22
16265.404.05−0.22
17148.673.85−0.24
18312.593.19−0.28
1973.132.99−0.30
20810.112.93−0.30
21228.652.93−0.30
223538.482.45−0.33
2331120.992.07−0.34
24234.042.44−0.34
25197.242.26−0.35
2632100.170.83−0.43
27243.940.91−0.44
281823.050.71−0.46
2938221.980.01−0.46
303386.130.45−0.46
312423.990.37−0.48
322554.590.26−0.48
33305.450.41−0.48
341731.100.31−0.48
352731.640.31−0.49
361619.410.33−0.49
37284.440.36−0.49
38924.550.28−0.49
Average834.226.571.43
Table 4. Characteristics of landscape types with the ecological corridor width of 30 m.
Table 4. Characteristics of landscape types with the ecological corridor width of 30 m.
Landscape TypeArea (km2)Proportion of Ecological Corridor Area (%)
Core12.6758.61
Islet0.642.95
Perforation0.572.63
Edge4.7722.08
Loop0.160.76
Bridge1.165.38
Branch1.647.58
Total area of green land21.6179.36
Total area of ecological corridor27.23100
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

Tao, L.; Chen, Y.; Chen, F.; Li, H. Construction of Green Ecological Network in Qingdao (Shandong, China) Based on the Combination of Morphological Spatial Pattern Analysis and Biodiversity Conservation Function Assessment. Sustainability 2023, 15, 16579. https://doi.org/10.3390/su152416579

AMA Style

Tao L, Chen Y, Chen F, Li H. Construction of Green Ecological Network in Qingdao (Shandong, China) Based on the Combination of Morphological Spatial Pattern Analysis and Biodiversity Conservation Function Assessment. Sustainability. 2023; 15(24):16579. https://doi.org/10.3390/su152416579

Chicago/Turabian Style

Tao, Ling, Yanni Chen, Fang Chen, and Haifang Li. 2023. "Construction of Green Ecological Network in Qingdao (Shandong, China) Based on the Combination of Morphological Spatial Pattern Analysis and Biodiversity Conservation Function Assessment" Sustainability 15, no. 24: 16579. https://doi.org/10.3390/su152416579

APA Style

Tao, L., Chen, Y., Chen, F., & Li, H. (2023). Construction of Green Ecological Network in Qingdao (Shandong, China) Based on the Combination of Morphological Spatial Pattern Analysis and Biodiversity Conservation Function Assessment. Sustainability, 15(24), 16579. https://doi.org/10.3390/su152416579

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