Mapping the Influence of Land Use / Land Cover Changes on the Urban Heat Island Effect — A Case Study of Changchun , China

The spatio-temporal patterns of land use/land cover changes (LUCC) can significantly affect the distribution and intensity of the urban heat island (UHI) effect. However, few studies have mapped a clear picture of the influence of LUCC on UHI. In this study, both qualitative and quantitative models are employed to explore the effect of LUCC on UHI. UHI and LUCC maps were retrieved from Landsat data acquired from 1984, 1992, 2000, 2007, and 2014 to show their spatio-temporal patterns. The results showed that: (1) both the patterns of LUCC and UHI have had dramatic changes in the past 30 years. The urban area of Changchun increased more than four times, from 143.15 km2 in 1984 to 577.45 km2 in 2014, and the proportion of UHI regions has increased from 15.27% in 1984 to 29.62% in 2014; (2) the spatiotemporal changes in thermal environment were consistent with the process of urbanization. The average LST of the study area has been continuously increasing as many other land use types have been transformed to urban regions. The mean temperatures were higher in urban regions than rural areas over all of the periods, but the UHI intensity varied based on different measurements; and (3) the thermal environment inside the city varied widely even within a small area. The LST possesses a very strong positive relationship with impervious surface area (ISA), and the relationship has become stronger in recent years. The UHI we employ, specifically in this study, is SUHI (surface urban heat island).


Introduction
UHI (urban heat island), which refers to the phenomenon where urban areas tend to have higher temperatures than those of surrounding rural areas, is one of the major problems the 21st century poses to human beings as a result of the unprecedented growth in urbanization [1][2][3].During the process of urbanization natural landscapes, like cropland and forests, are transformed into impervious surfaces consisting of chemical materials that effectively store short-wave radiation, such as cement, metal, and asphalt [4][5][6].The UHI has profound effects on social, economic, and environmental problems, such as human health and well-being, mortality and risk of violence, higher energy costs, air quality, and urban runoff [7][8][9][10].In terms of statistics, the UHI effects are found in more than 1000 cities of different sizes through all latitudes in both hemispheres, and more cities will suffer from it in the future [11].In past years, many efforts have been made to mitigate the negative impacts caused by UHI, such as new roofing and paving materials, and more green space areas [12][13][14][15].Thus, it has become a major research focus in urban climatology and urban ecology since first reported in 1818 [16].
To have a better understanding of the UHI effects is very important for providing a scientific basis for future climate mitigation and urban planning [17,18].
Traditionally, UHI studies have been conducted by ground-based observations [19].In recent years, with the development of remote sensing technology, thermal infrared (TIR) remote sensing has become popular in urban climate and environment studies.It has become a powerful research tool for UHI studies [20][21][22][23].Land surface temperature (LST), one of the most important variables measured by TIR, is considered the skin temperature of the Earth.LST plays an important role in the field of surface energy, environmental, and climate changes [24,25].Compared to the CLHI (canopy layer heat island) and BLHI (boundary layer heat island), the UHI characterized by LST is defined as SUHI (surface urban heat island) and it relates directly to land surface features [26].The SUHI is dealt with in this study.
There has been quite a significant number of publications to investigate the impacts of LUCC (land use/land cover changes) on the thermal environment [27][28][29][30][31]. Rao is the first one to attempt the use of remote sensing technology in the study of UHI [32].In the early stages remote sensing applications, NOAA AVHRR and MODIS, were the main data sources for retrieving LST [33,34].However, the coarse spatial resolution of these data are found suitable only for small-scale urban temperature mapping and cannot be used to establish more accurate relationships between LST and surface biophysical parameters in urban environments [24].Over 40 years, the Landsat series of satellites with higher spatial resolution provide the longest temporal record of space-based surface observations since 1972 and have become more popular among SUHI studies [35,36].Land use types with different thermal properties and radiation features make different contributions to the UHI.For example, normalized differential vegetation index (NDVI ) is traditionally used as the indicator of vegetation abundance to estimate the LST-vegetation relationship, and the scatter diagram proved that there was a negative correlation between NDVI and LST [37,38].Artificial buildings directly change the ratio of surface sensible heat flux and latent heat flux.Impervious surface area (ISA) and normalized different building index (NDBI), which are derived from images, are often used as the indicators to characterize the UHI effect.Goward found that there was a very close relationship between urban land cover and urban temperature patterns, but did not explore the urban inner structure in detail [39].Li et al. selected Shanghai as a case for quantifying the impact of LUCC changes on the patterns of surface UHI and found that the dramatic changes in LUCC have drastically altered the spatiotemporal patterns of UHI [40].At present, high resolution images, such as IKONOS and Quickbird, are widely used in urban LUCC studies, making it possible to explore the relationship between one single object and LST at finer scales [41].
The UHI studies are comprehensive studies which are influenced by the atmospheric environment, regional climate, hydrology, ecology, and many other factors.The remote sensing provides new perspectives for authors to study UHI and allow us to have a better understanding of UHI.However, most of the studies just focus on downtown urban areas over short periods and neglect the spatiotemporal changes of UHI over long periods considering the urban and suburban regions.Thus, this paper aims to: (1) map a clear picture of the spatio-temporal changes in the patterns of LUCC and the thermal environment of Changchun over the past 30 years; (2) analyze the impacts of LUCC on SUHI through both qualitative and quantitative models; and (3) measure UHI intensity and explore its relationship with ISA.

Study Site
Changchun (125 • 06 -125 • 36 E, 43 • 43 -44 • 04 N), known as "China's Detroit", is the capital city of Jilin province in Northeast China and serves as the province's financial, political, and cultural center with an estimated population of 3.66 million in 2014 (Figure 1).Changchun has a sub-humid continental climate with long extremely cold winters and short hot summers.According to the local meteorological records, the average mean air temperature of Changchun (1984-2014) is 5.5 • C (Figure 2) and there has been an obvious increasing trend in the mean air temperature.The highest monthly average temperature is 23.

