Diurnal and Seasonal Variation of Surface Urban Cool and Heat Islands in the Semi-arid City of Erbil, Iraq

The influence of land surface temperature (LST) makes the near-surface layer of the troposphere a key driver of urban climate. This paper assesses the temporal formation of the daytime Surface Urban Cool Island (SUCI) and night-time Surface Urban Heat Island (SUHI) effect in Erbil, Iraq, situated in a semi-arid climate region. LST retrievals from the Moderate Resolution Imaging Spectroradiometer (MODIS) Aqua and Terra and MODIS Normalized Difference Vegetation Index (NDVI) from January 2003 to December 2014 are analysed. The relationships of LST with NDVI and the Normalized Multi-band Drought Index (NMDI) are investigated in order to assess the influence of vegetation and moisture on the observed patterns of LST and the SUCI/SUHI. The results indicate that during the daytime, in summer, autumn and winter, densely built-up areas had lower LST acting as a SUCI compared to the non-urbanised area around the city. In contrast, at night-time, Erbil experienced higher LST and demonstrated a significant SUHI effect. The relationship between LST and NDVI is affected by seasonality and is strongly inverted during spring (r 2 = 0.73; p < 0.01). Contrary to previous studies of semi-arid cities, a SUCI was detected, not only in the morning, but also during the afternoon.


Introduction
Since land surface temperature moderates the lowest layer of the atmosphere, it plays an important role in the evolution of the urban climate.The comfort and health of urban residents is affected by surface temperature, which can influence internal building temperatures and energy exchange processes [1].The variation between the characteristics of surface cover in urban and rural areas, in terms of three-dimensional geometry of the built environment, heat absorption, construction materials, surface albedo and vegetation abundance cause different air and surface temperatures in a city relative to the surrounding less urbanised area.
The term "Surface Urban Cool Island" (SUCI) is defined as an urban area where lower surface temperatures are prevailing compared to the non-urbanised dry surroundings and is commonly found in arid and semi-arid regions [2,3].In contrast, the Urban Heat Island (UHI) effect refers to areas of higher air and surface temperatures experienced in a city compared to the rural surroundings [4][5][6].Satellite sensors can be used to investigate the surface temperature of cities and study the growth of SUHIs and their spatial variation [7].Rao [8] used satellite data of ITOS-1 to study the SUHI of the mid-Atlantic coast of the USA.Matson et al. [9] measured the UHI of the Midwestern and North-western United States during night-time with NOAA5 satellite data.Since then,

•
the determination of the seasonal and diurnal variation of the patterns of LST and SUCI/SUHI in a city with semi-arid climate; • the quantitative analysis of relationships between LST and NDVI and how they change between seasons

Study Area
Erbil is the capital of the Kurdistan Region and the main city in the north of Iraq.The city is situated 412 m above mean sea level [44] between 43 • 57'E to 44 • 03'E and 36 • 08'N to 36 • 14'N (Figure 1).Over the past two decades the infrastructure and population of the city have experienced extensive growth.In 2015, the population of the Erbil Governorate was estimated to be 1,530,722 inhabitants [45].Over the past 2 decades the annual growth rate of the population in Erbil is 2.9% [46].
The study area has a semi-arid continental climate and is classified as subtropical semi-arid (BSh) in the Koppen climate classification.It has a rainy and cool winter whilst the summer is dry and hot.The average annual precipitation is 380.26 ± 108.88 mm and the majority of it falls from December to March whilst the summer months are dry.The annual air temperature is 21.85 • C and July and August are the hottest months of the year; air temperature during these months can reach 49 • C [47].
Within the urban area, residential land use is the most dominant land use type comprising buildings made of concrete blocks.In the rural area, winter grain cultivation such as wheat and barley is most common and in general, agriculture in the study area is rain-fed rather than irrigated.
The citadel of Erbil is at the heart of the city and is the centre around which the surrounding urban pattern has grown in concentric rings [48].In the second half of the twentieth century, the third and fourth stages of the city expansion occurred.These newly developed areas are shaped in circular patterns around the city centre [49].Since 2003, the city changed rapidly and experienced abundant new architectural structures and large housing projects.As an example of the rapid growth, the number of quarters in Erbil reached 82 in 2012 while it was only eight in the 1950s.More tall buildings and small houses characterise this latest stage of city expansion.
new architectural structures and large housing projects.As an example of the rapid growth, the number of quarters in Erbil reached 82 in 2012 while it was only eight in the 1950s.More tall buildings and small houses characterise this latest stage of city expansion.The elevation of the city ranges from 319 m in the southwest to 878 m in the northeast of the study area.The mean elevation of the study area is 454 m.A general slope is present from south to northeast.Erbil is built-up on the Erbil Plain and has dark brown soils.Having good structure, large depth and rich in organic materials makes it one of the most fertile soils.The average soil depth of Erbil Plain is 140 cm.Soil texture in Erbil Governorate contains sand, loamy silt, loamy clay and silty clay [50].

