Next Article in Journal
Multiscale Market Integration and Nonlinear Granger Causality between Natural Gas Futures and Physical Markets
Next Article in Special Issue
Role of Urban Public Space and the Surrounding Environment in Promoting Sustainable Development from the Lens of Social Media
Previous Article in Journal
Management Scale Assessment of Practices to Mitigate Cattle Microbial Water Quality Impairments of Coastal Waters
Previous Article in Special Issue
Remote Sensing-Based Analysis of Landscape Pattern Evolution in Industrial Rural Areas: A Case of Southern Jiangsu, China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Impacts of Landscape Changes on Annual Mean Land Surface Temperature in the Tropical Mountain City of Sri Lanka: A Case Study of Nuwara Eliya (1996–2017)

1
Faculty of Life and Environmental Sciences, University of Tsukuba, 1-1-1, Tennodai, Tsukuba, Ibaraki 305-8572, Japan
2
Department of Environmental Management, Faculty of Social Sciences and Humanities, Rajarata University of Sri Lanka, Mihintale 50300, Sri Lanka
3
Graduate School of Life and Environmental Sciences, University of Tsukuba, 1-1-1, Tennodai, Tsukuba, Ibaraki 305-8572, Japan
4
Department of Plant and Environmental Sciences, School of Natural Resources, Copperbelt University, P.O. Box 21692, Kitwe 10101, Zambia
*
Authors to whom correspondence should be addressed.
Sustainability 2019, 11(19), 5517; https://0-doi-org.brum.beds.ac.uk/10.3390/su11195517
Submission received: 23 August 2019 / Revised: 1 October 2019 / Accepted: 1 October 2019 / Published: 6 October 2019

Abstract

:
Although urbanization has contributed to improving living conditions, it has had negative impacts on the natural environment in urbanized areas. Urbanization has changed the urban landscape and resulted in increasing land surface temperature (LST). Thus, studies related to LST in various urban environments have become popular. However, there are few LST studies focusing on mountain landscapes (i.e., hill stations). Therefore, this study investigated the changes in the landscape and their impacts on LST intensity (LSTI) in the tropical mountain city of Nuwara Eliya, Sri Lanka. The study utilized annual median temperatures extracted from Landsat data collected from 1996 to 2017 based on the Google Earth Engine (GEE) interface. The fractions of built-up (BL), forested (FL) and agricultural (AL) land, were calculated using land use and cover maps based on urban–rural zone (URZ) analysis. The urban–rural margin was demarcated based on the fractions of BL (<10%), and LSTI that were measured using the mean LST difference in the urban–rural zone. Besides, the mixture of land-use types was calculated using the AL/FL and BL/FL fraction ratios, and grid-based density analysis. The results revealed that the BL in all URZs rapidly developed, while AL decreased during the period 1996 to 2017. There was a minimal change in the forest area of the Nuwara Eliya owing to the government’s forest preservation policies. The mean temperature of the study area increased by 2.1 °C from 1996 to 2017. The magnitude of mean LST between urban–rural zones also increased from 1.0 °C (1996) to 3.5 °C (2017). The results also showed that mean LST was positively correlated with the increase and decrease of the BL/FL and AL/FL fraction ratios, respectively. The grid-based analysis showed an increasing, positive relationship between mean LST and density of BL. This indicated that BL density had been a crucial element in increasing LST in the study area. The results of this study will be a useful indicator to introduce improved landscape and urban planning in the future to minimize the negative impact of LST on urban sustainability.

1. Introduction

In recent decades, population growth and economic development have directly affected landscapes’ transformations in developing countries [1]. Rapid changes in the landscape have resulted in the conversion of natural vegetation and agricultural land into built-up (impervious) land, such as buildings, parking lots, roads, and other constructions [1,2,3,4]. This has caused several environmental problems at local, regional, and global scales [5], such as decreases in agricultural land [6]; habitat destruction [5,7,8]; air, soil, and water contamination [9,10,11]; increases in vector-borne diseases, such as malaria and dengue [12]; decreases in green space [13,14]; and increases in land surface temperature (LST) [15,16,17].
Increasing temperatures in urban environments are largely an outcome of rapid urbanization and anthropogenic activities [18]. Built-up areas have been formed using several materials, such as concrete, flooring, pebbles, stone, and gravel, which decrease evapotranspiration, increase the sensitivity of the city and notably affect its local climate [16,19,20]. However, studying temperature changes based on the air temperature is challenging, as there is a lack of meteorological stations, especially in developing countries [21]. Thus, satellite remote-sensing data provides vital information for observing the temperature patterns in urban areas [21]. Several studies have been conducted by using vast range of remote sensing data, such as Modis data [22,23,24,25], Synthetic Aperture Radar (SAR) data [26], Nightlights data [27,28], Landsat data [17,29,30], Land Scan data [31,32,33], and fossil fuel CO2 emission data [33,34], to understand LST patterns. Many urban landscapes from small to large scales have been studied worldwide [16,35,36], including coastal cities [16,33,37,38], desert cities [20,39], and mountain cities [13,14]. Mountain cities are attractive for rich people, as they have a cold climate and comfortable living conditions, which have resulted in rapid urban development [14]. Thus, studies focusing on mountain landscapes are vital for understanding the changing pattern of LST to introduce mitigation measures for comfortable living conditions.
A large number of studies have used two or more satellite images from different time points to analyze an LST pattern due to the unavailability of cloud-free images. However, the difference in the acquisition times might influence any resulting LST pattern due to varying environmental factors (wind speed, the Sun’s radiation, surface moisture, and humidity) [4,32]. The use of more satellite images captured in multiple time points can potentially provide more specific information to understand the changing pattern of LST [4]. Still, it is not easy to analyze extensive earth observation data sets, because of issues of spatial and temporal resolution [40,41]. Thus, big data analysis platforms can be used as an alternative for conducting accurate results [40]. The Google Earth Engine (GEE) provides the potential to process a large number of satellite images and researchers can easily access free public data archives for more than thirty years of historical data [40]. Hence, in this study, we used GEE to extract the annual median LST for three-time points such as 1996, 2006, and 2017 based on several images captured during the selected years. We hypothesized that the use of many images captured in multiple time points would provide a more precise pattern of LST in the study area.
The spatial distribution of the LST intensity (LSTI) provides essential environmental information for understanding a temperature pattern in detail. The literature provides two methods that can be used to study LSTI. The first one is the categorization of land use and cover as the local climate zone and then following the cross cover comparison method to calculate the LSTI [13,14,42,43,44]. The second method involves determining the difference in the mean LST between urban and rural zones based on the urban–rural gradient analysis [13,14,45]. In this study, we used the second method. Other studies utilizing the urban–rural gradient analysis have shown that the urban–rural demarcation is essential for determining the temperature differences [13,14].
Mountain cities in Asia have been developing since the colonization period of the 19th and 20th centuries [46,47]. The cool climate and natural landscapes became the most prominent factors driving the development of mountain cities. During the colonial period, cool climates were preferred, as they allowed the colonials to maintain their “western lifestyle.” Mountain cities were also selected to avoid wasting illnesses, episodic pestilence, sunstroke, and depression [46,48]. As a result, most mountain cities have been experiencing rapid urbanization since the end of World War II until now. Notwithstanding the vast literature on the impacts of urbanization on landscape changes and other associated consequences, like increasing LST, mountain cities remain one of the landscapes that have not been extensively studied.
Nuwara Eliya is a historical city that developed during the British colonial period, as reflected in the city’s architecture [49]. The British people preferred a cold climate and Nuwara Eliya was referred to as “Little England” [49]. They used Nuwara Eliya as a meeting place for wealthy British families [50]. Currently, several recreational sites, such as Lake Gregory, golf links, racecourses, and a large number of clubs, remain in Nuwara Eliya [50]. After gaining independence in 1948, Nuwara Eliya became renowned worldwide as a tourist destination in Sri Lanka, and thousands of tourists from Sri Lanka and overseas visited this area to enjoy the cold climate and natural beauty [49]. This has resulted in several environmental problems associated with urbanization in Nuwara Eliya. One of the critical environmental impacts of urbanization that has likely affected the mountain city Nuwara Eliya is the Urban Heat Island (UHI) resulting from variations in LST [13,14]. Still, previous studies related to LST have been focused on the Colombo City [3,4,37] and Kandy [13,51] in Sri Lanka.
Thus, in this study, we hypothesized that the rapid changes in the urban landscape of Nuwara Eliya influenced the patterns of LST. The study examined the changes in the urban landscape and their impacts on the spatiotemporal changes of LST, in order to contribute to efforts aimed at introducing proper landscape and urban planning in Nuwara Eliya. The study has three objectives to: (1) monitor urban landscape changes and their impacts on the spatial and temporal variations of LST intensity from 1996–2017 based on the urban–rural gradient analysis; (2) identify relationships between mean LST and the agricultural land (AL)/forested land (FL) and built-up land (BL)/FL fraction ratios; and (3) study the relationship between mean LST and densities of BL, FL, and AL based on the grid-based analysis method. The results of this study could be useful for enhancing capacity and knowledge to minimize the possible negative impacts of rapid urbanization in the study area.

2. Materials and Methods

2.1. Study Area: Nuwara Eliya, Sri Lanka