Data Sources
Landsat images provide the longest continuous records of earth's surface since 1972 and meet the requirements to capture the changes of Earth's land cover and conditions [36].Available Landsat-

Data Sources
Landsat images provide the longest continuous records of earth's surface since 1972 and meet the requirements to capture the changes of Earth's land cover and conditions [36].Available Landsat-

Data Sources
Landsat images provide the longest continuous records of earth's surface since 1972 and meet the requirements to capture the changes of Earth's land cover and conditions [36].Available Landsat-5 and Landsat-8 images (Table 1), collected from US Geological Survey (USGS) in the years from 1984 to 2014, were used to construct time series data to obtain the information of LUCC and to characterize the thermal environment in this study.All of the images were of good quality with no cloud or less than 10% cloud coverage.Prior to interpretation and LST retrievals, both radiometric correction and co-registration were applied to these images, which were corrected to the UTM projection system based on the images of 2014 by using the nearest-neighbor algorithm through ENVI 5.3 (provided by Exelis Visual Information Solutions).GF-1 images with 2 m spatial resolution, acquired on 22 June 2015 from the China Center for Resources Satellite Data and Application (CRESDA), were georeferenced to Landsat images to assess the accuracy of the classification.

LST Retrieval and Its Normalization
Through the radiative transfer equation, Landsat thermal infrared (TIR) data can establish a direct contact with LST [25].The history of surface temperature (ST) retrieval by remote sensing dates back to 1980s [42].Since then, many retrieval algorithms have been proposed based on different assumptions for different data sources [43][44][45][46].Among these algorithms, the mono window algorithm (MWA), developed for Landsat TM6 by Qin [45], was used to map LST in which only three parameters were required: emissivity, transmittance, and effective mean atmospheric temperature, and its basic form can be written as: . T s is LST; T i is the at-sensor brightness temperature of TIR bands i; ε is the land surface emissivity (LSE) of band i; τ is the atmospheric transmittance of band i which can be obtained on the USGS website based on the time of image acquisition and the central latitude and longitude of image; T a is the effective mean atmospheric temperature which can be calculated by using air temperature.T i can be calculated as follows: where the values of constants K 1 and K 2 can be gained from the metadata of Landsat images (USGS).L λ is the top of atmosphere (TOA) spectral radiance (W/(m 2 •sr•µm)) and can be calculated by the pixel values (DN, digital number) of the images.LSE that can change substantially over short distances is one of the most important inputs in LST retrieval algorithms.The NDVI threshold method, which uses certain NDVI values to distinguish different land use types, is chosen to estimate LSE in this study [47].According to the previous studies [48] and LSE libraries, the LSE for vegetation, buildings, soil, and water are 0.9844, 0.9700, 0.9731, and 0.9950, respectively.The NDVI of vegetation (NDVI V ) and buildings (NDVI B ) can be extracted from the NDVI histograms.For those pixels composed of vegetation and buildings (NDVI B ≤ NDVI ≤ NDVIv), LSE can be estimated as follows [49]: where ε mix is the LSE of mixed pixel; ε b is the LSE of buildings; ε v is the LSE of vegetation; P V is the fractional vegetation cover.P V can be obtained from the NDVI according to Carlson and Ripley [50]: According to the location of Changchun and the time of image acquisition, the following equation can be used for the estimation of the mean atmospheric temperature, as shown below [45]: where T a is mean atmospheric temperature (K), and T o is the near-surface air temperature (K) that can be obtained from meteorological station.
It is not proper to directly analyze the differences of LST of multiple year periods due to the seasonality and the inter-annual variability.As a result, LST maps need to be normalized before quantitative analysis, as follows [51]: where NLST is the normalized LST value of pixel i; T i is the initial LST of pixel i; T max and T min are the maximum and minimum LSTs of the whole area, respectively.Then the NLST maps were classified into seven zones: lowest, lower, low, medium, high, higher, and highest using ArcGIS 10.2 (provided by Environmental Systems Research Institute) (Table 2) [18].Areas where NLST belong to high, higher, or highest zones were defined as the UHI regions (UHIR) in this study.NLST mean is the mean value of NLST of the study area, and S is the standard deviation of NLST.

Land Use Classification
Based on the environment of the study area, prior knowledge, and previous studies [52], eight land use/land cover types, as defined in Table 3, were identified from false-color Landsat images in 2014 by visual interpretation through ArcGIS 10.2.Compared to the computer-based interpretation of remote sensing images, visual interpretation usually takes more time and effort, but yields a more accurate result.The land use type of each patch was identified according to the color, size, shape, texture, and its spatial relationships with neighbors.During this process, when it comes to an unknown type, Google Earth with high spatial resolution images and other auxiliary data including soil and vegetation type, DEM (digital elevation model), and climate data were combined to help identify the land use type.Three hundred samples of reference data acquired from GF-1 were used for accuracy assessment.Overall, the total accuracy of the land use map achieves 91.6%, meeting the recommended value [53].For other years, an interpreter with rich experience drew the land-use changes by comparing the Landsat images.According to the theoretical model "vegetation-impervious surface-soil" (V-I-S) and the characteristics of the urban environment, a linear spectral mixture analysis (LSMA) with an improved method to extract and synthesize the "most representative" endmembers was selected to estimate the urban impervious surface area [54,55].LSMA assumes that there is no intersection between endmembers, and it can be described as a linear combination of these endmembers: where R b is the reflectance for band b; N is the number of endmembers; f i is the fraction of endmember; R (i,b) is the reflectance of endmember i in band b; and e b is the residual.Additionally, another two restrictive conditions, where: ∑ N i=1 f i = 1 and o ≤ f i ≤ 1, were implemented to guarantee the real meaning of impervious surface and avoid unreasonable values.One hundred fifty samples with 5 × 5 pixels selected from the estimated impervious surface images were used to conduct accuracy assessment.At the same time, the impervious surface areas of the identical samples on GF-1 images were interpreted visually.The accuracy assessment indicates that more than 80% of the sample's absolute errors are less than 0.1 and the correlation coefficient of estimated and interpreted values is 0.94 which is significant at the 0.05 level of confidence.

