Comparative Analysis of Responses of Land Surface Temperature to Long-Term Land Use / Cover Changes between a Coastal and Inland City : A Case of Freetown and Bo Town in Sierra Leone

Urban growth and its associated expansion of built-up areas are expected to continue through to the twenty second century and at a faster pace in developing countries. This has the potential to increase thermal discomfort and heat-related distress. There is thus a need to monitor growth patterns, especially in resource constrained countries such as Africa, where few studies have so far been conducted. In view of this, this study compares urban growth and temperature response patterns in Freetown and Bo town in Sierra Leone. Multispectral Landsat images obtained in 1998, 2000, 2007, and 2015 are used to quantify growth and land surface temperature responses. The contribution index (CI) is used to explain how changes per land use and land cover class (LULC) contributed to average city surface temperatures. The population size of Freetown was about eight times greater than in Bo town. Landsat data mapped urban growth patterns with a high accuracy (Overall Accuracy > 80%) for both cities. Significant changes in LULC were noted in Freetown, characterized by a 114 km2 decrease in agriculture area, 23 km2 increase in dense vegetation, and 77 km2 increase in built-up area. Between 1998 and 2015, built-up area increased by 16 km2, while dense vegetation area decreased by 14 km2 in Bo town. Average surface temperature increased from 23.7 to 25.5 ◦C in Freetown and from 24.9 to 28.2 ◦C in Bo town during the same period. Despite the larger population size and greater built-up extent, as well as expansion rate, Freetown was 2 ◦C cooler than Bo town in all periods. The low temperatures are attributed to proximity to sea and the very large proportion of vegetation surrounding the city. Even close to the sea and abundant vegetation, the built-up area had an elevated temperature compared to the surroundings. The findings are important for formulating heat mitigation strategies for both inland and coastal cities in developing countries.


