Next Article in Journal
Sustainable Financing through Crowdfunding
Next Article in Special Issue
Long-Term Regional Environmental Risk Assessment and Future Scenario Projection at Ningbo, China Coupling the Impact of Sea Level Rise
Previous Article in Journal
The Research on Temporal–Spatial Distribution and Morphological Characteristics of Ancient Settlements in the Songhua River Basin
Previous Article in Special Issue
Water Bodies’ Cooling Effects on Urban Land Daytime Surface Temperature: Ecosystem Service Reducing Heat Island Effect
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Study on Urban Expansion Using the Spatial and Temporal Dynamic Changes in the Impervious Surface in Nanjing

1
School of Architecture, Nanjing Tech University, Nanjing 210000, China
2
College of Landscape Architecture, Nanjing Forestry University, Nanjing 210000, China
*
Author to whom correspondence should be addressed.
Sustainability 2019, 11(3), 933; https://0-doi-org.brum.beds.ac.uk/10.3390/su11030933
Submission received: 15 January 2019 / Revised: 6 February 2019 / Accepted: 8 February 2019 / Published: 12 February 2019

Abstract

:
Impervious surface area is a key factor affecting urbanization and urban environmental quality. It is of great significance to analysis timely and accurately the dynamic changes of impervious surface for urban development planning. In this study, we use a comprehensive method to extract the time series data on the impervious surface area (ISA) from the multi-temporal Landsat remote sensing images with a high overall accuracy of 90%. The processes and mechanisms of urban expansion at different political administration and direction level in the Nanjing metropolitan area are investigated by using the comprehensive classification method consisting of minimum noise fraction, linear spectral mixture analysis, spectral index, and decision tree classifiers. The expansion of Nanjing is examined by using various ISA indexes and concentric regression analyses. Results indicate that the overall classification accuracy of ISA is higher than 90%. The ISA in Nanjing has dramatically increased in the past three decades from 427.36 km2 to 1780.21 km2 and with a high expansion rate of 0.48 from 2000 to 2005. The city sprawls from monocentric to urban core with multiple subcenters in a concentric structure, and the geometric gravity center of construction land moves southward annually. The stages of urbanization in different district levels and the dynamic changes in different direction levels are influenced by the topographic and economic factors.

1. Introduction

Development of a city during rapid urbanization is often at the expense of natural landscape pervious surfaces which are replaced with an impervious surface composed of mixed clay, asphalt, stone, and metal [1]. The transformation of the earth’s surface from non-urban land to urban land becomes an irreversible process [2,3], thereby generating serious problems, such as farmland loss, heat island effect, flood hazard, and habitat fragmentation [4,5,6,7,8,9,10,11,12,13].
In recent years, many studies have focused on examining the urban land cover and its dynamic changes using the remote sensing (RS) method [14,15,16,17]. However, urban land cover extraction from RS images results in poor accuracy because an urban land is usually a mosaic of different land covers, such as plants, soils, buildings, and water [18,19]. Impervious surface area (ISA) is an essential component of urban landscapes and is the artificial surface that cannot be penetrated by water because of the composition of man-made construction materials in structures, such as road and building roofs [18]. ISA can truly reflect the expansion of urban construction land and can be extracted from RS images [13].
Ridd proposed the vegetation–impervious surface–soil conceptual model to explain urban composition [20]. He regarded urban surface as a combination of vegetation, impervious surface, and soil [20]. Wu and Murray applied the conceptual model of fixed end-member to the Landsat Thematic Mapper Plus (ETM+) image, and they developed the real and operable linear spectral mixture analysis (LSMA) [21]. Various methods, such as spectral mixture analysis, multiple end-member spectral mixture analysis, and ISA index, have been used to improve the precision of ISA extraction [21,22,23,24,25,26]. RS data, such as QuickBird, DMSP-OLS, Terra MODIS, and so on, have also been used to improve accuracy [27,28,29,30,31].
Developing countries, especially China, have become the main driving force of global urbanization in recent years [32,33,34]. The number of large-sized and medium-sized cities are increasing rapidly, and the scale of urban construction is growing. China has experienced an acceleration of urban expansion from the late 1980s [33]. In addition, 70% population of Chinese will live in urban areas by 2030, and 200 million rural immigrants will be received by Chinese cities [35]. As a result, a sharp increase in ISA will occur.
This study uses a large city, Nanjing, as an example. The sprawl of the city and each district is measured using impervious surface, and the consequences of city expansion are analyzed. Nanjing is the old capital of six dynasties in Chinese history and the capital of Jiangsu Province [36]. The urban expansion and development trajectories of Nanjing differ from those of globalizing first-tier metropolitan cities in China in recent years, such as Beijing and Shanghai [37]. Nanjing is a second-tier city with a rapidly growing global economy and is a representative site of research [38].
Some studies have been conducted on the spatiotemporal patterns of urban expansion by monitoring the dynamic changes in ISA [37,39,40,41,42,43]. The current research not only analyzes the status of ISA development in the entire city but also quantifies and compares spatiotemporal patterns of ISA and urban expansion details in different directions and districts. The impacts of the economy on different district development phases are also analyzed. Therefore, this study mainly aims to
(1)
Use a comprehensive method to extract ISA from Landsat images with an overall accuracy of 90%;
(2)
Map locations, density, the geometric center of gravity, and extents of ISA in Nanjing;
(3)
Quantify spatiotemporal patterns of urban expansion in the entire city, and
(4)
Comprehensively compare similarities and differences in the development phases of urbanization under different direction and district levels.
This type of research provides a method to improve the accuracy of urban expansion research data extracted from RS images and provide new understanding on the interactions of socioeconomic systems, different development phases, expansion mechanisms, and urbanization patterns from the perspective of direction, district levels, and the entire city. The other value of this study is to better understand the details of urban expansion in Nanjing, and to provide valuable insights for improving urban planning, management, and sustainability.

2. Materials and Methods

2.1. Study Area

Nanjing is a city with a 2000 years history and was an ancient Chinese capital [36]. Nanjing is the current capital of Jiangsu Province [38]. Thus, Nanjing is the center of politics, economy, and culture of the province and is the core area of the Yangtze River economic belt [44]. Its level of economic development is at the forefront in China [44].
Nanjing metropolis includes 11 districts, namely, Xuanwu, Qinhuai, Gulou, Jianye, Yuhuatai, Qixia, Pukou, Luhe, Jiangning, Lishui, and Gaochun, with a total area of 6587.02 km2 [45]. Its economic growth accelerated after the economic reform, especially in the 1990s. During this period, it was transformed from a medium-sized city to a megacity. Its Gross Domestic Product (GDP) and the municipal population has grown from RMB 16 billion and 5 million people in 1990 to RMB 9720 billion and 8 million people in 2015 [45] (Figure 1).

