Predicting of Land Surface Temperature Distribution in Freetown City, Sierra Leone by Using Polynomial Curve Fitting Model

Global warming has attracted much concern about the worldwide organiza-tion, civil society groups, researchers, and so forth because the worldwide surface temperature has been expanding. This investigation intends to assess and compare the ability of a combination of land cover indices to predict the future distribution of land surface temperatures in Freetown using the Polynomial model analysis. Landsat satellite images of 1988, 1998, 2000, 2010, and 2018 of the Freetown Metropolitan zone were utilized for analysis. The investigation had adopted two land covers indices, Modification of normalized difference water index and Urban Index (UI) (e.g., MNDWI and UI) and applied a multi regression equation for forecasting the future LST. The stimula-tion results propose that the development will be accompanied by surface temperature increases, especially in Freetown’s western urban area. The temperature prevailing in the west of the metropolitan area may increase in the city somewhere in the range from 1988 to 2018. Additionally, the res ults of the LST prediction show that the model is perfect. Our discoveries can be represented as a helpful device for policymakers and community awareness by giving a scientific basis for sustainable urban planning and management.


Introduction
The urban heat island (UHI) impact shows a higher air and land surface tem-Journal of Geographic Information System perature (LST) in urban regions in contrast with the encompassing rural area, generated by high levels of near-surface energy emission, solar radiation absorption of ground objects, and low rates of evapotranspiration [1]. This is as a result of urbanization and overpopulation. These have resulted in most important human activities, creating enormous impacts on the ecosystem at the local, regional, and global scales. It includes converting natural land surfaces into anthropogenic impervious surfaces due to the introduction of low and high albedo materials, which may lead to dramatic climate changes; consequently, an urban climate is warming systematically since the last decades. Such heat rise and warming can be embraced by numerous biophysical well-being intricacies like heat stress, air contamination, and other well-being-related issues [2]. In recent years, studies have proved that the developed area and bare land accelerate UHI's impact, though green space and water diminish the UHI intensity [3] [4]. Therefore, it is essential to study the prediction of land surface distribution [5]. This is a very vital in urban areas of developing countries that are experiencing an increase in land surface temperature, which may result in adverse climate and related factors. Studies have been confined to few cities in Africa, mainly in Ikom city in Nigeria, Harare metropolitan city in Zimbabwe. For example, Muduako et al., 2016 reported that LST in an urban area is one factor associated with urban heat rise and microclimatic warming within a city [6]. Mushore et al., 2017 documented that urban growth will increase warming and result in future high temperatures unless mitigation efforts are strengthened [7]. To the best of our knowledge, little work has been implemented to predict land surface distribution in Freetown, Sierra Leone, and quantify the implication of land surface temperature. Freetown is a coastal city with a rapid growth in urbanization, and such factor is expected to affect the increase in surface temperature. Such a study is vital both given quantifying the effects of expansion on lower tropospheric temperature and in developing knowledge on growth implications on future climate and thermal comfort on urban residents. Space-borne remote sensing has the benefit of simultaneously observing over large areas enabling spatial analysis and direct comparison in space. Remote sensing satellites such as Landsat have many benefits, including large stores of archival data and continuous improvements in data quality over time. For example, Landsat has archival data dating as far back as 1972, making it an important data source for detecting land surface changes over long periods. This paper used space-borne multi-spectral Landsat imagery to predict the future land surface temperature in Freetown city in Sierra Leone from 2018 to 2026.

Study Area
Study area Freetown is the first capital of Sierra Leone (Figure 1) that lies between latitude 8˚05'N and 8˚30'N and extends between longitudes 12˚50'W and 13˚20'W. It is the major port city on the Atlantic Ocean. It is located in the  western area of Sierra Leone. Freetown experiences a tropical climate with a rainy season from May to October and a dry season from November to April.
The yearly average minimum temperature for Freetown is around 23.8˚C, while the average maximum temperature is 29.9˚C. annual mean minimum temperature [8]. The topography of Freetown is undulated. Elevation ranges between 100 m and 700 m, with slopes exceeding 50 m above sea level. The prevailing winds are the southwest monsoon during the wet season and the northeastern harmattan, a dust-laden wind from the Sahara Desert during the dry season. The red-bordered areas (study areas) are a rapidly developing area.