Data
In this study, MODIS data from NASA's Terra and Aqua satellites are used.The resolution of LST data from MODIS is 1 km.Because it has four overpasses per day, it is better suited for investigating the diurnal variation of LST and the UHI effect than higher resolution sensors such as Landsat 8. Using MODIS Aqua (MYD11A2) and MODIS Terra (MOD11A2), an 8-day composite for a 12-year period was selected, from January 2003 to December 2014.Several categories of MODIS products from years 2003 to 2014 are used: The 8-day LST composite (MOD11A2 and MYD11A2), and 16-day NDVI composite (MOD13Q1) (Table 2).

Land Surface Temperature
LST is the radiative temperature of the land surface [51].It is influenced by albedo, vegetation cover and soil moisture [52].The "surface" can include snow and ice, bare soil, grass, or the roofs of The elevation of the city ranges from 319 m in the southwest to 878 m in the northeast of the study area.The mean elevation of the study area is 454 m.A general slope is present from south to northeast.Erbil is built-up on the Erbil Plain and has dark brown soils.Having good structure, large depth and rich in organic materials makes it one of the most fertile soils.The average soil depth of Erbil Plain is 140 cm.Soil texture in Erbil Governorate contains sand, loamy silt, loamy clay and silty clay [50].

Data
In this study, MODIS data from NASA's Terra and Aqua satellites are used.The resolution of LST data from MODIS is 1 km.Because it has four overpasses per day, it is better suited for investigating the diurnal variation of LST and the UHI effect than higher resolution sensors such as Landsat 8. Using MODIS Aqua (MYD11A2) and MODIS Terra (MOD11A2), an 8-day composite for a 12-year period was selected, from January 2003 to December 2014.Several categories of MODIS products from years 2003 to 2014 are used: The 8-day LST composite (MOD11A2 and MYD11A2), and 16-day NDVI composite (MOD13Q1) (Table 2).

Land Surface Temperature
LST is the radiative temperature of the land surface [51].It is influenced by albedo, vegetation cover and soil moisture [52].The "surface" can include snow and ice, bare soil, grass, or the roofs of buildings [53].Near-surface air temperature "is a measurement of the average kinetic energy of the air near the surface of the Earth" [54].Usually LST is measured by remote sensing whereas air Climate 2016, 4, 42 5 of 16 temperature is measured 1-2 m above the ground.Near-surface air temperature is a consequence of complex effects of the turbulent heat transports produced by nearby heated surfaces [55].
An 8-day composite product with quality assessment (QA) information from Terra (MOD11A2) and Aqua (MYD11A2) was used.These products were derived from the split-window method [56].Only LST images, of which 80% or higher pixels passed quality filter are used and lower-accuracy images are excluded.Pixels with zero values that represent NA values are ignored in the statistical analysis.Terra passes over the study area at approximately 11:00 and 22:00 local time, while Aqua images are taken around 02:00 and 13:00 local time.

NDVI
MOD13Q1 is a product from MODIS Terra, which contains a 16-day composite NDVI with 250 m spatial resolution and QA information.A total of 274 MOD13Q1 collection of 5 images from January 2003 to December 2014 were obtained.

