Next Article in Journal
Mineral and Citrate Concentrations in Milk Are Affected by Seasons, Stage of Lactation and Management Practices
Previous Article in Journal
The Past, Present, and Future of Barley Yellow Dwarf Management
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Influence of Land Use and Land Cover on Hydraulic and Physical Soil Properties at the Cerrado Agricultural Frontier

by
Emily Ane Dionizio
* and
Marcos Heil Costa
Department of Agricultural Engineering, Federal University of Viçosa (UFV), Viçosa, MG 36570-900, Brazil
*
Author to whom correspondence should be addressed.
Submission received: 31 December 2018 / Revised: 11 January 2019 / Accepted: 15 January 2019 / Published: 21 January 2019

Abstract

:
Western Bahia is one of the most active agricultural frontiers in the world, which raises concern about its natural resources conservation, especially regarding water availability. This study evaluated the influence of five different land uses and land covers on physical and hydraulic soil properties, and developed pedotransfer functions to derive regional hydraulic properties. Significant changes between physical and hydraulic soil properties under agricultural areas and under natural vegetation cover were found, reinforcing that agricultural activity may influence the soil water balance. Cerrado and Forest formation areas have higher infiltration rates ( K s a t ) compared to managed areas, with average values of 16.29 cm h−1, and 14.47 cm h−1, while irrigated croplands, rainfed croplands and pasture areas have much smaller infiltration rates, with K s a t equal to 3.01 cm h−1, 6.22 cm h−1 and 5.01 cm h−1, respectively. Our results suggest that the agriculture practices do not directly affect the vertical nature of hydrological flowpath, except in the case of intensive irrigated agriculture areas, where K s a t reduction can lead to erosive processes favoring organic matter losses, and decreases in productivity and soil quality. Impacts of land use change on hydraulic and physical soil properties are a reality in the Cerrado agriculture frontier and there is an urgent need to monitor how these changes occur over time to develop effective mitigation strategies of soil and water conservation.

1. Introduction

In the last decade, the rapid expansion of agribusiness in the Cerrado led to a new Brazilian agricultural frontier known as MATOPIBA (acronym for the states Maranhão, Tocantins, Piauí and Bahia), which have raised concerns about the natural resources conservation [1,2], especially regarding to the water availability [3,4,5].
In MATOPIBA, Western Bahia stands out by the agricultural expansion, representing 49% of the total agricultural area [6] with 1.8 million hectares in 2015, which is equivalent to an increase of 352% since 1985 [7,8]. The irrigated croplands alone had an increase of 526% between 1985 and 2002, with a simultaneous decrease of Cerrado, seasonal forest and transition vegetation areas by 881,483 ha, 66,417 ha and 269,592 ha, respectively [9]. In fact, the region also stands out for its high productivity records, reaching to 7.4 million tons of soybean, cotton and maize crops in the 2016/2017 harvest [10], confirming the potential of the agribusiness and the systematic dependence of the water availability and natural resources. Although Western Bahia is located above the Urucuia Aquifer and is drained by the Grande, Corrente and Carinhanha basins, there is a major concern about the regional water availability and the impacts of the agricultural activities on Cerrado biome [1,2,5].
The effects of land use change on soil physical proprieties are broadly known, especially when considered conversion of tropical forest to pasture or croplands [11,12,13]. Changes in soil bulk density, penetration resistance, porosity, near-surface hydraulic conductivity [14,15,16], infiltrability and saturated hydraulic conductivity [13,17] are described as possible consequences of the land use change.
In Cerrado, these effects have been studied in the last decade [18,19,20], and received a greater focus due to the advancement of the agricultural frontier from Southern Amazonia to MATOPIBA [5,20]. Soil compaction [21,22], erosion [23], and decreased permeability and water infiltration rates [22,24] are some effects of land use change that affects the physical and hydraulic soil properties causing damage on soil aeration, and soil water dynamics in Cerrado. However, most of these studies are based on few data and do not allow spatial generalizations. Remote sensing and modeling techniques have been used to broadly study the influence of land use on the water balance [5,25,26,27], in the dynamic of vegetation [1,7] and climate [28,29,30,31], and show that the replacement of native vegetation by croplands alters in Cerrado the amount of water recycled to the atmosphere at a large-scale [5], affecting the regional climate dynamics [28,29,30,31]. The results of these simulations, however, are very sensitive to the soil physical and hydraulic parameterizations in the model.
Thus, the knowledge of the hydraulic and physical soil properties contribute to development of powerful tools as hydrological and dynamic vegetation models, and remote sensing techniques to estimate the water recharge, the soil water availability, the cropland productivity, or the influence of extreme precipitation scenarios in the Cerrado agricultural frontier. Moreover, the hydraulic and physical soil properties are fundamental to make inferences about the soil quality and sustainability, allowing the development of alternatives that can prioritize the water and soil conservancy and preserve the remaining native vegetation while ensuring the increase of agricultural production.
The main goal of this study was to evaluate how land use and land cover affect the hydraulic and physical soil properties at the Cerrado agricultural frontier. In addition, we developed pedotransfer functions to derive hydraulic properties for use in dynamic vegetation models.

2. Methods

2.1. Study Area

The study area of this work is the Western Bahia located on the geological formation of the Urucuia Group (Upper Cretaceous), which is one of the main areas of agricultural expansion in the Cerrado biome. This region is drained by three important rivers, Grande, Correntes and Carinhanha, with an area of 131,168.59 km2 (Figure 1). Marked by the arrival of migrant rural producers from the southern region in 1990, Western Bahia is today one of the largest producers of soybean in MATOPIBA and the largest producer of cotton in Brazil.
In this region, the physiognomy vegetation is predominantly Cerrado stricto sensu, with predominance of the tree–shrub stratum [32]. The soils are acidic with low fertility, but located at flat or mildly hilly areas. Thus, the use of high level technology and the chemical inputs allows the correction of fertility providing favorable conditions to expansion and intensification of agriculture. The soil granulometry is predominantly sandy and has medium texture, classified as Latossolos (57%), Neossolos (29.6%) and Cambissolos (7%), according to the Brazilian soil classification system [33].
The regional climate is tropical humid according to Köppen [32], and presents two well-defined seasons, dry (October to April) and rainy (May to September). The average annual temperature is 24 C and the average annual rainfall is 1400 mm in the extreme west, gradually decreasing to 800 mm in the east.

