Integrating Ecological Assessments to Target Priority Restoration Areas: A Case Study in the Pearl River Delta Urban Agglomeration, China

: The identiﬁcation and management of ecological restoration areas play important roles in promoting sustainable urban development. However, current research lacks a scientiﬁc basis for the scope and scale of ecological restoration. Further, the absence of a framework to assess policy goals and public preferences that leads to identiﬁcation of ecological restoration areas across the science-policy interface is difﬁcult, and the existing frameworks’ performance has little applicability. We proposed a transdisciplinary framework to combine ecological quality, ecological health, and ecosystem services as an assessment endpoint to identify priority restoration areas. Further, we classiﬁed the ecological restoration areas on a township scale by K-means. Based upon policy goals and public preferences of the Pearl River Delta urban agglomeration, we chose air quality, biodiversity, soil fragility, recreation quality, ecosystem vigor, landscape metrics, and the water supply ecosystem service as elements of the evaluation system. This study showed that priority restoration areas accounted for 10.8% of the urban agglomeration area and classiﬁed township, largely in the difference between natural and semi-natural ecosystems and the human environment. Policymakers can use this framework comprehensively and ﬂexibly to identify and classify ecological restoration areas to achieve policy goals and fulﬁl public preferences.


Introduction
Human activities have significant adverse effects on the regional ecological environment that result in climate change, loss of biodiversity, and ecosystem degradation [1]. In developing countries in particular, urbanization has resulted in changes in land use, encroachment on ecological areas, and destruction of primitive habitats [2], all of which have led to the decline in regional ecosystems' quality and health, and varying degrees of degeneration. Approximately 2 billion hectares of the world's terrestrial ecosystems are considered degraded and in need of ecological restoration [3,4]. These ecosystems' decrease in area and widespread degradation has reduced ecosystem services and increased ecological risk [5]. With the promotion of sustainable urban development, the ecosystem's protection and restoration have become prominent research topics that focus on assessing ecosystem services for urban planning [6][7][8]. Assessing ecological space is dominated largely by social and economic activities, particularly policy-oriented and decision-making needs [9,10]. Thus, comprehensive research on the application of scientific assessment of ecological space to help make decisions and policies should be strengthened [11,12].

