Abstract

The savanna of Central Brazil (locally known as cerrado) has a long history of land cover change due to human activity. These changes have led to the degradation of cerrado forests and woodlands, leading to the expansion of grass-dominated cerrados and pastures. Thus, the aim of this study was to evaluate the temporal variation in energy flux in areas of degraded, grass-dominated cerrado (locally known as campo sujo) in Central Brazil. The amount of partitioned into H declined as monthly rainfall increased and reached a level of approximately 30% during the wet season, while the amount of partitioned into increased as monthly rainfall increased and reached a level of approximately 60% during the wet season. As a result, H was significantly higher than during the dry season, resulting in a Bowen ratio (β = H/) of 3-5, while Le was higher than H during the wet season, resulting in a . These data indicate that the energy partitioning of grass-dominated cerrado is relatively more sensitive to water availability than cerrado woodlands and forests, and have important implications for local and regional energy balance.

1. Introduction

Tropical savannas cover about 12% of the global land surface [1] and are characterized by high plant species diversity [2]. In Brazil, savanna (locally known as cerrado) covers about 24% of the territory, mostly in the central portion, and is the dominant vegetation in areas where the dry season causes prolonged periods of plant water stress [35]. The biodiversity of cerrado is extremely high and is estimated to be 160,000 species, including known plants, animals, and fungi [6].

Over the last few decades, cerrado has been converted to cattle pasture, and more recently, soybean and sugar cane agriculture [2, 7], and has experienced deforestation rates much higher than in the Amazon rainforest [8]. Another major threat to the remaining areas of cerrado is the decrease of the woody component due to the increase of anthropogenic fire frequency [9], converting the vegetation to a more open and shallow-rooted ecosystem. These changes have the potential to cause multiple changes in the structure and function of cerrado [7, 10, 11]; however, this biome has received relatively little attention from researchers in comparison with tropical rainforests [12, 13].

The change in land cover has the potential to change energy partitioning, by affecting the seasonal pattern and magnitude of radiation balance and albedo [12], and energy partitioning in the form of latent [1416], sensible [17, 18], and soil heat flux [19], which, in turn, will feedback on local, and perhaps regional, climate [20, 21]. For example, tropical forest conversion to pasture can cause a 1.5–2.0 kPa increase in vapor pressure deficit and a 5–10°C increase in soil surface temperature relative to intact forest [22]. Land cover change can also lead to an increase in the duration of the dry season, cause more rainfall to be partitioned into runoff, affect the development of the nocturnal and convective boundary layer, and destabilize regional rainfall regimes and surface water availability [12, 20].

Given the potential for land cover change to alter surface energy balance in cerrado, we evaluated the seasonal and interannual variations of energy partitioning in a degraded grass-dominated cerrado (locally known as campo sujo) of Central Brazil. We used Bowen ratio energy balance (BREB) methods over two consecutive years to characterize the seasonal and interannual variations in energy flux dynamics. We hypothesized that campo sujo cerrado would exhibit higher rates of sensible and ground heat flux than latent heat flux, especially during the dry season when surface water availability required for sustaining latent heat flux would be minimal. Little is known about the seasonal and interannual variations in energy balance for campo sujo cerrado, and testing these hypotheses is important for understanding how land degradation will affect energy balance.

2. Materials and Methods

2.1. Site Description

The experimental site was located in Santo Antonio de Leverger, MT, Brazil, which is 15 km south of Cuiabá (15°43′S : 56°04′W). The study site is within a grass-dominated cerrado that was degraded approximately 35 years ago after the partial clearing of cerrado woodland vegetation. According to Koppen, the climate of region is characterized as Aw, tropical semihumid, with dry winters and wet summers. Mean annual rainfall and temperature are 1420 mm and 26.5°C, respectively, and rainfall is seasonal with a dry season extending from May to September [11]. The research area is on flat terrain at an elevation of 157 m above sea level. The regional soil type is a rocky, dystrophic red-yellow latosol locally known as a Solo Concrecionário Distrófico [23].

2.2. Data Collection

