Spatiotemporal Dynamics of Urban Green Space Inﬂuenced by Rapid Urbanization and Land Use Policies in Shanghai

: Urbanization has led to the continuous expansion of built-up areas and the ever-growing urban population, threatening the quantity and quality of urban green space (UGS). Exploring the spatiotemporal variations of UGS is substantially conducive to the formulation of land-use policies to protect the ecosystems. As one of the largest megacities all around the world, Shanghai has experienced rapid urbanization in the past three decades. Insights into how UGS changes in response to urbanization and greening policies are essential for guiding sustainable urban development. This paper employed integrated approaches to characterize the changing patterns and intensities of green space in Shanghai, China from 1990 to 2015. The spatiotemporal dynamics of the UGS pattern were derived through four main methods: green space ratio, dynamic change degree (DCD), transition matrix and landscape metrics. The results showed that Shanghai’s green space decreased from 84.8% in 1990 to 61.9% in 2015 while the built-up areas increased from 15.0% to 36.5%. Among the green space sub-types, farmland was largely encroached and fragmented by urban sprawl, especially in the Outer Ring Expressway and Suburban Ring Expressway belts of the city. About 1522 km 2 of the green space has transferred into built-up areas, followed by farmland, waterbody, forest, and grassland in descending order. The 2000–2010 period witnessed the strong urban expansion and dramatic changes in UGS, but then the change around 2015 turned down and stable. The landscape pattern metrics showed that the entire green space in Shanghai was growingly fragmented and isolated during the past 25 years. Combined with the green space-related planning and policies issued in 1990–2015, the results revealed that both rapid urbanization and greening policies accounted for the spatiotemporal dynamics of UGS. Based on the results, some implicants to new urban planning and policies of Shanghai were highlighted.


Introduction
The global urbanization rate was predicted to reach 70% by 2050 [1]. The accelerated urbanization in the world has led to the ever-expanding built-up area and ever-shrinking urban green space (UGS) [2,3], causing ecosystem and environmental problems [4][5][6]. UGS is an important component of urban ecosystems, and its degradation disrupts the structure and process of urban ecosystems and affects urban sustainability. Increasing attention has been paid to the ecosystem service of UGS, including the improvement of microclimate regulation [7], mitigation of urban heat island effects [8], urban pollution controlling [9], providing aesthetic enjoyment and entertainment opportunities and maintaining wellbeing physically and psychologically [10].
Many countries and cities formulated greening policies and UGS planning actively to protect urban ecosystems and enhance the living quality of citizens [11,12]. For example, the United Kingdom and Brazil incorporated UGS planning into the urban planning system and issued corresponding planning policies since the 1990s [13,14]. The Maryland State Government of US established land protection programs to protect natural resources and developed Green Infrastructure Assessment to identify important ecological areas [15]. Great efforts also have been made in China in recent years to protect ecosystems. Especially the principle of ecological priority was insisted on in all land-use policies. For example, the ecological red lines controlling has become one of three important contents of the national spatial planning.
To assess the effects of rapid urbanization and land use policies on changes in UGS, many have studied the spatiotemporal dynamics and changing trend of UGS to help optimize their pattern and develop proper decisions to protect the urban ecosystem [16][17][18]. Satellite remote sensing techniques and geographical information systems (GIS) have been widely applied in recognizing the land-use changes over time [19,20]. A land-use transition matrix to reflect the transition among different land-use types in a specific period [21,22], and landscape metrics to measure the spatial changes of the landscape composition and configuration, were commonly used. In addition, a few studies have integrated gradient analysis with landscape metrics to quantitatively characterize green space patterns [23][24][25].
In this study, we applied the land use data of Shanghai to quantify the spatiotemporal dynamics of UGS from 1990 to 2015. We identified the spatiotemporal pattern and trend of UGS change in terms of rate, direction, gradient, land use transition and landscape pattern. This study aims to address three points: (1) the spatiotemporal pattern and changes of UGS in Shanghai during the past 25 years, (2) the discussion on the possible planning policies affecting changes in UGS, and (3) the suggestions to the Shanghai's official urban plan to optimize the layout of UGS. Insights into this change and its driving factors therefore could benefit sustainable urban development and UGS protection in Shanghai [34].

Study Area
Shanghai, as the leading city of the Yangtze River Delta and the economic and financial center of China, is located in China's eastern coast, covering an area of about 6340 km 2 [35]. It is flat with the Huangpu River running through its heartland and is bounded by three ring roads, i.e., the Middle Ring Expressway, the Outer Ring Expressway and the Suburban Ring Expressway (Figure 1a). which isolate urban clusters and realize interconnection with urban and rural ecological space, with the forest as the main body. In addition, the ecological network planning also includes multiple ecological intervals (vertical green belts that communicate and connect the city center and the surrounding green space, and restrict the continuous development of the city inside the Outer Ring Expressway) [37] and green wedges (wedge-shaped green space embedded in the city from the periphery of the city, like a city park or square) [38].