2.2. Data Collection and Preprocessing

RS data used in this research was mainly from the Landsat Thematic Mapper (TM), Enhanced ETM+, and Operational Land Imager from 1990, 1995, 2000, 2005, 2010, and 2017 (Path: 120, Row: 38) (Table 1). The digital elevation model (DEM) data came from the ASTER Global DEM (GDEM) in Nanjing. The Environmental Systems Research Institute (ESRI) Shapefile of Nanjing administrative units was used. The images were georeferenced in the Universal Transverse Mercator system (WGS84 datum, Zone51 North), and Google Earth images from 1990, 1995, 2000, 2005, 2010, and 2017 were used as the reference data for evaluating ISA accuracy.
Population and GDP data of Nanjing megacity were collected from Nanjing Statistical Yearbook, which is published by China Statistics Press. Ancillary data such as urban planning and management was also collected for examining the different influences of urban development policies on urbanization patterns and rates.

2.3. Mapping ISA Distribution

The proposed comprehensive classification framework for improved LSMA consisted of minimum noise fraction (MNF), LSMA, spectral index, and decision tree classifiers.
Figure 2 illustrates the major steps of the framework. The major steps are explained in detail below:
(1)
Basic RS image operations, such as atmospheric correction and radiometric calibration, were used. ISO unsupervised classification and spectral index of MNDWI were calculated through ENVI 5.2 [46,47].
(2)
MNF transform is a tool for determining the number of bands in the image data, separating the noise in the data, and reducing the computing demand [21]. This method can effectively eliminate noise and reduce the dimension of an image. We used the MNF method to extract good-quality end-members. High albedo, low albedo, and vegetation end-members were extracted on the basis of the geometric vertices in a two-dimensional scatter graph [48].
(3)
Three end-members were used as the ROIs to produce fractional images by using the LSMA approach [48].
(4)
Three fractional images, namely, high-albedo, low-albedo, and green vegetation fraction images, with the Iterative Selforganizing Data Analysis Techniques Algorithm (ISODATA) were used in the decision tree classifier to produce ISA and other data.
(5)
Post processing was conducted, and the results of ISA accuracy were evaluated [49].

2.4. Dynamic Changes in ISA

2.4.1. Density of Dynamic Change in ISA

We mapped the density of ISA to reflect the characteristics of its expansion in time and space. This can intuitively explain the time series of ISA distribution characteristic in the entire city and can illustrate the distribution of human activities. We transformed the original 30 m-sized ISA raster images into the 1 km-sized ISA raster images by using the aggregate method on ArcGIS software [50]. As a result, the scale and direction of the ISA transformation from the rural areas to the urban areas could be determined.
The following formula can be used to calculate the density of ISA:
Density   of   ISA   D O I = I S A T A .
where DOI is the density of ISA in the study area, ISA is the area of impervious surface in the study area, and TA is the total area in the study area.

2.4.2. Center of Gravity of Dynamic Change in ISA

The geometric center of gravity of ISA distribution was measured by ArcGIS to reflect the overall situation of urban expansion in Nanjing. The center of gravity model is composed of 11 district units, and the center of gravity model of ISA is as follows:
Gravity   center   of   ISA   I S A i ( x , y ) = i n ( I i Q ( x i , y i ) ) i n I i ,
where ISAi(x,y) are the geographic coordinates of the center of gravity of ISA in Nanjing, Ii is the area of the impervious surface in different districts, and Q(xi,yi) are the geographic coordinates of the center of gravity of ISA in different districts in Nanjing.

2.4.3. Expansion Rate and Area of Dynamic Change in ISA

Some indexes of ISA change were used to analyze the entire city and its districts. By analyzing the ISA change data in time and space, the expansion of a city can be described. The following variables are calculated to analyze the dynamic change in ISA [51]. The first two formulas are mainly analyzed from the area of ISA, and the latter formulas are mainly analyzed from the growth efficiency of ISA.
Overall ISA changed area OICA = ISA(t2) − ISA(t1)
Annual expansion area AEA = [ISA(t2) − ISA(t1)]/(t2 − t1)
Expansion rate ER = [ISA(t2) − ISA(t1)]/ISA(t1)
where ISA(t) is the area of ISA in different time, t1 is the prior time of t2.

2.4.4. Dynamic Change in ISA in Different Directions

One urban expansion model assumes that urban expansion is a concentric circle structure [44]. This model indicates that a city takes the Central Business District (CBD) as its expansion circle center and expands from the center to suburbs in the form of concentric circles [52,53].
On the basis of this theory, we created the concentric buffer zone structure of the ISA in ArcGIS to analyze the city expansion in different directions.
We defined Xinjiekou CBD as the concentric circle center of Nanjing on the basis of previous research experience [52]. This CBD is located at the intersection of Hanzhong Road and Zhongshan Road. Meanwhile, 5 km buffer zone intervals were created to cover the entire city [44]. An equal interval of 22.5 degrees for each direction was used, and the study area and the concentric circles were divided into eight parts: east, southeast, south, southwest, west, northwest, north, and northeast (Figure 1). The density and expansion rate of ISA were calculated using this method to investigate the spatial expansion of the city in different directions.

2.4.5. Impacts of Topography on Dynamic Change in ISA

Nanjing is in a hilly region [36]. Thus, the impact of terrain on its urban expansion should be investigated. In explaining the relationship between the impact of topography and city expansion, different interval ranges of elevation and slope were used as statistical indicators to describe the distribution and dynamic change in ISA in different topographies. In this case, the elevation was divided into six levels, and the slope was divided into four levels.

2.4.6. Relationship between ISA Dynamics and GDP in Different Districts

Economic development has directly promoted the process of urbanization and determined the level of urbanization [3,54]. However, the relationship between economic and ISA changes at the district level has not been investigated. In this study, we used the regression models to explore the relationship between the GDP data and the ISA change data from 11 districts in Nanjing. We aimed to find the different development phases of five urban core districts and six suburban districts in Nanjing and to develop effective strategies of urban planning on different urban expansion patterns in district level.

3. Results

3.1. Dynamic Change in ISA in the Entire City

