Optimal Satellite Sensor Selection Utilized to Monitor the Impact of Urban Sprawl on the Thermal Environment in Doha City, Qatar

In the present study, the most appropriate Landsat and ASTER Image spatial and spectral resolutions were selected to extract thermal infrared data in Doha City, Qatar from 1990-2015. Results indicated Urban Heat Island (UHI) calculated from Landsat sensor data was more consistent with ground truthed temperatures. Landsat TM thermal infrared data with low spatial resolution (60-120 m) was more appropriate for large-scale thermal studies, but was not suitable to assess complex urban thermal environments and establish Land Surface Temperatures (LSTs) for individual buildings. However, ASTER sensor, with a Thermal Infrared (TIR) subsystem of 90 m spatial resolution, showed more accurate determination of thermal patterns and LSTs. Landsat data revealed relatively high temperatures, which exhibited increased congruence with ground truthed measures. In Doha City, LST showed an increased trend from 1990 to 2015; low in 1990, moderate in 2000, to very high in 2015. Furthermore, our results confirmed accelerated urban sprawl during 1990-2015. The evaluated areas increased from 80 km2 in 1990 to 179 km2 in 2015. A strong positive relationship was observed between UHI and urban area intensity. Urban planning absence significantly affected the rural environment. Urbanized areas exhibited high thermal conductivity and radiation heat budgets, characterized by increased atmospheric and surface temperatures compared to surrounding rural areas. Therefore, rapid urbanization rates must undergo proper planning and management regimes to insure positive impacts on the local climate and environment. Journal of Earth Science & Climatic Change J o u r n a l o f E art hS cience &Cmatic C h a n g e


Introduction
Doha is a boom city on the Arabian littoral and the singular preeminent urban agglomeration on the low Qatari peninsula [1]. It is Qatar's largest city, with over 60% of the nation's population residing there. The city has grown from a port of less than thirty thousand inhabitants in the mid-twentieth century to today's cosmopolitan extravaganza having over a million people [1]. Doha is characterized by the striking modernity of the urban built space, and perhaps also the pace of social change for the inhabitants therein. Those inhabitants include Qataris, who number less than three hundred thousand, but also more than a million other people who come and go as guests of a sort, mostly to work. Like many other cities on the coast of the Arabian Peninsula, Doha has long been entwined in the maritime trade routes of the Indian Ocean world [1]. In recent years, it has become one of the major and most important producers and exporters of Liquefied Natural Gas (LNG) in the world. The wealth produced by Qatar's oil and gas exports has generated an unprecedented construction boom in the city and its surrounding peripheries. This has resulted in significant growth at all levels, from urban development and infrastructure provision to the building of office towers and sophisticated architecturally innovative cultural and educational facilities. Doha grew around the native pearling industry. In addition, hydrocarbons were discovered in 1939, but their exploitation and a sequence of further discoveries would escalate only in the waning decades of the twentieth century [1]. It is the economic center of the country located on the coast of the Arabian Gulf. The wealth generated by the resources was substantial, and much of it has been expended urban development. Recently, increasingly large geographical portions of the urban landscape are encompassed in singular planned ventures. Zones, partitions, walls, enclaves, and compounds are familiar aspects of everyday life in Doha. These are the "zoning technologies" that demarcate the increasingly large units of urban planned space in Doha. During the late 1990s and particularly during the first decade of the new millennium the structure of Qatar's governance was made new due to the establishment of stateowned and semi-privatized institutions intended to implement the liberalization strategy. This resulted in rapid urban growth in Doha and its metropolitan area, which grew from around 500,000 inhabitants in the late 1990s to 1.5 million in 2011 [1]. The everchanging and increasing population not only poses challenges to planners and decision makers, but architects also have to envisage and design sustainable and pleasing environments for the needs of a future. The first revenues made from oil were mainly invested in the development of infrastructure such as the construction of a regional road connecting the Dukhan field with Doha and the first airport, which was built in the east of the city [2,3].
Remote sensing and geographic information system techniques are considered powerful and effective tools and are widely applied to detect the spatio-temporal dynamics of land use/cover changes [4]. LST is derived from satellite thermal infrared (TIR) imagery [5,6]. This study aimed to select the most appropriate remote sensing sensors to monitor changes in the thermal environment of Doha City due to urban sprawl from 1990-2015.

