Economic Driven Patterns of the Urban Environment

This study investigates the land surface temperature (LST) distribution from thermal infrared data for analyzing the characteristics of surface coverage using the Vegetation-ImperviousSoil (VIS) approach. A set of ten images, obtained from Landsat-5 Thematic Mapper, between 2001 and 2010, were used to study the urban environmental conditions of 47 neighborhoods of Porto Alegre city, Brazil. Porto Alegre has had the smallest population growth rate of all 27 state capitals in the last two decades in Brazil, with an increase of 11.55% in inhabitants from 1,263 million in 1991 to 1,409 million in 2010. We applied the environmental Kuznets curve (EKC) theory in order to test the influence of the economically-related scenario on the spatial nature of social-environmental arrangement of the city at neighborhood scale. Our results suggest that the economically-related scenario exerts a non-negligible influence on the physically driven characteristics of the urban environmental conditions as predicted by EKC theory. The linear inverse correlation R2 between household income (HI) and LST is 0.36 and has shown to be comparable to all other studied variables. Future research may investigate the relation between other economically-related indicators to specific land surface characteristics.


Introduction
The powerful transformation initiated by the economic development in the 19th century is leading civilization to a continuous and rising population migration from rural to urban areas.At the same time, excessive local demands on environmental systems, due to urban growth demand, have become global in scope [1].Urban areas exist because living in concentrated spaces brings about economic and political advantages [2] in view of the fact that this creates interdependency among people due to the supply and demand for products of physical and spatial perspective [3].The agglomeration in space, through social and functional organizations, generates constructive densification, which, in turn, makes it possible for people to access more rapidly the activities in a shorter time and distance.Nevertheless, the high constructive and population density, allied with the diversity of activities, can cause the urban areas to enter a degenerative trend, because while they generate an agglomeration economy, maximizing urban equipment and infrastructure, they equally generate anti-economies [4], like holdups, soaring localization costs and alterations to the urban microclimate in terms of higher temperature.Since then, competition in soil occupation processes mainly due to urbanization and industrialization, has been demonstrated to be in disagreement with sustainable development concepts.Nowadays, although there are a number of diverse definitions of urban sustainability, the core value of urban sustainability always lies in the balance of environmental, economic and social development [5,6,7].
In highly urbanized areas, vegetation and impervious surface are two key urban components [8] and combined they are important indicators of environmental quality [9].Usually, Impervious Surface Areas (ISA) are considered as the complement of vegetation coverage fraction [10].
Cities use construction materials such as concrete and asphalt, which do not allow water to penetrate into the soil, leading to most of the incident shortwave radiation being absorbed and transformed into sensible heat [11].So, under certain imbalanced conditions between percent of vegetation coverage, impervious surface types and soil in the urban environment, with a large expanse of non-evaporating impervious materials, there is a consequent increase in sensible heat flux at the expense of latent heat flux [12].As a result, surface and atmospheric modifications caused by urbanization generally lead to a modified thermal climate, which is warmer than the surrounding non-urbanized areas [13].This phenomenon, which modulates the air temperature of the lowest layers in the urban atmosphere is called an Urban Heat Island (UHI).The knowledge about UHI is very important in earth sciences for urban climatology, global environmental change and management practices.UHIs are mainly due to a high percentage of non-reflective, water-resistant surfaces coverage and a low percentage of vegetated and moisture trapping surfaces [14] and are a particular public policy concern during the summer time due to their association with increases in energy, air pollution, and heat stress-related human health.Due to this phenomenon, the distribution of LST has a close relation to the spatial configuration of land cover characteristics in the urban environment [15].Under these conditions, it is quite conceivable that presence of urban green spaces, moisture trapping surfaces within other physical characteristics of urban design are always beneficial to mitigating UHI.
In surface areas characterized by fully or fractionally vegetated cover, thermal processes related to direct sunlight interception can significantly influence the surrounding measurements of LST.Furthermore, urban green spaces and trees have been recognized as an important component of urban landscapes [6,7].Several benefits can be associated with green spaces in the urban environment caused by their capacity for cooling the air [16] and reducing the available surfaces for successively radiating energy reflectance [17].The evapotranspiration process generated by trees and vegetation in cities can cool the air by using the surrounding urban environment heat from the air to evaporate water [18].A study from [19] concluded that a quality forest or green space has a positive economic ripple effect on nearby properties, because urban reforestation in quantity, size and appropriate form makes the cities comfortable and more livable.
In some countries, recent planning and research on the thermal characteristics related to its economic growth have been gaining interest.Studies have shown that urban trees increase the value of residential properties and also make neighborhoods aesthetically more appealing and add to the value of property, regardless of the tree species [20].The main assumption is that planting trees in cities has high costs and demands planning from policymakers and population in a way that lot size and tree presence reflect, to some extent, the market forces determined by the welfare of the citizens and their preferences [6].
In the 1980s, initial studies relating environmental quality indicators to income were based on the concept of the environmental Kuznets curve (EKC).The EKC effect is an empirical phenomenon and the main assumption of this concept is that economic growth through time is necessary in order for environmental quality to be maintained.This is an essential part of the sustainable development argument [21].
Most of the EKC studies limited the theory to principally problems related to air pollutant concentrations [21].Recent evidence on the effective mitigation of CO2 emissions per capita, and cheaper energy use intensity after 2005 [22], has reinforced it in highly developed economies in the world, such as Australia, France, Germany, United Kingdom and United States.Nevertheless, a more recent study from [6] has observed that urban forests were highly related to income for cities with a population over 100,000 in the United States.A detailed study of [23] in the city of Santiago, in Chile, also found that socioeconomic strata were related to different tree cover densities.In Brazil, preliminary analysis from [24] has shown that high LST averages at neighborhood scale were better spatially linked to those ones with lower sustainability conditions than vegetation index alone.In Australia, the Urban Climates Research (UCR) of the Cooperative Research Centre for Low Carbon Living (CRC-LCL), has been conducting several scale analyses of the urban environment conditions to UHI understanding and mitigation.
Considering this research issue, a literature review has indicated that thermal remote sensing data obtained from satellites, with adequate temporal and spatial resolution data distribution can improve the understanding of the main surface physical properties in the urban ecosystem conditions [25][26][27][28][29][30][31].
Currently, the NASA Landsat Data Collection (NLDC) provides historical Landsat data products with medium spatial resolution of 30 m, with a geometric quality of the images that allows the imagery composition for time series analysis and guarantee improved pixel geolocation.Historically, Landsat data has been widely used for urban development and UHI studies, mainly due to its long-term record provided with almost 40 years of data archive [31,32,33] and freely available for research.Landsat data series has an adequate spatial scale for detecting changes over many urban environmental conditions that allows long-term studies of regional and global land cover change [32].In this context, the launch of the LDCM (Landsat Data Continuity Mission) on February 11, 2013, can also be used for continued urban environment research and management by providing reliable spatial information.
In this study, we demonstrate the influence of socio-economic development on local microclimate conditions in the urban environmental scenario.We go further showing that LST can be used as an integrative indicator of general urban environmental conditions on several levels of social organization in Porto Alegre city in Brazil.Within this approach, we also applied the EKC theory in order to test the influence of the economically-related scenario on the spatial nature of social-environmental arrangement of the city, at neighborhood scale.

