Next Article in Journal
Blockchain Technology and Related Security Risks: Towards a Seven-Layer Perspective and Taxonomy
Next Article in Special Issue
Research on the Spatiotemporal Dynamic Relationship between Human Activity Intensity and Ecosystem Service Value in the Three Gorges Reservoir Area
Previous Article in Journal
Investigation of Meteorological Effects on Çivril Lake, Turkey, with Sentinel-2 Data on Google Earth Engine Platform
Previous Article in Special Issue
Soil Organic Carbon as Response to Reforestation Age and Land Use Changes: A Qualitative Approach to Ecosystem Services
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Impacts of Land Use Changes on Soil Functions and Water Security: Insights from a Three-Year-Long Study in the Cantareira System, Southeast of Brazil

by
Monna Lysa Teixeira Santana
1,
Vanêssa Lopes de Faria
1,
Samara Martins Barbosa
1,
Milson Evaldo Serafim
1,2,
Alexandre Uezu
3,
Bruno Montoani Silva
1 and
Junior Cesar Avanzi
1,*
1
Department of Soil Science, Federal University of Lavras (UFLA), Trevo Rotatório Professor Edmir Sá Santos, Lavras 37203-202, Brazil
2
Federal Institute of Education, Science and Technology of Mato Grosso (IFMT), Av. dos Ramires, Distrito Industrial, Cáceres 78200-000, Brazil
3
Faculty for Environmental Conservation and Sustainability (ESCAS), IPÊ—Institute for Ecological Research, 47 km Dom Pedro I hwy, Nazare Paulista 12960-000, Brazil
*
Author to whom correspondence should be addressed.
Sustainability 2023, 15(18), 13395; https://0-doi-org.brum.beds.ac.uk/10.3390/su151813395
Submission received: 31 July 2023 / Revised: 28 August 2023 / Accepted: 4 September 2023 / Published: 7 September 2023

Abstract

:
Maintaining soil functions is crucial for human well-being, but there is a lack of integration between soil, water security, ecosystem services, and climate change. To bridge this knowledge gap and address erosion-induced soil and water losses and considering intrinsic impacts of soil structure, a three-year-long study was conducted focused on three dominant soil types (Typic Hapludult, Typic Dystrudept, and Typic Usthortent) combined with different land uses (native forest, eucalyptus plantation, rotational grazing, and extensive grazing) in a critical water supply region for the São Paulo metropolitan area in Southeastern Brazil. Surface runoff, evaluated for erosion resistance, was measured using the Cornell infiltrometer, and soil electrical resistivity tomography estimated soil water content to a depth of 1.5 m for groundwater recharge analysis. Soil hydraulic properties were also measured. The results revealed that native forest soils had higher hydraulic conductivity, particularly in the surface layer, compared to eucalyptus and pastures. Native forests in Typic Hapludult showed a higher runoff rate (200 to 250 mm h−1) due to a naturally dense subsoil layer that negatively impacted water infiltration and recharge with a high erosion potential, therefore reducing the amount of water stored. Typic Usthortent maintained a higher soil water content in pastures than in other land uses and also showed a low rate of water infiltration, resulting in perched water in the surface layer. In Typic Dystrudept, the native forest presented higher hydraulic conductivity (0–5 cm: 115.9 cm h−1) than eucalyptus (0–5 cm: 36.4 cm h−1), rotational grazing (0–5 cm: 19.4 cm h−1), and extensive grazing (0–5 cm: 2.6 cm h−1), but there were no significant differences in soil water content among land uses. This work illustrates the crucial role of native forests in affecting deep water recharge, reducing the soil surface erosion, mainly in soils without naturally subsoil layer, maintaining recharge potential. For Ultisols, pastures preserved soil structure and are therefore less impacted by soil management. With these results, a contribution is made to soil and water conservation, providing support for sustainable management practices in erosion-prone areas.

1. Introduction

The soil is an essential component for human well-being, and its maintenance requires attention to ensure that it provides its ecosystem functions. Although the literature on the role of soils in providing ecosystem services (ESs) has grown recently, there is a lack of integration between soil, water security, and climate change [1,2]. Land use change, such as the conversion of native forests into extensive pastures, has led to the degradation of soil properties and functions, which can reduce its capacity to provide various ESs [3,4,5], including water supply regulation, erosion control, and climate regulation.
Native forests play a critical role in preventing erosion due to their dense cover vegetation and root system which helps in the soil aggregation. Lense et al. [6] reported that a 12% rise in deforestation in the Xingu River watershed resulted in an astonishing 312% surge in erosion rates. This significant correlation underscores the critical role of native vegetation in erosion prevention and highlights the urgent need for its preservation. In addition, native forest is indispensable for global groundwater storage due to their deep-rooted trees which facilitate water transportation to deeper soil layers, enabling more extended storage periods before evaporation or absorption by plants [7,8]. However, several studies have reported that in temperate regions, soil covered with native forest has a lower overall water content during the dry season [9,10,11]. Well-managed pastures can also serve as soil cover and extract substantial amounts of water. Nevertheless, it should be noted that pasture water uptake is primarily limited to the superficial soil layer, specifically within the 0–50 cm depth range [12,13,14].
In addition to land use, intrinsic soil factors, such as texture, structure, and hydraulic conductivity, also play a critical role in determining water permeability variation. Soils that have a subsurface layer of high bulk density may have limitations in water infiltration and, consequently, have more water runoff and lower water storage along the profile. Furthermore, clayey soils retain more water than sandy soils due to the latter presenting higher macropore continuity and greater water transport capacity throughout the entire profile, despite its increased susceptibility to erosion. Understanding these factors is essential for adequately managing soil water resources and promoting water security [15].
The Cantareira Water Supply System is a crucial water production area that provides several ecosystem functions to the metropolitan area of São Paulo, with emphasis on water supply, covering an area of approximately 2300 km2 and extending to 75 cities in the states of Minas Gerais and São Paulo in Southeastern Brazil [16]. Livestock accounts for a large portion of land use within the Cantareira System, comprising 46% of the landscape. Native forests cover 35% of the area, plantation forests cover 16%, and reservoirs and water bodies represent 3%. In the region, there are widespread degraded pastures [16] that are incapable of generating sufficient biomass to ensure effective land preservation. The 2014/2015 water scarcity in the Cantareira System was an event of prolonged drought, which caused serious supply problems for a large population. Despite the regional scope of the event, its impacts were felt throughout the territory [17]. The dry season in the region typically occurs during the winter months from May to September, with June, July, and August being the driest months. As a result of the water shortage, there was water rationing for the population not only in the catchment area but also in the metropolitan area of São Paulo. Therefore, it is essential to constantly monitor the provision of ecosystem services and the mitigation of degrading actions in these areas for water recharge.
Despite the importance of the Cantareira System, the lack of methods for assessing ecosystem services has prevented a clear view of the real impacts of transitions between planted forest and pasture on the soil’s ability to provide such services. Understanding the modification in soil functions due to land use change and intrinsic factors will be fundamental for formulating policies that promote sustainable development, soil conservation, maintenance of quality water supply, and payment schemes for water-related environmental services, as implemented in other parts of Brazil [5,18,19,20].
In this context, soil water infiltration and runoff surface are critical physical processes that respond to changes in soil structure as affected by management practices. Enhancing water infiltration conditions in the soil leads to reduced erosion and promotes underground aquifer recharge [21]. As it is linked to soil water dynamics, it serves as a vital indicator of agronomic performance and potential environmental impacts in a specific area. Once infiltrated, the water content can be monitored to facilitate decision-making and assist in understanding environmental recharge processes. For that, geophysical methods such as soil electrical resistivity (ρ) have been used in soil science to collect non-destructive information on soil heterogeneity, including soil water content [22]. ρ is a function of various soil properties, such as the nature of solid constituents (texture and mineralogy), pore arrangement (porosity, pore size distribution, and connectivity), degree of saturation (water content), solute concentration in the soil solution, and temperature [23,24,25,26]. Furthermore, the acquisition of indicators related to soil hydrology has always been limited by methodological constraints [27]. In this regard, the Cornell sprinkler infiltrometer [28] has emerged as an alternative to simplify the process of measuring infiltration and surface runoff, without compromising the reliability of the results. This infiltrometer enables the monitoring and assessment of soil structural quality, providing valuable insights for soil and water conservation studies.
Therefore, the main hypotheses of this study are the following: (i) surface runoff may vary among different land uses and soil types with lower rates in the native forest; (ii) accurate estimation of soil water content and its spatial–temporal changes through electrical resistivity tomography might require specific calibration tailored to individual soil types; (iii) the spatial–temporal patterns of soil water content might potentially be affected by soil uses, and these patterns may also show discrepancies depending on the soil type. This study aimed (i) to evaluate surface runoff in soils under different uses; (ii) to model ρ–θ relationship in soils for θ estimation; (iii) to evaluate the temporal dynamics in contrasting soil types: Inceptisol, Ultisol, and Entisol; and (iv) to verify the land use effect on maintaining water content throughout the seasons and identify priority areas for preservation. In addition, the information obtained will contribute to the planning and proper management of water resources, promoting water security and ensuring the availability of quality water for the population and for the region’s ecosystems.