The Intensity of UHI
According to the definition of UHI, in this study, the intensity of SUHI was calculated from: where NLST urban is the mean NLST of the urban area, and NLST rural is the mean NLST out of urban area.Additionally, in order to make a comparison with SUHI, ISUHI (Inner SUHI) was defined as follows: where NLST inner_urban is the mean NLST of the first traffic ring line in this study.Compared to the NLST urban , the area of NLST inner_urban was identical during the study periods.

Spatiotemporal Distribution of LST and Statistics
Based on the LST retrieval algorithm mentioned above, the 30-year Landsat LST maps were generated and are shown in Figure 3 to not only measure the magnitude of LST, but also to make it possible to quantify LST spatially explicit over the whole study area.In order to display the LST map clearly, the density slice function in ArcGIS 10.2 was used to distinguish the LST zones by different colors.From the perspective of the spatial characteristics, it showed that the temperature of urban areas was consistently higher than that of rural areas for all the years.In urban areas, the apparent hot spots with the highest LST-the China FAW Group in the southwest, the commercial district in the central area, and the industrial district in the east-could be identified easily.Although the LST in urban areas tended to have higher values, there were still relatively low LST zones, especially in 2014, indicating that both the water and vegetation areas were among the lowest or lower LST zones.By comparison, in rural areas, the forest in the southeast region of the study area and the lake areas had the lowest LST values.Similarly, there were also some high or higher LST regions in the rural areas where there were rural settlements and not much else.
The spatial dynamics of LST showed dramatic changes during 1984-2014.The most remarkable phenomenon was that the area of UHIR had increased tremendously.This was mainly because dramatic changes had taken place in the LUCC patterns which had great influence on the thermal environment.A large amount of cropland, including paddy land and dry land, have been translated to artificial buildings consisting of asphalt, cement, metal, and other chemical materials during this period.In the early years, the highest zones were found inside the urban area, and the higher or high zones had a scattered distribution that looked like an isolated island.As time went by, the isolated hot island had been merged into the central area gradually, creating a larger heat island.
Different from the early years, there were several hot spots in 2014 and the distribution of higher and highest zones were not as concentrated as before.It may be due to the fact that the city had been divided into several different function zones.The average NLST values of the whole study area were 0.51, 0.55, 0.63, 0.62, and 0.68 in 1984, 1992, 2000, 2007, and 2014, respectively.Table 4 showed the area proportion of different NLST zones and UHIR of the whole study area in different years.The UHIR in 2014 accounted for 29.62% of the whole area, which was 15.27% in 1984, indicating that more area would suffer from the effects of UHI.Both the area of higher zone and highest zone have increased, generally, from 1984 to 2014, and after 2007, there was a more obvious increase rate.The From the perspective of the spatial characteristics, it showed that the temperature of urban areas was consistently higher than that of rural areas for all the years.In urban areas, the apparent hot spots with the highest LST-the China FAW Group in the southwest, the commercial district in the central area, and the industrial district in the east-could be identified easily.Although the LST in urban areas tended to have higher values, there were still relatively low LST zones, especially in 2014, indicating that both the water and vegetation areas were among the lowest or lower LST zones.By comparison, in rural areas, the forest in the southeast region of the study area and the lake areas had the lowest LST values.Similarly, there were also some high or higher LST regions in the rural areas where there were rural settlements and not much else.
The spatial dynamics of LST showed dramatic changes during 1984-2014.The most remarkable phenomenon was that the area of UHIR had increased tremendously.This was mainly because dramatic changes had taken place in the LUCC patterns which had great influence on the thermal environment.A large amount of cropland, including paddy land and dry land, have been translated to artificial buildings consisting of asphalt, cement, metal, and other chemical materials during this period.In the early years, the highest zones were found inside the urban area, and the higher or high zones had a scattered distribution that looked like an isolated island.As time went by, the isolated hot island had been merged into the central area gradually, creating a larger heat island.
Different from the early years, there were several hot spots in 2014 and the distribution of higher and highest zones were not as concentrated as before.It may be due to the fact that the city had been divided into several different function zones.The average NLST values of the whole study area were 0.51, 0.55, 0.63, 0.62, and 0.68 in 1984, 1992, 2000, 2007, and 2014, respectively.Table 4 showed the area proportion of different NLST zones and UHIR of the whole study area in different years.The UHIR in 2014 accounted for 29.62% of the whole area, which was 15.27% in 1984, indicating that more area would suffer from the effects of UHI.Both the area of higher zone and highest zone have increased, generally, from 1984 to 2014, and after 2007, there was a more obvious increase rate.The area of the low zone had decreased from 45.89% in 1984 to 26.38% in 2014, indicating that a large amount of natural surfaces had been transformed into artificial surfaces.Near-ground air temperature collected from local meteorological stations were used to validate the precision of the retrieved LST from Landsat images.The average LST of the 9 × 9 pixels near the meteorological station was compared to the air temperature, and the average satellite-based LST was about 3.8 • C higher than the average air temperature in 2014, which was reasonable in terms of previous studies [56].It is well known that performing LST accuracy assessments is a complex and difficult task due to the scale effect and the difficulty to in situ LST data [25].After the process of normalization, in summary, the retrieved LST can reflect the evolution process of the thermal environment.Near-ground air temperature collected from local meteorological stations were used to validate the precision of the retrieved LST from Landsat images.The average LST of the 9 × 9 pixels near the meteorological station was compared to the air temperature, and the average satellite-based LST was about 3.8 °C higher than the average air temperature in 2014, which was reasonable in terms of previous studies [56].It is well known that performing LST accuracy assessments is a complex and difficult task due to the scale effect and the difficulty to obtain in situ LST data [25].After the process of normalization, in summary, the retrieved LST can reflect the evolution process of the thermal environment.Thus, we focused on the changes of the urban area, rural settlements and dry land (Figure 5 and Table 5).The urban area had increased from 143.51 km 2 in 1984 to 577.45 km 2 in 2014 with a mean annual expansion area of 14.46 km 2 .Especially, from 2007 to 2014, the urban area increased by 185.87 km 2 with the largest expansion rate.During the whole study period, dry land was the major contributor (76.82%) followed by rural settlements (17.81%) to the new emerging urban area.For the latter transformation, small towns expanded gradually by merging surrounding smaller settlements and becoming small urban areas.These small urban areas were ultimately merged into the main urban areas.The area of dry land has decreased from 1130.12 km 2 in 1984 to 737.70 km 2 in 2014, a decrease of 34.72%.

