Next Article in Journal
Spatiotemporal Degradation of Abandoned Farmland and Associated Eco-Environmental Risks in the High Mountains of the Nepalese Himalayas
Next Article in Special Issue
Quantifying the Effect of Land Use Change Model Coupling
Previous Article in Journal
Amman (City of Waters); Policy, Land Use, and Character Changes
Previous Article in Special Issue
Projecting Urbanization and Landscape Change at Large Scale Using the FUTURES Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Simulation of an Urban-Rural Spatial Structure on the Basis of Green Infrastructure Assessment: The Case of Harbin, China

School of Architecture, Harbin Institute of Technology, Key Laboratory of Cold Region Urban and Rural Human Settlement Environment Science and Technology, Ministry of Industry and Information Technology, 66 West Dazhi Street, Nangang District, Harbin 150006, China
*
Author to whom correspondence should be addressed.
Submission received: 30 October 2019 / Revised: 12 December 2019 / Accepted: 13 December 2019 / Published: 15 December 2019
(This article belongs to the Special Issue Land Change Modelling)

Abstract

:
Due to their long-term dual structures and rapid urbanization, cities and villages in developing countries are undergoing the challenges of urban-rural integration and ecological security. This study aims to determine the pattern of urban-rural spatial structures under the circumstances of ecological security in the future to promote the integrated, coordinated, green, and sustainable development of urban-rural spaces. Using a quantitative evaluation method, the logistic-CA model, the LCP (least cost path) model, and a classification of ecological importance, this study constructed an integrated simulation model based on a green infrastructure assessment and applied the model to simulate and predict the urban-rural spatial structure of the Harbin city territory (Harbin) in 2035. The results indicate that the urban-rural hierarchical scale structure of Harbin comprises a central city, sub-central city, central town, major town, common town, central village, and general village. The urban-rural traffic network structure forms a pattern of “radiation + grid”, with Harbin city at the center of the structure. The urban-rural land use zoning structure consists of eco-spaces, agricultural spaces, and construction spaces. It can be concluded that in 2035, the urban-rural spatial structure of Harbin will show an increasing development tendency, where single-center, medium, and small cities in will Harbin develop, and traffic systems above the county level will also improve.

1. Introduction

The cities and villages of developed countries, like Britain, France, the USA, and Japan, have broadly evolved into a state of balance [1,2]. However, influenced by the extreme differences between urban and rural areas [3,4], the contradiction of urban-rural development in China remains prominent. Urban and rural areas are equally essential and should be combined organically [5]. Therefore, how to intensively make use of urban-rural resources and spaces, rationally optimize the overall layout of urban and rural areas, protect and rectify the urban-rural ecological environment, and coordinate the relationship between urban and rural development to promote the sustainable and healthy development of urban-rural society and economy, is one of main tasks confronting China at the moment. The urban-rural spatial structure, a regional reflection of humans’ social, economic, and cultural activities under specific circumstances, illustrates the spatial dichotomy between urban and rural areas. In addition, with the rapid development of urbanization, urban sprawl has brought great challenges to urban and rural eco-environments [6,7,8], especially in China [9,10]. How to protect and maintain ecosystem security during the process of urban growth, to resolve the contradiction between ecological protection and urban—rural development, is another major problem that is necessary to solve. Against this background, studies on urban-rural spatial structure face new challenges, that is, using the bottom-line thinking of taking “three areas (urban, agricultural, and ecological spaces) and three lines (the red line of the ecological areas, the red line of permanent basic farmlands, and the urban development boundary)” as impassable boundaries of spatial development while ensuring the national layout of urban-rural spaces. This new challenge would first require the land use in urban and rural areas to be arranged from the perspective of ecology, so that the urban-rural spatial structure can evolve toward sustainability and ecological security.
It is an historical idea that the urban spatial structure can be arranged from the perspective of ecology. The ancient Chinese “Feng Shui [11,12]” model and Western “garden construction [13,14]” model, Howard’s “Garden Cities”, Saarinen’s “Theory of Organic Decentralization”, and Soleri’s “Bionic City” in the late nineteenth and early twentieth centuries are representations of this idea. After the 1950s, the outbreak of the global environmental crisis made human beings rethink their relationships with nature. McHarg suggested that urban design should be connected with nature and described how to conduct land development via natural processes in his book [15] where his ideas initiated the exploration of urban ecological space. Thereafter, words like “Ecological Priority”, “Ecological Infrastructure”, and “Green Infrastructure” emerged and are now used in ecological urban studies [16,17,18]. Green infrastructure is the network of natural and semi-natural green and blue areas [19,20,21,22], which play an important role in the process of urban and rural sustainable development. Green infrastructure provides support for ecological processes and can contribute to better human health and wellbeing [18,23]. The application and protection of green infrastructure is beneficial to mitigate urban sprawl and promote the rational development of the urban spatial structure [24]. However, previous studies have mainly concentrated on discussing the optimization or reconstruction of urban spatial structures from an ecological perspective. It is a complex and long process to change or shape the spatial structure of a city, which is why some researchers have developed simulation studies on urban development in the future using Cellular Automata (CA) models [25,26,27]. The research objects in these studies are actually urban land expansion or urban interior spatial structures. The Lowry model can simulate the urban spatial structure in terms of land use and transport [28]. However, this model is limited to explaining the formation mechanism of current urban spatial structures. With reference to existing research [29,30], the urban-rural spatial structure of this study refers to the distribution and combination of the point-line-surface in space, where points are settlements, lines are traffic lines, and surfaces are the land use division corresponding to the hierarchical size structure, the traffic network structure, and the land use zoning structure, respectively.
In this study, in order to promote the sustainable development of urban and rural space, an integrated simulation model from an ecological perspective was proposed. This model primarily includes four parts: (1) a green infrastructure assessment; (2) hierarchical size structure extraction; (3) traffic network structure simulation; and (4) land use zoning structure determination. Using a quantitative evaluation method, the CA model, the Least Cost Path (LCP) model, and a classification of ecological importance, coupled with the spatial analysis technology of Geographic Information System (GIS), this study simulated the spatial structure of the Harbin city territory in 2035. The quantitative evaluation method was used to assess green infrastructure. The results of this green infrastructure will be an ecological constraint on the model, which is the fundamental basis of the integrated model. The CA model is used to simulate urban-rural growth, from which the hierarchical size structure can be extracted. The LCP model can simulate potential traffic paths, from which the traffic network structure can be obtained, combined with the present and planned traffic situations. The classification of ecological importance is the basis of the land use zoning, aiming to give priority to determining the scope of ecological spaces. It can be concluded that this integrated model has limitations, but can reflect the pattern of urban-rural spatial structures in the future to some extent. Furthermore, in 2035, the urban-rural spatial structure of Harbin will show an increasing development tendency toward a single-center, and traffic systems above the county level will also be improved.
This study offers two contributions to knowledge: (1) This study attempts to construct an integrated simulation model to simulate urban-rural spatial structures; and (2) this study offers an analysis of Harbin’s urban-rural spatial structure on a regional scale from an ecological perspective. Previous research on urban spatial simulations have focused on central Harbin city [31] or Harbin city [32]. The scope of this study will be the administrative area of Harbin, that is the Harbin city territory. This territory contains Harbin city and nine counties (or county-level cities). Research on the spatial structure of Harbin aims to reveal the spatial expansion and evolution of the construction land in central Harbin city [33,34,35] without considering rural areas and ecological problems. Additionally, some studies have placed Harbin city into a bigger scale and broader context, analyzing the urban agglomeration’s spatial structure evolution [36] and its patterns [37].
The main purpose of the study was to understand the patterns of the urban-rural spatial structure under ecological security via the constructed simulation model to propose optimization and implementation strategies in the next step. This study also aims to accelerate harmonious urban and rural development through overall consideration of urban-rural settlements, traffic systems, and land use zoning, thus further promoting the achievement of the final target of urban-rural integration. This is a significant and key aspect of Chinese national spatial planning.