The accuracy assessment of the confusion matrix showed that the overall accuracy of the ISA data in Nanjing was more than 90% in 2010 and 2017 (Table 2). The producer’s and user’s accuracies in 2010 were 89.0% and 94.5%, respectively. The producer’s and user’s accuracies in 2017 were 92.0% and 95.5%, respectively. According to previous researchers, the results in other time series were in similar accuracy.
As shown in Figure 3, the density of ISA distribution in Nanjing has drastically changed in the past 30 years. The density of ISA increased from the urban core to the suburban, and the ISA expansion was separated into two parts by the Yangtze River. After 2000, the ISA located in the northern part of the Yangtze River accelerated and expanded. The urban expansion was marked by two axes in the east–west and in the north–south.
The transitions in the geometric gravity center of ISA change can explain the overall situation of urban expansion. From 1990 to 2017, the geometric center of gravity of ISA generally moved southward annually. This transition was divided into two stages. In the first stage, from 1990 to 1995, the center of gravity of ISA gradually moved to southwest; in the second stage, from 1995 to 2017, the center of gravity of ISA gradually shifted to the southeast. The southern expansion of Nanjing was obvious from 2005 to 2010. However, in the past years from 2010 to 2017, the distance between the geometric center of gravity of ISA of the two years became small. Therefore, the construction of ISA increased in the northwest, and the expansion became balanced between the north and south directions (Figure 4).
Overall, the ISA in Nanjing metropolis increased rapidly from 427.36 km2 in 1990 to 1780.21 km2 in 2017 (Table 3).
With regard to area change, the increasing amount of ISA and the annual expansion area reached the maximum in the period of 2010 to 2017; the period of 2000 to 2005 was in second place; the area change in other periods was relatively low. In terms of rate change, the expansion rate reached 0.48 in the period 2000 to 2005, respectively; on the contrary, the rates only reached 0.23 in the period 1995 to 2000. The calculation results of urban expansion showed that, before 2000, the urban expansion of Nanjing was relatively slow; after 2000, Nanjing rapidly developed, especially in the periods 2000 to 2005 and 2010 to 2017 (Table 4).

3.2. Dynamic Change in ISA in District Level

Planners usually distinguish the 11 districts of Nanjing as two types: urban core and suburban districts (Figure 1). Urban core districts include Xuanwu, Qinhuai, Gulou, Jianye, and Yuhuatai, while suburban districts include Qixia, Pukou, Luhe, Jiangning, Lishui, and Gaochun.
Figure 5(a1,a2) indicates obvious differences between the urban core and suburban districts in terms of ISA distribution. At the total tendency of ISA amount, the amount of ISA in all urban core districts was below 60 km2, and the trend of ISA growth gradually slowed down because of the limited space in the urban districts. The ISA amount was over 150 km2, and the trend in ISA growth steadily increased in the suburban districts compared with those in the urban core districts. Although the trend in ISA growth in the urban core districts slowed down, the amount of ISA in Jianye and Yuhuatai in the urban core increased relatively more than that in others. Among the suburban districts, the amount of ISA in Jiangning and Luhe was relatively large owing to the urban development and the large area of their own administration.
Figure 5(b1,b2) indicate that the urban core districts presented higher ISA densities than the suburban districts. The density of ISA in the urban core districts was more than 40%. Among them, the ISA density of Gulou and Qinhuai were the largest. By contrast, Xuanwu exhibited a low density owing to the Zijin mountains and the Xuanwu Lake, which are ecologically protected areas. In the future development of the urban core districts, Jianye and Yuhuatai will further develop to be the main areas of city construction.
In the suburbs, the density of ISA was less than 40%, except in Qixia. The development trend of Qixia was gradually near that of Jianye and Yuhuatai because its geographical location is near the core of the city. Gaochun, Luhe, and Lishui presented a low density of ISA because of their large agricultural and mountainous areas.
The density of ISA in other suburban districts was around 20%–30%. However, the development trends of the ISA density showed an accelerated increase in the period 2005 to 2017. The density will grow rapidly in the future. The results also coincided with the actual situation. In particular, three suburban cores were formed: Xianling suburban core in Qixia District, Dongshan suburban core in Jiangning District, and Jiangbei suburban core in Pukou District.
Figure 6 indicates that the increase of ISA in suburban districts was more than that in urban core districts. In recent years, the increase in urban area is below 1.5 km2/year, and the increase in the suburb is over 5 km2/year. The annual expansion area in Jianye and Yuhuatai increased more than that in three other districts; the annual expansion area development of the three districts in the city core was consistent. Jiangning exhibited the highest value of the annual expansion area. The construction of Jiangning and Lishui has been growing rapidly since 1990, and the construction of Qixia has slowed down since 2005.
After 2005, obvious differences were found in the development of the urban core and suburban districts. The growth point of ISA shifted from the urban area to the suburban, and the suburb entered the main period of development.
Figure 6(b2) shows that, in the period 1990 to 2005, Jianye and Yuhuatai were rapidly developed. However, the annual expansion growth rate of all the urban core districts slowed down after 2005, and the trend of the urban core districts consistently stabilized. In the period 2000 to 2005, the annual expansion rate of all the suburban areas grew fast, whereas the development of Qixia region became slow.

3.3. Dynamic Change in ISA in Different Directions

Figure 7 shows that Nanjing expansion showed a concentric circle characteristic. Compared with that in other cities, the overall density of ISA in Nanjing was lower in the same period. The reason was that Nanjing is a landscape city. The Yangtze River and Zijin Mountain are also in the middle of the city, thereby reducing the total density of ISA. From 1990 to 2017 within the distance of 0 to 5 km, the density of ISA was the highest in each period. The total ISA density in this area increased from 0.63 to 0.83. However, the density of ISA over 40 km was less than 0.3 in all the periods. At around 15 km, the ISA density drastically changed and increased from 0.10 in 1990 to 0.55 in 2017. This result also reflects that, before 2000, the main expansion mostly occurred in the city center in the range of 10 km. After 2000, especially after 2005, the city expansion concentrated on the space between 15 and 20 km.
Different direction figures show that a slight “U” appeared in the east, north, and northwest directions because of the topographic features, such as the Yangtze River and Zijin Mountain, which block the city expansion. However, another wave peak appeared after the U shape. This peak indicates that the urban subcenter was built around the 15 km area in each direction. The figure shows that the west side was not changing dramatically outside the 10 km area. The reason is that the west area belongs to the hilly area of Jianghuai and the urban construction is greatly influenced by the terrain. Figure 7 shows that the development of ISA in southwest Nanjing increased obviously around the 10 km area. This phenomenon is attributed to the formation of the new city center called Hexi. After 2000, the area between 10 and 25 km was influenced by the formation of the new city center. As the economy developed, the density of the 45 km area in southwest increased beyond the city concentric structure. The main reason is that the area between the two cities of the Yangtze River formed the economic zone. The Banqiao city belt is the economic belt between Nanjing and Maanshan. Three peaks in the south direction at 20, 55, and 80 km formed the southern New Town, Lukou Airport, and construction center of Gaochun District, respectively.
Fifteen and 50 km areas of the southeast formed the subcenter of Dongshan City and the construction center of Lishui District.
Analysis of the entire change of ISA density at different distances in Nanjing can reflect the focus area of city construction in different periods. Figure 8 indicates that, at the center of the city around 5 km, the change of ISA density was the largest in the period 1990 to 1995. Then, the rate gradually reduced to approximately 0.01 in 2010 to 2017. This finding shows that the city core formed in the period 1990 to 1995. After 1995, the focus of its construction shifted. Around the 10 km area, the ISA presented the greatest change rate in the period 2000 to 2005. After this period, the change rate of ISA decreased. Therefore, the region focused on its construction before 2005. Around the 15 km area, the change rate of ISA increased from 1990 to 2005 and reached the peak stage in the period 2000 to 2005. In the past years from 2010 to 2017, the change in ISA density was over 0.05 in the distance of 35 km area. Therefore, this trend is the focus of development and construction in recent years.
From Figure 8, the focus area of different directions can be analyzed. After 2000, the strength of construction increased in the 15 km area of the east. The high change value of the 30 to 35 km area in 2005 was due to the rapid construction of Tangshan according to the actual situation. In the northeast, the construction accelerated in the 15 km area in the period 1995 to 2005. At present, the construction focuses on the 25 km area. In the north, northwest, and west directions, the city grew very fast from 2010 to 2017. The main reason is the construction of Jiangbei national development zone and Hexi city core. After 2000, the change rate of ISA increased in the south, southeast, and southwest directions. The change rate of ISA has been the highest in recent years. Therefore, the south area is the main area of urban expansion after 2000.

