Next Article in Journal
Numerical Simulation of an Indirect Contact Mobilized Thermal Energy Storage Container with Different Tube Bundle Layout and Fin Structure
Next Article in Special Issue
Study on the Spatial and Temporal Evolution of NDVI and Its Driving Mechanism Based on Geodetector and Hurst Indexes: A Case Study of the Tibet Autonomous Region
Previous Article in Journal
Nature-Based Solutions for Cooling in High-Density Neighbourhoods in Shenzhen: A Case Study of Baishizhou
Previous Article in Special Issue
Study on Rationality of Public Fitness Service Facilities in Beijing Based on GIS
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Changes in Impervious Surfaces in Lhasa City, a Historical City on the Qinghai–Tibet Plateau

1
School of Environmental Science and Engineering, Tianjin University, Tianjin 300350, China
2
School of Science, Tibet University, Lhasa 850000, China
3
School of Geography, South China Normal University, Guangzhou 510631, China
4
The Joint Laboratory of Plateau Surface Remote Sensing, Tibet University, Lhasa 850000, China
*
Author to whom correspondence should be addressed.
Sustainability 2023, 15(6), 5510; https://doi.org/10.3390/su15065510
Submission received: 20 February 2023 / Revised: 8 March 2023 / Accepted: 18 March 2023 / Published: 21 March 2023
(This article belongs to the Special Issue Advances in Applications of Remote Sensing for Urban Sustainability)

Abstract

:
Impervious surface cover reflects the urban environment and urban expansion. Lhasa City is a historical city and one of the most populous on the Qinghai–Tibetan Plateau, and has been experiencing rapid urbanization in recent years. Analyzing the impervious surface distribution can reveal urban development characteristics and provide data for sustainable urban planning to protect the heritage. This study explored the spatial and temporal changes and expansion patterns of impervious surfaces in different zones of Lhasa City. Impervious surface maps (2014 and 2021) were extracted from Gaofen-1 images with a high spatial resolution (2 m) using an object-based image analysis method. Next, a gravity center, standard deviational ellipses and landscape indices were used to characterize impervious surface expansions in different zones. The result indicated that the impervious surface in Lhasa expanded from 51.149 km2 in 2014 to 63.299 km2 in 2021. The growth rates of impervious surfaces inside the Environmental Coordination zone were lower than in the zones outside. From 2014 to 2021, the impervious surface of Lhasa expanded in the southeast direction. Infilling and consolidation were the primary impervious surface development patterns. The expansion of the impervious surface was related to topography, population, and economic and policy factors.

1. Introduction

Impervious surfaces are defined as materials that prevent water from infiltrating the soil, including parking lots, rooftops, roads, and sidewalks [1]. Impervious surface changes have been considered an efficient indicator for monitoring urban changes and play an important role in natural environment assessment [2,3,4]. The urbanization process of an area and the direction of urban expansion are reflected in the geometry, location, magnitude, and spatial pattern of impervious surfaces, as well as the ratio of pervious to impervious surfaces [5,6]. Increases in the area of impervious surfaces cause a variety of environmental problems, including water quality degradation [7], urban heat islands [8,9], and urban waterlogging [10]. As a result, the estimation and mapping of impervious surface areas is extremely important for diverse issues in environmental science that are central to global environmental changes and human–environment interactions, environmental management, and urban planning, including the construction of infrastructure and achieving sustainable development [5].
Impervious surfaces have already been widely used to explore spatial and temporal changes in urban development and to analyze urban environmental effects [11,12,13]. A variety of approaches have been employed to identify urban land use changes, including concentric rings [14], Moran’s I index [15], bands [16], transects [17], gravity center analysis and standard deviational ellipses (SDE) [4], as well as landscape indices [18,19]. Single methods can be used to evaluate different aspects of urban dynamics, such as orientation and patterns. However, conducting these analyses alone cannot achieve a comprehensive understanding of the spatiotemporal heterogeneity of impervious surfaces in complex urban environments over extended time periods, such as at the decade scale. It is therefore important to investigate the expansion of urban impervious surfaces at various scales using methods in combination.
Much effort has recently been made to map impervious surfaces using remote sensing techniques. Major approaches over the past decade to impervious surface extraction include an impervious surface index, a support vector machine, a spectral mixture analysis, a decision tree model, a regression tree, a random forest (RF), an object-based image analysis (OBIA), and an artificial neural network. Among these approaches, many researchers have used the OBIA method to extract impervious surfaces from high-spatial-resolution images [20,21,22]. Rather than relying on individual pixels, this method uses nearby pixel groups (objects). These pixel groups can be segmented based on geometric, spectral, contextual, and textural information. This method can effectively reduce the spectral variation present in a variety of land cover types, thereby improving the precision of classification [23,24,25]. Numerous investigations have indicated that the OBIA method exhibits high efficiency and precision in revealing the distribution of impervious surfaces [26,27,28]. The nearest neighbor (NN) classifier is popular when performing the OBIA method; however, the NN classifier may be less effective with high-dimensional data due to issues related to feature correlation [29,30]. To reduce these limitations and enhance object-based classification accuracy, machine learning classifiers have attracted increasing research interests. Machine learning classifiers including the support vector machine, the classification and regression tree, and the RF have been reported to outperform most classifiers in many studies [23,31,32].
Since the beginning of the 21st century, Chinese cities have developed rapidly, and the urban geographical structure and spatial pattern have been continuously optimized and upgraded. The implementation of the “Western Development” strategy has significantly driven the urbanization development of the Qinghai–Tibet Plateau. Over the past 30 years, the Tibetan Autonomous Region has achieved rapid economic development with the support of national policies. Lhasa City is the capital of the Tibet Autonomous Region, with a history of over 1300 years. Lhasa City is a multifunctional urban center defined by the intersection between history and modernity. It is among the most populous cities on the Qinghai–Tibet Plateau and a center of Tibetan national and regional culture, with extremely elevated historical and cultural value and unique social and economic functions. The development status of Lhasa City is an important window to characterize the development level of the Tibet region, and the impervious surface can well represent the urban expansion and development trends. However, the ecological environment of the Qinghai–Tibet Plateau is extremely fragile and any human disturbance can easily change the ecological functions and ecosystem structures [33]. As a major human activity zone in the plateau, the land cover of Lhasa City experienced significant changes accompanied by urbanization [34,35]. The expansion of impervious surfaces in Lhasa City has changed the land surface temperature [36,37], eco-environmental quality [38], and riverine hydrological processes [39]. Furthermore, rapid urbanization has led to the rapid expansion of urban space, and the overall pattern and traditional features of many historical and cultural cities have been damaged during large-scale urban and rural construction [40,41,42]. Heritage conservation is an indispensable factor for sustainable development [43]. Balancing heritage conservation and urban development is very important in sustainable urban development plans [44]. This study is the first to analyze the impervious distribution of Lhasa City and heritage conservation areas with high spatial resolution images, which can reveal urban development characteristics and provide more precise fundamental data for sustainable urban planning to protect the heritage.
In this study, the Chengguan District of Lhasa City was used as the study area, and a variety of methods were adopted to perform a spatial and temporal analysis of its impervious surfaces in different zones from 2014 to 2021. The specific objectives of the study are to (1) extract impervious surface maps using the OBIA method and analyze the spatial and temporal dynamics of impervious surfaces over the study period; (2) identify the expansion direction of impervious surfaces by gravity centers and the SDE method; and (3) explore the landscape pattern of impervious surfaces using the landscape indices of the city and heritage conservation areas.