Datasets
We used the land-use datasets interpreted from Landsat TM/ETM (National Aeronautics and Space Administration, Washington, USA) and Landsat OLI (National Aeronautics and Space Administration, Washington, USA) images for 1990,1995,2000,2005,2010, and 2015. The land-use datasets used in the research had already been supervised classified and widely used in published research [39][40][41]. The satellite images are classified using the supervised maximum likelihood classification method in ENVI 5.2 (https://www.l3harrisgeospatial.com/docs/whatsnew_envi52.html) to generate land-use patterns and the overall accuracies of the classified land use data are above 97.3% [41]. With the development of the Pudong New Area, Shanghai's economy has been growing rapidly since 1990. The population also increased from 13.3 million in 1990 to 24.2 million in 2015 [36], which promoted urban growth together with economic prosperity. Shanghai has become one of the fastest-growing and most densely populated cities in China and the world. The urban expansion has resulted in the fragmentation and reduction of UGS.
To prevent the loss of UGS caused by urbanization and maintain ecological security, Shanghai has recently proposed multi-level, networked and functional urban ecological space systems, which are composed of "a double-ring, nine ecological corridors, ten ecological conservation districts" as issued by the "Shanghai Urban Master Plan (2017-2035)" (Figure 1b). The double-ring belts indicate the green belt around the central area (outer ring green belt) and that around the far suburbs (suburban green belt). Green belts refer to green space to prevent the blind expansion of the city. This kind of belt can be scenic area, forests, etc. Ecological corridors refer to radial and unobstructed corridors in the city, which isolate urban clusters and realize interconnection with urban and rural ecological space, with the forest as the main body. In addition, the ecological network planning also includes multiple ecological intervals (vertical green belts that communicate and connect the city center and the surrounding green space, and restrict the continuous development of the city inside the Outer Ring Expressway) [37] and green wedges (wedge-shaped green space embedded in the city from the periphery of the city, like a city park or square) [38].

Datasets
We used the land-use datasets interpreted from Landsat TM/ETM (National Aeronautics and Space Administration, Washington, DC, USA) and Landsat OLI (National Aeronautics and Space Administration, Washington, DC, USA) images for 1990,1995,2000,2005,2010, and 2015. The land-use datasets used in the research had already been supervised classified and widely used in published research [39][40][41]. The satellite images are classified using the supervised maximum likelihood classification method in  24 February 2021) to generate land-use patterns and the overall accuracies of the classified land use data are above 97.3% [41]. Additionally, we reclassify the land-use into green space and non-green space (Table 1) according to the land use classification criteria of the Chinese Academy of Sciences [42]. The green space includes farmland (paddy field, dry land), forest (forest land, shrub, woods, others), grassland (dense grass, moderate grass, sparse grass) and waterbody (stream and rivers, lakes, reservoir and ponds, permanent ice and snow, beach and shore, swampland, bottomland), and the non-green space includes the built-up area (urban built-up, rural settlements, others), and unused land (Sandy land, Gobi, bare soil, bare rock, others). We then mapped the land use in Shanghai for 1990Shanghai for , 1995Shanghai for , 2000Shanghai for , 2005Shanghai for , 2010, and 2015 using ArcGIS 10.2 ( Figure 2). Significant loss of green space and expanding built-up area can be clearly seen in Figure 2. Table 1. The land-use classification criteria applied in this study [42]. Additionally, we reclassify the land-use into green space and non-green space (Table 1) according to the land use classification criteria of the Chinese Academy of Sciences [42]. The green space includes farmland (paddy field, dry land), forest (forest land, shrub, woods, others), grassland (dense grass, moderate grass, sparse grass) and waterbody (stream and rivers, lakes, reservoir and ponds, permanent ice and snow, beach and shore, swampland, bottomland), and the non-green space includes the built-up area (urban built-up, rural settlements, others), and unused land (Sandy land, Gobi, bare soil, bare rock, others). We then mapped the land use in Shanghai for 1990Shanghai for , 1995Shanghai for , 2000Shanghai for , 2005Shanghai for , 2010, and 2015 using ArcGIS 10.2 ( Figure 2). Significant loss of green space and expanding builtup area can be clearly seen in Figure 2. Table 1. The land-use classification criteria applied in this study [42].

New Category Initial Category Description
Green space (GS)

GS1: Farmland Cropland
Cultivated lands for crops, which including mature cultivated land, newly cultivated land, fallow, shifting cultivated land.
GS2: Forest Woodland Lands growing trees including arbor, shrub, bamboo and forest.

GS3: Grassland Grassland
Lands covered by the herbaceous plants with coverage greater than 5%, including shrub rangeland and mixed rangeland with the coverage of shrub canopies less than 10%.

GS4: Waterbody Waterbody
Lands covered by natural water bodies or lands with facilities for irrigation and water reservation.

Non-green space
Built-up area Built-up area Lands used for urban and rural settlements, factories and transportation facilities.
Unused land Unused land Lands that are not put into practical use or difficult to use.

