# A Quantitative Assessment of Surface Urban Heat Islands Using Satellite Multitemporal Data over Abeokuta, Nigeria

Authors: K. A. Ishola; E. C. Okogbue; O. E. Adeyeri
Journal: International Journal of Atmospheric Sciences
 (2016)

---

## Abstract

The fast urban expansion has led to the transformation of the natural landscape into anthropogenic surfaces. The city of Abeokuta, for instance, is located in a region experiencing rapid urbanization, which has produced a remarkable effect on the surface thermal response. This effect significantly influences urban internal microclimatology on a regional scale. In this study, the surface temperatures and land cover types retrieved from Landsat TM and ETM+ images of Abeokuta city for 1984, 2003, and 2014 were analyzed. A quantitative approach was used to assess surface urban heat islands through the relationships among surface temperature and land cover types. Results showed that impervious surface areas were found to be correlated positively with high temperatures. Conversely, vegetated areas and bare surfaces correlated positively with mid temperature zones. This study found that areas with increasing impervious surfaces will accelerate LST rise and consequently lead to increasing effect of surface urban heat islands. These findings pose a major challenge to urban planners. However, the study would help to quantify the impacts of different scenarios (e.g., vegetation loss to accommodate urban growth) on LST and consequently to devise appropriate policy measures.


---

## Body

## 1. Introduction

Rapid rural-urban migration and population growth in the cities over the past decades have led to an unprecedented rate of urbanization globally, and it is likely to continue in the subsequent decades [1]. During the process of urbanization, a direct environmental consequence is the modification of land surfaces which results in the alteration of its physical properties including soil moisture, material heat capacity, and surface reflectivity and emissivity. This leads to change in land surface temperature and the decrease of evapotranspiration [2]. Consequently, the impacts result in surface urban microclimate variations locally. Over the last three decades, the rapid development of remote sensing techniques provided appropriate data for studying the regional surface characteristics of an area. Land surface parameters such as land surface temperature (LST), the normalized difference vegetation index (NDVI), and surface albedo dynamics were derived from data obtained from these sensors [3]. A quantitative relationship between LST and land use/land cover using both empirical and dynamic numerical methods has been documented [4, 5].The consistent expansion of the urban areas is evident in the transformation of natural landscape into anthropogenic impervious surfaces. In turn, increased impervious surface areas also lead to greater temperatures in urban areas than in nonurban areas resulting in urban heat effect [6]. The proportion of impervious surfaces has been reported to be a good indicator for the monitoring of the urban thermal response. A positive correlation between the proportion of impervious surfaces and land surface temperatures was identified by Yuan and Bauer [7].However, a consistent and quantitative detailed study to assess the changes in surface thermal response due to urban sprawl in Abeokuta city is very rare. This is probably due to the challenges of lack of consistent and repetitive data coverage on the changes in land surface physical conditions over this area. Hence, this study is aimed at assessing quantitatively the changes in surface urban thermal response due to various heterogeneous surface areas. The study will evaluate the dynamics of LST with respect to land cover changes and also determine the separate contributions of land cover types to surface thermal response change. Further, it will assist urban planners and policy makers in developing measures to curbing the environmental risks posed by this impact.


## 2. Materials and Methods

This study looks at the city of Abeokuta in southwest Nigeria (Figure1) due to its rapid human and industrial growth and development. It lies between latitudes 7°10′0′′N and 7°15′0′′N and longitudes 3°17′0′′E and 3°26′0′′E. Annual rainfall is about 963 mm and the temperature is usually between 26 and 28°C. The study employs cloudless Landsat 5 TM and 7 ETM+ images from 18 December 1984, 29 January 2003, and 29 December 2014. All bands from the images were applied in the analysis. Bands 1–5 and 7 have a spatial resolution of 30 m, and the thermal infrared band (band 6) has a spatial resolution of 120 m.Figure 1
Map showing the study area.


### 2.1. Land Cover Classifications

The maximum likelihood scheme was adopted for this work. Based on the prior knowledge of the land cover types in the study area, a chosen color composite (RGB = 432) was used for digitizing polygons around each training site for similar land cover. Then a unique identifier was assigned to each known land cover type. Afterwards, the statistical characterizations (i.e., signatures) of each land cover class were developed. Each composed image was ordered into 4 area classes: water, vegetation, impervious surface, and bare soil as shown in Figure2.Figure 2
Maps of land cover types of Abeokuta city for (a) 1984, (b) 2003, and (c) 2014.


### 2.2. Image Processing