3.4. Impacts of Topography on Dynamic Change in ISA

Table 5 shows the ISA at different altitude levels from 1990 to 2017. The results show that most of the distribution of ISA was less than 40 m in height. In the 0 to 20 m and 20 to 40 m levels, the total amount increased from 198.71 (176.45 km2) to 987.86 (629.37 km2) from 1990 to 2017 as the altitude increased. The increase of ISA in the same period was small. That is, ISA increased in different altitudes. However, as altitude increased, the expansion rate became small. This result indicates that urban expansion occurs mainly in relatively low-altitude areas and that elevation is a constraint on the urbanization of Nanjing. A similar situation is observed for the slope groups. Table 6 shows that the majority of ISA was located in areas of slope values less than 5 degrees, and these areas were suitable for construction. This result indicates that the relationship between Nanjing expansion and terrain was obvious. The expansion of Nanjing was mainly distributed in the flat terrain, and the elevation and slope were important constraints for the expansion of the city (Figure 9).

3.5. Relationship between ISA Dynamics and GDP in Different Districts

Economic condition is an important driver of urbanization. In this study, we examined the relationship between GDP and ISA by exploring regression models of 11 districts in Nanjing. The result shows a strong positive relation between ISA and GDP at different periods; the coefficients of determination (R2) reached more than 0.88 for all districts (Table 7).
Figure 10 shows that GDP was more sensitive to ISA in the suburban districts than the urban core districts in Nanjing. In the urban core districts, the changes in ISA in the Jianye and Yuhuatai areas were sensitive to the changes in GDP. Compared with other mature urban core districts, the changes in ISA in Jianye and Yuhuatai presented a great development space in the next few years. In the suburb districts, the driving force of Qixia’s economic development to ISA was the smallest, followed by that of Jiangning. The ISA of Pukou, Luhe, and Gaochun was significantly affected by GDP, and the state of the total ISA increased significantly.

4. Discussion

4.1. Improving the Accuracy of Urban Expansion Research Data

Urban expansion analysis has been usually conducted by extracting the data of urban construction land use [14,15,16,17]. However, the relatively coarse spatial resolution of Landsat images, mixed pixel, and complex urban land cover composition reduces the accuracy of urban construction land classification [19]. In this study, the comprehensive method of a pixel-based classification was used to extract ISA. This method is more accurate than the construction land that is directly extracted from Landsat files. We improved LSMA, which has been validated by scholars as having a valid ISA extraction accuracy, through a variety of methods, such as Modified Normalized Difference Water Index (MNDWI) and decision tree. The improved method can effectively extract ISA with an overall accuracy of 90%.
This method is effective in producing accurate distributions of ISA, which is especially necessary for testing urban expansion. Given that this study explores the expansion of land use in each district and different directions in Nanjing, the improvement in the accuracy of ISA plays a very crucial role in the entire study. Accurate extraction of ISA allows a clear analysis of ISA changes in each area and direction of the urban expansion.

4.2. Analyzing City Expansions at Different Spatial Scales and Direction Levels

Rational control of the pace and direction of the expansion of construction land in metropolitan areas has been the focus of Chinese urbanization and land management policies. Research on the expansion mechanism of construction land through different directions and spatial scales can provide a scientific basis for urban development management.
This study established a relatively comprehensive analysis system of urban expansion. From the macroscopic view of the overall urban expansion, the density changes, and center of gravity of ISA were judged, and, thus, the overall trend of urban expansion was analyzed. From the perspective of micro-regional and directional levels, we studied the expansion process and change trend of ISA in every district and direction of Nanjing and analyzed the key areas and directions of urban expansion in the future. The date base of multi-center development in Nanjing can be provided for urban planning and construction needs.
On the basis of this systematic research model and institutional analysis of the expansion of construction land use in Nanjing, additional research experience and resources can be provided. The results can be used to summarize a set of characteristics suitable for Chinese urbanization and land expansion construction.

4.3. Effects of Topographic and GDP Factors on Urban Expansion

Topography and GDP are two important factors affecting urban ISA expansion in Nanjing. This research confirmed this fact. Examining the relationships between topography and ISA at different periods shows that topography is still an important constraint factor affecting the expansion of the city even with developed engineering technology.
Using regression analysis parsing the relationship between ISA and GDP in different districts in Nanjing indicates that ISA is positively related to GDP, especially in districts at the development stage. GDP in most of these districts significantly influences urban expansion. Those districts are usually located on the urban fringe area. By contrast, GDP in other types of districts insignificantly influences the ISA expansion.
On the basis of the analysis of topography and economy, the future development of different districts in Nanjing can be analyzed and the current and future development focuses of different areas can be summarized. Effective strategies for urban planning, policy and improved service for urban citizens can be provided.

4.4. Social and Political Factors on Urban Expansion