Methods
To analyze the spatiotemporal dynamics of UGS in Shanghai, we proposed an analysis method composed of four main parts ( Figure 3).
12, x FOR PEER REVIEW 5 of 21

Methods
To analyze the spatiotemporal dynamics of UGS in Shanghai, we proposed an analysis method composed of four main parts ( Figure 3).

Green Space Ratio
To deeply capture the spatiotemporal changes of UGS in Shanghai from 1990 to 2015, we analyzed the UGS in three aspects: quantitative changes of the whole city, eight directional changes and concentric buffer-based changes.
The directional and concentric analysis methods are effective for characterizing the quantity and spatial distribution of various types of UGS in terms of orientation and distance with respect to a pre-determined urban center [31,43]-in this case the People's Square, the center of the circle. The directional analysis method can help to study the impacts of the urban construction in different directions on the green space and the concentric analysis can reflect the gradient changes of UGS from the dense urban construction area to vast suburban. We divided the study area into eight sectors, i.e., north, northeast, east, southeast, south, southwest, west and northwest (Figure 4a), for the directional analysis. While in the concentric analysis, the study area was divided into 15 concentric buffers from the center point with a width of 5 km (Figure 4b). Based on these two methods, we analyzed UGS changes and the green space ratio of each type. These were displayed on radar graphs and growth curves to show the spatiotemporal changes of the overall green space and its different subtypes in the eight directions or from the city center point to the boundary.

Green Space Ratio
To deeply capture the spatiotemporal changes of UGS in Shanghai from 1990 to 2015, we analyzed the UGS in three aspects: quantitative changes of the whole city, eight directional changes and concentric buffer-based changes.
The directional and concentric analysis methods are effective for characterizing the quantity and spatial distribution of various types of UGS in terms of orientation and distance with respect to a pre-determined urban center [31,43]-in this case the People's Square, the center of the circle. The directional analysis method can help to study the impacts of the urban construction in different directions on the green space and the concentric analysis can reflect the gradient changes of UGS from the dense urban construction area to vast suburban. We divided the study area into eight sectors, i.e., north, northeast, east, southeast, south, southwest, west and northwest (Figure 4a), for the directional analysis. While in the concentric analysis, the study area was divided into 15 concentric buffers from the center point with a width of 5 km (Figure 4b). Based on these two methods, we analyzed UGS changes and the green space ratio of each type. These were displayed on radar graphs and growth curves to show the spatiotemporal changes of the overall green space and its different subtypes in the eight directions or from the city center point to the boundary.

Land Use Transition Matrix
The land-use transition matrix reflects the structure and state of transitions of landuse types between every two periods [44]. By comparing the land-use patterns of two different periods pixel by pixel, we can identify the type and area transition of each landuse type. The land-use transition matrix can be given by [45]: where Pij indicates the transitional area from land-use type i to type j, where each element in the matrix is characterized assuming that Pij is non-negative and ∑ =1.

Dynamic Change Degree (DCD)
The DCD of UGS is a kind of land-use change rate index, which reflects the change density and rate of a certain green space type in a period [46]. The DCD of UGS can be calculated as [46]: where K indicates the DCD of UGS, Ua and Ub, respectively indicate the area of UGS at the start and end of a period, and T indicates the period in years.

Landscape Metrics
Landscape metrics are often used to detect and quantify the spatiotemporal changes in landscape composition and configuration [34,47,48]. In this study, we selected eight landscape metrics including five class-level metrics and three landscape-level metrics by referencing earlier work reported in many research (Table 2) [44,49]. They are used to reflect fragmentation, connectivity and diversity of landscape patterns [50]. The five classlevel metrics include: percentage of landscape (PLAND) depicts the relative abundance of each land-use type; patch density (PD) represents the landscape fragmentation; large patch index (LPI), landscape shape index (LSI) illustrates landscape dominance and patch cohesion index (COHESION) reflects the connectivity in patches. Three landscape-level metrics are: contagion index (CONTAG) describes the aggregation degree of different patch types in the landscape, Shannon's diversity index (SHDI) shows landscape diversity

Land Use Transition Matrix
The land-use transition matrix reflects the structure and state of transitions of land-use types between every two periods [44]. By comparing the land-use patterns of two different periods pixel by pixel, we can identify the type and area transition of each land-use type. The land-use transition matrix can be given by [45]: where P ij indicates the transitional area from land-use type i to type j, where each element in the matrix is characterized assuming that P ij is non-negative and ∑ n j=1 P ij =1.

Dynamic Change Degree (DCD)
The DCD of UGS is a kind of land-use change rate index, which reflects the change density and rate of a certain green space type in a period [46]. The DCD of UGS can be calculated as [46]: where K indicates the DCD of UGS, U a and U b , respectively indicate the area of UGS at the start and end of a period, and T indicates the period in years.

