Next Article in Journal
Does Land Fragmentation Affect the Effectiveness of Fiscal Subsidies for Agriculture: Evidence from China
Next Article in Special Issue
Assessing Spatial Heterogeneity in Urban Park Vitality for a Sustainable Built Environment: A Case Study of Changsha
Previous Article in Journal
Risk Assessment in Sustainable Infrastructure Development Projects: A Tool for Mitigating Cost Overruns
Previous Article in Special Issue
Spatial Trade-Offs and Synergies between Ecosystem Services in Guangdong Province, China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evolution Model, Mechanism, and Performance of Urban Park Green Areas in the Grand Canal of China

1
School of Architecture & Art Design, Hebei University of Technology, Tianjin 300401, China
2
School of Architecture, Southeast University, Nanjing 210096, China
3
Spatial Planning Center, Yangtze Delta Region Institute of Tsinghua University, Zhejiang, Jiaxing 314006, China
*
Authors to whom correspondence should be addressed.
Land 2024, 13(1), 42; https://doi.org/10.3390/land13010042
Submission received: 13 November 2023 / Revised: 26 December 2023 / Accepted: 27 December 2023 / Published: 30 December 2023

Abstract

:
Urban park green areas are part of territorial space planning, shouldering the mission of providing residents with high-quality ecological products and public space. Using a combination of several measurement models such as the BCG (Boston Consulting Group) matrix, ESDA (Exploratory Spatial Data Analysis), MLR (Machine Learning Regression), GWR (Geographically Weighted Regression), and GeoDetector, this paper presents an empirical study on the changes in Urban Park Green Areas (UPGAs) in the Grand Canal of China. By quantitatively measuring the spatio–temporal evolution patterns of UPGAs, this study reveals the driving mechanisms behind them and proposes policy recommendations for planning and management based on performance evaluation. The spatio–temporal evolution of UPGAs and their performance in China’s Grand Canal are characterized by significant spatial heterogeneity and correlation, with diversified development patterns such as HH (High-scale–High-growth), HL (High-scale–Low-growth), LH (Low-scale–High-growth), and LL (Low-scale–Low-growth) emerging. The evolution performance is dominated by positive oversupply and positive equilibrium, where undersupply coexists with oversupply. Therefore, this paper recommends the implementation of a zoning strategy in the future spatial planning of ecological green areas, urban parks, and green infrastructure. It is also recommended to design differentiated construction strategies and management policies for each zoning area, while promoting inter-city mutual cooperation in the joint preparation and implementation of integrated symbiosis planning. Furthermore, the spatio–temporal evolution of the UPGAs in the Grand Canal of China is influenced by many factors with very complex dynamic mechanisms, and there are significant differences in the nature, intensity, spatial effects, and interaction effects between different factors. Therefore, in the future management of ecological green areas, urban parks, and green infrastructure, it is necessary to interconnect policies to enhance their synergies in population, aging, industry and economy, and ecological civilization to maximize the policy performance.

1. Introduction

1.1. Background

In the national spatial planning system, the planning for urban green areas is one of the most important projects, aimed at providing high-quality ecosystem services and public spaces for residents in China. Scientific planning of the quantity, quality, and spatial structure of the supply of urban green areas, especially Urban Park Green Areas (UPGAs), will significantly foster a more sustainable and livable urban environment [1,2]. In addition, as the construction of the Yangtze River Economic Belt, the protection and utilization of the Grand Canal and the construction of its cultural belt, and the ecological protection and high-quality development of the Yellow River Basin have been successively upgraded to national strategies, watershed spatial governance has become a key task of territorial spatial planning and a hotspot in academic research [3]. Therefore, quantitatively analyzing the spatio–temporal evolution patterns and driving mechanisms of UPGAs from a regional holistic perspective and evaluating the performance of land supply and demand will provide a basis for green infrastructure planning practice in watersheds and will help establish a technological system adapted to territorial spatial planning in these areas.

1.2. Literature Review

1.2.1. Urban Green Areas and Parkland

Current research on urban green areas mainly deals with planning methods and management policies [4,5], accessibility and satisfaction assessment [6,7], spatio–temporal dynamics and their influencing factors [8,9], value for ecosystem services and health [10,11], and other areas [12], while research on urban parkland focuses on areas such as spatial siting and configuration [13,14], needs assessment [15], planning methods [16], and willingness to pay [17]. Overall, there has been a large body of research on urban green areas and parkland, but still less attention has been paid to UPGAs [18]. For example, given that the unbalanced distribution of UPGAs has a significant impact on the well-being of residents, Li [19] and Xu [20] proposed an optimal spatial division plan for the service levels of UPGAs from the perspective of opportunity equity and spatial scale, based on the case studies of Taiyuan and Xuchang. Doll [21] assessed the greenness of UPGAs in Australia from the perspective of landscape preference and water consumption. Wang [22] quantitatively measured the equity of UPGAs in the central city of Beijing, analyzing the space to reveal a serious mismatch between the supply and demand of green areas. Yin [23] quantitatively examined the retention capacity of 176 urban park green areas within the Fifth Ring Road of Beijing for PM2.5 and endeavored to provide a basis for the design and construction of UPGAs to improve air quality. Biernacka [24] mapped and analyzed the dynamics of UPGAs in Poland, and they suggested the inclusion of informal green areas in urban planning. Engstrom [25] analyzed the advantages and disadvantages of using a hedonic price approach to capture the values of UPGAs in urban planning, and Ayele [26] studied the management model of UPGAs in Addis Ababa during rapid urbanization in Ethiopia.

1.2.2. Basin Planning and the Grand Canal

As Molle [27] and Antwi [28] put it, watershed planning has its origins in regional water resource management, but the embeddedness of other natural and ecological resource management systems has contributed to its gradual transformation into a new concept of socio–political life and spatial governance systems. Suhardiman [29] argued that watershed planning has evolved in Nepal as an arena for power operations and struggles, with its cross administrative boundaries jointly created by different government agencies. Essentially, watershed planning is a debate between multiple interests on development opportunities around the two perspectives of conservation and utilization, during which a large number of integrated modeling methods [30,31] and planning tools are created [32]. At present, the study of artificially excavated canals plays an important role in watershed planning, especially in China, the United States, Egypt, and India [33]. For the Grand Canal of China, the current research mainly focuses on the fields of cultural heritage, heritage and ecological protection, tourism development, land use, urban and rural spatial changes, and human habitat analysis. It has moved beyond the construction of water facilities and the development of cultural belts to green belts and economic belts [34] (Table 1).

1.2.3. Research Gaps and Questions

There are three shortcomings in the current research. First, there is a wealth of research on urban green areas separated from parkland, but fewer studies combining the two, and such studies mainly focus on single-city case studies, lacking an analysis of the whole area and not matching the needs of watershed planning. Second, studies on the Grand Canal are mainly concerned with culture, ecological protection, and tourism development. Few scholars have focused on land use and change, except for Xia [65], who analyzed the impact of green areas on the well-being of residents in the Hangzhou section of the Grand Canal. Third, the Grand Canal is essentially a regional cultural, economic, and green belt, but current research lacks spatial correlation analysis from a regional perspective, and discussion of the driving mechanism ignores the influence of spatial and interactive effects.
To address the aforementioned shortcomings, this paper introduces a combination of spatial measurement models to study the whole area of the Grand Canal and analyze the spatio–temporal evolution patterns of UPGAs to reveal the driving mechanisms behind them, and it proposes suggestions and strategies for green areas or green infrastructure planning based on the change performance evaluation. This study aims to (1) quantitatively measure the evolutionary patterns of the UPGAs in the Grand Canal and reveal their spatial effects through the BCG (Boston Consulting Group) matrix and ESDA (Exploratory Spatial Data Analysis) in both temporal and spatial dimensions; (2) quantitatively measure the direct influence of different influencing factors on their spatio–temporal evolution patterns, as well as the spatial and interactive effects of the influencing factors using the MLR (Machine Learning Regression) method, the GWR (Geographically Weighted Regression) method, and GeoDetector; (3) evaluate the performance of their changes according to the indicators defined in the United Nations Sustainable Development Agenda, which dynamically analyzed the match between land supply and population demand; and (4) provide suggestions for their management policy and planning design based on the analysis results.

2. Materials and Methods

2.1. Study Area

The Grand Canal of China, with a history of more than 1000 years and spanning thousands of kilometers from north to south, crosses a number of water systems from north to south, such as the Haihe River, the Yellow River, the Huaihe River, the Yangtze River, Taihu Lake, and the Qiantang River, and has served as an important water transportation channel as well as an artery for the economic and cultural exchanges between the north and south of China since ancient times. The Grand Canal of China, consisting of the Beijing–Hangzhou Grand Canal, the Sui–Tangshan Grand Canal, and the Zhedong Canal, was officially approved by UNESCO on 22 June 2014 for inclusion in the World Heritage List. The State Council issued the Outline of the Plan for the Protection, Inheritance, and Utilization of the Grand Canal Culture in February 2019, which defines the spatial scope of the Grand Canal basin as the two municipalities of Beijing and Tianjin and the six provinces of Hebei, Shandong, Jiangsu, Zhejiang, Henan, and Anhui, comprising a total of 86 cities. The study area in this paper is the same as the planned scope, covering all 86 cities (Figure 1).

2.2. Research Steps and Technical Route

This study is based on a variety of measurement models, and it is performed in four steps (Figure 2):
The first step is to analyze the spatio–temporal evolution patterns of the supply scale and per-capita area of UPGAs in the Grand Canal using the BCG (Boston Consulting Group) matrix, and to analyze their spatial heterogeneity and autocorrelation through ESDA (Exploratory Spatial Data Analysis).
The second step is to quantitatively analyze the driving mechanisms of changes in the scale of their supply and per-capita area, including the importance of factors, spatial effects, and interaction effects, based on a combination of MLR (Machine Learning Regression), GWR (Geographically Weighted Regression), and GeoDetector.
The third step is to analyze the match between their supply and demand from the perspective of sustainable development based on the LCRPGR (Ratio of Land Consumption Rate to Population Growth Rate).
The fourth step is to propose planning and management suggestions as guidance and a basis for green area planning and green infrastructure policy design based on the analysis results of the first three steps.

2.3. Research Methods and Indicator Selection

2.3.1. Boston Consulting Group (BCG) Matrix

The BCG matrix is often used in enterprise strategy management to classify the development status of businesses or products into four types, Star, Question, Cow, and Dog, and propose differentiated development strategies based on the combined analysis of the relative market share and growth rate of business departments or products. A development strategy is needed for Star products, and further investment is required to support their rapid expansion in the market and to make them the leading products of the enterprise. Cow products require a profit strategy, but the core of future development strategies is not to increase business investment, enterprise output, or market supply, but rsther to quickly recover funds through high product profits to support the development of Star products. For Question products, more research is needed and flexible strategies should be developed according to the actuality. Specifically, it is necessary to, based on the survey and research results, select some potential subcategories for key investments, while abandoning others. Dog products should be stopped, either by abandonment or recycling, in order to decisively stop loss, as they are the products that have been in a period of decline, unable to create more earnings for the enterprise.
In this paper, a concept is introduced to analyze the spatio–temporal evolutionary patterns of UPGAs in the Grand Canal of China, where Relative Share ( R S ) and Growth Rate ( G R ) represent the regional status in the spatial dimension and the growth capacity in the temporal dimension, respectively. With their median as the threshold, the spatio–temporal evolutionary patterns of UPGAs of the 86 cities in the study area can be categorized into four quadrants—HH (High-scale–High-growth), HL (High-scale–Low-growth), LH (Low-scale–High-growth), LL (Low-scale–Low-growth)—by the Cartesian coordinate system. With U P G i and U P G i ' as the indicators for UPGAs in the i -th city in 2020 and 2010, respectively, and U P G m a x as the maximum value for 86 cities in the study area, R S and G R are calculated as follows [66]:
R S = U P G i U P G m a x × 100 %
G R = U P G i U P G i ' U P G i ' 1 × 100 %

2.3.2. Exploratory Spatial Data Analysis (ESDA)

This paper employs ESDA to quantitatively judge and visualize the spatial features of UPGAs, including spatial autocorrelation and spatial heterogeneity. Moran’s I index is used to measure the overall strength of spatial autocorrelation, and a value greater or less than zero represents the positive or negative spatial autocorrelation in the spatial distribution of UPGAs, respectively, otherwise it represents a random distribution [67,68]. To further measure the localized characteristics of spatial associations, the Getis–Ord G i * index is introduced to classify the cities in the study area into four types: hot, sub-hot, sub-cold, and cold. To measure regional differences in UPGAs, coefficient of variation and spatial clustering methods are introduced to represent and visually demonstrate spatial heterogeneity. A larger coefficient of variation indicates a greater regional difference in UPGAs, with 0.36 and 0.16 being the thresholds to determine high and low levels of spatial heterogeneity [69]. With U P G ¯ being the mean of UPGA metrics, S being their standard deviation, and W i j being the spatial weight, Moran’s I, Getis–Ord G i * , and CV are calculated as follows [70]:
I = n i = 1 n j = 1 n W i j ( U P G i U P G ¯ ) ( U P G j U P G ¯ ) ( i = 1 n j = 1 n W i j ) i = 1 n ( U P G i U P G ¯ ) 2
G i * = j = 1 n W i j U P G i U P G ¯ j = 1 n W i j S [ n j = 1 n W i j 2 ( j = 1 n W i j ) 2 ] n 1
C V = S / U P G ¯ ,     S = i = 1 n ( U P G i i = 1 n U P G i n ) 2 n ,     U P G ¯ i = 1 n U P G i n

2.3.3. Machine Learning Regression (MLR)

Given the nonlinear characteristics of distribution planning for UPGAs and machine learning regression methods that do not rely on a priori subjective human experience, a nonlinear econometric model is adopted to analyze the importance of influencing factors [71]. In this paper, the decision tree, random forest, adaboost, and ExtraTrees algorithms in machine learning regression models are used to analyze the influence of different factors on the spatio–temporal evolution patterns of UPGAs. Decision tree, a tree-like structure, tests the data sample from the root node, divides the data sample into different data sample subsets according to different results, and calculates the data through a series of rules [72]. Random forest is a supervised machine learning algorithm constructed by integrating decision tree-based learners, which introduces randomness into the training process of decision tree to make it excellent in overfitting and noise resistance [73]. The AdaBoost model is an iterative algorithm that adds a new weak classifier in each round until a predetermined sufficiently small error rate is reached [74]. Extra-trees are derived from the traditional decision tree algorithm, characterized by the direct use of random features and thresholds in the node division of the decision tree, resulting in a larger and more random shape and difference in each decision tree [75]. The purpose of introducing machine learning regression in this paper is to measure the importance of factors, not to make predictions. Therefore, all data are included in the calculation, and the final result is determined according to goodness of fit and the comparative analysis of different algorithms.
Dependent variables include the supply scale of UPGAs and per-capita area of UPGAs, labeled with Y 1 and Y 2 , respectively. For independent variables, the combined influence of social, economic, and natural factors should be considered (Table 2). For social factors, population density represents the overall impact of a population aged 60 and above, and the proportion of population aged 60 and above represents the impact of special aging groups [76]. Outflow population indicates the impact of a semi-urbanized, transient population [77]. Of the economic factors, GDP represents the impact of the size of the economy [78]. Per-capita GDP represents the impact of the stage and quality of economic development [79,80] and fiscal self-sufficiency rate represents the government’s ability to intervene in the economy [81,82]. In terms of natural factors, topography represents the impact of topographic complexity [83], average temperature represents the impact of climate change, especially the urban heat island effect [84], and ventilation coefficient represents the impact of regional wind environment and urban air quality [85,86]. According to the measurement of covariance between factors using the least squares linear regression model, the maximum value of VIF for per-capita GDP among the nine independent variables reaches 8.25, but is still less than 10, indicating that the covariance of the independent variables is weak and can be almost ignored. The dependent variable data came from the China Urban Construction Statistical Yearbook; the population data came from the population census; the economic data came from the China City Statistical Yearbook and the statistical yearbooks of eight provinces/municipalities directly under the central government; the topographic relief data come from the Relief Degree of Land Surface Dataset of China (1 km) [87,88]; the average temperature data came from data.cma.cn; and ventilation coefficients were calculated from ECMWF re-analysis-interim data by the methods of Broner [89], Hering [90], and Chen [91]. Equations (6) and (7) are used for the positive and negative indicators in the standardization of dependent and independent variables, where D i + / is a standardized value, D i is the original value, and D M a x   a n d   D M i n are the maximum and minimum values of the original data, respectively.
D i + = D i D M i n D M a x D M i n + 0.001
D i = D M a x D i D M a x D M i n + 0.001

2.3.4. Geographically Weighted Regression (GWR)

In this study, GWR is used to analyze the impact of each factor on the spatio–temporal evolution patterns of UPGAs. GWR improves the computational accuracy of the regression model by creating localized regression equations for each city and incorporating the spatial autocorrelation and heterogeneity of UPGA changes into the regression process [92]. With Y i representing the spatio–temporal evolution pattern of UPGAs of the i -th city (HH, HL, LH, LL are assigned values of 4, 3, 2, and 1, respectively, in the calculation), X i k being the k -th independent variable (influencing factor), β 0 being a constant term, ( μ i , v i ) being the spatial location of the i -th city (geographic center of gravity coordinate), β k ( μ i , v i ) being the correlation between the variables of the i -th city, and ϵ i being the error of the regression equation, GWR is calculated as follows [93]:
Y i = β 0 ( μ i ,   v i ) + k β k ( μ i ,   v i ) X i k + ϵ i

2.3.5. GeoDetector

In this study, GeoDetector is used to measure the interaction between different factors. Different factors interact with each other when they act together in the UPGA planning, and GeoDetector measures the interaction effect of factor pairs using the q-index. It calculates the spatial pattern of the dependent variable Y i and the similarity of independent variables X i k and X i l to obtain the single-factor and dual-factor influences q ( X i ), q ( X j ), and q ( X i X j ); furthermore, it compares q ( X i X j ) and other parameters to select and identify the final result—nonlinear weaken (q ( X i X j ) < Min q ( X i ), q ( X j )), single weaken (Min (q ( X i ), q ( X j ) < q ( X i X j ) < Max q ( X i ), q ( X j )), double enhance (q ( X i X j ) > Max q ( X i ), q ( X j )), independent (q ( X i X j ) = q ( X i ) + q ( X j )), and nonlinear enhance (q ( X i X j ) > q ( X i ) + q ( X j )) [94,95]. With h = 1, 2, 3, … l , where l is the number of partitions of spatial clustering, σ 2 is the total variance of dependent variables, σ h 2 is the variance of dependent variables of the h -th partition, and S S W and S S T are the sums of variances within the partition and the study area, the index q is calculated as follows [96]:
q = 1 h = 1 l n h σ h 2 n σ 2 = 1 S S W S S T ,    S S W = h = 1 l n h σ h 2 ,   S S T = n σ 2

2.3.6. Ratio of Land Consumption Rate to Population Growth Rate (LCRPGR)

The Transforming Our World—the 2030 Agenda for Sustainable Development proposes 17 SDGs (sustainable development goals). Indicator SDG 11.3.1 is defined as the ratio of the Land Consumption Rate (LCR) to the Population Growth Rate (PGR) and is used to represent the relationship between urban expansion and population change [97]. This study chooses to use this method to evaluate the performance of UPGAs from a sustainable development perspective. LCR is a reflection of the growth rate of land used for urban park purposes and represents the efficiency of changes in the supply of urban green areas. PGR reflects the change rate of urban population and is used to measure the change rate of the green area demand of the population in an area over a period of time. LCRPGR measures the relationship between the change rates of two variables, LCR and PGR, and is used to represent the match between supply and demand in UPGAs. Theoretically, an LCRPGR equal to 1 is the most desirable result, and in view of the elasticity in practice development, 0.75 and 1.25 are set as thresholds to classify the analysis results into eight categories (Table 3). LCRPGR is calculated as follows [98]:
L C R P G R = L C R P G R   = l n ( U P G i / U P G i ' ) n l n ( P D i / P D i ' ) n

3. Results

3.1. Spatiotemporal Evolution Model

3.1.1. Supply Scale

According to the relative share of supply scale of UPGAs, the coefficient of variation and Moran’s I are 1.47 and 0.05, respectively (Z = 1.90, p < 0.05), indicating huge inter-city differences and significant positive spatial autocorrelation. Most of the high-value cities are concentrated in Shandong, Beijing, Tianjin, and provincial capital metropolitan areas such as Zhengzhou, Hefei, Hangzhou, Shijiazhuang, and Nanjing. The Nanjing metropolitan area extends to cover the region of Southern Jiangsu (Suzhou, Wuxi, Changzhou, etc.). In addition, Ningbo, Linyi, Yantai, Wuxi, Nantong, Luoyang, and Zibo also have a leading edge, with a relative share of over 10%. Most of the low-value cities are concentrated in the north and south of Hebei and Henan and in the west of Anhui and Zhejiang, especially Zhoushan, Jinhua, Hengshui, Tongling, Lu’an, Xuchang, Xinxiang, Anyang, Xinyang, Puyang, Suzhou, Zhumadian, Hebi, Cangzhou, Sanmenxia, Zhoukou, Bozhou, Quzhou, Xuancheng, Huangshan, Chizhou, and Lishui, and have a large disadvantage, with a relative share of not more than 3%. The hotspots are concentrated in Beijing, Tianjin, and the north of Hebei Province; the sub-hotspots are in the Shandong Peninsula and the densely populated urban areas of Southern Jiangsu; and most of the coldspots are located in the border areas of Henan, Anhui, Shandong, and Jiangsu Provinces (Figure 3).
According to the growth of UPGA supply scale, the CV and Moran’s I are 0.70 and 0.01 (Z = 0.61, p > 0.05), respectively, indicating that neither the inter-urban differences nor the spatial correlations are significant. High-value cities are not spatially clustered geographically, and include Nantong, Fuyang, Zhengzhou, Wenzhou, Jining, Kaifeng, Ningbo, Chuzhou, Shangqiu, Yancheng, Taizhou, Qingdao, Suqian, and Liaocheng. Most of the low-value areas are clustered in northern Hebei and southeastern Shandong, with Pingdingshan, Luohe, Xinxiang, Qinhuangdao, Rizhao, Huaibei, Wuxi, Handan, Chizhou, Zhenjiang, Zhangjiakou, and Tangshan lagging behind in development, and Chengde in particular showing negative growth. The hotspot cities are mainly in the border areas of Henan, Anhui, and Shandong Provinces and extend to central Jiangsu. The sub-hotspot cities are distributed in the periphery of the hotspot cities in a “center-edge” structure. In southeastern Zhejiang, a small “center-edge” structure is developing. The coldspots are concentrated in Beijing, Tianjin, northern Hebei Province, central Shandong Province, and the border areas of Anhui and Zhejiang Provinces.
According to the spatiotemporal evolution model of the scale of UPGA supply, the median relative share and growth rates are 4.66% and 68.92%, respectively. A high proportion is found in HH and LL cities, both at approximately 30%. HH cities are scattered in distribution and are only in western Shandong; coastal Jiangsu; and the provincial capital metropolitan areas of Henan, Anhui, and Zhejiang, including Beijing, Tianjin, Xingtai, Changzhou, and others. LL cities are relatively clustered in the north of Hebei, the south of Anhui, the west of Zhejiang, and the northeast corner of Henan, including Zhangjiakou, Chengde, Cangzhou, Langfang, Zhenjiang, Jiaxing, Jinhua, Quzhou, Zhoushan, Lishui, Ma’anshan, Huaibei, Huangshan, Lu’an, Chizhou, Xuancheng, Zaozhuang, Rizhao, Binzhou, Pingdingshan, Anyang, Hebi, Xinxiang, Puyang, Luohe, and Xinyang. Most of the HL cities are concentrated in Shandong, Hebei, and Jiangsu, including Shijiazhuang, Tangshan, Qinhuangdao, Handan, Baoding, Nanjing, Wuxi, Xuzhou, Suzhou, Yangzhou, Huzhou, Taizhou, Huainan, Zibo, Yantai, Weifang, and Linyi. Most LH cities are concentrated in Henan and Anhui, especially in the border areas of the two provinces, including Hengshui, Taizhou, Suqian, Bengbu, Tongling, Anqing, Chuzhou, Suzhou, Bozhou, Heze, Kaifeng, Jiaozuo, Xuchang, Sanmenxia, Shangqiu, Zhoukou, and Zhumadian. Overall, most of the hotspot cities are clustered in the west of Shandong and extend to Jiangsu and Hebei, while most of the coldspot cities are in the central part of Henan and the border area of Anhui and Zhejiang, both geographically distributed in a band (Figure 4).

3.1.2. Per-Capita Area

According to the relative share of the per-capita area of UPGAs, the CV and Moran’s I are 0.19 and 0.17, respectively (Z = 3.93, p < 0.01), indicating moderate inter-city differences but significant positive spatial autocorrelation. Most of the high-value cities are concentrated in Shandong, with a few in the northern end of Hebei and southern Anhui. In addition, Chuzhou, Bozhou, Fuyang, Nantong, Yangzhou, and other cities also have a leading edge, with a relative share close to 70%. Most of the low-value cities are concentrated in Zhejiang, the western part of Hebei, the border area of Henan and Shandong, and the central part of Anhui, especially Taizhou, Wuhu, Pingdingshan, Hefei, Anyang, Suzhou, Jinhua, Hangzhou, Xinxiang, Jinan, Lishui, Cangzhou, Zhangjiakou, and Tianjin, and they have a significant disadvantage, with a relative share of less than 50%. The hotspot cities are all clustered in Shandong and the sub-hotspot cities are in its periphery and extend to Hebei and Jiangsu, forming a “center-edge” structure. There are three clusters of coldspot cities in Zhejiang (except Hangzhou), Beijing–Tianjin, and the east of Henan, while all other cities are sub-hotspots (Figure 5).
According to the growth rate of the per-capita area of UPGAs, the CV and Moran’s I are 1.17 and 0.20 (Z = 4.53, p < 0.01), respectively, indicating significant spatial heterogeneity and spatial autocorrelation. Most of the high-value cities are concentrated in Henan and extend in a continuous belt towards Anhui and a necklace (stepping stone) towards Shandong. Low-value cities are clustered in northern Hebei, central Zhejiang, Shandong, and southern Jiangsu, with Zhangjiakou, Shaoxing, Hefei, Qinhuangdao, Yantai, Handan, Rizhao, Hangzhou, Chengde, and Suzhou in particular showing negative growth. The hotspot cities are mainly concentrated in the border area of Henan and Anhui Provinces, while the sub-hotspot cities are in its periphery, forming a “center-edge” structure. There are three coldspot urban clusters in northern Zhejiang, the peninsula and south of Shandong, and the eastern part of the Beijing–Tianjin–Hebei metropolitan area (Beijing–Tangshan–Qinhuangdao region).
According to the spatiotemporal evolution model of per-capita area of UPGAs, the median relative share and growth rate are 59.33% and 22.86%, respectively. A high proportion is found in HL and LH cities, both at approximately 25%. HH cities are dispersed in a geographical distribution, including Xingtai, Langfang, Nantong, Taizhou, Suqian, Quzhou, Huaibei, Tongling, Chuzhou, Fuyang, Lu’an, Bozhou, Qingdao, Beijing, Zibo, Dongying, Jining, Binzhou, Luoyang, and Xuchang. Most of the LL cities are concentrated in Hebei and Zhejiang, a small number are located in the densely populated urban areas of southern Jiangsu, and very few are randomly distributed. LL cities include Tianjin, Shijiazhuang, Baoding, Zhangjiakou, Cangzhou, Wuxi, Changzhou, Suzhou, Hangzhou, Jiaxing, Shaoxing, Jinhua, Taizhou, Lishui, Hefei, Huainan, Zaozhuang, Liaocheng, Puyang, and Xinyang. Most of the HL cities are concentrated in Shandong, and there are two small clusters in the southern end of Anhui and the northeast corner of Hebei. These include Tangshan, Qinhuangdao, Handan, Chengde, Nanjing, Xuzhou, Yangzhou, Zhenjiang, Huzhou, Zhoushan, Huangshan, Chizhou, Xuancheng, Yantai, Weifang, Tai’an, Weihai, Rizhao, Linyi, Dezhou, Hebi, Luohe, and Sanmenxia. Most of the LH cities are concentrated in Henan and extend to northern Anhui and northern Jiangsu, including Hengshui, Lianyungang, Huai’an, Yancheng, Ningbo, Wenzhou, Wuhu, Bengbu, Ma’anshan, Anqing, Suzhou, Jinan, Heze, Zhengzhou, Kaifeng, Pingdingshan, Anyang, Xinxiang, Jiaozuo, Nanyang, Shangqiu, Zhoukou, and Zhumadian. Overall, most of the hotspot cities are clustered in Shandong; sub-hotspot cities are mostly found in Henan, Anhui, Jiangsu, Hebei, Beijing, and Tianjin; most of the coldspot cities are in Zhejiang; and most of the sub-coldspot cities are in the border areas of Hebei, Henan, and Shandong, showing significant clustering characteristics (Figure 6).
The comparative analysis in the above two dimensions shows that Anyang, Pingdingshan, Puyang, Shangqiu, Xinyang, and Zhengzhou are always at a low level in the supply and per-capita area of urban green areas, with the construction of urban green infrastructure far behind that of other cities along the Grand Canal. Therefore, more investment and support are needed in future planning, construction, and management. On the contrary, Qingdao, Linyi, Nantong, Zibo, Weifang, Jining, Weihai, and Dongying are consistently at a high level, with the urban green infrastructure development ahead of the other cities in the Grand Canal, making them of exemplary value in the region. It should be noted that historical and cultural cities and livable cities such as Hangzhou, Beijing, Nanjing, Tianjin, and Suzhou are regional leaders in the scale of urban green areas supply, but they have no competitive advantage per capita. They should attach special attention and focus to green infrastructure planning, construction, and management in the future.

3.2. Driving Mechanism

3.2.1. Importance and Nature of Factors

The analysis results from the machine learning regression algorithms show a high goodness, generally greater than 0.85, with decision tree being optimal in both schemes. The results of the four algorithms are in general similar. Although there are differences in the coefficients of factors with higher and lower importance, their ranking remains relatively stable (Table 4). To provide full play to the advantages of all algorithms and eliminate the defects of a single algorithm, this paper uses the average values of four algorithms to determine the importance of factors. For the supply scale of UPGAs, topography, outflow population, and population density have much higher importance than other factors, and they are defined as key factors; fiscal self-sufficiency rate and average temperature have much lower importance than other factors, and they are defined as auxiliary factors with direct influences that can be ignored; and ventilation coefficient, proportion of population aged 60 and above, GDP, and per-capita GDP are not more or less important and are defined as important factors. For per-capita area of UPGAs, per-capita GDP is a key factor; fiscal self-sufficiency rate, GDP, and average temperature are auxiliary factors; while topography, ventilation coefficient, population density, outflow population, and proportion of population aged 60 and above are important factors (Figure 7).

3.2.2. Spatial Effect of Factors

From the supply scale of UPGAs, the minimum values of proportion of population aged 60 and above ( X 2 ), GDP ( X 4 ), and fiscal self-sufficiency rate ( X 6 ) are all greater than zero, suggesting that they all play a positive role as a whole. The maximum values of topography ( X 7 ), average temperature ( X 8 ), and ventilation coefficient ( X 9 ) are all less than zero, suggesting that they all act as negative obstacles overall. The maximum values of population density ( X 1 ), output population ( X 3 ), and per-capita GDP ( X 5 ) are greater than zero, while the minimum values are less than zero, indicating that they both have positive driving and negative blocking effects with a complex impact mechanism (Table 5). The influence of population density shows a “dumbbell” pattern geographically, decreasing from the north and south to the middle. Beijing, Tianjin, and Hebei in the north are highlands of positive effects, while Henan is a depression of negative effects and Zhejiang is a new highland of positive effects. The proportion of population aged 60 and above is characterized by coastal highs and inland lows, with Bohai Bay (Beijing, Tianjin, Hebei, Shandong) as the high ground, the Yangtze River Delta as the second high ground, and Henan and Anhui as the depressions, with the weakest in Henan. The influence of the outflow population also presents a “dumbbell” pattern geographically, with the Yangtze River Delta, especially Zhejiang and southern Anhui Provinces, being the highlands of positive effects and Beijing–Tianjin–Hebei being the highlands of negative effects. The depressions are distributed in the border areas of Henan, Shandong, Anhui, and Jiangsu Provinces. The highland of GDP influence is in the west of Anhui and the south of Henan, and it is the origin of the gradient to the north and the coast, reaching the lowest in the Beijing–Tianjin–Hebei region. Per-capita GDP is characterized by a decreasing gradient from south to north, with highlands in Zhejiang and southern Anhui and depressions in Bohai Bay and the Shandong Peninsula. The fiscal self-sufficiency rate has a high impact on the northern region, with Beijing and Hebei as the highlands, and low for the center, with Anhui and the northern border region of Jiangsu as the depressions. The influence increases in a gradient from the depression to the south and reaches the highest in Zhejiang. The influence of topography, average temperature, and ventilation coefficient is characterized by clustering, and the depressions are all located in the Yangtze River Delta, especially in the southern part of Zhejiang and Anhui Provinces. Their highlands are all clustered in bands, but with differences in geographic location, where topography is located in the Shandong Peninsula, the average temperature is located in the junction area of Shandong, Henan, Anhui, and Jiangsu Provinces, and the ventilation coefficient is in the north of Henan and the west of Hebei (Figure 8).
From the per-capita area of UPGAs, only the minimum value of per-capita GDP ( X 5 ) is greater than zero, indicating that it plays a positive driving role. The maximum values of outflow population ( X 3 ), GDP ( X 4 ), topography ( X 7 ), average temperature ( X 8 ), and ventilation coefficient ( X 9 ) are all less than zero, suggesting that they all act as negative obstacles overall. The maximum values of population density ( X 1 ), proportion of population aged 60 and above ( X 2 ), and fiscal self-sufficiency rate ( X 6 ) are greater than zero, while the minimum values are less than zero, indicating that they both have positive driving and negative blocking effects with a complex impact mechanism (Table 6). The influence of population density is characterized by a geographic gradient of “high in the coastal area and low in the inland area”, with the two clusters in the Shandong Peninsula and the northern end of Hebei being the highlands, and the contiguous band-like areas in Henan, Anhui, and western Zhejiang being the depressions. The influence of the proportion of population aged 60 and above presents a “dumbbell” pattern geographically, with the Yangtze River Delta as the highland of positive influence and Beijing–Tianjin–Hebei as the highland of negative influence, divided by the border area of Henan, Shandong, and Hebei. The influence of the outflow population is geographically characterized by “high in the north and low in the south”, with Beijing–Tianjin–Hebei being the highland and the Yangtze River Delta being the depression. The influence of GDP is characterized by a geographic gradient of “high in the inland area and low in the coastal area”, with the highlands located in the border area between Anhui and Henan and the depressions in the Shandong Peninsula and the east coast of Zhejiang. Per-capita GDP and outflow population are highly similar in their geographic patterns of influence, except that the nature of the influence shifts from negative to positive. The influence of the fiscal self-sufficiency rate, the proportion of population aged 60 and above, and topography is geographically characterized by a “low in the center and high at the edge”, with depressions located in the border areas of Anhui, Henan, Shandong, and Jiangsu and highlands in the Shandong Peninsula, both with a small spatial scale. The average temperature influence highland is located in the border areas of Zhejiang, Anhui, and Jiangsu Provinces, and the depression is in Henan. The geographic pattern of influence of ventilation coefficient and outflow population is similar; however, the coverage of the latter uplands and depressions is larger than that of the latter, with highlands extending to the north and west of Henan and depressions expanding to the south of Jiangsu and the south of Anhui (Figure 9).

3.2.3. Interactive Effect of Factors

Different factors show a significant synergistic enhancement effect, mainly in the form of nonlinear enhancement, with only a few factor pairs in double enhancement. For the supply scale of UPGAs, the factor pair of topography ∩ average temperature ( X 7 X 8 ) is in double enhancement. For the per-capita area of UPGAs, population density ∩ proportion of population aged 60 and above ∩ proportion of population aged 60 and above ( X 1 X 2 ), population density ∩ GDP ( X 1 X 4 ), proportion of population aged 60 and above ∩ GDP ( X 2 X 4 ), and proportion of population aged 60 and above ∩ average temperature ( X 2 X 8 ) are in double enhancement. It is worth noting that a large number of super-factor pairs arise from factor interaction, and their interaction forces are much higher than those of other factor pairs and single factors. Population density ∩ GDP ( X 1 X 4 ), proportion of population aged 60 and above ∩ GDP ( X 2 X 4 ), outflow population ∩ GDP ( X 3 X 4 ), per-capita GDP ∩ GDP ( X 5 X 4 ), and ventilation coefficient ∩ GDP ( X 9 X 4 ) are super factor pairs of the supply scale of UPGAs, with an interaction force of more than 0.80. Population density ∩ per-capita GDP ( X 1 X 5 ), outflow population ∩ Per-capita GDP ( X 3 X 5 ), outflow population ∩ average temperature ( X 3 X 8 ), outflow population ∩ ventilation coefficient ( X 3 X 9 ), per-capita GDP ∩ average temperature ( X 5 X 8 ), and per-capita GDP ∩ ventilation coefficient ( X 5 X 9 ) are super factor pairs of the per-capita area of UPGAs, with an interaction force of more than 0.5 (Table 7).

3.3. Performance Evaluation

No cities in the Grand Canal region fall into the categories of negative oversupply, negative undersupply, or negative equilibrium. Positive oversupply has the highest proportion, more than 50%, followed by positive equilibrium at approximately 30%. Ma’anshan, Huaibei, Binzhou, Luoyang, Pingdingshan, and Shangqiu are super oversupply members. The only member of the super undersupply group is Chengde, while the only member of the positive undersupply group is Handan. The positive oversupply members include Beijing, Tangshan, Xingtai, Cangzhou, Langfang, Hengshui, Nanjing, Nantong, Lianyungang, Huai’an, and others. Most of these are clustered in Henan and Anhui and extend in a band towards Jiangsu, Zhejiang, and Hebei. The positive equilibrium members include Tianjin, Shijiazhuang, Qinhuangdao, Baoding, Zhangjiakou, Wuxi, Xuzhou, and others. These are distributed in bands in Shandong and Hebei and in clusters in Jiangsu and Zhejiang. Overall, there are three clusters of hotspot cities in Shandong, Jiangsu, and Zhejiang and three clusters of coldspot cities in Henan, the northern end of Hebei, and the border area between Henan and Anhui, with sub-hotspots and sub-coldspots distributed on their periphery, forming a “center-edge” spatial structure (Figure 10).

4. Discussion

4.1. Differentiated Zoning Planning

The study shows significant spatial inequalities and differences in the geographic distribution of UPGAs in the Grand Canal, with huge disparities in the supply scale and per-capita area across cities, including relative shares, growth rates, and spatial and temporal evolution patterns. In addition, most of the UPGA supply is in a mismatch with the population demand, with co-existence of both oversupply and undersupply. The high similarity of these analytical results to the findings of other scholars suggests that spatial imbalances and inequalities, supply–demand imbalances, and mismatches are regular features. For the former, Rigolon [99] concluded that there are significant spatial inequalities in urban parkland and quality and Ren [100] found geographic and social inequalities in the distribution of urban parks in Shanghai based on Gini coefficient analysis. For the latter, Tan [101], Gao [102], and Zhu [103] in their case studies of Wuhan, Shenzhen, and Beijing found that there is a serious mismatch between UPGA supply and resident demand, which is a great challenge for future planning and management. Unfortunately, current green areas planning, park planning, and green infrastructure planning focuses more on the design of spatial layout schemes for intra-city parks and green area systems, with less attention and planning response to inter-city imbalances and inequalities in general [104].
In summary, we recommend the adoption of differentiated zoning planning strategies in green areas and park planning for the Grand Canal. This takes the evaluation of the performance of matching supply and demand as the core basis to adjust the direction of the control of the supply scale and per-capita area in cities according to the classification results of the spatio–temporal evolution mode of urban parkland. For cities such as Tianjin, Shijiazhuang, Qinhuangdao, Baoding, Zhangjiakou, and Wuxi, as their supply and demand are in positive balance, the focus of their future planning will be to keep their management policies stable and to maintain and contribute to the city’s long-term balance of supply and demand. Handan and Chengde face a serious supply shortage, with UPGA supply size and per-capita area in HL and LL states, resulting in weak growth. Therefore, in the future, provincial governments should increase the quota of their UPGAs, and city governments should adopt speed control-oriented planning and policies to accelerate incremental supply and promote a balance between supply and demand. For Tangshan, Xingtai, Cangzhou, Langfang, Hengshui, Nanjing, Nantong, Lianyungang, and other positive oversupply cities, subdivision planning based on spatio–temporal evolution patterns is required. Cities with spatio–temporal evolution patterns in the HH and LH categories, such as Xingtai, Nantong, Dongying, Fuyang, and Anqing, should strictly limit the growth of quotas in the future and promote a balance between supply and demand by reducing land supply. They should focus UPGA planning and management on quality improvement rather than quantity growth in the future, and they could sell their surplus land index to undersupply cities via the regional inter-city trading platform. Cities with spatio–temporal evolution in LL and HL patterns, such as Cangzhou, Tangshan, Jiaxing, Huzhou, Lishui, and Xuancheng, should increase or purchase land quotas and strengthen the planning, design, and management of urban green areas to promote a balance between supply and demand through the growth of land supply and the strict protection of the stock.

4.2. Integrated Symbiotic Planning

We found significant spatial correlations of UPGAs in the Grand Canal in this study, with hotspot and coldspot cities clustering together. It should be noted that Choumert [105] and Kim [106] also noted such correlations, and in common with this paper they both found significant positive spatial autocorrelations in UPGAs between cities in France and South Korea, with local spatial clustering features prominent, although the global correlation is weak or insignificant. However, unfortunately, they did not propose a response strategy for spatial planning and green area management policies based on spatial correlation characteristics. We believe that clustered hotspot and coldspot cities face similar development challenges, and they can maximize their performance with minimum cost through inter-city collaboration in UPGA planning and management. Highly interconnected cities may seek to rapidly increase the quantity and quality of their own and regional green infrastructure through the construction of point-like regional parks and linear inter-city greenways.
For the construction of regional parks, cities can learn from the European experience of building large-scale green parks across administrative districts based on the characteristics of natural resources such as the Grand Canal’s green areas and water system, as well as planning and constructing open spaces, high-quality landscapes, and recreational and activity facilities to promote the symbiosis of spatial functions and social needs [107,108]. In urban green area system planning, park system planning, and green infrastructure planning in China, it is common to keep the planning scope in line with the scope of the city’s overall planning, i.e., the planning is more confined to the central urban area, which has seriously weakened the overall allocation of green resources in counties and towns within the municipal area and has led to no overall planning for the regional green area system. The construction of regional parks that are not confined to the central city and encouraged coordination between neighboring cities will contribute to the integration of regional ecological, spatial, tourism, and social resources to enhance regional sustainability and competitiveness. The regional greenway is a linear green open space. Natural and artificial landscapes such as mountains, lakes, fields, scenic spots, ancient cities, cultural heritages, traditional villages, and tree-lined roads along the main stream and tributaries of the Grand Canal should be connected together, and landscape recreation routes and service facilities serving pedestrians and cyclists should be built within them [109]. Hotspot cities can link and showcase quality parks and green areas in the united cities through regional greenways to enhance the overall image of the regional habitat. For coldspot cities, they can promote the sharing of parks and green area resources among different cities through regional greenways to enhance utilization efficiency and alleviate, to a certain extent, the imbalance between supply and demand.

4.3. Multi-Policies Mix Management

Both the supply scale and per-capita areas of UPGAs have very complex driving mechanisms. In terms of nature, the ventilation coefficient always plays a negative blocking role, while population density always has both positive driving and negative blocking effects. The roles of other factors are always in a complex state of change and may shift from positive to negative (e.g., GDP), mixed (e.g., proportion of population aged 60 and above), or from negative to mixed (e.g., average temperature and topography) or from mixed to positive (e.g., per-capita GDP) and negative (e.g., outflow population). In terms of intensity, the key factors of supply scales of UPGAs are completely different from those of their per-capita area. Proportion of population aged 60 and above and ventilation coefficient are always important factors for both, while average temperature and fiscal self-sufficiency rate are always common auxiliary factors. The intensity ratings of the other factors differ widely. For example, GDP has been reduced from an important factor of the supply scale of UPGAs to an auxiliary factor of their per-capita area. Per-capita GDP, by contrast, has been upgraded from an important factor to a key factor. In terms of interaction, the super factor pairs of the supply scale of UPGAs are completely different from those of their per-capita area. In terms of spatial effects, all factors show a significant spatial clustering of influence on geographic patterns, but the geographic distribution and spatial extent of highlands and depressions vary widely, with a variety of change patterns such as gradient rise and fall, high in the center and low at the edge, low in the center and high at the edge, dumbbell (high at the ends and low in the middle), contiguous belt, and clustering emerging (Figure 11).
A growing body of research reveals that the planning and management of UPGAs are complex systematic projects subject to the influence of many factors. Many scholars have discussed the factors influencing changes in urban parks. For example, Cheng [110], Nam [111], and Smith [112] argued that government funding plays a key role in the management of urban parks in China and the United Kingdom. Luo [113], Feng [114], and Kim [115] argued that both population density and size had significant spatial correlations with the level of service of urban parks. Guo [116,117] found that house prices, transportation accessibility, and the status of the surrounding commercial facility package are important factors influencing the accessibility of urban parks. Their findings are corroborated with the conclusions of this paper. Different from them, this paper fully explores the spatial effects and interaction effects of different factors, which has a great inspirational value for the design of policy combination. In view of the intensity, nature, and spatial and interaction effects of different factors, future planners and the government should better interconnect the policies for UPGAs. Since a single policy is limited and difficult to operationalize in practice, multiple policies should be designed and implemented in the management. In addition, attention should be paid to the mode of combining policies in policy design and implementation, integrating policies on population, ageing, social mobility, economic and industrial development, financial investment, and the building of an ecological civilization. While ensuring the precision of each policy, it is important to maximize the synergy of different policies and to maximize policy performance by means of the interaction effect of factor pairs.

5. Conclusions

This paper arrives at the following findings: (1) The spatio–temporal evolution patterns of the UPGAs in the Grand Canal of China are diversified, with many types emerging, such as HH (High-scale–High-growth), HL (High-scale–Low-growth), LH (Low-scale–High-growth), and LL (Low-scale–Low-growth). (2) The evolutionary performance of UPGAs in the Grand Canal is mainly characterized by positive oversupply and positive equilibrium, with super oversupply, super undersupply, and positive undersupply found in a small number of cities. (3) The spatio–temporal evolution patterns and performance of UPGAs are characterized by significant spatial heterogeneity and positive spatial autocorrelation, with huge inter-city differences, and both hotspot and coldspot cities are clustered geographically. (4) The spatio–temporal evolution of UPGAs is driven by a complex mechanism, and different factors vary greatly in nature, intensity, spatial effect, and interaction effect. (5) The planning and management of UPGAs in the Grand Canal should be implemented by classifying and zoning, and zoning planning and symbiosis planning should be prepared and implemented based on the results of the analysis. In addition, it is necessary to design differentiated and diversified policies for each planning zone in the future, and to focus on enhancing the synergy of multiple policies in the management, so as to maximize the benefits based on the “combination of policies”.
This paper presents innovations in the following areas: (1) It pushes UPGA research to shift from case studies of a single city to systematic studies of regional urban agglomerations. It is a big step forward as a single system is more than the sum of its parts. Against the backdrop of China’s urban development entering a new era of regional integration dominated by urban agglomerations, urban belts, urban contiguous areas, and metropolitan areas, the research conclusions of individual cities reached in the past are not entirely applicable to current regional green infrastructure planning, although they have provided good guidance for urban scale green space system planning. (2) By integrating the BCG (Boston Consulting Group) matrix, ESDA (Exploratory Spatial Data Analysis), MLR (Machine Learning Regression), GWR (Geographically Weighted Regression), GeoDetector, and LCRPGR (Ratio of Land Consumption Rate to Population Growth Rate), it comprehensively and systematically reveals the driving mechanisms behind UPGAs while quantitatively evaluating their spatio–temporal evolution patterns and performance, especially analyzing in detail the spatial and interactive effects of different factors. It is a brand-new exploration and discovery. (3) Instead of being limited to or stagnated in the analysis of the change characteristics of UPGAs, this study proposes differentiated zoning planning, integrated symbiotic planning, and multi-policies mix management based on the design of spatial planning and management policies for green areas, parks, and green infrastructure in the Grand Canal of China. It is a remarkable fact that canals are common in all countries of the world, and the technical approach, analytical methods, and results of this paper are not only applicable to China, but can also be used as the basis and reference for canal planning and management in Egypt, India, Indonesia, America, and other countries.
There are still some shortcomings in this study: (1) Due to data and information limitations, this study only took into account the effect of population size in the performance evaluation, while it did not include the heterogeneity of the needs of different populations in the analytical model, which may affect the accuracy of the analysis results. (2) This paper is based on a regional inter-city comparative study with no in-depth analysis of parkland in a single city within the Grand Canal, especially within the key cities, which somewhat constrains the breadth of application of the results. The canal is not only a navigation channel and cultural heritage zone for human beings, but also a natural ecological zone and a source of green well-being for the people living along the route. Beginning with the construction of ecological civilization and ending with the sustainable development of the Grand Canal of China, this study focuses on the regional analysis and planning of UPGAs, expanding the field of canal research from the traditional dimensions of shipping and transportation, water resources, and history and culture to the dimension of green areas. It demonstrates certain theoretical innovations while responding to the practical needs of ecological protection and high-quality development of the Grand Canal.

Author Contributions

Conceptualization, S.Z. and Z.C.; methodology, S.Z. and M.H.; software, M.H. and S.Z.; validation, Z.C., C.Z. and S.Z.; formal analysis, Z.C. and C.Z.; investigation, Z.C. and M.H; resources, Z.C. and S.Z.; data curation, S.Z. and Z.C.; writing—original draft preparation, Z.C. and S.Z.; writing—review and editing, C.Z. and S.Z.; visualization, Z.C. and C.Z.; supervision, M.H. and S.Z.; project administration, Z.C.; funding acquisition, Z.C. and C.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Social Science Foundation Project of Hebei Province (HB22YS031).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data used in this paper mainly came from the Ministry of Housing and Urban-Rural Development of the People’s Republic of China (https://www.mohurd.gov.cn/index.html, accessed on 12 March 2023), the National Bureau of Statistics (http://www.stats.gov.cn/sj/ndsj/, accessed on 9 February 2023), the Global Change Research Data Publishing & Repository (https://www.geodoi.ac.cn/WebCn/doi.aspx?Id=887, accessed on 13 April 2023), and the Atmospheric Composition Analysis Group of Dalhousie University in Canada (https://www.heywhale.com/mw/dataset/641cfe42a001a17c4784f212, accessed on 18 April 2023).

Acknowledgments

Thank you to Menghan Guo and Mushuang Sun for their assistance in data collection and literature review. All the authors are grateful to the reviewers and editors.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mwanzu, A.; Nguyu, W.; Nato, J.; Mwangi, J. Promoting Sustainable Environments through Urban Green Spaces: Insights from Kenya. Sustainability 2023, 15, 11873. [Google Scholar] [CrossRef]
  2. Ji, S.W.; Ma, R.F.; Ren, L.Y.; Wang, C.J. How to Find Vacant Green Space in the Process of Urban Park Planning: Case Study in Ningbo (China). Int. J. Environ. Res. Public Health 2020, 17, 8282. [Google Scholar] [CrossRef] [PubMed]
  3. Evers, M. Integrative River Basin Management: Challenges and Methodologies Within the German Planning System. Environ. Earth Sci. 2016, 75, 1085. [Google Scholar] [CrossRef]
  4. Eshetu, S.B.; Yeshitela, K.; Sieber, S. Urban Green Space Planning, Policy Implementation, and Challenges: The Case of Addis Ababa. Sustainability 2021, 13, 11344. [Google Scholar] [CrossRef]
  5. Feltynowski, M. Urban Green Spaces in Land-Use Policy—Types of Data, Sources of Data and Staff—The Case of Poland. Land Use Policy 2023, 127, 106570. [Google Scholar] [CrossRef]
  6. Buckland, M.; Pojani, D. Green Space Accessibility in Europe: A Comparative Study of Five Major Cities. Eur. Plan. Stud. 2022, 31, 146–167. [Google Scholar] [CrossRef]
  7. Cernicova-Buca, M.; Gherhes, V.; Obrad, C. Residents’ Satisfaction with Green Spaces and Daily Life in Small Urban Settings: Romanian Perspectives. Land 2023, 12, 689. [Google Scholar] [CrossRef]
  8. Sathyakumar, V.; Raaj, R.; Bardhan, R. Geospatial Approach for Assessing Spatiotemporal Dynamics of Urban Green Space Distribution Among Neighbourhoods: A Demonstration in Mumbai. Urban For. Urban Green. 2020, 48, 126585. [Google Scholar] [CrossRef]
  9. Chen, C.Y.; Bi, L.L.; Zhu, K.F. Study on Spatial-Temporal Change of Urban Green Space in Yangtze River Economic Belt and Its Driving Mechanism. Int. J. Environ. Res. Public Health 2021, 18, 12498. [Google Scholar] [CrossRef]
  10. Pinto, L.V.; Inacio, M.; Ferreira, C.S.S.; Ferreira, A.D.; Pereira, P. Ecosystem Services and Well-Being Dimensions Related to Urban Green Spaces-A Systematic Review. Sustain. Cities Soc. 2022, 85, 104072. [Google Scholar] [CrossRef]
  11. Lee, A.C.K.; Maheswaran, R. The Health Benefits of Urban Green Spaces: A Review of the Evidence. J. Public Health 2011, 33, 212–222. [Google Scholar] [CrossRef] [PubMed]
  12. Farkas, J.Z.; Hoyk, E.; de Morais, M.B.; Csomos, G. A Systematic Review of Urban Green Space Research over the Last 30 Years: A Bibliometric Analysis. Heliyon 2023, 9, e13406. [Google Scholar] [CrossRef] [PubMed]
  13. Li, C.Y.; Zhang, T.T.; Wang, X.; Lian, Z.F. Site Selection of Urban Parks Based on Fuzzy-Analytic Hierarchy Process (F-AHP): A Case Study of Nanjing, China. Int. J. Environ. Res. Public Health 2022, 19, 13159. [Google Scholar] [CrossRef] [PubMed]
  14. Fasihi, H.; Parizadi, T. Analysis of Spatial Equity and Access to Urban Parks in Ilam, Iran. J. Environ. Manag. 2020, 260, 110122. [Google Scholar] [CrossRef] [PubMed]
  15. Zhang, J.G.; Cheng, Y.Y.; Zhao, B. How To Accurately Identify the Underserved Areas of Peri-Urban Parks? An Integrated Accessibility Indicator. Ecol. Indic. 2021, 122, 107263. [Google Scholar] [CrossRef]
  16. Senik, B.; Uzun, O. A Process Approach to the Open Green Space System Planning. Landsc. Ecol. Eng. 2022, 18, 203–219. [Google Scholar] [CrossRef]
  17. Gelo, D.; Turpie, J. Bayesian Analysis of Demand for Urban Green Space: A Contingent Valuation of Developing a New Urban Park. Land Use Policy 2021, 109, 105623. [Google Scholar] [CrossRef]
  18. Halecki, W.; Stachura, T.; Fudala, W.; Stec, A.; Kubon, S. Assessment and Planning of Green Spaces in Urban Parks: A Review. Sustain. Cities Soc. 2023, 88, 104280. [Google Scholar] [CrossRef]
  19. Li, Z.; Bai, X.; Xu, Z.J.; Ma, H.Q.; Xu, Y.A.; Wang, N.; Yue, X. The Optimal Spatial Delineation Method for the Service Level of Urban Park Green Space from the Perspective of Opportunity Equity. Environ. Sci. Pollut. Res. 2023, 30, 85520–85533. [Google Scholar] [CrossRef]
  20. Xu, S.A.; Wang, Y.Z. Influence of Spatial Scale on the Study of Access Fairness of Urban Park Green Space. Front. Environ. Sci. 2023, 10, 1030796. [Google Scholar] [CrossRef]
  21. Doll, C.A.; Burton, M.P.; Pannell, D.J.; Rollins, C.L. Are Greenspaces Too Green? Landscape Preferences and Water Use in Urban Parks. Ecol. Econ. 2023, 211, 107896. [Google Scholar] [CrossRef]
  22. Wang, X.K.; Meng, Q.Y.; Liu, X.Z.; Allam, M.; Zhang, L.L.; Hu, X.L.; Bi, Y.X.; Jancso, T. Evaluation of Fairness of Urban Park Green Space Based on an Improved Supply Model of Green Space: A Case Study of Beijing Central City. Remote Sens. 2023, 15, 244. [Google Scholar] [CrossRef]
  23. Yin, Z.; Zhang, Y.X.; Ma, K.M. Evaluation of PM2.5 Retention Capacity and Structural Optimization of Urban Park Green Spaces in Beijing. Forests 2022, 13, 415. [Google Scholar] [CrossRef]
  24. Biernacka, M.; Kronenberg, J.; Laszkiewicz, E.; Czembrowski, P.; Parsa, V.A.; Sikorska, D. Beyond Urban Parks: Mapping Informal Green Spaces in an Urban-Peri-Urban Gradient. Land Use Policy 2023, 131, 106746. [Google Scholar] [CrossRef]
  25. Engstrom, G.; Gren, A. Capturing the Value of Green Space in Urban Parks in A Sustainable Urban Planning and Design Context: Pros and Cons of Hedonic Pricing. Ecol. Soc. 2017, 22, 21. [Google Scholar] [CrossRef]
  26. Ayele, B.Y.; Megento, T.L.; Habetemariam, K.Y. The Governance and Management of Green Spaces in Addis Ababa, Ethiopia. Heliyon 2022, 8, e09413. [Google Scholar] [CrossRef] [PubMed]
  27. Molle, F. River-Basin Planning and Management: The Social Life of a Concept. Geoforum 2009, 40, 484–494. [Google Scholar] [CrossRef]
  28. Antwi, S.H.; Linnane, S.; Getty, D.; Rolston, A. River Basin Management Planning in the Republic of Ireland: Past, Present and the Future. Water 2021, 13, 2074. [Google Scholar] [CrossRef]
  29. Suhardiman, D.; Bastakoti, R.C.; Karki, E.; Bharati, L. The Politics of River Basin Planning and State Transformation Processes in Nepal. Geoforum 2018, 96, 70–76. [Google Scholar] [CrossRef]
  30. Safiolea, E.; Baki, S.; Makropoulos, C.; Deliège, J.F.; Magermans, P.; Everbecq, E.; Gkesouli, A.; Stamou, A.; Mimikou, M. Integrated Modelling for River Basin Management Planning. Proc. Inst. Civ. Eng.-Water Manag. 2011, 164, 405–419. [Google Scholar] [CrossRef]
  31. Merriam, E.R.; Petty, J.T.; Strager, M.P. Watershed Planning within a Quantitative Scenario Analysis Framework. JOVE-J. Vis. Exp. 2016, 113, e54095. [Google Scholar] [CrossRef]
  32. Matt, J.E.; Underwood, K.L.; Diehl, R.M.; Lawson, K.S.; Worley, L.C.; Rizzo, D.M. Terrain-Derived Measures for Basin Conservation and Restoration Planning. River Res. Appl. 2023, 39, 1795–1811. [Google Scholar] [CrossRef]
  33. Dutta, S.; Sarkar, S. Canal-Oriented Development: Integrating an Urban Canal Front with the City. Territ. Ric. Insediamenti Ambiente 2020, 13, 47–65. [Google Scholar] [CrossRef]
  34. Chen, Y. The Hot Spots and Frontiers of Research on the Grand Canal Culture Belt in China: Literature and Academic Trends. Humanit. Soc. Sci. Commun. 2022, 9, 453. [Google Scholar] [CrossRef]
  35. Yang, J.; Wang, L.; Wei, S. Spatial Variation and Its Local Influencing Factors of Intangible Cultural Heritage Development along the Grand Canal in China. Int. J. Environ. Res. Public Health 2023, 20, 662. [Google Scholar] [CrossRef]
  36. Li, L. Cultural Communication and Diversity along the Grand Canal of China: A Case Study of Folk Songs in Intangible Cultural Heritage. Herit. Sci. 2023, 11, 66. [Google Scholar] [CrossRef]
  37. Liu, J.G.; Tan, X.M.; Deng, J.; Zhu, Y.F. Protection and Utilization Strategies for Grand Canal Heritage. In Proceedings of the 35th IAHR World Congress, Chengdu, China, 8–13 September 2013; Volume I–II, pp. 771–779. [Google Scholar]
  38. Cai, J.D.; Peng, J. Introduction of Beijing-Hangzhou Grand Canal and Analysis of Its Heritage Values. In Water Projects and Technologies in Asia; CRC Press: Boca Raton, FL, USA, 2019; Volume 26, pp. 2–7. [Google Scholar] [CrossRef]
  39. Li, S.H.; Wang, C.S.; Fu, X.X. Exploring the Cultural Heritage Space Adaptability of the Beijing-Hangzhou Grand Canal Based on Point of Interest Data. Landsc. Res. 2023, 48, 33–44. [Google Scholar] [CrossRef]
  40. Xu, H. Heritage Landscape Structure Analysis in Surrounding Environment of The Grand Canal Yangzhou Section. In Proceedings of the 3rd International Conference on Energy Equipment Science and Engineering, Beijing, China, 28–31 December 2017; Volume 128, p. 012042. [Google Scholar] [CrossRef]
  41. Zhang, Y.X.; Zhang, C.Y.; Zhang, X.D.; Wang, X.G.; Liu, T.; Li, Z.; Lin, Q.Y.; Jing, Z.H.; Wang, X.Y.; Huang, Q.Y.; et al. Habitat Quality Assessment and Ecological Risks Prediction: An Analysis in the Beijing-Hangzhou Grand Canal (Suzhou Section). Water 2022, 14, 2602. [Google Scholar] [CrossRef]
  42. Yang, D.Z.; Song, W. Ecological Function Regionalization of the Core Area of the Beijing-Hangzhou Grand Canal Based on the Leading Ecological Function Perspective. Ecol. Indic. 2022, 142, 109247. [Google Scholar] [CrossRef]
  43. Zhang, Q.; Wang, X.M.; Zhu, J.Q.; Li, Z.; Wang, Y. Occurrence and Risk Assessment of Persistent Organic Pollutants in A Branch of the Grand Canal in Hangzhou, China. Environ. Monit. Assess. 2018, 190, 211. [Google Scholar] [CrossRef]
  44. Zhou, S.B.; Ye, J.Y.; Li, J.X.; Zhang, G.Q.; Duan, Y.Q. Identifying Intrinsic Drivers to Changes in Riparian Ecosystem Services By Using Psr Framework: A Case Study of the Grand Canal in Jiangsu, China. Environ. Dev. 2022, 43, 100728. [Google Scholar] [CrossRef]
  45. Lu, L.; Jiao, M.; Weng, L.S. Influence of First-Time Visitors’ Perceptions of Destination Image on Perceived Value and Destination Loyalty: A Case Study of Grand Canal Forest Park, Beijing. Forests 2023, 14, 504. [Google Scholar] [CrossRef]
  46. Zhang, S.Y.; Liu, J.M.; Pei, T.; Chan, C.S.; Gao, C.X.; Meng, B. Perception in Cultural Heritage Tourism: An Analysis of Tourists to the Beijing-Hangzhou Grand Canal, China. J. Tour. Cult. Chang. 2023, 21, 569–591. [Google Scholar] [CrossRef]
  47. Zhang, S.Y.; Liu, J.M.; Pei, T.; Chan, C.S.; Wang, M.D.; Meng, B. Tourism Value Assessment of Linear Cultural Heritage: The Case of the Beijing-Hangzhou Grand Canal in China. Curr. Issues Tour. 2023, 26, 47–69. [Google Scholar] [CrossRef]
  48. Zhang, F.; Yang, L.S.; Leo, X. A Study on the Distribution and Utilization of Recreational Resources along the Grand Canal Culture Belt. Chin. J. Urban Environ. Stud. 2019, 7, 1950015. [Google Scholar] [CrossRef]
  49. He, Y.H.; Wu, L. Analysis on Spatial Development Mode of Eco-Sports Tourism in Grand Canal Landscape Environment Culture Belt. Environ. Monit. Assess. 2022, 194, 925. [Google Scholar] [CrossRef] [PubMed]
  50. Geng, L.P. The Construction of Public Leisure Space of Beijing-Hangzhou Grand Canal Hangzhou Section. Adv. Build. Mater. Sustain. Archit. 2012, 174–177, 2289–2292. [Google Scholar] [CrossRef]
  51. Wang, J.Y.; Wang, M.H.; Dou, H.H.; Su, M.M.; Dong, H.Y.; Liu, Z.H. Research on Climate Change and Water Heritage Tourism Based on the Adaptation Theory-A Case Study of the Grand Canal (Beijing Section). Sustainability 2023, 15, 7630. [Google Scholar] [CrossRef]
  52. Zhang, Y.; Tian, Q.; Wu, J. Coupling Coordination Degree and Obstacle Factors Between the Tourism Industry and Ecological Environment in the Beijing-Hangzhou Grand Canal Basin, China. Environ. Dev. Sustain. 2023. [Google Scholar] [CrossRef]
  53. Chen, M.; Wang, J.C.; Sun, J.; Ye, F.; Zhang, H.Y. Spatio-Temporal Distribution Characteristics of Intangible Cultural Heritage and Tourism Response in the Beijing-Hangzhou Grand Canal Basin in China. Sustainability 2023, 15, 10348. [Google Scholar] [CrossRef]
  54. Shi, F.; Lu, Y.Y.; Chen, L.G.; Hsu, W.L. Evaluation of the Sustainable Use of Land Resources in the Cities along the Jiangsu Section of the Beijing-Hangzhou Grand Canal. Land 2023, 12, 1173. [Google Scholar] [CrossRef]
  55. Mukherjee, S.; Shashtri, S.; Singh, C.K.; Srivastava, P.K.; Gupta, M. Effect of Canal on Land Use/Land Cover Using Remote Sensing and GIS. J. Indian Soc. Remote Sens. 2009, 37, 527–537. [Google Scholar] [CrossRef]
  56. Dale, V.H.; Brown, S.; Calderon, M.O.; Montoya, A.S.; Martinez, R.E. Projected Land-Use Change for the Eastern Panama Canal Watershed and Its Potential Impact. In The Río Chagres, Panama. A Multidisciplinary Profile of a Tropical Watershed; Springer: Dordrecht, The Netherlands, 2005; Volume 52, pp. 337–345. [Google Scholar]
  57. Huang, W.J.; Xi, M.W.; Lu, S.B.; Taghizadeh-Hesary, F. Rise and Fall of the Grand Canal in the Ancient Kaifeng City of China: Role of the Grand Canal and Water Supply in Urban and Regional Development. Water 2021, 13, 1932. [Google Scholar] [CrossRef]
  58. Zhu, C.H.; Nie, Y.P. Study on the Effects of Grand Canal on City Pattern Change of Hangzhou Based on Remote Sensing. In Proceedings of the 2009 Joint Urban Remote Sensing Event, Shanghai, China, 20–22 May 2009; Volume 1–3, pp. 1229–1234. [Google Scholar]
  59. Zhao, Y.; Yan, J.W.; Li, Y.; Bian, G.M.; Du, Y.Z. In-Site Phenotype of the Settlement Space along China’s Grand Canal Tianjin Section: GIS-sDNA-Based Model Analysis. Buildings 2022, 12, 394. [Google Scholar] [CrossRef]
  60. Huo, X.L.; Xu, X.W.; Tang, Y.; Zhang, Z. An Analysis of the Spatial Evolution and Influencing Factors of Rural Settlements along the Shandong Section of the Grand Canal of China. River Res. Appl. 2021, 39, 1283–1299. [Google Scholar] [CrossRef]
  61. Portela, J.F.; Roda, E.M.M. Rural Development Programs in the Riverside Municipalities of the Canal De Castilla: Their Impact on the Evolution of the Population (2000–2020). Bol. Asoc. Geogr. Esp. 2022, 93. [Google Scholar] [CrossRef]
  62. Tsung, N.; Corotis, R.; Chinowsky, P.; Amadei, B. A Retrospective Approach to Assessing the Sustainability of the Grand Canal of China. Struct. Infrastruct. Eng. 2013, 9, 297–316. [Google Scholar] [CrossRef]
  63. Gangopadhyay, A.; Patra, P. The Historical Background of the Canal System in Calcutta, India, and Its Contribution to Development. Proc. Inst. Civ. Eng.-Eng. Hist. Herit. 2020, 173, 80–91. [Google Scholar] [CrossRef]
  64. Cheng, Z.; He, J.L.; Li, Y.; Zhu, Y.X.; Dai, J.C. Coupling-Coordinated Development of the Water-Economy-Innovation Nexus: A Case Study of the Grand Canal Area in China. J. Water Resour. Plan. Manag. 2023, 149, 05023016. [Google Scholar] [CrossRef]
  65. Xia, Y.K. Impact of Green Space on Residents’ Wellbeing: A Case Study of the Grand Canal (Hangzhou Section). Front. Public Health 2023, 11, 1146892. [Google Scholar] [CrossRef]
  66. Chen, H.; Zhao, S.; Zhang, P.; Zhou, Y.; Li, K. Dynamics and Driving Mechanism of Real Estate in China’s Small Cities: A Case Study of Gansu Province. Buildings 2022, 12, 1512. [Google Scholar] [CrossRef]
  67. Watson, S.I. Efficient design of geographically-defined clusters with spatial autocorrelation. J. Appl. Stat. 2021, 49, 3300–3318. [Google Scholar] [CrossRef] [PubMed]
  68. Zhang, P.; Li, W.; Zhao, K.; Zhao, S. Spatial Pattern and Driving Mechanism of Urban–Rural Income Gap in Gansu Province of China. Land 2021, 10, 1002. [Google Scholar] [CrossRef]
  69. Zhang, P.; Chen, H.; Zhao, K.; Zhao, S.; Li, W. Dynamics, Risk and Management Performance of Urban Real Estate Inventory in Yangtze River Delta. Buildings 2022, 12, 2140. [Google Scholar] [CrossRef]
  70. Zhao, S.; Zhao, K.; Yan, Y.; Zhu, K.; Guan, C. Spatio-Temporal Evolution Characteristics and Influencing Factors of Urban Service-Industry Land in China. Land 2022, 11, 13. [Google Scholar] [CrossRef]
  71. Rustamov, J.; Rustamov, Z.; Zaki, N. Green Space Quality Analysis Using Machine Learning Approaches. Sustainability 2023, 15, 7782. [Google Scholar] [CrossRef]
  72. Yang, Z.Q.; Fang, C.L.; Li, G.D.; Mu, X.F. Integrating Multiple Semantics Data to Assess the Dynamic Change of Urban Green Space in Beijing, China. Int. J. Appl. Earth Obs. Geoinf. 2021, 103, 102479. [Google Scholar] [CrossRef]
  73. Wu, C.; Du, Y.H.; Li, S.; Liu, P.Y.; Ye, X.Y. Does Visual Contact with Green Space Impact Housing Prices? An Integrated Approach of Machine Learning and Hedonic Modeling Based on the Perception of Green Space. Land Use Policy 2022, 115, 106048. [Google Scholar] [CrossRef]
  74. Walker, K.W.; Jiang, Z.H. Application of adaptive boosting (AdaBoost) in demand-driven acquisition (DDA) prediction: A machine-learning approach. J. Acad. Libr. 2019, 45, 203–212. [Google Scholar] [CrossRef]
  75. Benocci, R.; Afify, A.; Potenza, A.; Roman, H.E.; Zambon, G. Toward the Definition of a Soundscape Ranking Index (Sri) in an Urban Park Using Machine Learning Techniques. Sensors 2023, 23, 4797. [Google Scholar] [CrossRef]
  76. Kim, D.; Jin, J. Does Happiness Data Say Urban Parks Are Worth It? Landsc. Urban Plan. 2018, 178, 1–11. [Google Scholar] [CrossRef]
  77. Zhang, W.J.; Li, J.K. A Quasi-Experimental Analysis on the Causal Effects of Covid-19 on Urban Park Visits: The Role of Park Features and the Surrounding Built Environment. Urban For. Urban Green. 2023, 82, 127898. [Google Scholar] [CrossRef] [PubMed]
  78. Perez-Martinez, J.; Hernandez-Gil, F.; San Miguel, G.; Ruiz, D.; Arredondo, M.T. Analysing Associations Between Digitalization and the Accomplishment of the Sustainable Development Goals. Sci. Total Environ. 2023, 857, 159700. [Google Scholar] [CrossRef] [PubMed]
  79. Zhang, P.; Li, W.; Zhao, K.; Zhao, Y.; Chen, H.; Zhao, S. The Impact Factors and Management Policy of Digital Village Development: A Case Study of Gansu Province, China. Land 2023, 12, 616. [Google Scholar] [CrossRef]
  80. Yang, L.; Li, Y.T. Analysis on the Spatial Differentiation of National Wetland Park and Influencing Factors in Hunan Province. J. Cent. South Univ. For. Technol. 2019, 13, 99–106. [Google Scholar] [CrossRef]
  81. Li, L.; Zhao, K.; Wang, X.; Zhao, S.; Liu, X.; Li, W. Spatio-Temporal Evolution and Driving Mechanism of Urbanization in Small Cities: Case Study from Guangxi. Land 2022, 11, 415. [Google Scholar] [CrossRef]
  82. Zhao, S.; Yan, Y.; Han, J. Industrial Land Change in Chinese Silk Road Cities and Its Influence on Environments. Land 2021, 10, 806. [Google Scholar] [CrossRef]
  83. Dang, Y.; Wang, C.J.; Chen, P.R. Identification and Optimization Strategy of Urban Park Service Areas Based on Accessibility by Public Transport: Beijing as a Case Study. Sustainability 2022, 14, 7112. [Google Scholar] [CrossRef]
  84. Hasan, M.; Hassan, L.; Al Mamun, A.; Abualreesh, M.H.; Idris, M.H.; Kamal, A.M. Urban Green Space Mediates Spatiotemporal Variation in Land Surface Temperature: A Case Study of an Urbanized City, Bangladesh. Environ. Sci. Pollut. Res. 2022, 29, 36376–36391. [Google Scholar] [CrossRef]
  85. Mu, B.; Liu, C.; Tian, G.H.; Xu, Y.Q.; Zhang, Y.L.; Mayer, A.L.; Lv, R.; He, R.Z.; Kim, G. Conceptual Planning of Urban-Rural Green Space from a Multidimensional Perspective: A Case Study of Zhengzhou, China. Sustainability 2020, 12, 2863. [Google Scholar] [CrossRef]
  86. Heshani, A.L.S.; Winijkul, E. Numerical Simulations of the Effects of Green Infrastructure on PM2.5 Dispersion in an Urban Park in Bangkok, Thailand. Heliyon 2022, 8, e10475. [Google Scholar] [CrossRef] [PubMed]
  87. Xing, Z.; Zhao, S.; Li, K. Evolution Pattern and Spatial Mismatch of Urban Greenspace and Its Impact Mechanism: Evidence from Parkland of Hunan Province. Land 2023, 12, 2071. [Google Scholar] [CrossRef]
  88. You, Z.; Feng, Z.; Yang, Y. Relief Degree of Land Surface Dataset of China (1km). Digital Journal of Global Change Data Repository. J. Glob. Chang. Data Discov. 2018, 2, 151–155. [Google Scholar] [CrossRef]
  89. Broner, F.; Bustos, P.; Carvalho, V.M. Sources of Comparative Advantage in Polluting Industries; NBER Working Paper; National Bureau of Economic Research: Cambridge, MA, USA, 2012; p. 8337. [Google Scholar] [CrossRef]
  90. Hering, L.; Poncet, S. Environmental policy and exports: Evidence from Chinese cities. J. Environ. Econ. Manag. 2014, 68, 296–318. [Google Scholar] [CrossRef]
  91. Chen, H.; Guo, W.; Feng, X.; Wei, W.; Liu, H.; Feng, Y.; Gong, W. The impact of low-carbon city pilot policy on the total factor productivity of listed enterprises in China. Resour. Conserv. Recycl. 2021, 169, 105457. [Google Scholar] [CrossRef]
  92. Fortheringham, S.; Brunsdon, C.H.; Charlton, M. Quantitative Geography; SAGE Publications: London, UK, 2000. [Google Scholar]
  93. Shrestha, A.; Luo, W. Analysis of Groundwater Nitrate Contamination in the Central Valley: Comparison of the Geodetector Method, Principal Component Analysis and Geographically Weighted Regression. ISPRS Int. J. Geo-Inf. 2017, 6, 297. [Google Scholar] [CrossRef]
  94. Wang, J.F.; Li, X.H.; Christakos, G.; Liao, Y.L.; Zhang, T.; Gu, X.; Zheng, X.Y. Geographical detectors-based health risk assessment and its application in the neural tube defects study of the Heshun region, China. Int. J. Geogr. Inf. Sci. 2010, 24, 107–127. [Google Scholar] [CrossRef]
  95. Zhao, S.; Li, W.; Zhao, K.; Zhang, P. Change Characteristics and Multilevel Influencing Factors of Real Estate Inventory—Case Studies from 35 Key Cities in China. Land 2021, 10, 928. [Google Scholar] [CrossRef]
  96. Wang, J.F.; Hu, Y. Environmental health risk detection with GeogDetector. Environ. Model. Softw. 2012, 33, 114–115. [Google Scholar] [CrossRef]
  97. Melchiorri, M.; Pesaresi, M.; Florczyk, A.J.; Corbane, C.; Kemper, T. Principles and Applications of the Global Human Settlement Layer as Baseline for the Land Use Efficiency Indicator-SDG 11.3.1. ISPRS Int. J. Geo-Inf. 2019, 8, 96. [Google Scholar] [CrossRef]
  98. Laituri, M.; Davis, D.; Sternlieb, F.; Galvin, K. SDG Indicator 11.3.1 and Secondary Cities: An Analysis and Assessment. ISPRS Int. J. Geo-Inf. 2021, 10, 713. [Google Scholar] [CrossRef]
  99. Rigolon, A. A Complex Landscape of Inequity in Access to Urban Parks: A Literature Review. Landsc. Urban Plan. 2016, 153, 160–169. [Google Scholar] [CrossRef]
  100. Ren, X.Y.; Guan, C.H. Evaluating Geographic and Social Inequity of Urban Parks in Shanghai Through Mobile Phone-Derived Human Activities. Urban For. Urban Green. 2022, 76, 127709. [Google Scholar] [CrossRef]
  101. Tan, C.D.; Tang, Y.H.; Wu, X.F. Evaluation of the Equity of Urban Park Green Space Based on Population Data Spatialization: A Case Study of a Central Area of Wuhan, China. Sensors 2019, 19, 2929. [Google Scholar] [CrossRef] [PubMed]
  102. Gao, W.X.; Lyu, Q.; Fan, X.; Yang, X.C.; Liu, J.T.; Zhang, X.R. Building-Based Analysis of the Spatial Provision of Urban Parks in Shenzhen, China. Int. J. Environ. Res. Public Health 2017, 14, 1521. [Google Scholar] [CrossRef]
  103. Zhu, J.Y.; Lu, H.T.; Zheng, T.C.; Rong, Y.J.; Wang, C.X.; Zhang, W.; Yan, Y.; Tang, L. Vitality of Urban Parks and Its Influencing Factors from the Perspective of Recreational Service Supply, Demand, and Spatial Links. Int. J. Environ. Res. Public Health 2020, 17, 1615. [Google Scholar] [CrossRef] [PubMed]
  104. Wu, L.F.; Kim, S.K. Exploring the Equality of Accessing Urban Green Spaces: A Comparative Study of 341 Chinese Cities. Ecol. Indic. 2021, 121, 107080. [Google Scholar] [CrossRef]
  105. Choumert, J.; Cormier, L. The Provision of Urban Parks: An Empirical Test of Spatial Spillovers in an Urban Area Using Geographic Information Systems. Ann. Reg. Sci. 2011, 47, 437–450. [Google Scholar] [CrossRef]
  106. Kim, J.; Choi, N.; Lee, D.K. Spatial Preference Heterogeneity in Policies for Improving Urban Green Spaces. Urban For. Urban Green. 2022, 78, 127781. [Google Scholar] [CrossRef]
  107. Cruickshank, J. Local Hegemonies Resisting a Green Shift and What to Do About It: The Introduction of a Regional Park in Southern Norway. J. Environ. Policy Plan. 2018, 20, 313–327. [Google Scholar] [CrossRef]
  108. Branca, D.; Haller, A.; Mossa, M. From the mountains to the sea: The Tepilora Natural Regional Park, Sardinia. eco.mont-J. Prot. Mt. Areas Res. Manag. 2023, 15, 28–34. [Google Scholar] [CrossRef]
  109. Chen, J.X.; Han, S.S.; Chen, S.Q. Understanding the Structure and Complexity of Regional Greenway Governance in China. Int. Dev. Plan. Rev. 2022, 44, 241–264. [Google Scholar] [CrossRef]
  110. Cheng, Y.; Shi, Y.; Andrew, S. Exploring the Link Between Fiscal Arrangements and the Quality of Public Services: Evidence from Major Us Urban Park Systems. Public Perform. Manag. Rev. 2020, 43, 1445–1470. [Google Scholar] [CrossRef]
  111. Nam, J.; Dempsey, N. Place-Keeping for Health? Charting the Challenges for Urban Park Management in Practice. Sustainability 2019, 11, 4383. [Google Scholar] [CrossRef]
  112. Smith, A.; Whitten, M.; Ernwein, M. De-municipalisation? Legacies of Austerity for England’s Urban Parks. Geogr. J. 2023. [Google Scholar] [CrossRef]
  113. Luo, T.; Yang, F.M.; Wu, L.L.; Gao, X.H. Equity Evaluation of Urban Park System: A Case Study of Xiamen, China. J. Environ. Eng. Landsc. Manag. 2020, 28, 125–136. [Google Scholar] [CrossRef]
  114. Feng, S.; Chen, L.D.; Sun, R.H.; Feng, Z.Q.; Li, J.R.; Khan, M.S.; Jing, Y.C. The Distribution and Accessibility of Urban Parks in Beijing, China: Implications of Social Equity. Int. J. Environ. Res. Public Health 2019, 16, 4894. [Google Scholar] [CrossRef]
  115. Kim, K.; Lee, C.K.; Kim, H.W. Understanding the Accessibility of Urban Parks and Connectivity of Green Spaces in Single-Person Household Distribution: Case Study of Incheon, South Korea. Land 2022, 11, 1441. [Google Scholar] [CrossRef]
  116. Guo, S.H.; Song, C.; Pei, T.; Liu, Y.X.; Ma, T.; Du, Y.Y.; Chen, J.; Fan, Z.D.; Tang, X.L.; Peng, Y.; et al. Accessibility to Urban Parks for Elderly Residents: Perspectives from Mobile Phone Data. Landsc. Urban Plan. 2019, 191, 103642. [Google Scholar] [CrossRef]
  117. Guo, S.H.; Yang, G.G.; Pei, T.; Ma, T.; Song, C.; Shu, H.; Du, Y.Y.; Zhou, C.H. Analysis of Factors Affecting Urban Park Service Area in Beijing: Perspectives from Multi-Source Geographic Data. Landsc. Urban Plan. 2019, 181, 103–117. [Google Scholar] [CrossRef]
Figure 1. The Grand Canal and its location in China.
Figure 1. The Grand Canal and its location in China.
Land 13 00042 g001
Figure 2. Technical route for the study of urban park green areas within the Grand Canal.
Figure 2. Technical route for the study of urban park green areas within the Grand Canal.
Land 13 00042 g002
Figure 3. Relative share and growth rate spatial analysis of the supply scale of UPGAs in the Grand Canal of China.
Figure 3. Relative share and growth rate spatial analysis of the supply scale of UPGAs in the Grand Canal of China.
Land 13 00042 g003
Figure 4. Spatiotemporal evolution model of the supply scale of urban park green areas on the Grand Canal of China.
Figure 4. Spatiotemporal evolution model of the supply scale of urban park green areas on the Grand Canal of China.
Land 13 00042 g004
Figure 5. Relative share and growth rate spatial analysis of the per-capita area of UPGAs in the Grand Canal of China.
Figure 5. Relative share and growth rate spatial analysis of the per-capita area of UPGAs in the Grand Canal of China.
Land 13 00042 g005
Figure 6. Spatiotemporal evolution model of the Ppr-capita area of urban park green areas on the Grand Canal of China.
Figure 6. Spatiotemporal evolution model of the Ppr-capita area of urban park green areas on the Grand Canal of China.
Land 13 00042 g006
Figure 7. Factor importance analysis of UPGAs in the Grand Canal of China.
Figure 7. Factor importance analysis of UPGAs in the Grand Canal of China.
Land 13 00042 g007
Figure 8. Factor spatial effect analysis of the supply scale of urban park green areas on the Grand Canal of China.
Figure 8. Factor spatial effect analysis of the supply scale of urban park green areas on the Grand Canal of China.
Land 13 00042 g008
Figure 9. Factor spatial effect analysis of the per-capita area of urban park green areas on the Grand Canal of China.
Figure 9. Factor spatial effect analysis of the per-capita area of urban park green areas on the Grand Canal of China.
Land 13 00042 g009
Figure 10. Performance analysis of urban park green areas on the Grand Canal of China.
Figure 10. Performance analysis of urban park green areas on the Grand Canal of China.
Land 13 00042 g010
Figure 11. Factor spatial effect analysis of the urban park green areas on the Grand Canal of China.
Figure 11. Factor spatial effect analysis of the urban park green areas on the Grand Canal of China.
Land 13 00042 g011
Table 1. Literature review of the Grand Canal of China.
Table 1. Literature review of the Grand Canal of China.
AreasViewpoints
Cultural Heritage and Ecological ProtectionI. Study the distribution, characteristics, and influencing factors of historical relics and intangible cultural heritage along the canal [35,36] and further propose strategies for protection and utilization [37]. II. Assess the value of cultural heritage and relics along the canal and determine the adaptive landscape development methods [38,39,40]. III. Emphasize the evaluation of canal habitat quality [41], ecological functions [42], and pollution risks [43] and analyze their impact on ecosystem services [44].
Tourism developmentI. Value the construction of tourism destinations and the construction of the tourism industry system, including the image perception of tourism destinations and its impact on tourism loyalty [45,46], tourism value assessment and resource utilization [47,48], tourism spatial development models [49], and regional tourism openness and cooperation [50]. II. Analyze the coupling relationship between tourism and ecology, heritage, and climate, including the impact of climate change on the development of canal tourism [51], the collaboration between tourism and ecosystems and their development obstacles [52], and the correlation between the spatio–temporal distribution of cultural heritage and tourism response [53].
Land use and Urban-rural changesI. Analyze the level of sustainable and healthy land use along the canal [54] and land use/cover changes [55] and their impact on regional development [56]. II. Analyze the rise and fall of cities along the canal and spatial pattern and structural changes and their influencing factors, especially the role of canal logistics and flooding [57,58]. III. Analyze the geographical evolution of rural spatial settlements along the canal and its influencing factors, especially traditional villages and historical and cultural ancient villages [59,60,61].
Sustainable DevelopmentI. Assess the spatial sustainable development of the canal basin [62] and its contribution to regional development [63]. II. Analyze the spatio–temporal characteristics of urbanization and the socio–economic benefits of canal land using the coupled coordination degree model to identify the synergistic development model of water–economy–innovation [64].
Table 2. Indicator selection of independent variables.
Table 2. Indicator selection of independent variables.
IndicatorCodeVIF
Supply scale of UPGAs Y 1 --
Per-capita area of UPGAs Y 2 --
SocietyPopulation density X 1 1.14
Proportion of population aged 60 and above X 2 1.93
Outflow population X 3 3.49
EconomicGDP X 4 2.84
Per-capita GDP X 5 8.25
Fiscal self-sufficiency rate X 6 5.41
NaturalTopography X 7 1.27
Average temperature X 8 1.58
Ventilation coefficient X 9 1.63
Table 3. Matching relationship between urban park green space supply and population demand based on LCRPGR measurements.
Table 3. Matching relationship between urban park green space supply and population demand based on LCRPGR measurements.
TypeLCRPGRLCRPGR
Super oversupply>0<0<0
Super undersupply<0>0<0
Negative oversupply<0<0>0 and <0.75
Negative undersupply<0<0≥1.25
Negative equilibrium<0<0>0.75 and <1.25
Positive oversupply>0>0≥1.25
Positive undersupply>0>0>0.75 and <1.25
Positive equilibrium>0>0>0 and <0.75
Table 4. Machine learning regression results based on different algorithm analysis of urban park green areas on the Grand Canal of China. The asterisk represents the result after rounding.
Table 4. Machine learning regression results based on different algorithm analysis of urban park green areas on the Grand Canal of China. The asterisk represents the result after rounding.
FactorsSupply Scale of UPGAsPer-Capita Area of UPGAs
Decision TreeRandom ForestAdaboostExtra TreesDecision TreeRandom ForestAdaboostExtra Trees
X 1 13.00%5.90%6.70%7.60%24.70%13.50%12.00%12.80%
X 2 4.80%7.30%8.80%5.80%7.10%6.70%7.90%10.10%
X 3 3.30%7.40%11.40%8.70%19.50%18.00%16.90%13.50%
X 4 3.10%5.50%6.70%6.00%3.00%7.70%10.20%8.40%
X 5 57.40%47.50%29.60%28.10%4.60%4.90%8.20%10.10%
X 6 0.00%5.00%6.80%10.50%0.90%6.30%7.70%10.80%
X 7 7.40%10.50%13.70%15.00%33.00%22.00%18.60%13.10%
X 8 0.00%3.60%6.80%7.60%4.20%6.00%5.50%6.70%
X 9 10.80%7.20%9.50%10.70%3.00%14.80%13.00%14.60%
R21 *0.891 *0.901 *0.870.990.86
Table 5. Descriptive statistical analysis of GWR parameters for the supply scale of urban park green areas on the Grand Canal of China.
Table 5. Descriptive statistical analysis of GWR parameters for the supply scale of urban park green areas on the Grand Canal of China.
FactorsMin25% QuantileMedian75% QuantileMax
Population density X 1 −0.00480.00430.01050.02410.0394
Proportion of population aged 60 and above X 2 0.01930.04530.05410.05990.0801
Outflow population X 3 −0.0272−0.0220−0.00490.02650.0320
GDP X 4 0.08530.11410.19140.21620.2379
Per-capita GDP X 5 −0.1189−0.0948−0.0622−0.02530.0030
Fiscal self-sufficiency rate X 6 0.02670.05340.06420.11810.1432
Topography X 7 −0.1802−0.1270−0.0838−0.0478−0.0358
Average temperature X 8 −0.1571−0.1136−0.1051−0.0777−0.0576
Ventilation coefficient X 9 −0.1224−0.0918−0.0458−0.0093−0.0043
Table 6. Descriptive statistical analysis of GWR parameters for the per-capita area of urban park green areas on the Grand Canal of China.
Table 6. Descriptive statistical analysis of GWR parameters for the per-capita area of urban park green areas on the Grand Canal of China.
FactorsMin25% QuantileMedian75% QuantileMax
Population density X 1 0.48710.51630.51950.52130.5392
Proportion of population aged 60 and above X 2 −0.0669−0.0463−0.0303−0.01640.0017
Outflow population X 3 −0.01880.00080.03350.06710.0815
GDP X 4 −0.1529−0.1385−0.1140−0.0711−0.0506
Per-capita GDP X 5 −0.0924−0.0607−0.0462−0.0397−0.0226
Fiscal self-sufficiency rate X 6 0.06670.09510.13200.15790.1921
Topography X 7 −0.0645−0.03690.01470.06830.0847
Average temperature X 8 −0.0728−0.0352−0.0251−0.00990.0374
Ventilation coefficient X 9 −0.1124−0.1106−0.0939−0.0724−0.0481
Table 7. Factor interactive effect analysis of urban park green areas on the Grand Canal of China.
Table 7. Factor interactive effect analysis of urban park green areas on the Grand Canal of China.
Code X 1 X 2 X 3 X 4 X 5 X 6 X 7 X 8 X 9
Supply scale of UPGAs X 1 0.08
X 2 0.540.12
X 3 0.630.630.28
X 4 0.840.800.840.45
X 5 0.610.670.560.870.26
X 6 0.620.630.650.780.570.25
X 7 0.150.190.430.520.340.330.00
X 8 0.270.350.420.690.380.410.080.05
X 9 0.530.560.690.820.640.760.160.310.12
Per-capita area of UPGAs X 1 0.03
X 2 0.090.05
X 3 0.490.360.19
X 4 0.060.070.250.01
X 5 0.520.340.690.330.20
X 6 0.280.140.380.120.330.06
X 7 0.110.110.470.070.380.140.02
X 8 0.300.240.560.250.620.380.300.18
X 9 0.210.210.550.110.520.170.080.350.03
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

Cai, Z.; Zhao, S.; Huang, M.; Zhang, C. Evolution Model, Mechanism, and Performance of Urban Park Green Areas in the Grand Canal of China. Land 2024, 13, 42. https://doi.org/10.3390/land13010042

AMA Style

Cai Z, Zhao S, Huang M, Zhang C. Evolution Model, Mechanism, and Performance of Urban Park Green Areas in the Grand Canal of China. Land. 2024; 13(1):42. https://doi.org/10.3390/land13010042

Chicago/Turabian Style

Cai, Zihan, Sidong Zhao, Mengshi Huang, and Congguo Zhang. 2024. "Evolution Model, Mechanism, and Performance of Urban Park Green Areas in the Grand Canal of China" Land 13, no. 1: 42. https://doi.org/10.3390/land13010042

APA Style

Cai, Z., Zhao, S., Huang, M., & Zhang, C. (2024). Evolution Model, Mechanism, and Performance of Urban Park Green Areas in the Grand Canal of China. Land, 13(1), 42. https://doi.org/10.3390/land13010042

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