Data Used
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, map projection of the collected satellite images is Universal Transverse Mercator (UTM) within Zone 28 N-Datum World Geodetic System (WGS) 84 (Table 1). Various software packages were utilized because each one has strength in some operations needed for this study.

Image Preprocessing
Landsat satellite images (Landsat 5 TM and Landsat 8 OLI/TIRS) with metadata (MTL) file for the study area were utilized for assessing land surface temperature (LST). The satellite images obtained are processed by a set of pre-processing procedures. The pre-processing in clouded radiometric calibration, atmospheric Journal of Geographic Information System correction (dark-object subtraction) using the Envi 5.3 and geometrical distortions correction. After the pre-processing procedure, the satellite images are again re-sampled with pixel sizes of 30 × 30 m for all bands which also includes the thermal band. As referenced before 30 m goals resolution reflective bands of Land-sat images obtained from various seasons were piled up into a single multi-band record for the categorization.

LST Recoupment from the Information in Thermal Infrared Data
In the proposed research, a mean thermal infrared image containing digital numbers was analyzed every year through computing 2 images collected in each month using LST based procedures that were proposed by [9]. The procedure suggested involves 1) transformation of digital number data contained inside Landsat images to radiance, 2) transformation of radiance into temperature of blackbody and 3) blackbody temperature due to brightness into temperature of land surface (emissivity correction). Adopted from [10] emissivity mapping of land surface was calculated using NDVI.

Transformation of the Digital Number (DN) to Spectral Radiance
As shown in equation no. 1 which is considered by [11], data gathered from digital numbers are transformed to spectral radiance. Data from digital numbers are of the TIR bands of ETM+ and TM5 images for individual years. Then the use of Equation no. (2) thermal infrared images of Landsat 8were transformed considering the standard of USGS.

Calculation Brightness Temperature of Blackbody Using Spectral Radiance
With the use of Equation (3) the radiant images were transformed into temperature detected in blackbody [12].
In the above Equation brightness temperature in Kelvin unit is measured using sensor which is denoted with T b , spectral radiance in W/(m 2 srμm) is denoted with L λ and through image MTL file prelaunch calibration constants in Kelvin unit is obtained. These are represented as K 1 and K 2 . Error in measurement of surface temperature can arise during the execution process due to detection of earth as a blackbody because of its brightness [12]. To minimize these error values, emissivity correction is performed to attain land surface temperature (LST) from Tb using Equation (4), [13].

Retrieval of Surface Emissivity (ε)
A threshold method was suggested by [14] to achieve land surface emissivity from Normalized Difference Vegetation Index (NDVI) threshold method. If NDVI < 0.2 then pixels are acknowledged as barren land. The emissivity of these barren lands was derived from red spectral region. If NDVI > 0.5 then pixels are acknowledged as lands thoroughly covered with plants or agricultural fields and its emissivity value was estimated to be 0.99 [10]. When NDVI value lies between the range of 0.2 to 0.5 then pixels are acknowledged as partially covered with vegetation. Using Equation (4) the emissivity is obtained.
where v ε the emissivity of vegetation coverage, s ε is the soil surface emissivity and, v P is computed to obtain the proportion of vegetation using Equation (5).
where NDVI s is the NDVI value of pure soil and NDVI v is the NDVI values of pure vegetation derived from NDVI image. ε ∆ in Equation (7) indicates distribution of the land surfaces according to its geometry in addition to its the internal reflection. The value of internal reflection is very low for the ordinary and uniform surfaces [10]. The value of internal reflection is however found to be 2% for rough and heterogeneous surface. ε ∆ is calculated using Equation (6) where shape factor is denoted by F. The mean value for different land surfaces according to its distributions based on its geometry is assumed as 0.5. Journal of Geographic Information System Summarizing Equation (5) and Equation (6), the final equation to calculate emissivity estimation is given as Equation (7). v mP n ε = + (7) where m and n co-efficient are calculated as below in Equation (8):

Brightness Temperature to LST
In the above Equation (9), wavelength of emitted radiance (11.5 μm) [11] is represented by λ, ρ = hc/σ [15] Table 2 involves urban indices and indices of vegetation. These were estimated through a digital number of measured bands. As recently referenced, the ultimate objective is to investigate the superiority in a relationship's characteristics with the surface temperature. In addition to this to perceive indices with highest grounded capacity for assessing temperature of surfaces of urban areas, a couple of indices were tried. Simulating LST using polynomial curve fitting was shown in     entire year [5]. Land cover indices and LST values were mined from each pixel in the designated area of study for individual point data type in order to achieve this objective. These points were utilized to calibrate the linear regression model [20].