Introduction
There has been an increase in the number of urban dwellers, together with an accompanying expansion of built-up area globally [1].Urban areas are strategic areas economically, as well as from an administrative perspective.They are important for issues such as the improvement of education and health delivery of a nation.Despite their socio-economic importance, urban areas and characteristic complex land use and land cover (LULC) spatial structure also pose a variety of environmental changes [2][3][4][5].According to Acharya et al. [5], the benefits of urban growth in developing countries include opportunities for employment, specialization, and the better production of goods and services.The challenges, however, include air pollution and water pollution in industrialized areas, while flash flooding is prevalent in highly impervious areas.Another notable challenge of urban development is temperature elevation, especially in densely built-up areas [6,7].Studies have shown that urban areas are comparatively warmer than undisturbed surroundings such as rural areas; a phenomenon called Urban Heat Island (UHI) [8][9][10][11][12].According to Gusso et al. [8], cities use construction materials such as concrete and asphalt, which do not allow water to penetrate and absorb a large amount of heat, thereby increasing urban temperatures.Elevated temperature results in increased outdoor and indoor human thermal discomfort, as well as increased heat-related health risk [13][14][15][16].Urban heat islands have maximized the number of heat wave days and tropical-like night conditions in several main cities, including Paris, Baltimore, Washington D.C., and Shanghai, during the summer [17][18][19].Furthermore, the Intergovernmental Panel on Climate Change (IPCC) [20] stressed that land cover changes have the potential to raise air temperatures of urbanized areas by 4 • C by 2100.The changes and associated adverse impacts seriously threaten the sustainable development of urban areas [21].Urban land use and land cover heterogeneity, as well as changes, result in the complex and varied spatial structure of heat intensities which also vary from city to city.It is thus important to establish city specific land surface temperature patterns in order to derive relevant mitigation and response strategies.
Remote sensing offers a variety of options for monitoring both LULC and LST spatial structure.Unfortunately, space-borne sensors detect thermal infra-red at either a low (e.g., above 500 m such as METEOSAT) or medium (e.g., 30-500 m such Landsat, ASTER and MODIS), but not high, spatial resolution (e.g., below 30 m such as SPOT).This results in mismatch in the resolution between retrieved LULC and LST maps.High resolution thermal data is often obtained from air-borne missions.Generally, high spatial resolution datasets are expensive to gather, have a low temporal resolution, usually lack a thermal infra-red component, and have very limited historical archives not sufficient for long term analysis [22].Medium resolution multi-spectral datasets are often reliable for urban LULC and LST analysis.For example, Landsat has large stores of visible, infra-red, and thermal data archives spanning from as early as 1972 to present [6,22].Recently, studies showed that Landsat data are effective and very accurate in mapping urban LULC distribution, as well as changes thereof [22][23][24][25].For example, using Landsat data, Mushore et al. [9] retrieved LULC spatial and temporal patterns in Harare between 1984 and 2015 at overall accuracies greater than 80%.Studies have also proved the effectiveness of Landsat thermal data in mapping land surface temperature variations, including those in complex urban settings [26][27][28].Recently, multi-temporal Landsat data was used to develop a model to predict future urban surface temperatures in Harare [29].Mushore et al. [29], showed that if historical growth patterns will persist, land surface temperatures will increase by as much as 5 • C by 2045.Therefore, the utility of medium resolution datasets in quantifying the impact of urban growth on LST patterns needs to be continually exploited.This is necessary in cities of low Gross Domestic Product countries such as in Africa, especially where similar studies have not yet been done; for example, in Sierra Leone.
In Africa, the studies have been confined to a few cities mainly in South Africa, Zimbabwe, and Nigeria.For example, Odindi et al. [7] investigated the impact of seasonality of urban greenery on heat island patterns in the Ethkwini municipality in South Africa.However, although they used 30 m multispectral Landsat 7 data for LULC mapping, surface temperatures were retrieved from course resolution (1 km) MODIS thermal data.Other studies in Africa were also confined to a single city; for example, Mushore et al. [9] only focused on Harare in Zimbabwe, while in West Africa, Abegunde and Adedeji [30] focused on Ibadan in Nigeria.Given the projected urban growth which must be faster in developing countries, there is thus a need to understand the implications in other parts of Africa [19].While Odindi et al. [31] compared LST patterns in coastal cities of South Africa, there is a general paucity of literature on comparing LST patterns between two cities of an African country.Precisely, there is a lack of literature comparing LST patterns of two cities, especially with one being inland and the other being coastal, such as Freetown and Bo town in Sierra Leone.As such, there is the need for a novel study to understand urban growth patterns, as well as responses of LST, in Sierra Leone, in West Africa.Such analysis is important for understanding both the differential effect of urban growth and of global warming between a coastal and an inland city in West Africa.Adaptation and mitigation strategies derived from such an analysis will take into account the position of a city relative to the ocean.Furthermore, the Contribution Index (CI) has not yet been used to compare growth patterns of two cities, as well as to explain the impacts of growth on surface temperatures in West Africa.To the best of our knowledge, the index has only been successfully tested on the African continent in South Africa [7,23,31] and in Zimbabwe [9].Odindi et al. [31] used CI to compare LULC and LST patterns between coastal cities of South Africa, but did not compare a coastal city with an inland city.Although Odindi et al. [31] compared LST variations in two cities; they used course resolution MODIS data, leaving a gap on comparison analysis using Landsat data in Africa.Liu and Weng [32] also found the 30 m visible and infrared, as well as the 90 to 120 m resolution thermal infra-red, Landsat data to be optimal in the analysis of the relationship between LULC and LST patterns.
The objectives of this study are thus to (1) use remote sensing to determine urban growth patterns in Sierra Leone; (2) quantify the effect of urban growth on spatial and temporal LST patterns in two major cities of Sierra Leone using the CI; and (3) understand the differences in responses of LST to urban growth and global warming between a coastal city (Freetown) and an inland city (Bo town) in Sierra Leone.The study hypothesizes that urban growth patterns should differ between Freetown and Bo town and thus influence LST spatial and temporal changes to differ between the two cities.

Study Area
The study was conducted in the two major cities of Sierra Lone; Freetown and Bo town (Figure 1).Freetown is the major port city on the Atlantic Ocean and is located in the western area of Sierra Leone.Bo town is the second largest city in Sierra Leone (after Freetown) and the biggest city in the Southern Province.Bo town serves as the capital and administrative focus of Bo District in the Southern Province.Freetown has a total area of 357 km 2 and a population of 772,873, constituting 15.53% of the total Sierra Leonean population [33].From the projected population of local administrative data from 2005 to 2014 (http://statistics.sl),out of 6,348,350 populations in Sierra Leone, 27.14% lives in urban areas, with 16.4% living in Freetown and 4% living in Bo [33].In Sierra Leone, the national census should be done once every 10 years.To date, five censuses have been conducted in 1963, 1974, 1985, 2004, and 2015.Another census was supposed to be done between 1994 and 1995, but was postponed due to the civil war which commenced in 1991 in the country.In order to ascertain that the population was growing in the study area, we used all the available data from the five censuses, although focus was on the period between 1998 and 2015.Census statistics are obtainable from Statistics Sierra Leone at national, town, and chiefdom levels.Therefore, population statistics for Freetown and Bo town were obtained at the town level.
Climate summaries were obtained from the Sierra Leone Meteorological Department under the Ministry of Transport and Aviation (http://mta.sl/meteorological-department).Freetown and Bo town experience a tropical climate, with a rainy season from May to October and a hot dry season from November to April.Freetown has an average annual precipitation of more than 3500 mm.It receives the highest amount of rainfall in the country due to its proximity to the Peninsula Mountains and Atlantic Ocean.The average annual precipitation of Bo town is around 2616.6 mm.The annual average minimum temperature for Freetown is around 23.8 • C, while the average maximum temperature  The topography of Freetown is undulated.Elevation ranges between 100 m and 700 m, with slopes exceeding 50 m and Bo town is 104 m above sea level.The prevailing winds are the south west monsoon during the wet season and the northeastern harmattan, which is a dust laden wind from the Sahara Desert, during the dry season.In this research, we have taken the most developed parts of Bo town (10,808.57ha or 108.08 km 2 ) and the most developed parts of Freetown (51,896.79ha or 518.97 km 2 ), as illustrated in Figure 1.The red bordered area of Freetown and Bo town (study areas) are the rapidly developing areas.The period between November and April was selected for temperature analysis because it is hot and dry, hence posing a threat to human thermal comfort., administrative spatial datasets from the National Tourist Board and Environmental Protection Agency, and ancillary secondary maps were also used as ground truth data for accuracy assessment.

Image Preprocessing
The remote sensing images are re-projected to the UTM WGS 84 N (UTM zone 29-North) following a third order polynomial fit and nearest neighbor resampling techniques.Digital numbers (DN) of TM5, ETM+, and OLI images are stored as 8 bit and 16 bit, respectively [34][35][36].These DNs of each image are converted to the top of atmospheric (TOA) spectral radiance using sensor specific calibration parameters directly obtained from the image MTL (metadata) file following the standard spectral radiance (Equation ( 1)).
where, ρ is the pixel surface reflectance, ρ e is an average surface reflectance for the pixel and a surrounding region, S is the spherical albedo of the atmosphere, L a is the radiance back scattered by the atmosphere, A and B are coefficients that depend on atmospheric and geometric conditions but not on the surface, and L is the spectral radiance.
The radiance of the reflective bands is then converted to a band interleaved by line (BIL) format to make them efficient for the atmospheric correction process in order to reduce atmospheric effects like water content, dust particles, aerosols, cloud, and varying sun angles, etc., which could significantly influence optical images and thereby degrade their spectral information.Hence, these are subjected to an atmospheric correction process to be applied to minimize those effects and produce corrected surface reflectance.The Fast Line-of-sight Atmospheric Analysis of Hypercube (FLAASH) is applied for the atmospheric correction process [37].FLAASH is a first principle of atmospheric correction tool which generally corrects wavelengths of visible, near-infrared, and shortwave infrared data.It uses the MODTRAN radiation transfer code [38] for retrieving atmospheric noises like aerosols, dusts, and water vapor content, etc., from dark land pixels in the scene based on a nearly fixed ratio between reflectance from pixels at 660 nm and 2100 nm [39].The overall FLAASH method takes input from the radiance and provides an atmospherically corrected surface reflectance image output using Equation (2).

Urban Growth Assessment Using Remote Sensing and Census Data in Freetown and Bo Town
Land use and land cover (LULC) maps for 1998, 2000, 2007, and 2015 were obtained using supervised image classification of multispectral Landsat data described in Section 2.2 above.Supervised image classification involves the use of ground control points obtained from field surveys or high resolution imagery to assist remote sensing software to assign LULC classes to pixels based on multi-spectral images.In each classification procedure, thermal data were left out since the objective was then to link LULC dynamics with LST derived from these data.The Support Vector Machine (SVM) algorithm was used because it was found to perform better than other common classifiers such as ANN, maximum likelihood, and Mahalanobis distance [22,40,41].SVM also comparatively requires very little training data.In each year, the ground truth LULC data collected from field work and auxiliary data were split into 70% (for classification) and 30% (for accuracy assessment) following the recommendation of Adelabu et al. [40].The area is classified into built-up, dense vegetation, sparse vegetation, water/wetlands, and agriculture land.A post classification change detection approach was used to determine the effect of growth on the spatial distribution and areal coverage of LULC types.According to Yu et al. [42], post classification is the most widely used change detection method.Due to simplicity and ease of interpretation, in this study, we detect changes in area per class, as was done by Salvati and Sabbi [43].
In order to link remotely sensed spatial and temporal patterns in LULC with population growth, census data for 1963,1974,1985,2004, and 2015 were used.Although the study focuses on the time interval from 1998 to 2015, the analysis of population dynamics includes time as far back as 1963 in order to take advantage of data availability, as well as to obtain a clearly convincing description of the population trends in the area.

LST Retrieval from Thermal Infrared Data
The steps as summarised by Weng et al. [26] and described in detail by Weng et al. [44] are followed to retrieve the land surface temperature from Landsat's thermal infrared data.The procedure involved (i) conversion of digital numbers (DN) to spectral radiance; (ii) computation of satellite brightness temperature from spectral radiance; and (iii) retrieval of land surface temperature from brightness temperature (emissivity correction).Full details of the steps are described in the Sections 2.5.1 and 2.5.2 below.

Conversion from Digital Numbers to Brightness Temperature
The DNs of the TIR bands of each year's ETM+ and TM5 images are converted to spectral radiance using the formula adopted by Chander and Markham [45] (Equation ( 3)) and Landsat 8's thermal infrared images were converted using the USGS standard (Equation ( 4)).
In the above equations, L λ is the spectral radiance in W/(m 2 srµm) received by the sensor from each pixel of the image.M L and A L are band specific multiplicative and additive rescaling factors obtained from the image MTL file, Qcal is the DN of each image, and QCAL max is the maximum DN (65535 for the 16-bit Landsat 8 and 255 for other Landsat missions).L max and L min are the maximum and minimum top of atmospheric (TOA) radiances in W/(m 2 srµm), respectively.
After the conversion of the DNs to the spectral radiance, the radiant images are converted to the blackbody temperature using (Equation ( 5)).
where T b is the effective at-sensor brightness temperature in Kelvin unit, L λ is the spectral radiance in W/(m 2 srµm), and K 1 and K 2 are prelaunch calibration constants in Kelvin unit obtained from the image MTL file.

Surface Emissivity (ε) Retrieval
The land surface emissivity is retrieved using the Normalized Difference Vegetation Index (NDVI) threshold method [45,46].According to the method, when NDVI < 0.2, the pixels are considered as bare lands and the emissivity is retrieved from the red spectral region.When NDVI > 0.5, the pixels are considered as fully vegetation coverage and the emissivity value is assumed to be 0.99.When NDVI ranges between 0.2 and 0.5, the pixels are considered as a mixture use of soil and vegetation.In this case, emissivity is retrieved using Equation ( 6), as follows: where ε v is the emissivity of vegetation coverage, ε s is the emissivity of soil surface, and, P v is the proportion of vegetation calculated from Equation ( 7), where NDV I s is the NDVI value of pure soil and NDV I v is the NDVI value of pure vegetation extracted from the NDVI image.
In Equation ( 6), the term ∆ε is the indication of the geometrical distribution of the natural surface, as well as the internal reflection whose value is considered as negligible for the plain and homogenous surfaces.However, in the case of a rough and heterogeneous surface, the value is assumed to be 2% Sobrino et al. [46] and is expressed by the following (Equation ( 8)): where F is the shape factor whose mean value for different geometrical distributions is assumed to be 0.55 [45,46].By summarizing Equations ( 6) and ( 8), the final equation for emissivity estimation is obtained by Equation ( 9), as follows: where m and n coefficients are calculated as: Brightness temperatures assume that the earth is a blackbody, which it is not, and this can result in errors in surface temperature.In order to minimize these errors, emissivity correction is necessary and this is done to finally obtain the land surface temperature (LST) from T b using Equation (11) [44].

Linking Urban Growth to LST
The effect of LULC in the warming or cooling of an area depends on the LULC type and the proportion of the total area occupied by each type.For example, vegetation cover and water/wetlands have a surface cooling effect due to latent heat transfer.However, even though they have a cooling effect, the overall value depends on the proportion of the total area they occupy [49].The warming or cooling extent of an LULC type taking into account the proportion of the total area it occupies is quantified using the Contribution Index (CI).The CI is used to link spatial structure, as well as long term changes in LULC, to LST intensities.The CI for each LULC type is computed for both cities using Equation ( 12) for all the periods mentioned in Section 2.2 [7,31,48].
D t is the difference between the average temperature of the entire study area and the average of the LULC class type.Variable S is the proportional area of the LULC type, which is the ratio of the area covered by the class to the total area of the study area.Positive values of CI indicate how much the LULC type contributes to raising the surface temperatures of an area, while negative values indicate a heat mitigation value.

Remote Sensing Based Urban Growth Assessment in Freetown and Bo Town
Visual inspection of Figure 2A-D indicates the expansion of built-up area in Freetown.This is notable in the northern, eastern, and western parts of the city.Since 1998, the city has been characterized by a tongue of dense vegetation occupying most of the central part of the city.This dense vegetation area is not diminishing, even as built-up area is expanding.The growth of Freetown concentrated along the northwestern and eastern margins is influenced by the ocean (Figure 2A-D).On the other hand, the growth of Bo town since 1998 has been largely characterized by expansion from the central to the southwestern areas of the city (Figure 2E-H).The growth of Bo town also infiltrated into densely vegetated areas between 1998 and 2015.Between 1985 and 2015, the agriculture area has decreased by about 84 km 2 , while the built-up area increased by almost 80 km 2 in Freetown (Figure 3).Dense vegetation areas increased by 22 km 2 , while sparse vegetation areas increased by 40 km 2 .The increase in vegetation areas could be part of an explanation of why bare areas reduced in area by 28 km 2 .A difference was observed in Freetown because, here, growth occurs along the coast away from the central zone of dense vegetation.Land use and land cover changes in Bo town were not as marked as in Freetown.For example, built-up areas increased by 15 km 2 , while areas with sparse vegetation increased by 7 km 2 in Bo town.During the same period, the dense vegetation and agriculture areas decreased by 14 km 2 and 10 km 2 , respectively.

Census Based Urban Growth Patterns in Freetown and Bo town
The population increased by almost ten times in both Freetown and Bo town between 1963 and 2015 (Table 2).Population densities also changed from 246.5 to 2023.8 (people/km 2 ) and from 246.2 to 1609.0 (people/km 2 ) in Freetown and Bo town, respectively.The population size of Freetown has always far exceeded that of Bo town, such that in 2015, the sizes were 1,050,301 and 173,905, respectively.

Responses of LST to Growth Patterns in Freetown and Bo Town
High surface temperatures (above 30 • C) are most notable in the northern and western parts of Freetown in 1998 (Figure 4A).Over the years, high surface temperatures have also been spreading southward along the western margin of the city (Figure 4B-D).Low surface temperatures below 22 • C have remained characteristic of the central and southwestern parts of the city.On the contrary, since 1998, the high surface temperature has spread from the central parts of the city of Bo town, especially towards the southwest (Figure 4E-H).Low temperature areas (below 20 • C) surround this expanding hot spot and are shrinking in size.The shape of high surface temperature areas in both Freetown and Bo town closely mimics that of the built-up area, indicating their strong warming influence.Conversely, low surface temperature patterns also track areas with vegetation cover, being low in dense vegetation areas in both cities.In both cities, average temperatures are rising with time (Table 3).

Link between Long Term Changes in LULC and LST Dynamics
The agriculture area has a positive contribution index (CI) in Freetown, indicating that such places increase heat in the city during the dry season (Table 3).Although the area under agriculture has reduced between 1998 and 2015, the CI has remained positive and increased.The cooling contribution of dense vegetation is increased as indicated by a CI of −0.47 in 1998 followed by −0.85 in 2015.Sparse vegetation also has a significant cooling effect in Freetown, although its value has decreased slightly between 1998 (CI = −0.34)and 2015 (CI = −0.24).The heat mitigation value of vegetation was also noted in Texas, where woodlands were 1.5-3.9• C cooler than neighboring areas.The cooling effect of dense vegetation was more than that of sparse vegetation indicated for an example by a CI of −0.67 compared to −0.13 in Freetown in 2015 for dense and sparse greenery, respectively.The built-up area in Freetown increased in terms of the warming effect by almost three times, as indicated by the CI of 0.26 in 1998 and 0.72 in 2015.
In Bo town, the surface cooling effect of dense vegetation is increasing significantly (CI = −0.55 in 2000 and −0.85 in 2015).At the same time, bare areas and areas with sparse vegetation are decreasing in terms of their warming and cooling effect, respectively (Table 4).The hot spot area expanding from the centre of the city, mainly to the southeast, can thus be explained by the increasing warming effect of the built-up area between 1998 (CI = 0.26) and 2015 (CI = 0.41).Although water bodies have a cooling effect, their contribution has remained minimal over the years due to the low proportion of the cities they occupy.In both cities, the CI for water has remained less than −0.2 in all the years.Urban growth patterns in Freetown are unique, in that they are characterized by the expansion of built-up and dense vegetation areas.Although Freetown is larger in size and growing faster, it was about 2 • C cooler than Bo town in all periods.

Discussion
The study obtained a high classification accuracy both in a coastal city (Freetown) and an inland city (Bo town).The overall classification accuracy reached the 85% recommendation by Anderson [49], because even at a 30 m resolution of Landsat optical data, the mixed pixel problem did not significantly affect the quality of the LULC maps produced.Despite the complexity of classification in urban areas due to surface heterogeneity, the mapping accuracies are also higher than the 80% overall accuracy recommended by Omran [50].The high level of accuracy can be justified by Voogt and Oke [51], who noticed that improvements that have occurred in satellite sensors over the years provide detailed and accurate land surface representation at a low cost.The high classification accuracy could also be attributed to the renowned performance of the Support Vector Machine algorithm [22,40,41].According to Jia et al. [41], the Support Vector Machine (SVM) algorithm was found to outperform other common classifiers such as ANN, maximum likelihood, and Mahalanobis distance.The algorithm was also used for multi-temporal Landsat-based classification in an urban setting in Harare, where overall accuracies above 80% were also obtained.These findings show the value of freely available medium resolution space-borne remotely sensed datasets for monitoring urban extent and growth, especially in resource-constrained nations.
The population increased by almost ten-fold in both Freetown and Bo town between 1963 and 2015, while the population densities also increased.In all the periods considered, the population size of Freetown has always far exceeded that of Bo town.Most of the economic and administrative activities of Sierra Leone are concentrated in Freetown, hence the larger population size and faster growth than Bo town.Furthermore, the beauty of the sea seems to make residents prefer to concentrate along the coastal margins of Freetown than to spread further inland towards the dense vegetation area.Besides increasing population sizes, built-up areas are also expanding in both cities.Growth patterns observed in both cities agree with earlier observations and predictions that urban population is growing, globally [19,52].Expansion of the built-up area in Freetown has been mainly concentrated along the coast and is most notable in the northern, eastern, and western parts of the city.This growth along the northern margins of Freetown explains why the dense vegetation area in the central part of the city is not diminishing even as the built-up area is expanding.A different pattern is observed in Bo town, where the built-up area is expanding from central locations outwards.Unlike in Freetown, the growth of Bo town has led to a reduction in the area of the densely vegetated LULC category between 1985 and 2015.As observed in Bo town, in most studies, the proportion of total area occupied by dense vegetation decreases with continuous urban expansion [26,27,53,54].Kamusoko et al. [54] observed that the expansion of built-up areas in Harare Zimbabwe pushed most dense vegetation locations outwards to the peripheries of the city.
As expected, temperature responded strongly to spatiotemporal dynamics of LULC in both Freetown and Bo town.High temperatures in both cities were observed in built-up areas and their extent increased with time as the cities were expanding.The influence of buildings explains why high surface temperatures (above 30 • C) were recorded in northern and eastern parts of Freetown.Over the years, surface temperatures in this regime have also been spreading southward along the western margin of the city following the expansion of the built-up area.The shape of high surface temperature areas in both Freetown and Bo town closely mimics that of the built-up area, indicating their strong warming influence.This concurs with Sha and Ghauri [28], who observed that surface urban heat island expands with expansion in a built-up area.Buildings reduce heat removal by advection and reduce the sky view factor, thus limiting heat escape to space, while walls and pavements absorb and emit heat [28,53,55,56].This results in large amounts of stagnant heat and high temperatures, especially in closely packed and high rise buildings.The warming in both cities could also be explained by increased anthropogenic activities supported by an increasing population size in both cities over time, which increases long wave radiation in the lower atmosphere.Nayak and Mandal [3] and Grimmond [57] also attributed urban warming to both LULC changes and other anthropogenic effects such as greenhouse gas emissions.The rising temperature in response to the growth of both cities can be captured by the explanation that, as population grows, urbanization increases and the magnitude of the urban heat island also expands [58].Similar findings were obtained in Australia between 1951 and 2003, where land cover changes produced statistically significant warming [59].
Vegetation cover has been indicated to be a strong mitigation measure against the elevation of surface temperatures in both cities.For example, in Freetown, low surface temperatures (below 22 • C) remained characteristic of the central and southwestern parts of the city where buildings have not yet replaced vegetation cover.Similarly, low temperature areas (below 20 • C) surround an expanding hot spot in the central parts of the city of Bo town.The heat mitigation value of vegetation was also captured by a strong negative Contribution Index (between −0.5 and −1) in areas with dense and sparse vegetation.This concurs with Odindi et al. [7] who in the EThekwini municipality, South Africa, showed that the temperature reduction effect of vegetation increases with the percentage of total area covered.Although water bodies also have a cooling effect (negative Contribution Index [CI]), their contribution has remained minimal over the years due to the low proportion of the cities they occupy in both cities.Based on CI, the cooling effect of dense vegetation was more than of sparse vegetation, which echoes the suggestion by Zhang et al. [60] that not only vegetation types but also spatial structure affects LST distribution.Vegetation cover promotes surface cooling due to latent heat transfer.
In both Freetown and Bo town, agriculture areas were causing warming of the city, as indicated by a positive Contribution Index (CI) in all periods.This could be because, during the dry seasons, agriculture areas will either be covered by drying crop residue or will be semi-bare/bare, thus absorbing a considerable amount of heat.This is in agreement with the findings of Mushore et al. [61] in Harare, which showed that, during the hot dry season, croplands act as a heat source as they absorb and release large amounts of heat due to negligible evaporation.Although areas under agriculture have reduced between 1998 and 2015, the CI has remained positive and increased, implying an increased warming contribution to the city.This could be because the temperature of these areas has increased over the years with the changes attributed to global warming.Early planting of crops means that by the dry season the residues will be completely dry, resulting in high heat absorption, which could also be another explanation.However, the decrease in area under agriculture may indicate a shift of agriculture to the secondary industry and services in both cities.In other cities such as Harare [9], growth is also characterized by the major replacement of dense vegetation and agriculture areas with building and impervious surfaces, resulting in warming.Therefore, the surface warming mostly of Freetown between 1998 and 2015 can be attributed to global warming, the warming effect of dry agricultural land, and increase in the built-up area which absorbs a significant amount of heat.This agrees with Jiang and Tian [62], who demonstrated that the construction of buildings leads to the transition of an area from a dense vegetation low temperature to sparse vegetation high temperature zone.
Even in coastal cities where the water table is presumed to be high and sea breezes cool the atmosphere, a high density of buildings still causes warming.Although Freetown is larger in population size as well as built-up extent and also growing faster, it was cooler than Bo town in all periods (by about 2 • C).The difference could be a result of surface moisture and cold air advection due to proximity to the sea.Surface wetness reduces the temperature of a surface due to increased evaporation and latent heat transfer [56].According to Rasul et al. [56], green areas and water bodies act as urban cool islands, hence the low temperature of Freetown despite being larger in size than Bo town.Besides being close to the sea, the proportion of dense vegetation cover is greater in Freetown than Bo town, which reduces the average temperature of the city.According to Sithole and Odindi, green spaces act as heat sinks, tend to be porous, and assimilate heat.Due to the influence of the sea, dense buildings and high surface temperature are found along the coast in Freetown.This has also led to the sustenance and expansion of a tongue of dense green area and low temperature in the central part.Water and vegetation which surround the built-up area of Freetown act as a sink to these gases, which may also explain the lower temperature there than in Bo town.According to Odindi et al. [7], the heat contribution of dense vegetation is similar to that of water, hence Freetown is surrounded by cool areas resulting a in lower mean surface temperature than Bo town.

Conclusions
We have compared urban growth and land surface temperature patterns between a coastal city (Freetown) and an inland city (Bo town) in Sierra Leone in this paper.Multi spectral Landsat data are used to quantify land use and land cover, as well as surface temperature, changes between 1998 and 2015.Based on the findings of the study, we conclude that multi-spectral Landsat data and the Support Vector Machine algorithm retrieve LULC spatial patterns and urban growth with a high accuracy.The growth patterns of Freetown are concentrated along city margins at the coast, while Bo town expanded from the center outwards.The abundance of dense vegetation and proximity to ocean makes Freetown cooler, although it is larger in population and is expanding in terms of the built-up area faster than Bo town.However, even in cool areas such as at the coast, built-up areas have warmer surface temperatures than non-built-up areas such as dense vegetation areas.Expansion of the built-up area from the city core pushes out vegetation towards the margin, resulting in a high temperature towards the center, as in Bo town.Overall, the built-up area expansion increases urban temperature, in addition to the effect of global warming, while vegetation has a strong heat mitigation effect.The Freetown-Bo town scenario has indicated that it is possible for a small city to be warmer than larger and faster growing cities within the same country.Temperature patterns depend heavily on position relative to ocean, as well as the size and spatial structure of dense vegetation area.Therefore, even vegetation and water patches around a built-up area (not only those within) have an influence on its temperature.Although the study managed to convincingly link urban growth induced LULC changes with LST dynamics, future efforts must address some limitations which could hamper the reliability of the findings.The study depended on medium spatial resolution Landsat datasets, whose temporal resolution of 16 days is low.This, together with the cloud free image requirement for surface analysis, resulted in a limited amount of data available for the study.
In the presence of sufficient data, averages could have been computed to eliminate the effects of randomness associated with the use of single date images to represent an entire month.Due to the low temporal resolution of Landsat data, it is difficult to obtain in-situ meteorological data at the exact time of satellite overpass for a comparison of temperatures obtained from remote sensing with in-situ observations of air temperature in Sierra Leone.Meteorological operations in Sierra Leone are still manned; taking observations at World Meteorological Organization (WMO) prescribed synoptic hours which do not coincide with the overpass times of Landsat missions.Limited access to in-situ meteorological data inhibited the analysis to test the validity of the findings of this study, although they agreed with global trends.Reflective bands of Landsat are at a higher spatial resolution than the thermal dataset (for example 30 m versus 100 m for Landsat 8).This mismatch has the potential to increase the mixed pixel problem on LST retrievals, thus compromising the link between LULC (30 m resolution) and LST (100 m), even though thermal data is downloaded at a resolution of 30 m after resampling.Other factors which affect thermal properties such as differences in building material and roof types between Freetown and Bo town were not investigated in this study.

is 29 . 9 •
C. The annual mean minimum temperature for Bo town is 21.2 • C and the average maximum temperature is 31.3• C.

Figure 1 .
Figure 1.Location of Freetown and Bo town in Sierra Leone, West Africa.
This study uses cloud free and geometrically corrected Landsat imagery from the Earth Resources Observation and Science (EROS) center through the United States Geological Survey (USGS) Global Visualization Viewer.The path/row was 202/54 for Freetown and 201/054 for Bo town.The image scenes dated to 27 February 1998 (TM5), 3 February 2000 (ETM), and 23 February 2007 and 28 January 2015 (OLI) for Bo town, and 28 February 1998 (TM5), 3 February 2000 (ETM), 27 February 2007 (ETM), and 4 February 2015 for Freetown.Apart from satellite imagery, several referenced datasets like ground GPS data of different LULC categories, Population census data from the Statistics Sierra Leone, and mean temperature and mean humidity from the Sierra Leone metrological Department were used.High resolution contemporary satellite imagery (GEOEYE-1 and Google Earth historical image of 2015)

Figure 2 .
Figure 2. Urban growth induced LULC changes in Freetown (A-D) and Bo town (E-H) between 1998 and 2015.

Table 1
indicates the overall accuracy (OA) and kappa coefficient (k) obtained in LULC classification for different years in Freetown and Bo town.The overall accuracies were greater than 85% for both cities in all years.Accuracies per individual LULC class (i.e., user accuracy (UA) and producer accuracy (PA)) are shown in Appendix A.

Figure 3 .
Figure 3. Urban growth induced LULC changes in Freetown and Bo town (1998 to 2015).

Figure 4 .
Figure 4. Land surface temperature change in Freetown (A-D) and Bo town (E-H) between 1998 and 2015.

Table 1 .
Accuracy of multi-temporal LULC classifications in Freetown and Bo town.

Table 2 .
Census-based population growth in Freetown and Bo town.

Table 3 .
Changes in the heat source/sink role of land use and land cover types in Freetown between 1998 and 2015.Green means vegetation.

Table 4 .
Changes in the heat source/sink role of land use and land cover types in Bo town between 1998 and 2015.