-In the wake of climate change, cities need to adapt to global warming. In this context, the use of afforestation to improve the microclimate may assist in raising the quality of life for population. This objective requires research that analyzes how the variations in parameters related to canopy dynamics, such as the leaf area index (LAI) and photosynthetically active radiation (PAR) can influence thermal comfort indices. To contribute to this research, this study measured the air temperature, relative air humidity, PAR, and LAI on a monthly basis from July, 2017, to June, 2018, in an urban park in a tropical region of Brazil. Kriging maps were created for the heat index (HI), and multiple polynomial regression models were adjusted to estimate the HI using PAR and LAI data. After defining the models, positive and negative variations of LAI were tested to observe if any changes in HI occurred. The simulated results showed greater sensitivity to negative variations in LAI, in which a 50% reduction in LAI decreased the HI by 28%, particularly during the dry period. As the area is a region that experiences considerable variability in terms of humidity and high temperatures throughout the year, conditions conducive to thermal discomfort had often occurred in the study area. This indicates that for tropical regions, even in urban parks, it is important that users pay close attention to their hydration needs and the duration and intensity of their physical activities.
Keywords: afforestation; air temperature; mobile transect; photosynthetically active radiation; relative humidity; thermal comfort. limate change and its impacts associated with anthropogenic actions including vegetation suppression, land use, and occupation, and the dense and compact set of constructions from constant urban expansion have become the focus of widespread discussion amongst the scientific community (Grimmond et al., 2009;Middel et al., 2014;Petralli et al., 2014). This is due to the changes in the thermodynamic field caused by the decrease in shortwave reflection and increased emission of long waves; this causes an increase in air temperature even during the periods of shorter heat stroke duration, intensifying the formation of heat islands (Ayoade, 2003;Peng et al., 2012).
Artificialized areas, mainly in the central regions of cities, produce greater changes in the local climate. Consequently, vegetated urban spaces contribute to better thermal comfort and the reduction of heat islands. This is because of the interception of a part of the solar radiation incident by the tree canopy; the extent of interception varies based on the species, resulting in differing attenuation of solar radiation. (Abreu et al., 2012). Herb et al. (2008) affirmed that the canopy of trees affects the heat transfer of the surface and the temperature of the soil below it, providing better thermal conditions for pedestrian movement (Souza et al., 2020).
Thus, the variations in leaf area indices (LAI) may influence thermal comfort, especially in tropical regions with seasonal precipitation. According to Llandert (1982), thermal comfort associated with the presence of afforestation is mainly related to the canopy density of each tree. Through its leaves, trees are able to absorb 15% to 35% of the received light energy, pass between 30% and 50% of the energy, and reflects the remainder of the energy (~30% to 40%), during the daytime.
The constant concern on promoting the quality of life in cities and the health of the population has fostered the use of research and different methodologies, to demonstrate the effective action of vegetation on the urban microclimate (Bueno- Bartholomei, 2003 Ren et al. (2013) confirmed that the population recognizes urban parks as an oasis amid the dense urban construction model, acting to alleviate the thermal discomfort of heat islands.
One way to measure a user's comfort index is the heat index (HI) method proposed by Steadman (1979aSteadman ( , 1979bSteadman ( , 1984)). This method is one of the most popular environmental health indices, providing the basis for heat warnings in many communities of the United States of America (NOAA, 2009). The heat index has also been applied in several regions and climates in
This research was conducted at the Conservation Unit Ilto Ferreira Coutinho Park, located in the central region of the municipality of Tangará da Serra, state of Mato Grosso, Brazil (Fig. 1). The study area is located at 14°37'08"S and 57°29'09"W, with an altitude of 452 m and spans approximately 12 ha (Melz and Tiago, 2009). The original growth of semideciduous forest is characterized as Cerrado-Amazonian ecotone. The most abundant families are Anacardiaceae, Malvaceae, Bignoniaceae, Annonaceae, Apocynaceae, Meliaceae, Miristicaceae, and Rubiaceae, where 81% of the species are native to the region (Rodrigues, et al., 2015). There are two types of soils in the study area: red dystrophic latosol (red latosol) and hydromorphic quartzenic neosol. Based on the Köppen climate classification, the climate is Aw which is characterized as hot and humid with rain in the summer and drought in the winter with wet and dry seasons (Alvares et al., 2013). The surroundings are composed of commercial areas and single-story residences, with very little shade generated by buildings.
Period Air temperature and relative humidity were collected for one year, on a monthly basis, Data collection commenced in July 2017 and ceased in June 2018, from 8 am to 5 pm, with an hourly collection, at a central park point; the averages were subsequently calculated. Collection days were preferentially chosen when there was no rain or cloud, based on the method described in Oke (1982). Measurements were obtained using a portable microclimate station (Kestrel 4500 Weather Tracker, NK Company, Boothwyn, PA, USA).
The incident photosynthetically active radiation (PAR) data and LAI were collected at 12 pm at 15 points distributed throughout the park; the averages of these values were calculated following the measurement. This data was measured using a linear ceptometer (AccuPar -LP 80, Decagon Devices, Washington, USA).
The mobile transect method was used for geostatistical analysis, with air temperature and relative humidity data collected every minute during a trip that covered all regions of the park. The collections were made at 8am, 12am, and 5pm , in January and June 2018, to cover the seasonality of the region. The HI was calculated based on these results, and the semivariograms and kriging map were obtained. Because the sensor is open path, not aspirated, it was also important that measurements occurred when the air was not completely stagnant during the study to minimize lags in sensor response (Sun, 2011).
The HI was calculated for thermal comfort analysis, often used for hot regions with low-intensity winds and where an individual is in the shade. This results in the body thermal sensation based on the air temperature and relative humidity (Steadman, 1979b).
The HI was proposed by Steadman (1979a), according to Equation ( 1):
HI = -42,379 + 2,049015230 * T air +10,14333127 * R h -0,22475541 * T air * R h -6,83783 * 10 -3 * (T air ) 2 -5,481717 * 10 -2 * (R h ) 2 + 1,22874 * 10 -3 * (T air ) 2 * R h +8,5282 * 10 -4 * T air * (R h ) 2 -1,99 * 10 -6 * (T air ) 2 * (R h ) 2 (1)where HI is the heat index (°F), T air is the actual air temperature (dry bulb temperature) (°F), and R h is relative humidity (%). The data in ºF was subsequently converted to ºC. The HI results were analyzed according to alert levels detailed in Table 1.
The T-test test was used to verify possible statistically significant differences to air temperature, relative humidity, and HI, at a significance level of 5%. This was carried out to verify possible differences between the periods, and the Pearson's correlation was used to identify dependencies between meteorological variables.
For spatial analysis, semi-variograms were adjusted to the results of the HI. The semi-variance measures the degree of dependence between two samples. It increases as the distance between points increases, until it stabilizes at the point known as the threshold (Co + C1), with half the hope of variance between the pairs of points separated by a distance "h" represented by the classic model, according to Equation (2) as follows:
??(?) = 1 2??(?) ? [??(?? ?? ) ? ??(?? ?? + ?)]² ??(?) ??=1(2)Where ? (h) is the estimator of semi-variance for each distance, h; N(h) is the number of pairs of points separated by the distance, h; Z(x) is the regionalized variable at point, x; and Z(x+h) is the value of point x+h (Burrough and Macdonnell, 1998). The semi-variogram is represented by the plot of h versus h. The theoretical semi-variogram generated by this function must be adjusted with a theoretical model that provides the parameters nugget effect (Co), sill (Co + C1), and range (Ao). The degree of spatial dependence of variables was classified according to Cambardella et al. (1994); there was strong spatial dependence when semi-variograms had a nugget effect that was 25% from the threshold, moderate when between 25 and 75%, and weak when greater than 75%. The theoretical semi-variogram models consisted of spherical, exponential, and Gaussian models described by Andriotti (2003) and Yakamoto and Landim (2013). These models were used to estimate the semi-variance at any distance between samples using the developed GS + (GS +, 2000) software by Gamma Design Software and spreadsheets. The selection of most optimal semivariogram adjustment method was important as this is the reference point from which the spatial correlation structure to be used in the inferential kriging procedure is interpreted (Dias et al., 2015).
The evaluation and selection of most optimal adjustments of experimental semi-variograms were based on the smallest sum of the squares of weighted deviations (SQDP) and the highest coefficient of determination (R²; Andriotti, 2003; Yakamoto and Landim, 2013). The interpolation and spatialization of the variables were then carried out using the punctual ordinary kriging; their spatial distribution was subsequently analyzed.
For models estimating the HI as a function of PAR and LAI, the multiple polynomial regression method was applied, using the Solver tool of Microsoft Excel®. This tool uses the generalized reduced gradient (GRG) as the solution method, adjusting five parameters (A, B, C, D, and E) as per Equation (3):
???? ?????? = (?? * ?????? ?? ) + (?? * ?????? ?? ) + ?? (3)After determining the adjusted equation, a t-test was conducted to compare the calculated and modeled data. The LAI was varied positively and negatively by 10%, 20%, 30%, 40%, and 50%, to enable projections on how LAI variation may interfere with HI.
Hourly averages were determined for air temperature, relative humidity, and HI data for mobile transects conducted in dry and rainy periods, as shown in Fig. 2. For the wet season, 12 h and 17 h were classified as "Extreme Caution", according Table 1, and 8 h as "Caution". The greatest discomfort observed occurred at sunset at 17 h, reaching 34.52 ±1.09 ºC.
The wet season is the period has the highest incidence of solar radiation, as it is summer in the Southern Hemisphere; this equates to higher temperatures and relative humidity. During the dry period, the average Possible edge effects may contribute to the greatest standard deviations observed in during the mornings; this was a variation of 2.62 ºC for the wet season and 1.99 °C for the dry season. The regions outside the park were observed to heat up relatively rapidly due to the specific heat of the building materials. This was an occurrence that had been attenuated during other times as the park heats up.
For seasonality analysis, the difference in median values between the two groups for air temperature, relative humidity, and HI was greater than expected, a statistically significant difference was observed (P < 0.001). Thus, the seasonality results follow local weather patterns.
Spatial analyses were conducted to observe the distribution of the HI in the park and possible edge effects. For this purpose, semi-variograms were calculated according to Table 2.
Table 2: Adjustment of semivariograms for heat index: model, nugget effect (Co), sill (Co + C), range (Ao), determination coefficient (R²), and sum of squares of residuals (SSR) and spatial relationship.
The largest nugget effect found (Co) was 0.57, and relatively small values of the nugget effect indicated minor errors in measurements (Dafonte et al., 2010). With the exception the 12 am measurement during the wet season in which the adjustment was exponential, all other adjustments were spherical. Of the six schedules analyzed, there was a strong spatial dependence in four; according to the classification of Cambardella et al. (1994), there is only strong spatial dependence when the nugget effect is up to 25% of the sill. Moderate dependency occurred twice, between 5 pm during the wet season and 12 am in the dry season.
The greatest range (A 0 ) occurred at midday in the dry period, with 501.5 m, generating the observed similarity in the special HI, as shown in Fig. 3. The lack of humidity and diminished leaf area accounted for the standardization observed in the dry season, as shown in Table 2 and Fig. 3. The range is of fundamental importance to interpret the semi-variograms, indicating the distance to where sampling points correlate with each other (Carvalho et al., 2002). This ensures that all neighboring points are so similar that they may be used to estimate values for any point between them (Machado et al., 2007).
The lowest coefficient of determination was 0.78, at 12am during the wet season, suggesting good model adjustments. Regarding the sum of the squares of residues (SSR), only the 8 am measurement during the wet season had values exceeding. As the SSR is a measure of discrepancy between the actual and modeled data, a small SSR indicates a tight fit between the model and the measured data (Draper and Smith, 1998).
The spatial patterns follow the hourly patterns shown in Fig. 2, in which in the spatial analyses of the wet season showed that thermal discomfort was greater in this season than the dry season. Dense vegetation cools the air as it prevents radiation from reaching the forest floor (Napoli et al., 2016) and provides a larger surface area for evaporative cooling. However, the incidence of radiation in the summer season was able to elevate air temperatures, resulting in higher heat levels during the wet season. According to the Pearson's correlation test (r), air temperature and relative humidity were negatively correlated (r = -0.666 and p > 0.01). This means that an increase or decrease in air temperature results in an increase or decrease in relative air humidity. Varejão-Silva (2006) explained that such behavior is caused by the inverse proportionality of the relative air humidity and the vapor saturation point. Therefore, the relative air humidity will also be inversely proportional to the air temperature.
The annual average LAI was approximately 6 m²m -2 , demonstrating regional seasonality. These values are close to those reported by Sanches et al. (2008) who used the Moderate Resolution Imaging Spectroradiometer (MODIS) satellite and reported values between 5.25 to 5.54 m²m -2 for another area of the Cerrado-Amazonian ecotone forest. As the urban park vegetation is a seasonal semi-deciduous forest, leaf loss occurs from August onwards, decreasing the LAI. Another factor that contributes to lower LAI is the solar declination angle; during the dry season this angle is higher (Spolador et al., 2006;Novai et al., 2018). The greater the zenith angle, the greater the path traveled by radiation within the canopy, increasing the chance of absorption by leaves and branches (Senna et al., 2005), this increase is associated due to the measurement technique, which uses the linear ceptometer. PAR data followed the radiation pattern of the Southern Hemisphere, according to Novais et al. (2016), and LAI and PAR were strongly correlated (r = 0.804 and p > 0.01), to the maximum incidence of PAR in January, approximately 2006 mol m²s -1 .
The accumulated rainfall throughout the analyzed period was 1968 mm; in July 2017, there was
no rain whereas February was the rainiest with 431 mm, which is approximately 22% of total annual precipitation. The seasonal pattern of precipitation in the study area is in line with the regional forecast as rainfall is concentrated between October and April in the Cerrado region (Dallacort et al., 2011).
???? ?????? = (51.63 * ?????? ?1. 19 ) + (107.57 * ?????? 0.11 ) ± 216.29 (4) The difference between the mean values of HI real and HI mod was insufficient to reject the possibility that this difference was due to random sampling variability. There was no statistically significant difference between the groups (P = 0.996), justifying the use of the model to simulate LAI variations.
Table 3 presents the annual averages of the HI for each model variations and their respective differences from the original value. The 50% variation of LAI showed the most concerning results, where three months were within the "Danger" classification, and the remaining months were in the "Extreme Caution" classification. The positive variations in LAI contributed to a greater number os months in the "No Warning" classification, ranging from zero months in this classification to four months for the LAI +50% model.
Based on the analysis of hourly average air temperature and relative humidity generated by the spatialized data, the wet season was more uncomfortable than the dry season. This result differs from the monthly analyses undertaken at a central point of the park. Such differences in results suggest that the spatial patterns of data collection as well as the possible influence of clearings and edges, influenced these The addition of 50% LAI improved the annual average HI by approximately 8.32%, whereas the decrease resulted in a 27.74% decrease in the average annual HI. Fig. 5 presents the monthly HI results from the model variations according to the LAI.
Using PAR and LAI data, the following model was obtained for HI estimates: results. These results corroborate the scientific findings from Silva Júnior et al. (2012) for the city of Belém (Pará, Brazil), whereby for the majority of the afternoon the city is thermally uncomfortable. Neighborhoods with a higher percentage of soil sealing and lower vegetation cover exhibited greater thermal discomfort, indicating that urban elements (paved streets, houses, buildings, and vehicle routes), influence the city's climate.
In terms of the influence of the edge effect, Vasconcelos and Zamparoni (2011) and Morakinyo et al. (2017) found that the effectiveness of trees to improve daytime thermal comfort reduces with increasing urban density whereas the opposite was true for the nighttime. Therefore, as the park is not influenced by external shading due to low construction density and people do not visit the park at night, there is no problem using species that have higher LAI, as it is expected that only daytime comfort occurs. This highlights the importance of urban planning to present solutions to promote quality of life within a healthy environment.
The most uncomfortable time was at 5 pm during the wet season, suggesting that despite the high LAI of this period, which mitigates the incidence of solar radiation, air exchange is reduced, making heat dissipation difficult. The cooling effect of vegetation is likely to be lower at night when there is very little transpiration. This means less evaporative cooling (Richards et al., 2020) because of the trapped heat and humidity within the urban canopy layer, compared with the rapid nocturnal cooling of open areas (Fahmy et al. 2010). Boone (2008) explains that humidity often creates a temperature that feels hotter than reality. This is because the body cools with the evaporation of sweat due to the consumption of latent heat on the surface of the skin. However, when air humidity rises, there is a decrease in the rate of sweat evaporation, causing greater heat retention and leading to discomfort and stress (Delworth et al., 1999).
In terms of the analysis of monthly averages, the temperature decreases in September (the dry season), may be related to cold fronts, where there are sharp drops in solar incidence due to the presence of cloudiness (Biudes et al., 2015).
For model analysis, the annual average LAI was approximately 6 m²m -2 ; this is a relatively high value compared to cerrado forests, where Hoffman et al. Thus, whilst dense canopy does contribution to thermal comfort, its contribution is not as great as the decrease in LAI, which leads to greater variations in HI, making the HI reach 40.25 ºC (extreme Caution Alert Level); this constitutes a ~28% increase in relation to measured values.
Despite the general consensus that urban parks are the places of leisure and exercise as well as oases within the urban area, in tropical areas, greater attention should be given to vulnerable users of the park, for instance, those susceptible to cramps, sunstroke and likely exhaustion, and possible brain damage due to prolonged exposure to more intense physical activities. This is because the heat indices during the dry season was largely unsatisfactory.
V.
The LAI influences the HI, and this influence is more pronounced in models that decrease the amount of leaves. A 50% reduction in LAI caused a 28% increase in the HI, being classified as "Extreme Caution" alert level during the dry season which was found to be more susceptible to the changes in the LAI. Hourly analysis showed that the evening, especially in the wet season, was of greater concern, as there was a lack of thermal comfort even in shaded areas. The possible nighttime discomfort generated by the heat retained in the park is insufficient to cause major concern, as park visit is typically limited to the daytime. Although urban parks are recognized as a place for leisure and physical activities, greater care is suggested during the use of urban parks in tropical regions. It is recommended that park users equip themselves with light clothing, hydration, and are careful with strenuous activities, avoiding the periods of greatest thermal discomfort.
Alert level | Heat index (HI) °F °C | Symptoms | |
No warning | < 80 | < 27 | No problem. |
Caution | 80-90 | 27-32 | Possible fatigue in cases of prolonged exposure and physical activity. |
Extreme caution | 90-105 | 32-41 | Possible cramps, sunstroke, and exhaustion due to prolonged exposure and physical activity. |
Danger | 105-130 | 41-54 | Cramps, sunstroke, and likely exhaustion. Possible brain damage due to prolonged exposure to physical activity. |
Extreme danger | > 130 | > 54 | Stroke or action and risk of imminent cerebral vascular accident. |
0% | 10% | 20% | 30% | 40% | 50% | ||
HI | LAI + LAI - | 31.504 | 30.77 32.42 | 30.17 33.58 | 29.67 35.11 | 29.25 37.21 | 28.89 40.25 |
Variation | LAI + LAI - | - | -2.34% 2.88% | -4.24% 6.58% | -5.83% 11.44% | -7.17% 18.11% | -8.32% 27.74% |
This work was supported by the Mato Grosso Research Support Foundation-FAPEMAT [grant number 0194288/2017].
Impacto da morfologia de parque urbano no microclima e no conforto térmico de Cuiabá-Brasil. Revista Brasileira de 2019.
Impact of urban form and design on mid-afternoon microclimate in Phoenix Local Climate Zones. 10.1016/j.landurbplan.2013.11.004. Landscape and Urban Planning 2014. 122 p. .
Field-scale variability of soil properties in central Iowa soils. 10.2136/sssaj1994.03615995005800050033x. Soil science society of America journal 1994. 58 (5) p. .
An evaporation-interception model for urban areas. 10.1029/91WR00557. Water Resour. Res 1991. 27 p. .
Urban surface energy balance models: model characteristics and methodology for a comparison study. 10.1007/978-3-642-00298-4_11. Meteorological and Air Quality Models for Urban Areas, (Berlin, Heidelberg
A street thermal environment study in summer by the mobile transect technique. 10.1007/s00704-011-0444-6. Theoretical and applied climatology 2011. 106 (3-4) p. .
Propriedades físico-químicas e microbiológicas do solo de um Parque em Tangará da Serra, MT, uma área de transição entre Amazônia e Cerrado. 10.1590/S0044-59672009000400011. Acta Amazonica 2009. 39 (4) p. .
Análise da distribuição espacial do conforto térmico da cidade de Belém, PA no período menos chuvoso. 10.26848/rbgf.v5.2.p218-232. Revista Brasileira de Geografia Física 2012. 5 (2) p. .
Mapping of soil micronutrients in an European Atlantic agricultural landscape using ordinary kriging and indicator approach. 10.1590/S0006-87052010000500018. Doi:doi.org/10.1590/S0006-87052010000500018 Bragantia 2010. 69 p. .
Geoestatística na determinação da variabilidade espacial de características químicas do solo sob diferentes preparos. Doi: 10. 1590/S0100-204X2002000800013. Pesquisa Agropecuária Brasileira 2002. 37 (8) p. .
Radiação fotossinteticamente ativa em uma floresta de transição Cerrado-Amazônica. Revista Brasileira de Meteorologia 2006. 21 (3b) p. .
Albedo do Solo abaixo do Dossel em Área de Vochysia divergens Pohl no Norte do Pantanal. 10.1590/0102-778631220150001. Revista Brasileira de Meteorologia 2016. 31 (2) p. .
Space-temporal variation of PAR reflected by the soil and transmitted by the canopy in a floodplain forest of Mato Grosso State. Doi: 10.5902/ 1980509835097. Brazil. Ciencia Florestal 2018. 28 (4) p. .
Áreas verdes urbanas: uma contribuição aos estudos das ilhas de frescor. Brazilian Geographical Journal: geosciences and humanities research medium 2015. 6 (2) p. .
Atmospheric Pollution and Meteorological Parameters in the City of Cuiabá-MT. 10.15341/mese(2333-2581)/06.05.2019/003. Modern Environmental Science and Engineering 2019. 5 (6) p. .
Índice de área foliar em floresta de transição amazonia cerrado em diferentes médotos de estimativa. Doi: 10. 5902/2179460X9750. Ciência e Natura 2008. 30 (1) p. .
Bioclima térmico como fator de planejamento urbano e arquitetônico em climas tropicais -o caso de Campinas. Doi: 10.1007/ s00704- 013-0886-0. Brasil. Urban Ecosystems 2014. 17 (2) p. .
Fração de radiação fotossinteticamente ativa absorvida pela floresta tropical amazônica: Uma comparação de medições de campo, modelagem e sensoriamento remoto. Journal of Geophysical Research 2005. (G1) p. 110. (Biogeosciences)
LAI based trees selection for mid latitude urban developments: A microclimatic study in Cairo. 10.1016/j.buildenv.2009.06.014. https://doi.org/10.1016/j.buildenv.2009.06.014 Egypt. Building and Environment 2010. 45 (2) p. .
Modeling tree shade effect on urban ground surface temperature. Doi: 10.2134/ jeq2015.02.0097. Journal of environmental quality 2016. 45 (1) p. .
Urban planning indicators: useful tools to measure the effect of urbanization and vegetation on summer air temperatures. 10.1002/joc.3760. International Journal of Climatology 2014. 34 (4) p. .
Patterns of energy exchange for tropical ecosystems across a climate gradient in Mato Grosso. Doi: 10.1016 / j.agrformet.2014.12.008. Brazil. Agricultural and Forest Meteorology 2015. 202 p. . (J.)
Estimativa do conforto térmico na cidade de Cuiabá/ MT. 10.18316/1981-8858.16.22. http://dx.doi.org/10.18316/1981-8858.16.22 Revista de Ciências Ambientais 2016. 10 (1) p. .
Métodos de ajuste de semivariogramas para modelagem espacial de íons de precipitação pluvial em Cuiabá. 10.5902/2179460X12101. Brasil. Ciência e Natura 2015. 37 (2) p. .
Soil nutrients and vegetation structure in a neotropical seasonal wetland. Doi: 10. 15666/aeer/1302_289305. Applied ecology and environmental research 2015. 13 (1) p. .
Inventário florístico no entorno das trilhas do parque natural municipal "Ilto Ferreira Coutinho. Anais: Resumos Expandidos. Vi Simpósio Da Amazônia Meridional Em Ciências Ambientais Scientific Electronic Archives, 2015. 8.
Distribuição das chuvas no município de Tangará da Serra, médio norte do Estado de Mato Grosso. Doi: 10. 4025/actasciagron.v33i2.5838. Brasil. Acta Scientiarum. Agronomy 2011. 33 (2) p. .
Ground cover and its influence on temperature of urban microclimate of Cuiabá-MT. Caminhos de Geografia 2011. (39) p. 12.
The assessment of sultriness. Part I: A temperature-humidity index based on human physiology and clothing science. 10.1175/1520-0450(1979)018<0861:TAOSPI>2.0.CO;2. Journal of applied meteorology 1979. 18 (7) p. .
The assessment of sultriness. Part II: effects of wind, extra radiation and barometric pressure on apparent temperature. 10.1175/1520-0450(1979)018<0874:TAOSPI>2.0.CO;2. Journal of Applied Meteorology 1979. 18 (7) p. .
O microclima e o (des) conforto térmico em ambientes abertos na cidade do recife. Revista de Geografia (Recife) 2011. 28 (1) p. .
Urban microclimate in vegetated and non-vegetated areas in rainy and sunny conditions. 10.5380/raega. Raega-O Espaço Geográfico em Análise 2020. 47 (1) p. .
Tendências das séries históricas do índice de calor no município de Santa Maria-RS. 10.1590/0103-8478cr20131345. Ciência Rural 2014. 44 (8) p. .
Surface urban heat island across 419 global big cities. 10.1021/es2030438. Environmental science & technology 2012. 46 (2) p. .
Temporal dynamic of CO 2 efflux in a savanna fragment in Mato Grosso. 10.6008/CBPC2179-6858.2018.002.0004. Revista Ibero-Americana de Ciências Ambientais 2018. 9 (2) p. .
A study on the impact of shadow-cast and tree species on in-canyon and neighborhood's thermal comfort. Doi: 10.1016 / j.buildenv.2017.01.005. Building and Environment 2017. 115 p. .
The energetic basis of the urban heat island. Doi: 10. 1002/qj.49710845502. Quarterly Journal of the Royal Meteorological Society 1982. 108 (455) p. .
Seasonal leaf dynamics across a tree density gradient in a Brazilian savanna. 10.1007/s00442-005-0129-x. Oecologia 2005. 145 (2) p. .
Simulação da temperatura da superfície do solo para diferentes coberturas do solo. 10.1016/j.jhydrol.2008.04.020. Journal of Hydrology 2008. 356 (3-4) p. .
Estimation of the relationship between urban park characteristics and park cool island intensity by remote sensing data and field measurement. Doi: 10.3390/ f4040868. Forests 2013. 4 (4) p. .