Urban Heat Island Formation in Greater Cairo: Spatio-Temporal Analysis of Daytime and Nighttime Land Surface Temperatures along the Urban–Rural Gradient

: An urban heat island (UHI) is a signiﬁcant anthropogenic modiﬁcation of urban land surfaces, and its geospatial pattern can increase the intensity of the heatwave effects. The complex mechanisms and interactivity of the land surface temperature in urban areas are still being examined. The urban–rural gradient analysis serves as a unique natural opportunity to identify and mitigate ecological worsening. Using Landsat Thematic Mapper (TM), Operational Land Imager/Thermal Infrared Sensor (OLI/TIRS) and Moderate Resolution Imaging Spectroradiometer (MODIS), Land Surface Temperature (LST) data in 2000, 2010, and 2019, we examined the spatial difference in daytime and nighttime LST trends along the urban–rural gradient in Greater Cairo, Egypt. Google Earth Engine (GEE) and machine learning techniques were employed to conduct the spatio-temporal analysis. The analysis results revealed that impervious surfaces (ISs) increased signiﬁcantly from 564.14 km 2 in 2000 to 869.35 km 2 in 2019 in Greater Cairo. The size, aggregation, and complexity of patches of ISs, green space (GS), and bare land (BL) showed a strong correlation with the mean LST. The average urban–rural difference in mean LST was − 3.59 ◦ C in the daytime and 2.33 ◦ C in the nighttime. In the daytime, Greater Cairo displayed the cool island effect, but in the nighttime, it showed the urban heat island effect. We estimated that dynamic human activities based on the urban structure are causing the spatial difference in the LST distribution between the day and night. The urban–rural gradient analysis indicated that this phenomenon became stronger from 2000 to 2019. Considering the drastic changes in the spatial patterns and the density of IS, GS, and BL, urban planners are urged to take immediate steps to mitigate increasing surface UHI; otherwise, urban dwellers might suffer from the severe effects of heatwaves.


Introduction
Urbanization impairs the natural landscape as it produces impervious surfaces. Rapid urbanization has become one of the most critical global issues in the 21st century [1][2][3][4]. Various socio-environmental problems, including climate change [5][6][7], energy systems [8], deforestation [9], water and air quality [10], and environmental health [11], have been attributed to large regions being urbanized too rapidly and without proper planning.
The urban heat island (UHI) phenomenon [12][13][14][15], which refers to higher temperatures in urban areas relative to the surroundings, has been studied in many cities around the world [15][16][17][18][19]. The distribution of impervious surfaces covered by cement, asphalt, and concrete raises the land's radiative surface temperature [20] and changes the humidity of urban areas [21,22]. The rise in heat in urban areas has caused various social problems such as increasing water and energy consumption [23], air pollution [24], discomfort, and

LUC Classification
We downloaded the atmospherically corrected pre-processed (level 2) Landsat data from the GEE platform in 2000, 2010, and 2019 [79]. First, we demarcated the 50 × 50 km boundary for the study area and prepared the boundary shapefile. The city center is located in the central business district (CBD). Second, the boundary shapefile was imported on Asset in GEE. Finally, we ran the script for Greater Cairo to download Landsat imageries for summer (July and August). In this stage, we ran the script three times for three time points. Using the Image Collection tool in GEE, we prepared the three final Landsat imageries, including four images for 2000 (Landsat 5), two images for 2010 (Landsat 5), and four images for 2019 (Landsat 8) (Appendix A Table A1) for Greater Cairo. All downloaded images were projected onto the WGS84/UTM 36N projection system before further processing.
The LUC maps of the study area were constructed by applying four machine learning methods: K-nearest neighbor (KNN), artificial neural network (ANN), support vector machine (SVM), and random forest (RF), facilitated by the R software [80]. In the LUC classification, bands 4, 3, and 2 for Landsat 5 and bands 5, 4, and 3 for Landsat 8 were used. The spatial resolution of the prepared LUC maps was 30 × 30 m. Four LUC types were derived from this classification: (i) bare land (BL; associated classes: desert area and stone land); (ii) green space (GS; associated classes: forest, cropland, grassland, and shrub); (iii) impervious surface (IS; associated classes: all kind of impervious surface areas including buildings, roads, and airports); and (iv) water (W; associated classes: all kinds of water bodies, e.g., rivers, and ponds). Google Earth's historical images were used as reference data for accuracy assessment. Accurate results were generated by automatic sampling in the algorithm using 400 points each year. We constructed four LUC maps for each year based on the four LUC types.
In the next step, we ranked the classified LUC maps based on the highest overall accuracy value and determined the LUC patterns generated by K-nearest neighbor (the overall accuracy was over 90% each year) (Appendix A Tables A2-A4). Subsequently, we applied the post-classification corrections as the majority filter and hybrid classification method to avoid misclassification errors and salt and pepper noise [81][82][83].

