Link Ecological and Social Composite Systems to Construct Sustainable Landscape Patterns: A New Framework Based on Ecosystem Service Flows

: Integrating the ﬂow of supply and demand of ecosystem services (ESs) into the ecological security pattern (ESP) of coastal ecosystems with extremely fragile ecological backgrounds and contradictory human–land relationships is beneﬁcial to the coordinated development of human–land systems. However, existing studies ignore the issue of scales of supply–demand linkages, making the ESP not properly guide sustainable development. Based on ESs delivery chain theory and landscape ecology approaches, we developed a sustainable development framework consisting of coupled microscopic natural–social systems. The method was tested using data from the Liao River Delta. In this study area, the natural supply potential and demand mapping distribution of key ESs were assessed to identify ecological sources in the Liao River Delta, a typical coastal zone in northern China. The resistance surface based on land use type assignment was modiﬁed using hydrological connectivity frequency and nighttime light intensity. Ecological corridors were extracted and optimized using a minimum cumulative resistance model and connectivity evaluation. The study found that the high supply area and the high demand reﬂection area are not consistent in location and supply level. Ecological source areas are evenly distributed, accounting for 12% of the total area. The ecological corridors are mainly concentrated in the west and southeast and do not cross the built-up areas in the east. This ESP framework safeguards the local demand for natural products and the natural potential to maintain services over the longer term and to a larger scale while informing the development of environmental management measures.


Introduction
Intensive human activities have sharply compressed urban ecological space, threatening the sustainable development of human society. For coastal cities, the intensive use of coastal wetlands and other ecological resources by reclamation activities is typical [1]. At the same time, human society has to rely on natural ecosystems for survival and development. They provide essential services to people, including a fresh water supply, carbon sequestration, and the provision of fish and wildlife habitat [2,3]. Under the combined pressures of urban sprawl and declining natural services, intervention strategies that prioritize natural spaces, considering the potential for natural services and human needs, are vital [4,5]. There is a critical need for methods that reflect supply and demand dynamics in oil field in the country. Under the background of severe ecological and environmental problems in the study area, existing small-scale restoration activities ignore ecosystem supply and demand and large-scale spatial integration. In contrast, ecosystem security pattern construction can compensate for the lack of research objects and study scales [26]. Therefore, the objectives of this study are to (1) deconstruct and quantify the process of spatial flow transport and delivery of ES so that ES demand is reflected in the natural spaces that provide local services at a clear spatial extent and size; (2) propose a new method for identifying ecological source sites by combining ES demand allocation mapping and natural supply potential; (3) use a minimum cumulative resistance (MCR) model, combined with hydrological connectivity assessment and nighttime light intensity (NLI) to extract ecological corridors and further optimize them; and (4) propose relevant recommendations for regional ecological management.

Study Area and Data Sources
Located in Northeast China, the Liao River Delta has one of the best-preserved wetlands globally ( Figure 1). We have delineated the study boundary based on the principle of wetland and administrative integrity. The total area is 1388.28 km 2 , of which wetlands account for 51.44%. The total population is about 480,000, mainly in the east and northwest. Several rivers, including the Liao River, flow into the sea here. The whole region has a temperate semi-humid monsoon climate with an average annual rainfall of 650 mm. Considering the study area's characteristics and the core ecological risks, we selected the following four key ESs: (1) habitat quality, where 446 species of various wildlife, mainly birds, are distributed; (2) carbon sequestration, where coastal wetlands are an important carbon sink in the Earth's surface system [27]; (3) water purification, related to surface water quality management at the national level; and (4) landscape aesthetic services, in the knowledge economy era, one of the urban dwellers' urgent needs and development priorities [21]. March 2022)). (8) The inundation probability was derived from a surface water remotesensing dataset (http://global-surface-water.appspot.com/ (accessed on 5 March 2022)). (9) Traffic points of interest came from map data service providers (http://www.rivermap.cn/ (accessed on 5 March 2022)). (10) Socioeconomic statistics: average grain production, average market price of grain, national per capita tourism income; per capita tourism income of the study area, etc., were from the Panjin City Statistical Yearbook (http://www.panjin.gov.cn/ (accessed on 5 March 2022)) and the China Statistical Yearbook (http://www.stats.gov.cn/ (accessed on 5 March 2022)). All spatial data were interpolated and refined to a 30 m resolution raster using ArcGIS.