The model offers a general idea about the correlation between LST and LULC in-
dices. This realization is relatively coherent with those narrated by other earlier researches [21].

Polynomial Curve Fitting for Land Use Land Cover Indices Prediction
The Land Use Land Cover Indices for 26. 12.1988, 9.3.1998 Curve fitting over m pairs of data ( ) ( ) ( ) x y x y x y  is a process of obtaining a polynomial regression between the two pairs as follows: ( )

Prediction Accuracy Assessment for LST
where: the modeled surface temperature and the actual land surface temperature recorded from Landsat data for the ith pixel is denoted by T predictedand T observed . Root Mean Square Error and ration RMSE/std can also be regarded as the measure for calculating accuracy of the prediction model in forecasting temperature. Land surface temperature distribution for the period from 2026 is then found using the prediction model after the accuracy assessment.

Monitoring and Assessing LST Changes
The   of the study area. When surface temperature increases, the amount of NDVI value decreases. This is the most common behavior found in grassland and sparse vegetation areas.

Retrieval of Surface Temperature from the Land Use Land Cover Indices
This technique of using temperature forecast to estimate LST pattern to model and estimate UHI is one such contribution of this study. As presented in Table 4 and Figure 4, For the Polynomial curve fitting analysis of land surface temperature distribution the Urban Index (UI) and Modification of normalized difference water index (MNDWI) were chosen as the predictors [5]. The regression model was tested on an independent Landsat information acquired in February 2010 and February 2018 and the model nearly resembled like the temperature patterns ( Figure 5, Figure 6). Based on the Information acquired from independent Landsat in 2018 Regression model was tested. Temperature redeemed from UI and MNDWI based on their direct relation to thermal infrared information (Band 10) of Landsat 8. Based on 161 points samples over the inspected region, the UI and MNDWI forecasted surface temperature with advanced correctness (mean relative rate percentage 5.88%, and 4.41% with root mean square error 1.61˚C, 1.31˚C and ration RMSE/std are 0.5 and 0.44 for 2010 and 2018 respectively.

Predicted Temperature Distribution in Freetown up to the Year 2026
The rising temperature patterns observed between 2000 and 2018 may proceed through 2026 (Table 5) and Figure 7. The scope of high-temperature classification (greater than 35˚C) was predicted to increase to the detriment of low-temperature classes. According to LST's values, the predicted values in 2026 were more than the simulated value in 2018. However, in the predictions for South-eastern regions where low-density residential areas are found were moderately cooler than Northern-eastern regions where high-density residential areas are found. Predictions demonstrate that land surface temperatures underneath 32˚C will possibly stay common in the northern half, where low and medium density residential is found. Besides, expecting that development patterns observed between 2000 and 2018 endure, the extension would high-density built-up areas would bring about high surface temperatures (over 35˚C) in eastern regions, Journal of Geographic Information System such as in high-density residential areas. It is worth noting that the LST expanded for all land cover types.

Conclusions and Outlook
The investigation aims to forecast future distribution of land surface temperatures in Freetown utilizing the Polynomial curve fitting analysis. The research used two land cover indices (e.g., MNDWI and UI) and implemented a multi-regression equation to predict future LST. We declare that the Urban Index (UI) and (MNDWI) predicted surface temperature with high accuracy (mean relative rate percentage 5.88%, 4.41%, and root mean square error 1.63˚C, 1.31˚C, and ration RMSE/std are 0.5 and 0.44 for 2010 and 2018 respectively). Urban expansion development will be accompanied by surface temperature increments, especially in Freetown's western metropolitan area. The temperature that dominates the west urban region may increase in the city between 2000 and 2018 may proceed through 2026. For example, different factors, such as effective mitigation procedures and changes in city development policies, can impact surface temperature patterns. By and large, this study's discoveries underscore the significance of medium resolution satellite information in foreseeing future surface temperatures in urban settings. However, there is a need for future studies to explore the feasibility of these methods and techniques at national or regional spatial levels.