A micrometeorological tower enabled the collection of data on air temperature , relative humidity , wind speed , precipitation , soil temperature , soil heat flux , net radiation , global solar radiation , and soil moisture . and RH were measured 5 m and 18 m above the ground level using thermohygrometers (HMP45AC, Vaisala Inc., Woburn, MA, USA). was measured using two heat flux plates (HFP01-L20, Hukseflux Thermal Sensors BV, Delft, The Netherlands) installed 1.0 cm below the soil surface, with one placed in a sandy soil type and the other placed in a laterite soil type, which were typical of the local soil. and were measured 5 m aboveground using a net radiometer (NR-LITE-L25, Kipp & Zonen, Delft, The Netherlands) and a pyranometer (LI200X, LI-COR Biosciences, Inc., Lincoln, NE, USA), respectively. Precipitation was measured using a tipping-bucket rainfall gauge (TR-525M; Texas Electronics, Inc., Dallas, TX, USA). The sensors are connected to a datalogger (CR1000, Campbell Scientific, Inc., Logan, UT, USA) that scanned each sensor every 30 seconds and stored average and, in the case of P, total quantities every 30 minutes.

2.3. Data Processing

Data were collected between May 2009 and April 2011. Fluxes of latent and sensible heat were calculated over 30-minute intervals using Bowen ratio and energy balance (BREB) techniques [24] following the guidelines and modifications described by Perez et al. [25]. Bowen ratio methods have been used for decades, and while other methods, such as eddy covariance, may be more direct and amenable to analysis of measurement error, there are objective methods that are available for minimizing errors associated with resolving small gradients in vapor pressure or temperature caused by poor instrument performance and/or atmospheric conditions [25, 26].

The balance of energy was calculated as where (J m−2 s−1) was measured by net radiometer and (J m−2 s−1) was the mean heat flux in the soil measured by the soil heat flux plates installed in the sandy and laterite soils. and were calculated as a function of the Bowen ratio which, in turn, can be calculated as a function of the air temperature and vapor pressure ) gradients and the psychrometric constant [24] For (3), actual vapor pressure was calculated as a function of saturation vapor pressure and RH using (4) and (5), respectively: while psychrometric constant was calculated as a function of the specific heat at constant pressure  J kg−1 C−1 according to [27]), the local atmospheric pressure ( at the research site), and the latent heat of vaporization , which varies as a function of temperature [28], With estimates of (3), can be calculated as and can be calculated as the difference between , , and using (1) [25].

The criteria for accepting data collected from the Bowen ratio method were based on those described by Perez et al. [25]. Briefly, the Bowen ratio method fails when sensor resolution is inadequate to resolve gradients in and , stable atmospheric conditions, such as during the dawn and dusk, cause , and conditions change abruptly leading to errors in measurement [25, 26]. Using this filtering method, physically realistic values of β can be obtained in an objective, quantitative manner which limits the potential for bias and error in estimating energy balance terms [25, 29]. Gaps in estimates of and were filled by using linear relationships between retained values of and/or and measured values of .

The percentage of available energy partitioned into and was determined using linear regression, where diel (24 h) average or (dependent variables) was regressed against over monthly intervals. The slope of these regressions indicates the relative partitioning of into or . Seasonal and annual differences between energy balance terms were statistically analyzed using bootstrap randomization techniques, where the mean and the 95% confidence interval were calculated by randomly resampling each energy flux variable time series over 1000 iterations [30].

3. Results and Discussion

3.1. System Performance

Approximately 63% of all possible β values were retained after filtering for inadequate resolution, stable atmospheric conditions, and abrupt changes in measurement conditions [25]. Independent measurements of were obtained from eddy covariance in April 2011 to assess the performance in the Bowen ratio energy balance estimates. Using linear regression with the Bowen ratio estimates of as the dependent variable, the mean (±95% confidence interval) intercept and slope were −13.54 ± 4.89 W/m2 and , respectively (; observations). These data indicate that estimates of derived from two independent measurement systems were comparable and provide confidence in the Bowen ratio time series reported here.

3.2. Seasonal and Interannual Variations in Micrometeorology

