Farmland expansion in the Brazilian Cerrado, considered one of the largest
agricultural frontiers in the world, has the potential to alter water fluxes
on different spatial scales. Despite some large-scale studies being
developed, there are still few investigations in experimental sites in this
region. Here, we investigate the water balance components in experimental
plots and the groundwater table fluctuation in different land covers: wooded
Cerrado, sugarcane, pasture and bare soil. Furthermore, we identify possible
water balance trade-offs due to the different land covers. This study was
developed between 2012 and 2016 in the central region of the state of São
Paulo in southern Brazil. Hydrometeorological variables, groundwater table,
surface runoff and other water balance components were monitored inside
experimental plots containing different land covers; the datasets were
analyzed using statistical parameters; and the water balance components
uncertainties were computed. Replacing wooded Cerrado by pastureland and
sugarcane shifts the overland flow (up to 42 mm yr
Brazil has significant areas used for extensive grazing over pasture and farmland (mainly soybeans and sugarcane), and parts of these areas were mostly occupied by the native Cerrado, which decreased significantly in the last century (Marris, 2005). It is estimated that 52.2 % of the original Cerrado area is now occupied by pasturelands and croplands, 0.8 % is occupied by other land uses, and 47 % remains undisturbed (Beuchle et al., 2015). The conversion from native Cerrado to pastureland, and afterwards, to sugarcane, can be considered a potential land cover and land use change (LCLUC) in southeastern Brazil (Alkimim et al., 2015). However, few studies investigate the effects of the LCLUC dynamics (Bonan, 2008; Loarie et al., 2011; Grecchi et al., 2014).
In the state of São Paulo, the minimum area required by the Brazilian forest code to maintain native forests in the Cerrado biome (20 % of the total area of a state) is not reached in much of the state (Soares-Filho et al., 2014). These areas were primary resources for expansion of the agricultural frontier in the tropics (Gibbs et al., 2010; Lapola et al., 2013), and these changes may cause significant disturbances in the hydrological processes (Loarie et al., 2011; Oliveira et al., 2014, 2015; Nobrega et al., 2017). These processes in tropical and subtropical zones are different from other regions across the world due to the increased solar energy availability and rainfall, calling for the need of basic field studies, long-term data acquisition and availability, and the development and application of mathematical models (Wohl et al., 2012; Burt and McDonnell, 2015).
In the context of the Cerrado biome, the conversion of undisturbed vegetation and pasturelands to mechanized crop systems (e.g., sugarcane, corn and soybeans) indicates that this region in Brazil has a dynamic LCLUC situation (Lapola et al., 2013). The sugarcane is the Brazilian backbone for energy security, as the ethanol production is the third-most cultivated crop after soybeans and corn, reflecting the increasing demand for automotive fuels over the years (Leal et al., 2013; Rodrigues et al., 2018). Thus, the country is the world's second largest ethanol producer, and the Cerrado comprises the sugarcane expansion frontier due to the availability of water and pasturelands for the crop expansion (Bellezoni et al., 2018).
Native forests help in maintaining the water cycle (Krishnaswamy et al., 2013; Ghimire et al., 2014). The evapotranspiration appears to be a key component in the aquifer recharge control in the Guarani aquifer system (GAS) outcrop zone (Wendland et al., 2007; Lucas and Wendland, 2015; Lucas et al., 2015). Thus, water partitioning in these areas should be well understood in different land cover and land use (LCLU) types, allowing the evaluation of possible water balance trade-offs due to the LCLUC (Krishnaswamy et al., 2013; Ghimire et al., 2014; Frank et al., 2014) caused by environmental and economic needs (Barretto et al., 2013).
In Brazil, many large-scale studies on water balance were developed in some hydrographic regions in the country (Paiva et al., 2013; da Paz et al., 2014; Oliveira et al., 2014; Melo et al., 2016). Nevertheless, experiment-scale studies are still rare due to the local heterogeneities and uncertainties from hydrological measurements and estimates (Beven, 2006; Graham et al., 2010). Basic field-hydrology studies are important for improving the agricultural production efficiency while promoting sustainable development. Therefore, these studies are important for promoting new solutions and techniques to maintain the water balance in spite of the rapid LCLUC (Dotterweich, 2013; Nobrega et al., 2017). Research of this kind can be done using experimental plots, which are delimitated hillslopes (control volume) where the overland flow is directed to one outlet (Sadeghi et al., 2013; Oliveira et al., 2015, 2016; Mwango et al., 2016; Strohmeier et al., 2016; Youlton et al., 2016b; Anache et al., 2017, 2018; Zhao et al., 2017).
The objective of this study is to evaluate the water fluxes in different LCLU types: wooded Cerrado (also known as Cerrado sensu stricto), sugarcane, pasture and bare soil. In addition, we identify possible water balance trade-offs due to the different LCLU types, compute the uncertainties for each water balance component and propagate the error separately for each LCLU.
Location of study sites, Cerrado biome borders, Guarani aquifer system (GAS) outcrop zone distribution in Brazil, and experimental design, where site 1 contains the plots with agricultural land uses (pasture, sugarcane and bare soil) and site 2 contains the plots with wooded Cerrado.
The study sites are located in the Arruda Botelho Institute (IAB) in
Itirapina, São Paulo, Brazil (latitude 22
The experiment began in October 2011 (Oliveira et al., 2015, 2016; Youlton et al., 2016b, a), and it is still in operation. This study considered the data collected during 5 years (2012–2016). It contains manual and automatic instruments, as well as permanent structures, all assembled in the two sites (Fig. 2). Two monitoring wells, meteorological stations (tripod and tower), and 12 steel-made bounded plots are the permanent structure. The instruments are described according to their positioning and function (Table 1). All instruments recorded data every 10 min, except the pressure transducers, which logged the groundwater table twice a day. The distance between sites 1 and 2 is 1.7 km. Site 1 has nine plots placed side by side (approximately 2.5 m of distance between each plot). Also in site 1, there is a meteorological station that concentrates almost all sensors placed at that site at one point, except for the soil moisture probes. They are connected to the meteorological station, but they are placed inside the plots (there is one placed inside the first sugarcane plot, and 20 m to the left, there is another one placed inside the first pasture plot). Site 2 has three plots inside a tropical woodland (wooded Cerrado), and due to the tree density and topography, the plots are approximately 5 to 10 m distant from each other. Approximately 50 m to the north of the plots, there is a meteorological tower (11 m height) containing all the sensors placed at site 2.
Hydrological monitoring performed on the four treatments: (1) relative humidity and air temperature probes at 2 m (site 1) and 11 m (site 2), (2) rainfall gauges, (3) soil moisture sensors, (4) solar radiation sensor, (5) wind speed and direction (anemometers) at 2 m (site 1) and at 11 m (site 2), (6) surface runoff collectors, (7) monitoring wells equipped with water table pressure transducers, (8) net radiation sensor, and (9) soil heat flux plate.
Monitoring site instrumentation characteristics.
The plots were designed to adequately represent the process heterogeneities
(Sadeghi et al., 2013): three replicates for each LCLU, with 5 m
width, 20 m length and a 9 % slope gradient. The plots located at site 1
contain three different LCLU types, including (i) pasture (
The water balance components (Eq. 1) were monitored over 5 years (2012–2016)
in a control volume defined by the bounded plots placed within the
experimental sites. We considered different techniques for monitoring these
components according to the LCLU conditions. The water balance residual
(d
The rainfall was monitored using a tipping bucket (Model TB4, Hydrological Services) with a gauging resolution of 0.254 mm (Table 1). The rainfall data were registered every 10 min in order to obtain rainfall intensity and duration.
Rectangular experimental plots directed the overland flow to tanks at the end of the slope, where the volume was measured at times after rainfall events. We also calculated the runoff coefficient for each LCLU using a genetic algorithm to minimize the squared errors between observed and estimated runoff values using the rational method (Wang, 1991).
The reference evapotranspiration (
The water stress coefficient (Eq. 4) was calculated for each day
All the necessary data and coefficients for calculating the water stress
coefficient (
Variables used to calculate evapotranspiration for sugarcane and pasture.
Approximate sugarcane age (days).
The wooded Cerrado evapotranspiration (
The bare soil condition has no vegetation, and, consequently, there is no transpiration. Thus, we applied the method developed by Ritchie (1972). This method has two phases: firstly, the soil evaporation is equal to the potential soil evaporation estimated using the Priestley and Taylor (1972) method adapted to free surfaces. During this phase, there is no water restriction (precipitation higher than evaporation) and the evaporation is governed by the available energy; secondly, the accumulated soil evaporation exceeds the precipitation, and the soil evaporation is currently given as a function of the dry days that followed the last wet day. The evaporation cycle is interrupted and returns to the first phase whenever the precipitation exceeds the accumulated evaporation during the second phase.
The water table was registered twice a day (at 06:00 and 18:00 BRT – Brasilia Time) using pressure transducers (Diver, Schlumberger) placed inside two monitoring wells (well 1 is located in the pasture area, and well 2 is located inside the wooded Cerrado area). In the study site, both wells presented similar hydraulic conductivity according to the slug test (Bouwer and Rice, 1976) previously performed. We evaluated the aquifer hydraulic conductivity from both wells in order to validate the water table comparison among each other, as if the aquifer condition in the wells were different, such comparison would not be fair. Both wells reach the water table at approximately 40 m depth in an unconfined sandstone formation (Botucatu Formation), which belongs to São Bento Group of the Mesozoic age. Furthermore, the soils above the aquifer that appears throughout the unsaturated zone are Cenozoic sediments weathered from the sandstone (Wendland et al., 2007). Additionally, despite the limited number of wells, the experimental design allowed a first look into the groundwater table behavior under different LCLU types (pasture and wooded Cerrado) and a cross-check with the surface water balance outcomes.
The normality assumption was tested using the Shapiro–Wilk test considering a 95 % confidence interval for rainfall, evapotranspiration, surface runoff and water balance residual datasets. The one-way analysis of variance (ANOVA) was applied to test the null and alternative hypothesis, that is, equality of surface runoff, evapotranspiration and water balance residual distribution functions between the four treatments (LCLU) versus the difference in distribution functions between at least two treatments. Additionally, the multiple comparisons between treatments were performed using the Tukey test (Montgomery, 2008). The rainfall, evapotranspiration, surface runoff, water balance residual and soil moisture graphs were plotted using a daily basis timescale. The groundwater table fluctuation was plotted using a monthly timescale due to the noise typically found in this kind of measurement. In order to present the order of magnitude over the years, the data were also resumed annually in tables and figures.
Data uncertainties are flaws found in the available information used to
represent the reality and basically depend on the knowledge of the observed
data (Refsgaard et al., 2007). The water balance uncertainty
can be given by the standard error propagated from its components. Thus, the
water balance residual (d
The absolute error of the rainfall estimations was calculated using Eq. (9)
based on the instrument accuracy informed by the manufacturer. Considering
the used tipping bucket rain gauge (TB4), the error may reach up to
The surface runoff may vary due to the heterogeneities found between plots'
replicates (Wendt et al., 1986; Nearing et al., 1999; Gómez et al.,
2001; Sadeghi et al., 2013). Thus, the standard error for surface runoff is
given by the standard deviation of the replicated plots:
Water balance components for different LCLU types, namely rainfall,
As previously mentioned, the evapotranspiration was estimated using the FAO
56 methodology (Allen et al., 1998) for pasture and sugarcane
and the Priestley and Taylor (1972) method for wooded Cerrado and
bare soil. However, this study does not have evapotranspiration observations
to evaluate how well this variable was estimated in the study sites.
Consequently, the uncertainties for evapotranspiration estimates were
calculated by combining all input variable uncertainties (Eq. 11) which were
measured in the field for the FAO 56 method (temperature, relative humidity,
solar radiation, barometric pressure and soil moisture) and the Priestley
and Taylor method (temperature, net radiation and soil heat flux). The
equation is as follows:
The water balance components show different patterns according to the LCLU
(Table 3 and Fig. 3). We verified this using a multiple comparison test
(Tukey) in which the water balance residual variation (d
Mean and standard deviation of the annual water balance components for 2012–2016 period.
d
The average annual rainfall in the study site was 1388 mm yr
The evapotranspiration estimates (Fig. 3b) were different between the considered land covers (Table 3). We observed that the wooded Cerrado evapotranspiration presented the highest rates among all analyzed LCLU types and also the smallest variability throughout the year. However, there is no agreement between the measurements and estimates performed in wooded Cerrado areas and the present study (Table 4), due to the diverse rainfall patterns among the study sites and the different methods used to measure or estimate the evapotranspiration. Thus, this study shows evidence of the need for reference values of evapotranspiration in undisturbed areas for a better understanding of their role in the water cycle. In addition, the rainfall is different for the studies compared in Table 4.
Evapotranspiration and rainfall from studies developed in Brazil under similar LCLU to the present study.
Lat: latitude. Long: longitude.
The sugarcane evapotranspiration rates were higher than pasture due to the higher crop water demand during the initial phase of its annual cycle, and both evapotranspiration estimates agree with measurements performed by previous studies (Sakai et al., 2004; Cabral et al., 2012; Nobrega et al., 2017). Both sugarcane and pastureland ceased the evapotranspiration during the dry season due to the decoupling condition, when the radiation is the only contributor to the evapotranspiration process (water stress condition; Pereira, 2004). This happens due to the lack of water readily available for the plants along the root zone. This condition does not repeat in the wooded Cerrado, as its root system is deeper (Canadell et al., 1996), and the plants may reach water in depths where sugarcane and pasture cannot do this. However, the root zone depth of undisturbed vegetation such as the wooded Cerrado is uncertain and may vary according to the soil characteristics (Canadell et al., 1996) and the groundwater table level (Leite et al., 2018). It may influence the plants' transpiration (Rawitscher, 1948; Oliveira et al., 2005) and, consequently, the water balance residual.
The surface runoff (Fig. 3c) presented the lowest values among all water balance components, and it was the most influenced variable by the land use. There was no surface runoff generation during August due to the non-occurrence of rainfall events. The bare soil was the most sensitive to the runoff generation throughout the dry season, when the occurrence of rainfall events was lower. The months of January, February and March registered the highest averages for rainfall due to the wet season. The runoff coefficients for wooded Cerrado, sugarcane, pasture and bare soil were 0.001, 0.007, 0.029 and 0.063, respectively. These values are in compliance with a previous study performed by Oliveira et al. (2016) for similar conditions. Although a higher runoff coefficient was expected for sugarcane in comparison with the pasture, sugarcane presented runoff values significantly lower than pasture as the soil tillage increased water infiltration and, consequently, reduced the surface runoff. Additionally, there is an apparent inverse linear relationship between culture size and runoff coefficient.
The surface runoff in sugarcane plantations is not well understood by the scientific community despite the economic importance to the country. The monitored values agreed with previous studies in the same study area (Oliveira et al., 2016; Youlton et al., 2016b; Anache et al., 2018). The highest runoff rates were registered after planting and harvesting events, when the soil is more exposed to the rainfall.
In addition, the surface runoff in the pasture presented significantly higher rates in comparison with sugarcane and wooded Cerrado, due to the top soil compaction caused by grazing. The order of magnitude of the measured runoff in the pasture plots is similar to that previously published (Saraiva et al., 1981; Silva et al., 2011; Dedecek, 1989). The high runoff variability observed in this land use was due to the variable presence of animals which were managed in an extensive approach, and, consequently, heterogeneities may happen in soil compaction and vegetation conditions (Nacinovic et al., 2014). Additionally, the deforestation and agricultural land uses may increase soil compaction, as the LCLUC influences the hydrological patterns along the soil profile by evident modifications in the soil characteristics (bulk density, infiltration capacity, etc.; Lamparter et al., 2016; Meister et al., 2017; de Almeida et al., 2018).
The wooded Cerrado presented the lowest runoff rates among all LCLU types due to the higher soil protection, which avoided the overland flow generation. The order of magnitude of the surface runoff in the wooded Cerrado was similar to those found in shrublands and forested areas (Dedecek, 1989; Silva et al., 2011; Oliveira, 2012, 2015; Nacinovic et al., 2014). These reduced surface runoff rates in the wooded Cerrado increase soil water infiltration in comparison with pasture, sugarcane and bare soil. Thus, higher infiltration rates increase plant water availability (Krishnaswamy et al., 2013). Additionally, the organic matter layer above the soil (wooded Cerrado's forest floor litter) reduces the overland flow. In some cases, the litter removal may increase the surface runoff up to 50 % (Gomyo and Kuraji, 2016).
Water balance uncertainties.
d
The water balance residuals (Fig. 3d) showed the water surplus (positive values) and deficits (negative values), evidencing the consumption of the soil water storage in the dry season (June–September). Such behavior is also evidenced by the average soil water content along the first meter of soil (Fig. 3e), as the soil moisture became lower in the wooded Cerrado in comparison with pasture and sugarcane observations during the dry season. This is because only the wooded Cerrado condition had negative values during the dry season, as the other land uses (sugarcane and pasture) had no conditions to remove water from deeper regions due to the shallow root system and their physiological characteristics. Thus, wooded Cerrado vegetation adapted to dry weather conditions due to the plant water demand for evapotranspiration, which reduced gradually as long as the water balance residual accumulated during the wet season was consumed (Oishi et al., 2010; Christoffersen et al., 2014; Cabral et al., 2015; Oliveira et al., 2015). In addition, the structural quality of soil from undisturbed woodlands leads to a higher capacity for retaining water than that of pasture soil (Tseng et al., 2018).
The groundwater table fluctuation inside the wooded Cerrado area was lower
compared to the pasture area considering the observations of the monitoring
wells (Fig. 3f). This variation in the wooded Cerrado was less than
1 m yr
There is clear evidence that a time lapse between the water infiltration through the soil and the aquifer recharge exists due to the huge non-saturated zone thickness (around 40 m; Fig. 3f). For this reason, we cannot ignore that the evapotranspiration influences the aquifer recharge, as previously reported by other studies under conditions similar to and different from the present study (Scanlon et al., 2005; Scott et al., 2014; Lucas et al., 2015; Oliveira et al., 2017; Lucas and Wendland, 2015). Changes in the LCLU, such as the potential conversion from wooded Cerrado to an agricultural LCLU (here we tested the pasture), may affect groundwater recharge, processes and availability. These effects tend to be non-linear and difficult to analyze as they result from complex interactions between LCLU and hydrological processes (Han et al., 2017). Thus, mathematical approaches (Archer and Fowler, 2018; Gómez et al., 2018) or natural tracers (Su et al., 2018) are useful tools for verifying the response time of the groundwater table to the water balance from different LCLU types. In further studies, such techniques may be part of a solution for investigating how responsive the aquifer is to the surface water partitioning in the conditions considered here over time.
The evapotranspiration and root zone depth controlled the water balance
residual and, consequently, the water percolation throughout the
non-saturated zone and aquifer recharge (Finch, 1998; Gouvêa and
Wendland, 2011; Krishnaswamy et al., 2013; Lucas and Wendland, 2015;
Domínguez et al., 2016; Manzione et al., 2017). In the pasture area
(site 1), the soil water that was not consumed by the plants did not
evaporate, flowing down along the unsaturated zone. Consequently, the water
uptake by the plants became unfeasible, as the roots were more shallow (see
Table 2) than in the wooded Cerrado. Furthermore, the aquifer recharge
decreases as the vegetation density increases in undisturbed Cerrado areas
(Oliveira et al., 2017), following the principle that the increased
canopy cover of the wooded Cerrado found in the study area may occur due to
the deep groundwater level (Leite et al., 2018; Villalobos-Vega et al.,
2014). However, the water balance analysis performed here focused on
hillslope hydrology, and the monitoring well depths reflect the aquifer
behavior in a broader area covered by pasture and wooded Cerrado in
comparison with the 100 m
Annual water balance residual (d
The water balance component (
The water balance residual (d
Annual observed means (2012–2016) of hydrological trade-offs
related to evapotranspiration (
The Brazilian Cerrado is very important economically, as it is responsible for most of the agricultural production that supplies both external and internal markets (Klink and Machado, 2005). A better understanding of the trade-offs between the ecosystem and economical needs that govern the land cover and land use dynamics in the Cerrado biome is necessary (Marris, 2005). The undisturbed Cerrado area was reduced to 50 % of its original extension due to intense land use (Lapola et al., 2013; Alkimim et al., 2015). The undisturbed vegetation helps in maintaining the water cycle, with low surface runoff and high evapotranspiration rates.
This study shows evidence that the suppression of the wooded Cerrado and the
conversion to agricultural land uses, such as sugarcane and pasture,
increased the surface runoff and decreased evapotranspiration, even
considering measurements and estimation method uncertainties. Consequently,
the water balance residual (d
Previous studies show that native forests help in maintaining aquifer recharge by the high infiltration rates, similar to observations performed in tropical forests in southern India and in mountainous areas in the Himalayas (Krishnaswamy et al., 2013; Ghimire et al., 2014). Nevertheless, evapotranspiration appears to be a key component in the aquifer recharge control in the study site condition, which is located in a Guarani aquifer system outcrop zone (Lucas et al., 2015; Lucas and Wendland, 2015). Therefore, the aquifer recharge rates, evidenced here by the groundwater table fluctuation (Fig. 3f), may be reduced in forested areas in comparison with agricultural landscapes due to the atmospheric and vegetation water demands and the increased soil water retention capacity (Adane et al., 2018; Dias et al., 2015; Wang et al., 2018; Tseng et al., 2018). This validates the information stating that the LCLU significantly impacts groundwater recharge (Scanlon et al., 2005; Scott et al., 2014; Lucas et al., 2015; Dawes et al., 2012). This fact may suggest that the ecosystem service of native forests, such as the wooded Cerrado, is not the aquifer recharge maintenance but rather the constant return of water to the atmosphere throughout the year.
Significant and non-significant changes in the water partitioning may be
observed in the case of substituting the wooded Cerrado with bare soil,
pasture or sugarcane (Fig. 5). We highlight the
Sugarcane and pasture presented trade-offs that were equal, as non-significant changes occurred among each other's water balance components. However, the sugarcane plantation presented a higher potential for maintaining the evapotranspiration rates that are closer to the undisturbed Cerrado conditions than the pasture, agreeing with previous estimates (Loarie et al., 2011). In addition, it is suggested that sugarcane has increased evapotranspiration rates in comparison with other annual crops (Guarenghi and Walter, 2016; Hernandes et al., 2018b, a).
This paper presented an experimental approach at the hillslope scale,
concerning the possible water partitioning trade-offs due to the LCLUC
dynamics. We monitored the water balance components over 5 years in
different land covers: wooded Cerrado, pasture, sugarcane and bare soil.
These land covers are subjected to the current LCLUC dynamics in
southeastern Brazil. The water partitioning observations in different LCLU
types confirm that modifications in the land surface conditions may significantly
change the water balance residual (up to 584 mm yr
Higher water consumption by dense and native vegetation, such as the wooded Cerrado, also happens due to the higher infiltration rates, increasing the plant water availability. In general, the root systems are deeper than pasturelands and sugarcane plantations and have the capacity of reaching water for transpiration deeper in the soil profile. This maintains evapotranspiration throughout the year, even during the dry season. Hence, less water becomes available for the aquifer recharge in areas where the canopy layer is higher and denser. However, reference values for evapotranspiration in undisturbed land covers, such as the wooded Cerrado, are still needed in order to reduce uncertainties from the current approximations and validate the water balance hypothesis. Particularly the undisturbed Cerrado and the agricultural LCLU should be investigated concerning their potential for maintaining the aquifer recharge and groundwater availability as well as answering how responsive the aquifer is to the surface water partitioning in different LCLU types over time.
The datasets underlying this research are available in the Supplement (S1 and
S2) of this paper, and they are accessible from the following data repository
link:
The supplement related to this article is available online at:
JAAA, EW, PTSO and CY designed the experiments, and JAAA, LMPR, PTSO and CY carried them out. JAAA, LMPR and EW analyzed the experiments' outcomes and discussed the results. JAAA and EW prepared the paper with contributions from all co-authors.
The authors declare that they have no conflict of interest.
This study was supported by grants from the Ministry of Science, Technology, Innovation and Communication – MCTIC – and the National Council for Scientific and Technological Development – CNPq (grant numbers 201109/2015-8, 142393/2015-0, 150057/2018-0, 441289/2017-7, 306830/2017-5, and 302991/2017-4); the São Paulo Research Foundation – FAPESP (grant number 2015/03806-1); and the Coordination of Improvement of Higher Education Personnel – CAPES (finance code 001). The authors acknowledge the graduate program in Hydraulics and Sanitary Engineering – PPGSHS (USP-EESC) – for the scientific support and the Arruda Botelho Institute – IAB – for allowing the development of this study on its private land. The authors would like also to thank the editor and the anonymous referees for their useful comments, which substantially improved the paper. Edited by: Martijn Westhoff Reviewed by: three anonymous referees