Next Article in Journal
Landslide Susceptibility Assessment in Nepal’s Chure Region: A Geospatial Analysis
Previous Article in Journal
Deciphering the Relationship between Regional Production Factors, Governance, and the Economic Performance of Metropolitan Areas in China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Study on the Spatial and Temporal Evolution of Regional Green Space Morphology Outside Built-Up Areas based on the Google Earth Engine and Biophysical Component Modeling

1
Art School, Anhui Jianzhu University, Hefei 230022, China
2
Anhui Cultural Tourism Innovative Development Research Institute, Hefei 230022, China
3
Key Laboratory of National Forestry and Grassland Administration on Ecological Landscaping of Challenging Urban Sites, Shanghai Academy of Landscape Architecture Science and Planning, Shanghai 200232, China
4
College of Architecture and Urban Planning, Tongji University, Shanghai 200092, China
5
State Key Laboratory of Hydraulic Engineering Simulation and Safety, School of Civil Engineering, Tianjin University, Tianjin 300072, China
6
China Institute of Water Resources and Hydropower Research, Beijing 100044, China
7
School of Architecture and Planning, Anhui Jianzhu University, Hefei 230022, China
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Land 2023, 12(12), 2184; https://doi.org/10.3390/land12122184
Submission received: 3 October 2023 / Revised: 13 December 2023 / Accepted: 15 December 2023 / Published: 18 December 2023
(This article belongs to the Topic Low Carbon Economy and Sustainable Development)

Abstract

:
The spatial pattern of regional green space is an important dimension to describe and quantitatively express the characteristics of regional green spaces outside the built-up area of a city. With the expansion of urban and rural construction land, regional green space has been continuously encroached upon. This leads to a decline in regional ecological well-being and the loss of biodiversity. Based on the remote sensing data of Shanghai city from 2000 to 2020, we quantitatively studied the spatial morphological change characteristics of regional green space outside the built-up area of Shanghai city. Firstly, with the help of the GEE platform, the optimal decoding accuracy classification method was selected through machine learning (random forest, support vector machine, classification regression tree); then, based on the biophysical component (BCI) and CA binarization, the built-up area ranges for up to five time nodes were obtained; finally, through GIS spatial data analysis and processing technology, the regional green space dynamic data of Shanghai for five time nodes were extracted. Based on the above data, an analysis index system was constructed to quantitatively analyze the spatial morphology characteristics of the regional green space outside the built-up area of Shanghai. The results show that (1) the area of regional green space outside the built-up area of Shanghai had a fluctuating growth pattern of “decreasing and then increasing”. The arable land and water areas in Shanghai decreased, and the woodland area increased steadily, while the wetland and grassland areas showed a trend of first decreasing and then increasing. (2) The regional green patch fragmentation shows a fluctuating development trend of increasing, decreasing, and increasing. (3) The change in the spatial center of gravity of the regional green space in Shanghai had a high degree of consistency with the overall green space change. The center of gravity of the grasslands in the regional green space moved substantially to the northwest, while the center of gravity of the other types remained basically unchanged. This study reveals the spatial morphology characteristics of regional green spaces and provides a research method to study the dynamic changes in regional ecological resources. The results of this study can provide a scientific basis for the identification, protection, and development of regional ecological resources.

1. Introduction

Rapid urban population agglomeration has led to the continuous low-quality spread of urban construction land to the suburbs and the continuous erosion of green resources outside the built-up area of Shanghai [1]. This results in the destruction of the original relatively stable and intact green resources, thus reducing the regional ecosystem’s service function and seriously affecting the health and well-being of urban and rural residents [2]. Therefore, this study focused on the regional green space outside the urban built-up area, which is defined in the Standard for classification of urban green space (CJJ/T 85-2017) of the People’s Republic of China [3] as non-constructed land located outside the built-up area, and is divided into four main types: woodland, grassland, wetland, and all three combined. Regional green space is a collection of important ecological resources that support the operation of regional ecosystem services and directly or indirectly affect human well-being, including climate regulation, the production and protection of regional facilities, seedling production, recreation and leisure spaces, and the provision of habitats for species [4]. Currently, the regional green space outside the built-up area suffers from the impact of urban expansion, especially the loss of regional green space in economically developed areas [5]. Shanghai, as the most economically developed area in China, has experienced dramatic changes in the surrounding area of regional green space after rapid urban expansion due to the dual role of economic development and the rise of the concept of ecological civilization [6].
The concept of regional green space appears less frequently in theory and practice outside of China, but land with the characteristics of regional green space is more represented in the relevant green concepts of various countries, such as the greenway network in North America [7], the pan-European ecological network planning [8], the UK’s idyllic city [9], the UK Green Belt Act (Green Belt Act) [10], etc. These concepts are more reflective of regionally integrated conservation features [11]. In different countries and at different research scales, the form of regional green space also varies greatly, and the scope of research includes both micro, single green space resources, such as the United States Central Park [12], Hong Kong country parks [13], the United Kingdom Green Belt [14], etc., and macro-scale combinations of a variety of regional green space resources across several administrative jurisdictions, such as the greenway network in North America, the sapphire necklace in Europe, etc. [15]. These pathways link regional green spaces such as nature reserves and national parks to maintain the regional ecological balance and conserve wildlife habitats [16]. These concepts do not distinguish between inside and outside of cities, and they include all green resources, emphasizing the natural functions of regional green spaces such as species protection, habitat maintenance, and ecological environmental protection [17]. The concept of regional green space in China was put forward late and is now still in the early stage; in addition, the concept of regional green space varied greatly across different disciplines before the promulgation of the concept of regional green space in the Standard for classification of urban green space (CJJ/T 85-2017). Scholars [18] have discussed the classification of regional green space and the significance of generating it, which has greatly promoted the development of regional green spaces. After the concept of regional green space was clarified, Scoot discussed a control strategy for regional green spaces in the discipline of landscape architecture from the perspective of a park city [19]. Zhou analyzed the connotation and characteristics of regional green space from the categorization adjustment of “other green space” to “regional green space” [20]. Through this analysis, it was found that most of the existing regional green space studies start from the perspectives of classification, connotation, and control strategies, and there are few studies involving dynamically changing regional green space data and the identification of spatial pattern characteristics. Although there have been studies on urban green spaces or green space patterns [21] (although only within the built-up area or the entire region), there are no studies on the dynamically changing regional green space outside the built-up area and its spatial characteristics [22]. The lack of accurate identification techniques and the lack of spatial characteristics of regional green space data are not conducive to the development of regional ecological environments and regional integration [1,23,24].
Thus, this study took the regional green space outside the built-up area of Shanghai as the research object based on three points: (1) Shanghai’s rapid social and economic development in the past three decades attaches great importance to the construction of an ecological civilization, and takes the lead in the preparation of the relevant ecological planning within the country, making outstanding contributions to the protection and development of ecological resources outside the scope of the built-up area of the city. In addition, the development of regional green space outside the urban built-up area is protected and developed accordingly, which has value as a reference for other cities [25]. (2) Shanghai is ranked as the leading city in the whole country in terms of its gross social and economic development value, and its early social and economic development cannot be separated from its urban expansion, which has led to dramatic changes in the ecological resources of the green space outside the built-up area of the city. These dramatic changes in the characteristics of and laws regarding this regional green space are of greater reference value for the regional environment of Shanghai as well as for other developing cities. (3) The sustainable development of Shanghai under the new situation will inevitably put forward higher requirements for the development and management of regional ecological resources. At the regional level, Shanghai is a core city in promoting the integrated development of the Yangtze River Delta and, at the international level, Shanghai is developing into a global city of excellence [26]. Based on development needs at the regional and international levels, the identification of Shanghai’s regional spatial green resources and their characteristics not only can provide important practical support for the reconstruction of Shanghai’s regional ecological environmental protection and optimization, but also plays an important role in the construction of Shanghai as a city of global excellence and the development of the Yangtze River Delta Integration Demonstration Area.
Accordingly, the scientific question of this paper is: how can we scientifically collect dynamic data on the regional green space in Shanghai, which is constantly changing with the expansion of the built-up area? Further, how can we scientifically identify the multi-stage evolution characteristics of the spatial pattern of Shanghai’s regional green space? To address these scientific problems, we used the Google Earth Engine (GEE) platform [27], machine learning [28], biophysical composition (biophysical composition index) [29], CA binarization, and other technical methods to construct a technical method for the identification of regional green spaces. We extracted dynamic data information for the regional green space in Shanghai and adopted a Sankey diagram and spatial pattern index [30] to discriminate the multi-stage evolution characteristics of regional green space pattern in Shanghai. Therefore, this paper reveals the spatial pattern evolution characteristics during the period of 2000–2020 from the three steps of “dynamic data extraction—characterization of spatial pattern evolution—discussion and summary of the problem”, which will provide a reference for the future construction of green space in Shanghai as well as other urban areas.