The images acquired from the United States Geological Survey (USGS) database were radiometrically corrected before the calculation to avoid a dataset-specific result. The radiometric correction employed the algorithm of Chander and Markham [8] with the addition of an atmospheric correction.This converted the digital number (DN) of the raw image to at-satellite reflectance. Thermal band 6 was converted to at-satellite temperature and then corrected in an emissive way to LST. The LST was retrieved using the following models:(1)Tb=K2ln⁡K1/Lλ+1,LST=Tbε00.25,where Lλ equals spectral radiance at the sensor’s aperture (w/m2∗ster∗μm); Tb equals brightness temperature (K); K1 and K2 equal calibration constants. K1 equals 607.76/666.09 and K2 equals 1260.56/1282.71 (Wm−2 sr−1 μm−1) for Landsat 5/7; ε0 equals surface broadband emissivity; LST equals land surface temperature (K).


## 2.1. Land Cover Classifications

The maximum likelihood scheme was adopted for this work. Based on the prior knowledge of the land cover types in the study area, a chosen color composite (RGB = 432) was used for digitizing polygons around each training site for similar land cover. Then a unique identifier was assigned to each known land cover type. Afterwards, the statistical characterizations (i.e., signatures) of each land cover class were developed. Each composed image was ordered into 4 area classes: water, vegetation, impervious surface, and bare soil as shown in Figure2.Figure 2
Maps of land cover types of Abeokuta city for (a) 1984, (b) 2003, and (c) 2014.


## 2.2. Image Processing

The images acquired from the United States Geological Survey (USGS) database were radiometrically corrected before the calculation to avoid a dataset-specific result. The radiometric correction employed the algorithm of Chander and Markham [8] with the addition of an atmospheric correction.This converted the digital number (DN) of the raw image to at-satellite reflectance. Thermal band 6 was converted to at-satellite temperature and then corrected in an emissive way to LST. The LST was retrieved using the following models:(1)Tb=K2ln⁡K1/Lλ+1,LST=Tbε00.25,where Lλ equals spectral radiance at the sensor’s aperture (w/m2∗ster∗μm); Tb equals brightness temperature (K); K1 and K2 equal calibration constants. K1 equals 607.76/666.09 and K2 equals 1260.56/1282.71 (Wm−2 sr−1 μm−1) for Landsat 5/7; ε0 equals surface broadband emissivity; LST equals land surface temperature (K).


## 3. Results and Discussion

Surface temperature has been established as a major indicator of the presence of surface urban heat island in cities and urban areas. The relationship between land cover and surface temperature was analyzed for each year by producing images that provided a visual resource for analyzing the intensity and spatiality of LST change.


### 3.1. Spatiotemporal Dynamics of LST

The spatial-temporal distributions of LST changes over Abeokuta in the three examined dates (Figure3) show that the LSTs ranged from 25.0°C to 37.9°C, 18.6°C to 34.4°C, and 23.7°C to 40.2°C in 1984, 2003, and 2014, respectively.Figure 3
Maps of Abeokuta city showing the thermal zones for (a) 1984, (b) 2003, and (c) 2014.Figure4 shows that no areas in Abeokuta experienced an extreme temperature (≥41°C) as at the time satellite overpassed in all the periods. In 1984, a large part of the area (67.7% in all) fell within the higher temperature zones (≥30°C) and other areas (32.3%) fell into mid temperature zones (25°C–30°C). This high surface temperature might be due to high insolation and few atmospheric compositions (like the cloud cover) as at the time satellite overpassed, large distribution of imperviousness, or geographical relief of the area.Figure 4
Changing pattern of surface temperature zones.But in 2003, about 66.3% proportions of the surface area were dominated by mid-surface temperature zones (25°C–30°C). Other surface areas (29.1% and 4.6%) fell in the lower (<25°C) and higher (≥30°C) temperature zones, respectively. Moreover, a higher LST zone (≥30°C) was also observed in majority (81.6%) of the area in 2014. A noticeable surface area proportion (0.3% and 18.2%) were found in the lower (<25°C) and mid (25°C–30°C) temperature zones, respectively (Figure4).The mean surface temperature variation over different land cover types (Figure5) showed that the majority of the bare soil and impervious surface areas of Abeokuta had temperatures between 32°C and 35°C in both 1984 and 2014. In 2003, the temperature was cooler between 26°C and 29°C in both land cover types.Figure 5
Mean LST variations over different land cover types.However, vegetation and water bodies have similar temperature variations between 28 and 30°C in 1984 and 2014. This could be attributed to less urbanization process in the former year and increase in population and infrastructural development which has led to the extension of the urban core (higher urbanization) into the adjoining rural lands in the later year. Moreover, the mean LSTs of vegetation cover and water bodies were taken to between 24°C and 26°C in 2003 which might be due to the presence of cloud cover and aerosol loads that reduced the incoming shortwave radiation. The high surface temperature even over vegetation cover in 2014 was as a result of more surface modifications, anthropogenic activities, little vegetation, and consequently little/no evapotranspiration. In addition, the magnitude of LST change is almost the same (5-6°C/16–19 years) over both impervious surface and bare soil (Figure5).


### 3.2. Quantitative Relationship between LST and Land Covers

To show a quantitative relationship between land covers and temperature change, the study used the Pearson Correlation to calculate the correlation matrix (Table1). The area size data extracted from the supervised classification for each year was analyzed using the R software environment for statistical computing [4]. The higher positive correlations were found between the ranges of temperature <25°C and 25–30°C and high temperature of 30–33°C and ≥33°C. The highest negative correlations were found between the temperature range of <25°C and 30–33°C and the higher temperature ranges (25–30°C and ≥33°C).Table 1
Correlations among LSTs and land cover types.
ISA
Bare
Veg.
Water
<25°C
25–30°C
30–33°C
≥33°C
ISA
1.00
−0.99
−0.99
0.92
−0.22
−0.50
−0.10
0.78
Bare
1.00
0.97
−0.96
0.34
0.60
−0.02
−0.85
Veg.
1.00
−0.86
0.10
0.39
0.23
−0.70
Water
1.00
−0.59
−0.80
0.30
0.97
<25°C
1.00
0.96
−0.95
−0.78
25–30°C
1.00
−0.81
−0.93
30–33°C
1.00
0.54
≥33°C
1.00A high negative correlation was also found between the 25–30°C and 30–33°C temperature range. The higher positive correlation between land covers was found between areas of vegetation and bare soil and high correlation between impervious surface area (ISA) and water bodies. The highest negative correlation was between ISA and vegetated and ISA and bare areas.To relate the land cover patterns to surface urban heat, the correlations between the classes of land cover and LST were analyzed (Table1). The best highest positive correlations were found between water, ISA areas, and temperatures ≥33°C.These findings were consistent with [6, 7]. Other positive correlations were found between bare areas and temperatures between 25 and 30°C. However, these correlations were not as high compared to what was found in ISA and water areas.The highest positive correlations with low temperatures were found in the classes of bare soil and vegetation. Conversely, vegetation and bare areas had a high negative correlation with high temperatures.As shown in Figures2 and 4 and Table 1, changes in land cover patterns are highly correlated to changes in LSTs and the surface urban heat island (SUHI) effect.The impervious surface and water areas were found to contribute positively to LST rise, averagely through the three examined dates, with percentage contributions of moderate 60% to high 90% (Figure6). Consequently, this might lead to SUHI effect because the growth of impervious surfaces coincided with an increase in temperatures. However, vegetation and bare areas contributed negatively to SUHI effect.Figure 6
Separate contributions of land cover types to surface temperature change.


## 3.1. Spatiotemporal Dynamics of LST

The spatial-temporal distributions of LST changes over Abeokuta in the three examined dates (Figure3) show that the LSTs ranged from 25.0°C to 37.9°C, 18.6°C to 34.4°C, and 23.7°C to 40.2°C in 1984, 2003, and 2014, respectively.Figure 3
Maps of Abeokuta city showing the thermal zones for (a) 1984, (b) 2003, and (c) 2014.Figure4 shows that no areas in Abeokuta experienced an extreme temperature (≥41°C) as at the time satellite overpassed in all the periods. In 1984, a large part of the area (67.7% in all) fell within the higher temperature zones (≥30°C) and other areas (32.3%) fell into mid temperature zones (25°C–30°C). This high surface temperature might be due to high insolation and few atmospheric compositions (like the cloud cover) as at the time satellite overpassed, large distribution of imperviousness, or geographical relief of the area.Figure 4
Changing pattern of surface temperature zones.But in 2003, about 66.3% proportions of the surface area were dominated by mid-surface temperature zones (25°C–30°C). Other surface areas (29.1% and 4.6%) fell in the lower (<25°C) and higher (≥30°C) temperature zones, respectively. Moreover, a higher LST zone (≥30°C) was also observed in majority (81.6%) of the area in 2014. A noticeable surface area proportion (0.3% and 18.2%) were found in the lower (<25°C) and mid (25°C–30°C) temperature zones, respectively (Figure4).The mean surface temperature variation over different land cover types (Figure5) showed that the majority of the bare soil and impervious surface areas of Abeokuta had temperatures between 32°C and 35°C in both 1984 and 2014. In 2003, the temperature was cooler between 26°C and 29°C in both land cover types.Figure 5
Mean LST variations over different land cover types.However, vegetation and water bodies have similar temperature variations between 28 and 30°C in 1984 and 2014. This could be attributed to less urbanization process in the former year and increase in population and infrastructural development which has led to the extension of the urban core (higher urbanization) into the adjoining rural lands in the later year. Moreover, the mean LSTs of vegetation cover and water bodies were taken to between 24°C and 26°C in 2003 which might be due to the presence of cloud cover and aerosol loads that reduced the incoming shortwave radiation. The high surface temperature even over vegetation cover in 2014 was as a result of more surface modifications, anthropogenic activities, little vegetation, and consequently little/no evapotranspiration. In addition, the magnitude of LST change is almost the same (5-6°C/16–19 years) over both impervious surface and bare soil (Figure5).


