Seasonal and Spatial Variation of Surface Urban Heat Island Intensity in a Small Urban Agglomerate in Brazil

In recent years, SUHIs (surface urban heat islands) have been greatly emphasized in urban climate studies, since it is one of the climate phenomena most influenced by human action. In this study, temporal and spatial variations of SUHIs in the cities of Ceres and Rialma (Brazil) were investigated; satellite Landsat 8 TIRS/OLI images from 2013 to 2016 were used for this purpose. The results showed that in all seasons, two relationships were observed, one positive and one negative. An NDVI (Normalized Difference Vegetation Index) of 0.2 is the divider of this relationship: up to this value, the relationship is positive, that is, the higher the NDVI value, the higher the surface temperature, while the relationship is negative at an NDVI greater than 0.2. There was high seasonal variation in the SUHIs, with the highest intensities recorded in the spring and summer (±12 ◦C), and the lowest in the winter. These temporal variations were attributed to the annual cycle of precipitation, which directly involves the robustness of the Cerrado vegetation. SUHIs occupied, on average, an area three times larger than the area of SUCIs (surface urban cool islands). The highest values of SUCIs were observed in water bodies and in valley bottoms. Overall, SUHIs showed high intensities; however, a more intense core area, such as in large cities, was not observed.


Introduction
Rapid and unplanned urban growth results in a continuous increase in land use, defined as the conversion of open spaces into artificial surfaces as well as green urban areas and sporting and leisure facilities [1].
The land surface temperature (LST) is an important indicator for quantifying urban heat islands (UHIs) and surface urban heat islands (SUHIs).UHIs have been regarded as the most well-documented example of anthropogenic climate modification within the field of urban climate [2,3].
The UHI phenomenon describes the excess warmth of the urban atmosphere and surfaces compared to non-urbanized rural surroundings.In general, three types of heat islands are recognized: (A) the canopy layer heat island; (B) the boundary layer heat island, and (C) the surface urban heat island.The first two types are atmospheric heat islands produced by urbanization.Both refer to a warming of the urban atmosphere [4].A surface urban heat island (SUHI) refers to the relative warmth of the urban surfaces compared to their non-urbanized surroundings [4][5][6].In contrast, the term surface urban cool island (SUCI) is defined as an urban area where lower surface temperatures prevail compared to non-urbanized dry surroundings.
SUHIs can be controlled by numerous factors, the most important of which are the modification in radiation balance, the emission of anthropogenic heat, the reduction in evapotranspiration from vegetation and soils, thermal accumulation in buildings and pavements and the reduction of average wind speed [2,7].The adverse effects of UHIs include: (A) increasing thermal discomfort; a UHI will increase the duration and the degree of thermal discomfort [8]; (B) air pollution; a UHI increases the production of ozone near the ground [9] and, in the form of mesoscale wind, dispersed air pollution [10]; (C) increase in energy consumption [11]; (D) reduction in water quality; rapid temperature changes in water ecosystems caused by surface runoff from hot pavements and roofs covered with asphalt, mostly stressful and fatal for aquatic life [12]; and (E) increase in per capita water consumption in summer time [13].
Thus, greater knowledge of the thermal consequences of sealed soils could be very useful for urban planners and land-use decision makers for promoting an efficient soil sealing management approach in urban environments.Therefore, the objective of this work was to analyze the seasonal and spatial variation of the surface urban heat island intensity in a small urban agglomerate in Brazil.

Study Area
The population in the urban core of the cities of Ceres and Rialma was 31,245 inhabitants in 2014 [14].The municipality of Ceres has twice the population of Rialma, but Rialma has a municipal area greater than Ceres.The relief of the region is uneven with steep slopes, surrounded by valleys and hills.As seen in Figure 1, the two cities are separated by the Rio das Almas, which features two well-defined regimes (a full regime and an ebb regime).Ceres streets are wide, while in Rialma, they are narrower.
Ceres is a city known regionally for its services in health and education.Ceres and Rialma have an interdependent relationship of working-housing.Therefore, the landscape of Ceres/Rialma is a dimension of the social relationships that are established in the city and region, i.e., it is the result of the social relationships of production and domination that are established there.vegetation and soils, thermal accumulation in buildings and pavements and the reduction of average wind speed [2,7].
The adverse effects of UHIs include: (A) increasing thermal discomfort; a UHI will increase the duration and the degree of thermal discomfort [8]; (B) air pollution; a UHI increases the production of ozone near the ground [9] and, in the form of mesoscale wind, dispersed air pollution [10]; (C) increase in energy consumption [11]; (D) reduction in water quality; rapid temperature changes in water ecosystems caused by surface runoff from hot pavements and roofs covered with asphalt, mostly stressful and fatal for aquatic life [12]; and (E) increase in per capita water consumption in summer time [13].
Thus, greater knowledge of the thermal consequences of sealed soils could be very useful for urban planners and land-use decision makers for promoting an efficient soil sealing management approach in urban environments.Therefore, the objective of this work was to analyze the seasonal and spatial variation of the surface urban heat island intensity in a small urban agglomerate in Brazil.