2. Materials and Methods

2.1. Overview of the Study Area

Shanghai is a province-level administrative region, municipality directly under the central government, national central city, mega city, and core city of the Shanghai metropolitan area of the People’s Republic of China (Figure 1). It has been approved by the State Council as China’s international economic, financial, trade, shipping, and technological innovation center (1). As of 2019, the city has 16 districts under its jurisdiction, with a total area of 6340.5 km2 and a built-up area of 1237.85 km2. At the end of 2021, the permanent population of the city was 24.8943 million. In 2021, the gross domestic product of Shanghai was CNY 4321.485 billion, with a year-on-year increase of 8.1%.
The land area of the Shanghai coastal zone, except for a few exposed volcanic hills in the west and south, is covered by extremely thick Quaternary loose sediment, which thickens from southwest to northeast. The thickness of the Quaternary system on the north bank of Hangzhou Bay is about 100–220 m, the thickness on the eastern coastal plain is 280–300 m, and the thickness in the northern estuarine sand island area is 300–350 m. The coastal zone of Shanghai belongs to a subtropical monsoon climate, exhibiting monsoon and oceanic climate characteristics, with four distinct seasons, sufficient sunlight, and abundant rainfall. The climate in the coastal zone of Shanghai is mild and humid, with a shorter spring and autumn, and longer winter and summer. The annual rainfall is moderate, with around 60% of the year’s rainfall concentrated during the flood season from May to September. The annual average temperature is 15.8 °C, with an average of 3.6 °C in January and 27.8 °C in July.

2.2. Data Sources and Processing

Firstly, by comparing data sources from different reports in the literature, research data sources with a higher accuracy were selected to obtain socio-economic and administrative boundary data.
The administrative boundary data were from the Resource and Environmental Science and Data Center of the Chinese Academy of Sciences. The GDP 1 km × 1 km spatial grids for 2000, 2005, 2010, 2015, and 2020 were generated by Xu Xinliang et al. based on the national GDP statistical data of each county, taking into account the spatial interaction between land use types, nighttime light brightness, residential density data, and GDP closely related to human activities, and through spatial interpolation. The multi-factor weight distribution method was used to segregate the population data using administrative regions as the basic statistical unit on the spatial grid, thus achieving spatialization of the population. The night light data are a type of satellite remote sensing data mainly used for research work such as monitoring urbanization progress. These data are jointly produced by the Digital China Research Institute (Fujian) and the Key Laboratory of Spatial Data Mining and Information Sharing at Fuzhou University, with a resolution of 500 m. Through cross-sensor correction, DMSP-OLS and NPP-VIIRS can be continuously used.

2.3. Subsection

The regional green space is located outside the built-up area, and its land use will change dynamically with the expansion of the built-up area, resulting in changes in coverage. Finding a method to collect dynamic regional green space data is the goal of this research. This article defines the coverage of regional green space as non-construction land outside the built-up area, and defines the composition of regional green space as forest land, grassland, and wetland. With the help of the GEE platform [31], a dynamic recognition technology system for regional green space was constructed through machine learning (random forest, support vector machine, decision tree, and object-oriented classification), biophysical components (BCI), and GIS methods to extract dynamic data information on the scope of built-up areas and regional green spaces at five time nodes in Shanghai (2000, 2005, 2010, 2015, and 2020). Based on the above data, this study analyzed the evolution characteristics and main problems of regional green space spatial pattern identification from three perspectives: area changes, spatial structure changes, and spatial center of gravity changes over time (Figure 2).

2.3.1. Land Use Interpretation

In order to study the spatiotemporal changes in the regional green space in Shanghai, it is necessary to interpret the land cover and the performance of different remote sensing classification algorithms to map land cover based on the Google Earth Engine (GEE) cloud platform. High-spatial-resolution remote sensing images from Landsat and digital elevation data were used. Support vector machine (SVM), classification regression tree (CART), and random forest (RF) algorithms, as well as object-oriented methods were used in the supervised classification of Landsat 5/8 data for Shanghai [32]. In addition, different remote sensing classification algorithms and feature variable combinations were used to compare and analyze the results obtained from the classification process, and the classification results with the highest classification accuracy and the best spatiotemporal continuous change were used as the research data. Remote sensing images were used to extract the urban built-up areas. Currently, the most popular choice for interpreting urban cover at long-term nodes is the Landsat series remote sensing images.
This study utilized the GEE cloud platform and its remote sensing classification algorithms and remote sensing data to classify the land cover situation in the Shanghai area. The research methods used include the collection and optimization of sample data, the construction of a feature space and machine learning classification algorithms, and the evaluation of the classification results. The specific classification methods and results are provided in the supporting materials. The classification results are shown in Figure 3.

2.3.2. Extraction of Urban Built-Up Areas

This study used a combination of the BCI and NDVI (normalized difference vegetation index) to obtain the initial impermeable spatial distribution (Figure 4). Compared with commonly used indices such as the NDBI (normalized difference built-up index), NDVI, and NDSI (normalized difference snow index), the BCI can better distinguish between bare soil, vegetation, and impermeable surfaces, making it widely used for extracting impermeable surfaces. To obtain the urban impermeable surface (Figure 4A), the CA binarization method (Figure 4B) was used to extract the urban built-up area range (Figure 4C) of Shanghai in 2000, 2005, 2010, 2015, and 2020. The specific research steps are provided in the supporting materials.

2.3.3. Measurement of Land Use Quantity and Spatial Change

  • Land use transfer matrix. The land use transfer matrix can reflect the direction and quantity of land use transfer during a certain period of time. We used the reclassification and calculation functions in ArcGIS10.6 to obtain the transfer area of land use functions in each time period [33].
S i j = S 11 S 1 n S n 1 S n n
Here, S is the land area; N is the number of types of land use; and i and j are the land use types at the beginning and end of the research period.
2.
Shifts in land use functions. The transfer of the center of gravity of land use function can reflect the transfer of the center of gravity of various types of land in different years. Using the drawing function of ArcGIS 10.8 software, we obtain the standard deviation ellipses of various types of land use in different periods, and can generate the coordinates of the distribution of the center of gravity of various types of land use.
X = i = 1 n X i W i i = 1 n W i , Y = i = 1 n Y i W i i = 1 n W i
Here, X and Y represent the longitude and latitude coordinates of the distribution of the center of gravity of various types of land use; W i represents the area of the i-th land patch; X i   and   Y i represent the longitude and latitude, respectively, of the geometric center of gravity of the ith land patch; and N represents the total number of patches of various types of land use in that year.

2.3.4. Landscape Pattern Index Screening