Social and political factors are also important factors affecting ISA expansion in Nanjing. The household registration system and the reform and opening up policies are the main factors affecting the expansion of Chinese cities in the 21st century.
In 1958, China launched the household registration system, which divided the identity of urban residents and rural residents, imposed strict restrictions on the free movement of people [55]. Since the 1990s, with the reform of the household registration system, the urban and rural residents have flowed with each other, but the inequality of resources between urban and rural areas has caused the population concentrated in the central urban area which can provide the relatively-perfected educational and medical conditions. Many rural residents have flocked to the cities, and it indirectly led to the emergence of a basic pattern of city expansion in China, and Nanjing is no exception. In recent years, to attract talents and improve the competitiveness of Nanjing, the local government has formulated a new policy that people with the bachelor’s degree or above can be directly settled down, the possibility of sustainable expansion of the city’s scale has been provided.
Since the reform and opening up, Chinese society has entered a new stage of transformation and development [56]. Under the guidance of market economic policy, Chinese enterprises have begun to change from government-led enterprises to market-oriented enterprises. Meanwhile, the hot spots of foreign investment in China are the Yangtze River Delta and the Pearl River Delta, and Nanjing, as a hot investment city in the Yangtze River Delta, has accelerated its urban expansion due to the influx of large quantities of local and foreign capital, and it is embodied in the agglomeration development mode with industrial parks. Those industrial parks have brought the agglomeration of industry and population, accelerated the development of cities, and promoted the change of urban spatial structure.

5. Conclusions

Compared with other European and American megacities, Chinese megacities exhibit different urbanization processes, patterns, and mechanisms. This study used Nanjing as a case study to assess and analyze urban expansion by using a comprehensive method to extract and analyze ISA.
From the results, the following conclusions were drawn:
(1)
The proposed comprehensive method can effectively improve the accuracy of ISA extraction. In this case study, the overall accuracy of ISA in Nanjing was over 90%. The ISA of the city was comprehensively analyzed by macro analysis of urban expansion and microcosmic analysis of ISA in different districts and directions.
(2)
Urban development of Nanjing has been expanding rapidly, slowing down, and expanding rapidly again. The annual expansion area from 23.64 km2/year in 1990 to 1995 increased to 63.77 km2/year in 2000 to 2005, slowed down in 2005 to 2010, and then increased again in 2010 to 2017 with the annual expansion area of 72.61 km2/year. The region of high ISA density continues to grow. The center of gravity of ISA obviously moves southeast, and the overall direction of expansion is relatively balanced in recent years.
(3)
Types and trends of ISA changes in urban core and suburban districts obviously differ. In recent years, the increase in the area and density of ISA in the urban core districts is relatively slow, and the change in ISA in the suburban districts is the opposite. At the same time, the different areas and rates of ISA change from the different directions of urban CBD and are greatly influenced by the urban development policy and geographical factors.
(4)
The case of Nanjing indicates that topography and economy are important influencing factors of urbanization. In different districts at dissimilar stages of development, the changes and trends of ISA differ. The increase of the ISA in the districts at the stage of development is positively proportional to GDP. Based on the regression analysis of the economy and ISA, the future expansion trends of each district can be approximately estimated.
(5)
Urban expansion in Nanjing represents a common phenomenon in the development of the second-tier cities in China. The ISA analysis of Nanjing’s expansion indicates that Nanjing is still in a period of rapid urban expansion, the ISA density of the original urban central district continues to increase, and the overall ISA area of the city continues to expand. These are also the basic characteristics of the expansion of Nanjing and the other second-tier cities in China. Meanwhile, Nanjing’s urban expansion can provide some development experience for the same type of second-tier cities and the third-tier cities in China, and the government can formulate corresponding planning strategies according to the development of the ISA distribution and the future ISA expansion trends, analyze the urban spatial form, improve the level of urban management, and make rational use of urban space.

Author Contributions

Z.W. conceived and designed the experiments; Y.Q. made contributions to the design, data processing; Z.W. wrote the paper; All authors had read and approved the manuscript.

Acknowledgments