Study Area
The population in the urban core of the cities of Ceres and Rialma was 31,245 inhabitants in 2014 [14].The municipality of Ceres has twice the population of Rialma, but Rialma has a municipal area greater than Ceres.The relief of the region is uneven with steep slopes, surrounded by valleys and hills.As seen in Figure 1, the two cities are separated by the Rio das Almas, which features two welldefined regimes (a full regime and an ebb regime).Ceres streets are wide, while in Rialma, they are narrower.
Ceres is a city known regionally for its services in health and education.Ceres and Rialma have an interdependent relationship of working-housing.Therefore, the landscape of Ceres/Rialma is a dimension of the social relationships that are established in the city and region, i.e., it is the result of the social relationships of production and domination that are established there.

Data
The Landsat sensors have been monitoring the Earth for more than four decades, providing a continuity of data thoughout their lifetime [15].The first series of these satellites was launched in 1972 and it was named the Earth Resources Technology Satellite; it was later renamed to Landsat 1.Since then, there has been a total of 8 Landsat satellites.Among the eight, Landsat 6 failed to attain

Data
The Landsat sensors have been monitoring the Earth for more than four decades, providing a continuity of data thoughout their lifetime [15].The first series of these satellites was launched in 1972 and it was named the Earth Resources Technology Satellite; it was later renamed to Landsat 1.Since then, there has been a total of 8 Landsat satellites.Among the eight, Landsat 6 failed to attain orbit and fell to Earth in 1993.The remaining satellites have proven to be successful and have provided researchers with massive volumes of data which have been used in many studies [15,16].Landsat data are freely available to download from the United States Geological Survey (USGS) website [17].Table 1 shows the Landsat scenes used in the study.The methodology in this study is divided into sub-algorithms, which are sequentially presented in the next subsections.

Estimation of SUHIs
The calculation of the surface urban heat islands was carried out in six steps: The first was the conversion of digital numbers into radiance; the second was the calculation of the brightness temperature; the third was the calculation of the NDV I; the fourth consisted of estimating the surface soil emissivity from the NDV I values obtained; and the fifth step, after the emissivity was estimated, consisted of determining the corrected surface temperature, which was subsequently converted into the surface urban heat island in the sixth step, as suggested by [15].

Step 1-Conversion of Digital Numbers (DN) to Radiance
The thermal data in satellite imagery of Landsat sensors are stored in digital numbers (DN).Digital numbers are used as a way of representing pixels that have not yet been calibrated.They are a representation of the different levels of radiance in the raster image.After obtaining the satellite images, the first step was the conversion of digital numbers into radiance.Equation (1) shows the equation used to convert DN into spectral radiance [18] of the Landsat 8 satellite TIRS sensor.
where L λ is the spectral radiance in W/(m 2 × sr × µm) at the top of the atmosphere.M L is the band-specific multiplicative rescaling factor from the metadata (radiance_mult_band_10). Q cal is the quantized and calibrated standard product pixel values (DN).A L is the band-specific additive rescaling factor from the metadata (radiance_add_band_10).O i are the corrections published by USGS for calibration of the TIRS bands.

Step 2-Computation of Brightness Temperature (BT)
The brightness temperature was calculated by Equation (2).
Climate 2016, 4, 61 where BT is the brightness temperature in (K), L λ is the spectral radiance at the top of the atmosphere, K 1 and K 2 is a specific constant for conversion into the thermal band.

