Temporal Patterns of Energy Balance for a Brazilian Tropical Savanna under Contrasting Seasonal Conditions

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 ofR n partitioned intoH declined asmonthly rainfall increased and reached a level of approximately 30% during the wet season, while the amount of R n partitioned into L e increased as monthly rainfall increased and reached a level of approximately 60% during the wet season. As a result, H was significantly higher than L e


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 [3][4][5].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 [14][15][16], 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.0kPa 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].

International Journal of Atmospheric Sciences
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.

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].

Data Collection.
A micrometeorological tower enabled the collection of data on air temperature (  ), relative humidity (RH), 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.

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]  =  Δ Δ .
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 (1) sensor resolution is inadequate to resolve gradients in  and   , (2) stable atmospheric conditions, such as during the dawn and dusk, cause  ≈ −1, and (3) 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].

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/m 2 and 0.96 ± 0.05, respectively ( 2 = 0.80;  = 1326 observations).These data indicate that estimates of  e derived from two independent measurement systems were comparable and provide confidence in the Bowen ratio time series reported here.

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 longterm (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 ( = 0.67;  < 0.05) 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 ( = 0.65;  < 0.05); 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 ( = −47;  < 0.05; Figure 1(b)), due to frequent cloud cover during the wet season [31].However,   was negatively correlated with   ( = −0.67),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].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 [31][32][33].

Average Diel
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 ( 2 ) for linear regressions between   and  were typically >0.94, while  2 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 International Journal of Atmospheric Sciences Jun  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 [33][34][35][36] but are much more variable compared to cerrado woodlands and tropical forests [11,14,[37][38][39][40].
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.

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 southcentral Mato Grosso [14].Mean (±95% confidence interval (CI)) values of  were 4.81 ± 0.37 and 6.00 ± 0.28 MJ m −2 d −1 during the 2009 and 2010 dry seasons, respectively, resulting in a mean (±95% CI)  of 3.04±0.81 in 2009 and 5.17±1.00 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], International Journal of Atmospheric Sciences   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].

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 () > 1, while   was higher than  during the wet season, resulting in a  ≈ 1.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 grassdominated cerrado has important implications to local and regional energy balance.
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)).

Figure 1 :
Figure 1: (a) Average monthly air temperature (  ; closed symbols), relative humidity (RH; open symbols), and total monthly precipitation (P; closed bars); (b) average monthly solar radiation (  ; closed symbols) and net radiation (  ; open symbols) for the campo sujo cerrado at Fazenda Miranda during the study period.The shaded areas define the climatological dry season defined as the consecutive months when precipitation < 100 mm/month.

Figure 2 :
Figure 2: Mean monthly diel (24 hour) variation in (a) solar radiation (  ; closed symbols) and net radiation (  ; open symbols); (b) latent (  ; open symbols), sensible (; closed symbols), and ground (; crosses) heat fluxes for the campo sujo cerrado at Fazenda Miranda during the study period.The shaded areas define the climatological dry season defined as the consecutive months when precipitation< 100 mm/month.

Figure 3 :
Figure3: The proportion of net radiation (  ) portioned into (a) latent heat flux (  ) and (b) sensible heat flux (), calculated as the slope of the slope of the linear regression between   (independent variable) and  or   (dependent variables), as a function of total monthly rainfall.Linear regression statistics are from Table1.

Figure 4 :
Figure 4: Average monthly values of latent (  ; open symbols, dashed lines), sensible (; closed symbols, solid lines), and ground (; inverted triangles; dotted lines) heat fluxes, and the Bowen ratio ( ; cross symbols, long-dashed lines) for the campo sujo cerrado at Fazenda Miranda during the study period.The shaded areas define the climatological dry season defined as the consecutive months when precipitation < 100 mm/month.
(independent variable) and  or   (dependent variables), as a function of total monthly rainfall.Linear regression statistics are from Table1.

Table 1 :
Linear regression statistics for mean diel latent (  ) and sensible () fluxes (dependent variables) versus net radiation for each month during the 2009-2011 study period at Fazenda Miranda.

Table 2 :
Mean (+95% confidence interval) net radiation (  ), sensible heat flux (), latent heat flux (  ), ground heat flux (), and the Bowen ratio for the dry and wet seasons and annual cycle in 2009-2010 and 2010-2011.The number of days for each season and year is shown in parentheses.Confidence intervals were calculated using bootstrap randomization techniques ( = 1000 iterations).