2. Materials and Methods

2.1. Study Area and Data

2.1.1. Study Area

Lhasa has an average altitude of 3656 m and a population of 867,891, making it the one of the most populous cities on the Qinghai–Tibet Plateau. Lhasa is among the highest cities in the world. Since the mid-17th century, the city has been both the administrative and religious capital of Tibet and is home to numerous culturally significant Tibetan Buddhist sites, including the Norbulingka Palaces, Jokhang Temple, and Potala Palace. Chengguan District is found in the valley plain area where the Lhasa River flows. It has a plateau mountain climate. The rainy season occurs from July to September every year and it rains frequently at night. The unique geographical location gives the area an average annual sunshine duration of over 3000 h and it is known as the “Sunlight City”. The urbanization rate is 42.40% and the GDP is CNY 47.925 billion. In 2020, the resident population of Chengguan District was 473,586, accounting for 54.57% of the total population of Lhasa.
According to the “Regulations on Protection of Famous Historical and Cultural Cities of Lhasa,” the study area (Figure 1) was divided into five zones: the Central City zone (Chengguan District), the Environmental Coordination zone, the Historic zone, and two historical and cultural districts—Bakuojie zone and Caiyicun zone. Caiyicun is a historical and cultural district established in 2019.

2.1.2. Data Recourses and Preprocessing