Seasonal and interannual variations in micrometeorology were large over the study period (Figure 1). For example, precipitation varied from 0 mm in June 2010 to 380 mm in March 2011, in general, the months of June–August were the driest, and the months of January–March were the wettest (Figure 1(a)). Both years had a pronounced dry season (defined as the number of months when precipitation was <100 mm/month); however, the dry season in 2010 (April–December, 9 months) was substantially longer than the dry season in 2009 (April–August, 5 months) (Figure 1(a)). The dry season in 2009 was about 1 month shorter than the long-term (30 year) average, and 2010 dry season was about 4 months longer than the long-term average [11, 31]. Normally, the rainy season commences in October–November in the Cuiaba Basin [31], indicating an early transition to the rainy season in 2009 and a late transition in 2010 (Figure 1(a)). Total rainfall was 1415 m in 2009-2010 (May 1–31 April) and 1353 mm in 2010-2011, even with large interannual differences in dry season length. These rainfall totals are similar to the long-term average of 1420 mm for the region [11, 31].

Temporal trends in relative humidity were positively correlated (; ) with trends in precipitation, with the highest values (70–80%) during the peak of the wet season and the lowest values (ca. 60% in 2009 and 45% in 2010) in August-September of the dry season (Figure 1(a)). Air temperature trends were positively correlated with RH (; ); however, cross-correlation analysis indicated that RH lagged behind air temperature by approximately two months (Figure 1(a)). Air temperature typically reached a minimum in June and increased consistently at the end of the dry season (August-September) to a peak in the wet season. The increase in temperature toward the end of the dry season, coupled with corresponding increases in convection and humidity, serves as an important trigger in the transition to the wet season [31].

Seasonal patterns in solar radiation were negatively correlated with rainfall (; ; Figure 1(b)), due to frequent cloud cover during the wet season [31]. However, was negatively correlated with   , presumably because low leaf area index (LAI) during the dry season [32], when was at a seasonal maximum, causes a higher proportion of to be reflected [33].

3.3. Average Diel Patterns in Energy Partitioning

Average diel (24 h) patterns of and were consistent from month to month with maximum values observed during the midday (1200 h) local time (Figure 2(a)). Diel peaks in ranged from a maximum of 823 J m−2 s−1 in April 2009 to a minimum of 568 J m−2 s−1 in June 2009, while peaks in ranged from a maximum of 619 J m−2 s−1 in April 2009 to a minimum of 407 J m−2 s−1 in June 2009 (Figure 2(a)). However, peaks in average diel were typically lower during the wet season, while peaks in average diel were higher in the wet season, which is consistent with patterns observed in average monthly and described above (Figure 1(b)). These seasonal dynamics were likely due to variations in cloud cover, and LAI was described above [3133].

Average diel trends in and followed average diel trends in radiation closely (Table 1); however, seasonal variations in energy fluxes were large (Figure 2(b)). For example, coefficients of determination for linear regressions between and were typically >0.94, while values for linear regressions between and were typically >0.87, expect during the peak of the dry season (August and/or September) when water limitation caused reductions in (Table 1; Figure 2(b)). While diel variations in controlled the diel variations in and (Table 1), the proportion of energy partitioned into or varied depending on rainfall. Substantially more was partitioned into during the wet season, while substantially more was dissipated by during the dry season (Table 1). Midday peaks in exceeded peaks in , at times by as much as 2-fold, during the dry season, but, during the wet season, the midday peak in exceeded the peak in by the same amount (Figure 2(b)). The relative difference in the amount of partitioned into and was found to be in part controlled by monthly rainfall (Figure 3). For example, the proportion of partitioned into increased as rainfall increased up to 100 mm/month and then leveled off at approximately 0.50 at higher monthly rainfall rates (Figure 3(a)). Similarly, the proportion of partitioned into declined as rainfall increased to approximately 100 mm/month and then stabilized at 0.30 at higher monthly levels of rainfall (Figure 3(b)). These seasonal dynamics in energy partitioning are consistent with data reported from topical pastures, grass-dominated savanna, and semiarid temperate ecosystems [3336] but are much more variable compared to cerrado woodlands and tropical forests [11, 14, 3740].