2. Materials

Study Area

The study was established in the Cantareira water supply system in Southeastern Brazil, with sampling performed in the municipalities of Joanópolis (22°56′16″ S 46°05′50″ W, altitude of 1200 m), Nazaré Paulista (23°12′20″ S 46°21′12″ W, altitude of 800 m), and Piracaia (23°01′39″ S 46°19′35″ W, altitude of 840 m) (Figure 1). In each area, the representative land uses adopted by local farmers—native forest (NF), eucalyptus (E), rotational grazing (RG), and extensive grazing (EG)—were evaluated. The predominant climate in the region is classified as Cwb (Köppen) with cold and dry winters and hot and humid summers [29]. The average annual rainfall is 1570 mm, and yearly temperatures range from 18 to 20 °C [16]. The soils were classified as Argissolo Vermelho Amarelo, Cambissolo Háplico, and Neossolo Regolítico, according to the Brazilian System of Soil Classification [30], and as Ultisol (Typic Hapludult), Inceptisol (Typic Dystrudept) and Entisol (Typic Usthortent), according to Soil Taxonomy [31] (Figure 1).
In February 2019, disturbed and undisturbed soil were randomly sampled in trenches (40 × 40 × 40 cm) for each land use and soil type for physical parameters (Table 1). In each layer, eight samples were taken at the surface (0–5 cm) and the subsurface (30–35 cm).
Particle size distribution in disturbed samples was determined by the hydrometer method [32]. Undisturbed samples were collected using steel cylinders (2.5 cm in height × 6.3 cm diameter) for the determination of soil bulk density (BD), total porosity (TP), macroporosity (MAC), and microporosity (MIP) [33,34]. The saturated hydraulic conductivity (KS) was measured with a constant-charge permeameter [35] in steel cylinders (8.0 cm in height and 6.4 cm diameter). Additional soil characterization and land use management information for each site is available in Santana et al. [36].
In order to characterize the weather patterns during the experimental campaign, agro-meteorological data were obtained from the Meteorological Database for Teaching and Research (BDMEP) of the National Institute of Meteorology [37]. Figure 2 depicts the weekly precipitation and the average monthly precipitation, which were measured at the two closest meteorological stations to the study areas: Monte Verde, situated in the state of Minas Gerais in proximity to the city of Joanópolis, and the station in Bragança Paulista located in the state of São Paulo, close to the cities of Nazaré Paulista and Piracaia.

3. Methods

3.1. Surface Runoff

Around sampling points, surface runoff (SR) tests were performed, with two replicates for each land use. The SR was performed through a sprinkler infiltrometer (Cornell Sprinkler infiltrometer) described by van Es et al. [28] using harvested rain. The infiltrometer is a portable rainfall simulator consisting of 69 drippers at the bottom with a diameter of 0.00063 m and a length of 0.19 m each, providing a capacity of 20.6 L. It is assembled on a cylinder with a diameter of 0.24 m and is capable of simulating various rain intensities through an air inlet regulation system. A high intensity of rainfall was simulated in order to quantify the runoff.
The surface runoff was determined by Equation (1):
SR = [(Vt × 1000)/(45,730 × t)],
where SR is surface runoff (mm h−1); Vt is water volume collected (mL); 45,730 is ring area (mm2); and t is the time interval (hours) between the runoff collections (utilizing 0.05 h in this study).

3.2. Soil Electrical Resistivity Surveys

3.2.1. Description of ρ Measurements

The measurement of the soil electrical resistivity (ρ) was carried out in the study areas in the period of 2019 (summer, fall, winter, and spring), 2020 (winter and spring), and 2021 (summer and fall). Due to the limitations of COVID-19, it was not possible to carry out the measurements in the summer and fall of 2021 and in the eucalyptus and extensive pasture areas in the spring of 2021 for the Typic Hapludult. The measurement aims to evaluate the volumetric moisture in the soil profile and thus the water storage as a function of time, for each land use and soil type.
The ρ was directly measured with the X5xtal 250 Resistivity Meter (Auto Energia, Alvinópolis, MG, Brazil), with two multimeters in a dipole–dipole arrangement (A–B–M–N), where A–B are current electrodes and M–N are the potential electrodes [22]. The measurements were obtained from a horizontal transect of 7.78 m in length, with a spacing of 0.38 m between electrodes. The evaluated depth ranged from 0 to 1.50 m, with a total of 179 measurements being carried out in each land use. Hence, ρ was calculated by Equations (2) and (3):
ρ = (K × ∆V)/I,
K = 2πa,
where ρ is electrical resistivity (Ω m); K is geometric coefficient; ∆V is potential difference (mV); I is injected electrical current (mA); and a is spacing between electrodes (m).