2.2. Soil Sampling Design

Soil samples were collected in Western Bahia in two field campaigns in 2017 for five land use and land cover classes (LULCC): (1) Cerrado formations (CDO); (2) forest formations (FOR); (3) rainfed agriculture (RAG); (4) irrigated agriculture (IRR); and (5) pasture (PAST) (Figure 1).
Each class was sampled at 20 sampling points following the criteria: opening year of the agriculture or pasture area between 1990 and 2017, and logistical access to farms, including road conditions and permission to enter the farm by the farmer. In some cases, CDO and the FOR areas were sampled along the road between farms in natural vegetation areas rather than inside the farms. In total, 700 samples were collected considering 5 LULCCs × 20 sample points × 7 depth levels. At each sampling point, undisturbed samples were collected for 0–5, 5–10, 10–15, 15–20, 30, 50 and 70 cm of depth using an Uhland soil sampler, with metal cylinders of 100 cm3 volume. The 100 samples of 0–5 cm of depth were used to calculate the saturated hydraulic conductivity ( K s a t ), soil water content at 10 kPa and 1500 kPa ( θ f c , θ w p ), texture (coarse sand fraction ( f c s a ), fine sand ( f f s a ), and silt ( f s ), total porosity ( ϕ t o t ), microporosity ( ϕ m ), macroporosity ( ϕ M ), and soil particle density ( ρ s ), while all layers (0–5, 5–10, 10–15, 15–20, 30, 50 and 70 cm) were used to calculate the soil bulk density ( ρ a ).

2.3. Measuring Methods

The K s a t of soils was determined via the constant head permeameter method [34]. The constant water flow that flowed through the soil sample was measured and applied into the equation for direct calculation of K s a t . For textural analysis and particle density, samples were dried in the open air, and sieved in a 100-mesh sieve. Soil fractions ( f c s a , f f s a , and f s ) were separated according to Ruiz [35], using the sieve method for the sand fraction (0.05–2 mm) and the pipette method to determine the silt (0.002–0.05 mm) and clay (<0.002 mm) fraction.The soil water retention curve (SWRC) was measured using a sand tension table at the matric potentials of 1 and −6 kPa, and a Richards chamber with porous plates, at matric potentials of 10, 50, 100, 500 and −1500 kPa [36]. The centrifuge method was used to measure the soil water content at matric potential of 30 kPa. The SWRC was measured for 0–5 cm layer for all LULCC, and also for 5–20 cm layer in IRR areas. ρ a and ϕ t o t analyses were carried out using methods described by CLAESSEN et al. [37]. The ϕ m was determined by water content in a volumetric ring under 0.6 m of water column tension, and the soil ϕ M was estimated by the difference between total soil porosity and soil microporosity.

2.4. Pedotransfer Model

2.4.1. Adjustment of Soil Water retention Curve

The Campbell and Norman [38] soil model was fitted to soil water content data using Equation (1).
ψ m = ψ e θ θ s b
where ψ m is the soil water matric potential in kPa, ψ e is the air entry matric potential in kPa, θ is the volumetric water content cm3 cm−3, θ s is the saturated volumetric water content cm3 cm−3, and b is the empirical Campbell parameter, related to the particle size distribution. It is strongly dependent on soil texture [39] and is considered an index for soil pore-size distribution [40]. This model was chosen due to its minimal set of parameters necessary to describe the soil hydraulic properties, favoring its implementation in regional and global scales, and it has been widely used in modeling studies [41,42].
A linear regression with log-transformed data was used in the equation log y = a + x b , to determine the ψ e and the empirical value of b Campbell parameter, where a is the intercept of the soil water retention in the log-log system and b represents the slope.
The linear least squares method was used to adjust the ψ e and b Campbell constant for each sample. Then, the raw data were also used to develop the pedotransfer functions.

2.4.2. Development of Pedotransfer Functions

The pedotransfer functions (PTFs) were developed using the soil physic proprieties and the SWRC for 100 undisturbed samples collected in the 0–5 cm soil layer. Initially, we tested the normal distribution for f c s a , f f s a , f s , ϕ t o t , ϕ M , ϕ m , ρ s and ρ a properties using Shapiro–Wilk statistic test considering α = 0.05. The properties that are normally distributed were used to develop the PTF through multiple regression, while the predictors not normally-distributed, such as ρ s , f f s a , f s , ϕ M and ϕ m , were excluded from analysis.
To develop and validate the pedotransfer function, the observed data for water retention curve data and physical properties were separated into two groups, where the first one contained 75% of the data (calibration group) and the second 25% of the data (validation group). The samples of calibration and validation groups were randomly chosen. Between these two groups the average of soil physical properties were not significantly different at α = 0.05 according to the Student t test, except for f c s a in the forest formations areas FOR. However, for hydraulic proprieties ( ψ e and b), a significant difference was found between b Campbell parameter used in Calibration and Validation Group for CDO and FOR classes.
A multiple linear regression model was adjusted for each hydraulic Campbell parameters ( log ψ e , b, K s a t , θ f c , and θ w p ) using the soil physical properties measurements as predictors (Equation (2)):
y = β 0 + β 1 ρ a + β 2 ϕ t o t + β 3 f c s a
where y is one of the five hydraulic parameters ( log ψ e , b, K s a t , θ f c , and θ w p ), ρ a is the soil bulk density in g cm−3, ϕ t o t is the total soil porosity in cm3 cm−3, and f c s a represents the coarse sand fraction in percent. The stepwise method with 5% of significance was chosen to select the most important variables for determination of y through the backward and forward mechanism. This stepwise method uses the Akaike criterion to eliminate collinear variables, excluding non-informative variables of the final model [43].

2.5. Data Analysis

We tested the homoscedastic and normality of residuals applying the Shapiro–Wilk test [44] for all soil physical properties, and used the Tukey Cramer method to compare the soil physical properties among LULCC classes. The soil physical properties that were normally distributed (i.e., f c s a , ϕ t o t , and ρ a , with p-values greater than 0.05 in the Shapiro–Wilk test) were applied to the analysis of variance and were used for the development of pedotransfer function.The t-test was used to compare the average of the 75% of soil properties used for the development and 25% used for the validation of pedotransfer function (0–5 cm layer), and for the comparisons of the differences for each soil layers among LULCC. All data analyses were carried out using R software [43].

3. Results

3.1. Soil Physical and Hydraulic Properties

Most of the soil samples were sandy clay loams (40%), followed by 23.3% of sandy loams and 23.3% of loamy sands. The 13.4% of remaining samples were distributed among sandy clay, sand and silt clay classes (Figure 2). In the Brazilian soil classification system (SiBCS), these samples presented a similar distribution along the soil classes, with predominance of 40% of the samples in MeA (Médio argilosa), and 56% divided equally between MeAr (Médio Arenosa) and ArMe (Arenosa Média) classes, highlighting the predominance of sandy soils in the region (Figure 2). The average sand content ranged between 69.47% and 85.79% within the predominant soil classes, while the average clay content ranged between 11.29% and 26.87%.
The results of average soil bulk density ( ρ a ) suggest a compacted soil surface for agriculture land use, ranging between 1.52 and 1.61 g cm−3 (Figure 3a and Table 1). Under natural land cover, FOR and CDO, the ρ a was lower and statistically different from agricultural areas with average ρ a ≤ 1.36 g cm−3 (Figure 3a and Table 1). Along the soil profile, the soil ρ a showed an increase trend from 0–5 cm to 30 cm for all LULCCs (Figure 4). This increase of ρ a in the subsurface layers is higher in areas under managed soil, than in FOR and CDO soils with values above 1.65 g cm−3 (Figure 4).
The particles density ( ρ s ) ranged between 2.5 g cm−3 and 2.65 g cm−3 for all LULCCs, showing no statistical differences according to the Tukey Cramer test. The soil total porosity ( ϕ t o t ) ranged between 26% and 60%, with average of 43% for all samples collected (Table 1 and Figure 3b). The compaction pattern found in ρ a was also observed in ϕ t o t , with reduction of the ϕ M in agriculture land use compared to CDO and FOR areas. In natural ecosystems, ϕ t o t ranged between 26% and 57% for CDO and between 38% and 60% for FOR, while in agriculture systems the total soil porosity ranged between 30% and 53% (Figure 3e). Average CDO and FOR ϕ t o t was greater than 45%, while under agriculture land use ϕ t o t were smaller than 41% (Table 1).
The soil compaction decreased ϕ M in PAST and RAG areas, reducing the soil capacity to drain excess water after a heavy precipitation—lower K s a t values. The reduction of ϕ M contributed to an increase of ϕ m , which may alter the soil aeration and roots growing conditions. The ϕ M in managed areas ranged between 9% and 12%, while ϕ m were above 25%. Although the values of ϕ M and ϕ m of the IRR class were not measured, significant differences were found for total soil porosity in relation to CDO and FOR (Figure 3e).
Soils under natural vegetation cover had a higher infiltration rates ( K s a t ) compared to managed areas, with mean values of 16.29 cm h−1 for CDO, and 14.47 cm h−1 for FOR, while among the agriculture land uses, average infiltration rates were much smaller, ranging from 3.01 cm h−1 in irrigated croplands to 6.22 cm h−1 in RAG (Table 1 and Figure 3f).
All LULCCs showed low volumetric water content at the field capacity (−10 kPa), varying from 0.17 cm3 cm−3 to 0.27 cm3 cm−3. Although LULCCs present different values of field capacity ( θ f c ) and wilt point ( θ w p ), the average difference between θ f c and θ w p for each LULCC is typically around 0.06 cm3 cm−3, highlighting the low water retention capacity for these soils (Table 1).
Table 1 shows all hydraulic and physical properties for the Campbell and Norman model. Soils under natural vegetation presented lower air entry potential than soils under agriculture or pasture, with average ψ e values equal to 0.71 kPa for CDO and 0.87 kPa for FOR. The highest values of air entry matric potential were found in soils under pasture ( ψ e = 1.80 kPa) and irrigated agriculture ( ψ e = 1.46 kPa), which are associated to higher compaction (Table 1 and Figure 3). The Campbell b parameter average values ranged between −4.06 and −5.30, which are similar to literature values for sandy, loamy sand and sandy loam soils.

3.2. Pedotransfer Functions

All PTFs used to estimate ψ e , b, K s a t , θ f c , and θ p a w for soils under PAST and RAG were significant at α = 0.05. The PTFs for IRR, CDO, and FOR showed significance only for a few parameters (Table 2). ψ e was significant in all LULCC, while b was significant only in CDO, PAST and RAG. The PTFs for K s a t were not significant for any LULCCs, except PAST and RAG. θ f c PTFs were significant only for CDO, PAST and RAG, while θ w p , PTFs were significant for FOR, PAST and RAG (Table 2).
The PTF performance to estimate of the ψ e was generally good, showing the maximum adjustment in FOR areas with r = 0.98, followed by PAST (r = 0.93), RAG (r = 0.84) and IRR (r = 0.84) (Table 2). The lowest correlation was found for CDO areas with r = 0.55. For this LUC, the average log ψ e used for calibration and validation were significantly different, which may have contributed to the lower performance in ψ e estimate.
For the Campbell b parameter, the performance of the PTF estimative showed poor correlation for CDO, IRR, and PAST, with correlations 0.36, −0.69 and 0.24, respectively. However, for RAG and FOR areas, the PTFs showed r > 0.70 agreement between the estimated and observed Campbell b parameter (Table 2).
The soil water content represented by θ f c , and θ w p also had a poor correlation for IRR and RAG agriculture LUs with r < 0.20. In CDO, FOR and PAST, the estimated parameters had a correlation >0.45. The worst correlations were for K s a t for all LULCCs with negative r values.

4. Discussion

The sand fraction is predominant in the granulometric composition of the soils, the clay content varies between 10% and 41% and the silt fraction presents the lower values, typically below 4% (Figure 2). This granulometric composition is the result of the high rate of weathering during the genesis of these soils and the soil source material—in this case, the sandstones of the Urucuia Formation [45,46,47,48].
The general trends showed higher values of ρ a and lower values of ϕ t o t in agriculture areas compared to native vegetation, revealing that soils under agricultural land use are slightly compacted in Western Bahia. This is in line with Fontana et al. [48], who found ρ a in agriculture areas in the region of Luis Eduardo Magalhães equal or above 1.52 g cm−3. Cunha et al. [46] also found increased ρ a in different crop areas with different periods of cultivation, demonstrating that duration of land use also has an influence on soil structure. Naturally, the Cerrado Western Bahia soils present a cohesive sub-surface horizon [46,48,49,50], making it more susceptible to compaction by grazing, mechanization and management applied at the soil surface. This transitional or subsurface horizon was observed and characterized by a slight increase in ρ a , generally within 10–30 cm [48]. Indeed, our results show the natural cohesive sub-surface in the CDO and FOR areas (Figure 4b), and the intensification of the increase of ρ a in agriculture land uses. The ρ a in RAG was twice as high in the 10–15 cm layer in relation to the 0–5 cm surface layer. In IRR areas, there was also an increase in ρ a in these layers, although with less intensity than observed in RAG areas (Figure 4).
In Western Bahia, where the irrigated and rainfed croplands have a key role in the development of agriculture, this natural sub-surface horizon cohesion might be a concern factor for the maintenance of the rates of infiltration, the permeability and the water availability to the croplands and recharge of the aquifer. In general, a reduction of ϕ t o t , and an increase of ρ a may be occurring due to a combination of agricultural implements used to remove the natural vegetation cover and the applied agriculture management, consequently reducing K s a t values.
During the fieldwork, the farmers reported that, after removing the Cerrado, it is mandatory to use the conventional tillage system, with the plowing, complemented by the subsoilers and scarifiers, to break any physical impediments. According to our interviews, after the conversion to agriculture, most farmers adopt a planting cycle system, where for the four years after the removal of Cerrado, the no-till system is applied and in the fifth year the conventional tillage system is used, characterizing a mixed management. However, we do not have enough long-term data to claim that this management is the best practice to adopt in order to avoid cohesive sub-surfaces [51,52], but we evaluated the effects of surface and subsurface compaction in hydraulic properties in IRR areas (Table 3). In these areas, where there are two or more crops planted per year, and the traffic of machinery is more intense, ρ a at the 15–20 cm layer is 1.67 g cm−3, significantly higher than at the surface (0–5 cm) layer (1.57 g cm−3) (Table 3), which is a vertical pattern similar to all the other LULCCs (Figure 4). While the difference in ρ a between surface and subsurface layers is significant, it is not sufficient to influence the hydraulic parameters, such as K s a t , ϕ t o t , θ w p and θ p a w (Table 3), although ψ e and θ f c are significantly different in the vertical.
The spatial variability of Ksat in soils under natural vegetation was also extremely high. In CDO, Ksat varied between 224.35 cm h−1 and 1.14 cm h−1, while, in FOR, it varied between 376.69 cm h−1 and 4.86 cm h−1. For Cerrado areas in the MATOPIBA, other studies have found K s a t in the range between 403.8 cm h−1 [48] and 5.26 cm h−1 [20].
In this study, the variability of K s a t in the natural areas can be explained by the higher heterogeneity of the soils, and, consequently, different accumulation of organic matter, litter, tree density, the soil fauna, source material and root systems acting on the soils structures. Despite the soil compaction, the average K s a t for agriculture areas can be considered higher when compared to K s a t for latossolos under conventional tillage system in Goiás, K s a t = 0.535 cm h−1 [53]. For other agriculture areas, in the Cerrado biome, the K s a t values presented in the literature range between 5.41 cm h−1 [20], and 15.47 cm h−1 [48]. Likely, the use of rotation among maize, soybean, cotton and other croplands, in addition to the mixed management in Western Bahia, contribute to the maintenance of the high rates of K s a t in agriculture land uses areas even with the presence of soil compaction.

5. Conclusions

This study analyzed hydraulic and physical soil properties in Western Bahia at a local scale and in different land uses and land covers. Significant changes were found between some soil properties in agricultural areas and natural vegetation cover, indicating that agricultural activity can influence the soil properties.
The agriculture land use increased soil bulk density at soil surface and subsurface, reducing the K s a t by an order of magnitude in relation to Cerrado and Forest areas, and also decreasing the soil porosity. Despite the reduction, K s a t still ranges between 30 mm h−1 and 62 mm h−1, which are still considered high hydrological infiltration rates. Thus, our results reveal that in Western Bahia the agriculture land use areas do not affect directly the vertical nature of hydrological flowpath for visited areas, but, in the case of very intense precipitation events, the K s a t reduction may lead to erosive processes favring nutrient and soil losses. In Western Bahia, however, the farmers are very interested in adopting sustainable practices that preserve the soil quality, investing in state-of-the-art technology, increasing the intervals of soil revolving and implementing crop rotation system.
Nonetheless, one must be careful to not extrapolate these results for all Western Bahia region, since we do not know how irrigated agriculture, rainfed agriculture, and pasture influence the hydraulic and physical properties in the long-term, and how the management (conventional tillage or no-tillage) can influence these properties in loamy/sand soils.
In the literature, several studies monitoring physical soil properties show that reduced tillage practices have promising results for soil moisture conservation and for crops growth. However, the number of scientific studies remains low for the Cerrado agricultural frontier. We emphasize that the impacts of land use change on hydraulic and physic soil properties are a reality in MATOPIBA and there is an urgent need to monitor how these changes occur over time in order to develop effective mitigation strategies of soil and water conservation.

Author Contributions

Experiment design, M.H.C.; data collection, E.A.D.; formal analysis, E.A.D.; data curation, M.H.C.; writing—original draft preparation, E.A.D.; writing—review and editing, E.A.D. and M.H.C.; visualization, E.A.D.; supervision, M.H.C.; project administration, M.H.C.; funding acquisition, M.H.C.

Funding

This study was financed in part by the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior—Brasil (CAPES)—Finance Code 001 and PRODEAGRO (grant 011/2016).

Acknowledgments

We thank the Associação de Agricultores e Irrigantes da Bahia (AIBA-team) for all logistical support and the opportunity to know the Western Bahia region. For dedicated help during fieldwork, we would like to thank Glauciana, Jonatas, Jales, Samuel, Anderson, Sérgio, Eneas, Liliane, Lucas Lima and José Domingos and several other short-term helpers. In particular, we would to thank all farmers who welcomed our team on to their farms.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study, in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. De Oliveira, S.N.; de Carvalho Júnior, O.A.; Gomes, R.A.T.; Guimarães, R.F.; McManus, C.M. Landscape-fragmentation change due to recent agricultural expansion in the Brazilian Savanna, Western Bahia, Brazil. Reg. Environ. Chang. 2017, 17, 411–423. [Google Scholar] [CrossRef]
  2. Nunes de Oliveira, S.; Abílio de Carvalho Júnior, O.; Trancoso Gomes, R.A.; Fontes Guimarães, R.; McManus, C.M. Deforestation analysis in protected areas and scenario simulation for structural corridors in the agricultural frontier of Western Bahia, Brazil. Land Use Policy 2017, 61, 40–52. [Google Scholar] [CrossRef]
  3. Brannstrom, C.; Jepson, W.; Filippi, A.M.; Redo, D.; Xu, Z.; Ganesh, S. Land change in the Brazilian Savanna (Cerrado), 1986-2002: Comparative analysis and implications for land-use policy. Land Use Policy 2008. [Google Scholar] [CrossRef]
  4. Batistella, M.; Valladares, G.S. Farming expansion and land degradation in Western Bahia, Brazil. Biota Neotrop. 2009. [Google Scholar] [CrossRef]
  5. Spera, S.A.; Galford, G.L.; Coe, M.T.; Macedo, M.N.; Mustard, J.F. Land-use change affects water recycling in Brazil’s last agricultural frontier. Global Chang. Biol. 2016, 22, 3405–3413. [Google Scholar] [CrossRef] [PubMed]
  6. Filho, A.C.; Costa, K. A Expansão da Soja no Cerrado—Caminhos Para a Ocupação Territorial, Uso do Solo e Produção Sustentável; Technical Report; AGROICONE: São Paulo, Brazil, 2016. [Google Scholar]
  7. Sano, E.E.; dos Santos, C.C.M.; da Silva, E.M.; Chaves, J.M. Fronteira agrícola do oeste baiano: Considerações sobre os aspectos temporais e ambientais. Geociencias 2011, 30, 479–489. [Google Scholar]
  8. AIBA. Anuário Agropecuário do Oeste da Bahia Safra 2014/2015; AIBA: Barreiras, Brazil, 2015; p. 64. [Google Scholar]
  9. Batistella, M.; Guimarães, M.; de Miranda, E.E.; Vieira, H.R.; Valladares, G.S.; Mangabeira, J.A.d.C.; de Assis, M.C. Monitoramento da expansão agropecuária na reião Oeste da Bahia. Embr. Monit. Satél. 2002, 20, 41. [Google Scholar]
  10. AIBA. Anuário Agropecuário do Oeste da Bahia Safra 2016/2017; AIBA: Barreiras, Brazil, 2017; p. 53. [Google Scholar]
  11. Elsenbeer, H.; Newton, B.E.; Dunne, T.; De Moraes, J.M. Soil hydraulic conductivities of latosols under pasture, forest and teak in Rondonia, Brazil. Hydrol. Process. 1999. [Google Scholar] [CrossRef]
  12. Muller, M.M.L.; Guimarães, M.D.F.; Desjardins, T.; Martins, P.F.D.S. Degradação de pastagens na Região Amazônica: Propriedades físicas do solo e crescimento de raízes. Pesquisa Agropec. Bras. 2001, 36, 1409–1418. [Google Scholar] [CrossRef]
  13. Scheffler, R.; Neill, C.; Krusche, A.V.; Elsenbeer, H. Soil hydraulic response to land-use change associated with the recent soybean expansion at the Amazon agricultural frontier. Agric. Ecosyst. Environ. 2011, 144, 281–289. [Google Scholar] [CrossRef] [Green Version]
  14. Lal, R. Deforestation and land-use effects on soil degradation and rehabilitation in western Nigeria 1. Soil physical and hydrological properties. Land Degrad. Dev. 1996, 7, 19–45. [Google Scholar] [CrossRef]
  15. De Moraes, J.F.; Volkoff, B.; Cerri, C.C.; Bernoux, M. Soil properties under Amazon forest and changes due to pasture installation in Rondônia, Brazil. Geoderma 1996, 70, 63–81. [Google Scholar] [CrossRef]
  16. Zimmermann, B.; Elsenbeer, H.; De Moraes, J.M. The influence of land-use changes on soil hydraulic properties: Implications for runoff generation. For. Ecol. Manag. 2006, 222, 29–38. [Google Scholar] [CrossRef]
  17. Horel, Á.; Tóth, E.; Gelybó, G.; Kása, I.; Bakacsi, Z.; Farkas, C. Effects of Land Use and Management on Soil Hydraulic Properties. Open Geosci. 2015, 7. [Google Scholar] [CrossRef] [Green Version]
  18. De Almeida, V.P.; Alves, M.C.; Da Silva, E.C.; De Oliveira, S.A. Rotação de culturas e propriedades físicas e químicas em latossolo vermelho de cerrado sob preparo convencional e semeadura direta em adoção. Rev. Bras. Cienc. Solo 2008, 32, 1227–1237. [Google Scholar] [CrossRef]
  19. De Faria, A.F.G.; dos Santos, A.C.; dos Santos, T.M.; Batistella Filho, F. Influência do manejo do solo nas propriedades químicas e físicas em topossequência na bacia do rio Araguaia, Estado do Tocantins. Rev. Bras. Cienc. Solo 2010, 34, 517–524. [Google Scholar] [CrossRef] [Green Version]
  20. Filizola, H.F.; Fontana, A.; Donagemma, G.K.; de Souza, M.D.; Bortolon, E.S.O.; Bortolon, L. Qualidade física de solos influenciada pelo uso e manejo na região de Guaraí-TO. Bol. Pesquisa Desenvol. 2017, 72, 35. [Google Scholar]
  21. Bonetti, J.d.A.; Paulino, H.B.; de Souza, E.D.; Carneiro, M.A.C.; Caetano, J.O. Soil physical and biological properties in an integrated crop-livestock system in the Brazilian Cerrado. Pesquisa Agropec. Bras. 2018, 53, 1239–1247. [Google Scholar] [CrossRef]
  22. Nóbrega, R.L.B.; Guzha, A.C.; Torres, G.N.; Kovacs, K.; Lamparter, G.; Amorim, R.S.S.; Couto, E.; Gerold, G. Effects of conversion of native cerrado vegetation to pasture on soil hydro-physical properties, evapotranspiration and streamflow on the Amazonian agricultural frontier. PLoS ONE 2017, 12, e0179414. [Google Scholar] [CrossRef] [PubMed]
  23. Anache, J.A.; Flanagan, D.C.; Srivastava, A.; Wendland, E.C. Land use and climate change impacts on runoff and soil erosion at the hillslope scale in the Brazilian Cerrado. Sci. Total Environ. 2018, 622–623, 140–151. [Google Scholar] [CrossRef]
  24. Hunke, P.; Mueller, E.N.; Schröder, B.; Zeilhofer, P. The Brazilian Cerrado: assessment of water and soil degradation in catchments under intensive agricultural use. Ecohydrology 2015, 8, 1154–1180. [Google Scholar] [CrossRef]
  25. Dias, L.C.P.; Macedo, M.N.; Costa, M.H.; Coe, M.T.; Neill, C. Effects of land cover change on evapotranspiration and streamflow of small catchments in the Upper Xingu River Basin, Central Brazil. J. Hydrol. Reg. Stud. 2015, 4, 108–122. [Google Scholar] [CrossRef] [Green Version]
  26. Moreira, V.S.; Candido, A.; Roberti, D.R.; Webler, G.; Diaz, M.B.; Gon, G.; Gonçalves, D.; Pousa, R.; Degrazia, A. Influence of Soil Properties in Different Management Systems: Estimating Soybean Water Changes in the Agro-IBIS Model. Earth Interact. 2018, 22. [Google Scholar] [CrossRef]
  27. Przeździecki, K.; Zawadzki, J.; Cieszewski, C.; Bettinger, P. Georgia and South Carolina using the triangle method. Silva Fennica 2017, 51, 1–19. [Google Scholar]
  28. Wang, J.; Chagnon, F.J.F.; Williams, E.R.; Betts, A.K.; Renno, N.O.; Machado, L.A.T.; Bisht, G.; Knox, R.; Bras, R.L. Impact of deforestation in the Amazon basin on cloud climatology. Proc. Natl. Acad. Sci. USA 2009, 106, 3670–3674. [Google Scholar] [CrossRef] [Green Version]
  29. Malhado, A.C.M.; Pires, G.F.; Costa, M.H. Cerrado conservation is essential to protect the Amazon rainforest. Ambio 2010, 39, 580–584. [Google Scholar] [CrossRef]
  30. Coe, M.T.; Marthews, T.R.; Costa, M.H.; Galbraith, D.R.; Greenglass, N.L.; Imbuzeiro, H.M.; Levine, N.M.; Malhi, Y.; Moorcroft, P.R.; Muza, M.N.; et al. Deforestation and climate feedbacks threaten the ecological integrity of south-southeastern Amazonia. Philos. Trans. R. Soc. Lond. Ser. B Biol. Sci. 2013, 368, 20120155. [Google Scholar] [CrossRef]
  31. Pires, G.F.; Costa, M.H. Deforestation causes different subregional effects on the Amazon bioclimatic equilibrium. Geophys. Res. Lett. 2013, 40, 3618–3623. [Google Scholar] [CrossRef] [Green Version]
  32. IBGE. Manual Técnico da Vegetação Brasileira (Manuais Técnicos em Geociências n. 1); Fundação Instituto Brasileiro de Geografia e Estatística (IBGE): Rio de Janeiro, RJ, Brazil, 1992.
  33. Dos Santos, H.G.; Jacomine, P.K.T.; dos Anjos, L.H.C.; de Oliveira, V.A.; Lumbreras, J.F.; Coelho, M.R.; de Almedia, J.A.; de Araujo Filho, J.C.; de Oliveria, J.B.; Cunha, T.J.F. Brazilian Soil Classification System, 5th re ed.; EMBRAPA Solos: Brasília, DF, Brazil, 2018. [Google Scholar]
  34. Youngs, E. Hydraulic Conductivity of Saturated Soils. In Soil Analysis: Physical Methods; Marcel Dekker Inc.: New York, NY, USA, 1991; pp. 161–207. [Google Scholar]
  35. Ruiz, H.A. Incremento da exatidão da análise granulométrica do solo por meio da coleta da suspensão (silte + argila). Rev. Bras. Cien. Solo 2005, 29, 297–300. [Google Scholar] [CrossRef]
  36. Richards, L.A.; Fireman, M. Pressure-plate apparatus for measuring moisture sorption and transmission by soils. Soil Sci. 1943, 56, 395–404. [Google Scholar] [CrossRef]
  37. ClassenL, M.E.C.; Barreto, W.D.O.; Paula, J.L.D.; Duarte, M.N. Manual de Métodos de Análise de Solo; EMBRAPA; Centro Nacional de Pesquisa de Solos: Rio de Janeiro, Brazil, 1997; Volume 2, pp. 7–9. [Google Scholar]
  38. Campbell, G.S.; Norman, J.M. An Introduction to Environmental Biophysics, 2nd ed.; Springer: New York, NY, USA, 1998; Volume 6, p. 474. [Google Scholar]
  39. Clapp, R.B.; Hornberger, G.M. Empirical equations for some soil hydraulic properties. Water Resour. Res. 1978, 14, 601–604. [Google Scholar] [CrossRef]
  40. Moldrup, P.; Olesen, T.; Komatsu, T.; Schjønning, P.; Rolston, D. Tortuosity, Diffusivity, and Permeability in the Soil Liquid and Gaseous Phases. Soil Sci. Soc. Am. J. 2001, 65, 613. [Google Scholar] [CrossRef]
  41. Van Looy, K.; Bouma, J.; Herbst, M.; Koestel, J.; Minasny, B.; Mishra, U.; Montzka, C.; Nemes, A.; Pachepsky, Y.A.; Padarian, J.; et al. Pedotransfer Functions in Earth System Science: Challenges and Perspectives. Rev. Geophys. 2017, 55, 1199–1256. [Google Scholar] [CrossRef] [Green Version]
  42. Pittaki-Chrysodonta, Z.; Moldrup, P.; Knadel, M.; Iversen, B.V.; Hermansen, C.; Greve, M.H.; de Jonge, L.W. Predicting the Campbell Soil Water Retention Function: Comparing Visible–Near-Infrared Spectroscopy with Classical Pedotransfer Function. Vadose Zone J. 2018, 17. [Google Scholar] [CrossRef]
  43. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2018. [Google Scholar]
  44. Wilk, M.B.; Shapiro, S.S. An Analysis of Variance Test for Normality. Biometrika 1965, 52, 591–611. [Google Scholar]
  45. Jacomine, P.K.T.; Cavalcanti, A.C.; Ribeiro, M.R.; Montenegro, J.O.; Burgos, N.; Filho, T.F.R.d.M.; Formiga, R.A. Levantamento exploratório - Reconhecimento de solos da margem esquerda do Rio São Francisco da Bahia. Bol. Técnico EMBRAPA Rec. 1976, 38, 1–439. [Google Scholar]
  46. Cunha, T.J.F.; Macedo, J.R.; Ribeiro, L.P.; Palmieri, F.; de Freitas, P.L.; Aguiar, A.d.C. Impacto do manejo convencional sobre propriedades físicas e substâncias húmicas de solos sob cerrado. Ciênc. Rural 2001. [Google Scholar] [CrossRef]
  47. Fontana, A.; de Oliveira, A.P. Levantamento e caracterização de solos frágeis no Município de Luís Eduardo Magalhães, BA. In Solos Frágeis: Caracterização, Manejo e Sustentabilidade; de Castor, S.S., Hernani, L.C., Eds.; EMBRAPA Solos: Brasília, DF, Brazil, 2015; Chapter 3. [Google Scholar]
  48. Fontana, A.; Teixeira, W.G.; Balieiro, F.d.C.; de Moura, T.P.A.; de Menezes, A.R.; Santana, C.I. Características e atributos de Latossolos sob diferentes usos na região Oeste do Estado da Bahia. Pesquisa Agropec. Bras. 2016. [Google Scholar] [CrossRef]
  49. Giarola, N.F.B.; Silva, A.P.; Tormena, C.; Souza, L.S.; Ribeiro, L.P. Similaridades entre o caráter coeso dos solos e o comportamento hardsetting: Estudo de caso. Revista Bras. Ciênc. Solo 2001, 25, 239–247. [Google Scholar] [CrossRef]
  50. Giarola, N.F.B.; Silva, A.P.D. Conceitos sobre solos coesos e hardsetting. Sci. Agric. 2002, 59, 613–620. [Google Scholar] [CrossRef] [Green Version]
  51. Stone, L.F.; Silveira, P.M. Efeitos do sistema de preparo e da rotação de culturas na porosidade e densidade do solo. Rev. Bras. Ciênc. Solo 2001, 25, 395–401. [Google Scholar] [CrossRef] [Green Version]
  52. Santos, G.G.; Marchão, R.L.; da Silva, E.M.; da Silveira, P.M.; Becquer, T. Qualidade física do solo sob sistemas de integração lavoura-pecuária. Pesquisa Agropec. Bras. 2011, 46, 1339–1348. [Google Scholar] [CrossRef] [Green Version]
  53. Da Silva, N.F.; Nobre, F.; de Carvalho, R.C.; Filho, F.R.C.; Teixeira, M.B.; Joaquim, J. Características físico-hídricas de um latossolo sob diffrentes sistemas de manejo. Rev. Bras. Agric. Irrigada 2014, 8, 375–390. [Google Scholar]
Figure 1. Study area and location of the soil samples collected considering different land use and land cover in Western Bahia.
Figure 1. Study area and location of the soil samples collected considering different land use and land cover in Western Bahia.
Agriculture 09 00024 g001
Figure 2. Average soil texture fractions for 0–5 cm layer in Western Bahia according to the USDA soil classification and to the Brazilian soil classification system—SiBCS.
Figure 2. Average soil texture fractions for 0–5 cm layer in Western Bahia according to the USDA soil classification and to the Brazilian soil classification system—SiBCS.
Agriculture 09 00024 g002
Figure 3. Soil hydraulic and physical properties considering 0–5 cm depth for different LULCCs in Western Bahia.In the box plots, the lower limit of the box indicates the 25th percentile, the black line within the box marks the median, the red point within the box marks the mean, and the upper limit of the box indicates the 75th percentile. Bars above and below the box indicate the confidence interval. The samples are distributed in Cerrado areas (CDO), n = 22; Forest formations (FOR), n = 19; Irrigated agriculture (IRR), n= 20; Pasture (PAST), n = 21; and Rainfed agriculture (RAG), n = 20. Different letters means that averages are statistically different according to Tukey Cramer test at α = 0.05.
Figure 3. Soil hydraulic and physical properties considering 0–5 cm depth for different LULCCs in Western Bahia.In the box plots, the lower limit of the box indicates the 25th percentile, the black line within the box marks the median, the red point within the box marks the mean, and the upper limit of the box indicates the 75th percentile. Bars above and below the box indicate the confidence interval. The samples are distributed in Cerrado areas (CDO), n = 22; Forest formations (FOR), n = 19; Irrigated agriculture (IRR), n= 20; Pasture (PAST), n = 21; and Rainfed agriculture (RAG), n = 20. Different letters means that averages are statistically different according to Tukey Cramer test at α = 0.05.
Agriculture 09 00024 g003
Figure 4. Profile average of soil bulk density standard deviation for different land use and land cover classes. The total samples are 714 distributed in: Cerrado areas (CDO), n = 22; Forest formations (FOR), n = 19; Irrigated agriculture (IRR), n= 20; Pasture (PAST), n = 21; and Rainfed agriculture areas (RAG), n = 20, multiplied by 7 depths.
Figure 4. Profile average of soil bulk density standard deviation for different land use and land cover classes. The total samples are 714 distributed in: Cerrado areas (CDO), n = 22; Forest formations (FOR), n = 19; Irrigated agriculture (IRR), n= 20; Pasture (PAST), n = 21; and Rainfed agriculture areas (RAG), n = 20, multiplied by 7 depths.
Agriculture 09 00024 g004
Table 1. Average soil physical parameters for 0–5 cm layer under different land use cover in the Western of Bahia used in this study.
Table 1. Average soil physical parameters for 0–5 cm layer under different land use cover in the Western of Bahia used in this study.
LULCC ψ e b ρ a K sat ρ s ϕ tot θ fc θ wp θ paw
CDO0.71−4.181.3616.292.530.48130.12740.08480.0425
FOR0.87−5.101.3514.472.520.46190.17890.11090.0680
IRR1.46−4.261.573.012.610.39910.14940.08580.0636
PAST1.80−5.001.615.102.580.37620.16630.10520.0611
RAG1.20−5.301.526.222.590.41080.17190.11770.0542
ψ e : soil air potential entry, kPa; b: Campbell parameter; ρ a : soil bulk density, g cm−3; K s a t : saturated hydraulic conductivity, cm h−1; ρ s : soil particle density, g cm−3; ϕ t o t : soil total porosity, fraction; θ f c : volumetric moisture at 10 kPa, cm3 cm−3; θ w p : volumetric moisture at 1500 kPa, cm3 cm−3; θ p a w : volumetric moisture available to plants, cm3 cm−3
Table 2. Pedotransfer functions for Western Bahia.
Table 2. Pedotransfer functions for Western Bahia.
Pedotransfer Function—PTFValidation
yEquationR2Fp-Valuer
IRR—Irrigated Agriculture
log ψ e 3.5824 − 1.2283 f c s a − 7.7754 ϕ t o t 0.598.600.0050.84
b−3.7421 − 1.6219 f c s a 0.172.710.123−0.69
K s a t −0.8328 + 3.2159 ϕ t o t 0.061.950.186−0.008
θ f c 0.86625 − 0.01727 f c s a − 0.21859 ρ a − 0.85690 ϕ t o t 0.080.310.8160.29
θ w p 0.5624 − 0.1561 ρ a − 0.5720 ϕ t o t 0.110.670.5280.20
CDO—Cerrado formations
log ψ e 1.6879 − 0.5922 f c s a − 3.6247 ϕ t o t 0.6310.910.0020.55
b−6.204 + 3.530 f c s a 0.162.710.1220.36
K s a t 1.7913 − 1.2343 f c s a 0.163.580.081−0.15
θ f c 0.46082 − 0.18014 ρ a 0.275.090.0410.80
θ w p 0.24049 − 0.10555 ρ a 0.183.080.1010.68
PAST—Pasture
log ψ e 4.4995 − 1.5082 ρ a − 5.1395 ϕ t o t 0.516.330.0130.93
b−12.938 + 6.370 f c s a + 12.937 ϕ t o t 0.598.530.050.24
K s a t 4.4227 + 1.5955 f c s a − 2.7828 ρ a 0.507.960.0060.58
θ f c 0.59324 − 0.14310 f c s a + 0.18177 ρ a 0.7315.644.54 × 10−40.93
θ w p 0.17058 − 0.14580 f c s a 0.5616.470.0010.45
FOR—Forest Formations
log ψ e 2.5457 − 0.8746 f c s a − 5.3052 ϕ t o t 0.536.190.0160.98
b−36.460 + 7.915 f c s a + 11.956 ρ a + 26.145 ϕ t o t 0.594.810.0250.74
K s a t −0.6224 + 1.2973 f c s a + 2.7719 ϕ t o t 0.182.480.130−0.76
θ f c 1.42937 − 0.28125 f c s a − 0.46640 ρ a − 1.00496 ϕ t o t 0.533.670.0510.79
θ w p 1.24529 − 0.21743 f c s a − 0.43746 ρ a − 0.99032 ϕ t o t 0.687.300.0070.84
RAG—Rainfed Agriculture
log ψ e 11.351 − 4.110 f c s a − 12.448 ϕ t o t 0.475.300.0220.84
b−36.157 + 10.086 f c s a + 7.275 ρ a + 35.835 ϕ t o t 0.8623.034.80 × 10−50.73
K s a t −0.1303 + 1.6849 f c s a 0.214.820.047−0.28
θ f c 2.84267 − 0.26601 f c s a − 0.97267 ρ a −2.40257 ϕ t o t 0.8014.6037.5 × 10−5−0.33
θ w p 1.30028 − 0.20449 f c s a − 0.40620 ρ a − 1.12476 ϕ t o t 0.7611.760.0010.42
r: Pearson correlation coefficient.
Table 3. Hydraulic and soil physical properties for irrigated land use in different depths.
Table 3. Hydraulic and soil physical properties for irrigated land use in different depths.
Depth ψ e b ρ a K sat ρ s ϕ tot θ fc θ wp θ paw
0–51.46 a−4.26 a1.574 a3.00 a2.612 a0.397 a0.169 a0.0858 a0.0840 a
15–200.03 b−4.16 a1.665 b3.07 a2.662 a0.374 a0.172 b0.0902 a0.0818 a
ab Values significantly different according to the t Student test at α = 0.05 are followed by different letters. ψ e : soil air potential entry, kPa; b: Campbell parameter; ρ a : soil bulk density, g cm−3; K s a t : saturated hydraulic conductivity, cm h−1; ρ s : soil particle density, g cm−3; ϕ t o t : soil total porosity; fraction; θ f c : volumetric moisture at 10 kPa, cm3 cm−3; θ w p : volumetric moisture at 1500 kPa, cm3 cm−3; θ p a w : volumetric moisture available to plants, cm3 cm−3.

Share and Cite

MDPI and ACS Style

Dionizio, E.A.; Costa, M.H. Influence of Land Use and Land Cover on Hydraulic and Physical Soil Properties at the Cerrado Agricultural Frontier. Agriculture 2019, 9, 24. https://0-doi-org.brum.beds.ac.uk/10.3390/agriculture9010024

AMA Style

Dionizio EA, Costa MH. Influence of Land Use and Land Cover on Hydraulic and Physical Soil Properties at the Cerrado Agricultural Frontier. Agriculture. 2019; 9(1):24. https://0-doi-org.brum.beds.ac.uk/10.3390/agriculture9010024

Chicago/Turabian Style

Dionizio, Emily Ane, and Marcos Heil Costa. 2019. "Influence of Land Use and Land Cover on Hydraulic and Physical Soil Properties at the Cerrado Agricultural Frontier" Agriculture 9, no. 1: 24. https://0-doi-org.brum.beds.ac.uk/10.3390/agriculture9010024

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop