Spatio-Temporal Features of Urban Heat Island and Its Relationship with Land Use/Cover in Mountainous City: A Case Study in Chongqing

The urban heat island (UHI) becomes more and more serious with the acceleration of urbanization. Many researchers have shown interest in studying the UHI by using remote sensing data, but these studies rarely examine the mountainous cities. Studies on UHI in mountainous cities often used empirical parameters to estimate the land surface temperature (LST), and lacked satellite-ground synchronous experiments to test the accuracy. This paper revised the parameters in the mono-window algorithm used to retrieve the LST according to the characteristics of mountainous cities. This study examined the spatial and temporal patterns of the UHI intensity in Chongqing, a typical mountainous city, and its relationship with land cover from 2007 to 2016 based on the Landsat 5 TM and Landsat 8 TIRS data and the improved method. The accuracy of the LST derivation increased by about 1 ◦C compared to the traditional method. The high LST areas increased and extended from the downtown to suburban area each year, but the rate of change decreased. The UHI is dramatically impacted by the rivers. There is a good relationship between the urban sprawl and the UHI. The LST was reduced by about 1 ◦C within a 300 m distance from large urban fringe green spaces. The urban landscape parks had a strong effect relieving the UHI at a 100 m distance. The LST was reduced by about 0.5 ◦C. This study greatly improves the accuracy of LST derivation, and provides reliable parameters for the UHI researched in mountainous cities.


Introduction
Urbanization is one important attribute of regional social and economic development. It results in land use/cover changes and increased building density, both horizontally and vertically, in urban areas. These changes associated with urbanization impact on urban climate considerably [1,2]. The consistent and significant increase in surface air temperature in urban areas relative to in surrounding rural areas, a phenomenon termed as urban heat island (UHI) effect, has been exacerbated [3]. At the same time, the eco-environment problems caused by UHI become more and more serious with the acceleration of urbanization [4,5]. In recent decades, UHI effects have been the focus of study of sciences, and large UHI effects have been measured and reported for most regions of the world. The UHI of New York City is performed using a mesoscale network of weather stations by Gedzelman et al. [6]. It averages about 4 • C in summer and autumn and 3 • C in winter and spring. According to Weng, urban land development raised surface radiant temperature by 13.01 K between 1989 and 1997 in Zhujiang Delta of South China. The spatial pattern of radiant temperature increase was correlated with the pattern of of the air circulation caused by the terrain barrier, and the high population and building density caused by rapid urbanization [33].

Data and Image Pro-Processing
Satellite remote sensing offers a great opportunity to acquire continuous LST data without direct physical contact with the surface, with sufficient spatial resolution to distinguish between urban and surrounding rural areas [34]. The Landsat 5 TM and Landsat 8 TIRS are the most widely used remote sensing data for UHI studies. Five Landsat 5 TM images (i.e., 20 September 2007;20 July 2008;24 September 2009;11 August 2010 and 30 August 2011) and four Landsat 8 TIRS images (i.e., 19 August 2013; 6 August 2014; 8 July 2015 and 10 July 2016) were used for retrieval of LST in this study. The data acquisition dates had highly clear atmospheric conditions, and the images were acquired from the Institute of Remote Sensing and Digital Earth Chinese Academic of Sciences Data Center, which corrected the radiometric and geometrical distortions of the images before delivery. The images were further rectified to a common Universal Transverse Mercator coordinate system based on 1:10,000 scale topographic maps, and were resampled using the nearest neighbor algorithm with a pixel size of 30 by 30 m for all bands, including the thermal bands (with the resolution of 120 by 120 m and 100 by 100 m for Landsat 5 TM and Landsat 8 TIRS respectively). The resultant RMSE was found to be less than 0.5 pixels. In addition, the meteorological data from 4 standard weather stations, 11 intensive automatic weather stations, and 19 field observing sites, the urban planning data, the land use/cover, and some other auxiliary data were used to assist in the analysis.

Data and Image Pro-Processing
Satellite remote sensing offers a great opportunity to acquire continuous LST data without direct physical contact with the surface, with sufficient spatial resolution to distinguish between urban and surrounding rural areas [34]. The Landsat 5 TM and Landsat 8 TIRS are the most widely used remote sensing data for UHI studies. Five Landsat 5 TM images (i.e., 20 September 2007; 20 July 2008; 24 September 2009; 11 August 2010 and 30 August 2011) and four Landsat 8 TIRS images (i.e., 19 August 2013; 6 August 2014; 8 July 2015 and 10 July 2016) were used for retrieval of LST in this study. The data acquisition dates had highly clear atmospheric conditions, and the images were acquired from the Institute of Remote Sensing and Digital Earth Chinese Academic of Sciences Data Center, which corrected the radiometric and geometrical distortions of the images before delivery. The images were further rectified to a common Universal Transverse Mercator coordinate system based on 1:10,000 scale topographic maps, and were resampled using the nearest neighbor algorithm with a pixel size of 30 by 30 m for all bands, including the thermal bands (with the resolution of 120 by 120 m and 100 by 100 m for Landsat 5 TM and Landsat 8 TIRS respectively). The resultant RMSE was found to be less than 0.5 pixels. In addition, the meteorological data from 4 standard weather stations,  11 intensive automatic weather stations, and 19 field observing sites, the urban planning data, the land use/cover, and some other auxiliary data were used to assist in the analysis.
The sensor TM band 6 and TIRS band 1 on Landsat 5 and Landsat 8 measures the radiances at the top of the atmosphere, and its brightness temperatures can be derived using Plank's law [35]. The approach to the retrieval of temperature was described in the User's Handbook. The following equation was used to convert the digital number (DN) of TIR bands into spectral radiance: where L λ is spectral radiance; DN represents the digital number of band 6 of TM and band 1 of TIRS; L max and L min (given in the header files of the images) are the maximum and minimum spectral radiance of band 6 of TM and band 1 of TIRS respectively. The next step is to convert the spectral radiance to at-satellite brightness temperature under the assumption of uniform emissivity. The conversion formula is: where T is the effective at-satellite brightness temperature in K of band 6 of TM and band 1 of TIRS; L λ is spectral radiance; and K 1 and K 2 are pre-launch calibration constants.

Derivation of LST and Improvement of the Parameters
The methods for the retrieval of LST from satellite TIR (Thermal Infrared) data can be broadly classified into three categories: single-channel methods, multi-channel methods (split-window algorithm, SWA), and multi-angle methods [36][37][38]. In comparison with other methods, the mono-window algorithm is simpler and only three parameters (the average atmospheric operative temperature, land surface emissivity and atmospheric transmittance) are needed; therefore, the mono-window algorithm was used to calculate the LST for the Landsat TM and TIRS images. The mono-window algorithm's calculation formula is as follows [39]: where T s is the emissivity-corrected land surface temperatures, and a and b are regression coefficients. In most studies, the values of the a and b were given roughly as −67.35535, 0.458608, respectively, within the temperature range of 0~70 • C. The images used in this study are all at noon in summer, and the range of temperature variation is small. In this case, the a and b are equal to -67.9542 and 0.45987, respectively (within the temperature range of 20~50 • C); C = τε and D = (1 − τ)[1 + τ(1 − ε)] are the intermediate variables; τ is the atmospheric transmittance; ε is the land surface emissivity; T is the effective at-satellite brightness temperature of band 6 (TM) and band 1 (TIRS); and T a is the effective mean atmospheric temperature. τ, ε and T a are three essential parameters in the mono-window algorithm. Many researchers used these empirical values in their studies. But, in fact, the mono-window algorithm is sensitive to parameters. So it is obvious that the empirical values are not suitable to the retrieval of LST in Chongqing, a typical mountainous city with complex terrain features. In this study, we improved the method to estimate the three essential parameters, and the retrieval accuracy of LST was also improved.

Determination of the Land Surface Emissivity (ε)
The urban area can be divided into water, natural surface, and built-up. If land over is water then ε w = 0.995. If land cover is natural surface, and NDVI < 0.05, then it is bare land, and ε s = 0.972. If NDVI > 0.7, it is fully covered by vegetation, and ε v = 0.986. When NDVI [0.05, 0.7], the land surface emissivity (ε) was computed as follows [19]: When land cover is built-up, the land surface emissivity (ε) was computed as follows [19]: where P v is vegetation coverage; R v = 0.9332 + 0.0585P v is vegetation temperature ratio; ε v = 0.986 is land surface emissivity with very high vegetation coverage; R m = 0.9886 + 0.1287P v is the building surface temperature ratio; and ε m = 0.972 is the building surface emissivity. The effect of d ε is negligible in many studies. They all assume that the land surface is flat. As mentioned before, the terrain of Chongqing is complex. In this paper, we estimated the value of d ε according to the vegetation coverage (when P v ≤ 0.05, d ε = 0.0038; when P v > 0.05, d ε = 0.0038(1 − P v )).

Determination of Atmospheric Transmittance (τ)
The atmospheric transmittance (τ) can be estimated by the following equation [39]: where ω is atmospheric water vapor content in g·cm −2 . In many studies, ω was often estimated by experience. In our study, we calculate it as follows [40]: where a 0 and a 1 are empirical constants.
where T 0 is air temperature, and RH is relative humidity. Due to a lack of meteorological data, most of the studies use the meteorological data from standard weather stations to estimate e. Unfortunately, urban areas often have few weather stations, making the accuracy of e unsatisfactory. In this study, the meteorological data from 4 standard weather stations, 11 intensive automatic weather stations, and 19 field observing sites are combined to determine the value of e. We therefore generated satisfactory results.

Determination of Effective Mean Atmospheric Temperature (T a )
According to Qin et al. for mid-latitude summers, T a can be approximated as [39]: where T 0 is near-surface air temperature. In this paper, we improved the accuracy of T a by using the integrated meteorological data from standard weather stations, intensive automatic weather stations, and field observing sites.

Retrieval Accuracy Validation of LST Based on Satellite-Ground Synchronous Experiment
The best way to validate retrieval accuracy is to compare the in-situ ground truth measurements of LST with the retrieved ones from the remote sensing images data of a specific region. On 30 August 2011, we carried out a satellite-ground synchronous experiment at a Landsat TM pass over the study area. We recorded the air temperature, LST, wind speed, and air humidity at 46 field sites each 15 min Sustainability 2018, 10, 1943 6 of 15 from 9:30 a.m. to 3 p.m. In actual operation, we used the data from 19 sites with high reliability to validate the retrieval accuracy of LST (Table 1). Table 1. Retrieval accuracy validation of LST based on satellite-ground synchronous experiment.

Method of Improved Parameters
Method of Empirical Parameters The results in Table 1 indicate that the LSTs retrieved after improving the parameters are generally higher than the in-situ ground truth measurements of LSTs. The mean values are 34.30 • C and 33.65 • C, respectively, with a difference of 0.65 • C. The maximum, minimum, and mean absolute difference between T true and T s1 are 4.22 • C, 0.17 • C, and 1.86 • C. We also compared the retrieved LST by comparing the empirical parameters with the ground truth LST. The former is generally lower than the latter. The mean value is lower by 1.94 • C. The maximum, minimum, and mean absolute difference between T true and T s2 are 5.08 • C, 0.05 • C, and 2.73 • C. The retrieval accuracy of LST in different land use/cover types are shown in Figure 2. The results suggest that the method of improved parameters has higher accuracy in wood, shrub, residential area, square, and river beach than in the method of empirical parameters. In farmland and pavement, the method of improved parameters shows a lower LST retrieval accuracy. On the whole, the experimental results show that the LSTs retrieved by using the improved parameters are more accurate compared to the LSTs retrieved based on the previous parameters. This good matching of the LSTs retrieved to the actual ones confirms the applicability of the improved parameters to the UHI studies in mountainous cities.
LST images from nine periods Landsat 5 TM6 and Landsat 8 TIRS1 data of the study area were calculated using the improved parameters mono-window algorithm ( Figure 3, Table 2). In order to eliminate the effects of the environmental factors, we normalized the LSTs to 0~1, and divided them to seven levels (i.e., very high temperature area, high temperature area, sub-high temperature area, medium temperature area, sub-medium temperature area, low temperature area, and very low temperature area) by using the natural breaks grading method ( Figure 4, Table 3). method of empirical parameters. In farmland and pavement, the method of improved parameters shows a lower LST retrieval accuracy. On the whole, the experimental results show that the LSTs retrieved by using the improved parameters are more accurate compared to the LSTs retrieved based on the previous parameters. This good matching of the LSTs retrieved to the actual ones confirms the applicability of the improved parameters to the UHI studies in mountainous cities. LST images from nine periods Landsat 5 TM6 and Landsat 8 TIRS1 data of the study area were calculated using the improved parameters mono-window algorithm ( Figure 3, Table 2). In order to eliminate the effects of the environmental factors, we normalized the LSTs to 0~1, and divided them to seven levels (i.e., very high temperature area, high temperature area, sub-high temperature area, medium temperature area, sub-medium temperature area, low temperature area, and very low temperature area) by using the natural breaks grading method ( Figure 4, Table 3).

Spatio-Temporal Features of UHI
The temporal variation and spatial distribution of annual UHI over the study area are shown in

Spatio-Temporal Features of UHI
The temporal variation and spatial distribution of annual UHI over the study area are shown in Figures 3 and 4, Tables 2 and 3. The results suggest that The LSTs show a gradually increasing trend from 2007 to 2011, and a decreasing trend from 2011 to 2016. The spatial pattern of the UHI changed from a scattered pattern in 2007 to a more contiguous pattern in 2011 and 2016, along with the expansion of the regional urban system. The centers of high temperature were consistent with built-up and bare land areas, and the low-temperature centers are located in the area of two rivers (i.e., Yangtze River and Jialing River) and two mountains area (i.e., Mount Zhongliang and Mount Tongluo). The relatively high temperature areas (i.e., very high LST area and high LST area) shifted from south to north and expand from the central city to the outskirts. On 20 August 2011, the mean LST was 37 • C. In most of the study area, the LST was above 35 • C, and the maximum LST was even up to 51.26 • C. Compared to 2007, the mean LST increased by 0.29 • C, the distribution of the relatively high temperature areas became more even. The minimum LST was 28.01 • C, and the maximum LST was 46.52 • C on 10 July 2016. Compared to 2011, the mean LST decreased by 2.39 • C. The distribution of high-temperature centers was more disperse from 2007 to 2016, and their inner structure and composition obviously changed. These distributions and changes of the UHI effect are closely related to the types and functions of the different underlying surface and their changes.
In the study area, the very low LST areas were mainly the large water bodies, such as ponds, reservoirs and rivers, and changed with the water level. There were marked changes from 2007 to 2016, and the area ratio showed an obvious increasing trend from 2007 to 2016. (Table 3) expansion of the regional urban system. The centers of high temperature were consistent with builtup and bare land areas, and the low-temperature centers are located in the area of two rivers (i.e., Yangtze River and Jialing River) and two mountains area (i.e., Mount Zhongliang and Mount Tongluo). The relatively high temperature areas (i.e., very high LST area and high LST area) shifted from south to north and expand from the central city to the outskirts. On 20 August 2011, the mean LST was 37 °C. In most of the study area, the LST was above 35 °C, and the maximum LST was even up to 51.26 °C. Compared to 2007, the mean LST increased by 0.29 °C, the distribution of the relatively high temperature areas became more even. The minimum LST was 28.01 °C, and the maximum LST was 46.52 °C on 10 July 2016. Compared to 2011, the mean LST decreased by 2.39 °C. The distribution of high-temperature centers was more disperse from 2007 to 2016, and their inner structure and composition obviously changed. These distributions and changes of the UHI effect are closely related to the types and functions of the different underlying surface and their changes. In the study area, the very low LST areas were mainly the large water bodies, such as ponds, reservoirs and rivers, and changed with the water level. There were marked changes from 2007 to 2016, and the area ratio showed an obvious increasing trend from 2007 to 2016. (Table 3). The low LST areas concentrated on the high vegetation coverage areas and some small water bodies. After 2009, the area ratio sharply decreased, while it increased quickly from 2011 to 2013. Then, the increasing trend slowed down. The sub-medium LST areas were located in the sub-high vegetation coverage areas and farmlands. They changed in volatility from 2007 to 2016. The medium LST areas  centers, and the expansion of the relatively high temperature areas (from 0.76 km 2 in 2009 to 57.07 km 2 in 2010). In 2011, the relatively high temperature areas decreased compared to 2010 (from 2.35 km 2 in 2010 to 1.31 km 2 in 2011), but it was still far more than it was in 2007, 2008, and 2009. Although the relatively high LST areas were largest in 2010, apparently there was a general trend of increasing UHI effect from 2007 to 2011. From 2011 to 2015, the relatively high temperature areas continually increased. The area ratio increased from 0.17% to 0.80. In 2016, the area ratio decreased to 0.25%.

Relationship between LST and NDVI
The Normalized Difference Vegetation Index (NDVI) was calculated for each image. The negative correlations were found between NDVI and LST by comparing both images. In order to further study their quantitative relationship, the NDVI and LST images on 30 August 2011 are selected as samples to establish the linear regression. Its regression equation was given as follows.
With the increasing of NDVI, the LST decreased from approximately 45 • C to 33 • C. When NDVI was approximately 0.18, some pixels had a higher LST. After analyzing the NDVI, LST and high-resolution images, we find that most of these pixels were in the new construction areas.

Relationship between UHI and Urban Expansion
In the paper, the LST changes were detected from 2007 to 2016, and the amplitude of changes are classified on Table 4, Figure 5.
In Figure 5, we find that the LST decrease areas are mainly the water bodies and high vegetation coverage areas, the no-change areas are distributed in the old urban districts, and the LST increase areas are located at the north and south of the study area.

Relationship between UHI and Land Use/Cover
In this paper, we also explored the effect of different land use/cover types for UHI in a mountainous city. Five land use/cover types were classified in the study area (i.e., built-up, bare land, vegetation, water and road). Vegetation is the area containing forests, shrubs, grasslands and crops. Water includes reservoirs, ponds and rivers. The area under development is similar to bare land. The maps of land use/cover were retrieved from a Landsat images on 20 September 2007, 30 August 2011, and 10 July 2016. To study the temperature relationship of different land use/cover types, the maximum, minimum, and mean LST of different land use/cover types were derived by averaging all corresponding pixel values (Table 5). It can be found from Table 5 that the maximum and mean LST in bare land are all highest followed by build-up, and the maximum and mean LST in water are lowest due to its high heat capacity.
To further study the relationship between LST and land use/cover types, the areas of each LST level in different land use/cover types in 2007, 2011, and 2016 were calculated (Table 6). Using the information, how the land use/cover types may have contributed to UHI can be estimated. From our analysis, it can be seen that almost all land use/cover types except water were distributed over the medium and sub-medium LST areas. The difference is that the areas of bare land and road located in medium LST areas were larger than those in sub-medium LST areas, but the vegetation was the

Relationship between UHI and Land Use/Cover
In this paper, we also explored the effect of different land use/cover types for UHI in a mountainous city. Five land use/cover types were classified in the study area (i.e., built-up, bare land, vegetation, water and road). Vegetation is the area containing forests, shrubs, grasslands and crops. Water includes reservoirs, ponds and rivers. The area under development is similar to bare land. The maps of land use/cover were retrieved from a Landsat images on 20 September 2007, 30 August 2011, and 10 July 2016. To study the temperature relationship of different land use/cover types, the maximum, minimum, and mean LST of different land use/cover types were derived by averaging all corresponding pixel values (Table 5). It can be found from Table 5 that the maximum and mean LST in bare land are all highest followed by build-up, and the maximum and mean LST in water are lowest due to its high heat capacity. To further study the relationship between LST and land use/cover types, the areas of each LST level in different land use/cover types in 2007, 2011, and 2016 were calculated (Table 6). Using the information, how the land use/cover types may have contributed to UHI can be estimated. From our analysis, it can be seen that almost all land use/cover types except water were distributed over the medium and sub-medium LST areas. The difference is that the areas of bare land and road located in medium LST areas were larger than those in sub-medium LST areas, but the vegetation was the opposite. The vegetation areas mainly showed the characteristics of low and sub-medium LST. The relatively high LST areas concentrated in the built-up (0.9 km 2 , 14.11 km 2 , 16.01 km 2 ) and bare land (2.54 km 2 , 8.95 km 2 , 5.95 km 2 ) in 2007, 2011, and 2016. They nearly matched the total contribution. Most of the areas of built-up were medium and sub-medium LST. However, the low and very low LST areas were highly concentrated in water and vegetation. Both of them contributed more than 90% to the total area of low and very low LST areas in the study area. So, it is easy to conclude that water and high coverage vegetation areas play an important role mitigating the UHI effect.

Mitigating Effects of Different Urban Green Spaces on UHI
Two types of urban green spaces (i.e., urban fringe forests and urban parks) were selected in the study. The urban fringe forests with large areas were retrieved from Landsat TM on 30 August 2011 ( Figure 6), and the urban parks data was provided by Chongqing Park Management and Landscape Planning Bureau. The urban fringe forests include three primary parts: Mount Tongluo located in the east, Mount Zhongliang located in the west, and some hills located in the north. Nine urban parks were selected as samples in the paper (i.e., Bolin park, Dongbu park, Eling park, Huahui park, Zhongyang park, Pingdingshan park, Shaping park, Shanhu park, Shimen park). We created six buffers per 50 m for each urban fringe forest and urban park in this study and estimated the mitigating effects of these urban green spaces on UHI (Tables 7 and 8).
From Tables 7 and 8, it can be seen that there are obvious mitigating effects of the urban green spaces on UHI, but the mitigating effects decreased over distance. The urban fringe forests can achieve better performance than urban parks in mitigating UHI effects. For urban fringe forests, the two mountains (Mount Zhongliang in the west and Mount Tongluo in the east) showed a similar trend because of their similar features in area and spatial distribution. The decrease of LST changed from 1.85 • C and 1.93 • C to 1.47 • C and 1.32 • C, respectively, from the first buffer to the sixth buffer. The LSTs of buffers in the north showed a fast decline over distance, from 2.01 • C to 0.53 • C within 200 m. This may be caused by the small area and different spatial pattern of the north urban fringe forest. In general, the urban fringe forests contributed to the decrease of LST by more than 1 • C within Two types of urban green spaces (i.e., urban fringe forests and urban parks) were selected in the study. The urban fringe forests with large areas were retrieved from Landsat TM on 30 August 2011 (Figure 6), and the urban parks data was provided by Chongqing Park Management and Landscape Planning Bureau. The urban fringe forests include three primary parts: Mount Tongluo located in the east, Mount Zhongliang located in the west, and some hills located in the north. Nine urban parks were selected as samples in the paper (i.e., Bolin park, Dongbu park, Eling park, Huahui park, Zhongyang park, Pingdingshan park, Shaping park, Shanhu park, Shimen park). We created six buffers per 50 m for each urban fringe forest and urban park in this study and estimated the mitigating effects of these urban green spaces on UHI (Tables 7 and 8).  Tables 7 and 8, it can be seen that there are obvious mitigating effects of the urban green spaces on UHI, but the mitigating effects decreased over distance. The urban fringe forests can achieve better performance than urban parks in mitigating UHI effects. For urban fringe forests, the two mountains (Mount Zhongliang in the west and Mount Tongluo in the east) showed a similar trend because of their similar features in area and spatial distribution. The decrease of LST changed from 1.85 °C and 1.93 °C to 1.47 °C and 1.32 °C, respectively, from the first buffer to the sixth buffer. The LSTs of buffers in the north showed a fast decline over distance, from 2.01 °C to 0.53 °C within 200 m. This may be caused by the small area and different spatial pattern of the north urban fringe forest. In general, the urban fringe forests contributed to the decrease of LST by more than 1 °C within 300 m. The mitigating effects of the shorter distance to the forests the mitigating effects on UHI are more obvious.   Urban parks are natural landscapes for citizens to rest and can also, to some extent, play an important role in mitigating UHI effects. Table 8 shows the differences between these parks, specifically the different areas and patterns mitigating UHI effects. In general, the LST in the first and second buffer of the parks, except Pingdingshan park, obviously decreased. The average contribution to the LST decrease in the first buffer (0-50 m) was 0.51 • C. The maximum and minimum decrease of LST were 1.73 • C (Zhongyang park) and 0.06 • C (Shaping park), respectively, within the 0-50 m distance. The parks had only slight effects on the LST when the distance exceeded 150 m. Therefore, we should design and plan urban parks (area, pattern, location and quantity) scientifically and reasonably to maximize their effects on the UHI.

