Next Article in Journal
Can Tourism Development Help Improve Urban Liveability? An Examination of the Chinese Case
Previous Article in Journal
What’s Happening with the Patent Box Regimes? A Systematic Review
Previous Article in Special Issue
Land Comprehensive Carrying Capacity of Major Grain-Producing Areas in Northeast China: Spatial–Temporal Evolution, Obstacle Factors and Regulatory Policies
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Demarcation of Urban Development Boundary Based on the Maxent-CA Model: A Case Study of Wuxi in China

School of Geography and Ocean Science, Nanjing University, Nanjing 210023, China
*
Author to whom correspondence should be addressed.
Sustainability 2022, 14(18), 11426; https://doi.org/10.3390/su141811426
Submission received: 11 August 2022 / Revised: 3 September 2022 / Accepted: 7 September 2022 / Published: 12 September 2022
(This article belongs to the Special Issue Sustainable Urban Development and Land Use Policy)

Abstract

:
With the rapid development of urbanization, the demarcation of the urban development boundary (UDB) is of great practical significance to curb the disorderly spread of urban land, avoid losing control of urban development space, and build a barrier to green development space. In this paper, we propose a method to support the demarcation of the UDB by combining the Maxent model and the cellular automata (CA) model. This approach comprehensively considers the relationship between urban construction suitability, neighborhood effect, spatial constraint, and random interference based on a spatio-temporal dynamic simulation. This contributes to the analysis of the driving mechanism and distribution pattern of urban expansion. According to the principle of scale expansion and centralization, the simulation result is modified to demarcate the UDB. The following conclusions are drawn: the Maxent-CA model can intuitively reflect the driving mechanism and accurately simulate urban expansion in specific cities, which contributes to demarcating the UDB. Considering that this method fully embodies the principle of combining top-down and bottom-up approaches in the demarcation of UDB, we argue that the Maxent-CA model is of vital importance for the sustainable development of the living environment and is of great reference value for territorial spatial planning.

1. Introduction

Since the reform and opening up in 1978, social economy development and urbanization in China have experienced rapid growth [1,2]. According to the seventh population census in China, the urbanization rate of the population reached 63.89% in 2020. From 2010 to 2020, the built-up area in Chinese cities increased from 40,058.01 km2 to 60,721.32 km2, an increase of 1.52 times. Comparatively, the urban population increased from 669.78 million to 902.20 million in a decade, which was only an increase of 0.34 times. The rate of urban expansion is much higher than that of population growth. This pattern of urban sprawl is prone to causing serious urban and rural problems, such as the occupation of high-quality cultivated land resources, tough livelihoods for land-lost farmers, a declining trend of industrial development, low land use efficiency, aggravated environmental pollution, and population hollowing issues [3,4,5,6]. Controlling the excessive urban expansion and ensuring the sustainable intensification of cultivated land [7] to achieve a green, low-carbon, efficient and intensive spatial pattern of land use have become one of the hot issues of academic research [8]. However, as the largest developing country, China is still in a period of rapid urbanization, which means that urban expansion is inevitable. Therefore, facing the contradiction between economic development and environmental pressure, China has put forward the concept of the “urban development boundary (UDB), ecological protection red line, and permanent basic farmland” in its territorial space planning to achieve a balance between urban development, ecological protection, and food security [9,10]. Among the three control lines, the UDB is an effective means to curb the disorderly spread of urban land, avoid losing control of urban development space, and build a barrier to green development space. It has become a key task in China’s territorial space planning and a primary measure used to control urban expansion.
The UDB, originally called the urban growth boundary (UGB), is a regional planning tool to protect rural areas and promote the compactness of cities [11,12]. After World War II, the rapid development of car traffic in Western countries resulted in problems such as suburban urbanization, central city recession, and disordered urban sprawl [13]. Thus, various countries began to pay attention to the coordinated development of cities and peripheral space. In 1976, the USA first applied the UGB theory to define the urban development scope in Salem [13]. Subsequently, the USA incorporated the UDB theory as a policy control tool [14] into the smart growth movement [15]. The demarcation of the UDB in the Portland metropolitan area is a typical form of integration of policy and theory [16]. Zhang [17] first introduced the UDB theory into China, but it was not until 2014 that Chinese research on the UDB turned from academic inquiry to practical implementation. The report of the 19th National Congress of the Communist Party of China in 2017 stressed the need to complete the demarcation of three control lines. In the Guiding Opinions regarding Coordinating the Delimitation and Implementation of Three Control Lines in National Spatial Planning in 2019, the General Office of the State Council, once again, emphasized the strict demarcation and implementation of three control lines, hoping to form a more scientific and orderly territorial space management system by 2035. As the maximum rigid boundary for the expansion of urban construction land, the UDB is the demarcation line between developable construction areas and forbidden construction areas. Its application is gradually being recognized and supported by the Chinese government.
The existing demarcation methods for UDB are mainly divided into two types: the reverse method and the forward method [18]. Based on the urban development suitability, the reverse method draws the forbidden development boundary to frame the suitable construction area of urban land [19,20]. This method can visually distinguish between suitable and unsuitable areas for urban construction, but it has a low degree of differentiation for the location conditions. Based on the current distribution of urban construction land, the forward method simulates urban expansion through an analysis of the driving mechanism and then demarcates the UDB. The method places demand on the training and accuracy of the model. Studies have combined the two to form a better model system, which is called the synthesis method [21,22,23]. Based on the forbidden development boundary, the synthesis method establishes a mathematical model to simulate the expansion pattern so as to demarcate the UDB. The key to simulating urban expansion lies in the selection of models, which includes the logistic regression model, artificial neural networks (ANN), the SLEUTH (slope, land use, exclusion, urban extent, transportation, and hillshade) model, the future land use simulation (FLUS) model, ant colony optimization model algorithm [24,25,26,27,28,29], and the cellular automaton (CA) model.
The CA model assigns a value of 0 or 1 to each unitary cell according to its binary state after spatial cellularization. Based on certain state transition rules, the spatial evolution of the next time series is carried out considering the states of other cells that form a specific neighborhood relationship with a single cell [30,31]. This model has the advantages of spatio-temporal dispersion, spatial interaction, and local temporal causality. It can achieve spatial interaction in the adjacent hierarchical structure [32] and effectively integrate the geographical driving factors of urban expansion [33]. Therefore, it has been widely applied in the study of urban expansion and UDB demarcation [34,35,36,37]. However, the CA model has the disadvantage of being unable to obtain model parameters from images of a single period and a single class, so it cannot explain the variable mechanism well. In addition, if only the neighborhood effect is considered in the spatio-temporal dynamic simulation, the simulation result will be far different from the complex process of real urban expansion. Therefore, existing studies on the CA model usually aim to explore state transition rules to improve the simulation accuracy [38]. Combining the CA model with other intelligent algorithms, such as SD-CA [34], SLEUTH-CA [27], Markov-CA [39], and other models, tends to give simulation results that are more similar to real urban expansion.
This paper attempted to combine the CA model and Maxent model to study the spatial simulation of urban expansion. Based on the partial constraint information of the unknown distribution, the Maxent model can select the probability distribution with the maximum information entropy under the constraint condition, which has the most uncertainty from random events [40]. The Maxent model has been widely applied in the field of species habitat assessment [41,42]. This model can automatically obtain the weight coefficients of variables from images of a single period and a single class, which contributes to reflecting the driving mechanism and formulating state transition rules. Combined with the CA model, on the one hand, it can make up for the disadvantage that the CA model cannot obtain the model parameters from images of a single period and a single class. On the other hand, it can give full play to the advantages of both models and greatly improve the simulation accuracy. At present, only a few studies on the Maxent-CA model have proved its practical value. Lin combines the Maxent-CA model and multiperiod dynamic programming algorithm to propose a novel method for the scientific optimization of nature reserves [43]. Tan constructs the SD-Maxent-CA model for dynamic simulation of natural forests to present its potentiality in future land use simulation [44]. However, only a few scholars have used the Maxent-CA model to make empirical studies in urban expansion simulation. Zhang firstly applied the Maxent-CA model to simulate urban expansion in China [45]. The results show that the model performs well in simulating the land use of a single class. However, this study only simulates urban expansion at a large regional scale (China) from a macro perspective without conducting model verification for a specific small regional scale. Additionally, the resolution of this study is 1 km, so the variability and uniqueness of urban expansion in small regions tend to be missed. Chen verified that the Maxent-CA model can simulate a more compact land pattern [46] but did not explore more practical value in the field of UDB demarcation. In view of this, this paper simulated urban land expansion based on the Maxent-CA model to propose a new approach to demarcate the UDB. It can provide a scientific basis and decision reference for the demarcation of the three zones and three lines in China.
This paper selected Wuxi city as the study area and applied the Maxent-CA model to simulate urban expansion in a small region. From the perspective of both natural and socioeconomic conditions, 11 environmental variables were selected. The maximum entropy probability was calculated based on the new range of urban land in Wuxi from 2005 to 2015. Under dual control of the constraint condition and the maximum expansion area in the CA model, this study simulated the range of urban expansion in 2035 and demarcated the UDB of Wuxi using the synthesis method. In this way, we were able to delve into the mechanisms by which relevant environmental variables drive urban expansion as well as explore a more scientific method of demarcating the UDB. This study aimed to provide a decision-making reference for improving the urban construction layout, coordinating urban and rural development, and promoting economical and intensive land use.

2. Materials and Methods

2.1. Study Area

Wuxi (119°33′–120°38′ E, 31°07′–32°02′ N) is located in Jiangsu Province in Eastern China. It is a vital part of the Yangtze River Delta city cluster and contributes to regional development by undertaking economic radiation from Shanghai and driving the industrial upgrade of the hinterland. Wuxi has a total area of 4627.47 km2, including five districts (Liangxi, Binhu, Huishan, Xishan, and Xinwu) and two county-level cities under its jurisdiction (Jiangyin and Yixing) (Figure 1). The topography of Wuxi is mainly plain. The city has a dense network of rivers, with the Yangtze River to the north and Taihu Lake to the south. The humid subtropical monsoon climate in this area gives Wuxi four distinct seasons and sufficient precipitation. The excellent location and climate conditions have led to Wuxi becoming one of the most economically developed cities in China. In 2020, Wuxi’s GDP totaled CNY 123,748 billion with a per capita GDP of CNY 165,777, ranking first in Jiangsu Province and second in China. The urbanization rate of the permanent resident population reached 82.8%.
The rapid development of the economic level and urbanization led to the rapid expansion of construction land in Wuxi. Between 2010 and 2020, the built-up area expanded from 348.50 km2 to 563.41 km2 (an increase of 61.67%), while the permanent resident population grew from 6.38 million to 7.46 million (an increase of only 17.07%), which meant that land urbanization occurred much faster than population urbanization in Wuxi. In terms of the pattern of construction land expansion, although the urban expansion in Wuxi between 1990 and 2005 was dominated by edge expansion and infilling expansion, its small-area patches were obviously reduced [47]. Between 2005 and 2015, the density of urban land patches decreased from 0.44 to 0.22, which means that the pattern of urban expansion became more fragmented. The efficiency of land development needs to be improved urgently.

2.2. Data

(1)
Land use data
Land use data from Wuxi were vector data obtained from the Resource and Environment Science and Data Center (RESDC) of the Chinese Academy of Sciences (CAS) (http://www.resdc.cn/, accessed on 23 December 2021), and the years 2005, 2015, and 2020 were selected as the time nodes. The classification accuracy of the original data was as high as 95% [48], allowing the extraction of the secondary land class of “the built-up area above the county” and “other construction land such as industrial and mining, roads, airports, special land” as urban land in Wuxi. Meanwhile, the “water area” and “forest land” in the primary land class can be extracted as vital components of the urban land-restricted area.
(2)
Resource and environmental data
SRTM digital elevation model (DEM) data from Jiangsu Province (90 m × 90 m) and the spatial distribution dataset of the normalized difference vegetation index (NDVI) (1 km × 1 km, 2015) were obtained from RESDC. The DEM data were used to calculate the index of the surface relief and aspect. A 1:4,000,000 map of soil types in the Yangtze River Delta region (2013) was obtained from the Cold and Arid Regions Environmental and Engineering Research Institute (CAREERI) of the CAS. The soil types were classified into seven categories: iron bauxite, alfisols, primary soil, semi-hydromorphic soil, anthropogenic soil, lake–reservoir, and river. The ecological protection red line and permanent basic farmland data were sourced from the Wuxi Natural Resources Bureau as essential components of the urban land restriction area.
(3)
POI and network data
Using web crawling, we obtained the distribution data for points of interest (POI) and road network data from Wuxi on the Baidu map. The POI data included transportation facilities, medical care, science and culture education, enterprise companies, government agencies, and other categories. Data cleaning was needed to remove irrelevant information to ensure authenticity and unification, and spatial processing was achieved through coordinate transformation. The road network data were in shapefile format and included county roads, provincial roads, national roads, expressways, and high-speed railways. Due to the limitation of multivariate data collection, the year of network big data was unified as 2019.
(4)
Socioeconomic data
Spatial data related to China’s gross domestic product (GDP) (1 km × 1 km, 2015) and population (1 km × 1 km, 2010) were obtained from the RESDC. Permanent resident population data from 2006 to 2020 were taken from the Wuxi Statistical Yearbook.
The coordinate system was set as a Krasovsky ellipsoid. Tools for resampling and extracting by mask were used to unify the resolution (30 m × 30 m) and range size (3234 rows × 3553 columns) of the raster data.

2.3. Method

2.3.1. Index System

Urban expansion is a spatiotemporal evolution process involving the joint action of various natural and human factors [49,50]. This study selected 11 environmental variables as influencing factors of the urban expansion mechanism, considering both natural conditions and socioeconomic conditions. In terms of natural conditions, the attributes of the urban underlying surface affect the construction cost of urban land and the geographical locations of buildings, while the vegetation coverage and soil type affect the selection of urban development land [25]. All of these factors play a role in guiding the spatial layout of urban expansion. Thus, the surface relief, aspect, NDVI, and soil type were selected as evaluation indexes. In terms of socioeconomic conditions, due to the impact of traffic accessibility, urban expansion near major traffic arteries is significantly faster than that in other areas [51,52,53]. The closer the area is to large-scale transport stations, the more prominent the direct demand for urban land is. This indirectly affects infrastructure construction, industrial development, and population agglomeration through the multiplier effect, further expanding the demand for urban land [54]. Therefore, distances to roads (including county roads, provincial roads, national roads, highways, and high-speed railways), transport stations (including airports, ports, railway stations, metro stations, and bus stations), infrastructure (including hospitals, schools, parks, and squares), industrial enterprises, and government departments at all levels were selected. Additionally, the economic level and population are positive independent factors for most types of urban expansion [45,54]. Economic development and population growth can accelerate the transformation of non-urban land, which was characterized by the gross domestic product (GDP) and population distribution in this paper (Table 1, Figure 2).

2.3.2. Suitability Evaluation of Urban Construction Based on the Maxent Model

Maxent (version 3.4.1) (CBC, New York, NY, USA) [55], under the Java runtime environment, is a type of software based on the principle of maximum entropy. According to the known distribution information ( X ) (Figure 2) of the new points of urban land, by taking the environmental variable ( Y ) as the constraint condition, the influence probability of Y on the distribution of X can be calculated. The probability distribution ( X ) with the largest information entropy and the most uniform distribution in the probability model will be obtained [56], which means that the urban land distribution results will be the most stochastic.
H X , Y = i = 1 n p x i , y i log p x i , y i ,
X = a r g m a x H X , Y ,
The 11 environmental variables were set to ASCII format, while 516,542 new urban land points were set to CSV data in the form of latitude and longitude coordinates. Considering the possibility of overfitting, a random sample of the new points was taken before importing them into Maxent software. In order to explore the extent to which different numbers of samples affect the accuracy of the Maxent model, 500, 600, 800, 1000, 2000, 5000, and 10,000 sample points were selected randomly. Subject to the same environmental constraints, a comparative analysis was carried out to determine the optimal number of samples.
To improve the accuracy of the results, the sample data were divided into two parts, 75% of the points were randomly selected for model training and the remaining 25% were selected for the resulting test. The software’s built-in jackknife test can assist with interpreting the contribution and importance of each environmental variable to urban expansion, of which the output file is in logistic format. The software also creates the Receiver Operating Characteristic Curve (ROC) to verify the model, and the area under ROC (AUC) can be used to evaluate the accuracy of the Maxent model. The AUC usually ranges from 0.5 to 1.0, with larger values indicating that the model is more accurate [57], while an AUC of 0.5 is judged to be of no reference value, and an AUC of 1.0 is judged to indicate a perfect prediction. Accordingly, the evaluation criteria can be set as follows: 0.5–0.6 for fail, 0.6–0.7 for poor, 0.7–0.8 for fair, 0.8–0.9 for good, and 0.9–1.0 for excellent. Only when the AUC values of both training data and test data reach good or excellent levels can the results of the Maxent model be incorporated into the CA model for further analysis.

2.3.3. Simulation of Urban Expansion Based on the CA Model

Permanent Resident Population Forecast

In the study of urban expansion simulation, the maximum expansion area is generally regarded as the termination threshold of circular iterations for the CA model and is one of the basic building blocks of the model. However, the latest plan for Wuxi is the Master Plan for Land Use of Wuxi (2006–2020), developed in 2011, which is not applicable to the simulation of Wuxi’s urban expansion in 2035. Moreover, the new round of territorial spatial planning has not yet been completed, so relevant data on the future urban land planning of Wuxi are not available. Therefore, from the perspective of population and land use, this study used the permanent resident population from 2006 to 2020 as the database to forecast the population size and urban land demand of Wuxi in 2035.
According to the existing urban land area and permanent resident population data, the per capita urban land area of Wuxi in 2015 was 118.55 m2. Referring to the Standard of Climatic Regionalization for Architecture (GB50178-93) and the Code for Classification of Urban Land Use and Planning Standards of Development Land (GB50137-2011), Wuxi is located in climate zone III, so its ideal urban land area per capita was set as 110 m2 based on the principle of regulated development and reasonable expansion. By using this value, the maximum urban land area required to accommodate the above population in 2035 can be obtained. Combined with the setting of 900 m2 per cell, the urban land area can be converted into a corresponding number of cells, which is used as the termination threshold of circular iterations.
In order to improve the accuracy of population forecast, three population forecasting methods, namely the comprehensive growth rate method, the unitary linear regression, and the GM (1, 1) grey prediction model, were used in this study to calculate the average result. In particular, the comprehensive growth rate method follows the expected indicators of the Wuxi 14th Five-year Plan for Population Development to set the annual growth rate of the permanent resident population as 1.4%. The grey prediction passed the quasi-smooth test ρ t and the quasi-exponential law test σ t , indicating that the population forecast is feasible. Detailed results are shown in Table 2. The average urban land raster obtained was 1,177,403, and this was used in the subsequent CA model.

CA Model

The CA model, with distinct spatio-temporal coupling characteristics, is especially suitable for the dynamic simulation of complex geographical systems, such as urban expansion. There are three important components to be considered in the construction of the CA model: the neighborhood effect, spatial constraint, and random interference. The neighborhood effect can control the self-organizing behavior of the model and thus influence the formation of a complex system [58,59]. The spatial constraint describes the areas where urban land development is prohibited, while random interference is used to simulate random and unpredictable factors. The state transition rule serves as the core that determines the relationship between the transition probability of urban land and the neighborhood effect, spatial constraint, and random interference [33].
The simulation process used in this study is based on the assumption that the CA model contains only a binary transformation; that is, it only contains two cellular attributes (urban land and non-urban land). Only the possibility of the conversion from non-urban land to urban land was considered. After the study area was rasterized with a resolution of 30 m, urban land was defined as 1 and non-urban land as 0, giving each cell its state S x t at a specific time step t . Its state S x t + 1 at the next time step can be affected by multiple factors, which can be expressed as follows.
S x t + 1 = f S x t ,   M x t ,   N x t ,   L x t ,   R x t ,
where M x t represents the transition potential of the cell driven by environmental variables, namely the urban construction suitability obtained by the Maxent model; N x t represents the neighborhood coefficient, namely the urban development intensity of the cellular neighborhood; L x t represents the spatial constraint of the restricted area on the cell; and R x t represents random interference.
For the neighborhood structure applied by the CA model, the Moore type is one of the most frequently used neighborhood structures. It consists of 8 surrounding cells, which, except for the central cell, are arranged in a 3 × 3 grid structure. The neighborhood coefficient N x t can be expressed as:
N x t = c o n S x t = 1 8 ,
Combined with the urban construction suitability and neighborhood effect, the global transition probability of each cell can be calculated.
p x t = M x t + w N N x t ,
where w N represents the weight coefficient of neighborhood effect.
Thus, the final probability of cellular transition P x t is shown as:
P x t = e x p α p x t p x . m a x t 1 R x t L x t ,
where α represents the diffusion coefficient, and p x . m a x t represents the maximum global transition probability in each iteration loop, which is continuously updated in each subloop.
The spatial constraint of the CA model mainly includes four restricted areas: the water area, forest land, ecological protection red line, and permanent basic farmland in Wuxi. The water area and forest land data were extracted from the Wuxi land use data. The four data layers were merged and set as the forbidden construction area of urban land in the CA model, as shown by the following formula:
L x t = 1 ,   not   subject   to   restricted   zones 0 ,   subject   to   restricted   zones ,
According to P x t in each iteration loop, the model updated the cellular state as well as counting the number s t e p N U M of cells whose state changed to urban land and the total number s u m N U M of cells whose state was urban land in this iteration. When the s u m N U M was equal to or greater than the iteration threshold, that is, when the total size of urban land exceeded the maximum urban land area calculated by permanent resident population forecast, the model stopped. The distribution result of urban land at this time was taken as the prediction result for urban expansion in Wuxi in 2035.

2.3.4. Identification of Three Urban Expansion Types

In order to quantitatively describe the characteristics and distribution pattern of urban expansion in Wuxi in 2035, this study converted the new urban land data into vector format and extracted the common edge according to the spatial attributes. The ratio between the length of the common edge and the total perimeter of each patch was used as the basis for classification [52,60]. The new urban land patches were divided into three urban expansion types, infilling, edge expansion, and outlying, for visual display [61].
P = l c l ,
where l c represents the length of the common edge between a new urban patch and an old urban patch, and l represents the total perimeter of the new urban patch. The value of P ranges from 0 to 1. The infilling expansion has 0.5 < P 1 , which is characterized as old urban land surrounding more than half of the new urban area. The edge expansion has 0 < P 0.5 , which refers to the development of new urban land spreading outward from the edge of an old urban area and surrounded by less than 50% old urban area. The outlying expansion has P = 0 , which indicates that there is no common edge between the new urban land and old urban land.
The above methods can be integrated into the following research process (Figure 3).

3. Results

3.1. Analysis of the Maxent Model

3.1.1. Influence of the Sample Size on the Model Accuracy

Under the same environmental constraints, different numbers of sample points were randomly selected to predict the maximum entropy distribution. The AUC values of the training data set and test data set are shown in Table 3. The results show that when the sample size exceeded 1000, the model accuracy decreased as the sample size increased, and when the sample size was between 800 and 1000, the AUC values of both were greater than 0.8, indicating that the Maxent model can predict urban expansion well. The follow-up results on the mechanism of environmental variables show that the best result was obtained with a sample size of 800. This value was applied to subsequent simulations to ensure the high accuracy of the Maxent-CA model.

3.1.2. The Contribution and Importance of Variables

Eight hundred random sample points were used for the trial operation of the Maxent model. The contributions and importance of 11 environmental variables to the maximum entropy distribution of urban land in Wuxi are shown in Table 4. The contribution covers the variable interaction, while the importance reflects the dependent influence of a single environmental variable on the maximum entropy distribution.
The results show that the contributions and importance of the aspect, infra, and GDP are relatively low. It is speculated this is due to the predominantly plain topography, the generally developed economy, and the perfect infrastructure of Wuxi. The homogenization of regional nature leads to homogeneity in the construction conditions. After excluding these three low-value variables, the new index system was applied to the model. The AUC values of the training data set and test data set were 0.817 and 0.806, respectively, indicating that the model is credible. The final results are shown in Table 5.
Among the eight environmental variables, the more prominent contributors are road (30.23%), NDVI (22.69%) and indus (16.07%), with a total contribution of 68.99%, indicating that these three variables play key roles in the systematic model of variable interaction. While the single variables with large influences are the NDVI (22.91%), gov (17.52%), and road (17.24%), the importance of traffic, indus, and pop fluctuate around 10–15%, which means that their influences on the model are basically equal. Overall, the contributions and importance of socioeconomic conditions are higher than those of natural conditions, indicating that the urban expansion process in Wuxi is more likely to be affected by human socioeconomic activities. In addition, the contributions of gov and pop are obviously lower than their importance levels, while the contributions of road and indus are significantly higher than their importance levels, indicating that a system with variable interaction will reduce the roles of gov and pop in the model, and conversely, the roles of road and indus will increase.
The model also outputs the jackknife test of the variables’ regularized training gain (Figure 4), where the dark bar represents the extent to which the environment variable affects the model when used in isolation. The traffic is found to have the highest gain and therefore appears to represent the most useful information. The light bar represents the extent to which the environment variable decreases the gain the most when it is omitted, where the shorter the bar is, the more indispensable the presence of the environment variable is to the gain in results. The NDVI has the shortest bar, indicating that the NDVI contains the most information that is not present in the other variables and has the greatest impact on the prediction results when it changes alone.

3.1.3. The Relationship between the Probability Distribution of Urban Land and Variables

The response curve can reflect the relationship between the existence probability and environmental variables, that is, the quantitative change between the probability distribution of urban land and environmental variables (Figure 5). The interpretation of the response curve provides a deeper understanding of the relationship between the dependent and independent variables. A probability of greater than 0.5 indicates that the numerical interval of the environmental variable contributes to urban expansion and vice versa.
The response curve of soil shows that the patches of anthropogenic soil (code 5) are most likely to be covered by urban expansion, while the possibility of the patches of lake–reservoir (code 6) being converted into urban land decreases sharply, which is consistent with general cognition. The NDVI response curve shows an increasing and then decreasing trend, with the optimum range for facilitating urban expansion being between 0.2 and 0.6. When the vegetation coverage is higher than a certain level, the patches are not suitable for urban development. The relief curve shows that areas with a topographic difference of more than 10 degrees are also not suitable. In terms of socioeconomic conditions, 1500 to 9000 people per square kilometer is a relatively optimal range for urban expansion, and 2000 people per square kilometer is the optimal density. The response curves of gov, indus, road, and traffic all show non-linear negative correlations with the probability of urban expansion decreasing with an increasing distance.

3.1.4. Maximum Entropy Distribution of Urban Land Based on the Maxent Model

High values of information entropy in the maximum entropy distribution of urban land in Wuxi (Figure 6) are mainly found in the central urban areas of Wuxi city, Jiangyin city, and Yixing city, indicating that future urban expansion will tend to occur in these places which already have the foundation for urban development and where the dense populations and convenient transportation are conducive to urban expansion. The low information entropy in the scattered areas of Jiangyin city and most areas of Yixing city, on the other hand, make it difficult for these areas to be reached by urban economic radiation due to their distance from the existing urban land, resulting in a lower rate of urban development. Additionally, the information entropy of water areas, such as Taihu Lake, is extremely low, which reflects the constraint of the water environment on urban expansion in Wuxi.

3.2. Analysis of the Maxent-CA Model

3.2.1. The Simulation of Urban Expansion Based on the Maxent-CA Model

A total of 252 iterations of the Maxent-CA model were used to generate 1,184,450 urban cells to simulate the urban expansion of Wuxi in 2035. From the perspective of the overall distribution pattern of urban land, the total scale of urban land is expected to reach 1066.01 km2 in 2035. The new urban land is predicted to be mainly distributed in the center of Wuxi city and the subcentral urban area of Jiangyin city and scattered in the subcenter of Yixing city. The northern parts surrounded by the Yangtze River and the eastern parts adjacent to the Taihu Lake Basin are expected to be the key areas for future development (Figure 7) while adhering to the bottom line of ecological protection. From the perspective of urban expansion types (Table 6), 57.37% of the new urban land is predicted to be infilling expansion concentrated in the Jiangyin High-tech Industrial Development Zone (mainly in Chengdong Street) and the Wuxi Airport Development Zone (mainly in Shuofang Street and Hongshan Street), mostly clustered in blocks; 37.83% of the new urban land is predicted to be edge-expansion, mainly distributed in the Huishan Economic Development Zone, Anzhen Street (the eastern part of the Xishan Economic Development Zone) and the corridor of the Huyi Highway between Qianqiao Street and Hu Dai Town; and only 4.80% is predicted to be outlying expansion, with the typical areas such as Zhouzhuang Town and Changjing Town developing as node towns on the main traffic routes.
In general, urban expansion basically spreads outwards around the distribution range of the old urban land and focuses on the transition of non-urban land surrounded by original urban land to meet the requirements of vacant landfilling, patch integration, and concentrated contiguity. New municipalities such as the Kegang–Huangtu Cluster, Donggang–Xibei Cluster and Zhangzhu–Hufu Cluster are predicted to develop in integrated clusters, integrating into the overall urban development and connecting three main urban areas through the Taihu–Yangtze River Development Axis and the Taihu Bay Science and Technology Innovation Belt. This simulation may allow the co-construction and sharing of public services and infrastructure, promoting more compact and intensive urban spatial development.

3.2.2. The Demarcation of the UDB

Based on the urban expansion simulation results obtained by the Maxent-CA model, this study used the Guidelines for the Compilation of the Municipal Territorial Spatial Master Plan (Trial Implementation) issued by the Ministry of Natural Resources as a reference to conduct vacant landfilling, patch deletion, and plot merging on the simulation results. Originally, the enclosed area within a single closure line was more than 30 hectares in size. Restricted areas, such as water areas, forest land, ecological protection red line areas, and permanent basic farmland were not filled in as necessary avoidance areas. Finally, the patches of node towns or special towns with small areas were kept. Based on the above three points, the demarcation results of the UDB are shown in Figure 8. The area of urban land is 1066.01 km2, slightly exceeding the threshold of 1059.66 km2 obtained from the forecast of the permanent resident population by 0.60%. After strict correction, the area of urban land within the rigid UDB is 1053.46 km2, which is approximately 99.41% of the maximum expansion area. On the whole, the results are in line with the study’s expectations.
In general, the predicted UDB of Wuxi in 2035 meets the requirements of the urban spatial development strategy proposed in the Wuxi Territorial Spatial Master Plan (2021–2035). This UDB mainly covers the Xinwu New District in the southeast of the main urban area of Wuxi, the Xidong urban area in the east, the Huishan urban area in the north, and the new urban area on the periphery of the Jiangyin sub-city, along with the old urban area of Yixing city. It focuses on guiding the construction activities of industrial development zones, such as high-tech industries, advanced manufacturing, and productive services. Meanwhile, the demarcation of the UDB also strictly abides by the constraint rules of restricted areas, preserving a certain amount of ecological space and agricultural space for urban development. By reserving ecological nodes such as Taihu Lake Basin, Yangtze River Wetland in Jiangyin, and National Forest Park in Yixing and avoiding permanent basic farmland, the pressure on Wuxi’s ecosystem is alleviated, providing a guarantee for the city’s sustainable development.

4. Discussion and Conclusions

4.1. Discussion

4.1.1. The Accuracy of the Maxent-CA Model

In order to ensure the feasibility of the Maxent-CA model, we used land use data from 2005 and 2015 to simulate the urban expansion of Wuxi in 2020, using the Maxent-CA model and logistic-CA model, respectively. In particular, the land use data from Wuxi in 2020 were sourced from the RESDC, from which the actual area of the extracted urban land was determined to be 874.69 km2. According to the setting of 900 m2 per cell, the actual number of cells corresponding to the urban land of Wuxi in 2020 was 971,873, and this was used as the termination threshold of the circular iteration. The results showed the Kappa index of the Maxent-CA model was 93.87% and that of the logistic-CA model was 90.56%. In terms of spatial distribution (Figure 9), the logistic-CA model showed a more fragmented distribution of urban land, with more small urban patches in remote areas. The Maxent-CA model showed a more concentrated and contiguous distribution of urban land, effectively controlling the disorderly spread of urban land. This could provide technical support for the demarcation of a complete, scientifical, concise, and easily implemented UDB. The comparison proved that the Maxent-CA model was applicable to the prediction of the land use pattern. The results of the model were credible for the urban expansion simulation and UDB demarcation of Wuxi in 2035.

4.1.2. Further Analysis

Through an in-depth analysis of the Maxnet-CA model, the original urban land in 2015 was found to occupy 928,127 cells whose states remained as urban land in 2020, accounting for an absolute proportion of the simulation results and therefore bringing a high Kappa index. In contrast, 48,634 cells were shown to be incorrectly simulated and 49,993 cells exceeded the actual scope of urban land. Combined with Figure 10, further spatial analysis of the simulation results revealed that the “actual, not simulated” area (Figure 10a) is mainly distributed in the peripheral area separated from the original urban land, and the “simulated, not actual” area (Figure 10b) is mainly concentrated in the central area surrounded by the original urban land. This indicates that the model has a degree of deviation when simulating the outlying expansion of urban land. It is actually more applicable to the small study area where urban expansion is mainly characterized by a cluster-like distribution and central infilling. The reason is that urban development is a dynamic and open evolution process and outlying expansion tends to follow policy-oriented guidance [52], which is more uncertain and complex. The Maxent-CA model summarizes the land transition rules based on the distribution characteristics of existing urban land, so it cannot fully predict future urban development.
On the other hand, the inadequate protection of key ecological restricted areas is also a major reason affecting the accuracy. The model result shows that 25,550 cells in the restricted area (52.53% of “actual, not simulated” area) were occupied by urban land in 2015, including about 22.99 km2 of water area, forest land, ecological protection red line area, and permanent basic farmland area that is not actually protected (such as Huishan National Forest Park shown in Figure 10c). As a result, the simulation process has an error of 25,550 cells in advance, which impacts the accuracy of the results. The Maxent-CA model in this paper strictly prohibits urban expansion within the restricted areas to protect the key ecological reserves, which ensures the scientifical and reasonable demarcation of UDB.
Additionally, while the Maxent-CA model took full advantage of the Maxent model’s advantage of automatically obtaining weight coefficients from images of a single period and a single class, it also exposed the model’s shortcomings. Actually, the land use planning and management tend to involve the interactive transformation of multiple land use types [45,62]. Thus, the Maxent-CA model would have limitations in practical application. Meanwhile, it was difficult to obtain all the data at 30 m resolution directly. The non-spatial data used in this paper, such as GDP and population distribution, were generally only counted to county-level administrative regions, which had a large resolution. There was a possibility of errors caused by loss of details in the process of resampling to 30 m.
In general, the future study needs to consider the influence of political decisions and human activities on urban expansion more comprehensively. Using an improved model and data with higher resolution, more accurate simulation results are expected to be obtained to improve the practical value of the Maxent-CA model in UDB demarcation. Through the demarcation of UDB, we hope to effectively protect the ecological land, rationalize the layout of functional urban land, and control urban expansion.

4.2. Conclusions

In this study, a Maxent-CA model was constructed to demarcate the UDB. The model not only integrates the advantage of the Maxent model by automatically obtaining the weight coefficients of independent variables from images of a single period and a single class but also the advantage of the CA model by realizing spatio-temporal dynamic simulation in neighborhood structure. Based on the Maxent-CA model, this study took a typical rapidly urbanizing city (Wuxi) as the study area, land use data from 2005 and 2015 as the dependent variable, and selected environmental variables as the index system. The model was used to simulate the urban expansion of Wuxi in 2035 to demarcate the UDB, demonstrating the value of the Maxent-CA model in urban master planning. The results show that:
(1)
The Maxent-CA model proposed in this paper can be applied to demarcate the UDB of specific cities. This study comprehensively considered the relationship between the urban construction suitability, neighborhood effect, spatial constraint, and random interference of urban expansion and fully embodied the principle of combining top-down and bottom-up UDB demarcation approaches in the land use master plan.
(2)
The Maxent-CA model can intuitively reflect the driving mechanism of urban expansion. From the perspective of the importance of a single variable, NDVI and gov were found to have higher influences. However, road and indus were shown to contribute more to the model influenced by the interactions between variables. Moreover, the response relationship between urban expansion and environmental variables is complex and non-linear.
(3)
The expansion of new urban land is dominated by the infilling expansion type, followed by edge expansion, with outlying expansion being the least common. The accuracy of the model was shown to decrease when simulating outlying expansion and in situations where the key ecological restricted areas were not adequately protected. Future study needs to consider the influence of political decisions and human activities on urban expansion more comprehensively. The relevant departments should pay attention to the phenomenon of urban land occupation and coordinate the scale of urban and rural land.
Under the current urban spatial planning system in China, the policy criteria, methods, and management systems of boundary demarcation for urban and rural planning, land use planning, development priority zones, and so on vary among different administrative departments due to their different requirements. The method of demarcating the UDB based on the Maxent-CA model proposed in this paper has great potential. More in-depth research should be conducted in the future to explore the reference value of the Maxent-CA model in territorial spatial planning.

Author Contributions

Formal analysis and investigation, J.Z.; writing—original draft preparation, J.Z.; writing—review and editing, J.Z. and Y.C.; methodology, X.Y.; resources, J.Z., W.Q. and D.W.; data curation, J.Z.; supervision, Y.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the project “Assessment of Resource and Environmental Carrying Capacity and Territory Spatial Development Suitability in Wuxi” (0209151950).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bai, X.; Shi, P.; Liu, Y. Society: Realizing China’s urban dream. Nature 2014, 509, 158–160. [Google Scholar] [CrossRef] [PubMed]
  2. Kuang, W. 70 years of urban expansion across China: Trajectory, pattern, and national policies. Sci. Bull. 2020, 65, 1970–1974. [Google Scholar] [CrossRef]
  3. Guo, B.; Bian, Y.; Pei, L.; Zhu, X.; Zhang, D.; Zhang, W.; Guo, X.; Chen, Q. Identifying Population Hollowing Out Regions and Their Dynamic Characteristics across Central China. Sustainability 2022, 14, 9815. [Google Scholar] [CrossRef]
  4. Lyu, X.; Peng, W.; Niu, S.; Qu, Y.; Xin, Z. Evaluation of sustainable intensification of cultivated land use according to farming households’ livelihood types. Ecol. Indic. 2022, 138, 108848. [Google Scholar] [CrossRef]
  5. Liu, Y.; Li, J.; Yang, Y. Strategic adjustment of land use policy under the economic transformation. Land Use Policy 2018, 74, 5–14. [Google Scholar] [CrossRef]
  6. Long, H.; Liu, Y.; Hou, X.; Li, T.; Li, Y. Effects of land use transitions due to rapid urbanization on ecosystem services: Implications for urban planning in the new developing area of China. Habitat Int. 2014, 44, 536–544. [Google Scholar] [CrossRef]
  7. Lyu, X.; Peng, W.; Yu, W.; Xin, Z.; Niu, S.; Qu, Y. Sustainable intensification to coordinate agricultural efficiency and environmental protection: A systematic review based on metrological visualization. J. Land Use Sci. 2021, 16, 313–338. [Google Scholar] [CrossRef]
  8. Niu, S.; Lyu, X.; Gu, G. A New Framework of Green Transition of Cultivated Land-Use for the Coordination among the Water-Land-Food-Carbon Nexus in China. Land 2022, 11, 933. [Google Scholar] [CrossRef]
  9. Liu, Y.; Zhou, Y. Territory spatial planning and national governance system in China. Land Use Policy 2021, 102, 105288. [Google Scholar] [CrossRef]
  10. Lu, X.; Qu, Y.; Sun, P.; Yu, W.; Peng, W. Green Transition of Cultivated Land Use in the Yellow River Basin: A Perspective of Green Utilization Efficiency Evaluation. Land 2020, 9, 475. [Google Scholar] [CrossRef]
  11. Xie, X.; Kang, H.; Behnisch, M.; Baildon, M.; Krüger, T. To What Extent Can the Green Belts Prevent Urban Sprawl?—A Comparative Study of Frankfurt am Main, London and Seoul. Sustainability 2020, 12, 679. [Google Scholar] [CrossRef]
  12. Ball, M.; Cigdem, M.; Taylor, E.; Wood, G. Urban growth boundaries and their impact on land prices. Environ. Plan. A 2014, 46, 3010–3026. [Google Scholar] [CrossRef]
  13. Wang, Y.; Gu, C.; Li, X. Research progress of urban growth boundary at home and abroad. Urban Plan. Int. 2014, 29, 1–11. Available online: http://lib.cqvip.com/Qikan/Article/Detail?id=661962783 (accessed on 8 December 2020). (In Chinese).
  14. Ding, C.; Knaap, G.J.; Hopkins, L.D. Managing urban growth with urban growth boundaries: A theoretical analysis. J. Urban Econ. 1999, 46, 53–68. [Google Scholar] [CrossRef]
  15. Hepinstall-Cymerman, J.; Coe, S.; Hutyra, L.R. Urban growth patterns and growth management boundaries in the Central Puget Sound, Washington, 1986–2007. Urban Ecosyst. 2013, 16, 109–129. [Google Scholar] [CrossRef]
  16. Nelson, A.C.; Moore, T. Assessing urban growth management: The case of Portland, Oregon, the USA’s largest urban growth boundary. Land Use Policy 1993, 10, 293–302. [Google Scholar] [CrossRef]
  17. Zhang, J. Urban growth management in the United States. Urban Plan. Overseas 2002, 2, 37–40. (In Chinese) [Google Scholar] [CrossRef]
  18. Lin, J.; Qiao, Z.; Ye, Z. “Delimitation” and “implementation” of urban growth boundary: Analysis and thoughts on the practice in 14 pilot cities in China. Urban Plan. Forum 2017, 4, 37–43. (In Chinese) [Google Scholar] [CrossRef]
  19. Wang, W.; Jiao, L.; Zhang, W.; Jia, Q.; Su, F.; Xu, G.; Ma, S. Delineating urban growth boundaries under multi-objective and constraints. Sustain. Cities Soc. 2020, 61, 102279. [Google Scholar] [CrossRef]
  20. Wang, Z.; Zhang, Q.; Zhang, X.; Guan, X. Urban growth boundary delimitation of Hefei City based on the resources and environment carrying capability. Geogr. Res. 2013, 32, 2302–2311. (In Chinese) [Google Scholar] [CrossRef]
  21. Li, X.; Yeh, A.G.O. Cellular automata for simulating complex land use systems using neural networks. Geogr. Res. 2005, 24, 19–27. (In Chinese) [Google Scholar] [CrossRef]
  22. He, Q.; Tan, R.; Gao, Y.; Zhang, M.; Xie, P.; Liu, Y. Modeling urban growth boundary based on the evaluation of the extension potential: A case study of Wuhan city in China. Habitat Int. 2018, 72, 57–65. [Google Scholar] [CrossRef]
  23. Liu, X.; Wei, M.; Li, Z.; Zeng, J. Multi-scenario simulation of urban growth boundaries with an ESP-FLUS model: A case study of the Min Delta region, China. Ecol. Indic. 2022, 135, 108538. [Google Scholar] [CrossRef]
  24. Salem, M.; Bose, A.; Bashir, B.; Basak, D.; Roy, S.; Chowdhury, I.R.; Alsalman, A.; Tsurusaki, N. Urban Expansion Simulation Based on Various Driving Factors Using a Logistic Regression Model: Delhi as a Case Study. Sustainability 2021, 13, 10805. [Google Scholar] [CrossRef]
  25. Tayyebi, A.; Pijanowski, B.C.; Tayyebi, A.H. An urban growth boundary model using neural networks, GIS and radial parameterization: An application to Tehran, Iran. Landsc. Urban Plan. 2011, 100, 35–44. [Google Scholar] [CrossRef]
  26. Silva, E.A.; Clarke, K.C. Calibration of the SLEUTH urban growth model for Lisbon and Porto, Portugal. Comput. Environ. Urban Syst. 2002, 26, 525–552. [Google Scholar] [CrossRef]
  27. Yin, H.; Kong, F.; Yang, X.; James, P.; Dronova, I. Exploring zoning scenario impacts upon urban growth simulations using a dynamic spatial model. Cities 2018, 81, 214–229. [Google Scholar] [CrossRef]
  28. Liu, X.; Liang, X.; Li, X.; Xu, X.; Ou, J.; Chen, Y.; Li, S.; Wang, S.; Pei, F. A future land use simulation model (FLUS) for simulating multiple land use scenarios by coupling human and natural effects. Landsc. Urban Plan. 2017, 168, 94–116. [Google Scholar] [CrossRef]
  29. Ma, S.; Li, X.; Cai, Y. Delimiting the urban growth boundaries with a modified ant colony optimization model. Comput. Environ. Urban Syst. 2017, 62, 146–155. [Google Scholar] [CrossRef]
  30. Liu, J.; Chen, Y. GIS-based cellular automata models and researches on spatial complexity of man-land relationship. Geogr. Res. 2002, 21, 155–162. (In Chinese) [Google Scholar] [CrossRef]
  31. White, R.; Engelen, G.; Uljee, I. The use of constrained cellular automata for high-resolution modelling of urban land-use dynamics. Environ. Plan. B 1997, 24, 323–343. [Google Scholar] [CrossRef]
  32. Li, X.; Chen, G.; Liu, X.; Liang, X.; Wang, S.; Chen, Y.; Pei, F.; Xu, X. A new global land-use and land-cover change product at a 1-km resolution for 2010 to 2100 based on human-environment interactions. Ann. Am. Assoc. Geogr. 2017, 107, 1040–1059. [Google Scholar] [CrossRef]
  33. Clarke, K.C.; Gaydos, L.J. Loose-coupling a cellular automaton model and GIS: Long-term urban growth prediction for San Francisco and Washington/Baltimore. Int. J. Geogr. Inf. Sci. 1998, 12, 699–714. [Google Scholar] [CrossRef] [PubMed]
  34. Liang, X.; Liu, X.; Li, X.; Chen, Y.; Tian, H.; Yao, Y. Delineating multi-scenario urban growth boundaries with a CA-based FLUS model and morphological method. Landsc. Urban Plan. 2018, 177, 47–63. [Google Scholar] [CrossRef]
  35. Lv, J.; Wang, Y.; Liang, X.; Yao, Y.; Ma, T.; Guan, Q. Simulating urban expansion by incorporating an integrated gravitational field model into a demand-driven random forest-cellular automata model. Cities 2021, 109, 103044. [Google Scholar] [CrossRef]
  36. Long, Y.; Han, H.; Lai, S.; Mao, Q. Urban growth boundaries of the Beijing Metropolitan Area: Comparison of simulation and artwork. Cities 2013, 31, 337–348. [Google Scholar] [CrossRef]
  37. Sun, W.; Liu, C.; Wang, S. Simulation research of urban development boundary based on ecological constraints: A case study of Nanjing. J. Nat. Resour. 2021, 36, 2913–2925. [Google Scholar] [CrossRef]
  38. Barredo, J.I.; Demicheli, L. Urban sustainability in developing countries’ megacities: Modelling and predicting future urban growth in Lagos. Cities 2003, 20, 297–310. [Google Scholar] [CrossRef]
  39. Yi, D.; Zhao, X.; Guo, X.; Zhao, L.; Zhang, H.; Han, Y.; Roshan, S.; Luo, Z. Delimitation of urban development boundary based on ecological sensitivity evaluation and CA-Markov simulation in plain city: A case of Nanchang, Jiangxi, China. Chin. J. Appl. Ecol. 2020, 31, 208–218. (In Chinese) [Google Scholar] [CrossRef]
  40. Phillips, S.J.; Anderson, R.P.; Dudík, M.; Schapire, R.E.; Blair, M.E. Opening the black box: An open-source release of Maxent. Ecography 2017, 40, 887–893. [Google Scholar] [CrossRef]
  41. Kaky, E.; Nolan, V.; Alatawi, A.; Gilbert, F. A comparison between Ensemble and MaxEnt species distribution modelling approaches for conservation: A case study with Egyptian medicinal plants. Ecol. Inform. 2020, 60, 101150. [Google Scholar] [CrossRef]
  42. Wang, Y.; Chao, B.; Dong, P.; Zhang, D.; Yu, W.; Hu, W.; Ma, Z.; Chen, G.; Liu, Z.; Chen, B. Simulating spatial change of mangrove habitat under the impact of coastal land use: Coupling MaxEnt and Dyna-CLUE models. Sci. Total Environ. 2021, 788, 147914. [Google Scholar] [CrossRef] [PubMed]
  43. Lin, C.; Hong, Y.; Tu, W.; Liu, J. Multiperiod Dynamic Programming Algorithm for Optimizing a Nature Reserve. Sustainability 2022, 14, 3266. [Google Scholar] [CrossRef]
  44. Tan, J.; Li, A.; Lei, G.; Xie, X. A SD-MaxEnt-CA model for simulating the landscape dynamic of natural ecosystem by considering socio-economic and natural impacts. Ecol. Model. 2019, 410, 108783. [Google Scholar] [CrossRef]
  45. Zhang, Y.; Liu, X.; Chen, G.; Hu, G. Simulation of urban expansion based on cellular automata and maximum entropy model. Sci. Sin. 2020, 50, 339–352. [Google Scholar] [CrossRef]
  46. Chen, Z. Simulation and Prediction Analysis of Land Use Change in Changzhou Based on Maxent-CA. Geomat. World 2020, 27, 29–35,41. (In Chinese) [Google Scholar] [CrossRef]
  47. Tian, G.; Jiang, J.; Yang, Z.; Zhang, Y. The urban growth, size distribution and spatio-temporal dynamic pattern of the Yangtze River Delta megalopolitan region, China. Ecol. Model. 2011, 222, 865–878. [Google Scholar] [CrossRef]
  48. Liu, J.; Zhang, Z.; Zhuang, D.; Wang, Y.; Zhou, W.; Zhang, S.; Li, R.; Jiang, N.; Wu, S. A study on the spatial-temporal dynamic changes of land-use and driving forces analyses of China in the 1990s. Geogr. Res. 2003, 22, 1–12. (In Chinese) [Google Scholar] [CrossRef]
  49. Tripathy, P.; Kumar, A. Monitoring and modelling spatio-temporal urban growth of Delhi using Cellular Automata and geoinformatics. Cities 2019, 90, 52–63. [Google Scholar] [CrossRef]
  50. Maimaitijiang, M.; Ghulam, A.; Sandoval, J.S.O.; Maimaitiyiming, M. Drivers of land cover and land use changes in St. Louis metropolitan area over the past 40 years characterized by remote sensing and census population data. Int. J. Appl. Earth Obs. 2015, 35, 161–174. [Google Scholar] [CrossRef]
  51. Zhang, R.; Pu, L.; Zhu, M. Impacts of transportation arteries on land use patterns in urban-rural fringe: A comparative gradient analysis of Qixia District, Nanjing City, China. Chin. Geogr. Sci. 2012, 23, 378–388. [Google Scholar] [CrossRef]
  52. Sun, C.; Wu, Z.; Lv, Z.; Yao, N.; Wei, J. Quantifying different types of urban growth and the change dynamic in Guangzhou using multi-temporal remote sensing data. Int. J. Appl. Earth Obs. 2013, 21, 409–417. [Google Scholar] [CrossRef]
  53. Yang, J.; Xi, J.; Xie, P.; Ge, Q.; Li, X.; Gao, H. Spatio-temporal simulation and application of the tourist town growth based on cellular automaton model: A case study of Sanpo Town in Hebei Province. Geogr. Res. 2013, 32, 915–923. (In Chinese) [Google Scholar] [CrossRef]
  54. Zhang, Q.; Su, S. Determinants of urban expansion and their relative importance: A comparative analysis of 30 major metropolitans in China. Habitat Int. 2016, 58, 89–107. [Google Scholar] [CrossRef]
  55. Phillips, S.J.; Dudík, M.; Schapire, R.E. Maxent Software for Modeling Species Niches and Distributions (version 3.4.1). Available online: http://biodiversityinformatics.amnh.org/open_source/maxent/ (accessed on 25 November 2020).
  56. Phillips, S.J.; Anderson, R.P.; Schapire, R.E. Maximum entropy modeling of species geographic distributions. Ecol. Model. 2006, 190, 231–259. [Google Scholar] [CrossRef]
  57. Swets, J.A. Measuring the accuracy of diagnostic systems. Science 1988, 240, 1285–1293. [Google Scholar] [CrossRef]
  58. Feng, Y.; Tong, X. Dynamic land use change simulation using cellular automata with spatially nonstationary transition rules. GISci. Remote Sens. 2018, 55, 678–698. [Google Scholar] [CrossRef]
  59. Feng, Y.; Han, Z. Impact of neighbor configurations on spatially-explicit modeling results. Geogr. Res. 2011, 30, 1055–1065. (In Chinese) [Google Scholar] [CrossRef]
  60. Xu, C.; Liu, M.; Zhang, C.; An, S.; Yu, W.; Chen, J. The spatiotemporal dynamics of rapid urban growth in the Nanjing metropolitan region of China. Landsc. Ecol. 2007, 22, 925–937. [Google Scholar] [CrossRef]
  61. Wilson, E.H.; Hurd, J.D.; Civco, D.L.; Prisloe, M.P.; Arnold, C. Development of a geospatial model to quantify, describe and map urban growth. Remote Sens. Environ. 2003, 86, 275–285. [Google Scholar] [CrossRef]
  62. Dadashpoor, H.; Azizi, P.; Moghadasi, M. Land use change, urbanization, and change in landscape pattern in a metropolitan area. Sci. Total Environ. 2019, 655, 707–719. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Study area.
Figure 1. Study area.
Sustainability 14 11426 g001
Figure 2. Dependent variable and environmental variables (ak) in the urban expansion model. (a) Surface relief, (b) Aspect, (c) NDVI, (d) Soil type, (e) Distance to road, (f) Distance to transport station, (g) Distance to infrastructure, (h) Distance to industrial enterprise, (i) Distance to government department, (j) GDP density and (k) Population density.
Figure 2. Dependent variable and environmental variables (ak) in the urban expansion model. (a) Surface relief, (b) Aspect, (c) NDVI, (d) Soil type, (e) Distance to road, (f) Distance to transport station, (g) Distance to infrastructure, (h) Distance to industrial enterprise, (i) Distance to government department, (j) GDP density and (k) Population density.
Sustainability 14 11426 g002
Figure 3. Framework of the Maxent-CA model for the demarcation of the UDB.
Figure 3. Framework of the Maxent-CA model for the demarcation of the UDB.
Sustainability 14 11426 g003
Figure 4. Jackknife test of the variables’ regularized training gain.
Figure 4. Jackknife test of the variables’ regularized training gain.
Sustainability 14 11426 g004
Figure 5. Response curves illustrating the relationship between the probability distribution of urban land and variables.
Figure 5. Response curves illustrating the relationship between the probability distribution of urban land and variables.
Sustainability 14 11426 g005aSustainability 14 11426 g005b
Figure 6. Maximum entropy distribution of urban land in Wuxi.
Figure 6. Maximum entropy distribution of urban land in Wuxi.
Sustainability 14 11426 g006
Figure 7. Spatial distribution of urban expansion types.
Figure 7. Spatial distribution of urban expansion types.
Sustainability 14 11426 g007
Figure 8. The demarcation of the UDB in 2035.
Figure 8. The demarcation of the UDB in 2035.
Sustainability 14 11426 g008
Figure 9. Comparison between the results of the Maxent-CA model and logistic-CA model in 2020.
Figure 9. Comparison between the results of the Maxent-CA model and logistic-CA model in 2020.
Sustainability 14 11426 g009
Figure 10. Comparison of the simulated and actual scope of urban land of Wuxi in 2020. (a) The typical “actual, not simulated” area, (b) the typical “simulated, not actual” area, and (c) the typical unprotected area.
Figure 10. Comparison of the simulated and actual scope of urban land of Wuxi in 2020. (a) The typical “actual, not simulated” area, (b) the typical “simulated, not actual” area, and (c) the typical unprotected area.
Sustainability 14 11426 g010
Table 1. Specific details of environmental variables in the Maxent-CA model.
Table 1. Specific details of environmental variables in the Maxent-CA model.
TypeVariableAbbreviationMinMaxUnit
Natural
conditions
Surface reliefrelief0175m
Aspectaspect0360°
NDVINDVI0.0080.9
Soil typesoilIron bauxite: 1; alfisols: 2; primary soil: 3; semi-hydromorphic soil: 4; anthropogenic soil: 5; lake–reservoir: 6; river: 7
Socio-economic conditionsDistance to roadroad013,872.1m
Distance to
transport station
traffic013,754.7m
Distance to infrastructureinfra015,171.5m
Distance to
industrial enterprise
indus015,540.8m
Distance to
government department
gov013,738.7m
GDP densityGDP834248,673yuan/km2
Population densitypop501.310,949.5pop/km2
Table 2. Prediction of the urban expansion scale for Wuxi in 2035.
Table 2. Prediction of the urban expansion scale for Wuxi in 2035.
MethodPopulation/Thousand PeopleArea/m2Raster NumberAccuracy
Comprehensive growth rate9194.771011.421,123,805
Unitary linear regression9640.301060.431,178,259R2 = 0.9371
GM (1, 1) model10,064.831107.131,230,146δ = 2.14%
Average9633.301059.661,177,403
Table 3. Accuracy of the Maxent model with different sample sizes.
Table 3. Accuracy of the Maxent model with different sample sizes.
Sample Size60080010002000500010,000
11 environment variablesAUC of the training data set0.8250.8220.8160.7860.7390.685
AUC of the test data set0.8090.8110.8130.7790.7340.684
8 environment variablesAUC of the training data set0.8250.8170.8140.7850.7380.685
AUC of the test data set0.7840.8060.8050.7770.7340.692
Table 4. Preliminary results for the contribution and importance of each environmental variable.
Table 4. Preliminary results for the contribution and importance of each environmental variable.
VariableContribution (%)Importance (%)
road27.7420.56
NDVI23.5122.86
pop13.9911.37
indus10.809.99
traffic9.0411.06
relief8.752.77
gov3.6311.32
soil1.142.67
GDP0.583.34
infra0.553.03
aspect0.281.04
Table 5. Final results for the contribution and importance of environmental variables.
Table 5. Final results for the contribution and importance of environmental variables.
TypeVariableContribution (%)Importance (%)
Natural conditionsNDVI22.6922.91
relief6.243.11
soil0.553.64
Socio-economic conditionsroad30.2317.24
indus16.079.53
traffic11.3814.63
gov7.3817.52
pop5.4611.43
Table 6. Urban expansion types of Wuxi in 2035.
Table 6. Urban expansion types of Wuxi in 2035.
TypeOutlying ExpansionEdge ExpansionInfilling ExpansionTotal
Area (km2)11.0887.30132.40230.78
Percent (%)4.8037.8357.37100.00
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhang, J.; Chen, Y.; Yang, X.; Qiao, W.; Wang, D. The Demarcation of Urban Development Boundary Based on the Maxent-CA Model: A Case Study of Wuxi in China. Sustainability 2022, 14, 11426. https://doi.org/10.3390/su141811426

AMA Style

Zhang J, Chen Y, Yang X, Qiao W, Wang D. The Demarcation of Urban Development Boundary Based on the Maxent-CA Model: A Case Study of Wuxi in China. Sustainability. 2022; 14(18):11426. https://doi.org/10.3390/su141811426

Chicago/Turabian Style

Zhang, Jiaying, Yi Chen, Xuhong Yang, Wenyi Qiao, and Danyang Wang. 2022. "The Demarcation of Urban Development Boundary Based on the Maxent-CA Model: A Case Study of Wuxi in China" Sustainability 14, no. 18: 11426. https://doi.org/10.3390/su141811426

APA Style

Zhang, J., Chen, Y., Yang, X., Qiao, W., & Wang, D. (2022). The Demarcation of Urban Development Boundary Based on the Maxent-CA Model: A Case Study of Wuxi in China. Sustainability, 14(18), 11426. https://doi.org/10.3390/su141811426

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