2. Materials and Methods

2.1. Study Area and Data

Harbin city territory (i.e., the administrative area of Harbin), lies in the northeastern part of China (Figure 1) and is the capital city of Heilongjiang Province, which contains Harbin city and nine counties (or county-level cities). It has a total administrative area of 53,100 km2, a city area of 10,198 km2, and a central city area of 4187 km2. The total population in the administrative area was above 9.61 million by the end of 2015 (data from the ‘Harbin Statistics Yearbook 2016’). The urban-rural spatial structure for 2015 is shown in Figure 2. The Harbin urban-rural residential system comprises the central city (1), the sub-central city (2), the central town (14), the major town (58), the common town (111), the central village (1689), and the general village. The traffic network structure consists of a railway network and a highway network. The density of the traffic network above the county level is 0.39 km/km2.
The urban-rural development of Harbin city territory has a long history dating back to the first appearance of the five Ancient Cities during the Jin Dynasty (1115–1234 BC). Compared with other cities in Heilongjiang Province, the level of urban-rural integration in Harbin city is high [38], but compared to the whole of China, its level is relatively low [39]. Furthermore, the urban-rural integration between Harbin city and the other counties shows a development imbalance [40]. In order to enhance the overall urban-rural integration, this study, taking the Harbin city territory as the object of the study, simulates the urban-rural spatial structure of ecological priority to develop optimization strategies. A flowchart of the study structure is shown in Figure 3. For convenience, the abbreviation “Harbin” hereafter represents “Harbin city territory”.
Land use data used in this study were offered by Beijing Digital Technology Co. Ltd. Digital Elevation Model (DEM) and Moderate-resolution Imaging Spectroradiometer (MODIS) data were downloaded separately from the Geospatial Data Cloud [41] and National Aeronautics and Space Administration (NASA) [42]. Traffic, demographic, and Gross Domestic Product (GDP) data were taken from the ‘Harbin Yearbook’, ‘Harbin Statistics Yearbook’, and relevant government departments. All data were processed using ArcGIS.

2.2. Quantitative Evaluation Method

Green infrastructure refers to a network of natural and semi-natural green and blue areas that can maintain and improve ecosystem services, which comprises different kinds of green vegetation and water bodies [20,21,22]. Assessing green infrastructure with a quantitative evaluation method aims to understand the distribution of key ecological resources and sequence the ecological importance of these resources. The results will be used as restrictive conditions in the integrated simulation model. Therefore, an evaluation index system was built from resource protection and ecological services. According to the data availability, the scientific principle of index selection, and previous research [43,44,45,46,47], this study selected 12 indices to assess Harbin’s green infrastructure. Specific indicators can be seen in Table 1. The distance to each resource refers to the distance of the resource boundaries from a certain point, which can be measured by the buffer tool in ArcGIS. Generally, the closer the distance to each resource, the higher ecological sensitivity the point has during the process of land development. Thus, the ecological importance of this point is also higher. The distance of the buffer in this study was determined by multiple empirical studies [47,48,49]. The indices (e.g., the capability of support, the ability for temperature reduction, etc.) were measured by the product of the equivalents of the ecosystem service value supplied by the per unit area of the ecosystem [50] and the areas of various green infrastructures. The weights of the indicators were determined by expert marking and literature research on the network construction, planning, and evaluation of green infrastructure [45,46,47,51,52].
The expert marking method refers to consulting experts in the fields of landscape, ecology, and planning in the form of a questionnaire to obtain the importance value (Xij) of each index. Table 2 is an example of the expert marking table. A new score ( X i j ) is reassigned to the average score ( X i j ¯ ) of 10 experts according to Table 3. This new value is considered to be the final importance value for each indicator. Then, the weight ( w i j ) of each indicator can be calculated by Equation (1). Meanwhile, these 10 experts are invited to score the relative importance between the resource protection factors and the ecological service factors to obtain the combined weight ( w i ). Consequently, the final weight ( w i j ) of each indicator is the result of Equation (2). The overlay process is fulfilled in ArcGIS.
w i j = X i j / j = 1 n X i j
w i j = w i × w i j
where n is the number of indicators in the i category.

2.3. Prediction of the Urban-Rural Construction Land Scale

In general, the increase in construction land originates from the requirements for residential, public service, and industrial land due to population growth [53,54]. Therefore, this study built a linear regression equation between population, GDP, and construction land, but the results showed that this equation could not pass the DW (Durbin–Watson) test. According to the ‘Harbin Statistics Yearbook 2018’, the population of Harbin has experienced a downward trend since 2013. As a result, this study supposes that industrial production brought about the expansion of construction land, so the correlation between construction land and per capita GDP was analyzed and the correlation coefficient reached 0.9978. Meanwhile, the per capita GDP increased over time, and the two factors showed a linear correlation. Thus, using data from 2000 to 2016, first, the per capita GDP for the year 2035 was projected by linear regression Equation (3), and then the scale of the construction land was calculated in accordance with linear regression Equation (4). In addition, the linear regression equation of the construction land and time was also established [55] (Equation (5)). All equations passed the R test and the DW test.
y 1 = 3506.1 t 7008271.1
y 2 = 0.006489 y 1 + 1519.2
y 3 = 21.57806 t + 41590.6
where y 1 is GDP; y2 and y3 are the scale of urban-rural construction land; and t is time (year). Finally, the average of y2 and y3 was determined as the scale of urban-rural construction land in 2035.

2.4. Logistic-Celluar Aotomata Model

Geographical CA is a commonly used model in simulating urban land development [56,57] and comprises four basic parameters: cells, state, neighborhood, and transformation rules. The state of each cell at time t + 1 is decided by the state of the cell itself and its surroundings and the transformation rules at time t. This can be expressed by Equation (6).
S i j t + 1 = f ( S i j t , N i j t , R )
where S i j t + 1 and S i j t are the land use states of cell ij at time t + 1 and t, respectively; N i j t is the state set of cells around cell ij; and R represents the transformation rule. With the development of the CA model, the constrained urban CA model is able to mimic changes in the urban patterns of different developing targets or scenarios [25,27,58]. For instance, if ecological conditions are added to a CA model, then this CA model can simulate the urban development pattern of ecological security. This would provide a basis for making sustainable development policies. Equation (6) can be changed into
S i j t + 1 = f ( S i j t , N i j t , c o n ( S i j t ) , R )
where c o n ( S i j t ) expresses a constrained condition set that decides whether the state of cell ij at t moment can be transformed; its value is 0~1.
According to Equation (7), adopting an 8-neighborhood rule and transformation rule obtained by the logistic regression equation method [59], taking the results of green infrastructure assessment as an ecological constrained condition, this study simulates the evolution of urban-rural construction land. On the basis of the simulation results, the size of the urban and rural construction land was computed. In accordance with the computed size, urban and rural settlements in a specific region can be divided into the central city, sub-central city, central town, major town, common town, etc., which were combined into the hierarchical size structure [29,60]. Generally, the determination of a city’s categories is based on population size, which was decided by the size of the construction land in this study due to the limitations of data acquisition. Therefore, the “size” in the hierarchical size structure of this study was the size of the construction land.

2.5. Least Cost Path Model