Nuwara Eliya is located in the central province of Sri Lanka and resides at latitudes of 6°54′21.94″ to 7°2′28.53″ N and longitudes of 80°50′5.89″ to 80°41′57.23″ E (Figure 1). The study area includes the landscape within a 7.5 km2 radius of the center of Nuwara Eliya. Nuwara Eliya is surrounded by one of the tallest mountains (~2524 m high) in the study area, known as Pidurutalagala. The average elevation of Nuwara Eliya is about 1800 m. Monthly rainfall in Nuwara Eliya ranges from approximately 71.5 to 226.8 mm, with an annual average of 1905 mm. The lowest average monthly rainfall is recorded in January, February, and March, resulting in a relatively dry period. According to the Meteorology Department of Sri Lanka, the daily mean temperature is about 15.9 °C, with average minimum and maximum temperatures of 11.6 °C and 20.2 °C, respectively. Nuwara Eliya is the most popular site in Sri Lanka among local and foreign tourists due to the cold climate that persists throughout the year. In terms of the urban development pattern, Nuwara Eliya exhibits a single locus.

2.2. Overall Workflow

Figure 2 shows the overall workflow of the study to achieve the objectives described above. The study workflow included five major steps: (i) extraction of median upper-atmosphere brightness temperature from the thermal bands of Landsat images (Figure 2a); (ii) retrieval of LST (Figure 2b); (iii) land use/land cover (LULC) classification into five classes, forest land (FL), built-up (BL), agriculture land (AL), other land, and water using machine learning algorithms, based on supervised classifications (Figure 2c); (iv) urban–rural gradient, statistical, and intensity analyses, based on the BL, FL, and AL fractions, AL/FL and BL/FL fraction ratios, and mean LST; and finally (v) grid analysis based on the mean LSTs and the densities of BL, FL, and AL.

2.3. Calculating Annual Median At-Satellite Brightness Temperature Using the Google Earth Engine (GEE)

The study employed the GEE to calculate annual median at-satellite brightness temperatures using atmospherically corrected pre-processed Landsat datasets (Level 2) [52]. We hypothesized that the use of more images would produce clearer output images to understand the LST patterns in the study area. Past studies have shown the aptness of this method to generate median LSTs in different areas [40,41]. In this process, several steps were performed as follows.
(i) The study area was defined and imported as an “Asset” in GEE, and subsequently used as the primary geometry source throughout the process. Masking was then conducted due to cloud disturbance in the available Landsat imageries. The cloud disturbance could be attributed to Nuwara Eliya being located in a tropical area.
(ii) The Image Collection tool in GEE was used to prepare the imagery for the study, including 15 images from 1996 (Landsat 5), 17 images from 2006 (Landsat 5), and 20 images from 2017 (Landsat 8) (Table A1).
(iii) Afterwards, median temperatures of each pixel were extracted for 1996, 2006, and 2017 by using “ee.reducer” methods available in GEE [53]. Extraction of the pixel based median temperatures was done to take care of inaccurate extraction caused by disturbance from clouds in the study area. Figure 3 shows the graphical illustration of creating image collection and the code used to generate the annual, median, at-satellite brightness temperature for each pixel of each year (provided in Annex Code A1 (Landsat 5TM) and A2 (Landsat 8)).
(iv) The extracted median temperatures based on the upper-atmosphere brightness temperatures from the thermal bands (in Kelvin) (band 6 for Landsat 5 and band 10 for Landsat 8) were used to calculate the LSTs, as described in Section 2.4.
(v) All extracted images were georectified using the WGS84/UTM 44N projection system before further processing.

2.4. LST Calculations

