Surface Urban Heat Island in Middle City: Spatial and Temporal Characteristics

: Currently, cities have aroused the interest of researchers due the local climate change caused by the surface urban heat island (SUHI) e ﬀ ect. The impact of anthropogenic land use and cover changes has led to more frequent intense SUHI, with direct consequences on urban quality of life. Therefore, this research aims at analyzing the inﬂuences of natural and anthropogenic variables on the seasonality and spatial SUHI intensity in a Brazilian city, using remote sensing data and analysis of several physical parameters. Results show that the city of S ã o Carlos has an SUHI mosaic and surface urban freshness island (SUFI). On average, 86% of the urban area presented a SUHI, whilst most SUFIs are located near watercourses, parks, slopes and valley bottoms, revealing the e ﬀ ects of green areas and relief on creation of microclimates. The SUHI showed signiﬁcant seasonal variability.


Introduction
Currently, cities have increasingly aroused the interest of researchers for the quality of life in the urban environment. Concerns about climatic conditions of metropolises have led several scholars to conduct research on this theme [1][2][3][4][5][6][7][8][9].
Studies conducted in the urban atmosphere have shown that cities generate a climate of their own, as result of the interference of various factors that occur within urban boundaries, that alters the climate at a local or regional scale [10]. In cities, changes in land surface temperatures are mainly due to building concentration, waterproofing by the pavement of streets and avenues, increases in the flow of people and automobiles and decreases in the vegetation cover. Such changes reverberate both local and regional climatic effects, depending on the size of the respective urban area [6,[11][12][13].
The most prominent feature of urban climate is the urban heat island (UHI) phenomenon. The surface urban heat island (SUHI) [14] is a type of UHI and occurs over the surface envelope mediated by the changes in the surface temperature between urban and rural areas [15][16][17]. The UHI has been assessed by using meteorological observations from the mobile transects, near surface networks [7] [18,19], and climate classifications [20]. With the advance of new geospatial technologies, researchers have increasingly used remote sensing data (e.g., satellite imagery) to study SUHI [6,[21][22][23]. Satellite remote sensing covers large areas and provides a direct and consistent way to examine the thermal environment of cities large [7,8,[24][25][26].
The SUHI causes high summer temperatures and leads to serious health problems associated with heat, such as human thermal discomfort [27][28][29][30]. Therefore, mitigating SUHIs is critical for the management of urban space. The spatialization of differences between intra-urban and rural temperatures may offer subsidies to combat the magnitude of heat islands [31]. The SUHI simulations have been used in numerous studies [32][33][34][35], so that such simulations may provide prospects for reordering a territory compared with real data or based on observations of other studies, The SUHI is largely influenced by a number of natural and anthropogenic factors, which includes city size and population impermeable surfaces (low albedo, decreased evapotranspiration), topography, anthropogenic heat and weather conditions [13,[35][36][37][38]. Quantifying SUHI intensity and their independent variables is very important to improve the quality of life and achieve sustainable development of the city [22]. While most studies on UHI and SUHIs focus on large cities, the overall objective of this research was to investigate the SUHI intensity in a southeastern Brazilian city. We aim to (i) verify the seasonal and spatial variability of SUHI and (ii) discuss the relations between natural and anthropogenic variables that favor the development of the SUHI.

Study Area
The municipality of São Carlos, located in the state of São Paulo (Figure 1), is a transition region between the Cerrado biome and the Atlantic Forest. The climate is controlled by equatorial, tropical and polar masses [39]. Seasonality is well defined, with dry periods from April to September and humid periods from October to March. The monthly average temperature is 20.5 • C, the highest monthly average temperature occurs in February (23 • C) and January (22.9 • C) and the lowest average temperature occurs in July (16.9 • C) and June (17 • C) [40]. Rainfalls have a complex inter-annual variability. The average rainfall is 1528.1 mm, with a standard deviation of 288.8 mm. However, there may be exceptional rainfalls: the minimum rainfall recorded is 1143.4 and the maximum is 2486.7 mm [41].
Urban Sci. 2020, 4, x FOR PEER REVIEW 2 of 19 management of urban space. The spatialization of differences between intra-urban and rural temperatures may offer subsidies to combat the magnitude of heat islands [31]. The SUHI simulations have been used in numerous studies [32][33][34][35], so that such simulations may provide prospects for reordering a territory compared with real data or based on observations of other studies, The SUHI is largely influenced by a number of natural and anthropogenic factors, which includes city size and population impermeable surfaces (low albedo, decreased evapotranspiration), topography, anthropogenic heat and weather conditions [13,[35][36][37][38]. Quantifying SUHI intensity and their independent variables is very important to improve the quality of life and achieve sustainable development of the city [22]. While most studies on UHI and SUHIs focus on large cities, the overall objective of this research was to investigate the SUHI intensity in a southeastern Brazilian city. We aim to (i) verify the seasonal and spatial variability of SUHI and (ii) discuss the relations between natural and anthropogenic variables that favor the development of the SUHI.

Study Area
The municipality of São Carlos, located in the state of São Paulo (Figure 1), is a transition region between the Cerrado biome and the Atlantic Forest. The climate is controlled by equatorial, tropical and polar masses [39]. Seasonality is well defined, with dry periods from April to September and humid periods from October to March. The monthly average temperature is 20.5 °C, the highest monthly average temperature occurs in February (23 °C) and January (22.9 °C) and the lowest average temperature occurs in July (16.9 °C) and June (17 °C) [40]. Rainfalls have a complex inter-annual variability. The average rainfall is 1528.1 mm, with a standard deviation of 288.8 mm. However, there may be exceptional rainfalls: the minimum rainfall recorded is 1143.4 and the maximum is 2486.7 mm [41]. São Carlos is geomorphologically located on the western edge of the Basaltic Cuestas, near the Western Plateau and the Peripheral Depression. Thus, the relief of the city comprises isolated tables and hills, which are the result of an intense erosion caused by rivers that deeply carved the plateau. São Carlos is geomorphologically located on the western edge of the Basaltic Cuestas, near the Western Plateau and the Peripheral Depression. Thus, the relief of the city comprises isolated tables and hills, which are the result of an intense erosion caused by rivers that deeply carved the plateau. The altitude has a range of 480 m. The lowest altitudes are around 775 m in valleys to the west of the main urban area, whilst the highest altitudes of 910 m are located to the east in a plateau [42].
The soils of the São Carlos are highly correlated with the relief and lithology shapes. There are predominantly clayey soils at the highest levels; shallow soils and rocky outcrops on cliffs; and medium textured soils (clay content between 15% and 35%) in lower altitudes areas. In the latter, there are also extensive areas of sandy and deep soils. Therefore, the city is characterized by high permeability soils consisting mostly of Red-Yellow Latosol and deep quartz sand with sandy characteristics [42].
The original vegetation of the region is the Cerrado, characterized by a herbaceous-woody and herbaceous-shrub forest formation with evergreen trees [42]. This type of vegetation is associated with an almost flat topographic sandy terrain extending to the north to the banks of the Mogi-Guaçu River and to the south to the edges of the Jacaré-Guaçu River valley. However, there is currently little original vegetation due to most non-urbanized areas being occupied by pastures, plantations, reforestation and mining [42].
In 2018, the population of São Carlos was estimated at 249,415 inhabitants [43]. However, the urban drainage, the road system and leisure areas did not keep up with the city's expansion, as 98.4% of the households have adequate sanitation, 89.2% of urban households are on public roads with forestation, and 34.2% of urban households are on public roads with an adequate urbanization (manholes, sidewalks, pavement and curbs). The city has also expanded without taking into account environmental, climatic, geological and topographical factors [44,45], leading to numerous social and environmental problems, such as floods, landslides, disorderly construction and high temperatures associated with UHI. Figure 2 shows succinctly the steps to be performed in this research. This methodological procedure was elaborated based on several researches, seeking the most appropriate method for each step. The altitude has a range of 480 m. The lowest altitudes are around 775 m in valleys to the west of the main urban area, whilst the highest altitudes of 910 m are located to the east in a plateau [42]. The soils of the São Carlos are highly correlated with the relief and lithology shapes. There are predominantly clayey soils at the highest levels; shallow soils and rocky outcrops on cliffs; and medium textured soils (clay content between 15% and 35%) in lower altitudes areas. In the latter, there are also extensive areas of sandy and deep soils. Therefore, the city is characterized by high permeability soils consisting mostly of Red-Yellow Latosol and deep quartz sand with sandy characteristics [42].

Materials and Methods
The original vegetation of the region is the Cerrado, characterized by a herbaceous-woody and herbaceous-shrub forest formation with evergreen trees [42]. This type of vegetation is associated with an almost flat topographic sandy terrain extending to the north to the banks of the Mogi-Guaçu River and to the south to the edges of the Jacaré-Guaçu River valley. However, there is currently little original vegetation due to most non-urbanized areas being occupied by pastures, plantations, reforestation and mining [42].
In 2018, the population of São Carlos was estimated at 249,415 inhabitants [43]. However, the urban drainage, the road system and leisure areas did not keep up with the city's expansion, as 98.4% of the households have adequate sanitation, 89.2% of urban households are on public roads with forestation, and 34.2% of urban households are on public roads with an adequate urbanization (manholes, sidewalks, pavement and curbs). The city has also expanded without taking into account environmental, climatic, geological and topographical factors [44,45], leading to numerous social and environmental problems, such as floods, landslides, disorderly construction and high temperatures associated with UHI. Figure 2 shows succinctly the steps to be performed in this research. This methodological procedure was elaborated based on several researches, seeking the most appropriate method for each step.  (Tables 1 and 2). All available images had reduced cloud cover with a maximum of 10.45% (Table 1). The images were generated at the same time with a discrepancy of less than two minutes, which ensures comparison over the years. All images had excellent quality, with a value of nine on a scale from zero to ten, so that zero means poor quality and ten means excellent  (Tables 1 and 2). All available images had reduced cloud cover with a maximum of 10.45% (Table 1). The images were generated at the same time with a discrepancy of less than two minutes, which ensures comparison over the years. All images had excellent quality, with a value of nine on a scale from zero to ten, so that zero means poor quality and ten means excellent quality. The weather conditions for each image acquisition campaign are observed in Table 2. Relative humidity was high (above 86%) on 23 April 2015 and 14 July 2016, while the average air temperature was below 30 • C on all occasions. The highest wind speed was 7.23 m/s.

Dependent Variable
Procedures for Calculation of the SUHI LANDSAT data are freely available for download on the website of the United States Geological Survey (USGS) [46]. LANDSAT satellite sensors have been monitoring the Earth for over four decades, providing data continuity during this period [47,48]. The first series of this satellite was launched in 1972, named LANDSAT 1. Since then, there have been eight LANDSAT satellites. Among all eight, LANDSAT 6 failed to reach orbit and fell to the Earth in 1993. The other satellites were successful, and provided the researchers with a large amount of data used in many studies [47,49].
The calculation of the SUHI was performed in six steps ( Figure 2), as follows: (i) was to convert the digital numbers into radiance; (ii) to calculate brightness temperature; (iii) to calculate the NDVI (Normalized Difference Vegetation Index); (iv) to estimate the surface emissivity using the values Urban Sci. 2020, 4, 54 5 of 19 obtained by the NDVI; (v) after estimating emissivity, to determine the corrected surface temperature, which was later transformed into SUHI in the sixth step, as suggested by Ndossi and Avdan [47,48].
Step 1: Digital Numbers (DN) for Radiance Thermal data from LANDSAT satellite sensor images are stored in digital numbers (DN). DNs are pixels that have not yet been calibrated and represent different levels of radiance in the raster image. After obtaining satellite images, the first step converted DN into spectral radiance using data from the LANDSAT 8 satellite TIRS sensor by applying the Equation (1) [50]: where L λ is the spectral radiance at W/ m 2 × sr × µm at the top of the atmosphere. M L is the multiplicative scaling factor of band 10. Q cal is the quantized value calibrated by the pixel in DN. O i is the USGS published corrections for TIRS band calibration.
where BT is the brightness temperature in (K), L λ is the spectral radiance at the top of the atmosphere, K 1 is a constant value for conversion to thermal band, as is K 2 .
Step 3: NDVI Calculation Following the above calculations, the NDVI was calculated by obtaining the ratio between the difference in reflectance measurement in the near (ρ 5 ) and red (ρ 4 ) infrared channels, and the sum of these channels [51] (Equation (3)): NDVI is an indicator of the quantity and condition of green vegetation. The values range from −1 to +1. Positive values are obtained some type of vegetation areas and vary according to the photosynthetic activity of the plants and its density. In the water surfaces, the NDVI is generally lower than zero [23,52,53].
Step 4: Determination of Earth's Surface Emissivity Terrestrial surface emissivity was estimated using the NDVI values. According to Zhang, Wang and Li [54], when the NDVI of an area is known, emissivity can be estimated. Thus, the emissivity of a pixel was estimated by classifying the pixels according to the class into which they fit (Table 3). When a pixel obtained a NDVI value below −0.185, the emissivity value was 0.955 and the NDVI value was higher than or equal to −0.185 and lower than 0.157, the emissivity of 0.985 was attributed to the pixel. When NDVI was higher or equal to 0.157 and lower than or equal to 0.727, a logarithmic relation between NDVI and emissivity was used [47]. Finally, when the NDVI was higher than 0.727, a value of 0.990 was assigned to the pixel (Table 3). Table 3. NDVI (Normalized Difference Vegetation Index)-based algorithm [47].

NDVI
Land Surface Emissivity (LSE) Step 5: Surface Temperature Correction After estimating the emissivity of the soil surface, the brightness temperature correction was performed, using the Planck's function, that corrects the emissivity of a surface compared to a black body [55,56], according to the Equation (4): and ε is the surface emissivity.
Step 6: Surface Temperature Transformation into SUHI Usually, the intensity of the urban heat island (UHII) is measured using air temperature observations along transects or at fixed stations in the surrounding urban and rural areas [22,57]. However, the SUHI requires another meaning because data for it are obtained by remote sensing [23]. Thus, the SUHI was calculated as the difference in surface temperature of each pixel in the urban area by the average surface temperature of the rural area [23,58,59] according to Equation (5): where BT U is the surface temperature in each pixel of the urban area, and BT R is the average of the rural surface temperature obtained from a 1-km buffer from the urban boundaries ( Figure 3). A positive value means the occurrence of a SUHI and a negative value represents surface urban freshness island (SUFI). There is no standardized method for delimiting the rural area buffer and measuring SUHIs. Alves [23] analyzed the seasonal and spatial variation of SUHI intensity in a small urban cluster in Brazil using a 500-m buffer; then he used 200-m and 500-m buffers to quantify the effects of watercourses and green areas on surface temperature. Cheval

Independent Variables
The independent variables altitude and topographic position index were calculated using SRTM data. The variables UI, NDVI (obtained in Step 3) and NDWI were obtained using images of the LANDSAT 8 satellite.

Altitude
The altitude or height above mean sea level is a key factor that controls air temperature. Areas with different altitudes may cause variations in air temperature and atmospheric emissivity, which may cause variations in surface temperature [62]. We used data from the TOPODATA project [63], which offers free access to geomorphometric variables derived from SRTM data (Shuttle Radar Topographic Mission) for the entire Brazilian territory. Subsequently, these data were refined from the original spatial resolution from 3 arc-seconds (~90 m) to 1 arc-second (~30 m) per kriging.

Topographical Position Index
The Topographic Position Index (TPI), an algorithm created by Weiss [64], which calculates differences in the elevation between a central pixel (Z0) and the average elevation of its surroundings (Zm). The value of TPI expresses the contrast intensity [64,65].
The TPI allows for determining characteristics of the landscape taking into account the topographic attributes. Thus, the units defined are closely related to water flow conditions and surface dynamics [64].
The TPI is the difference between a cell elevation value and the average neighborhood altitude around these cells. Positive values mean that the cell is higher than its surroundings, while negative values mean it is lower. If it is significantly higher than the surrounding area, then it is likely to be near or on top of a hill or a ridge. Significantly low values suggest that the cell is near or at the bottom of a valley [64]. TPI values close to zero may mean either a flat area or a medium slope area.
The variation in the inclination and orientation of relief slopes, associated with variations in altitude, plays an important role in the distribution of calorific-luminous energy [62,66]. The radiative flow of an inclined slope may be more intense than others at a certain time of day, but the incidence of sunlight may be significantly lower in other periods. Likewise, the shading of buildings and trees, among others, are more expressive in flat areas [62,65]. For modeling, we used data from 1000 random points within the buffers. This reduction in sampling was necessary due to the complexity of GWR calculation, which requires high memory and computational processing.

Independent Variables
The independent variables altitude and topographic position index were calculated using SRTM data. The variables UI, NDVI (obtained in Step 3) and NDWI were obtained using images of the LANDSAT 8 satellite.

Altitude
The altitude or height above mean sea level is a key factor that controls air temperature. Areas with different altitudes may cause variations in air temperature and atmospheric emissivity, which may cause variations in surface temperature [62]. We used data from the TOPODATA project [63], which offers free access to geomorphometric variables derived from SRTM data (Shuttle Radar Topographic Mission) for the entire Brazilian territory. Subsequently, these data were refined from the original spatial resolution from 3 arc-seconds (~90 m) to 1 arc-second (~30 m) per kriging.

Topographical Position Index
The Topographic Position Index (TPI), an algorithm created by Weiss [64], which calculates differences in the elevation between a central pixel (Z0) and the average elevation of its surroundings (Zm). The value of TPI expresses the contrast intensity [64,65].
The TPI allows for determining characteristics of the landscape taking into account the topographic attributes. Thus, the units defined are closely related to water flow conditions and surface dynamics [64].
The TPI is the difference between a cell elevation value and the average neighborhood altitude around these cells. Positive values mean that the cell is higher than its surroundings, while negative values mean it is lower. If it is significantly higher than the surrounding area, then it is likely to be near or on top of a hill or a ridge. Significantly low values suggest that the cell is near or at the bottom of a valley [64]. TPI values close to zero may mean either a flat area or a medium slope area.
The variation in the inclination and orientation of relief slopes, associated with variations in altitude, plays an important role in the distribution of calorific-luminous energy [62,66]. The radiative flow of an inclined slope may be more intense than others at a certain time of day, but the incidence of sunlight may be significantly lower in other periods. Likewise, the shading of buildings and trees, among others, are more expressive in flat areas [62,65].

Urbanization Index
The Urbanization Index (UI) is calculated using the bands 5 and 7 of the LANDSAT 8 satellite. This index was proposed by Kawamura, Jayamanna and Tsujiko [67] to estimate urbanization. It may vary from −1 to 1: the higher the value, the greater the degree of urbanization, following the Equation (6): where UI is the Urbanization index, ρ 7 is the reflectance of the band 7 and ρ 5 is the reflectance of the band 5 of the LANDSAT 8 satellite.

Normalized Difference Water Index
The Normalized Difference Water Index (NDWI) was developed by McFeeters [68], who purposes to delineate features in the aquatic environment and highlight their presence in the images. Water is important for the thermal balance because of its ability to absorb energy and release it slowly. The NDWI can be defined according to Equation (7): where ρ 3 is the reflectance of the band 3 and ρ 5 is the reflectance of the band 5 of the LANDSAT 8 satellite. NDWI values range from −1 to 1. The value zero is a threshold, i.e., the type of cover is considered water if NDWI ≥ 0 and not water if NDWI ≤ 0.

SUHI in São Carlos
Cities have a complex thermal system. Surface varieties are largely responsible for creating a mosaic of microclimates. In addition to the artificialization of the urban environment, natural issues, such as topography, also establish a direct relationship in the creation of microclimates [6,9,12,35,69,70].
It is noteworthy that the mean value of any phenomenon should be carefully considered; it also masks some episodic specificities and makes the temporal patterns more evident. In the city of São Carlos, houses, buildings, pavements, roofs, streams, riparian forests and parks create a mosaic of microclimates influenced by both natural and anthropogenic characteristics ( Figure 4).
The areas corresponding to Miguel Petroni Street (2) (14) presented, on average, the greatest SUHIs (4 • C). The areas 2, 5, 7, 10, 12 and 13 ( Figure 4) consist of residential buildings, mostly with ceramic roofs, while the areas 3 and 8 correspond to industrial areas with fiber cement roofs. There is a shopping center in the area 15, which is a shopping area with ample parking space with an asphalt pavement. Therefore, in São Carlos, the greatest SUHIs were observed in areas built with materials of a high thermal conductivity and heat capacity (concrete and asphalt). According to Gartland [13], materials with a high thermal conductivity are better able to direct the absorbed thermal energy into their interior, and materials with a high heat capacity can store more energy without losing it to the environment: the more energy is stored, the greater the material temperature will be. areas on the surface temperature, and consequently on SUHIs. [9,69,72] found a similar effect exerted by green areas. Urban parks have been considered an important part of urban vegetation because they are cooler than their surrounding areas, creating SUFI effect [57,73]. There was a predominance of SUHI over SUFI in the urban area of São Carlos ( Figure 5). The SUHI, on average, occupied an area of 52.6 km 2 , which means six times the area occupied by the SUFI (8.5 km 2 ). This finding highlights a problem observed in several Brazilian cities: the lack of urban planning considering microclimate as a factor for spatial planning [19,23,[74][75][76]. Most SUFIs are located near watercourses, parks, slopes and valley bottoms ( Figure 5), revealing the effects of green areas and relief on SUHI intensity. A similar result was observed by Alves [77] and Lopes [78]. However, in urban areas, some of them may have surface urban freshness islands (SUFI), also called negative SUHIs [23,71], due to the presence of green areas (parks, woods, riparian forest etc.) or the presence of surfaces with water. In São Carlos, the green areas located between the Quebec and Fehr Parks (1), São Carlos Ecological Park (4), a remnant green area west of Hermano Bomfim da Silva Street (6), riparian forests west of Vila Conceição (11), Bicão Park (14) and the riparian forest east of the University of São Paulo-campus 2 (16) are examples of the influence of green and humid areas on the surface temperature, and consequently on SUHIs. Refs. [9,69,72] found a similar effect exerted by green areas. Urban parks have been considered an important part of urban vegetation because they are cooler than their surrounding areas, creating SUFI effect [57,73].
There was a predominance of SUHI over SUFI in the urban area of São Carlos ( Figure 5). The SUHI, on average, occupied an area of 52.6 km 2 , which means six times the area occupied by the SUFI (8.5 km 2 ). This finding highlights a problem observed in several Brazilian cities: the lack of urban planning considering microclimate as a factor for spatial planning [19,23,[74][75][76]. Most SUFIs are located near watercourses, parks, slopes and valley bottoms ( Figure 5), revealing the effects of green areas and relief on SUHI intensity. A similar result was observed by Alves [77] and Lopes [78]. Urban Sci. 2020, 4, x FOR PEER REVIEW 10 of 19 Figure 5. Area occupied by SUHI and surface urban freshness island (SUFI).
In the city of São Carlos, the seasonality of SUHIs was evident ( Figure 6). Although, the development of UHIs and SUHIs is generally associated with urban expansion, increased human and industrial activities, geographical location, climate, urban materials, no-presence of vegetation, seasons and times of day and night may exert different effects [80]. Lazzarini et al. [81] reported the SUHI in Abu Dhabi, a city in the desert, and showed a reversal of what is considered a standard for SUHIs during the day: central areas had lower surface temperatures compared to the suburbs, with a daily difference of 5-6 K in summer and 2-3 K in winter. In contrast, the highest intensities of negative SUHIs were verified within the urban area during the winter in São Carlos.
The seasonal variation of SUHIs ( Figure 6) results from changes in solar radiation intensity due to the earth's translational movement and the variation of Cerrado and pasture leaf areas, which follow the rhythm of rainfalls [23,82]. The highest SUHI intensities were observed in the summer, with SUHIs above 8 °C. According to Buyantuyev and Wu [79], surface urban heat islands (SUHIs) exist at any time of the day. However, they are more intense closer to noon and during the summer.
The summer is the rainy season in São Carlos, where vegetation and pasture become denser due to a higher metabolic activity, as emphasized by the NDVI values. As a result, the urban and rural contrast is sharper, and the SUHI intensities tend to be higher. In late autumn and winter, there was little difference in the SUHIs intensity from urban to rural areas. In autumn, the amplitude of SUHI intensity was lower than in other seasons, ranging from −2.6 to 4.4 °C because, during this period, the Cerrado and pastures dry or lose their leaves to reduce the adverse effects of the lack of rainfalls. Besides that, green areas are known to decrease surface temperature due to shading and evaporative cooling [83,84].
In the city of São Carlos, the seasonality of SUHIs was evident ( Figure 6). Although, the development of UHIs and SUHIs is generally associated with urban expansion, increased human and industrial activities, geographical location, climate, urban materials, no-presence of vegetation, seasons and times of day and night may exert different effects [80]. Lazzarini et al. [81] reported the SUHI in Abu Dhabi, a city in the desert, and showed a reversal of what is considered a standard for SUHIs during the day: central areas had lower surface temperatures compared to the suburbs, with a daily difference of 5-6 K in summer and 2-3 K in winter. In contrast, the highest intensities of negative SUHIs were verified within the urban area during the winter in São Carlos.
The seasonal variation of SUHIs ( Figure 6) results from changes in solar radiation intensity due to the earth's translational movement and the variation of Cerrado and pasture leaf areas, which follow the rhythm of rainfalls [23,82]. The highest SUHI intensities were observed in the summer, with SUHIs above 8 • C. According to Buyantuyev and Wu [79], surface urban heat islands (SUHIs) exist at any time of the day. However, they are more intense closer to noon and during the summer.
The summer is the rainy season in São Carlos, where vegetation and pasture become denser due to a higher metabolic activity, as emphasized by the NDVI values. As a result, the urban and rural contrast is sharper, and the SUHI intensities tend to be higher. In late autumn and winter, there was little difference in the SUHIs intensity from urban to rural areas. In autumn, the amplitude of SUHI intensity was lower than in other seasons, ranging from −2.6 to 4.4 • C because, during this period, the Cerrado and pastures dry or lose their leaves to reduce the adverse effects of the lack of rainfalls. Besides that, green areas are known to decrease surface temperature due to shading and evaporative cooling [83,84]. Urban Sci. 2020, 4, x FOR PEER REVIEW 11 of 19 Figure 6. Urban surface heat islands, NDVI and RGB (Red, Green, Blue) images obtained in the summer, autumn, winter and spring. The sites A, B, C, D and E, represented by dashed ellipses, refer to areas that show different patterns in SUHI and NDVI intensity.
The sites A, B, C, D and E represent areas with different surface types. A, B and D are high urbanization sites, and the C represents a green area with riparian forest, while the E eucalyptus plantation represents an area outside the urban limits. Regarding the sites A, B and D, all seasons presented SUHIs with a higher intensity in the summer. The site D had the greatest SUHI of 4.8 °C (Figure 7). The winter was the season with the lowest SUHI intensities. The lowest NDVI values were observed at these sites. The site C showed only SUFIs with a higher intensity in the summer (−2.4 °C) The sites A, B, C, D and E represent areas with different surface types. A, B and D are high urbanization sites, and the C represents a green area with riparian forest, while the E eucalyptus plantation represents an area outside the urban limits. Regarding the sites A, B and D, all seasons presented SUHIs with a higher intensity in the summer. The site D had the greatest SUHI of 4.8 • C (Figure 7). The winter was the season with the lowest SUHI intensities. The lowest NDVI values were observed at these sites. The site C showed only SUFIs with a higher intensity in the summer (−2.4 • C) because of green area with riparian forests and a high NDVI, located in a valley bottom, which favors the absorption of solar radiation throughout the day, especially in the morning (Landsat 8 satellite passing time). Site E had the highest SUFIs, reaching −3.4 • C in the winter and −3 • C in the summer. Unlike other areas, site E showed a change in its thermal pattern from SUFI to SUHI in spring. It nevertheless obtained the highest intensity among the sites (5.5 • C). This site was covered by eucalyptus. However, trees were cut in the spring, leaving the soil exposed and substantially changing the NDVI (from 0.73 to 0.3), confirming the effects of green areas on surface temperature and consequently on SUHI/SUFI intensity.
Urban Sci. 2020, 4, x FOR PEER REVIEW 12 of 19 because of green area with riparian forests and a high NDVI, located in a valley bottom, which favors the absorption of solar radiation throughout the day, especially in the morning (Landsat 8 satellite passing time). Site E had the highest SUFIs, reaching −3.4 °C in the winter and −3 °C in the summer. Unlike other areas, site E showed a change in its thermal pattern from SUFI to SUHI in spring. It nevertheless obtained the highest intensity among the sites (5.5 °C). This site was covered by eucalyptus. However, trees were cut in the spring, leaving the soil exposed and substantially changing the NDVI (from 0.73 to 0.3), confirming the effects of green areas on surface temperature and consequently on SUHI/SUFI intensity.

Relationship Between SUHI, UI, NDVI, NDWI and TPI
In order to determine the correlation between SUHIs and natural and anthropogenic variables, a Pearson's correlation coefficient was performed. One thousand samples were randomly selected to investigate the relationship between SUHI and altitude, UI, NDVI, NDWI and TPI during the seasons. The results of the correlation analyses are summarized in Table 4 and Figure 8. The correlation between the SUHI and the variables has a seasonal factor. The greatest variations occurred among the variables SUHI, UI and NDWI, i.e., variations of 0.55 and 0.23 in the coefficients, respectively. Seasonal variation may cause an inverse relationship: one season may show a positive correlation and the other a negative one. Haashemi et al. [80] observed that the relationship between NST (Normalized Surface Temperature) and the variables SVF (Sky Vision Factor), albedo, WS (Waterproof Surface) and altitude showed positive and negative correlations depending on the season and time. This change in the direction of significant correlations (p-value < 0.05) was not observed in this study because the variables used in one research study do not always have the same pattern as the SUHIs in other studies. The local settings of the urban site and the rural area create a complex system, varying their patterns in space and time.

Relationship between SUHI, UI, NDVI, NDWI and TPI
In order to determine the correlation between SUHIs and natural and anthropogenic variables, a Pearson's correlation coefficient was performed. One thousand samples were randomly selected to investigate the relationship between SUHI and altitude, UI, NDVI, NDWI and TPI during the seasons. The results of the correlation analyses are summarized in Table 4 and Figure 8. The correlation between the SUHI and the variables has a seasonal factor. The greatest variations occurred among the variables SUHI, UI and NDWI, i.e., variations of 0.55 and 0.23 in the coefficients, respectively. Seasonal variation may cause an inverse relationship: one season may show a positive correlation and the other a negative one. Haashemi et al. [80] observed that the relationship between NST (Normalized Surface Temperature) and the variables SVF (Sky Vision Factor), albedo, WS (Waterproof Surface) and altitude showed positive and negative correlations depending on the season and time. This change in the direction of significant correlations (p-value < 0.05) was not observed in this study because the variables used in one research study do not always have the same pattern as the SUHIs in other studies. The local settings of the urban site and the rural area create a complex system, varying their patterns in space and time. The correlation coefficient of SUHI with altitude was 0.08 in the fall, denoting a low correlation. However, the p-value was below the significance level of 0.05. Great caution is required when interpreting p-values, as this measurement is greatly influenced by sample size. Large samples tend to result in low p-values, although the observed effect is not of great practical importance [85]. Altitude is an important geographical factor for air temperature variation. However, the relationship with SUHI variation was not equally observed (Table 4), since it is a surface phenomenon. Li et al. [62] showed that this relationship could occur due to the variation in atmospheric emissivity that occurs with the variation of altitude, which directly influences the radiation that reaches the surface. However, the differences in altitude in the urban area were not expressive at an atmospheric scale sufficient to establish a significant correlation. The correlation coefficient of SUHI with altitude was 0.08 in the fall, denoting a low correlation. However, the p-value was below the significance level of 0.05. Great caution is required when interpreting p-values, as this measurement is greatly influenced by sample size. Large samples tend to result in low p-values, although the observed effect is not of great practical importance [85]. Altitude is an important geographical factor for air temperature variation. However, the relationship with SUHI variation was not equally observed (Table 4), since it is a surface phenomenon. Li et al. [62] showed that this relationship could occur due to the variation in atmospheric emissivity that occurs with the variation of altitude, which directly influences the radiation that reaches the surface. However, the differences in altitude in the urban area were not expressive at an atmospheric scale sufficient to establish a significant correlation.
A study conducted in Ljutomer (Slovenia) by Ivajnšič et al. [65] found a slightly positive relationship between mean air temperature and TPI, which means that convex areas are slightly warmer than concave areas of that city. In São Carlos, TPI showed no significant correlation in any season, as well as altitude. Therefore, although air temperature, surface temperature, UHI and SUHI are directly related, and the correlations of altitude and TPI did not extend to other phenomena.
The NDWI showed a positive and significant correlation with SUHI in all seasons. The highest correlation was observed in the summer and autumn (0.82 and 0.83) ( Table 4 and Figure 8), which disagrees with the results of the Deilami et al. [86], Lee et al. [87] and Tyubee et al. [88], who found a negative association of SUHI and surface temperature with NDWI.
Green areas play an important role in radiation balance. NDVI is the most used index to diagnose the condition of the vegetation. Its relationship with NST and SUHI has been observed in several studies. The negative correlation is the most observed result [79,89,90], corroborating the results found in this research. There were negative correlations in all seasons, which were more significant in the summer ( Figure 8). However, it is possible to notice different results, such as in Lee et al. [87], who observed a negative relationship between NST and NDVI when NDVI was greater than zero and a positive relation when NDVI was lower than zero. Alves [23] found a similar result, noting that an NDVI equal to 0.2 was a divisor of the positive and negative relationship with NST; up to this value, the relation was positive, i.e., the higher the value of the NDVI, the higher the value of NST, while for a NDVI greater than 0.2, the relation was negative. Rasul et al. [91] observed negative relationships in the autumn and spring. However, the relationship was positive in the winter and non-significant in the summer. In a study by Sun and Kafatos [92], the correlations between NDVI and NST were generally positive in the winter and early spring. However, the NDVI-NST relationship was negative in the summer and autumn. Bala et al. [17] did not report a significant correlation.
All stations had significant correlations with UI, NDVI and NDWI ( Table 4). The summer had the highest correlations (Table 4 and Figure 8), in which the effects of the artificial variable (UI) and natural variables (NDVI and NDWI) were more pronounced because of major contrasts between urban and rural areas.

Conclusions
The urbanization process has generated cities with diverse land use and covers according to the spatial, political, economic and social configuration of the region in which they are inserted. However, cities have developed a common feature: SUHI, a phenomenon that alters the energy balance and directly influences the local climate.
In the city of São Carlos, there was a predominance of SUHIs over SUFIs (Surface Urban Freshness Island). SUHIs, on average, occupied 86% of the urban area, while SUFIs occupied only 14%. This finding confirms need for urban planning that takes into account the microclimate as an inseparable factor for sustainable development. Most SUFIs are located near watercourses, parks, slopes and valley bottoms, revealing the effects of green areas and relief on the SUFIs intensity.
The associations between SUHIs and all five independent variables (Altitude, UI, NDVI, NDWI and TPI) show significant seasonal variations. Linear relationships were observed for the variables NDVI, NDWI and UI.
Several limitations deserve to be discussed. First, this study did not include variables related to urban functions due to a lack of data in the models. The urban climate is an integrated product of a complex interactivity of urban form and urban function with the overlying atmosphere. Therefore, it is also important to address urban function, such as industrial and anthropogenic heat emission in modeling.
Despite the constant frequency of satellite sensor data sampling, in the case of Landsat 8, periods of heavy cloudiness prevent further spatiotemporal analysis of the urban thermal environment. In addition, the Landsat 8 satellite passes over study area, usually at 10 a.m., making intra-daily analysis impossible. Therefore, the results obtained must be interpreted carefully.
However, the relationship between SUHIs and urban landscape has a great seasonal variability. This seasonal effect is generally created through the seasonal variation of solar radiation and the landscape's biogeophysical properties. For example, evergreen trees are able to reduce exposure to cold air and wind. Deciduous trees are unable to do so. Conversely, evergreen trees intercept solar radiation and prevent it from reaching the surface, reducing the surface temperature.