Landscape Metrics
Landscape metrics are often used to detect and quantify the spatiotemporal changes in landscape composition and configuration [34,47,48]. In this study, we selected eight landscape metrics including five class-level metrics and three landscape-level metrics by referencing earlier work reported in many research (Table 2) [44,49]. They are used to reflect fragmentation, connectivity and diversity of landscape patterns [50]. The five classlevel metrics include: percentage of landscape (PLAND) depicts the relative abundance of each land-use type; patch density (PD) represents the landscape fragmentation; large patch index (LPI), landscape shape index (LSI) illustrates landscape dominance and patch cohesion index (COHESION) reflects the connectivity in patches. Three landscape-level metrics are: contagion index (CONTAG) describes the aggregation degree of different patch types in the landscape, Shannon's diversity index (SHDI) shows landscape diversity and heterogeneity, and Shannon's evenness index (SHEI) describes the distribution of various patches in the landscape. Among these, PLAND, PD, LPI and LSI can represent fragmentation; COHESION and CONTAG reflect the connectivity in the green space sub-types and landscape, respectively; SHDI and SHEI can show the diversity. Table 2. Landscape metrics used in this study [48].

Class-level
Percentage of landscape (PLAND)

Percent
Fragmentation and dominance index Landscape shape index (LSI) ik total length (m) of edge in the landscape between patch types (classes) i and k; includes the entire landscape boundary and some or all background edge segments involving class i A = total landscape area (m 2 ).

None
Shape index Patch cohesion index (COHESION) (p * i ln p i ) P i = proportion of the landscape occupied by patch type (class) i .

None
Diversity index Shannon's evenness index (SHEI)  Table 3 shows Shanghai's green space area has been decreasing from 1990 to 2015, while non-green space has been increasing. The ratio of green space declined from 84.8% in 1990 to 61.9% in 2015 while the built-up area increased from 15.0% to 36.5%. Among the various types of green space, the loss of farmland was the most intense, the proportion of forest remained relatively stable, the grassland increased, and the waterbody showed an overall decreasing trend. Farmland experienced a substantial reduction from 71.3% to 52.9%, but still shared the largest proportion than any other type of Shanghai's green space. Therefore, over 80% of the lost green space is contributed by the farmland.   Figure 5 shows that the overall green space in Shanghai was evenly distributed in eight directions in 1990, but the distribution of subtypes of green space in eight directions in the city was quite different. The overall green space in all directions declined over time, and the changes of the four green space types in eight directions had significant differences.

Green Space Changes along the Concentric Buffers
Since the mainland of Shanghai is enclosed in buffer 12 and green space changed significantly (Figure 4b), we focused on the UGS between Buffer-1 and Buffer-12 where built-up areas of Shanghai are mostly distributed in this study. The UGS and built-up area show an obvious gradient change with a reversal trend. Figure 6 shows that the green space ratio of each buffer increased rapidly from the city center to the suburbs, especially between Buffer-2 and Buffer-4, but tended to be stable outside of Buffer-9 because there In 1990, Shanghai's green space accounted for similar areas in all directions, suggesting that the UGS was roughly evenly distributed, which may be greatly attributed to the farmland. With the rapid urbanization, the proportion of green space decreased in all directions from 1990 to 2015, especially in the east-west direction. The distribution and change of farmland were similar to the overall green space with a major loss in the eastwest direction. In contrast, the built-up area shows substantial expansion in all directions, particularly in the east and west, indicating that encroachment on suburban farmland is the main way of Shanghai's urbanization.
Forest changed slightly in all directions but was mainly distributed in the north (Chongming Island), the southwest (Sheshan National Forest Park) and southeast (Fengxian Forest park). The grassland is the only type that experienced an increase in all directions, and a notable increase was particularly observed in the southeast after 2005. Probably because Lingang New Town began to be developed. Water bodies are mainly located in the north to the east, where the Yangtze River runs through Shanghai in the north and Shanghai has a long coast in the east. However, the water bodies in the southeast and east reduced greatly, which may be attributed to the tidal flat reclamation, converting many sea areas to useable land. The water bodies in the northeast and north also showed an evident loss due to sediment deposition, leading to the inland expansion in the Chongming Island.

Green Space Changes along the Concentric Buffers
Since the mainland of Shanghai is enclosed in buffer 12 and green space changed significantly (Figure 4b), we focused on the UGS between Buffer-1 and Buffer-12 where built-up areas of Shanghai are mostly distributed in this study. The UGS and built-up area show an obvious gradient change with a reversal trend. Figure 6 shows that the green space ratio of each buffer increased rapidly from the city center to the suburbs, especially between Buffer-2 and Buffer-4, but tended to be stable outside of Buffer-9 because there are many satellite towns and new towns beyond Suburban Ring Expressway. From the time dimension of view, the green space ratio of each buffer continued to decrease from 1990 to 2015 with obvious growth of built-up area occurred, especially between Buffer-2 and Buffer-9. Grassland and water bodies have a higher rate in Buffer-6 because upstream of the Huangpu River and its Dazhi river runs along with this buffer with many wetlands. However, both of them increased quickly outside of Buffer-9, which means the typical water towns are well kept in this region. During 1990-2000, the grassland and water bodies changed minor in each buffer zones; whereas, after 2000, the proportion of grassland increased significantly from Buffer-10 to Buffer-12 while waterbody declined dramatically.
Forest is mainly concentrated around Buffer-6, where the over 500 m width green belt was developed since 1990 and many ecological corridors are distributed. However, it then declined away from the center city. It seemed that forest development had not been taken seriously.

Transitions between Green Space and Built-Up Area
The land-use transition matrixes (Table 4) reflect the spatiotemporal transition between the sub-types of green space and the built-up area. From 1990 to 2015, farmland was the sub-type of green space that transferred out most, of which 1428.3 km 2 were converted into built-up area, accounting for 95.3% of the total loss of farmland. The transition from water bodies to the built-up area was also relatively large (73.0 km 2 ), accounting for  The concentric buffers roughly parallel to the three ring-expressways in Shanghai. The Middle Ring Expressway is located between Buffer-2 and Buffer-3, the Outer Ring Expressway is located between Buffer-3 and Buffer-4, and the Suburban Ring Expressway spans from Buffer-5 to Buffer-9. The outer and suburban ring expressway zoned Shanghai into 3 belts. The center part inner the Outer Ring Expressway is the central city where the green space ratio is low and rises slowly. The green space especially farmland between the Outer and Suburban Ring Expressways was encroached by the built-up area. While the belt out of the Suburban Ring Expressway, waterbody and grassland are the main green space types, which shows the characteristic of "water village" in the Middle-Lower Yangtze Plain.
Among the green space sub-types, farmland dominated all zones outside of Buffer-3. Figure 6 shows that the spatiotemporal dynamics and distribution of farmland are similar to the total green space, because farmland takes up the highest proportion of green space. The farmland's proportion shows an overall upward trend from Buffer-1 to Buffer-9, except from Buffer-6 where farmland decreases slightly. While outside of Buffer-9, the proportion of farmland shows a downward trend, and from 1990 to 2015 the farmland density of each buffer zones declines with time, especially in the regions between Buffer-3 and Buffer-7.
Grassland and water bodies have a higher rate in Buffer-6 because upstream of the Huangpu River and its Dazhi river runs along with this buffer with many wetlands. However, both of them increased quickly outside of Buffer-9, which means the typical water towns are well kept in this region. During 1990-2000, the grassland and water bodies changed minor in each buffer zones; whereas, after 2000, the proportion of grassland increased significantly from Buffer-10 to Buffer-12 while waterbody declined dramatically.
Forest is mainly concentrated around Buffer-6, where the over 500 m width green belt was developed since 1990 and many ecological corridors are distributed. However, it then declined away from the center city. It seemed that forest development had not been taken seriously.

Transitions between Green Space and Built-Up Area
The land-use transition matrixes (Table 4) reflect the spatiotemporal transition between the sub-types of green space and the built-up area. From 1990 to 2015, farmland was the sub-type of green space that transferred out most, of which 1428.3 km 2 were converted into built-up area, accounting for 95.3% of the total loss of farmland. The transition from water bodies to the built-up area was also relatively large (73.0 km 2 ), accounting for 25.4% of the transition of water bodies. The built-up area transferred from the forest was fewer (16.2 km 2 ), but accounted for 81.0% of the losses of forest. The area transferred from grassland to built-up area was the least (4.6 km 2 ). The transition was also found from built-up area to green space (26.3 km 2 ), such as farmland and forest, and the quantity was generally increasing from 1990 to 2015.
The transition among different green space types was predominantly from water bodies to farmland. The past 25 years have witnessed 286.93 km 2 reduction in water body, of which 59.8% were turned into farmland. With the increasing demand for urban space, tidal flats were considered the most important land reserves in Shanghai. Considerable tidal flats along the Pudong New Area's coasts were developed and transferred into builtup area. Another feature was the increasing transitional area from farmland to forest and grassland after 2000. The transition area during 2010-2015 increased by more than 70 times of those during 2005-2010, which might be attributed to the implementation of the greening policy in Shanghai.  Figure 7 shows the spatial transition from the green space to the built-up area during 1990 and 2015, where the farmland around the urban fringes intensively was transferred into built-up area. The area under transition expanded outward, which revealed the encroachment of urbanization on farmland over time. The decrease of water bodies mainly occurred in the eastern and southern coastal area. The same result can also be drawn from Figure 5. However, the grassland and forest transferred to built-up area was distributed dispersedly. The built-up area exhibits a radial growth pattern from the city center to the fringes and small towns in far suburbs. The UGS increased greatly with the distance from the city center.    Figure 8 shows that the DCD was relatively small from 1995 to 2000 for the overall green space, but the reduction rate of green space rose during 2000-2010 while the builtup area always maintained a high growth rate. In the period of 2005-2010, DCD of green space, built-up area and sub-types of green space reached the highest point which witnessed strong urban expansion and green space change rate, and after 2010 the change rate of both the green space and built-up area turned down to be stable again. The result is a reflection of the urbanization process of shanghai. Shanghai experienced a quick development during 2000-2010, and then began to slow down the construction rate and emphasize high-quality development after 2010.

Dynamic Change Degree (DCD) of Green Space
The DCD of the sub-types of UGS is quite different. Most importantly, farmland shows a similar trend as the overall green space. The dynamic degree of the forest was kept between −2% and 2% per year. Grassland got the fast increase rate to 14% per year during 2005-2010. However, every sub-type turned to be stable in the last period (2010-2015) except that the water bodies decreased sharply at 7.5%.

Green Space Changes in Landscape Patterns
We applied Fragstats 4.2 to calculate the landscape-level metrics of the entire landscape, that is, all land-use types (forest, grassland, farmland, waterbody, build-up area, unused land). The contagion index (CONTAG) shows a consistent decreasing trend in Figure 9a. The negative trend of CONTAG suggested high landscape fragmentation and low landscape connectivity in Shanghai. As shown in Figure 9b,c, Shannon's diversity index (SHDI) and Shannon's evenness index (SHEI) have been increasing during 1995-2015. The increasing SHDI and SHEI of green space indicate the increase of landscape fragmentation because of the rapid urban expansion and the reduction of the originally dominant land-use type (farmland). These in turn strengthened the relative balance of various landscape patches. The DCD of the sub-types of UGS is quite different. Most importantly, farmland shows a similar trend as the overall green space. The dynamic degree of the forest was kept between −2% and 2% per year. Grassland got the fast increase rate to 14% per year during 2005-2010. However, every sub-type turned to be stable in the last period (2010-2015) except that the water bodies decreased sharply at 7.5%.

Green Space Changes in Landscape Patterns
We applied Fragstats 4.2 to calculate the landscape-level metrics of the entire landscape, that is, all land-use types (forest, grassland, farmland, waterbody, build-up area, unused land). The contagion index (CONTAG) shows a consistent decreasing trend in Figure 9a. The negative trend of CONTAG suggested high landscape fragmentation and low landscape connectivity in Shanghai. As shown in Figure 9b,c, Shannon's diversity index (SHDI) and Shannon's evenness index (SHEI) have been increasing during 1995-2015. The increasing SHDI and SHEI of green space indicate the increase of landscape fragmentation because of the rapid urban expansion and the reduction of the originally dominant land-use type (farmland). These in turn strengthened the relative balance of various landscape patches. We used Fragstats 4.2 to calculate every class-level metric of UGS (forest, grassland, farmland, waterbody) and built-up area. Figure 10 shows the changes in the class-level metrics for all sub-types of the UGS during 1990-2015, where the farmland remained the dominant according to percentage of landscape (PLAND) and large patch index (LPI). However, the urban land use pattern changed dramatically in response to the continuous interference of human activities (e.g., urban development and construction). The increasing patch density (PD) for farmland indicates that the urban development accompanying by rapid urbanization has resulted in the fragmentation of farmland. The patch cohesion index (COHESION) of farmland kept steady and relatively higher than other sub-type, which means good connectivity between the farmland patches. The increasing Landscape shape index (LSI) indicates the gradual complexity of farmland patches and a decrease in their aggregation. During the study period, the built-up area became another dominant landscape type in Shanghai as revealed by higher LPI of built-up area compared to farmland since 2010. PD and LSI of forest decreased while COHESION increased, which might be attributed to the artificial forests (e.g., the Shanghai Bay National Forest Park) that promoted We used Fragstats 4.2 to calculate every class-level metric of UGS (forest, grassland, farmland, waterbody) and built-up area. Figure 10 shows the changes in the class-level metrics for all sub-types of the UGS during 1990-2015, where the farmland remained the dominant according to percentage of landscape (PLAND) and large patch index (LPI). However, the urban land use pattern changed dramatically in response to the continuous interference of human activities (e.g., urban development and construction). The increasing patch density (PD) for farmland indicates that the urban development accompanying by rapid urbanization has resulted in the fragmentation of farmland. The patch cohesion index (COHESION) of farmland kept steady and relatively higher than other sub-type, which means good connectivity between the farmland patches. The increasing Landscape shape index (LSI) indicates the gradual complexity of farmland patches and a decrease in their aggregation. During the study period, the built-up area became another dominant landscape type in Shanghai as revealed by higher LPI of built-up area compared to farmland since 2010. We used Fragstats 4.2 to calculate every class-level metric of UGS (forest, grassland, farmland, waterbody) and built-up area. Figure 10 shows the changes in the class-level metrics for all sub-types of the UGS during 1990-2015, where the farmland remained the dominant according to percentage of landscape (PLAND) and large patch index (LPI). However, the urban land use pattern changed dramatically in response to the continuous interference of human activities (e.g., urban development and construction). The increasing patch density (PD) for farmland indicates that the urban development accompanying by rapid urbanization has resulted in the fragmentation of farmland. The patch cohesion index (COHESION) of farmland kept steady and relatively higher than other sub-type, which means good connectivity between the farmland patches. The increasing Landscape shape index (LSI) indicates the gradual complexity of farmland patches and a decrease in their aggregation. During the study period, the built-up area became another dominant landscape type in Shanghai as revealed by higher LPI of built-up area compared to farmland since 2010. PD and LSI of forest decreased while COHESION increased, which might be attributed to the artificial forests (e.g., the Shanghai Bay National Forest Park) that promoted PD and LSI of forest decreased while COHESION increased, which might be attributed to the artificial forests (e.g., the Shanghai Bay National Forest Park) that promoted the aggregation and connectivity of the forests. The increasing PD and LSI of grassland suggest that grassland was interfered out due to intensive human activities. For water bodies, the PD and LSI show an upward trend while the LPI and COHESION display a downward trend, indicating more scattered, more complicated, increasingly fragmented, and poorly connected water patches.