Study area
The study area is the city of Porto Alegre, which is the capital of Rio Grande do Sul State (RS) and was founded on March 26, 1772.The city is the capital of the most southern state in Brazil and it is located at latitude of 30°01′59′′ South and longitude of 51º13′48′′ West.Porto Alegre is well suited to this study because it has passed through a period of slow growth in the last two decades, with an increase of 11.55% in inhabitants from 1,263 million in 1991 to 1,409 million in 2010.This represents the smallest population growth rate of all 27 capitals in Brazil and far lower than the overall average growth rate of capitals with 43.75% and also from overall Brazilian average of 29.92% [34].
The region predominantly experiences a subtropical mid-latitude climate (Cfb) with four welldefined seasons [36].The mean annual temperature is 19.5 o C, with monthly averages varying between 12 o C and 26 o C. The monthly mean precipitation is 114 mm and the yearly-accumulated average precipitation is 1373 mm with no dry period.Figure 1 shows the study area.The city is completely covered with a Landsat scene (Path 221 and Row 081).The neighborhood scale elevation average of the study area is 39 m above the sea level with a standard deviation (SD) of 14 m.The maximum average altitude is 103 m and minimum is 7 m.The study area comprises 47 densely populated neighborhoods surrounding the downtown, covering a total area of 9077 ha.

