Recognition of Thermal Hot and Cold Spots in Urban Areas in Support of Mitigation Plans to Counteract Overheating : Application for Athens

Mitigation plans to counteract overheating in urban areas need to be based on a thorough knowledge of the state of the thermal environment, most importantly on the presence of areas which consistently demonstrate higher or lower urban land surface temperatures (hereinafter referred to as “hot spots” or “cold spots”, respectively). The main objective of this research study is to develop a methodological approach for the recognition of thermal “hot spots” and “cold spots” in urban areas during summer; this is accomplished with (a) the combined use of high and medium spatial resolution satellite data (Landsat 8 and Terra-MODIS, respectively); (b) the downscaling of the Terra-MODIS satellite data so as to acquire spatial resolution similar to the Landsat one and at the same time take advantage of the high revisit time as compared to the respective one of Landsat (16 days); and (c) the application of a statistical clustering technique to recognize “hot spots” and “cold spots”. The methodological approach was applied as a case study for the urban area of Athens, Greece for a summer period. Results demonstrated the capacity of the methodological approach to recognize “hot spots” and “cold spots”, revealed a strong relationship between land use and “hot spots” and “cold spots”, and showed that the average land surface temperature (LST) difference between the “hot spots” and “cold spots” can reach 9.1 ◦K.


Introduction
The thermal environment in several urban areas worldwide has experienced a rapid deterioration, leading to the overheating of cities.The increase of temperature may be attributed to the influence of the radiative properties of urban surfaces and the three-dimensional configuration and heat capacity of buildings [1][2][3][4][5][6].In response to the increased urban (air and land surface) temperatures and the strong relationship between high urban air temperatures and adverse heat-related health outcomes (e.g., thermal discomfort, mortality, and heat-related illness) [7][8][9][10][11][12], urban planners and decision makers are recognizing the growing need to create more attractive, thermally comfortable, and sustainable cities, especially as urban populations expand and climatic variability and extremes increase [13][14][15].Results of specific building adaptation studies for Athens, Greece showed that this is possible when appropriate and energy-efficient building technologies are used [16].
In support of mitigation plans, satellite imagery in the thermal infrared (TIR) can be used as a tool to identify thermal "hot spots" and "cold spots" and correlate the state of the thermal environment to land use composition [17][18][19][20]."Hot spots" and "cold spots" are defined as areas which consistently demonstrate higher or lower land surface temperature (LST), respectively, as compared to their surrounding areas.Taken that areas of high air temperature coincide with areas of LST when coarse spatial resolution data is used (with a resolution >30 m) [13], the recognition of "hot spots" and "cold spots" can also support the extraction of information on the spatial distribution of air temperatures.As a matter of fact, several studies have successfully correlated air and land surface temperatures, while several researchers have successfully estimated air temperature from remote sensing data [21][22][23][24][25][26].Coutts et al. [13] used very high resolution airborne thermal infrared data in conjunction with Landsat and MODIS data in order to assess the adequacy of satellite TIR data to "hot spot" recognition.They found that Landsat data are suitable for "hot spot" recognition and supportive for mitigation plans to counteract overheating due to the adequate spatial coverage, the free data accessibility, the routine capture, and the well-documented processes and corrections.A disadvantage for the use of Landsat TIR data is the low temporal resolution, as the period between two successive passes of the satellite is 16 days [27,28].This may result to an unavailability of cloud free images for a long period if only Landsat data are used.To overcome this problem, TIR data from the MODIS instrument on board the Terra satellite are used; such data have a daily temporal resolution, but fall back in spatial resolution (1 km) [29][30][31][32].To this end, disaggregation methods (i.e., the improvement of the coarser spatial resolution of TIR images by introducing spatial detail from a finer resolution image) for downscaling MODIS TIR need to be applied, so as to match the respective spatial resolution of the Landsat data.
Several disaggregation methods have been suggested in order to enhance the spatial resolution of the TIR data by linking TIR and visible and near infrared (VNIR) data [33][34][35].These methods utilize the relation between the land cover and the LST; at first the relation between the coarser resolution TIR and VNIR bands is estimated and is thereafter applied at the finer resolution VNIR data in order to obtain finer resolution TIR data [33].Different disaggregation methods in support of the MODIS LST product have been proposed in literature [36,37], while Bisquert et al. [33] applied and compared many of these methods to MODIS and Landsat sensors.Results showed that the linear regression between the normalized difference vegetation index (NDVI) and LST lead to better results, with RMSE values around ±2 K, compared to other regressions or to more sophisticated approaches (Neural Networks and Data Mining).
The main objective of the current research is to develop a methodological approach for the recognition of "hot spots" and "cold spots" in urban areas during summer; this is accomplished with the combined use of high and medium spatial resolution satellite data in the thermal infrared, from the Landsat 8 and Terra satellites, respectively.A secondary objective is to investigate the relationship of "hot spots" and "cold spots" with land use/land cover.Both objectives support the drafting of mitigation plans for counteracting the overheating of urban areas.The methodological approach for "hot spots" and "cold spots" recognition is presented in Section 2. In Section 3, the methodological approach is applied for the urban area of Athens, and results are discussed accordingly.