Processing Data
MODIS products were downloaded from daac.ornl.govusing the WGS 84 ellipsoid model in UTM projection (zone 38N); associated with QA data for each image.The value of each pixel was multiplied by the scale factors for k LST = 0.02 and k NDVI = 0.0001.All spatial data were stacked and cropped to the study area.Statistics for subsets of the study area, stratified by the city and surrounding area were extracted and the images categorized into different seasons and fill gaps by the linear interpolation method.

LST
For determining the seasonal variation, the images are categorised into four seasons based on astronomical seasons of Erbil.Winter was defined between 22 December and 20 March, spring from 21 March to 20 June, summer from 21 June to 22 September and autumn from 23 September to 21 December.Average LST and standard deviation were calculated for twelve years for urban and the rural areas.The significance of the regression model was tested with Welsh's t-test.

NDVI
Due to the amount of vegetation varying according to the season and the amount of precipitation, the relationship between LST and NDVI depends on the season.The NDVI data for 12 years is divided into four seasons.Then NDVI and LST data are displayed as scatterplots and linear regression are used to quantify this relationship.

Soil Moisture
Soil moisture is one of the main drivers of LST especially in dry climates.NMDI uses multi near infrared and short wave infrared channels to estimate both soil and vegetation moisture.To calculate NMDI, Equation (1) was employed [28]: where DT is the diurnal temperature range (LSTday-LSTnight or Tmax-Tmin), DTmax is the maximum diurnal temperature (constant: 60 • C), fsm is scaled between zero and one.

Identification of the Surface UCI and UHI Intensity
There are different methods of measuring UHII.Here, a 10-km buffer zone surrounding the city of Erbil is selected as a "rural" (non-urban) reference area.The city limit used to define urban areas and the shapefile of the city was created in ArcGIS based on high-resolution satellite imagery.The city boundary delimits predominantly urban areas, which contain some pixels of other LULC classes while the buffer principally represents rural areas (92%), which contain only a small number (1.3%) of urban pixels.The intensity of SUCI and SUHI is calculated using Equation (3) [58].
where SUHI = Surface Urban Heat Island; SUCI = Surface Urban Cool Island; LST u = mean LST of urban area; LST r = mean LST of rural buffer area.A positive value represents a SUHI situation, while a negative result represents SUCI.

Variation between Urban and Rural Areas
The means and standard deviations of LST, NDVI and NMDI were calculated for both the city and rural areas.The percentage of low quality LST images was low at night (Terra 2.5% and Aqua 2.4%) while it was higher for daytime data (Terra 12% and Aqua 20%).

Relationship of LST with NDVI, fsm and NMDI
The relationships between LST and the various influencing factors are presented as scatterplots.In addition, assuming linear relationships, linear regression is used.The significance of the regression model was tested using the Welch method.

Diurnal Variation of LST and SUCI/SUHI
There is considerable variation of LST and SUCI/SUHI between the city and rural areas due to the variation of heat absorption and radiation during the day and night, respectively.Figure 2 indicates that generally during the day the city experiences a SUCI effect.The magnitude of the SUCI intensity tends to be greatest (-1.78 ± 1.68 • C) in the morning (11:00) although there is some seasonal variation.However, at night the situation is reversed and LST in the city is higher than the rural surroundings and a SUHI is experienced.The SUHI intensity attains its highest value, around 2.67 ± 0.72 • C, in the early night (~22:00 pass); this drops off during the late night.In terms of the most extreme SUCI and SUHI during these 12 years, in the morning (11:00) SUCI was −5.37 • C, and in the afternoon it was −5.32 • C. The highest SUHI of early night was 4.59 • C and it was 4.28 • C in the late night.
Figure 3 demonstrates the variation between the average LST and SUCI/SUHI over all seasons between the city of Erbil and the rural surroundings for the four overpasses of the Terra and Aqua satellite during 2014.At ~11:00, the LST of the city was 30.92 ± 0.98 • C, whilst the LST of the rural area was 32.44 ± 1.46 • C. The city experienced a daytime SUCI averaging −1.5 • C, decreasing to −1.2 • C by the afternoon.In contrast, at night the surface of the rural environment radiates more heat than the surface of the city.Consequently, the average LST of the city in 2014 at ~22:00 was 17.86 ± 1.12 • C whilst in the rural area it was 14.6 • C ± 1.09 • C. The SUHI intensity was 3.

Seasonal Variation of LST and SUCI/SUHI
The city exhibits a night-time SUHI with varying intensity depending on the season (Figures 4 and 5).In the spring and summer, the intensity of the SUHI in the study area is higher than during autumn and winter.In contrast, the intensity of the daytime SUCI (negative SUHI) is higher in autumn and summer; the incidence of a daytime SUCI decreases in winter and spring.Spring is the growing season for grasses and crops in the rural area, which cools the local environment via reduced heat absorption and greater evapotranspiration and latent heat flux during this season.
In autumn, during the daytime LST in the rural area is higher than the LST of Erbil and it exhibits a SUCI.In contrast, at night the urban areas experienced a SUHI (Table 3).In spring during the daytime, the city typically exhibited a mild SUCI (at ~13:00 the SUCII was −0.05 ± 0.64 • C).During the night-time SUHI was 2.82 • C ± 1.05.In summer, the average of LST in the city at ~11:00 was 46.98 ± 1.43 • C and the LST of the rural area 49.76 ± 1.30 • C. The SUCI intensity was −2.78 • C at 11:00, whilst during the night-time the city exhibited a SUHI intensity of 3.05 • C.This is particularly problematic from a public health perspective; as it is reducing thermal comfort while increasing the risk of heatstroke, exhaustion, and heat-related mortality.In addition, it increases the effect of heat waves on inhabitants who are sensitive to high temperatures, such as the elderly, children and the sick [59][60][61][62].

Seasonal Variation of LST and SUCI/SUHI
The city exhibits a night-time SUHI with varying intensity depending on the season (Figures 4  and 5).In the spring and summer, the intensity of the SUHI in the study area is higher than during autumn and winter.In contrast, the intensity of the daytime SUCI (negative SUHI) is higher in autumn and summer; the incidence of a daytime SUCI decreases in winter and spring.Spring is the growing season for grasses and crops in the rural area, which cools the local environment via reduced heat absorption and greater evapotranspiration and latent heat flux during this season.
In autumn, during the daytime LST in the rural area is higher than the LST of Erbil and it exhibits a SUCI.In contrast, at night the urban areas experienced a SUHI (Table 3).In spring during the daytime, the city typically exhibited a mild SUCI (at ~13:00 the SUCII was −0.05 ± 0.64 °C).During the night-time SUHI was 2.82 °C ± 1.05.In summer, the average of LST in the city at ~11:00 was 46.98 ± 1.43 °C and the LST of the rural area 49.76 ± 1.30 °C.The SUCI intensity was −2.78 °C at 11:00, whilst during the night-time the city exhibited a SUHI intensity of 3.05 °C.This is particularly problematic from a public health perspective; as it is reducing thermal comfort while increasing the risk of heatstroke, exhaustion, and heat-related mortality.In addition, it increases the effect of heat waves on inhabitants who are sensitive to high temperatures, such as the elderly, children and the sick [59][60][61][62].In general, during the last 12 years the NDVI of the rural area was higher than the NDVI of Erbil.The average NDVI of the rural area was 0.22 ± 0.11 whereas the average NDVI of the city was 0.16 ± 0.05 (Table 4).Since vegetation in the study area is strongly dependent upon precipitation, the value of NDVI varies between years depending on the amount of rain.The peak NDVI value was recorded in 2003 when it was 0.57 in the rural area and 0.37 in the city.The total precipitation in that year was high (496.3mm) compared to the annual average (386.2mm).The lowest values of NDVI (0.12 ± 0.02) coincided with reduced precipitation in 2008 (297.5 mm).Moreover, the value of NDVI demonstrates seasonal variation.High values (0.31 ± 0.07 in the rural and 0.19 ± 0.05 in the city) occur in spring when precipitation is high and temperatures start to increase.In contrast, during the dry season (summer) the seasonal average of NDVI dropped to the lowest value with 0.14 ± 0.01 in the rural area and 0.13 ± 0.01 in the city.

NMDI
In general, the NMDI in the city is higher than the rural surroundings and it displays annual and seasonal variation.The average NMDI in the city during 2003 to 2014 was 0.69 ± 0.03 and it was 0.60 ± 0.12 in the rural area (Table 4).The high moisture in the city compare to the rural areas is a result of irrigation activity of city dwellers and low precipitation in the area.

Soil Moisture Constraint
In general, the fsm in the city is higher than the rural surroundings and it displays annual and seasonal variation.The average fsm in the city during 2003 to 2014 was 0.35 ± 0.11 and it was 0.27 ± 0.11 in the rural area.The peak fsm of the city was 0.77 in December 2013 and the lowest value was 0.18 in June 2006.Moreover, the highest value of the rural area was 0.82 in December 2013 and the lowest value (0.12) was recorded in June 2006.In general, fsm is high in the winter and autumn when it is between 0.4 to 0.7, while in the summer it reduced to between 0.1 and 0.3.

Relationships between LST and NDVI, fsm and NMDI
A weak but statistically significant inverse relationship (r 2 = 0.18; p < 0.01) between surface temperature and NDVI was found (Figure 6).As a result of seasonal vegetation phenology, the relationship of LST with NDVI varies by season.In the wet spring, vegetation in the study area increases significantly.During this season the highest inverse relationship (r 2 = 0.73, p < 0.01) is found between LST and NDVI (Figure 7).In the dry summer, NDVI values drop to the lowest level and the relationship with LST becomes very weak (r 2 = 0.02, p = 0.22).At the start of the wet season in autumn, NDVI increases slightly and the relationship with LST is stronger (r 2 = 0.34, p < 0.01).However, in the winter the relationship between NDVI and LST during daytime switches to a positive relationship (r 2 = 0.37, p < 0.01), especially when LST is low.As an extension of our previous results [20], applying MODIS products to establish the relationship of LST and NDVI for each season in this research adds further insight into the nature of this relationship.
A moderate inverse relationship (r 2 = 0.38; p < 0.01) is found between LST and NMDI.Inside the city the NMDI is 0.35 ± 0.11 and the mean LST at ~11:00 is 30.92 ± 0.98 • C, whilst in the 10 km buffer of the city the NMDI decreases to 0.27 ± 0.11.Moreover, the mean LST increased to 32.44 ± 1.46 A moderate inverse relationship (r 2 = 0.38; p < 0.01) is found between LST and NMDI.Inside the city the NMDI is 0.35 ± 0.11 and the mean LST at ~11:00 is 30.92 ± 0.98 °C, whilst in the 10 km buffer of the city the NMDI decreases to 0.27 ± 0.11.Moreover, the mean LST increased to 32.44 ± 1.46 °C.A strong statistically significant inverse relationship (r 2 = 0.81; p < 0.01) is ascertained between LST and fsm.Areas with high fsm have low LST; for inside the city border the fsm is 0.35 ± 0.11 and the mean LST at ~11:00 is 30.92 ± 0.98 °C, whilst in the 10 km buffer of the city the fsm decreases to 0.27 ± 0.11.Moreover, the mean LST increased to 32.44 ± 1.46 °C.A moderate inverse relationship (r 2 = 0.38; p < 0.01) is found between LST and NMDI.Inside the city the NMDI is 0.35 ± 0.11 and the mean LST at ~11:00 is 30.92 ± 0.98 °C, whilst in the 10 km buffer of the city the NMDI decreases to 0.27 ± 0.11.Moreover, the mean LST increased to 32.44 ± 1.46 °C.A strong statistically significant inverse relationship (r 2 = 0.81; p < 0.01) is ascertained between LST and fsm.Areas with high fsm have low LST; for inside the city border the fsm is 0.35 ± 0.11 and the mean LST at ~11:00 is 30.92 ± 0.98 °C, whilst in the 10 km buffer of the city the fsm decreases to 0.27 ± 0.11.Moreover, the mean LST increased to 32.44 ± 1.46 °C.A strong statistically significant inverse relationship (r 2 = 0.81; p < 0.01) is ascertained between LST and fsm.Areas with high fsm have low LST; for inside the city border the fsm is 0.35 ± 0.11 and the mean LST at ~11:00 is 30.92 ± 0.98 • C, whilst in the 10 km buffer of the city the fsm decreases to 0.27 ± 0.11.Moreover, the mean LST increased to 32.44 ± 1.46 • C.

Discussion
The results presented here demonstrate a different dynamic behaviour in the semi-arid climate of Erbil, where the strongest SUHI occurs at night, confirming earlier results by Cheval et al. [63], and turning into a SUCI during the day.The results indicate the existence of SUCI during different times of the day and not only in the morning as stated in the literature [64][65][66][67].
The seasonal variation of SUCI/SUHI indicates that the night-time SUHI is stronger in the spring and summer than the autumn and winter.This is in agreement with et al. [68], who found that SUHI is higher in summer and spring and weaker in the winter and autumn.Similarly, Chow and Roth [38] observed the highest UHI intensity during the dry season and lowest intensities during the wet season.On the other hand, daytime SUCI intensity in the autumn and summer is higher than in winter and turns neutral in the spring.In general, during summer the UHI effect is more important to the urban residential areas.This is due to the more extreme thermal conditions, which increase illnesses and deaths related to heat [69].
During the study period, NDVI of the rural areas was on average 0.06 higher than the NDVI of the urban areas.The result for Erbil as a semi-arid area disagrees with results from arid studies, for instance for Phoenix in the USA [11,12].Nevertheless, because the vegetation in the area mainly depends on precipitation, this difference of NDVI between the city and rural areas increases significantly in spring and decreases during the dry season.In March 2010 and 2013, this difference reached 0.31 while in August and September 2013 no differences were observed.
This study also examined the relationship of LST and NDVI by season for twelve years of MODIS data.The general relationship for the 12-year period was weakly negative.In general vegetation in the study area depends on the precipitation, rather than irrigation.Therefore, the majority of croplands are dry in the summer time.Therefore, this seasonality in vegetation lead to seasonal variation in relationship between LST and NDVI.The results indicate that this relationship changes significantly with season, confirming the results of Yuan and Bauer [27] and Sun and Kafatos [21].In the summer, it is very weak.The peak of the inverse relationship is located in the spring within the period of vegetation growth, meaning that in the area of high vegetation greenness LST decreases.On the contrary, during the winter, which is a relatively cold and wet season in the study area, this relationship turns positive and LST rises with increasing NDVI.This positive relationship is consistent with Sun and Kafatos [21] during the winter, Haashemi et al. [42] and the results of Schultz and Halpert [70] in the high and mid-latitudes.The reason for this positive relationship is potentially related to vegetation water content which retains heat during the cold period; the decrease of photosynthesis may also have an influence [42].

Conclusions
During the twelve years of MODIS observations, the greatest magnitude of SUCI occurs in the morning.However, the greatest magnitude of SUHI emerged at early night.Diurnal variation shows that in the morning, the average LST of the city was lower than the rural and the city experienced a SUCI during the day with higher intensity in the morning than the afternoon.Average LST of the city in 2014 at night was higher than the LST of the rural area.
In terms of seasonal variation, during the night in all seasons the city exhibits a SUHI with different intensity.In the spring and summer, the intensity of the SUHI in the study area is higher than in autumn and winter.During the day-time the SUCII is higher in autumn and summer and decreases in winter, changing to a weak SUCI in the spring, which is the growing season of grasses and crops in rural areas.For this reason, the rural area absorbs heat more slowly in spring.In summer, the hottest and driest season, at 13:00 the average LST of the city is 50.36 ± 1.35 The examination of the key factors that influence LST in the study area shows that during the last twelve years NDVI in rural areas was higher than in Erbil city.In contrast, NMDI in the city is higher than in the rural area and shows annual and seasonal variation.The reason is that in the urban areas people wet the surface by irrigation of urban green spaces, mainly in residential areas and urban parks [71], while in the rural areas agriculture is rain-fed rather than irrigated.In the study area, 93% of croplands are rain-fed [72], therefore, the vast majority of croplands are very dry in the summertime.
A strong significant inverse relationship is ascertained between LST and fsm.A weak but significant inverse relationship between LST and NDVI was found.This relationship varies in magnitude with season.During spring, the strongest inverse relationship is found between LST and NDVI.This relationship weakened in summer and was moderately strong (inverse) in autumn.However, in winter this relationship turns into a moderate positive relationship.A moderate (significant) inverse relationship is found between surface temperature and NMDI.
MODIS LST was effective for studying the temporal variation of SUHI/SUHI (e.g., diurnal and seasonal).In addition, using QA data in the quality control process helped to achieve results that are more accurate.Using R scripts in addition to ENVI Software for processing MODIS products in this study was time-saving.Applying MODIS products for the relationship of LST and NDVI investigated for each season in this research adds details of the nature of the relationship between them, and explains the reason of argument on this relationship in literature.
A SUCI was detected in the semi-arid city of Erbil during different times of the day and not only in the morning as stated in previous literature.This study has focused on the temporal variation of the SUCI/SUHI.Spatial variations, based on land cover and land use in Erbil, have been investigated previously [20].Future studies on the surface temperature of cities in arid and semi-arid environments are required to confirm the results presented here.Other research should focus on long-term air temperature measurements for UCI/UHI exploration.This study was limited to assessing a subset of all possible factors influencing the UHI effect while future studies should address this limitation and explore more factors that affect LST, such as urban geometry, albedo and building components.