Satellite Imagery and Data Set
Absolute LST measures are commonly used in physiological studies because they is related to vegetation canopy [37].Multi-temporal imagery analysis is essential for UHI studies [38].LST obtained from satellite imagery is, strictly speaking, a measure of the "skin temperature" or surface radiometric energy (kinetic) emitted from the land surface and is related to the thermal infrared (TIR) radiation from surface rather than air temperature [39].
The Landsat data series satellites have collected several years of data in the thermal spectral band and have been used in several studies for UHI [40].However, the challenge in urban environmental LST estimation is the numerous factors that need to be quantified in order to assess accurate LST retrieval from thermal satellite data.Those factors include sensor radiometric calibration [41,42], atmospheric correction [43,44,45], simplified atmosphere approaches [39] but mainly surface emissivity correction [46] and changing physically driven dynamics of land coverage through time [47,31].Due to radiometric and atmospheric conditions only ten images were selected.The generated image product is composed of cloud-free images covering the 47 neighborhoods studied.
The information levels and data sources included in this study are the following: i. Landsat-5 TM imagery from 2000 to 2010, Path and Row 221-081 [48]; ii.Shuttle Radar Topography Mission (SRTM) data [49], version 3, product GL1, was used to analyze the relationship between elevation and distribution of other variables studied on a 30-meter spatial resolution [48]; iii.Monthly rainfall data obtained from the Database for Meteorological Research [50], which covers the period from October to December.This data was used to identify possible drought periods in the warmest season; iv.Soil types map in a 1:5,000,000 scale [51]; v. Census data of Neighborhoods limits and urban sectors [34]; vi.Census data of Population by neighborhoods in 2010 [34]; vii.Climatological norm of accumulated rainfall from October to December [50;52]; viii.Census data of Household Incomes (Brazilian salary units) at neighborhood scale, from 2000 and 2010 [34]; ix.Digital geospatial reference from National Aeronautics and Space Administration-Global Land Survey [48].

Reflectance Data Generation
Timely analysis based on remote sensing data, such as the Landsat series, requires accurate input parameters from meta data, ground based meteorological measurements and look up tables based on surface-related parameters [32,53,54].Full absolute image correction includes atmospheric and sensor-related parameters and then the derivation of physical units, such as reflectance [17].Planetary reflectance (or apparent) is limited by atmospheric characteristics.To compare imaged targets at different times and locations and even to compare targets imaged by different sensors, required transforming the images to reflectance of the surface, since this is a target property, and which can be compared for checking differences and surface modeling parameters.In this way, a complete rescaling calibration and correction of digital numbers (DN) was performed by NASA for the conversion of DN into Top-of-Atmosphere (TOA) reflectance data according to a time-dependent calibration look-up tables (LUT), generated from the lifetime gain model equations for all bands.This methodology substitutes the previous National Landsat Archive Production System (NLAPS) based on internal instrument calibration procedure [54].
In order to overcome future difficulties related to complex vicarious calibration processes and methodologies it is very important to develop an independent technical source for satellite instrument calibration to ensure confidence in the accuracy, interoperability, and quality of satellite observations reducing analysis uncertainties [55].Independent algorithms allow continuous monitoring and provide transparency to the methodology.Transparency is needed not only to permit other investigators to scrutinize the quality of the data, but also so they can pursue alternative research inquiries because parameters, literature and adjacent supporting methodologies are not always available and published [56].Thus, the calibration and validation have been primary to the measurements obtained from satellites and surface parameters determinations through satellite images [53].With this purpose in mind, the radiometric calibration has the purpose of ensuring the update of the quality of data generated by a remotely located sensor and ensuring the possibility of converting the data recorded by the sensors to physical quantities that can be correlated to geophysical parameters [57].
Regarding the calibration of reflectance channels, for the accurate transformation of radiance into TOA and then into surface reflectance data, the images were also atmospherically corrected according to [43], as variations of atmospheric conditions are spatially and temporally significant [54].
First of all, DN were converted into radiance and next to reflectance, according to the calibration parameters of [54] and LST [41].After conversion to at-satellite radiance, each image was converted to at-satellite reflectance (assuming a uniform Lambertian surface under cloudless conditions) according to the general mathematical function for reflectance calculation shown in Equation 1. Radiometric calibration coefficients and parameters were obtained from [54].
The TOA calculation can be described by applying the simple radiance-reflectance conversion for each band of LANDSAT-5 TM according to Equation 1.
where, Lλ is the spectral radiance at the sensor in each band (W m −2 sr −1 μm −1 ); d is the distance between the Earth and the Sun in astronomic units (UA); ESUNλ is the average solar exoatmospheric irradiance (W m −2 μm −1 ); and θz is the solar zenith angle (degrees).
The atmosphere correction process is applied by substituting for the term L obtained during the DN conversion to radiance by the factor (Lsat -L p )/Tz when correcting TOA to canopy-surface reflectance in the Equation 1, according to Equations 2, 3 and 4. The atmosphere correction process is based on the Dark-Object Subtraction (DOS) concept which depends on the wavelength of visible bands and assumes near-zero percent reflectance of specific features [17].The DOS approach is related to a weighted displacement effect to each band on a per-pixel-basis, according to the processes described in [58] and [43].By using this approach, the data necessary to perform atmospheric correction in the visible bands 1-5 from Landsat 5 TM can be obtained from the image itself [46] by restricting the atmospheric influence, as the radiation transfer from the Earth's surface to on-orbit sensors.As an additional control for radiometric quality and a reduction in the amount of atmosphere effects, Landsat-5 TM imagery was submitted to a threshold criterion for the image processing protocol from [47], whereby no image was used when the amount of atmospheric contamination effects exceeded 70 DN in the blue band (452-518μm), of Landsat-5 TM.