Four images from 2014 to 2021 were acquired from the GF-1 satellite at the China Centre for Resources Satellite Data and Application (CRESDA) (https://data.cresda.cn/, accessed on 18 March 2023) and were selected to extract impervious surfaces. Table 1 lists the images used in this study. The images selected in the study had 2 m panchromatic and 8 m multispectral bands. The selected images were almost cloud-free, with good atmospheric conditions and high image quality.
Radiometric correction, geometric correction, image registration, and image mosaicking were used to preprocess the original remote sensing images in order to obtain remote sensing image data for the study area.
The GF-1 satellite imagery has four multispectral bands including blue, green, red, and near-infrared spectra [45]. Radiometric correction was applied to convert raw digital image data from satellite to apparent reflectance [46]. The FLAASH Atmospheric Correction Model was applied to the multispectral images to minimize the atmospheric effects. Next, geometric correction was applied using an RPC sensor model and the images were reprojected into the WGS84/UTM coordinate systems. The NNDiffusion pan-sharpening method was implemented to obtain a fused image with a spatial resolution from 8 to 2 m for multispectral bands.

2.2. Methods

The OBIA method for the extraction of impervious surfaces is composed of two main steps: image segmentation and object-based classification. The OBIA method was performed with eCognition Developer 9.0 software. Furthermore, an accuracy assessment was performed to evaluate the classification results. Figure 2 shows the workflow of the study.

2.2.1. Mapping Impervious Surfaces

Image Segmentation

Image segmentation refers to the process of dividing an image into a group of adjacent pixels (objects) that exhibit similar spatial features, spectral information, and textual features [20]. The multiscale segmentation algorithm was used to generate segmented images. In the algorithm, scale is a key parameter which determines spectral variation within objects and their size. Greater scale values result in larger and heterogeneous objects. Furthermore, the shape and compactness parameter values are also important because they determine the weight of the color and the smoothness of segmented objects. Using different segmentation parameters leads to different results; thus, the generation of discernible image objects on the greatest possible scale can be accomplished through trial and error. In the present work, the three parameters (scale = 70, shape = 0.3, and compactness = 0.4) were chosen using an iterative experimental approach at different levels.

Feature Selection

The feature selection plays a key role in classification and can enhance the computational performance of classifiers by eliminating redundant information [47]. The Trimble eCognition software provides various sets of feature options; besides image statistics, many remote sensing indices, including the normalized difference water index (NDWI) and normalized difference vegetation index (NDVI), are an option that can improve the classification results. Moreover, the perpendicular impervious surface index (PISI) is used as a feature in the OBIA method because of its good performance in mapping impervious surfaces in high spatial resolution images with great accuracy [28]. Experiments were performed during the feature selection step to determine the appropriate statistics, and the most widely used features including the shape and texture (derived from the grey-level co-occurrence matrices (GLCM)) were also extracted, including the homogeneity, dissimilarity, entropy, and correlation, while spectral statistics (mean value and brightness) were defined for the classifier.

RF Classifier

RF is a widely used ensemble learning method that has been used effectively for classification to solve a range of pixel- and object-based tasks owing to its processing speed, accuracy, and robustness [48,49,50,51]. The main concept underlying the RF algorithm is the creation of a decision forest constructed from a multitude of decision trees (DT) classifiers trained based on various bootstrap sub-samples from the original dataset [52]. For unknown samples, class labels are predicted based on majority voting after combining the outputs of all individual trees [52,53,54]. Each DT utilizes approximately two-thirds of the training samples selected using the bootstrapping technique, while the remaining one-third samples are used for out-of-bag error estimation.

Accuracy Assessment

High-resolution images obtained from Google Earth imagery for 2014 and 2021 were utilized for the evaluation of the accuracy of the impervious surface map. Five hundred verification points were randomly selected for each year. The confusion matrix method was applied for the extraction of the following four accuracy indicators: overall accuracy (OA), user’s accuracy (UA), producer’s accuracy (PA), and overall kappa coefficient [55].
The OA is the proportion of correctly classified pixels from the total accuracy assessment pixels. The PA refers to the probability that a certain land cover of an area on the ground is correctly classified and represents the omission error. The UA refers to the proportion of correctly classified pixels within the image and represents the commission error. The kappa coefficient is a statistical measure used to test interrater reliability [56].

2.2.2. Impervious Surface Mean Center Analysis and Standard Deviational Ellipse

The gravity center [57], assigned as the impervious surface mean center (ISMC), was applied to determine the direction of urban expansion based on impervious surface maps from 2014 to 2020. The coordinates of the ISMC are calculated as the following equation:
x ¯ = i = 1 n x i n y ¯ = i = 1 n y i n
where x and y are the longitude and latitude coordinates of the ISMC, respectively; xi and yi are the coordinates of the ith grid in the impervious surface map; and n is the grid number of the impervious surface map.
Determining the spatial patterns of a city facilitates the comprehension of the development direction and degree of urbanization [58]. Therefore, SDE were applied to illustrate the spatial expansion of urbanization. Based on the impervious surface images, four SDE parameters were calculated: the mean center, long and short axis, and azimuth angle. The mean center and the azimuth angle were determined to capture the urbanization change direction of the entire region [50,51]. The mean center of the SDE was the previously calculated ISMC. The azimuth of the SDE is calculated as the following equations:
tan θ = i = 1 n x ˜ i 2 i = 1 n y ˜ i 2 + i = 1 n x ˜ i 2 i = 1 n y ˜ i 2 2 4 i = 1 n x ˜ i y ˜ i 2 2 i = 1 n x ˜ i y ˜ i
x ˜ i = x i x ¯ y ˜ i = y i y ¯
where xi and yi are the deviations of the i-th grid center and the ISMC for the x- and y-axis, respectively. The standard deviations, σx and σy, for the x- and y-axis, respectively, are calculated as follows:
σ x = 2 i = 1 n x ˜ i cos θ y ˜ i sin θ 2 n σ y = 2 i = 1 n x ˜ i sin θ + y ˜ i cos θ 2 n
The long and short axis and azimuth angle describe the direction and range of impervious surface distributions. The ratio of the long axis to the short axis reflects the clustering or dispersion degree of the impervious surface. A greater value of the ratio demonstrates stronger orientational effects of the impervious surface. By contrast, a ratio of 1 reflects no directional trend of the data. In this work, variation of the SDE of the impervious surface helped to identify whether the spatial distribution of the impervious surface was elongated and showed a particular direction of the impervious surface.

2.2.3. Analysis of the Changes in Landscape Indices

Landscape indices quantitatively reflect the configuration and composition of urban landscapes [13,59]. In the present study, the impervious surfaces expansion degree was calculated based on these indices. The analysis was performed with FRAGSTATS 4.2 software using an eight-neighbor rule. Five landscape indices were adopted to characterize the impervious surface changes which are commonly used in the study of urban expansion and its environmental impacts. These landscape indices were the number of patches (NP), patch density (PD), mean patch size (MPS), largest patch index (LPI), landscape shape index (LSI), and aggregation index (AI) (Table 2).

3. Results

3.1. Impervious Surface Changes and Accuracy Assessment

The result of the accuracy assessment is presented in Table 3. In 2014, the OA and kappa coefficient of the impervious surface map were 96.60% and 89.41%, respectively. In 2021, these values were 95.80% and 89.37%, respectively. As such, the results were considered reliable and could be used for further analyses.
The impervious surface maps for 2014 and 2021 are shown in Figure 3. Table 4 displays the impervious surface areas in different zones. In the Central City zone, the impervious surface area rose from 51.15 km2 in 2014 to 63.30 km2 in 2021, for an overall increase rate of 23.75% and a change rate of 1.74 km2/a. This indicates a clear expansion of impervious surfaces in Lhasa from 2014 to 2021. The added impervious surface was most concentrated in the southeast of the Central City zone. The impervious surfaces increasing rates in the Central City and Caiyicun were higher than in other zones. Furthermore, the Bakuojie zone had the highest impervious coverage in both 2014 and 2021. Moreover, impervious surface coverage of the Historic and Caiyicun zones were similar in 2014, while impervious surface coverage of the Caiyicun zone were much higher than the Historic zone in 2021.

3.2. Spatial and Temporal Analysis of ISMC and SDE

To illustrate the barycenter changes during urban development, the ISMCs of the entire area were identified based on an impervious surface map. The ISMC of the entire area in 2014 and 2021 is displayed in Figure 4 and Figure 5. In 2014, the ISMC was located at 91.146° E, 29.662° N, and moved to 91.148° E, 29.660° N in 2021. The changes in ISMC revealed the direction expansion. Over time, the impervious surfaces throughout the entire area expanded in a southeast direction. This suggested that the direction for urban expansion was towards the southeast.
The SDE were applied to determine the orientation and spatial and temporal development trends of impervious surfaces. The SDE of impervious surfaces in 2014 and 2021 are shown in Figure 4. According to Figure 4, the impervious surface presents an obvious orientation and direction. Between 2014 and 2021, the azimuth increased from 93.634° to 94.053° (Table 5), showing that the impervious surface expansion was oriented from west to east. The long and short axis, and their ratio of the SDE represented the clustering or dispersion degree of impervious surfaces. If the ratio is close to 1, the expansion of the impervious surface does not have a certain principal. The ratio was 2.172 and 2.119 in 2014 and 2021, respectively, indicating that the impervious surface had a significant distinct direction. The long axis increased from 8108.17 to 8493.82 m and the short axis increased from 3732.33 to 4009.02 m. This indicates a dispersion trend of impervious surfaces in both directions.

3.3. Landscape Indices of Different Zones

The landscape indices for the five zones in 2014 and 2021 are shown in Figure 5.
NP increased in the Central City zone but decreased in the other four zones. The Bakuojie and Caiyicun zones had the lowest NP in 2014 and 2021. The MSP in all four zones increased, with no obvious change in the Central City zone. The Central City zone had the lowest MSP. The LPI increased in the Central City zone, while no obvious changes were observed in other zones. The Central City zone had the lowest LPI and the Bakuojie zone had the highest LPI. This indicated that the dominance of impervious surfaces clearly increased only in the Central City zone. The LSI declined in the five zones, demonstrating that the spatial pattern of the impervious surfaces was simplified and regularized. AI increased in the five zones, indicating that as the impervious surface expanded, the impervious surface agglomerated. The Bakuojie and Caiyicun zones had the lowest LSI and relatively high AI compared to the other zones. This indicates that the impervious surfaces in the Bakuojie and Caiyicun zones have strong connectivity and simple shapes.

4. Discussion

4.1. Expansion of Impervious Surface

The results indicated that the impervious surface experienced clear expansion in Lhasa. However, the analysis results demonstrated that the impervious surface in different zones of Lhasa underwent different growth patterns.
From 2014 to 2021, the growth rates of impervious surfaces in the Central City and Caiyicun zones were 23.755% and 23.157%, respectively. The growth rates in the Historic, Environmental Coordination, and Bakuojie zones, which are inside the Environmental Coordination zone, were 3.017%, 4.812%, and 2.233%, respectively. This indicates that the expansion area was mainly outside the Environmental Coordination zone. The impervious surface covered most areas in the Historic zone, Environmental Coordination zone, and Bakuojie zone in 2014. With the continuous expansion of the urban areas, the available lands for construction reduced; therefore, the expansion area and expansion rate were limited. Moreover, because Lhasa is one of the National Famous Historical and Cultural Cities in China, new construction and expansion activities in conservation zones are limited. Therefore, the expansion rate and areas in the Historic, Environmental Coordination, and Bakuojie zones are lower than those in the Central City and Caiyicun zones.
The ratio of the long axis to the short axis of the impervious surface SDE showed that there was a clear direction in the distribution of impervious surfaces. Changes in the long and short axes showed that the impervious surfaces exhibited a dispersion trend. Between 2014 and 2021, the shift of the ISMC showed that the impervious surfaces expanded toward the southeast throughout the entire region. Chengguan District is a valley plain distributed from west to east. The topography makes urban development have a clear direction. The impervious surface distribution also showed a west-to-east trend. The newly added impervious surfaces were mainly located in the Najin zone, which is southeast of Lhasa. The expansion direction was consistent with the urban planning requirements of the city. With the expansion of the city, according to the Overall Urban Plan for Lhasa (2009–2020) [60], the government has actively developed the new urban area to ease the pressure coming from the central district and balance the preservation of historic towns and the rapid development. As a result, the ISMC has shifted toward the southeast.

4.2. Landscape Pattern Changes of Impervious Surface

In the Central City zone, the values of NP and LPI showed a clear increasing trend, indicating that the impervious surface had an obvious expansion and the dominance of impervious surfaces became more prominent. This is caused by economic development, population growth, and urban expansion. In the other four conservation zones, the values of NP showed a decreasing trend and the LPI showed no obvious changes. Owing to the limited impervious surface expansion in these zones, the dominance of the impervious surface did not exhibit an obvious change. In addition, the decreased values in mean patch size indicated that most of the added patches were connected to previous impervious surface patches, meaning that the expansion of the impervious surface was organized. This indicates that infilling and consolidation was the impervious surface development pattern inside these areas. With the increase in impervious surfaces, landscape fragmentation decreased. The low LSI value suggests that the spatial shape of the impervious surface became increasingly regular, and the expansion of the impervious surface was orderly. The high AI value means that with the increase in impervious surfaces, the distribution of impervious surfaces became increasingly concentrated, and the degree of aggregation increased. Since 2013, a series of regulations have been released to protect the old city and promote the sustainable development of Lhasa [61,62,63]. In the conservation zone, the regions are well planned, especially in the Bakuojie and Caiyicun zones. The Bakuojie and Caiyicun zones had the lowest NP and strongest connectivity of impervious surfaces. In these zones, the spatial shape of impervious surfaces was regular, and the landscape inside the areas became more stable, regular, and aggregated.

4.3. Driving Forces of Impervious Surface Expansion

Previous studies have shown that topography is an important factor restricting urban expansion and influencing development patterns [12,64]. The distribution of the impervious surfaces indicated that the urban space of Lhasa is restricted by the surrounding mountains.
Population and economic growth are the driving forces of urban development [65,66,67]. Generally, rapid population and economic growth are accompanied by significant ISA expansion. With the implementation of the “Western Development” strategy and the “One Belt and One Road” initiative, the economy of Lhasa has developed rapidly. In addition, urbanization accelerates population and economic growth [68].
The population of Lhasa increased from 6.56 × 105 in 2014 to 8.68 × 105 in 2021, with an increasing rate of 32.32%. Population has always been an active factor in urban formation and has an important impact on urban development [69]. The population growth in cities leads to increasing demands for housing, infrastructure such as roads, water conservancy, and electricity, as well as public service facilities for education, employment, medical care, and entertainment. The construction of new buildings and facilities occupied a large amount of land and resulted in the increase in impervious surfaces.
The gross domestic product (GDP) of Lhasa increased from CNY 34.74 billion in 2014 to CNY 74.18 billion in 2021. As a historical city and a popular tourism destination, Lhasa City has developed tourism and service industries. In 2014, the total GDP of Lhasa was CNY 34.745 billion, and the primary, secondary and tertiary industries were CNY 1.294 billion, 12.775 billion, and 20.676 billion, respectively. In 2021, the total GDP of Lhasa was CNY 74.184 billion, and the primary, secondary and tertiary industries were CNY 2.445 billion, 27.808 billion, and 43.931 billion, respectively [70,71]. The proportion of secondary and tertiary industries increased from 96.27% in 2014 to 96.70% in 2021. The rapid development of the secondary and tertiary industries led to an expansion in industrial land and service facilities, resulting in an increase in impervious surfaces. In addition, the development of secondary and tertiary industries and better services attracted more people to migrate to the city and further increase the area of impervious surfaces.
Policy is another important factor influencing urban development patterns. According to the Overall Urban Plan for Lhasa [60], the functions of the Central City should be optimized and the historic buildings should be well protected. This policy proposed the future expansion layouts of Lhasa: renovate the old city zones to improve the efficiency of space utilization; establish new urban cores and promote the development of the surrounding areas of the old city zones; and increase the connection between the new urban areas and the old city zones, to form a polycentric urban structure. Furthermore, with the implementation of the “Regulations on the Protection of the Old Town of Lhasa” [61], the “Regulations on the Protection of the Ancient Village of Lhasa” [62], and the “Historic and Cultural City Protection Plan for Lhasa” [63], the government put much effort into protecting the old city. At the same time, the impervious surfaces in the Central City zone tended to be saturated; thus, the increasing rates of saturation within the Environmental Coordination zone were lower than in the Central City and Caiyicun zones. Moreover, the government has made environmental protection a priority together with pursuing high-quality development for the city [72] and balancing rapid urbanization and ecological protection. Farmland and ecological land were protected by policies, which also resulted in an infilling and consolidation growth pattern of impervious surfaces in Lhasa City.

5. Conclusions

In this study, the impervious surface maps of 2014 and 2021 were extracted by an OBIA method using Gaofen-1 images of Lhasa. This work utilized a variety of complementary methods to characterize and analyze the spatial and temporal changes in the impervious surfaces in different zones of Lhasa. The main findings are as follows:
  • The impervious surface area in Lhasa grew significantly from 51.149 km2 in 2014 to 63.299 km2 in 2021, with an overall increase rate of 23.75% and change rate of 1.74 km2/a. The Environmental Coordination, Historic, and Bakuojie zones, which are inside the Environmental Coordination zone, had increasing rates of 3.017%, 4.812%, and 4.246%, respectively. The Caiyicun zone, which is outside the Environmental Coordination zone, had an increasing rate of 23.16%.
  • Between 2014 and 2021, the shift of the ISMC revealed that the impervious surfaces of Lhasa expanded toward the southeast. The changes in SDE showed that the impervious surface distribution had a clear direction; the impervious surface expansion was oriented from west to east, and the impervious surfaces had an observable dispersion trend.
  • In general, the landscape pattern in Lhasa City tends to be orderly, aggregated, and regularized. In the Central City zone, the impervious surface with an obvious expansion and dominance became more prominent. In the other four conservation zones, the dominance of impervious surfaces did not change significantly. Infilling and consolidation were the primary impervious surface development patterns in these areas. With the expansion of the impervious surface, the landscape in Lhasa City tended to be stable, regular, and aggregated, especially in the Bakuojie and Caiyicun zones.
The findings of this study will improve the understanding of the spatial and temporal variation in urban expansion and provide a foundation for the analysis of its environmental effects as well as an exploration of its driving factors. In addition, these results will help decision-makers to formulate urban planning policies and provide a reference for further research in this area.

Author Contributions

Methodology, software, validation, visualization, writing, S.W.; resources, data curation, X.T.; conceptualization, supervision, F.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by “The National Social Science Fund of China” grant number “21AMZ011” And “The Major Cultivation Fund for Philosophy and Social Sciences of South China Normal University” grant number “ZDPY2206”.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Arnold, 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]
  2. Yuan, F.; Bauer, M.E. Comparison of Impervious Surface Area and Normalized Difference Vegetation Index as Indicators of Surface Urban Heat Island Effects in Landsat Imagery. Remote Sens. Environ. 2007, 106, 375–386. [Google Scholar] [CrossRef]
  3. Sekertekin, A.; Zadbagher, E. Simulation of Future Land Surface Temperature Distribution and Evaluating Surface Urban Heat Island Based on Impervious Surface Area. Ecol. Indic. 2021, 122, 107230. [Google Scholar] [CrossRef]
  4. Xu, J.; Zhao, Y.; Zhong, K.; Zhang, F.; Liu, X.; Sun, C. Measuring Spatio-Temporal Dynamics of Impervious Surface in Guangzhou, China, from 1988 to 2015, Using Time-Series Landsat Imagery. Sci. Total Environ. 2018, 627, 264–281. [Google Scholar] [CrossRef]
  5. 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]
  6. Xu, H.; Wang, M. Remote Sensing-Based Retrieval of Ground Impervious Surfaces. J. Remote Sens. 2016, 20, 1270–1289. [Google Scholar] [CrossRef]
  7. Alberti, M.; Booth, D.; Hill, K.; Coburn, B.; Avolio, C.; Coe, S.; Spirandelli, D. The Impact of Urban Patterns on Aquatic Ecosystems: An Empirical Analysis in Puget Lowland Sub-Basins. Landsc. Urban Plan. 2007, 80, 345–361. [Google Scholar] [CrossRef]
  8. Bektaş Balçik, F. Determining the Impact of Urban Components on Land Surface Temperature of Istanbul by Using Remote Sensing Indices. Environ. Monit. Assess. 2014, 186, 859–872. [Google Scholar] [CrossRef]
  9. Hua, L.; Zhang, X.; Nie, Q.; Sun, F.; Tang, L. The Impacts of the Expansion of Urban Impervious Surfaces on Urban Heat Islands in a Coastal City in China. Sustainability 2020, 12, 475. [Google Scholar] [CrossRef] [Green Version]
  10. Siddik, M.S.; Tulip, S.S.; Rahman, A.; Islam, M.N.; Haghighi, A.T.; Mustafa, S.M.T. The Impact of Land Use and Land Cover Change on Groundwater Recharge in Northwestern Bangladesh. J. Environ. Manag. 2022, 315, 115130. [Google Scholar] [CrossRef]
  11. Xiao, R.; Ouyang, Z.; Zheng, H.; Li, W.; Schienke, E.W.; Wang, X. Spatial Pattern of Impervious Surfaces and Their Impacts on Land Surface Temperature in Beijing, China. J. Environ. Sci. 2007, 19, 250–256. [Google Scholar] [CrossRef]
  12. Wu, W.; Li, C.; Liu, M.; Hu, Y.; Xiu, C. Change of Impervious Surface Area and Its Impacts on Urban Landscape: An Example of Shenyang between 2010 and 2017. Ecosyst. Health Sustain. 2020, 6, 1767511. [Google Scholar] [CrossRef]
  13. Zhang, S.; Yang, K.; Ma, Y.; Li, M. The Expansion Dynamics and Modes of Impervious Surfaces in the Guangdong-Hong Kong-Macau Bay Area, China. Land 2021, 10, 1167. [Google Scholar] [CrossRef]
  14. Seto, K.C.; Fragkias, M. Quantifying Spatiotemporal Patterns of Urban Land-Use Change in Four Cities of China with Time Series Landscape Metrics. Landsc. Ecol. 2005, 20, 871–888. [Google Scholar] [CrossRef]
  15. Peng, J.; Liu, Y.; Shen, H.; Xie, P.; Hu, X.; Wang, Y. Using Impervious Surfaces to Detect Urban Expansion in Beijing of China in 2000s. Chin. Geogr. Sci. 2016, 26, 229–243. [Google Scholar] [CrossRef] [Green Version]
  16. Zhou, W.; Huang, G.; Pickett, S.T.A.; Cadenasso, M.L. 90 Years of Forest Cover Change in an Urbanizing Watershed: Spatial and Temporal Dynamics. Landsc. Ecol. 2011, 26, 645–659. [Google Scholar] [CrossRef]
  17. Luck, M.; Wu, J. A Gradient Analysis of Urban Landscape Pattern: A Case Study from the Phoenix Metropolitan Region, Arizona, USA. Landsc. Ecol. 2002, 17, 327–339. [Google Scholar] [CrossRef]
  18. Kowe, P.; Pedzisai, E.; Gumindoga, W.; Rwasoka, D.T. An Analysis of Changes in the Urban Landscape Composition and Configuration in the Sancaktepe District of Istanbul Metropolitan City, Turkey Using Landscape Metrics and Satellite Data. Geocarto Int. 2015, 30, 506–519. [Google Scholar] [CrossRef]
  19. Feng, S.; Fan, F. Spatiotemporal Changes of Landscape Pattern Using Impervious Surface in Guangdong-Hong Kong-Macao Greater Bay Area, China. Chin. J. Appl. Ecol. 2018, 29, 2907–2914. [Google Scholar] [CrossRef]
  20. Hu, X.; Weng, Q. Impervious Surface Area Extraction from IKONOS Imagery Using an Object-Based Fuzzy Method. Geocarto Int. 2011, 26, 3–20. [Google Scholar] [CrossRef]
  21. Sebari, I.; He, D.-C. Automatic Fuzzy Object-Based Analysis of VHSR Images for Urban Objects Extraction. ISPRS J. Photogramm. Remote Sens. 2013, 79, 171–184. [Google Scholar] [CrossRef]
  22. Zhang, X.; Xiao, P.; Feng, X. Impervious Surface Extraction from High-Resolution Satellite Image Using Pixel- and Object-Based Hybrid Analysis. Int. J. Remote Sens. 2013, 34, 4449–4465. [Google Scholar] [CrossRef]
  23. Ma, L.; Li, M.; Ma, X.; Cheng, L.; Du, P.; Liu, Y. A Review of Supervised Object-Based Land-Cover Image Classification. ISPRS J. Photogramm. Remote Sens. 2017, 130, 277–293. [Google Scholar] [CrossRef]
  24. Mallinis, G.; Koutsias, N.; Tsakiri-Strati, M.; Karteris, M. Object-Based Classification Using Quickbird Imagery for Delineating Forest Vegetation Polygons in a Mediterranean Test Site. ISPRS J. Photogramm. Remote Sens. 2008, 63, 237–250. [Google Scholar] [CrossRef]
  25. Laliberte, A.S.; Rango, A.; Havstad, K.M.; Paris, J.F.; Beck, R.F.; McNeely, R.; Gonzalez, A.L. Object-Oriented Image Analysis for Mapping Shrub Encroachment from 1937 to 2003 in Southern New Mexico. Remote Sens. Environ. 2004, 93, 198–210. [Google Scholar] [CrossRef]
  26. Myint, S.W.; Gober, P.; Brazel, A.; Grossman-Clarke, S.; Weng, Q. Per-Pixel vs. Object-Based Classification of Urban Land Cover Extraction Using High Spatial Resolution Imagery. Remote Sens. Environ. 2011, 115, 1145–1161. [Google Scholar] [CrossRef]
  27. Lichtblau, E.; Oswald, C.J. Classification of Impervious Land-Use Features Using Object-Based Image Analysis and Data Fusion. Comput. Environ. Urban Syst. 2019, 75, 103–116. [Google Scholar] [CrossRef]
  28. Feng, S.; Fan, F. Impervious Surface Extraction Based on Different Methods from Multiple Spatial Resolution Images: A Comprehensive Comparison. Int. J. Digit. Earth 2021, 14, 1148–1174. [Google Scholar] [CrossRef]
  29. Maxwell, A.E.; Warner, T.A.; Strager, M.P.; Conley, J.F.; Sharp, A.L. Assessing Machine-Learning Algorithms and Image- and Lidar-Derived Variables for GEOBIA Classification of Mining and Mine Reclamation. Int. J. Remote Sens. 2015, 36, 954–978. [Google Scholar] [CrossRef]
  30. Colkesen, I.; Kavzoglu, T. Comparative Evaluation of Decision-Forest Algorithms in Object-Based Land Use and Land Cover Mapping. In Spatial Modeling in GIS and R for Earth and Environmental Sciences; Elsevier: Amsterdam, The Netherlands, 2019; pp. 499–517. [Google Scholar]
  31. Li, M.; Ma, L.; Blaschke, T.; Cheng, L.; Tiede, D. A Systematic Comparison of Different Object-Based Classification Techniques Using High Spatial Resolution Imagery in Agricultural Environments. Int. J. Appl. Earth Obs. Geoinf. 2016, 49, 87–98. [Google Scholar] [CrossRef]
  32. Talukdar, S.; Singha, P.; Mahato, S.; Shahfahad; Pal, S.; Liou, Y.-A.; Rahman, A. Land-Use Land-Cover Classification by Machine Learning Classifiers for Satellite Observations—A Review. Remote Sens. 2020, 12, 1135. [Google Scholar] [CrossRef] [Green Version]
  33. Gao, Q.; Guo, Y.; Xu, H.; Ganjurjav, H.; Li, Y.; Wan, Y.; Qin, X.; Ma, X.; Liu, S. Climate Change and Its Impacts on Vegetation Distribution and Net Primary Productivity of the Alpine Ecosystem in the Qinghai-Tibetan Plateau. Sci. Total Environ. 2016, 554, 34–41. [Google Scholar] [CrossRef] [PubMed]
  34. Fu, S.; Zhang, X.; Kuang, W.; Guo, C. Characteristics of Changes in Urban Land Use and Efficiency Evaluation in the Qinghai–Tibet Plateau from 1990 to 2020. Land 2022, 11, 757. [Google Scholar] [CrossRef]
  35. Chen, T.; Jiao, J.; Wei, W.; Li, J.; Zhang, Z.; Yang, H.; Ma, H. Spatiotemporal Variation in the Land Use/Cover of Alluvial Fans in Lhasa River Basin, Qinghai–Tibet Plateau. Agriculture 2023, 13, 312. [Google Scholar] [CrossRef]
  36. Wu, J.; Sun, C.; Fan, F. An Analysis of the Spatiotemporal Distribution of Impervious Surface and Land Surface Temperature in Plateau Areas. J. South China Norm. Univ. Sci. Ed. 2020, 52, 92–100. [Google Scholar]
  37. Li, D.; Tian, P.; Luo, H.; Hu, T.; Dong, B.; Cui, Y.; Khan, S.; Luo, Y. Impacts of Land Use and Land Cover Changes on Regional Climate in the Lhasa River Basin, Tibetan Plateau. Sci. Total Environ. 2020, 742, 140570. [Google Scholar] [CrossRef]
  38. Han, W.; Lu, H.; Liu, G.; Wang, J.; Su, X. Quantifying Degradation Classifications on Alpine Grassland in the Lhasa River Basin, Qinghai-Tibetan Plateau. Sustainability 2019, 11, 7067. [Google Scholar] [CrossRef] [Green Version]
  39. Chen, J.; Kong, Y.; Mei, Y. Riverine Health Assessment Using Coordinated Development Degree Model Based on Natural and Social Functions in the Lhasa River, China. Int. J. Environ. Res. Public Health 2022, 19, 7182. [Google Scholar] [CrossRef] [PubMed]
  40. Li, X.; Hou, W.; Liu, M.; Yu, Z. Traditional Thoughts and Modern Development of the Historical Urban Landscape in China: Lessons Learned from the Example of Pingyao Historical City. Land 2022, 11, 247. [Google Scholar] [CrossRef]
  41. Arkaraprasertkul, N. Gentrifying Heritage: How Historic Preservation Drives Gentrification in Urban Shanghai. Int. J. Herit. Stud. 2019, 25, 882–896. [Google Scholar] [CrossRef]
  42. Chen, H.; Wang, L.; Waley, P. The Right to Envision the City? The Emerging Vision Conflicts in Redeveloping Historic Nanjing, China. Urban Aff. Rev. 2020, 56, 1746–1778. [Google Scholar] [CrossRef]
  43. Appendino, F. Balancing Heritage Conservation and Sustainable Development—The Case of Bordeaux. IOP Conf. Ser. Mater. Sci. Eng. 2017, 245, 062002. [Google Scholar] [CrossRef] [Green Version]
  44. Ragheb, A.; Aly, R.; Ahmed, G. Toward Sustainable Urban Development of Historical Cities: Case Study of Fouh City, Egypt. Ain Shams Eng. J. 2022, 13, 101520. [Google Scholar] [CrossRef]
  45. Zhou, Q.; Yu, Q.; Liu, J.; Wu, W.; Tang, H. Perspective of Chinese GF-1 High-Resolution Satellite Data in Agricultural Remote Sensing Monitoring. J. Integr. Agric. 2017, 16, 242–251. [Google Scholar] [CrossRef]
  46. Chang, S.; Wu, B.; Yan, N.; Zhu, J.; Wen, Q.; Xu, F. A Refined Crop Drought Monitoring Method Based on the Chinese GF-1 Wide Field View Data. Sensors 2018, 18, 1297. [Google Scholar] [CrossRef] [Green Version]
  47. Ma, L.; Fu, T.; Blaschke, T.; Li, M.; Tiede, D.; Zhou, Z.; Ma, X.; Chen, D. Evaluation of Feature Selection Methods for Object-Based Land Cover Mapping of Unmanned Aerial Vehicle Imagery Using Random Forest and Support Vector Machine Classifiers. ISPRS Int. J. Geo-Inf. 2017, 6, 51. [Google Scholar] [CrossRef] [Green Version]
  48. Guan, H.; Li, J.; Chapman, M.; Deng, F.; Ji, Z.; Yang, X. Integration of Orthoimagery and Lidar Data for Object-Based Urban Thematic Mapping Using Random Forests. Int. J. Remote Sens. 2013, 34, 5166–5186. [Google Scholar] [CrossRef]
  49. Lebourgeois, V.; Dupuy, S.; Vintrou, É.; Ameline, M.; Butler, S.; Bégué, A. A Combined Random Forest and OBIA Classification Scheme for Mapping Smallholder Agriculture at Different Nomenclature Levels Using Multisource Data (Simulated Sentinel-2 Time Series, VHRS and DEM). Remote Sens. 2017, 9, 259. [Google Scholar] [CrossRef] [Green Version]
  50. Song, Q.; Hu, Q.; Zhou, Q.; Hovis, C.; Xiang, M.; Tang, H.; Wu, W. In-Season Crop Mapping with GF-1/WFV Data by Combining Object-Based Image Analysis and Random Forest. Remote Sens. 2017, 9, 1184. [Google Scholar] [CrossRef] [Green Version]
  51. Phiri, D.; Morgenroth, J.; Xu, C.; Hermosilla, T. Effects of Pre-Processing Methods on Landsat OLI-8 Land Cover Classification Using OBIA and Random Forests Classifier. Int. J. Appl. Earth Obs. Geoinf. 2018, 73, 170–178. [Google Scholar] [CrossRef]
  52. Breiman, L. Random Forests. Mach. Learn. 2001, 45, 5–32. [Google Scholar] [CrossRef] [Green Version]
  53. Pal, M. Random Forest Classifier for Remote Sensing Classification. Int. J. Remote Sens. 2005, 26, 217–222. [Google Scholar] [CrossRef]
  54. Rokach, L. Decision Forest: Twenty Years of Research. Inf. Fusion 2016, 27, 111–125. [Google Scholar] [CrossRef]
  55. Olofsson, P.; Foody, G.M.; Herold, M.; Stehman, S.V.; Woodcock, C.E.; Wulder, M.A. Good Practices for Estimating Area and Assessing Accuracy of Land Change. Remote Sens. Environ. 2014, 148, 42–57. [Google Scholar] [CrossRef]
  56. Udali, A.; Lingua, E.; Persson, H.J. Assessing Forest Type and Tree Species Classification Using Sentinel-1 C-Band SAR Data in Southern Sweden. Remote Sens. 2021, 13, 3237. [Google Scholar] [CrossRef]
  57. Kendall, M.G. Rank Correlation Methods; Griffin: London, UK, 1948. [Google Scholar]
  58. Xiang, K.; Zhao, A.; Liu, H.; Zhang, X.; Zhang, A.; Tian, X.; Jin, Z. Spatiotemporal Evolution and Coupling Pattern Analysis of Urbanization and Ecological Environmental Quality of the Chinese Loess Plateau. Sustainability 2022, 14, 7236. [Google Scholar] [CrossRef]
  59. Estoque, R.C.; Murayama, Y.; Myint, S.W. Effects of Landscape Composition and Pattern on Land Surface Temperature: An Urban Heat Island Study in the Megacities of Southeast Asia. Sci. Total Environ. 2017, 577, 349–359. [Google Scholar] [CrossRef]
  60. People’s Government of Lhasa Municipality. Overall Urban Plan for Lhasa(2009–2020). 2017. Available online: https://www.docin.com/p-1894397395.html (accessed on 18 March 2023).
  61. People’s Government of Lhasa Municipality. Regulations on the Protection of Old Town of Lhasa. 2013. Available online: https://flk.npc.gov.cn/detail2.html?NDAyOGFiY2M2MTI3Nzc5MzAxNjEyODFjZGZhZTVlMTM%3D (accessed on 18 March 2023).
  62. People’s Government of Lhasa Municipality. Regulations on the Protection of Ancient Village of Lhasa. 2016. Available online: https://flk.npc.gov.cn/detail2.html?ZmY4MDgwODE3MWU5ZTE4MTAxNzIwMTljODU1MjJkMWE (accessed on 18 March 2023).
  63. People’s Government of Lhasa Municipality. Historic and Cultural City Protection Plan for Lhasa. 2019. Available online: http://www.xizang.gov.cn/rsxz/msgj/201912/t20191218_126510.html (accessed on 18 March 2023).
  64. Li, G.; Sun, S.; Fang, C. The Varying Driving Forces of Urban Expansion in China: Insights from a Spatial-Temporal Analysis. Landsc. Urban Plan. 2018, 174, 63–77. [Google Scholar] [CrossRef]
  65. Shahtahmassebi, A.R.; Song, J.; Zheng, Q.; Blackburn, G.A.; Wang, K.; Huang, L.Y.; Pan, Y.; Moore, N.; Shahtahmassebi, G.; Sadrabadi Haghighi, R.; et al. Remote Sensing of Impervious Surface Growth: A Framework for Quantifying Urban Expansion and Re-Densification Mechanisms. Int. J. Appl. Earth Obs. Geoinf. 2016, 46, 94–112. [Google Scholar] [CrossRef] [Green Version]
  66. Ma, Q.; He, C.; Wu, J. Behind the Rapid Expansion of Urban Impervious Surfaces in China: Major Influencing Factors Revealed by a Hierarchical Multiscale Analysis. Land Use Policy 2016, 59, 434–445. [Google Scholar] [CrossRef]
  67. Ma, Y.; Xu, R. Remote Sensing Monitoring and Driving Force Analysis of Urban Expansion in Guangzhou City, China. Habitat Int. 2010, 34, 228–235. [Google Scholar] [CrossRef]
  68. Chen, T.; Lang, W.; Chan, E.; Philipp, C.H. Lhasa: Urbanising China in the Frontier Regions. Cities 2018, 74, 343–353. [Google Scholar] [CrossRef]
  69. Zhou, P.; Xie, Y.; Gao, Z.; Jiang, G. Changes and Driving Forces of Impervious Surfaces in Xi’an City. Bull. Soil Water Conserv. 2020, 40, 274–281. [Google Scholar] [CrossRef]
  70. Tibet Statistical Bureau. Tibet Statistical Yearbook; China Statistics Press: Beijing, China, 2022; ISBN 9787503796791.
  71. Tibet Statistical Bureau. Tibet Statistical Yearbook; China Statistics Press: Beijing, China, 2015; ISBN 9787503774768.
  72. The State Council Information Office of the People’s Republic of China. Ecological Progress on the Qinghai-Tibet Plateau; The State Council Information Office of the People’s Republic of China: Beijing, China, 2018. Available online: http://english.www.gov.cn/archive/white_paper/2018/07/18/content_281476227186598.htm (accessed on 18 March 2023).