Based on the characteristics of Shanghai’s urban development and environmental factors, a reasonable and effective use of the moving window method was to deeply analyze the spatial structure and internal differences in the fragmentation of green spaces outside the built-up areas in Shanghai. The moving windows method can clearly determine the dynamic changes in the spatial structure of regional green spaces, and can provide a more in-depth analysis of the changes in regional green land use and related land use, as well as the resulting urban pattern evolution [33]. The moving window method mainly outputs a new grid map corresponding to the selected landscape index by statistically analyzing the selected features within the window. The specific steps are as follows: (1) pixel verification was conducted on the current land use situation in Shanghai from 2000 to 2020, ensuring that the spatial position and size of pixels in each period were equal and, (2) after proofreading the grid layers from 2000 to 2020, the land use in Shanghai was analyzed using moving windows at the landscape level and type level, and the corresponding landscape indicator spatial grid maps were obtained. The landscape pattern calculation based on the moving window method was implemented using the Fragstats 4.2 platform. This study mainly analyzed the landscape pattern of and changes in green spaces outside the built-up areas in Shanghai. At the same time, in order to analyze the changes in landscape patterns in different regions, this study also calculated the landscape patterns on a 1 km × 1 km scale and used their spatial visualization to analyze the changes in landscape indices in the study area.
A single landscape pattern index cannot easily express the structure of the entire green space system’s landscape pattern, so multiple landscape pattern indices need to be used to jointly reflect the characteristics of the landscape pattern. To comprehensively express the characteristics of a landscape structure, scholars usually divide the landscape pattern index into three levels: the patch level, patch type level, and landscape level. Fargstats 4.2 is the most commonly used software for landscape pattern analysis and the calculation of landscape indices, with up to 277 indicators describing the three levels of landscape patterns in the software. The patch level index expresses the spatial structure characteristics of individual patches, and has little significance in explaining various types of patches, i.e., the entire landscape structure, at the urban scale. Many landscape pattern indices at the other two levels express a similar content and meaning, and many different indices have similarities in calculation and meaning, and their results also have a certain correlation. Therefore, it is necessary to choose the most suitable simple and persuasive landscape pattern index to illustrate the problem, in order to quantitatively describe the changes in landscape structural characteristics over different periods [34].
The main principles for selecting landscape pattern indices in this study were as follows: (1) the calculation formula for landscape pattern indices is relatively simple and will not be affected by changes in the other indices; (2) the ecological or spatial significance of the index is clear, which is conducive to describing the characteristics of the spatial structure; and (3) the selected index can comprehensively reflect the heterogeneity or spatial structure characteristics of the landscape from various perspectives, which means that suitable landscape pattern indices can be selected from multiple perspectives to comprehensively and systematically describe the structure, characteristics, and dynamics of the landscape pattern. This study ultimately selected indicators that can describe the patch type and landscape, including the patch type area (CA), fragmentation index (PD), landscape shape index (LSI), connectivity index (cohesion), and Shannon diversity index (SHDI). The calculation formula for the landscape index is shown in Table 1.

3. Results

3.1. Dynamic Degree of Green Space Area in Shanghai

3.1.1. Change in Green Area

The land use type results based on spatial remote sensing images of Shanghai’s regional green spaces in the five periods of 2000, 2005, 2010, 2015, and 2020 were used to produce a land use data statistical table (Table 2) and Sangji map through ArcGIS 10.6 software. From the temporal changes in regional green space in Shanghai, a growth trend of “first decreasing and then increasing” can be seen in the green space outside the built-up areas of Shanghai during the research period. From 2000 to 2005, there was overall decrease in regional green space, with the area of regional green space decreasing from 580.556 km2 to 506.426 km2. However, due to urban expansion, the proportion of regional green space in the research area increased from 7.13% to 8.77%. Subsequently, from 2005 to 2020, the total area increased rapidly from 506.426 km2 to 880.99 km2, with the growth rate reaching its maximum in 2010 and slowing down after 2010. The trend of the increases and decreases in arable land, water system, and urban built-up areas is consistent with previous research.
Between 2000 and 2020, the large-scale shrinkage of arable land became the main source for the increase in regional green space area, with forest land continuing to increase and wetlands and grasslands fluctuating, resulting in an overall increase in area (Figure 5). From the perspective of the spatiotemporal changes in the conversion of green spaces outside the built-up areas in Shanghai, the land-type changes in the regional green spaces in Shanghai mainly involved the conversion of arable land to construction land and forest land, and the conversion of forest land and grasslands to construction land. Overall, there was a conversion of regional green spaces to construction land. In the early stage, regional green space was mainly formed naturally. During this period, the expansion of urban built-up areas resulted in the export of a large amount of arable land, and the conversion of regional green space to other land uses was not significant in terms of the total amount. With the development of the economy and society, construction projects such as land reclamation, as well as various ecological security and protection projects, led to two spatial transformation trends: firstly, human spontaneous protection measures led to a large number of water systems being converted into construction land and wetlands, making the spatial pattern of the regional green spaces more stable. The period from 2000 to 2015 was one of intense changes, and, from 2015 to 2020, the overall transformation of regional green spaces and other types of land became stable.

3.1.2. Dynamic Changes in Regional Green Space and Related Land Use

Overall, according to Table 3, during the land use change in the study area from 2000 to 2020, the cultivated land area showed a downward trend, with an initial rapidly declining rate and then a gradual decline. From 2000 to 2005, the overall area of the regional green space decreased. From 2005 to 2020, the overall area of regional green space showed an upward trend, but the upward trend gradually slowed down.
From 2000 to 2005, due to the expansion of urban construction land, the area of various types of land decreased on a large scale. The cultivated land area decreased by 186.226 km2, the grassland area decreased by 85.088 km2, the water area decreased by 62.702 km2, and the wetland area decreased by 70.548 km2. The changes in the six types of land use were −4.4%, 37.8%, −47.9%, −3.7%, 46.4%, and −36.7%, respectively.
From 2005 to 2010, due to the expansion of urban construction land, the area of various types of land decreased, which was the largest change in land use types during the entire study period. The cultivated land area decreased by 186.226 km2, the grassland area increased by 85.088 km2, the water area decreased by 62.702 km2, and the wetland area increased by 21.775 km2. The changes in the six types of land use were 13.1%, 49.3%, −47.9%, −70.3%, −3.4%, and 21.77%, respectively.
From 2010 to 2015, due to the expansion of urban construction land, the area of cultivated land continued to decrease, with an increase of 260.39 km2. The grassland area decreased by 85.088 km2, the water area increased by 62.702 km2, and the wetland area increased by 70.548 km2. The changes in the six types of land use were −4.4%, 37.8%, −47.9%, −3.7%, 46.4%, and −36.7%, respectively.
From 2015 to 2020, due to the expansion of urban construction land, there were significant differences in the areas of various types of land, with an increase in the areas of forest land, grasslands, and wetlands. The cultivated land area still decreased, with a decrease of 186.226 km2; meanwhile, there was an increase of 3.706 km2 in grassland area, a decrease of 118.607 km2 in water area, and an increase of 21.389 km2 in wetland area. The changes in the six types of land use were −4.4%, 37.8%, −47.9%, −8.2%, 9.9%, and 8.9%, respectively.

3.2. Spatial Center of Gravity of Regional Green Spaces

3.2.1. LUCC Center of Gravity and Ellipse Direction