The extracted annual, median, at-satellite brightness temperatures, as mentioned in Section 2.3, were scaled using the land surface emissivity derived from Equation (1) [54].
ε = { mPV + n }
where m = (ε − ε) − (1 − εσ); Fεv and n = εs + (1 − εs) Fεv. εs and εv are the soil emissivity and vegetation emissivity, respectively. In this study, we used the result of [54] for m = 0.004 and n = 0.986. The proportion of vegetation (Pv) was calculated using the normalized difference vegetation index (NDVI) based on Equation (2) [55].
NDVI = ρ N I R ρ R e d ρ N I R + ρ R e d
where ρNIR refers to the surface reflectance values of bands 4 (Landsat-5) and 5 (Landsat-8); and ρRed refers to the surface reflectance values of bands 3 (Landsat-5), and 4 (Landsat-8 OLI).
Pv was extracted using Equation (3).
P v = ( ( NDVI NDVI min ) / ( NDVI max   NDVI min ) ) 2
where, Pv is the proportion of vegetation, NDVI is original NDVI value calculated using Equation (2), and NDVImin and NDVImax are the minimum and maximum values of the NDVI dataset, respectively.
The emissivity corrected images were used to extract LSTs using Equation (4) [3,30].
LST = T b / 1 + ( λ   × T b / ρ ) ln ε
where Tb is the at-satellite brightness temperature in degrees Kelvin; λ is the central band wavelength of emitted radiance (11.5 μm for band 6 [30] and 10.8 μm for band 10 [13]; ρ is h × c/σ (1.438 × 10–2 m K), with σ as the Boltzmann constant (1.38 × 10–23 J/K), h as Planck’s constant (6.626 × 10–34 J *s), and c as the velocity of light (2.998 × 108 m/s); and ε is the land-surface emissivity estimated using the NDVI method [54]. The calculated LST value (Kelvin) is then converted to °C.

2.5. Land Use/Land Cover (LULC) Mapping

Machine learning methods, such as support vector machines (SVM), K-nearest neighbor, random forest, and neural networks, have been widely used to classify LULC [13,14]. Among these methods, the SVM has provided higher overall accuracy [56,57,58,59]. Thus, in this study, we used the SVM algorithm to conduct LULC mapping. The classification scheme used in the study included five LULC categories; namely, built-up land (BL), forested land (FL), agricultural land (AL), other lands, and water. The other lands category was comprised of a combination of grasslands and bare lands.
We produced three LULC maps containing the five LULC categories (BL, FL, AL, other land and water) for the years 1996, 2006, and 2017 with overall accuracies of 85%, 93%, and 92%, respectively (see Table A2, Table A3 and Table A4). The accuracy was assessed by using 500 reference points generated by a stratified random sampling method [14] for all the LULC categories. Google Earth imageries were used to assess the accuracy of the classified LULC maps for 2006 and 2017. The accuracy assessment for the 1996 LULC map was conducted with the aid of available topographic maps and different band combinations of Landsat imageries.

2.6. LST Intensity (LSTI) Measurement

The LSTI was calculated based on the urban–rural gradient analysis approach, which involves the creation of concentric rings or buffer zones around the city center with standard distance intervals extending to the rural areas [13,14,42]. Urban–rural gradient analysis has been conducted to identify the spatial and temporal variations in environmental variables in many previous environmental studies [13,14,16,20,60,61]. In this study, thirty-five 210 m buffer zones (hereafter referred to as urban–rural zones (URZs)) were created for the study area. The mean LST of each URZ was then extracted using zonal statistics. The fractions of BL, FL, and AL were determined by calculating their respective proportions in each URZ. Of note is that previous studies combined AL and FL into one LULC category defined as green space [13,14,16]. However, AL has been one of the key drivers of LULC changes in the study area. Therefore, AL in Nuwara Eliya and its surroundings was considered separately.
The magnitude of the LSTI along the urban–rural gradient (LSTI U–R) was determined based on the ∆ mean LST, ∆ fraction of BL, ∆ fraction of FL, and ∆ fraction of AL, following methodology proposed by Estoque and Murayama [14]. To calculate LSTIU–R, we first determined the ∆ mean LST by finding the difference between the mean LST in the URZ with the highest fractions of BL, FL and AL (defined as URZ1) and other URZs (i.e., URZ1 − URZ2…35). We then applied the same method to determine the ∆ fraction of BL, ∆ fraction of FL, and ∆ fraction of AL along the urban–rural gradient. We used the same threshold for delineating the urban and rural zones as the authors of [13,14]; i.e., from the city center, all URZs with an >10% fraction of BL were considered as urban and those beyond the first URZ with <10% fraction of BL were considered as rural.

2.7. AL/FL and BL/FL Fraction Ratios and Their Intensities

The AL/FL and BL/FL fraction ratios were calculated using the URZs created in Section 2.6. Ranagalage et al., 2018 [13] proposed the green space (GS)/impervious surface (IS) fraction ratio and its intensity, and we followed their methodology to extract AL/FL and BL/FL in each URZ using Equations (5) and (6), respectively.
AL / FL   fraction   ratio = ALz FLz
BL / FL   fraction   ratio = BLz FLz

2.8. Grid-Based Analysis

The grid-based method was used to analyze the spatial distributions of the densities of BL, FL, AL, and mean LSTs in 1996, 2006, and 2017. In this analysis, we used 210m × 210 m grids (7 × 7 pixels) to demarcate the relationships between BL, FL, and AL density with mean LST. The 210 × 210 m grid size was used in previous studies and showed a high correlation with mean LST [16,39,44,62]. After the creation of the set of grids, the mean LST and the densities of BL, FL, and AL in each grid were calculated. Then, scatter plots were created and linear regression analysis was performed to identify the relationships between mean LST and densities of BL, FL, and AL.

3. Results

3.1. Landscape Changes and LST Distribution of Nuwara Eliya

The results revealed that Nuwara Eliya experienced rapid urbanization in the last 21 years. The built-up area increased from 289.9 ha to 2,080.4 ha from 1996 to 2017 with an annual growth rate of 85.3 ha per year (Figure 4 and Table 1 and Table 2). The forest area did not change significantly due to the implementation of forest reserves by the government [63]. It was observed that the rapid changes in the built-up land negatively affected the agricultural sector. The area of agricultural land decreased from 8503.2 ha to 6583.9 ha from 1996 to 2017 (Table 1).
Figure 4d–f shows the spatial distribution of LST in Nuwara Eliya for the years 1996, 2006, and 2017, respectively. In 1996, LST ranged between 9.1 and 29.2 °C with a mean value of 18.9 °C. In 2006, LST ranged between 5.9 and 30.2 °C with a mean value of 17.9 °C. In 2017, LST ranged between 14.3 and 31.0 °C with a mean value of 21.0 °C. Therefore, the mean LST exhibited an increasing trend from 1996 to 2017. The mean LST increased by 2.1 °C during the past 21 years. Figure 4d–f shows increasing temperatures in the city center and in the southwestern and northeastern areas. The observed pattern of increasing temperatures was similar to the increase in BL.

3.2. Magnitude and Trend of LSTI

3.2.1. LSTIU–R along the Urban–Rural Gradient

Figure 5a–c shows the spatial pattern of mean LSTs and the fractions of BL, FL, and AL along the urban–rural gradient. The highest mean LST value was recorded in the URZ1 near the city center. In 1996, the mean LST of URZ1 was 19.8 °C which decreased to 18.9 °C in 2006 and increased to 24.4 °C in 2017. Conversely, the mean LST of all other URZs increased from 1996 to 2017. The results revealed that the mean LST of all other URZs ranged from 18.9 °C to 21.2 °C between 1996 and 2017. The results further revealed that the lowest temperatures were recorded in URZ12 in 1996, URZ10 in 2006, and URZ19 in 2017 (Figure 5).
The fraction of BL increased rapidly between 1996 and 2017. In URZ1, the fraction of BL increased by 32.4%, 48.7%, and 58.7%, in 1996, 2006, and 2017, respectively. The fraction of FL did not significantly change during the last 21 years. However, rapid changes in AL occurred in all URZs throughout the study’s temporal extent. The lowest fraction of AL was continually recorded in URZ1 and exhibited a decreasing pattern. The AL fraction in URZ1 was 43.6% in 1996, 29.6% in 2006, and 15.8% in 2017.
Figure 5b shows the results of the statistical analysis of mean LST with the fractions of BL, FL, and AL. The mean LST had a strongly significant, positive correlation (p < 0.001) with the fraction of BL at all the three time-points. The fraction of FL exhibited a strongly significant, negative correlation (p < 0.001) with mean LST. On the other hand, mean LSTs exhibited dissimilar relationships with the fractions of AL in 1996, 2006, and 2017. In 1996, the mean LST had a significant, positive correlation (p < 0.001) with the fraction of AL. In 2006, the mean LST and the fraction of AL had a weak relationship (p < 0.5), and in 2017, the mean LST had a significant, negative correlation with the fraction of AL (p < 0.1).

3.2.2. The Magnitude of LSTIU–R along the Urban–Rural Gradient

The magnitudes of the mean LST and fractions of BL, FL, and AL were calculated based on URZ1 as described in Section 2.6 above. The ∆ mean LST in URZ13 showed a considerable decrease in 1996 and 2006. However, the decrease shifted to URZ20 in 2017 (Figure 6a). The ∆ fractions of BL, FL, and AL also showed a similar pattern. The ∆ mean LST had a significant relationship (positive) with the ∆ fraction of BL, and a strong relationship (negative) with the ∆ fraction of FL. The ∆ mean LST exhibited a significant, positive relationship in 1996 (p < 0.001), and in 2006 (p < 0.5), with the ∆ fraction of AL. However, the relationship changed in 2017, showing a weak, negative correlation (p < 0.5) with the ∆ fraction of AL.

3.3. AL/FL and BL/FL Ratios

An annual decrease in AL of 91.6 ha was observed from 1996 to 2017 (Table 2). FL was stable due to the protection measures implemented by the Forest Department of Sri Lanka (Table 1). A rapid decreasing trend in the AL/FL fraction ratio was observed from 1996 to 2017. The AL/FL fraction ratio decreased in all zones from URZ1 to URZ35 with average values of 0.84, 0.7, and 0.53 in 1996, 2006, and 2017, respectively (Figure 7a). The mean LST values had a significant (p < 0.001), positive relationship with the AL/FL fraction ratio. Figure 7c shows the statistical relationship between the change in mean LST and Δ(AL/FL) fraction ratio along the urban–rural gradient, and the result shows a significant positive relationship in all the three time-points in this study.
The BL/FL ratio increased from 1996 to 2017, with average values of 0.13, 0.22, and 0.42 in 1996, 2006, and 2017, respectively (Figure 8a). The mean LST exhibited a significant (p < 0.001), positive relationship with the BL/FL fraction ratio (Figure 8b). The relationship between the change in mean LST and Δ(BL/FL) fraction ratio is presented in Figure 8c. The results revealed that there was a strongly significant, positive relationship between the change in mean LST and Δ(BL/FL) fraction ratio for all the three time-points.

3.4. The Density of BL, FL, and AL Compared to Mean LST

Figure 9a shows the spatial distribution of mean LST and the density of BL, FL, and AL in 1996, 2006, and 2017. High temperatures were observed around the city center and the northwestern and southwestern parts of the study area. The BL density was also high around the city center area in 1996 and it spread in the southern, southwestern and southeastern directions within the Nuwara Eliya area in 2006 and 2017. The FL density had an opposite pattern, showing low FL densities around all areas of high BL densities in all the three-time points. However, AL density showed a decreasing trend, especially in the city center area and southwestern direction of the study area.
Figure 9b shows the relationship between the densities of BL, FL, AL, and mean LST based on the grid-based analysis. The mean LST was statistically significant with the densities of BL, FL, and AL in all the three-time points (p < 0.001). There was an increasing trend of a positive correlation between BL density and mean LST from 1996 to 2017. The correlation between mean LST and the density of FL showed a significant negative relationship throughout the study period. In addition to that, there was a decreasing trend of a positive correlation between AL density and mean LST from 1996 to 2017.

4. Discussion

In this study, we examined the landscape changes and their impacts on LST in Nuwara Eliya, Sri Lanka over a 21-year period (1996 to 2017). The results revealed that the urbanization pattern of Nuwara Eliya was concentrated around the city center areas in 1996 which then expanded towards the southern, southwestern, northwestern, and eastern parts by 2017. The development pattern of Nuwara Eliya is similar to that of other mountain cities (e.g., Kandy) in Sri Lanka [13]. The built-up land has expanded rapidly in the past 21 years (Table 1 and Table 2, and Figure 4) with an overall increase of 1790.6 ha and an annual increase of 85.3 ha. The tourism industry has been reported to be the main driving force of the observed rapid development in the study area, as thousands of local and foreign tourists have been visiting Nuwara Eliya [49]. Previous studies, such as Baguio in the Philippines [14,46], Kandy in Sri Lanka [13,51], Bengaluru in India [64], Bogor in Indonesia [46], and Dalat in Vietnam [46], revealed that green spaces, such as forests, declined due to the urban development pressure in their respective study areas. However, our study revealed that the forest areas were intact throughout the study period, owing to strong government policies on forest preservation. Most of the forest areas located in Nuwara Eliya are protected. That has resulted in other LULC types other than forests being consumed. Our results demonstrate that BL expanded mainly at the expense of agricultural areas (Figure 5a).
The results revealed that the fraction of BL greatly influenced the LSTI in Nuwara Eliya. During the past 21 years, annual mean LST increased by 2.1 °C (18.9 °C in 1996 and 21.0 °C in 2017). Additionally, the mean annual temperature difference between the urban (URZ1) and rural (<10% BL) zones showed an increasing pattern, recording values of 1.0 °C (19.8 °C in URZ1 and 18.8 °C in URZ5) in 1996, 1.3 °C in 2006 (18.9 °C in URZ1 and 17.6 °C in URZ6), and 3.5 °C in 2017 (24.4 °C in URZ1 and 21.0 °C in URZ18) (Figure 5a). These results show an interesting pattern of LST in the study area. Urbanization has considerably influenced the increasing annual mean LST of the study area. The highest mean LST values were recorded around the city center (URZ1) in all the three-time points. The fractions of BL increased by 32.4%, 48.7%, and 58.7%, while the fraction of AL decreased by 43.6%, 29.6%, and 15.8% in 1996, 2006, and 2017, respectively. The mean LST had a significant, positive relationship with the fraction of BL and significant negative relationship with the fraction of FL (p < 0.001) in all the three time-periods (Figure 5b). These results are similar to previous research findings in other mountain cities in Asia, such as Kandy in Sri Lanka [13], and Baguio in the Philippines [14]. The mean LST of all URZs (i.e., urban and rural combined) also showed an increasing pattern throughout the temporal extent of the study (1996–2017). The patterns observed in Nuwara Eliya not only showed changes in the landscape, but also agreed with other researchers that have reported an increasing trend of global temperatures at high-elevation regions [65,66,67]. It is thus plausible to speculate that these global temperatures in high-elevation regions could have also contributed to the observed LST increases in the study area and other urban environments worldwide.
The statistical relationship between mean LST and fraction of AL was different in all the three time-points. The mean LST showed a strongly significant, positive relationship with the fraction of AL in 1996, while in 2006, the relationship was significant, but weakly positive. However, in 2017, the relationship became weakly negative. We observed that poorly managed agricultural lands were converted into built-up lands during the 21 years of this study. The average fraction of AL exhibited a decreasing pattern, with values of 38.7%, 35.8%, and 27.2% in 1996, 2006, and 2017, respectively. This result is similar to the observations in other mountain cities in Sri Lanka [13]. The observed pattern is critical evidence for policymakers to introduce proper mechanisms for minimizing the adverse impacts associated with the rapid changes in the urban landscape of the study area. The negative impacts of urban development on agricultural land could cause several environmental issues [59], such as increased soil erosion [68], and landslides [69,70,71]. The rapid development of built-up land, resulting in increased levels of LST, could also negatively affect the tourism industry. Policies should be implemented to maintain the natural beauty of Nuwara Eliya for the sustainable tourism industry. Most tourists wish to enjoy the landscape and comfortable climate of the area [49,59].
The change in the urban–rural pattern provides a clear picture of the urbanization pattern in the study area. The urban–rural zones were demarcated using the fraction of BL (<10%). The rural zone exhibited a shifting pattern from 1996 to 2017. The rural zone located in URZ5 (1 km from the city center) in 1996 and the rural zone located in URZ6 (1.2 km from the city center) in 2006, had shifted to URZ18 by 2017 (3.8 km from the city center) (Figure 6a). Rapid changes to built-up land resulted in a shift in the rural zone away from the city center. Previous researchers have shown a similar trend in the shifting rural zone in mountain cities [13,14]. Results showed that the temperature of URZ1 was higher than that of the rural zones. The change in mean LST had a strong, positive relationship with the ∆ fraction of BL and a strong negative relationship with the ∆ fraction FL. However, the relationship between the ∆ mean LST and ∆ fraction of AL shifted from positive (1996 and 2006) to negative (2017). Policymakers should pay much attention to the shifting of the urban zone.
The mixture of land-use types, which is important for controlling LST [72], was investigated based on two indices: the AL/FL and BL/FL fraction ratios. The results showed that the AL/FL fraction ratio decreased from 1996 to 2017 (Figure 7a). The maximum AL/FL fraction ratio declined with values of 3.3, 2.3, and 1.3 in 1996, 2006, and 2017, respectively. The mean LST exhibited a significantly strong positive relationship with the AL/FL fraction ratio in 1996 and 2006. However, the degree of the relationship became weaker in 2017 due to the changes in the AL of the study area. The change in mean LST exhibited a positive relationship with the ∆(AL/FL) fraction ratio over the three time-points (Figure 7c). On the other hand, the result of the BL/FL fraction ratio provided an indicator of the LSTI pattern in the study area. Changes in built-up land have been proven to be vital in understanding the LST changing pattern in urban areas [20,33,73]. The maximum BL/FL fraction ratio increased with values of 1.7, 2.3, and 3.4 (URZ1) in 1996, 2006, and 2017, respectively. This result could be the reason for the observed increase in the mean LST over the 21-year study period (Figure 8a). The increasing BL/FL fraction ratio directly affected the LSTI in the study area, as the mean LST had a strongly significant relationship with the BL/FL fraction ratio over the 21-year study period.
The spatial distribution pattern of the BL, FL, and AL densities provided vital information for understanding the distribution pattern of the study area. The linear regression analysis results also revealed the vital formation of capturing the land use mixture of the study area. The average density of BL increased with values of 1.3%, 3.5%, and 9.2% while the density of AL decreased with values 37.6%, 35.9%, and 29.2% in 1996, 2006, and 2017, respectively (Figure 9a). This pattern indicates that the BL density has been increasing in the last 21 years. The linear regression analysis results showed that the coefficient of determination (R2) between mean LST and the density of BL increased from 0.04 to 0.29 from 1996 to 2017. This shows that the density of BL had a significant influence on the increase of LST in the study area. The R2 between the mean LST and the density of AL declined from 0.59 to 0. 47 from 1996 to 2017. This also shows that AL density had a positive influence on the increase of LST. The density of FL, on the other hand, influenced the control of LST in the study area, showing a negative, strongly significant relationship between the two (Figure 9b). Most of the previous studies have combined AL as green spaces, but AL of the Nuwara Eliya area has been positively influencing to the increase of LST in the three-time points. The policymakers and planners need to consider maintaining a good mixture of land use to control the LST pattern in the future.
On the whole, the rapid increase of built-up areas has been the main force driving the high temperatures recorded in 2017 compared to 2006 and 1996. The ∆ mean LST exhibited a positive relationship with the ∆(BL/FL) fraction ratio over the three time-points (Figure 8c). Previous research has demonstrated that having a mixture of land-use types can reduce temperatures in urban areas [72]. Based on the results from this study, maintenance of forest coverage will continue at a similar level due to governmental preservation policies. However, rapidly developing urban areas require proper urban planning to maintain sustainability. We have noticed that the present urban development pattern is unplanned. The walls and roofs of most of the buildings are also green in color; something which could reduce the indoor temperatures [74,75]. To achieve sustainable development, we recommend that urban planners and policymakers should consider constructing green belts along the main road, schools, hospitals, and other government buildings. The development pattern observed was mainly horizontal rather than vertical. Vertical development could reduce pressure on land and leave a lot of free land for improving greenspace coverage. Vertical development can also be used as a proxy indicator to understand urban-development intensity [76,77]. In addition, greenspace cover plays a vital role in reducing the impacts of LST in urban areas [78]. We believe that the current urban planning must focus on protecting this popular tourist destination in Sri Lanka’s sustainably. The findings of this study can be used as an indicator for reorienting current urban planning policies to improve the natural and social environment in Nuwara Eliya.
Overall, this study successfully used annual mean LST to quantify landscape changes and their impacts on LSTI from 1996 to 2017. However, there are two main limitations of the study that should be considered. Firstly, comparing the mean LST change magnitudes with changes in observational air temperature measured from metrological stations could have helped in verifying our findings. However, that was challenging due to the unavailability of air temperature observational data from meteorological stations. Nonetheless, we used previous similar studies to validate our results, and we recommend that more studies should be conducted to confirm the increasing temperature in mountain cities of Sri Lanka. The lack of observational air temperature data in most developing countries has been the basis for the proliferation of UHI studies based on LST retrieved from remote sensing satellite imagery, compared to atmospheric UHI studies based on measured air temperature. Secondly, we believe that horizontal and vertical development of buildings both play a vital role in controlling the annual mean LST of the study area. However, vertical data on various types of buildings and built-structures were not used in this study due to unavailability in Sri Lanka. Thus, the results should be interpreted in the context of those limitations.

5. Conclusions

The spatial distribution of LST provides vital information for understanding the local climatic variation of cities and can be used as a proxy indicator to introduce sufficient measures to minimize the negative impacts associated with high LST. This study revealed that Nuwara Eliya had undergone rapid urban development over the 21 years studied. The study used multiple images to extract the median annual LSTs for 1996, 2006, and 2017. The results revealed that the annual mean LST increased by 2.1 °C over the past 21 years. The annual mean LST difference between the urban and rural zones increased from 1.0 °C to 3.5 °C from 1996 to 2017. The study shows that rapid development was spreading towards rural zones, and the fraction of BL influenced the increase in annual mean LST. The government policies have been protecting the FL, resulting in minimal changes which have been the major contributor to the control of the LST pattern. This study recommends that having a mixture of land-use types would considerably control the increasing LST in the study area. Thus, policymakers need to consider the importance of the land use mixture to reduce the impacts of a high LST in the study area. The future urban policy must focus on maintaining the natural splendor of the Nuwara Eliya to promote tourism in the study area and Sri Lanka. Nuwara Eliya is renowned as a tourist destination in Sri Lanka. Thus, the findings of this study can be used as an indicator to introduce sustainable future landscape and urban planning to protect the world-renowned tourist hub of Sri Lanka.

Author Contributions

The author to receive correspondence, M.R., proposed the topic and spearheaded the data processing and analysis, as well as the writing of the manuscript. Y.M., D.D., and M.S. helped in the design, research implementation, analysis, and the writing of the manuscript.

Funding

This study was supported by the Japan Society for the Promotion of Science (JSPS) through Grant-in-Aid for Scientific Research (B) 18H00763 (2018-20).

Acknowledgments

The authors are grateful to the anonymous reviewers for their helpful comments and suggestions to improve the quality of this paper.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Properties of the Landsat images (Level 2) used in this study.
Table A1. Properties of the Landsat images (Level 2) used in this study.
SensorScene IDAcquisition
Date
Time
(GMT)
Cloud Cover (%) in Landsat Title
Landsat 5 TMLT05_L1TP_141055_19960221_20170106_01_T11996-02-2103:59:4519
LT05_L1TP_141055_19960308_20170106_01_T11996-03-0804:00:496
LT05_L1TP_141055_19960324_20170105_01_T11996-03-2404:01:5111
LT05_L1TP_141055_19960409_20170105_01_T11996-04-0904:02:5146
LT05_L1TP_141055_19960425_20170104_01_T11996-04-2504:03:499
LT05_L1TP_141055_19960511_20170104_01_T11996-05-1104:04:469
LT05_L1TP_141055_19960527_20170104_01_T11996-05-2704:05:4129
LT05_L1TP_141055_19960730_20170103_01_T11996-07-3004:09:0678
LT05_L1TP_141055_19960815_20170103_01_T11996-08-1504:09:5666
LT05_L1TP_141055_19960831_20170103_01_T11996-08-3104:10:4841
LT05_L1TP_141055_19960916_20170102_01_T11996-09-1604:11:4170
LT05_L1TP_141055_19961002_20170102_01_T11996-10-0204:12:3348
LT05_L1TP_141055_19961103_20170102_01_T11996-11-0304:14:0938
LT05_L1TP_141055_19961119_20170101_01_T11996-11-1904:14:5319
LT05_L1TP_141055_19961205_20170101_01_T11996-12-0504:15:4084
Landsat 5 TMLT05_L1TP_141055_20060131_20161123_01_T12006-01-3104:44:0646
LT05_L1TP_141055_20060216_20161123_01_T12006-02-1604:44:2828
LT05_L1TP_141055_20060304_20161122_01_T12006-03-0404:44:4941
LT05_L1TP_141055_20060405_20161123_01_T12006-04-0504:45:2610
LT05_L1TP_141055_20060421_20161122_01_T12006-04-2104:45:4168
LT05_L1TP_141055_20060507_20161122_01_T12006-05-0704:45:5533
LT05_L1TP_141055_20060523_20161121_01_T12006-05-2304:46:0731
LT05_L1TP_141055_20060608_20161121_01_T12006-06-0804:46:2327
LT05_L1TP_141055_20060624_20161121_01_T12006-06-2404:46:3972
LT05_L1TP_141055_20060710_20161120_01_T12006-07-1004:46:5330
LT05_L1TP_141055_20060811_20161119_01_T12006-08-1104:47:1779
LT05_L1TP_141055_20060827_20161119_01_T12006-08-2704:47:2945
LT05_L1TP_141055_20060912_20161119_01_T12006-09-1204:47:4157
LT05_L1TP_141055_20060928_20161119_01_T12006-09-2804:47:5255
LT05_L1TP_141055_20061014_20161118_01_T12006-10-1404:48:0322
LT05_L1TP_141055_20061030_20161118_01_T12006-10-3004:48:1330
LT05_L1TP_141055_20061115_20161118_01_T12006-11-1504:48:2160
Landsat 8 OLI/TIRSLC08_L1TP_141055_20170113_20170311_01_T12017-01-1304:54:053
LC08_L1TP_141055_20170129_20170214_01_T12017-01-2904:53:5940
LC08_L1TP_141055_20170214_20170228_01_T12017-02-1404:53:5268
LC08_L1TP_141055_20170302_20170316_01_T12017-03-0204:53:4664
LC08_L1TP_141055_20170318_20170328_01_T12017-03-1804:53:3613
LC08_L1TP_141055_20170403_20170414_01_T12017-04-0304:53:2916
LC08_L1TP_141055_20170419_20170501_01_T12017-04-1904:53:2015
LC08_L1TP_141055_20170505_20170515_01_T12017-05-0504:53:1332
LC08_L1TP_141055_20170606_20170616_01_T12017-06-0604:53:3454
LC08_L1TP_141055_20170622_20170630_01_T12017-06-2204:53:4028
LC08_L1TP_141055_20170708_20170716_01_T12017-07-0804:53:4344
LC08_L1TP_141055_20170724_20170809_01_T12017-07-2404:53:4938
LC08_L1TP_141055_20170809_20170824_01_T12017-08-0904:53:5647
LC08_L1TP_141055_20170825_20170913_01_T12017-08-2504:54:0038
LC08_L1TP_141055_20170910_20170927_01_T12017-09-1004:54:0261
LC08_L1TP_141055_20170926_20171013_01_T12017-09-2604:54:0782
LC08_L1TP_141055_20171012_20171024_01_T12017-10-1204:54:1276
LC08_L1TP_141055_20171028_20171108_01_T12017-10-2804:54:1319
LC08_L1TP_141055_20171113_20171121_01_T12017-11-1304:54:1038
LC08_L1TP_141055_20171215_20171223_01_T12017-12-1504:54:0532
Code A1. The code used to generate median temperature of Landsat 5 TM images for 1996 and 2006.
//Export Landsat 5 SR data
/**
* Function to mask clouds based on the pixel_qa band of Landsat SR data.
* @param {ee.Image} image Input Landsat SR image
* @return {ee.Image} Cloudmasked Landsat image
*/
var cloudMaskL457 = function(image) {
var qa = image.select(‘pixel_qa’);
//If the cloud bit (5) is set and the cloud confidence (7) is high
//or the cloud shadow bit is set (3), then it’s a bad pixel.
var cloud = qa.bitwiseAnd(1 << 5)
.and(qa.bitwiseAnd(1 << 7))
.or(qa.bitwiseAnd(1 << 3));
//Remove edge pixels that don’t occur in all bands
var mask2 = image.mask().reduce(ee.Reducer.min());
return image.updateMask(cloud.not()).updateMask(mask2);
};
//Load the Landsat 5 data collection - example for B6
{
var col = ee.ImageCollection(‘LANDSAT/LT05/C01/T1_SR’)
.map(cloudMaskL457)
.filterDate(‘1996-01-01’,‘1996-12-30’)
.select([‘B6’])
.filterBounds(geometry);
}
//Compute the median of the collection -example for B6
{
var image = col.median().clip(geometry);
print(image, ‘Selected band’);
Map.addLayer(image);
}
//Export the B6, specifying scale and region.
Export.image.toDrive({
image: image,
description: ‘B6’,
scale: 30,
region: geometry,
maxPixels:1e13,
folder: ‘Landsat 5 data collection’,
skipEmptyTiles: true
});
//Export the B6 information as table.
Export.table.toDrive({
collection: col,
description: ‘B6_information’,
fileFormat: ‘CSV’,
folder:’Landsat 5 data collection’
});
//Get the number of collections.
var count = col.size();
print(‘Count: ‘, count);
//Get the date range of images in the collection.
var range = col.reduceColumns(ee.Reducer.minMax(), [“system:time_start”])
print(‘Date range: ‘, ee.Date(range.get(‘min’)), ee.Date(range.get(‘max’)))
Code A2. The code used to generate median temperature of Landsat 8 images for 2017.
//Export Landsat 8 SR data
/**
* Function to mask clouds based on the pixel_qa band of Landsat 8 SR data.
* @param {ee.Image} image input Landsat 8 SR image
* @return {ee.Image} cloudmasked Landsat 8 image
*/
function maskL8sr(image) {
//Bits 3 and 5 are cloud shadow and cloud, respectively.
var cloudShadowBitMask = (1 << 3);
var cloudsBitMask = (1 << 5);
//Get the pixel QA band.
var qa = image.select(‘pixel_qa’);
//Both flags should be set to zero, indicating clear conditions.
var mask = qa.bitwiseAnd(cloudShadowBitMask).eq(0)
.and(qa.bitwiseAnd(cloudsBitMask).eq(0));
return image.updateMask(mask);
}
//Load the Landsat 8 data collection - example for B10
{
var col = ee.ImageCollection(‘LANDSAT/LC08/C01/T1_SR’)
.map(maskL8sr)
.filterDate(‘2018-01-01’,‘2018-12-30’)
.select([‘B10’])
.filterBounds(geometry);
}
//Compute the median of the collection - example for B10
{
var image = col.median().clip(geometry);
print(image, ‘Selected band’);
Map.addLayer(image);
}
//Export the B10, specifying scale and region.
Export.image.toDrive({
image: image,
description: ‘B10’,
scale: 30,
region: geometry,
maxPixels:1e13,
folder: ‘Landsat 8 data collection’,
skipEmptyTiles: true
});
//Export the B10 information as table.
Export.table.toDrive({
collection: col,
description: ‘B10_information’,
fileFormat: ‘CSV’,
folder:’Landsat 8 data collection’
});
//Get the number of collections.
var count = col.size();
print(‘Count: ‘, count);
//Get the date range of images in the collection.
var range = col.reduceColumns(ee.Reducer.minMax(), [“system:time_start”])
print(‘Date range: ‘, ee.Date(range.get(‘min’)), ee.Date(range.get(‘max’)))
Table A2. Error matrix for the classified 1996 land use/cover map, classified.
Table A2. Error matrix for the classified 1996 land use/cover map, classified.
Classified DataReference Data
Built-upForestAgricultural LandOtherWaterTotalUser’s Accuracy (%)
Built-up7286308980.9
Forest2165103318390.2
Agricultural Land5171402216684.3
Other Land2233414281.0
Water0220162080.0
Total811941614222500
Producer’s Accuracy (%)88.985.187.081.072.7
Overall accuracy (%) = 85%.
Table A3. Error matrix for the classified 2006 land use/cover map, classified.
Table A3. Error matrix for the classified 2006 land use/cover map, classified.
Classified DataReference Data
Built-upForestAgricultural LandOtherWaterTotalUser’s Accuracy (%)
Built-up98250010593.3
Forest318280019394.3
Agricultural Land451551116693.4
Other Land00151771.4
Water0111262989.7
Total105190170728500
Producer’s Accuracy (%)93.395.891.271.492.9
Overall accuracy (%) = 93%.
Table A4. Error matrix for the classified 1996 land use/cover map, classified.
Table A4. Error matrix for the classified 1996 land use/cover map, classified.
Classified DataReference Data
Built-upForestAgricultural LandOtherWaterTotalUser’s Accuracy (%)
Built-up110331011794.0
Forest516071017392.5
Agricultural Land321352114394.4
Other Land2323314180.5
Water0110242692.3
Total1201691483726500
Producer’s Accuracy (%)91.794.791.289.292.3
Overall accuracy (%) = 92%.

References

  1. Chen, X.; Zhang, Y. Impacts of urban surface characteristics on spatiotemporal pattern of land surface temperature in Kunming of China. Sustain. Cities Soc. 2017, 32, 87–99. [Google Scholar] [CrossRef] [Green Version]
  2. Thi Van, T.; Bao, H.D.X. Study of the impact of urban development on surface temperature using remote sensing in Ho Chi Minh City, Northern Vietnam. Geogr. Res. 2010, 48, 86–96. [Google Scholar] [CrossRef]
  3. Ranagalage, M.; Estoque, R.C.; Murayama, Y. An urban heat island study of the Colombo Metropolitan Area, Sri Lanka, based on Landsat data (1997–2017). ISPRS Int. J. Geo-Inf. 2017, 6, 189. [Google Scholar] [CrossRef]
  4. Ranagalage, M.; Estoque, R.C.; Zhang, X.; Murayama, Y. Spatial changes of urban heat island formation in the Colombo District, Sri Lanka: Implications for sustainability planning. Sustainability 2018, 10, 1367. [Google Scholar] [CrossRef]
  5. Son, N.-T.; Chen, C.-F.; Chen, C.-R.; Thanh, B.-X.; Vuong, T.-H. Assessment of urbanization and urban heat islands in Ho Chi Minh City, Vietnam using Landsat data. Sustain. Cities Soc. 2017, 30, 150–161. [Google Scholar] [CrossRef]
  6. Hou, H.; Wang, R.; Murayama, Y. Scenario-based modelling for urban sustainability focusing on changes in cropland under rapid urbanization: A case study of Hangzhou from 1990 to 2035. Sci. Total Environ. 2019, 661, 422–431. [Google Scholar] [CrossRef] [PubMed]
  7. Alphan, H. Land-use change and urbanization of Adana, Turkey. Land Degrad. Dev. 2003, 14, 575–586. [Google Scholar] [CrossRef]
  8. Shalaby, A.; Ghar, M.A.; Tateishi, R. Desertification impact assessment in Egypt using low resolution satellite data and GIS. Int. J. Environ. Stud. 2004, 61, 375–383. [Google Scholar] [CrossRef]
  9. Simwanda, M.; Murayama, Y. Spatiotemporal patterns of urban land use change in the rapidly growing city of Lusaka, Zambia: Implications for sustainable urban development. Sustain. Cities Soc. 2018, 39, 262–274. [Google Scholar] [CrossRef]
  10. El Araby, M. Urban growth and environmental degradation. Cities 2002, 19, 389–400. [Google Scholar] [CrossRef]
  11. Debela, T.H.; Beyene, A.; Tesfahun, E.; Getaneh, A.; Gize, A.; Mekonnen, Z. Fecal contamination of soil and water in sub-Saharan Africa cities: The case of Addis Ababa, Ethiopia. Ecohydrol. Hydrobiol. 2018, 18, 225–230. [Google Scholar] [CrossRef]
  12. Knudsen, A.B.; Slooff, R. Vector-borne disease problems in rapid urbanization: New approaches to vector control. Bull. World Health Organ. 1992, 70, 1–6. [Google Scholar] [PubMed]
  13. Ranagalage, M.; Dissanayake, D.; Murayama, Y.; Zhang, X.; Estoque, R.C.; Perera, E.; Morimoto, T. Quantifying surface urban heat island formation in the World Heritage Tropical Mountain City of Sri Lanka. ISPRS Int. J. Geo-Inf. 2018, 7, 341. [Google Scholar] [CrossRef]
  14. Estoque, R.C.; Murayama, Y. Monitoring surface urban heat island formation in a tropical mountain city using Landsat data (1987–2015). ISPRS J. Photogramm. Remote Sens. 2017, 133, 18–29. [Google Scholar] [CrossRef]
  15. Voogt, J.A.; Oke, T.R. Thermal remote sensing of urban climates. Remote Sens. Environ. 2003, 86, 370–384. [Google Scholar] [CrossRef]
  16. 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] [PubMed]
  17. Weng, Q. A remote sensing—GIS evaluation of urban expansion and its impact on surface temperature in the Zhujiang Delta, China. Int. J. Remote Sens. 2014, 22, 1999–2014. [Google Scholar]
  18. Singh, P.; Kikon, N.; Verma, P. Impact of land use change and urbanization on urban heat island in Lucknow city, Central India. A remote sensing based estimate. Sustain. Cities Soc. 2017, 32, 100–114. [Google Scholar] [CrossRef]
  19. Mirzaei, P.A. Recent challenges in modeling of urban heat island. Sustain. Cities Soc. 2015, 19, 200–206. [Google Scholar] [CrossRef] [Green Version]
  20. Rousta, I.; Sarif, M.O.; Gupta, R.D.; Olafsson, H.; Ranagalage, M.; Murayama, Y.; Zhang, H.; Mushore, T.D. Spatiotemporal analysis of land use/land cover and its effects on surface urban heat island using Landsat data: A case study of Metropolitan City Tehran (1988–2018). Sustainability 2018, 10, 4433. [Google Scholar] [CrossRef]
  21. Teferi, E.; Abraha, H. Urban heat island effect of Addis Ababa City: Implications of urban green spaces for climate change adaptation. In Climate Change Adaptation in Africa; Springer: Berlin, Germany, 2017. [Google Scholar]
  22. Tran, H.; Uchihama, D.; Ochi, S.; Yasuoka, Y. Assessment with satellite data of the urban heat island effects in Asian mega cities. Int. J. Appl. Earth Obs. Geoinf. 2006, 8, 34–48. [Google Scholar] [CrossRef]
  23. Van Nguyen, O.; Kawamura, K.; Trong, D.P.; Gong, Z.; Suwandana, E. Temporal change and its spatial variety on land surface temperature and land use changes in the Red River Delta, Vietnam, using MODIS time-series imagery. Environ. Monit. Assess. 2015, 187, 464. [Google Scholar] [CrossRef] [PubMed]
  24. Ayanlade, A. Seasonality in the daytime and night-time intensity of land surface temperature in a tropical city area. Sci. Total Environ. 2016, 557–558, 415–424. [Google Scholar] [CrossRef] [PubMed]
  25. Zhou, D.; Zhang, L.; Hao, L.; Sun, G.; Liu, Y.; Zhu, C. Spatiotemporal trends of urban heat island effect along the urban development intensity gradient in China. Sci. Total Environ. 2016, 544, 617–626. [Google Scholar] [CrossRef] [PubMed]
  26. Marconcini, M.; Metz, A.; Esch, T.; Zeidler, J. Global urban growth monitoring by means of SAR data. In Proceedings of the 2014 IEEE Geoscience and Remote Sensing Symposium, Quebec, QC, Canada, 13–18 July 2014; pp. 1477–1480. [Google Scholar]
  27. Paranunzio, R.; Ceola, S.; Laio, F.; Montanari, A. Evaluating the effects of urbanization evolution on air temperature trends using nightlight satellite data. Atmosphere 2019, 10, 117. [Google Scholar] [CrossRef]
  28. Zhou, Y.; Smith, S.J.; Zhao, K.; Imhoff, M.; Thomson, A.; Bond-Lamberty, B.; Asrar, G.R.; Zhang, X.; He, C.; Elvidge, C.D. A global map of urban extent from nightlights. Environ. Res. Lett. 2015, 10, 2000–2010. [Google Scholar] [CrossRef]
  29. Weng, Q. Thermal infrared remote sensing for urban climate and environmental studies: Methods, applications, and trends. ISPRS J. Photogramm. Remote Sens. 2009, 64, 335–344. [Google Scholar] [CrossRef]
  30. Weng, Q.; Lu, D.; Schubring, J. Estimation of land surface temperature-vegetation abundance relationship for urban heat island studies. Remote Sens. Environ. 2004, 89, 467–483. [Google Scholar] [CrossRef]
  31. Bhaduri, B.; Bright, E.; Coleman, P.; Urban, M.L. LandScan USA: A high-resolution geospatial and temporal modeling approach for population distribution and dynamics. GeoJournal 2007, 69, 103–117. [Google Scholar] [CrossRef]
  32. Zhang, X.; Estoque, R.C.; Murayama, Y. An urban heat island study in Nanchang City, China based on land surface temperature and social-ecological variables. Sustain. Cities Soc. 2017, 32, 557–568. [Google Scholar] [CrossRef]
  33. Dissanayake, D.; Morimoto, T.; Ranagalage, M.; Murayama, Y. Impact of urban surface characteristics and socio-economic variables on the spatial variation of land surface temperature in Lagos City, Nigeria. Sustainability 2019, 11, 25. [Google Scholar] [CrossRef]
  34. Bhargava, A.; Lakmini, S.; Bhargava, S. Urban Heat Island Effect: It’s Relevance in Urban Planning. J. Biodivers. Endanger. Species 2017, 5, 5–187. [Google Scholar]
  35. Xu, S. An approach to analyzing the intensity of the daytime surface urban heat island effect at a local scale. Environ. Monit. Assess. 2009, 151, 289–300. [Google Scholar] [CrossRef] [PubMed]
  36. Li, Y.Y.; Zhang, H.; Kainz, W. Monitoring patterns of urban heat islands of the fast-growing Shanghai metropolis, China: Using time-series of Landsat TM/ETM+ data. Int. J. Appl. Earth Obs. Geoinf. 2012, 19, 127–138. [Google Scholar] [CrossRef]
  37. Senanayake, I.P.; Welivitiya, W.D.D.P.; Nadeeka, P.M. Remote sensing based analysis of urban heat islands with vegetation cover in Colombo city, Sri Lanka using Landsat-7 ETM+ data. Urban Clim. 2013, 5, 19–35. [Google Scholar] [CrossRef]
  38. Sakakibara, Y.; Owa, K. Urban–rural temperature differences in coastal cities: Influence of rural sites. Int. J. Climatol. 2005, 25, 811–820. [Google Scholar] [CrossRef]
  39. Myint, S.W.; Brazel, A.; Okin, G.; Buyantuyev, A. Combined effects of impervious surface and vegetation cover on air temperature variations in a rapidly expanding desert city. GIScience Remote Sens. 2010, 47, 301–320. [Google Scholar] [CrossRef]
  40. Ravanelli, R.; Nascetti, A.; Cirigliano, R.V.; Di Rico, C.; Leuzzi, G.; Monti, P.; Crespi, M. Monitoring the impact of land cover change on surface urban heat island through Google Earth Engine: Proposal of a global methodology, first applications and problems. Remote Sens. 2018, 10, 1488. [Google Scholar] [CrossRef]
  41. Parastatidis, D.; Mitraka, Z.; Chrysoulakis, N.; Abrams, M. Online global land surface temperature estimation from landsat. Remote Sens. 2017, 9, 1208. [Google Scholar] [CrossRef]
  42. Stewart, I.; Oke, T. Classifying urban climate field sites by “local climate zones”: The case of nagano, japan. In Proceedings of the Seventh International Conference on Urban Climate, Yokohama, Japan, 29–30 June 2009. [Google Scholar]
  43. Oke, T.R. Initial Guidance to Obtain Representative Meteorological Observations at Urban Sites; World Meteorological Organization: Geneva, Switzerland, 2004. [Google Scholar]
  44. Ranagalage, M.; Estoque, R.C.; Handayani, H.H.; Zhang, X.; Morimoto, T.; Tadono, T.; Murayama, Y. Relation between urban volume and land surface temperature: A comparative study of planned and traditional cities in Japan. Sustainability 2018, 10, 2366. [Google Scholar] [CrossRef]
  45. Priyankara, P.; Ranagalage, M.; Dissanayake, D.; Morimoto, T.; Murayama, Y. Spatial process of surface urban heat island in rapidly growing Seoul metropolitan area for sustainable urban planning using Landsat Data. Climate 2019, 7, 110. [Google Scholar] [CrossRef]
  46. Estoque, R.C.; Murayama, Y. Quantifying landscape pattern and ecosystem service value changes in four rapidly urbanizing hill stations of Southeast Asia. Landsc. Ecol. 2016, 31, 1481–1507. [Google Scholar] [CrossRef]
  47. Estoque, R.C.; Murayama, Y. City Profile: Baguio. Cities 2013, 30, 240–251. [Google Scholar] [CrossRef]
  48. Crossette, B. The Great Hill Stations of Asia; Basic Books: New York, NY, USA, 1999; Volume 73. [Google Scholar]
  49. Jayasinghe, M.K.D.; Gnanapala, W.K.A.; Sandaruwani, J.A.R. Factors affecting tourists’ perception and satisfaction in Nuwara Eliya, Sri Lanka. Ilorin J. Econ. Policy 2015, 2, 1–15. [Google Scholar]
  50. Weerasinghe, W.W.K. Transformation of the landscape of Nuwara-Eliya. Ph.D. Thesis, University of Moratuwa, Colombo, Sri Lanka, 2003. [Google Scholar]
  51. Dissanayake, D.; Morimoto, T.; Ranagalage, M.; Murayama, Y. Land-use/land-cover changes and their impact on surface urban heat islands: Case study of Kandy City, Sri Lanka. Climate 2019, 7, 99. [Google Scholar] [CrossRef]
  52. Google Erath Engine, Landsat Collection Structure. Available online: https://developers.google.com/earth-engine/landsat (accessed on 25 May 2019).
  53. Google Erath Engine, Image Collection Reductions. Available online: https://developers.google.com/earth-engine/reducers_image_collection (accessed on 25 May 2019).
  54. Sobrino, J.A.; Jiménez-Muñoz, J.C.; Paolini, L. Land surface temperature retrieval from LANDSAT TM 5. Remote Sens. Environ. 2004, 90, 434–440. [Google Scholar] [CrossRef]
  55. Zhang, Y.; Odeh, I.O.A.; 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. Geoinf. 2009, 11, 256–264. [Google Scholar] [CrossRef]
  56. Huang, C.; Davis, L.S.; Townshend, J.R.G. An assessment of support vector machines for land cover classification. Int. J. Remote Sens. 2002, 23, 725–749. [Google Scholar] [CrossRef]
  57. Yang, X. Parameterizing Support Vector Machines for land cover classification. Photogramm. Eng. Remote Sens. 2011, 77, 27–37. [Google Scholar] [CrossRef]
  58. Shi, D.; Yang, X. Support vector machines for land cover mapping from remote sensor imagery. In Monitoring and Modeling of Global Changes: A Geomatics Perspective; Springer: Dordrecht, The Netherlands, 2015; pp. 265–279. [Google Scholar]
  59. Ranagalage, M.; Wang, R.; Gunarathna, M.H.J.P.; Dissanayake, D.; Murayama, Y.; Simwanda, M. Spatial forecasting of the landscape in rapidly urbanizing hill stations of South Asia: A case study of Nuwara Eliya, Sri Lanka (1996–2037). Remote Sens. 2019, 11, 1743. [Google Scholar] [CrossRef]
  60. Gunaalan, K.; Ranagalage, M.; Gunarathna, M.H.J.P.; Kumari, M.K.N.; Vithanage, M.; Srivaratharasan, T.; Saravanan, S.; Warnasuriya, T.W.S. Application of geospatial techniques for groundwater quality and availability assessment: A case study in Jaffna Peninsula, Sri Lanka. ISPRS Int. J. Geo-Inf. 2018, 7, 20. [Google Scholar] [CrossRef]
  61. Simwanda, M.; Ranagalage, M.; Estoque, R.C.; Murayama, Y. Spatial analysis of surface urban heat islands in four rapidly growing African Cities. Remote Sens. 2019, 11, 1645. [Google Scholar] [CrossRef]
  62. Dissanayake, D.; Morimoto, T.; Murayama, Y.; Ranagalage, M. Impact of landscape structure on the variation of land surface temperature in Sub-Saharan Region: A case study of Addis Ababa using Landsat data (1986–2016). Sustainability 2019, 11, 2257. [Google Scholar] [CrossRef]
  63. Department of Forest. Forest Ordinance; Department of Forest: Sri Jayawardenepura Kotte, Sri Lanka, 1908.
  64. Sultana, S.; Satyanarayana, A.N.V. Urban heat island intensity during winter over metropolitan cities of India using remote-sensing techniques: Impact of urbanization. Int. J. Remote Sens. 2018, 39, 6692–6730. [Google Scholar] [CrossRef]
  65. Pepin, N.; Bradley, R.S.; Diaz, H.F.; Baraer, M.; Caceres, E.B.; Forsythe, N.; Fowler, H.; Greenwood, G.; Hashmi, M.Z.; Liu, X.D.; et al. Elevation-dependent warming in mountain regions of the world. Nat. Clim. Chang. 2015, 5, 424–430. [Google Scholar] [Green Version]
  66. Diaz, H.F.; Bradley, R.S.; Ning, L. Climatic changes in mountain regions of the American Cordillera and the tropics: Historical changes and future outlook. Arct. Antarct Alp. Res. 2014, 46, 735–743. [Google Scholar] [CrossRef]
  67. Bradley, R.S.; Keimig, F.T.; Diaz, H.F.; Hardy, D.R. Recent changes in freezing level heights in the Tropics with implications for the deglacierization of high mountain regions. Geophys. Res. Lett. 2009, 36, 2–5. [Google Scholar] [CrossRef]
  68. Dissanayake, D.; Morimoto, T.; Ranagalage, M. Accessing the soil erosion rate based on RUSLE model for sustainable land use management: A case study of the Kotmale watershed, Sri Lanka. Model. Earth Syst. Environ. 2018, 4, 291–306. [Google Scholar] [CrossRef]
  69. Ranagalage, M. Landslide hazards assessment in Nuwara Eliya District in Sri Lanka. In Proceedings of the Japanese Geographical Meeting, Tsukuba, Japan, 28–30 March 2017; p. 100336. [Google Scholar]
  70. Perera, E.N.C.; Jayawardana, D.T.; Jayasinghe, P.; Ranagalage, M. Landslide vulnerability assessment based on entropy method: A case study from Kegalle district, Sri Lanka. Model. Earth Syst. Environ. 2019, 1–15. [Google Scholar] [CrossRef]
  71. Perera, E.N.C.; Jayawardana, D.T.; Ranagalage, M.; Jayasinghe, P. Spatial multi criteria evaluation (SMCE) model for landslide hazard zonation in tropical hilly environment: A case study from Kegalle. Geoinform. Geostat. Overv. 2018, S3. [Google Scholar] [CrossRef]
  72. Thanh Hoan, N.; Liou, Y.-A.; Nguyen, K.-A.; Sharma, R.; Tran, D.-P.; Liou, C.-L.; Cham, D. Assessing the effects of land-use types in surface urban heat islands for developing comfortable living in Hanoi City. Remote Sens. 2018, 10, 1965. [Google Scholar] [CrossRef]
  73. 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]
  74. Galagoda, R.U.; Jayasinghe, G.Y.; Halwatura, R.U.; Rupasinghe, H.T. The impact of urban green infrastructure as a sustainable approach towards tropical micro-climatic changes and human thermal comfort. Urban For. Urban Green. 2018, 34, 1–9. [Google Scholar] [CrossRef]
  75. Manawadu, L.; Ranagalage, M. Urban heat islands and vegetation cover as a controlling factor. In Proceedings of the International Forestry and Environment Symposium 2013 of the Department of Forestry and Environmental Science, University of Sri Jayewardenepura, Thulhiriya, Sri Lanka, 10–11 January 2014; p. 125. [Google Scholar]
  76. Estoque, R.C.; Murayama, Y.; Ranagalage, M.; Hou, H.; Subasinghe, S. Validating ALOS PRISM DSM-derived surface feature height: Implications for urban volume estimation. Tsukuba Geoenviron. Sci. 2017, 13, 13–22. [Google Scholar]
  77. Ranagalage, M.; Murayama, Y. Measurement of urban built-up volume using remote sensing data and geospatial techniques. Tsukuba Geoenviron. Sci. 2018, 14, 19–29. [Google Scholar]
  78. Handayani, H.H.; Murayama, Y.; Ranagalage, M.; Liu, F.; Dissanayake, D. Geospatial analysis of horizontal and vertical urban expansion using multi-spatial resolution data: A case study of Surabaya, Indonesia. Remote Sens. 2018, 10, 1599. [Google Scholar] [CrossRef]