3.2.2. Description of ρ–θ Calibration Test

To model the ρ–θ relationship, a calibration test was performed on undisturbed soil sampled in situ with PVC (polyvinyl chloride) rigid plastic cylinders (both diameter and height of 0.1 m), as described by Melo et al. [38]. Undisturbed soil was sampled from two depths (0–10 and 35–45 cm depths), which represent the surface and subsurface layer, respectively. Duplicates were collected in each soil type (Typic Hapludult, Typic Dystrudept, and Typic Usthortent).
At the laboratory, a Wenner arrangement (A–M–N–B) of electrodes [22] was carried out with aluminum electrodes spaced 0.019 m apart and inserted at a depth of 0.05 m to determine ρ. Each soil sample was saturated by capillarity for 72 h, and the test was conducted by acquiring a set of ρ measurements and θ—obtained from the PVC soil sample weight and determination of gravimetric water content and the bulk density—as the soil dried by natural evaporation. Therefore, ρ and θ were measured in gradual degrees of water saturation acquired for 29 days with an X5xtal 250 Resistivity Meter (Auto Energia, Alvinópolis, MG, Brazil). The electrical conductivity of water was assumed constant during the measurement process.
After this period, the soil samples were oven-dried at a temperature of 105–110 °C for 48 h to calculate the gravimetric water content (w) and the soil bulk density (BD) [33]. The value of θ was calculated by Equation (4):
θ = w (BD/Dw),
where θ is volumetric soil water content (m3 m−3); w is the gravimetric water content (kg kg−1); BD is the soil bulk density (Mg m−3); and Dw is the density of water (considered as 1 Mg m−3).
To estimate the θ through ρ data, the ρ–θ relationship was modeled for each soil type using the power Equation (5):
θ = a ρb,
where a and b are fitting parameters (dimensionless).

3.3. Statistical Analysis

For the analysis of spatial and in-depth variability of electrical resistivity, the ρ was derived by performing an inversion process executed in RES2Dinvx64 software (Geotomo Software, v.4.08), and 2D images were generated using Surfer 13.5.583 (Golden Software). The triangulation method was applied to interpolate the ρ points and obtain 2D tomograms at each treatment.
In order to verify the accuracy of the ρ–θ models for each treatment, the determination coefficient (R2) and the root mean square error (RMSE) calculated between the observed and predicted θ values were used. The calibration models were evaluated using the R (3.6.1) statistical environment [39].

4. Results and Discussion

4.1. Runoff with Cornell Infiltrometer

Figure 3 illustrates the surface runoff obtained from the Cornell Sprinkle infiltrometer test. For Typic Hapludult (Figure 3A), it was observed that the native forest began producing runoff after 6 min, whereas the other land uses exhibited runoff in the initial minute. Upon reaching a plateau, the land uses showed runoff rates ranging between 200 and 250 mm h−1 for NF and RG, respectively. However, there were no substantial differences in the runoff rates observed among the different land uses. Chittolina et al. [40] found a runoff coefficient ranging from 25 to 37%, while studies conducted in other water supply areas of the Cantareira System found a runoff coefficient of 39% for the Jaguari river sub-basin and ranging from 29 to 42% for the Piracicaba river sub-basin [17].
The native forest for Typic Dystrudept showed low runoff rates (Figure 3B), with values of approximately 5 mm h−1, and the eucalyptus-initiated runoff reached a plateau at 11 mm h−1 after 3 min. In contrast, the grazing areas exhibited considerably higher runoff rates (average of RG: 215 mm h−1 and EG: 370 mm h−1). Studies conducted in Dystrudept in the Atlantic Forest—the Mantiqueira Range [41], near the Cantareira System, predominantly covered by forests and pastures—have shown that the native forest has a significant impact on base flow. The authors suggested that native forests offer more favorable conditions for water infiltration, resulting in erosion reduction, greater groundwater recharge, and, consequently, a higher base flow. Conversely, anthropogenic activities have a negative impact on soil permeability. In contrast, preserved areas such as natural forests have the potential to enhance infiltration, which can lead to a decrease in overland flow and sediment transport. As a result, the high runoff rate observed in pasture areas can be attributed to the replacement of native forests with this type of land use. This phenomenon is associated with the specific features of mountainous regions of the Atlantic Forest, which possesses a thick layer of litter, low soil bulk density, high organic carbon content, and greater biological activity when compared to pasture lands. Grazing lands do not facilitate water infiltration, as reported by de Menezes et al. [42] in the Mantiqueira Range Region, southern, Brazil.
There was no surface runoff observed in the native forest and eucalyptus for Typic Ustorthent (Figure 3C). However, the EG showed the highest runoff rates observed as early as the first minute of the test, reaching a maximum runoff rate of 288 mm h−1, and maintaining a plateau at 245 mm h−1. Nevertheless, RG began producing runoff after 9 min and sustained a runoff rate of approximately 160 mm h−1. It is evident that the adoption of pasture management practices has the potential to reduce surface runoff when compared to conventional grazing. In studies conducted by Oliveira et al. [43], groundwater recharge decreased with increasing vegetation density in Brazilian Cerrado Entisols, where grasslands showed higher groundwater recharge rates when compared to areas with higher vegetation density, above 350 mm year−1. This emphasizes the relevance of incorporating soil type and land use information in hydrological and climatic modeling. It is important to highlight those younger soils, such as Inceptisol (Typic Dystrudept) and Entisol (Typic Usthortent), in which soil structure is not well-developed through a field evaluation, can still play an important role in water infiltration compared to soil that presents a dense subsoil horizon like Ultisol.
Generally, in the forest formation areas, the occurrence of a low surface runoff rate contributes to water recharge and soil preservation, resulting in soil losses below the tolerance limit in the majority of the Cantareira System [44].

4.2. ρ–θ Modeling

The ρ–θ modeling, parameters, and accuracies are shown in Figure 4. The model showed good accuracy with R2 values between 0.79 and 0.91 and low RMSE values (between 0.034 and 0.051). It is known that soil moisture and electrical resistivity have non-linear and inverse mathematical relationships, which can be applied to various model types: power, logarithmic, or even exponential [22,45,46]. However, soil texture is a differential for a better fit of this relationship. Archie’s power equation [47] is often applied to saturated sandy soils due to the relationship between resistivity and the number of pores generated by the arrangement of particles. In studies conducted by Melo et al. [48], the hypothesis that Archie’s Law also applies to tropical clayey soils was accepted, showing a good correlation when using power calibration models. The authors evaluated clayey Oxisols, but the results of this study also showed good results for Entisols, Inceptisols, and Ultisols.