The LCP model reflects the cost consumed or the work done when organisms move from the sources to the destinations. As it features a concise data structure, a fast algorithm, and intuitive results, LCP is widely used to study route simulations and planning problems [61,62,63,64,65]. This study aims to simulate the optimal travel path between urban-rural settlements on the basis of ecological protection. Consequently, an LCP model was selected. In order to apply the LCP model, three elements (sources, resistance factors, and resistance) need to be determined in advance. Sources are the origins of species dispersal and maintenance, which are usually composed of protected species or existing habitats. In this study, sources refer to urban-rural settlements that support the agglomeration and diffusion of the population and economic elements. The resistance factors vary based on the research objectives. Construction suitability, ecological condition, and land use types constitute the resistance factors of this study. The resistance values of the factors are different in terms of the subjects, which are usually decided upon by experts and empirical values [49,63], but these values only represent the relative resistance for organisms to migrate and disperse. Taking the setting of the resistance value for land use types as an example, in the field of ecology, the resistance of construction land is normally the highest, while the resistances of natural reserves and forests are relatively lower in terms of the migration of species, which is the opposite in this study because the resistance value of this LCP model refers to the cost of the movement of residents or material flows. Finally, the resistance factors and their resistance values were overlaid to obtain the resistance surface, which represents the difficulty of flowing.

2.6. Classification of Ecological Importance

The purpose of dividing land by its level of ecological importance is to provide instruction for making land protection and development measures through understanding the geographical distribution of land with its eco-value and, ultimately, to realize urban-rural sustainable development. In this study, the boundaries and ranks of ecological spaces can be determined according to the classification of ecological importance, which is the basis of dividing urban-rural functional spaces and determining controlling measures.
This study estimates and sequences the ecological importance of urban-rural land on the basis of green infrastructure assessment and adopts a natural break method to rank importance. Based on the obtained results of ecological importance zoning, coupled with the scope of the urban-rural construction land and farmland, the urban-rural land use zoning structure can be determined.

3. Results

3.1. Harbin Green Infrastructure Assessment

According to the definition of green infrastructure used in this study, Harbin’s green infrastructure consists of several land use types including farmland, forest, grassland, water bodies, and other land [66]. Applying the established evaluation index system to the green infrastructure assessment of Harbin in 2015, Figure 4 describes the sequencing results of the ecological importance of Harbin’s green infrastructure. In general, the green infrastructure with high ecological importance is distributed in the middle region of Harbin, while those with lower ecological importance are primarily concentrated in Harbin city. Natural reserves, great parks, etc. that need to be protected, gain the highest value, which means that strategies of priority protection should be applied to these areas. It should be noted that the purple parts in Figure 4 are urban-rural construction areas, which do not belong to the green infrastructure. Consequently, in the subsequent study on urban-rural spatial structure simulation, areas with higher ecological importance are regarded as limited development areas, while urban-rural construction areas are considered suitable development areas, where the value of importance is 0. This would benefit eco-system security and sustainable land development in Harbin.

3.2. Hierarchical Size Structure Extraction

This study adopted a land scale instead of a population scale to discuss the urban-rural hierarchical size structure of Harbin in 2035. Therefore, a logistic-CA model was used to simulate the expansion of urban-rural construction land. On the basis of this simulation, the urban-rural hierarchical size structure can be extracted.

3.2.1. Parameter Settings and Accuracy Test

Before applying a logistic-CA model, it is necessary to set parameters and examine the accuracy of the model in advance. The resolution of the raster data used in the model in this study was 90 m x 90 m. The parameter settings of the model are as follows:
  • Spatial variables of the transformation rule (R) based on logistic regression: distance to the highway, distance to the main road, distance to the urban center, distance to the town center, per capita GDP, population density, elevation, and slope;
  • Neighborhood rule ( N i j t ): a neighborhood space of 3 x 3 grid cells;
  • Constrained condition ( c o n ( S i j t ) ): the results of green infrastructure assessment.
Using the land use data of Harbin from 2005, 2013, and 2015, the transformation rule and the accuracy of the logistic-CA model was obtained and tested. In SPSS, taking two land use maps as the inputs of the dependent variable, and the distance to the highway (X1), distance to the main road (X2), distance to the urban center (X3), distance to the town center (X4), per capita GDP (X5), population density (X6), elevation (X7), and slope (X8) as the independent variables, the logistic regression Equation (8) was established:
R = 1 1 + e ^ [ ( 0.812 0.091 X 1 + 0.227 X 2 + 1.128 X 3 + 0.623 X 4 + 0.531 X 5 + 0.856 X 6 4.497 X 7 7.687 X 8 ) ]
The slope and elevation factors have the most obvious influence on the development of urban-rural construction land. Comparing the simulation results of 2015 (Figure 5a) with the real urban-rural construction land development in 2015 (Figure 5b), and adopting a pixel by pixel comparison method, the results indicate that the overall accuracy value of the model reached 96.88%. The Kappa simulation, which refers to the coefficient of agreement between the simulated land use transitions and the observed land use transitions [67], was 0.82. In addition, the Moran I index, which can reflect the extent of the spatial concentration or dispersion [68], was also calculated (Table 4). The results show that the simulated Moran I index value was lower than the real value of 2015 as well as the value of 2005. Thus, it can be seen that the land use pattern of Harbin tends to disperse evolution. On the whole, the simulated results from the model were close to the observed development of construction land. Consequently, this model can be used to predict land use changes in the future.

3.2.2. The Hierarchical Size Structure of Harbin in 2035

According to Equation (3), the per capita GDP of Harbin increased up to 126,655 Yuan in 2035. Thus, the scale of the construction land of Harbin for 2035, calculated by Equations (4) and (5), is 2341.06 km2 and 2320.75 km2, respectively. Therefore, the mean of 2330.91 km2 is considered as the urban-rural construction land scale for 2035.
Based on the predicted scale of 2035, the results of the green infrastructure assessment, and the tested logistic-CA model, the land use change of Harbin in 2035 was obtained (Figure 6). It was found that the Moran I index for 2035 declines, which means that the land use pattern in Harbin is more dispersed. Appropriate land use adjustment measures should be taken to enhance the intensive utilization of urban-rural land. Subsequently, the urban-rural hierarchical size structure for Harbin in 2035 was extracted based on the land use of 2035. As shown in Figure 7, Harbin’s urban-rural residential system comprises a central city (1), sub-central city (5), central town (12), major town (72), common town (94), central village (2301), and general village (28,665). Out of these categories, because Harbin city has a land area of more than 400 km2, it was classed as the central city; cities with land areas between 10 and 20 km2 including Wuchang, Shangzhi, Yanshou, Tonghe, and Binzhou would be classified as sub-central cities; the land areas of Yilan, Bayan, Mulan, Fangzheng, Dalianhe, Yimianpo, Weihe, Yabuli, Lalin, Binxi, Xinglong, and Shanhe towns range from 5 to 10 km2, which would make them central towns; towns with a scale of construction land from 1 to 5 km2 are major towns; other towns whose areas are smaller than 1 km2 are common towns; villages with a land area between 0.2 and 1 km2 are considered central villages; other villages are general villages. By comparing the hierarchical size structure for 2035 with that of 2015, the number of sub-central cities and major towns increases, but there is a decrease in the number of central towns, which occurs due to the upgrading of the Wuchang, Yanshou, and Binzhou towns, which indicates that the medium and small cities of Harbin will develop.

3.3. Traffic Network Structure Simulation

With the help of ArcGIS, based on the extracted hierarchical size structure, an LCP model was applied to simulate potential traffic routes. Harbin’s urban-rural traffic network system, based on the green infrastructure assessment used in this study, was obtained via the extraction of key traffic corridors and integration with current and forecasted traffic.

3.3.1. Simulation of Potential Traffic Routes