The standard deviation ellipse and center of gravity of different land use functions in each period were plotted to derive the range and directionality of the main distribution of different land use/cover in Shanghai, as well as the direction of transfer and expansion. In the generated ellipses, the long axis indicates the main distribution direction of the type of land use, and the short axis indicates the main distribution range. A smaller short axis indicates a more obvious centripetal force, and a flatter ellipse indicates a stronger directionality.
From Figure 6A, it can be seen that the patterns of construction land, forest land, grasslands, and wetlands in the overall land area of Shanghai changed significantly from 2000 to 2020, while the spatial pattern of the water area changed less. As the overall spatial center of gravity of urban development and construction increased centripetally, the length of the short axis of the ellipse became smaller, but because of the urban expansion in the “big pie” style, the overall directionality of the construction space was weakened, and the direction of the overall center of gravity gradually moved to the west. The cultivated land space gradually shifted from east to west, with increasing directionality, and the pattern of cultivated land basically remained stable after 2015. The center of gravity of the wetlands moved in the northwest direction, which was mainly due to the release of the “Shanghai Wetland Protection and Restoration Plan” after 2006, which restored the original damaged wetlands to the status quo. In addition, the reclamation project of Shanghai Municipality made the area of the Chongming wetland and other wetlands increase dramatically, which caused the center of gravity of the wetlands to shift in the northwest direction on a large scale. From Figure 6B, it can be seen that the overall movement direction of the regional green space in the study area was basically the same as that of the whole area of Shanghai, because the green space within the built-up area of the city is mainly an artificial landscape, and generally speaking, there will not be a wide range of changes over a long period of urban development.

3.2.2. Regional Green Space Center of Gravity Shift

From Figure 7, it can be seen that the transfer direction of the center of gravity of the regional green space was highly consistent with the development direction of green space in Shanghai, and the spatial locations basically coincided, indicating that changes in the regional green space were the main type of change in the evolution of the spatial pattern of green space in Shanghai. According to the analysis of the center of gravity transfer distance and orientation of the woodlands, grasslands, and wetlands, except for the center of gravity of the wetlands which basically remained unchanged during this 20-year period, the transfer direction of the grasslands mainly shifted in the southwestern direction during 2000–2005, and shifted in the southeastern direction during the period of 2005–2010 for a short distance. Between 2010 and 2015, the most drastic change in the grasslands occurred, with the center of gravity of the grasslands shifting to the northwest; between 2015 and 2020, the center of gravity of the grasslands shifted to the east; and during the 20-year period, the overall center of gravity of the grasslands shifted to the north. The center of gravity of the woodlands shifted to the north in 2000–2005, and to the south-east during 2005–2010, with small changes in the distance shifted; the center of gravity of the woodlands tended to stabilize during the 10-year period from 2010 to 2020, and there was relatively little fluctuation in the distance shifted.

3.3. Evolution of Regional Green Space Landscape Pattern

3.3.1. Spatial and Temporal Evolution of Patch Area

The area of a patch type mainly indicates the total area of each type of land and regional green space as a whole. From Table 4 and Figure 8, it can be seen that the grassland area firstly decreased and then increased, and between 2000 and 2010, the area gradually decreased by nearly half; between 2010 and 2020, the area increased again, and the overall area increased relative to 2000. The area of forested land increased dramatically during the 20-year period, from 215.59 km2 in 2000 to 454.14 km2 in 2020, and the wetland area first increased and then decreased, from 192.31 km2 to 143.53 km2 during 2000–2010, then began to increase from 143.53 km2 to 261.90 km2 after 2010. Overall, the grassland and wetland areas were closer in size at the beginning of the study in 2000, and then the gap began to widen considerably. The areas of the two types of regional greenspace became gradually closer after 2010, but the gap between the areas had increased relative to that at the beginning of the study. The woodlands always led the way, with the highest proportion of the total regional green space.

3.3.2. Degree of Regional Green Space Fragmentation

The patch fragmentation index mainly indicates the size of the number of patches in each type of regional green space and the regional green space as a whole, which can reflect the degree of fragmentation and completeness of a certain type of regional green space or the regional green space as a whole. From Table 5 and Figure 9, it can be seen that the fragmentation degree is increasing in the region as a whole, which is the result of the development of urban construction and is similar to the findings of many scholars. The patch fragmentation of woodland increased from 0.82 to 1.91, that of wetland increased from 0.02 to 0.04, and the patch fragmentation (PD) of grassland decreased from 1.17 to 1.01 during 2000–2020 in the regional green space of Shanghai, whereas the increase in the patch area of grassland indicates that the distribution of grassland is still tending to be intact and the degree of fragmentation is decreasing on the whole.

3.3.3. Regional Green Space Connectivity

The connectivity index, also known as the aggregation index, mainly indicates the degree of aggregation of each type of regional green space and the overall patches of regional green space, with a larger value indicating a better spatial connectivity. As shown in Table 6 and Figure 10, the connectivity index of the wetlands in the study area showed a continuous decreasing trend during the study period, while the connectivity index of the woodlands showed a continuous increasing trend, and the connectivity index of the grasslands fluctuated. The overall connectivity index of the green space in the study area also showed a fluctuating decrease from 94.76 to 94.65. Combined with the increase in the area of the woodlands, grasslands, and wetlands, the increase in the woodlands was a piecemeal increase; for example, the protection forest on both sides of the watershed is in the form of a band, which will make the connectivity of the woodland in the area increase significantly. The connectivity of grassland slightly decreased, but the area increased. This was related to the distribution of grasslands in the form of a point distribution.

3.3.4. Shannon Diversity of Regional Green Space

The Shannon diversity index (SHDI), which mainly indicates the diversity of regional green space types, is expressed by the uniformity of the distribution of each type of regional green space patch and the overall patches, and a larger value indicates that types of patch are more diverse, and that they are more uniformly or evenly distributed in the overall landscape. As can be seen from Table 7, the SHDI of regional green spaces in the study area increased, indicating that the diversity of regional green spaces in the study area gradually increased. However, spatially (Figure 11), with urban expansion, the SHDI around it decreased substantially; over the past 20 years, the SHDI along the Jinshan–Fengxian–Pudong line decreased and then increased, and the SHDI near Dianshan Lake continued to increase.

4. Discussion

