Articles | Volume 13, issue 1
Earth Syst. Sci. Data, 13, 63–82, 2021
Earth Syst. Sci. Data, 13, 63–82, 2021

Data description paper 14 Jan 2021

Data description paper | 14 Jan 2021

A 30 m resolution dataset of China's urban impervious surface area and green space, 2000–2018

A 30 m resolution dataset of China's urban impervious surface area and green space, 2000–2018
Wenhui Kuang1, Shu Zhang1,2, Xiaoyong Li2,3, and Dengsheng Lu4,5 Wenhui Kuang et al.
  • 1Key Laboratory of Land Surface Pattern and Simulation, Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing 100101, China
  • 2College of Resources and Environment, University of Chinese Academy of Sciences, Beijing 10049, China
  • 3State Key Laboratory of Urban and Regional Ecology, Research Center for Eco-Environmental Sciences, Chinese Academy of Sciences, Beijing 100085, China
  • 4School of Geographical Sciences, Fujian Normal University, Fuzhou 350007, China
  • 5Fujian Provincial Key Laboratory of Subtropical Resources and Environment, Fujian Normal University, Fuzhou 350007, China

Correspondence: Wenhui Kuang (


Accurate and timely maps of urban underlying land properties at the national scale are of significance in improving habitat environment and achieving sustainable development goals. Urban impervious surface (UIS) and urban green space (UGS) are two core components for characterizing urban underlying environments. However, the UIS and UGS are often mosaicked in the urban landscape with complex structures and composites. The “hard classification” or binary single type cannot be used effectively to delineate spatially explicit urban land surface property. Although six mainstream datasets on global or national urban land use and land cover products with a 30 m spatial resolution have been developed, they only provide the binary pattern or dynamic of a single urban land type, which cannot effectively delineate the quantitative components or structure of intra-urban land cover. Here we propose a new mapping strategy to acquire the multitemporal and fractional information of the essential urban land cover types at a national scale through synergizing the advantage of both big data processing and human interpretation with the aid of geoknowledge. Firstly, the vector polygons of urban boundaries in 2000, 2005, 2010, 2015 and 2018 were extracted from China's Land Use/cover Dataset (CLUD) derived from Landsat images. Secondly, the national settlement and vegetation percentages were retrieved using a sub-pixel decomposition method through a random forest algorithm using the Google Earth Engine (GEE) platform. Finally, the products of China's UIS and UGS fractions (CLUD-Urban) at a 30 m resolution were developed in 2000, 2005, 2010, 2015 and 2018. We also compared our products with six existing mainstream datasets in terms of quality and accuracy. The assessment results showed that the CLUD-Urban product has higher accuracies in urban-boundary and urban-expansion detection than other products and in addition that the accurate UIS and UGS fractions were developed in each period. The overall accuracy of urban boundaries in 2000–2018 are over 92.65 %; and the correlation coefficient (R) and root mean square errors (RMSEs) of UIS and UGS fractions are 0.91 and 0.10 (UIS) and 0.89 and 0.11 (UGS), respectively. Our result indicates that 71 % of pixels of urban land were mosaicked by the UIS and UGS within cities in 2018; a single UIS classification may highly increase the mapping uncertainty. The high spatial heterogeneity of urban underlying covers was exhibited with average fractions of 68.21 % for UIS and 22.30 % for UGS in 2018 at a national scale. The UIS and UGS increased unprecedentedly with annual rates of 1605.56 and 627.78 km2 yr−1 in 2000–2018, driven by fast urbanization. The CLUD-Urban mapping can fill the knowledge gap in understanding impacts of the UIS and UGS patterns on ecosystem services and habitat environments and is valuable for detecting the hotspots of waterlogging and improving urban greening for planning and management practices. The datasets can be downloaded from (Kuang et al., 2020a).

1 Introduction

The effects of rapid urbanization on environments have been witnessed around the world (Seto et al., 2012; Bai et al., 2018; Kuang et al., 2020b, Zhang et al., 2021) and profoundly contribute to changes in biosphere, hydrosphere and atmosphere (Wu et al., 2014; Kuang et al., 2018). In China, a rapid urbanization process emerged in the 21st century (Xu and Min, 2013; Ma et al., 2014; Bai et al., 2014; Kuang, 2012; Kuang et al., 2016), resulting in a rapid increase in urban impervious surface area (UIS) (Kuang et al., 2013; Kuang and Dou, 2020; Lu et al., 2008, Kuang and Yan, 2018). This process further triggered various urban environmental problems such as urban heat islands and urban flooding (Haase et al., 2014; Hamdi and Schayes; 2007; Kuang, 2011; Kuang et al., 2015, 2017; Xu, 2006; Zhang et al., 2017). Although many green areas have been constructed in Chinese cities recently, China has a relatively lower percentage of urban green space (UGS) than developed countries such as the United States (Nowak and Greenfield, 2012; Kuang et al., 2014). These urban environmental problems triggered the urgency of developing accurate urban land cover datasets with high spatial resolutions for delineating the underlying urban environments. Along with the development of Earth observation technologies, remote sensing has become the mainstream method for mapping UIS and UGS and monitoring their changes (Weng, 2012; Wang et al., 2013; Lu et al., 2014, 2018; Zhang et al., 2009).

Various land use products such as the Global Land Cover product (GlobaLand30) (Chen et al., 2015), the University of Maryland (UMD) Land Cover Classification (Hansen et al., 2000), Moderate Resolution Imaging Spectroradiometer (MODIS)-based land use and land cover products (Friedl et al., 2010), GlobCover (Bontemps et al., 2011), and Finer Resolution Observation and Monitoring of Global Land Cover (FROM-GLC) (Gong et al., 2013) are freely available worldwide (Grekousis et al., 2015; Dong et al., 2017). These products have different definitions of urban areas or settlements due to their different classification systems, such as the International Geosphere-Biosphere Programme (IGBP) (Belward, 1996). Some urban land datasets, such as the Normalized Urban Areas Composite Index (NUACI), which were constructed by supervised learning approaches, have been released at a national or global scale with spatial resolutions from 30 m to 1 km (Liu et al., 2018; He et al., 2019; Gong et al., 2019). Others such as the built-up grid of the Global Human Settlement Layer (GHSL Built) (Pesaresi et al., 2013) and Global Urban Footprint (GUF) (Esch et al., 2017, 2018) have been published too. Most urban land products have focused on built-up land or urban area classification but cannot delineate urban land as a heterogeneous unit consisting of UIS, UGS and other fractions (Chen et al., 2015). Therefore, few urban land products have provided intra-urban UIS and UGS fractions at the sub-pixel level.

A detailed UIS dataset inside a city is required as a primary urban environmental index. Numerous studies on impervious surface mapping at the national scale mainly rely on medium–low-spatial-resolution remotely sensed data such as MODIS and the Defense Meteorological Satellite Program's Operational Linescan System (DMSP OLS) (Gong et al., 2013; Zhou et al., 2014; Grekousis et al., 2015; Zhou et al., 2015; Kuang et al., 2016; Zhou et al., 2018). Recently, more research has shifted to employ medium–high-spatial-resolution data (e.g., Landsat) to improve data products (Li et al., 2018; Liu et al., 2018; Gong et al., 2019, 2020a; Li et al., 2020; Lin et al., 2020). The US Geological Survey have developed the National Land Cover Database (NLCD) and provided impervious surface fraction, percent tree canopy, land cover classes and their changes with a spatial resolution of 30 m (Falcone and Homer, 2012; Yang et al., 2018). However, a detailed intra-urban UIS and UGS dataset with a 30 m spatial resolution for China at the national scale is not available yet, making it difficult to conduct detailed analysis of aspects such as urban living environments.

A systematic assessment on urban land mapping algorithms indicates that previous research mainly classified urban land into a single type with “urban area” or impervious surface area (ISA), which limits recognition of the urban environment (Reba and Seto, 2020). There are two critical challenges in mapping urban land cover composites. Firstly, the conceptual definition of urban land or ISA in previous research is unclear; thus, the spatial extent is inconsistent, resulting in a large divergence in the statistical area of urban land. Meanwhile, the segmentation on urban–rural boundaries has not been accurate from moderate-resolution satellite images using computer-based automatic classification owing to differences in geographic conditions, social economic conditions and land policies. Therefore, accurate mapping of urban–rural boundaries is pivotal in detecting urban land cover change. Secondly, the spatial heterogeneity of urban surface properties has resulted in difficulty in decomposing urban land cover types with complex surface materials at the pixel scale, which has been limited by the huge amounts of data processing and storage capacities required for a 30 m resolution.

In reality, the urban land cover is composed of UIS, UGS and other fractions. UIS refers to the urban impervious surface features caused by artificial land use activities, like building roofs, asphalt or cement roads, and parking lots. UGS is an important component of the green infrastructure of cities and provides a range of ecosystem services, including parks, trees and grass. Previous studies have proven that spectral mixture analysis (SMA) provides an effective tool to retrieve the UIS and UGS fractions from Landsat multispectral imagery (Lu and Weng, 2004, 2006; Peng et al., 2016; Kuang et al., 2018). However, this method needs local knowledge for problem-specific analysis such as intra-urban land cover analysis of a single city or a single urban agglomeration (Zhang and Weng, 2016; Xu et al., 2018). Although the globally standardized SMA can effectively extract substrate, dark areas and vegetation (Small and Milesi, 2013), the UIS cannot be accurately and directly extracted from multispectral images without post-processing considering its widely spectral variation and different meanings of UIS and substrate (Lu et al., 2014). Because of the high correlation between UIS and vegetation indices in the urban landscape (Weng et al., 2004), a fractional UIS dataset can be estimated from vegetation indices using a regression-based approach (Sexton et al., 2013; Wang et al., 2017).

To address the above issues, we propose a synthetical strategy to utilize the advantage of both accurate urban boundary information from China's Land Use/cover Dataset (CLUD) extracted by human–computer digitalization and the retrieval of UIS and UGS fractions through big-data processing from the GEE platform. Based on this strategy, we developed the product of a national UIS and UGS fraction dataset at a 30 m spatial resolution in 2000, 2005, 2010, 2015 and 2018 across China. This dataset provides a foundation for understanding urban dwellers' environments and enhances our understanding of the impacts of urbanization on ecological services and functions, and it will also be helpful in future research and practices in urban planning and urban environmental sustainability.

2 The strategy of developing the CLUD-Urban product

To acquire the accurate CLUD-Urban product, three steps were generally implemented according to our mapping strategy. Firstly, national urban boundaries in 2000–2018 were extracted from CLUD which was generated using a uniform technological flow and classification system in a human–computer digitalization environment. Time series of urban boundaries and their expansions have good performance in accuracy and data quality. The national urban vector boundaries in 2000, 2005, 2010, 2015 and 2018 were converted to raster data with a 30 m resolution for further processing (Fig. 1). Secondly, the settlement and vegetation fractions with a 30 m resolution were retrieved using a random forest algorithm in the GEE platform. Thirdly, the UIS and UGS fractions with a 30 m resolution were mapped through overlaying the urban boundaries of CLUD with settlement and vegetation fractions, respectively (Fig. 1). Accuracy assessment both of urban boundaries and of UIS and UGS fractions was implemented using samples from Google Earth images. Quality control was conducted throughout the data processing in mapping the CLUD-Urban product. A detailed description is given in the following sections.

Figure 1The technological flowchart of generating the CLUD-Urban product.


3 Data sources and pre-processing

Landsat is the longest-running satellite series for Earth observation. Landsat Thematic Mapper (TM), Enhanced Thematic Mapper Plus (ETM+) and Operational Land Imager (OLI) data with path ranges of 118–149 and row ranges of 23–43 in China were selected (Table 1). In mapping CLUD, Landsat TM, ETM+ and OLI in each period, China–Brazil Earth Resources Satellite program (CBERS) and Huan Jing (HJ-1A and HJ-1B) satellite images in 2010 were used to generate false-color composite images with near-infrared, red and green spectral bands as red, green and blue. Image enhancement was performed to improve the visual interpretation quality. Image-to-image registration was conducted to control the image rectification errors of less than 2 pixels (60 m). CBERS-1 and Huan Jing (HJ-1A and HJ-1B) satellite images were only used in extracting the vector polygons of CLUD in 2010, which was conducted using uniform data processing with Landsat images.

Table 1The multitemporal data series used in this research.

Download Print Version | Download XLSX

In the retrieval of settlement and vegetation fractions, Landsat TM, ETM+, and OLI data in each period from January to December were collected using the GEE platform. Shuttle Radar Topography Mission (SRTM) digital elevation model data and the normalized difference vegetation index (NDVI) with a 30 m resolution were acquired as input parameters to retrieve settlement and vegetation fractions. Google Earth images in selected cities with a 0.6 m resolution were used to assess the accuracy of the CLUD-Urban product.

4 Extraction of urban boundaries from CLUD

4.1 The classification system and interpretation symbols

CLUD with 30 m resolution was developed by the Chinese Academy of Sciences and has been updated from 2000 to 2018 every 5 or 3 years. This dataset can delineate land use or land cover change associated with human activities, including urbanization at a scale of 1:100 000 (Liu et al., 2005a, b, 2010). This product adopted a hierarchical classification system covering the 6 first-level classes and the 25 second-level classes. Here the 6 first-level classes comprise cropland, woodland, grassland, water body, construction land and unused land. A detailed description of each class can be found in previous publications (Liu et al., 2005b; Zhang et al., 2014). The construction land was divided into three second-level classes, including urban land, rural settlements, and industrial and mining lands beyond cities. Urban land was defined as a built-up area of the concentrated construction, i.e., buildings, roads, squares, green infrastructure and other lands for providing a living, industrial production, and ecosystem services for the dwellers of cities or towns. According to this definition, urban land can be megacities (more than 10 million population), megalopolises (5–10 million population), large cities (1–5 million population), medium cities (0.5–1 million population), small cities (0.2–0.5 million population) and towns (less than 0.2 million population) (Kuang, 2020a). The industrial and traffic lands outside cities are excluded in the urban land. Based on the designed classification system, the interpretation symbols from the second-level classes were built for the false-color composite images as a reference to aid the human–computer interpretation (Fig. 2) (Zhang et al., 2014).

Figure 2The interpretation symbols and extracted urban boundaries from Landsat images in Beijing city. The images were provided by the Geospatial Data Cloud site, Computer Network Information Center, Chinese Academy of Sciences (, last access: 27 September 2017).

4.2 Land use and dynamic polygon interpretation

According to the CLUD classification system, the vector polygons of land use classes in 2000 were digitalized through overlying the false-color composite images with the aid of interpretation symbols and geoknowledge from each zone (Fig. 3). The polygons of urban lands were identified through using detailed image interpretation symbols for each second-level land use class based on Landsat or similar-resolution images. Usually, the polygons of urban lands exhibit larger sizes than rural settlements and others (e.g., industrial and traffic lands) in a cinerous color ornamented with white. Digitalization personnel differentiated the urban land from rural settlements and others based on interpretation symbols and geoknowledge from field investigation (Fig. 2). In the digitalization environment, each vector polygon was assigned a code of the second-level classes. The vector polygons of land use classes in 2000 were double-checked to ensure the correct type in interpretation. The dynamic polygons were extracted through comparing the difference in two differently dated images and assigned the codes including the types before and after changes (Fig. 3). The land use changes within 5 or 3 years were detected using the uniform method. Finally, the land use maps in 2000, 2005, 2010, 2015 and 2018 and their changes at 5- or 3-year intervals were generated for CLUD. The detailed technological flow can be found in previous publications (Liu et al., 2005b; Zhang et al., 2014). An example of a land use map in 2010 in the Conghua district of Guangzhou and dynamic land use changes in 2010–2015 is illustrated in Fig. 3.

Figure 3Land use classification and extracted vector polygons as an example with the Conghua district of Guangzhou.


4.3 Retrieval of multitemporal urban boundaries

The vector boundaries of urban extents were extracted from the CLUD land use maps in each period (Kuang et al., 2016). We also examined 10 732 urban vector polygons in 2000. The number of polygons increase to 50 061 in 2018. The urban vector boundaries were acquired from Landsat images or similar-resolution images. The vector polygons of urban boundaries were converted to raster data with a 30 m×30 m cell size. The datasets on urban land across China in 2000, 2005, 2010, 2015 and 2018 were generated with 30 m resolution. Here we showed urban boundaries and the expansion process with a 30 m resolution in the cities of Xi'an, Wuhan, Guangzhou and Ürümqi (Fig. 4).

Figure 4The urban boundaries extracted from CLUD with a 30 m resolution in selected cities. The administrative boundaries were provided by the National Geomatics Center of China (, last access: 15 September 2020); DEM dataset was downloaded from SRTM 90 m digital elevation data (, last access: 15 September 2020).

5 Mapping UIS and UGS fractions using GEE platform

5.1 Collection of training samples

The training samples of UIS and UGS fractions are a pivotal input parameter in the random forest model for mapping national settlement and vegetation fractions. In light of large discrepancies among UIS and UGS composites in different climate zones with various geographical and social economic conditions, we collected a total of 2570 samples from randomly selected cities in different climate zones (Schneider et al., 2010) (Fig. 5). Here we also refer to the existing UIS dataset to acquire samples with 10 % intervals of the ISA fraction, and those samples are primarily distributed in the homogeneous UIS or UGS areas, which might provide more effective samples and decrease the impact of imagery mismatch. The samples of UIS and UGS covered with diversified types, including buildings, roads and squares and grass and trees from parks, roads and residential green spaces. The UIS and UGS percentages were interpreted within each sample using Google Earth images (Fig. 5b1–b4). Finally, the training samples in 2000, 2005, 2010, 2015 and 2018 were used for training the random forest model.

Figure 5Distribution of sampling cities in China and training samples in selected cities. The images were provided by the Geospatial Data Cloud site, Computer Network Information Center, Chinese Academy of Sciences (, last access: 13 August 2015 for Shijiazhuang, 1 July 2015 for Xi'an, 2 April 2015 for Chengdu, and 20 November 2015 for Kunming). The administrative boundaries were provided by the National Geomatics Center of China (, last access: 15 September 2020).

5.2 Retrieval of settlement and vegetation fractions using random forest model

Many previous studies have indicated that random forest is more effective and accurate in classifying urban land types than other machine learning approaches such as support vector machines (SVMs) and artificial neural networks (ANNs) (Zhang et al., 2020). Random forest exhibits a strong capacity in processing high-dimensional datasets and has been successfully applied to mapping global ISA at a 30 m resolution (Zhang et al., 2020). In this research, we proposed a strategy to acquire the settlement and vegetation percentage at the pixel scale using the advantage of random forest and big-data processing based on the GEE platform.

According to 16 global urban ecoregions based on temperature, precipitation, topographic conditions and social economic factors (Schneider et al., 2010), China has three urban ecoregions. In each urban ecoregion, the annual maximum NDVI; spectral bands in Landsat TM, ETM+ and OLI; and the slope index derived from the SRTM DEM with a 30 m resolution were selected as the input parameters to run the random forest model. The Landsat images were from 1 January to 31 December of each baseline year. The annual maximum NDVI (NDVImax) was retrieved using Eq. (1):

(1) NDVI max = max ( NDVI 1 , NDVI 2 , , NDVI i ) ,

where NDVIi is the NDVI value of the ith image. The individual NDVI was calculated from Landsat images in the period between 1 January to 31 December, and all images were collected using GEE (Gorelick et al., 2017).

In the GEE platform, the settlement and vegetation fractions were calculated for each urban ecoregion through using the training parametrizations. The lawn, forest or their mosaicked areas were selected as input samples in mapping UGS. Post-processing was implemented to remove the pixels with NDVI values greater than 0.5 or DEM slope values greater than 15. In arid and semi-arid areas, the enhanced bare-soil index (EBSI) was utilized to separate UIS from bare soils (As-syakur et al., 2012; Li et al., 2019). As a result, the settlement and vegetation fractions with 30 m×30 m in 2000, 2005, 2010, 2015 and 2018 were generated for developing the CLUD-Urban product (Fig. 6).

Figure 6Distribution of sampling cities in China and training samples in selected cities. The administrative boundaries and residential point information were provided by the National Geomatics Center of China (, last access: 15 September 2020).

5.3 Mapping of UIS and UGS fractions

The settlement and vegetation fractions with a 1× 1 grid of each period were downloaded from GEE platform. In ArcGIS 10.0 software, the settlement and vegetation layers were merged at a provincial scale with a 30 m×30 m. The national UIS and UGS fractions with 30 m×30 m resolution in 2000, 2005, 2010, 2015 and 2018 were produced through overlaying the urban boundaries of CLUD with settlement and vegetation fractions, respectively (Figs. 7, 8 and 9).

Figure 7Spatial distribution of urban impervious surface (UIS) in 2000–2018 across China. The administrative boundaries were provided by the National Geomatics Center of China (, last access: 15 September 2020).

Figure 8Spatial distribution of urban green space (UGS) in 2000–2018 across China. The administrative boundaries were provided by the National Geomatics Center of China (, last access: 15 September 2020).

Figure 9The spatial distribution of urban impervious surface (UIS) in selected cities from 2000 to 2018. DEM dataset was downloaded from SRTM 90 m digital elevation data (, last access: 15 September 2020).

6 Accuracy assessment of the CLUD-Urban product

The national urban boundaries and UIS and UGS fractions were assessed through qualitative and quantitative indices, respectively. Firstly, on the accuracy of CLUD in 2000, 2005 and 2010 we referred to our previous publications (Liu et al., 2010, 2014; Zhang et al., 2014). The accuracy of the six first-level classes – cropland, forest, grassland, built-up area, water body and unused – and of the second-level land use and land cover types, including urban land, rural settlements, industrial and traffic lands, was assessed using the field investigation data and the Google Earth images (Liu et al., 2010, 2014; Zhang et al., 2014). We also implemented accuracy assessment on urban boundaries of CLUD from 2000 to 2018 using overall accuracy, producer's accuracy and user's accuracy (Fig. 10) (Kuang et al., 2016; Kuang, 2020a).

Figure 10Validation samples on the CLUD-Urban product. The administrative boundaries were provided by the National Geomatics Center of China (, last access: 15 September 2020).

The validation samples for assessing the accuracy of UIS and UGS fractions were collected within urban boundaries using a stratified random sampling method with the ISA fraction at 10 % intervals. Those samples with a window size of 3 pixels × 3 pixels (90 m×90 m) were used to digitalize the UIS and UGS polygons through the human–computer interaction based on Google Earth images (Kuang et al., 2014; Kuang, 2020b). A total of 1869 validation samples were randomly acquired in different regions in China in 2000–2018, including 1070 samples located in the changed UIS and UGS areas (Fig. 10). Mean UIS and UGS fractions in each grid were calculated. The comparison between estimated values and validation values was conducted using the correlation coefficient (R) and root mean square error (RMSE) (Kuang et al., 2014; Kuang, 2020b). We also evaluated the changed UIS and UGS areas using R and RMSE based on 1070 validation samples.

7 Results

7.1 The accuracy of CLUD-Urban

A quality check and data integration were performed for the years 2000, 2005, 2010, 2015 and 2018 to ensure the quality and consistency of the interpretation results. Our assessment results indicated the overall accuracy of the first-level land use and land cover types is 98.04 % in 2000, 94.3 % in 2010, 91.64 % in 2015 and 91.12 % in 2018 (Liu et al., 2014; Zhang et al., 2014; Kuang et al., 2016; Ning et al., 2018). The built-up area has the highest accuracy among the six land use types owing to its clear urban boundaries, and the accuracy reached 98.92 % in 2000 and 97.01 % in 2005 according to previous assessment (Zhang et al., 2014). The user's accuracy of urban land type is relatively high with 93.67 % in 2010, 92.65 % in 2015 and 91.32 % in 2018 (Table 2). Overall, the urban land accuracy shows a decreasing trend, which resulted from the fuzzy and unidentifiable urban–rural boundaries owing to the continuous pattern of urban–rural land driven by China's fast urban development since the 21st century. In CLUD, the change polygons were identified based on human interpretation. The validation of UIS and UGS fractions in each period showed that the RMSEs were 0.11–0.12 and 0.11–0.12 respectively, and the R values were 0.89–0.91 and 0.87–0.90, respectively (Table 3). The R and RMSE values for the changed UIS areas in 2000–2018 are 0.88 and 0.12, and those for the changed UGS areas in the same period are 0.85 and 0.12, respectively.

Table 2Confusion matrix of China's Land Use/cover Dataset.

Download Print Version | Download XLSX

Table 3Accuracy assessments for the UIS and UGS fractions.

Download Print Version | Download XLSX

7.2 Patterns and dynamics of UIS and UGS since the beginning of the 21st century

Our result indicated that China's UIS increased from 2.46×104km2 in 2000 to 5.35×104km2 in 2018 (Fig. 7). From the perspective of the quality of dwellers' habitat environments, the percentage of UIS in China's urban areas in 2018 is 74.42 %, showing a higher UIS density than developed countries like the USA (Kuang et al., 2014). However, the UIS percentage in urban areas decreased from 74.42 % in 2000 to 68.21 % in 2018 owing to the improvement of urban greening conditions. As shown in Fig. 7, the UIS across China is mainly distributed in the coastal and central regions and is relatively discrete in the western regions. The pattern of “high in east and low in west” national UIS remained unchanged between 2000 and 2018 (Fig. 7). China's UGS shows an increasing trend. The total UGS area increased from 1.00×104km2 in 2000 to 1.83×104km2 in 2018 (Fig. 8). Looking at both UIS and UGS in urban areas, our results indicate a slight increase in UGS density and decrease in UIS density, which has resulted from strengthening urban greening since the start of the 21st century. The UGS percentage rose from 18.91 % in 2000 to 22.30 % in 2018. As shown in Fig. 9, UIS and UGS of cities from coastal, northeastern and southwestern China have high spatial heterogeneity and showed the different urban expansion rate in the past 28 years.

Figure 11The fractions of urban impervious surface (UIS) and urban green space (UGS) at national and regional scales (coastal, central, western and northeastern zones) in 2000 and 2018.


The large discrepancies of the UIS and UGS percentage in urban areas were exhibited among eastern, central, western and coastal zones. The coastal zone showed a remarkable increasing trend from 16.50 % in 2000 to 21.66 % in 2018 (Figs. 9 and 11). We also found that urban greening conditions were positively improved in Beijing in the same period, which resulted in the increase in UGS percentage and decrease in UIS percentage in urban areas (Fig. 9). This means that urban habitat environment in coastal zones has become more liveable and comfortable, which is associated with the greening of parks and forests. We also found that the western cities have a relatively low UGS percentage in urban areas, which is 0.86 % lower than the average of China owing to the low-greening conditions (Figs. 9 and 11).

Table 4A summary of existing urban land products.

Download Print Version | Download XLSX

7.3 Comparisons of the CLUD-Urban product with other datasets

We compared the vector boundaries of urban areas with the existing land use products and found obvious discrepancies because of the differences in data production, data source, resolution and definition of urban land use types. The spatial resolutions of land cover products range from 30 to 1000 m. Figure 12 provides a comparison of urban land datasets (see Table 4 for these datasets), showing that our product has better performance in delineating the detailed spatial patterns of intra-urban land cover, i.e. the composite of UIS and UGS (note both the GHSL Built and GlobaLand30 products cover only 2 years). The accuracy of urban boundaries from CLUD-Urban is over 92 % and is basically inconsistent with that of the impervious surface map (Zhang et al., 2020). Our dataset has a higher classification accuracy in urban boundaries than that of GHSL with 90.3 %, FROM-GLC with 89.6 %, Human Built-up and Settlement Extent (HBASE) with 88.0 %, GlobaLand30 with 88.4 % and NUACI with 85.6 %. Furthermore, our CLUD-Urban product can accurately delineate the spatial heterogeneity of UIS and UGS composites, which showed the R with 0.90 and 0.89 and RMSE with 0.11 and 0.11, respectively. In those existing datasets, the UIS and UGS composites can't be effectively decomposed at the pixel scale (Fig. 12).

Figure 12A comparison of urban land cover between this product and other datasets in Beijing. The Landsat images were provided by the Geospatial Data Cloud site, Computer Network Information Center, Chinese Academy of Sciences (, last access: 20 August 2000 for 2000, 8 September 2005 for 2005, 8 June 2010 for 2010, 23 September 2015 for 2015, and 8 April 2018 for 2018).

8 Discussions

8.1 The mapping advantages integrated with human–computer interpretation and GEE platform

In mapping urban land use and land cover change at national scale, two pivotal steps were required to segment the urban land, rural settlements, and industrial and traffic lands in periphery of cities for accurately acquiring the urban boundaries and to retrieve the UIS and UGS fractions at pixel scale. The urban boundaries are generally mapped using classification methods such as unsupervised classifiers, supervised classifiers, human–computer interpretation and other advanced algorithms (i.e., ANN, SVM and random forest) (Wu and Murray, 2003; Zhang et al., 2020). Among these methods, human–computer interpretation is widely regarded as a most accurate method in classifying urban land use and land cover changes, especially in both detecting changing information and extracting vector polygons as whole geofeatures. However, this method takes more time and manual labor to digitalize a large number of polygons. CLUD has an advantage for providing the accurate urban boundaries and is updated at an interval of every 5 or 3 years from 2000 to 2018.

Cities or towns were classified as a homogeneous unit in CLUD. We developed the UIS and UGS fractions to fill the data gap for the requirement of urban environmental management. Here we adopted the advantage of high accuracy and long-time series in mapping urban land from CLUD. Meanwhile we also utilized the highly efficient computation and large storage capacities on the GEE platform. In mapping the CLUD-Urban product, we proposed quantitatively retrieving the UIS and UGS fractions using random forest. Because we used advantages of manual interpretation and intelligent computation, CLUD-Urban exhibits high accuracy and reliability in delineating urban land surface properties.

8.2 The potential implications of promoting habitat environment and urban sustainability

The CLUD-Urban product may effectively delineate the “built-up environment” of Chinese cities, especially maps on surface imperviousness and greening conditions (Kuang, 2020b). CLUD-Urban can be applied to such fields as enhancing the quality of the urban habitat environment, reducing urban heat islands, and improving prevention of rainstorm and flood disasters (Huang et al., 2018). Our previous study indicated that the thermal dissipation strength of forest canopy or lawns in cities may be assessed at the pixel scale and that the greening projects are more effective in alleviating urban-heat-island intensity (Kuang et al., 2015). The CLUD-Urban product also helps identify urban flood regulation priority areas based on ecosystem service approaches (Li et al., 2020).

The analysis of CLUD-Urban indicates an unprecedented rate and magnitude of urban expansion since the start of the 21st century. The low UGS of cities in western zones indicates the need to promote their greening level (Kuang and Dou, 2020). The CLUD-Urban product can also be used to assess sustainable-development-goal (SDG) targets such as the ratio of land consumption to population growth or average share of the built-up area that is open space for public use. Therefore, CLUD-Urban can have many potential applications in the development of sustainable, liveable and resilient cities.

8.3 Limitations of the method and dataset

Although state-of-the-art technologies and methodologies were applied to the development of CLUD-Urban (Dong et al., 2017; Kuang et al., 2020b), improvement to the mapping CLUD-Urban quality can still be made. For example, the retrieval of UIS and UGS was conducted as a prerequisite of CLUD, which focused on the pixel decomposition of UIS and UGS in urban areas. If the UIS and UGS fractions are parameterized to be input into a hydrological process model or urban climate, the settlements or impervious surfaces located on the outskirts of a city or in rural areas are removed from CLUD. To address this issue, the first-level classification or second-level classification of CLUD should be utilized to merge with UIS and UGS using the method in our previous publication (Kuang et al., 2020b). Mapping CLUD requires a large amount of labor and time as many interpreters are involved in this work. The extraction of urban boundaries might be subjective, and there is a time lag in mapping UIS and UGS. It is needed to develop some advanced tools to extract urban boundaries using automatic algorithms.

Fine urban land use and land cover change mapping at a national scale with high-resolution multi-source data may be developed with the aid of big data and cloud platforms (Gong et al., 2020a). The development of a series of new algorithms and models are pivotal for improving the accuracy of datasets in retrieving urban boundaries and land cover composites. However, geoknowledge is still essential for retrieving a high-quality dataset (Kuang et al., 2018). CLUD-Urban can contribute to the development of sustainable cities, such as with Global Ecosystems and Environment Observation (GEO) and UN-Habitat, in future.

9 Data availability

All data presented in this paper are available at (Kuang et al., 2020a). These new-version datasets include the UIS and UGS fractions with a 30 m spatial resolution in 2000, 2005, 2010 2015 and 2018. A detailed metadata description is provided, including contact information.

10 Conclusion

CLUD-Urban – China's UIS and UGS fraction datasets with 30 m spatial resolution – was generated using multiple data sources. CLUD-Urban provides detailed delineation of UIS and UGS components in China for the years of 2000, 2005, 2010, 2015 and 2018. Comparing to other products, the novelty of this dataset is to take cities as heterogeneous units at the pixel level, which consist of UIS, UGS and others. The accuracy of the CLUD-Urban dataset is over 92.65 % using the integrated approach of visual interpretation and prior knowledge. The RMSEs of UIS and UGS fractions are 0.10 and 0.14, respectively. Results from the analysis of urban areas, including UIS and UGS, show large regional differences in China. CLUD-Urban provides fundamental data sources for examining urban environment issues and for delineating intra-urban structure or urban landscape at the national scale.

Author contributions

WK, SZ and XL designed the research; SZ and XL implemented the research; WK, SZ and DL wrote the paper.

Competing interests

The authors declare that they have no conflict of interest.


We thank Yali Hou and Changqing Guo for processing the data and Fengyun Sun and Rafiq Hamdi for their help in manuscript editing.

Financial support

This research has been supported by the National Natural Science Foundation of China (grant no. 41871343) and the Strategic Priority Research Program A of the Chinese Academy of Sciences (grant no. XDA23100201).

Review statement

This paper was edited by Yuyu Zhou and reviewed by four anonymous referees.


As-syakur, A. R., Adnyana, I. W. S., Arthana, I. W., and Nuarsa, I. W.: Enhanced built-up and bareness index (EBBI) for mapping built-up and bare land in an urban area, Remote Sens., 4, 2957–2970,, 2012. 

Bai, X., Shi, P., and Liu, Y.: Society: realizing China's urban dream, Nature, 509, 158–160,, 2014. 

Bai, X., Dawson, R. J., Urge-Vorsatz, D., Delgado, G. C., Barau, A. S., Dhakal, S., Dodman, D., Leonardsen, L., Masson-Delmotte, V., Roberts, D., and Schultz, S.: Six research priorities for cities and climate change, Nature, 555, 19–21,, 2018. 

Belward, A. (Ed.): The IGBP-DIS global 1 km land cover data set “DISCover”: proposal and implementation plans. Report of the Land Cover Working Group of the IGBP-DIS, IGBP-DIS Working Paper, No. 13, Stockholm, 1996. 

Bontemps, S., Defourny, P., Bogaert, E., Arino, O., Kalogirou, V., and Perez, J.: Globcover 2009. Products Description and Validation Reports, available at: (last access: 18 April 2019), 2011. 

Chen, J., Chen, J., Liao, A., Cao, X., Chen, L., Chen, X., He, C., Han, G., Peng, S., Lu, M., Zhang, W., Tong, X., and Mills, J.: Global land cover mapping at 30 m resolution: a POK-based operational approach, ISPRS J. Photogramm., 103, 7–27,, 2015. 

Dong, J., Kuang, W., and Liu, J.: Continuous land cover change monitoring in the remote sensing big data era, Sci. China Earth Sci., 60, 2223–2224,, 2017. 

Esch, T., Heldens, W., Hirner, A., Keil, M., Marconcini, M., Roth, A., Zeidler, J., Dech, S., and Strano, E.: Breaking new ground in mapping human settlements from space – the Global Urban Footprint, Isprs J. Photogramm., 134, 30–42,, 2017. 

Esch, T., Bachofer, F., Heldens, W., Hirner, A., Marconcini, M., Palacios-Lopez, D., Roth, A., Üreyen, S., Zeidler, J., Dech, S., and Gorelick, N.: Where we live – a summary of the achievements and planned evolution of the global urban footprint, Remote Sens., 10, 895,, 2018. 

Falcone, J. A. and Homer, C. G.: Generation of a U.S. national urban land-use product, Photogramm. Eng. Rem. S., 78, 1057–1068,, 2012. 

Friedl, M. A., Sulla-Menashe, D., Tan, B., Schneider, A., Ramankutty, N., Sibley, A., and Huang, X.: MODIS Collection 5 global land cover: algorithm refinements and characterization of new datasets, Remote Sens. Environ., 114, 168–182,, 2010. 

Gong, P., Wang, J., Yu, L., Zhao, Y., Zhao, Y., Liang, L., Niu, Z., Huang, X., Fu, H., Liu, S., Li, C., Li, X., Fu, W., Liu, C., Xu, Y., Wang, X., Cheng, Q., Hu, L., Yao, W., Zhang, H., Zhu, P., Zhao, Z., Zhang, H., Zheng, Y., Ji, L., Zhang, Y., Chen, H., Yan, A., Guo, J., Yu, L., Wang, L., Liu, X., Shi, T., Zhu, M., Chen, Y., Yang, G., Tang, P., Xu, B., Giri, C., Clinton, N., Zhu, Z., Chen, J., and Chen, J.: Finer resolution observation and monitoring of global land cover: first mapping results with Landsat TM and ETM+ data, Int. J. Remote Sens., 34, 2607–2654,, 2013. 

Gong, P., Li, X., and Zhang, W.: 40-Year (1978–2017) human settlement changes in China reflected by impervious surfaces from satellite remote sensing, Sci. Bull., 64, 756–763,, 2019. 

Gong, P., Chen, B., Li, X., Liu, H., Wang, J., Bai, Y., Chen, J., Chen, X., Fang, L., and Feng, S.: Mapping essential urban land use categories in China (EULUC-China): preliminary results for 2018, Sci. Bull., 65, 182–187,, 2020a. 

Gong, P., Li, X., Wang, J., Bai, Y., Chen, B., Hu, T., Liu, X., Xu, B., Yang, J., Zhang, W., and Zhou, Y.: Annual maps of global artificial impervious area (GAIA) between 1985 and 2018, Remote Sens. Environ., 236, 111510,, 2020b. 

Gorelick, N., Hancher, M., Dixon, M., Ilyushchenko, S., Thau, D., and Moore, R.: Google Earth Engine: planetary-scale geospatial analysis for everyone, Remote Sens. Environ., 202, 18–27,, 2017. 

Grekousis, G., Mountrakis, G., and Kavouras, M.: An overview of 21 global and 43 regional land-cover mapping products, Int. J. Remote Sens., 36, 5309–5335,, 2015. 

Haase, D., Larondelle, N., Andersson, E., Artmann, M., Borgström, S., Breuste, J., Gomez-Baggethun, E., Gren, Å., Hamstead, Z., Hansen, R., Kabisch, N., Kremer, P., Langemeyer, J., Rall, E. L., McPhearson, T., Pauleit, S., Qureshi, S., Schwarz, N., Voigt, A., Wurster, D., and Elmqvist, T.: A quantitative review of urban ecosystem service assessments: concepts, models, and implementation, Ambio, 43, 413–433,, 2014. 

Hamdi, R. and Schayes, G.: Sensitivity study of the urban heat island intensity to urban characteristics, Int. J. Climatol., 28, 973–982,, 2007. 

Hansen, M. C., Defries, R. S., Townshend, J. R. G., and Sohlberg, R.: Global land cover classification at 1 km spatial resolution using a classification tree approach, Int. J. Remote Sens., 21, 1331–1364,, 2000. 

He, C., Liu, Z., Gou, S., Zhang, Q., Zhang, J., and Xu, L.: Detecting global urban expansion over the last three decades using a fully convolutional network, Environ. Res. Lett., 14, 34008,, 2019. 

Huang, C., Yang, J., and Jiang, P.: Assessing impacts of urban form on landscape structure of urban green spaces in China using Landsat images based on Google Earth Engine, Environ. Res. Lett., 10, 054011,, 2018. 

Kuang, W.: Simulating dynamic urban expansion at regional scale in Beijing-Tianjin-Tangshan Metropolitan Area, J. Geogr. Sci., 21, 317–330,, 2011. 

Kuang, W.: Evaluating impervious surface growth and its impacts on water environment in Beijing-Tianjin-Tangshan metropolitan area, J. Geogr. Sci., 22, 535–547,, 2012. 

Kuang, W.: 70 years of urban expansion across China: Trajectory, pattern, and national policies. Sci. Bull., 65, 1970–1974,, 2020a. 

Kuang, W.: National urban land-use/cover change since the beginning of the 21st century and its policy implications in China, Land Use Pol., 97, 104747,, 2020b. 

Kuang, W. and Dou, Y.: Investigating the patterns and dynamics of urban green space in China's 70 major cities using satellite remote sensing, Remote Sens., 12, 1929,, 2020. 

Kuang, W. and Yan, F.: Urban structural evolution over a century in Changchun city, Northeast China, J. Geogr. Sci., 28, 1877–1895,, 2018. 

Kuang, W., Liu, J., Zhang, Z., Lu, D., and Xiang, B.: Spatiotemporal dynamics of impervious surface areas across China during the early 21st century, Chin. Sci. Bull., 58, 1691–1701,, 2013. 

Kuang, W., Chi, W., Lu, D., and Dou, Y.: A comparative analysis of megacity expansions in China and the U.S.: patterns, rates and driving forces, Landscape Urban Plan., 132, 121–135,, 2014. 

Kuang, W., Dou, Y., Zhang, C., Chi, W., Liu, A., Liu, Y., Zhang, R., and Liu, J.: Quantifying the heat flux regulation of metropolitan land use/land cover components by coupling remote sensing modeling with in situ measurement, J. Geophys. Res.-Atmos., 120, 113–130,, 2015. 

Kuang, W., Liu, J., Dong, J., Chi, W., and Zhang, C.: The rapid and massive urban and industrial land expansions in China between 1990 and 2010: a CLUD-based analysis of their trajectories, patterns, and drivers, Landscape Urban Plan., 145, 21–33,, 2016. 

Kuang, W., Yang, T., Liu, A., Zhang, C., Lu, D., and Chi, W.: An EcoCity model for regulating urban land cover structure and thermal environment: taking Beijing as an example, Sci. China Ser. D-Earth Sci., 60, 1098–1109,, 2017. 

Kuang, W., Yang, T., and Yan, F.: Examining urban land-cover characteristics and ecological regulation during the construction of Xiong'an New District, Hebei Province, China, J. Geogr. Sci., 28, 109–123,, 2018. 

Kuang, W., Zhang, S., Li, X., and Lu, D.: A 30-meter resolution dataset of impervious surface area and green space fractions of China's cities, 2000–2018, Zenodo,, 2020a. 

Kuang, W., Du, G., and Lu, D.: Global observation of urban expansion and land-cover dynamics using satellite big-data, Sci. Bull., accepted,, 2020b. 

Li, H., Wang, C., Zhong, C., Su, A., Xiong, C., Wang, J., and Liu, J.: Mapping urban bare land automatically from Landsat imagery with a simple index, Remote Sens., 9, 249,, 2019. 

Li, X., Zhou, Y., Zhu, Z., Liang, L., Yu, B., and Cao, W.: Mapping annual urban dynamics (1985–2015) using time series of Landsat data, Remote Sens. Environ., 216, 674–683,, 2018. 

Li, X., Zhou, Y., Zhu, Z., and Cao, W.: A national dataset of 30 m annual urban extent dynamics (1985–2015) in the conterminous United States, Earth Syst. Sci. Data, 12, 357–371,, 2020. 

Lin, Y., Zhang, H., Lin, H., Gamba, P. E., and Liu, X.: Incorporating synthetic aperture radar and optical images to investigate the annual dynamics of anthropogenic impervious surface at large scale, Remote Sens. Environ., 242, 111757,, 2020. 

Liu, J., Liu, M., Tian, H., Zhuang, D., Zhang, Z., Zhang, W., Tang, X., and Deng, X.: Spatial and temporal patterns of China's cropland during 1990–2000: an analysis based on Landsat TM data, Remote Sens. Environ., 98, 442–456,, 2005a. 

Liu, J., Liu, M., Zhuang, D., Zhang, Z., and Deng, X.: Study on spatial pattern of land-use change in China during 1995–2000, Sci. China Ser. D-Earth Sci., 46, 3732003,, 2005b. 

Liu, J., Zhang, Z., Xu, X., Kuang, W., Zhou, W., Zhang, S., Li, R., Yan, C., Yu, D., Wu, S., and Jiang, N.: Spatial patterns and driving forces of land use change in China during the early 21st century, J. Geogr. Sci., 20, 483–494,, 2010. 

Liu, J., Kuang, W., Zhang, Z., Xu, X., Qin, Y., Ning, J., Zhou, W., Zhang, S., Li, R., Yan, C., Wu, S., Shi, X., Jiang, N., Yu, D., Pan, X., and Chi, W.: Spatiotemporal characteristics, patterns, and causes of land-use changes in China since the late 1980s, J. Geogr. Sci., 24, 195–210,, 2014. 

Liu, X., Hu, G., Chen, Y., Li, X., Xu, X., Li, S., Pei, F., and Wang, S.: High-resolution multi-temporal mapping of global urban land using Landsat images based on the Google Earth Engine Platform, Remote Sens. Environ., 209, 227–239,, 2018. 

Lu, D. and Weng, Q.: Spectral mixture analysis of the urban landscape in Indianapolis with Landsat ETM plus imagery, Photogramm. Eng. Rem. S., 70, 1053–1062,, 2004. 

Lu, D. and Weng, Q.: Use of impervious surface in urban land-use classification, Remote Sens. Environ., 102, 146–160,, 2006. 

Lu, D., Tian, H., Zhou, G., and Ge, H.: Regional mapping of human settlements in southeastern China with multisensor remotely sensed data, Remote Sens. Environ., 112, 3668–3679,, 2008. 

Lu, D., Li, G., Kuang, W., and Moran, E.: Methods to extract impervious surface areas from satellite images, Int. J. Digit. Earth, 7, 93–112,, 2014. 

Lu, D., Li, L., Li, G., Fan, P., Ouyang, Z., and Moran, E.: Examining spatial patterns of urban distribution and impacts of physical conditions on urbanization in coastal and inland metropoles, Remote Sens., 10, 1101,, 2018. 

Ma, Q., He, C., Wu, J., Liu, Z., Zhang, Q., and Sun, Z.: Quantifying spatiotemporal patterns of urban impervious surfaces in China: an improved assessment using nighttime light data, Landscape Urban Plan., 130, 36–49,, 2014. 

Ning, J., Liu, J., Kuang, W., Xu, X., Zhang, S., Yan, C., Li, R., Wu, S., Hu, Y., Du, G., Chi, W., Pan, T., and Ning, J.: Spatiotemporal patterns and characteristics of land-use change in China during 2010–2015, J. Geogr. Sci., 28, 547–662,, 2018. 

Nowak, D. J. and Greenfield, E. J.: Tree and impervious cover in the United States, Landscape Urban Plan., 107, 21–30,, 2012. 

Peng, J., Shen, H., Wu, W., Liu, Y., and Wang, Y.: Net primary productivity (NPP) dynamics and associated urbanization driving forces in metropolitan areas: a case study in Beijing City, China, Landscape Ecol., 31, 1077–1092,, 2016. 

Pesaresi, M., Huadong, G., Blaes, X., Ehrlich, D., Ferri, S., Gueguen, L., Halkia, M., Kauffmann, M., Kemper, T., Lu, L., Marin-Herrera, M. A., Ouzounis, G. K., Scavazzon, M., Soille, P., Syrris, V., and Zanchetta, L.: A Global Human Settlement Layer from optical hr/vhrrs data: concept and first results, Ieee J-Stars., 6, 2102–2131,, 2013. 

Reba, M. and Seto, K. C.: A systematic review and assessment of algorithms to detect, characterize, and monitor urban land change, Remote Sens. Environ., 242, 111739,, 2020. 

Schneider, A., Friedl, M. A., and Potere, D.: Mapping global urban areas using MODIS 500 m data: new methods and datasets based on “urban ecoregions”, Remote Sens. Environ., 114, 1733–1746,, 2010. 

Seto, K. C., Guneralp, B., and Hutyra, L. R.: Global forecasts of urban expansion to 2030 and direct impacts on biodiversity and carbon pools, P. Natl. Acad. Sci. USA, 109, 16083–16088,, 2012. 

Sexton, J. O., Song, X., Huang, C., Channan, S., Baker, M. E., and Townshend, J. R.: Urban growth of the Washington, D.C.–Baltimore, MD metropolitan region from 1984 to 2010 by annual, Landsat-based estimates of impervious cover, Remote Sens. Environ., 129, 42–53,, 2013. 

Small, C. and Milesi, C.: Multi-scale standardized spectral mixture models, Remote Sens. Environ., 136, 442–454,, 2013. 

Wang, H., Lu, S., Wu, B., and Li, X.: Advances in remote sensing of impervious surfaces extraction and its applications, Adv. Earth Sci., 28, 327–336, 2013. 

Wang, P., Huang, C., Brown de Colstoun, E. C., Tilton, J. C., and Tan, B.: Global Human Built-up And Settlement Extent (HBASE) dataset from Landsat, NASA Socioeconomic Data and Applications Center (SEDAC), Palisades, NY,, 2017. 

Weng, Q.: Remote sensing of impervious surfaces in the urban areas: requirements, methods, and trends, Remote Sens. Environ., 117, 34–49,, 2012. 

Weng, Q., Lu, D., and Schubring, J.: Estimation of land surface temperature–vegetation abundance relationship for urban heat island studies, Remote Sens. Environ., 89, 467–483,, 2004. 

Wu, C. and Murray, A. T.: Estimating impervious surface distribution by spectral mixture analysis, Remote Sens. Environ., 84, 493–505,, 2003. 

Wu, J., Xiang, W., and Zhao, J.: Urban ecology in China: historical developments and future directions, Landscape Urban Plan., 125, 222–233,, 2014. 

Xu, H.: Modification of normalised difference water index (NDWI) to enhance open water features in remotely sensed imagery, Int. J. Remote Sens., 27, 3025–3033,, 2006. 

Xu, J., Zhao, Y., Zhong, K., Zhang, F., Liu, X., and Sun, C.: Measuring spatio-temporal dynamics of impervious surface in Guangzhou, China, from 1988 to 2015, using time-series Landsat imagery, Sci. Total Environ., 627, 264–281,, 2018. 

Xu, X. and Min, X.: Quantifying spatiotemporal patterns of urban expansion in China using remote sensing data, Cities, 35, 104–113,, 2013. 

Yang, L., Jin, S., Danielson, P., Homer, C., Gass, L., Bender, S. M., Case, A., Costello, C., Dewitz, J., Fry, J., Funk, M., Granneman, B., Liknes, G. C., Rigge, M., and Xian, G.: A new generation of the United States National Land Cover Database: Requirements, research priorities, design, and implementation strategies, Isprs J. Photogramm., 146, 108–123,, 2018. 

Zhang, C., Kuang, W., Wu, J., Liu, J., and Tian, H.: Industrial land expansion in rural China threatens food and environmental securities, Front. Env. Sci. Eng., 15, 29,, 2021.  

Zhang, L. and Weng, Q.: Annual dynamics of impervious surface in the Pearl River delta, China, from 1988 to 2013, using time series Landsat imagery, ISPRS J. Photogramm., 113, 86–96,, 2016. 

Zhang, X., Liu, L., Wu, C., Chen, X., Gao, Y., Xie, S., and Zhang, B.: Development of a global 30 m impervious surface map using multisource and multitemporal remote sensing datasets with the Google Earth Engine platform, Earth Syst. Sci. Data, 12, 1625–1648,, 2020. 

Zhang, Y., Odeh, I. O. A., and Han, C.: Bi-temporal characterization of land surface temperature in relation to impervious surface area, NDVI and NDBI, using a sub-pixel image analysis, Int. J. Appl. Earth. Obs., 11, 256–264,, 2009. 

Zhang, Z., Wang, X., Zhao, X., Liu, B., Liu, F., Yi, L., Zuo, L., Wen, Q., Xu, J., and Hu, S.: A 2010 update of National Land Use/Cover Database of China at 1:100 000 scale using medium spatial resolution satellite images, Remote Sens. Environ., 149, 142–154,, 2014. 

Zhang, Z., Wang, W., Cheng, M., Liu, S., Xu, J., He, Y., and Meng, F.: The contribution of residential coal combustion to PM2.5 pollution over China's Beijing-Tianjin-Hebei region in winter, Atmos. Environ., 159, 147–161,, 2017. 

Zhou, Y., Smith, J. S., Elvidge, C. D., Zhao, K., Thomson, A. M., and Imhoff, L. M.: A cluster-based method to map urban area from DMSP/OLS nightlights, Remote Sens. Environ., 147, 173–185,, 2014. 

Zhou, Y., Smith, J. S., Zhao, K., Imhoff, L. M., Thomson, A. M., Bondlamberty, B., Asrar, G., Zhang, X., He, C., and Elvidge, C. D.: A global map of urban extent from nightlights, Environ. Res. Lett., 10, 054011,, 2015. 

Zhou, Y., Li, X., Asrar, G. R., Smith, S. J., and Imhoff, M.: A global record of annual urban dynamics (1992–2013) from nighttime lights, Remote Sens. Environ., 219, 206–220,, 2018. 

Short summary
We propose a hierarchical principle for remotely sensed urban land use and land cover change for mapping intra-urban structure and component dynamics. China’s Land Use/cover Dataset (CLUD) is updated, delineating the imperviousness and green surface conditions in cities from 2000 to 2018. The newly developed datasets can be used to enhance our understanding of urbanization impacts on ecological and regional climatic conditions and on urban dwellers' environments.