Classification of surface cover types
Once the surface reflectance by correcting TOA reflectance was calculated, we applied an empirical reflectance-based method.The VIS approach from [8] conceptualizes an urban environment in terms of three primary physical components: vegetation (V), impervious surfaces (I), and soil (S) coverage [10].The modified VIS model considered the spectral signature of pixel covering urban areas as a linear combination of those three components.
Impervious surface areas (ISA) are surfaces where water cannot infiltrate and are primarily associated with transportation (streets, highways, parking lots and sidewalks) and building rooftops.ISA extractions based on an image classification is a complex task and is often underestimated due to the heterogeneity of urban landscapes due its complexity [59].ISA calculation is based mainly on identifying and categorizing the spectral profile of highly reflective surfaces such as pavements, rooftops and buildings, non-vegetated or even sparsely vegetated areas.At pixel and sub-pixel level impervious surfaces may be mixed with other land cover types, such as trees, grasses, and soils [59].
The non-blue band versions of Enhanced Vegetation Index (EVI-2) can be computed and then can be used on the LST-emissivity dependence calculation.The adequate correction of atmospheric effects is very important because it will allow correct EVI-2 values to be generated, which will also be used for the classification of surface coverage types.EVI-2 is a two-band version of EVI that has been developed for sensors without a blue band [60].Reflectance values from bands 3 and 4 were used to generate the EVI-2 data according to Equation 2 from [60].EVI-2 retains sensitivity and the same as EVI for high leaf area index (LAI) canopies, but does not rely on the usually poor-quality blue band [61].In this way, EVI-2 calculation (Equation 9) can be used as an acceptable substitute of EVI over atmospherically corrected pixels [60].
Soil classification was obtained by considering those pixels that were not classified either urban neither densely vegetated areas when in accordance to Equation 7. The main assumptions of this mathematical approach are summarized in equations 5 to 9, as follows: = 0.19 < < 0.29; 0.12 < < 0.19; 0.07 < < 0.12; 0.02 < < 0.07 = 0.29 < 2 (8) where, %ISAN is the percent of impervious surface area of a given neighborhood; NA is the total area of the neighborhood (ha); ISA is the impervious surface area of a given neighborhood (ha); VegetationA is the vegetated surface area (ha); SoilA is the soil area of a given neighborhood (ha); Soil represents the obtained reflectance of dry bare soil coverage of Landsat-5 TM at a given pixel, empirically obtained for red clay soils and red-yellow clay soils from Porto Alegre region; NIR, Red, Green and Blue represent the obtained reflectance in the near-infrared, red, green and blue bands of Landsat-5 TM at a given pixel, respectively.
Vegetation represents the EVI-2 of green spaces and vegetated areas covering the surface background obtained with Landsat-5 TM image at a given pixel; EVI-2 is the no-blue band Enhanced Vegetation Index (EVI).