MODIS Data
The MODIS sensor is onboard on the Terra satellite, which was launched by the National Aeronautics and Space Administration (NASA) Earth Observing System (EOS) in December 1999 [84]. The Terra-MODIS satellite captures images for 10:30 a.m. (local solar time) in descending mode and 10:30 p.m. in ascending mode [85]. However, the Terra-MODIS acquisition time is a nominal time that varies with the location [86]. The swath width of MODIS instruments is 2330 km, and they observe the entire planet every one to two days. Swath acquisition produces outputs every 5 min. MODIS sensors include 36 spectral bands of the electromagnetic spectrum, with visible light and infrared radiation. MODIS instruments capture data at three spatial resolutions and various temporal resolutions (i.e., spatial resolutions: bands from 1 to 2250 m; bands from 3 to 7, 500 m; bands from 8 to 36, 1000 m; and temporal resolutions: daily, 4 day, 8 day, 16 day, monthly, quarterly, and annually) [87].
MODIS LST products obtain LST at high temporal (daily) and low spatial resolutions (1 × 1 km), gridded at intervals of sinusoidal projection. This study used Terra-MODIS LST version 6 product data to investigate the surface UHI in Greater Cairo [88]. Terra-MODIS LST data were downloaded from the GEE platform [89] for three time points to characterize LST along the urban-rural gradient. In this stage, we used ee.reducer to extract each pixel's mean temperature values in 2000, 2010, and 2019, available on the GEE platform. The temperature value was obtained from the MOD11_L2 product [88]. MODIS bands 31 and 32 and six quality indicator layers provided land surface temperatures for both daytime and nighttime. Bit flags were used to manage the quality control parameters, i.e., mandatory quality assessment flags, data quality flags, emissivity quality flags, and LST error flags [88].
Greater Cairo summer was defined as the period from July to August. We used the same Cairo boundary shapefile to download MODIS data. Subsequently, all downloaded MODIS data were projected onto the WGS84/UTM 36N projection system. We used ArcGIS MODIS Python Toolbox to pre-process the data [90].

MODIS LST and Density of IS, GS, and BL
The LUC density is apparent in the cross-sectional analysis of the urban-rural gradient of a typical surface UHI profile. Based on the summer daytime/nighttime mean LST, we created a surface UHI intensity profiling for this study. First, kilometer 0, the city center, was allocated ( Figure 1c). Second, raster daytime and nighttime LST maps were created by snapping together polygon grids. These polygon grids were 1 × 1 km, the same as in previous studies [16]. The densities of IS, GS, and BL were defined as the percentage of the total area within a 1 × 1 km grid, which is equal to the MODIS LST data's spatial resolution. Finally, the relationship between the mean LST and the urban intensity measured on each grid was examined using bivariate correlation analysis and scatter plot diagrams. We excluded the water category in this stage.

