Regional Ecological Security Pattern Construction Based on Ecological Barriers: A Case Study of the Bohai Bay Terrestrial Ecosystem

: The construction of ecological barriers and ecological security patterns is an important way of maintaining regional ecological security in landscape ecology. However, there is still no consensus on the concept and connotation of ecological barriers, and the zoning and adaptive management of ecological sources are rarely considered in the construction of ecological security patterns. This study uses the terrestrial ecosystem of Bohai Bay, China as a study area, and the identiﬁcation and zoning of ecological sources in the ecological security pattern are achieved by combining an ecosystem service assessment with an ecological risk assessment, and on this basis, ecological barriers are identiﬁed to optimize the structure and function of ecological sources. The minimum cumulative resistance model is used to identify ecological corridors and ecological strategic nodes and to construct an ecological security pattern based on the modiﬁed ecological sources. The results demonstrate that ﬁrstly, 2873.25 km 2 was identiﬁed as the ecological source, accounting for 14.28% of the total. Secondly, there are three large ecological barrier zones and nine ecological barrier cells with a total area of 1173.06 km 2 , accounting for 40.83% of the ecological sources. Thirdly, a total of 35 ecological corridors were extracted, and 32 ecological strategic nodes were marked, mainly distributed at the intersection and branches of important ecological corridors. An ecological security pattern construction system was formed with the collection of ecological source selection, ecological barrier identiﬁcation, ecological resistance surface construction, and ecological corridor extraction. Fourthly, the concept and connotation of ecological barriers was analyzed, and the complementary relationship between ecological barriers and ecological security patterns in terms of structure and function is discussed. This study enriches the deﬁnition and connotation of ecological barriers, provides a new framework for identifying the ecological security patterns, and provides scientiﬁc guidance for ecological protection and management in coastal areas. data J.Z.;


Introduction
As urbanization accelerates, the ecological services (ESs) that ecosystems can provide are gradually becoming more and more homogeneous, and ecological security is under serious threat [1]. Building ecological security patterns (ESPs) is an important way to maintain the structure and function of ecosystems, ensure ecological security, and balance the relationship between economic and social development and ecological protection [2][3][4].
The traditional methodology of constructing ESPs can be divided into three steps: identifying ecological sources, constructing ecological resistance surface, and determining elements, such as ecological corridors and ecological strategic nodes [5]. Ecological sources

Study Area
The study area is the Bohai Bay terrestrial ecosystem, which covers 14 coastal districts and counties belonging to Hebei, Shandong, and Tianjin, respectively. Bohai Bay is located at 117°30′-118°51′E, 38°0′-39°15′N, with a width of 111 km from east to west and a length of 130 km from north to south, covering an area of 19,500 km 2 . The continental monsoon climate is marked by cold winters and hot summers, with four distinct seasons. The area has many rivers, lakes, ponds, reservoirs, estuaries, and shallow mudflats, forming a rich and diverse wetland landscape. The unique geographical location and good wetland environment make it an important distribution area for wetland waterbirds in eastern China, an important feeding and breeding site on the migration route of migratory birds from East Asia to Australasia, and an important refueling station for migratory waterbirds to replenish their energy. Bohai Bay, as an important maritime corridor in the Bohai Economic Circle, is rapidly urbanizing, with exceptionally active human activities and huge pressure on the ecosystem. The study area is shown in Figure 1.

Data Sources
The land use data of Bohai Bay for 1990, 2000, 2010, and 2018 used in this study were obtained from the Resource and Environment Science and Data Centre of the Chinese Academy of Sciences (https://www.resdc.cn, 13 March 2021. The land use data for 1990, 2000, and 2010 are based on Landsat TM/ETM+ remote sensing images, and the land use data for 2018 are based on Landsat 8 remote sensing images, which were generated by manual visual interpretation. The spatial resolution of the data has a certain impact on the model [43]; in this study, the resolution of the raster data is all unified at 30 m and ensures consistent spatial reference. The land use types include 6 primary types of arable land, forest land, grassland, water, residential land, and unused land and 25 secondary types.
In the habitat quality model, the types of land use maps in 1990, 2000, 2010, and 2018 with frequent human activities-arable land, urban land, rural settlements, other construction land, and unused land-were extracted as threat source layers and ecological land-for-

Data Sources
The land use data of Bohai Bay for 1990, 2000, 2010, and 2018 used in this study were obtained from the Resource and Environment Science and Data Centre of the Chinese Academy of Sciences (https://www.resdc.cn, 13 March 2021. The land use data for 1990, 2000, and 2010 are based on Landsat TM/ETM+ remote sensing images, and the land use data for 2018 are based on Landsat 8 remote sensing images, which were generated by manual visual interpretation. The spatial resolution of the data has a certain impact on the model [43]; in this study, the resolution of the raster data is all unified at 30 m and ensures consistent spatial reference. The land use types include 6 primary types of arable land, forest land, grassland, water, residential land, and unused land and 25 secondary types. In the habitat quality model, the types of land use maps in 1990, 2000, 2010, and 2018 with frequent human activities-arable land, urban land, rural settlements, other construction land, and unused land-were extracted as threat source layers and ecological land-forests, grasslands, rivers, and lakes-were extracted as habitat layers. In the habitat risk assessment of the Boai Bay ecological sources, based on the frequent reclamation activities in the Bohai Bay area, we identified the Bhai Bay reclamation risk layer on the basis of the original threat layers and the etraction function of ArcGIS. Arable land, urban land, rural settlements, other construction land, unused land and reclamation areas were all etracted as the risk source layer; the reclamation areas were extracted as the risk source layer, and the terrestrial ecosystems with high habitat quality in Bohai Bay were taken as Sustainability 2022, 14, 5384 5 of 23 the ecological service supply area layer. The software packages used were ArcGIS 10.5 and InVEST 3.8.9.

Integrated Framework
This study devised an idea for delineating the scope of the regional ecological barriers and proposes a specific method for optimizing the ESPs. Firstly, the habitat quality of terrestrial ecosystems in Bohai Bay was assessed through the habitat quality module in the InVEST model, and areas with high habitat quality were selected as ecological sources. Using the habitat risk assessment model, the ecological sources were assessed for the risk of degradation caused by human activities. Ecological sources with high risk ratings were defined as areas where ecological barriers should be constructed. The ecological sources and ecological barriers were together referred to as ecological service supply areas. The ESPs of the Bohai Bay land area were then constructed based on the paradigm of "ecological sources-ecological barriers-resistance surface construction-ecological corridor extraction". This resulted in ESPs based on ESs and ecosystem risks. The research idea is shown in Figure 2 below. ests, grasslands, rivers, and lakes-were extracted as habitat layers. In the habitat risk assessment of the Boai Bay ecological sources, based on the frequent reclamation activities in the Bohai Bay area, we identified the Bhai Bay reclamation risk layer on the basis of the original threat layers and the etraction function of ArcGIS. Arable land, urban land, rural settlements, other construction land, unused land and reclamation areas were all etracted as the risk source layer; the reclamation areas were extracted as the risk source layer, and the terrestrial ecosystems with high habitat quality in Bohai Bay were taken as the ecological service supply area layer. The software packages used were ArcGIS 10.5 and InVEST 3.8.9.

Integrated Framework
This study devised an idea for delineating the scope of the regional ecological barriers and proposes a specific method for optimizing the ESPs. Firstly, the habitat quality of terrestrial ecosystems in Bohai Bay was assessed through the habitat quality module in the InVEST model, and areas with high habitat quality were selected as ecological sources. Using the habitat risk assessment model, the ecological sources were assessed for the risk of degradation caused by human activities. Ecological sources with high risk ratings were defined as areas where ecological barriers should be constructed. The ecological sources and ecological barriers were together referred to as ecological service supply areas. The ESPs of the Bohai Bay land area were then constructed based on the paradigm of "ecological sources-ecological barriers-resistance surface construction-ecological corridor extraction". This resulted in ESPs based on ESs and ecosystem risks. The research idea is shown in Figure 2 below.

Identifying Ecological Barriers
The Habitat Quality (HQ) model assesses the ESs of an ecosystem and is expressed as the magnitude of the Habitat Quality Index. Habitat quality is calculated as follows: where Qxj is the habitat quality of raster x in land use type j; H is the habitat suitability of land use type j; Dxj is the degree of habitat degradation of raster x in land use type j; z is the normalization constant and takes the value of 2.5, k is the half-saturation constant, and the default value is 0.5.

Identifying Ecological Barriers
The Habitat Quality (HQ) model assesses the ESs of an ecosystem and is expressed as the magnitude of the Habitat Quality Index. Habitat quality is calculated as follows: where Q xj is the habitat quality of raster x in land use type j; H is the habitat suitability of land use type j; D xj is the degree of habitat degradation of raster x in land use type j; z is the normalization constant and takes the value of 2.5, k is the half-saturation constant, and the default value is 0.5.
where R denotes the number of threat factors; r is the threat factor of the habitat; Y r denotes the number of rasters on the r threat factor layer; w r denotes the weight of the threat factor; r y denotes the level of influence of the threat factor on habitat y, taking values between 0 and 1; i rxy denotes the influence of the habitat with threat factor r in raster x on raster y; β x Sustainability 2022, 14, 5384 6 of 23 denotes the approachable level of raster x; S jr denotes the land use sensitivity of type j to threat factor r.
where d xy is the linear distance between raster x and y; d rmax is the maximum action distance of threat r.
In this paper, on the basis of the InVEST model manual and related studies [44,45], the threat factors and threat factor sensitivities were assigned with reference to the economic and social development of the study area, as shown in Tables 1 and 2.
The degree of risk to ecosystems due to human activities is another important basis for identifying ecological barriers in this study. The Habitat Risk Assessment (HRA) model can be used to assess the risk to ecosystems from human activities and defines risk as the potential for human activities to reduce habitat quality and thus damage ecosystem services [46].
In this study, two indicators, exposure and consequence, were selected for the ecological risk assessment of the terrestrial ecosystem in Bohai Bay, with reference to the index requirements of the HRA model and related studies [16,19,47]. Exposure refers to the extent to which the ecological service supply area is exposed to human activities and is divided Sustainability 2022, 14, 5384 7 of 23 into three indicators: the spatial superposition of the ecological service supply area and human activities, the frequency of human activities, and the intensity of disturbance. Consequence refers to the sensitivity of the ecological service supply area to human activities and includes three indicators: habitat heterogeneity, management effectiveness, and habitat loss. Habitat heterogeneity refers to the strength of ESs provided by the ecological service areas, as shown in Table 3. The specific calculation process for habitat risk is as follows: where R ij denotes the risk to habitat i caused by pressure j, E denotes exposure, and C denotes consequence. The scores for total exposure E and total consequence C are the weighted average of the exposure value e i and consequence value c i for each evaluation indicator i. Each evaluation indicator in this study is equally weighted.
where d i represents the data quality score for indicator i, w i represents the importance weight of indicator i, and N represents the number of indicators used to evaluate each habitat.
Finally, the model calculates the cumulative risk to the habitat for all pressures, with the cumulative risk to habitat i being the sum of the risk scores under the influence of each pressure. where R i represents the cumulative risk to habitat i, J represents the number of pressures, and R ij represents the risk to habitat i caused by pressure j. The identification of ecological barriers was mainly done with the help of ArcGIS extraction analysis tools, and the areas with high habitat quality were independently extracted and defined as ecological sources. Afterward, the areas with median habitat risk and above in the study area were overlaid with the ecological service supply areas to obtain the specific extent of the ecological barriers.

ESP Construction
The study used the Minimum Cumulative Resistance (MCR) model to construct a resistance surface for terrestrial ecosystems in Bohai Bay. The model requires consideration of three factors: the source of ecological service supply, distance, and landscape basal characteristics [48]. The source is the starting point for the outward diffusion of ecological services in the MCR model, with internal homogeneity and the ability to expand or attract in all directions, which in this study refers to the ES supply areas. The basic equation of the MCR model is as follows: where f represents the positive correlation between the minimum resistance at any point in space and its distance to all sources and the basal characteristics of the landscape. D ij is the spatial distance from source j across a landscape i to a point; R i is the cumulative resistance that needs to be overcome for the landscape i to extend to the source. By assigning a value to the degree of vegetation cover for different land use types, the higher the value, the higher the cost of crossing the patch. Referring to previous studies and taking into account the actual situation in Bohai Bay [49][50][51], the resistance values were set as shown in Table 4. Ecological corridors ensure the exchange of information, energy, and organisms between ecological sources and are an effective way to maintain the continuity of ecological functions and guarantee regional ESPs [52]. The two ends of the ecological corridor are connected to different ecological sources. To make the ecological corridor extraction more accurate, we used the ArcGIS Find Geometric Center tool to convert the faceted ecological sources into point elements (geometric centers of ecological sources). Based on the construction of the minimum cumulative resistance surface, the ecological corridor was extracted by using the ArcGIS minimum cost path analysis tool. This tool identified a potential path of ecosystem service delivery between ecological sources by calculating the cumulative minimum value of the cumulative cost generated by the different resistance of ecological sources through the landscape. The important ecological corridor transmits the ecological services of more than three ecological sources, and the common ecological corridor transmits the ecological services between two ecological sources. Based on the ecological corridor network, the intersections of the ecological corridors were defined as ecological strategic nodes using the element editing function of ArcGIS.

Habitat Quality in the Terrestrial Ecosystem of Bohai Bay
The habitat quality of terrestrial ecosystems in Bohai Bay was assessed based on the habitat quality module of the InVEST model, and the spatial and temporal evolution patterns of habitat quality and habitat degradation in the study area were investigated, as shown in Figure 3a,b. In general, from 1990 to 2018, the mean value of habitat quality in the study area was 0.2901, and the overall habitat quality was at a low level [53][54][55]. In the three time periods of 1990-2000, 2000-2010, and 2010-2018, habitat quality showed a small increase in the 1990-2000 period, a small increase in the 2000-2010, and 2010-2018 showed a large decrease and a large increase, with the opposite changes in the degree of habitat degradation. The change between habitat quality classes was particularly pronounced from 2000-2010, with large areas of the class IV habitat transforming into class V habitat during this decade. Over the past three decades, the degree of habitat degradation of the terrestrial ecosystem in Bohai Bay has increased by 14.59%, and the increase rate from 1990 to 2010 is low. The high-value area of habitat degradation of the terrestrial ecosystem in Bohai Bay is scattered, mostly concentrated in the south of the study area. The high degradation area is mainly concentrated in the middle of the study area, distributed around urban construction land, and the medium degradation area is distributed in a belt around Bohai Bay, mainly between cultivated land and offshore coast. The reason for such distribution is related to the distribution of land use types. Urban construction land and other construction land are mainly distributed in the middle of the study area; rural residential areas are scattered in the whole study area, and cultivated land is mainly distributed around rural residential areas. These areas have frequent human activities and strong environmental disturbance. Therefore, the degree of habitat degradation around these areas is relatively high, and the habitat quality is relatively low. In terms of time scale, the mean values of habitat quality in the Bohai Bay terrestrial ecosystem in 1990, 2000, 2010, and 2018 were 0.2931, 0.2943, 0.2685, and 0.3045, respectively, indicating that the habitat quality of the terrestrial ecosystem in Bohai Bay is improving. On a spatial scale, the overall habitat quality of the terrestrial ecosystem in Bohai Bay is roughly distributed as high in the east and low in the west. Habitat quality in Bohai Bay was classified from high to low in categories I-V, with category I habitat quality being the highest and category V habitat quality the lowest, as shown in Figure 3a and Table 5. indicating that the habitat quality of the terrestrial ecosystem in Bohai Bay is improving. On a spatial scale, the overall habitat quality of the terrestrial ecosystem in Bohai Bay is roughly distributed as high in the east and low in the west. Habitat quality in Bohai Bay was classified from high to low in categories I-V, with category I habitat quality being the highest and category V habitat quality the lowest, as shown in Figure 3a and Table 5. Specifically, the area of habitat quality category I only accounts for 0.4% of the total area of the region and is scattered in forested areas, such as forest nature reserves. Habitat quality category II areas are mainly located in the eastern part of the study area, along the coast of Bohai Bay, where the land use types are mainly near-coastal wetlands and lakes and reservoirs, including many wetland nature reserves, such as the Yellow River Delta, which are very rich in species diversity. Habitat quality category III areas are scattered throughout the study area, with concentrations in the northeast and central parts, where the land use types are mainly lakes, rivers, and grasslands, with high species diversity. Areas in habitat quality category IV, the largest of all categories, are scattered throughout the study area, with land use types dominated by arable land and rural residential land and heavily influenced by human disturbance, with serious damage to ecological modification and poor species diversity. Habitat quality category V is mainly concentrated in the central and northeastern parts of the study area and is also scattered throughout the study area. The land use types in this area are mainly urban construction land and other construction land, with frequent human activities, strong disturbance, low species diversity, and poor habitat quality. During the period 1990-2018, there were certain trends in the changes of habitat quality classes in the terrestrial ecosystem of Bohai Bay, as shown in Table 5. In general, the area of categories V and category II habitats in the study area increased by 333.28% and 103.88%, respectively, while the area of category I, category III, and category IV habitats decreased by 56.48%, 33.90%, and 83.09%, respectively. In terms of trends, the area of category I, category III, and category IV habitats showed an overall decreasing trend, while the area of category V and category II habitats showed an overall increasing trend, and the area of category V habitat increased significantly. This indicates that the quality of habitats in the study area is decreasing and polarizing. Figure 4 expresses the transfer into and out of the habitat area of each category, and the horizontal part of the figure represents the transfer out of each category of habitat, and the vertical part represents the transfer into each category of the habitat area. The red and green colors indicate the size of the transfer out and transfer into, respectively, and the larger the amount, the darker the color. The area transferred out of each category of habitat quality in Bohai Bay was ranked from largest to smallest: V > IV > II > III > I. The transferred area of each category of habitat quality in Bohai Bay was ranked from largest to smallest: III > V > IV > II > I. Figure 5 shows the spatial distribution of the transformation of habitat quality category. In the north of Bohai Bay and most coastal areas, the transformation of habitat quality category is relatively frequent. The habitat quality in the northern part of Bohai Bay has deteriorated significantly, while the habitat quality in the southern coastal area has improved.
The area transferred out of each category of habitat quality in Bohai Bay was ranked from largest to smallest: V > IV > II > III > I. The transferred area of each category of habitat quality in Bohai Bay was ranked from largest to smallest: III > V > IV > II > I. Figure 5 shows the spatial distribution of the transformation of habitat quality category. In the north of Bohai Bay and most coastal areas, the transformation of habitat quality category is relatively frequent. The habitat quality in the northern part of Bohai Bay has deteriorated significantly, while the habitat quality in the southern coastal area has improved.  In this study, Global Moran's I cluster analysis and hotspot analysis were used to further investigate the spatial distribution characteristics and patterns of habitat quality in the terrestrial ecosystems of Bohai Bay. The Global Moran's I was used to describe the average degree of association of all spatial units with the surrounding area over the entire region. Moran's I > 0 indicates a spatial positive correlation. The larger its value, the more obvious the spatial correlation. Moran's I < 0 indicates spatial negative correlation. The smaller its value, the greater the spatial difference. Otherwise, Moran's I = 0, and the space is random. This paper used hotspot analysis to identify the spatial location and area changes of high and low habitat quality in Bohai Bay. In Figure 3c, +3 and −3 represent high-value aggregation and low-value aggregation with 99% confidence in this region, respectively, and +2 and −2 represent high-value aggregation and low-value aggregation In this study, Global Moran's I cluster analysis and hotspot analysis were used to further investigate the spatial distribution characteristics and patterns of habitat quality in the terrestrial ecosystems of Bohai Bay. The Global Moran's I was used to describe the average degree of association of all spatial units with the surrounding area over the entire region. Moran's I > 0 indicates a spatial positive correlation. The larger its value, the more obvious the spatial correlation. Moran's I < 0 indicates spatial negative correlation. The smaller its value, the greater the spatial difference. Otherwise, Moran's I = 0, and the space is random. This paper used hotspot analysis to identify the spatial location and area changes of high and low habitat quality in Bohai Bay. In Figure 3c, +3 and −3 represent high-value aggregation and low-value aggregation with 99% confidence in this region, respectively, and +2 and −2 represent high-value aggregation and low-value aggregation with 95% confidence in this area, respectively, and +1 and −1 represent high-value aggregation and low-value aggregation with 90% confidence in this area, respectively; 0 represents the regional habitat quality, and the aggregation phenomenon of high value or low value is not obvious. Cluster analysis can be used to analyze the spatial aggregation characteristics of a set of elements. In this study, it was used to analyze the spatial aggregation of high-value and low-value areas of habitat quality in Bohai Bay. In Figure 3d, HH represents habitat quality high-value aggregation areas, and the habitat quality in this area is stable; LL represents habitat quality low-value aggregation areas, and the habitat quality in this area is stable; HL represents the area with high values of habitat quality mainly surrounded by low values, and habitat quality in this region varies significantly; LH represents the area where the low value of habitat quality is mainly surrounded by high values, and habitat quality in this region varies significantly.
According to the spatial autocorrelation analysis of habitat quality, there is some spatial correlation in the distribution of habitat quality in the study area. The Moran's I increased from 0.0816 in 1990 to 0.1663 in 2018, an increase of 103.80%, which indicates a certain spatial correlation and an increasing trend in the correlation of habitat quality in Bohai Bay, as shown in Table 6. Based on the cluster analysis, the spatial distribution of high habitat quality areas and low habitat quality areas in Bohai Bay can be understood. The proportion of the area of HL and LH areas of terrestrial ecosystem habitat quality in Bohai Bay to the total study area increased significantly in 1990, 2000, 2010, and 2018, with 18.98%, 19.78%, 45.28%, and 45.49%, respectively, which indicates that the mutual conversion between high and low habitat quality areas in Bohai Bay has increased in the last three decades. Based on the habitat quality hotspot analysis from 1990 to 2018, the proportion of cold spot area to the total area of the study area increased from 14.91% to 19.57%, and the proportion of hotspot area to the total area of the study area decreased from 14.48% to 13.31%, which indicates that the area of high habitat quality in the study area is shrinking, while the area of low habitat quality is continuously spreading showing aggregated distribution. With the extraction and analysis function of ArcGIS, areas with habitat quality ratings of high in Bohai Bay in 2018 were extracted, and these areas were defined as ecological sources, as shown in Figure 6. the land side; the medium-risk area is mainly located in the northeast of Wudi County in Binzhou City in Shandong Province, in the Hekou District in Dongying City, and in most of the ecological sources of Huanghua County in Hebei Province. The medium-risk areas are mainly located in the northeastern part of Wudi County, Binzhou City, Shandong Province, the northeastern part of Hekou District, Dongying City, and most of Huanghua County in Hebei Province, in the ecological sources.  The total area of the ecological sources in Bohai Bay is 2873.25 km 2 , accounting for 14.28% of the total area of the study area. Among them, the area at risk is 2821.25 km 2 , accounting for 98.19% of the total area of the ecological sources, which is mainly medium and low risk. Among the risky areas, the medium-risk area is 1191.5 km 2 , accounting for 41.47% of the ecological sources; the low-risk area is 1128 km 2 , accounting for 39.26% of the ecological sources; and the high-risk area is 501.75 km 2 , accounting for 17.46% of the ecological sources. The risk-free area is very small and is mainly located in the center of the Yellow River Delta Natural Reserve in Shandong Province; the low-risk area is mainly located in the ecological sources of Dongying City and Binzhou City in Shandong Province near the sea side and the ecological sources of Caofeidian City in Hebei Province near the land side; the medium-risk area is mainly located in the northeast of Wudi County in Binzhou City in Shandong Province, in the Hekou District in Dongying City, and in most of the ecological sources of Huanghua County in Hebei Province. The medium-risk areas are mainly located in the northeastern part of Wudi County, Binzhou City, Shandong Province, the northeastern part of Hekou District, Dongying City, and most of Huanghua County in Hebei Province, in the ecological sources. The risk situation of ecological sources in different provinces of Bohai Bay is as follows ( Figure 8): the mean values of risk to ecological sources from human activities in Tianjin, Hebei, and Shandong provinces are very similar, with risk values of 2.69, 2.73, and 2.63, respectively. In Tianjin and Shandong provinces, the three main sources of risk to ecological sources are industrial and construction sites, urban sites, and reclamation areas. In Hebei Province, the three largest sources of risk to ecological sources are industrial and construction sites, urban sites, and rural settlements. The difference in the risk of regions from the perspective of the contribution of different risk sources to the risk value is shown in Figure 9. The proportion of risk caused by paddy fields to the ecological sources of Hebei in Bohai Bay is the largest at 10.47% of the total risk value, which is much larger than the proportion of risk caused by paddy fields The risk situation of ecological sources in different provinces of Bohai Bay is as follows (Figure 8): the mean values of risk to ecological sources from human activities in Tianjin, Hebei, and Shandong provinces are very similar, with risk values of 2.69, 2.73, and 2.63, respectively. In Tianjin and Shandong provinces, the three main sources of risk to ecological sources are industrial and construction sites, urban sites, and reclamation areas. In Hebei Province, the three largest sources of risk to ecological sources are industrial and construction sites, urban sites, and rural settlements. The risk situation of ecological sources in different provinces of Bohai Bay is as follows ( Figure 8): the mean values of risk to ecological sources from human activities in Tianjin, Hebei, and Shandong provinces are very similar, with risk values of 2.69, 2.73, and 2.63, respectively. In Tianjin and Shandong provinces, the three main sources of risk to ecological sources are industrial and construction sites, urban sites, and reclamation areas. In Hebei Province, the three largest sources of risk to ecological sources are industrial and construction sites, urban sites, and rural settlements. The difference in the risk of regions from the perspective of the contribution of different risk sources to the risk value is shown in Figure 9. The proportion of risk caused by paddy fields to the ecological sources of Hebei in Bohai Bay is the largest at 10.47% of the total risk value, which is much larger than the proportion of risk caused by paddy fields  The difference in the risk of regions from the perspective of the contribution of different risk sources to the risk value is shown in Figure 9. The proportion of risk caused by paddy fields to the ecological sources of Hebei in Bohai Bay is the largest at 10.47% of the total risk value, which is much larger than the proportion of risk caused by paddy fields to Shandong and Tianjin at 4.79%. The proportion of risk posed by drylands to the ecological sources in Tianjin is the smallest at 11.47%, which is lower than its risk to Hebei at 15.83%. The proportion of risk posed by urban land to the ecological sources in Tianjin is the largest at 20.78%, while its risk to Hebei is the lowest at 18.06%. The proportion of risk posed by rural settlements to the ecological sources in Hebei is the largest at 17.16%, while its risk to Shandong is the lowest at 14.16%. The proportion of risk posed by industrial and construction land to the ecological sources in Tianjin is very large at 23.43% of the total risk value, while its risk to Hebei is the lowest at 18.11%. The proportion of risk posed by unused land to the ecological sources in Shandong in Bohai Bay is the largest at 10.37% of the total risk value, much larger than its risk to Hebei at 5.00%. The proportion of risk to the ecological sources in Tianjin, Shandong, and Hebei is relatively similar, at 15.87%, 16.83%, and 15.36%, respectively. to Shandong and Tianjin at 4.79%. The proportion of risk posed by drylands to the ecological sources in Tianjin is the smallest at 11.47%, which is lower than its risk to Hebei at 15.83%. The proportion of risk posed by urban land to the ecological sources in Tianjin is the largest at 20.78%, while its risk to Hebei is the lowest at 18.06%. The proportion of risk posed by rural settlements to the ecological sources in Hebei is the largest at 17.16%, while its risk to Shandong is the lowest at 14.16%. The proportion of risk posed by industrial and construction land to the ecological sources in Tianjin is very large at 23.43% of the total risk value, while its risk to Hebei is the lowest at 18.11%. The proportion of risk posed by unused land to the ecological sources in Shandong in Bohai Bay is the largest at 10.37% of the total risk value, much larger than its risk to Hebei at 5.00%. The proportion of risk to the ecological sources in Tianjin, Shandong, and Hebei is relatively similar, at 15.87%, 16.83%, and 15.36%, respectively.

Ecological Barriers Identification
Ecological barriers should be constructed in ecological sources that are at high risk of being affected by human activities. The barrier areas cannot only reduce the interference of human activities with ecological sources, but also help to amplify or transmit ecological services. This study identifies areas suitable for ecological barrier construction with the help of habitat quality and habitat risk models and uses the extraction tool of ArcGIS to eliminate discrete areas to obtain nine ecological barrier construction areas, which are divided into three major zones according to their provinces. Hebei is the first ecological barrier zone with three ecological barriers; Tianjin is the second ecological barrier zone with one ecological barrier, and Shandong is the third ecological barrier zone with five ecological barriers, as shown in Figure 10.

Ecological Barriers Identification
Ecological barriers should be constructed in ecological sources that are at high risk of being affected by human activities. The barrier areas cannot only reduce the interference of human activities with ecological sources, but also help to amplify or transmit ecological services. This study identifies areas suitable for ecological barrier construction with the help of habitat quality and habitat risk models and uses the extraction tool of ArcGIS to eliminate discrete areas to obtain nine ecological barrier construction areas, which are divided into three major zones according to their provinces. Hebei is the first ecological barrier zone with three ecological barriers; Tianjin is the second ecological barrier zone with one ecological barrier, and Shandong is the third ecological barrier zone with five ecological barriers, as shown in Figure 10.

Ecological Security Patterns
The core elements of the ESPs of the Bohai Bay terrestrial ecosystem include the ecological sources, ecological barriers, ecological corridors, and ecological strategic nodes. All the least costly paths were combined and filtered to obtain a total of 35 ecological corridors, which together form the Bohai Bay Ecological Corridor Network, including 16 important ecological corridors and 19 common ecological corridors. Most of the 19 common corridors are located in Binzhou and Dongying in Shandong Province, where the ecological barriers are relatively dense and the corridors are interspersed. Ecological strategic nodes are located mainly at the intersections and branches of important ecological corridors. By combining the ecological sources, ecological barrier construction areas, minimum cumulative resistance surface, ecological corridors, and ecological strategic nodes, the ESPs of Bohai Bay was obtained.
The ecological strategic node is an area that plays a decisive role in regional landscape connectivity. It includes not only the areas that play a key role in connecting ESs and functions, but also the areas with weak functions on the ecological corridor and the areas most vulnerable to interference on the ecological corridor. Based on the ecological corridor network, the intersection of the ecological corridor is defined as the ecological strategy node by using the element editing function of ArcGIS, and 32 nodes were finally determined, as shown in Figure 11.

Ecological Security Patterns
The core elements of the ESPs of the Bohai Bay terrestrial ecosystem include the ecological sources, ecological barriers, ecological corridors, and ecological strategic nodes. All the least costly paths were combined and filtered to obtain a total of 35 ecological corridors, which together form the Bohai Bay Ecological Corridor Network, including 16 important ecological corridors and 19 common ecological corridors. Most of the 19 common corridors are located in Binzhou and Dongying in Shandong Province, where the ecological barriers are relatively dense and the corridors are interspersed. Ecological strategic nodes are located mainly at the intersections and branches of important ecological corridors. By combining the ecological sources, ecological barrier construction areas, minimum cumulative resistance surface, ecological corridors, and ecological strategic nodes, the ESPs of Bohai Bay was obtained.
The ecological strategic node is an area that plays a decisive role in regional landscape connectivity. It includes not only the areas that play a key role in connecting ESs and functions, but also the areas with weak functions on the ecological corridor and the areas most vulnerable to interference on the ecological corridor. Based on the ecological corridor network, the intersection of the ecological corridor is defined as the ecological strategy node by using the element editing function of ArcGIS, and 32 nodes were finally determined, as shown in Figure 11.

Using Ecological Risk Assessment for Ecological Source Zoning
At present, in the research on the identification of ecological sources, one takes the ecosystem importance assessment as the only basis for the screening of ecological sources, and the other combines the ecosystem importance assessment with the risk assessment of ecosystem degradation [13,56]. It only assesses the importance of regional ecosystems and evaluates the current structure and function of ecosystems, ignoring the interaction of human activities with ecosystems and the long-term potential of the ecological source to provide ESs. Taking the risk of ecosystem degradation as another indicator for the screening of ecological sources, although the screened ecological sources have long-term potential to provide ESs, this source identification method directly ignores the ecological sources that can provide ESs but bears high risks. Taking the importance of ecosystem as the basis for the identification of ecological sources and conducting risk analysis on the identified ecological sources can comprehensively show the current state of ecological sources and the long-term potential of providing ESs. The risk assessment identifies internal differences in ecological sources and reflects them in spatial distribution, allowing for finegrained and adaptive ecological management. In this study, risk assessment was used to grade the risk of human activities in the terrestrial ecosystems of Bohai Bay. Most of the ecological sources in Bohai Bay are at varying degrees of risk, mainly medium and low risk, while about 20% of the ecological source areas are at high risk. Most of the sporadically distributed ecological sources are at high risk due to their high degree of habitat fragmentation and complex sources of risk, while most of the ecological sources near industrial and construction sites and towns are also at high risk due to the intensity of disturbance and frequent disturbance. Medium-risk areas are mainly located on the marine side of the ecological sources in Caofeidian District, Tangshan City, and in most of the ecological sources in Huanghua County in Hebei Province. Low-risk areas are mainly located in Dongying City and Binzhou City in Shandong Province, near the marine side of

Using Ecological Risk Assessment for Ecological Source Zoning
At present, in the research on the identification of ecological sources, one takes the ecosystem importance assessment as the only basis for the screening of ecological sources, and the other combines the ecosystem importance assessment with the risk assessment of ecosystem degradation [13,56]. It only assesses the importance of regional ecosystems and evaluates the current structure and function of ecosystems, ignoring the interaction of human activities with ecosystems and the long-term potential of the ecological source to provide ESs. Taking the risk of ecosystem degradation as another indicator for the screening of ecological sources, although the screened ecological sources have long-term potential to provide ESs, this source identification method directly ignores the ecological sources that can provide ESs but bears high risks. Taking the importance of ecosystem as the basis for the identification of ecological sources and conducting risk analysis on the identified ecological sources can comprehensively show the current state of ecological sources and the long-term potential of providing ESs. The risk assessment identifies internal differences in ecological sources and reflects them in spatial distribution, allowing for fine-grained and adaptive ecological management. In this study, risk assessment was used to grade the risk of human activities in the terrestrial ecosystems of Bohai Bay. Most of the ecological sources in Bohai Bay are at varying degrees of risk, mainly medium and low risk, while about 20% of the ecological source areas are at high risk. Most of the sporadically distributed ecological sources are at high risk due to their high degree of habitat fragmentation and complex sources of risk, while most of the ecological sources near industrial and construction sites and towns are also at high risk due to the intensity of disturbance and frequent disturbance. Medium-risk areas are mainly located on the marine side of the ecological sources in Caofeidian District, Tangshan City, and in most of the ecological sources in Huanghua County in Hebei Province. Low-risk areas are mainly located in Dongying City and Binzhou City in Shandong Province, near the marine side of the ecological sources. The area of no-risk areas is very small, with only a small area within the Yellow River Delta Nature Reserve. Using the HRA model for risk assessment can evaluate the risk level of ecological sources, identify the risk causes of ecological sources from the perspective of exposure and sensitivity (Figure 12), and provide a basis for formulating corresponding risk mitigation measures. The sources of risk in the ecological sources of Bohai Bay are mainly industrial and construction land, rural settlements, and urban land, which are closely related to the type of land use after the reclamation. Generally speaking, urban construction land poses a greater risk to ecological sources than rural settlements, but in Bohai Bay, as part of the ecological sources are scattered in various areas, they show a staggered distribution with rural settlements, resulting in a wider range of disturbance from rural settlements, while urban construction land is mainly concentrated. Although the extreme value of the risk caused by rural settlements to the ecological sources is not as large as that caused by urban construction land, the average value of the risk caused by rural settlements is larger than urban construction land. The risk to the ecological sources in Bohai Bay is more due to the vulnerability of the ecosystem itself, and the ecological restoration work in Bohai Bay can focus on improving the ability of the area to resist the risk sources. This paper introduces ecological risk assessment to demonstrate the structural and functional differentiation within ecological sources. The ecological sources are divided into high, medium, and low zones according to risk levels, and the exposure and sensitivity analysis of the HRA model provides the basis for adaptive ecological management. the ecological sources. The area of no-risk areas is very small, with only a small area within the Yellow River Delta Nature Reserve. Using the HRA model for risk assessment can evaluate the risk level of ecological sources, identify the risk causes of ecological sources from the perspective of exposure and sensitivity (Figure 12), and provide a basis for formulating corresponding risk mitigation measures. The sources of risk in the ecological sources of Bohai Bay are mainly industrial and construction land, rural settlements, and urban land, which are closely related to the type of land use after the reclamation. Generally speaking, urban construction land poses a greater risk to ecological sources than rural settlements, but in Bohai Bay, as part of the ecological sources are scattered in various areas, they show a staggered distribution with rural settlements, resulting in a wider range of disturbance from rural settlements, while urban construction land is mainly concentrated. Although the extreme value of the risk caused by rural settlements to the ecological sources is not as large as that caused by urban construction land, the average value of the risk caused by rural settlements is larger than urban construction land. The risk to the ecological sources in Bohai Bay is more due to the vulnerability of the ecosystem itself, and the ecological restoration work in Bohai Bay can focus on improving the ability of the area to resist the risk sources. This paper introduces ecological risk assessment to demonstrate the structural and functional differentiation within ecological sources. The ecological sources are divided into high, medium, and low zones according to risk levels, and the exposure and sensitivity analysis of the HRA model provides the basis for adaptive ecological management.

Identification of Ecological Barrier Construction Areas
Based on the risk zoning of ecological sources and the analysis of the concept and connotation of ecological barriers, this paper presents a differentiated ecological management of ecological sources by identifying ecological barrier zones. The main function of ecological barriers is to resist the adverse impact of risk sources on the ecological service supply area and transfer ESs to the demand area. In Figure 13a, the embedded ecological barrier shows that the beneficiary and supply areas of ESs are in the same location, and the ecological barrier is at the periphery of the beneficiary or supply area, acting as a barrier to undesirable factors. In Figure 13b, a guiding ecological barrier shows that there is a directionality of ESs from the supply area to the beneficiary area, and a connection area is needed so that ESs reach the beneficiary area along the connection area. When there is a spatial mismatch between the ecosystem service supply area and the beneficiary area,

Identification of Ecological Barrier Construction Areas
Based on the risk zoning of ecological sources and the analysis of the concept and connotation of ecological barriers, this paper presents a differentiated ecological management of ecological sources by identifying ecological barrier zones. The main function of ecological barriers is to resist the adverse impact of risk sources on the ecological service supply area and transfer ESs to the demand area. In Figure 13a, the embedded ecological barrier shows that the beneficiary and supply areas of ESs are in the same location, and the ecological barrier is at the periphery of the beneficiary or supply area, acting as a barrier to undesirable factors. In Figure 13b, a guiding ecological barrier shows that there is a directionality of ESs from the supply area to the beneficiary area, and a connection area is needed so that ESs reach the beneficiary area along the connection area. When there is a spatial mismatch between the ecosystem service supply area and the beneficiary area, there is a connecting area between them, which generates ecosystem service flows, and the connecting area influences the ES flow process. In Figure 13c, a strong transmission ecological barrier indicates that there is an exchange of ESs within the ecosystem service supply area and the ecological barriers, and the ecological barriers have a diffusion effect on the transfer of ESs. In Figure 13d, an external diffusion-type ecological barrier indicates that ESs spread in all directions without direction; the supply and benefit areas are adjacent to each other, and the ecological barriers are around the supply area, acting as a barrier to undesirable factors and facilitating the flow of ESs. there is a connecting area between them, which generates ecosystem service flows, and the connecting area influences the ES flow process. In Figure 13c, a strong transmission ecological barrier indicates that there is an exchange of ESs within the ecosystem service supply area and the ecological barriers, and the ecological barriers have a diffusion effect on the transfer of ESs. In Figure 13d, an external diffusion-type ecological barrier indicates that ESs spread in all directions without direction; the supply and benefit areas are adjacent to each other, and the ecological barriers are around the supply area, acting as a barrier to undesirable factors and facilitating the flow of ESs. However, existing studies treat ecological barriers more as macroscopic concepts with a huge research scale, making it difficult to analyze their structure and function after identifying the extent of the ecological barriers. In regional or city-scale studies, this way of constructing ecological barriers is similar, that is, based on methods, such as environmental carrying capacity assessment and analysis of ecosystem service functions, and areas with the excellent ecological environment in the region are identified as ecological barrier areas. This research approach neither fulfils the function of ecological barriers in maintaining ecological service supply areas, nor does it reflect the role of ecological barrier areas in facilitating the delivery of ecosystem services. The areas identified in this study for constructing ecological barriers through ecosystem service provisioning area identification and risk assessment can maintain the stability of the structure and function of the ecological service provisioning areas, while facilitating the delivery of ESs.

The Relationship between Ecological Barriers and Ecological Security Patterns
Both ecological barriers and ESPs are landscape ecological plans to maintain regional ecological security [5,31]. The construction of ecological barriers is to protect the structure and function of ecological sources, and the construction of ESPs is to optimize land spatial However, existing studies treat ecological barriers more as macroscopic concepts with a huge research scale, making it difficult to analyze their structure and function after identifying the extent of the ecological barriers. In regional or city-scale studies, this way of constructing ecological barriers is similar, that is, based on methods, such as environmental carrying capacity assessment and analysis of ecosystem service functions, and areas with the excellent ecological environment in the region are identified as ecological barrier areas. This research approach neither fulfils the function of ecological barriers in maintaining ecological service supply areas, nor does it reflect the role of ecological barrier areas in facilitating the delivery of ecosystem services. The areas identified in this study for constructing ecological barriers through ecosystem service provisioning area identification and risk assessment can maintain the stability of the structure and function of the ecological service provisioning areas, while facilitating the delivery of ESs.

The Relationship between Ecological Barriers and Ecological Security Patterns
Both ecological barriers and ESPs are landscape ecological plans to maintain regional ecological security [5,31]. The construction of ecological barriers is to protect the structure and function of ecological sources, and the construction of ESPs is to optimize land spatial planning by using the ESs provided by ecological sources. From the perspective of ecological restoration, ecological source is not only the protection object of ecological barrier, but also the basis of ESPs [6,7]. Ecological barriers and ESPs are complementary in both structure and function. The protection of ecological sources by ecological barriers safeguards the long-term potential of the ESPs to provide ESs, while ecological corridors linking multiple ecological sources in the ESPs complete the ecological barriers to facilitate the exchange and transfer of ESs to ecological service beneficiary areas and other ecological service supply areas. In Figure 14, (c) is the basic ESP construction paradigm, focusing on the identification of ecological elements and the correlation among them. (d) is the optimized ESP construction methods, where the interrelationship between human activities and ecological elements is taken into account through the construction of ecological barriers. In Figure 14, (c) ecological sources not only transmit ESs to ecosystem service demand areas, but also exchange ESs among multiple ecological sources. Ecological sources are not only exporters of ESs, but also receivers of adverse impacts from risk sources, and there is instability in the structure and function of ecological sources when disturbed by risk sources. In Figure 14, (d) ecological barriers are constructed in areas of ecological sources disturbed by risk sources. The use of ecological barriers improves the stability of the structure and function of ecological sources, thus safeguarding the long-term potential of the ESPs to maintain regional ecological security. planning by using the ESs provided by ecological sources. From the perspective of ecological restoration, ecological source is not only the protection object of ecological barrier, but also the basis of ESPs [6,7]. Ecological barriers and ESPs are complementary in both structure and function. The protection of ecological sources by ecological barriers safeguards the long-term potential of the ESPs to provide ESs, while ecological corridors linking multiple ecological sources in the ESPs complete the ecological barriers to facilitate the exchange and transfer of ESs to ecological service beneficiary areas and other ecological service supply areas. In Figure 14, (c) is the basic ESP construction paradigm, focusing on the identification of ecological elements and the correlation among them. (d) is the optimized ESP construction methods, where the interrelationship between human activities and ecological elements is taken into account through the construction of ecological barriers. In Figure 14, (c) ecological sources not only transmit ESs to ecosystem service demand areas, but also exchange ESs among multiple ecological sources. Ecological sources are not only exporters of ESs, but also receivers of adverse impacts from risk sources, and there is instability in the structure and function of ecological sources when disturbed by risk sources. In Figure 14, (d) ecological barriers are constructed in areas of ecological sources disturbed by risk sources. The use of ecological barriers improves the stability of the structure and function of ecological sources, thus safeguarding the longterm potential of the ESPs to maintain regional ecological security.

Conclusions
In rapidly urbanizing developing countries, where economic development and resource conservation are in conflict, ecological barriers and ESPs in regional landscape

Conclusions
In rapidly urbanizing developing countries, where economic development and resource conservation are in conflict, ecological barriers and ESPs in regional landscape planning are both practical ways of guarding regional ecological resources and maintaining regional ecological security using ecological bottom-line thinking. In China, the construction of ecological barriers as a vague concept has played a remarkable role in the protection and restoration of ecosystem structure and function. However, most of the existing ecological barrier construction takes megaregions as the research scale, equating areas with the excellent ecological environments as ecological barriers and focusing on the maintenance of ecological security of the whole region. The overly large research object, combined with the complexity of the ecosystem itself, makes the analysis of ecological barrier structures and functions complicated. This study identifies ecological barriers at the scale of ecological patches using habitat quality and habitat risk assessment, focusing on the maintenance of ecological patches and the transmission of ESs by ecological barriers. At the same time, ecological barriers are integrated into the ESPs, and the framework of "ecological sources-ecological barriers-ecological resistance surface-ecological corridors" was used to construct the ESPs.
The results of the habitat quality analysis show that the overall habitat quality of the Bohai Bay land area is at a low level, and the spatial distribution is polarized. The area of ecological sources is decreasing in the Bohai Bay land area. The results of the habitat risk analysis show that the ecological sources in the Bohai Bay land area are mainly of medium-risk and low-risk, and the risk mainly comes from the vulnerability of the ecosystem itself. The ecological sources in the Bohai Bay land area are mainly located in the land-sea interface, and the land use types are mainly woodlands, rivers, and grasslands, accounting for 14.3% of the entire area. A total of nine ecological barrier construction cells were identified, accounting for 40.8% of the ecological sources and were divided into three major zones according to the provinces they belonged to. A total of 35 ecological corridors and 32 ecological strategic nodes with rivers and woodlands as the main land types were obtained, including 16 important ecological corridors and 19 general ecological corridors.
However, there are still many problems in the study of ESPs. Although land use is a comprehensive indicator of the structure and function of regional ecosystems and has been used as the basis for most studies on ESs, the link between land cover and ESs behind land use has not been fully resolved. This study focuses on the specific scope and relative capacity of the study area to provide ESs, but does not consider the types of ESs that can be provided in the ecosystem service supply area. Moreover, there is a lack of evaluation on the effectiveness of ecological barriers and ESPs. Therefore, more research is needed to deal with these problems.