Figure 1. Study area.
Figure 1. Study area.
Sustainability 15 05510 g001
Figure 2. Workflow of this study.
Figure 2. Workflow of this study.
Sustainability 15 05510 g002
Figure 3. Impervious surface maps of (a) 2014 and (b) 2021.
Figure 3. Impervious surface maps of (a) 2014 and (b) 2021.
Sustainability 15 05510 g003
Figure 4. Standard deviational ellipses (SDEs) of impervious surfaces of 2014 and 2021.
Figure 4. Standard deviational ellipses (SDEs) of impervious surfaces of 2014 and 2021.
Sustainability 15 05510 g004
Figure 5. Landscape indices of impervious surface of different zones in 2014 and 2021.
Figure 5. Landscape indices of impervious surface of different zones in 2014 and 2021.
Sustainability 15 05510 g005
Table 1. Remote sensing images of Lhasa.
Table 1. Remote sensing images of Lhasa.
Acquisition DateCloud CoverCenter Latitude and LongitudeSensor
12 February 2014091.0 E, 29.7 NPMS
12 February 20141%91.3 E, 29.6 NPMS
5 February 20211%91.1 E, 29.7 NPMS
5 February 20211%91.1 E, 29.4 NPMS
Table 2. Landscape indices used in the study.
Table 2. Landscape indices used in the study.
Landscape IndicesCalculation and Description
Number of patches (NP) NP = n i
where ni is the number of patches in the landscape of patch type.
Mean patch size (MPS) MPS = i = 1 n a i / n
where n refers to the number of patches and ai is the area of patch i.
Largest patch index (LPI) LPI = max a i j / A × 100
where aij refers to the area (m2) of patch ij.
Landscape shape index (LSI) LSI = 0.25 E / A
where E is the total length (m) of the edge between patch types i and k in the landscape; this includes the entire landscape boundary and some or all background edge segments that involve patch types i.
Aggregation index (AI) AI = i = 1 n g i i / max g i i P i × 100
where gii is the number of like adjacencies (joins) between pixels of patch type i based on the single-count method and max gii refers to the maximum number of like adjacencies (joins) between pixels of patch type i based on the single-count method. Pi is the proportion of landscape comprising patch type i.
Table 3. Confusion matrix for the impervious surface mapping.
Table 3. Confusion matrix for the impervious surface mapping.
2014 Google Earth Images
Classified resultsClassificationImpervious surfaceNon-impervious surfaceSumUser’s accuracy
Impervious surface92910191.09%
Non-impervious surface839139997.99%
SUM100266500
Producer’s accuracy92.00%97.75%
OA96.60%
KAPPA89.41%
2021 Google Earth Images
Classified resultsClassificationImpervious surfaceNon-impervious surfaceSumUser’s accuracy
Impervious surface1251213791.24%
Non-impervious surface935436397.52%
SUM134366500
Producer’s accuracy93.28%96.72%
OA95.80%
KAPPA89.37%
Table 4. Impervious surface area changes in different zones of 2014 and 2021.
Table 4. Impervious surface area changes in different zones of 2014 and 2021.
ZoneImpervious Surface Area (km2)Impervious Surface CoverageIncrease Rate
2014202120142021
Central City zone51.14963.2999.975%12.345%23.755%
Historic zone4.5064.64265.575%67.554%3.017%
Environmental Coordination zone23.82724.97355.800%58.485%4.812%
Bakuojie zone1.1501.19983.967%87.532%4.246%
Caiyicun zone0.3280.40467.434%83.050%23.157%
Table 5. Parameters of the standard deviational ellipses (SDE) of 2014 and 2021.
Table 5. Parameters of the standard deviational ellipses (SDE) of 2014 and 2021.
YearLong Axis (m)Short Axis (m)Orientation (°)Long Axis/Short Axis
20148108.1763732.33393.6342.172
20218493.8194009.02094.0532.119
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wang, S.; Tan, X.; Fan, F. Changes in Impervious Surfaces in Lhasa City, a Historical City on the Qinghai–Tibet Plateau. Sustainability 2023, 15, 5510. https://0-doi-org.brum.beds.ac.uk/10.3390/su15065510

AMA Style

Wang S, Tan X, Fan F. Changes in Impervious Surfaces in Lhasa City, a Historical City on the Qinghai–Tibet Plateau. Sustainability. 2023; 15(6):5510. https://0-doi-org.brum.beds.ac.uk/10.3390/su15065510

Chicago/Turabian Style

Wang, Sishi, Xin Tan, and Fenglei Fan. 2023. "Changes in Impervious Surfaces in Lhasa City, a Historical City on the Qinghai–Tibet Plateau" Sustainability 15, no. 6: 5510. https://0-doi-org.brum.beds.ac.uk/10.3390/su15065510

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