Utilizing Remotely Sensed Observations to Estimate the Urban Heat Island Effect at a Local Scale: Case Study of a University Campus

: The urban heat island (UHI) effect has become a significant focus of research in today’s era of climate change, and a key consideration for the next generation of urban planning focused on green and livable cities. UHI has traditionally been measured using in situ data and ground-based measurements. However, with the increased availability of satellite-based thermal observations of the Earth, remotely sensed observations are increasingly being utilized to estimate surface urban heat island (SUHI), using land surface temperature (LST) as a critical indicator, due to its spatial coverage. In this study, we estimated LST based on Landsat-8 observations to demonstrate the relationship between LST and the characteristics of the land use and land cover on the campus of King Abdulaziz University (KAU), Jeddah, Saudi Arabia. We found a consistent variation of between 7 and 9 degrees Celsius for LST across campus, spanning all summer and winter seasons between 2014 and 2019. The LST correlates strongly with both green vegetation and built-up land cover, with a slightly stronger correlation with the latter. The relationship between LST and green vegetation has a notable seasonality, with higher correlation in the summer seasons compared to the winter seasons. Our study also found an overall increase in LST between 2014 and 2019, due to intentional changes in the built-up land cover, for example from the conversion of natural green surfaces to artificial surfaces. The findings of this study highlight the utility of the remotely sensed observation of LST to assess the SUHI phenomenon and can be used to inform future planning aimed at securing green and livable urban areas in the face of a changing climate. show a significant and positive correlation between LST, NDBI, and UI (up to r = 0.69 and r = 0.68, respectively, p < 0.001 for both), in parallel to a significant negative correlation between NDVI and LST (up to r = −0.43, p < 0.001). These results suggest that LST is negatively correlated with green vegetation and positively correlated with built-up land cover. Interestingly, while the correlation between LST, NDBI, and UI does not vary much by season, the correlation between LST and NDVI