LUCC Dynamics
Sustainability 2017, 9, 312 9 of 17 Thus, we focused on the changes of the urban area, rural settlements and dry land (Figure 5 and Table 5).The urban area had increased from 143.51 km 2 in 1984 to 577.45 km 2 in 2014 with a mean annual expansion area of 14.46 km 2 .Especially, from 2007 to 2014, the urban area increased by 185.87 km 2 with the largest expansion rate.During the whole study period, dry land was the major contributor (76.82%) followed by rural settlements (17.81%) to the new emerging urban area.For the latter transformation, small towns expanded gradually by merging surrounding smaller settlements and becoming small urban areas.These small urban areas were ultimately merged into the main urban areas.The area of dry land has decreased from 1130.12 km 2 in 1984 to 737.70 km 2 in 2014, a decrease of 34.72%.

The Difference of NLST among Land Use Classes
The distribution of LST is the reflection of the patterns of land use types with different thermal characteristics.Zonal Statistic in ArcGIS 10.2 was used to quantify the differences of NLST among different land use classes, and the results were shown in Figure 6.It showed that the NLST of different land use types had obvious differences.Based on the mean annual values of NLST, urban areas had the highest NLST (0.76) and was followed by other the built-up areas (0.69), unused land (0.67), rural settlements (0.64), dry land (0.53), paddy land (0.49), woodland (0.47), and water areas (0.42).Urban areas are always linked to buildings and denser population.The thermal characteristics of chemical materials abates surface evaporation and enhances the sensible heat flux, together with anthropogenic heat sources, making the urban area become the hottest regions.Other built-up areas and unused land usually have little vegetation cover.As a result, both of them have relatively higher NLST.Previous studies showed that the vegetation could reduce UHI through two major processes: shading and evapotranspiration [57,58].Thus, land use types with high vegetation cover tend to have lower NLST.The heat capacity of water is higher than other types, so the temperature of water bodies is lower in the daytime.

The Difference of NLST among Land Use Classes
The distribution of LST is the reflection of the patterns of land use types with different thermal characteristics.Zonal Statistic in ArcGIS 10.2 was used to quantify the differences of NLST among different land use classes, and the results were shown in Figure 6.It showed that the NLST of different land use types had obvious differences.Based on the mean annual values of NLST, urban areas had the highest NLST (0.76) and was followed by other the built-up areas (0.69), unused land (0.67), rural settlements (0.64), dry land (0.53), paddy land (0.49), woodland (0.47), and water areas (0.42).Urban areas are always linked to buildings and denser population.The thermal characteristics of chemical materials abates surface evaporation and enhances the sensible heat flux, together with anthropogenic heat sources, making the urban area become the hottest regions.Other built-up areas and unused land usually have little vegetation cover.As a result, both of them have relatively higher NLST.Previous studies showed that the vegetation could reduce UHI through two major processes: shading and evapotranspiration [57,58].Thus, land use types with high vegetation cover tend to have lower NLST.The heat capacity of water is higher than other types, so the temperature of water bodies is lower in the daytime.In order to investigate the impacts of LUCC on LST, and to show the explicit spatial influence, a profile from west to east was generated through the profiles tool in Envi5.3.Additionally, four samples with 100 × 100 pixels were selected on the profile.Two of them (a and c) were on the edge of the urban boundary and had dramatic changes in land use.The left side of sample (a) and the right side of sample (c) were the left and right boundaries of the city, respectively.Sample (b), which had little changes in land use, was located inside the city while (d) was in the rural area and also faced few land use changes.
As shown clearly in Figure 7, the NLST on the profile had considerable fluctuations from west to east.The phenomenon that the temperature of urban areas tended to be higher than that of rural areas was apparent in both 1984 and 2014.The NLST of urban areas also had big differences.This was because the environment of the urban area was very complicated and the surface features can change greatly in a very short distance.For example, sample (b) is the largest park, named "South Lake Park", in Changchun whose NLST is lower compared to surrounding areas.For samples (a) and (c), the changes of land use patterns lead to the significant changes of the NLST distribution.Based on the statistics in this study, on average, the NLST can increase by 0.21 when the natural landscape was transformed to impervious surfaces.For sample (b), there were almost no changes in NLST as the surface features were almost identical.For sample (d), compared to 1984, several small hot spots emerged as some industrial or mining districts appeared.There is no denying that urbanization can change the thermal environment considerably.Therefore, when mitigating the negative effects of UHI, it counts a significant amount to optimize the land use planning by reasonable allocations of different land use types.In order to investigate the impacts of LUCC on LST, and to show the explicit spatial influence, a profile from west to east was generated through the profiles tool in Envi5.3.Additionally, four samples with 100 × 100 pixels were selected on the profile.Two of them (a and c) were on the edge of the urban boundary and had dramatic changes in land use.The left side of sample (a) and the right side of sample (c) were the left and right boundaries of the city, respectively.Sample (b), which had little changes in land use, was located inside the city while (d) was in the rural area and also faced few land use changes.
As shown clearly in Figure 7, the NLST on the profile had considerable fluctuations from west to east.The phenomenon that the temperature of urban areas tended to be higher than that of rural areas was apparent in both 1984 and 2014.The NLST of urban areas also had big differences.This was because the environment of the urban area was very complicated and the surface features can change greatly in a very short distance.For example, sample (b) is the largest park, named "South Lake Park", in Changchun whose NLST is lower compared to surrounding areas.For samples (a) and (c), the changes of land use patterns lead to the significant changes of the NLST distribution.Based on the statistics in this study, on average, the NLST can increase by 0.21 when the natural landscape was transformed to impervious surfaces.For sample (b), there were almost no changes in NLST as the surface features were almost identical.For sample (d), compared to 1984, several small hot spots emerged as some industrial or mining districts appeared.There is no denying that urbanization can change the thermal environment considerably.Therefore, when mitigating the negative effects of UHI, it counts a significant amount to optimize the land use planning by reasonable allocations of different land use types.In order to investigate the impacts of LUCC on LST, and to show the explicit spatial influence, a profile from west to east was generated through the profiles tool in Envi5.3.Additionally, four samples with 100 × 100 pixels were selected on the profile.Two of them (a and c) were on the edge of the urban boundary and had dramatic changes in land use.The left side of sample (a) and the right side of sample (c) were the left and right boundaries of the city, respectively.Sample (b), which had little changes in land use, was located inside the city while (d) was in the rural area and also faced few land use changes.
As shown clearly in Figure 7, the NLST on the profile had considerable fluctuations from west to east.The phenomenon that the temperature of urban areas tended to be higher than that of rural areas was apparent in both 1984 and 2014.The NLST of urban areas also had big differences.This was because the environment of the urban area was very complicated and the surface features can change greatly in a very short distance.For example, sample (b) is the largest park, named "South Lake Park", in Changchun whose NLST is lower compared to surrounding areas.For samples (a) and (c), the changes of land use patterns lead to the significant changes of the NLST distribution.Based on the statistics in this study, on average, the NLST can increase by 0.21 when the natural landscape was transformed to impervious surfaces.For sample (b), there were almost no changes in NLST as the surface features were almost identical.For sample (d), compared to 1984, several small hot spots emerged as some industrial or mining districts appeared.There is no denying that urbanization can change the thermal environment considerably.Therefore, when mitigating the negative effects of UHI, it counts a significant amount to optimize the land use planning by reasonable allocations of different land use types.

The Changes of UHI Intensity
As mentioned before, two kinds of UHI intensity were used in this study, SUHI and ISUHI.The largest difference between them lies in how to represent the NSLT of urban areas.In the traditional ways, the mean NLST of the whole urban area is used to represent the NSLT of urban areas.While, the mean NLST of the first traffic ring inside the city is used in ISHUI.Figure 8 showed the boundary of the urban area in different years and the same first traffic ring boundary.

The Changes of UHI Intensity
As mentioned before, two kinds of UHI intensity were used in this study, SUHI and ISUHI.The largest difference between them lies in how to represent the NSLT of urban areas.In the traditional ways, the mean NLST of the whole urban area is used to represent the NSLT of urban areas.While, the mean NLST of the first traffic ring inside the city is used in ISHUI.Figure 8 showed the boundary of the urban area in different years and the same first traffic ring boundary.

The Changes of UHI Intensity
As mentioned before, two kinds of UHI intensity were used in this study, SUHI and ISUHI.The largest difference between them lies in how to represent the NSLT of urban areas.In the traditional ways, the mean NLST of the whole urban area is used to represent the NSLT of urban areas.While, the mean NLST of the first traffic ring inside the city is used in ISHUI.Figure 8 showed the boundary of the urban area in different years and the same first traffic ring boundary.Urban impervious surface area (ISA) is one of the most important indicators to characterize the degree of urbanization and environmental quality [56].Using the methods mentioned above, Figure 10 showed the statistics of the ISA of the whole urban area and first traffic ring line.It can be seen that the percentage of ISA in the first ring has continuously increased from 1984 to 2007 and approached saturation since 2007.Although the area of urban region has been increasing since 1984, on the contrary, the percentage of ISA of the whole urban area has decreased since 2000.This can be explained by two main reasons.First, the new emerging urban area added to the old ones may have a low ISA value; Second, more water and green space were placed in the old regions to decrease the ISA of the whole urban area.
Comparing Urban impervious surface area (ISA) is one of the most important indicators to characterize the degree of urbanization and environmental quality [56].Using the methods mentioned above, Figure 10 showed the statistics of the ISA of the whole urban area and first traffic ring line.It can be seen that the percentage of ISA in the first ring has continuously increased from 1984 to 2007 and approached saturation since 2007.Although the area of urban region has been increasing since 1984, on the contrary, the percentage of ISA of the whole urban area has decreased since 2000.This can be explained by two main reasons.First, the new emerging urban area added to the old ones may have a low ISA value.Second, more water and green space were placed in the old regions to decrease the ISA of the whole urban area.
Comparing  In order to quantify the relationships between NLST and ISA, 300 samples with 3 × 3 pixels were selected in the urban area in the year 1984 and 2014.A zonal statistics tool was used to calculate the mean values of NLST and ISA for each sample.Then regression models were further performed with the help of SPSS 19.0 (provided by International Business Machines Corporation).The two dimensional scatter plots between NLST and ISA were shown in Figure 11.We can see that there are more high values in 2014 than 1984, indicating that there are more artificial surfaces in these samples.In order to quantify the relationships between NLST and ISA, 300 samples with 3 × 3 pixels were selected in the urban area in the year 1984 and 2014.A zonal statistics tool was used to calculate the mean values of NLST and ISA for each sample.Then regression models were further performed with the help of SPSS 19.0 (provided by International Business Machines Corporation).The two dimensional scatter plots between NLST and ISA were shown in Figure 11.We can see that there are more high values in 2014 than 1984, indicating that there are more artificial surfaces in these samples.
Urban impervious surface area (ISA) is one of the most important indicators to characterize the degree of urbanization and environmental quality [56].Using the methods mentioned above, Figure 10 showed the statistics of the ISA of the whole urban area and first traffic ring line.It can be seen that the percentage of ISA in the first ring has continuously increased from 1984 to 2007 and approached saturation since 2007.Although the area of urban region has been increasing since 1984, on the contrary, the percentage of ISA of the whole urban area has decreased since 2000.This can be explained by two main reasons.First, the new emerging urban area added to the old ones may have a low ISA value.Second, more water and green space were placed in the old regions to decrease the ISA of the whole urban area.
Comparing  In order to quantify the relationships between NLST and ISA, 300 samples with 3 × 3 pixels were selected in the urban area in the year 1984 and 2014.A zonal statistics tool was used to calculate the mean values of NLST and ISA for each sample.Then regression models were further performed with the help of SPSS 19.0 (provided by International Business Machines Corporation).The two dimensional scatter plots between NLST and ISA were shown in Figure 11.We can see that there are more high values in 2014 than 1984, indicating that there are more artificial surfaces in these samples.Additionally, positive linear relationships between them were found in both years, but in 2014 the NLST had a stronger relationship with ISA than that of 1984.Why is the relationship stronger and the scope larger in 2014?In other words, why is the UHI intensity stronger in the sample with the same values of ISA in 2014 compared to that of 1984?One possible reason is that although the value of ISA is same, there are more buildings in the sample (the sample is more dense) or the buildings are taller, which would generate more heat.

Discussion
As an emerging science, remote sensing, with the advantages of wide coverage, good temporal synchronization and various data sources, has been widely used in UHI studies as a supplement to traditional approaches, which have been conducted by ground-based observation taken from standard height metrological stations [24,26,59].The definition of UHI is clear.However, the most fundamental issue are the ways to measure UHI from the perspective of remote sensing.In this study, five-stages of Landsat data were used to provide a clear picture of the spatio-temporal changes in the patterns of LUCC and the thermal environment.Mapping the UHI intensity is very important because it can provide more explicit spatial information about the thermal environment.Additionally, two methods were used in this study to quantify the changes of UHI intensity, but yielded different results.It is difficult to say which method is better than the other and difficult to conclude the change trend of UHI intensity of Changchun in the past 30 years; is it stronger or weaker?Both of the methods are reasonable.However, it will help to generalize a better understanding of the definition and measurement of UHI.
Based on the information provided by the maps of LUCC and NLST, we found several important characteristics of the process of urbanization and the impacts of LUCC on NLST.Firstly, in general, the urban regions have increased continuously since 1984, but after 2000 the expansion rate was much larger.This was because the central government of China proposed a policy to help promote the development of the northeast part of China in 2003 and 2006.With the stimulation policy, several small local towns around Changchun began to develop quickly and, ultimately, merged into the main urban area.Thus, the land use and development policy do make differences on urbanization; Secondly, the differences of mean NLST among land use types can be as high as 0.34.Therefore, the land use changes could lead to considerable changes on the thermal environment.In general, the most common transformation is from cropland to urban areas in our case.According to the mean values of NLST, when this transformation took place, the land would change from low or lower LST zones directly to higher or highest LST zones.During the analysis process, the profile method with samples made it very clear to show the impacts of LUCC on LST.This approach may provide some implications for other similar research.
This study also has some implications for urban and land use planning.Our results showed that inside the urban area, the LST could change dramatically in a short distance as a result of the complicated urban environment.We found that there was a very strong positive relationship between LST and ISA, which was similar to other studies [56,60].A higher percentage of ISA leads to higher LST, and the effect of ISA on LST has become more and more apparent in recent years.That is because urban expansion not only lies in the size, but also in the height and density, of those spaces.In our case, a larger urban area does not mean a stronger UHI intensity.It is our opinion that the negative effect of UHI could be mitigated by altering the ISA of the urban regions.Specific advice and measures include adding more cooling land use, such as water bodies and woodlands in urban areas.We found, in our study, that the new land merged into urban regions had the effect of lowering UHI intensity because of a relatively low value of the ISA.
There are also some limitations in this study.It is better to use year by year Landsat data to explore the impacts of LUCC on thermal environments in a finer time scale which may provide more effective information for both the readers and policy-makers.However, it is difficult to gain quality images because of the poor atmospheric conditions.Second, more attention should be paid to the validation of satellite-derived LST, since no results from remote sensing data without validation can be used with confidence [25].Finally, the characteristics of the thermal environment inside the urban area should be given more attention in future research.