Figure 1. Study area: (a) map from part of South Asia (http://www.maps-world.net); (b) location of Nuwara Eliya; and (c) population density of Nuwara Eliya (100 × 100 m).
Figure 1. Study area: (a) map from part of South Asia (http://www.maps-world.net); (b) location of Nuwara Eliya; and (c) population density of Nuwara Eliya (100 × 100 m).
Sustainability 11 05517 g001
Figure 2. The workflow of the study. Note BL, FL, and AL refer to built-up land, forested land, and agricultural land, respectively.
Figure 2. The workflow of the study. Note BL, FL, and AL refer to built-up land, forested land, and agricultural land, respectively.
Sustainability 11 05517 g002
Figure 3. Graphical illustration of image collection [53].
Figure 3. Graphical illustration of image collection [53].
Sustainability 11 05517 g003
Figure 4. Land use/cover maps and land surface temperature (LST) values in Nuwara Eliya: land use/cover in (a) 1996, (b) 2006, and (c) 2017; and LST in (d) 1996, (e) 2006, and (f) 2017.
Figure 4. Land use/cover maps and land surface temperature (LST) values in Nuwara Eliya: land use/cover in (a) 1996, (b) 2006, and (c) 2017; and LST in (d) 1996, (e) 2006, and (f) 2017.
Sustainability 11 05517 g004
Figure 5. (a) Spatial distribution of mean LST, and the fractions of BL, FL, and AL along the urban–rural gradient; and scatter plots of the mean LST and fractions of (b) BL, (c) FL, and (d) AL. Note: changes in BL, FL, and AL types of land are not independent of each other.
Figure 5. (a) Spatial distribution of mean LST, and the fractions of BL, FL, and AL along the urban–rural gradient; and scatter plots of the mean LST and fractions of (b) BL, (c) FL, and (d) AL. Note: changes in BL, FL, and AL types of land are not independent of each other.
Sustainability 11 05517 g005
Figure 6. (a) Magnitudes of the trends of LST intensity (LSTI)U–R (°C), and fractions of BL, FL, and AL along the URZ; and (b) statistical relationships between the magnitudes of mean LSTs and fractions of BL, FL, and AL in Nuwara Eliya during 1996, 2006, and 2017. Note: changes in BL, FL, and AL types of land are not independent of each other.
Figure 6. (a) Magnitudes of the trends of LST intensity (LSTI)U–R (°C), and fractions of BL, FL, and AL along the URZ; and (b) statistical relationships between the magnitudes of mean LSTs and fractions of BL, FL, and AL in Nuwara Eliya during 1996, 2006, and 2017. Note: changes in BL, FL, and AL types of land are not independent of each other.
Sustainability 11 05517 g006
Figure 7. (a) Spatial distribution of the mean LST, AL/FL fraction ratio; (b) scatter plots of the mean LST and AL/FL fraction ratio; and (c) the Δ(mean LST) (°C) and Δ(AL/FL) fraction ratio in Nuwara Eliya in 1996, 2006, and 2017.
Figure 7. (a) Spatial distribution of the mean LST, AL/FL fraction ratio; (b) scatter plots of the mean LST and AL/FL fraction ratio; and (c) the Δ(mean LST) (°C) and Δ(AL/FL) fraction ratio in Nuwara Eliya in 1996, 2006, and 2017.
Sustainability 11 05517 g007
Figure 8. (a) Spatial distribution of the mean LST and BL/FL fraction ratio; (b) scatter plots of the mean LST and BL/FL fraction ratio; and (c) Δ(mean LST) (°C) and Δ(BL/FL) fraction ratio in Nuwara Eliya in 1996, 2006, and 2017.
Figure 8. (a) Spatial distribution of the mean LST and BL/FL fraction ratio; (b) scatter plots of the mean LST and BL/FL fraction ratio; and (c) Δ(mean LST) (°C) and Δ(BL/FL) fraction ratio in Nuwara Eliya in 1996, 2006, and 2017.
Sustainability 11 05517 g008
Figure 9. (a) Graphical view of the 210 × 210 m polygon grid showing maps of mean LST, built-up density, forest density, and agricultural density in 1996, 2006, and 2017; (b) Scatterplots showing the relationship between the densities of BL, FL, AL, and mean LST in 1996, 2006, and 2017.
Figure 9. (a) Graphical view of the 210 × 210 m polygon grid showing maps of mean LST, built-up density, forest density, and agricultural density in 1996, 2006, and 2017; (b) Scatterplots showing the relationship between the densities of BL, FL, AL, and mean LST in 1996, 2006, and 2017.
Sustainability 11 05517 g009aSustainability 11 05517 g009b
Table 1. Details of the land use/land cover (LULC) changes in Nuwara Eliya (1996, 2006, and 2017).
Table 1. Details of the land use/land cover (LULC) changes in Nuwara Eliya (1996, 2006, and 2017).
Land Use/Cover1996 2006 2017
Area (ha)%Area (ha)%Area (ha)%
Built-up289.91.3785.53.52080.49.3
Forest13,076.758.213,502.160.113,234.358.9
Agricultural Land8503.237.98085.5366583.929.3
Other Land511.82.360481.82.1
Water73.40.375.90.374.60.3
Total22,45510022,45510022,455100
Table 2. LULC changes during 1996–2006, 2006–2017, and 1996–2017.
Table 2. LULC changes during 1996–2006, 2006–2017, and 1996–2017.
Land Use/Cover1996–20062006–20171996–2017
Land Use/Cover Changes (ha)Annual Growth Rate (ha per year)Land Use/Cover Changes (ha)Annual Growth Rate (ha per year)Land Use/Cover Changes (ha)Annual Growth Rate (ha per year)
Built-up495.649.61294.9117.71790.685.3
Forest425.342.5−267.8−24.3157.67.5
Agricultural Land−417.7−41.8−1501.7−136.5−1919.3−91.4
Other Land−505.8−50.6475.743.2−30.1−1.4
Water2.50.3−1.3−0.11.30.1

Share and Cite

MDPI and ACS Style

Ranagalage, M.; Murayama, Y.; Dissanayake, D.; Simwanda, M. The Impacts of Landscape Changes on Annual Mean Land Surface Temperature in the Tropical Mountain City of Sri Lanka: A Case Study of Nuwara Eliya (1996–2017). Sustainability 2019, 11, 5517. https://0-doi-org.brum.beds.ac.uk/10.3390/su11195517

AMA Style

Ranagalage M, Murayama Y, Dissanayake D, Simwanda M. The Impacts of Landscape Changes on Annual Mean Land Surface Temperature in the Tropical Mountain City of Sri Lanka: A Case Study of Nuwara Eliya (1996–2017). Sustainability. 2019; 11(19):5517. https://0-doi-org.brum.beds.ac.uk/10.3390/su11195517

Chicago/Turabian Style

Ranagalage, Manjula, Yuji Murayama, DMSLB Dissanayake, and Matamyo Simwanda. 2019. "The Impacts of Landscape Changes on Annual Mean Land Surface Temperature in the Tropical Mountain City of Sri Lanka: A Case Study of Nuwara Eliya (1996–2017)" Sustainability 11, no. 19: 5517. https://0-doi-org.brum.beds.ac.uk/10.3390/su11195517

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