Thermal Data Generation
Regarding the thermal analysis from Landsat 5 TM band 6, after converting DN into absolute radiance values, LST is computed from at-satellite brightness temperatures (i.e., blackbody temperature) using pre-launch calibration constants [41,54] and under the assumption of emissivity factor equal to unity [62].Under this assumption (εNB=1), the conversion of the detected thermal radiation into brightness temperature is given in Equation 10, [54], which is a special case of Plank's radiation law.In order for LST to become representative of realistic surface conditions at this point, the LST still needs to be corrected to more realistic values of non-unity surface emissivity [41,63] by means a robust technique of emissivity dependent LST retrieval.Surface emissivity is controlled by factors, such as water content, chemical composition, structure and roughness [64].
This formulation does not result on directly corrections to atmospheric effects (absorption and emissions along the path) because of the difficulty with estimating water vapor content from thermal detection in the mono-window band 6 [62].In the thermal channel, atmospherically corrected surface conditions can be obtained by means an indirect estimation of physical characteristics from atmospherically corrected reflectance of surface in the visible bands, according to the previously mentioned Equation 2. Emissivity calculation for non-black body surface with the narrow band adjustment (NB) for Landsat data series depends on the empirical equation of Soil Adjusted Vegetation Index (SAVI) from [65].Parameters SAVI and LAI are obtained from previously atmospherically corrected band 1 and 2. The value for L in SAVI formulation (Equation 11) can be derived from analysis of multiple images where vegetation does not change, but surface soil moisture does.According to [65] optimal value of L=0.5 (light soil < L=0.5 < dark soil) must be used to account for first-order soil background variations.Afterwards, SAVI (Equation 11) was used for calculation of the LAI which parameters were empirically obtained from Equation 12according to [66].Finally, the adjusted LST (hereafter LST) is obtained by correcting the radiating surface temperature to the surface emissivity (ε), at pixel scale as given by Equation 13, which is the ratio of the thermal energy radiated by the surface to the thermal energy radiated by a blackbody at the same temperature [66,67,68], as follows: where, LST is the emissivity-corrected surface temperature (K); K1 is the calibration constant 1 (607.76W m −2 sr −1 μm −1 ); K2 is the calibration constant 2 (1260.56W m −2 sr −1 μm −1 ); L is the blackbody radiance of the thermal band 6 (W m −2 sr −1 μm −1 ); and εNB is the emissivity factor, which depends on the type of surface coverage conditions when LAI > 3.0, εNB = 0.98 because the increased water content in vegetation actually increases emissivity capacity.where, L is the constant parameter in SAVI with value 0.5; SAVI is the Soil Adjusted Vegetation Index; LAI is the Leaf Area Index empirically adjusted.

LST and Biophysical Descriptors
Although the presence of trees in urban ecosystems is primarily conceptualized in terms of the aesthetic, environmental, and wildlife habitat functions [69], it is well known that shade trees and even small plants for land cover, help cool the urban environment [70].
In the natural environment, dry, bare and low-density soils, for example, have been linked to higher LST as a result of the relatively low thermal inertia [71].These thermal properties vary with the type of soil and moisture content [71].On the other hand, in the urban environment, nonevaporative and impervious surface areas not only play an important role in UHI formation [12] but also in heat generated by anthropogenic sources, such as traffic and industries [30].This is because the emissivity of soils or sparsely vegetated areas is a function of soil moisture conditions and soil density [30,47,72].The physical fundamentals for such a relationship rely on the fact that for any surface material, certain internal properties, such as heat capacity, thermal conductivity and inertia, play important roles in governing the temperature of a body at equilibrium with its surroundings [73].In the summer season, generally 10%-30% of the sun's energy reaches the area below a tree, with most of it being absorbed by leaves and used for photosynthesis and some being reflected back into the atmosphere [70].Even so, thermal responses for vegetation can also vary highly as a function of the biophysical properties of the vegetation itself [15,74].The inverse relationship between thermal characteristics of surface and vegetation indices has been extensively documented in the literature [30,72].This concept is based on the assumption that vegetation coverage mitigates high LST occurrence values [47].Given this and considering that temperature is closely related to physiological activities of vegetation cover, LST can be a useful measure of the physiological activity of the top canopy leaves when leaf cover is sufficiently high such that they are not affected by background temperature from soil [37].However, vegetation index measurements are subject to seasonal variations, which may influence the results on UHI studies, and thus, such vegetation index measurements alone may not serve as a good indicator for urban development [30].
Regarding the seasonal effects and fluctuation of vegetation cover conditions, it is well known that in drought-free years, well-developed vegetation reflects only a small portion of incident solar radiation in the visible band of the spectrum because of chlorophyll absorption properties and other plant pigments that absorb sunlight [31].In the NIR, plants reflect much more because of a scattering effect caused by the internal structure and water content of leaves [75].
Furthermore, the percentage of impervious surfaces coverage under green vegetated areas can also change with very small variations on vegetation index values.Based on studies from [6], which accounted for tree number distributed in the urban environment, we aggregated the public and private trees among the city environment into one single variable, which is representative of the amount of vegetation.