The simulation of potential traffic routes can be divided into three steps: source identification, building a resistant surface, and traffic network construction. In this study, urban-rural settlements were regarded as sources. The building of a resistant surface was implemented by using an overlay tool in ArcGIS, and the resistance factors and their resistance can be seen in Table 5. To set the range, the range of the green infrastructure assessment results was divided according to a manual classification method in ArcGIS, and the range of the slope was determined based on the suitability of the engineering construction [69]. The setting of resistance references was done according to previous studies [47,49,63]. The resistance value is a relative value, which only expresses the relative resistance of the landscape by a numerical value. Finally, with the aid of cost distance and cost path tools, the potential urban-rural traffic network (Figure 8) of Harbin for 2035, based on green infrastructure assessment, was mapped.

3.3.2. Extraction of Key Traffic Corridors

In accordance with the grade of the urban-rural settlements connected by traffic routes, key traffic corridors were extracted. Settlements with a higher grade in their hierarchical size structure are more important, and the traffic routes that link those settlements are more crucial. As a result, this study considered the potential routes that connected the central city, sub-central cities, and central towns as key traffic paths, and those that link major towns are deemed vice-key traffic paths (Figure 9).

3.3.3. Construction of Traffic Network

According to the importance of the highway, the grading of roads can be divided into five grades: motorway, national road, provincial road, county road, and rural road. On the basis of the ‘National Highway Network Planning’, and the ‘Heilongjiang Provincial Highway Network Planning’, ‘Comprehensive Transportation Planning of Harbin during the 13th Five-Year Plan’ proposed that the layout of the motorway system is “two rings, nine axes, and two links” (Figure 10a), the layout of the common trunk highway system is “one ring, thirteen axes, three vertical lines, and eleven links” (Figure 10b). In accordance with this transportation planning, coupled with Figure 8 and Figure 9, and the existing roads, Harbin’s urban-rural traffic network structure in 2035, based on a green infrastructure assessment, was constructed by adjusting the grade of some simulated roads. As shown in Figure 11, this structure consists of a railway network and a highway network. In this structure, the highway network contains a motorway system, a common trunk highway system composed of national roads, provincial roads, and county roads, and a rural highway system. The density of the urban-rural traffic network excluding the rural highway system will have increased by 0.08 km/km2 since 2015, when the density was 0.39 km/km2. This increase shows that traffic systems above the county level have improved.

3.4. Land Use Zoning Structure Determination

In urban-rural spatial planning, the division of urban-rural functional spaces on the basis of ecological importance zoning is an essential way to put the idea of ecological priority into practice. This study divides the urban-rural spaces of Harbin into urban-rural construction spaces, ecological spaces, and agricultural spaces based on the ecological importance zoning obtained via a classification of ecological importance. In this way, the land use zoning structure is determined.

3.4.1. Ecological Importance Zoning

According to the sequencing of ecological importance from the green infrastructure assessment, and taking a natural break classification approach [70,71] (natural break values in this approach is obtained by the Jenk optimization method of statistics, which can minimize the sum of internal variances at all levels [72]) in GIS, Harbin (Figure 12) can be divided into extremely important (5.335–7.28), very important (3.758–5.335), relatively important (2.826–3.758), generally important (1.135–2.826), and not important (0–1.135) zones. Extremely important and very important zones are regarded as excluded from development areas, where protection and restoration behavior are conducted. Relatively important and generally important zones are taken as restricted areas, in which measures of protection priority will be taken. Not important zones contain all the construction land that can be developed and can be classified as areas suitable for construction.

3.4.2. Delimitation of Urban-Rural Functional Spaces

In accordance with the theory of urban-rural integration and the definition of territorial functional space, this study classifies the functional spaces of Harbin as urban-rural construction spaces, eco-spaces, and agricultural spaces (Figure 13). Eco-spaces are divided into control regions of ecological red lines and other eco-spaces. The boundary lines of excluded from development areas with high eco-values are ecological red lines. Regions within the scope of the ecological red lines are classified as control regions that should be strictly protected to keep the urban-rural eco-system safe. The scopes of urban-rural construction land, as determined by the land use simulation for 2035, are identified as the urban-rural construction spaces of 2035. Comparing the current ecological importance zoning picture with the predicted land use change picture of 2035, agricultural spaces and other eco-spaces located in restricted areas were defined. In this study, agricultural spaces refer to agricultural production spaces, but do not include rural living spaces. Other eco-spaces mainly comprise grasslands and other types of land. The percentages of each kind of space in the total land area are as follows: eco-space, 51.4% (which controls regions of ecological red lines comprising 46.7%, while other eco-spaces comprise 4.7%); urban-rural space, 4.4%; and agricultural space, 44.0%. Construction spaces and agricultural spaces are primarily distributed in the smoother western areas, while eco-spaces are mostly concentrated in the middle area, where the terrain is higher, and the surface is covered by forests and rivers.

3.5. Integration of Urban-Rural Spatial Structures

The obtained urban-rural hierarchical size structure, traffic network structure, and land use zoning structure were integrated to obtain the overall urban-rural spatial structure of Harbin in 2035 based on the green infrastructure assessment. As shown in Figure 14, generally, the structure features a single center (Harbin city), a core distribution, a point-axis, and main cities and towns distributed in series along a trunk network that is made up of highways and railways. Measuring the extent of the single center by the urban primacy index showed that the centrality of Harbin city is predicted to keep increasing from 25.74 (2015) to 28 (2035). The urban-rural hierarchical size structure can be described as a seven-level settlement system composed of a central city, a sub-central city, a central town, a major town, a common town, a central village, and a general village. The urban-rural traffic network structure comprises railways and highways of different grades and forms a pattern of “radiation + grid”, with Harbin city as the center of the structure. The urban-rural land use zoning structure was classified into eco-spaces, agricultural spaces, and construction spaces.

4. Discussion

