Spatio-Temporal Variability of Simulated 2 m Air Temperature for Nairobi City, Kenya

Nairobi city in Kenya is an example of the many cities in Africa experiencing rapid growth and too many environmental impacts, including the urban heat islands (UHI) phenomenon. UHI develops and escalates under outrageous hot periods, such as during heatwave, which can affect outdoor microclimate, human thermal comfort, and health and also increase the energy demand for cooling. This study used the MUKLIMO_3 model to reveal the spatial heterogeneity and temporal variability of 2 m air temperature fields and ther-mally sensitive areas within Nairobi city. The model reproduced the potential summer day conditions through idealized simulations of temperature, wind, and relative humidity based on the urban topography and local climate zones (LCZs) data at a spatial resolution of 200 m. The aim was to expand knowledge of urban climate change based on the accompanying extensive modifi-cation of land use and land cover that are critical for the local-scale atmospheric circulation. The model results revealed intricate spatiotemporal patterns of 2 m air temperature fields, which accrued from terrain-induced flows and land surface heterogeneity as described by the LCZ parameters. The variation in canopy layer UHI was a joined impact of the unfavorable location of the town and urbanization process, which added


Introduction
Urban climate change is a severe problem for cities around the world (Hashim & Hashim, 2016), particularly in developing countries like Kenya, where urbanization is happening at a fast rate (Alavipanah, Wegmann et al., 2015). Urbanization prompts changes in the landscape (Oyugi, Odenyo, & Karanja, 2017), as infrastructure replace vegetation and open land making the surfaces that were once porous and sodden to end up impermeable and dry (Ali, Patnaik, & Madguni, 2017;Hashim & Hashim, 2016;Middel et al., 2014;Ren et al., 2016;See et al., 2015;Singh, Kikon, & Verma, 2017). Tropical cities such as Akure in Nigeria, Cairo in Egypt, and Nairobi city in Kenya are currently experiencing heat waves, and they are expected to face further the challenge of increased intensity, duration, and spatial distribution of heat waves . The rapid urbanization and increased global warming (Oyugi et al., 2017), led to a sequential of environmental problems (Stone & Rodgers, 2001), especially the formation of urban heat islands (UHI) (Bechtel & Daneke, 2012;Stewart & Oke, 2012) in Nairobi city, a situation whereby urban built-up zones experience hotter temperature than their rural surroundings areas (Santamouris et al., 2017;Stone, Hess, & Frumkin, 2010;Stone & Rodgers, 2001) and which is considered to be the most interpretive indication of urban climate change in the current time of increasing urbanization (Santamouris et al., 2007;Santamouris et al., 2017). The changing temperature brought about by the UHI impacts leads to different air qualities, which relate to outdoor microclimate, cooling demands in buildings as well as human thermal comfort levels and health (Bechtel & Daneke, 2012;Ren et al., 2016;Stone et al., 2010;Vargo, Stone, Habeeb, Liu, & Russell, 2016). Its (UHI) fundamental drivers are believed to identify with the three attributes of urban areas: increased roughness due to building structural geometries, drier and progressively impenetrable surfaces and anthropogenic warmth and moisture discharges Bechtel & Daneke, 2012;Geletič et al., 2018;Stevan et al., 2013).
Climate monitoring for planning and adaptation is gradually becoming an essential element of city planning and development of Nairobi city. However, urban climate change planning and adaptation actions in a city require recognition of the possible range of UHI increase, regarding both its magnitude and spatial distribution that captures the urban footprint analysis, internal structure and E. O. Matsaba et al. DOI: 10.4236/cus.2020.82011 207 Current Urban Studies texture of the urban landscape to answer some critical questions about a rapidly urbanizing city. However, given the diversity of urban structure Bechtel & Daneke, 2012;Lelovics, Unger, Gál, & Gál, 2014) qualities of relief (Bokva et al., 2014;Bokwa, Hajto, Walawender, & Szymanowski, 2015;Saaroni & Ziv, 2010) and the variability of synoptic conditions (Przybylak et al., 2017) it is a relatively complex and challenging task since UHI studies require high range, frequency, quality and density of the measurement network and stations. On the contrary, the contemporary techniques for remote sensing focusing on land surface temperature variability may not provide critical information about 2 m air temperature (Voogt & Oke, 2003).
The spatial and temporal variation of the UHI of local cities around the world has been studied frequently (Arnds, Böhner, & Bechtel, 2017;Arnfield, 2003), as a result of the new development of climate data sources and methodological approaches in recent years, which shifted urban climatology research from identifying UHIs and prediction of UHI intensity to simulating for the exact patterns of the spatiotemporal variability of UHIs and temperature fields in urban environments. Although a large number of UHI studies are available, the results for tropical cities like Nairobi are still missing. During the current progress in exploring the complexity of the process driving the climate system, numerical models designed for application on the city scale have gradually been developed to provide more detailed and more accurate information about the spatiotemporal variability of canopy layer UHI Ching et al., 2018Ching et al., , 2011See et al., 2015). The present state of the art of numerical modeling makes it possible to solve the thermodynamics of the atmosphere and complex relations between variables, such as the buildings and vegetation types on the scale of the urban environment Bechtel & Daneke, 2012;Sun, 2013 This study focused on the MUKLIMO_3 model to simulate and identify Spatio-temporal variability of 2 m air temperature for a typical summer (warmest) day conditions in Nairobi City and its surroundings. The study applied the concept of topography and local climate zones (LCZs) to obtain representative spatial units for analyses of 2 m air temperature differences between various types of urban neighborhoods (Alexander & Mills, 2014;Bokwa et al., 2015;. The aim was to expand knowledge of urban climate change and to document the canopy UHI phenomenon for Nairobi city. MUKLIMO_3 results can be used to assess human thermal stress utilizing the perceived temperature. The study was borne out of the need to inform city planners and policymakers to manage the heat load in the city.

Description of the Study Area
Nairobi city longitude 36.6˚ and latitude −1.1˚ (upper left corner) and longitude 37.1˚ and latitude −1.5˚ (lower right corner) and altitude of 1684 m above sea level (Oyugi et al., 2017), is the largest and capital city of Kenya, with a populace of 4.40 million in 2019 census ( Figure 1). The city started as a railway town in 1899 with an initial city plan which covered 18 km 2 and later extended to 25 km 2 in 1920. The city plan completely ignored Asian laborers and Africans, which resulted in informal housing, began being built on the perimeter of the cities. The city has grown in concentric zones, clearly-defined historical centres, industrial areas, old and new residential buildings, housing estates, modern shopping centres, and malls with allotments. With functionalism as the central principle, Nairobi was classified into zones that led to a neighborhood with irregular, as well as regular street plans. There is a presence of an exceptionally polluted stream (the Nairobi River) and a greater extent of open low plant vegetated spaces. The natural landscape outside the city centre is prevalently natural, with agriculture areas on the Western and sand and uncovered soils on the Eastern part. Patches of LCZs A-dense trees B-scattered trees, and C-bush scrub complete the scene mosaics of the city (Ochola et al., 2020).

Simulation of 2 m Air Temperature Using MUKLIMO_3 Model
The thermodynamic version of the microscale urban climate MUKLIMO_3 model included prognostic equations for atmospheric temperature, relative humidity, wind speed, and wind direction (Früh et al., 2011). MUKLIMO_3 model was employed on the urban scale to simulate the 2 m air temperature field, relative humidity field, and wind field for Nairobi city and its surroundings on a three-dimensional (3D) model The 2 m air temperature field for the potential warmest day (26th February 2018) condition was reproduced through idealized simulations of temperature, wind, and relative humidity in Nairobi city based on the city's topography and LCZs data (Ochola et al., 2020) at a spatial resolution of 200 m. The model used LCZ data (Land use file) to describe the urban land use classes and unique thermal surfaces within Nairobi city ( Figure 3).
Topographical data (elevation file) was used to describe the influence of city  The five layers of meteorological data (air temperature, relative humidity, and speed of wind) in the vertical profile of the atmosphere (up to 3142 m above ground level) for the operation of the MUKLIMO_3 model was derived from the sounding data of 63741 HKNC Nairobi (Nairobi Dagoreti weather station located at latitude −1.300587˚ and longitude 36.759368˚) observations at 00Z 28 February 2018 (Table 1).  , which solves the Reynolds-averaged Navier-Stokes equations to simulate atmospheric flow fields in the presence of unresolved buildings (Früh et al., 2011). The model considered the friction impacts of building structures and turbulence generated because of wind stream separation to calculate turbulent fluxes of heat, moisture and momentum (Hollosi & Zuvela-aloise, 2017;Žuvela-Aloise, 2017;Žuvela-Aloise, Koch, Buchholz, & Früh, 2016) at ground level, the walls and the roofs of buildings in a domain with unresolved build-up, is described in Früh et al., 2011(Geleti, Lehnert, Dobrovolný, & Bulletin, 2016. The model also considered the effects of cloud cover on radiation. However, it did not include cloud processes, precipitation, horizontal ru-

Results
The model simulations were represented with 24 temperature maps at a horizontal resolution of 200 m. The time step between two successive modeled fields was 60 minutes. The results indicated a complex spatial structure of the simulated 2 m air temperature and urban heat load in Nairobi City. The minimum day time simulated 2 m air temperature, and maximum nighttime simulated 2 m2 m air temperature was found in densely built-up LCZs in the city center and residential areas with flat terrain north-east of Nairobi city. Lower 2 m air temperature early in the morning was typical of Natural LCZ types, while the warmer 2 m air temperature was typical of built-up LCZ types. The places with the lowest air temperature in the early morning hours before sunrise was located in areas dominated with the featureless landscape of soil or sand cover with few or no plants (LCZ F) in the southeastern (Mlolongo, Syiokimau, Athi River and some parts of Emabakasi) parts of Nairobi ( Figure 5: 2 m air temperature for Nairobi at 06 am to 08 am). After the sunrise, lower 2 m air temperature is typical to natural landscapes with forests, scattered trees, bushes, low plants, and water bodies (LCZs A, B, and G).
In comparison, warmer 2 m air temperature are typical to both open and high-density built-up areas very large low rise and heavy industrial buildings (LCZs 1 -10), and open natural with bushes and scrubs, low plants (mostly grasslands, agricultural fields, and parks), and paved allotment areas, and bare soils and sand (LCZs C, D, E, and F) ( Figure 5 and Figure 6: 2 m air temperature for Nairobi from 09 am to 05 pm). As the air gets warmer (from 11 am to 05 pm), the model gradually generates areas with a higher proportion of LCZ 2, 3, 7, 8, 9, and 10 that are warmer than their surroundings, including areas located outside the compact urban development. However, the distinct UHI is not formed until noon to 03 pm. At 03 pm, the UHI is formed over most of the city center.
After the sunset, the distribution of the 2 m air temperature was predominantly influenced by the spatial distribution of wind fields within LCZs. Areas with lower temperature were located in high-density built-up areas very large low rise and heavy industrial buildings (LCZs 1 -10) ( Figure 6 and Figure 7: 2 m air temperature for Nairobi at 06 pm to 05 am). In comparison, the warmer air temperature was located in build-up LCZ types. As the air gets colder (at 06 pm to 05 pm), the model gradually generated areas with a lower proportion of hot spots ( build-up areas) that are warmer than their surroundings, including areas located outside the compact urban development. The distinct nocturnal UHI formed only in built-up areas of the city center and reduced gradually after the sunset. Hotspots within the city were formed in areas with a higher proportion of high-density built-up classes (LCZ 2, 3, 7 and 10), with 2 m air temperature range of 30˚C to 40˚C during daytime and 28˚C to 34˚C at night. On the contrary, relatively cooler spots within the city correspond to larger areas of LCZ  up to 15˚C between the warmest part of the city and the most cooling forests at 1200 noon to 03 pm ( Figure 5 and Figure 6). The model, however, forecasts a relatively higher intensity of canopy layer UHI during the evening and night hours. At 9 p.m., the warmest parts of the city are only about 12˚C warmer than the natural landscape around the city (predominantly LCZ D) and up to 10˚C warmer than the forested areas.
In summary, the MUKLIMO_3 model output revealed intricate spatial and temporal patterns of 2 m air temperature fields across Nairobi city and its environment. High-density built-up areas of the city recorded possible maximum urban cooling during the daytime and maximum urban heating patterns at night compared to areas that were sparsely built or covered with vegetations. The spatial pattern of simulated 2 m air temperature confirmed the presence of canopy layer UHI across Nairobi city. The variation in simulated warmer 2 m air temperature in parts of the town accrued from terrain-induced flows and land surface heterogeneity as described by the LCZ parameters.

Discussion
February 26, 2018, was considered in this study because it was a typical summer day having recorded the hottest temperature with maximum thermal stress within the city. The model results confirmed the presence of hotspots in the city both daytime and night time, especially in high-density built-up areas of the city with a high percentage of water-resistant, non-reflective surfaces and low vegetation compared to the surrounding rural areas with urban green spaces. Spatial analysis of simulated of 2 m temperature fields and the increase of heat load in the high-density built-up areas demonstrated canopy layer UHI in Nairobi city as a joined impact of both terrain (for example, terrain sheltering and the unfavorable location of the city), or local land surface heating (due to a high percentage of water-resistant, non-reflective surfaces, building structure type and density, a high fraction of asphalt and a small amount of vegetation cover) which add to the formation of excessive urban heat load.
During the day, the heat load in high altitude vegetated areas was locally reduced compared to the lowland high built-up areas, which are visible in the vertical profile of temperature (Table 1) due to the temperature inversion, which makes the elevated and natural areas warmer. The location of the city and the building type, consisting of high density, a high fraction of pavement, and a small amount of vegetation, contributed to the formation of excessive heat load at night. During the night, robust and reliable cooling of the free surface results in a temperature inversion, and therefore MUKLIMO_3 model assumed a temperature that was just 1˚C higher in the high-density built-up areas than in the rural natural surroundings and temperature that was about 10˚C higher in the built-up areas than in the coldest forests.
MUKLIMO_3 model primarily reflected the effect of elevation and land cover on UHI distribution evident during the daytime and at night. The model pre- Matsaba et al. dicts globally higher temperature at lower elevations. It can be considered as a simplification (Bokva et al., 2014;Bokwa et al., 2015). On the other hand, the model did not reflect the extent of the variability of building density (i.e., the amount and effect of accumulated heat). It is anticipated that for more accurate simulation of the spatiotemporal temperature field, it is necessary to focus attention on the LCZ table.
Moreover, it is possible that the concept of LCZ is too general for modeling on a detailed level and may have caused incorrect settings of the thermal capacity of individual surfaces. However, the use of LCZ classes to describe urban features decreased model uncertainties related to the choice of selection by limiting and standardizing the parameterization of the unresolved buildings in the MUKLIMO_3 model. Generally, the MUKLIMO_3 model considered not only the physical features of the unresolved urban structure but also the elevation and advection of the air masses between the surrounding surfaces computing atmospheric processes in the urban condition. The heat load for an entire individual city relied upon, relief, the surrounding land use types, and location in the city.
The comparison of the real station measurements and MUKLIMO_3 simulations in Nairobi city and its surroundings was not possible because of a lack of access to the hourly 2 m air temperature data from the metrological stations. However, the model prediction was consistent and, therefore, could correspond best with the situation on 28 February 2018 (see Figures 5-7); this may be related to more stable atmospheric conditions very close to climatological warm days. However, Individual problems with the accuracy of the model simulation could be primarily related to the starting phase of the modeling (at 10.00 am, and 11.00 am; see Figure 5). Additionally, the MUKLIMO_3 model should be used to simulate future projections of the UHI phenomenon and create spatial configurations of planting designs with more cooling and fewer energy demands and their application to the whole in Nairobi city.

Conclusion
MUKLIMO_3 model was used to reproduce a simplified urban structure and indicate both spatial and temporal variability in 2 m air temperature for a potential summer day in Nairobi city. Simulation results showed complex spatial and temporal patterns of 2 m air temperature fields, which confirmed the presence of canopy layer UHI. The variations in simulated 2 m air temperature patterns accrued from terrain-induced flows and land surface heterogeneity, as described by the urban parameters. Therefore, the enhancement of urban heat load in Nairobi city could be linked to the concept of rapid urbanization process and its historical development. Hence urban planning strategies such as optimization of mixed building heights and the introduction of green & blue infrastructure were critical to mitigating heat-stress across Nairobi city. The results were the first to provide the relationship between LCZs and the spatial distribution of UHI for Nairobi city using the MUKLIMO_3 model. The model outputs could be used to study the development of urban microclimates within Nairobi city. The results could also be useful for evaluating the influence of location, urban parameters, and weather conditions on urban microclimatic within Nairobi city. However, the MUKLIMO_3 model must be validated for Nairobi city to reach a better performance. Therefore, there will be further studies to understand the model settings and optimal input data preparations.