Introduction
Land surface temperature (LST), which is also referred to as the 'surface radiometric temperature' [1] or 'radiative "skin" temperature' [2], can be measured remotely, for example by satellite or airborne remote sensing. LST can be used to estimate radiation budgets and heat balance, and so can be used as a control for global climatological models [3]. LST is also a key factor when estimating the urban heat island (UHI) effect, i.e., when an urban area exhibits higher atmospheric and surface temperatures compared to adjacent rural areas [4], with a significant impact on the ecology, environment, and society [5]. In general, the UHI effect can be categorized into three main types: boundary layer urban heat island (BLUHI), canopy layer urban heat island (CLUHI) and surface urban heat island (SUHI) [6]. In this study, we examine the SUHI, as measured on the surface. Ensuring that cities are green and livable in the face of a changing climate requires careful consideration of the implications of SUHIs in urban planning [7]. Accurate and timely information regarding spatial and temporal variations in LST are essential for climatological, ecological, hydrological, and biogeochemical applications [8].
Patterns of land use and land cover (LULC) on Earth are essential factors affecting LST. Different LULC patterns are distinguished by unique reflectance and roughness, thereby leading to variations in LST [9,10]. Other factors, such as seasonality, climatological conditions [11], and the spatial clustering of the land cover [12], including its composition, configuration, complexity, intensity, morphology, and heterogeneity [3,[13][14][15][16] are also key variables affecting LST.
Because LST is impacted by the characteristics of the land [17], LULC changes within urban areas and their fringes have a direct impact on LST [18,19]. In general, urban artificial impervious surfaces (roads, industrial areas, and other buildings) are associated with higher LST [20,21], compared to water bodies and vegetation land cover, which mitigate the microclimate conditions and lower the temperature in urban areas [22]. As a result, an increase in the area of asphalt roads and commercial and industrial land use, along with a decrease in the area of urban parks and gardens, other green spaces, and agricultural land is often associated with increased surface temperature [23], especially during the daytime [21]. Other factors affecting LST include, for example, building materials [10], the density of built-up structures, building heights, and population density [24]. Although the relationship between LST and the distribution of built-up land cover is not always linear and can vary by geography and seasonality [25], previous studies have found an association between LST and vegetation or built-up land cover, as measured utilizing remotely sensed observations; for example, using the normalized difference vegetation index (NDVI) and the normalized difference built-up index (NDBI) [26][27][28][29][30].
While previous studies have evaluated SUHI on the scale of a city or neighborhood, very few have measured and examined the SUHI effect on a local or micro scale by means of satellite observations. Many of the studies that examine micro-scale variations of UHI have used either in situ air measurements or airborne measurements at the surface level (e.g. [31,32]), with less emphasis on satellite-based observations.
Various micro-scale factors can potentially reduce the SUHI effect for local climate mitigation. For example, according to [33], SUHI can be eliminated when wind speed reaches a certain velocity, based on a "city-size-dependent threshold." Similarly, an open low-rise gridiron precinct can potentially impact ventilation performance and, as a result, outdoor thermal comfort. According to [34], precinct ventilation performance has a clear positive impact in diminishing the SUHI effect by reducing humidity. When considered and implemented in urban planning, natural ventilation can often reduce the air temperature and enhance the subjective comfort experience [35,36]. A local climate zone (LCZ) relates to a classification system that aims to unify heat island measurement; LCZ was one of the first climate-based classifications for heat island studies [37]. The classification of the zones is not related to time, culture, or location, but it is quantifiable, based on land cover, building morphology, and thermal climate.
In light of the scarcity of studies that evaluate the relationship between LST and SUHI on a local scale, the main objective of this study was to examine the SUHI effect on a local scale, i.e., a university campus in Saudi Arabia. The study derived LST from publicly available Landsat-8 observations in the period 2013-2019, to identify and explain the relationship between LST and LULC. The analysis was implemented using the Google Earth Engine (GEE) (we provide a description about GEE in Section 2). This paper contributes to the existing body of literature by demonstrating the potential of Landsat-based observations to be used to estimate temporal trends in LST (annual and seasonal) at a local scale, and the relationship between these trends and changes in LULC characteristics. The study aimed to link LST and LULC to micro-scale planning decisions and show that Landsat-based measurements can be used to capture local and micro-scale changes in LST patterns.

Materials and Methods
In this study, we analyzed the relationship between annual and seasonal variations in LST and changes in LULC in the King Abdulaziz University campus (KAU), between 2013 and 2019. We estimated and derived LST based on Landsat-8 observations (Tier 1 raw scenes). Furthermore, we estimate the characteristics of the land cover based on spectral indices. This is because the spectral indices prove to be an effective of mapping land cover [38] and provide better classification results [39]. In addition, it was derived from Landsat-8, as a proxy for vegetation land cover (NDVI), and for soil and built-up land cover ((NDBI) and Urban Index (UI)). Currently, there are two operating Landsat satellites, Landsat-7 (launched in 1999; carries the enhanced thematic mapper (ETM+) sensor) and Landsat-8 (launched in 2013; carries the operational land imager and TIRS camera). While both satellites can be used to derive LST, Landsat-8 has narrower spectral bands compared to Landsat 7, and TIRS provides additional capabilities for the estimation of LST compared to previous versions of Landsat. While previous studies have shown the potential of estimating LST from Landsat 7 data, some methods, such as the single-channel method which we used in this study, rely on Landsat-8's B10 or B11 thermal bands.
We performed the analysis in Google Earth Engine (GEE), which is a platform that is frequently used for planetary-scale analysis in a wide range of applications [40], including mapping populations [41,42] and urban areas [43], as well as monitoring LULC changes and their relationship to SUHI [44,45].

Study Area
We evaluated the relationship between LST and land cover changes in KAU. Established in 1967, KAU occupies an area of 2224 acres (9 km²) in Jeddah, Saudi Arabia (between 21°29′N and 21°30′N latitude and 39°14′E and 39°16′E longitude). When KAU was established, the university buildings were constructed, primarily from natural stone, bricks, and concrete. The university campus has since grown, and new buildings have been constructed to expand faculties and accommodate more students. The newer buildings on campus often use granitic rock for decoration of the walls and floors [46]. The campus area includes faculty buildings, classrooms, the university hospital, sports facilities and the equestrian club, labs, central library, and administrative buildings and centers, as well as staff and female student housing. The campus has five public open spaces in the female campus and three in the male campus ( Figure 1). These spaces account for 0.62% of the campus area. Furthermore, due to limited public transportation to serve the staff and students, there are many areas designated for car parking, which accounts for 10.42% of the campus area. (Figure 2 presents air photos of selected sites on campus. A detailed description of these sites and their characteristics is provided in Appendix A (Table A1)).

Estimating LST
As with the methodology proposed by [47], we estimated LST based on Landsat-8 observations using a single-channel method, which has previously been shown to provide relatively accurate LST retrieval [48] and has been widely used due to its simplicity [49,50].
We relied on Band 10 (thermal infrared 1, resampled from 100m to 30m, 10.60-11.19 µm). Previous studies suggest that stray light artifacts due to out-of-field radiance that scatters onto the detectors are often found in TIRS data, adding a signal across the field-of-view that is non-uniform, and approximately twice as large in Band 11 as it is in Band 10 [48]. We estimated the average LST during the summer (21 June to 22 Sept) and the winter (22 Dec to 20 March) of each year between 2013 to 2019 (totaling 14-time intervals).
The procedure involves the following steps: Identify Landsat scenes: We obtained Landsat 8 Collection 1 Tier 1 DN values, representing scaled, calibrated at-sensor radiance, which are available in GEE. Tier 1 includes Level-1 Precision Terrain (L1TP) processed data with well-defined radiometry, inter-calibrated across different Landsat sensors. We retrieved all Landsat-8 scenes that were captured during each time interval (we selected only scenes with less than 5% cloud coverage). TOA spectral radiance: We converted the raw reflectance value in Band 10 (B10) of each scene into a TOA spectral radiance (radiance measured by the sensor). Two parameters are required to convert the DN value into TOA: the radiance multiplicative and the radiance add, which we extracted from the metadata of each scene using this expression: where TOA (L) is the top of atmosphere radiance of B10, ML is the multiplicative rescaling factor of B10 in a given image, DN is the reflectance value of B10, and AL is the additive rescaling factor of B10 in a given image. LST can be determined using the following expression: where ψ1, ψ2, and ψ3 can be derived as a function of total atmospheric water vapor content (w) using the following matrix approximation: where w is the water vapor, and coefficients cij are derived from [47] as follows: We retrieved water vapor content information from NCEP/NCAR. The NCEP/NCAR Reanalysis Project is a joint initiative between the National Centers for Environmental Prediction (NCEP, formerly NMC) and the National Center for Atmospheric Research (NCAR). The goal of the project is to produce atmospheric analyses based on historical data, along with analyses of the current atmospheric state (climate data assimilation system, CDAS) (also available in GEE). We retrieved surface water vapor information for the corresponding date of each Landsat scene (four observations per day) and calculated the average water vapor over the four observations per day. We calculated emissivity (ε) according to the following expression [47]: Where fractional vegetation cover (FVC) is calculated as: Moreover, is the surface emissivity of vegetated ( ) and soil ( ) surfaces. Like [51], we assumed soil and vegetation emissivities of 0.97 and 0.99, respectively.
Finally, we converted the LST values in Kelvin to degrees Celsius. We estimated the LST for each pixel in each Landsat scene over each interval (season) and calculated per pixel, per season LST values as the average value of all overlapping pixels in an interval. As explained above, the main objective of this study was to estimate seasonal variations in LST and their relationship to LULC changes. In this respect, it should be mentioned that, while LST, by its nature, may fluctuate due to microclimate conditions, the characteristics of LULC do not tend to change much monthly. While some previous studies have examined the relationship between temperature and LULC by using single daily images [52], other studies have captured long-term variations using monthly, seasonal [53] or annual composites [54]. In this research, we examined seasonal and annual variations in LST, rather than daily or bi-weekly trends. This method aims to capture broad seasonal trends in LST and LULC (estimated as the average across an entire season), rather than to examine daily climatological trends that are less related to LULC patterns.

LULC Characteristics
To establish the relationship between LST and the characteristics of land cover, we relied on several remote-sensor-derived spectral indices, which we obtained from Landsat-8 Collection 1 Tier 1 calibrated TOA reflectance composites (available in GEE). NDVI: NDVI [55] conveys the relationship between visible red light (which is absorbed by chlorophyll) and near-infrared wavelength (which is scattered by leaf mesophyll cells). It is calculated as: where NIR is the near-infrared wavelength, and the RED is the red wavelength (Bands 5 and 4, respectively, in Landsat 8). NDVI ranges from (−1) to (+1). An NDVI value was calculated for each pixel per season (average value for all overlapping pixels in scenes with less the 25% cloud coverage during a season). NDBI: NDBI [56] conveys the relationship between medium infrared and near-infrared wavelengths. It is calculated as: Where SWIR1 is the reflectance in the shortwave infrared wavelength (1.57-1.65 µm), and NIR is the reflectance in the near-infrared wavelength (Bands 6 and 5, respectively, in Landsat 8). NDBI ranges from (−1) to (+1). An NDBI value was calculated for each pixel per season (average value for all overlapping pixels in all scenes with less the 25% cloud coverage during a season).
UI: UI [57] normalizes the NIR and SWIR 2 bands (Bands 7 and 5 in Landsat 8) and relies on the inverse relationship between the brightness of the NIR and SWIR in built-up areas. Previous studies have used UI to measure urbanization and urban sprawl [58] and to extract the extent of built-up land cover from remotely sensed observations [59]. The index is calculated as follows: A UI value is calculated for each pixel per season (average value for all overlapping pixels in scenes with less than 25% cloud coverage during a season).

Aggregation to the Hexagon Level and Analysis
 Aggregation to the hexagon unit of analysis (hexagonal tessellation) For this study, we used a hexagon tessellation mapping procedure to estimate the relationship between patterns of LST and LULC. Aggregating pixel-level data to a given unit of analysis is necessary when analyzing LULC characteristics at different geographical scales. The approach of relying on hexagonal tessellation for LULC research, which can also be referred to as 'zone mapping' [60], has previously been used in the remote sensing domain [61,62]. It allows the creation of areas that are relatively homogenous in terms of their land cover and land use spectral reflectance, as well as their LST. Here, we created a grid of cells (hexagonal tessellation) of 2546 hexagons that span the campus; each hexagon is around 3100 m 2 (compared to an approximate area of a Landsat pixel being 900 m 2 , i.e., around 3-4 Landsat pixels overlapping with one hexagon). For each hexagon, we calculated the average LST, NDVI, NDBI, and UI values of each overlapping pixel.
 Correlation and multiple regression analysis We relied on the hexagons in the hexagonal tessellation as the units of analysis. We evaluated the correlation (Spearman correlation test) between LST and each of the indices. In addition, we calculated the per-hexagon slope of change (SOC) of LST, NDVI, NDBI, and UI values over the period 2014 to 2019 (we define SOC as the rate of change or slope of a linear function between LST and the spectral indices). Furthermore, we performed a multiple regression analysis to estimate the potential of the explanatory variables to explain, together, variations in LST. Multiple regression analysis has previously been used in the field of remote sensing research; for example, by [63], who estimated the potential of 28 remotely sensed derived variables to explain variations in air temperature, or by [64], who performed a thermal unmixing approach, calculated using multiple linear regression to evaluate the relation between LULC patterns and LST.
 Getis-Ord Gi* for hot spot analysis We also evaluated changes in patterns of LST and patterns of LULC using the Getis-Ord Gi*method for hot spot analysis [65]. This method is ideal for identifying local patterns in the spatial distribution of data, which could not have been revealed using global statistical measures [65]. In other words, because global statistics conducted over an entire area of interest are often too coarse and are insensitive to local "nuances," the Getis-Ord Gi* method was used to locate and identify clusters of concentration of high or low values [66]. Given a set of weighted features, Getis-Ord Gi* statistics identify statistically hot and cold spots. Here, we calculated the Getis-Ord Gi* statistic for each hexagon in the dataset, which allowed us to identify hexagons that exhibit high and low values that are spatially clustered. The Getis-Ord Gi* statistic is calculated as [67]: is an attribute value of feature j, , represents a spatial weight between feature i and j, while n is the number of features, and:

Results
This section is divided into four sub-sections. First, we present the spatial and temporal patterns of LST. Then, we present a detailed analysis of the relationship between spatial and temporal patterns of LST and LULC, along with the results of correlation and regression tests to estimate the relationship between changes in LST and changes in LULC. Finally, we present the results of the Getis-Ord Gi* hot spot analysis, to identify temporal changes in hot and cold spots of LST and LULC.

Spatial and Temporal Patterns of LST
The analysis shows clear spatial and temporal variations in the distribution of the LST across the KAU campus. As can be seen in Figure 3 (seasonal trends in average LST in the period between 2013 and 2019), the average remotely sensed, derived seasonal LST ranged between 45.96 °C (Std = 2.57 °C) in the summer, and 35.39 °C (Std = 1.53 °C) in the winter. The highest summer temperature was recorded in 2014, and the highest winter temperature was recorded in 2018. In addition, we find a clear variation in the spatial distribution of LST across campus. Figure 4 presents, as an illustration, the distribution of LST in the summer and winter of 2019. In summer 2019, the LST ranged between 48.5 °C and 40.4 °C. In the winter of 2019, LST ranged between 37.87 °C and 28.7 °C. The results suggest some similarities between summer and winter in the distribution of the LST.
In general, we find that open spaces that are not covered with green vegetation (typically, bare land or other types of surface, such as asphalt and compressed soil) are characterized by a relatively high LST. This trend is observed, for example, over the open areas in the northern part of campus (area number 1, Figure 2), and east to the staff housing (area number 2, Figure 2), as well as over several paved parking spaces (e.g., areas number 3, Figure 2). For example, in summer 2019, the temperature over the main campus parking lot (area number 3, Figure 2) was 6.5 °C higher than over the green space south of the parking lot (area number 4, Figure 2) (i.e., 47.18 °C compared to 40.72 °C, respectively) (the distance between these areas is approximately 250 meters). Due to a lack of public transportation within the campus and between the campus and Jeddah and surrounding cities, there are many parking spaces (e.g., areas number 3 and 5, Figure 2), most of which are not proximal to trees or shading structures. The LST over these parking spaces is up to 6 °C higher, compared to the temperature over adjacent green spaces across campus. Some of the cooler areas on campus include the space near the sports tent (area number 4, Figure 2), the park in the central-eastern part of the campus (which is characterized by mature trees and various planting types) (area number 3, Figure 2), and even the shaded parking space west of the staff housing (area number 5, Figure 2). These trends are also evident in the winter. In winter 2019, the temperature over the parking spaces across campus was around 34.5 °C (areas number 3 and 5, Figure 2), compared to the temperature over the green spaces, which was approximately 30.7 °C (areas number 6 and 4, Figure 2). Study area number 7, which is new staff housing at the University, shows a consistently low LST, both in winter and in summer. A potential explanation for this result could be the white-painted roof of staff housing, which is a rarity on the KAU campus and the orientation of the buildings. These buildings are structures on a 45°-degree angle to the direction of the wind, which enhances natural ventilation. In addition, the form of the buildings and the staggered layout help in accentuating the wind movement between blocks. Hence, the orientation of the buildings, which allows for wind movement and increased velocity, together with the painted roof, may explain, at least in part, the lower LST observed over Area 7 (staff housing).

The Relation Between Patterns of LST and LULC
The spatial variations in LST can be explained, at least partly, by the distribution of green vegetation across the campus. To identify the relationship between the distribution of green vegetation and built-up land cover with LST, we overlaid hexagons with the highest NDVI and NDBI values (i.e., hexagons with an average NDVI and NDBI values that exceeded a set threshold, which we defined as the 80th percentile of the NDVI and NDBI value of all hexagons in the dataset). Figure  5 presents, as an illustration, the hexagons that exhibit the highest NDVI and NDBI values in the winter and summer of 2019, against LST (these hexagons are characterized by an NDVI value that is greater than 0.115 and 0.136 and with an NDBI value greater than 0.049 and 0.043 for the winter and summer, respectively). As clearly illustrated in Figure 5, hexagons that exhibit the highest NDVI value are also characterized by relatively low temperatures, while hexagons that exhibit high NDBI are characterized by higher temperatures. This trend is observed in the winter and the summer. In other words, higher LST corresponds with less green vegetation and more built-up land cover.

Getis-Ord Gi* Hot Spot Analysis
To identify temporal trends in LST, we performed hot spot analysis using Getis-Ord Gi* [65], which allows the estimation of hot and cold spots by examining features and their neighbors. For example, a feature (hexagon) is characterized as a statistically hot spot if both it and other surrounding features (hexagons) have high values [25]. Because the correlation and regression tests are performed at the geographical unit of the hexagon, we also perform the Getis-Ord Gi* at the hexagonal unit. This allows us to estimate the extent to which the observed changes in the correlations between the variables are also reflected in their spatial patterns. We calculated the average temperature in the summers and winters of 2014-2015, 2016-2017, and 2018-2019, and identified significantly hot and cold spots based on Getis-Ord Gi* analysis (we choose to calculate the average over two-year intervals to identify changes in the patterns over the time visually). We found changes in the spatial patterns of LST across the campus, which we related to changes in LULC. For example, the campus stadium (area number 8, Figure 2) changed from a cold spot (in 2014-2015 and 2016-2017) to a hot spot (in 2018-2019). We linked this to changes in the stadium's field ground from natural to artificial turf (in 2018, all sports grounds on campus were changed to artificial turf to reduce maintenance and irrigation costs), which has probably resulted in a significant increase in the LST. Another example can be found in the southwestern part of the campus (the driving school, area number 13, Figure 2 To evaluate the relationship between changes in LST and changes in LULC, we calculated hot and cold spots for built-up land cover (NDBI and UI) and green vegetation (NDVI) (we used Getis-Ord Gi*). As can be seen in Figure 6 (hot and cold spots of live vegetation), several changes can be observed across time in the patterns of green vegetation (i.e., areas that exhibit significantly high NDVI values). In 2018-2019, for example, there were some new hot spots of green vegetation along the main street and entrance to the Higher Management Building (area number 10, Figure 2), which are likely to be the result of incorporating new green vegetation and natural turf. Several hot spots of built-up land cover (measured according to NDBI) can be seen, especially along the eastern and southern boundaries of the campus. A significant hot spot of built-up land cover can be seen around the academic square, or the new male campus, which is located at the center area of the campus (area number 11, Figure 2). This area incorporates large faculty and administration buildings and has little green vegetation. This area is also a hot spot with high LST (Figure 7). Changes in the observed spatial patterns of built-up areas (according to NDBI) in the period between 2014 and 2019 can be seen in the area surrounding the staff housing (area number 12, Figure 2), which was a hot spot of high NDBI in 2014-2015 (the construction of this area was completed at the end of 2018) (Figure 8). In parallel with the change of this area from a hot spot of NDBI to a non-hot spot, this area became a cold spot for LST. Additionally, in 2014-2015, the area in the southern part of the campus (area number 13, Figure  2) used to be a cold spot of NDBI. This area was constructed in 2018. As a result, this area showed an increase in LST between 2014-2015 and 2018-2019, causing it to change gradually from a cold spot of LST to a non-cold spot.

LST and LULC Correlation and Regression Tests
To further evaluate the relationship between LST and the patterns of LULC, we performed a correlation test to estimate the correlation between green vegetation (NDVI) and built-up land cover (NDBI and UI) and LST. We provide, in Appendix A ( Figure A1), as an illustration, the distributions (histograms) of LST, NDVI and NDBI values across all the hexagons in the study area. As can be observed, the values are not normally distributed and are skewed either to the right (NDVI) or left (NDBI). A Shapiro-Wilk normality test confirmed that the values are not normally distributed ( Table  1 presents the results of the test, indicating a p-value of 0, i.e., the test rejects the hypothesis of a normal distribution). While the Pearson Correlation test measures the linear relationship between two variables, the Spearman rank-order test measures the rank order of two variables and is a nonparametric version of the Pearson test. Spearman's test may be more appropriate than the Pearson test, especially when the variables are characterized by a high kurtosis and can be used to determine if there is an association between two nonparametric variables. Here, we performed a Spearman's rank-difference coefficient test to estimate the relationship (correlation) between the variables. does vary by season, with a stronger (negative) correlation in the summer compared to the winter (indicated by a correlation of up to r = −0.42 in the summer and up to r = −0.33 in the winter; p < 0.001 for both). These results imply that, first, LST is more impacted by the distribution of built-up land cover on campus than by the distribution of vegetation, and second, that the effect of vegetation on moderating the temperature is more significant in the summer than in the winter. Figure 9. Spearman correlation test between LST and NDVI, NDBI, and UI. P < 0.01 for all.
Next, to evaluate the explanatory power of the spectral indices to predict LST, we performed a multiple regression analysis. The results (presented in Table 2) show that, together, NDVI, NDBI, and UI explain up to 49% of the variation in LST (in winter, 2014). In all seasons examined, we found that the explanatory power of NDVI is relatively low (up to R 2 = 24% in summer, 2013), and that it is lower in the summer than in the winter. We also found that the addition of NDBI and UI as additional explanatory variables improves the prediction, in which there is much less variation between winter and summer.
Finally, to evaluate the relationship between temporal changes in LST and LULC (namely, green vegetation and built-up land cover), we performed a Spearman correlation test to (Table 3)  Correlation Coefficients (r) NDVI NDBI UI

Discussion
The SUHI effect, i.e., when an urban area exhibits a higher atmospheric and surface temperature than adjacent rural areas [4], has been the subject of a large body of studies. Many factors affect the extent and nature of SUHI, including LULC characteristics [68] and their spatial patterns and configuration [17]. Changes in the characteristics of LULC, including the extent of impervious surfaces and built-up structures, have implications for the local and regional-scale climate, including on the SUHI [17,25,69,70]. Because green spaces and vegetation can mitigate the impact of the SUHI [71] effect and decrease the temperature [72,73], urban university campuses, which often have large areas of vegetation and green space, play a valuable role in mitigating the SUHI effect and helping to reduce its negative impact [74].
Until recently, the SUHI studies mainly use three type of data: in situ temperature observation, remote sensing data and climate models simulation [75]. However, with the ever-growing availability, accuracy, and sensitivity of satellite observations, an increasing number of studies are utilizing thermal remotely sensed data to measure LST as a proxy for SUHI [76,77].
This study utilized the LST measurement derived from remotely sensed observations to understand the SUHI characteristics at a local microscale, and the relation with LULC patterns, which is based on many studies, becomes highly essential in SUHI studies [78].
In this study, we showed the potential utility of LST measurements derived from remotely sensed observations to understand the SUHI characteristics at a local micro-scale, and the relationship with patterns of LULC. We relied on Landsat-8 observations to estimate LST and its relationship with the characteristics and spatial patterns of LULC at the scale of a university campus (KAU in Saudi Arabia).
We found that the distribution of hot and cold spots of LST is closely linked to the characteristics of the land cover on campus, namely, vegetation and built-up land cover. Green spaces (i.e., areas characterized by high NDVI values) exhibit significantly lower temperatures compared to their surroundings, while areas characterized by built-up land cover (i.e., high NDBI and UI values) exhibit higher temperatures relative to their surroundings. The characteristics of the LULC on campus have a significant impact on the variation of LST. As with previous studies that show that green urban parks can be cooler by 5-7 °C relative to their surroundings [79] and other types of built-up land cover [80], we found that the LST varies greatly between green and artificial surfaces. For example, the estimated temperature over the main campus parking lot (in Summer 2019) was 6.5 °C higher than over the adjacent green space.
The relationship we found between LST and the LULC characteristics are also indicated by a significantly negative correlation between LST and NDVI (r = −0.49) and a significantly positive correlation between LST and NDBI and UI (up to r = 0.66 and r = 0.64, respectively). Thus, in agreement with the findings of [29], these results suggest a weaker correlation between LST and NDVI compared to the correlations between LST, NDBI, and UI, indicating that the built-up indices (NDBI and UI) may be better predictors/descriptors of LST than NDVI. Furthermore, the strong positive correlation between LST, NDBI, and UI demonstrates the potential impact of rapid urbanization processes on the SUHI effect and illustrates how the built environment alters the temperature, even at a local scale. For example, policies related to the conversion of natural green vegetation to artificial turf may-as observed, for example, in the case of the University sports stadium-have immense implications for the microclimate by increasing LST significantly.
We report that while the correlation between LST, NDBI, and UI does not vary significantly between seasons, the correlation between LST and NDVI does vary by season, with a higher (negative) correlation in the summer than in the winter. These findings, which indicate the potential of green spaces to mitigate the temperature across campus, especially in the summer, align with the findings of [81], who reported that the relationship between LST and NDVI is much weaker than with NDBI, and varies by season. Similar results were suggested by [82], who showed a strong negative correlation between LST and NDVI, especially during the warm months.
The results suggest that, in the period from 2014 to 2019, changes occurred in the spatial patterns of LST hot and cold spots across campus. These changes can be explained by differences in the characteristics of the surface; for example, the conversion of green vegetation (including native vegetation) to built-up and paved land cover; for example, the campus sports stadium was converted from green vegetation to artificial turf (AT) in 2018. While AT has become increasingly popular for sports, residential and commercial uses [83], it often results in significantly higher temperatures compared to natural green surfaces. According to [84], due to the lack of evaporation, the temperatures of AT surfaces can exceed those of natural grass by as much as 21 °C and of air temperature by 17 °C.
Lastly, we investigated the relationship between the rate of change of LST over time and the rate of change in the characteristics of land cover and land use. We did this by calculating a SOC measure, which indicates the slope of a regression line between LST, NDVI, NDBI, and UI and time. As with the results described above, we found a positive correlation between the SOC of LST and the SOCs of NDBI and UI, and a negative correlation with the SOC of NDVI, again with stronger correlations in summer than winter.
To summarize, the findings we report in this study demonstrate how remotely sensed, derived LST measures can be used to assess temporal changes in surface temperatures at local scales (the scale of a university campus). Our findings also illustrate the impact that patterns of land cover and land use have on micro-scale SUHI. Understanding the impact of LULC on the SUHI is fundamental for urban planning, from a regional to a micro-scale, and should be a top priority for any planning that aims to promote not only sustainable and livable cities, but also sustainable, livable and inviting university campuses. We note several limitations to this study and the potential extension for future evaluation. First, in this study, we relied on seasonal composites based on Landsat-8 observations, which were calculated as the average value of all scenes in an entire season. These measures do not capture daily and weekly fluctuation, due to microclimate conditions. Second, because Landsat-8 data are available only from 2013, this study examined the relationship between changes in LST and LULC in the period from 2013 to 2019. A future study could evaluate this relationship also with Landsat-7 observations (which are available from 1999), although they do not include the TIR sensor. Third, the main objective of this study was to estimate LST measures based on remotely sensed data (Landsat). Although previous studies have evaluated the accuracy of the remotely sensed estimations using in situ data, an extension to this study would collect additional in situ measurements in the university campus to validate the accuracy of the predictions. Finally, we related the observed changes in LST to changes in patterns of LULC (for example, the conversion of the sports stadium from natural to artificial turf). An extension to this study could validate these assumptions by including additional control and comparison studies, including through in situ measurements.

Conclusions
The year 2019 is likely to have been the second warmest recorded by instruments, and the past five years and the decade 2010-2019 were the warmest on record [85]. While these trends have been observed globally, arid and semi-arid countries, such as Saudi Arabia, are especially vulnerable to the challenges imposed by climate change. Goal 11 of the 2030 Agenda for Sustainable Development ("sustainable cities and communities") seeks to "make cities and human settlements inclusive, safe, resilient, and sustainable." In addition, one of the programs of the Saudi Arabian National Transformation Program is the Quality of Life Program, which aims to enhance livability in Saudi cities. While achieving this target is challenging, sustainable urban planning can help mitigate some of the impacts of climate change.
Evidence of the impact of climate change on the living environment includes the intensified SUHI phenomenon. The SUHI effect characterizes not only cities, but also smaller geographical scales, including university campuses. Remotely sensed data is effective for the estimation of LSTa fundamental proxy for SUHI. Previous studies have used remotely sensed derived LST, primarily for analysis at large scales. Here, we highlight its effectiveness for LST estimation on a micro-scale (university campus). The results of this study demonstrate (1) how thermal remote sensing can be used for localized LST estimations; (2) the significance of green spaces in moderating the surface temperature in a university campus-especially during the hot summer of an arid region; (3) that NDBI and UI spectral indices are superior to NDVI in explaining LST variations.
Developing methods to determine the heat island effect on campuses and understand the relationship with land use, land cover, and planning, can facilitate improvements in the outdoor experience of students and staff on campus. These methods can also be used to provide recommendations for future sustainable planning that would ultimately improve the quality of life on campus and help universities make energy consumption more efficient and become more sustainable.

Conflicts of Interest:
The authors declare no conflict of interest.
Appendix A: Table A1. Description of KAU sites, as presented in Figure 2.


It is located in the north of the University.  Open area.  Land cover is asphalt (parking area and the old driving school). Some parts include old and unused buildings and warehouses.  It has some vegetation around the buildings.
2  It is located east of the campus.  Open area.  Land cover is compacted soil, and in the past, it was more natural, with natural vegetation.  Currently, it has no vegetation or buildings.
3  It is located in the central part of the campus, toward the north.  Parking space-not shaded.


Only the parking close to the buildings is shaded with a tensile structure.  Land cover is asphalt.  From the east and the west, surrounded by administrative granite buildings and classrooms.  There is no vegetation inside the area, but there are street trees in the north. It has open areas around mosques and pocket gardens.  The building materials used for elevations are paint, due to the budget of the project.  The building's construction does not include green building approaches.
8  It is located in the northern part of the campus.  Campus stadium and sports facilities.  It includes the main stadium, small football fields and main gym buildings, made of granite and cladding.  At the end of 2018, the natural turf was changed to artificial in the stadium and the football fields.


It is located in the southern part of the campus.  Open area.  Land cover is compacted land.  It has scattered natural vegetation.  The south of the area is Abdullah AlSuliman Street (main road).  In the northwest border is the old staff housing made by precast concrete.
10  It is located in the central part of the campus, toward the northwest  It is the main street to the higher administrivia buildings.  The road edges and the roundabout have various types of vegetation. The current land cover is asphalt with small buildings.  It has no vegetation inside it, but there are mature trees from the west and the north.