4.3. Seasonal Water Content Distribution

The soil water content estimated by ρ for each monitoring period is presented in Figure 5, Figure 6 and Figure 7. In 2019, for Typic Hapludult (Figure 5), it is possible to observe that the native forest was the land use that kept the water in profile up to 1.50 m deep throughout the summer until the winter seasons (from February to August, respectively), especially during the driest period (i.e., winter). Typic Hapludult was characterized by a dense pedogenetic horizon (starting near 30 cm depth), as demonstrated by a high bulk density (1.59 Mg m−3) in the subsurface layer of the native forest, which may limit surface water infiltration (KS) (Table 1) and a surface runoff rate of 200 mm h−1 (Figure 3A). However, soil preparation for eucalyptus plantation may have contributed to the relief of subsurface consolidation, promoting root development below this restrictive layer, consuming water at shallow depths, and drying out the subsurface soil layer throughout the seasons, indicated by the pale colors in the Figure 5E. Among the pastures, it was possible to visualize that the rotated grazing maintained more water throughout the seasons remaining stationary in the soil up to a depth of 1.0 m, unlike the extensive pasture, which showed greater soil dryness from a depth of 13 cm. This result is consistent with the high bulk density, greater impermeability (<KS) (Table 1), and high runoff rates (Figure 3A) of the extensive pasture, which may have impeded the water’s reach to this depth, indicating drying during the evaluation period. Despite the aggressive and penetrating nature of grass roots, achieving this effect requires the implementation of appropriate management practices. In previous studies, low groundwater recharge potentials were observed for these pastures in this soil, directly related to the common pasture management practice of rotating heavier animals in the rotational area and heifers in the extensive area, leading to an animal stocking rate beyond the ideal and causing soil compaction from animal trampling, resulting in water accumulation indicating impaired drainage and inefficient root uptake with a low rate of evapotranspiration [49]. According to a study conducted by Bispo et al. [50] in a Typic Hapludult in Atlantic Forest in Brazil, it was reported that surface runoff was reduced by 41% in well-managed pastures. Sone et al. [51] suggest integrating crops and livestock as an alternative to extensive pastures, noting a 60% increase in water infiltration and a 50% reduction in soil loss rates.
In the spring of 2019 (October), there was a significant precipitation event (Figure 2), which can be mainly observed by the water recharge in the native forest. The monitoring conducted in August 2020 showed a drier winter with no precipitation when compared to the previous year, in which eucalyptus was the only use that was able to maintain a water content greater than 0.24 m3 m−3 up to a depth of approximately 0.50 m. During the dry period of the previous year, the extensive pasture exhibited the same drying pattern as before and dried out more compared to the rotated pasture, particularly deeper in the soil. In the following summer (March 2021), high precipitation events were observed (Figure 2) after a long dry period, in which eucalyptus stored water along the soil profile up to 1.50 m. Nonetheless, the estimated percolation cannot be considered as water recharge since the water balance analysis only considered the 1.5 m soil depth, while certain tree species are adapted to extracting water from deeper layers [52]. Yu et al. [11] conducted a study on the seasonal variation of soil moisture in medium-textured soil under different land uses to a depth of 5.0 m. They were able to infer that moisture levels stabilized with minor fluctuations below a depth of 2.0 m.
Next monitoring, in fall (June 2021), the beginning of the dry period was observed, with the native forest presenting a water content below 0.10 m3 m−3, maintaining the same behavior during dry seasons. In this year, the rotated grazing showed greater variability compared to the extensive pasture. Despite all climatic anomalies, and due to it being a drier year, it is possible to observe a pattern in the soil water content, the same as in 2019, where the native forest was the driest land use; this may be due to evapotranspiration, as a result of voluminous and dense root system tropical forest [53] or groundwater recharge, due to the good soil water–air relationship in this use. This is indicated in Table 1 by the good porous and KS conditions compared to other uses.
The water content for Typic Dystrudept (Figure 6) showed greater stability in maintaining water within the soil system throughout seasons when compared to Typic Hapludult. The native forest exhibited higher water content in the soil profile and a more homogeneous behavior among the other land uses with soil moisture decreasing with depth, mainly for the second year. The pastures produced greater drainage, indicated by the low water content, particularly during the dry winter season (August), regardless of the monitoring year. Rotated grazing became drier compared to the EG condition. Extensive grazing exhibited an atypical behavior in October 2019. In the summer (March 2021), after a rainy season (Figure 2), it was observed that water remained stagnant in the soil surface of the EG land use. The extensive grazing exhibits low hydraulic conductivity—KS (0–5 cm: 2.61 cm h−1; 30–35 cm: 4.18 cm h−1; Table 1) resulting in surface water retention and inefficient water recharge; meanwhile, the native forest presented the highest values (0–5 cm: 115.90 cm h−1; 30–35 cm: 15.58 cm h−1) and the lowest surface runoff rate. Furthermore, in regions where rainforests are predominant and forest conversion takes place, the substitution of deep-rooted native forest with shallow-rooted pasture can disrupt the hydrological cycle in various ways, like the increase in streamflow [54].
Salemi et al. [55] reported a significantly higher hydraulic conductivity at a depth of 15 cm in studies conducted on young soils (Entisols and Inceptisols) in the Atlantic Forest, when comparing soil under the forest with soils under eucalyptus and pasture areas. The observation made by the authors is consistent with the forest soil exhibiting the lowest recorded soil bulk density (1.19 Mg m−3) at the same depth, as well as the highest degree of aggregation and pore space because it has the highest organic matter content. According to Centeno et al. [56], the enhanced aggregation and pore space in the forest soil have a positive effect on water conductivity.
The variation in water storage for Typic Usthortent is presented in Figure 7. It was possible to observe that pastures maintained a high soil moisture content, likely associated with its low infiltration rate, visualized by the high runoff rate (Figure 3C). Typic Usthortent is a very young soil with stoniness characteristics and a poorly developed structure, which favored the maintenance of water in pastures throughout the monitoring period. In contrast, the deep and aggressive roots of the native forest and eucalyptus trees explored water deeply, resulting in drier soil throughout the seasons, indicating that actual evapotranspiration, potential percolation, and soil water storage variation were the most significant hydrological factors in the water balance.
In addition to root water uptake, planting density is a critical factor in determining the water use efficiency of eucalyptus trees. The influence of planting density on water use efficiency in forestry has been the subject of extensive research, and studies have demonstrated that higher planting densities can lead to increased water use [55]. Nevertheless, the impact of planting density on the ecosystem’s water balance is complex and influenced by multiple factors. In this regard, Hakamada et al. [57] have made a valuable contribution to the understanding of this topic by studying different eucalyptus clones and plant densities in Brazil. Their study reported that planting fast-growing trees in tighter spacings might have an adverse effect on the ecosystem’s water balance, potentially depleting stored soil water. Moreover, the results indicated that forests consume more water than pasture, but the maintenance of forest cover in the watershed is required for soil and water conservation, reducing the impacts of soil erosion on water yield and its quality from springs [49].

