1. Introduction
Knowledge of land surface temperature aids understanding of the temporal and spatial variations in global land surfaces [1] . This is because LST and emissivity are vital parameters in earth surface characteristics for energy budget estimation and land cover assessment [2] . LST data are typically retrieved from raw Landsat datasets by converting the digital number values of the thermal bands into absolute radiance values [3] [4] . These radiance values are then used to find satellite brightness temperatures, calculated under the assumption of unity emissivity and using pre-launch calibration constants [3] [4] [5] . Emissivity-based methods have been applied to urban areas [6] and vegetation cover-based methods to agricultural areas [7] . Current understanding of spatial thermal patterns and their association with surface characteristics is based on numerous studies of urban surface temperatures over the past two decades [8] . Based on [9] , urban areas are generally 3.5˚C - 4.5˚C warmer than surrounding rural areas. [10] projected that more than 5.87 million km2 of land globally will probably be converted into urban areas by 2030, and 20% of this area has a probability of less than 75% of further urban expansion. [8] used detailed land-cover maps with reasonably similar surfaces to help in measuring temperature gradients across a region. This is because surface temperature is uniquely related to surface properties. Estimating thermal conditions of land surface using satellite images requires the relationship between surface temperature, surrounding topography, and land use and land cover to be determined [4] [11] .
Cities experience direct environmental impacts, such as uneven effects from temperature increases, variations in glacial melt, and sea-level rise. Possible future in-creases in temperature are expected to worsen the urban heat island (UHI) effects [12] , where cities become warmer than their surrounding area [13] [14] . This warming is of particular concern for urban areas due to the UHI effects [15] . The UHI effect results from the retention of solar energy in building materials during the day and lower rates of radiant cooling at night. Urban areas commonly experience lower wind speeds, less convective heat loss and evapotranspiration, and more energy from surface warming [16] . As a product of urbanization and industrialization, UHI is believed to be one of the main dilemmas of the century [17] . Rapid urbanization imposes dramatic land surface alteration and replacement of natural landscapes, which can significantly change the surface radiation, thermal properties, and humidity [18] [19] . Transformation of natural land cover and the need for urban materials lead to climatic changes [11] .
According to [8] , the UHI effect occurs when land cover is transformed, mostly via the replacement of natural vegetation and agricultural land by concrete, asphalt, and building walls related to urban land uses. Extreme heat in cities is constantly generated and released into the atmosphere by building infrastructure and transportation [20] . Urban are-as have consequently become vital drivers of local and regional climatic and environmental changes, with harmful consequences for social and ecological processes [21] . Commercial and residential buildings, roads, and parking lots covered by impervious surfaces worsen the UHI effect, whereas vegetation, such as trees, grass, and shrubs, can alleviate it [22] . In the past few decades, the land use land cover (LULC) of the UAE has significantly changed. Where the vegetation has in-creased and the built up areas have expanded [23] . [24] found that large parks in Dubai City had significantly lower temperatures than the rest of the city during summer. Increasing the proportion of green spaces, such as public parks, gardens, trees, and green roofs, is an effective strategy for urban areas to adapt both climate change and the UHI [25] [22] [26] . The preliminary hypothesis of this study is that urban centers would portray higher temperature levels than the surrounding desert surface. The objective is to detect and study the change in LST with respect to the change in Land cover in the past 30 years in Al-Ain City.
2. Study Area and Satellite Data
Al-Ain City is located in the eastern part of the Emirate of Abu Dhabi, United Arab Emirates (UAE). The area is characterized by an arid climate of extremely hot summers and warm winters, with an average annual temperature of 27˚C and a mean annual rainfall of 96.4 mm [27] . The city is considered one of the largest oases in the Arabian Peninsula. The study area covered the airport area in the northwest, Buraimi border (Oman) in the northeast, central district in the east, Sanaeya industrial area in the southeast, Maqam in the west, and Jabal Hafeet and cement factory in the south. The central residential area includes some palm plantation farms and parks (see Figure 1). Landsat images captured during
Figure 1. Landsat false color images of the study area in Al-Ain City, United Arab Emirates taken from original Landsat images.
the summer (June 18, 1988; June 12, 1994; August 2, 2001; June 16, 2013; and June 11, 2017) (http://earthexplorer.usgs.gov/) were used in this study (Table 1). Images were selected during the summer to portray the peak LST of the study area.
3. Methodology
The image processing workflow shown in Figure 2 was followed:
1) Image normalization was applied to eliminate the effect of different sensors and dates. The image normalization techniques developed by [28] were used to reduce artificial LST changes introduced by sensor technology, illumination, and viewing geometry changes on the thermal bands of the images.
2) Using the equation (mb1 + b). m = (Br Dr)/(Bs Ds), b = (Dr × Bs) (Ds × Br)/(Bs Ds), b1 is the thermal band of the sub image, Br is the mean temperature of the bright pixels of the reference image, Dr is the mean temperature of the dark pixels of the reference image, Bs is the mean temperature of the bright pixels of the sub image, Ds is the mean temperature of the dark pixels of the sub image. Corrected values are shown in Table 2.
Atmospheric correction and LST retrieval was done. LST retrieval from satellite data is considered an important but challenging topic in thermal studies of urban environments [29] . A single-channel algorithm developed by [30] was employed, which includes an updated fit using MODTRAN4 radiative transfer code [31] . This algorithm has proven its effectiveness for Landsat 8 thermal
Table 1. Details of Landsat images used in the study.
Figure 2. Image processing workflow followed in this study.
Table 2. Al-Ain City normalization equation variables.
infrared sensors (band 10) [32] . In this algorithm, the emissivity effects were corrected using an approach developed by [33] . LST values were then converted from Kelvin to degrees Celsius. Table 3 shows the corrected mean temperatures.
3) Image warping was performed to locate and match the tie points in a base and warp image to ensure that they were in the same location that was changed after normalization. The thermal differences are shown in Figure 3.
4) Raster color slicing was applied by setting the temperature change value based on multiple trails and chose the most likely proper values for each difference image. Blue color indicates the cooler areas and red indicate the warmer areas (Figure 4).
5) The green space areas in km2 of the study area was found using NDVI by applying the equation (NIR − R)/(NIR + R), taking into consideration that for images taking from Landsat 4 - 7, NIR is band 4 and R is band 3 while in Landsat 8 the NIR is band 5 and R is band 5 (see Table 4).
4. Results and Discussion
Land cover areas were manually identified in the thermal images based on our
Table 3. Al-Ain City corrected mean temperatures in study period.
Figure 3. Thermal difference images of Al-Ain City throughout the study period.
Table 4. Green spaces in the study area.
Figure 4. Difference in land surface temperature in Al-Ain City throughout the study period.
knowledge of the area and available land use and land cover maps from AL Ain city Municipality. Our results showed that, the mean LST in Al-Ain city has increased by 1.31˚C from 53.37˚C in 1988 to 54.68˚C in 2017. The difference in LST from 1988 to 2017 was 3.5 degrees the largest difference is from 2001 to 1994 by 4.5 degrees (Figure 4). Table 5 shows the changes in mean temperature in four areas in Al-Ain City, namely the airport area, Falaj Hazzaa, Sanaeya industrial area, and Maqam rugby club.
The temperature in the airport area decreased by 3.49˚C from 1988 to 2017. In Falaj Hazzaa, the temperature decreased by 5.1˚C between 2017 and 1988. In contrast, the temperature in Sanaeya has increased by 4˚C since 1988. The temperature in the rugby club area decreased by 4.39˚C during the study period. It was noticed that the green space in the study area over the study period decreased in 2001, was 3.3 km2 in 1994 and dropped to 1.1 km2 in 2001 (Table 4).
[34] , reported that the hottest urban surfaces are paved areas during both day and night because of their high thermal inertia, which allows them to absorb and store more sunlight. In addition to the presence of factories, metal workshops, and other related services in the industrial area (Sanaeya), the high density and high traffic as well as the absence of plantations and green areas, Sanaeya experienced increased temperature compared to other areas in Al-Ain. The rugby club in Maqam was established in 2007, and the plot was previously bare land. Changes in temperature could have been due to the vegetation cover in the rugby course and surrounding urbanized areas. Figure 5 shows the temperature trends in Al-Ain City.
In a study by [35] on Abu Dhabi city, they found that in desert area the UHI phenomenon is commonly inverted because of the high amount of vegetation in urban areas compared to suburban areas where bare ground and sand are found. They also reported that, during the night, the temperature of the different surfaces drops and the suburbs become cooler compared to downtown areas; consequently, the typical UHI effect can be detected only at night [35] . According to [36] , commonly the UHI effect makes the city center warmer than the surrounding areas. Yet this phenomenon is inverted in desert environments. This inversion is due to the vegetation found in urban areas with comparison to suburban areas, where bare lands and hot sand dunes are spreading. [36] ’ results showed that the UHI has been affected by the change in vegetation fraction and a negative relationship between vegetation fraction and temperature was reported.
Table 5. Mean land surface temperature change (˚C) found in Al-Ain City.
Figure 5. Temperature trends in degree Celsius over the study period in Al-Ain City.
The findings of [37] study on Dubai City, showed that LST has greatly reduced because of the increase in urban vegetation cover. In this study the difference in LST from 1988 to 2017 was 3.5 degrees the largest difference is from 2001 to 1994 by 4.5 degrees. This period experienced a drop in the green area which can explain the big variance in LST. The seasonal variation can be noticed clearly in downtown vegetated areas. Vegetation is stressed due to high temperatures in summer months, accordingly NDVI values drop [38] . In this study, the green space in the study area over the study period decreased in 2001 (Table 4). This decrease can be due to seasonal changes. The 2001 image was taken in August, whereas the other images were taken in June. It was not possible to use an image of June 2001 because of the defects in the images from that month. The findings of a study by [29] on Dubai showed that LST decreases when the urban areas are planted and the amount of green space influences the urban areas. In exception of the decrease in 2001, the green area was increasing with time which can explain the decrease in LST.
5. Conclusion
Urbanization is generally thought to lead to the replacement of natural landscapes, resulting in changes in surface radiation, thermal properties, and humidity. Urban areas can be 3.5˚C - 4.5˚C warmer than neighboring rural areas due to the UHI effect. Contrary to previous studies, urbanization of desert surfaces in this study did not show a pronounced UHI effect. Generally, a decrease in LST was observed rather than an increase in the studied urban areas compared to surrounding desert surfaces. LST decreased by 3.49 in the airport area and 5.1˚C in Falaj Hazzaa from 1988 to 2017. Green areas, such as the rugby club, also experienced a temperature decrease of 4.39˚C over the study period. In contrast, Sanaeya industrial area has experienced an increased LST of about 4˚C since 1988. The overall decrease was likely primarily due to a significant increase in urban vegetation, such as trees, public parks, and lawns in residential areas. The Sanaeya industrial area showed a classical UHI due to the absence of such urban vegetation
Acknowledgements
The authors would like to acknowledge the UAE University Graduate School for funding this project.