Spatial Variation of Land Use / Cover Composition and Impact on Surface Urban Heat Island in a Tropical Sub-Saharan City of Accra, Ghana

: Rapid urbanization is one of the most crucial issues in the world of the 21st century. Notably, the urban heat island phenomenon is becoming more prominent in megacities and their hinterlands in temperate and subtropical climatic regions. In the daytime in summer, there exists a high possibility of accelerating the land surface temperature (LST) in desert cities, due to the alterations made by human beings in the natural environment. In this study, we investigate the spatial formation of LST in a tropical sub-Saharan city of Accra, a gateway to West Africa, using Landsat data in 2003 and 2017. Machine learning techniques and the di ﬀ erent spatial and statistical methods such as tasseled cap transformation (TCT), urban-rural gradient, and multiresolution grid-based and landscape metrics were employed to examine procured land use / cover (LUC) and LST maps. LUC was classiﬁed into ﬁve categories: Built up, Green 1, Green 2, Bare land, and Water. The results of the analysis indicate that Built up, Green 2, and Bare land had caused the highest heating e ﬀ ect while Green 1 and Water had caused the considerable cooling e ﬀ ect during the daytime in Accra. The urban-rural di ﬀ erence in LST recorded 1.4 ◦ C in 2003 and 0.28 ◦ C in 2017. The mean size, mean shape, largest patch, and aggregation of Built up, Green 1, and Green 2 had a strong relationship with the mean LST. It is essential for urban planners to carefully examine the formation and e ﬀ ect of the urban heat island (UHI) for sustainable urban development and landscape policy toward mitigation and adaptation planning in Accra.

The UHI phenomenon results from the relative rise of land surface temperature in the urban areas compared to the surrounding rural environs [11][12][13]. There are two different concepts, namely, atmospheric UHI and surface UHI. Canopy layer UHI and boundary layer UHI are measured to calculate atmospheric UHI, while land surface temperature (LST) is used to determine surface  [26]. Ghana boundary layer was downloaded from DIVA-GIS [27].

Satellite Data Used and Pre-Processing
Landsat is the most long-lasting satellite program launched by the National Aeronautics and Space Administration (NASA) and United States Geological Survey (USGS) with a history of more than 48 years, formerly named ERTS (Earth Resources Technology Satellite) [28]. Landsat satellites provide earth surface images for LUC detection. This data has been widely used by LST researchers [11,29,30]. Captured Landsat images in 2003 Landsat 7 ETM (12 February; 10: 04 GMT) and Landsat OLI/TRIS 8 (25 January; 10: 15 GMT) in 2017 were employed for our analysis [26]. The pre-processing  [26]. Ghana boundary layer was downloaded from DIVA-GIS [27].

Satellite Data Used and Pre-Processing
Landsat is the most long-lasting satellite program launched by the National Aeronautics and Space Administration (NASA) and United States Geological Survey (USGS) with a history of more than 48 years, formerly named ERTS (Earth Resources Technology Satellite) [28]. Landsat satellites provide earth surface images for LUC detection. This data has been widely used by LST researchers [11,29,30]. Captured Landsat images in 2003 Landsat 7 ETM (12 February; 10: 04 GMT) and Landsat OLI/TRIS 8 (25 January; 10: 15 GMT) in 2017 were employed for our analysis [26]. The pre-processing was conducted before the retrieval of LST. Atmospheric correction and radiometric calibration through the conversion of digital number values were taken from the satellite data: (1) multispectral bands' value of surface reflectance and (2) thermal bands' value of satellite brightness temperature. L λ was derived for Landsat 7 ETM and Landsat 8 OLI/TRIS from the quantized calibrated pixel value using Equations (1) and (2).
The brightness temperature was calculated using Equation (3). (The DN values of the thermal bands of Landsat ETM, Landsat-8 OLI/TIRS were transformed).