Methodological Approach
The adopted methodology for the identification of "hot spot" and "cold spot" areas is based on the use of satellite data.These areas are characterized by very high ("hot spots") and very low ("cold spots") LST values in multi-temporal satellite images, assuming stable meteorological conditions and very limited, if any, land cover changes.Consequently, the satellite-based approach needs to be replicated for more than one day in order to assess the spatial consistency of the hot and cold spot areas.To this end, the analysis cannot be based merely on Landsat, taken that it overpasses every 16 days.MODIS provides LST data on a daily basis, whereas the application of downscaling allows for LST products at the same spatial resolution as Landsat.The multi-temporal analysis is even more necessary for the implementation of the methodology in other regions that may be characterized by unstable weather conditions, in order to investigate any large spatial variations of the hot/cold spot areas.In the following sections the required steps are described in detail.
where ε is the surface emissivity, L sen is the radiance of the thermal channel TIRS-10, and γ, δ are two parameters given by where T sen is the brightness temperature of the thermal channel TIRS 10, b γ = 1324 for TIRS-10, and ψ 1 , ψ 2 , and ψ 3 are atmospheric functions given by where τ is the atmospheric transmission, L d is the downwelling radiance, and L u is the upwelling radiance.These parameters are calculated from the Atmospheric Correction Parameter Calculator provided by NASA (http://atmcorr.gsfc.nasa.gov)[40].Land surface emissivity (ε) was estimated using the NDVI Thresholds Method as proposed by Sobrino et al. [41].

Downscaling Modis Data
In this study, the adopted disaggregation method for downscaling MODIS TIR data is based on the linear regression approach that leads to better results, as stated in [33].Linear regression between Terra MODIS and Landsat 8 TIRS is applied in the event that no Landsat overpass coincides with a MODIS one, with low viewing angle, thus a close Landsat NDVI image is used under the assumption that NDVI remains constant.
Initially, a normalization process needs to be applied to Landsat data in order to minimize the discrepancies that result from the differences in spectral resolution, atmospheric correction, viewing angle, pixel footprint, time acquisition, etc., between the two sensors [33].Subsequently, according to the disaggregation procedure followed in this study and presented schematically in Figure 1, the relationship between the coarser MODIS NDVI and LST at 960 m is obtained and subsequently applied to the MODIS NDVI at 960 m, resulting in the LST' 960 m.The residuals as calculated by subtracting the LST' 960 m from the LST 960 m, were added to the finer LST' 30 m, which resulted from the application of the initial relationship to the finer Landsat NDVI 30 m, giving eventually the downscaled MODIS LST 30 m.
where ε is the surface emissivity, is the radiance of the thermal channel TIRS-10, and γ, δ are two parameters given by where is the brightness temperature of the thermal channel TIRS 10, = 1324 for TIRS-10, and , , and are atmospheric functions given by where τ is the atmospheric transmission, is the downwelling radiance, and is the upwelling radiance.These parameters are calculated from the Atmospheric Correction Parameter Calculator provided by NASA (http://atmcorr.gsfc.nasa.gov)[40].Land surface emissivity (ε) was estimated using the NDVI Thresholds Method as proposed by Sobrino et al. [41].

Downscaling Modis Data
In this study, the adopted disaggregation method for downscaling MODIS TIR data is based on the linear regression approach that leads to better results, as stated in [33].Linear regression between Terra MODIS and Landsat 8 TIRS is applied in the event that no Landsat overpass coincides with a MODIS one, with low viewing angle, thus a close Landsat NDVI image is used under the assumption that NDVI remains constant.
Initially, a normalization process needs to be applied to Landsat data in order to minimize the discrepancies that result from the differences in spectral resolution, atmospheric correction, viewing angle, pixel footprint, time acquisition, etc., between the two sensors [33].Subsequently, according to the disaggregation procedure followed in this study and presented schematically in Figure 1, the relationship between the coarser MODIS NDVI and LST at 960 m is obtained and subsequently applied to the MODIS NDVI at 960 m, resulting in the LST' 960 m.The residuals as calculated by subtracting the LST' 960 m from the LST 960 m, were added to the finer LST' 30 m, which resulted from the application of the initial relationship to the finer Landsat NDVI 30 m, giving eventually the downscaled MODIS LST 30 m.

"Hot Spot" Recognition
A statistical approach is applied so as to detect-and cluster-the presence of "hot spot" or "cold spot" areas over a study area by calculating G statistics [42] for each pixel in the image, as shown in Equations ( 2)-( 4).The resulting z-scores show where the pixels with either high or low LST values cluster spatially, while they are also measures of statistical significance (Table 1).It is important to note that a pixel with a high LST value may not be a statistically significant "hot spot" if its neighbor pixels have much lower LST values.To be a statistically significant "hot spot", a pixel must have a high value and be surrounded by other pixels with high values as well.The methodological approach suggests the use of all eight neighbor pixels in addition to the central one, in order to calculate the local sum and subsequently obtain the local average (i.e., local sum/9) of the main pixel (Figure 2).This procedure was selected in order to ensure that LST "hot spots" are related to air temperature "hot spots" given that a larger spatial domain is expected to relate better with the near-surface air temperature source area.This pixel value is then compared to the respective average value of the image (sum of all pixels in the study area/number of pixels in the study area).When the local average is very different from the image average, and the difference is too large to be the result of random chance, a statistically significant z-score is found.For statistically significant positive z-scores, the larger the z-score, the more intense is the clustering of high values ("hot spot").For statistically significant negative z-scores, the smaller the z-score, the more intense is the clustering of low values ("cold spot").
where G i * is the resultant G statistics (z-scores and p-values) for pixel i, x j is the LST value for pixel j, w i,j is the spatial weight between pixel i and neighboring pixel j, n is equal to the total number of pixels, X is the mean LST of all pixels, and S the variance.

"Hot Spot" Recognition
A statistical approach is applied so as to detect-and cluster-the presence of "hot spot" or "cold spot" areas over a study area by calculating G statistics [42] for each pixel in the image, as shown in Equations ( 2)-( 4).The resulting z-scores show where the pixels with either high or low LST values cluster spatially, while they are also measures of statistical significance (Table 1).It is important to note that a pixel with a high LST value may not be a statistically significant "hot spot" if its neighbor pixels have much lower LST values.To be a statistically significant "hot spot", a pixel must have a high value and be surrounded by other pixels with high values as well.The methodological approach suggests the use of all eight neighbor pixels in addition to the central one, in order to calculate the local sum and subsequently obtain the local average (i.e., local sum/9) of the main pixel (Figure 2).This procedure was selected in order to ensure that LST "hot spots" are related to air temperature "hot spots" given that a larger spatial domain is expected to relate better with the near-surface air temperature source area.This pixel value is then compared to the respective average value of the image (sum of all pixels in the study area/number of pixels in the study area).When the local average is very different from the image average, and the difference is too large to be the result of random chance, a statistically significant z-score is found.For statistically significant positive z-scores, the larger the z-score, the more intense is the clustering of high values ("hot spot").For statistically significant negative z-scores, the smaller the z-score, the more intense is the clustering of low values ("cold spot").
where Gi* is the resultant G statistics (z-scores and p-values) for pixel i, xj is the LST value for pixel j, wi,j is the spatial weight between pixel i and neighboring pixel j, n is equal to the total number of pixels, is the mean LST of all pixels, and S the variance.

Application Area
In this section, the methodological approach is applied for the area denoted in Figure 3 with the red polygon representing the Municipality of Athens.The area has coverage of about 39 km 2 and includes high-density urban residential areas, commercial and industrial areas, transport modes, and the associated road network.Urban space is characterized in particular by a high degree of mixed land-use, limited green spaces or open public spaces, and limited pathways that allow for the influx of airflow from the surrounding countryside [43,44].

Application Area
In this section, the methodological approach is applied for the area denoted in Figure 3 with the red polygon representing the Municipality of Athens.The area has coverage of about 39 km 2 and includes high-density urban residential areas, commercial and industrial areas, transport modes, and the associated road network.Urban space is characterized in particular by a high degree of mixed land-use, limited green spaces or open public spaces, and limited pathways that allow for the influx of airflow from the surrounding countryside [43,44].

Data
Landsat 8 images (morning pass) were used in order to calculate the LST of the urban agglomeration of Athens.The acquisition dates were 27 July 2016 and 12 August 2016, and the selection of these dates were in compliance with the need for similar meteorological conditions.As a matter of fact, during this period no precipitation occurred over the Athens basin, while the wind and the temperature regime were similar for all days under study.In particular, clear sky conditions were observed from 26 July 2016 until 12 August 2016, the mean temperature was 30.5 °C, the maximum air temperature was 36.9 °C, the minimum air temperature was 23.4 °C, and the mean wind speed was 5.5 km/h southwards-typical Mediterranean summer conditions.
MODIS LST data from Terra satellite (MOD11A1 product [45]) were selected so as to refer to the acquisition time of the Landsat 8 images.They were used for four days with low MODIS viewing angle (<±30°), as high angles lead to higher errors, between the two Landsat 8 overpasses (27 July and 12 August).These days are the 26th and 28th of July as well as the 2nd and 11th of August.The selected time period along with the acquisition time of the satellite data and the satellite viewing angle are presented in Table 2.In addition, the MODIS Terra VNIR data (MOD13Q1) were used for the NDVI calculation, while the Urban Atlas dataset [46] was employed for the land use analysis.

Data
Landsat 8 images (morning pass) were used in order to calculate the LST of the urban agglomeration of Athens.The acquisition dates were 27 July 2016 and 12 August 2016, and the selection of these dates were in compliance with the need for similar meteorological conditions.As a matter of fact, during this period no precipitation occurred over the Athens basin, while the wind and the temperature regime were similar for all days under study.In particular, clear sky conditions were observed from 26 July 2016 until 12 August 2016, the mean temperature was 30.5 • C, the maximum air temperature was 36.9 • C, the minimum air temperature was 23.4 • C, and the mean wind speed was 5.5 km/h southwards-typical Mediterranean summer conditions.
MODIS LST data from Terra satellite (MOD11A1 product [45]) were selected so as to refer to the acquisition time of the Landsat 8 images.They were used for four days with low MODIS viewing angle (<±30 • ), as high angles lead to higher errors, between the two Landsat 8 overpasses (27 July and 12 August).These days are the 26th and 28th of July as well as the 2nd and 11th of August.The selected time period along with the acquisition time of the satellite data and the satellite viewing angle are presented in Table 2.In addition, the MODIS Terra VNIR data (MOD13Q1) were used for the NDVI calculation, while the Urban Atlas dataset [46] was employed for the land use analysis.

Land Surface Temperature
The LST evolution during the period under study is presented in Figure 4.The surface temperature patterns are similar in all images, whilst the mean LST values range from 312.505 • K to 314.076 • K, indicating rather stable conditions from the 26th of July to the 12th of August.The latter is in accordance with the stable meteorological conditions characterising the period under study.Furthermore, all images clearly depict that the western part of the Municipality of Athens exhibits higher LST than the eastern part, mainly due to the presence of industrial buildings and the high-density urban areas.On the contrary, the southeastern part shows lower LST as it is less densely built and consists of many urban green areas.Interestingly, small areas inside the boundaries of the municipality of Athens exhibit very high LST values that exceed 319 • K in all cases.This temperature value is much higher than the average LST of the municipality, indicating the presence of "hot spots".Intra-urban variations of LST are clearly depicted in the Landsat images (Figure 4b,f) and in the downscaled images too (Figure 4a,c-e).
It should be mentioned that the Landsat LST error is typically less than 1.5 K [39] and the accuracy of the MODIS LST product is better than 1 K [47].Furthermore, a discrepancy of ±1.9 K is observed between MODIS downscaled LST and Landsat LST [33], when Landsat overpass coincided with a MODIS one, the latter with low viewing angle.The summation of the above leads to the worst case scenario error of the downscaled MODIS LST.The above errors do not affect the recognition of the "hot spots" and "cold spots" patterns in the urban area under investigation; this is due to the fact that the recognition is based on differences between local and image averages, rather than to absolute LST values.

Recognition of "Hot Spots"
The calculation of the G statistics for each LST image resulted in the assignment of a z-value to each pixel.In order to define the "hot spot" areas, the suitable threshold for the z-values (Table 1) was determined as follows.Usually a z-value over 2.58 (99% confidence level) indicates a statistically significant result, but a more rigorous threshold of a z-value over 3.3 (99.9% confidence level) can be used in order to narrow down the number and the extent of the "hot spots".Comparing the results using both thresholds, it was found out that the less rigorous threshold of 2.58 increases the extent of the "hot spot" areas by 29% and the number of "hot spots" by 27%.Additionally, this threshold increases the extent and the number of "cold spots" by 68% and 109% respectively.Taken the above, the threshold 3.3 of z-value was chosen as the most suitable for this study.
Subsequently, "hot spots" and "cold spots" were determined when a pixel was characterized as a "hot spot" or "cold spot", respectively, in at least four images (i.e., over 50% of the available data).Based on the above, three categories of "hot spots" and "cold spots" emerged according to the number of times they were recognized.Table 3 provides the statistical analysis for each "hot spot" and "cold spot" category.A gradual change in the minimum, maximum, and mean values is evident.
The average LST of a 6-day "hot spot" is 4.43 • K higher than the average LST of the areas that are neither "hot spots" nor "cold spots".Even the 4-day "hot spots" have a substantial difference from these areas, reaching 2.16 • K.In addition, the average LST difference between the 6-day "hot spots" and the 6-day "cold spots" reaches 9.1 • K.
Climate 2018, 6, x FOR PEER REVIEW 7 of 13 difference from these areas, reaching 2.16 °K.In addition, the average LST difference between the 6-day "hot spots" and the 6-day "cold spots" reaches 9.1 °K.   Figure 5 depicts the "hot spots" and "cold spots" of the Municipality of Athens.The majority of the "hot spots" are located in the western part of the city, which has greater LST values.The old industrial area (denoted as A in Figure 5) in the western boundaries of the municipality appears as a large "hot spot" due to the homogeneous land cover of impervious materials and the lack of vegetation.Some of the "hot spots" are located in the historic center of Athens (denoted as B in Figure 5), where major sightseeing places are located.On the contrary, the majority of "cold spots" are located at the eastern part of the municipality in the vicinity of urban green areas.The above results in conjunction with the LSΤ differences as derived between the 6-day and the 4-day "hot spots" (2.27 °K) highlight the need to prioritize the interventions to the 6-day "hot spot" areas.Furthermore, land use analysis (Figure 6) lead to the recognition of various land use patterns in each "hot spot" or "cold spot" category, indicating the need for customized interventions.The extracted map of "hot spots" and "cold spots" in the Municipality of Athens (Figure 5) can be easily processed in Geographic Information Systems (GIS) environment by urban planners and decision makers as it contains geospatial information.Using GIS software, it was calculated that the "hot spots" occupy an area of 2.75 km 2 (7% of the total municipality area) and the "cold spots" occupy an area of 2.78 km 2 .Moreover, various data, such as census, energy consumption, greenery, age of buildings, etc., can be overlaid, thus enabling the design of mitigation plans to counteract overheating.
The above results indicate that further investigation of the "hot spot" and "cold spot" areas is needed in order to examine the relation of land use/land cover with the presence of "hot spots" or "cold spots".To this end, land use data from the Urban Atlas were used; Table 4 provides the land use information for the municipality of Athens.High-density urban fabric covers more than 45% of the area, while proportions of land are covered by roads (26.73%), industrial and commercial structures (19.83%), and urban parks (15.53%).The above results in conjunction with the LST differences as derived between the 6-day and the 4-day "hot spots" (2.27 • K) highlight the need to prioritize the interventions to the 6-day "hot spot" areas.Furthermore, land use analysis (Figure 6) lead to the recognition of various land use patterns in each "hot spot" or "cold spot" category, indicating the need for customized interventions.The above results in conjunction with the LSΤ differences as derived between the 6-day and the 4-day "hot spots" (2.27 °K) highlight the need to prioritize the interventions to the 6-day "hot spot" areas.Furthermore, land use analysis (Figure 6) lead to the recognition of various land use patterns in each "hot spot" or "cold spot" category, indicating the need for customized interventions.

Conclusions
In this study, a methodological approach was developed in order to support mitigation plans for counteracting overheating in urban areas.The methodological approach recognizes-with the use of satellite remote sensing data in the thermal infrared-areas that consistently demonstrate high or low land surface temperatures during summer time.A critical criterion in the methodological approach is that a pixel (an area) with a high LST value may not be a "hot spot" if its neighbor pixels (areas) have much lower LST values.To be considered as a "hot spot", a pixel (area) must have high LST value and must be surrounded by pixels (areas) with high LST values as well.This criterion overcomes trivial approaches to represent the thermal environment of urban areas, based merely on the LST value on a pixel by pixel basis.
The methodological approach is based on: (a) the combined use of various satellite data, in order to circumvent-through a downscaling technique-the inherent difficulty which characterizes such data, namely the inverse relationship between spatial and temporal resolutions; (b) the estimation of LST on a daily basis and at a spatial resolution of 30 meters; (c) the statistical analysis of the land surface temperatures so as to recognize and cluster "hot spots" and "cold spots"; and (d) the correlation of the "hot spots" and "cold spots" temporal endurance and spatial extent to land use classes.
The methodological approach was applied for the Municipality of Athens on the basis of six-either original or downscaled-high-resolution LST images during a 16-day summer period; "hot spots" were recognized when a pixel was characterized as a "hot spot" in at least in four images (i.e., over the 50% of the available data) using G statistics.
The average LST of a 6-day "hot spot" is 4.43 • K higher than the average LST of the areas that are neither "hot spots" nor "cold spots".Even the 4-day "hot spots" have a substantial difference from these areas, reaching 2.16 • K.In addition, the average LST of the 6-day "hot spot" is 2.27 • C higher than the 4-day "hot spots", whereas the average LST difference between the 6-day "hot spots" and the 6-day "cold spots" reaches 9.1 • K.
Furthermore, land use analysis revealed the strong relationship between land use on the one hand and "hot spots" and "cold spots" on the other."Hot spots" are almost exclusively covered by industrial/commercial, high-density urban fabric, and street type surface materials.Conversely, "cold spots" have, as expected, an increased presence of urban greenery.
The methodological approach, as developed in this study, is replicable to all urban areas and can be highly supportive to mitigation plans to counteract overheating as it can detect "hot spot" or "cold spot" areas in fine resolution (30 × 30 m).Moreover, various data, such as census, energy consumption, greenery, age of buildings, etc., can be overlaid, thus enabling a thorough design of mitigation plans to counteract overheating.Further work to be done refers to the application of the methodological approach for night time satellite images from MODIS as Landsat night time images are very limited.The successful identification of nighttime "hot spots" will lead to an improved understanding of the thermal environment of the city.

Figure 2 .
Figure 2. The calculation of the local LST sum for a pixel under consideration (red) includes all of its neighbors (orange).

Figure 2 .
Figure 2. The calculation of the local LST sum for a pixel under consideration (red) includes all of its neighbors (orange).

Figure 3 .
Figure 3.The urban agglomeration of Athens and the administrative boundaries of the municipality of Athens (red polygon).

Figure 3 .
Figure 3.The urban agglomeration of Athens and the administrative boundaries of the municipality of Athens (red polygon).

13 Figure 5 .
Figure 5. Location of the hot/cold spots within Athens Municipality.(A: old industrial area of Athens, B: the historic center of Athens).

Figure 5 .
Figure 5. Location of the hot/cold spots within Athens Municipality.(A: old industrial area of Athens, B: the historic center of Athens).

Figure 5 .
Figure 5. Location of the hot/cold spots within Athens Municipality.(A: old industrial area of Athens, B: the historic center of Athens).

Figure 6 .
Figure 6.Land use percentages of the "hot spot" and "cold spot" categories.Figure 6. Land use percentages of the "hot spot" and "cold spot" categories.

Figure 6 .
Figure 6.Land use percentages of the "hot spot" and "cold spot" categories.Figure 6. Land use percentages of the "hot spot" and "cold spot" categories.

Table 1 .
Critical p-values and z-scores for different confidence levels.

Table 1 .
Critical p-values and z-scores for different confidence levels.Significance Level (p Value) Critical Value (z Score) Confidence Level

Table 2 .
Selected dates and characteristics of the LST data.

Table 3 .
LST statistics for the "hot spots" and "cold spots" categories (in • K).

Table 4 .
Land use percentages for the municipality of Athens based on the Urban Atlas data.