5. Conclusions

In a broader context, this study holds significant universal implications, shedding light on the consequences of converting native forests to pastures as a localized agricultural practice. This conversion exhibited a notable surge in potential surface runoff within the soil, amplifying sediment transport towards water bodies. Notably, the Typic Hapludult, eucalyptus roots demonstrated the ability to extract water from deep soil layers by penetrating through the compacted layer. In the less weathered soil, Typic Usthortent, water infiltration was impeded and water remained in the studied soil layer.
The investigation underscored the pivotal role of inherent soil properties in governing water dynamics. Consequently, a soil-specific calibration model was devised to establish a correlation between soil electric resistivity and water content. This model exhibited remarkable accuracy, enabling robust water content estimation.
Furthermore, the research emphasized the impact of land use changes, highlighting the critical role of native forests in facilitating deep water recharge. This is attributed to their greater root penetration capacity inducing biopores, evapotranspiration potential, and lower surface runoff rates. In contrast, pastures with shallow roots can only retain water in surface layers, making them inefficient in recharging groundwater. Therefore, the quality of pastures, determined by grass management and animal stocking rate, is paramount in ensuring water security. Nevertheless, further research beyond the 2.0 m is necessary to enhance the understanding of groundwater recharge.
By comprehensively addressing these interrelated factors, this study not only contributes to regional water resource management strategies but also offers a universal framework for understanding the intricate dynamics between land use, soil properties, and groundwater recharge.

Author Contributions

M.L.T.S.: Conceptualization, Methodology, Formal analysis, Writing—Original Draft, Writing—Review & Editing. V.L.d.F.: Conceptualization, Writing—Review & Editing. S.M.B.: Conceptualization, Writing—Review & Editing. M.E.S.: Conceptualization, Writing—Review & Editing. A.U.: Conceptualization, Methodology, Writing—Review & Editing. B.M.S.: Conceptualization, Methodology, Formal analysis, Writing—Original Draft, Writing—Review & Editing, Supervision. J.C.A.: Conceptualization, Methodology, Formal analysis, Writing—Review & Editing, Supervision. All authors have read and agreed to the published version of the manuscript.

Funding