The regional green space outside the built-up area of Shanghai, which is the most seriously disturbed area in the process of rapid urbanization [35], is important in the study of its landscape spatial pattern characteristics. Previous studies on the spatial and temporal evolution of the green space landscape pattern in Shanghai were mostly based on the built-up area or the overall green space identified using remote sensing image data [36]. The green space data cover the built-up area inside and outside of the city, and it is difficult to determine the unique characteristics of the changes in the spatial pattern of the green space landscape in the areas outside the built-up area. For example, Zhang Li-Quan et al. quantitatively revealed the spatial distribution characteristics of the landscape pattern in Shanghai. With the help of GIS spatial analysis technology, Zhang et al. [37] studied the typical regional green corridor space in Shanghai, constructed an evaluation system, and provided a scientific and effective basis for the development of regional green corridors. The overall lack of spatial characterization of regional green landscapes outside the built-up area prevents regional ecological environment construction [38]. Liu et al. made a relevant analysis of the overall green space pattern changes in Shanghai from 1990 to 2015 in terms of the degree of change and landscape indicators, mainly exploring the impact of urban expansion on arable land, and did not provide an in-depth discussion of the phase of regional green space [39]. Based on the kilometer grid method, Da et al. analyzed the change patterns of the green space landscape in the central urban area of Shanghai from 2005 to 2015, and quantified the obvious “spatial characteristics” of different types of green space along the gradient of landscape patterns; however, the study was based on the forest soil and did not investigate the regional green space. The rapid urbanization of Shanghai since the turn of the century has led to the rapid development of the city. Since the 21st century, Shanghai has experienced a period of rapid urbanization, with the regional green space landscape spatial pattern change rate being significantly accelerated as a result. The regional environmental coverage of ecological resources has become much larger than the built-up area [1]. The existing research is mostly focused on the scope of the built-up area, but the green space within the built-up area and the regional green space outside the built-up area are very different [40,41]. The collection of regional green space dynamic data and their spatial evolution characteristics need to be further strengthened to fill the gap in the relevant research. This study analyzed the spatial pattern of the regional green space landscape outside the built-up area of Shanghai using remote sensing data from 2000 to 2020, effectively making up for the lack of research on regional green spaces outside built-up areas, to help ensure the healthy and stable development of the regional environmental system.
The relationships between the area growth and landscape pattern indices (fragmentation, connectivity, diversity) at different stages did not necessarily show positive correlations, but there was a positive correlation between the area reduction and landscape pattern indices. For example, Zhang et al. [42] found that the area of green space in Jinan City continued to increase during 1996–2004, but the diversity did not change much, while Deng Jinsong et al. [43] showed that the total area of green space in Hangzhou City continued to decrease slowly during 1996–2004, and that the fragmentation was still high.
The results of this study show that, between 2000 and 2020, the total area of green space in the Shanghai Municipality showed a change characteristic of “decreasing first and then increasing”. For example, in the four phases of 2000–2005, 2005–2010, 2010–2015, and 2015–2020, the changes in the area of regional green space were a rapid decrease, slow increase, rapid increase, and then slow increase. The trend of the regional green space growth was closely related to the socio-economic development of Shanghai. Socio-economic development leads to urban expansion, and also promotes investments in landscaping and greening. In this stage, from 2001 to 2005, the regional green space area decreased by 74.1 km2. However, the cognition of ecological and environmental protection tends to lag behind socio-economic development [44], and there was a spatial and temporal divergence of green space and fragmentation of the urban area of Guangzhou in the last 35 years. Real estate investment and the number of factories have a negative effect on regional green space. The amount of regional green space occupied by urban expansion is much larger than the amount of new construction, and there is a “game” relationship between regional green spaces and socio-economic development, during which the fragmentation, connectivity, and diversity indices are negatively affected. In 2005–2010, socio-economic development increased to a certain degree, and began to focus on green development, resulting in the area of regional green space increasing by 106.7 km2. To a certain extent, the government began to pay attention to green development, and the government began to increase the construction of regional green spaces outside built-up areas to curb urban sprawl. The regional green space area began to show good growth, which led to an increase in the connectivity and biodiversity in the landscape pattern index, but there was still severe fragmentation, and the phase of synchronous growth of the regional green space and the expansion of the city was realized. In 2010–2015, the regional green area increased by 210.1 km2, and the regional green area in this phase was more fragmented than in the previous one. The regional green space growth in this stage was more obvious than in the previous stage, as the positive role of government policy continued to strengthen, which leaned toward ecological environmental protection and green development by increasing the investment in ecological environment construction and greatly promoting the construction of regional green spaces. From 2015 to 2020, urban expansion was basically stable, with none of the rapid expansion which was seen in the previous ten years, which led to a stage where the regional green space area increased by only 82.4 km2, only half of the growth seen in the previous stage. This stage of urban development was the era where the regional green space slowly increased, and the regional green space and region entered the stage of high-quality development.
The rapid development of regional integration in Shanghai and the construction of regional facilities accelerated the fragmentation of natural ecological resources such as landscapes and forests in the region, resulting in a reduction in regional ecosystem stability and deterioration of connectivity. By accurately analyzing the regional green space data, we can scientifically discern the characteristics of the regional green space pattern evolution and identify the main problems. Under this new situation, the integrated development of Shanghai’s regional environment will inevitably put forward higher requirements for the development and management of the regional green space, and targeted research on the evolution of its spatial pattern can not only provide important theoretical support for Shanghai’s regional ecological environmental protection work and optimization of environmental reconstruction, but also has great significance for the city of Shanghai’s goals of global excellence and serving as the model for the integration of the Yangtze River Delta [23]. From a governing perspective, these research results can provide a decision-making basis for regional planning, ecological spatial planning, regional integration management, etc., which can effectively support national land spatial planning [22], enhance the stability of the regional ecology, and provide urban and rural residents with open spaces for recreation [45]. From a market perspective, the results of this research can provide a direct reference for the planning, construction, and management of regional green spaces, and promote the comprehensive function of regional green spaces’ inherent ecological resources [34]. From a social perspective, the results of this research can help to curb regional ecological environment fragmentation, promote the integration of gray infrastructure and the blue-green ecological resources in the region, enhance the image of the regional environment, and strengthen the quality of life for residents.
However, due to the limitation of space, there are some shortcomings in this study, which only focused on the evolutionary characteristics of green space patterns in Shanghai, a region with a higher level of socio-economic development, and did not study these characteristics in regions with a lower level of socio-economic development. A comparative study of the differences in the characteristics of regional green space evolution between Shanghai and underdeveloped cities, and further investigation into these characteristics, will play an important role in the decision-making regarding regional green space planning and the construction practices of cities.

5. Conclusions

Based on the moving window method, landscape pattern indices, transfer matrices, and other methods, this study comprehensively analyzed the evolution process and characteristics of regional green space patterns outside the built-up area of Shanghai between 2000 and 2020 from the perspectives of changes in regional green space, changes in spatial structure, and changes in the center of gravity of development. The specific conclusions are as follows:
(1)
Constructing technical methods for acquiring scientific regional green space data.
Regional green space is located outside the built-up area, which not only has a large number of patches, large area, and wide coverage, but also involves different management departments such as water conservancy, transportation, agriculture, environmental protection, etc., and its scope and land use change dynamically with the expansion of the built-up area. A method to accurately analyze the dynamically changing regional green space data is needed in practice and research. In this study, through machine learning, biophysical components (BCI), CA binarization, and other technical methods, we constructed a regional green space identification technology system to identify regional green space dynamic data. This method was used to obtain a more comprehensive and relatively reasonable urban built-up area scope as well as regional green space data, and also provides a reference for the development of regional green space in Shanghai or other cities, as well as regional environment-related research.
(2)
The evolution of Shanghai’s regional green area shows a fluctuating growth trend of “first decline and then rise”.
The cultivated land area in Shanghai showed a large area of reduction, and the water area continuously decreased, but the regional green space in terms of wetlands and grasslands showed a development trend of “first decline then rise”. In addition, the forested land area maintained continuous and steady growth. The changes in the land type of the regional green space in Shanghai were divided into several stages: the initial stage of regional green space was mainly formed naturally; the expansion of urban built-up areas eroded a large amount of arable land during this period, when the conversion of regional green space into other land types was not obvious. With economic and social development, reclamation and other construction projects, as well as various types of ecological security and ecological protection projects, caused the regional green space to change dramatically; with the socio-economic level of development at a higher level, the overall area of the regional green space and other types of land tended to stabilize. Therefore, the development of similar urban environments can learn from the Shanghai Municipality in terms of the construction of a large number of woodlands to promote the regional ecological environment, but also to avoid eroding the ecological environment after the destruction process has begun.
(3)
The number of regional green spaces in Shanghai increased, but the spatial pattern at the overall level was fragmented and heterogeneous; at the local level, the woodlands, and wetlands in some areas showed continuous development, with the development trend of point-line synergy.
Large patches in the regional green space were mostly scenic recreational green spaces and ecological conservation green spaces, and the number of individual patches was large and relatively dispersed, such as the large patches of regional green space in the Sheshan area and Chongming Island area due to the concentration of original natural resources, which are relatively isolated and dispersed. The number of small patches was large, and they were aggregated; they were ecological conservation green spaces and regional facilities protection green spaces, and were distributed in water systems or around regional facilities in a strip-like aggregated distribution, with a strong spatial connectivity value. It can be seen that the coordination and integration between the regional green space was poor, presenting the spatial layout characteristics of a small number of extra-large patches that were dispersed, and a large number of small patches that were aggregated. Therefore, in the future planning and construction of regional green spaces in Shanghai or other cities, it is not only necessary to strengthen the protection of large regional green spaces, but also to consider the interconnectivity between them to avoid fragmentation.
(4)
The spatial center of gravity of Shanghai’s regional green space has a high degree of consistency with the overall green space changes.
The difficulty of increasing the urban green space area leads to changes in the regional green space outside the built-up area of the city, which directly affects the change in the overall green space pattern of Shanghai. According to our analysis of the center of gravity transfer distance and orientation of the woodlands, grasslands, and wetlands in Shanghai’s regional green space, the center of gravity of the wetlands basically remained unchanged during the 20-year period, while that of the grasslands moved significantly to the northwest and that of the woodlands was transferred by a smaller distance and tended to be stabilized. Therefore, in the future, regional ecological environment protection and construction in Shanghai or other cities can not only focus on the green space resources within the urban built-up area, but also need to strengthen the regional green space resources outside of the built-up area.

Supplementary Materials

The supporting information can be downloaded at https://www.mdpi.com/article/10.3390/land12122184/s1.

Author Contributions

Conceptualization, Y.J. and L.Z. (Lang Zhang); methodology, Y.J.; software, X.G.; validation, X.G.; formal analysis, L.Z. (Lei Zhang); investigation, L.Z. (Lei Zhang); resources, L.Z. (Lei Zhang); data curation, X.G.; writing—original draft preparation, Y.J.; writing—review and editing, L.Z. (Lang Zhang); visualization, Y.J.; supervision, Y.J.; project administration, L.Z. (Lang Zhang); funding acquisition, L.Z. (Lang Zhang). All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Doctoral Research Initiation Fund Project of Anhui Jianzhu University (grant number: 2023QDR03), the National Natural Science Foundation of China (grant number: 32171569), the Anhui Cultural Tourism Innovative Development Research Institute (grant number: ACTZ2022YB01), and the Open Fund of Sichuan County Economic Development Research Center (grant number: XY2023030).

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lai, X.; Yu, H.; Liu, G.; Zhang, X.; Feng, Y.; Ji, Y.; Zhao, Q.; Jiang, J.; Gu, X. Construction and Analysis of Ecological Security Patterns in the Southern Anhui Region of China from a Circuit Theory Perspective. Remote Sens. 2023, 15, 1385. [Google Scholar] [CrossRef]
  2. Liu, K.; Xue, Y.; Chen, Z.; Miao, Y. The Spatiotemporal Evolution and Influencing Factors of Urban Green Innovation in China. Sci. Total Environ. 2023, 857, 159426. [Google Scholar] [CrossRef] [PubMed]
  3. Vigneshwaran, S.; Vasantha Kumar, S. Comparison of Classification Methods for Urban Green Space Extraction Using Very High Resolution Worldview-3 Imagery. Geocarto Int. 2021, 36, 1429–1442. [Google Scholar] [CrossRef]
  4. Mohan, E.; Rajendran, K. Effect of Plant Growth-Promoting Microorganisms on Quality Seedling Production of Feronia elephantum (Corr.) in Semi-Arid Region of Southern India. Int. J. Curr. Microbiol. Appl. Sci. 2014, 3, 103–116. [Google Scholar]
  5. Wang, Y.; Chang, Q.; Fan, P.; Shi, X. From Urban Greenspace to Health Behaviors: An Ecosystem Services-Mediated Perspective. Environ. Res. 2022, 213, 113664. [Google Scholar] [CrossRef] [PubMed]
  6. Gougeon, S.; Kemp, K.M.; Blicher, M.E.; Yesson, C. Mapping and Classifying the Seabed of the West Greenland Continental Shelf. Estuar. Coast. Shelf Sci. 2017, 187, 231–240. [Google Scholar] [CrossRef]
  7. Erickson, D.L. Connecting Corridors: Implementing Metropolitan Greenway Networks in North America. In Ecological Networks and Greenways: Concept, Design, Implementation; Cambridge University Press: Cambridge, UK, 2004; pp. 200–221. [Google Scholar]
  8. Biondi, E.; Casavecchia, S.; Pesaresi, S.; Zivkovic, L. Natura 2000 and the Pan-European Ecological Network: A New Methodology for Data Integration. Biodivers. Conserv. 2012, 21, 1741–1754. [Google Scholar] [CrossRef]
  9. Briquel, V.; Collicard, J.-J. Diversity in the Rural Hinterlands of European Cities. In The City’s Hinterland; Routledge: London, UK, 2016; pp. 19–40. [Google Scholar]
  10. Evans, A.W. The Myths and Reality of Urban Constraint in United Kingdom: Changing Circumstances and Unchanged Policies; European Real Estate Society (ERES): Amsterdam, The Netherlands, 2013. [Google Scholar]
  11. Wolch, J.R.; Byrne, J.; Newell, J.P. Urban Green Space, Public Health, and Environmental Justice: The Challenge of Making Cities “Just Green Enough”. Landsc. Urban Plan. 2014, 125, 234–244. [Google Scholar] [CrossRef]
  12. ElDeeb, H. Sustainable Development and Green Growth: Study of Lafayette Central Park Proposed Master Plan, USA. Procedia Eng. 2017, 181, 44–51. [Google Scholar] [CrossRef]
  13. Vohora, V.K.; Donoghue, S.L. Application of Remote Sensing Data to Landslide Mapping in Hong Kong. Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci. 2004. Available online: https://www.researchgate.net/profile/Vijay-Vohora/publication/228606641_Application_of_remote_sensing_data_to_landslide_mapping_in_Hong_Kong/links/0deec52babd4c25a4e000000/Application-of-remote-sensing-data-to-landslide-mapping-in-Hong-Kong.pdf (accessed on 1 October 2023).
  14. Thomas, K.; Littlewood, S. From Green Belts to Green Infrastructure? The Evolution of a New Concept in the Emerging Soft Governance of Spatial Strategies. Plan. Pract. Res. 2010, 25, 203–222. [Google Scholar] [CrossRef]
  15. Zhang, Y.B.; Wu, R.W. The Theory and Practice of Greenway Construction in Europe. Chin. Landsc. Archit. 2007, 8, 33–38. [Google Scholar]
  16. Yang, R.; Ren, F.; Xu, W.; Ma, X.; Zhang, H.; He, W. China’s Ecosystem Service Value in 1992–2018: Pattern and Anthropogenic Driving Factors Detection Using Bayesian Spatiotemporal Hierarchy Model. J. Environ. Manag. 2022, 302, 114089. [Google Scholar] [CrossRef] [PubMed]
  17. Yang, G.; Ge, Y.; Xue, H.; Yang, W.; Shi, Y.; Peng, C.; Du, Y.; Fan, X.; Ren, Y.; Chang, J. Using Ecosystem Service Bundles to Detect Trade-Offs and Synergies across Urban–Rural Complexes. Landsc. Urban Plan. 2015, 136, 110–121. [Google Scholar] [CrossRef]
  18. Zhao, W.; Wang, Y.; Chen, D.; Wang, L.; Tang, X. Exploring the Influencing Factors of the Recreational Utilization and Evaluation of Urban Ecological Protection Green Belts for Urban Renewal: A Case Study in Shanghai. Int. J. Environ. Res. Public Health 2021, 18, 10244. [Google Scholar] [CrossRef] [PubMed]
  19. Scott, A.J.; Carter, C.; Reed, M.R.; Larkham, P.; Adams, D.; Morton, N.; Waters, R.; Collier, D.; Crean, C.; Curzon, R. Disintegrated Development at the Rural–Urban Fringe: Re-Connecting Spatial Planning Theory and Practice. Prog. Plan. 2013, 83, 1–52. [Google Scholar] [CrossRef]
  20. Zhou, S.; Song, Y.; Li, Y.; Wang, J.; Zhang, L. Construction of Ecological Security Pattern for Plateau Lake Based on MSPA–MCR Model: A Case Study of Dianchi Lake Area. Sustainability 2022, 14, 14532. [Google Scholar] [CrossRef]
  21. Shen, J.; Li, S.; Wang, H.; Wu, S.; Liang, Z.; Zhang, Y.; Wei, F.; Li, S.; Ma, L.; Wang, Y. Understanding the Spatial Relationships and Drivers of Ecosystem Service Supply-Demand Mismatches towards Spatially-Targeted Management of Social-Ecological System. J. Clean. Prod. 2023, 406, 136882. [Google Scholar] [CrossRef]
  22. Lyu, R.; Zhao, W.; Pang, J.; Tian, X.; Zhang, J.; Wang, N. Towards a Sustainable Nature Reserve Management: Using Bayesian Network to Quantify the Threat of Disturbance to Ecosystem Services. Ecosyst. Serv. 2022, 58, 101483. [Google Scholar] [CrossRef]
  23. Wang, Y.; Gu, X.; Yu, H. Spatiotemporal Variation in the Yangtze River Delta Urban Agglomeration from 1980 to 2020 and Future Trends in Ecosystem Services. Land 2023, 12, 929. [Google Scholar] [CrossRef]
  24. Wang, Y.; Wang, Y.; Xia, T.; Li, Y.; Li, Z. Land-Use Function Evolution and Eco-Environmental Effects in the Tarim River Basin from the Perspective of Production–Living–Ecological Space. Front. Environ. Sci. 2022, 10, 1004274. [Google Scholar] [CrossRef]
  25. Liu, S.; Wu, P. The Impact of Urban Sprawl on Green Total Factor Productivity: A Spatial Econometric Analysis in China. Front. Environ. Sci. 2023, 11, 1095349. [Google Scholar] [CrossRef]
  26. Li, N.; Song, Y.; Xia, W.; Fu, S.-N. Regional Transportation Integration and High-Quality Economic Development, Coupling Coordination Analysis, in the Yangtze River Delta, China. Systems 2023, 11, 279. [Google Scholar] [CrossRef]
  27. Yunfeng, H.; Lingjie, S.; Qianli, Z.; Zhaohai, W. Land Change Patterns and Driving Mechanism in Beijing since 1990 Based on GEE Platform. Remote Sens. Technol. Appl. 2018, 33, 573–583. [Google Scholar]
  28. Rui, Q.; Cheng, H. Quantifying the Spatial Quality of Urban Streets with Open Street View Images: A Case Study of the Main Urban Area of Fuzhou. Ecol. Indic. 2023, 156, 111204. [Google Scholar] [CrossRef]
  29. Yang, F.; Pan, S.; Yao, X. Regional Convergence and Sustainable Development in China. Sustainability 2016, 8, 121. [Google Scholar] [CrossRef]
  30. Chang, Y.; Zhang, G.; Zhang, T.; Xie, Z.; Wang, J. Vegetation Dynamics and Their Response to the Urbanization of the Beijing–Tianjin–Hebei Region, China. Sustainability 2020, 12, 8550. [Google Scholar] [CrossRef]
  31. Zhang, F.; Gao, Y. Composite Extraction Index to Enhance Impervious Surface Information in Remotely Sensed Imagery. Egypt. J. Remote Sens. Space Sci. 2023, 26, 141–150. [Google Scholar] [CrossRef]
  32. Chan, F.K.S.; Griffiths, J.A.; Higgitt, D.; Xu, S.; Zhu, F.; Tang, Y.-T.; Xu, Y.; Thorne, C.R. “Sponge City” in China—A Breakthrough of Planning and Flood Risk Management in the Urban Context. Land Use Policy 2018, 76, 772–778. [Google Scholar] [CrossRef]
  33. Hong, W.; Guo, R.; Su, M.; Tang, H.; Chen, L.; Hu, W. Sensitivity Evaluation and Land-Use Control of Urban Ecological Corridors: A Case Study of Shenzhen, China. Land Use Policy 2017, 62, 316–325. [Google Scholar] [CrossRef]
  34. Fan, X.; Gu, X.; Yu, H.; Long, A.; Tiando, D.S.; Ou, S.; Li, J.; Rong, Y.; Tang, G.; Zheng, Y.; et al. The Spatial and Temporal Evolution and Drivers of Habitat Quality in the Hung River Valley. Land 2021, 10, 1369. [Google Scholar] [CrossRef]
  35. Wang, Z.; Shi, P.; Shi, J.; Zhang, X.; Yao, L. Research on Land Use Pattern and Ecological Risk of Lanzhou–Xining Urban Agglomeration from the Perspective of Terrain Gradient. Land 2023, 12, 996. [Google Scholar] [CrossRef]
  36. Zhou, Z.C.; Wang, J. Evolution of Landscape Dynamics in the Yangtze River Delta from 2000 to 2020. J. Water Clim. Chang. 2022, 13, 1241–1256. [Google Scholar] [CrossRef]
  37. Zhang, B.; Shao, D.; Zhang, Z. Spatio-Temporal Evolution Dynamic, Effect and Governance Policy of Construction Land Use in Urban Agglomeration: Case Study of Yangtze River Delta, China. Sustainability 2022, 14, 6204. [Google Scholar] [CrossRef]
  38. Pan, N.; Guan, Q.; Wang, Q.; Sun, Y.; Li, H.; Ma, Y. Spatial Differentiation and Driving Mechanisms in Ecosystem Service Value of Arid Region: A Case Study in the Middle and Lower Reaches of Shule River Basin, NW China. J. Clean. Prod. 2021, 319, 128718. [Google Scholar] [CrossRef]
  39. Liu, S.; Zhang, X.; Feng, Y.; Xie, H.; Jiang, L.; Lei, Z. Spatiotemporal Dynamics of Urban Green Space Influenced by Rapid Urbanization and Land Use Policies in Shanghai. Forests 2021, 12, 476. [Google Scholar] [CrossRef]
  40. Wu, J.; Fan, X.; Li, K.; Wu, Y. Assessment of Ecosystem Service Flow and Optimization of Spatial Pattern of Supply and Demand Matching in Pearl River Delta, China. Ecol. Indic. 2023, 153, 110452. [Google Scholar] [CrossRef]
  41. Gu, X.; Long, A.; Liu, G.; Yu, J.; Wang, H.; Yang, Y.; Zhang, P. Changes in Ecosystem Service Value in the 1 Km Lakeshore Zone of Poyang Lake from 1980 to 2020. Land 2021, 10, 951. [Google Scholar] [CrossRef]
  42. Zhang, L.; Kong, F.; Sun, Z.; Zhuang, Y.; Yin, H. Evaluation of the Land-Use Suitability of Jinan City Based on GIS and Landscape Metrics. In Proceedings of the 2010 18th International Conference on Geoinformatics, IEEE, Beijing, China, 18–20 June 2010; pp. 1–6. [Google Scholar]
  43. Deng, J.-S.; Li, J.; Yu, L.; Wang, K. Dynamics of Land Use Landscape Pattern in Hangzhou City during Its Rapid Urbanization. Yingyong Shengtai Xuebao 2008, 19, 2003–2008. [Google Scholar]
  44. Zhao, X.; Zhang, M.; Li, Y.; Huang, X.; Wang, B.; Zhang, L. Urban Residential Land Expansion and Agglomeration in China: A Spatial Analysis Approach. Environ. Dev. Sustain. 2020, 22, 5317–5335. [Google Scholar] [CrossRef]
  45. Sun, X.; Yang, P.; Tao, Y.; Bian, H. Improving Ecosystem Services Supply Provides Insights for Sustainable Landscape Planning: A Case Study in Beijing, China. Sci. Total Environ. 2022, 802, 149849. [Google Scholar] [CrossRef]
Figure 1. Overview of the research area. (a) Elevation and administrative divisions of Shanghai; (b) geographical location of Shanghai.
Figure 1. Overview of the research area. (a) Elevation and administrative divisions of Shanghai; (b) geographical location of Shanghai.
Land 12 02184 g001
Figure 2. Research procedure.
Figure 2. Research procedure.
Land 12 02184 g002
Figure 3. Land use results.
Figure 3. Land use results.
Land 12 02184 g003
Figure 4. Method for extracting built-up areas in Shanghai. (A) Results after extraction; (B) CA binarization diagram; (C) CA binarization results; (D) biophysical composition (BCI) calculation process.
Figure 4. Method for extracting built-up areas in Shanghai. (A) Results after extraction; (B) CA binarization diagram; (C) CA binarization results; (D) biophysical composition (BCI) calculation process.
Land 12 02184 g004
Figure 5. Sangji map of Shanghai land use changes.
Figure 5. Sangji map of Shanghai land use changes.
Land 12 02184 g005
Figure 6. Spatial distribution map of land use centers and ellipses of different land use types in Shanghai from 2000 to 2020. (A) Shanghai as a whole. (B) Outside urban built-up area (study area).
Figure 6. Spatial distribution map of land use centers and ellipses of different land use types in Shanghai from 2000 to 2020. (A) Shanghai as a whole. (B) Outside urban built-up area (study area).
Land 12 02184 g006
Figure 7. Spatial distribution map of land use centers and ellipses in Shanghai from 2000 to 2020.
Figure 7. Spatial distribution map of land use centers and ellipses in Shanghai from 2000 to 2020.
Land 12 02184 g007
Figure 8. Spatial distribution of green space fragmentation in Shanghai from 2000 to 2020.
Figure 8. Spatial distribution of green space fragmentation in Shanghai from 2000 to 2020.
Land 12 02184 g008
Figure 9. Spatial distribution of regional green space connectivity in Shanghai from 2000 to 2020.
Figure 9. Spatial distribution of regional green space connectivity in Shanghai from 2000 to 2020.
Land 12 02184 g009
Figure 10. Spatial distribution of regional green space diversity in Shanghai from 2000 to 2020.
Figure 10. Spatial distribution of regional green space diversity in Shanghai from 2000 to 2020.
Land 12 02184 g010
Figure 11. Spatial distribution of green space diversity in Shanghai from 2000 to 2020.
Figure 11. Spatial distribution of green space diversity in Shanghai from 2000 to 2020.
Land 12 02184 g011
Table 1. Landscape pattern index.
Table 1. Landscape pattern index.
Index NameFormulaDescription
Patch Type Area C A = j = 1 n a i j ( 1 10000 ) a i j is the area of ij patches (m2) divided by 10,000 and converted into CA.
Landscape Shape Index L S I = 0.25 E A E is the total length of all patch boundaries in the landscape, and A is the total area of the landscape.
Patch Fragmentation Index P D = N P 1 N C NP is the sum of various patches; NC is the total number of cells in the grid of the landscape data matrix. FN describes the difference between the total number of various patches in the landscape and 1 divided by the total number of grids in the landscape data matrix. The more patches in the landscape, the greater the PD value and fragmentation.
Shannon Diversity Index S H D I = i = 1 m ( P i × l n P i ) The Shannon diversity index reflects landscape heterogeneity, with larger values indicating a balanced distribution of various patch types in the landscape. This is manifested by a richer land use in the landscape, a higher the degree of fragmentation, and a greater uncertain information content.
Table 2. Statistical table of land use area in Shanghai.
Table 2. Statistical table of land use area in Shanghai.
Land Use Type20002005201020152020
Area (km2)RatioArea (km2)RatioArea (km2)RatioArea (km2)RatioArea (km2)Ratio
Agriculture4243.080.604056.850.573526.990.513266.670.473001.500.44
Forests215.590.030297.090.041443.580.064396.910.057454.140.067
Grasslands172.670.02487.580.01226.010.0037185.950.027189.650.028
Water1693.180.231630.480.231574.030.231441.990.211323.390.20
Constructed616.430.085902.180.131269.380.181401.280.201541.390.23
Wetlands192.310.027121.760.017143.530.021240.520.035261.910.039
Table 3. Statistical table of land use area in Shanghai.
Table 3. Statistical table of land use area in Shanghai.
Land Use Type20002005201020152020
Area (km2)ChangeArea (km2)ChangeArea (km2)ChangeArea (km2)ChangeArea (km2)Change
Agriculture−186.23−0.044−529.86−0.13−260.32−0.074−265.17−0.081−186.23−0.044
Forests81.510.38146.480.4946.670.1157.230.1481.510.38
Grasslands−85.088−0.4861.570.70159.946.153.710.019−85.09−0.48
Water−62.70−0.037−56.44−0.034−132.04−0.084−118.60−0.082−62.70−0.037
Constructed285.750.46367.210.41367.220.29140.110.099285.750.46
Wetlands−70.55−0.3721.770.1896.980.6821.390.089−70.55−0.37
Table 4. Shanghai regional green patch type area statistical table (km2), 2000–2020.
Table 4. Shanghai regional green patch type area statistical table (km2), 2000–2020.
Land Use Type20002005201020152020
Grasslands172.6787.5886.01185.95189.65
Forests215.59297.09443.58396.91454.14
Wetlands192.31121.76143.53240.52261.90
Total regional green patches580.56506.43613.12823.37905.70
Table 5. Shanghai regional green patch fragmentation statistical table (km2), 2000–2020.
Table 5. Shanghai regional green patch fragmentation statistical table (km2), 2000–2020.
Land Use Type20002005201020152020
Grasslands1.170.610.111.171.01
Forests0.821.151.751.541.91
Wetlands0.0160.0540.0280.0400.04
Total regional green patches2.011.811.892.752.97
Table 6. Statistical table of regional green space connectivity of Shanghai from 2000 to 2020.
Table 6. Statistical table of regional green space connectivity of Shanghai from 2000 to 2020.
Land Use Type20002005201020152020
Grasslands85.396886.505793.124285.986484.6414
Forests89.952890.898790.324290.980390.1874
Wetlands99.221798.76198.654699.255199.1913
Total regional green patches94.760694.269594.576895.115294.6495
Table 7. Statistical table of Shanghai green space diversity from 2000 to 2020.
Table 7. Statistical table of Shanghai green space diversity from 2000 to 2020.
Year20002005201020152020
SHDI0.7869260.7869260.7869260.7869260.786926
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ji, Y.; Zhang, L.; Gu, X.; Zhang, L. Study on the Spatial and Temporal Evolution of Regional Green Space Morphology Outside Built-Up Areas based on the Google Earth Engine and Biophysical Component Modeling. Land 2023, 12, 2184. https://doi.org/10.3390/land12122184

AMA Style

Ji Y, Zhang L, Gu X, Zhang L. Study on the Spatial and Temporal Evolution of Regional Green Space Morphology Outside Built-Up Areas based on the Google Earth Engine and Biophysical Component Modeling. Land. 2023; 12(12):2184. https://doi.org/10.3390/land12122184

Chicago/Turabian Style

Ji, Yiwen, Lang Zhang, Xinchen Gu, and Lei Zhang. 2023. "Study on the Spatial and Temporal Evolution of Regional Green Space Morphology Outside Built-Up Areas based on the Google Earth Engine and Biophysical Component Modeling" Land 12, no. 12: 2184. https://doi.org/10.3390/land12122184

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop