Next Article in Journal
Motivations, Behaviors and Expectancies of Sexting: The Role of Defensive Strategies and Social Media Addiction in a Sample of Adolescents
Previous Article in Journal
Differences between Self-Report and Biomarkers in Smoking Indicators: The Necessity of Biomonitoring in Global Surveillance
Previous Article in Special Issue
Evaluation of Human Settlement Environment and Identification of Development Barriers Based on the Ecological Niche Theory: A Case Study of Northern Shaanxi, China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimization of the Territorial Spatial Patterns Based on MOP and PLUS Models: A Case Study from Hefei City, China

1
School of Economics and Management, Anhui Agricultural University, Hefei 230036, China
2
Institute of Land and Resources, Anhui Agricultural University, Hefei 230036, China
3
Wuxi Forestry Station, Wuxi 214000, China
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2023, 20(3), 1804; https://doi.org/10.3390/ijerph20031804
Submission received: 6 December 2022 / Revised: 13 January 2023 / Accepted: 16 January 2023 / Published: 18 January 2023

Abstract

:
Optimization of the territorial spatial patterns can promote the functional balance and utilization efficiency of space, which is influenced by economic, social, ecological, and environmental factors. Consequently, the final implementation of spatial planning should address the issue of sustainable optimization of territorial spatial patterns, driven by multiple objectives. It has two components—the territorial spatial scale prediction and its layout simulation. Because a one-sided study of scale or layout is divisive, it is necessary to combine the two to form complete territorial spatial patterns. This paper took Hefei city as an example and optimized its territorial spatial scale using the multiple objective programming (MOP) model, with four objective functions. A computer simulation of the territorial spatial layout was created, using the patch-generating land use simulation (PLUS) model, with spatial driving factors, conversion rules, and the scale optimization result. To do this, statistical, empirical, land utilization, and spatially driven data were used. The function results showed that carbon accumulation and economic and ecological benefits would be ever-increasing, and carbon emissions would reach their peak in 2030. The year 2030 was a vital node for the two most important land use types in the spatial scale—construction land and farmland. It was projected that construction land would commence its transition from reduced to negative growth after that time, and farmland would start to rebound. The simulation results indicated that construction land in the main urban area would expand primarily to the west, with supplemental expansion to the east and north. In contrast, construction land in the counties would experience a nominal increase, and a future ecological corridor would develop along the route south of Chaohu County–Chaohu Waters–Lujiang County–south of Feixi County.

1. Introduction

Economic and social development has always affected territorial spatial patterns. It is a major strategic deployment of China to establish a scientific and efficient territorial spatial planning system. In 2019, Several Opinions on Establishing and Supervising the Implementation of the Territorial Spatial Planning System was issued by the Central Committee of the Communist Party of China and the State Council. It has stressed the need to integrate the various plannings such as the main functional area planning, urban and rural planning, and land use planning into a unified territorial spatial planning to realize the “multiple-planning integration”. According to research on various types of planning, such as regional planning [1], urban planning [2], land use planning [3], landscape planning [4], habitat management planning [5], etc., the spatial planning system exhibits a good development history in the international arena. However, these spatial plannings are relatively independent, so there are many works of literature to discuss their convergence or integration. For example, Lopes et al. thought that city planning encompasses disciplines related to socio-economic, land-use, transport, environment, and others, but these disciplines face communication difficulties and objectives divergence due to contradictory interests and isolated evolution [6]. Sangawongse et al. took Wat Ket, Chiang Mai, Thailand as an example to study its transition between planning modes, from centralized planning to collaborative urban land use planning [7]. Kaczmarek et al. developed a machine learning approach for the integration of spatial development plans based on natural language processing [8]. In China, territorial spatial planning represents a significant innovation in terms of spatial planning. Due to the relative independence and even contradictions among various spatial planning strategies, their articulation or integration has been explored widely. Thus, the focus of territorial spatial planning is on dealing with the problems of excessive planning types, overlapping and even conflicting planning contents. Consequently, the final implementation of territorial spatial planning should be embodied in the sustainable optimization of territorial spatial patterns driven by multiple objectives [9].
Studies concerning the territorial spatial pattern are mainly categorized into “production–living–ecological space” and “ecological–agricultural–urban zone”, but the focus is still on the scale [10] and layout [11] in the spatial expression. In other words, the scale structure and spatial layout of territorial utilization serve as the basis for further implementation of “production–living–ecological space” [12] and “ecological–agricultural–urban zone” [13]. The optimization of the territorial spatial scale is primarily oriented to economic, social, ecological, and environmental objectives. The scale structure of territorial utilization is deduced by using regression prediction, Markov prediction, neural network, system dynamics, and other methods. While the single-objective oriented settings are mainly low carbon [14], ecosystem service [10], land suitability [15], etc., the dual-objective oriented settings include land capability and suitability [16], economic benefits and ecosystem services [17], eco-environmental constraints [18], etc. Multi-objective settings with more than three objectives are rarely found. A variety of spatial models are used to conduct a computer simulation in time and space to optimize the territorial spatial layout by setting space conversion rules usually categorized as two types. One is parallelism, which uses many juxtaposition rules based on the current spatial configuration [19]. The other is the progressive rules that simulate in steps, and the former step is the necessary condition for the latter [20]. Currently, the prevailing spatial models include Agent [21], CLUE-S (the Conversation of Land Use and its Effects at Small regional extent) [22], CA (Cellular Automata) [11], and other improved models based on them. There are many articles that address scale optimization based on objective realization, which is relatively easy to operate when the spatial layout is not considered. In recent years, increasing numbers of papers on spatial simulation have appeared, but most set transformation rules based on the spatial status quo. In spatial planning, however, the scale and layout of any one moment are interrelated, so a one-sided study of scale or layout is divisive to some extent.
Planning provides vision and scientific guidance for regional future development, and the adjustment of spatial scale and layout should meet more comprehensive developmental objectives as far as possible. At the same time, good planning should be combined closely with current natural, economic, and social conditions and avoid extensive demolition and construction. This paper, therefore, carried out a multi-objective prediction for the territorial spatial scale and then combined the objective prediction results with spatial driving factors based on the GIS (Geographic Information Systems) platform, to optimize the territorial spatial layout. In the setting of multiple-objective functions, minimum emission and maximum accumulation of carbon were added to the conventional goals of maximum economic and ecological benefits, considering the “double carbon” target of peak and neutrality. In the setting of spatial driving factors, POI datasets that could characterize geographical entities were added to the conventional nature and location traffic conditions. The research results can therefore provide an integrated concept for the study of scale and layout in spatial planning and are a reference for planners to adjust regional scale and layout.
For the core organization of the article, Section 3.1. is the optimization of the territorial spatial scale, and Section 3.2. is the optimization of the territorial spatial layout. Section 3.1.1. is the optimization schemes of four functions. Section 3.1.2. gives the setting of constraint conditions in the optimization schemes. Section 3.1.3. is the scale optimization result obtained by using LINGO12 software based on MOP model. Section 3.2.1. is spatial driving factors, including three types of basic driving factors and nine POI datasets. Section 3.2.2. is the accuracy test of PLUS model. Section 3.2.3 is the simulation result of space layout obtained by using ArcGIS10.2 software based on PLUS model.

2. Materials and Methods

2.1. Study Area

As the capital of Anhui Province, Hefei is one of the four major cities of science and education in China and a sub-center city of the Yangtze River Delta. Hefei reached an urbanization rate of 82.28%, with GDP exceeding 1 trillion yuan in 2020, ranking among the top 20 cities in China in terms of the two indicators for the first time. According to the seventh national census, the resident population of urban areas in Hefei surpassed 5 million in 2020, advancing into the ranks of megacities in China. However, a previous study suggested that the construction land in Hefei grew from 113,000 to 207,100 ha from 1997 to 2012, with total carbon emissions increasing from 4,462,100 to 16,567,300 tons [23]. With its rapid economic and social development, the city faces the serious challenge of optimizing land use structure, maintaining ecosystem balance, and protecting resources and the environment. In particular, the situation raises a great challenge for achieving the “double carbon” goals (peak carbon emissions before 2030 and reach carbon neutrality before 2060) and territorial spatial planning. Hefei City not only has to keep up the strong impetus of rapid economic development but must also consider ecological benefits and strive to achieve carbon peak and carbon neutrality. As a city facing all these challenges, Hefei can be considered a strong candidate for requiring optimization of territorial spatial patterns, which makes it the subject of our study. The location map of Hefei is shown in Figure 1.

2.2. Data Source