## 3.2. Quantitative Relationship between LST and Land Covers

To show a quantitative relationship between land covers and temperature change, the study used the Pearson Correlation to calculate the correlation matrix (Table1). The area size data extracted from the supervised classification for each year was analyzed using the R software environment for statistical computing [4]. The higher positive correlations were found between the ranges of temperature <25°C and 25–30°C and high temperature of 30–33°C and ≥33°C. The highest negative correlations were found between the temperature range of <25°C and 30–33°C and the higher temperature ranges (25–30°C and ≥33°C).Table 1
Correlations among LSTs and land cover types.
ISA
Bare
Veg.
Water
<25°C
25–30°C
30–33°C
≥33°C
ISA
1.00
−0.99
−0.99
0.92
−0.22
−0.50
−0.10
0.78
Bare
1.00
0.97
−0.96
0.34
0.60
−0.02
−0.85
Veg.
1.00
−0.86
0.10
0.39
0.23
−0.70
Water
1.00
−0.59
−0.80
0.30
0.97
<25°C
1.00
0.96
−0.95
−0.78
25–30°C
1.00
−0.81
−0.93
30–33°C
1.00
0.54
≥33°C
1.00A high negative correlation was also found between the 25–30°C and 30–33°C temperature range. The higher positive correlation between land covers was found between areas of vegetation and bare soil and high correlation between impervious surface area (ISA) and water bodies. The highest negative correlation was between ISA and vegetated and ISA and bare areas.To relate the land cover patterns to surface urban heat, the correlations between the classes of land cover and LST were analyzed (Table1). The best highest positive correlations were found between water, ISA areas, and temperatures ≥33°C.These findings were consistent with [6, 7]. Other positive correlations were found between bare areas and temperatures between 25 and 30°C. However, these correlations were not as high compared to what was found in ISA and water areas.The highest positive correlations with low temperatures were found in the classes of bare soil and vegetation. Conversely, vegetation and bare areas had a high negative correlation with high temperatures.As shown in Figures2 and 4 and Table 1, changes in land cover patterns are highly correlated to changes in LSTs and the surface urban heat island (SUHI) effect.The impervious surface and water areas were found to contribute positively to LST rise, averagely through the three examined dates, with percentage contributions of moderate 60% to high 90% (Figure6). Consequently, this might lead to SUHI effect because the growth of impervious surfaces coincided with an increase in temperatures. However, vegetation and bare areas contributed negatively to SUHI effect.Figure 6
Separate contributions of land cover types to surface temperature change.


## 4. Conclusion

In this paper, qualitative and quantitative analyses were performed to study the relationship between land cover changes and SUHIs. The local-scale analyses focused on the city of Abeokuta, and the following conclusions were made: (1) the distribution of surface heat in the urban area of Abeokuta has changed from a mixed pattern to urban heat as the impervious areas expanded rapidly in the city since almost 13% of the area of the city changed to high LST (≥33°C); (2) land cover patterns and changes contributed to the variations in the microclimate and affected the SUHI intensity primarily through an extension of urban core to the rural areas, soil impaction, and deforestation; (3) SUHIs were proportionally related to water bodies (r=0.97) because the slight expansion of the major water body (River Lafenwa) located in the city centre brings about increase in the surface temperature and ISA (r=0.78) which includes buildings, roads, rock surfaces, parking lots, and other paved surfaces; (4) quantitative analysis among temperature and land cover patterns showed that increasing impervious surface areas will result in LST rise and consequently SUHIs. However, vegetation reduces the LST. Thus, future city planning should focus on urban greening. (5) It is suggested that representative land cover indices can be used together as independent factors in a multiple regression equation model in order to evaluate the effectiveness of land covers as indicators to LST change in a more robust way. Moreover, the change in LST does not fully explain the surface urban heat islands phenomenon as it has been used in this study, but LST is considered one of the important land surface components contributing to surface urban heat. Thus, the calculation of surface urban heat requires data from a comparative geography (e.g., rural versus urban), which was beyond the scope of the study.



---

*Source: 3170789-2016-05-03.xml*