Integrated Approaches to Describe Spatiotemporal Dynamics of UGS
In this study, we used time-series of land-use datasets between 1990 and 2015, which help to find the green space changes in a long time. We also integrated several methods of spatiotemporal dynamics analysis in order to understand well about UGS spatial pattern and change trend in detail.
Green space ratio is a general index in reflecting the quantity of UGS. The study demonstrated the usefulness of the concentric and directional green space ratio analysis in characterizing the spatiotemporal variations of UGS [51], which help us to find the UGS of Shanghai lost mainly in the east-west direction and changed along the urbanrural gradient. Unfortunately, the directional analysis showed that the built-up area also increased in the east-west direction substantially. This may threaten the Dianshan Lake area in the west of Shanghai where many small-area water bodies and wetlands serve as very important ecological reserves in Shanghai. The expressway became an important restricted boundary to the city development and showed different characteristics of UGS by the concentric analysis. The classic concentric zone urban development theory seems to account adequately for the land-use pattern of the Shanghai metropolitan region [32].
DCD reflected that green space types and built-up area development were extremely uneven during the rapid urbanization, but their ups and downs matched the urbanization process of Shanghai precisely.
The landscape pattern analysis and transition matrix has been widely applied to analyze the urban land use pattern quantitatively [22]. Since the landscape metrics have scale effects at both the landscape and class levels [27] and many metrics relate with each other, minimal indicators were encouraged to be used to characterize the UGS [29]. We selected three landscape-level metrics and five class-level metrics to reflect fragmentation, connectivity and diversity of landscape pattern. Similar to numerous cities around the world reported in previous studies, urbanization increased the diversity, fragmentation, and configurational complexity of the urban landscape of Shanghai [29,32]. The green space in Shanghai was sacrificed to meet the demand for urbanization and industrialization accompanied by population growth and the increase in demand for their substance. The initial farmland occupied by built-up area suggests the major reason for green space loss, especially between Outer and Suburban Ring Expressway. The grassland and forest in plain were also turned into the built-up area. However, the magnitude was relatively small, and the positions were relatively scattered. This result is the same as other researches [34].

The Impacts of Policies on UGS
Many factors drive the changes of UGS [52]. Greening and land-use policies are the major driving factors that can restrict the negative impacts of urbanization on the UGS. Some strategic land-use policies can even lead to a substantial increase in UGS [53][54][55]. The Shanghai government has proposed many policies since the 1980s [33], which played a substantial role on green space.
(1) The protection policies of farmland restricted its massive losses that may be caused by urban expansion. For instance, as shown in Table 3, about 327 km 2 of farmland was lost to built-up area in 1990-1995, but in 1995-2000 the losses of farmland were only 23.3% of those in 1990-1995. This could be attributed to the strict implementation of the "Basic Farmland Protection Policy" issued by the Shanghai local government in 1996. In 1998, a proposal encouraging land consolidation was put forward, leading to the growth of the opposite transition from built-up area to farmland in the coming period 2000-2015 according to Table 4. Supported by this proposal, the reclamation of abandoned and inefficient built-up areas along with the comprehensive rehabilitation of farmlands, water bodies, roads, forests, and villages effectively increased the farmlands. In the meanwhile, the State Council of the PRC has promulgated the "Regulations on the Protection of Basic Farmland" in 1999, and put forward the requirements for delimiting basic farmland protection areas to protect the cultivated land from the institutional point of view. It has also been repeatedly emphasized in the "Shanghai Land-use Master Plan (2006-2020)" that in order to ensure urban food security and the health of ecological systems, it is necessary to strictly protect the basic farmland, which may slow down the reduction of the farmland in 2010-2015 (Table 3).
(2) The implementation of greening policies and master plan also had positive impacts on afforestation and planting grass. In the early 1990s, Shanghai began to pay attention to the protection of green space across the jurisdictions and implemented the "Shanghai UGS System Plan (1994-2010)". It can be clearly seen in Table 3 that after 2000, both forest and grassland showed an increasing trend. Although forest began to decrease after 2005, the overall area remained basically unchanged. This may be due to the "Shanghai Urban Master Plan (1999-2020)" issued in 1999. It proposed to prevent and control the possible adverse effects of urban development on the natural environments. The local government also promoted policies for increasing the urban public greening land and suburban forests. After 2000, the Shanghai local government emphasized more on the improvement of urban and rural ecological environments by implementing greening policies, including the "Shanghai UGS System Plan (2002-2020)" implemented in 1999 and the "Shanghai Basic Ecological Network Space Plan" issued in 2009. Due to this reasonable green space planning, the aggregation and connectivity of the grassland patches have been significantly improved after 2005 according to Figure 10a. The ecological structure, as well as the ecological benefits of UGS, has been further improved. In the meanwhile, the grassland had a notable increase in the southeast according to Figure 5d, and increased significantly in Buffer-10 to Buffer-12 according to Figure 6d, because 2005 was the year that the Lingang New Town began to be developed. The sea reclamation has increased part of the grassland area. In summary, all these planning policies have an influence on the area of forest and grassland.
(3) The composition of green space and the conversion of its sub-types were also influenced by the policies. After 2000, especially in 2000-2005, more farmlands transferred to forests (Table 4), which might be related to the decisions made by the State Council of the PRC to expand reforestation on consideration of the nation's grain surplus in 1998. A strategy named "Regulations on Conversion of Farmland to Forests" issued in 2002 might also contribute to the farmland transferring outward. Shanghai's government has seized this opportunity and made the "Shanghai Forest Planning" in 2003 to accelerate artificial forestry construction and make the proper spatial layout of the forest network, which accordingly improved the aggregation and connectivity of the forest landscape patches.
In summary, these policies are important driving factors for the dynamics of Shanghai's green space, where the balance between green space and urbanization can be maintained under the regulations of proper local and national policies.