TSDS Approach
The study of LST is closely related to the distribution of LULC characteristics.Models that describe LULC can provide a link between LST and physical spatial distribution in a way that absolute values of LST are not key fundamental variables for the understanding of its distribution is.This is because LST is highly variable through time and does not allow the comparison of absolute values in a simplified conceptual approach.
Detailed studies from [30], conducted in the metropolitan area of Minnesota (US), found that the vegetation index-LST relationship is more suitable for the analysis of UHI during summer and early autumn.In this way, for a better understanding of the physical descriptors of heating effects over the urban environment, as a representative of their development through time, we only compared the summer season imagery from late December to February, in a way to detect the most pronounced LST values associated with surface coverage.In this paper the absolute temperature fluctuations through time have thus not been evaluated.Instead, the focus is on the evaluation of the trend of LST-emissivity dependent in neighborhoods related and its distribution.
A previous analysis of the LST distribution was done by applying the Thermal Spatial Distribution Signature (TSDS) procedure.The TSDS procedure [24,47] using TIR data obtained from Landsat-5 TM was applied as a criterion for multi-temporal imagery combination of data, which is based on the relation between the thermal characteristics of surface and coverage and consists of the standard deviation (SD) analysis of LST to compare images by means similarity characteristics.So, in order to evaluate the most precise combination data for multi-temporal analysis, the key question is what is the time window and their extension for performing an adequate imagery combination in order to perform a reliable LULC diagnosis [47].
According to [47], when studying a histogram of LST distribution through time, it is expected that the highly urbanized or sparsely vegetated areas are unable to cool the surrounding sensible heat by evapotranspiration, as they typically do.Such fractional vegetation and soil cover conditions leads to a spread of the LST range by increasing statistical indicators in terms of the variance of distribution towards higher LST values [47].

Data Retrieval
After choosing a selection of image groups by means of the TSDS approach, we can compare overall data retrieved between neighborhoods for surface coverage types and LST (Table 1 and Table  2).By adequately combining all imagery data, the major patterns of LST distribution could be identified.The overall average of the LST, from all 10 images selected, was 300.36 K (27.21 o C) with a standard deviation of 1.34 K.This low value of standard deviation compared to the mean is also indicative of a robust technique.Note: SD: Standard Deviation.

Spatially distributed patterns of UHI
On the pixel scale, differences of 10.71 K, were found in the LST distribution.Maximum absolute LST with 305.80 K (32.65 o C) was found at São João and minimum of 295.09K (21.94 o C) was found at Jardim Botânico.On the neighborhood scale, differences of 2.83 K were found between neighborhoods averages.Maximum LST average of neighborhood with 301.67 K (28.52 o C) was found at São Geraldo which is characterized by low vegetated areas coverage with high density of impervious coverage.The minimum LST average of 298.65 K (25.50 o C) was found at Farroupilha which is characterized by highly vegetated coverage area and where ISA do not prevail.
Regarding the analysis of LST dependence to EVI-2 we have found that the main pattern is a clear core distribution of lower LST values located at the center of the study area (dark red and orange) which is surrounded by higher LST values (dark blue and purple).The lower LST configuration is associated with high maintenance residential area and the higher LST values are representative of main arterial avenues and streets of the city.Most of the higher LST averages were observed in the northern neighborhoods, also called as the 4th district, especially São Geraldo (301.67K),Floresta (301.66K),Cristo Redentor (301.37K) and Vila João Pessoa (301.35K).There are two main factors that stand out for the higher LST in Porto Alegre, as also observed by [24]: 1) physically, they are related to the installation of surface train facilities in April 1985 that goes out to the north, spreading to the adjacent neighborhoods, with a further enhancement of impervious surfaces and avenues; 2) economically, there is a concentration of most industrial plants, corporate services and companies located in these regions.Furthermore, UHI mitigation strategies involving vegetation tend to be more expensive per unit area than strategies involving high-albedo surfaces [76] as for roof tops and industrial facilities.These two factors are interrelated to the increase of pollution, population decline and imbalance in land use between impervious surfaces and vegetation coverage.
Regarding the lower LST cluster at the middle of the study area, we found that are mostly covering residential areas.Lower LST values were observed in the center-southern neighborhoods (orange and dark red), especially the Farroupilha (298.65K),Bela Vista (298.68K),Moinhos de Vento (299.17K) and Petrópolis (299.40K).There are two main factors that stand out for the lower LST in Porto Alegre: 1) physically, they are distributed as a concentration of residential ones; 2) economically, there is a concentration of parks and malls located within these neighborhoods.
Figure 2 shows the LST averages covering the urban area of Porto Alegre city obtained from all selected images in the study area.