We appreciate the financial support of Coordination for the Improvement of Higher Education Personnel (CAPES), National Council for Scientific and Technological Development (CNPq—Process 441244/2017-3), Foundation for Research Support of the State of Minas Gerais (FAPEMIG), São Paulo Research Foundation (FAPESP—Process 2019/19429-3), Department of Soil Science at Federal University of Lavras (DCS—UFLA), and Institute for Ecological Research (IPÊ). B.M.S. and J.C.A. thank CNPq for grants number 311743/2021-8 and 307059/2022-7, respectively.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. de Melo, M.C.; Fernandes, L.F.S.; Pissarra, T.C.T.; Valera, C.A.; da Costa, A.M.; Pacheco, F.A.L. The COP27 screened through the lens of global water security. Sci. Total Environ. 2023, 873, 162303. [Google Scholar] [CrossRef]
  2. Pokhrel, Y.; Felfelani, F.; Satoh, Y.; Boulange, J.; Burek, P.; Gädeke, A.; Gerten, D.; Gosling, S.N.; Grillakis, M.; Gudmundsson, L.; et al. Global terrestrial water storage and drought severity under climate change. Nat. Clim. Chang. 2021, 11, 226–233. [Google Scholar] [CrossRef]
  3. De Groot, R.S.; Wilson, M.A.; Boumans, R.M.J. A typology for the classification, description and valuation of ecosystem functions, goods and services. Ecol. Econ. 2002, 41, 393–408. [Google Scholar] [CrossRef]
  4. MEA Millenium Ecosystem Assessment; Island Press: Washington, DC, USA, 2005; p. 1597260401.
  5. Alves, M.A.B.; de Souza, A.P.; de Almeida, F.T.; Hoshide, A.K.; Araújo, H.B.; da Silva, A.F.; de Carvalho, D.F. Effects of Land Use and Cropping on Soil Erosion in Agricultural Frontier Areas in the Cerrado-Amazon Ecotone, Brazil, Using a Rainfall Simulator Experiment. Sustainability 2023, 15, 4954. [Google Scholar] [CrossRef]
  6. Lense, G.H.E.; Avanzi, J.C.; Parreiras, T.C.; Mincato, R.L. Effects of deforestation on water erosion rates in the Amazon region. Rev. Bras. Cienc. Agrar. 2020, 15, 1–7. [Google Scholar] [CrossRef]
  7. Bassiouni, M.; Manzoni, S.; Vico, G. Optimal plant water use strategies explain soil moisture variability. Adv. Water Resour. 2023, 173, 104405. [Google Scholar] [CrossRef]
  8. Rodríguez-Iturbe, I.; Porporato, A. Ecohydrology of Water-Controlled Ecosystems: Soil Moisture and Plant Dynamics; Cambridge University Press: Cambridge, UK, 2007; Volume I–VI. [Google Scholar] [CrossRef]
  9. Jayawickreme, D.H.; Jobbágy, E.G.; Jackson, R.B. Geophysical subsurface imaging for ecological applications. New Phytol. 2014, 201, 1170–1175. [Google Scholar] [CrossRef]
  10. Kim, J.H.; Fourcaud, T.; Jourdan, C.; Maeght, J.-L.; Mao, Z.; Metayer, J.; Meylan, L.; Pierret, A.; Rapidel, B.; Roupsard, O.; et al. Vegetation as a driver of temporal variations in slope stability: The impact of hydrological processes. Geophys. Res. Lett. 2017, 44, 4897–4907. [Google Scholar] [CrossRef]
  11. Yu, B.; Liu, G.; Liu, Q.; Huang, C.; Li, H.; Zhao, Z. Seasonal variation of deep soil moisture under different land uses on the semi-arid Loess Plateau of China. J. Soils Sediments 2019, 19, 1179–1189. [Google Scholar] [CrossRef]
  12. Bengtsson, J.; Bullock, J.M.; Egoh, B.; Everson, C.; Everson, T.; O’Connor, T.; O’Farrell, P.J.; Smith, H.G.; Lindborg, R. Grasslands-more important for ecosystem services than you might think. Ecosphere 2019, 10, e02582. [Google Scholar] [CrossRef]
  13. Milazzo, F.; Francksen, R.M.; Zavattaro, L.; Abdalla, M.; Hejduk, S.; Enri, S.R.; Pittarello, M.; Price, P.N.; Schils, R.L.; Smith, P.; et al. The role of grassland for erosion and flood mitigation in Europe: A meta-analysis. Agric. Ecosyst. Environ. 2023, 348, 108443. [Google Scholar] [CrossRef]
  14. Sala, O.; Paruelo, J. Ecosystem services in grasslands. In Nature’s Services: Societal Dependence on Natural Ecosystems; Daily, G., Ed.; Island Press: Washington, DC, USA, 1997; pp. 237–252. [Google Scholar]
  15. Leul, Y.; Assen, M.; Damene, S.; Legass, A. Effects of land use types on soil quality dynamics in a tropical sub-humid ecosystem, western Ethiopia. Ecol. Indic. 2023, 147, 110024. [Google Scholar] [CrossRef]
  16. Uezu, A.; Sarcinelli, O.; Chiodi, R.; Jenkins, C.; Martins, C. Atlas Dos Serviços Ambientais Do Sistema Cantareira, 1st ed.; IPÊ: São Paulo, Brazil, 2017; Volume 1, ISSN 9788579541131. [Google Scholar]
  17. Domingues, L.M.; da Rocha, H.R. Serial droughts and loss of hydrologic resilience in a subtropical basin: The case of water inflow into the Cantareira reservoir system in Brazil during 2013–2021. J. Hydrol. Reg. Stud. 2022, 44, 101235. [Google Scholar] [CrossRef]
  18. Mamedes, I.; Guerra, A.; Rodrigues, D.B.; Garcia, L.C.; Godoi, R.d.F.; Oliveira, P.T.S. Brazilian payment for environmental services programs emphasize water-related services. Int. Soil Water Conserv. Res. 2023, 11, 276–289. [Google Scholar] [CrossRef]
  19. Ruggiero, P.G.; Metzger, J.P.; Tambosi, L.R.; Nichols, E. Payment for ecosystem services programs in the Brazilian Atlantic Forest: Effective but not enough. Land Use Policy 2018, 82, 283–291. [Google Scholar] [CrossRef]
  20. Taffarello, D.; Srinivasan, R.; Mohor, G.S.; Guimarães, J.L.B.; Calijuri, M.D.C.; Mendiondo, E.M. Modeling freshwater quality scenarios with ecosystem-based adaptation in the headwaters of the Cantareira system, Brazil. Hydrol. Earth Syst. Sci. 2018, 22, 4699–4723. [Google Scholar] [CrossRef]
  21. Luz, C.C.d.S.; de Almeida, W.S.; de Souza, A.P.; Schultz, N.; Anache, J.A.A.; de Carvalho, D.F. Simulated rainfall in Brazil: An alternative for assesment of soil surface processes and an opportunity for technological development. Int. Soil Water Conserv. Res. 2023; in press. [Google Scholar] [CrossRef]
  22. Samouëlian, A.; Cousin, I.; Tabbagh, A.; Bruand, A.; Richard, G. Electrical resistivity survey in soil science: A review. Soil Tillage Res. 2005, 83, 173–193. [Google Scholar] [CrossRef]
  23. Garcia-Montiel, D.C.; Coe, M.T.; Cruz, M.P.; Ferreira, J.N.; da Silva, E.M.; Davidson, E.A. Estimating Seasonal Changes in Volumetric Soil Water Content at Landscape Scales in a Savanna Ecosystem Using Two-Dimensional Resistivity Profiling. Earth Interact. 2008, 12, 1–25. [Google Scholar] [CrossRef]
  24. Grubbs, R.A.; Straw, C.M.; Bowling, W.J.; Radcliffe, D.E.; Taylor, Z.; Henry, G.M. Predicting spatial structure of soil physical and chemical properties of golf course fairways using an apparent electrical conductivity sensor. Precis. Agric. 2019, 20, 496–519. [Google Scholar] [CrossRef]
  25. Jeřábek, J.; Zumr, D.; Dostál, T. Identifying the plough pan position on cultivated soils by measurements of electrical resistivity and penetration resistance. Soil Tillage Res. 2017, 174, 231–240. [Google Scholar] [CrossRef]
  26. Roodposhti, H.R.; Hafizi, M.K.; Kermani, M.R.S.; Nik, M.R.G. Electrical resistivity method for water content and compaction evaluation, a laboratory test on construction material. J. Appl. Geophys. 2019, 168, 49–58. [Google Scholar] [CrossRef]
  27. Cheng, Q.; Chen, X.; Chen, X.; Zhang, Z.; Ling, M. Water infiltration underneath single-ring permeameters and hydraulic conductivity determination. J. Hydrol. 2011, 398, 135–143. [Google Scholar] [CrossRef]
  28. van Es, H.M.; Schindelbeck, R. Field Procedures and Data Analysis for the Cornell Sprinkle Infiltrometer; Cornell University: Ithaca, NY, USA, 2003. [Google Scholar]
  29. Alvares, C.A.; Stape, J.L.; Sentelhas, P.C.; Moraes, G.J.L.; Sparovek, G. Köppen’s climate classification map for Brazil. Meteorol. Z. 2013, 22, 711–728. [Google Scholar] [CrossRef]
  30. Santos, H.G. Sistema Brasileiro de Classificação de Solos, 5th ed.; Embrapa: Brasília, Brazil, 2018. [Google Scholar]
  31. Staff Keys to Soil Taxonomy, 12th ed.; USDA: Washington, DC, USA, 2014; Volume 11, ISBN 0926487221.
  32. Gee, G.W.; Or, D.; Dane, J.H.; Topp, C. Particle-Size Analysis. In Methods of Soil Analysis, Part 4—Physical Methods; John Wiley & Sons, Ltd.: Hoboken, NJ, USA, 2018; pp. 255–293. ISBN 9780891188933. [Google Scholar]
  33. Grossman, R.B.; Reinsch, T.G. Bulk density and linear extensibility. In Methods of Soil Analysis, Part 4—Physical Methods; John Wiley & Sons, Ltd.: Hoboken, NJ, USA, 2002; pp. 201–228. ISBN 2002109389. [Google Scholar]
  34. Reynolds, W.; Drury, C.; Yang, X.; Tan, C. Optimal soil physical quality inferred through structural regression and parameter interactions. Geoderma 2008, 146, 466–474. [Google Scholar] [CrossRef]
  35. Klute, A. Laboratory Measurement of Hydraulic Conductivity of Saturated Soil. In Methods of Soil Analysis, Part 1—Physical and Mineralogical Properties, Including Statistics of Measurement and Sampling; American Society of Agronomy, Inc.: Madison, WI, USA, 1965; pp. 210–221. [Google Scholar] [CrossRef]
  36. Santana, M.L.T.; dos Santos, F.F.; de Carvalho, K.M.; Peixoto, D.S.; Uezu, A.; Avanzi, J.C.; Serafim, M.E.; Nunes, M.R.; van Es, H.M.; Curi, N.; et al. Interactions between land use and soil type drive soil functions, highlighting water recharge potential, in the Cantareira System, Southeast of Brazil. Sci. Total Environ. 2023, 903, 166125. [Google Scholar] [CrossRef] [PubMed]
  37. INMET. Banco de Dados Meteorológicos Para Ensino e Pesquisa. Available online: http://www.inmet.gov.br/portal/index.php?r=bdmep/bdmep (accessed on 11 December 2021).
  38. Melo, L.B.B.; Benevenute, P.A.N.; Barbosa, S.M.; Chiarini, T.P.A.; Oliveira, G.C.; Lima, J.M.; Vanella, D.; Consoli, S.; Ferreira, E.A.; Silva, B.M. Spatial and temporal electrical resistivity dynamics in a dense Ultisol under deep tillage and different citrus root-stocks. Soil Tillage Res. 2023, 228, 105629. [Google Scholar] [CrossRef]
  39. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2022; Available online: https://www.r-project.org (accessed on 15 November 2022).
  40. Chittolina, M.; da Rocha, H.R.; Domingues, L.M.; Lobo, G.d.A. Hydrological response of a headwater catchment in Southeast Brazil—Threshold patterns of stormflow response. Hydrol. Process. 2023, 37, e14879. [Google Scholar] [CrossRef]
  41. Pinto, L.C.; de Mello, C.R.; Norton, L.D.; Silva, S.H.G.; Taveira, L.R.S.; Curi, N. Land-use effect on hydropedology in a mountainous region of Southeastern Brazil. Ciência Agrotecnol. 2017, 41, 413–427. [Google Scholar] [CrossRef]
  42. de Menezes, M.D.; Silva, S.H.G.; de Mello, C.R.; Owens, P.R.; Curi, N. Spatial prediction of soil properties in two contrasting physiographic regions in Brazil. Sci. Agric. 2016, 73, 274–285. [Google Scholar] [CrossRef]
  43. Oliveira, P.T.S.; Leite, M.B.; Mattos, T.; Nearing, M.A.; Scott, R.L.; Xavier, R.d.O.; Matos, D.M.d.S.; Wendland, E. Groundwater recharge decrease with increased vegetation density in the Brazilian cerrado. Ecohydrology 2017, 10, e1759. [Google Scholar] [CrossRef]
  44. Lense, G.H.E.; Lämmle, L.; Ayer, J.E.B.; Lama, G.F.C.; Rubira, F.G.; Mincato, R.L. Modeling of Soil Loss by Water Erosion and Its Impacts on the Cantareira System, Brazil. Water 2023, 15, 1490. [Google Scholar] [CrossRef]
  45. Cosenza, P.; Marmet, E.; Rejiba, F.; Cui, Y.J.; Tabbagh, A.; Charlery, Y. Correlations between geotechnical and electrical data: A case study at Garchy in France. J. Appl. Geophys. 2006, 60, 165–178. [Google Scholar] [CrossRef]
  46. McCarter, W.J. The electrical resistivity characteristics of compacted clays. Geotechnique 1984, 34, 263–267. [Google Scholar] [CrossRef]
  47. Archie, G.E. The Electrical Resistivity Log as an Aid in Determining Some Reservoir Characteristics. Trans. AIME 1942, 146, 54–62. [Google Scholar] [CrossRef]
  48. de Melo, L.B.B.; Silva, B.M.; Peixoto, D.S.; Chiarini, T.P.A.; de Oliveira, G.C.; Curi, N. Effect of compaction on the relationship between electrical resistivity and soil water content in Oxisol. Soil Tillage Res. 2021, 208, 104876. [Google Scholar] [CrossRef]
  49. Alvarenga, L.; de Mello, C.; Colombo, A.; Cuartas, L.; Bowling, L. Assessment of land cover change on the hydrology of a Brazilian headwater watershed using the Distributed Hydrology-Soil-Vegetation Model. CATENA 2016, 143, 7–17. [Google Scholar] [CrossRef]
  50. Bispo, D.F.A.; Guimarães, D.V.; Marques, J.J.G.d.S.e.M.; Beniaich, A.; Acuña-Guzman, S.F.; Silva, M.L.N.; Curi, N. Soil Organic Carbon as Response to Reforestation Age and Land Use Changes: A Qualitative Approach to Ecosystem Services. Sustainability 2023, 15, 6863. [Google Scholar] [CrossRef]
  51. Sone, J.S.; de Oliveira, P.T.S.; Zamboni, P.A.P.; Vieira, N.O.M.; Carvalho, G.A.; Macedo, M.C.M.; de Araujo, A.R.; Montagner, D.B.; Sobrinho, T.A. Effects of Long-Term Crop-Livestock-Forestry Systems on Soil Erosion and Water Infiltration in a Brazilian Cerrado Site. Sustainability 2019, 11, 5339. [Google Scholar] [CrossRef]
  52. Rodrigues, A.F.; de Mello, C.R.; Terra, M.d.C.N.S.; Beskow, S. Water balance of an Atlantic forest remnant under a prolonged drought period. Ciência Agrotecnol. 2021, 45, e008421. [Google Scholar] [CrossRef]
  53. O’Connor, J.; Santos, M.J.; Rebel, K.T.; Dekker, S.C. The influence of water table depth on evapotranspiration in the Amazon arc of deforestation. Hydrol. Earth Syst. Sci. 2019, 23, 3917–3931. [Google Scholar] [CrossRef]
  54. Bruijnzeel, L. Hydrological functions of tropical forests: Not seeing the soil for the trees? Agric. Ecosyst. Environ. 2004, 104, 185–228. [Google Scholar] [CrossRef]
  55. Salemi, L.F.; Groppo, J.D.; Trevisan, R.; de Moraes, J.M.; Ferraz, S.F.d.B.; Villani, J.P.; Duarte-Neto, P.J.; Martinelli, L.A. Land-use change in the Atlantic rainforest region: Consequences for the hydrology of small catchments. J. Hydrol. 2013, 499, 100–109. [Google Scholar] [CrossRef]
  56. Centeno, L.N.; Hu, W.; Timm, L.C.; She, D.; Ferreira, A.d.S.; Barros, W.S.; Beskow, S.; Caldeira, T.L. Dominant Control of Macroporosity on Saturated Soil Hydraulic Conductivity at Multiple Scales and Locations Revealed by Wavelet Analyses. J. Soil Sci. Plant Nutr. 2020, 20, 1686–1702. [Google Scholar] [CrossRef]
  57. Hakamada, R.E.; Hubbard, R.M.; Stape, J.L.; Lima, W.d.P.; Moreira, G.G.; Ferraz, S.F.d.B. Stocking effects on seasonal tree transpiration and ecosystem water balance in a fast-growing Eucalyptus plantation in Brazil. For. Ecol. Manag. 2020, 466, 118149. [Google Scholar] [CrossRef]
