Identifying Spatial Priority of Ecological Restoration Dependent on Landscape Quality Trends in Metropolitan Areas

Ecological restoration has become an important tool for mitigating and adapting to environmental degradation caused by global urbanization. However, current research has focused on single indicators and qualitative analysis, meaning that ecological restoration has not been effectively and comprehensively addressed. This study constructed a spatial priority identification system for ecological restoration, with landscape area, landscape structure and landscape function as the core indicators. The system has wide adaptability. In this work, the spatial classification of ecological degradation was performed by overlay analysis. The results showed the following: (1) In the Shanghai metropolitan area, the landscape quality showed a trend of degradation, with built-up areas encroaching on forests and cropland. (2) Ecological degradation in the suburbs was more severe than that in the urban center. Forests had the highest landscape area indicator (LAI) stability. Significant degradation of landscape structure indicators (LSIs) occurred when built-up area and cropland were transformed into forests. (3) Different types of ecological restoration had significant spatial distribution patterns. Through this identification system, this study aimed to help planners/managers of ecological restoration to recognize the changing patterns of regional landscape quality and its relationship with land cover. It ultimately provides a basis for the formulation of regional ecological objectives and spatial strategies.


Introduction
In recent years, cities, as human settlements, have seen continuous expansion of their built-up areas in efforts to meet the growing needs of mankind [1]. From 2001 to 2018, the global built-up area increased from 7.47 × 10 5 km 2 to 8.0 × 10 5 km 2 [2]. With the population gathering, industrial upgrading, and the spread of rapid transportation, metropolitan areas are emerging globally as large, highly integrated regions [3]. However, humanity's growing ecological footprint has altered the original vegetation, climate systems, and material and energy cycles [4]. These changes have led to numerous urban risks, such as the urban heat island (UHI) effect, air pollution, and urban flooding, which severely limit the sustainable development of cities [5][6][7]. Therefore, in metropolitan areas, the study of landscape status, change processes, and ecological restoration has become a worldwide research hotspot [8].
Landscape quality as a term to measure spatial quality has an extensive definition that includes many factors, such as environmental pollution and cleanliness, visual and aesthetic quality, and regional ecological benefits [9]. Thus, the evaluation of landscape quality should include a broader analysis of the structure and relationships between the different landscape elements. In this study, landscape quality was defined and understood as the state in which the landscape area, structure, and function is found at a given time. This quality is an effect superimposed on a set of environmental components, processes, and phenomena that is the direct result of human activity (e.g., urbanization and spatial design) or indirect influence (e.g., agriculture and industry) [10]. In other words, the change in urban landscape quality is more influenced by human disturbance than natural factors [11].
Under the influence of long-term human activities, the regional landscape quality will be degraded, and the resilience to external disturbances and ecological stability and other natural regulation functions will be reduced [12]. When spaces with high-level landscape quality are degraded to low-level ones, it is mainly reflected in the reduction in vegetation area and quantity and in the decrease in structural connectivity and stability [13,14]. Currently, the most common and generally accepted indicators are those related to landscape ecology, which also have the advantage of being applicable to any geographical environment [15].
On the one hand, many researchers have studied the objective evaluation of landscapes. Landscape pattern indicators are the most common evaluation metrics. Such indicators have been used by scholars to evaluate the variability of study areas from different scales, times, and regional perspectives. For instance, existing studies have proven keen to explore the spatial and temporal variability of landscape quality across different regions [16], to analyze the scale effects of individual landscape elements (green spaces and water bodies) [17], and to model the influence of land use on landscape quality [18]. In particular, landscape connectivity has been explored in many studies as an important part of evaluating landscape quality. Some studies have quantified landscape connectivity through spatial metrics, such as COHESION metrics or dispersal range [19,20]. Some scholars have applied morphological spatial pattern analysis (MSPA) to analyze the fragmentation of ecological networks and to monitor ecological changes in different periods [21].
On the other hand, research on ecological restoration focuses on the influential mechanisms of biotic and abiotic factors, which means analyzing the transitions of plant and animal communities, the quality of abiotic elements, and the differences between geographical environments during the degradation of landscape quality [22][23][24]. In addition, some studies have used landscape pattern indicators and statistical models to extract ecological information from large-scale areas in order to explore the characteristics and processes of landscape quality degradation. These data are used to restore ecological networks and to identify the priorities of ecosystem services [25][26][27][28][29].
With advancements in remote sensing technology and the accumulation of data, it has become possible to detect and evaluate regional landscape quality changes in the long term [30][31][32]. For instance, in the Guangdong-Hong Kong-Macau Greater Bay Area, Yang et al. integrated vegetative cover, the vegetation health index, the normalized differential build-up and bare soil index, land surface moisture, and land surface temperature to construct a comprehensive index to describe the spatial and temporal variation characteristics of landscape quality [33]. Figueira Branco et al. evaluated the temporal trends in vegetation patterns within the Sooretama Biological Reserve and its surroundings, located in the north of Espírito Santo state, Brazil, and found that the Vegetation Condition Index had a potential for drought occurrence analysis in regions and areas with different vegetation densities [34].
However, most studies have focused on the qualitative analysis of ecological status [35]. Some were biased toward individual indicators, and few conducted comprehensive evaluations from multiple dimensions [36,37]. Fewer researchers still have further analyzed the priority of regional ecological restoration. As a complex, dynamic ecosystem, the regional landscape cannot rely on a single feature to summarize its transformation; various factors, such as area, structure, and function, should also be taken into account [38]. A framework for intuitive and effective ecological restoration, in conjunction with the gradual development of the region, is required, especially in a global first-tier city with a high population density. To fill the gaps in existing studies, this study proposed an integrated approach to evaluate ecological restoration priorities based on regional landscape quality tendencies, taking the Shanghai metropolitan area as the study site. The results should help urban planners/managers to recognize the objective rules of landscape quality in different regions. More importantly, spatial prioritization of landscape quality can be identified to point to regions that need optimal protection. Furthermore, targeted restoration measures can be implemented according to different types of degradation. The specific aims of this study were as follows: (1) to develop a comprehensive index system with wide applicability to quantitatively analyze the regional landscape area, landscape structure, and landscape function; (2) to reveal trends in landscape quality over a long time series and explore the relationship between land cover and landscape quality; and (3) to identify the most degraded areas of landscape quality and classify the types of ecological restoration under different degradation levels.

