Identiﬁcation of SUHI in Urban Areas by Remote Sensing Data and Mitigation Hypothesis through Solar Reﬂective Materials

: The urban heat island (UHI) is an increasingly widespread phenomenon of concern to the wellbeing and the health of populations living in urban environments. The SUHI (Surface UHI) is directly related to UHI and inﬂuences its extension and intensity. Satellite images in the thermal infrared spectral region can be used to identify and study the SUHI. In this work, Landsat 8 TIR images were acquired to study the SUHI of a medium-sized municipality of the Po valley in the northern part of Italy. An additional Worldview 3 satellite image was used to classify the study area and retrieve the surface albedo of building roofs. Using the Local Climate Zone approach, existing roof materials were virtually replaced by solar reﬂective materials, and the mitigation potential of the SUHI and the UHI was quantiﬁed. This virtual scenario shows a decrease in the overheating of building roofs with respect to the ambient temperature of up to 33% compared to the current situation in the industrial areas. Focusing on UHI intensity, the air temperature decrease could be up to 0.5 ◦ C.


Introduction
The global population is increasingly moving from small rural areas to cities. In 2018, 55.3% of the world's population lived in cities, and the threshold of 60% should be reached by 2030. There were 548 cities with at least one million inhabitants in 2018, but that could rise to 706 in 2030 [1,2]. Over-urbanization, combined with the growing problems related to climate change and global warming, will lead in the near future to the worsening of many of the social, health, and environmental negative effects in urban areas [3,4].
A main problem of cities is the Urban Heat Island (UHI), a phenomenon which is well-known and widely reported in the literature by the scientific community [5][6][7][8]. The UHI is a microclimatic phenomenon occurring in metropolitan areas and consists of a significant temperature increase in urban areas compared to the neighboring peripheral and rural areas [9]. The UHI is mainly due to manmade surfaces, such as concrete and asphalt, which adsorb more solar radiation than natural surfaces such as vegetation and bare soils. Furthermore, urban impervious surfaces, such as asphalt and concrete, show very low evapotranspiration, resulting in a low dissipation of energy as latent heat [10,11]. This effect is more pronounced in summer and results in a significant difference in mean temperature between the city and the suburbs [12]. Heat stored during the day is slowly dissipated by radiation, starting in the late afternoon and during the night, while in rural areas, bare soils and vegetation dissipate heat faster [13].
In the scientific literature, a distinction is made between the Surface Urban Heat Island (SUHI) and the (atmospheric) Urban Heat Island (UHI). The latter represents the difference in air temperature between the urban and rural areas while the former refers to ground temperatures. The SUHI is usually at a maximum during the day, when the sun heats the

Study Area
The study area of this work is the municipality of Reggio Emilia (Figure 1). Reggio Emilia has 160,000 inhabitants and covers an area of 231 km 2 . The city is located in the Po Valley, halfway between the Po river and the Apennine mountains [46]. High population density, industrialization processes, and intensive agriculture are characteristics of the Po Valley of which Reggio Emilia is a part [47]. This is an area characterized by temperate continental climatic conditions. The summers are hot and muggy, and the winters are cold and humid. The rainfall is distributed throughout the year, but the rainiest seasons are autumn and spring, while the driest ones are winter and summer [48][49][50]. The winters are typically characterized by the thermal inversion phenomenon, with an accumulation of pollutants near to the soil surface; the autumns and winters are frequently foggy, mainly in the plains [50]. decrease in the building roof surface temperature (mitigation action for the SUHI) and the air temperature decrease for the UHI.

Study Area
The study area of this work is the municipality of Reggio Emilia (Figure 1). Reggio Emilia has 160,000 inhabitants and covers an area of 231 km 2 . The city is located in the Po Valley, halfway between the Po river and the Apennine mountains [46]. High population density, industrialization processes, and intensive agriculture are characteristics of the Po Valley of which Reggio Emilia is a part [47]. This is an area characterized by temperate continental climatic conditions. The summers are hot and muggy, and the winters are cold and humid. The rainfall is distributed throughout the year, but the rainiest seasons are autumn and spring, while the driest ones are winter and summer [48][49][50]. The winters are typically characterized by the thermal inversion phenomenon, with an accumulation of pollutants near to the soil surface; the autumns and winters are frequently foggy, mainly in the plains [50].
The meteorological data of the area were measured by the "Reggio Emilia San Lazzaro" meteorological station (Lat. 44.69, Lon. 10.67), located inside the San Lazzaro university campus and part of the network of measurement stations of the Geophysical Observatory of the University of Modena and Reggio Emilia [51].
Reggio Emilia has been divided into 63 municipal districts by the local administration. These areas have very different sizes: there are many districts in the historic center with small dimensions, while in the rural areas they cover large portions of the territory.

Local Climate Zones (LCZs) Approach
The LCZs were defined by Stewart and Oke as "Regions of uniform surfaces cover, structure, material, and human activity that span hundreds of meters to several kilometers in horizontal scale" [22]. The LCZs were introduced here to assess the potentiality of mitigation actions for SUHIs and, consequently, for UHIs at a local scale. In Italian historical cities such as Reggio Emilia, the urbanization process throughout the years caused an irregular growth of the city that has resulted in commercial/industrial areas close to the city center or generally mixed with residential areas. Mitigation strategies must consider the The meteorological data of the area were measured by the "Reggio Emilia San Lazzaro" meteorological station (Lat. 44.69, Lon. 10.67), located inside the San Lazzaro university campus and part of the network of measurement stations of the Geophysical Observatory of the University of Modena and Reggio Emilia [51].
Reggio Emilia has been divided into 63 municipal districts by the local administration. These areas have very different sizes: there are many districts in the historic center with small dimensions, while in the rural areas they cover large portions of the territory.