Figure 1. Location of the Cantareira System and study sites, soil types, and land uses. The gray color represents Brazil.
Figure 1. Location of the Cantareira System and study sites, soil types, and land uses. The gray color represents Brazil.
Sustainability 15 13395 g001
Figure 2. Weekly rainfall during the studied period in (a) Bragança Paulista (monthly average rainfall from 2004 to 2021); and (b) Monte Verde (monthly average rainfall from 2017 to 2021). The black arrows indicate the timing of the electrical resistivity tomography acquisitions.
Figure 2. Weekly rainfall during the studied period in (a) Bragança Paulista (monthly average rainfall from 2004 to 2021); and (b) Monte Verde (monthly average rainfall from 2017 to 2021). The black arrows indicate the timing of the electrical resistivity tomography acquisitions.
Sustainability 15 13395 g002
Figure 3. Surface runoff determined by Cornell Sprinkle infiltrometer for Typic Hapludult (a), Typic Dystrudept (b) and, Typic Usthortent (c), under native forest (NF), eucalyptus (E), rotational grazing (RG), and extensive grazing (EG).
Figure 3. Surface runoff determined by Cornell Sprinkle infiltrometer for Typic Hapludult (a), Typic Dystrudept (b) and, Typic Usthortent (c), under native forest (NF), eucalyptus (E), rotational grazing (RG), and extensive grazing (EG).
Sustainability 15 13395 g003
Figure 4. ρ–θ modeling calibration curves (lines) for soil type: Typic Hapludult, Typic Dystrudept, and Typic Usthortent and its observed data.
Figure 4. ρ–θ modeling calibration curves (lines) for soil type: Typic Hapludult, Typic Dystrudept, and Typic Usthortent and its observed data.
Sustainability 15 13395 g004
Figure 5. Estimated soil water content (m3 m−3) for the Typic Hapludult under native forest (NF), eucalyptus (E), rotational grazing (RG), and extensive grazing (EG) during the monitoring periods. The absolute error is indicated by the percentages on the lower left side.
Figure 5. Estimated soil water content (m3 m−3) for the Typic Hapludult under native forest (NF), eucalyptus (E), rotational grazing (RG), and extensive grazing (EG) during the monitoring periods. The absolute error is indicated by the percentages on the lower left side.
Sustainability 15 13395 g005
Figure 6. Estimated soil water content (m3 m−3) for the Typic Dystrudept under native forest (NF), eucalyptus (E), rotational grazing (RG), and extensive grazing (EG) during the monitoring periods. The absolute error is indicated by the percentages on the lower left side.
Figure 6. Estimated soil water content (m3 m−3) for the Typic Dystrudept under native forest (NF), eucalyptus (E), rotational grazing (RG), and extensive grazing (EG) during the monitoring periods. The absolute error is indicated by the percentages on the lower left side.
Sustainability 15 13395 g006
Figure 7. Estimated soil water content (m3 m−3) for the Typic Usthortent under native forest (NF), eucalyptus (E), rotational grazing (RG), and extensive grazing (EG) during the monitoring periods. The absolute error is indicated by the percentages on the lower left side.
Figure 7. Estimated soil water content (m3 m−3) for the Typic Usthortent under native forest (NF), eucalyptus (E), rotational grazing (RG), and extensive grazing (EG) during the monitoring periods. The absolute error is indicated by the percentages on the lower left side.
Sustainability 15 13395 g007
Table 1. Soil classification and characteristics in the soil profile in the native forest (NF), eucalyptus (E), rotation grazing (RG), and extensive grazing (EG).
Table 1. Soil classification and characteristics in the soil profile in the native forest (NF), eucalyptus (E), rotation grazing (RG), and extensive grazing (EG).
Soil/CityLand UseDepthClaySiltSandBDKSTPMACMIP
cm-----------------%---------------Mg m−3cm h−1------------m3 m−3------------
Typic Hapludult/NF0–53822401.14117.210.510.150.36
Piracaia/ 30–354518371.5954.530.480.070.41
E0–53318491.266.130.500.170.33
30–355615291.3934.870.480.140.34
RG0–52316611.3214.600.510.150.36
30–354013471.582.010.470.100.37
EG0–52612621.441.800.510.120.39
30–352815571.520.740.450.080.37
Typic Dystrudept/NF0–53615490.88115.900.460.180.29
Nazaré Paulista 30–3541542193915.580.560.160.40
E0–53314531.1936.390.500.190.31
30–354312451.482.970.460.110.35
RG0–54017431.3119.450.480.090.39
30–354614401.408.520.510.140.37
EG0–53419471.452.610.460.090.37
30–354216421.524.180.470.130.35
Typic Usthortent/NF0–55512330.67221.810.620.240.37
Joanópolis 30–352920511.11150.540.560.200.35
E0–53823390.774.060.650.180.46
30–354532231.20234.760.530.190.35
RG0–54316410.8744.470.660.240.41
30–354422341.2482.090.640.290.34
EG0–55513320.967.330.670.190.48
30–352717561.2795.070.520.240.27
BD: soil bulk density; KS: hydraulic conductivity; TP: total porosity; MAC: macroporosity (>0.05 mm); MIP: microporosity.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Santana, M.L.T.; Faria, V.L.d.; Barbosa, S.M.; Serafim, M.E.; Uezu, A.; Silva, B.M.; Avanzi, J.C. Impacts of Land Use Changes on Soil Functions and Water Security: Insights from a Three-Year-Long Study in the Cantareira System, Southeast of Brazil. Sustainability 2023, 15, 13395. https://0-doi-org.brum.beds.ac.uk/10.3390/su151813395

AMA Style

Santana MLT, Faria VLd, Barbosa SM, Serafim ME, Uezu A, Silva BM, Avanzi JC. Impacts of Land Use Changes on Soil Functions and Water Security: Insights from a Three-Year-Long Study in the Cantareira System, Southeast of Brazil. Sustainability. 2023; 15(18):13395. https://0-doi-org.brum.beds.ac.uk/10.3390/su151813395

Chicago/Turabian Style

Santana, Monna Lysa Teixeira, Vanêssa Lopes de Faria, Samara Martins Barbosa, Milson Evaldo Serafim, Alexandre Uezu, Bruno Montoani Silva, and Junior Cesar Avanzi. 2023. "Impacts of Land Use Changes on Soil Functions and Water Security: Insights from a Three-Year-Long Study in the Cantareira System, Southeast of Brazil" Sustainability 15, no. 18: 13395. https://0-doi-org.brum.beds.ac.uk/10.3390/su151813395

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