Methods and Techniques
Four sets of remote sensing data were obtained, from the USGS website, to conduct this research, including Landsat 5 Thematic Mapper (TM) during 1990, Enhanced Landsat Thematic Mapper Images Plus (ETM+) during 2000, ASTER images during 2013, and Landsat 8 during 2015. All sensors represented the winter season, i.e. January and February. Spatial, spectral, and temporal resolutions of all sensors are summarized in Table 1 [7][8][9][10].
Remote sensing image processing and Geographic Information Systems (GIS) techniques were utilized to classify land cover and extract the UHI in Doha City from 1990-2015.

Image classification
Image classification was applied to automatically categorize all pixels in an image into land cover classes or themes. Unsupervised classification was performed separately on the remotely-sensed images using a histogram peak cluster technique to identify dense or frequently occurring pixelated areas [11][12][13][14][15]. In the unsupervised approach, spectrally separable classes were determined to form a supervised classification scheme. ENVI 4.8 was employed to process satellite image classifications [16]. Authors applied supervised maximum likelihood classification (MLC) to detect land cover classes. Once a sufficient number of spectral subclasses were obtained for all information classes, MLC was performed with the full set of refined spectral classes [15]. MLC is a statistically complex approach to classify images for change detection; however several studies reported the methodology was most effective for supervised classification [17,18]. Here, supervised classification was verified for each date using field investigation information collected from more than 34 ground checkpoints. ENVI 4.8 and ERDAS Imagine 9.3 were used to retrieve UHI for the study area. Landsat calibration equations were obtained and corrected from Landsat 7 and Landsat 8 handbooks [7,8]. ASTER calibration equations and required parameters were obtained from the ASTER user handbook [9,10].

Calculation of urban heat island
In order to retrieve UHI of the study area, the following processing was performed: Convert DN values into spectral radiance: Landsat DN values were converted into spectral radiance (TOA) in units (W m -2 µm -1 sr -1 ) using the Landsat calibration equation following author [19]. All specific calibration coefficients and other related parameters (such as sun elevation, gain, and offset) were obtained from a level 1-product header or ancillary data record. For Landsat 8, the following equation was used: where L λ represents spectral radiance in units (W m -2 µm -1 sr -1 ), M L indicates the band specific multiplication rescaling factor obtained from metadata, A L represents the band specific rescaling additive factor obtained from metadata, and Q cal the quantified and calibrated standard product pixel values (DN) [8].
ASTER data radiance can be obtained from DN values as follows: Where the unit conversion coefficient of each band was obtained from the ASTER user handbook [9,10].
Convert the spectral radiance into effective radiance at sensor brightness temperature: For Landsat 5,7, and 8 and ASTER data, the spectral radiance of thermal infrared bands were converted into effective radiance at sensor brightness temperature (the temperature values of black body) using Planck's function equation as follows: where T i represents sensor brightness temperature in kelvin; C 1 and C 2 are constants; λ is emitted radiance wavelength for peak response and average limiting wavelengths (λ=11.5 µm for Landsat band 6 data, λ=10.9 µm and 12 µm for Landsat 8 bands 10 and 11, respectively; λ=10.6 µm and 11.3 µm for ASTER bands 13 and 14, respectively); and L i is spectral radiance, [20,21].
Landsat data were converted using the following equation in which Planck's function is similar to the two free parameters K 1 and K 2 [19,22]: Where T i indicates effective at-satellite brightness temperature in kelvin; K 2 represents the calibration constant 2 in kelvin; K 1 is the calibration constant 1 in W/ (m 2 .sr.µm); and L λ represents spectral radiance at the sensor's aperture ( Table 2).