Diel variations in also followed average diel trends in radiation closely (Figure 2(b)), but, as with H and , there were large seasonal variations. exhibited wider diel variations during the dry season, ranging from −60 J m−2 s−1 at night to as much as 200 J m−2 s−1 during the day, while during the wet season ranged from approximately −30 J m−2 s−1 at night to on average 100 J m−2 s−1 during the day (Figure 2(b)). Such large seasonal fluctuations in reflect the seasonal variations in soil thermal conductivity, which are affected by variations in rainfall and soil moisture, and seasonal variation in vegetation coverage, which influences exposure of soil to [19]. In general, soil thermal conductivity increases with soil moisture [41], which should result in higher during the wet season; however, increased plant growth and cover during the wet season cause shading of the soil surface, which counteracts the increase in soil thermal conductivity.

3.4. Seasonal Patterns in Energy Balance

Temporal variations in energy fluxes were large over daily, seasonal, and annual time scales (Figure 4; Table 2). During both years, was significantly larger than during the dry season, reflecting the larger surface-to-air temperature gradient [31] and the lower water availability that are typical of the dry season in south-central Mato Grosso [14]. Mean (±95% confidence interval (CI)) values of were and  MJ m−2 d−1 during the 2009 and 2010 dry seasons, respectively, resulting in a mean (±95% CI) β of in 2009 and in 2010 (Table 2). These large β values are comparable to those observed for grass-dominated cerrado [2, 4, 42] and semiarid temperate ecosystems [34] but substantially higher than those observed in cerrado woodlands and forests [11, 14, 43]. These variations reflect a decline in the as the density of woody vegetation declines [2].

The significantly higher in 2010 presumably reflected the 4-month longer dry season experienced that year (Figure 4). However, it is interesting to note that was statistically similar in the wet and dry seasons of 2009, but, in 2010, was significantly lower during the wet season (Table 2). The lower wet season in 2010 appeared to be due to heavy rainfall that occurred during the 2010-2011 wet season. For example, after a long (9 months) dry season, approximately 1005 mm of rain was recorded for January–March 2011 (Figure 1), accounting for nearly 75% of all of the rainfall for the 2010-2011 measurement year. In reality, began to decline relative to as early as November 2010, when rainfall increased but was still below the 100 mm/month threshold for the dry season (Figure 1); however, such a high amount of rainfall during the 2011 dry season would act to increase surface water availability and hence energy partitioning to and decrease surface-air temperature gradients that drive [31].

In contrast, exhibited the largest and most consistent seasonal and interannual variations that were coincident with seasonal and interannual variations in rainfall (Table 2; Figure 4). The decline in dry season was presumably due to declines in both transpiration and evaporation during the long dry season. For example, declines in surface water availability lead directly to a decline in surface evaporation and transpiration in shallow-rooted grasses [42]. Similarly, while cerrado trees are thought to be deeply rooted [9], declines in root and stem hydraulic conductance with soil drying often lead to a concomitant decline in stomatal conductance, transpiration, and leaf area during the dry season [44, 45].

4. Conclusions

Energy fluxes were measured over a two-year period in a Brazilian savanna (campo sujo cerrado) using Bowen ratio energy balance techniques. Our data indicate that rainfall was the primary control on the partitioning of available energy into sensible , latent , and ground heat fluxes. The amount of partitioned into declined as monthly rainfall increased and reached a level of approximately 30% during the wet season, while the amount of partitioned into increased as monthly rainfall increased and reached a level of approximately 60% during the wet season. As a result, was significantly higher than during the dry season, resulting in a Bowen ratio , while was higher than during the wet season, resulting in a . Our data are comparable to data collected from other grass-dominated cerrado ecosystems, but seasonal variations in energy fluxes are much higher in our system compared to tree-dominated cerrado and tropical forest because of the importance of water limitation to grasses and evaporation. Given that land cover and climate changes are expected to lead to an increase in the dry season duration and a decrease in rainfall, the high sensitivity energy partitioning to water availability in grass-dominated cerrado has important implications to local and regional energy balance.

Acknowledgments

The research was supported by Universidade Federal de Mato Grosso (UFMT) Programa de Pós Graduação em Física Ambiental (PPGFA), UFMT-Grupo de Ecofisiologia vegetal (GPEV), and Coordenação de Aperfeiçoamento de Pessoal do Ensino Superior (CAPES). Special thanks are due to Dr. Clovis Miranda and his family for allowing this work to be conducted at Fazenda Miranda.