The present study is supported by the National Natural Science Foundation of China (Grant No. 31570703). The authors are grateful for the technical support of Professor Shiguang Shen.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Strahler, A.H.; Woodcock, C.E.; Smith, J.A. On the nature of models in remote sensing. Remote Sens. Environ. 1986, 20, 121–139. [Google Scholar] [CrossRef]
  2. Carrion-Flores, C.; Irwin, E.G. Determinants of residential land-use conversion and sprawl at the rural-urban fringe. Am. J. Agric. Econ. 2004, 86, 889–904. [Google Scholar] [CrossRef]
  3. Seto, K.C.; Fragkias, M.; Güneralp, B.; Reilly, M.K. A meta-analysis of global urban land expansion. PLoS ONE 2011, 6, e23777. [Google Scholar] [CrossRef]
  4. Bhatta, B.; Saraswati, S.; Bandyopadhyay, D. Urban sprawl measurement from remote sensing data. Appl. Geogr. 2010, 30, 731–740. [Google Scholar] [CrossRef]
  5. Alberti, M. Maintaining ecological integrity and sustaining ecosystem function in urban areas. Curr. Opin. Environ. Sustain. 2010, 2, 178–184. [Google Scholar] [CrossRef]
  6. Grimm, N.B.; Faeth, S.H.; Golubiewski, N.E.; Redman, C.L.; Wu, J.; Bai, X.; Briggs, J.M. Global change and the ecology of cities. Science 2008, 319, 756–760. [Google Scholar] [CrossRef]
  7. Foley, J.A.; DeFries, R.; Asner, G.P.; Barford, C.; Bonan, G.; Carpenter, S.R.; Chapin, F.S.; Coe, M.T.; Daily, G.C.; Gibbs, H.K. Global consequences of land use. Science 2005, 309, 570–574. [Google Scholar] [CrossRef]
  8. Dupras, J.; Marull, J.; Parcerisas, L.; Coll, F.; Gonzalez, A.; Girard, M.; Tello, E. The impacts of urban sprawl on ecological connectivity in the montreal metropolitan region. Environ. Sci. Policy 2016, 58, 61–73. [Google Scholar] [CrossRef]
  9. Chan, K.M.; Vu, T.T. A landscape ecological perspective of the impacts of urbanization on urban green spaces in the klang valley. Appl. Geogr. 2017, 85, 89–100. [Google Scholar]
  10. Dupras, J.; Alam, M. Urban sprawl and ecosystem services: A half century perspective in the montreal area (quebec, canada). J. Environ. Policy Plan. 2015, 17, 180–200. [Google Scholar] [CrossRef]
  11. Arnold Jr, C.L.; Gibbons, C.J. Impervious surface coverage: The emergence of a key environmental indicator. J. Am. Plan. Assoc. 1996, 62, 243–258. [Google Scholar] [CrossRef]
  12. Myronidis, D.; Loannou, K. Forecasting the urban expansion effects on the design storm hydrograph and sediment yield using artificial neural networks. Water 2019, 11, 31. [Google Scholar] [CrossRef]
  13. Myronidis, D.; Stathis, D.; Sapountzis, M. Post-Evaluation of flood hazards induced by former artificial interventions along a coastal mediterranean settlement. J. Hydrol. Eng. 2016, 21, 05016022. [Google Scholar] [CrossRef]
  14. El Garouani, A.; Mulla, D.J.; El Garouani, S.; Knight, J. Analysis of urban growth and sprawl from remote sensing data: Case of fez, morocco. Int. J. Sustain. Built Environ. 2017, 6, 160–169. [Google Scholar] [CrossRef]
  15. Cabral, A.I.; Costa, F.L. Land cover changes and landscape pattern dynamics in senegal and guinea bissau borderland. Appl. Geogr. 2017, 82, 115–128. [Google Scholar]
  16. Hassan, M.M. Monitoring land use/land cover change, urban growth dynamics and landscape pattern analysis in five fastest urbanized cities in bangladesh. Remote Sens. Appl. Soc. Environ. 2017, 7, 69–83. [Google Scholar] [CrossRef]
  17. Dahal, K.R.; Benner, S.; Lindquist, E. Urban hypotheses and spatiotemporal characterization of urban growth in the treasure valley of idaho, USA. Appl. Geogr. 2017, 79, 11–25. [Google Scholar]
  18. Lu, D.; Weng, Q. Spectral mixture analysis of the urban landscape in indianapolis with landsat etm+ imagery. Photogram. Eng. Remote Sens. 2004, 70, 1053–1062. [Google Scholar] [CrossRef]
  19. Schneider, A. Monitoring land cover change in urban and peri-urban areas using dense time stacks of landsat satellite data and a data mining approach. Remote Sens. Environ. 2012, 124, 689–704. [Google Scholar] [CrossRef]
  20. Ridd, M.K. Exploring a vis (vegetation-impervious surface-soil) model for urban ecosystem analysis through remote sensing: Comparative anatomy for cities. Int. J. Remote Sens. 1995, 16, 2165–2185. [Google Scholar] [CrossRef]
  21. Wu, C.; Murray, A.T. Estimating impervious surface distribution by spectral mixture analysis. Remote Sens. Environ. 2003, 84, 493–505. [Google Scholar] [CrossRef]
  22. Lu, D.; Moran, E.; Hetrick, S. Detection of impervious surface change with multitemporal landsat images in an urban–rural frontier. ISPRS J. Photogram. Remote Sens. 2011, 66, 298–306. [Google Scholar] [CrossRef]
  23. Xu, H. Rule-based impervious surface mapping using high spatial resolution imagery. Int. J. Remote Sens. 2013, 34, 27–44. [Google Scholar] [CrossRef]
  24. Okujeni, A.; van der Linden, S.; Hostert, P. Extending the vegetation–impervious–soil model using simulated enmap data and machine learning. Remote Sens. Environ. 2015, 158, 69–80. [Google Scholar] [CrossRef]
  25. Carlson, T.N.; Arthur, S.T. The impact of land use—land cover changes due to urbanization on surface microclimate and hydrology: A satellite perspective. Glob. Planet. Chang. 2000, 25, 49–65. [Google Scholar] [CrossRef]
  26. Liu, S.; Gu, G. Improving the impervious surface estimation from hyperspectral images using a spectral-spatial feature sparse representation and post-processing approach. Remote Sens. 2017, 9, 456. [Google Scholar]
  27. Yang, L.; Huang, C.; Homer, C.G.; Wylie, B.K.; Coan, M.J. An approach for mapping large-area impervious surfaces: Synergistic use of landsat-7 etm+ and high spatial resolution imagery. Can. J. Remote Sens. 2003, 29, 230–240. [Google Scholar] [CrossRef]
  28. Alahmadi, M.; Atkinson, P.; Martin, D. Fine spatial resolution residential land-use data for small-area population mapping: A case study in riyadh, saudi arabia. Int. J. Remote Sens. 2015, 36, 4315–4331. [Google Scholar] [CrossRef]
  29. Friedl, M.A.; McIver, D.K.; Hodges, J.C.; Zhang, X.; Muchoney, D.; Strahler, A.H.; Woodcock, C.E.; Gopal, S.; Schneider, A.; Cooper, A. Global land cover mapping from modis: Algorithms and early results. Remote Sens. Environ. 2002, 83, 287–302. [Google Scholar] [CrossRef]
  30. Guo, W.; Lu, D.; Wu, Y.; Zhang, J. Mapping impervious surface distribution with integration of snnp viirs-dnb and modis ndvi data. Remote Sens. 2015, 7, 12459–12477. [Google Scholar] [CrossRef]
  31. Weng, Q. Remote sensing of impervious surfaces in the urban areas: Requirements, methods, and trends. Remote Sens. Environ. 2012, 117, 34–49. [Google Scholar] [CrossRef]
  32. Li, H.; Wei, Y.D.; Liao, F.H.; Huang, Z. Administrative hierarchy and urban land expansion in transitional china. Applied Geography 2015, 56, 177–186. [Google Scholar] [CrossRef]
  33. Schneider, A.; Mertes, C. Expansion and growth in chinese cities, 1978–2010. Environ. Res. Lett. 2014, 9, 024008. [Google Scholar] [CrossRef]
  34. Liu, Z.; He, C.; Zhang, Q.; Huang, Q.; Yang, Y. Extracting the dynamics of urban expansion in china using dmsp-ols nighttime light data from 1992 to 2008. Landsc. Urban Plan. 2012, 106, 62–72. [Google Scholar] [CrossRef]
  35. Chinese national population development plan (2016–2030). Available online: http://www.gov.cn/zhengce/content/2017-01/25/content_5163309.htm (accessed on 28 November 2017).
  36. Yuan, F.; Gao, J.; Wu, J. Nanjing-an ancient city rising in transitional china. Cities 2016, 50, 82–92. [Google Scholar] [CrossRef]
  37. Yin, J.; Yin, Z.; Zhong, H.; Xu, S.; Hu, X.; Wang, J.; Wu, J. Monitoring urban expansion and land use/land cover changes of shanghai metropolitan area during the transitional economy (1979–2009) in china. Environ. Monit. Assess. 2011, 177, 609–621. [Google Scholar] [CrossRef]
  38. Xu, C.; Liu, M.; Zhang, C.; An, S.; Yu, W.; Chen, J.M. The spatiotemporal dynamics of rapid urban growth in the nanjing metropolitan region of china. Landsc. Ecol. 2007, 22, 925–937. [Google Scholar] [CrossRef]
  39. Shi, L.; Ling, F.; Ge, Y.; Foody, G.M.; Li, X.; Wang, L.; Zhang, Y.; Du, Y. Impervious surface change mapping with an uncertainty-based spatial-temporal consistency model: A case study in wuhan city using landsat time-series datasets from 1987 to 2016. Remote Sens. 2017, 9, 1148. [Google Scholar] [CrossRef]
  40. Gao, F.; de Colstoun, E.B.; Ma, R.; Weng, Q.; Masek, J.G.; Chen, J.; Pan, Y.; Song, C. Mapping impervious surface expansion using medium-resolution satellite image time series: A case study in the yangtze river delta, china. International Journal of Remote Sensing 2012, 33, 7609–7628. [Google Scholar] [CrossRef]
  41. Zhang, L.; Weng, Q. Annual dynamics of impervious surface in the pearl river delta, china, from 1988 to 2013, using time series landsat imagery. ISPRS J. Photogram. Remote Sens. 2016, 113, 86–96. [Google Scholar] [CrossRef]
  42. Van de Voorde, T.; Jacquet, W.; Canters, F. Mapping form and function in urban areas: An approach based on urban metrics and continuous impervious surface data. Landsc. Urban Plan. 2011, 102, 143–155. [Google Scholar] [CrossRef]
  43. Scott, D.; Petropoulos, G.P.; Moxley, J.; Malcolm, H. Quantifying the physical composition of urban morphology throughout wales based on the time series (1989–2011) analysis of landsat tm/etm+ images and supporting gis data. Remote Sens. 2014, 6, 11731–11752. [Google Scholar] [CrossRef]
  44. Tian, G.; Jiang, J.; Yang, Z.; Zhang, Y. The urban growth, size distribution and spatio-temporal dynamic pattern of the yangtze river delta megalopolitan region, china. Ecol. Model. 2011, 222, 865–878. [Google Scholar] [CrossRef]
  45. Nanjing Statistical Bureau. Nanjing Statistical Yearbook; China Statistics Press: Beijing, China, 2017.
  46. Han-Qiu, X. A study on information extraction of water body with the modified normalized difference water index (mndwi). J. Remote Sens. 2005, 5, 589–595. [Google Scholar]
  47. Xu, X.; Zhou, Y.; Ning, Y. Urban Geography; High Education Publisher: Beijing, China, 2009. [Google Scholar]
  48. Lu, D.; Weng, Q. Use of impervious surface in urban land-use classification. Remote Sens. Environ. 2006, 102, 146–160. [Google Scholar] [CrossRef]
  49. Congalton, R.G.; Green, K. Assessing the Accuracy of Remotely Sensed Data: Principles and Practices; CRC Press: Boca Raton, FL, USA, 2008. [Google Scholar]
  50. Li, L.; Lu, D.; Kuang, W. Examining urban impervious surface distribution and its dynamic change in hangzhou metropolis. Remote Sens. 2016, 8, 265. [Google Scholar] [CrossRef]
  51. Kuang, W.; Chi, W.; Lu, D.; Dou, Y. A comparative analysis of megacity expansions in china and the us: Patterns, rates and driving forces. Landsc. Urban Plan. 2014, 132, 121–135. [Google Scholar] [CrossRef]
  52. Chen, J.; Gao, J.; Chen, W. Urban land expansion and the transitional mechanisms in Nanjing, China. Habitat Int. 2016, 53, 274–283. [Google Scholar] [CrossRef]
  53. Wei, Y.D.; Ye, X. Urbanization, urban land expansion and environmental change in china. Stoch. Environ. Res. Risk Assess. 2014, 28, 757–765. [Google Scholar] [CrossRef]
  54. Bai, X.; Chen, J.; Shi, P. Landscape urbanization and economic growth in china: Positive feedbacks and sustainability dilemmas. Environ. Sci. Technol. 2011, 46, 132–139. [Google Scholar] [CrossRef]
  55. Wu, X.G.; Zheng, B.D. Household registration, urban status attainment, and social stratification in China. Res. Soc. Stratific. Mobil. 2018, 53, 40–49. [Google Scholar] [CrossRef]
  56. Lv, J.; Yang, B.D.; Yang, Y.J.; Zhang, Z.H.; Liu, G.J. Spatial patterns of China’s major cities and their evolution mechanisms during the past decades of Reform and Opening Up. Procedia Eng. 2017, 198, 915–925. [Google Scholar] [CrossRef]