Local Climate Zones (LCZs) Approach
The LCZs were defined by Stewart and Oke as "Regions of uniform surfaces cover, structure, material, and human activity that span hundreds of meters to several kilometers in horizontal scale" [22]. The LCZs were introduced here to assess the potentiality of mitigation actions for SUHIs and, consequently, for UHIs at a local scale. In Italian historical cities such as Reggio Emilia, the urbanization process throughout the years caused an irregular growth of the city that has resulted in commercial/industrial areas close to the city center or generally mixed with residential areas. Mitigation strategies must consider Atmosphere 2022, 13, 70 4 of 24 the replacement of roofs with different solar reflective materials depending on the area of the city. City center roofs are almost exclusively made of clay tiles because they are part of the cultural heritage and therefore it is very difficult to consider a change in the color of the materials. On the other hand, flat roofs on residential or commercial/industrial buildings are suitable for replacement with cool roofs. As these zones are often contiguous to the city center, mitigation effects provided by new roofs can lend benefits to the surrounding areas. For these reasons, it is preferable to act considering single areas rather than the entire urban area. Thus, after the SUHI identification by the Landsat 8 TIR images, some LCZs in which it is possible and advantageous to consider the replacement of roofs will be considered. Some rural areas will also be considered for comparison.

Data Set
For this study, both satellite images and vector files were used. In particular, the satellite data from Landsat 8 were used to retrieve the LST maps, while an image from the Worldview3 satellite allowed us to classify urban surfaces and to calculate the surface albedo. The municipality of Reggio Emilia and the topographic geodatabase of the cartographic archive of Emilia Romagna Region (Geoportale-Emilia Romagna) provided vector shapefiles and high resolution orthoimages, which were used as a support for the image processing and classification.

Landsat 8 Images
Landsat 8 is an American satellite for Earth observation, in operation since 2013 and developed and managed by NASA and USGS (United States Geological Survey) [42]. The spacecraft voyages at an altitude of 705 km, completing one orbit every 99 min and repeating the cycle every 16 days. It moves along a near polar orbit, flying over every point of the Earth at 10:00 a.m. +/− 15 min (helio synchronous orbit).
The satellite is equipped with two sensors: the OLI (Operational Land Imager) collects image data from bands in the VIS-NIR spectral region; the TIRS (Thermal Infrared Sensor) operates with two bands in the TIR spectral region. The OLI bands have a spatial resolution of 30 m; the TIRS bands have a spatial resolution of 100 m (but the images are resampled at 30 m). Bands features are shown in Table 1. Both the OLI and the TIRS bands were used to calculate the LST and the vegetation information of the study area. A time series of TIR satellite images is required to achieve a full assessment of the time and space variability of the SUHI because the latter is driven by climate and meteorological conditions. The years from 2016 to 2021 were considered for this study and, because the UHI shows its effects particularly during spring and summer, only the months from May to September of every year were analyzed. During the cold months, the UHI has a lower intensity and it results in a free heat source [46,52,53]. First of all, the Landsat 8 images available in the period of interest were considered. The first selection rejected all the cloudy images and the ones collected after a storm; in fact, clouds obstruct satellite visibility, and wet surfaces show lower values of LST due to rainfall. Then, for the remaining images, the daily weather conditions were compared to the mean weather conditions of the respective months in order to identify the most suitable image for every investigated month. The daily and monthly weather data were retrieved from the meteorological station at the San Lazzaro university campus of Reggio Emilia. Table 2 lists all the processed images, with the relative daily weather conditions, and the monthly average meteorological data. Table 2. Daily and monthly meteorological data measured at San Lazzaro university campus, Reggio Emilia.

Daily Weather Conditions Monthly Averaged Weather Conditions
Year Day T min

Worldview Image
A Worldview3 (WV3) image was acquired to classify the urban surfaces (especially building roofs) and to compute the surface albedo. This commercial satellite has eight bands in the visible and near infrared, with a spatial resolution of 1.6 m, allowing the assessment of many physical properties of the surfaces [54][55][56]; contrariwise, aerial photos return only "color" information. Figure 2 and Table 3 show the image used and the sensor characteristics. The image was acquired by the satellite on 2 June 2017 at 10:58:08 CET. Although the surface albedo depends on the seasons for vegetation, it is practically constant for the artificial surfaces, with some slight changes over the years due to surface deterioration (weathering, soiling, biological growth, etc.), especially in the first month after construction [57]. The impact of ageing, solar elevation, rain, and cloud cover on the albedo for conventional and non-conventional materials has been deeply assessed in the scientific literature. The impervious surfaces showed no significative variations in the surface albedo over time, if measured across midday (with a delta of ± 2 h), and thus no significant seasonal variation [57][58][59].
As this study is focused on building roofs, the surface albedo could be calculated with only one WV3 image. Several works are available in the scientific literature to prove this behavior of impervious surfaces [57][58][59]; for example, [58] found that the measured albedo of impervious surfaces is constant over time and has no significant seasonal variation. The WV3 image was acquired in 2017; thus, the surface albedo could be slightly different in 2021 due to the aging of the roofs.
However, one of the purposes of this work is the assessment of the benefits (in terms of the albedo and the subsequent decrease in surface temperature) that could be achieved  Although the surface albedo depends on the seasons for vegetation, it is practically constant for the artificial surfaces, with some slight changes over the years due to surface deterioration (weathering, soiling, biological growth, etc.), especially in the first month after construction [57]. The impact of ageing, solar elevation, rain, and cloud cover on the albedo for conventional and non-conventional materials has been deeply assessed in the scientific literature. The impervious surfaces showed no significative variations in the surface albedo over time, if measured across midday (with a delta of ± 2 h), and thus no significant seasonal variation [57][58][59].
As this study is focused on building roofs, the surface albedo could be calculated with only one WV3 image. Several works are available in the scientific literature to prove this behavior of impervious surfaces [57][58][59]; for example, [58] found that the measured albedo of impervious surfaces is constant over time and has no significant seasonal variation. The WV3 image was acquired in 2017; thus, the surface albedo could be slightly different in 2021 due to the aging of the roofs.
However, one of the purposes of this work is the assessment of the benefits (in terms of the albedo and the subsequent decrease in surface temperature) that could be achieved if the roofs were replaced with solar reflective materials. In this framework, an overestimation of the albedo of the current situation (that does not take into account the aging) could be associated with a conservative hypothesis.