This study simulated and predicted the overall arrangement of the urban-rural spatial structure under the basis of ecological constraints in order to maintain the lowest effect on ecosystem security when urban-rural development is in progress. Controlling urban growth is a common approach for promoting urban-rural sustainable development [73,74]. Green infrastructures in this study were assessed in terms of resource protection and ecological services, instead of only from the perspective of ecological services [51,75]. The results showed that key areas, like natural reserves, gained a higher value of importance. This would be especially helpful for the establishment of subsequent development strategies.
This study considered settlements, traffic systems, and land use zoning as the three determinants of urban-rural spatial structures. In the process of hierarchical size structure extraction, although the Moran I index, which can reflect the extent of spatial concentration or dispersion, is lower for 2035 than for 2015, medium and small cities of Harbin will be further developed by 2035. These results are in accord with the national strategy in China. The rapid development of medium and small cities could have positive effects on the development of surrounding rural areas [76], which is of great significance for realizing balanced regional development and urban-rural integration, especially for developing countries. Meanwhile, the urbanization progress of the central city (Harbin city) is also continuing, whose speed is faster than that of medium and small cities. Consequently, the urban primacy index increased, which indicates that the centrality of the central city will still increase. The traffic network supports urban-rural material flow and is the necessary basis for urban-rural integration. Many studies have predicted traffic flow and explained spatial-temporal correlations [77,78], or forecasted traffic volume through a significant number of surveys [79], but their aims were not to explore the spatial layout of the traffic network. According to the simulation result of this study, the traffic network structure of Harbin in 2035 tends to be more complex. The grade and quality of the roads will also be enhanced. Drawing on the experience of urban-rural integration in developed countries [80,81], the current service capacity including the grade and quality of highways and railways will be improved to connect all settlements through the traffic network. It should be noted that because the Harbin railway development will mainly be fulfilled by upgrading existing railway facilities, the layout of the railway in this study was maintained the same as the current situation. This could be further studied when more data are available. The land in this study was zoned in terms of living, productive, and ecological functions responding to national spatial planning in China. Finally, the overlaid result showed that the urban-rural spatial structure still featured a significant single center. This is currently a common phenomenon for the development of regions in China. A review of urbanization in developed countries has found that the tendency to develop a single center could be an inevitable process preceding urban-rural integration [82].
In accordance with the equivalents of the ecosystem service values supplied by the per unit area of the ecosystem [50], water bodies can offer a superior water resources supply service, hydrological adjustment service, biodiversity support service, and cultural service; forests can provide a better gas regulation service, climate regulation service, and biodiversity support service; while agricultural land can offer a better food supply service. The application of the integrated model in this study can lead to policy recommendations that may support more sustainable development in Harbin. For example, through the policy of prohibiting development activities in water bodies and forest areas and the restricting misuse of agricultural land, multiple ecological services including water resource supply services, hydrological adjustment services, biodiversity support services, and so on, can be preserved to support sustainable development in Harbin.
Previous research on spatial structure only simulated land use structures from the perspective of the surface [83,84,85], while the integrated model developed in this study comprehensively predicts the urban-rural spatial structure in the future from a point-line-surface perspective. In addition, the model of this study can determine the sequence of land development in advance on the basis of the results of a green infrastructure assessment. Then, the urban-rural spatial structure can be simulated. Therefore, the integrated model in this study can be considered as a model of sustainable development. The simulation results based on this model could effectively protect important ecological resources. To further explain the strength of the integrated model in this study, a comparison between the ecological constraint scenario that we developed and the scenario of urbanization as usual should be undertaken. This work is being done at present. Furthermore, to realize the results of the simulation in this study, some optimization measurements must be undertaken, for example, how to promote the fast development of medium and small settlements in Harbin, how to improve the service capacity of the traffic network in Harbin, and how to control and manage developing activities in ecological, agricultural, and construction spaces for Harbin. All of these questions may be addressed in future research in this area.
This study focused on problems at the macro scale and considers urban areas and rural areas in a comprehensive way. To some extent, by giving priority to the protection of areas with higher ecological importance to restrict the location of construction land development, the problem of mixed urban-rural land uses (which entails more detailed and complicated issues such as fragmentation, leap-frog development, and low-density urbanization) can be modified. However, this study is not concerned about problems of large-scale development within the city. For these problems, some restoration and compulsory measures will be necessary. This will also be studied more deeply in the future.

5. Conclusions

From an ecological perspective, this study proposed an integrated simulation model to simulate the urban-rural spatial structure of Harbin for 2035, the result of which is expected to further promote the sustainable development of urban-rural integration through the limitation of green infrastructure on land use transitions induced by urbanization. This study, using a quantitative evaluation method, the logistic-CA model, the LCP model, and a classification of ecological importance, simulated the hierarchical size structure, traffic network structure, and land use zoning structure, and integrated them into an urban-rural spatial structure for Harbin. The conclusions based on the simulated results from this integrated model are as follows:
  • The integrated simulation model based on green infrastructure assessment in this study can be used to simulate urban-rural spatial structures.
  • According to the simulated results of this study, although land use patterns in Harbin in 2035 will be more dispersed, the medium and small cities of Harbin will become developed.
  • The simulated results indicate that in 2035, the traffic network of Harbin will still consist of a railway network and highway network, and the traffic systems above the county level will be improved.
  • For Harbin, based on the simulated results, most of the land in 2035 will be covered by ecological land, followed by agricultural land, while construction spaces will be primarily distributed in the smoother western area.
  • Through the integrated simulation model, the urban-rural spatial structure of Harbin shows an increasing development tendency for a single center. The main cities and towns are distributed in series along with a trunk network.

Author Contributions

Conceptualization, R.G. and Y.B.; Methodology, Y.B.; Validation, Y.B.; Formal analysis, Y.B.; Investigation, Y.B.; Data curation, R.G.; Writing—original draft preparation, Y.B.; Writing—review and editing, Y.B.; Visualization, Y.B.; Supervision, R.G.; Project administration, R.G.; Funding acquisition, R.G.

Funding

This research was funded by the Later Funded Program of Chinese Ministry of Education, grant number 17JHQ036.

Acknowledgments

We thank Tong Wu for her technical support in the paper’s revisions.

Conflicts of Interest

The authors declare no conflicts of interest. The funders had roles in the design of the study; in the collection of data; in the decision to publish the results.