Step 3-Calculation of NDV I
The normalized difference vegetation index (NDV I) is a simple graphical indicator that can be used to analyse remote sensing measurements, typically from satellite image data, and assess whether the target being observed contains live green vegetation.Consequently, we also explored the transformation of NDV I into values associated with cover fraction using empirical relationships with vegetation indices, as a possible basis function.The Normalized Difference Vegetation Index (NDV I) is given by: where NDV I is the Normalized Difference Vegetation Index.N IR is the near infrared band and R is the red band.

Step 4-Determination of Land Surface Emissivity (LSE)
The Land Surface Emissivity was estimated from the NDV I values.According to [19], when the NDV I of an area is known, the LSE can be estimated.The LSE of a pixel is estimated by classifying the pixels according to the class that they fall into.When a pixel has an NDV I value that is below −0.185, then the LSE value of 0.995 is assigned to the pixel; when the NDV I value is greater than or equal to −0.185 and is less than 0.157, an LSE value of 0.985 is assigned to the pixel, when the NDV I value is greater than or equal to 0.157 and is less than or equal to 0.727, a logarithmic relationship between NDV I and LSE is used [15], and finally, when the NDV I value is greater than 0.727, the pixel is assigned a value of 0.990, as shown in Table 2.After estimating the soil surface emissivity, brightness temperature correction was performed.For this, Planck's function was used.Equation (4) shows Planck's function [20,21].This function corrects the emissivity of a surface, in comparison to a black body.
where LST is the surface temperature (K), BT is the brightness temperature (K), λ is the wavelength of the emitted radiation, ρ is (h To calculate the SUHIs, two buffers were produced, one rural, one urban.The rural buffer is located 500 m from the urban buffer (Figure 2).Usually, the urban heat island (UHI) intensity is measured from observations of the air temperature along transects or in fixed stations in the urban and surrounding rural areas [22,23].Nonetheless, the SUHI has another meaning; therefore, SUHI was calculated as the difference between the surface temperature of each pixel of the buffer in the urban area and the average surface temperature of the buffer in the rural area [5,24], according to Equation ( 5), shown in Figure 2.
where LST U is the surface temperature of the buffer urban and LST is the surface temperature of the buffer rural area.A positive value represents an SUHI situation, while a negative result represents a SUCI.
Climate 2016, 4, 61 5 of 11 Usually, the urban heat island (UHI) intensity is measured from observations of the air temperature along transects or in fixed stations in the urban and surrounding rural areas [22,23].Nonetheless, the SUHI has another meaning; therefore, SUHI was calculated as the difference between the surface temperature of each pixel of the buffer in the urban area and the average surface temperature of the buffer in the rural area [5,24], according to Equation ( 5), shown in Figure 2.
where is the surface temperature of the buffer urban and LST is the surface temperature of the buffer rural area.A positive value represents an SUHI situation, while a negative result represents a SUCI.
In Ceres-Rialma, the seasonality was evident, as shown in Figures 3 and 4. In the winter, the intensity of SUHI was lower than in other seasons.The highest intensities of the SUHIs occurred in the summer, i.e., 12 °C.In autumn and spring, there was a significant time difference between the SUHI values.
The highest amount and intensity of the negative SHUIs, also called the surface urban cool islands (SUCIs), was observed in winter.SUHI bloxplots (Figure 3) contained all data presented in Figure 4, and therefore not only represent the urbanized area but also the rural area.
In Ceres-Rialma, the seasonality was evident, as shown in Figures 3 and 4. In the winter, the intensity of SUHI was lower than in other seasons.The highest intensities of the SUHIs occurred in the summer, i.e., 12 • C. In autumn and spring, there was a significant time difference between the SUHI values.The seasonal variation in the SUHIs (Figure 4) is likely to result from the change in the intensity of solar radiation due to the movement of the Earth and the variation in the Cerrado vegetation and pastures, according to the rainfall [28].
The highest intensities in the SUHIs were observed in the summer (31 January 2015), a time that corresponds to the rainy season in the region, where the vegetation and pasture become more dense with increased metabolic activity.At this time, the contrast between urban and rural areas is higher, and therefore, the intensities of the SUHIs tend to be higher.By contrast, in the winter (11 August  corresponds to the rainy season in the region, where the vegetation and pasture become more dense with increased metabolic activity.At this time, the contrast between urban and rural areas is higher, and therefore, the intensities of the SUHIs tend to be higher.By contrast, in the winter (11 August 2015), there was little difference in the intensity of the SUHIs of the urban area compared to the SUHIs of the rural area, as in this period, the Cerrado vegetation and pastures dry out or lose their leaves to minimize the adverse effects of the lack of precipitation.The highest amount and intensity of the negative SHUIs, also called the surface urban cool islands (SUCIs), was observed in winter.SUHI bloxplots (Figure 3) contained all data presented in Figure 4, and therefore not only represent the urbanized area but also the rural area.
The seasonal variation in the SUHIs (Figure 4) is likely to result from the change in the intensity of solar radiation due to the movement of the Earth and the variation in the Cerrado vegetation and pastures, according to the rainfall [28].
The highest intensities in the SUHIs were observed in the summer (31 January 2015), a time that corresponds to the rainy season in the region, where the vegetation and pasture become more dense with increased metabolic activity.At this time, the contrast between urban and rural areas is higher, and therefore, the intensities of the SUHIs tend to be higher.By contrast, in the winter (11 August 2015), there was little difference in the intensity of the SUHIs of the urban area compared to the SUHIs of the rural area, as in this period, the Cerrado vegetation and pastures dry out or lose their leaves to minimize the adverse effects of the lack of precipitation.

Spatial Characteristics of the SUHIs
In order to analyze the temporal and spatial variation of the average intensity of the SUHIs, a transect from point A to point B (Figure 5) was established, with a regular distance of 30 m, totaling 225 data points.The profile of Figure 5 shows the influence of different types of land use and occupation on the intensity of the SUHIs; initially, in the first few meters of the transect, SUHIs close to 8 • C were observed in the summer.Between 1200 and 1900 m, lower intensities were measured.From the point of 3600 to 4600 m, lower intensities of the SUHIs were observed in all seasons.This occurred because this space is outside the urban boundary, encompassing vegetation sites and the Rio das Almas.

Rio das Almas.
It can be noted in Figure 5 that in almost all transects in winter, the intensities of the SUHIs were lower than in other seasons, while in the summer, most SHUIs had the highest intensities.The spatial variation of the SUHIs in the autumn and spring had similar patterns.Therefore, there is a sharp contrast, especially in the summer and winter seasons, in the intensity of SUHIs in the urban agglomerate of Ceres-Rialma.The average SUHIs superimposed on the relief (Figure 6) enable one to check that valley bottom areas or steep slope areas showed the highest intensities of the SUCIs.Rio das Almas, which divides the cities of Ceres and Rialma, has lower intensities due to the heat capacity of water and the lower intensity of solar radiation that reaches the inclined surfaces.It can be verified that the urban area of Ceres and Rialma had, on average, the highest intensities of SUHIs (Figure 6), which was verified in numerous studies [1,5,26,27,[29][30][31][32].
As shown in Figure 7 and Table 3, in the study area, there was predominance of SUHIs at the expense of SUCIs.SUHIs occupied an area three times larger than SUCIs; this finding reflects a typically urban problem, the lack of green areas and trees on sidewalks, as it is noted that most of the existing SUCIs are located next to watercourses and steep slopes.SUHIs of 1 °C to 2 °C occurred in most of the study area (21.75% of the total area).The intensity of 5 °C to 6 °C covered a smaller area (0.9 km 2 ), and in the SUCIs with intensity of −6 °C to −5 °C, there was an even smaller area (0.03 km 2 ).The largest area of SUCIs was found from −1 °C to 0 °C (3.93 km 2 ).It can be noted in Figure 5 that in almost all transects in winter, the intensities of the SUHIs were lower than in other seasons, while in the summer, most SHUIs had the highest intensities.The spatial variation of the SUHIs in the autumn and spring had similar patterns.Therefore, there is a sharp contrast, especially in the summer and winter seasons, in the intensity of SUHIs in the urban agglomerate of Ceres-Rialma.
The average SUHIs superimposed on the relief (Figure 6) enable one to check that valley bottom areas or steep slope areas showed the highest intensities of the SUCIs.Rio das Almas, which divides the cities of Ceres and Rialma, has lower intensities due to the heat capacity of water and the lower intensity of solar radiation that reaches the inclined surfaces.It can be verified that the urban area of Ceres and Rialma had, on average, the highest intensities of SUHIs (Figure 6), which was verified in numerous studies [1,5,26,27,[29][30][31][32].As shown in Figure 7 and Table 3, in the study area, there was predominance of SUHIs at the expense of SUCIs.SUHIs occupied an area three times larger than SUCIs; this finding reflects a typically urban problem, the lack of green areas and trees on sidewalks, as it is noted that most of the existing SUCIs are located next to watercourses and steep slopes.SUHIs of 1 • C to 2 • C occurred in most of the study area (21.75% of the total area).The intensity of 5 • C to 6 • C covered a smaller area (0.9 km 2 ), and in the SUCIs with intensity of −6 • C to −5 • C, there was an even smaller area (0.03 km 2 ).The largest area of SUCIs was found from −1 • C to 0 • C (3.93 km 2 ).

Relationships between LST and
In the summer and spring, periods of high solar radiation and insolation, the highest LST (±40 °C) and (0.57) values were observed.In Erbil, Iraq, a maximum LST of 53 was observed [27].The lowest LST values were observed in the autumn (minimum of 18.87 °C), followed by winter.The standard deviations of LST were similar in all seasons; in relation to the , the smallest deviations deviations occurred during autumn and winter, as shown in Table 4.
The relationship between LST and is not always negative [25,27,33] and may vary seasonally, depending on the type of vegetation.In Ceres and Rialma, there is an interesting observation (Figure 8): in all seasons, two relationships were observed, one positive and one negative.An of 0.2 is the divider of this relationship; up to this point, the relationship is positive, that

Relationships between LST and NDV I
In the summer and spring, periods of high solar radiation and insolation, the highest LST (±40 • C) and NDV I (0.57) values were observed.In Erbil, Iraq, a maximum LST of 53 was observed [27].The lowest LST values were observed in the autumn (minimum of 18.87 • C), followed by winter.The standard deviations of LST were similar in all seasons; in relation to the NDV I, the smallest deviations occurred during autumn and winter, as shown in Table 4.
The relationship between LST and NDV I is not always negative [25,27,33] and may vary seasonally, depending on the type of vegetation.In Ceres and Rialma, there is an interesting observation (Figure 8): in all seasons, two relationships were observed, one positive and one negative.An NDV I of 0.2 is the divider of this relationship; up to this point, the relationship is positive, that is, the higher the NDV I value, the higher the surface temperature, while above an NDV I of 0.2, the relationship is negative, i.e., the higher the NDV I, the lower the surface temperature.This relationship is interesting because it shows that vegetation with poor physiological conditions cannot minimize the surface temperature in the study area.
is not observed: the center consists of commercial buildings with a maximum of three floors and households; the area surrounding the center has mainly, households; the peripheral region has spaced homes.This configuration does not allow such a difference as in large cities but creates a relatively homogeneous area of SUHIs.

Figure 1 .
Figure 1.Location of the cities of Ceres and Rialma.

Figure 1 .
Figure 1.Location of the cities of Ceres and Rialma.

Figure 2 .
Figure 2. Delimitation of the urban and rural buffer.

igure 2 .
Delimitation of the urban and rural buffer.

Figure 4 .
Figure 4. Seasonal variation in the intensity of SUHIs.Figure 4. Seasonal variation in the intensity of SUHIs.

Figure 4 .
Figure 4. Seasonal variation in the intensity of SUHIs.Figure 4. Seasonal variation in the intensity of SUHIs.

Figure 5 .
Figure 5. Spatial and seasonal variation of the average SUHIs along line A-B.

Table 3 .
Area occupied by the classes of average SUHIs and SUCIs in the Ceres and Rialma cities and surroundings.

Figure 5 .
Figure 5. Spatial and seasonal variation of the average SUHIs along line A-B.

Figure 6 .
Figure 6.Average SUHIs of the cities of Ceres and Rialma, superimposed on the relief.

Figure 6 .
Figure 6.Average SUHIs of the cities of Ceres and Rialma, superimposed on the relief.

Figure 6 .
Figure 6.Average SUHIs of the cities of Ceres and Rialma, superimposed on the relief.

Figure 7 .
Figure 7. Area occupied by the average SUHIs and SUCIs.

Table 2 .
Algorithm based on the NDV I image.