Conclusions
It is important to understand the mechanism of the impacts of LUCC on the thermal environment for better urban planning and UHI mitigation.In this paper, taking Changchun, a cold temperature zone city, as a case study, both qualitative and quantitative analyses were used to investigate a long term LUCC and UHI dynamics based on the information extracted form Landsat series data from 1984 to 2014.Particularly, we used NLST to characterize the patterns of the thermal environment and analyzed the differences of NLST among land use types, and then quantified the effects of urbanization on UHI intensity.
Three main conclusions can be made from the results of this study.First, both of the patterns of LUCC and the thermal environment have had dramatic changes in the past 30 years.The urban area of Changchun increased more than four times from 143.15 km 2 in 1984 to 577.45 km 2 in 2014.The proportion of UHI regions in the whole study area has increased from 15.27% in 1984 to 29.62% in 2014; Second, the changes of LUCC do have a great influence on the thermal environment.In other words, the spatiotemporal changes in thermal environments were consistent with the process of urbanization.Based on the mean NLST values of different land use types, urban areas have the highest NLST.Thus, the average NLST of the study area has been continuously increasing as many other land use types have been transformed to urban regions during the process.Additionally, the mean temperatures were higher in urban regions than rural ones over all the periods, but the UHI intensity varied based on different measurements.Third, the thermal environment inside the city changes considerably in a very short distance.The NLST possessed a very strong positive relationship with ISA, and the positive relationship has become stronger in recent years.In our case, a larger urban area does not mean a stronger UHI intensity.
The UHI effect will be a continuing challenge as the world continues to urbanize, as cities grow in size and density.With the development of remote sensing TIR sensors and LST retrieval algorithms, we can gather more accurate and comprehensive LST data to characterize UHI.On this basis, remote sensing with the help of other discipline's methods will allow us to conduct UHI studies at wide temporal and spatial scales with different themes.The development of UHI studies is a great challenge that requires the concerted effort of multiple disciplines.However, such an effort would contribute substantially to advance our analytical capacity for understanding UHI and, thus, ultimately provide a basis for forging sustainable urban strategies.
1 • C in July and the lowest is −15.6 • C in January.The average annual precipitation of Changchun is 561.6 mm, most of which occurs from June to August when the temperature and humidity are highly persistent.The topography of the study area is flat with minimal elevation changes from 250 m to 350 m.Like many other megacities, Changchun has been undergoing an accelerated rate of urbanization and population explosion.In the past 30 years, dramatic changes have occurred in land use and landscape patterns, which have altered the thermal environment significantly.Sustainability 2017, 9, 312 3 of 17 and there has been an obvious increasing trend in the mean air temperature.The highest monthly average temperature is 23.1 °C in July and the lowest is −15.6 °C in January.The average annual precipitation of Changchun is 561.6 mm, most of which occurs from June to August when the temperature and humidity are highly persistent.The topography of the study area is flat with minimal elevation changes from 250 m to 350 m.Like many other megacities, Changchun has been undergoing an accelerated rate of urbanization and population explosion.In the past 30 years, dramatic changes have occurred in land use and landscape patterns, which have altered the thermal environment significantly.

Figure 1 .
Figure 1.Location of the study site.

Figure 2 .
Figure 2. The average mean air temperature of Changchun from 1984-2014.

Figure 1 .
Figure 1.Location of the study site.

Figure 1 .
Figure 1.Location of the study site.

Figure 2 .
Figure 2. The average mean air temperature of Changchun from 1984-2014.

Figure 2 .
Figure 2. The average mean air temperature of Changchun from 1984-2014.

Figure 3 .
Figure 3.The NLST maps of the study area in different years.

Figure 3 .
Figure 3.The NLST maps of the study area in different years.

Figure 4
Figure 4 demonstrated the dynamics of LUCC in the study area in the past 30 years.It is clear that Changchun has experienced a rapid process of urbanization and the LUCC patterns have changed considerably.Among the eight land use classes, the area of woodland, water, and unused land changed little; the area of paddy land had fluctuations depending on the availability of water; the area of dry land and rural settlements had a continuously decreasing trend, while the urban area and other built-up land had been expanding continuously.

Figure 4
Figure 4 demonstrated the dynamics of LUCC in the study area in the past 30 years.It is clear that Changchun has experienced a rapid process of urbanization and the LUCC patterns have changed considerably.Among the eight land use classes, the area of woodland, water, and unused land changed little; the area of paddy land had fluctuations depending on the availability of water; the area of dry land and rural settlements had a continuously decreasing trend, while the urban area and other built-up land had been expanding continuously.

Figure 5 .
Figure 5.The LUCC statistics of Changchun from 1984 to 2014.

Figure 5 .
Figure 5.The LUCC statistics of Changchun from 1984 to 2014.

Figure 6 .
Figure 6.The statistics of NLST for different land use types.

Figure 7 .
Figure 7.The profile of NLST in 1984 and 2014 from west to east with samples (a-d) in 100 × 100 pixels.

Figure 6 .
Figure 6.The statistics of NLST for different land use types.

Figure 7 .
Figure 7.The profile of NLST in 1984 and 2014 from west to east with samples (a-d) in 100 × 100 pixels.

Figure 7 .
Figure 7.The profile of NLST in 1984 and 2014 from west to east with samples (a-d) in 100 × 100 pixels.

Figure 9
Figure 9 shows the changes of the intensity of UHI during the study periods.From Figure 9, a very interesting phenomenon was that the SUHI and ISUHI showed different change trends.The values of ISUHI in 1984, 1992, 2000, 2007, and 2014 were 0.341, 0.348, 0.346, 0.359, and 0.374, respectively.The ISUHI generally has an upward trend.In comparison, the values of SUHI in 1984, 1992, 2000, 2007, and 2014 were 0.225, 0.249, 0.252, 0.243, and 0.239.SUHI increased before 2000, but decreased from 2000 to 2014.According to the information collected from the National Economy and Society statistical bulletin of Changchun from 2000 to 2014, the area of green space in the built-up area of Changchun has increased from 6058 ha (38.1%) in 2000 to 18,244 ha (41.5%) in 2014.The increase of the area proportion of the green space may explain the decrease of the SUHI of Changchun since 2000, to some extent.As shown in Figure 8, the area of urban regions has increased a great deal since 1984.Therefore, it does not mean that the larger area of urban regions can generate stronger UHI intensity.This interesting question motivated us to conduct further research.

Figure 9 .
Figure 9.The changes of SUHI and ISUHI in Changchun.

Figure 9
Figure 9 shows the changes of the intensity of UHI during the study periods.From Figure 9, a very interesting phenomenon was that the SUHI and ISUHI showed different change trends.The values of ISUHI in 1984, 1992, 2000, 2007, and 2014 were 0.341, 0.348, 0.346, 0.359, and 0.374, respectively.The ISUHI generally has an upward trend.In comparison, the values of SUHI in 1984, 1992, 2000, 2007, and 2014 were 0.225, 0.249, 0.252, 0.243, and 0.239.SUHI increased before 2000, but decreased from 2000 to 2014.According to the information collected from the National Economy and Society statistical bulletin of Changchun from 2000 to 2014, the area of green space in the built-up area of Changchun has increased from 6058 ha (38.1%) in 2000 to 18,244 ha (41.5%) in 2014.The increase of the area proportion of the green space may explain the decrease of the SUHI of Changchun since 2000, to some extent.As shown in Figure8, the area of urban regions has increased a great deal since 1984.Therefore, it does not mean that the larger area of urban regions can generate stronger UHI intensity.This interesting question motivated us to conduct further research.

Figure 9
Figure 9 shows the changes of the intensity of UHI during the study periods.From Figure 9, a very interesting phenomenon was that the SUHI and ISUHI showed different change trends.The values of ISUHI in 1984, 1992, 2000, 2007, and 2014 were 0.341, 0.348, 0.346, 0.359, and 0.374, respectively.The ISUHI generally has an upward trend.In comparison, the values of SUHI in 1984, 1992, 2000, 2007, and 2014 were 0.225, 0.249, 0.252, 0.243, and 0.239.SUHI increased before 2000, but decreased from 2000 to 2014.According to the information collected from the National Economy and Society statistical bulletin of Changchun from 2000 to 2014, the area of green space in the built-up area of Changchun has increased from 6058 ha (38.1%) in 2000 to 18,244 ha (41.5%) in 2014.The increase of the area proportion of the green space may explain the decrease of the SUHI of Changchun since 2000, to some extent.As shown in Figure 8, the area of urban regions has increased a great deal since 1984.Therefore, it does not mean that the larger area of urban regions can generate stronger UHI intensity.This interesting question motivated us to conduct further research.

Figure 9 .
Figure 9.The changes of SUHI and ISUHI in Changchun.Figure 9.The changes of SUHI and ISUHI in Changchun.

Figure 9 .
Figure 9.The changes of SUHI and ISUHI in Changchun.Figure 9.The changes of SUHI and ISUHI in Changchun.
Figures 9 and 10, the intensity of UHI corresponds well with the ISA values.The correlation coefficient of UHI intensity and ISA is 0.97 which indicates that the UHI intensity possess a strong correlation with the ISA.The intensity of UHI is influenced by many factors, such as the definition of UHI, the components and structure of urban area and the proportion of cooling land use types like water bodies and woodland.
Figures 9 and 10, the intensity of UHI corresponds well with the ISA values.The correlation coefficient of UHI intensity and ISA is 0.97 which indicates that the UHI intensity possess a strong correlation with the ISA.The intensity of UHI is influenced by many factors, such as the definition of UHI, the components and structure of urban area and the proportion of cooling land use types like water bodies and woodland.

Figure 10 .
Figure 10.The percentage ISA of the whole urban area and first traffic ring.

Figure 11 .
Figure 11.The regression results between NLST and ISA in (a) 1984 and (b) 2014.

Figure 10 .
Figure 10.The percentage ISA of the whole urban area and first traffic ring.
Figures 9 and 10, the intensity of UHI corresponds well with the ISA values.The correlation coefficient of UHI intensity and ISA is 0.97 which indicates that the UHI intensity possess a strong correlation with the ISA.The intensity of UHI is influenced by many factors, such as the definition of UHI, the components and structure of urban area and the proportion of cooling land use types like water bodies and woodland.

Figure 10 .
Figure 10.The percentage ISA of the whole urban area and first traffic ring.

Figure 11 .
Figure 11.The regression results between NLST and ISA in (a) 1984 and (b) 2014.Figure 11.The regression results between NLST and ISA in (a) 1984 and (b) 2014.

Figure 11 .
Figure 11.The regression results between NLST and ISA in (a) 1984 and (b) 2014.Figure 11.The regression results between NLST and ISA in (a) 1984 and (b) 2014.

Table 1 .
The satellite data used in the study.

Table 2 .
The classification of NLST zones based on the NLST range.

Table 3 .
The land use classification system and the details can be seen in[52].

Table 4 .
The area proportion (%) of different NLST zones and UHIR in different years.
Sustainability 2017, 9, 312 8 of 17 area of the low zone had decreased from 45.89% in 1984 to 26.38% in 2014, indicating that a large amount of natural surfaces had been transformed into artificial surfaces.

Table 4 .
The area proportion (%) of different NLST zones and UHIR in different years.