1. Introduction
Urbanization is a complicated phenomenon that is accompanied by population movements, changes in the rural-urban divide, and flows of matter and energy. In recent years, urbanization has expanded into a significant spatial phenomenon on a global scale. In 1950, 25% of the world’s population lived in cities [
1] and currently, urban areas are home to 50% of the world’s population [
2]. Especially in less developed countries, the urban population will increase to 56.4% by 2030 [
3]. Urbanization is one of the most significant factors in human-induced land cover and land-use change, even though urban areas make up only a small portion of the entire Earth's surface [
4]. Urbanization and its evolution have been studied using a range of spatial-temporal scales, using remote sensing techniques as the primary source of consistent and continuous data [
4,
5,
6]. The urban system is intricate, with many interconnected components. A key component of the urban system is the impervious surface, which is made up mostly by buildings, roads, and other infrastructure and serves as population and economic hubs [
6].
Urban entity expansion needs to be continuously monitored due to worsening conditions of overcrowding, housing shortages, inadequate infrastructure, and growing ecological issues [
2]. Therefore, it is necessary to implement effective spatial planning and monitoring measures to reduce the negative effects of urban entity expansion. It is especially well suited to anticipate those processes in urban environments using multi-temporal and multi-spectral data. Different size cities can be seen in the geographical space as a result of urbanization, especially major cities that have been rapidly urbanizing during the last two decades in China [
7]. Sustainable urban development in China requires maintaining a balance between human activities, the environment, population, and urbanization by monitoring changes in urban entities in order to inform policy decisions on urban environmental management and planning. Major changes brought about by urbanization in China include population and economic transition from the countryside to urban areas. As a country of the global south, China experiences the largest population concentration around urban areas and as a result of these developments, it also bought about more socioeconomic benefits and adverse environmental impacts. Following the implementation of economic reforms, the country underwent a significant rural-to-urban shift, with an increase in the concentration of inhabitants in and surrounding urban regions from 1978 to 2012, from 17.9% to 52.6%, respectively [
8]. Recent decades have seen a marked increase in the rate of urbanization due to migration and the expansion of urban boundaries. In contrast, the 2010 population census in China showed that about 50% of the country’s population was categorized as urban. The turning point for China came after 2011 when for the first time in its history, the majority of the country’s settlements were classified as urban. Today’s new stage of development in China incorporates urbanization as a key component and direction. The second stage of China’s urbanization, or the stage of steady expansion, has been reached, according to the theory of urbanization development [
9].
Identifying urban areas based on distinct urban characteristics might be a reliable entry point for analyzing and predicting urban expansion to address future sustainable development issues. Urban areas are frequently thought of as impervious areas filled by roads, buildings, and other constructions [
4,
10]. Typically, the qualities of man-made structures are used to demarcate urban entities [
11]. However, it is ineffective in evaluating how human activities are concentrated together [
12]. Additionally, it is essential to illustrate the other socioeconomic functions of urban dwellers using various urban structures [
13]. However, it is challenging to apply this to wider spatial-temporal coverage. Urban entities, which have a concentration of population and human activities, are the areas that spread positive externality. The concentration of human socioeconomic activities, continuity, and physical structures across a larger spatial coverage are the three main features of urban entities. Applying economic-based variables to characterize urban entities, such as the vacancy rate of urban housing, will increase the accuracy of the urban entity definition. Upon this, some researchers [
14,
15] advocated using an integrated definition in conjunction with socioeconomic factors and physical structures as a more effective way to define urban boundaries and their expansion. Government institutes normally use population density or size as the primary criterion for designating urban entities. The two most popular datasets that utilize for the identification of urban entities are LandScan, and World Pop datasets [
16]. Since different countries have diverse population and density characteristics, it is challenging to extract cohesive urban entities using population size variances or density. Urban entity mapping on a wider scale has plenty many options when using remotely sensed data. Nowadays, Landsat, Sentinel, and Moderate Resolution Imaging Spectrometer (MODIS) data are typically used to produce urban entity maps. These datasets are primarily intended to reveal the temporal and spatial distributions of urban impervious areas. But prevailing urban laws, suggest that urban expansion is a complicated process and urban entities shouldn't be primarily identified upon their impervious surfaces [
17]. Although daylight satellite imagery can be used to track land use changes, their spectrally varied nature makes it difficult to use them to study urbanization dynamics [
18,
19]. However, Defense Meteorological Satellite Program’s Operational Line-scan System (DMSP/OLS), and Suomi National Polar Partnership-Visible Infrared Imaging Radiometer Suite instrument onboard (SNPP-VIIRS) nighttime light (NTL) images offer a unique proxy to extract urban dynamics. As a result, since the 1990s DMSP/OLS images have been utilized for different urban studies, including demographic dynamics, city lights, electricity consumption, urbanization, gross domestic products in urban areas, and greenhouse gas emissions [
18,
19,
20,
21].
Urban centers will become difficult, and inconvenient places for commuters, service seekers, and urban dwellers as a result of a lack of space and amenities. Given that it allows urban geographers and planners to map urban areas, and has a significant impact on resource management and socioeconomic development in cities, the evaluation of urban entities is a main urban research direction. As a rapidly urbanizing country, numerous studies have been conducted to extract urban entities from NTL data in China. While few scholars [
22] tried to introduce correcting techniques for multi-satellite NTL stable light data to identify urban expansion in China others [
23] quantified the dynamics of urban entities. Most researchers have used new methods, and new-generation NTL data for urban entity extraction as well as compared their findings with DMSP-OLS data [
7,
18,
24,
25]. Urban extraction using DMSP/OLS data and NTL intensity were also key focusing areas in previous researchers [
26,
27]. Previous researchers have also used neighborhood statistics, NDVI, and local-optimized threshold analysis using nighttime stable light data and VIIRS data to trace urban growth in China [
28,
29,
30,
31]. Few scholars explored the characteristics and trends of urban development in urban agglomerations using DMSP/OLS data [
32]. Predicting and quantifying urban growth in coastal cities using DMSP-OLS was another research focus of some scholars [
33,
34,
35,
36]. In most recent research few scholars used the K-means algorithm to pinpoint sub-urbanization patterns using SNPP-VIIRS data comparing finding with LandScan population, road network, and impervious extent data products [
37,
38,
39].
Due to the long coverage of DMSP/OLS satellite imagery, it is possible to identify multi-stage urban dynamics in both spatial and temporal dimensions. Thus, using these NTL data, the rapid, sluggish, saturated, and de-urbanization stages may be nicely observed. Due to the spectral and geographical complexity of the land cover within cities, fine-resolution remote sensing of urban areas is also difficult. Fortunately, night light sensors offer the unique ability to track human activities from a distance, in contrast to most conventional sensors [
26]. DMSP/OLS nighttime light composite data is one of the most efficient urban information sources. As well as, SNPP-VIIRS has a significant amount of enhancements in terms of spatial resolution, and calibrations [
40]. Artificial lights emitted from cities both human settlements and functions can be captured at night using DMSP/OLS and VIIRS. Since NTL data have a positive correlation with socioeconomic indicators and the density and height of built-up areas, it has been widely used to detect dynamics in urban land use, particularly in impervious surface and urban land use intensification. In particular, DMSP-OLS and SNPP-VIIRS data, are frequently used to extract urban entities, population, and gross domestic product, and analyze carbon emissions [
27,
41,
42,
43]. There are two advantages to using NTL images to detect urban entities. These data have the capability to map urban areas on large spatial-temporal scales and can precisely pinpoint where human activities are concentrated with low spatial resolution [
44]. According to scholars in the field of night light remote sensing, there are two main problems with urban mapping utilizing NTL data. Most studies only take into account impervious surfaces as urban entities. Thus, impervious surfaces were considered as the urban entities in the majority of urban entity mapping research [
4,
45,
46,
47,
48] although few studies recognized the relevance of urban entity extraction and the ability of NTL to extract urban entities. Though, DMSP-OLS and SNPP-VIIRS act as valuable resources for urban growth monitoring those studies lack prefecture-level cities in China. Thus, the applicability and validity of urban entities derived through NTL are to be analyzed further. Urban mapping in China using NTL is solely upon the short temporal scales based on DMSP and SNPP-VIIRS data except for a few studies [
48,
49]. But, the accurate, long-term SNPP-VIIRS-like data have the ability to identify, and extract urban growth at prefecture levels.
The study tried to extract and map urban entities in prefecture cities during 2000-2020 from an urban entity perspective using SNPP-VIIRS-Like data from the Harvard Dataverse. The research was organized under five sub-sections listed below. Whereas
Section 1 is devoted to explaining the research background and literature survey and
Section 2 describes the study area, materials, and methods. Results and findings were thoroughly explained in
Section 3 where compression was made between NTL data and other data products like LandScan and HE. In the discussion Section, it was presented the similarities and differences of the key findings with similar research works along with limitations.
Section 5 gives the conclusions and future research directions.
Figure 1.
(a) Provinces in PR of China; (b) Prefecture divisions; (c) selected prefecture divisions.
Figure 1.
(a) Provinces in PR of China; (b) Prefecture divisions; (c) selected prefecture divisions.
Figure 2.
Methodological flowchart of the study.
Figure 2.
Methodological flowchart of the study.
Figure 3.
Growth of UE, HE, and MODIS urban areas during 2000-2020.
Figure 3.
Growth of UE, HE, and MODIS urban areas during 2000-2020.
Figure 4.
Expansion of urban entities of capitals in northern, eastern, and southern provinces, 2000-2020.
Figure 4.
Expansion of urban entities of capitals in northern, eastern, and southern provinces, 2000-2020.
Figure 5.
Expansion of urban entities of capitals in central, northeast, northwest, and southwest provinces, 2000-2020.
Figure 5.
Expansion of urban entities of capitals in central, northeast, northwest, and southwest provinces, 2000-2020.
Figure 6.
Extracted urban entities (a) Shanghai; (b) Beijing; (c) Shenzhen; (d) Guangzhou with LandScan product and Landsat8 images, 2015.
Figure 6.
Extracted urban entities (a) Shanghai; (b) Beijing; (c) Shenzhen; (d) Guangzhou with LandScan product and Landsat8 images, 2015.
Figure 7.
Extracted urban entities of six provincial capitals with OSM road networks, 2015.
Figure 7.
Extracted urban entities of six provincial capitals with OSM road networks, 2015.
Figure 8.
Regression results of total urban areas of UE, HE, and MODIS from 2000-2020.
Figure 8.
Regression results of total urban areas of UE, HE, and MODIS from 2000-2020.
Figure 9.
Comparison results of urban areas for UE, HE, MODIS on provincial, and, national levels, 2000-2020.
Figure 9.
Comparison results of urban areas for UE, HE, MODIS on provincial, and, national levels, 2000-2020.
Figure 10.
Extracted urban entities of selected prefectures with HE, MODIS, and LandScan population product, 2015.
Figure 10.
Extracted urban entities of selected prefectures with HE, MODIS, and LandScan population product, 2015.
Figure 11.
Visual differences of urban entities of five capitals in Globeland, Landsat8, and SNPP-VIIRS-like images in 2020.
Figure 11.
Visual differences of urban entities of five capitals in Globeland, Landsat8, and SNPP-VIIRS-like images in 2020.
Figure 12.
Regression results of total urban areas for UE, HE, MODIS, and urban GDP 2000-2020.
Figure 12.
Regression results of total urban areas for UE, HE, MODIS, and urban GDP 2000-2020.
Table 1.
Spatial data sources of the study.
Table 1.
Spatial data sources of the study.
Data |
Year |
Format |
Resolution/scale |
Source |
SNPP-VIIRS-like |
2015-2020 |
Raster |
500 m |
https://dataverse.harvard.edu/dataset.xhtml |
NPP-VIIRS-like |
2000,2005, 2010 |
Raster |
500 m |
https://dataverse.harvard.edu/dataset.xhtml |
LandScan |
2015 |
Raster |
1000 m |
https://www.un-spider.org/links-and-resources/data-sources/landscan |
HE |
2015 |
Raster |
1000 m |
http://data.tpdc.ac.cn/zh-hans/data/3100de5c-ac8d-4091-9bbf-6a02de100c88/ |
MODIS |
2015 |
Raster |
500 m |
https://ladsweb.modaps.eosdis.nasa.gov/search/order/1/MCD12Q1--6 |
GlobeLand30 |
2020 |
Raster |
30 m |
http://www.globallandcover.com/defaults_en.html? |
OSM |
2015 |
Vector |
1:5000 |
https://www.openstreetmap.org |
LandSat8 |
2015 |
Raster |
30 m |
https://earthexplorer.usgs.gov/ |
Prefecture boundaries |
2019 |
Vector |
1:50,000,000 |
http://ngcc.sbsm.gov.cn/article/en/ |