Support Information
Aerial images were acquired at a spatial resolution of 30 cm by CGR S.p.a. in 2018. Very high resolution has been used in support of the surface classification, for the choice of algorithm training area, and also for validation. These images are not useful for the surface classification as they do not provide any information about the spectral signature of the surfaces (simple RGB photographs). On the contrary, the WV3 image provides 8 spectral bands in the visible-near infrared spectra, expressed in surface reflectance with a spatial resolution of 1.6 m, which is suitable for the study of building roofs.
The data set used also includes some urban information in vector format provided by the municipality of Reggio Emilia and Geoportale-Emilia Romagna: (1) Buildings: shapefile containing geometries and information on building roofs; (2) Streets: layer containing street names, used for localizing buildings; (3) Municipal districts: shapefile of the municipal districts, used to identify LCZs. Meteorological data have been retrieved from the meteorological station at the San Lazzaro university campus, Reggio Emilia. This station belongs to the network of meteorological stations of the Geophysical Observatory of Modena [60]. if the roofs were replaced with solar reflective materials. In this framework, an overestimation of the albedo of the current situation (that does not take into account the aging) could be associated with a conservative hypothesis.

Support Information
Aerial images were acquired at a spatial resolution of 30 cm by CGR S.p.a. in 2018. Very high resolution has been used in support of the surface classification, for the choice of algorithm training area, and also for validation. These images are not useful for the surface classification as they do not provide any information about the spectral signature of the surfaces (simple RGB photographs). On the contrary, the WV3 image provides 8 spectral bands in the visible-near infrared spectra, expressed in surface reflectance with a spatial resolution of 1.6 m, which is suitable for the study of building roofs.
The data set used also includes some urban information in vector format provided by the municipality of Reggio Emilia and Geoportale-Emilia Romagna: (1) Buildings: shapefile containing geometries and information on building roofs; (2) Streets: layer containing street names, used for localizing buildings; (3) Municipal districts: shapefile of the municipal districts, used to identify LCZs. Meteorological data have been retrieved from the meteorological station at the San Lazzaro university campus, Reggio Emilia. This station belongs to the network of meteorological stations of the Geophysical Observatory of Modena [60].

Land Surface Temperature
The Landsat 8 images were obtained from the Earth Explorer portal as level 1 images; this means that the bands are in digital numbers (DN). The Semi-Automatic Classification Plugin (SCP) [61] was used to convert the OLI bands from DN data to surface reflectance