Figure 1 .
Figure 1.Map of Iraq, the Kurdistan Region and true colour Landsat image of the study area.(a) Iraq's map; (b) Kurdistan Regional Government; (c) the study area.

Figure 1 .
Figure 1.Map of Iraq, the Kurdistan Region and true colour Landsat image of the study area.(a) Iraq's map; (b) Kurdistan Regional Government; (c) the study area.

Figure 3 .
Figure 3.Diurnal variation of average 2014 Land Surface Temperature and Surface Urban Cool Island/Surface Urban Heat Island during the day and night.

Figure 3 .
Figure 3.Diurnal variation of average 2014 Land Surface Temperature and Surface Urban Cool Island/Surface Urban Heat Island during the day and night.

Figure 3 .
Figure 3.Diurnal variation of average 2014 Land Surface Temperature and Surface Urban Cool Island/Surface Urban Heat Island during the day and night.

Figure 4 .
Figure 4. Diurnal variation of average Land Surface Temperature in different seasons over 12 years from Aqua and Terra MODIS and Surface Urban Cool Island/Surface Urban Heat Island.

Figure 4 .
Figure 4. Diurnal variation of average Land Surface Temperature in different seasons over 12 years from Aqua and Terra MODIS and Surface Urban Cool Island/Surface Urban Heat Island.

Figure 5 .
Figure 5. Seasonal variation of average Land Surface Temperature over 12 years, based on MODIS Aqua, the green line is the boundary of Erbil city.

Figure 5 .
Figure 5. Seasonal variation of average Land Surface Temperature over 12 years, based on MODIS Aqua, the green line is the boundary of Erbil city.

Figure 6 .
Figure 6.Relationship of Land Surface Temperature with NDVI, fsm and NMDI, from MODIS data.

Figure 6 .
Figure 6.Relationship of Land Surface Temperature with NDVI, fsm and NMDI, from MODIS data.

Figure 6 .
Figure 6.Relationship of Land Surface Temperature with NDVI, fsm and NMDI, from MODIS data.

Table 1 .
A summary of published studies of Urban Heat Island Intensity (UHII)/Urban Cool Island Intensity (UCII) around the world.

Table 2 .
Summary of MODIS image products used in this study.

Table 2 .
Summary of MODIS image products used in this study.

Table 3 .
Diurnal variation of average Land Surface Temperature and ∆ LST in different seasons over 12 years from Aqua and Terra MODIS.

Table 3 .
Diurnal variation of average Land Surface Temperature and ∆ LST in different seasons over 12 years from Aqua and Terra MODIS.

Table 4 .
Min, max and mean of NDVI, NMDI and fsm of Erbil city and rural area of 12 years from MODIS data.
• C and LST of the rural areas is 52.84 ± 1.01 • C and SUCII is −2.48 • C.During the night-time, the city exhibits a SUHI with SUHII around 3.05 • C.