Land surface temperature calculation
Satellite thermal infrared sensors measure Top of the Atmosphere (TOA) radiances, from which brightness temperature (known as blackbody temperatures) can be derived based on Plank's law. TOA radiances are the result of mixing three parts of energy. The first is the emitted radiance from earth's surface, the second is the upwelling radiance from the atmosphere, and the third is the down welling radiance from the sky [23,24]. The difference among TOA Band-specific thermal conversion constant from the metadata associated with the images(scene dependent) and land surface brightness temperature is subject to the influence of atmospheric conditions. Therefore, to obtain an actual land surface brightness temperature, atmospheric effects, including absorption upward emission and downward irradiance reflected from the surface, should be corrected first. This correction was conducted by calculating spectral emissivity (ɛ), [23,24]. LSTs were obtained by recovering satellite temperature T i by applying the correction for emissivity. The following steps were performed:

Calculation of land surface emissivity ɛ
The effects of land surface emissivity on satellite measurements can be generalized into three categories [23]: i) emissivity causes a reduction of surface-emitted radiance; ii) non-black surfaces reflect radiance; and iii) anisotropy of reflectivity and emissivity might reduce or increase total surface radiance. It is known that emissivity as a function of wavelength is controlled by surface water content, chemical composition, structure, and roughness. For vegetated areas, emissivity can vary significantly with plant species, areal densities, and growth rates. In fact, land surface emissivity is closely related to NDVI. Therefore, emissivity can be estimated using the following steps: Computing NDVI: Normalized Difference Vegetation Index (NDVI) was used to transform multispectral data into a single image band representing vegetation distribution. NDVI values between -1 and 1 indicate the amount of green vegetation present in a pixel. The higher the NDVI value, the more green vegetation exists. A standard algorithm was computed using ENVI (4.8) software following the equation in [25]. For Aster data, this equation was applied using band 2; (0.63-0.69 µm) and band 3N; (0.78-0.86 µm). The equation was also applied to Landsat 8 using bands 4 (0.64-0.67) and 5 (0.85-0.88).

LST normalizing and obtaining urban heat island
Although we chose to sample the study area during approximately the same season, it was not appropriate to compare multiple data images from different years. This was due to inter annual variation and different atmospheric conditions within the same time period among years. Therefore, to compare LST from different dates, a normalization method was performed using the following equation [26]: where U i represents the LST of pixel i after normalizing, T si is LST of date i in degree kelvin as obtained in step 3, T min and T max indicate the minimum and maximum LST retrieved from all data sets, extremes excluded.
The output value U i varied from 0 to 1 (0 ≤ U i ≤ 1.0), which was classified as follows: Class 4 and class 5 were merged and defined as the UHI zone.

Results and Discussion
Image classification was applied to determine changes among the available images. Unsupervised/supervised classification was employed to obtain detailed land use classification data. The classes in the unsupervised classification scheme were very similar to those generated in the supervised classification. Maximum likelihood supervised classification was applied to Landsat images acquired from 1990, 2000, and 2015 for land cover class assessment. Four land use categories were identified, including sea-water, urban areas, and sand and vegetation cover (Figure 1). Each class was ground truthed using a Garmin 38 GPS unit; more than 34 ground data sites were visited and confirmed with GPS coordinates. The results were consistent with urban area acceleration from 1990-2015, with an increase from 80 km 2 in 1990 to 179 km 2 in 2015.Remote sensing data were used to monitor discrete land cover types and estimate biophysical characteristics of land surfaces [28]. Potential parameters describing urban construction materials, composition, and urban canopy structure can be derived from remote sensing approaches. Consequently, UHI effects produced from urban surface energy budgets can be recognized. High resolution remotely sensed thermal imagery has the advantage of providing a time-synchronized dense grid of temperature data over an entire city and distinctive temperatures for individual buildings [27]. Remote sensing technology was used to monitor the impacts of surface water and vegetation cover on the thermal environment. Abutaleb et al., assessed UBI over greater Cairo using remotely sensed imagery during 2002-2012 [29]. Results indicated expansion of urban areas increased land surface thermal radiation in highly populated areas. The growing interest in thermal remote sensing methodologies to extract LST and urban energy budgets is due to the importance of this information to study urban climatology, global environmental changes, human-environment interactions, and planning and management practices. A strong positive relationship was observed between UHI intensity and urban areas. Urban planning significantly affected the thermal environment. Urbanized areas exhibited high thermal conductivity and radiation heat budgets, characterized by increased atmospheric and surface temperatures compared to surrounding rural areas [30].
Ahmed and kaiser reported local thermal changes occurred during rapid urbanization [31]. Cultivated land expansion and relocation of industrial facilities from cities to suburban and rural areas achieved further reductions in UHI. In the present study, remote sensing technology was used to monitor urban growth impacts on the thermal environment from 1990-2015 in Doha City. These changes  were analyzed using thermal infrared. Thermal bands were processed to obtain radiant surface temperatures to spatially and temporally investigate UHI effects associated with increased urban areas over the 25-year period ( Figure 2). LST increased from low in 1990 to moderate in 2000 to very high in 2015. UHI and the energy budget extracted from Landsat 5 and 7 thermal infrared data were significantly correlated with urban planning and development activities. Vegetation cover reduced the amount of solar radiation reaching man-made surfaces; in addition, vegetation caused cooling of the surrounding air by releasing moisture through a process known as evapotranspiration [32]. Therefore, proper planning and management regimes must be put in place as rapid urbanization proceeds to mitigate negative impacts on the local climate and environment. All processed data were acquired in the early morning (approximately 7:00) in the winter season; i.e. January and February. The actual mean near-surface air temperature was kindly supplied from the Qatar International airport. Measured data obtained from the field were used to verify calculated LST extracted from remote sensing images. UHI was also extracted from ASTER images for 2013 ( Figure 3). The results were compared with those extracted from the Landsat 8 image. Results showed UHI calculated from the Landsat sensor were more consistent with ground truthed temperatures.

Conclusions
Analysis of four remote sensing datasets, including Landsat 5 Thematic Mapper (TM) from 1990, ETM+ from 2000, ASTER images from 2013, and Landsat 8 from 2015 established the most appropriate remotely sensed approach to monitor thermal environment change in Doha City due to urbanization sprawl during the 25-year period from 1990-2015. All sensors were chosen to represent the same season, which was during the Qatar January and February winter season. Remote sensing image processing and GIS techniques were utilized to classify land cover and extract UHI from data obtained for Doha City from 1990-2015. Unsupervised/supervised classification was employed to obtain detailed land use classification data. Four land use categories were identified, including seawater, urban areas, sand and vegetation covers. Each class was verified in the field using a Garmin 38 GPS unit; more than 34 ground data sites were visited and confirmed with GPS coordinates. The results were consistent with urban area acceleration from 1990-2015, with an increase from 80 km 2 in 1990 to 179 km 2 in 2015. A strong positive correlation was observed between UHI intensity and urban areas. Urban planning significantly affected the thermal environment. Urbanized areas exhibited high thermal conductivity and radiation heat budgets, characterized by increased atmospheric and surface temperatures compared to surrounding rural areas. In Doha City, LST increased from low in 1990 to moderate in 2000 to very high in 2015. UHI and the energy budget extracted from Landsat 5 and 7 thermal infrared data were significantly correlated with urban planning and development activities. Vegetation cover reduced the amount of solar radiation reaching man-made surfaces; in addition, vegetation caused cooling of the surrounding air by releasing moisture through a process known as evapotranspiration. Therefore, proper planning and management regimes must be put in place as rapid urbanization proceeds to mitigate negative impact on the local climate and environment. Landsat TM thermal infrared data with low spatial resolution (60-120 m) was more appropriate for large-area thermal studies, however it might have utility to investigate complex urban thermal environments and determine LSTs for individual buildings. However, ASTER sensors with a TIR subsystem of 90 m spatial resolution were more accurate in resolving thermal patterns and LSTs. Landsat data show relatively high temperatures, which were more reliable relative to the actual measured data. ATLAS data (10 m resolution) or TABI (2 m resolution) performed well in community level studies in terms of thermal exchange [33,34], while ASTER and Landsat data were more appropriate for large area thermal studies. In this study, UHI calculated from the Landsat sensor was more consistent with ground truthed temperatures.