Retrieval of LST
The commonly used method of acquiring unprocessed Landsat data requires the DN value of thermal bands (Landsat thematic mapper (TM) and ETM = band 6, and Landsat 8 OLI/TIRS = bands 10 and 11) [11,31,34]. First, we derived absolute radiance values and conducted the derivation of satellite brightness temperature. We used the pre-processed bands (band 6 for Landsat 7 ETM and bands 10 and 11 of Landsat 8 OLI/TRIS). Subsequently, Kelvin values of the top of atmosphere brightness temperature were transferred to Celsius ( • C). The proportion of vegetation and then emissivity corrected LST maps were derived using Equation (6): where, T b = Landsat 7 ETM band 6, Landsat 8 OLI/TRIS band 10 brightness temperature; λ = wavelength of emitted radiance (for band 6, 11.5 µm); and ( [11]. Land surface emissivity (ε) was calculated using Equation (7): where, ε s and ε v are the soil emissivity and vegetation emissivity; F is a form factor, the mean value of which assumes that various geometrical distributions is 0.55, respectively [34]. P v , the vegetation proportion, was derived using Equation (8): (8) where, NDVI is the Normalized Difference Vegetation Index; the NDVI min and NDVI max = the minimum value of NDVI and the maximum value of the NDVI, respectively [35]. NDVI was derived using Equation (9): where, NDVI was derived using the surface reflectance of Landsat 7 ETM band 3 (ρRED) and band 4 (near-infrared (ρNIR)), Landsat 8 OLI/TRIS band 4 (ρRED), and band 5 (ρNIR) [11].

LUC Classifications
The LUC of the study area was classified using three machine learning methods facilitating R software: K-nearest neighbor, support vector machines, and random forest [36]. We employed supervised learning for the classification. Five types of LCCs were considered in this classification, including, (i) Built up (consists of the impervious surface such as buildings, roads, airports, schools, and industrialized areas.); (ii) Green 1 (consists of forest areas, woodland, mangrove, and bushland with big tree crowns); (iii) Green 2 (includes grassland, and cropland scattered green space); (iv) Bare land (includes exposure areas without any natural and human-made areas); and (v) Water (consists of rivers, lake reservoirs, lagoon, and wetlands). We applied a random sampling method to cover all LUC categories in the study area and produced 400 points each year. Google Earth's historical images were utilized as reference data for accuracy assessment. Then, we made three LUC maps for each year, based on five LCCs.
We ranked the classified LUC maps based on the highest value of overall accuracy and made the LUC maps generated by K-nearest neighbor (overall accuracy was over 88% each year). Afterward, we conducted the post-classification corrections as the majority filter and hybrid classification method to avoid misclassification errors and salt and pepper noises [37][38][39].

Tasseled Cap Transformation (TCT)
Only LUC classification is not enough to understand the spatial formation of LST in Accra. Therefore, we use TCT in this study because TCT gives useful information on brightness, greenness, and wetness in Accra. Using the TCT concept, which was first introduced by Kauth, R.J. and Thomas G.S. [40], we can measure ecological footprints and urban development with high-resolution multispectral sensors (multispectral scanner system (MSS), Landsat TM, ETM+, Landsat OLI/TIRS, and Quick bird). The first component is brightness (orthogonal to the brightness), which is comprised of bare soil and cover up to some extent, built up areas, and natural environments such as asphalt, rock, gravel, and other bare places. Greenness, the second component, is related to healthy vegetation. The third is wetness, being orthogonal to brightness and greenness. The wetness category is associated with hydro-related features such as soil moisture, water, and wet elements. Brightness, wetness, and greenness of the TCT consist of 97% important information of the image. The Brightness Index (BI), Greenness Index (GI), and Wetness Index (WI) of Landsat 7 ETM and Landsat 8 OLI/TRIS are derived using Tables 1 and 2. Finally, the TCT components of BI, GI, and WI are rescaled to 1-100 using a rescale function tool available in ArcGIS 10.6.

Spatial Analysis of Accra
Many studies related to the UHI have examined the usefulness of the satellite remote sensing data to explore the relationship between landscape pattern and LST [30,[43][44][45][46][47][48][49]. Recently, urban-rural gradient analysis, multiresolution grid-based analysis, and landscape configuration analysis have been conducted in UHI studies. Each technique gives some insights into a better understanding of the UHI in the study area, which provides useful information for urban designers and urban planners to mitigate the UHI effect. In the spatial analysis, we conducted three methods, i.e., urban-rural gradient analysis, multiresolution grid-based analysis, and landscape metrics-based analysis. We excluded water and bare land categories from this section because (i) we already calculated the WI in the study area in Section 2.5, and (ii) bare land showed a very low amount compared to the other LCCs in the study area.

Urban-rural Gradient Analysis
We examined spatial differences between LST and LCCs along the urban-rural gradient. First, we prepared the multiple ring buffer zones from the city center, and 50 ring buffer zones were created based on a 500 m distance interval. Second, we extracted the mean LST and the density of LCCs. The same 50 ring buffer zones were applied to extract TCT of the study area.

Multiresolution Grid-Based Analysis
This analysis enables us to explore the LCCs' influence on the mean LST using polygon grids. The statistical relationship between LCCs and mean LST was calculated to identify the impact of LCCs on the mean LST. However, a suitable grid size is not yet identified in the literature, and they are still discussing what the proper grid size is. Myint  After creating these polygon grids, we measured the mean LST and the density of LCCs in Accra. The set of multiresolution polygon grids was used to examine the TCT's influence on the mean LST of the study area. An attempt was made to investigate the most suitable grid size to simulate and predict meteorological and physical attributes using statistical nexus between atmospheric and surface temperature and LCCs in Accra.

Landscape Metrics-Based Analysis
This analysis enables us to examine the spatial traits of patches of LCCs influencing on the LST. Based on FRAGSTATS v4.2 software, we calculated the mean patch area of LCCs. We got Built up of 3.94 ha and 7.49 ha, Green 1 of 2.59 ha and 2.48 ha, Green 2 of 3.92 ha, and 3.39 ha in 2003 and 2017, respectively. Considering these results, multiresolution polygon grids analysis in Section 2.6.2 was not suitable for landscape metrics-based analysis. Based on previous studies [11,[51][52][53], we created a 3 km × 3 km fishnet to divide the study area. Forty out of 116 polygon grids were selected randomly for further investigation. We examined the four-landscape metrics, namely, mean patch area (AREA_MN), mean shape index (SHAPE_MN), largest patch index (LPI), and aggregation index (AI) ( Table 3) by employing FRAGSTATS v4.2 software. The analytical parameter was 8-cell neighborhood in 2003 and 2017, respectively. Bivariate correlation analysis was applied to determine the statistical relationships between LCCs and mean LST in 3 km × 3 km grid size. Table 3. Four spatial metrics used in this analysis [54].

Spatio-Temporal Change of LUC and LST (from 2003 to 2017)
The LUC maps show that Accra experienced rapid urbanization from 2003 to 2017 ( Figure 2). Built up increased from 30.51% to 51.90%, and Green 1 decreased from 22.45% to 11.26% during the study period. Green 2 decreased from 46.71% to 32.49%. Overall, the rapid expansion of Built up has led to changes in the landscape pattern in the study area. It can be identified that the spatial dispersion of land cover categories (Figure 2) determines the spatial variation of the land surface temperature of Accra ( Figure 3).
The LUC maps show that Accra experienced rapid urbanization from 2003 to 2017 ( Figure 2). Built up increased from 30.51% to 51.90%, and Green 1 decreased from 22.45% to 11.26% during the study period. Green 2 decreased from 46.71% to 32.49%. Overall, the rapid expansion of Built up has led to changes in the landscape pattern in the study area. It can be identified that the spatial dispersion of land cover categories (Figure 2) determines the spatial variation of the land surface temperature of Accra ( Figure 3).
The    The LUC maps show that Accra experienced rapid urbanization from 2003 to 2017 ( Figure 2). Built up increased from 30.51% to 51.90%, and Green 1 decreased from 22.45% to 11.26% during the study period. Green 2 decreased from 46.71% to 32.49%. Overall, the rapid expansion of Built up has led to changes in the landscape pattern in the study area. It can be identified that the spatial dispersion of land cover categories (Figure 2) determines the spatial variation of the land surface temperature of Accra ( Figure 3). The

TCT; LCCs vs. LST
The TCT maps and the result of correlation analysis are shown in Figure 5 and Appendix A1. In 2003, BI was mostly accumulated in the eastern areas; GI was mainly accumulated in the northeast direction from the city center, and WI was concentrated in the surrounding zones of the city center. However, in 2017, BI was spread out widely ( The results indicate that the density of Built up and the mean LST had a consistent trend. The ring zones near the city center had the highest mean LST, which was decreased along the urban-rural gradient both in 2003 and 2017 (Figure 6a). In contrast, the density of Green 1 and Green 2 shows the lowest mean LST closest to the city center. The density increased away from the city center showing the majority of the rural area both in 2003 and 2017 (Figure 6a).
The BI did not have the highest LST in the city center. This trend gradually raised radially, and again fell near the 18km zone (Figure 6b). The WI and the mean LST had a consistent pattern, and both fell away from the city center. It is characteristic that WI was dominant in urban areas. However, GI showed the opposite trend along the urban-rural gradient. In 2017, the mean LST was low near the city center and continuously increased until the 5 km zone along the urban-rural gradient ( Figure  6b).

TCT; LCCs vs. LST
The TCT maps and the result of correlation analysis are shown in Figure 5 and Appendix A. In 2003, BI was mostly accumulated in the eastern areas; GI was mainly accumulated in the northeast direction from the city center, and WI was concentrated in the surrounding zones of the city center. However, in 2017, BI was spread out widely ( The results indicate that the density of Built up and the mean LST had a consistent trend. The ring zones near the city center had the highest mean LST, which was decreased along the urban-rural gradient both in 2003 and 2017 (Figure 6a). In contrast, the density of Green 1 and Green 2 shows the lowest mean LST closest to the city center. The density increased away from the city center showing the majority of the rural area both in 2003 and 2017 (Figure 6a).
The BI did not have the highest LST in the city center. This trend gradually raised radially, and again fell near the 18km zone (Figure 6b). The WI and the mean LST had a consistent pattern, and both fell away from the city center. It is characteristic that WI was dominant in urban areas. However, GI showed the opposite trend along the urban-rural gradient. In 2017, the mean LST was low near the city center and continuously increased until the 5 km zone along the urban-rural gradient (Figure 6b).          mean LST and WI and the mean LST as grid size increased.

Land Cover Density and TCT vs. Multiresolution Grid-Based Analysis
Conversely, the correlation results were decreased between the mean LST and GI as the grid size decreased. BI can make a heating effect in closer vicinities. The trend of coefficient value between the BI and the mean LST increases when the grid size becomes smaller. In contrast, GI can make a cooling environment in closer proximity. This explains that the coefficient value between the mean LST and GI lowers when the grid size becomes larger.  Notably, the density of Built up and the mean LST had a strong relationship when grid size increased both in 2003 and 2017 (Figure 7). However, the decreasing trend at the coefficient value between the mean LST and the Green 1 shows that the impact of the Green 1 on the mean LST became stronger when the grid size became smaller. Green 1 has a cooling effect in closer contiguities. It indicates that the coefficient value between the mean LST and Green 1 reduces when the grid size becomes larger. Green 2 and the mean LST coefficient values also decrease when the grid size becomes larger in the two-time points.

Landscape Metrics vs. LST
Among the grid-based analysis of the TCT, the mean LST and BI had a strong and positive correlation as grid size became larger. Mean LST and GI and WI had a strong negative relationship as grid size became smaller both in 2003 and 2017. In general, BI can retain more solar radiation to enhance the surface temperature. In contrast, GI and WI have a cooling effect of supporting evaporation. These mechanisms explain why the correlation became stronger between BI and the mean LST and WI and the mean LST as grid size increased.
Conversely, the correlation results were decreased between the mean LST and GI as the grid size decreased. BI can make a heating effect in closer vicinities. The trend of coefficient value between the BI and the mean LST increases when the grid size becomes smaller. In contrast, GI can make a cooling environment in closer proximity. This explains that the coefficient value between the mean LST and GI lowers when the grid size becomes larger.    Table 4). The large patches of Green 1 can make a robust cooling environment than the small patches of Green 1 [55][56][57][58]. Our study also reveals that the large patches of Built up produce a significant heat effect, while the smaller patches of Built up make the relatively lower UHI effect.

Impacts of Urbanization
The findings show that urbanization took place mainly towards east and west directions (Figure 2). It also spread toward the north direction. Some studies mention that Accra has experienced rapid urbanization [24,59]. The number of recorded inhabitants was 1.9 million in 1985 [60], and the population of the region had become 4.6 million in 2016 [61]. The annual population growth rate was 3.5% between 2000 and 2017, and the projected population will be doubled by 2040 [62]. This situation has had not only positive impacts on their inhabitants, socially and economically but also considerable adverse effects on the natural environs. The urbanization effect of Accra has been weakened by the existence of Green 1 and Green 2 ( Figure 2). The different push and pull factors are influencing the urban growth pattern in Accra [62]. For example, push factors are the high rent price in the downtown of Accra and flood occurrences in the Accra core area. The government projects on the redevelopment plan are extending from the urban core area to the urban periphery [23,63]. There are a few pull factors. The availability of free lands to build houses is one of them [23,62]. The study results prove other impacts of rapid urban growth, including the UHI phenomenon.

LST Pattern Along the Urban-Rural Gradient
Based on the urban-rural gradient analysis, the zones near the city center did not show the highest mean LST in 2003 and 2017. According to that, the density of Built up and Green 1 was strongly correlated with each other, both in 2003 and 2017. The highest WI values showed closer to the city center along the urban-rural gradient. The neighborhood water bodies, other water resources, and humidity can result in the cooling environment around the city center. Therefore, the city center and surrounding environment did not have the highest mean LST value in the given study period ( Figure 6). However, there are three peaks of the mean LST, which can be observed from 2 km to 7 km, from 13 km to 18 km, and from 22 km to 23 km from the city center. The population growth and the development of multicore urban structures from 2003 to 2017 accounted for enhancing the LST in these zones. The visual interpretation shows more fragmentation of the Built up area in the classified 2003 map and less fragmentation of the Built up area in the 2017 map. However, the real situation indicates that the Built up area in 2003 was less fragmented and the Built up area in 2017 was more fragmented due to high rise buildings and the expansion of concentric urban centers away from the city center. These kinds of human activities enhanced the spatial formation of LST in the given zones.
The lowest mean LST can be observed from 20 km to 21 km, with a contrast of 1.4 • C in 2003 and 0.28 • C in 2017. The formation of the mean LST along the urban-rural gradient of Accra shows typical characteristics of the UHI, with LST steadily declining along the urban-rural gradient with physical features [30]. Similar findings were reported by   [11] in the UHI study in three megacities of Southeast Asia. The Built up density was high around these zones, and the Green 1 density was very low. However, from 22 km to 23 km, the peak area had a high density of Green 2. The Green 2 category is vital for Accra because it has a higher mean LST than Green 1. Therefore, the temperature of Green 2 continuously influences the LST in the region. City planners of Accra should pay much attention to the existence of Green 2 areas in their urban design. As a result, GI of the TCT had the lowest value from 1 km to 17 km, and high BI value in these zones accounted for increased LST in the area. The WI gradually fell along the urban-rural gradient.
The decreasing trend of the mean LST was evident from 7-8 km to 12-13 km ( Figure 6) and "Densu Delta Protected Area", "Pambros Salt Basin", and "Achimota Forest Reserve" were positioned in corresponding zones (Figure 1d). These natural environments have plenty of vegetation, including subtropical trees, grassland, and bushland. This is a cause of the cooling effect near the city center. Another decreasing trend of the mean LST was prominent from 18 km to 21 km from the city center. "Weija Reservoir," and "Sakumo Lagoon Protected Area" were located in these areas (Figure 1d). The vast water reservoir and plenty of mangrove forests in the lagoon area have caused the cooling effect of the city. These results were similar to other studies that analyzed the influence of natural environs on the LST cooling [18,64,65]. Urban planners need to pay much attention to the central area of Accra to mitigate the UHI effect. The massive size of water areas and surrounding natural zones should be preserved. We identified that the LST value of Accra in 2003 and 2017 was dependent on not only LCCs, but also on different environmental factors such as significant air masses, dust, humidity, solar radiation, and air moisture when the Landsat thermal band imageries were investigated.

Density of Land Lover Categories and TCT vs. LST at Multiresolution Grid-Based Analysis
The multiresolution grid-based analysis explored that Built up density and Green 2 density were most associated with the raising of the LST (Figure 7). However, Green 1 was the most relevant cooling effect on the LST.   [11] found similar results using the data on an impervious surface and green space in three megacities in Southeast Asia (a suitable grid size was 210 × 210, and 240 × 240 for Bangkok, Jakarta, and Manila and all results were statistically significant at the (ρ < 0.001). Xiao et al. (2007) [47] made clear that the relationship between impervious surface and mean LST became stronger while grid size enlarged around 30 m to 960 m. All of the correlations were significant (ρ ≤ 0.05).
Testimonies regarding the mean LST in each TCT analysis with multiresolution grids indicate that BI has improved the UHI of Accra (Figure 8). The BI of the TCT map shows that the higher BI had emitted of Built up category compared to the LUC maps in 2003 and 2017 (Figures 3 and 5). By contrast, the GI and WI have visible evidence of the cooling effect on the LST. The GI of the TCT map has the highest GI value in the northeast part of the study area. Compared to the LUC maps, these areas correspond to Green 1 of the regions (Figures 3 and 5). However, in comparison with the Green 2 category of LUC maps and the GI of TCT, we observe that Green 2 areas in the GI maps have a low degree effect (1-40) for cooling down the LST. Rasul et al. (2015) [66] focused on the urban cool island effect in Erbil, Iraqi Kurdistan, using TCT. They found a strong positive correlation with the mean LST and BI (r = 0.75) and the negative relationship with the mean LST and GI and WI (r = −0.39 and r = −0.9) using Landsat 8 OLI/TRIS data.
It is essential to mention that the correlation between all LCCs with the mean LST, and the correlation between TCT and the mean LST did not increase the equal quantity around the grid size of 210 m × 210 m or 240 m × 240 m. The slope of the correlation was flat after passing the mentioned grid size (Figures 7 and 8). The 210 m × 210 m and 240 m × 240 m were proper grids of LCCs that affected land temperature. On the other hand, these grid sizes are the most suitable for forecasting methodological parameters and ecological features employing statistical relationships in land temperature and atmospheric temperature and LCCs, as well as spatial collocation. The results were similar to other studies: Myint et al. (2010) [50] have studied the UHI effect in Phoenix in the USA, and Hao Hou et al. (2020) [51] have discussed the cooling effect by using landscape composition and configuration in Hangzhou, China.  [70] revealed that a patch shape with highly aggregate edges produced more substantial mitigating of the LST effect in Isfahan, Iran. Masoudi et al. (2019) [71] examined the correlation between the LST pattern and green areas in four Asian cities, namely Kuala Lumpur, Hong Kong, Jakarta, and Singapore, using 240 m, revealing that landscape configuration in Hong Kong and Kuala Lumpur had not affected the LST and the patches of the low fragmented green area were created in LST cooling in Singapore and Jakarta.
It is important to note that the rapid urbanization transforms natural environments into the Built up, which is retaining more solar radiation and releasing less solar radiation reflected [3]. According to that, LST influences the land temperature, but its surrounding air temperature and other environmental factors (wind/humidity) impact human health and thermal comfort. Vegetation and shadow help to reduce surface temperature. One solution to minimize the UHI is to grow plenty of trees and grass in the study area. Based on SHAPE_MN, we found that Green 2 could act to provide a lower surface temperature (Table 4). On the other hand, tree covers can create an urban cool island condition by evapotranspiration. In general, evapotranspiration by plants absorbs ambient heat and release water vapor to the atmosphere [11,13].
Furthermore, urban vegetation cover serves to enhance the air quality of urban environs, which help reduce air pollution and mitigate the UHI effect of vital for the city dwellers. The public parks and recreational facilities are needed to minimize the UHI effect. In this context, city planners should play a significant role when designing the urban spatial structure by paying attention to the UHI phenomenon. The scattered and fragmented vegetation has a smaller cooling effect than the gathered vegetation pattern. Large clusters of vegetation increase the urban cool island effect. Therefore, it is crucial for city planners to attain the suitable spatial allocation of the urban landscape by designing large clusters of vegetation and scattering building environs.

Conclusions
In this study, we examined the spatial process of UHI's formation based on the LST in Accra by facilitating various geospatial techniques of TCT, urban-rural gradient, multiresolution grid-based analysis and spatial metrics in 2003 and 2017. We found that the distribution of UHI has dispersed and fragmented over time because of suburbanization. The Built up area increased from 30.51% to 51.09%, and the Green 1 area decreased from 22.45% to 11.26%, and the Green 2 area decreased from 46.71% to 32.49% from 2003 to 2017.
Daytime mean LST is high in Accra because it is located in a tropical sub-Sahara. We revealed that Built up and Green 2 caused the heating effect, and Green 1 and Water caused the cooling effect in the study area. The results indicated that the Accra s heat core was located 1-2 km far from the city center inside the central business district. The WI of TCT was gradually decreased along the urban-rural gradient, indicating considerable heat stress around the suburb. The density of Built up and the mean LST had a strong relationship when grid size became larger, and the density of Green 1 and mean LST had a more robust relationship when the grid size became smaller. It is concluded that 210 m × 210 m and 240 m × 240 m grid sizes are the most suitable for forecasting parameters and ecological features employing statistical relationships in land temperature and atmospheric temperature and LCCs in Accra. The urban cool island effect is dominant in Densu Delta Protected Area, Achimota Forest Reserve, Weija Reservoir, and Sakumo Lagoon Protected in Accra.
Urban planners and policymakers should pay attention to Green 2 because this category causes a higher mean LST. These findings give some insights for understanding spatial variations of the urban thermal environment and their spatial characteristics in a coastal city in a tropical sub-Saharan climate. Also, we found that multiresolution grid-based and landscape configuration analysis is useful to predict the impacts of the LUC pattern on the future UHI. People are exposed to the increase in surface temperature in their life. Their daily socio-economic activities are significantly related to built up environs. In this case, it is better to control the UHI effect. Expanding and maintaining the large cluster of green covers is a critical time factor. We should consider the increase of dense vegetation areas, green roof, and scattered building environs to mitigate the UHI effect. Urban environments can have better-living conditions for city dwellers.