ESP Identification Framework
In this paper, a sustainable landscape pattern of coastal wetland ecosystem is proposed. This ESP is an ecological network based on the theory of landscape pattern optimization and ES delivery chain consisting of core patches that provide actual ESs, highvalue ecological functional areas, and ecological corridors. Through the protection and management of important ecological elements, it has the following benefits. Firstly, it safeguards the demand of the socioeconomic system for natural patches to provide indirect Ten data types, including the following for 2020, were collected for this study. (1) Land use data and Normalized difference vegetation index (NDVI) data were interpreted and . Land use types are divided into 14 categories: tidal flat, tidal creek, intertidal salt marsh, reed, inland salt plant, river, inland water, paddy field, dryland, residential land, industrial land, estuarine water and sandbar, uncultivated land and aquaculture. Based on a comparison with data from high-resolution remote-sensing imagery and field validation sites, the Kappa index was 0.91, and the classification results were reliable. The accuracy of each category and the sample size are detailed in Supplementary Figure S1. (2) Soil data, including aboveground biomass and soil organic carbon content, were obtained from data from 121 sampling sites measured by the team. (3) Elevation data were from the National Aeronautics and Space Administration (NASA; http://reverb.echo.nasa.gov/ (accessed on 2 March 2022)) and had a spatial resolution of 30 m; the data were processed by cropping and filling depressions. (4) Population data with a spatial resolution of 100 m were from WorldPop (https://www.worldpop.org/ (accessed on 5 March 2022)). (5) Precipitation data were monthly surface climate data from the China Meteorological Data Service Center (http://data.cma.cn/ (accessed on 5 March 2022)). The data of seven meteorological stations in the study area and adjacent areas were selected and spatially interpolated by an ArcGIS inverse distance weighting method to obtain a meteorological raster. (6) Suomi NPP VIIRS nighttime light data were from NOAA (https://ngdc.noaa.gov/ (accessed on 5 March 2022)). (7) A total of 6467 migratory bird observation records from 2013 to the present were obtained from the open-source eBird website (www.ebird.org (accessed on 5 March 2022)). (8) The inundation probability was derived from a surface water remote-sensing dataset (http://global-surface-water.appspot.com/ (accessed on 5 March 2022)). (9) Traffic points of interest came from map data service providers (http://www.rivermap.cn/ (accessed on 5 March 2022)). (10) Socioeconomic statistics: average grain production, average market price of grain, national per capita tourism income; per capita tourism income of the study area, etc., were from the Panjin City Statistical Yearbook (http://www.panjin.gov.cn/ (accessed on 5 March 2022)) and the China Statistical Yearbook (http://www.stats.gov.cn/ (accessed on 5 March 2022)). All spatial data were interpolated and refined to a 30 m resolution raster using ArcGIS.

ESP Identification Framework
In this paper, a sustainable landscape pattern of coastal wetland ecosystem is proposed. This ESP is an ecological network based on the theory of landscape pattern optimization and ES delivery chain consisting of core patches that provide actual ESs, high-value ecological functional areas, and ecological corridors. Through the protection and management of important ecological elements, it has the following benefits. Firstly, it safeguards the demand of the socioeconomic system for natural patches to provide indirect or direct ESs; secondly, it enhances network connectivity and reduces ecological problems of habitat fragmentation and patch isolation; thirdly, it regulates green infrastructure configuration according to the matching relationship between supply and demand of natural patch services. This ESP construction process mainly includes the following two parts.
Mapping the supply and demand of ESs: Human well-being and economic activities depend heavily on ecosystem functions and processes. We constructed a multi-source database based on data availability and regional characteristics and mapped the spatial distribution and magnitude of supply capacity for habitat quality, water purification, carbon sequestration, and landscape aesthetics services. Considering the spatial separation of supply and benefit areas, we hypothesized that ES demand could also be allocated and projected to the corresponding natural spaces based on the movement of spatial flows from the sources that produce the services to the sinks that generate the demand. In this study area, ES demand originates not only from built-up settlements but also from people's pressure for surface water management and the requirement to protect migratory birds [28]. After identifying the demand subjects, we used an abstract spatial flow quantification method based on traffic flows. ES demands are reflected in natural space according to the flow direction, flow rate, boundaries, and distribution patterns along with distances of ecological flows [14,29]. Figure 2 illustrates the conceptual framework of this study.
projected to the corresponding natural spaces based on the movement of spatial flows from the sources that produce the services to the sinks that generate the demand. In this study area, ES demand originates not only from built-up settlements but also from people's pressure for surface water management and the requirement to protect migratory birds [28]. After identifying the demand subjects, we used an abstract spatial flow quantification method based on traffic flows. ES demands are reflected in natural space according to the flow direction, flow rate, boundaries, and distribution patterns along with distances of ecological flows [14,29]. Figure 2 illustrates the conceptual framework of this study. Ecological element identification: The supply and demand mapping of ESs was used to quantify the ecosystem background and the ecological pressure given by humans. Ecological sources are further identified by combining the mapping results of both. The resistance surface represents how different landscape units or habitat patches impede species migration [30]. It is mainly determined based on land use types and the influence of human facilities such as transportation and oil wells. Due to the influence of the regional hydrological environment on biological migration, this study incorporates hydrological connectivity to correct the resistance surface. As a feature of human activity intensity, NLI was also used to modify patch resistance [31]. Based on the resistance surface and connectivity evaluation, ecological corridors were extracted and optimized, facilitating more efficient ecological management [32]. Ecological sources and corridors form a complete ecological security pattern. Ecological element identification: The supply and demand mapping of ESs was used to quantify the ecosystem background and the ecological pressure given by humans. Ecological sources are further identified by combining the mapping results of both. The resistance surface represents how different landscape units or habitat patches impede species migration [30]. It is mainly determined based on land use types and the influence of human facilities such as transportation and oil wells. Due to the influence of the regional hydrological environment on biological migration, this study incorporates hydrological connectivity to correct the resistance surface. As a feature of human activity intensity, NLI was also used to modify patch resistance [31]. Based on the resistance surface and connectivity evaluation, ecological corridors were extracted and optimized, facilitating more efficient ecological management [32]. Ecological sources and corridors form a complete ecological security pattern.

Assessment of Ecosystem Services
These four key ESs encompass the many types of ESs specified in the CICES classification v4.3 (https://cices.eu (accessed on 10 March 2022)), such as regulating and supporting and cultural services, comprehensively reflecting the regional ecological security. Some of the services assessed represent relative rather than absolute sizes within the region, and we have standardized their classification to guide the identification of ecological elements. The detailed calculation process is shown in Tables 1 and 2.

ES Indicator Calculation Method Unit Explanation Data Source and Parameter Reference
Habitat Quality is the habitat quality of grid x in land use and land cover (habitat type) j; D xj is the stress level of grid x in land use and land cover (habitat type) j; k is the half-saturation constant, usually D xj half of the maximum value; H j is the habitat suitability of land use and land cover j; and z is a normalized constant, usually 2.5. [33] Carbon sequestration Sc is the total regional carbon storage, C_above is the aboveground carbon storage, C_below is the underground carbon storage, C_soil is the soil carbon storage, and C_dead is the carbon storage stored in dead organic matter, all in units (t/hm 2 ).
1, 2 [34,35] Water purification The vegetation purification ability index is constructed by using the vegetation phosphorus accumulation P i ; the left side of the multiplication sign represents the relative size of the i type in the vegetation accumulation of all land uses; the vegetation phosphorus content of the i type land use is averaged, and the NDVI is used as the correction coefficient to obtain each grid; and the vegetation purification capacity index of grid j.

1, 2
Landscape Aesthetics S L is the service supply of landscape aesthetics (CNY·hm −2 ); 1/7 is the ratio of the economic value of the equivalent factor of grain production per unit to the market value of food production per unit of arable land; Q is the average grain yield per unit of the study area (Kg·hm −2 ); P is the average market price of grain in the study area (CNY·kg −1 ); A i is the landscape aesthetic service equivalent factor of the i-th land use type; S is the total area of the study area (hm 2 ); S i is the area of the i-th land use type (hm 2 ); t is the national per capita tourism income; and T is the per capita tourism income of the study area.

1, 10 [36]
Note: The data source number is the data number in Section 2.1.

Habitat Quality
Indicated species protection Natural patches in the study area Q' j is the probability of empirical bird abundance records where k(s,h) is the Gaussian kernel function; h is the bandwidth; s is the distance from point x to the target point; X tj is the number of migratory birds based on eBirds uploaded by bird watchers over the years; and x ij is the distance from pixel j to i Distance (km). 7

Carbon sequestration
Carbon sink protection Carbon sink S j = S j × HM j S' j is the probability of grid j transitioning from a carbon sink to a carbon source. S j is the total carbon storage of grid j; HM j is the disturbance degree of human activities to grid j.

1, 2, 6
Water purification Environmental Department's Requirements for Surface Water Quality The catchment area of the water system in the study area A' i is the water purification capacity of natural patches in catchment i; load surf,j is the surface pollutant load of grid j in catchment i (kg/(ha·year)); NDR surf , f is grid j, the proportion of pollutant loads reaching the river; eff' j is the maximum retention efficiency of the land between the grid and the river; IC 0 and k are calibration parameters; and IC j is the terrain index, representing the probability of the nutrients in grid j reaching the river.
1, 3, 5 [37,38] Landscape Aesthetics Daily travel guarantee for residents S' i is the degree to which settlement j allocates people's travel needs to natural space i. POP j is the population of settlement j; DIST j is the distance between settlement j and natural space i.

Habitat Quality
Natural supply potential: Habitat quality is the ability of an ecosystem to provide the conditions necessary for survival and reproduction [39]. Habitat quality depends on the relative impact of each threat, the relative sensitivity of each habitat type to each threat, the distance between the habitat and the source of the threat, and the degree to which the land is legally protected. A decline in habitat quality is seen due to the increased intensity of nearby land use. The habitat quality module of the InVEST model generates a habitat quality index map by combining information on land cover and biodiversity threat factors [38]. The relevant parameters were set using an expert scoring method and adjusted concerning literature from similar study areas, with the participation of 10 experts, including ecologists, geographers, and planners (please see Supplementary Tables S1 and S2 in the Supplementary Materials).
Demand reflection: The study area is an important migration route for migratory birds in eastern China, the southernmost part of the natural breeding area and the northernmost part of the wintering area of the Grus japonensis, and one of only a few critical breeding sites for Larus saundersi in the world. In addition, compared to other waterbirds, the Grus japonensis and the Larus saundersi are more sensitive to changes in the external environment; hence, the Larus saundersi and the Grus japonensis were used as sensitive species and con- servation targets in this study [40,41]. From the perspective of human needs, conservation actions maintain or increase the abundance and diversity of other species used or enjoyed by people.
The main ranges of indicator species are reed and saltmarsh patches, and they may move arbitrarily through these patches. Kernel density estimation, a non-parametric method for estimating an unknown probability density function, was used to determine indicator species' spatial distribution [42]. The intersection of the estimated spatial distribution with the main habitat patches reflects conservation needs in natural space.

Carbon Sequestration
Natural supply potential: The study area has large areas of mature natural wetlands and artificial coastal wetlands, which play an essential role in absorbing greenhouse gases from the atmosphere and mitigating global warming. Normally, organic matter is released back into the atmosphere in the form of carbon dioxide when plants wither. Due to the particular environment of wetland soils, the decomposition process will be slowed down or terminated by the lack of oxygen and low pH so that carbon is sequestered in the peat soil. In the long term, wetlands are global carbon sinks [43]. This study modeled carbon stocks based on aboveground biomass, belowground biomass, soil carbon pools, and dead organic matter [35]. The carbon pool table and carbon sink accounting instructions for the region is detailed in Supplementary Table S3.
Demand reflection: Carbon flow has the characteristics of in situ, directional, and omnidirectional flow, and its range is global [44]. In large-scale studies, the proportion of natural cover to the area of the administrative district is often used to represent carbon sink capacity, and all natural patches within the administrative district are also considered SPAs for the absorption of human carbon emissions [45]. This approach does not apply to this study area due to the relatively small size and concentration of natural patches.
Mature natural wetland systems are enormous carbon reservoirs, with 50 times more peatland carbon storage than terrestrial ecosystems globally [43]. Coastal wetlands and their ESs in large river deltas show sensitive responses to anthropogenic disturbances such as hydrological engineering, seawall construction, and land reclamation [46,47]. When wetlands are disturbed, peatlands cease to be a carbon reservoir and act as a source of carbon [48]. Therefore, wetland patches that are subject to high levels of human disturbance should be identified and protected, both for large-scale climate protection and for the prioritization of ecological elements for management at the local scale. Two indicators, therefore, influence the mapping of demand for this service: human disturbance and natural capacity to sequester carbon.

Water Purification
Natural supply potential: water purification means improving water quality, including the deposition, filtration, adsorption, biological absorption, and biochemical transformation of pollutants in sewage. Pollutants in the upper reaches of estuaries and near shore mainly come from agricultural fertilizers, pesticides, and industrial effluents. Coastal wetlands can partially decompose and eliminate the nutrients of these pollutants and serve the function of water purification. The actual capacity of the ecosystem to remove nitrogen or pollutants cannot be measured [49]. Because the removal rate of nitrogen and phosphorus by wetland vegetation was significantly correlated with vegetation nitrogen and phosphorus accumulation [50], this study used the relative amount of nitrogen and phosphorus accumulation in each land use to measure the water purification capacity and used NDVI to correct the vegetation purification capacity index.
Demand reflection: Several rivers in the study area are downstream sections into the sea, and the pressure on surface water quality assessment is high. Total P was chosen as the pollution indicator for this study in combination with the main pollution indicators for surface water communicated by the Chinese Ministry of Ecology and Environment and the availability of nitrogen-and phosphorus-related data. Many studies have used the catchment area as a physiographic unit for statistics and analysis of water-related ecosystems [14,44,51]. Like carbon sequestration, water purification service flows are in situ, directional, and all-encompassing but are relatively limited in extent. The spatial flow area for water purification services in this study was identified as the catchment area where the river network outlet is located. This service flow allocation mechanism is a complex exponential function related to topographic indices and filtration efficiency [52,53], and the Nutrient Delivery Ratio (NDR) module of the Invest model can help to assess the nutrient retention by natural vegetation as it moves from nutrient sources to streams.

Landscape Aesthetics
Natural supply potential: Landscape aesthetics is mainly reflected in the recreational, cultural, and artistic aspects of landscape evaluation. This study used the ESs value equivalent scale method to assess the capacity of this service [36,54]. Its capacity values are related to land use type, food production, and value. In addition, a correction is needed considering that there is some error in characterizing regional features with a national parameter scale. As tourism can indirectly reflect the aesthetic value of the landscape, the aesthetic value was corrected for tourism income.
Demand reflection: According to Maslow's needs theory, people need access to natural spaces to satisfy their physical, cognitive, and aesthetic needs [55]. We used the population size of the settlement and the distance to the nearest natural space to express the outdoor recreation demand index. Based on the travel patterns of residents, 5 km is considered the furthest range that people are willing to travel for daily recreation and short trips. As demand is proportional to the likelihood of people's activities and the proximity of surrounding natural spaces [56], this study linearly assigned demand to natural spaces based on population size and distance from natural spaces within 5 km.

Ecological Sources
Ecological sources are one of the core ecological elements of an ESP and need to have high quality ecological functions and provide sustainable ESs for human well-being [57]. The former is expressed as the natural supply potential of ESs, and the latter is reflected in the mapping of demand. Due to conservation costs and efficiency, the top 20% of patches for each of the two indicators were extracted, and the ensemble of the two was defined as the ecological source (Equation (1)) [58]. In particular, the intersection of the two indicators was chosen for the ecological source of this service, as the natural potential of habitat quality cannot be distinguished in the calculation method from the extent to which reed habitats are adapted to biodiversity. Patches less than 20 ha in size were excluded, taking into account the natural condition of the study area and the cost of conservation.
where Supply i > 1/5 refers to the patches with the first 1/5 score of the layer of the i-th ES supply potential; Demand i > 1/5 refers to the natural spaces with the first 1/5 score reflected by the i-th ES demand; and Area means the ecological source area, with Area min for the minimum of the source area.

Ecological Corridors
Ecological corridors are located between different ecological sources and enhance the combined supply of biodiversity and ES, providing social and economic benefits to humans [59]. The resistance surface is based on a combination of the assignment of ecological resistance to land use types and artificial facilities and the extent of their influence (Table 3) [60,61]. It reflects the state and trends of physical movement, indicating the resistance that species need to overcome in their flow from source to other spaces [62].
Referring to [10], we further corrected the resistance values for land type using hydrological connectivity, which represents a natural attribute, and NLI, which is a social attribute (Equation (2)) [60]. The relationship between biodiversity and hydrological connectivity has been demonstrated in macrobenthos and aquatic organisms [63,64]. In this paper, hydrological connectivity is expressed as an index of inundation probability [65]. As a common practice for extracting corridors, we used the MCR model to identify least-cost pathways as potential ecological corridors based on a new resistance surface (Equation (3)). Next, potential corridors were optimized to determine how much they contribute to network connectivity by evaluating the overall connectivity change caused by removing them [66]. This is defined by the following expressions: where R is the original resistance surface layer; R' is the modified resistance surface; HC is the hydrological connectivity; NLI is the nighttime light intensity; n is the total number of core patches; MCR is the minimum cumulative resistance value; D ij is the spatial travel distance from the geometric centroid of source i to target j; R i denotes the degree of impediment to species migration to destination i or landscape dispersal; f denotes the positive correlation function between MCR and ecological processes; a i and a j are the areas of patches i and j, respectively; A L is the total area of the landscape; p ij is the maximum product probability of all paths between patches i and j; dPC i is the individual importance of a certain corridor; and PC i -remove is the PC index after its removal.

Supply Potential and Demand Distribution of Ecosystem Services
The spatial distribution of habitat quality services in the Liao River Delta is shown in Figure 3a. Overall, the regional habitat quality is good, with a mean value of 0.62. The high and low values are distinguished, with the high values concentrated in the natural patches in the central and southwestern parts of the study area, mainly in coastal wetlands, marshy wetlands, and water bodies. The low values are primarily located in the western part of the study area, mainly in residential, industrial, and agricultural land use. The current distribution of indicator species is estimated in Figure 3b. A higher probability of occurrence indicates a more substantial conservation need. Demand mapping was concentrated in the central part of the natural cover, mainly sandbars and reed wetlands near the estuary. In addition, bird record data indicated more indicator species sightings in the northeastern paddy fields, which may be caused by birds flying from their roosts to the fields to feed.

Supply Potential and Demand Distribution of Ecosystem Services
The spatial distribution of habitat quality services in the Liao River Delta is sho Figure 3a. Overall, the regional habitat quality is good, with a mean value of 0.6 high and low values are distinguished, with the high values concentrated in the n patches in the central and southwestern parts of the study area, mainly in coasta lands, marshy wetlands, and water bodies. The low values are primarily located western part of the study area, mainly in residential, industrial, and agricultural lan The current distribution of indicator species is estimated in Figure 3b. A higher proba of occurrence indicates a more substantial conservation need. Demand mapping wa centrated in the central part of the natural cover, mainly sandbars and reed wetland the estuary. In addition, bird record data indicated more indicator species sightings northeastern paddy fields, which may be caused by birds flying from their roosts fields to feed.  The total supply of carbon storage service was 1.93 × 10 6 t, with a mean value of 1.24 t/unit. Overall, it shows the characteristics of high in the west and low in the east and high in the north and low in the south. The high values are concentrated in the reed wetland in the west of the study area, where the realization of carbon storage service depends on the accumulation of plant organic matter, and the high carbon density in the marsh wetland is more favorable to the accumulation of organic matter. The low values are mainly located in the southern part of the study area, where the vegetation cover is low, such as wasteland, water areas, and construction land. This area has low vegetation cover, low carbon density, and poor carbon storage capacity. The higher conservation demand for carbon sequestration services is interpreted as the greater the degree of disturbance and the greater the possibility of patches becoming carbon sources. High values are mainly distributed in the northwest and southeast near human habitation or artificial facilities.
The natural services for water purification are concentrated in the central and southwestern riverine reed patches in a banded pattern. The low values are the built-up areas and agricultural land in the east and the natural land with low vegetation cover in the south. From the mapping of water purification demand, it can be seen that high external and low internal values characterize the overall distribution. The high values are evenly distributed in several catchments in the east and northwest, mainly related to the distance between the catchment and the pollution source, the pollution load, and the slope.
The total and mean landscape aesthetic values were CNY 7.36 × 10 9 and 4794.65 CNY/hm 2 , respectively. The high values were mainly concentrated in the reed area and southern watershed in the western part of the study area, where the water system is well developed and the vegetation cover is high, which can bring high landscape aesthetic values. The mapped distribution of demand for natural space is shown in Figure 3h. The high values are mainly located in the reed and inland salt vegetation patches in the southeast and northwest.

Spatial Distribution of Ecological Elements
The total number of ecological source sites was 47, with 172.66 km 2 , accounting for 12% of the study area. The distribution within the study area was relatively balanced, with large strip patches in the southern coastal and estuarine areas. The proportion of area in the core protection zone was almost the same as the study area at similar scales. This facilitates the development of prioritized interventions to reduce costs and improve conservation efficiency. The ecological source was a collection of high-value areas for ecological functions, with the four ES core areas of habitat quality, carbon sequestration, water purification, and landscape aesthetics accounting for 4%, 25%, 11%, and 60% of the area, respectively.
The spatial distribution of the initial resistance surface is shown in Figure 4, showing the characteristics of several large settlements in the north and southeast as the center of high values and slowly decreasing to the periphery. According to the related study, the values of NLI above 300.53 and below 0.24 were set to 300.53 and 0, respectively, to exclude the data singularity values. The high values of hydrological connectivity frequency were mainly located in the central and southern water systems, artificial water bodies, and flooding areas, which can provide a suitable living environment for birds, plankton, and benthic organisms. As shown in Figure 4, the resistance surface corrected by NLI and hydrological connectivity frequency improved in both global and local accuracy. Relatively, grids with higher hydrological connectivity frequencies have lower resistance values for the same level of human disturbance.  Based on the ecological sources and improved resistance surfaces, 128 potential e logical corridors were identified (Figure 5a). The potential corridors were in a tic-tac network, mainly concentrated in the west and southeast, and did not cross the built area in the east. Their total length was 716.31 km, accounting for 258.33% of the study a boundary length. Some corridors have a high degree of overlap and a high network d sity in the northwest, with a high degree of network redundancy. The total numbe corridor networks optimized by connectivity was 54, with a total length of 145.74 Based on the ecological sources and improved resistance surfaces, 128 potential ecological corridors were identified (Figure 5a). The potential corridors were in a tic-tac-toe network, mainly concentrated in the west and southeast, and did not cross the built-up area in the east. Their total length was 716.31 km, accounting for 258.33% of the study area boundary length. Some corridors have a high degree of overlap and a high network density in the northwest, with a high degree of network redundancy. The total number of corridor networks optimized by connectivity was 54, with a total length of 145.74 km, accounting for 52.56% of the study area boundary. The optimized ESP is shown in Figure 5b.

Discussion
In the case study of the Liao River Delta, we applied the spatial flow approac connect the service supply and demand areas effectively. The vital ecological element monitoring, control, and restoration were identified through the ESP framework to f an overall ecological protection network. The critical process of ecological sources ext tion, framework application, and uncertainty still needs to be discussed next.

The Way of Selecting Ecological Sources
A set of interconnected interregional human and natural systems constitute a motely coupled system through the flow connections [67]. In previous studies, the i of connectivity between supply and demand was often neglected, incorrectly narrow the coupled system's scope. In this study, the natural supply and demand mapping a of the ES together represent the natural system, which is remotely coupled to the hum system [14]. That is, the ability to deliver services to local demand is maintained, and potential to deliver services to remote and longer time scales is retained.
The consistency of the natural supply and demand mapping zones of the ES answ whether natural supply can meet human demand and whether human well-being ca met by protecting only ecologically important zones, i.e., high natural supply zones. determines how we should choose ecological sources. For example, in landscape aest ics and water purification services, there is little overlap between high supply and h demand mapping zones. This means that both zones need to be protected simultaneo if we are to achieve the pre-determined goals of coupling. Habitat quality services carbon sequestration services are another way of thinking about source site selection. high supply and high demand mapping areas for habitat quality services are exten due to the large wetland areas and signs of biological activity in the study area. In cont carbon sequestration services generate demand only when the supply area is distur due to the unique properties of wetland carbon sinks. Therefore, the core patches of th two services are the intersection of high natural supply area and high demand mapp

Discussion
In the case study of the Liao River Delta, we applied the spatial flow approach to connect the service supply and demand areas effectively. The vital ecological elements for monitoring, control, and restoration were identified through the ESP framework to form an overall ecological protection network. The critical process of ecological sources extraction, framework application, and uncertainty still needs to be discussed next.

The Way of Selecting Ecological Sources
A set of interconnected interregional human and natural systems constitute a remotely coupled system through the flow connections [67]. In previous studies, the issue of connectivity between supply and demand was often neglected, incorrectly narrowing the coupled system's scope. In this study, the natural supply and demand mapping areas of the ES together represent the natural system, which is remotely coupled to the human system [14]. That is, the ability to deliver services to local demand is maintained, and the potential to deliver services to remote and longer time scales is retained.
The consistency of the natural supply and demand mapping zones of the ES answers whether natural supply can meet human demand and whether human well-being can be met by protecting only ecologically important zones, i.e., high natural supply zones. This determines how we should choose ecological sources. For example, in landscape aesthetics and water purification services, there is little overlap between high supply and high demand mapping zones. This means that both zones need to be protected simultaneously if we are to achieve the pre-determined goals of coupling. Habitat quality services and carbon sequestration services are another way of thinking about source site selection. The high supply and high demand mapping areas for habitat quality services are extensive due to the large wetland areas and signs of biological activity in the study area. In contrast, carbon sequestration services generate demand only when the supply area is disturbed due to the unique properties of wetland carbon sinks. Therefore, the core patches of these two services are the intersection of high natural supply area and high demand mapping.

Ecological Management Logic and Application Based on This ESP Framework
Many ESP studies based on supply and demand often calculate the ratio or coordination between natural supply and macro demand (GDP, population) to determine the degree of supply and demand matching within each larger basic unit and thus make recommendations for ecological management [5,68]. They ignore the objective fact that supply and demand are remotely coupled, and it is difficult to provide concrete knowledge for environmental management. As shown in Figure 6(1-3), the ESP further integrates ecological elements representing demand by means of service flows that allocate demand to natural space, thus improving the coupled system. ESP is a multi-element network system consisting of high natural supply areas and high demand mapping areas and corridors that can provide an important reference point for spatial policy implementation to protect nature and natural resources [2]. In principle, all development and construction activities are prohibited in the high natural supply areas, and the original land functions are permanently preserved, with planned monitoring, protection, and restoration. The high demand mapping area needs to retain both a strong natural supply capacity and the high intensity of pressure exerted by humans. The area further needs to determine whether the natural capacity is meeting human demand, such as whether surface water quality is up to standard and whether people's recreational needs are being met qualitatively and quantitatively. It is also essential to confirm whether human demand disturbance in the region excessively affects the regular natural supply, such as whether the high demand reflecting area for carbon sequestration is undergoing a functional shift from a carbon sink to a carbon source. These findings will help us select appropriate prioritization interventions to regulate the supply capacity from the natural system and the demand from the social system, respectively.
For this study area, the Pearson coefficients of demand mapping and supply potential are moderately negatively correlated, indicating that the two are not coordinated in growth trends. This suggests that further action should be taken in the demand mapping area to ensure that natural supply can meet human demand. Combined with the results in Figure 6(4a), the natural patches of Liaodong Bay New Town in the southeast and the ecological economic zone in the northwest face high human demand and low supply capacity for landscape aesthetics, and the supply capacity of the SPA,. such as improving the connectivity of these areas to increase the distribution of service flows along the route, expanding flow boundaries, and conducting green infrastructure construction, needs to be further strengthened. An urban central park can be built in Liaodong Bay New Town by relying on the existing large-area natural patches. In addition to the natural patches of ecological sources, community parks will be built near large residential areas in the ecological economic zone, and traffic connections will be ensured. Figure 6(4b) shows that the natural supply of regional habitat quality is sufficient, and there is a pointwise mismatch between supply and demand in the central Red Beach Scenic Area. Ecological restoration measures such as reducing the hydrological impact caused by the construction of upstream agricultural canals, restricting oilfield exploitation activities in nature reserves, and providing habitats on the leeward side are needed. Figure 6(4c) shows that the natural supply of water purification is sufficient, and the natural patches in the southeastern port area face tremendous purification pressure, which can be alleviated by reducing pollution discharge or enhancing green infrastructure construction in conjunction with the Liaodong Bay New City Central Park. In addition, it should strictly implement the return of farmland to wetlands and accelerate the shutdown and withdrawal of oil well facilities in the ecological source of carbon sequestration services.
human behavioral tendencies and water-related models for pathway allocation. We simplify the complex service flow between the demand and supply areas by taking the conservation target as the demand subject for the demand mapping of habitat quality and carbon sequestration services. Using black box theory to map demand to natural space provides the possibility of incorporating some service types with complex demand allocation processes into the ESP framework [42,70]. Figure 6. Ecological management logic and application: (1a) protection of habitat patches indicative of the physical presence of birds; (1b) protection of carbon sinks that are undergoing severe disturbance; (1c) protection of patches that perform the function of purifying pollutants emitted by humans that flow into water systems; (1d) attention to natural patches within people's living circles; (2) Connecting natural-social systems through ecosystem service flows; (3) A part of the high supply area, the actual supply area and corridors together form an overall ecological safety network; (4) pairwise comparisons were made between the natural service potential and demand-reflecting values of different ecosystem services with SBA as the boundary; (4a) landscape aesthetics; (4b) habitat quality; (4c) water purification.

Limitations and Future Research
Although this study constructs an ESP based on landscape-scale supply and demand to strengthen the role played by ecological source sites in coupling human-land systems, Figure 6. Ecological management logic and application: (1a) protection of habitat patches indicative of the physical presence of birds; (1b) protection of carbon sinks that are undergoing severe disturbance; (1c) protection of patches that perform the function of purifying pollutants emitted by humans that flow into water systems; (1d) attention to natural patches within people's living circles; (2) Connecting natural-social systems through ecosystem service flows; (3) A part of the high supply area, the actual supply area and corridors together form an overall ecological safety network; (4) pairwise comparisons were made between the natural service potential and demand-reflecting values of different ecosystem services with SBA as the boundary; (4a) landscape aesthetics; (4b) habitat quality; (4c) water purification.
The demand mapping in this study is an extension and development of the existing studies, which can provide reference and inspiration for ESP studies in other study areas [4,14,42,69]. The technical paths regarding the allocation and mapping of demand for landscape aesthetic and water-related services are relatively clear. Their demand boundaries (e.g., travel willingness range and water catchment area) can be easily accessed through behavioral and spatial analyses. More established studies provide references to human behavioral tendencies and water-related models for pathway allocation. We simplify the complex service flow between the demand and supply areas by taking the conservation target as the demand subject for the demand mapping of habitat quality and carbon sequestration services. Using black box theory to map demand to natural space provides the possibility of incorporating some service types with complex demand allocation processes into the ESP framework [42,70].

Limitations and Future Research
Although this study constructs an ESP based on landscape-scale supply and demand to strengthen the role played by ecological source sites in coupling human-land systems, it still has limitations and uncertainties. First, this paper used bird observation data and precipitation data to assess the mapping of habitat quality and water purification demand. These data themselves may have uncertainties. For example, the number of bird observations is uploaded to the system by a large number of users, and the precipitation distribution is interpolated from the observed values of multiple stations, which may be somewhat different from the actual values in a specific space. Second, in the assessment of the natural supply of carbon sequestration, the assessment model relies on the carbon stock estimation for each land use type, so we have averaged the sampled data for different land use types. This may be more reasonable in assessing the overall carbon pool, but the highly generalized rank may lead to errors and loss of some spatial information, thus affecting the extraction of ecological sources. Datasets with a large number of sample points and long time spans could be used in future studies to compensate for this deficiency.
Third, in this paper, the threshold settings for corridor and source site extraction take the approach of empirical values and expert judgment, which is somewhat subjective. We did not put a limit on the maximum length of corridors because there is no research showing that ecological corridors have a distance limit in transporting materials and energy for ecosystems. The flow properties of key species at different scales should be included in future studies. Finally, we included ecological source sites jointly for high natural supply and high demand mapping areas and discussed their consistency to ensure the high value of ecological source sites. However, there are similar problems to other ESP studies. The performance and accuracy of ESP cannot be directly verified due to the limitations of controlled experiments [5]. This content will become necessary in the future when sufficient practice samples are available.

Conclusions
This study comprehensively considers the natural supply potential of the ecosystem and human needs and constructs an ecological security pattern with the characteristics of coastal wetlands to maintain the integrity of the ecosystem structure, function, and process. The service supply area and the service benefit area are connected by the method of spatial flow to identify the demand reflecting area and make up for the problem of remote coupling of supply and demand. This study proposes a new ESP framework in which various ecological elements can provide decision support on paths and locations for ecological management, safeguard the demand for natural products within the region, and maintain the natural potential for longer-term and wider-scale services. The frame is tried in the Liao River Delta. Based on understanding the processes by which multiple spatial streams of services transport and deliver ES, supply potential and service demand are quantified to identify ecological source sites consisting of high natural supply areas and high demand mapped areas. Based on the resistance surface corrected by hydrological connectivity and NLI, potential ecological corridors were extracted, and then network connectivity was optimized through connectivity evaluation. This study proposes ecological management logic and application according to the multi-level ESP and supply-demand relationship.