Analysis of surface physical descriptors
Each neighborhood in Porto Alegre city has different characteristics caused by the patterns of the prevailing urban development over the past decades.Different combinations of VIS play an important role in the LST distribution.So, an analysis of the connection of VIS surface cover types with LST distribution was performed.For simplicity, correlations were stablished by using a negative first order polynomial function as described by the mathematical equation y=−β×x+γ.
Our results show that the simple analysis of %Vegetation indicates a poor correlation to LST averages (R 2 =0.23) as expected, and %ISA (R 2 =0.31).However, %Soil coverage plays a major role possibly because it is indirectly related to vegetation coverage and the soil moisture associated to soil background in green areas.This result is probably linked to the fact that the empirical model used in this study for VIS surface type classification considered sparsely vegetation coverage and non- impervious surfaces as the same set of soil surface coverage type.Even so, although %Soil (or sparsely vegetation presence) does show a better correlation to LST averages distribution, it does not express any spatiality regarding its values as the UHI patterns recognized in Figure 2.
In general, we have found poor correlation of physical variables VIS, Altitude or Population density to LST as shown in Figure 3.Although studies from [30] identified that LST is better correlated to percent of impervious surfaces than vegetation index, the physical conditions and properties of the vegetation types and surrounding urban environment also have important influence on determination of LST distribution [77] as well for soil surface area.As it can be seen in Figure 3(B) and 3(C) results show that %ISA increases due to urbanization and subsequently leads to %Vegetation decrease according [78].Even so, it is possible to find some dense canopy trees in areas of almost completely impervious surfaces at surface level.Regarding low correlation of %Vegetation and Vegetation density (EVI-2) to LST averages, the main reason probably is linked to the fact that vegetation density by itself does not represent value on the real estate configuration.In this way, higher surface vegetation coverage (%Vegetation) are not necessarily representative of green areas and structured and humanized urban green spaces.Based on the emissivity dependence of LST to EVI-2, we applied the EKC concept in order to test the correlation between the welfare due to the social-environmental arrangement of the citizens (Household Income) and LST distribution at neighborhood scale.On the analysis of LST distribution associated to HI, a better correlation was found by using a negative first order polynomial function as described by the mathematical equation y=−β×x+γ in Figure 4. On the neighborhood scale, results show a linear inverse correlation (R²=0.36) between HI and LST averages.Physically, it is expected that there is a superior threshold to those higher LST averages as function of surface coverage conditions.

Urban development modelling
The different urban microclimates show alterations that stem from the loss of vegetation cover combined with soil imperviousness that come from the density and constructive compactness which are directly related to the rising socioeconomic development of these areas.So, in order to test the relation observed between the environmental variable to the economical one, we analyzed the spatially distributed dependence of LST averages to HI by means of a mathematically combination of the ratio (as Household Income/LST average) because they are inversely related.We obtained a rank which is adequately associated to the spatial distribution of UHI clusters in the city of Porto Alegre.Results in Figure 5(A) and 5(B) indicate that the ratio between HI and LST averages, increases the rank when LST is lower and HI is high, which can be considered as an indicator of better urban development conditions.The ratio decreases when LST averages are higher and HI is lower which, in the same way, can be considered as representative of not well developed urban areas.The top five higher ranked neighborhoods are: 1) Três Figueiras, 2) Bela Vista, 3) Moinhos de Vento, 4) Chácara das Pedras and 5) Mont Serrat, respectively.The top five lower ranked neighborhoods are: 6) Vila João Pessoa, 7) Bom Jesus, 8) Sarandi, 9) Cel.Aparício Borges and 10) São José, respectively.