Optimizing Green Space Patterns through Urban Planning
Accurately examining the spatiotemporal process of green space changes can provide essential information for improving the landscape and urban planning [22]. Particularly, it is urgent and necessary for the sustained development and management of Shanghai to explore the hotspot areas of disappeared green space and the unstable area during the last three decades. Considering the characteristics and change trend of the green space in Shanghai, something should be highlighted in new green space planning.
First, increase or maintain green space by reducing farmland loss and strengthening afforestation outside the Suburban Ring Expressway. The average decrease rate of the green space between 1990 and 2015 was about 1%. The continuous change would make Shanghai face the risk of ecological space loss. If Shanghai wants to maintain green space in more than 50% of all land-use by 2035, the average annual reduction of green space area should not exceed 0.4% in 2015-2035. Thus, improving the efficiency of inefficient construction land to meet urbanization demand later instead of occupying farmland is a rational way at present. Meanwhile, the construction of more coastal shelterbelt and suburb forest parks will improve the situation of poor connectivity in the forest.
Second, optimizing landscape patterns to cure the unhealthy structure and excessive fragmentation of the UGS. Only when green space has good connectivity and a certain area can it effectively play its important ecological service function. Therefore, optimizing and adjusting the composition type and landscape pattern of UGS will be an important content of UGS construction in the future. It should be stressed to promote the permanent basic farmland to be larger parcels, build a stronger network of ecological space, improve the connectivity of UGS, and control the urban expansion in the east-west direction.
Third, to control the decreasing intensity of water bodies and the serious reclamation on coastal wetlands. It is necessary to control the tidal flat reclamation and maintain the coastal resources to ensure the strategic space for long-term development. Considering the ecological systems, we should strengthen the protection of the coastal shelterbelt, wetland and aquatic plant, to improve the environmental quality of coastal areas.

Conclusions
This study observed the sharp decreases of green space in Shanghai from 1990 to 2015 based on land-use datasets, which was largely due to the rapid urbanization with the population explosion. The green space losses mainly were concentrated in the suburban areas between the Outer Ring Expressway and the Suburban Ring Expressway. The pattern of overall green space was increasingly fragmented. The sub-types of the green space present different changes in various directions along the urban-rural gradients. The dynamic changes are most significant in the period of 2000-2010. Urbanization and related policies are the main factors driving the spatiotemporal dynamics of UGS.
Shanghai's new urban planning is conducive to curbing the decline of green space. We suggest strictly implementation of conserving farmland and wetland, increasing forest land and optimizing green space patterns to improve the ecological benefits of UGS.
The methods used in this study to describe the spatiotemporal dynamics of the green space are useful to understand the dynamic changes of UGS and provides a basis for future green space and ecological space planning. The land use datasets interpreted from Landsat images have good temporality and can be used to analyze the dynamic changes of green space in long time series. Overall, this study provides references for analyzing and optimizing green space in other cities experiencing similar rapid urbanization.
Due to the limitations of data resolution and the classification method, a few small green space patches in the downtown area were not well identified to be included in the analysis, which could affect the result accuracy. Future work should focus on the analysis using high-resolution satellite images that may provide more details and multi-scale investigation of UGS.