Land Surface Temperature
The Landsat 8 images were obtained from the Earth Explorer portal as level 1 images; this means that the bands are in digital numbers (DN). The Semi-Automatic Classification Plugin (SCP) [61] was used to convert the OLI bands from DN data to surface reflectance and the TIRS bands from DN data to brightness temperature. The SCP uses the DOS1 atmospheric corrections to do that (based on the Dark Object Method) [61].
Surface temperature depends on surface emissivity (considering Planck's equation for the emittance of a real body). This parameter has been calculated using the Normalized Dif-ference Vegetation Index (NDVI). In turn, the NDVI has been corrected from the mixed pixel (vegetated/not vegetated) effect using the Fractional Vegetation Cover (FVC) [46,62,63].
The NDVI index provides an estimation of the vegetation presence on the pixel and is defined as: where ρ RED represents the red band reflectance (band 4), and ρ NIR is the near-infrared band reflectance (band 5). The FVC index estimates the portion of an area covered by vegetation and can be calculated as: where NDVIs is the NDVI value for bare soil, set as 0.1, and NDVIv is the NDVI value for fully vegetated soil, set as 0.65. Surface emissivity is computed using the equation [61]: where ε s is the typical soil emissivity, set as 0.93, and ε v is the typical vegetation emissivity, set as 0.98. Finally, the LST is computed by [61]: where T b is the brightness temperature referred to band 10 (K), and λ is the wavelength (for band 10, λ = 10.8 µm). c 2 = 1.4388 × 10 −2 m ·K is the second radiative constant.
The LST was first computed in degrees Kelvin and then converted to Celsius. Here, it is important to point out the limits of the Landsat 8 data: the spatial resolution of the TIRS sensor and the acquisition time of the images.
TIR bands have a spatial resolution of 30 m (resampled from a starting resolution of 100 m), and this leads to the mixed pixel problem: several kind of surfaces can be melted in the same pixel (such as the roof of an industrial building and the adjacent tree-lined road) [64][65][66]. The TIR images were acquired at 10:00 a.m. CET (± 15 min); thus, it is not possible to detect the highest daily surface temperature values. Furthermore, some sheet metal roofs presenting a lower thermal inertia than tile roofs can result in being hotter at 10:00 am than the other roofs. This factor needs to be reported because materials with medium-high thermal inertia have a stronger influence on UHIs given that they can release heat more slowly over time. However, the Landsat 8 is the satellite currently operating with the highest resolution in the TIR region; thus, it represents the better option for a multitemporal image acquisition. These limits highlight that the LST maps from Landsat 8 are not suitable for a surface temperature analysis at a detailed level, i.e., they cannot be used for the single-building roof temperature, but they can provide useful information for larger areas, such as an LCZ or a whole municipality [40,[67][68][69]. Therefore, in this study, the Landsat 8 images were used to create a time series of LSTs in order to assess the SUHI extension and intensity and identify critical areas for surface temperature. These areas were assimilated to the respective Local Climate Zones and will be deeply analyzed in the following sections.

Classification of Building Roofs
The Worldview3 image was used to identify the kinds of building roof. The high spatial resolution of the image (1.6 m) permits an accurate classification of these surfaces.
First of all, radiometric and atmospheric corrections were applied in order to convert the Digital Number raw data into TOA radiance (Top of the Atmosphere) and then into BOA reflectance (Bottom of the Atmosphere) or surface reflectance. Useful information for image pre-processing was obtained from Kuester et al. [35]. An object-based approach was used for the classification process. At first, the image was divided into significant objects with a multiresolution segmentation, using the shapefiles as support information. Then, the Nearest Neighbor algorithm [70] was trained with ground truth areas (retrieved from the RGB high-resolution images) and applied to the whole image. The building roofs were divided into three classes: -"Clay Tile Roofs", representing the typical coverages of the residential buildings; -"Bright Grey Roofs", as aluminum roofs, metal roofs, etc.; -"Medium/Dark Grey Roofs", as dark bituminous roofs, etc.
One more class was identified as "Shadows". The image was acquired at 10:58 and some tree/building shades are present on the roofs. These shadows could influence the albedo; thus, it is very important to identify and mask these areas.

Surface Albedo
The surface albedo (ALB current ) was computed using the WV3 image and the method suggested by Kuester et al. and Kaplan et al. [35,71]. The surface spectral reflectance of every band (ρ λ ) was multiplied for a weight coefficient (w λ ), calculated from the ratio between the spectral solar irradiance in each measurement band and the sum of the irradiances of all bands [38,55]: The w λ weights coefficients are the same used by Thuilier and suggested in Kuester et al. [35] for WV3.
As already mentioned, the shadows were excluded from the calculation in order to not underestimate the albedo values. A final map of the albedo values for each pixel has been retrieved.

Surface Temperature of Building Roofs
As mentioned above, Landsat 8 temperature maps cannot be used to study building roof surface temperature because of the low spatial resolution [40]. Thus, this surface temperature was computed using the albedo value, the surface's classification, and the meteorological data of the area. From Baldinelli et al. and UNI EN ISO 13,790 [40,41], in steady state conditions, the building's external surface temperature, T e,s , is computed as: where I is the solar irradiance [W/m 2 ], ALB is the surface albedo, T e,a is the external air temperature [K], T i,a is the temperature maintained in the inner environment [K], h e [W/m 2 ·K] is the external heat transfer surface coefficient, and U [W/m 2 ·K] is the thermal transmittance. The third term of the equation accounts for no more than 1 • C in the most common condition [40]; thus, it has been ruled out. The T e,a value and the I value were retrieved from the meteorological station of Reggio Emilia San Lazzaro [51]. The external heat transfer surface coefficient was computed using the model by McAdams [72] and following [73]. In particular, for a rough surface, while for smooth surfaces h e = 3.6V loc + 6.5 where V loc is the local wind speed. For horizontal surfaces and surfaces with slope angle (Φ) in the range 0 • ≤ Φ ≤ 45 • or 135 • ≤ Φ ≤ 180 • , V loc is equal to the wind velocity measured at 10 m (V 10 ). The surface external temperature equation was used to compute the surface temperature of building roofs in the five LCZs analyzed [73]. The hourly averaged values of the 10:00 am CET of 2 June 2017 were T e,a 300.95 K, W 814 W/m 2 , and V 10 4.4 m/s. After the building roof surface temperature calculations for current materials, this study hypothesizes a virtual substitution of all building surfaces with solar reflective materials applied on field. Subsequently, the same Equation (5) was applied to building roofs using the new value of albedo (ALB improved ). This new value depends on the solar reflective material used: cool colored materials were used for clay tile roofs, while white cool materials were used for bright/medium and dark grey roofs. The surface temperature was estimated using Equation (5) for building roofs, with the new albedo value (ALB improved ) and with the aged albedo value (ALB aged ). The latter is included in the study to understand the mitigation effect of solar reflective materials even after a few years from their application. Its value was estimated using the formula in [74]: where ALB improved is the albedo of the new material and β is an index of the durability of the material. It can range between 0.65 and 0.7, depending on whether the coating has been applied on site or with an industrial process. In this study, we consider coatings applied on site.
In the end, two building surface temperature values (T improved and T aged ) were obtained for every roof and compared to the current temperature value (T current ) in order to study the SUHI mitigation effect.

Estimation of UHI Mitigation for LCZs
The new albedo values were used to study the UHI mitigation action in the chosen LCZ. The mean albedo value for each LCZ was computed for the improved scenario and the aged scenario. Applying the equation from Santamouris [75], it was possible to estimate the air temperature decrease for each LCZ related to the albedo difference between the current conditions and the suggested scenarios. The UHI intensity variation is expressed (ATD − Air Temperature Decrease) as: where a is an empirical coefficient of correlation equal to 3.11, and ALBIN is the albedo increase after applying solar reflective materials (ALB improved − ALB current or ALB aged − ALB current ).

SUHI Analysis
The Landsat 8 data were processed, as explained in Section 2.1, to retrieve LST maps for all the images acquired from 2016 to 2021. These maps were analyzed firstly by photointerpretation and by using GIS spatial analysis algorithms. In this way, it was possible to identify several "critical" areas, i.e., areas characterized by high surface temperature values in all the considered images. Figure 4 shows an example of the daily LST maps obtained for the Reggio Emilia municipality.  Figure 4 shows that higher temperature values are generally detected in the city center and in the nearby commercial/industrial areas; in contrast, lower values are measured in the surrounding rural areas. This is especially true for the LST maps computed in June, but less noticeable for the LST maps in August. This different behavior can be attributed mainly to the state of the vegetation: in June, the vegetation is lush and the arable crops in the fields are close to ripening. In August, however, the fields and parks appear as bare soil dried by the high summer temperatures and the minimum rainfall that occurs at these latitudes between mid-to late July and the first weeks of August. As the purpose of the work is to develop mitigation actions on artificial surfaces, it is more efficient to focus on the images of June (for example considering Figure 4) where the vegetation gives its contribution in terms of temperature decrease. Here, it is simple to find that the temperature hotspots are mainly located in the proximity of the city center and the surrounding industrial district (especially in the northern part of the city). The SUHI intensity obviously decreases going from the urban center to the peripheral areas, where the building density is lower.

Local Climate Zones Identification
By adopting the concept of the Local Climate Zone (LCZ), five areas were chosen from the LST maps for the application of surface temperature equation (5) and for the investigation of SUHI and UHI mitigation actions. Critical areas with high surface temperatures were selected together with the rural areas for comparison. In particular, the five chosen areas (corresponding to the five municipal districts of Reggio Emilia) were: 1. Carrozzone; 2. Tondo; 3. Tribunale; 4. Cavazzoli; and 5. Coviolo ( Figure 5).
The first three areas (Carrozzone, Tondo, and Tribunale) are industrial areas, with a strong presence of warehouses, impervious surfaces, streets, parking, etc. These three areas were chosen because they are particularly close to the historic city center and because they are representative of the typical industrial areas of medium-sized municipalities in Italy. Here, many flat bituminous sheath roofs represent hotspots of temperature that can  Figure 4 shows that higher temperature values are generally detected in the city center and in the nearby commercial/industrial areas; in contrast, lower values are measured in the surrounding rural areas. This is especially true for the LST maps computed in June, but less noticeable for the LST maps in August. This different behavior can be attributed mainly to the state of the vegetation: in June, the vegetation is lush and the arable crops in the fields are close to ripening. In August, however, the fields and parks appear as bare soil dried by the high summer temperatures and the minimum rainfall that occurs at these latitudes between mid-to late July and the first weeks of August. As the purpose of the work is to develop mitigation actions on artificial surfaces, it is more efficient to focus on the images of June (for example considering Figure 4) where the vegetation gives its contribution in terms of temperature decrease. Here, it is simple to find that the temperature hotspots are mainly located in the proximity of the city center and the surrounding industrial district (especially in the northern part of the city). The SUHI intensity obviously decreases going from the urban center to the peripheral areas, where the building density is lower.

Local Climate Zones Identification
By adopting the concept of the Local Climate Zone (LCZ), five areas were chosen from the LST maps for the application of surface temperature Equation (5) and for the investigation of SUHI and UHI mitigation actions. Critical areas with high surface temperatures were selected together with the rural areas for comparison. In particular, the five chosen areas (corresponding to the five municipal districts of Reggio Emilia) were: 1. Carrozzone; 2. Tondo; 3. Tribunale; 4. Cavazzoli; and 5. Coviolo ( Figure 5). Atmosphere 2022, 13, x FOR PEER REVIEW 14 of 25 Subsequently, a quantitative assessment was made by computing the average temperatures for each of the 5 LCZs in the LST maps from 2016 to 2021 and reported in Figure  6.   The first three areas (Carrozzone, Tondo, and Tribunale) are industrial areas, with a strong presence of warehouses, impervious surfaces, streets, parking, etc. These three areas were chosen because they are particularly close to the historic city center and because they are representative of the typical industrial areas of medium-sized municipalities in Italy. Here, many flat bituminous sheath roofs represent hotspots of temperature that can be easily improved by cool roofs. While the temperature decrease affects only the considered roofs, the lower heat exchanged with the nearby atmosphere will provide benefits for the whole local city climate. The UHI mitigation actions in these areas could therefore positively affect the residential part of the city. The other two areas (Cavazzoli and Coviolo) are mainly rural areas, with few settlements, numerous crops, and few impervious surfaces.
The LCZs are described in [22,23] through a range of values for several parameters. These values have been computed for the five chosen areas in order to correctly classify them. The first three areas can be attributed to the LCZ 5 "Open midrise", while the last two can be attributed to the LCZ D "Low Plants". The LCZ 5 "Open midrise" is characterized by a building surface fraction between 20 and 40%, an impervious surface fraction between 30 and 50%, and a surface albedo between 0.12 and 0.25 [22]. The LCZ D "Low Plants" is characterized by a building surface fraction lower than 10%, an impervious surface fraction lower than 10%, and a surface albedo between 0.15 and 0.25 [23]. The choice of these areas made it possible to compare industrialized areas with rural areas and to show how the intervention on urban surfaces brings greater advantage in densely built areas. Table 4 shows the main parameters calculated for the five areas. LCZ 5 and LCZ D are the most common LCZs in the medium-sized municipalities, such as Reggio Emilia and a lot of Italian cities. Moreover, another LCZ is really common in the Italian cities and is called "LCZ3-Compact Low Rise". This LCZ could be identified as the typical Italian historical city center. This LCZ was not considered in this work because it is really difficult to replace current materials with solar reflective materials in this kind of area due to city cultural heritage. Figure 5 shows the five chosen LCZs. This LST map was obtained as a mean of all the LST maps from 2016 to 2021. The industrial areas are clearly visible as hotspots of the SUHI; this is true for the municipal districts of Carrozzone, Tondo, and Tribunale, but other hotspots, due to the commercial/industrial areas, are also present in other parts of the city. For example, Figure 4 shows a hotspot in the proximity of the industrial park of Mancasale, north of the city. However, this zone was not selected as a critical area because its municipal district contains several rural areas, and it leads to an average LST lower than the already-considered districts of Carrozzone, Tondo, and Tribunale.
Subsequently, a quantitative assessment was made by computing the average temperatures for each of the 5 LCZs in the LST maps from 2016 to 2021 and reported in Figure 6. Subsequently, a quantitative assessment was made by computing the average temperatures for each of the 5 LCZs in the LST maps from 2016 to 2021 and reported in Figure  6.    Figure 6 confirms the evaluations made by photointerpretation: the industrial areas have higher LST values than the surrounding rural areas. Higher temperatures are retrieved in the summer months, i.e., in June, July, and August, while in May and September the average temperatures are lower. In August 2016, 2017, and 2021, the average temperature of the Coviolo area is close to that of the industrial areas, probably due to the uncultivated agricultural fields full of dry scrub. The acquisition time of the Landsat 8 images, which was in the first half of the morning (10.00 am), must be recalled. This allows the users to have an idea of the temperature trend but not to monitor the peak temperature values of the day.

WV3 Image Classification
The WV3 image, after pre-processing elaboration and conversion into BOA reflectance, was classified using the Nearest Neighbor algorithm to identify three different building roof classes: "Clay Tile Roofs", "Bright Grey Roofs" and "Medium/Dark Grey Roofs". The classification was validated through a dataset of ground truth areas retrieved by the high-resolution RGB images. The overall accuracy of the classification was 90%.
Focusing on the LCZ analysis, 2626 buildings were classified. Table 5 shows the number of buildings of each class for each LCZ. Obviously the three industrial LCZs analyzed include several flat roofs, such as metal roofs and dark bituminous roofs, while the two rural LCZs include a greater number of tile roofs.

Surface Albedo Values and Building Roof Surface Temperature Calculation
Starting from the WV3 image, the surface albedo was then computed. The surface albedo was calculated for each pixel of the BOA reflectance image. The shaded areas were masked in order to avoid underestimations/overestimations of surface albedo. Zonal statistics and support vector files allowed the retrieval of the albedo average value on the surfaces of interest. In particular, the albedo statistics (mean, median, standard deviation, and percentiles) were extracted for each building within the five LCZs analyzed. The average albedo of the entire LCZ was also calculated. Table 6 shows the albedo statistics calculated for the five LCZs and the average albedo only of the buildings in each area. The albedo mean value of Carrozzone was rather higher than the ones in the other LCZs, but it also showed a higher variance. This fluctuation is mainly due to the sun elevation at the acquisition time of the WV3 image, which causes light reflections for some types of roof materials, i.e., metals. Therefore, the median value was considered more significant for representing the albedo of each LCZ. Consequently, observing the median albedo values for each LCZ, it can be seen that the rural areas have slightly higher values than the industrial areas. This low difference is due to the bare fields that maintain low albedo values. The albedo of the building roofs is higher in industrial areas due to the presence of numerous light roofs, such as bright aluminum roofs.
Then, the building roof surface temperature was computed using Equation (5), which requires as inputs the surface albedo and the meteorological data. This equation provides the peak surface temperature, i.e., the maximum surface temperature of the building roofs. In this study, the WV3 acquisition occurred on 2 June 2017 at 10:58 am; thus, a slight underestimation of the daily maximum temperature of the roof may be possible. This does not represent a problem because in the next sections Equation (5) has been used to retrieve the temperature of the roofs replaced with solar reflective materials. In this framework, the calculated temperature difference between the current scenario and the improved one was obtained by setting the same conditions (i.e., air temperature, solar irradiance, and wind conditions at the time of the acquisition). The albedo values of the considered LCZs and their roofs are plotted together with the temperature values in Figure 7. albedo values for each LCZ, it can be seen that the rural areas have slightly higher values than the industrial areas. This low difference is due to the bare fields that maintain low albedo values. The albedo of the building roofs is higher in industrial areas due to the presence of numerous light roofs, such as bright aluminum roofs. Then, the building roof surface temperature was computed using Equation (5), which requires as inputs the surface albedo and the meteorological data. This equation provides the peak surface temperature, i.e., the maximum surface temperature of the building roofs. In this study, the WV3 acquisition occurred on 2 June 2017 at 10:58 am; thus, a slight underestimation of the daily maximum temperature of the roof may be possible. This does not represent a problem because in the next sections equation 5 has been used to retrieve the temperature of the roofs replaced with solar reflective materials. In this framework, the calculated temperature difference between the current scenario and the improved one was obtained by setting the same conditions (i.e., air temperature, solar irradiance, and wind conditions at the time of the acquisition). The albedo values of the considered LCZs and their roofs are plotted together with the temperature values in Figure 7.   Figure 7, it is noticeable that the building surface temperature has comparable values for all the studied areas; slightly higher values were measured in rural areas because of the type of the roof materials (and the lack of bright roofs), but given the low number of buildings, they do not affect the average temperature of the area. In industrial LCZs, the overall temperature is high everywhere due to large impervious surfaces, such as parking lots, roads, etc. The albedo value is not particularly high in both rural and industrial areas. The cause of this value is the presence of numerous fields that do not have active crops and are similar to bare soil and therefore do not contribute to the increase in albedo. However, in rural areas the overall temperature is lower than the industrial areas because of other important phenomena implemented by vegetation, such as evapotranspiration [11].
The analysis then moved to the single building scale and the surface temperature of each roof was retrieved for each of the five LCZs. It is interesting to evaluate the surface albedo for individual buildings and analyze the albedo statistics for the different classes identified on the image. Figure 8 shows the temperature and surface albedo values for the different building classes.   Figure 7, it is noticeable that the building surface temperature has comparable values for all the studied areas; slightly higher values were measured in rural areas because of the type of the roof materials (and the lack of bright roofs), but given the low number of buildings, they do not affect the average temperature of the area. In industrial LCZs, the overall temperature is high everywhere due to large impervious surfaces, such as parking lots, roads, etc. The albedo value is not particularly high in both rural and industrial areas. The cause of this value is the presence of numerous fields that do not have active crops and are similar to bare soil and therefore do not contribute to the increase in albedo. However, in rural areas the overall temperature is lower than the industrial areas because of other important phenomena implemented by vegetation, such as evapotranspiration [11].
The analysis then moved to the single building scale and the surface temperature of each roof was retrieved for each of the five LCZs. It is interesting to evaluate the surface albedo for individual buildings and analyze the albedo statistics for the different classes identified on the image. Figure 8 shows the temperature and surface albedo values for the different building classes. Atmosphere 2022, 13, x FOR PEER REVIEW 17 of 25 The roofs classified as "Medium/Dark Grey Roofs" show, as is obvious, lower albedo values compared to other kinds of surfaces. On the other hand, the "Bright Grey Roofs" present the highest value of albedo and, as a consequence, the lowest value of temperature. Clay tile roofs and medium/dark grey roofs show similar values of albedo and surface temperature, but the last category can be significantly improved using solar reflective materials.

Solar Reflective Materials and SUHI Mitigation Actions
The last step of the analysis concerns the hypothesis of virtually replacing the current roof materials with solar reflective materials. The typical albedo values of these materials were found in the scientific literature. The values used for the different kinds of surfaces are presented in Table 7. Here, the aged albedo calculated by equation (8) is also shown, considering the application of the coating on field.  [26,[81][82][83] The surface albedo can be greatly increased only for some types of roofing (e.g., flat roofs or industrial roofs). For pitched tile roofs, the application of colored solar reflective material allows one to obtain albedo values lower than the white solar reflective material [77,79]. The current roof albedo values (ALBcurrent) were then replaced by the new and aged albedo values (ALBimproved and ALBaged), depending on the building's classification. Two new surface temperature maps were therefore obtained for the building roofs: one represents the scenario that immediately follows the application of solar reflective materials (Timproved) and the other one represents the aged scenario (Taged). The new temperature maps were computed using equation (5). Thus, two new values of the albedo were obtained for every LCZ: the first one referred to the improved scenario, and the second one referred to the The roofs classified as "Medium/Dark Grey Roofs" show, as is obvious, lower albedo values compared to other kinds of surfaces. On the other hand, the "Bright Grey Roofs" present the highest value of albedo and, as a consequence, the lowest value of temperature. Clay tile roofs and medium/dark grey roofs show similar values of albedo and surface temperature, but the last category can be significantly improved using solar reflective materials.

Solar Reflective Materials and SUHI Mitigation Actions
The last step of the analysis concerns the hypothesis of virtually replacing the current roof materials with solar reflective materials. The typical albedo values of these materials were found in the scientific literature. The values used for the different kinds of surfaces are presented in Table 7. Here, the aged albedo calculated by Equation (8) is also shown, considering the application of the coating on field.  [26,[81][82][83] The surface albedo can be greatly increased only for some types of roofing (e.g., flat roofs or industrial roofs). For pitched tile roofs, the application of colored solar reflective material allows one to obtain albedo values lower than the white solar reflective material [77,79]. The current roof albedo values (ALB current ) were then replaced by the new and aged albedo values (ALB improved and ALB aged ), depending on the building's classification. Two new surface temperature maps were therefore obtained for the building roofs: one represents the scenario that immediately follows the application of solar reflective materials (T improved ) and the other one represents the aged scenario (T aged ). The new temperature maps were computed using Equation (5). Thus, two new values of the albedo were obtained for every LCZ: the first one referred to the improved scenario, and the second one referred to the aged one. Even if it is not possible to analyze the LST variation of the entire LCZ, it is possible to make assessments on the building roof temperature changes in the different scenarios. Table 8 shows the albedo values for the current, improved, and aged scenario of the whole LCZ. Table 8 also presents the mean of the building roof temperatures for the current, improved, and aged scenarios. Compared to the current state, the albedo value grows more in the LCZs with more buildings and therefore in the industrial ones. In these areas the "improved" scenario allows an increase in the albedo from 0.125 up to 0.271, while the aged scenario is around 0.214. The rural areas Cavazzoli and Coviolo, which have a 3% percentage of buildings, have smaller albedo increases than the industrial LCZs as the influence of the buildings is minimal. Figure 9 shows the surface temperature of the building roofs for a portion of the Carrozzone area as an example. These maps were obtained using Equation (5) with the albedo and meteorological parameters of the current scenario (b), with the solar reflective materials improved values (c), and in the aged scenario (d).
Atmosphere 2022, 13, x FOR PEER REVIEW 18 of 25 aged one. Even if it is not possible to analyze the LST variation of the entire LCZ, it is possible to make assessments on the building roof temperature changes in the different scenarios. Table 8 shows the albedo values for the current, improved, and aged scenario of the whole LCZ. Table 8 also presents the mean of the building roof temperatures for the current, improved, and aged scenarios. 2 Compared to the current state, the albedo value grows more in the LCZs with more buildings and therefore in the industrial ones. In these areas the "improved" scenario allows an increase in the albedo from 0.125 up to 0.271, while the aged scenario is around 0.214. The rural areas Cavazzoli and Coviolo, which have a 3% percentage of buildings, have smaller albedo increases than the industrial LCZs as the influence of the buildings is minimal. Figure 9 shows the surface temperature of the building roofs for a portion of the Carrozzone area as an example. These maps were obtained using equation (5) with the albedo and meteorological parameters of the current scenario (b), with the solar reflective materials improved values (c), and in the aged scenario (d).  Observing Figure 9, the decrease in roof surface temperatures is clear both in the "improved" scenario and the "aged" scenario. Industrial buildings roofs take full advantage of the application of white solar reflective materials, showing noticeable drops in surface temperatures. Figure 10 reports for each LCZ: • the current values of albedo and surface temperature of building roofs; • the values of albedo and surface temperature of building roofs after the application on field of solar reflective materials ("improved" scenario); • the values of albedo and surface temperature of building roofs after a few years from the application of solar reflective materials ("aged" scenario).
the application of solar reflective materials ("aged" scenario).
The industrial areas present the highest gains of albedo value because of the massive presence of buildings and because of the kind of roofs that are flat and suitable for the application of white solar reflective materials. Rural areas include mainly clay tile roofs, on which only the application of colored solar reflective material is possible. Here, the albedo increase is limited. For all LCZs, the aging of the surface solar reflective materials leads to a loss of albedo capacity, but anyway, the "aged" scenario is still better than the current one.

UHI Mitigation
The Surface Urban Heat Island is surely related to the Urban Heat Island and actions to mitigate the first one also influence the second one [15,17]. Equation (9) was used to estimate the temperature decrease for the two scenarios considered (improved and aged), as compared to the current situation. This formula only provides an indicative temperature variation that has to be applied to current air temperature values. Thus, it provides a measure of the UHI loss of intensity. However, the equation has some limitations because it is not site specific [84,85]. Table 9 shows the albedo increase (ALBIN) obtained for the improved and aged scenario and the consequent air temperature decrease (ATD) that has to be applied to the UHI. The ALBIN is retrieved by subtracting the current value of albedo for each LCZ from the albedo value in the improved or aged scenario. The application of solar reflective materials improves the albedo of the LCZ by 0.146 for Carrozzone, 0.122 for Tondo, 0.1 for Tribunale, 0.013 for Cavazzoli, and 0.014 for Coviolo. Consequently, roof surface temperature decreases by 33.3% for Carrozzone, 28.6% for Tondo, 28.4% for Tribunale, 28.1% for Cavazzoli, and 28.3% for Coviolo. Aged performances compared to the current situation show limited benefits in terms of albedo (improvement of 0.09 for Carrozzone, 0.08 for Tondo, 0.07 for Tribunale, 0.008 for Cavazzoli, and 0.01 for Carrozzone) and cooling effects (temperature decrease of 20.7% for Carrozzone, 19% for Tondo, 18.8% for Tribunale, 28.1% for Cavazzoli, and 28.3% for Coviolo).
The industrial areas present the highest gains of albedo value because of the massive presence of buildings and because of the kind of roofs that are flat and suitable for the application of white solar reflective materials. Rural areas include mainly clay tile roofs, on which only the application of colored solar reflective material is possible. Here, the albedo increase is limited. For all LCZs, the aging of the surface solar reflective materials leads to a loss of albedo capacity, but anyway, the "aged" scenario is still better than the current one.

UHI Mitigation
The Surface Urban Heat Island is surely related to the Urban Heat Island and actions to mitigate the first one also influence the second one [15,17]. Equation (9) was used to estimate the temperature decrease for the two scenarios considered (improved and aged), as compared to the current situation. This formula only provides an indicative temperature variation that has to be applied to current air temperature values. Thus, it provides a measure of the UHI loss of intensity. However, the equation has some limitations because it is not site specific [84,85]. Table 9 shows the albedo increase (ALBIN) obtained for the improved and aged scenario and the consequent air temperature decrease (ATD) that has to be applied to the UHI. The ALBIN is retrieved by subtracting the current value of albedo for each LCZ from the albedo value in the improved or aged scenario. The use of solar reflective materials on building roofs brings a temperature decrease of up to 0.5 • C for the Carrozzone area, i.e., the area with the largest number of buildings. Slightly lower results are obtained for the other industrial areas. For the "aged" scenario, on the other hand, it is possible to obtain a reduction in the intensity of the UHI of up to 0.3 • C. For rural areas, the air temperature decrease is practically irrelevant given the low number of buildings. Solar reflective materials are therefore confirmed as suitable for UHI mitigation in the most densely built areas with dark, flat roofs [75,86,87]. This analysis has been focused only on building roofs, but in the future, it will also be extended to the other impervious surfaces present within the LCZs which can be replaced with solar reflective materials (for example cool pavements applied to streets, parking, etc. [81,[88][89][90]. Furthermore, new LCZs will be considered to evaluate the impacts of a wide application of solar reflective materials in urban areas.

Conclusions
This study considered the use of satellite images for the identification of SUHIs, the classification of building surfaces, and the calculation of surface albedo.
Landsat 8 images from 2016 to 2021 were processed to obtain Land Surface Temperature maps. These maps show the extension of the SUHI that usually decreases going from the urban center to the peripheral areas where building density is lower. The SUHI phenomenon is surely more noticeable when the vegetation is lush, and the presence of bare soil is limited.
Critical areas for the SUHI were first identified and categorized with the Local Climate Zones approach to studying local mitigation action on the SUHI. A Worldview 3 image allowed us to classify building roofs and to compute surface albedo values. Subsequently, empirical formulas were used to calculate the surface temperature of the building roofs, starting from the albedo value and the meteorological data of the study area.
Eventually, two mitigation scenarios of the SUHI, and consequently of the UHI, were constructed, assuming the replacement of building roofs with solar reflective materials that have a higher albedo value. Depending on the kind of roof, different materials were chosen: white solar reflective materials with an albedo value of up to 0.9 are suitable for industrial roofs, while colored solar reflective materials with an albedo of up to 0.55 could be used for clay tile roofs. The first scenario, called "improved", is the one immediately after the application of the solar reflective materials on field. After a few years from application, however, the albedo of these materials decreases due to atmospheric agents, dust deposition, or even biological growth; therefore, the "aged" scenario is also considered in this study. The industrial areas present the highest gains of albedo value (with an increase of up to 0.146 in the first scenario and 0.09 after aging for the LCZ of Carrozzone) because of the massive presence of buildings and because of the kind of roofs that are flat and suitable for the application of white solar reflective materials. Here, the benefits in terms of roof temperature decrease can reach 33.3% compared to the current situation. Rural areas include mainly clay tile roofs on which only the application of colored solar reflective material is possible.
Regarding the UHI mitigation effect, the "improved" scenario shows a significant increase in the albedo in densely built industrial areas and an air temperature decrease