Study Area
The PRD urban agglomeration located in southern Guangdong province, China, consisting mainly of nine cities: Guangzhou (GZ), Dongguan (DG), Shenzhen (SZ), Huizhou (HZ), Foshan (FS), Zhongshan (ZS), Zhuhai (ZH), Jiangmen (JM), and Zhaoqing (ZQ) (Figure 1a). The PRD covered an area of 54,770 km 2 inhabited by approximately 63 million people in 2018 (Statistical yearbook of Guangdong Province). With the rapid development of urbanization, the ecological problems and ecological risks were increasing significantly, and regional ecological security was facing certain threats [7]. The PRD is currently concerned with important ecological restoration planning issues, with a political leadership assuming the ambition to reach sustainable development by 2035. To reach this objective, the most important thing is to identify the ecological restoration area and implement ecological restoration zoning and hierarchical control. Currently, policymakers need more information about ecological restoration targets and ecological restoration locations to achieve policy goals. A major challenge is developing capacities on ecological restoration identification for urban planning, which is difficult because ecological assessment requires new interdisciplinary methodologies to assess multiple ecological space and social concerns.
Ecological space is a collection of natural, artificial, and semi-artificial ecological units [28]. Regions that play an important role in the ecosystem can be regarded as ecological spaces. They include not only forestland, grassland, wetland, water bodies, and other natural ecological spaces, but also cropland, parks, ecological corridors, and other areas that provide ecosystem services [13,29]. According to the definition of ecological space and land use data, the ecological space in this paper was divided into five broad types: forest land, cropland, grassland, water, and unused land [13] (Figure 1b). The ecological space data were derived from the Resource and Environment Science Data Center of the Chinese Academy of Sciences (http://www.resdc.cn/, accessed on 13 June 2020), with a spatial resolution of 30m. Landsat 8 remote sensing images were used as data sources to interpret ecological space data obtained in 2018. Through random sampling inspections, the accuracy rates of land-use classifications were greater than 94% (Table S1), which can be used to explain the land changes [10,27].
Remote Sens. 2021, 13, 2424 3 of 14 ecological assessment and ecological management to cope with the complexity and variability of the real world.

Study Area
The PRD urban agglomeration located in southern Guangdong province, China, consisting mainly of nine cities: Guangzhou (GZ), Dongguan (DG), Shenzhen (SZ), Huizhou (HZ), Foshan (FS), Zhongshan (ZS), Zhuhai (ZH), Jiangmen (JM), and Zhaoqing (ZQ) (Figure 1a). The PRD covered an area of 54,770 km 2 inhabited by approximately 63 million people in 2018 (Statistical yearbook of Guangdong Province). With the rapid development of urbanization, the ecological problems and ecological risks were increasing significantly, and regional ecological security was facing certain threats [7]. The PRD is currently concerned with important ecological restoration planning issues, with a political leadership assuming the ambition to reach sustainable development by 2035. To reach this objective, the most important thing is to identify the ecological restoration area and implement ecological restoration zoning and hierarchical control. Currently, policymakers need more information about ecological restoration targets and ecological restoration locations to achieve policy goals. A major challenge is developing capacities on ecological restoration identification for urban planning, which is difficult because ecological assessment requires new interdisciplinary methodologies to assess multiple ecological space and social concerns.
Ecological space is a collection of natural, artificial, and semi-artificial ecological units [28]. Regions that play an important role in the ecosystem can be regarded as ecological spaces. They include not only forestland, grassland, wetland, water bodies, and other natural ecological spaces, but also cropland, parks, ecological corridors, and other areas that provide ecosystem services [13,29]. According to the definition of ecological space and land use data, the ecological space in this paper was divided into five broad types: forest land, cropland, grassland, water, and unused land [13] (Figure 1b). The ecological space data were derived from the Resource and Environment Science Data Center of the Chinese Academy of Sciences (http://www.resdc.cn/, accessed on 13 June 2020), with a spatial resolution of 30m. Landsat 8 remote sensing images were used as data sources to interpret ecological space data obtained in 2018. Through random sampling inspections, the accuracy rates of land-use classifications were greater than 94% (Table S1), which can be used to explain the land changes [10,27].

Selection of Ecological Quality and Ecological Health Indicators
Ecological restoration is the result of multiple trade-offs. Ecological restoration areas have been defined as land that has experienced persistent decline or loss in biodiversity and ecosystem functions and services that cannot recover fully unaided within decadal time scales [30,31]. Considering the sustainable development of cities, this study defined ecological restoration areas as low-quality or high-risk areas, which represent lower values Remote Sens. 2021, 13, 2424 4 of 14 of regional biodiversity, ecosystem functions, and services that needed to be improved [32]. The assessment indicators should be selected based upon differing knowledge systems, such as policy goals, public preferences, expert knowledge, environmental quality, and successful restoration engineering [33,34]. The ecological goals of urban agglomeration planning and environmental protection planning are described in such documents as the "Planning for restoration of land in Guangdong (2020-2035)". This proposes to divide ecological space with important ecological functions or fragile ecosystems or natural ecological protection into an important ecological management space. "The Development planning of the Guangdong-Hong Kong-Macao Greater Bay Area" suggested further that urban planning should guarantee water resources' safety, and build ecological corridor and biodiversity protection networks. The policy goals were to maintain the natural ecosystems' integrity and health as a starting point, expand high-quality ecological areas, and safeguard the human environment's health and safety. Based upon this, the methodology adopted was based on participatory workshops involving policymakers and researchers. All of them have a good knowledge of the study area and the ecological restoration, through their involvement in field-based research or restoration activities. We found that policymakers and researchers hope to enhance biodiversity, air quality, water supply, public environmental health, and regional ecological security through the restoration of ecological space. To meet these requirement, ecological quality and health were integrated to represent the policy goals and public preferences to identify the ecological restoration areas ( Figure 2) and then identify the imbalance in regional ecological development and analyze the ecological areas that need restoration to meet policy goals and public preferences. These requirements may be met by selecting hotspot areas of ecological quality and health [14]. Incorporating the different types of hotspots into spatial planning can help policymakers draft more comprehensive strategies to improve the ways to balance environmental protection and urban development [14,35]. In this study, policy objectives, expert knowledge, and public preferences were combined to determine the choice of hotspots. Ecological areas that ranked as those that provided less than 10% of ecological quality and health were prioritized as ecological restoration areas.

Selection and Calculation of Ecological Quality Indicators
We identified a preliminary assessment index system based on participatory workshops. Then, questionnaires were administered to the residents of the PRD urban agglomeration to determine their preferences for the ecological space. We asked the respondents what kind of ecological qualities concern them generally (Select at least one item or choose three items at most): (1) Biodiversity; (2) vegetation coverage; (3) water network density; (4) river water quality; (5) soil fragility; (6) air quality; (7) outdoor recreation quality; (8) landscape connectivity, and (9) other. Random sampling was used in the survey. A total of 500 questionnaires were distributed, and 361 valid questionnaires were collected, including 204 questionnaires in the urban community and 157 questionnaires in the rural community (Table S2). Respondents mainly included farmers, employees, retirees, and whose preferences cover all or part of the urban agglomeration. More than 50% of respondents preferred biodiversity (BD), outdoor recreation quality (OR), air quality (AQ), and soil fragility (SF) ( Figure S1). This reflected residents and stakeholders focused largely on meeting their leisure and recreational needs, as well as maintaining a good living and working environment with forest coverage, high air quality, and park accessibility. Finally, these indices were selected as ecological space quality (EQ) assessment indicators, and the weight of each was determined according to the Analytic Hierarchy Process. The calculation formula was as follows:

Figure 2.
A multi-disciplinary framework to identify ecological restoration areas by linking policy and ecological information. Identification of ecological restoration areas depends upon their ecological quality and health (black dotted bordered rectangle). The assessment of ecological quality and health should incorporate policy goals, local stakeholder needs, and expert knowledge (red dotted bordered rectangle) to select the desired index to ground the analysis in the given decision-making context to achieve equal and full representation. NPP: Net primary productivity.

Selection and Calculation of Ecological Quality Indicators
We identified a preliminary assessment index system based on participatory workshops. Then, questionnaires were administered to the residents of the PRD urban agglomeration to determine their preferences for the ecological space. We asked the respondents what kind of ecological qualities concern them generally (Select at least one item or choose three items at most): (1) Biodiversity; (2) vegetation coverage; (3) water network density; (4) river water quality; (5) soil fragility; (6) air quality; (7) outdoor recreation quality; (8) landscape connectivity, and (9) other. Random sampling was used in the survey. A total of 500 questionnaires were distributed, and 361 valid questionnaires were collected, including 204 questionnaires in the urban community and 157 questionnaires in the rural community (Table S2). Respondents mainly included farmers, employees, retirees, and whose preferences cover all or part of the urban agglomeration. More than 50% of respondents preferred biodiversity (BD), outdoor recreation quality (OR), air quality (AQ), and soil fragility (SF) ( Figure S1). This reflected residents and stakeholders focused largely on meeting their leisure and recreational needs, as well as maintaining a good living and working environment with forest coverage, high air quality, and park accessibility. Finally, these indices were selected as ecological space quality (EQ) assessment indicators, and the weight of each was determined according to the Analytic Hierarchy Process. The calculation formula was as follows: Figure 2. A multi-disciplinary framework to identify ecological restoration areas by linking policy and ecological information. Identification of ecological restoration areas depends upon their ecological quality and health (black dotted bordered rectangle). The assessment of ecological quality and health should incorporate policy goals, local stakeholder needs, and expert knowledge (red dotted bordered rectangle) to select the desired index to ground the analysis in the given decision-making context to achieve equal and full representation. NPP: Net primary productivity.
In view of the experimental data's availability and spatial representation, this study referred to previous studies [13] that established the basic equivalent of different types of ecosystem service functions' value based upon expert experience to determine the land use types' biodiversity (Table S3). The air quality index derived from the Environmental Data Center of the Ministry of Environmental Protection of China (http://datacenter.mep. gov.cn, accessed on 5 April 2019). The ESTIMAP-recreation method was adopted for residents' outdoor recreation [36]. The score of each grid was calculated by summing three components (degree of naturalness, nature protection, and water) ( Table S4). The soil fragility index was weighted by soil erosion using the Ministry of Environmental Protection of China's habitat quality index as the standard. Soil erosion originated from the Resource and Environment Science Data Center of the Chinese Academy of Sciences (http://www.resdc.cn/, accessed on 17 March 2019) ( Table S5). The indices were all normalized to 0-100 according to their positive or negative effects.

Selection and Calculation of Ecological Health Indicators
The selection of ecological health indicators was based upon the ecosystem's ability to maintain a healthy structure, self-regulate, and recover under pressure, and was characterized synthetically by vigor, organization, and resilience [25]. Vigor represents the ecosystem's primary productivity. In this study, net primary productivity (NPP) was used to evaluate ecosystem vigor, and NPP data (MOD17A3) were acquired from http://ladswed.nascom.nasa.gov [37], accessed on 8 December 2019. The organization was determined primarily by landscape pattern [18], landscape connectivity, landscape heterogeneity, forest land, and wetland landscape characteristics that were characterized synthetically; their weight was determined according to previous studies [25].
in which O is ecosystem organization, AWMPFD is the index of the area-weighted mean patch fractal dimension, FN 1 is landscape fragmentation index, SHDI is Shannon's diversity index, MSIDI is modified Simpson's diversity index, CONT is the landscape contagion index, FN 2 is the fragmentation index of forest land, CONNECT 1 is the index of the patch connect of forest land, FN 3 is the fragmentation index of wetland, and CONNECT 2 is the index of the patch connect of wetland. These indicators were obtained using the moving window method in Fragstats v. 4.2 (the University of Massachusetts); the landscape index ranges from 0 to 100, and the unit is the percentage. Ecosystem resilience refers to the regional ecosystem's ability to restore itself to its original structure and function after external disturbance. This was determined largely by service indicators that regulate ecosystems' important ecological processes and functions [38]. In this study, the ecosystem services' regulatory services were used to represent the policy objectives of the water supply's resilience, which reflected the ecosystem's ability to conserve water, the most important ecosystem service in the PRD. The ecosystem service indicators above were calculated using the InVEST v. 3.3.3 model [14,39].

Identification of Management Bundles via a K-Means Cluster Analysis
Classifying ecological restoration areas helps identify different types of restrictive ecological elements in the areas and furthers the implementation of ecological management measures. Through a trade-off and synergistic analysis, air quality was found to be correlated significantly only with resilience. Based upon this, we eliminated the air quality index in the selection of classification indices of ecological restoration areas. To formulate the ecological policy and manage the ecological space for policymakers, we took the township as the statistical unit. To identify ecological restoration bundles, we selected the K-means clustering consistent with other similar studies [17,40]. Thereafter, we used the NbClust package in the R statistical software (v. 4.0.3) to choose the optimal number of K-means clusters [41,42].

Spatial Identification of Ecological Restoration Areas
Incorporating the different types of hotspots into spatial planning can help policymakers draft more comprehensive strategies to balance environmental protection and urban development. Ecological areas that demonstrated the lower 10% of ecological quality and health were identified as restoration hotspots (Figure 3). By choosing the hotspot areas according to ecological space quality, a land area of 3692.5 km 2 was selected as the ecological quality hotspot, which was distributed primarily in the cropland that accounted for 99.5% of the built-up land (Figure 4a). A land area of 2632.7 km 2 was selected as the ecological health hotspot, which was distributed largely in cropland and water bodies (mainly ponds) and accounted for 52.7% of the area (Figure 4b). Finally, we superimposed the ecological quality and health hotspots synthetically, excluding 436.9 km 2 of overlapping areas, and the potential ecological restoration areas totaled approximately 5888.3 km 2 that accounted for 10.8% of the PRD's urban agglomeration area (Figure 4c). With respect to the total amount, the restoration areas were concentrated primarily in ZQ, FS, and HZ. However, in terms of unit area, the restoration areas were concentrated in FS (Figure 4d).
(mainly ponds) and accounted for 52.7% of the area (Figure 4b). Finally, we superimposed the ecological quality and health hotspots synthetically, excluding 436.9 km 2 of overlapping areas, and the potential ecological restoration areas totaled approximately 5888.3 km 2 that accounted for 10.8% of the PRD's urban agglomeration area (Figure 4c). With respect to the total amount, the restoration areas were concentrated primarily in ZQ, FS, and HZ. However, in terms of unit area, the restoration areas were concentrated in FS (Figure 4d).

Classification of Township Management in Ecological Restoration Areas
This study adopted the K-means algorithm to classify the indicators into five categories ( Figure 5), and then explained them according to land classification characteristics, and the characteristics of leading ecological indicators (Table 1). Land use types can respond directly to policies and socioeconomic effects, which is the core of ecological management [6,43]. Bundle 1 (34.6% of the restoration area, 32.7% of the urban agglomeration area) consisted largely of forestland (59.6%) and cropland (25.1%). The regional ecological indicators were relatively balanced, but the values were low. The regional human activities were primarily agriculture-based and low-level, and were found in the primary agricultural economic region. Bundle 2 (11.4% of the restoration area, 31.2% of the urban agglomeration area) was composed mainly of large areas of natural forestland (84.6%). It had a high ecosystem service function and was an important area in which to construct a regional ecological security pattern. The ecological restoration area was distributed largely around cropland and forestland. Bundles 1 and 2 were largely natural ecosystems

Classification of Township Management in Ecological Restoration Areas
This study adopted the K-means algorithm to classify the indicators into five categories ( Figure 5), and then explained them according to land classification characteristics, and the characteristics of leading ecological indicators (Table 1). Land use types can respond directly to policies and socioeconomic effects, which is the core of ecological management [6,43]. Bundle 1 (34.6% of the restoration area, 32.7% of the urban agglomeration area) consisted largely of forestland (59.6%) and cropland (25.1%). The regional ecological indicators were relatively balanced, but the values were low. The regional human activities were primarily agriculture-based and low-level, and were found in the primary agricultural economic region. Bundle 2 (11.4% of the restoration area, 31.2% of the urban agglomeration area) was composed mainly of large areas of natural forestland (84.6%). It had a high ecosystem service function and was an important area in which to construct a regional ecological security pattern. The ecological restoration area was distributed largely around cropland and forestland. Bundles 1 and 2 were largely natural ecosystems with high quality ecological space and low health risk, and were the main areas of regional ecosystem services (supply and regulation). This study found that the areas that require ecological restoration were concentrated largely in the border areas between cropland and forestland. The main cause of regional degradation was the fact that, because of new agricultural activities, cropland had invaded natural ecological space. The ecological management policies should strengthen the protection of forestland, improve existing croplands' quality, convert cropland to forestland and grassland, and establish nature reserves and ecological recreational areas to enhance the regional outdoor recreational capacity. Bundle 3 (32.5% of the restoration area, 21.4% of the urban agglomeration area) consisted primarily of cropland (38.4%) and built-up land (20.8%). This area had little urbanization and a high level of agriculture. Except for organization, the other indicators in this bundle were lower than the average value, which was the key restoration area in the region. Bundle 3 was primarily a semi-natural ecosystem found in an area with little urbanization. It was the main area in which agriculture was distributed in the urban agglomeration and also the focus of regional ecological restoration. It is necessary to adjust the agricultural planting structure, transform "simple" crops into agroforestry systems, and integrate natural processes into agricultural production. Given the regional ecological characteristics, it is recommended to build fishponds with appropriate regional characteristics. Ponds function to regulate the amount of water and increase water quality [39]. The government can consider implementing ecological compensation to encourage this practice. Bundle 4 (3.5% of the restoration area, 1.6% of the urban agglomeration area) consisted largely of water bodies (37.8%) and cropland (29%). It was distributed primarily along the coast of the Pearl River and around the reservoir, and had high outdoor recreational potential. Bundle 4 was a special area of the PRD urban agglomeration that was composed primarily of rivers and reservoirs. It has a dense river network and is a good place for outdoor recreation. However, in identifying ecological restoration areas, it was found that this area also faces encroachment by cropland or built-up land. The government should implement measures to convert cropland to water bodies and expand the water surface area, as well as engage in greening to enhance the opportunities for recreation and leisure. Particularly in the case of ponds, ecological management should change their main service from aquaculture to leisure tourism, carry out greening and beautification, enhance biodiversity, and develop ecological tourism, such as fishing, moderately. Bundle 5 (18% of the restoration area, 13.1% of the agglomeration area) was composed primarily of built-up land (39.2%) and cropland (23%) distributed in the central part of the PRD urban agglomeration in a highly urbanized area. Bundle 5 is related closely to highly urbanized areas where the expansion of urban built-up land is the main threat to the ecological space and the intensive change in land use and increased population density has led to the decline in ecosystem services and ecological health. Hence, ecological management should rely upon the PRD's good river network structure and traffic trunk lines, select appropriate local, rather than exotic plant species, build ecological corridors, small urban green spaces, diversified natural ecological spaces, and enhance the urban ecological space's outdoor leisure capacity. the ecological space and the intensive change in land use and increased population density has led to the decline in ecosystem services and ecological health. Hence, ecological management should rely upon the PRD's good river network structure and traffic trunk lines, select appropriate local, rather than exotic plant species, build ecological corridors, small urban green spaces, diversified natural ecological spaces, and enhance the urban ecological space's outdoor leisure capacity.

Bundle 1 Forestland and Cropland Cropland and built-up land encroach on forest land
Strengthen the protection of forestland; Improve existing croplands' quality; Establish ecological recreational areas.

Bundle 2
Forestland Cropland encroaches on forest land Strengthen the protection and construction of forestland; Convert cropland to forestland and grassland; Establish nature reserves.

Bundle 3 Cropland and built-up land
Built-up land encroaches on cropland Ensure quantity and quality of essential cropland; Adjust the agricultural planting structure, transform 'simple' crops into agroforestry systems; Integrate natural processes into agricultural production.

Bundle 4 Water bodies and cropland Built-up land and cropland encroach on water body
Convert cropland to water bodies and expand the water surface area; Engage in greening to enhance the opportunities for recreation and leisure.

Bundle 5 Built-up land
Built-up land encroaches on other ecological space Build ecological corridors, small urban green spaces, diversified natural ecological spaces; Enhance the urban ecological space's outdoor leisure capacity.

Discussion
Our work began by acknowledging the necessary trade-off of multiple needs, identifying the priority ecological restoration areas with hotspots of ecological quality and health, and conducting a cluster analysis of urban agglomerations in township administrative units by machine learning. Our work provides a new method to identify and rank the management of ecological restoration areas scientifically. Previous studies have assessed ecological space largely from a single perspective and focused on its spatial and temporal distribution, as well as changes in ecological risk [8,44]. Previous studies have given more attention to the function of ecological space, including physical and biological composition and its benefits. Further, the research scale in the studies has been small (e.g., landscape, community scale), which may not capture the true magnitude of the damage that biodiversity and ecosystem functions and services have suffered, such that the application of the data to urban planning and policymaking had limitations [14]. Few studies have addressed the combination of ecological policy and ecological restoration areas' identification [45]. In China, ecological restoration is led primarily by the government, which is concerned principally with the construction of regional-scale ecological security, and has insisted that a holistic approach to protecting and restoring mountains, rivers, forests, croplands, lakes, and grasslands is necessary. Based upon this concept, our research filled a knowledge gap between ecological management, restoration, and assessment, which was based upon stakeholders' general concerns, policymakers' policy goals, and the public's needs. Our results may help policymakers formulate ecological policy and implement urban planning. The methods of this study will help reveal the public preferences on ecological space under urbanization in urban agglomeration. Our model can be applied to many scales, including county, city, urban agglomeration, province, and country. Through the spatial pattern distribution, we can explore their differences and internal relations, comprehensively analyze the results of ecological restoration spatial identification at different scales, and explore their internal correlations, as well as the ways and methods of mutual conversion between different scales. However, for some small scales, such as landscape patches, because the landscape evaluation index is focused on structure and function evaluation, this model cannot objectively reflect the actual results. The index system of this study can provide an initial choice for other urban agglomerations. However, the local actual situation should be considered, then establish the index system to suit local characteristics.
The combination of ecological quality and ecological health is an ideal endpoint in identifying ecological restoration areas. Considerable previous research had assessed ecological quality and health. However, there was still no way to assess ecological space comprehensively to identify ecological restoration areas. Hence, this study adopted the vigor-organization-resilience system assessment framework to assess ecological health based upon previous studies [25,46]. However, in selecting resilience indicators, this study did not adopt area-weighted ecosystem resilience coefficients to calculate land-use types; rather, we adopted the ecosystem services indicators with which policymakers and stakeholders are concerned generally. Existing studies had incorporated ecosystem services into the assessment of ecological health risk [25], as these have always been important indicators of ecological health from the human-nature coupled view [47]. Land use changes affect the ecosystem itself primarily through the services it provides, so the assessment of ecosystem services reflects the ecosystems' flexibility and sustainability better. Ecological quality and health and ecosystem services may be considered to have a reciprocal causal relation, which was an appropriate medium to evaluate ecological space comprehensively on the regional scale.
Ecological zoning is the core method of regional ecological management [48], and it has been recognized that the classification of ecological restoration areas through scientific and adaptive evaluation methods is one of policymakers' common concerns. Particularly for urban agglomerations, different cities' development level and ecological goals differ, and the failure to identify ecological zoning from the perspective of urban agglomerations has led to confusion and inconsistent management. Further, the mismatch between ecological restoration areas and administrative boundaries has led to the formulation of irrational policy and control implementation [5,49]. Based upon the assessment indicators of ecological quality and ecological health, this study used machine learning to classify the township scale according to the ecological restoration area's characteristics. The bundle could classify towns with similar ecological conditions clearly, identify ecological shortcomings in towns, and provide a further theoretical basis for policymakers to formulate ecological management policies on the urban agglomeration scale [17,40,50].
Our case study also had some limitations. Firstly, the research framework serves only as a stage in delineating the priority ecological restoration area. To determine the final ecological restoration area, policymakers need to consider more socioeconomic factors, particularly the interests of stakeholders in the restoration area. This research provided the basic framework with which to identify the priority ecological restoration area. Secondly, to determine the hotspot area's threshold, this study considered the scale of the PRD urban agglomeration, and combined the experts' knowledge and policy targets. Thus, the threshold was based upon a multiparty trade-off and societally-informed choice. However, the hotspot area's threshold may be updated over time depending upon socioeconomic development or the development of more advanced technology, and may generally vary in another region. Thirdly, although our cluster analysis divided the urban agglomeration into five bundles, the basis of classification was primarily the ecological indicators. Hence, the analysis of the bundles' determining socioeconomic forces is still incomplete. The zoning management of ecological restoration areas should be the focus of future research on socioeconomic determinants. Establish a suitable local spatial clustering method focused on social economy and natural factors may be a way for future additional investigations of ecological restoration to meet the complex and changeable urban ecosystem.

Conclusions
Based upon policy goals and public preferences, this study proposed a comprehensive method to identify priority ecological restoration areas by integrating the ecological space's quality and health, and performing a cluster analysis of a township using machine learning. In this way, we produced useful insights about the influence of policy decisions or actions to define ecological restoration areas at the scale of urban agglomeration, as well as basic assumptions for urban planning and ecological management. This study found that the priority restoration areas were distributed primarily in the areas surrounding built-up urban areas. This can be explained in part by the current state of the urbanization, urbanized for many years without real consideration of the regional ecological quality and ecological health, reinforced by the global context of the urbanization expansion. Thus, identification of ecological restoration areas requires a strong and assumed multi-actors' governance in favor of improving ecological quality and ecological health facing past and future urbanization pressures. To achieve this, we still need to develop frameworks outlining core socioecological criteria for determining ecological restoration targets. We need clear illustrations of the main methodological steps for identifying ecological restoration areas to communicate the importance of different socio-ecological components to stakeholders. The next step is to develop society-economic-nature index datasets to link ecosystem characteristics and socioeconomic characteristics to evaluate the socio-ecological outcomes from the multiple trade-offs to refine actual ecological restoration on the ground. Future ecological restoration should be unified planning according to the multiple trade-offs with the development trend of urbanization degree and should strengthen the restoration and construction of ecological space intra-city built-up area. It can be concluded that this approach permitted to fill a gap in the identity of ecological restoration in the spheres of science and policy interdependencies, with the aim of the consideration of ecological quality and ecological health in urban agglomeration ecological management.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/10 .3390/rs13122424/s1, Figure S1. The public preference of ecological space quality. Table S1. The error matrix of classified and reference data for 2018. Table S2. Attributes of respondents. Table S3. Weight of land use types of biodiversity index. Table S4. Weight of the outdoor recreation index. Table S5. Weight of the land fragile index.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.