Study Site
The Shanghai metropolitan area is located on the west coast of the Pacific Ocean, the eastern edge of the Asian continent (120 • 52 to 122 • 12 E, 30 • 40 to 31 • 53 N). The total area of the region is 6340 km 2 , and the landform type is mainly plain, part of the alluvial plain of the Yangtze River Delta. It has a subtropical humid monsoon climate with four distinct seasons, sufficient sunshine, and abundant rainfall. During July and August, the study site experiences abundant rainfall, as well as the densest vegetation foliage. The land cover type appears different, with the outer ring as the dividing line. The region within the outer ring is dominated by residential and commercial land, with a lack of natural vegetation cover [39]. The region beyond the outer ring, with large areas of secondary forest, has an inland lake (Taihu Lake) and numerous water nets ( Figure 1).
As an important economic, transportation, technology, financial, and shipping center in China, Shanghai is one of the world's largest metropolitan areas [40]. Compared to other metropolitan areas in the world, the Shanghai metropolitan area, on the one hand, is located on an impact plain, which is similar to many metropolitan areas, such as Tokyo, New York, and Jakarta [41]. The flat terrain is conducive to promoting urbanization and industrial gathering, while the concentration of population and industrial development also mean that the regional ecology suffers more disturbances.
On the other hand, metropolitan areas share a similar demographic background. For example, Tokyo and New York have population densities of 14,386 and 28,490 inhabitants/km 2 , respectively [42,43]. For the Shanghai metropolitan area, by the end of 2020, the population density was about 23,092 inhabitants/km 2 , which is lower than the 24,137 inhabitants/km 2 in 2010 (Shanghai Statistics Bureau, 2020). This phenomenon implies a population flow to regions outside the urban center, which is also reflected in the development of the rest of the metropolitan area. High-density habitation inevitably leads to the spread of grey infrastructure, meaning that the Shanghai metropolitan area is facing landscape quality degradation and land-use conflicts [44].

Data and Processing
Remote sensing images (RSIs), a digital elevation model (DEM), and Google HD imagery were used as the primary data to analyze landscape quality and land cover characteristics and changes in the study site. RSIs were downloaded from the USGS website (https://glovis.usgs.gov/app, accessed on 10 December 2021). DEM was collected from the Data Center for Resource and Environmental Sciences, Chinese Academy of Sciences (http://www.resdc.cn/, accessed on 10 December 2021), with a resolution of 30 m × 30 m.
Through field surveys and consultation with experts from the Shanghai Greening Bureau and Environmental Protection Bureau, it should be noted that around the year 2005, numerous policies were introduced in the Shanghai metropolitan area, resulting in large-scale landscape changes. The data used in this study are detailed in Table 1. For each year, two remote sensing images needed to be stitched together to extract the extent of the Shanghai metropolitan area (WRS_PATH = 118, WRS_ROW = 038 and WRS_PATH = 118, WRS_ROW = 039).

Data and Processing
Remote sensing images (RSIs), a digital elevation model (DEM), and Google HD imagery were used as the primary data to analyze landscape quality and land cover characteristics and changes in the study site. RSIs were downloaded from the USGS website (https://glovis.usgs.gov/app, accessed on 21 November 2021). DEM was collected from the Data Center for Resource and Environmental Sciences, Chinese Academy of Sciences (http://www.resdc.cn/, accessed on 21 November 2021), with a resolution of 30 m × 30 m.
Through field surveys and consultation with experts from the Shanghai Greening Bureau and Environmental Protection Bureau, it should be noted that around the year 2005, numerous policies were introduced in the Shanghai metropolitan area, resulting in large-scale landscape changes. The data used in this study are detailed in Table 1. For each year, two remote sensing images needed to be stitched together to extract the extent of the Shanghai metropolitan area (WRS_PATH = 118, WRS_ROW = 038 and WRS_PATH = 118, WRS_ROW = 039).
Data from all nine years were used to identify trends in landscape quality. The data of the starting year 2007 and the final year 2020 were used to classify the land cover. The land cover was classified into five categories based on the geography of Shanghai and the land classification in existing urban planning, which contains built-up area, water bodies, cropland, forest, grassland, and barren land. The Region of Interest (ROI) tool was used to create 100 samples for each category of land cover. We ensured that the separation value of ROI was greater than 1.9, which indicated that they were qualified samples. Finally, the maximum likelihood estimation was applied to supervise the classification of land cover in the study site [45]. In addition, the confusion matrix method was chosen for the accuracy validation of the classification results [46]. In total, 100 sample squares were selected separately for each land cover type using Google HD images and field surveys. The accuracy statements were obtained by selecting Confusion Matrix Using Ground Truth ROIs in ENVI 5.3. The results showed that the classification accuracy was above 85%, which satisfied the requirements. Meanwhile, the transition matrix was used to analyze the changes in land cover within the Shanghai metropolitan area, comparing the differences in land cover between 2007 and 2020 [47]. In this study, the proportion of land cover types in 2007 was set as Equation (1), and the proportion of land in category j in 2020 was set as Equation (2).
where C ij (i = j) represents the proportion of land cover shifted from type i to type j between 2007 and 2020, and n is the total number of land cover types. The diagonal entries, C jj , represent the proportion of land cover that shows the persistence of type j [48]. Loss i represents the reduction in land cover in type i, which is equal to the difference between C i+ and C ii . Loss j represents the total increase in land cover in type j from 2007 to 2020, which is equal to the difference between C +j and C jj [49].

Landscape Quality Index System
With landscape ecology as the theoretical basis of research, this study considered the spatial and temporal changes in the landscape in the regional ecosystem. By analyzing existing studies, indicators characterizing three dimensions of landscape area, landscape structure, and landscape function were selected to construct an index system of landscape quality in the Shanghai metropolitan area.

1.
Landscape area indicator (LAI) There are differences in the ability of various landscape types to maintain ecosystem structure and function and to resist external disturbances [30]. This is directly affected by the area occupied by the land cover type. The landscape area is a visual representation of the change in the regional ecosystem and an important indicator of the landscape quality. The enhanced vegetation index (EVI) is similar to the normalized difference vegetation index (NDVI). However, EVI corrects for some atmospheric conditions and canopy background noise and is more sensitive in regions with dense vegetation, compensating for the inability of NDVI to correct for changes in solar incidence angle [50]. Therefore, this study chose EVI to quantify the degree of variation in regional LAI change [51]. The equation for the LAI is as follows:

Landscape structure indicator (LSI)
To quantify the shape and configuration of the regional landscape, the LSI was constructed by fusing the fractal dimension index (FRAC) and the contiguity index (CONTAG). It is important to note that for LSI, the target was all land cover types, excluding barren land and built-up area. On the one hand, FRAC is often used to quantitatively describe the size of the core area and the clutter of the boundaries of a patch [52]. Compared to other landscape metrics that describe spatial patterns, such as the shape index, the mean perimeter-area ratio (PARA_MN), and MSPA, FRAC compensates for some shortcomings through algorithms. First, it overcomes the main limitation of PARA_MN in measuring shape complexity [53]. Second, FRAC solves the dimensional dependence problem, which is caused by the calculation of the perimeter area ratio as a shape complexity metric [54]. Thus, FRAC was chosen to describe the shape of landscape patches. The calculation formula is as follows: where P ij = perimeter (m) of the patch ij, and A ij = area (m 2 ) of the patch ij. The range of FRAC is 1 to 2. The closer the FRAC is to 2, the more fragmented the patch's shape is (i.e., the more it is disturbed by human activities). The units: none. The contiguity index (CONTAG) was used to assess the spatial connectivity and continuity within the raster cells, which can reflect the degree of fragmentation and the interference of human activities [55]. For example, landscapes with a few large, continuous patches (containing most internal units with similar proximity) may have higher values of CONTAG. Thus, compared to the interspersion and juxtaposition index (based on patch adjacencies), CONTAG can visualize the dispersion and interspersion of landscape patches by their numerical size [52]. The calculation formula is as follows: where c ijr = contiguity value for pixel r in patch ij, v = sum of the values in a 3-by-3 cell template, and a ij = area of patch ij in terms of the number of cells. The range of CONTIG is 0-1. The closer the value of CONTIG is to 0, the more it represents an increase in the continuity or connectivity of the patch. In addition, it can also represent a decrease in patch connectivity. The units: none. The LSI was composed of FRAC and CONTAG together. It should be noted that FRAC and CONTIG are both dimensionless indicators, but their values have different trends. When the value of FRAC and CONTIG is smaller, the shape of the patch is more regular, but the connectivity is worse. Thus, while normalizing the data, the value of FRAC needs to be inverted, which means that when the value of FRAC is smaller, the patch's shape is more fragmented. The first component, which is LSI, is then obtained by principal component analysis (PCA). The formula is as follows: The range of LSI is 0-1. The closer the value of LSI is to 0, the worse the connectivity of the patch and the more fragmented the shape. Moreover, it represents an increase in the morphological regularity and connectivity of the patch. The units: none.

Landscape function index (LFI)
Many studies have argued that heat, greenness, moisture, and dryness are closely related to the ecological functions corresponding to the ecosystem [26]. Any perturbation or change in these ecological factors will eventually act on the whole ecosystem [56,57]. The traditional evaluation method only incorporates a single type of remote sensing data into the index system. However, some indicators are difficult to obtain directly by remote sensing, such as surface air temperature, which can often have significant bias by inversion only [58]. Using a combination of multiple layers of remote sensing as the input data set for LFI is a good way to reduce error and attenuate the effect of a single type of dataset.
Land surface temperature (LST), normalized difference vegetation index (NDVI), WET, and normalized differential build-up and bare soil index (NDBSI), corresponding to heat, greenness, moisture, and dryness, respectively, have been widely used to evaluate landscape quality [59]. Therefore, this study obtained four indicators by band calculation [60]. They were then normalized from 0 to 1. Finally, the first component was taken to represent LFI by PCA [61]. The formula is as follows:

Tendency of Landscape Quality and Its Relationship to Land Cover Conversion
Tests for detecting trends in time-series data can be divided into parametric and non-parametric methods [62]. Non-parametric methods require only that the data be independent and insensitive to measurement errors and outlier data and to have been applied to the detection of time-series data in fields such as hydrology, meteorology, and vegetation [11,63]. The Mann-Kendall (MK) test, a widely used non-parametric statistical test, has the advantage that it does not require the measured values to follow a normal distribution and is not affected by missing values and outliers [64,65]. Therefore, this study used Mann-Kendall and Sen's slope estimation to identify and test the trend of landscape quality in the Shanghai metropolitan area.
The statistic S of the MK test is calculated as follows [66]: where n is the number of data points, and x i and x j are the data values in time series i and j (j > i), respectively. The detailed calculation of sgn x j − x i and the variance within Equation (10) can be found in the research of Milan Gocic et al. [62]. It should be noted that 8 years in the 2007-2020 time-series data were selected in this study, which means that n ≥ 8. Mann et al. explained that when n ≥ 8, S is almost normally distributed with consequent mean and variance [67,68].
The standard normal test statistic Z is expressed as A positive (negative) value of Z indicates that the data tend to increase (decrease) with time. Testing trends were carried out at the specific α significance level. When Z > Z 1−α/2 , the null hypothesis is rejected, and a significant trend exists. Significance levels α = 0.01 and α = 0.05 were used in this study. The null hypothesis of no trend was rejected for |Z| > 1.96 at a 5% significance level. Sen's slope estimator (known as the Theil-Sen median method) is a robust nonparametric statistical trend calculation method [69]. The calculation formula is as follows: where x j and x k are observations for years j and k (j > k), respectively; N is the number of time periods, N = n(n − 1)/2; and Q med is the median of Q i , whose value indicates the steepness of the data trend [62]. We applied the reclassification toolsets −1 for degradation, 0 for stability, and 1 for improvement. Based on the trend analysis, the Z values were classified into 2 categories (significant and insignificant), and the Sen's slope estimator results were classified into 3 categories (degradation, stable, and improvement) and multiplied by the raster calculator. The final trend was divided into 5 levels: significant degradation (−2), slight degradation (−1), stable (0), slight improvement (1), and significant improvement (2). Each of these 5 levels was assigned to each raster as a quantitative value of landscape quality. In other words, each type of land cover that changed had a corresponding trend level. Meanwhile, each raster had a corresponding land cover type in 2020 (transformed from the different land cover types in 2007). Furthermore, the spaces of different degradation levels and the land cover types to which they belong could establish a one-to-one relationship. Based on the land cover transfer matrix, the percentage of different degradation grading in the transformation type could be calculated.
All steps were calculated with the software R 4.1.0 and graphically represented in ArcGIS Pro 2.0. The specific steps and code data were referenced in the package "trend" [70]. The data package "raster" was used for pixel-by-pixel analysis [71]. Based on the above, the "terra" package was added to the original code to enable parallel computation and increase the computation speed [72]. Details of all packages and R scripts can be found on the related web page (https://cran.r-project.org/, accessed on 21 November 2021).

Ecological Restoration Priorities
The identification of ecological restoration priorities by combining all indicators is a difficult task. There may be different dimensions and ranges of values among indicators. Therefore, first, all metrics were selected and constructed in a dimensionless and standardized operation to ensure that their values all ranged from 0 to 1. According to the expert survey (Supplementary), the indicators of the three dimensions were assigned equal weights. The spatial distribution of ecological degradation was calculated by the spatial overlay. Second, it was further divided into five grades using the natural break method to reflect the priority level of ecological restoration, implying that the more severe the degradation of landscape quality, the higher the corresponding grading of ecological restoration.
Furthermore, the following steps were followed to achieve the purpose of targeted ecological restoration in regions with low landscape quality, with the identification of regions in need of ecological restoration and those whose values of landscape quality indicators (LAI, LSI, and LFI) of −2 and −1 were selected for multi-layer overlay. When an indicator or multiple indicators were −2, it signified that it was the dominant type of landscape quality degradation. Thus, the types of degradation dominated by different indicators were able to be classified into 8 types: Type A: degradation led by LFI; Type B: degradation led by LAI; Type C: degradation led by LSI; Type D: degradation co-led by LAI and LSI; Type E: degradation co-led by LAI and LFI; Type F: degradation co-led by LSI and LFI; Type G: no indicator as degradation led; and Type H: degradation co-led by LSI, LAI and LFI.

Changes in Land Cover
The land cover types in the Shanghai metropolitan area in 2007 and 2020 are shown in Figure 2. In 2007, a large proportion of forests were distributed in the periphery of the Shanghai metropolitan area. Small-and medium-sized ponds and wetlands formed a nested relationship with cropland and grassland. Affected by rapid urbanization, it was expected to see a significant outward expansion of the urban center by 2020. Compared with 2007, excluding built-up area (68.2% of the land area), the main land cover in 2020 was cropland and forest, at 20.9% and 9.8%, respectively. From the urban center to the suburbs, the vegetation coverage showed a gradual trend of decreasing, with the existing forest concentrated in large open green spaces and the Lin-Gang Special Area [73]. Meanwhile, former ponds and rivers were encroached upon by built-up area, mainly in the western and northern tracts.
Land 2022, 10, x FOR PEER REVIEW 9 of 22 and LFI; Type G: no indicator as degradation led; and Type H: degradation co-led by LSI, LAI and LFI.

Changes in Land Cover
The land cover types in the Shanghai metropolitan area in 2007 and 2020 are shown in Figure 2. In 2007, a large proportion of forests were distributed in the periphery of the Shanghai metropolitan area. Small-and medium-sized ponds and wetlands formed a nested relationship with cropland and grassland. Affected by rapid urbanization, it was expected to see a significant outward expansion of the urban center by 2020. Compared with 2007, excluding built-up area (68.2% of the land area), the main land cover in 2020 was cropland and forest, at 20.9% and 9.8%, respectively. From the urban center to the suburbs, the vegetation coverage showed a gradual trend of decreasing, with the existing forest concentrated in large open green spaces and the Lin-Gang Special Area [73]. Meanwhile, former ponds and rivers were encroached upon by built-up area, mainly in the western and northern tracts. The total area of land cover transition during this period was about 3082.18 km 2 . The land cover with the largest net increase in area was built-up area, mainly from the transition of cropland, forests, and water bodies. The land cover type with the first net reduction was cropland, which changed to forests and built-up area, at 14.14% and 46.30%, respectively (Table 2). Although 748.96 km 2 of cropland was turned into built-up area, the area of built-up area and forests that transformed into cropland was also sizeable (712.91 km 2 ). The second largest reduction in land cover was water bodies, which mainly changed to built-up area, cropland, and forests, accounting for 53.16%, 25.40%, and 19.18% of the total converted area of water bodies, respectively. The grassland and barren land took up a smaller proportion of the area. The total area of land cover transition during this period was about 3082.18 km 2 . The land cover with the largest net increase in area was built-up area, mainly from the transition of cropland, forests, and water bodies. The land cover type with the first net reduction was cropland, which changed to forests and built-up area, at 14.14% and 46.30%, respectively ( Table 2). Although 748.96 km 2 of cropland was turned into built-up area, the area of built-up area and forests that transformed into cropland was also sizeable (712.91 km 2 ). The second largest reduction in land cover was water bodies, which mainly changed to built-up area, cropland, and forests, accounting for 53.16%, 25.40%, and 19.18% of the total converted area of water bodies, respectively. The grassland and barren land took up a smaller proportion of the area. The rows represent which types of land cover this category will be converted to from 2007 to 2020. The numbers on columns represent the area of this land cover type converted from other types in 2020. 1 represents the percentage area of a land cover type converted to the remaining types compared to 2007. 2 represents the percentage area of land cover type that has not been converted compared to the area of that type in 2020.

Grading of Landscape Quality Trends
As expected, the regions with low LAI in 2007 were concentrated in the port terminals and river banks. The regions with low LAI in 2020 gradually expanded to the periphery of the region and suburban villages. Meanwhile, the median of LAI in 2007 was larger than that in 2020 (Med LAI.2007 is 0.707 and Med LAI.2020 is 0.395), indicating that LAI decreased during this period and that the distribution of areas with lower LAI grew more concentrated. As a result of the multi-year trend analysis (Figure 3c1), there was 4684.23 km 2 of patches with LAI degradation from 2007 to 2020, or about 57.37% of the study site area. Significant degradation regions occupied 5.16%. However, in terms of total area, the difference between the improved and degraded areas of LAI was not significant. Regardless of whether in the urban center or the suburbs, the area (LAI improvement) was much larger than the degraded area. Meanwhile, compared to the suburbs, the proportion of patches, showing LAI improvement, was relatively higher in the urban center ( Figure 3d1). Figure 3a2,b2 show that the landscape structure of the Shanghai metropolitan area underwent degradation between 2007 and 2020. The LSI in 2007 was, on the whole, less disturbed by human activities. Only in Jinshan District and the Lin-Gang Special Area, where factories were concentrated, did the landscape structure show fragmentation. After 13 years of urbanization, the LSI of the study site declined in 2020, which was especially obvious in the coastal area. The area of degradation reached 1729.67 km 2 , which is 68.2% of the total trend area. The rate of degradation was slower in the urban center compared to the suburbs (Figure 3d2). Moreover, the concentration of degraded LSI shifted from the southern part of Jinshan District to the periphery of this region (including Chongming Island), where the main land cover types were forests and cropland.
Comparing Figure 3a3,b3, the Shanghai metropolitan area also showed a significant degradation in LFI. Regions with low LFI in 2007 were restricted to the urban center and suburban villages, yet expanded to the entire metropolitan area by 2020. The number of regions exhibiting extremely low LFI increased. Only parts of the outer coastal areas, such as Chongming Island, showed green. In addition, we found that the trends in LFI were differentially distributed in the study site (Figure 3c3). On the one hand, regions with degradation covered more space (4560.44 km 2 , 67.81%) compared to those with improved and stable LFI. On the other hand, the percentage of area with stable LFI in the urban center was 13.8%, larger than that in the suburbs. In the suburbs, there were more areas where LFI improved. The regions with a significant improvement in LFI were concentrated in the northern part of Chongming Island-mainly cropland and forests. d 2022, 10, x FOR PEER REVIEW 11 LFI improved. The regions with a significant improvement in LFI were concentrated the northern part of Chongming Island-mainly cropland and forests.

Relationship between Land Cover and Landscape Quality Indicators
Based on the change in land cover and trends of landscape quality from 2007 to 20 compared to other land covers, the areas of water bodies, barren land, and grassland w rare, and their transition rates were low. Therefore, this study focused on the relations between built-up area, forest, and cropland with landscape quality indicators. A twotwo pairing of land cover and the type of its transfer was made through the grading landscape indicators.
As can be seen from Table 3, the landscape quality of the same area changed to ferent degrees from 2007 to 2020. The landscape quality of land cover was mainly conc trated in the −1 to 1 grade range, which was similar to the trend in 2007 and 2020 for same land cover type. Firstly, for LAI, the forests had the highest stability over long-te trends, whether from built-up area, cropland, or unconverted regions of the forest its However, water bodies had the second highest LAI variability, showing a stable slightly degraded trend (degradation grades of 0 and −1 were in the range of 10-20 Secondly, for LSI, the forests and water bodies showed a larger degree of degradation maximum percentage in −2 grade). It is noteworthy that the LSI of forests had an incr ing trend on the +1 grade. Finally, for LFI, built-up area showed a significant trend slight degradation (−1 grade, accounting for 20-70%). The process of mutual transit between forest and cropland also showed a slight degradation trend in the region.

Relationship between Land Cover and Landscape Quality Indicators
Based on the change in land cover and trends of landscape quality from 2007 to 2020, compared to other land covers, the areas of water bodies, barren land, and grassland were rare, and their transition rates were low. Therefore, this study focused on the relationship between built-up area, forest, and cropland with landscape quality indicators. A two-bytwo pairing of land cover and the type of its transfer was made through the grading of landscape indicators.
As can be seen from Table 3, the landscape quality of the same area changed to different degrees from 2007 to 2020. The landscape quality of land cover was mainly concentrated in the −1 to 1 grade range, which was similar to the trend in 2007 and 2020 for the same land cover type. Firstly, for LAI, the forests had the highest stability over long-term trends, whether from built-up area, cropland, or unconverted regions of the forest itself. However, water bodies had the second highest LAI variability, showing a stable or slightly degraded trend (degradation grades of 0 and −1 were in the range of 10-20%). Secondly, for LSI, the forests and water bodies showed a larger degree of degradation (the maximum percentage in −2 grade). It is noteworthy that the LSI of forests had an increasing trend on the +1 grade. Finally, for LFI, built-up area showed a significant trend of slight degradation (−1 grade, accounting for 20-70%). The process of mutual transition between forest and cropland also showed a slight degradation trend in the region.

Prioritization and Classification of Ecological Restoration
As shown in Figure 4a, the primary ecological restoration zone was 1160.62 km 2 , the secondary ecological restoration zone was 3324.20 km 2 , the ecological stable zone was 365.49 km 2 , the weak ecological improvement zone was 1799.37 km 2 , and the strong ecological improvement zone was 75.80 km 2 . The priority restoration zone showed a pattern wherein the suburbs had more area than the urban center; see, for example, the external region of the outer ring, the northern part of Chongming Island, and the village agglomeration regions in the suburbs. The sub-restoration zone was widely distributed, spreading outward with the urban center as the core. When land cover areas with grades −2 and −1 were considered to have undergone degradation of landscape quality, the land cover transition patterns that showed the degradation were as follows (transition from the former type to the latter type): (1) for LAI, built-up area-water bodies, forest-water bodies, and cropland-water bodies; (2) for LSI, built-up area-forest, forest-forest, and cropland-forest; (3) for LFI, built-up area-built-up area, forest-cropland/built-up area, and cropland-cropland/built-up area.

Prioritization and Classification of Ecological Restoration
As shown in Figure 4a, the primary ecological restoration zone was 1160.62 km 2 , the secondary ecological restoration zone was 3324.20 km 2 , the ecological stable zone was 365.49 km 2 , the weak ecological improvement zone was 1799.37 km 2 , and the strong ecological improvement zone was 75.80 km 2 . The priority restoration zone showed a pattern wherein the suburbs had more area than the urban center; see, for example, the external region of the outer ring, the northern part of Chongming Island, and the village agglomeration regions in the suburbs. The sub-restoration zone was widely distributed, spreading outward with the urban center as the core.

Prioritization and Classification of Ecological Restoration
As shown in Figure 4a, the primary ecological restoration zone was 1160.62 km 2 , the secondary ecological restoration zone was 3324.20 km 2 , the ecological stable zone was 365.49 km 2 , the weak ecological improvement zone was 1799.37 km 2 , and the strong ecological improvement zone was 75.80 km 2 . The priority restoration zone showed a pattern wherein the suburbs had more area than the urban center; see, for example, the external region of the outer ring, the northern part of Chongming Island, and the village agglomeration regions in the suburbs. The sub-restoration zone was widely distributed, spreading outward with the urban center as the core. Furthermore, the regions with grades −2 and −1 were selected, in which the ecological restoration zones dominated by different landscape quality indicators could be visually identified. Figure 4b shows that Types B, C, D, and G have larger areas compared to the remaining types. Type B was primarily a green space near the port along the Puyang River, including residential and factory districts. Type C was distributed around linear roads and water bodies, with a discrete spatial distribution. Type D was mostly found near the built-up area and water bodies, dominated by public green space and woodland. Type G was the type with the largest area percentage and was distributed around cropland in the suburbs, showing a belt-shaped distribution and fragmented space (Table  4). Furthermore, the regions with grades −2 and −1 were selected, in which the ecological restoration zones dominated by different landscape quality indicators could be visually identified. Figure 4b shows that Types B, C, D, and G have larger areas compared to the remaining types. Type B was primarily a green space near the port along the Puyang River, including residential and factory districts. Type C was distributed around linear roads and water bodies, with a discrete spatial distribution. Type D was mostly found near the built-up area and water bodies, dominated by public green space and woodland. Type G was the type with the largest area percentage and was distributed around cropland in the suburbs, showing a belt-shaped distribution and fragmented space (Table 4).

0.048
Type B The grading of LAI is −2, when the grading of LSI and LFI is −1.

97.357
Type F The grading of LSI and LFI is −2, when the grading of LAI is −1.

2.644
Type E The grading of LAI and LFI is −2, when the grading of LSI is −1.

0.048
Type B The grading of LAI is −2, when the grading of LSI and LFI is −1.

97.357
Type F The grading of LSI and LFI is −2, when the grading of LAI is −1.

0.048
Type B The grading of LAI is −2, when the grading of LSI and LFI is −1.

97.357
Type F The grading of LSI and LFI is −2, when the grading of LAI is −1.

2.644
Type C The grading of LSI is −2, when the grading of LAI and LFI is −1.

502.990
Type G The grading of LAI, LSI and LFI is −1.

962.478
Type D The grading of LAI and LSI is −2, when the grading of LFI is −1.

97.403
Type H The grading of LAI, LSI and LFI is −2. 0.219

Trends of Landscape Quality in Metropolitan Areas
With the urbanization in the last decade, the difference in landscape quality between the urban center and the outer-ring regions is particularly significant. The urban center has high population densities and large impervious surfaces, which means fragmented green space, UHI effects, and dry surface environments [74,75]. This study found that the landscape quality of the Shanghai metropolitan area has shown a trend of degradation over long-term urbanization. Regions of low landscape quality in 2007 existed mainly in the urban center. In 2020, significant degradation and slight degradation spread across the boundaries of the urban center (the outer-ring road) and toward the distant suburbs. It is notable that the variation in low values of landscape quality indicators also showed a gradual shift from the urban center to the regional boundaries (Figure 3c1-c3), consistent with the stage of suburbanization development in the van den Berg model [76].
The reason for such a result is that, on the one hand, urban development is inevitably accompanied by the adaptation of the original environment. As the urban center is the core of the initial development, the land-use efficiency is usually high [77]. In spatially fixed areas, such as green infrastructure and built-up area, the trend of landscape quality is milder compared to that of the suburbs. On the contrary, the suburbs are vast and sparsely populated, where the inefficient use of land and waste of resources are common [78]. The land cover in the area studied was mostly secondary forest and cropland, which are sensitive to human activities (Figure 2). On the other hand, in China, the governmentled land market is gradually improving, and the transformation of land properties depends on the implementation of urban planning. Notably, after the year 2000, the government gained control over the suburban districts and counties through administrative annexation [79], which led to the prevalence of land commoditization [80]. Suburbanization

Type F
The grading of LSI and LFI is −2, when the grading of LAI is −1.

0.048
Type B The grading of LAI is −2, when the grading of LSI and LFI is −1.

97.357
Type F The grading of LSI and LFI is −2, when the grading of LAI is −1.

2.644
Type C The grading of LSI is −2, when the grading of LAI and LFI is −1.

502.990
Type G The grading of LAI, LSI and LFI is −1.

962.478
Type D The grading of LAI and LSI is −2, when the grading of LFI is −1.

97.403
Type H The grading of LAI, LSI and LFI is −2. 0.219

Trends of Landscape Quality in Metropolitan Areas
With the urbanization in the last decade, the difference in landscape quality between the urban center and the outer-ring regions is particularly significant. The urban center has high population densities and large impervious surfaces, which means fragmented green space, UHI effects, and dry surface environments [74,75]. This study found that the landscape quality of the Shanghai metropolitan area has shown a trend of degradation over long-term urbanization. Regions of low landscape quality in 2007 existed mainly in the urban center. In 2020, significant degradation and slight degradation spread across the boundaries of the urban center (the outer-ring road) and toward the distant suburbs. It is notable that the variation in low values of landscape quality indicators also showed a gradual shift from the urban center to the regional boundaries (Figure 3c1-c3), consistent with the stage of suburbanization development in the van den Berg model [76].
The reason for such a result is that, on the one hand, urban development is inevitably accompanied by the adaptation of the original environment. As the urban center is the core of the initial development, the land-use efficiency is usually high [77]. In spatially fixed areas, such as green infrastructure and built-up area, the trend of landscape quality is milder compared to that of the suburbs. On the contrary, the suburbs are vast and sparsely populated, where the inefficient use of land and waste of resources are common [78]. The land cover in the area studied was mostly secondary forest and cropland, which are sensitive to human activities (Figure 2). On the other hand, in China, the governmentled land market is gradually improving, and the transformation of land properties depends on the implementation of urban planning. Notably, after the year 2000, the government gained control over the suburban districts and counties through administrative annexation [79], which led to the prevalence of land commoditization [80]. Suburbanization

Type C
The grading of LSI is −2, when the grading of LAI and LFI is −1. The grading of LFI is −2, when the grading of LAI and LSI is −1.

Type E
The grading of LAI and LFI is −2, when the grading of LSI is −1.

0.048
Type B The grading of LAI is −2, when the grading of LSI and LFI is −1.

97.357
Type F The grading of LSI and LFI is −2, when the grading of LAI is −1.

2.644
Type C The grading of LSI is −2, when the grading of LAI and LFI is −1.

502.990
Type G The grading of LAI, LSI and LFI is −1.

962.478
Type D The grading of LAI and LSI is −2, when the grading of LFI is −1.

97.403
Type H The grading of LAI, LSI and LFI is −2. 0.219

Trends of Landscape Quality in Metropolitan Areas
With the urbanization in the last decade, the difference in landscape quality between the urban center and the outer-ring regions is particularly significant. The urban center has high population densities and large impervious surfaces, which means fragmented green space, UHI effects, and dry surface environments [74,75]. This study found that the landscape quality of the Shanghai metropolitan area has shown a trend of degradation over long-term urbanization. Regions of low landscape quality in 2007 existed mainly in the urban center. In 2020, significant degradation and slight degradation spread across the boundaries of the urban center (the outer-ring road) and toward the distant suburbs. It is notable that the variation in low values of landscape quality indicators also showed a gradual shift from the urban center to the regional boundaries (Figure 3c1-c3), consistent with the stage of suburbanization development in the van den Berg model [76].
The reason for such a result is that, on the one hand, urban development is inevitably accompanied by the adaptation of the original environment. As the urban center is the core of the initial development, the land-use efficiency is usually high [77]. In spatially fixed areas, such as green infrastructure and built-up area, the trend of landscape quality is milder compared to that of the suburbs. On the contrary, the suburbs are vast and sparsely populated, where the inefficient use of land and waste of resources are common [78]. The land cover in the area studied was mostly secondary forest and cropland, which are sensitive to human activities (Figure 2). On the other hand, in China, the governmentled land market is gradually improving, and the transformation of land properties depends on the implementation of urban planning. Notably, after the year 2000, the government gained control over the suburban districts and counties through administrative annexation [79], which led to the prevalence of land commoditization [80]. Suburbanization

Type G
The grading of LAI, LSI and LFI is −1. The grading of LFI is −2, when the grading of LAI and LSI is −1.

Type E
The grading of LAI and LFI is −2, when the grading of LSI is −1.

0.048
Type B The grading of LAI is −2, when the grading of LSI and LFI is −1.

97.357
Type F The grading of LSI and LFI is −2, when the grading of LAI is −1.

2.644
Type C The grading of LSI is −2, when the grading of LAI and LFI is −1.

502.990
Type G The grading of LAI, LSI and LFI is −1.

962.478
Type D The grading of LAI and LSI is −2, when the grading of LFI is −1.

97.403
Type H The grading of LAI, LSI and LFI is −2. 0.219

Trends of Landscape Quality in Metropolitan Areas
With the urbanization in the last decade, the difference in landscape quality between the urban center and the outer-ring regions is particularly significant. The urban center has high population densities and large impervious surfaces, which means fragmented green space, UHI effects, and dry surface environments [74,75]. This study found that the landscape quality of the Shanghai metropolitan area has shown a trend of degradation over long-term urbanization. Regions of low landscape quality in 2007 existed mainly in the urban center. In 2020, significant degradation and slight degradation spread across the boundaries of the urban center (the outer-ring road) and toward the distant suburbs. It is notable that the variation in low values of landscape quality indicators also showed a gradual shift from the urban center to the regional boundaries (Figure 3c1-c3), consistent with the stage of suburbanization development in the van den Berg model [76].
The reason for such a result is that, on the one hand, urban development is inevitably accompanied by the adaptation of the original environment. As the urban center is the core of the initial development, the land-use efficiency is usually high [77]. In spatially fixed areas, such as green infrastructure and built-up area, the trend of landscape quality is milder compared to that of the suburbs. On the contrary, the suburbs are vast and sparsely populated, where the inefficient use of land and waste of resources are common [78]. The land cover in the area studied was mostly secondary forest and cropland, which are sensitive to human activities (Figure 2). On the other hand, in China, the governmentled land market is gradually improving, and the transformation of land properties depends on the implementation of urban planning. Notably, after the year 2000, the government gained control over the suburban districts and counties through administrative annexation [79], which led to the prevalence of land commoditization [80]. Suburbanization

Type D
The grading of LAI and LSI is −2, when the grading of LFI is −1. The grading of LFI is −2, when the grading of LAI and LSI is −1.

Type E
The grading of LAI and LFI is −2, when the grading of LSI is −1.

0.048
Type B The grading of LAI is −2, when the grading of LSI and LFI is −1.

97.357
Type F The grading of LSI and LFI is −2, when the grading of LAI is −1.

2.644
Type C The grading of LSI is −2, when the grading of LAI and LFI is −1.

502.990
Type G The grading of LAI, LSI and LFI is −1.

962.478
Type D The grading of LAI and LSI is −2, when the grading of LFI is −1.

97.403
Type H The grading of LAI, LSI and LFI is −2. 0.219

Trends of Landscape Quality in Metropolitan Areas
With the urbanization in the last decade, the difference in landscape quality between the urban center and the outer-ring regions is particularly significant. The urban center has high population densities and large impervious surfaces, which means fragmented green space, UHI effects, and dry surface environments [74,75]. This study found that the landscape quality of the Shanghai metropolitan area has shown a trend of degradation over long-term urbanization. Regions of low landscape quality in 2007 existed mainly in the urban center. In 2020, significant degradation and slight degradation spread across the boundaries of the urban center (the outer-ring road) and toward the distant suburbs. It is notable that the variation in low values of landscape quality indicators also showed a gradual shift from the urban center to the regional boundaries (Figure 3c1-c3), consistent with the stage of suburbanization development in the van den Berg model [76].
The reason for such a result is that, on the one hand, urban development is inevitably accompanied by the adaptation of the original environment. As the urban center is the core of the initial development, the land-use efficiency is usually high [77]. In spatially fixed areas, such as green infrastructure and built-up area, the trend of landscape quality is milder compared to that of the suburbs. On the contrary, the suburbs are vast and sparsely populated, where the inefficient use of land and waste of resources are common [78]. The land cover in the area studied was mostly secondary forest and cropland, which are sensitive to human activities ( Figure 2). On the other hand, in China, the governmentled land market is gradually improving, and the transformation of land properties depends on the implementation of urban planning. Notably, after the year 2000, the government gained control over the suburban districts and counties through administrative annexation [79], which led to the prevalence of land commoditization [80]. Suburbanization

0.048
Type B The grading of LAI is −2, when the grading of LSI and LFI is −1.

97.357
Type F The grading of LSI and LFI is −2, when the grading of LAI is −1.

2.644
Type C The grading of LSI is −2, when the grading of LAI and LFI is −1.

502.990
Type G The grading of LAI, LSI and LFI is −1.

962.478
Type D The grading of LAI and LSI is −2, when the grading of LFI is −1.

97.403
Type H The grading of LAI, LSI and LFI is −2. 0.219

Trends of Landscape Quality in Metropolitan Areas
With the urbanization in the last decade, the difference in landscape quality between the urban center and the outer-ring regions is particularly significant. The urban center has high population densities and large impervious surfaces, which means fragmented green space, UHI effects, and dry surface environments [74,75]. This study found that the landscape quality of the Shanghai metropolitan area has shown a trend of degradation over long-term urbanization. Regions of low landscape quality in 2007 existed mainly in the urban center. In 2020, significant degradation and slight degradation spread across the boundaries of the urban center (the outer-ring road) and toward the distant suburbs. It is notable that the variation in low values of landscape quality indicators also showed a gradual shift from the urban center to the regional boundaries (Figure 3c1-c3), consistent with the stage of suburbanization development in the van den Berg model [76].
The reason for such a result is that, on the one hand, urban development is inevitably accompanied by the adaptation of the original environment. As the urban center is the core of the initial development, the land-use efficiency is usually high [77]. In spatially fixed areas, such as green infrastructure and built-up area, the trend of landscape quality is milder compared to that of the suburbs. On the contrary, the suburbs are vast and sparsely populated, where the inefficient use of land and waste of resources are common [78]. The land cover in the area studied was mostly secondary forest and cropland, which are sensitive to human activities ( Figure 2). On the other hand, in China, the governmentled land market is gradually improving, and the transformation of land properties depends on the implementation of urban planning. Notably, after the year 2000, the government gained control over the suburban districts and counties through administrative annexation [79], which led to the prevalence of land commoditization [80]. Suburbanization

Trends of Landscape Quality in Metropolitan Areas
With the urbanization in the last decade, the difference in landscape quality between the urban center and the outer-ring regions is particularly significant. The urban center has high population densities and large impervious surfaces, which means fragmented green space, UHI effects, and dry surface environments [74,75]. This study found that the landscape quality of the Shanghai metropolitan area has shown a trend of degradation over long-term urbanization. Regions of low landscape quality in 2007 existed mainly in the urban center. In 2020, significant degradation and slight degradation spread across the boundaries of the urban center (the outer-ring road) and toward the distant suburbs. It is notable that the variation in low values of landscape quality indicators also showed a gradual shift from the urban center to the regional boundaries (Figure 3c1-c3), consistent with the stage of suburbanization development in the van den Berg model [76].
The reason for such a result is that, on the one hand, urban development is inevitably accompanied by the adaptation of the original environment. As the urban center is the core of the initial development, the land-use efficiency is usually high [77]. In spatially fixed areas, such as green infrastructure and built-up area, the trend of landscape quality is milder compared to that of the suburbs. On the contrary, the suburbs are vast and sparsely populated, where the inefficient use of land and waste of resources are common [78]. The land cover in the area studied was mostly secondary forest and cropland, which are sensitive to human activities ( Figure 2). On the other hand, in China, the government-led land market is gradually improving, and the transformation of land properties depends on the implementation of urban planning. Notably, after the year 2000, the government gained control over the suburban districts and counties through administrative annexation [79], which led to the prevalence of land commoditization [80]. Suburbanization serves to satisfy the need to slow down capital flows and is driven by the potential devaluation of capital [81,82]. Consequently, suburbanization is occurring in many metropolitan areas around the world [83], leading to a gradual proliferation of human interventions in the environment [84]. These may explain the phenomena of the emergence of landscape quality in this study.

Drivers of Land Cover Transition on Landscape Quality
Landscape quality is inextricably linked to land cover type [85]. With the expansion of built-up area, the degradation of landscape quality spreads to the edge of metropolitan areas. From 2007 to 2020, built-up area encroached on about 1255 km 2 of forests and cropland. However, due to China's long-standing basic cropland protection policy and the strictly enforced green space target [86,87], a large amount of cropland and forests remain in the suburbs of the Shanghai metropolitan area (Figure 2).
This study found that the trend of LAI was stable, whether the land was built-up area, cropland converted to forests, or forests that remained unchanged, indicating that the transition of other types to forests did not lead to significant area loss in the Shanghai metropolitan area. This may explain the trend of a slightly improved LSI for the transition of built-up area and cropland to forests. However, the transition of built-up area-water bodies, forest-water bodies, and cropland-water bodies in LAI showed a degradation trend of regional LAI. On the one hand, this may require more area for water bodies in land transition compared to other land cover types, such as flood protection levees for rivers and the hardening of water bodies. On the other hand, remote sensing was taken in September, when the rainy season is over in Shanghai and the relative water level is low, which may also have led to degradation in the final results.
The results of LSI revealed that the percentage of significantly degraded area of cropland-forest was 55.31%. This means that the landscape structure can be severely affected by the process of reforestation. The ecological value of forests is higher than that of cropland, but in a region with a high degree of human activities, such as a metropolitan area, forests converted from cropland are the lands with lowest productivity [88]. These lands are limited in size and tend to form isolated islands of forest, which, in turn, reduce land connectivity [89]. However, 11.55% of the areas slightly improved during the transformation process, indicating that the possibility of LSI improvement still exists. This is consistent with the finding that a small amount of forest growth is beneficial for landscape connectivity, as found by Mrinmay Mandal et al. [90].
However, cropland is an important vehicle for maintaining regional biological flows and material exchange [91], which indicates that, when cropland is transformed into builtup area, it is often accompanied by a degradation of landscape quality. Therefore, maintaining the natural state of cropland and solving the problem of soil degradation/residue contamination can effectively guarantee the landscape function of the metropolitan area [92].
The results of the LFI showed significant degradation of built-up area-built-up area, forest-cropland/built-up area, and cropland-cropland/built-up area. The essential difference between forest, cropland, and built-up area is vegetation cover and impermeable surface [26]. Thus, the transition from forests and cropland to built-up area inevitably causes a reduction in LFI. It is noteworthy that a significant degradation of LFI also occurs when other land cover types are converted to cropland. Cropland (as a special type of land cover), different seasons, and crop varieties can lead to LFI changes [93][94][95]. Crops planted on cropland within the Shanghai metropolitan area were dominated by green leafy vegetables, rice, and fruit trees, implying that changes in crop species may have caused the degradation of LFI.

The Trade-Off between Ecological Restoration and Urban Development
This study found that the development of the Shanghai metropolitan area has mainly come at the expense of invading cropland and forests in the past decades. This is similar to the development pattern of many cities around the world [96]. However, the percentage of the priority restoration zone in the suburbs was much higher than in the urban center (Figure 4a), which can be attributed to a series of urban renewal and ecological restoration projects implemented by the government in the urban center [97]. The renewal and ecological restoration of central urban areas has been a hot research topic. The environmental integrity and sustainability of the built environments were improved through extensive land renewal and the replacement of functional regions [98]. Currently, the urban center is  Table 4). The main vegetation units in these regions include artificial grassland, water bodies, and riparian woodland. As new green spaces emerge, the quality of the new ecosystems does not fully match that of the original ecology [99], implying that ecological restoration needs to identify the original habitats.
Suburban villages, as regions where built-up area is gathered, were surrounded by secondary forests and cropland (Figure 3b1-b3). Green spaces of a certain area are preserved and planned within the villages, which may explain why Type G boasts the largest area and the most extensive coverage. As rivers and lakes are prone to fragmentation after being affected by human activities [100], Type C is distributed near natural water bodies (502.99 km 2 ). Ecological restoration usually improves connectivity between fragmented patches by creating or restoring corridors and stepping stones [101]. Additionally, the implementation of ecological restoration also requires consideration of time costs due to the difference in the impact on landscape quality between the early and mid-late stages of urban development. Thus, depending on the stage, the capital investment for and implementation of ecological restoration must be adjusted [102], accelerating succession through integrated projects and restoring landscape quality to its pre-damaged state.
However, existing studies have caused controversy regarding the positive effects of ecological restoration projects on regional ecological quality. L-J Wang et al. found that the Grain for Green Program (GFGP) and the Natural Forests Conservation Program increased ecosystem values during 2000-2010 but decreased them during 2010-2015 [103]. Ecological restoration projects have a positive impact on the optimization of the structure of existing forests and the improvement in green cover [104,105]. The reasons for this could include, on the one hand, the climate and geography of the study site being highly differentiated, ensuring that the logic of ecological restoration can be applied in the face of geographically and climatically diverse regions. However, the results are prone to discrepancies [106]. On the other hand, the restoration of ecological quality has a time-lag effect [107]. The planning and implementation of ecological restoration projects are extremely complex from ecological, economic, and social perspectives [108]. In addition, there are many uncertainties in the evaluation of benefits after implementation [109].

Limitations and Future Research
In this study, the grading of ecologically degraded regions was achieved by analyzing trends in regional landscape quality and combining the transformation mechanism of land cover types. Additionally, eight types of ecological restoration requiring different strategies, were classified based on different factors that dominated the degradation. However, there were some limitations in this study.
First, Landsat series data with a resolution of 30m×30m were used in this study, and they have also been used in many other studies [110]. The differences in resolution of the research results need to be noted in subsequent studies. The remote sensing data used should be multi-sourced. For example, Sentinel-2 has higher accuracy multi-band remote sensing images, which are especially effective for monitoring vegetation health information [88]. Meanwhile, combining panchromatic bands with pan-sharpening of the original multispectral bands can effectively improve the resolution of the base data [111].
Second, this study comprehensively evaluated the trend changes in regional ecology through a multidimensional index system and identified regions with different levels of landscape quality. However, the essence of the study was still a characterization of regional ecology. Therefore, future research should explore the factors affecting ecological degradation retrospectively, such as economic factors and social development factors [112,113].
In addition, an analytical method based on geometric concepts, namely, the morphological spatial pattern analysis (MSPA), is widely used in the field of ecological sciences [16]. MSPA can visually reflect the spatial characteristics of regional ecology. In other words, the application of MSPA to analyze and identify trend changes in regional ecological spa-tial structure could directly guide the adjustment and implementation of urban planning, which should be further explored in future research.

Conclusions
Ecological degradation has always been an issue of concern for economic and social development [114]. Taking the Shanghai metropolitan area as an example, this study provided a comprehensive evaluation of regional landscape quality through indicators in three dimensions: landscape area, structure, and function. These indicators covered the main characteristics of spatial changes in the landscape while remaining independent of the scale of the study unit. Trends in landscape quality over long periods of time were analyzed using Mann-Kendall and Sen's slope estimator statistical tests. The results showed that, in the process of urban development, built-up areas in the Shanghai metropolitan area encroached on forests (506.18 km 2 ) and cropland (748.96 km 2 ) from 2007 to 2020. Furthermore, forests with stable LAI accounted for 47.71% of the whole study area and had the highest LAI stability compared to cropland and built-up area. When the built-up area and agricultural land were converted to forest, 18.63% and 50.31% of the land showed significant LSI degradation, respectively. The interconversion of forest and agricultural land affected about 10% of the land with LFI degradation, which meant that this behavior could lead to a slight degradation of landscape functions.
Based on the results of landscape quality, restoration priorities were identified, and different types of ecological restoration goals were set. We found that the suburbs warranted a much higher position of priority with respect to ecological restoration as compared to the urban center in terms of area and percentage. The most common types of ecological restoration in the Shanghai metropolitan area were Type C (covering 502.99 km 2 ), where the degradation of the landscape structure was dominant, and Type G (covering 962.478km 2 ), with a slight degradation of landscape quality. This study highlights ecological restorationoriented development planning in the Shanghai metropolitan area. The results illustrate the need to pay attention to the evolution process of landscape quality in long-term analysis under the premise of sustainable urban development. This can provide a more rapid scientific basis for ecological monitoring and urban planning.