Discussion
In this study, we rectify the parameters used in the mono-window algorithm according to the characteristics of Chongqing, a typical mountainous city. We use a combination of Landsat TM data, automatic weather station data, and satellite-ground synchronous experiment data to estimate the LST, test its accuracy, and examine the relationship between the impact factors and the UHI variations in Chongqing. The accuracy of the LST derivation increased by about 1 • C compared to the traditional method and provides reliable parameters for the UHI researched in mountainous city. Furthermore, all the analyses in this paper were based on the interpretation of remote sensing images, by which we analyzed not only the methods to retrieve the UHI, but also the phenomenon of UHI and its impact factors. The temporal and spatial variations of UHI were conducted through the analysis of multi-temporal remote sensing images. The remote sensing images are ideal for analyzing UHI, but it is difficult to select images with similar environmental conditions. In future studies, several additional topics need to be explored. Firstly, the retrieval of UHI in mountainous cities is more complex than the cities on the plain. Although we have improved the method to estimate the three essential parameters, and the retrieval accuracy of LST was also improved, the retrieval method of LST needs to be improved to accommodate the complex terrain features and reduce the influence of unhomogeneous atmospheric conditions in the mountainous cities. Secondly, although we have analyzed the effects of urban spraw, LUCC, and urban green spaces on the UHI, the impact mechanism of the shape, area, structure, and distribution pattern of different land use/cover types in the urbanized area on UHI needs to be further studied. Thirdly, the effect of human activities and other impact factors on UHI must be investigated.

Conclusions
Chongqing is a typical mountainous city. In this paper, qualitative and quantitative analyses were used to study the spatio-temporal features of UHI and the relationship between LUCC and UHI in the study area. Several conclusions were drawn: (1) The accuracy of the LST derivation was improved by about 1 • C compared to the traditional method; (2) The high LST areas increased and extended from the city center to suburban area each year, but the rate of change decreased. There were many minor changes inside high LST areas; (3) There were many rivers in the city and the UHI is dramatically impacted by the rivers. There is a good relationship between the urban expansion and the UHI; and (4) The urban green spaces reduced the effects of UHI, but their functions decreased with the increase of distance from the green spaces. The large urban fringe green spaces composed a relative low temperature area and relieved the UHI in the large area. The LST was reduced by about 1 • C within a 300 m distance from the large urban fringe green spaces. The urban landscape parks had strong effects relieving the UHI at a 100 m distance. The LST was reduced by about 0.5 • C. The study greatly improves the accuracy of LST derivation and provides reliable parameters for the UHI researched in mountainous city.
Author Contributions: Y.L. and C.L. designed the paper. C.L. collected the data and wrote the paper. Y.L. revised the paper. All authors have read and approved the final manuscript.