Figure 1. Location of the study area.
Figure 1. Location of the study area.
Sustainability 11 00933 g001
Figure 2. Flowchart for mapping impervious surface distribution.
Figure 2. Flowchart for mapping impervious surface distribution.
Sustainability 11 00933 g002
Figure 3. Density of dynamic change in impervious surface area (ISA) in Nanjing from 1990 to 2017.
Figure 3. Density of dynamic change in impervious surface area (ISA) in Nanjing from 1990 to 2017.
Sustainability 11 00933 g003
Figure 4. Center of gravity of dynamic change in ISA in Nanjing from 1990 to 2017.
Figure 4. Center of gravity of dynamic change in ISA in Nanjing from 1990 to 2017.
Sustainability 11 00933 g004
Figure 5. Amount and density of ISA in different districts in Nanjing from 1990 to 2017. (Note: (a1) Amount of ISA in urban core districts; (a2) Amount of ISA in suburban districts; (b1) Density of ISA in urban core districts; (b2) Density of ISA in suburban districts).
Figure 5. Amount and density of ISA in different districts in Nanjing from 1990 to 2017. (Note: (a1) Amount of ISA in urban core districts; (a2) Amount of ISA in suburban districts; (b1) Density of ISA in urban core districts; (b2) Density of ISA in suburban districts).
Sustainability 11 00933 g005
Figure 6. Annual growth amount and rates of ISA in different districts in Nanjing from 1990 to 2017. (Note: (a1) Annual growth amount of ISA in urban core districts; (a2) Annual growth amount of ISA in suburban districts; (b1) Annual growth rate of ISA in urban core districts; (b2) Annual growth rate of ISA in suburban districts).
Figure 6. Annual growth amount and rates of ISA in different districts in Nanjing from 1990 to 2017. (Note: (a1) Annual growth amount of ISA in urban core districts; (a2) Annual growth amount of ISA in suburban districts; (b1) Annual growth rate of ISA in urban core districts; (b2) Annual growth rate of ISA in suburban districts).
Sustainability 11 00933 g006
Figure 7. Density of ISA in different directions in Nanjing from 1990 to 2017. (Note: (a1) Entire area; (a2) East; (a3) Northeast; (a4) North; (a5) Northwest; (a6) West; (a7) Southwest; (a8) South; (a9) Southeast).
Figure 7. Density of ISA in different directions in Nanjing from 1990 to 2017. (Note: (a1) Entire area; (a2) East; (a3) Northeast; (a4) North; (a5) Northwest; (a6) West; (a7) Southwest; (a8) South; (a9) Southeast).
Sustainability 11 00933 g007aSustainability 11 00933 g007b
Figure 8. Density of change in ISA in different directions in Nanjing from 1990 to 2017. (Note: (a1) Entire area; (a2) East; (a3) Northeast; (a4) North; (a5) Northwest; (a6) West; (a7) Southwest; (a8) South; (a9) Southeast).
Figure 8. Density of change in ISA in different directions in Nanjing from 1990 to 2017. (Note: (a1) Entire area; (a2) East; (a3) Northeast; (a4) North; (a5) Northwest; (a6) West; (a7) Southwest; (a8) South; (a9) Southeast).
Sustainability 11 00933 g008aSustainability 11 00933 g008b
Figure 9. Impacts of Topography in ISA. (Note: (a) Elevation; (b) Slope).
Figure 9. Impacts of Topography in ISA. (Note: (a) Elevation; (b) Slope).
Sustainability 11 00933 g009
Figure 10. Relationships between ISA and GDP in different directions in Nanjing from 1990 to 2017. (Note: (a) Relationships between ISA and GDP in urban core districts; (b) Relationships between ISA and GDP in suburban districts).
Figure 10. Relationships between ISA and GDP in different directions in Nanjing from 1990 to 2017. (Note: (a) Relationships between ISA and GDP in urban core districts; (b) Relationships between ISA and GDP in suburban districts).
Sustainability 11 00933 g010
Table 1. Summary of data used in research.
Table 1. Summary of data used in research.
Data TypesDatasets
Remote sensing dataTimePath/rowImage type
1990-7-11Landsat 5-TM120/38
1995-10-13Landsat 5-TM120/38
2000-6-12Landsat 5-TM120/38
2005-4-7Landsat 5-TM120/38
2010-4-5Landsat 7-ETM120/38
2017-7-21Landsat 8-OLI120/38
DEM dataASTER GDEM with 30m spatial resolution
Boundary filesNanjing administration boundary shape data
GDP dataGDP data in 1995, 2000, 2005, 2010, and 2015 year at district level in Nanjing
Ancillary dataUrban planning and policies of Nanjing
Table 2. Accuracy assessment of impervious surface mapping result.
Table 2. Accuracy assessment of impervious surface mapping result.
20102017
ISAOthersProducer’s accuracy (%)User’s accuracy (%)ISAOthersProducer’s accuracy (%)User’s accuracy (%)
ISA891189.089.092891.192.0
Others1118994.594.5919196.095.5
Overall accuracy (%)92.792.6
Table 3. Impervious surface area (ISA) data in Nanjing from 1990 to 2017.
Table 3. Impervious surface area (ISA) data in Nanjing from 1990 to 2017.
Year199019952000200520102017
ISA (km2)427.36545.57668.72987.581271.931780.21
Table 4. Analysis of dynamic change in ISA in Nanjing from 1990 to 2017.
Table 4. Analysis of dynamic change in ISA in Nanjing from 1990 to 2017.
1990–19951995–20002000–20052005–20102010–2017
Overall ISA changed area (km2)118.22123.15318.86284.34508.29
Annual expansion area (km2/year)23.6424.6363.7756.8772.61
Expansion rate0.280.230.480.290.40
Table 5. Impacts of elevation on dynamic change in ISA in Nanjing from 1990 to 2017.
Table 5. Impacts of elevation on dynamic change in ISA in Nanjing from 1990 to 2017.
0–2020–4040–6060–8080–100100–150150–200>200
km2%km2%km2%km2%km2%km2%km2%km2%
1990198.7145.40176.4540.3144.1210.0811.132.543.800.872.980.680.460.110.050.01
1995251.7845.05227.9340.7854.599.7714.672.625.270.943.880.690.660.120.080.01
2000329.5948.14268.4239.2161.058.9215.082.205.420.794.030.590.920.130.140.02
2005518.7851.32371.4336.7483.848.2920.612.047.770.776.350.631.790.180.400.04
2010689.5852.96463.1835.57104.948.0624.421.889.120.707.850.602.340.180.590.05
2017987.8654.19629.3734.53147.608.1032.441.7811.620.649.860.543.130.171.020.06
Table 6. Impacts of slope on dynamic change in ISA in Nanjing from 1990 to 2017.
Table 6. Impacts of slope on dynamic change in ISA in Nanjing from 1990 to 2017.
0–55–1010–1515–20>20
km2%km2%km2%km2%km2%
1985121.7241.67119.4440.8938.1113.059.163.143.681.26
1990184.0943.01173.6840.5752.3912.2412.452.915.441.27
1995235.143.03222.2440.6766.7512.2215.592.856.741.23
2000294.1843.93270.8940.4579.1811.8218.032.697.431.11
2005446.1345.10395.239.95111.4811.2725.112.5411.211.13
2010585.9745.98504.6939.61138.1310.8430.952.4314.551.14
2017836.2546.88703.139.41185.6610.4140.272.2618.651.05
Table 7. Relationship between ISA dynamics and GDP in different districts in Nanjing.
Table 7. Relationship between ISA dynamics and GDP in different districts in Nanjing.
DistrictFormulaR2
Xuanwuy = 2.9239 ln(x) + 16.240.9672
Qinhuaiy = 2.8946 ln(x) + 20.1980.9747
Gulouy = 0.9895 ln(x) + 37.5880.9711
Yuhuataiy = 10.965 ln(x) - 10.4460.9864
Jianyey = 8.5131 ln(x) + 11.010.9683
Qixiay = 26.873 ln(x) − 16.0890.9636
Gaochuny = 38.893 ln(x) − 82.0510.9259
Pukouy = 44.196 ln(x) − 92.630.8891
Luhey = 54.399 ln(x) − 91.6790.9044
Lishuiy = 48.067 ln(x) − 119.110.9529
Jiangningy = 104.21 ln(x) − 312.020.9828

Share and Cite

MDPI and ACS Style

Qian, Y.; Wu, Z. Study on Urban Expansion Using the Spatial and Temporal Dynamic Changes in the Impervious Surface in Nanjing. Sustainability 2019, 11, 933. https://0-doi-org.brum.beds.ac.uk/10.3390/su11030933

AMA Style

Qian Y, Wu Z. Study on Urban Expansion Using the Spatial and Temporal Dynamic Changes in the Impervious Surface in Nanjing. Sustainability. 2019; 11(3):933. https://0-doi-org.brum.beds.ac.uk/10.3390/su11030933

Chicago/Turabian Style

Qian, Yanping, and Zhen Wu. 2019. "Study on Urban Expansion Using the Spatial and Temporal Dynamic Changes in the Impervious Surface in Nanjing" Sustainability 11, no. 3: 933. https://0-doi-org.brum.beds.ac.uk/10.3390/su11030933

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