Trend in the Daytime and Nighttime Surface UHI Intensity
We examined the mean LST difference (∆ mean LST) along the urban-rural gradient between the daytime and nighttime. All 1 km polygon grids in the same direction were targets of analysis. Four polygons of the surrounding city center (hereafter referred to as the central grid area = 4 km 2 ) were defined as Urban-Rural Zone 1 (URZ 1 ) ((Appendix A Figure A1). The 25 buffer areas were delineated as urban-rural zones (URZs), e.g., URZ 1 , URZ 2 , . . . , URZ 25 .
Zones derived the daytime and nighttime mean LST and densities of IS, GS, and BL in each URZ at 1 km intervals. The daytime/nighttime change in mean LST between URZ 1 (the zone with the highest urban intensity each year) and other URZs (i.e., URZ 2 , URZ 3 , ..., URZ 25 ) were defined by the surface UHI intensity. The same procedure was applied for the three time points.
We focused on the URZ with a high IS density as urban zones and URZs with <15% IS density as rural zones [16,91]. The trend in the daytime and nighttime surface UHI intensity was calculated using Equation (1). All the extracted values were plotted along the urbanrural gradient. Scatter plot diagrams were drawn to delineate statistical relationships.

Population Density Data
WorldPop is a database for estimating the world's population [92,93]. It provides various types of the gridded population data sets [94]. Numerous researchers have used this population data [95][96][97][98], including United Nations Development Programme (UNDP), World Health Organization (WHO), The World Bank, and the World Wildlife Fund (WWF) [99]. The large population or high population density (PD) of a city contributes to settlement expansion and vertical development of high-rise residential apartments. Indirectly, the PD has become one of the critical factors of urban heat island formation. In this study, 2000, 2010, and 2019 WorldPop data (at a 30 arc-second approximately 1 × 1 km resolution) were used to explore the relationship between mean LST and PD in both the daytime and the nighttime in Greater Cairo using scatter plots diagrams and linear regression. We used a 1 km polygon grid, produced in Section 2.4, to extract raster WorldPop PD values.

Landscape Configuration Analysis
Based on the landscape configuration analysis, we examined how IS, GS, and BL influenced the variability in the LST in Greater Cairo. The grid size used in Section 2.5 was insufficient for this configuration analysis due to the total patch's possible effect being allocated to 1 km grids. Moreover, instead of identifying a single patch's optimum size, we found the best configuration of patches in closer proximity. Considering previous studies [12,54,100] and the reasons mentioned above, we applied a relatively large grid size. We used a 5 × 5 km fishnet to divide the whole study area into 100 sub-parts. Based on other studies [12,54], 25 sub-parts were selected using a random sampling technique. The selected sub-parts represented 25% of the total population. The LUC map in each study area was clipped with fixed 25 polygon grids for further analysis. We determined three class-level spatial matrices: mean patch area (AREA_MN), largest patch index (LPI), and aggregation index (AI) ( Table 1). These spatial metrics have been widely used in previous UHI studies [12,101]. The 8-cell neighbor rule was applied to calculate the three metrics. To analyze the impact of landscape configuration on LST, the generated metric values were compared with the mean LST of the LUC classes of each sub-part. Table 1. Class-level metrics used in this study [102].

Spatial Metrics Formula Description Units
Mean patch area (AREA_MN) AREA_MN (ha) = The spatial pattern and heterogeneity of the area. ha LPI ability to detect the advantages of the LUC. 0-100 The calculation of class-level aggregation in the area. percentage Where n i = number of patches of land use/cover (LUC) class i; n = number of patches; j = total of the specific patch type; x ij = patch metrics value of patch ij; A = total area of LUC; a ij = total pixels of patch area ij; g ii = number of joins between pixels of class type; and max → g ii = maximum number of joins between pixels of class type.

LUC Changes and Magnitude and Trends of LST
The classified LUC maps in 2000, 2010, and 2019 show that Greater Cairo has undergone rapid urbanization over the 19-year period ( Figure 2). From 2000 to 2010, the IS area increased from 564.14 to 698.65 km 2 ; and from 2010 to 2019, the IS area increased by 171 km 2 , with a total gain of 305.21 km 2 ( Table 2). The increase in the IS area in the 2010-2019 period was faster than in the 2000-2009 period. The GS and BL areas showed a decrease, i.e., a total net loss of 71.83 km 2 and 229.19 km 2 from 2000 to 2019, respectively. Overall, the IS area has greatly enlarged by rapid urbanization in Greater Cairo.  Figure 3 shows the density of IS, GS, and BL with a 1 × 1 km grid size. Figure Figure 5 indicates the derived mean LST of four LUC types: IS, GS, BL, and W (water). In the daytime, BL had the highest mean LST in the three time points. The mean LST difference (daytime-nighttime) of IS was 15.23, 15.63, and 14.59 • C, while that of BL was 21.58, 18.45, 18.99 • C in 2000, 2010, and 2019, respectively. In the nighttime, the mean LST of the four LUC types was lower compared with the daytime mean LST. Though GS showed the lowest mean LST among the four types at the three time points, the difference in the four types was not as large in the nighttime. This means that the impacts of IS and BL on the LST increase were less influential at night.

Magnitude and Trend of the Surface UHI Intensity in the Daytime and Nighttime
The highest density of IS was in URZ 1 over time (>96.39% in 2000, >97.98% in 2010, and >99.21% in 2019) (urban zone), whereas the boundaries between the urban and rural zones with <15% density of IS were in URZ 19 in 2000, URZ 23 in 2010, and URZ 25 in 2019 ( Figures 6 and 7). The urban-rural gradient analysis showed an almost similar trend in the three periods. The composition ratio of the IS in 1-6 km zones decreased, while the composition ratio of the IS in 7-25 km zones increased from 2000 to 2019. In the daytime, the central grid area (URZ 1 ) had the lowest mean LST ( Figure 6a). Then, the mean LST gradually rose with an increase in distance. However, the mean LST tended to drop slowly from 6 to 9 km, and after 10 km, it began to rise again. Along the same line, the IS density and BL density revealed a consistent pattern with mean LST. GS density showed an inverse relationship with mean LST. The statistical analysis based on the 25 URZs indicated positive correlations between the mean LST and the IS density in the daytime. The GS density showed a high negative correlation, and the BL density had a high positive correlation with the mean LST (ρ < 0.001) in the daytime (Figure 6b).        Figure 7 shows the urban-rural gradient pattern in the nighttime. In 2000, 2010, and 2019, the mean LST increased gradually between URZ 1 and URZ 7 (i.e., approximately 8 km), and then, the 2000 mean LST decreased rapidly after URZ7, whereas the 2010 and 2019 decreases were more gradual (Figure 7a). After URZ 13 , the mean LST gradually increased until URZ 25 at the three time points. We conclude that the area of surface UHIs in the nighttime enlarged spatially over time. Overall, the correlation analysis indicated a high positive correlation between the mean LST and the density of IS in the nighttime. Conversely, the density of GS and the density of BL had a high negative correlation with the mean LST in the nighttime (ρ < 0.001) (Figure 7b). Notably, the correlations of mean LST with IS density in the nighttime were higher than in the daytime in the three periods. Figure 8 shows the actual situation of surface UHI intensity in the daytime and nighttime in greater Cairo. We use negative values of daytime (based on Equation (1) (1)) as negative values (Figure 8b) along the urban-rural gradient because, in the daytime, the urban zones in Greater Cairo (central grid area) experienced an urban cool island effect; in contrast, in the nighttime, urban heat islands were recognized, and the influence became stronger over time. The analysis revealed that the change in mean LST between URZ 1 to URZ 25 gradually increased with distance in the daytime, whereas the change in mean LST between them gradually decreased in the nighttime. The difference in daytime surface UHI intensity (between the urban zone and rural zone) was 3. of the four LUC types was lower compared with the daytime mean LST. Though showed the lowest mean LST among the four types at the three time points, the differe in the four types was not as large in the nighttime. This means that the impacts of IS a BL on the LST increase were less influential at night.

Magnitude and Trend of the Surface UHI Intensity in the Daytime and Nighttime
The highest density of IS was in URZ1 over time (>96.39% in 2000, >97.98% in 20 and >99.21% in 2019) (urban zone), whereas the boundaries between the urban and ru zones with <15% density of IS were in URZ19 in 2000, URZ23 in 2010, and URZ25 in 2 ( Figures 6 and 7). The urban-rural gradient analysis showed an almost similar trend   Figure 7 shows the urban-rural gradient pattern in the nighttime. In 2000, 2010, and 2019, the mean LST increased gradually between URZ1 and URZ7 (i.e., approximately 8 km), and then, the 2000 mean LST decreased rapidly after URZ7, whereas the 2010 and 2019 decreases were more gradual (Figure 7a). After URZ13, the mean LST gradually in- the nighttime enlarged spatially over time. Overall, the correlation analysis indicated a high positive correlation between the mean LST and the density of IS in the nighttime. Conversely, the density of GS and the density of BL had a high negative correlation with the mean LST in the nighttime ( < 0.001) (Figure 7b). Notably, the correlations of mean LST with IS density in the nighttime were higher than in the daytime in the three periods.   (Figure 8 b). In the daytime, we observed a correlation between the change in mean LST and the change in density of IS (positive), the change in density of GS (negative), and the change in density of BL (positive) ( Figure  8a). In the nighttime, we observed a strong relationship between the changes in mean LST and density of IS (positive), the change in density of GS (negative), and the change in density of BL (negative) ( < 0.001; Figure 8b).   ulation of Greater Cairo was shifting from the city core outward. The relationship between the mean LST and PD was not well-correlated in the daytime (coefficient of determination (R 2 ) = -0.0033 in 2000, R 2 = -0.0002 in 2010, and R 2 = 0.0177 in 2019). However, a positive correlation between mean LST and PD was indicated in the nighttime (R 2 = 0.2418 in 2000, R 2 = 0.1643 in 2010, and R 2 = 0.1681 in 2019). All the data analyses were statistically significant ( < 0.001). We estimated that energy use in the nighttime would cause the warming of the densely habited regions and industries.

Spatial-Metrics-Based Analysis vs. LST
The LUC types had a possible nexus with LST (Table 3), and all the results were statistically significant (p = 0.000). In the daytime, the three indices of AREA_MN, LPI, and AI were positively correlated with the mean LST for BL, but they were not as strongly correlated with IS in the three periods. Conversely, the results showed the opposite condition with IS in the nighttime from 2000 to 2019. All three indices showed negative relationships with mean LST for GS both in the daytime and nighttime (Table 3) The increase in LPI values of IS in both the daytime and the nighttime shows that mean patches of IS were less fragmented in 2010 than in 2000, and they became more fragmented by 2019. Our findings also indicated that large IS patches promoted a significant heat effect, whereas smaller IS patches generated a lower surface UHI effect. The LPI results in the daytime showed that the mean patches of BL were more fragmented in 2010 than in 2000. It became less fragmented by 2019, indicating the strong positive impacts of large patches with the mean LST. Conversely, these values indicated a low positive influence in the nighttime.

Rapid Urbanization and Its Impact on Greater Cairo
Previous studies showed that rapid urbanization increased informal settlements and environmental degradation in developing countries [103,104]. According to the World Urbanization Prospects Report, Greater Cairo's population is projected to rise from 5.7 million in 1970 to 14.7 million in 2025 [105]. In 1969, the President of Egypt proposed a master plan for establishing new towns on the fringes of the Cairo desert area [72]. As a result of this project, Greater Cairo has undergone scattered urban expansion. It is estimated that many surrounding desert areas will change to built-up environments in the not-too-distant future, primarily in the eastern side of Greater Cairo. Our results showed that Greater Cairo had experienced rapid urbanization in the 19-year period. Similar results reported by Mohamed and Worku [106] for Addis Ababa and its surrounding environment in Africa (built-up areas increased 3.7% in 2005, 5.7% in 2011, and 7% in 2015, whereas natural environment and agricultural lands were in continuous decline. Siddiqui et al. [107] found that the urban growth rate increased (4.6% in 1993 to 26% in 2013) in Uttar Pradesh of the Indian metropolitan, showing scattered and infilling urban expansion. Han and Jia [108] found that urban areas grew by 590 km 2 from 1995 to 2015, with a 4% annual growth rate, while agricultural areas declined to 397 km 2 by 2015 in Foshan, China. The urban development policies for Greater Cairo will enable the use of desert areas for urban sustainability through the advances in industrial zones and the transportation system with balanced control of urban expansion [109].
Greater Cairo's urbanization has transformed the natural landscape to IS areas, including buildings, roads, and other human-made surface materials, enhancing the surface UHI effect. Surface features such as buildings, roads, and other IS areas can absorb more solar radiation than natural surface areas [12]. Due to daytime solar radiation, these IS areas absorb more solar energy, but the absorbed solar energy is released during the night. As a result, the surface UHI effect over city areas is more pronounced than in the surrounding natural areas during the nighttime because, at night, there is no solar energy. Still, the urban core area of Greater Cairo shows an urban heat island effect (Figures 4 and 8), and the green spaces, which help to reduce urban heat at nighttime, of city area is relatively small compared with the surrounding area ( Figure 8). The outcome of the study analysis revealed the significant influence of urbanization on the spatial intensity of the surface UHI effect.

Surface UHI Nexus with LUC Classes and PD
In Greater Cairo, studying surface UHI in both daytime and nighttime is critical because the city is located in a hot desert region, and a large portion of the area is covered by desert sand (Figure 2). Another important factor is its LUC composition: the northern part of Greater Cairo is covered by green areas (Figure 2). Our findings revealed that the surface UHI gradually increased from the central grid (URZ 1 ) in the daytime and gradually decreased in the nighttime. There was a significant difference between LUC categories in the daytime (Figure 8). Conversely, there was no substantial difference between them in the nighttime ( Figure 5). The surface UHI effect in the nighttime did not fully correspond with LUC categories (Figures 4 and 5), indicating that another mechanism works to produce the nighttime surface UHI effect. Therefore, it is crucial to examine why these phenomena happen in Greater Cairo.
First, the urban cool island effect was prominent in the daytime CBD. The decline in IS density and increase in GS density in the central area promoted the cooling effect in the daytime ( Figure 6). Cairo's government announced a provisioning service in the urban green areas such as green roofs and urban agriculture at the rooftop [110], green corridors, urban parks, green pedestrian, green parking for buses and taxis, as well as urban planning mainly in the urban core area and its environs [111]. These projects have promoted the urban cool island effect in the daytime. Building shading also affects LST more significantly than tree shading because there are more high-rise buildings in the city core area than trees and many buildings have a light roof, producing a positive effect on the energy balance. The suburban area showed high-temperature values in the daytime because IS density gradually increased with distance. New urban development projects have been constructed on desert sand without proper green space planning.
Second, the LUC categories did not substantially promote the surface UHI effect during the nighttime. As mentioned above, due to the daytime solar radiation, such IS areas absorb more solar energy, but in the nighttime, the absorbed solar energy is released. Another factor is anthropogenic activities. Our result revealed a positive correlation between PD and mean LST during the nighttime (Figure 9), indicating that anthropogenic heat released from industry, traffic, and air conditioners from apartments has promoted the surface UHI effect in the nighttime. Light-color paints, urban materials, and cool building materials would help to decrease the temperature. Urban water is also one solution to adopt to mitigate the surface UHI phenomenon in the nighttime [112]. Home gardens should also be promoted during landscape and urban planning to minimize the surface UHI effect during the nighttime.

Trend in Surface UHI Intensity along the Urban-Rural Gradient
Based on MODIS surface temperature data, the surface UHI effect was identified in the study area (Figure 4). Along with the urban-rural gradient pattern, the lowest LST appeared in the CBD. In contrast, the highest LST appeared 17 to 25 km away from the city center, with a difference in surface UHI intensity between urban zone and rural zone −3.45 • C in 2000, −3.87 • C in 2010, and −3.45 • C in 2019 in the daytime ( Figure 6). Unlike previous studies in Africa [54,113], the CBD in Greater Cairo experienced the urban cool island phenomena in the daytime. Conversely, the surrounding rural zones experienced a relative rise in the surface temperature compared with the CBD. Several studies around the world have reported the same phenomenon. For instance, Rasul et al. (2015) [114] examined the daytime urban cool island effect in Erbil, Iraqi Kurdistan, using Landsat 8 data. They reported that the urban cool island intensity in the CBD differed from 3.5 to 4.6 • C in comparison with a 10 km buffer zone around the city. Haashemi et al. (2016) [115] explored the seasonal changes in the urban heat island in the semi-arid city of Tehran, Iran, using MODIS and Landsat 8 data. Their findings indicated that surface urban cool islands remained in the daytime (the maximal urban-rural contrast was −4 Kelvin in March). The maximal nighttime value of the urban-rural difference was 3.9 Kelvin in May. Lazzarini et al. (2015) [71] examined urban climate modification (urban heat island study) in eight hot and semi-arid cities in the world, including Abu Dhabi in the United Arab Emirates; Kuwait City in Kuwait; Riyadh in Saudi Arabia; Doha in Qatar; Las Vegas and Phoenix in the U.S.; Biskra, Algeria; and Bikaner, Rajasthan. They revealed that six cities, Abu Dhabi, Kuwait, Las Vegas, Phoenix, Biskra, and Bikaner, showed the cool UHI effect in the central area during the daytime.
The magnitude of the surface UHI intensity changed considerably with LUC types (Figures 5-7). However, we found that the LST values of daytime and nighttime were dependent not only on LUC composition but also on different environmental factors such as significant air masses, dust, humidity, and solar radiation when the MODIS thermal band imageries were taken. In the daytime, the zones 8 to 10 km away from the central grid had the lowest drop in the mean LST because these zones had an abundance of the natural environment (e.g., MAZHAR botanical garden area, a large portion of the Nile River, and the JazĪrat Warrāq al Hadar area) compared with the surrounding areas. The Nile River and expansive green spaces were the sources of cool air in this zone. It is important to note that the urban cool island effect was dominant in the central grid area in the daytime. According to the density analysis, this area was almost covered with IS (96.39% in 2000, 97.98% in 2010, and 99.21% in 2019) (Figures 6 and 7). This means that the IS s effect temperature rise was less than that of BL. Furthermore, the central areas consisted of mixed land use of IS and GS with the river environment. The zones stretching from 3 to 7 km did not show a significant cooling island effect along the urban-rural gradient.
In the nighttime, the density of IS was the critical factor in the increase in the urban heat effect in the CBD (Figure 7). Thermal energy absorbed in the daytime is released at nighttime, contributing to the rise in surface UHI [12,116]. In contrast, we observed that the highest mean LST dropped in areas 12 to 16 km away from the city center because green-related land use was dominant in these zones (Figure 7), e.g., cropland and different kinds of desert vegetation such as low canopy trees, shrubs, herbs, and water-based plants (cactus). The daytime and nighttime and urban-rural anomalies were related to the LUC composition. In the nighttime, the impervious surface enhanced surface UHI intensity, combined with the effect of anthropogenic activities.
Though most surface UHI intensity studies have used IS and GS as explanatory variables [56,117,118], a few studies have focused on spatial urban-rural gradient patterns [72]. The variation in the change in mean LST along the urban-rural gradient provides a clue for future urban design and urban thermal mitigation strategies. In our study, the change in surface UHI intensity was examined through surface temperature variations (the change in mean LST) along the urban-rural gradient, i.e., between the density of IS, GS, and BL and between urban-rural zones in the daytime and nighttime at three time points (Figure 8). We compared the results of surface UHI intensities over time. We examined the overall increase and decrease in the surface UHI intensity in the daytime and nighttime in the study area (Figure 8). Along the urban-rural gradient, the rise in surface UHI intensity from 2000 to 2019 based on the urban zone (the highest density of IS) and the rural zone (the lowest density of IS < 15%) was -3.41 • C lower in the daytime and was 2.1 • C higher in the nighttime. Notably, the surface UHI intensity fluctuated (either increasing or decreasing) at critical thresholds along the urban-rural gradient.
To the best of our knowledge, no previous studies identified the critical threshold of surface UHI intensity for LST change along the urban-rural gradient. However, most studies proved that the surface UHI effect is most critical in summer [34]. We determined the threshold value of using URZs to clearly explain the mean LST variations in the urban and rural zones in the daytime and nighttime. The urban design is less aggregated in urban settings, an increase in irregular vegetation areas, coupled with small and medium water areas, would be a practical approach to mitigate the surface UHI effect. There is an urgent need for designing a sustainable urban landscape to avoid the risk of strong surface UHI and heatwaves.

Landscape Configuration on Surface UHI Formation
The AREA_MN, LPI, and AI indices were significantly correlated with mean LST. AREA_MN, LPI, and AI of IS did not substantially influence the daytime LST variation (Table 3). Conversely, the three above-mentioned metrics of IS showed a considerable positive influence on the nighttime LST variation in Greater Cairo, indicating that urban planners must pay attention to the IS category because it promotes the surface UHI effect in the nighttime. However, except for AI, AREA_MN and LPI of GS did not produce large differences in daytime and nighttime (Table 3).
Our findings are similar to those of other studies. For instance, Zhou et al. (2011) [23] revealed a strong correlation of mean patch size and mean shape index with mean LST for IS (positive) and for GS (negative) at the Gwynns Fall watershed in the USA. Li et al. (2012) [119] reported a negative correlation between mean patch size, mean shape index of patches of GS, and mean LST in Beijing, China. Among the three spatial metrics of IS and GS, mean LST showed significant correlations, positive for IS and negative for GS, in both daytime and nighttime in the study area (Table 3). Large patches of GS provide cooling, thus lowering LST more than smaller patches of GS.
Continuous and larger patches of IS provided a stronger surface UHI effect than smaller patches of IS in both daytime and nighttime. The area experienced a stronger surface UHI effect in the nighttime than the daytime due to continuous and large patches of IS around the central grid area. Human activities in large patches of the IS area have supported the increase in the nighttime temperature due to anthropogenic heat release, the use of nighttime air conditioners, and traffic in the urban core area. Additionally, BL's continuous and large patches produced a strong surface UHI effect in the daytime in Greater Cairo. In contrast, the same condition was sufficiently influential in producing the cool island effect in the nighttime in Greater Cairo (Figures 4 and 8).
In general, fragmented green areas are less effective in mitigating the surface UHI effect in urban areas. Therefore, we suggest increasing vegetation and ponds in urban areas, introducing green roofs and green walls, and other such practical strategies to mitigate the strong surface UHI effect in desert cities like Greater Cairo.

Conclusions
We examined the spatial change in the local climate of Greater Cairo influenced by the surface UHI phenomenon in recent decades. MODIS LST data were used to study the daytime and nighttime temperature distribution using the urban-rural gradient and landscape-metrics-based analysis. The LUC classification showed that Greater Cairo had experienced rapid urbanization since 2000. The AREA_MN, LPI, and AI of the patches of IS (positive) and GS (negative) indicated strong correlations with mean LST at the three time points. The central area (URZ 1 ) in Greater Cairo has experienced the cool island effect in the daytime but the surface UHI effect in the nighttime. The IS density and mean LST had a week (positive) correlation in the daytime but a strong positive relationship in the nighttime. The surface UHI formation in the nighttime did not correspond with the LUC categories compared with the surface UHI formation in the daytime, indicating anthropogenic activities that strengthen the surface UHI effect.
The urban-rural gradient serves as a unique natural opportunity to identify and mitigate environmental distortion. The urban-rural gradient analysis and landscapemetrics-based analysis are beneficial for predicting the surface UHI surface and its nexus with spatial-temporal LUC changes in future surface UHIs.
We discussed the importance of increasing vegetative cover to reduce IS in Greater Cairo. Urban landscape planners must pay attention to the 1 to 11 km zones where the mean LST is high and GS density is low. It is time to mitigate the surface UHI effect in Greater Cairo by increasing green cover and dispersing the dense building distribution. In this situation, it is vital to introduce green roofs and green walls, grow clustering or irregular trees with large crowns, and construct small-and medium-sized water ponds to cope with the rising surface UHI.

Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.

Data Availability Statement:
The data that support the findings of this study are available from the corresponding author on reasonable request.

Acknowledgments:
We would like to thank Kamusoko Courage for his valuable comments and suggestions.

Conflicts of Interest:
The authors declare no conflict of interest.  Overall accuracy (%) = 87.5 Overall accuracy (%) = 81.  x FOR PEER REVIEW 24 of 28 Figure A1. Graphical illustration of the URZs used to estimate daytime/nighttime surface UHI intensity along the urban-rural gradient.

Abbreviations
, K.C. A systematic review and assessment of algorithms to detect, characterize, and monitor urban land change.