This paper uses statistical, empirical, land utilization, and spatially driven data for analysis, derived from the government yearbook, and the research of eminent institutions, platforms, and scholars. Socio-economic and energy data, obtained mainly from Anhui Statistical Yearbook, Hefei Statistical Yearbook, and China Energy Statistical Yearbook from 2011 to 2019, form the statistical data. Empirical data comprising coefficients of carbon emission, carbon accumulation, ecological benefits, and economic benefits are based on the published research results. The land utilization data consists of the remote sensing images of Hefei City in 2010 and 2018, acquired from the Resource and Environment Science and Data Center of the Chinese Academy of Sciences (http://www.resdc.cn (accessed on 30 November 2022)). The data with a resolution of 30 m were classified into six categories following the Chinese Academy of Sciences classification system: farmland, forest land, grassland, water bodies, construction land, and unused land. The 2010 remote sensing image was derived from Landsat-TM/ETM, and the 2018 image from Landsat 8, with less than 10% cloud cover. The accuracy of human–computer interactive interpretation was more than 90%, which met research requirements. The spatially driven data mainly covered the elevation and slope of the land, distances from the highways, town roads, railroads, district and county headquarters, and other points of interest (POI). Elevation data were obtained from Geospatial Data Cloud (http://www.gscloud.cn (accessed on 30 November 2022)); the slope was calculated from elevation; POI data were obtained from Amap; and other spatially driven data were obtained from OpenStreetMap.

2.3. Methods

Both the multiple objective programming (MOP) and patch-generating land use simulation (PLUS) models were used in the study. First, the territorial spatial scale was optimized using the MOP model with four objective functions—minimum carbon emission, maximum carbon accumulation, maximum economic benefits, and maximum ecological benefits. In the next step, computer simulation was carried out on the territorial spatial layout with the PLUS model, setting space conversion rules and using the scale optimization result as quantitative objectives. The MOP and PLUS models were implemented using LINGO12 and ArcGIS10.2, respectively.

2.3.1. MOP Model

The MOP model provides an optimal solution for multiple objectives within a given area. Generally, it comprises two or more objective functions, decision variables, and several constraining conditions, as shown in the following equation:
{ m i n   F ( X c e ) = i = 1 6 X i × c e i m a x   F ( X c a / X e e / X e b ) = i = 1 6 X i × c a i / e e i / e b i i , j = 1 n a i j X i = ( , )   Z j ,       X i 0
where F(Xce), F(Xca), F(Xee), and F(Xeb) denote the four objective functions representing carbon emission, carbon accumulation, and ecological and economic benefits, respectively. The decision variable Xi collectively denotes the area of farmland, forest land, grassland, water bodies, construction land, and unused land. Entries cei, cai, eei, and ebi represent the coefficients of carbon emission, carbon accumulation, ecological benefits, and economic benefits of each land use type. aij is the constraint coefficient corresponding to the jth variable in the ith constraint condition, and Zj is the constraint value. In terms of the objective function and constraint condition setting, the Markov chain [24] and GM (1,1) [25] were adopted to predict some indicators.
Since there is more than one function in multiple objective programming, it is usually impossible to find a point where all the objective functions reach their maximum/minimum values. Therefore, in the specific solution, it is necessary to adopt a compromise scheme to make each objective function as large/small as possible. The fuzzy deviation method of LINGO12 software [26] can be used to fuzzy the objective functions, transform the multiple objective linear problem into a single objective, and obtain the optimal solution to the problem.

2.3.2. PLUS Model

The PLUS model is a CA-based simulation model of spatial land utilization with a promising application in optimizing the territorial spatial layout [20]. It contains two modules. The first is a rule-mining framework based on a land expansion analysis strategy that explores each land type’s expansion and driving factors individually through a random forest algorithm. The second is a CA module based on multi-type random patch seeds, incorporating a decreasing threshold mechanism.
1.
Random Forest Classification for Dual-State Decisions
P i , k d ( x ) = n = 1 M I ( h n ( x ) = d ) M
where P i , k d denotes the growth probability of cell i to covert from the original land use type to the target type k; the value of d is either 0 or 1, and a value of 1 indicates that there were other land use types that changed to land use type k, while 0 represents other transitions; x is a vector that consists of multiple driving factors; I(∙) is the indicative function of the decision tree set; hn(x) is the prediction type of the nth decision tree for vector x; and M is the total count of decision trees.
2.
CA Model Based on Multi-Type Random Patch Seeds
O P i , k d = 1 , t = { P i , k d = 1 × ( r × μ k ) × D k t i f   Ω i , k t = 0   a n d   r < P i , k d = 1 P i , k d = 1 × Ω i , k t × D k t a l l   o t h e r s
Ω i , k t = c o n ( c i t 1 = k ) n × n 1 × w k
D k t = { D k t 1 i f   | G k t 1 | | G k t 2 | D k t 1 × G k t 2 G k t 1 i f   0 > G k t 2 > G k t 1 D k t 1 × G k t 1 G k t 2 i f   G k t 1 > G k t 2 > 0
where O P i , k d = 1 , t denotes the growth probability of cell i to covert from the original land use type to the target type k at iteration time t; Ω i , k t denotes the neighborhood effects of grid cell i at time t; D k t denotes the inertia coefficient for land use type k at time t; r is a random value ranging from 0 to 1; μ k is the threshold to generate the new land use patches for land use type k; c o n ( c i t 1 = k ) represents the total number of cells occupied by land use type k at the t − 1th iteration within the n × n window; w k is the variable parameter ranging from 0 to 1 among the different land use types, and the closer it is to 1, the stronger the expansion ability of the land use type; and G k t 1 and G k t 2 are the differences between the current amount of, and future demand for, land use type k at the t − 1th and t − 2th iteration.
To control the generation of multiple land use patches, a threshold descending rule for the competition process was proposed to restrict both the organic and spontaneous growth of all land use types.
If   k = 1 N | G c t 1 | k = 1 N | G c t | < S t e p   T h e n , l = l + 1
{ C h a n g e   P i , c d = 1 > τ     a n d     T M k , c = 1 N o   c h a n g e   P i , c d = 1 > τ   a n d   T M k , c = 0   τ = δ l × r 1
where Step is the step size of the PLUS model to approximate the land use demand; δ is the decay factor of decreasing threshold τ, which ranges from 0 to 1 and is set according to the expert; r1 is a normally distributed stochastic value with a mean of 1, ranging from 0 to 2; l is the number of decay steps; and T M k , c is the transition matrix that defines whether land use type k is allowed to convert to type c, and a value of 1 indicates permission, while 0 indicates no permission.

3. Results

3.1. Optimization of the Territorial Spatial Scale

3.1.1. Optimization Schemes

Following the timing of the current territorial spatial planning, optimization schemes for four objective functions were developed, choosing 2025, 2030, and 2035 as target years. Scenarios developed for carbon emission/accumulation as well as ecological/economic benefits are as follows. Their coefficients are shown in Table 1.

Minimum Carbon Emission

It is calculated by multiplying the area of each type of land by its carbon emission coefficient. The carbon emission/absorption ability of farmland, forest land, grassland, water bodies, and unused land, which have relatively stable carbon source/carbon sink capacities, is influenced mainly by changes in their areas. However, the carbon source capacity of the land for construction activities is significantly affected by human activities. As expected, the carbon emission coefficient changes with economic and social development. Thus, the minimum carbon emission level was decided based on results from Yu and Tian’s research [27] that were based on land use activities in Hefei City.

Maximum Carbon Accumulation

Based on Zhao et al. [28], the carbon accumulation coefficient of each type of land was obtained by calculating the total amount of accumulation and dividing it by the area of the land. The land’s relatively stable carbon accumulation coefficient without any construction activities was set as a fixed value, while the accumulation coefficient of construction land for each target year was obtained using the GM (1,1) model based on the calculated historical value.

Maximum Ecological Benefits

We used the method of Xie et al. [29] to calculate the ecological benefit coefficient of each land type for the target years. Based on the value of the world’s ecosystem services put forward by Costanza et al. [30], Xie et al. surveyed 700 Chinese professionals with ecological backgrounds. Based on their observations, they developed a set of evaluation methods for ecosystem services matching the conditions of China, which is what we followed.

Maximum Economic Benefits

The economic benefits of farmland, forest land, water bodies, and construction land were respectively represented by the output values of the planting industry, forestry, fisheries, and other secondary and tertiary industries. Given the difficulty to determine the economic benefits of grassland, the research results of Feng et al. [31] were used as a reference. Moreover, the minimum value of 0.001 ten thousand CNY/ha·year was taken for unused land that did not generate economic benefits. Then, the GM (1,1) model was used to obtain the economic benefit coefficients for the target years, based on the historical data of economic benefits per unit area in Hefei City from 2010 to 2018.

3.1.2. Constraints

The perspectives of area, development of society, ecology, and economy were considered for the optimization schemes to set the various constraints and are itemized and discussed below.
In the absence of administrative divisions for the land, the sum of the area of all land types should be equal to the total area. In other words, the total area of Hefei City should be S = X1 + X2 + X3 + X4 + X5 + X6, where X1 to X6 represents various land use types, same as Xi in Equation (1).
In the current period of rapid urbanization, the area of farmland will not exceed what is considered in the base period within a short time. Hence, the farmland area in 2018, considered the base period, was set as the upper limit. The base area in the target year predicted by Markov, based on historical data, was selected as the lower limit. Thus, the farmland area in the target year is predicted by Markov ≤ X1 ≤ the farmland area in 2018.
The forest land significantly contributes to reducing carbon emissions, increasing carbon accumulation and ecological benefits. Thus, its area in 2018, considered the base period, was set as the lower limit. There is an upper limit of a 10% increase of the current area because of its small increment, i.e., the forest land area in 2018 ≤ X2 ≤ current area × 110%.
Because of the low demand for grassland and its slow decline in the past years, the constraint of grassland area was set by referring to that of farmland area, i.e., the grassland area in the target year predicted by Markov ≤ X3 ≤ the grassland area in 2018.
The area covered by water shows stable fluctuations, within 0.5%. Thus, the constraint was formulated by increasing or decreasing the current area by 0.5%, i.e., the current area × 99.5% ≤ X4 ≤ the current area × 100.5%.
The construction land area for each target year was formulated based on previous research results by Yu and Huang [25]. Specifically, the construction land areas, X5, in 2025, 2030, and 2035 were designated as 227,297.70, 231,223.59, and 225,642.60 hectares, respectively.
Similar to grassland, the constraint was set as the unused land area in the target year predicted by Markov ≤ X6 ≤ the unused land area in 2018.
According to the population sizes and urbanization rates over the years, the total population, the population density of non-construction land, and the population density of construction land in each target year were predicted by the GM (1,1) model to develop the constraint on the population. This is defined as the population density of non-construction land × (X1 + X2 + X3 + X4 + X6) + population density of construction land × X5 ≤ total population.
For grain security, the grain demand per capita and the multiple-cropping index were set as 0.3925 tons/person [32] and 1.57 [33], respectively. Under the premise that no significant changes are expected in agricultural technology, the grain production capacity in each target year was predicted by the GM (1,1) model. The average values of 68.03% and 67.24% over the years were taken as the grain self-sufficiency rate and the ratio of the planting areas of food crops and cash crops, respectively. Thus, the constraint on grain demand was formulated as total population × 0.3925 × 68.03% ≤ X1 × grain production capacity × 67.24% × 1.57.
Considering the effectiveness of green equivalent in evaluating the quality of the regional ecological environment, green equivalent ≥ 1 was set as the constraint [34]. After the calculation of the ecological service values of farmland, forest land, and grassland, the green equivalents of farmland (GEfa), forest land (GEfo), and grassland (GEg) were calculated, respectively. Moreover, the target-year forest coverage rate (FCR) was predicted based on the target value and the increment in the National Economic and Social Development Planning of Hefei City. Thus, the constraint was set as (GEfa × X1+ GEfo× X2 + GEg × X3)/S × FCR ≥ 1.
For sustainable social and economic development, the ecological and economic benefits were set to follow a yearly growth trend.

3.1.3. Results of Scale Optimization

For the optimization scheme of a single objective function, different schemes will result in different territorial spatial scales (Figure 2).
Synthesizing these single objectives to make each function as large/small as possible, as mentioned in Section 2.3.1., the territorial spatial scale optimization results for the Hefei City for 2025, 2030, and 2035 were obtained using the fuzzy deviation method in the LINGO12 software (Table 2).
Overall, an upward trend in the economic benefits was predicted, consistent with the current strategic positioning and rate of development of the city. Ecological benefits and carbon accumulation would show steady growth, and the peak values of carbon emissions were estimated to appear in 2030. A positive development trend was found in all the land use structures, carbon sources/sinks, and ecological benefits, indicating that relevant national policies were actively implemented. To be more specific, the forest land area in the ecological land would keep increasing; the grassland area would basically stabilize after a temporary and gentle decline. The area covered by water was estimated to fluctuate and show relatively smaller growth. The unused land would experience a low base but an insignificant reduction. Correspondingly, carbon accumulation and ecological benefits would exhibit an upward trend, reflecting the contribution of Hefei City to the construction of ecological civilization in the new era. During the year 2030, an important node in the model, the construction land area was predicted to first grow in decreasing increments and later show a negative growth rate. However, in the meantime, the area of farmland would begin to rebound, and carbon emissions would reach their peak.
On the one hand, there would be a balance between land under construction and farming by ensuring the healthy and rapid development of the economy and society. The indices of growth would reflect the improvements, both in the efficient use of construction land and the increase in farmland productivity. On the other hand, it would witness the maximum carbon emissions in 2030 but will actively move towards the higher carbon neutrality target by 2060.

3.2. Optimization of the Territorial Spatial Layout

3.2.1. Spatial Driving Factors

Nature and location traffic conditions constitute the basic foundation of the scope, spatial distribution, and spatial expansion of a region. Based on previous studies, three types of basic driving factors have been generally identified: elevation, slope, and proximity factors [35,36,37,38]. In addition to the elevation and slope, distance to the highways, town roads, railroads, and district and county headquarters were used to characterize proximity factors (Figure 3).
A point of interest (POI) is a point location with name, address, and category attributes; all geographical entities that can be abstracted as points can be referred to in terms of POIs. POI datasets can thus provide useful socioeconomic information for significant locations, especially in relation to hospitals, schools, banks, shopping malls, and other facilities closely related to people’s lives. In addition, the introduction of dynamic POI data can effectively improve the accuracy of land use simulation [38]. Based on the development status and the changing trend of Hefei City, we select nine categories of POI data, including medical institutions, education institutions, financial institutions, commercial buildings, public service facilities, entertainment facilities, food and beverage outlets, transport sites, and living service facilities. Figure 4 reflects the density distribution of these POIs after normalization.

3.2.2. Accuracy Testing

The Kappa coefficient or the figure of merit (FOM) coefficient is commonly used to test simulation accuracy. The value of the former lies between 0 and 1, and the closer the value to 1, the higher the simulation accuracy. While a Kappa coefficient larger than 0.7 is acceptable, values larger than 0.8 are considered satisfactory [19]. The value of the FOM coefficient also ranges between 0 and 1, and a larger value suggests higher simulation accuracy. In general, the FOM coefficient of about 0.15 means that the test is reliable [39]. For improved validity and reliability of the simulation, we used both coefficients to test the accuracy of the results. According to the territorial spatial pattern of Hefei City in 2010, the simulation was conducted under the settings described above for 2018 and then compared with the situations in the same year. The Kappa coefficient of 0.9235 and the FOM coefficient of 0.2487 obtained in our calculations suggest high and acceptable simulation accuracy. Based on the territorial spatial pattern for the year 2018, a simulation was conducted to generate the territorial spatial layout of Hefei City for three years—2025, 2030, and 2035, and the results are discussed next.

3.2.3. Results of Simulation Optimization

Figure 5 is the visualization result of spatial layout optimization based on all previous studies. It meets the function equation, transformation rules, and other research goals set above: i.e., these equations and rules can be used for setting during the planning process, to ensure that implementation follows the optimal direction.
Our results show that during 2020–2025 and 2025–2030, the main urban area of Hefei City will continue to display a high demand for construction land. The models also predict that the city’s expansion would mainly be directed westward, with some additional activities to the east and north. This pattern matches the proposed long-term development plan of the city. It is consistent with the proposed expansion of the Xinqiao Science and Technology Innovation Demonstration Zone of Shushan District, which envisages the concentration of large scientific installations. It incorporates the New Canal City and the Feixi Industrial City Integration Demonstration Zone as the pivots. Growth in the northern wing would be around the Xiatang Industrial New City and Beicheng New District of Changfeng County. To the east, the city would expand, with the Feidong Industrial New District as its pivot.
The spatial coverage of construction land in Lujiang County, Chaohu County, and the central and northern Changfeng County would undergo a small incremental growth, mainly concentrated around the county headquarters. The incremental growth in the forest land areas was predicted to be mainly focused on Lujiang and Chaohu counties, with the future functional activities dominated by ecological space. Furthermore, the forest land areas of western Shushan District and eastern Feixi County would see a nominal rise because of the existence of Dashu Mountain National Forest Park and Zipeng Mountain National Forest Park. As mentioned previously, 2030 would be a prominent node due to various significant changes. During this period, carbon emissions are expected to peak. The land area growth under construction activities would change from a reduced state to negative values, while the farmland area is expected to show a nominal increase. Eventually, with all these changes, the territorial spatial layout of the city would enter a stage of steady consolidation. During 2030–2035, the land area under construction in the main urban area would basically remain stable, with a small reduction in the southern Lujiang and Chaohu counties and the nearby Zipeng Mountain in Feixi County, which is dominated by ecological space. Based on the above simulation results, the areas extending from the southern Chaohu County to the south and west, connecting with the Chaohu Waters and Lujiang County and southern Feixi County, would form an ecological corridor of Hefei City as a whole.

4. Discussion

Although the scope of government power at varying organizational levels of unitary, federal countries may differ in detail, they generally have similar planning logic and execution frameworks [1,2,3,4,5]. There are many types of spatial planning, however, and this can lead to contradictions, and even conflicts, in the implementation process [6,7,8], as some departments may not consider the interests of others in their plans. The innovation of China’s territorial spatial planning lies in the construction of a multiple-planning integration system, which would integrate economic and social development goals and balance the interests of different departments. However, the public’s greatest concern is how the plan will be implemented, and whether it will promote economic and social development. A specific spatial pattern supported by spatial scale [10] and layout [11] must therefore be in place, especially at the implementation stage. There are many research results on spatial scale optimization, but one or two factors, among economy, ecology, environment, and low carbon, are often chosen in the selection of optimization objectives. Although it reflects pertinence, it is difficult to reflect comprehensiveness. Spatial layout optimization is usually dominated by the setting of spatial factors, which is mostly a direct simulation based on the spatial status quo. It is difficult to fully reflect the statistical objectives of economic and social development through the spatial layout simulation. There are, however, few in-depth studies that combine the two, despite the fact that a one-sided study of scale or layout is divisive, because the two coexist and echo each other; it is, therefore, necessary to combine them to form complete territorial spatial patterns.
This paper put forward an overall optimization idea covering spatial scale and spatial layout. Specifically, compared to the common single-objective [10,14,15] or dual-objective [16,17,18] oriented settings, more comprehensive objective functions were scientifically set to optimize the scale of territorial use in a broad and effective manner. In the next step, the parallelism rule [19] and progressive rule [20] were integrated to provide a reference for the simulation of the territorial spatial layout based on multiple objective constraints, basic driving factors, and POI driving factors. In other words, various modules to be considered during territorial spatial planning and implementation were explored to achieve the various functions. These functions included minimum carbon emission, maximum carbon accumulation, maximum economic benefit, and maximum ecological benefit. Whether the objective functions are realized is a matter for post evaluation: i.e., the effects of the goals achieved when the planning reaches a certain stage of implementation. This implementation process will be directly reflected in the evolution of the spatial layout, but the evolution of this layout must be based on the current situation and is limited by factors that include nature, location of transportation, and built facilities. In addition to clarifying the conversion rules, spatial driving factors should therefore be included. Elevation and slope are natural geographical factors, while proximity factors are mainly location traffic factors, which measure the distance between a block and main roads and nodes. POI data can be obtained in real-time through Amap, and all geographical entities that can be abstracted as points can be referred to in terms of POIs. This paper selects facilities that are closely related to people’s lives, because planning is a reflection of human behavior, and both development and conservation should be centered on people. Perhaps, for this reason, the development of international spatial planning has continued to focus on human settlements.
As territorial spatial planning focuses on an all-encompassing structure that incorporates land use, urban and rural planning, as well as other developmental activities, it requires multiple levels of planning and integration. Based on the results of available research and with the prospects of increased data availability in the future, the introduction of more constraints or conversion rules is worth trying to improve the results of the optimization. For example, a criterion for the evaluation of resource and environmental carrying capacity [40] could be introduced in spatial scale optimization. Suitability evaluation of territory space [41] in the simulation of spatial layout could be another useful step that can be introduced. Spatial driving factor selection, the introduction of space management and control boundaries of permanent basic farmland [38], and the ecological protection red line [42], following the space conversion rules, are the other inputs that will contribute to more scientific and efficient planning and a clearer implementation path. In addition, once planning has been implemented to a certain stage, it can be evaluated according to the method used in this paper. Through the implementation evaluation, the phased results of planning implementation can be compared with prior optimization results, and the difference can be evaluated, so as to realize the mutual transformation of theory and practice.

5. Conclusions

In this paper, four objective functions were developed to realize multi-objective optimization of territorial spatial scale in Hefei City using the MOP model. On this basis, space conversion rules were developed to perform computer simulations of the territorial spatial layout via the PLUS model. The main conclusions are summarized below.
According to the results of scale optimization obtained by the MOP model, the area of forest land, a significant ecological land, would continue to grow with a large increment. However, the area of unused land and those covered by grassland and water bodies would reach a stable state with some fluctuations with small changes. Regarding the two critical land types, i.e., the farmland and construction land expanding to farmland, the year 2030 would be an important node. After 2030, the area of construction land would change from decreased to negative growth, while the net area of farmland would start to increase. These predictions demonstrate that while guaranteeing sustainable economic and social development, Hefei City will also meet other favorable changes in ecological parameters as well as improve the carbon footprint. With the creation of more carbon sinks and with higher carbon accumulation, the ecological benefits would increase steadily. The simulations show that carbon emissions will peak in 2030.
As revealed by the simulation results obtained by the PLUS model, from 2020 to 2030, construction land in the main urban area of Hefei City would mainly expand westward and then towards the east and the north. However, at the county level, the area of construction land would undergo only nominal growth. After 2030, the territorial spatial layout of the city would step into a period of steady consolidation. Lujiang and Chaohu counties would develop as important ecological spaces, and future years will see the growth of an ecological corridor to the south of Chaohu County–Chaohu Waters–Lujiang County–south of Feixi County.

Author Contributions

Conceptualization, R.Y.; methodology, R.Y. and H.C.; software, H.C.; validation, S.F., T.L., C.W. and Y.S.; formal analysis, H.C., Y.Y., and Q.W.; investigation, H.C., Y.Y. and Q.W.; data curation, H.C.; writing—original draft preparation, R.Y. and H.C.; writing—review and editing, R.Y.; supervision, R.Y.; project administration, R.Y. and H.C.; funding acquisition, H.C., S.F., T.L., C.W. and Y.S.; Resources, X.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Postgraduate Scientific Research Project of Education Department of Anhui Province (grant number YJS20210255). This research was supported by National Natural Science Foundation of China (grant number 71873003, 32201346, 42201281) and Natural Science Foundation of Anhui Province (grant number 1908085QG310, 2208085QD102).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available upon reasonable request from the corresponding author. The data are not publicly available due to consent provided by participants.

Acknowledgments

We are grateful for the anonymous reviewers and the editor for their constructive comments and suggestion for this paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cattivelli, V. Planning peri-urban areas at regional level: The experience of Lombardy and Emilia-Romagna (Italy). Land Use Policy 2021, 103, 105282. [Google Scholar] [CrossRef]
  2. Fox, D.M.; Carrega, P.; Ren, Y.; Caillouet, P.; Bouillon, C.; Robrt, S. How wildfire risk is related to urban planning and Fire Weather Index in SE France (1990–2013). Sci. Total Environ. 2018, 621, 120–129. [Google Scholar] [CrossRef] [PubMed]
  3. Maleki, J.; Masoumi, Z.; Hakimpourc, F.; Coello, C.A. A spatial land-use planning support system based on game theory. Land Use Policy 2020, 99, 105013. [Google Scholar] [CrossRef]
  4. Sánchez, M.L.; Cabrera, A.T.; Gómez Del Pulgar, M.L. Guidelines from the heritage field for the integration of landscape and heritage planning: A systematic literature review. Landsc. Urban Plan. 2020, 204, 103931. [Google Scholar] [CrossRef]
  5. Mukherjee, T.; Sharma, V.; Sharma, L.K.; Thakur, M.; Joshi, B.D.; Sharief, A.; Thapa, A.; Dutta, R.; Dolker, S.; Tripathy, B.; et al. Landscape-level habitat management plan through geometric reserve design for critically endangered Hangul (Cervus hanglu hanglu). Sci. Total Environ. 2021, 777, 146031. [Google Scholar] [CrossRef] [PubMed]
  6. Lopes, A.S.; Cavalcantea, C.B.; Vale, D.S.; Loureiro, C.F.G. Convergence of planning practices towards LUT integration: Seeking evidences in a developing country. Land Use Policy 2020, 99, 104842. [Google Scholar] [CrossRef]
  7. Sangawongse, S.; Fisher, R.; Prabudhanitisarn, S. From centralised planning to collaborative urban land use planning: The case of Wat Ket, Chiang Mai, Thailand. Soc. Sci. Humanit. Open 2021, 4, 100154. [Google Scholar] [CrossRef]
  8. Kaczmarek, I.; Iwaniak, A.; Świetlicka, A.; Piwowarczyk, M.; Nadolny, A. A machine learning approach for integration of spatial development plans based on natural language processing. Sustain. Cities Soc. 2022, 76, 103479. [Google Scholar] [CrossRef]
  9. Jia, K.J.; He, H.F.; Zhang, H.; Guo, J. Optimization of territorial space pattern based on resources and environment carrying capacity and land suitability assessment. China Land Sci. 2020, 34, 43–51. [Google Scholar]
  10. Wu, X.T.; Wang, S.; Fu, B.J.; Liu, Y.; Zhu, Y. Land use optimization based on ecosystem service assessment: A case study in the Yanhe watershed. Land Use Policy 2018, 72, 303–312. [Google Scholar] [CrossRef]
  11. Newland, C.P.; Maier, H.R.; Zecchin, A.C.; Newman, J.P.; van Delden, H. Multi-objective optimisation framework for calibration of Cellular Automata land-use models. Environ. Model. Softw. 2018, 100, 175–200. [Google Scholar] [CrossRef]
  12. Zhao, H.B.; Wei, J.C.; Sun, D.Q.; Liu, Y.X.; Wang, S.; Tan, J.T.; Miao, C.H. Recognition and spatio-temporal evolution analysis of production-living-ecological spaces based on the random forest model: A case study of Zhengzhou city, China. Geogr. Res. 2021, 40, 945–957. [Google Scholar]
  13. Ke, X.L.; Xiao, B.Y.; Zheng, W.W.; Ma, Y.C.; Li, Y.H. Urban-agricultural- ecological space zoning based on scenario simulation. J. Geo. Inf. Sci. 2020, 22, 580–591. [Google Scholar]
  14. Wang, G.Z.; Han, Q.; de Vries, B. The multi-objective spatial optimization of urban land use based on low-carbon city planning. Ecol. Indic. 2021, 125, 107540. [Google Scholar] [CrossRef]
  15. Akbaria, M.; Neamatollahib, E.; Neamatollahi, P. Evaluating land suitability for spatial planning in arid regions of eastern Iran using fuzzy logic and multi-criteria analysis. Ecol. Indic. 2019, 98, 587–598. [Google Scholar] [CrossRef]
  16. Feudis, M.D.; Falsone, G.; Gherardi, M.; Speranza, M.; Vianello, G.; Antisari, L.V. GIS-based soil maps as tools to evaluate land capability and suitability in a coastal reclaimed area (Ravenna, northern Italy). Int. Soil Water Conserv. Res. 2021, 9, 167–179. [Google Scholar] [CrossRef]
  17. Ma, S.H.; Wen, Z.Z. Optimization of land use structure to balance economic benefits and ecosystem services under uncertainties: A case study in Wuhan, China. J. Clean. Prod. 2021, 311, 127537. [Google Scholar] [CrossRef]
  18. Luo, X.; Lu, X.H.; Jin, G.; Wan, Q.; Zhou, M. Optimization of urban land-use structure in China’s rapidly developing regions with eco-environmental constraints. Phys. Chem. Earth 2019, 110, 8–13. [Google Scholar] [CrossRef]
  19. Liu, X.P.; Liang, X.; Li, X.; Xu, X.C.; Ou, J.P.; Chen, Y.M.; Li, S.Y.; Wang, S.J.; Pei, F.S. 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]
  20. Liang, X.; Guan, Q.F.; Clarke, K.C.; Liu, S.S.; Wang, B.Y.; Yao, Y. Understanding the drivers of sustainable land expansion using a patch-generating land use simulation (PLUS) model: A case study in Wuhan, China. Comput. Environ. Urban Syst. 2021, 85, 101569. [Google Scholar] [CrossRef]
  21. Chowdhury, S.; Peddle, D.R.; Wulder, M.A.; Heckbert, S.; Shipman, T.C.; Chao, D.K. Estimation of land-use/land-cover changes associated with energy footprints and other disturbance agents in the Upper Peace Region of Alberta Canada from 1985 to 2015 using Landsat data. Int. J. Appl. Earth Obs. Geoinf. 2021, 94, 102224. [Google Scholar] [CrossRef]
  22. Li, Y.J.; Ye, C.S.; Huang, X.L. Temporal-spatial evolution and scenario simulation of production-living-ecological space in Nanchang based on CLUE-S model. Res. Soil Water Conserv. 2021, 28, 325–332. [Google Scholar]
  23. Yu, R.; Huang, X.J.; Xu, G.L. Towards a Theory of Decoupling: The relationship of carbon emission and construction land expansion between 1997 and 2012 in Hefei City, China. Nat. Environ. Pollut. Technol. 2019, 18, 569–577. [Google Scholar]
  24. Feng, D.R.; Bao, W.K.; Fu, M.C.; Zhang, M.; Sun, Y.Y. Current and future land use characters of a national central city in eco-fragile region—A case study in Xi’an City based on FLUS model. Land 2021, 10, 286. [Google Scholar] [CrossRef]
  25. Yu, R.; Huang, X.J. Research on the expansion scale of construction land under the restriction of carbon emission peak. China Popul. Resour. Environ. 2019, 29, 66–72. [Google Scholar]
  26. Golfam, P.; Ashofteh, P.S.; Lo´aiciga, H.A. Modeling adaptation policies to increase the synergies of the water-climate-agriculture nexus under climate change. Environ. Dev. 2021, 37, 100612. [Google Scholar] [CrossRef]
  27. Yu, R.; Tian, S.M. An analysis of land use and carbon emission in Hefei city based on carrying relationship. J. Anhui Agric. Univ. 2016, 43, 939–945. [Google Scholar]
  28. Zhao, R.Q.; Huang, X.J.; Zhong, T.Y.; Chuai, X.W. Carbon effect evaluation and low-carbon optimization of regional land use. Trans. Chin. Soc. Agric. Eng. 2013, 29, 220–229. [Google Scholar]
  29. Xie, G.D.; Zhen, L.; Lu, C.X.; Xiao, Y.; Chen, C. Expert knowledge based valuation method of ecosystem services in China. J. Nat. Resour. 2008, 23, 911–919. [Google Scholar]
  30. Costanza, R.; d’Arge, R.; de Groot, R.; Farber, S.; Grasso, M.; Hannon, B.; Limburg, K.; Naeem, S.; O’Neill, R.V.; Paruelo, J.; et al. The value of the world’s ecosystem services and natural capital. Nature 1997, 387, 253–260. [Google Scholar] [CrossRef]
  31. Feng, C.C.; Cao, M.Z.; Xie, T.T. Optimization of land use structure in Tongling city based on different ecological conservation scales. Geogr. Res. 2014, 23, 2217–2227. [Google Scholar]
  32. Tang, H.J.; Li, M.Z. Study on per capita grain demand based on Chinese reasonable dietary pattern. Sci. Agric. Sin. 2012, 45, 2315–2327. [Google Scholar]
  33. Ding, M.J.; Chen, Q.; Xin, L.J.; Li, L.H.; Li, X.B. Spatial and temporal variations of multiple cropping index in China based on SPOT-NDVI during 1999–2013. Acta Geogr. Sin. 2015, 70, 1080–1090. [Google Scholar]
  34. Xu, J.; Gong, P.; Jin, L. Optimization evaluation and analysis on land use structure based on ecological green equivalent in Hohhot. Ecol. Econ. 2019, 35, 196–201. [Google Scholar]
  35. Dubovyk, O.; Sliuzas, R.; Flacke, J. Spatio-temporal modelling of informal settlement development in Sancaktepe district, Istanbul, Turkey. ISPRS J. Photogramm. Remote Sens. 2011, 66, 235–246. [Google Scholar] [CrossRef]
  36. Li, X.M.; Zhou, W.Q.; Ouyang, Z.Y. Forty years of urban expansion in Beijing: What is the relative importance of physical, socioeconomic, and neighborhood factors? Appl. Geogr. 2013, 38, 1–10. [Google Scholar] [CrossRef]
  37. Li, G.D.; Sun, S.A.; Fang, C.L. The varying driving forces of urban expansion in China: Insights from a spatial-temporal analysis. Landsc. Urban Plan. 2018, 174, 63–77. [Google Scholar] [CrossRef]
  38. Wu, R.; Wang, J.Y.; Zhang, D.C.; Wang, S.J. Identifying different types of urban land use dynamics using Point-of-interest (POI) and Random Forest algorithm: The case of Huizhou, China. Cities 2021, 114, 103202. [Google Scholar] [CrossRef]
  39. Pontius, R.G.; Boersma, W.; Castella, J.C.; Clarke, K.; de Nijs, T.; Dietzel, C.; Duan, Z.; Fotsing, E.; Goldstein, N.; Kok, K.; et al. Comparing the input, output, and validation maps for several models of land change. Ann. Reg. Sci. 2008, 42, 11–37. [Google Scholar] [CrossRef] [Green Version]
  40. Wu, Y.; Zong, T.; Shuai, C.Y.; Liao, S.J.; Jiao, L.D.; Shen, L.Y. Does resource environment carrying capacity have a coercive effect on urbanization quality? Evidence from the Yangtze River Economic Belt, China. J. Clean. Prod. 2022, 365, 132612. [Google Scholar] [CrossRef]
  41. Zhang, Z.; Li, J.M. Spatial suitability and multi-scenarios for land use: Simulation and policy insights from the production-living-ecological perspective. Land Use Policy 2022, 119, 106219. [Google Scholar] [CrossRef]
  42. Li, Z.J.; Liu, Y.M.; Zeng, H. Application of the MaxEnt model in improving the accuracy of ecological red line identification: A case study of Zhanjiang, China. Ecol. Indic. 2022, 137, 108767. [Google Scholar] [CrossRef]
Figure 1. Geographical location of Hefei City.
Figure 1. Geographical location of Hefei City.
Ijerph 20 01804 g001
Figure 2. Area of various types of land under single objective optimization in 2025, 2030, and 2035.
Figure 2. Area of various types of land under single objective optimization in 2025, 2030, and 2035.
Ijerph 20 01804 g002aIjerph 20 01804 g002b
Figure 3. Basic driving factors of nature and location traffic conditions.
Figure 3. Basic driving factors of nature and location traffic conditions.
Ijerph 20 01804 g003
Figure 4. The density of selected POI driving factors.
Figure 4. The density of selected POI driving factors.
Ijerph 20 01804 g004
Figure 5. Simulation results of territorial spatial layout.
Figure 5. Simulation results of territorial spatial layout.
Ijerph 20 01804 g005
Table 1. Coefficients of carbon emission/accumulation and ecological/economic benefits of various land use types in target years.
Table 1. Coefficients of carbon emission/accumulation and ecological/economic benefits of various land use types in target years.
Target YearFarmlandForest LandGrasslandWater BodyConstruction LandUnused Land
carbon emission coefficient
(t/ha)
20250.422−0.644−0.022−0.253169.900−0.005
2030333.191
2035332.136
carbon accumulation coefficient
(t/ha)
202563.015114.33143.96522.69089.43666.973
203098.693
2035103.198
ecological benefit coefficient
(ten thousand CNY/ha·year)
20250.5353.1782.7117.2880.0280.151
2030
2035
economic benefit coefficient
(ten thousand CNY/ha·year)
20254.3810.9821.9906.296377.0420.001
20307.1842.55311.6531362.647
20358.6953.82314.8692379.650
Table 2. Optimization results of territorial spatial scales and objective functions under multiple objective programming.
Table 2. Optimization results of territorial spatial scales and objective functions under multiple objective programming.
2018202520302035
farmland (ha)748,597.50723,702.87714,912.12717,800.85
forest land (ha)48,031.6549,539.1554,493.0257,181.50
grassland (ha)39,998.7039,593.2539,597.4839,592.62
water bodies (ha)105,536.34106,114.32106,021.08106,030.17
construction land (ha)204,083.01227,297.70231,223.59225,642.60
unused land (ha)24.4824.3924.3923.94
carbon emissions (ten thousand tons)3493.117578.117728.057518.24
carbon accumulation (ten thousand tons)7507.067694.087824.757920.27
ecological benefits (CNY one hundred million)249.18249.29250.23251.24
economic benefits (CNY one hundred million)8101.8718,278.6432,166.5454,506.56
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

Yu, R.; Cheng, H.; Ye, Y.; Wang, Q.; Fan, S.; Li, T.; Wang, C.; Su, Y.; Zhang, X. Optimization of the Territorial Spatial Patterns Based on MOP and PLUS Models: A Case Study from Hefei City, China. Int. J. Environ. Res. Public Health 2023, 20, 1804. https://doi.org/10.3390/ijerph20031804

AMA Style

Yu R, Cheng H, Ye Y, Wang Q, Fan S, Li T, Wang C, Su Y, Zhang X. Optimization of the Territorial Spatial Patterns Based on MOP and PLUS Models: A Case Study from Hefei City, China. International Journal of Environmental Research and Public Health. 2023; 20(3):1804. https://doi.org/10.3390/ijerph20031804

Chicago/Turabian Style

Yu, Ran, Hongsheng Cheng, Yun Ye, Qin Wang, Shuping Fan, Tan Li, Cheng Wang, Yue Su, and Xingyu Zhang. 2023. "Optimization of the Territorial Spatial Patterns Based on MOP and PLUS Models: A Case Study from Hefei City, China" International Journal of Environmental Research and Public Health 20, no. 3: 1804. https://doi.org/10.3390/ijerph20031804

APA Style

Yu, R., Cheng, H., Ye, Y., Wang, Q., Fan, S., Li, T., Wang, C., Su, Y., & Zhang, X. (2023). Optimization of the Territorial Spatial Patterns Based on MOP and PLUS Models: A Case Study from Hefei City, China. International Journal of Environmental Research and Public Health, 20(3), 1804. https://doi.org/10.3390/ijerph20031804

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