References

  1. Luo, L.; Wang, Y. The Operation “Mode” of the Practice of Urban-rural Integration: From the Prospective of Britain, France and the United States. Reform 2014, 3, 116–122. [Google Scholar]
  2. Xue, Q.; Sun, H. Examples of successful experience in development of urban and rural integration abroad. Agric. Econ. 2014, 1, 61–63. [Google Scholar]
  3. Wang, S.; Tan, S.; Yang, S.; Lin, Q.; Zhang, L. Urban-biased land development policy and the urban-rural income gap: Evidence from Hubei Province, China. Land Use Policy 2019, 87, 104066. [Google Scholar] [CrossRef]
  4. Wang, Y.; Li, Y.; Huang, Y.; Yi, C.; Ren, J. Housing wealth inequality in China: An urban-rural comparison. Cities 2020, 96, 102428. [Google Scholar] [CrossRef]
  5. Mumford, L. The City in History. Its Origins, Its Transformation, and Its Prospects, 1st ed.; Houghton Mifflin Harcourt: Boston, MA, USA, 1961. [Google Scholar]
  6. Wan, G.; Wang, C. Unprecedented Urbanisation in Asia and Its Impacts on the Environment. Aust. Econ. Rev. 2014, 47, 378–385. [Google Scholar] [CrossRef]
  7. Dame, J.; Schmidt, S.; Müller, J.; Nüsser, M. Urbanisation and socio-ecological challenges in high mountain towns: Insights from Leh (Ladakh), India. Landsc. Urban Plan. 2019, 189, 189–199. [Google Scholar] [CrossRef]
  8. Komugabe-Dixson, A.F.; de Ville, N.S.E.; Trundle, A.; McEvoy, D. Environmental change, urbanisation, and socio-ecological resilience in the Pacific: Community narratives from Port Vila, Vanuatu. Ecosyst. Serv. 2019, 39, 100973. [Google Scholar] [CrossRef]
  9. Fang, C.; Cui, X.; Li, G.; Bao, C.; Wang, Z.; Ma, H.; Sun, S.; Liu, H.; Luo, K.; Ren, Y. Modeling regional sustainable development scenarios using the Urbanization and Eco-environment Coupler: Case study of Beijing-Tianjin-Hebei urban agglomeration, China. Sci. Total Environ. 2019, 689, 820–830. [Google Scholar] [CrossRef]
  10. Chen, D.; Lu, X.; Liu, X.; Wang, X. Measurement of the eco-environmental effects of urban sprawl: Theoretical mechanism and spatiotemporal differentiation. Ecol. Indic. 2019, 105, 6–15. [Google Scholar] [CrossRef]
  11. Zeng, Z.; Xia, H.; Chen, H. Research of Chinese Ancient Urban Morphologies Based on Climate Adaptability. In LTLGB 2012; Springer: Berlin/Heidelberg, Germany, 2013; pp. 781–785. [Google Scholar]
  12. Xu, P. Feng-shui models structured traditional Beijing courtyard houses. J. Archit. Plann. Res. 1998, 15, 271–282. [Google Scholar]
  13. Dalley, S. Ancient Mesopotamian Gardens and the Identification of the Hanging Gardens of Babylon Resolved. Gard. Hist. 1993, 21, 1–13. [Google Scholar] [CrossRef]
  14. Bassin, J. The English Landscape Garden in the Eighteenth Century: The Cultural Importance of an English Institution. Albion 1979, 11, 15–32. [Google Scholar] [CrossRef]
  15. McHarg, I.L. Design with Nature; John Wiley & Sons: New York, NY, USA, 1992. [Google Scholar]
  16. Yang, Z. Eco-Cities: A Planning Guide, 1st ed.; CRC Press: Boca Raton, FL, USA, 2012; ISBN 9781439883228. [Google Scholar]
  17. Zhang, X.; Hes, D.; Wu, Y.; Hafkamp, W.; Lu, W.; Bayulken, B.; Schnitzer, H.; Li, F. Catalyzing sustainable urban transformations towards smarter, healthier cities through urban ecological infrastructure, regenerative development, eco towns and regional prosperity. J. Clean. Prod. 2016, 122, 2–4. [Google Scholar] [CrossRef] [Green Version]
  18. Lafortezza, R.; Davies, C.; Sanesi, G.; Konijnendijk, C.C. Green Infrastructure as a tool to support spatial planning in European urban regions. iForest-Biogeosciences For. 2013, 6, 102–108. [Google Scholar] [CrossRef] [Green Version]
  19. European Commission. The Multifunctionality of Green Infrastructure; European Commission: Brussels, Belgium, 2012. [Google Scholar]
  20. The North West Green Infrastructure Think Tank. North West Green Infrastructure Guide. The North West Green Infrastructure Think Tank: UK, 12 May 2008. [Google Scholar]
  21. Naumann, S.; Davis, M.; Kaphengst, T.; Pieterse, M.; Rayment, M. Design, Implementation and Cost Elements of Green Infrastructure Projects. Final report to the European Commission; DG Environment: Brussels, Belgium, 2011. [Google Scholar]
  22. Ely, M.; Pitman, S. Green Infrastructure: Life Support for Human Habitats; Botanic Gardens of Adelaide, Department of Environment, Water and Natural Resources: Adelade, Australia, 2014.
  23. Pitman, S.D.; Daniels, C.B.; Ely, M.E. Green infrastructure as life support: Urban nature and climate change. Trans. R. Soc. S. Aust. 2015, 139, 97–112. [Google Scholar] [CrossRef]
  24. Zhang, J. Green Infrastructure: A Systematic Solution to Urban Space and Environmental Issues. Mod. Urban Res. 2009, 24, 81–86. [Google Scholar]
  25. Zhu, J.; Yu, Y.; Zhou, S.; Wang, X.; Lv, L. Simulating Sustainable Urban Development by Incorporating Social-ecological Risks into a Constrained CA Model. Chin. Geogr. Sci. 2018, 28, 600–611. [Google Scholar] [CrossRef] [Green Version]
  26. Milad Aburas, M.; Ming Ho, Y.; Firuz Ramli, M.; Hanan Ash, Z. Improving the capability of an integrated CA-Markov model to simulate spatio-temporal urban growth trends using an Analytical Hierarchy Process and Frequency Ratio. Int. J. Appl. Earth Obs. Geoinf. 2017, 59, 65–78. [Google Scholar] [CrossRef]
  27. de Noronha Vaz, E.; Nijkamp, P.; Painho, M.; Caetano, M. A multi-scenario forecast of urban change: A study on urban growth in the Algarve. Landsc. Urban Plan. 2012, 104, 201–211. [Google Scholar] [CrossRef] [Green Version]
  28. Zhou, B.; Dai, T.; Liang, J.; Zhang, H. Simulation of urban spatial structure in Beijing based on Lowry model. Acta Georraphica Sin. 2013, 68, 491–505. [Google Scholar]
  29. Zou, H. Study on Rapid Urbanization Oriented Urban and Rural Spatial Structure’s Transformation and Optimization mode of Weinan. Master’s Thesis, Xi’an University of Architecture and Technology, Xi’ an, China, 2014. [Google Scholar]
  30. Peng, W. Study on Reconstructing Urban-Rural Spatial Structure of Panjin City. Master’s Thesis, Lanzhou University, Lanzhou, China, 2013. [Google Scholar]
  31. Huang, M.; Guo, R. Scenario simulation of Harbin urban expansion and its spatial pattern evaluation with GIS/RS. In Proceedings of the Earth Resources and Environmental Remote Sensing/GIS Applications IX, Berlin, Germany, 10–13 September 2018; Michel, U., Schulz, K., Eds.; SPIE: Bellingham, WA, USA, 2018; Volume 10790. 107901Y. [Google Scholar]
  32. Gong, W.; Yuan, L.; Fan, W.; Stott, P. Analysis and simulation of land use spatial pattern in Harbin prefecture based on trajectories and cellular automata—Markov modelling. Int. J. Appl. Earth Obs. Geoinf. 2015, 34, 207–216. [Google Scholar] [CrossRef]
  33. Zhang, X. Study on Urban Spatial Evolution of Harbin. Master’s Thesis, Northest Normal University, Changchun, China, 2002. [Google Scholar]
  34. Liu, J. Research on urban spatial structure of Harbin. J. Harbin Comm. Sch. CCP 2012, 68–74. [Google Scholar]
  35. Liu, Y. Effects of Social Changes in Modern Harbin on Evolvement of Urban Spatial Structure. Master’s Thesis, Northeast Normal University, Changchun, China, 2007. [Google Scholar]
  36. Ma, C.; Zhao, T. The Correlation between Industrial Organization and Sustainable Growth of City Cluster Spatial: A Case Study of Harbin-Changchun City Cluster. Energy Procedia 2019, 157, 773–781. [Google Scholar] [CrossRef]
  37. Wang, S.; Song, Y.; Feng, Z.; Jiang, L. Pattern and Progress of Large Urban Agglomerations and Urban Flows Intensity in Northeast China. Sci. Geogr. Sin. 2011, 31, 287–294. [Google Scholar]
  38. Liu, J.; Chen, D. Evaluation and Analysis on the Process of Urban-rural Integration in Heilongjiang Province. Chin. Agric. Sci. Bull. 2014, 30, 107–113. [Google Scholar]
  39. Xiao, X.; Luo, N. Measurement and Cluster Analysis of Regional Urban-Rural Integration Level: Based on the Interprovincial Panel Data of China from 2003 to 2012. Seeker 2015, 2, 48–52. [Google Scholar]
  40. Ma, H.; Zhang, Y.; Li, X. A Quantitative Analysis and Research for the Integration of Town and Country:Evidence from Harbin Region. Mod. Urban Res. 2009, 6, 90–95. [Google Scholar]
  41. Geospatial Data Cloud. Available online: http://www.gscloud.cn/ (accessed on 14 December 2019).
  42. Earthdata Search. Available online: https://search.earthdata.nasa.gov/search (accessed on 14 December 2019).
  43. Schäffler, A.; Swilling, M. Valuing green infrastructure in an urban environment under pressure—The Johannesburg case. Ecol. Econ. 2013, 86, 246–257. [Google Scholar] [CrossRef]
  44. Sandstrom, U.G. Green Infrastructure Planning in Urban Sweden. Plan. Pract. Res. 2002, 17, 373–385. [Google Scholar] [CrossRef]
  45. Liu, M.; Fan, J.; Yu, Z. Multifunctional Green Infrastructure Planning: A Case Study of Haidian District, Beijing. Chin. Landsc. Archit. 2013, 29, 61–66. [Google Scholar]
  46. Pakzad, P.; Osmond, P. Developing a Sustainability Indicator Set for Measuring Green Infrastructure Performance. Procedia-Soc. Behav. Sci. 2016, 216, 68–79. [Google Scholar] [CrossRef] [Green Version]
  47. An, C.; Shen, Q. Method of Constructing Network of Green Infrastructure based on Ecological Performance of Space Utilization. Landsc. Archit. 2013, 2, 22–31. [Google Scholar]
  48. Zhang, X.; Wang, Z.; Lin, J. GIS Based Measurement and Regulatory Zoning of Urban Ecological Vulnerability. Sustainability 2015, 7, 9924–9942. [Google Scholar] [CrossRef] [Green Version]
  49. Yin, H.; Kong, F. Lab Manual for Spatial Analysis in Urban and Rural Planning, 1st ed.; Southeast University Press: Nanjing, China, 2014. [Google Scholar]
  50. Xie, G.; Zhang, C.; Zhang, C.; Xiao, Y.; Lu, C. The value of ecosystem services in China. Resour. Sci. 2015, 37, 1740–1746. [Google Scholar]
  51. Weber, T.; Sloan, A.; Wolf, J. Maryland’s Green Infrastructure Assessment: Development of a comprehensive approach to land conservation. Landsc. Urban Plan. 2006, 77, 94–110. [Google Scholar] [CrossRef]
  52. Zong, M. A Study of Urban Green Infrastructure Network Construction and Planning Pattern. Shanghai Urban Plan. Rev. 2015, 3, 104–109. [Google Scholar]
  53. Han, J.; Hayashi, Y.; Cao, X.; Imura, H. Application of an integrated system dynamics and cellular automata model for urban growth assessment: A case study of Shanghai, China. Landsc. Urban Plan. 2009, 91, 133–141. [Google Scholar] [CrossRef]
  54. He, C.; Okada, N.; Zhang, Q.; Shi, P.; Zhang, J. Modeling urban expansion scenarios by coupling cellular automata model and system dynamic model in Beijing, China. Appl. Geogr. 2006, 26, 323–345. [Google Scholar] [CrossRef]
  55. Quan, Q.; Tian, G.; Sha, M. Dynamic simulation of Shanghai urban expansion based on multi-agent system and cellular automata models. Acta Ecol. Sin. 2011, 31, 2875–2887. [Google Scholar]
  56. Shen, T.-Y.; Wang, W.-D.; Hou, M.; Guo, Z.-C.; Xue, L.; Yang, K.-Z. Study on Spatio-Temporal System Dynamic Models of Urban Growth. Syst. Eng.-Theory Pract. 2007, 27, 10–17. [Google Scholar] [CrossRef]
  57. Shafizadeh Moghadam, H.; Helbich, M. Spatiotemporal urbanization processes in the megacity of Mumbai, India: A Markov chains-cellular automata urban growth model. Appl. Geogr. 2013, 40, 140–149. [Google Scholar] [CrossRef]
  58. Guan, Q.; Wang, L.; Clarke, K.C. An artificial-neural-network-based, constrained CA model for simulating urban growth. Cartogr. Geogr. Inf. Sci. 2005, 32, 369–380. [Google Scholar] [CrossRef] [Green Version]
  59. Shafizadeh-Moghadam, H. Improving spatial accuracy of urban growth simulation models using ensemble forecasting approaches. Comput. Environ. Urban Syst. 2019, 76, 91–100. [Google Scholar] [CrossRef]
  60. Wu, H. Study on Urban and Rural Spatial Structure’s Optimization of Pi County Based on the Itegration of Urban and Rural areas. Master’s Thesis, Xi’an University of Architecture and Technology, Xi’ an, China, 2011. [Google Scholar]
  61. Collischonn, W.; Pilar, J.V. A direction dependent least cost path algorithm for roads and canals. Int. J. Geogr. Inf. Sci. 2000, 14, 397–406. [Google Scholar] [CrossRef]
  62. Lee, J.; Munro-stasiuk, M.J. Extensions to least-cost path algorithms for roadway planning. Geogr. Inf. Sci. 2003, 17, 361–376. [Google Scholar]
  63. Effat, H.A.; Hassan, O.A. Designing and evaluation of three alternatives highway routes using the Analytical Hierarchy Process and the least-cost path analysis, application in Sinai Peninsula, Egypt. Egypt. J. Remote Sens. Sp. Sci. 2013, 16, 141–151. [Google Scholar] [CrossRef] [Green Version]
  64. Driezen, K.; Adriaensen, F.; Rondinini, C.; Doncaster, C.P.; Matthysen, E. Evaluating least-cost model predictions with empirical dispersal data: A case-study using radiotracking data of hedgehogs (Erinaceus europaeus). Ecol. Model. 2007, 209, 314–322. [Google Scholar] [CrossRef]
  65. LaRue, M.A.; Nielsen, C.K. Modelling potential dispersal corridors for cougars in midwestern North America using least-cost path methods. Ecol. Model. 2008, 212, 372–381. [Google Scholar] [CrossRef]
  66. Bai, Y.; Guo, R.; Xing, Y. Relationship between urban heat island and green infrastructure fraction in Harbin. In Proceedings of the Remote Sensing Technologies and Applications in Urban Environments III, Berlin, Germany, 10–13 September 2018; Chrysoulakis, N., Erbertseder, T., Zhang, Y., Eds.; SPIE: Bellingham, WA, USA, 2018; Volume 10793. 107930D. [Google Scholar]
  67. Van Vliet, J.; Bregt, A.K.; Hagen-Zanker, A. Revisiting Kappa to account for change in the accuracy assessment of land-use change models. Ecol. Model. 2011, 222, 1367–1375. [Google Scholar] [CrossRef]
  68. Wu, F. International Journal of Geographical Information Science Calibration of stochastic cellular automata: The application to rural-urban land conversions Calibration of stochastic cellular automata: The application to rural-urban land conversions. Int. J. Geogr. Inf. Sci. 2002, 16, 795–818. [Google Scholar] [CrossRef]
  69. Chen, Y.; Du, P.; Zheng, X.; Lin, J. Evaluation on ecological applicability of land construction in Nanning city based on GIS. J. Tsighua Univ. Technol. 2006, 46, 801–804. [Google Scholar]
  70. Febrianto, H.; Fariza, A.; Hasim, J.A.N. Urban flood risk mapping using analytic hierarchy process and natural break classification (Case study:Surabaya, East Java, Indonesia). In Proceedings of the International Conference on Knowledge Creation & Intelligent Computing, Manado, Indonesia, 15–17 November 2016; IEEE: Manado, Indonesia, 2017; pp. 148–154. [Google Scholar]
  71. Zhu, L.; Wang, H.; Tang, L. Importance evaluation and spatial distribution analysis of ecosystem services in Min triangle area. Acta Ecol. Sin. 2018, 38, 7254–7268. [Google Scholar]
  72. Li, Y.; Liu, C.; Min, J.; Wang, C.; Zhang, H.; Wang, Y. RS/GIS-based integrated evaluation of the ecosystem services of the Three Gorges Reservoir area. Acta Ecol. Sin. 2013, 33, 168–178. [Google Scholar]
  73. Siedentop, S.; Fina, S.; Krehl, A. Greenbelts in Germany’s regional plans—An effective growth management policy? Landsc. Urban Plan. 2016, 145, 71–82. [Google Scholar] [CrossRef]
  74. Yu, D.; Yanxu, L.; Bojie, F. Urban growth simulation guided by ecological constraints in Beijing city: Methods and implications for spatial planning. J. Environ. Manag. 2019, 243, 402–410. [Google Scholar] [CrossRef]
  75. Flynn, K.M.; Traver, R.G. Green infrastructure life cycle assessment: A bio-infiltration case study. Ecol. Eng. 2013, 55, 9–22. [Google Scholar] [CrossRef]
  76. Berdegué, J.A.; Soloaga, I. Small and medium cities and development of Mexican rural areas. World Dev. 2018, 107, 277–288. [Google Scholar] [CrossRef]
  77. Do, L.N.N.; Vu, H.L.; Vo, B.Q.; Liu, Z.; Phung, D. An effective spatial-temporal attention based neural network for traffic flow prediction. Transp. Res. Part C Emerg. Technol. 2019, 108, 12–28. [Google Scholar] [CrossRef]
  78. Polson, N.G.; Sokolov, V.O. Deep learning for short-term traffic flow prediction. Transp. Res. Part C Emerg. Technol. 2017, 79, 1–17. [Google Scholar] [CrossRef] [Green Version]
  79. Xu, H.-J.; Li, W.-Y.; Wang, T.; Yang, A.-L. Research on dynamic prediction method for traffic demand based on trip generation analysis. Adv. Mech. Eng. 2019, 11, 1687814019854561. [Google Scholar] [CrossRef] [Green Version]
  80. Kwŏn, A.; Kim, K. (Eds.) Urban Management in Seoul: Policy Issues Responses, 1st ed.; Seoul Development Institute: Seoul, Korea, 2001. [Google Scholar]
  81. Henderson, D.A. Urbanization of Rural America; Nova Science Publishers: New York, NY, USA, 1997. [Google Scholar]
  82. Shi, Y. From A Single Center City to A Multi-centers City—Spatial Organization Pattern of Metropolitan Development in China. Urban Plan. Forum 1999, 36–39, 26. [Google Scholar]
  83. Zhu, Z.; He, Q. Dynamic Simulation of Evolution of Urban Spatial Structure of Changsha City. Econ. Geogr. 2016, 36, 50–58. [Google Scholar]
  84. Zhang, L. Simulation of Urban Spatial Structure Dynamic Process in Nanjing Based on Cellular Automata. Geomat. Spat. Inf. Technol. 2010, 33, 169–172. [Google Scholar]
  85. Guan, Z. Study on Urban Spatial Structure Changes of Jiaozuo Based on SLEUTH Model. Master’s Thesis, Suzhou University of Science and Technology, Suzhou, China, 2010. [Google Scholar]