Discussion
These ranked results, which are essentially the same as HI alone, suggest that there is an economical influence on the historical formation of neighborhoods.The main lower ranked cluster and its adjacent areas are characterized by a particular and historical behavior in the urban environment development.Vila João Pessoa (between São José and Cel.Aparício Borges) is the first suburban neighborhood, formed around 1940, and developed along with the same sub-standard conditions through decades without adequate planning [79].Most of Vila João Pessoa population still live in a society with no fixed job, and survive on the collection of waste (industrial waste, plastic bottles).In 2010, the average income of residents of Vila João Pessoa was 3.23 salary units.The average for the entire Porto Alegre city, with the 83 neighborhoods was 5.3 salary units ranging between 1.84 to 18.24 salary units [79].In January, 4th 2010, the value of 1 salary unit (monthly payment) was on average equivalent to US$296.5 (BACEN) [80].
Considering the main higher ranked cluster and its adjacent areas, Moinhos de Vento (Windmills) is characterized by a particular historical behavior in the urban environment development.It has been formed right after arrival of the Açorianos (migrants from the Açores Islands region, of Portugal), in 1772.The windmills were brought by the Açorianos and had its apogee during the wheat cultivation period.Even so, after 1835 the cultivation in the region was abandoned.Neighborhood development was boosted in 1893 with tram line installation, implemented for public transportation by urban municipality company "Carris".Another driving factor for the Moinhos de Vento development was the construction of hydraulic facilities in the hills of the area, in 1904, which led to the opening of several nearby streets.In the southern border of the neighborhood area the German Hospital took place in 1927, which in 1942 was renamed Hospital Moinhos de Vento [79].Still today is one of the best and most modern city facilities.In 2010, the average income of residents of Moinhos de Vento was 16.05 salary units.Nowadays, it exhibits high-standard conditions and with adequate urban planning maintaining higher well-infrastructured conditions of green spaces and residential areas.Such results suggest that the economic resulting geolocation is a function of a veiled and underlying sectoring of city areas which, by its turn, promotes an environmentally unsustainable configuration of sub-urban areas but which is not directly related to imperviousness, vegetation coverage, soil exposure or lot size.
The spatial distribution of the districts and villages in the territory of the city, in general, tend to reflect the distribution of the various social segments.Peripheral squatter settlements, either legal or clandestine, devoid of infrastructure, in different parts of the city, can present higher environmental temperature rates due to denser soil occupation.Ultimately, such macro-scale results observed from satellite imagery indicate that there is need for environmental maintenance in the lower ranked neighborhoods.These results indicate that not only the prevailing physically driven conditions in the urban environment but also the social environmental arrangement associated mainly to infrastructure and aesthetically configuration of neighborhoods, is exerting a non-negligible influence on LST distribution.This also suggests that reduced quality of public administration and management policies thorough years has been contributing to the precarious environmental situation of some neighborhoods, which in its turn can also lead to the strengthening of low educational status of individuals [79].

Conclusions
Results suggest that the economically-related scenario exerts a non-negligible influence on the physically driven characteristics of the urban environmental conditions as predicted by EKC theory.
Different combinations of VIS surface cover types play an important role on the LST distribution.However, the characteristics of the urban environmental conditions that can be detected by thermal remote sensing are linked to the aspects which are associated to different levels of social organization in the urban contexts that cannot be objectively measured.Examples are: public administration efficiency, welfare, education or even sustainable urban environment concept which also depends on social and economic standards.In this way, the investigation on the association between the HI and LST averages at neighborhood scale has been demonstrated to be straightforward for the understanding of the development process that leads to UHI.This study will also be available as documentation of the existing findings on UHI in Porto Alegre city and thus can help further research.
Finally, we conclude that even at the neighborhood scale, there are variables regarding social organization, public management and city design that should be considered when studying urban development and environmental conditions.Future research may investigate the relationship between other economically-related variables to specific land surface characteristics.

Figure 1 .
Figure 1.Urban area of Porto Alegre in Brazil and the shaded study area with 47 neighborhoods.

Figure 2 .
Figure 2. Average of LST between 2001 and 2010 and estimated UHI in Porto Alegre municipality with the neighborhood limits highlighted.

Figure 3 .
Figure 3. Scattergram comparing the estimated LST distribution considering the 47 neighborhoods in the urban area of Porto Alegre municipality to A) average of EVI-2; B) %Vegetation; C) %ISA; D) %Soil; E) Altitude and F) Population density.

Figure 4 .
Figure 4. Scattergram comparing the estimated LST averages distribution to Household Income in the urban area of Porto Alegre municipality.

Figure 5 .
Figure 5. Average of LST in Porto Alegre municipality between 2001 and 2010 showing the location of development clusters obtained by ranking the ratio LST averages/HI with the streets design in the background, where: (A) Top 5 higher ranked neighborhoods with the main cluster highlighted and Moinhos de Vento limited by the dotted-red line, and (B) Top 5 lower ranked neighborhoods with the main cluster highlighted and Vila João Pessoa limited by dotted-red line.

Table 1 .
Selected Landsat-5 TM imagery for monitoring the studied area.

Table 2 .
Summary of retrieved averages obtained from the selected Landsat-5 TM.