Figure 1. The location and extent of study area.
Figure 1. The location and extent of study area.
Land 08 00196 g001
Figure 2. The urban-rural spatial structure in 2015. 哈尔滨市行政区划图: Administrative district map of Harbin.
Figure 2. The urban-rural spatial structure in 2015. 哈尔滨市行政区划图: Administrative district map of Harbin.
Land 08 00196 g002
Figure 3. A flowchart of the study structure.
Figure 3. A flowchart of the study structure.
Land 08 00196 g003
Figure 4. The sequencing of ecological importance for Harbin’s green infrastructure in 2015.
Figure 4. The sequencing of ecological importance for Harbin’s green infrastructure in 2015.
Land 08 00196 g004
Figure 5. A comparison of the urban-rural construction land development for Harbin in 2015. (a) Simulated result. (b) Observed development.
Figure 5. A comparison of the urban-rural construction land development for Harbin in 2015. (a) Simulated result. (b) Observed development.
Land 08 00196 g005
Figure 6. Land use change in 2035.
Figure 6. Land use change in 2035.
Land 08 00196 g006
Figure 7. Hierarchical size structure in 2035.
Figure 7. Hierarchical size structure in 2035.
Land 08 00196 g007
Figure 8. Potential traffic network (Harbin 2035).
Figure 8. Potential traffic network (Harbin 2035).
Land 08 00196 g008
Figure 9. Key traffic corridors (Harbin 2035).
Figure 9. Key traffic corridors (Harbin 2035).
Land 08 00196 g009
Figure 10. Comprehensive transportation planning of Harbin during the 13th Five-Year Plan. (a) The planned motorway system. (b) The planned common trunk highway system.
Figure 10. Comprehensive transportation planning of Harbin during the 13th Five-Year Plan. (a) The planned motorway system. (b) The planned common trunk highway system.
Land 08 00196 g010
Figure 11. The urban-rural traffic network structure of Harbin in 2035.
Figure 11. The urban-rural traffic network structure of Harbin in 2035.
Land 08 00196 g011
Figure 12. Ecological importance zoning.
Figure 12. Ecological importance zoning.
Land 08 00196 g012
Figure 13. Land use zoning structure.
Figure 13. Land use zoning structure.
Land 08 00196 g013
Figure 14. Urban-rural spatial structure of Harbin in 2035 based on the green infrastructure assessment.
Figure 14. Urban-rural spatial structure of Harbin in 2035 based on the green infrastructure assessment.
Land 08 00196 g014
Table 1. The evaluation index system of green infrastructure 1.
Table 1. The evaluation index system of green infrastructure 1.
CategoriesFactorsIndices (Xij)Weights
Resources protection (X1)Drinking water sourcesThe distance to water sources (X11)0.10
Great parks (>= 100 km2)The distance to parks (X12)0.07
Nature reservesThe distance to nature reserves (X13)0.09
Basic farmland protection zonesThe distance to basic farmland protection zones (X14)0.09
Vegetation coverageNDVI (X15)0.05
Ecological services (X2)Biodiversity supportThe capability of support (X21)0.15
Climate modificationThe ability of temperature reduction (X22)0.09
Carbon regulationNet capacity of carbon absorption (X23)0.09
Hydrological adjustmentAdjusting capacity of hydrology (X24)0.09
Food supplyThe ability of food production (X25)0.06
Water resources supplyWater supply ability (X26)0.06
Cultural entertainmentThe capacity of entertainment experience (X27)0.06
1 The units of the indicators associated with distance are in kilometers, and the units of the indicators related to capacity or ability are ten thousand yuan.
Table 2. An example of the expert marking table.
Table 2. An example of the expert marking table.
CategoriesIndicesEcological Importance Score (Xij)
97531
Resource protectionThe distance to parks (X1)Extremely importantvery importantrelatively importantgenerally importantnot important
Ecological servicesThe capability of support (X6)Extremely importantvery importantrelatively importantgenerally importantnot important
Table 3. The comparison table of the newly reassigned score and the average score of 10 experts.
Table 3. The comparison table of the newly reassigned score and the average score of 10 experts.
The newly reassigned score ( X i j ¯ ) < 2[2,4)[4,6)[6,8)≥ 8
the average score ( X i j )13579
Table 4. Moran I index of the model.
Table 4. Moran I index of the model.
YearValue of Moran I Index
2005real value0.55
2013real value0.37
2015real value0.45
simulated value0.37
Table 5. The setting of resistance factors and resistance value.
Table 5. The setting of resistance factors and resistance value.
Resistance FactorsRangeResistance Value
Green infrastructure assessment<1.511
1.51–2.8325
2.83–4.05100
4.05–5.33200
5.33–7.28500
Slope0–8%1
8–15%50
15–25%200
> 25%500
Land use typesConstruction landTraffic land0
Others1
Other land3
Grassland50
Farmland100
Water bodies300
Forest land500

Share and Cite

MDPI and ACS Style

Guo, R.; Bai, Y. Simulation of an Urban-Rural Spatial Structure on the Basis of Green Infrastructure Assessment: The Case of Harbin, China. Land 2019, 8, 196. https://0-doi-org.brum.beds.ac.uk/10.3390/land8120196

AMA Style

Guo R, Bai Y. Simulation of an Urban-Rural Spatial Structure on the Basis of Green Infrastructure Assessment: The Case of Harbin, China. Land. 2019; 8(12):196. https://0-doi-org.brum.beds.ac.uk/10.3390/land8120196

Chicago/Turabian Style

Guo, Rong, and Yujing Bai. 2019. "Simulation of an Urban-Rural Spatial Structure on the Basis of Green Infrastructure Assessment: The Case of Harbin, China" Land 8, no. 12: 196. https://0-doi-org.brum.beds.ac.uk/10.3390/land8120196

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