Next Article in Journal
Transferrable Principles to Revolutionize Drinking Water Governance in First Nation Communities in Canada
Previous Article in Journal
Long-Term Modelling of an Agricultural and Urban River Catchment with SWMM Upgraded by the Evapotranspiration Model UrbanEVA
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Abacus to Predict Groundwater Recharge at Non-Instrumented Hydrographic Basins

1
Department of Forestry Engineering, Federal Institute of the North of Minas Gerais—Campus Salinas, Salinas 39560-000, Brazil
2
Department of Civil and Environmental Engineering, University of Brasília, Brasília 70910-900, Brazil
*
Author to whom correspondence should be addressed.
Water 2020, 12(11), 3090; https://doi.org/10.3390/w12113090
Submission received: 23 September 2020 / Revised: 23 October 2020 / Accepted: 31 October 2020 / Published: 4 November 2020
(This article belongs to the Section Hydrology)

Abstract

:
One of the first steps to implement a policy for groundwater resources management is knowing the groundwater recharge. However, the unavailability of data and resources to execute field studies increase the uncertainty associated with the estimation of groundwater recharge. To fill this gap, the present work aimed to propose a method to predict groundwater recharge at non-instrumented hydrographic basins. The approach proposed is based on using an abacus to execute the transposition and/or regionalization of results generated in an experimental basin. The methodology comprised the estimation and mapping of recharge rates in the experimental basin using three distinct approaches—numerical modelling of the saturated zone, distributed hydrological modelling of the vadose zone, and the method of fluctuation of the water table elevation—and the following generation of the abacus, with average recharge values for combinations of soil class, land use/cover and slope using geographic information systems. The results indicate that the abacus is consistent for some Ferrasol areas, that the reliability of average regionalized values depends on the complexity of the physical environment—soil class, land use/cover, and slope—and that new studies, focusing on the hydro-physical characterization of soils, might produce more reliable estimations.

1. Introduction

Despite the importance of groundwater in volume, ecology, economy, and society, the lack of knowledge about hydrogeological systems contrasts with scenarios of contamination and depletion of these systems because of inordinate exploitation [1,2,3]. Considering that the remediation and/or recovery of aquifers is, at times, unfeasible [3], management is the best alternative to protect and maintain groundwater reserves [3,4,5]. One of the first steps to implement a groundwater management policy is the knowledge of water availability and for that it is essential to investigate aquifer’s renewal rate [6,7,8].
Rainfall is the main source for groundwater natural recharge, which is regulated by an array of physical and climatic factors, such as geological framework, soil type, slope, vegetation, land cover, and evapotranspiration [9,10,11]. According to Scanlon et al. [12] and Singh et al. [13], there is no universal approach to study groundwater recharge and that is why the amount of methods reported in scientific literature is high. Amongst these methods, it can be highlighted: (i) the ones based on measuring the water table level [14,15,16]; (ii) the ones that use chemical tracers [17,18,19,20,21]; (iii) the ones based on direct or indirect measures of hydrological flows, such as baseflow [15,22,23,24] and water balance [15,25,26,27]; and (iv) the ones based on hydrological or numerical modelling and/or simulations [27,28,29,30,31,32,33].
Quantifying groundwater recharge, however, is a difficult task [13,34,35]. Such difficulty is based on the fact that each method, among those available, may have applicability restricted to certain physical conditions of the study area, such as types of soil or climate, and it can valuate a specific type of recharge, such as potential recharge or effective recharge, or it can be indicated only for punctual estimates [15]. Some methods that use direct flow measures, for example, may generate good results when obtaining point values, but may be infeasible to characterize the spatial variability of the process. On the other hand, distribution methods are disregarded, at times, due to the amount and variety of required data, a demand that is limited precisely by the lack of such data [7,15]. This complexity represents a considerable challenge to quantify the recharge and the only solution to obtain reliable estimates, in many cases, is the combined application of several methods which in turn, would require resources and data not always available [15]. The problem is not restricted to groundwater recharge studies, and a way to overcome it is the transposition or regionalization of the results generated at experimental basins [36]. In such a case, the hydrological monitoring is intensified by adopting a smaller study area, and the selected environment may have physical and environmental characteristics that represent the larger basin or the biome to which the experimental basin belongs. Consequently, the data collected may be used as the foundation for studies that can generate results applicable also for non-instrumented areas with similar characteristics [37]. Lima et al. [38] and Creed et al. [39] used the solution when they adopted an experimental basin to study and comprehend the general hydrological and climatological processes, such as evapotranspiration, surface flow, generation of sediments, and global warming. In groundwater recharge research, Chang and Rubin [36] evaluated the performance of the statistical hierarchical analysis to understand the relationship between hydrological response watershed characteristics, for the regionalization of recharge in ungauged basins and Ibrahim et al. [40] used multiple linear regression to predict recharge in ungauged basins, generated from the relationship between simulated recharge in five gauged basins and their physical and climatic characteristics.
Considering that the groundwater recharge-regionalization remains a non-trivial procedure [36], transferring the results generated at experimental or representative basins still requires scientific efforts, mainly when observing estimations required for undocumented areas where the values adopted are arbitrary or based on studies executed at other areas with completely different characteristics than the study area. In Santos and Koide [41], for example, the simulated recharge rates are higher than 70% of the rainfall in a Cerrado area, while the management institutions and modelling studies in Brazil frequently adopt the average value of 20%, according to Araujo [42]. Recent works have shown the extent of the error that such general assumptions induce, especially when the objective is land use/cover management for water management, because even within the same area, biome, or hydrographic basin, the recharge rate can have wide spatial variations [43].In Bortolin et al. [1], the estimated recharges at a hydrographic basin located in the state of Rio Grande do Sul varied between 9% and 36% related to annual rainfall. The results of Ramires and Manzione [44] indicate a recharge between 20% and 36% at a Cerrado area in the center west of the state of São Paulo. For the same region, Gonçalves and Manzione [7] suggest that the percentage recharge rates do not suffer significative variations in function of the rainfall volume. Silva [45] estimated values with a higher range of variations, between 15% and 50%, also for a Cerrado biome area in São Paulo. In the state of Minas Gerais, Souza et al. [2] simulated and mapped the spatial distribution of recharge rates at the Doce river basin and, for the Ferrasol class, reached results very close to the ones estimated by Cambraia-Neto and Rodrigues [8], close to 30% of the annual rainfall for the same soil type at a basin in the Federal District.
The scientific efforts to comprehend the recharge process would be fruitless if not for the conversion of knowledge in usable data, without which the management stakeholders would hardly fulfil their role regarding the guarantee of universal and sustainable water availability. In this context, scientific works research the relation between urbanization and water quality, between land cover, flow, and sediments, and between vegetation and evapotranspiration, amongst others. Regardless of the issue, the intention, for the most part, is to provide tools for information gathering, in a way that water management may be more rational and efficient. According to Bortolin et al. [1], for surface water a lot has been accomplished in that sense, but there is still a huge gap to bridge when it comes to groundwater.
Charles et al. [46] estimated recharge rates using water balance for several combinations of climate, soil type and land cover. They proposed an electronic spreadsheet to estimate recharge in other sites using the mentioned environmental factors as input data. Although the study with the proposed spreadsheet did not have an experimental basin, it is a case of results transposition, which is useful for the process of integrated water resources management, but hardly used in Brazil. Santos [47] presented, preliminarily, an alternative, however, at the time, there were no published studies that could validate the proposition.
In that sense, the aim of the present work is to propose an alternative method to predict groundwater recharge at non-instrumented hydrographic basins based on the transposition and/or regionalization of the results generated at an experimental basin representative of the Cerrado biome.

2. Materials and Methods

The study area comprised, as the experimental basin, the Capão Comprido basin, located at the west portion of the Federal District, between the west longitude meridians 48°10′07″ and 48°06′13″, and the south latitude parallels 15°43′42″ and 15°45′41″ covering an area of approximately 16 km2. The Ribeirão Rodeador basin, with 113 km2, was used to validate the study (Figure 1).
The two basins are located in the Brazilian Cerrado biome and are under the same climatic domain—Aw humid/sub-humid—with well-defined rainy (spring/summer) and dry (autumn/winter) seasons. Rainfall occurs mainly between October and April, and the average temperatures of the hottest month and the coldest month are 24 °C and 10 °C, respectively. The air relative humidity can reach a critical level of 13%. The annual average rainfall between 1971 and 2016 was 1500 mm [42,47]. Geologically, the same formations of the Paranoá Group are present in both basins—Clay Metarritmito, Medium Quartzites and Sandy Metarritmito [48]. According to Reatto et al. [49], in the Capão Comprido basin the predominant types of soil are Red and Red-yellow Ferrasol, presenting also, in fewer areas, Cambisol, Plinthosol, and Gleysol. In the Ribeirão Rodeador basin, the same soil classes occur, with, however, more occurrence of Cambisols in relation to the Capão Comprido basin, and a small occurrence of Quartzarenic Neosol [42]. Regarding land use/cover, the Capão Comprido basin is predominantly rural, occupied by small farms, with remnants of riparian forests, typical formations of the Cerrado Biome—such as open cerrado grassland, shrub cerrado, and cerrado fields and cerrado—areas under anthropic uses, such as pastures, reforestation, olericulture and perenial/fruit crops, and bare-soil areas. In the Rodeador basin, the majority land use/cover are Annual Crops/olericulture, cerrado fields, open cerrado grassland and reforestation, also occurring small areas of pastures, perenial/fruit crops, riparian forest, bare-soil, and urbanization. Eight land use cover classes are common to both basins: Riparian Forest; Forest/Reforestation; Perenial Tree Culture, Cerrado Fields, Shrub Cerrado, Open Cerrado Grassland, Pasture and Bare Soil [42,47]. The main difference between the two basins, in terms of land use/cover, refers to the Annual Cultures/Olericulture class, with wide distribution in the Rodeador basin, but is not very significant in the Capão Comprido basin.
The materials used in the present work consisted of field data and theme plans of cartography information. The field data, continuously monitored at the study area, comprised historical series of pluviometry, fluviometry, meteorology, piezometry, and hydrogeological information, while the geographic data is constituted of maps of geology [48], soils [49], hydrogeology [48], altimetry from the SICAD-Federal District Cartography System and satellite imagery from the HRV/SPOT5 sensor. The elaboration of models and simulation were executed with the software Visual MODFLOW v.4.2, WETSPA (Water and Energy Transfer between Soils, Plants, and Atmosphere) and IDRISI, Andes version, while the maps were prepared in the software QGIS Desktop 2.18.23.
The methodology comprehended (i) mapping the spatial distribution of the recharge rates using three different approaches, (ii) generation of recharge codes and average recharge values for each possible combination of soil, land use/cover and slope, (iii) and verification of average abacus recharge rates, when regionalized in another area under the same characteristics. The steps are detailed in the following sections.

2.1. Mapping the Spatial Distribution of the Recharge Rates

Considering the lack of direct methods to quantify groundwater recharge and the specificities of the available approaches, the best way to study recharge is to apply at least two different methods, preferably ones that have no correlation [12]. Taking into account the aims of the present work, distributed methods were selected and applied. These methods are capable of generating not only estimations, but also spatially characterizing the recharge process. In that sense, three approaches were proposed: numerical modelling of the flow in saturated porous medium (MODFLOW); hydrological modelling of the vadose zone (WETSPA); and spatialization of point estimations calculated using the water table elevation (WTE).

2.1.1. Numerical Modelling of the Saturated Zone

This approach comprehended the elaboration of a numerical hydrogeological model of the phreatic aquifer at the study area, which was used to simulate groundwater flows and recharge via vertical and horizontal numerical simulation, using finite difference cells at the Visual MODFLOW v.4.2 software. The modelling contained the following reasoning: (a) recharge is one of the processes responsible for the alteration of the condition of the aquifer system over time and there are no measured values for the recharge processes at the study area; (b) however, the “fixed” characteristics of the system (structural configuration and properties of the aquifer’s material) and the temporal variations of its behavior (storage and phreatic level) are known, as well as other correlated measurable factors (evapotranspiration and surface water overflow); (c) thus, recharge is the unknown variable to be found. Since recharge is usually an input data, the problem can be solved with reverse modelling. Therefore, spatial distribution and recharge rates were estimated using calibration. The methodology developed, summarized in Figure 2, is more detailed in Santos [47] and Santos and Koide [41].

2.1.2. Distributed Hydrological Modelling of the Vadose Zone

This step was based on applying the model WETSPA—Water and Energy Transfer between Soils, Plants, and Atmosphere [50]. It is a physically based distributed model with routines to simulate water flow between the atmosphere, vegetation’s canopy, root zone, transmission zone, and saturation zone, as illustrated in Figure 3.
Using basic maps of slope, soil types, and land cover, this model executes the mapping of parameters spatial distribution in geographic information systems. The parameters modelled are related to water storage and movement through physical compartments, such as field capacity, permanent wilting point, hydraulic conductivity, surface water flow coefficients, leaf area index and Manning coefficient, amongst others. The mapping is based on linking spatial information in the basic maps with tables with parameters values, indicated in function of constant values or classes in the maps. Then, using the spatially distributed parameters, the measured data of average rainfall and the estimated values of potential evapotranspiration, the algorithm simulates the hydrological processes, including recharge, for each cell in the model. The daily total surface flow at the basin discharge point, estimated using the discharge rating curve, was used as a control observation in the calibration and verification steps. The complete formulation of the model is detailed in Santos [47] and Liu and Smedt [50].

2.1.3. Spatialization of Point Estimations Calculated Using Water Table Level Elevation

Having the monitoring data of the phreatic level, the piezometric elevations resulting from the recharge period during the rainy season were estimated. The elevations observed in each monitoring well were obtained graphically, measuring vertical distance between the maximum piezometric head and the prolongation of the recession line of the piezometric curve, a method indicated by Healy and Cook [51]. The point recharge, valid for the monitoring well’s influence area where the phreatic elevation was measured, was estimated by multiplying the total elevation, in mm/year, by the aquifer’s local value for the specific productivity coefficient, Sy. The values for Sy were obtained using a pumping test and a slug test, performed by Santos [47], in the monitoring wells at the Capão Comprido basin.
The spatial distribution of the recharge was obtained using a model for spatial multiple regression, of the type “Y = a + b1X1 + b2X2 + b3X3 + ⋯ + bnXn”, placing the point estimation in the dependent variable (Y) in relation to the spatial distribution of the predefined factors (X1, X2, X3, …Xn) that are determining to the process (terrain slope, land cover and use, hydrodynamic behavior of the aquifer, aquifer width, and phreatic depth, in this work). Since the order of magnitude and the unities of the mentioned factors are different, the respective normalized maps were used. More details about the method proposed in this step is described in Santos [47].

2.2. Generation of Recharge Codes and Average Recharge Values for Possible Combinations between Soil, Land Use/Cover, and Slope

Based on the principle which states that groundwater recharge is a result of the combination of “fixed” environmental factors, it can be expected that areas in the same biome or geographic region will have predictable hydrological behavior regarding the conversion of rainfall in groundwater at the aquifer. Identifying the possible combinations of such factors and quantifying the recharge for the mentioned conditions in representative sites or basins, it is supposed to be possible to estimate the unknown variable for non-instrumented areas that have the same combination of considered environmental factors. In this sense, the last step comprised of mapping the spatial distribution of the combinations of factors related to recharge, and the intersection of this map with the spatial distribution of the recharge rates generated in previous steps of the present work.
In the first procedure, based on the conceptual proposed by Santos [47], were considered the factors “soil class”, “land use/cover”, and “slope”. Unique identification numbers were attributed to the classes or variations from each factor in the following way:
  • soil class: 1 = cambisol; 2 = red-yellow ferrasol; 3 = gleysol; 4 = red ferrasol; 5 = plinthosol;
  • land use/cover: 1 = riparian forest; 2 = forest/reforestation; 3 = perennial tree culture; 4 = cerrado fields; 5 = shrub cerrado; 6 = open cerrado grassland; 7 = pasture; 8 = baresoil;
  • slope: 1 = slope between 0 and 5%; 2 = slope between 5 and 10%; 3 = slope between 10 and 15%; 4 = slope over 15%;
The maps related to each factor were converted to vector format and the data crossing was performed by an intersection operation, both spatial and of attribute, resulting in polygons with unique three digits identification numbers representing the physical conditions of the site. Lastly, it was executed by a spatial intersection, using a zone operation between the polygons with the combination of factors and the spatial distribution of recharges, thus obtaining an abacus, with which it is possible to estimate, for each combination of type of soil, land cover and slope, the average values of simulated recharge, considering the three aforementioned methods.

2.3. Verification of Average Abacus Recharge Rates

The validation of the average estimation of recharge in function of the combination of physical and environmental factors was performed using the abacus at another area in the Cerrado biome. Using the maps of type of soil, land cover/use, and slope of the validation basin, new maps were generated with the combination and the respective codes. Lastly, applying the abacus to the map and recharge codes of the validation basin, transposed/regionalized values for the recharge were obtained, which were confronted with the estimations generated by Araujo [42]. In this research, the spatial distribution of the recharge rates was obtained through integrated hydrological modeling of surface and saturated zone processes. The model used was the SWAT-MODFLOW [52], a tool that integrates the models SWAT (Soil and Water Assessment Tool) and MODFLOW. The model input was the time series of meteorological, precipitation, discharge and water table data, physical digital cartographic information of the basin (soil map, digital elevation model and land use/cover) and geological and hydrodynamic information for the saturated porous medium (spatial distribution of saturated hydraulic conductivity and specific yield and vertical stratification of the model domain). The discharge, base flow and water table level data refer to the hydrological year 2008/2009—the same time period adopted in the generation of the abacus—were used for model calibration and verification. The recharge simulated by this method is “potential”, as the abacus recharges obtained by the WETSPA model. The other results of the abacus—MN and WTF, of effective recharge—were evaluated by comparison with average values of effective recharge found in the literature.

3. Results and Discussion

The present work mapped 85 possible combinations of type of soil, land use and slope at the Ribeirão Rodeador basin. The average values for the simulated recharge in Araujo [42] are presented on the map in Figure 4a. Figure 4b presents the regionalized values from the abacus application. For some areas in the Rodeador basin it was not possible to map and compare the values due to the existence of combination between factors that are not present in the Capão Comprido basin, used to formulate the abacus. The recharge rates estimated by the methods “Numerical Modelling—NM” and “Water Table Elevation—WTE” were not included in the comparison, because they perform “effective recharge” estimations [8,41] and, consequently, give lower estimations than the ones obtained with the vadose zone methods, such as SWAT-MODFLOW and WETSPA.
According to the results simulated by Araujo [42], the average recharge rates varied between 20% and 60% of the total annual rainfall (316.2 and 948.6 mm/yr, based on the rainfall of the hydrological year 2008/2009, of 1581 mm), with higher values associated with types of soil Cambisol, Red-yellow Ferrasol, Plinthosol, and Gleysol. It was also observed the influence of the slope on the groundwater recharge process, that is, the recharge rates are low when the slope is high, in a same type of soil. For the regionalized values, the recharge rates were lower, in general, and the rates increased, presenting the higher values, when associated with deep soils, such as Ferrasols and Plinthosols, in low slope areas. This behavior in function of the individual variables of the factors “soil class”, “land use/cover”, and “slope” is synthetized on Figure 5.
Despite the difficulty of analyzing separately the influence of each factor [36], it was observed that the recharge varies according to the type of soil—between 35% and 48% in the simulated average rates per class and between 19% and 38% in the regionalized rates averages per class. However, depending on the class, land use/cover may have a major or minor influence than the soil, as, for example, in the case of Cambisols and Red Ferrasols, respectively. In the former, the type of land use/cover caused average rates between 30% and 60%, while in the latter, the variation observed between rates was lower, with values between 32% and 38%, for the simulated rates. Considering the differences between the approaches, a greater influence of the type of soil is observed, followed by the declivity and land use, which can be better visualized by mapping the relative error and spatially comparing it with each factor analyzed (Figure 6), and by observing the results from Table 1: a preliminary abacus, through which recharge rates can be regionalized in non-instrumented hydrographic basins for combinations of soil class, land use/cover, and slope factors. Estimations with low relative error can be considered more reliable [15,35].
The average relative error was −27.5, with the highest values occurring in steeper areas with Cambisol, and some areas with Red-yellow Ferrasol and Plinthosol covered with trees along the drainage network. In those areas, that correspond to circa 29% of the basin area, the relative errors modules of the regionalized rates were higher than 25%—a limit reference value of performance for the statistical index PBIAS [53]—reaching 96% (this extreme error value occurred in an area of about 3 hectares; less than 1% of the total area of the validation basin). It can be highlighted that the mentioned combinations of physical factors are spatially correlated and correspond to configurations of high physical complexity, to which the hydrological models are still limited. In the other areas, the average relative error was −19%.
The model SWAT is a respectable reference when it comes to hydrological modelling. Such status is due to the good performance of the model, mostly to simulate flow rates, exemplified by the countless works developed over the last 20 years in a variety of environments [32]. Nonetheless, Araujo [42] and Arroio Junior [54] emphasize that the model still presents some limitations that may compromise the consistency of the water balance, even if the simulated flow rates are well adhered to the measured flow rates. A good example of these limitations is the lack of reliable parameters for forested steep lands, which frequently results in inconsistencies in the simulation of subsurface flow [54], and another example is the underestimation of the evapotranspiration, that may lead to the overestimation of surface flow rates or recharge rates.
Analyzing the hydrological behavior of five hydrographic basins, Arroio Junior [54] faced the issue of identifying variations amongst the simulated actual evapotranspiration and the measured reference values. After modifying the routine to calculate evapotranspiration, there was a significant reduction in recharge in the five basins.
Despite the discussion about the limitation of the model SWAT-MODFLOW, the results obtained by Araujo [42] indicate that the model’s recharge estimations are more reliable. In this model, the simulation of surface and groundwater flows are integrated, exploring the strengths from each original model—SWAT simulating surface water flow and MODFLOW simulating groundwater flow—thus making the global modelling more robust to represent every flow in the hydrographic basin [55]. Using the integrated model, Araujo [42] compared simulated flow rates with measured flow rates, obtaining satisfactory results for the performance indexes, which indicate SWAT-MODFLOW model’s ability to represent the hydrological processes in a Cerrado environment. However, the model tends to overestimate baseflow, which could indicate the overestimation of the recharge. Notably, that could happen in steeper areas with Cambisol, where the rates simulated by SWAT-MODFLOW were higher than in other types of soil in the basin and overcame in more than 100% the regionalized rates generated by WETSPA model.
On the other hand, Santos [47] demonstrates the overestimation of recharges due to the underestimation of baseflows in the period used to generate the reference values. Even if comparative studies, such as Nossent et al. [56], reveal an equality in both model’s performance when applied to the simulation of hydrological processes in hydrographic basins, WETSPA model’s performance was not satisfactory on the case study that generated the reference values.
Then, which recharge rates are indicated? According to Polanco et al. [57], a common limitation in hydrological modelling is the calibration and assessment of results based only in outflow rates. In such cases, a satisfactory statistical performance does not guarantee the consistency of values simulated for water balance nor hydrological processes at the local scale. In that sense, it is cautiously advised to adopt intermediary recharge rate values, between the overestimated ones generated by SWAT-MODFLOW and the underestimated regionalized ones.
For the Cambisol class (number “1”, first digit), the great variation between simulated and regionalized values suggest caution when adopting a mean reference value since, according to Castro [58], this type of soil has a wide range of hydrological behavior, which would cause the difference to be actually from real physical variations in the basins. Nonetheless, it can be highlighted that Cambisol in the Ribeirão Rodeador basin is located, predominantly, in steeper areas with Cerrado fields and riparian forests, a combination that, according to Bortolin et al. [1] and Santos [47], is theoretically unfavorable to recharge groundwater.
The estimations obtained for the Red-yellow Ferrasol class (number “4”, first digit) were considered more reliable because, in general, they presented the lowest differences. Furthermore, the spatialized values simulated with SWAT-MODFLOW [42] and regionalized with WETSPA in the abacus, which varied between 32% and 39% and between 32% and 36%, respectively, were situated close to the ones obtained by Cambraia-Neto and Rodrigues [8], of 31%, related to annual rainfall for the same physical conditions of soil type, land use and slope in another hydrographic basin located in the east portion of the Federal District. Souza et al. [2] also attained percentage recharge rates of the same order of magnitude, of 36%, for Red-yellow Ferrasol in the Rio Doce basin, west of the state of Minas Gerais. However, in this case the rates of the same order of magnitude for Ferrasols were obtained in different environmental conditions (biome and climate), indicating a dominance of the “soil class” factor in the recharge process and that the role of this factor may be a good variable to use in recharge rate regionalization.
Regarding the effective recharge, the average rates for combinations of soil, land use/cover, and slope have been converted into weighted average rates, by the area factor, for the basin. The average values suggested by the abacus and the reference values obtained in the literature are shown in Table 2.
The average effective recharges suggested by the abacus are smaller than the potential simulated and regionalized recharges, corroborating Gonçalves and Manzione [7] and Healy and Cook [51] regarding the expected difference between the two types of recharge. In comparison to the reference values, the suggested effective rates were considered consistent, in the same order of magnitude as the general average of the rates reported by the mentioned works. It is important to highlight that the study area adopted by Cambraia-Neto and Rodrigues [8]—eastern region of the Federal District—is relatively close to the abacus generating basin, and that the effective average rate of Araújo [42]—19%, is equal to the average of the abacus, also obtained by numerical modeling—it was obtained in the validation basin. The high upper limit of the recharge rates obtained in Oliveira et al. [62] is because the soil in your study area is a highly permeable sandy Quartzarenic Neosol. Thus, the regionalization of the recharge rates proposed by the abacus would only be possible if the study basin has at least the three predominant classes of soil in the Capão Comprido basin—Red Ferrasol, Red-yellow Ferrasol and Cambisols—as occurs in Gaspar [59], Albuquerque and Chaves [60], Cunha [61], ANA [64], Gonçalves and Manzione [7], and Cambraia-Neto and Rodrigues [8].

4. Conclusions

The present work proposed an abacus to predict groundwater recharge at non-instrumented hydrographic basins elaborated from the results obtained at an experimental basin in the Cerrado biome. Confronting the proposed values to the ones available in scientific literature and estimated by Araujo [42], for a hydrographic basin with the same environmental conditions, it can be concluded that:
  • the potential recharge varies according to the soil class. However, in the same type of soil, land use/cover and slope can cause large variations in recharge rates;
  • for combinations between Cambisols and high slopes, the differences between simulated and regionalized potential recharge rates are high, indicating, in this case, caution on the adoption of values suggested by the abacus. As the calibration and verification of hydrological models are proceeds at the basin or sub-basin level, it is necessary to perform site specific studies to confirm the recharge rate actual order of magnitude;
  • for soils with more homogeneous porous matrix, such as Ferrasol, the difference between simulated and regionalized rates was smaller, indicating greater reliability on the values suggested by the abacus for such situations;
  • the regionalized values of potential recharge suggested by the abacus for combinations with Ferrasols were considered consistent, with order of magnitude confirmed by Cambraia-Neto and Rodrigues [8] and Araujo [42]. Thus, for these combinations, the suggested values can be transferred to other areas or basins in the Cerrado biome with similar characteristics. Similar values estimated by Souza et al. [2] for this type of soil, but with different climate and biome, suggest that the abacus values for Ferrasols may also be valid for other Brazilian biomes;
  • the divergences between simulated and regionalized estimations observed in areas with the same soil type, land cover and slope are a result of the conceptual and mathematical formulations specific in the models used, mostly regarding the calculation of actual evapotranspiration, which confirms Araujo [42] and Arroio Junior [54];
  • the average effective rates of recharge—NM and WTE—suggested by the abacus were considered consistent, with values in the same order of magnitude as the reference values in the literature. However, the regionalization of these recharge rates would only be possible if the study basin has at least the three predominant classes of soil in the Capão Comprido basin—conditions in which the abacus was generated;
  • new studies could produce a more robust abacus. In this sense, it is recommended to focus on the hydro-physical characterization of the soil types, which would enhance the range and reliability of future regionalization of proposed values;
  • for non-instrumented basins with physical characteristics similar to the study area, also used to validate the results, the abacus could be an useful tool for surface and groundwater resources integrated management, since the recharge rates can be obtained in an easier way, just by knowing annual average rainfall, land use, soil type, and slope of the new study area, which are data relatively easy to obtain.

Author Contributions

Conceptualization, R.M.d.S., S.K.; methodology, R.M.d.S., S.K.; validation, B.E.T. and D.L.d.A.; resources, S.K.; writing—original draft preparation, R.M.d.S.; writing—review and editing, R.M.d.S., S.K., B.E.T. and D.L.d.A.; supervision, S.K.; funding acquisition, S.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by ConselhoNacional de DesenvolvimentoCientífico e Tecnológico (CNPq), Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) and Fundação de Apoio à Pesquisa do Distrito Federal (FAP-DF), for the PhD scholarship and further financial support granted to execute the work.

Acknowledgments

The authors thank the public worker Irioman Junior for the support on field activities andEmpresaBrasileira de PesquisaAgropecuária (EMBRAPA)—Cerrados for providing the meteorological data.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bortolin, T.A.; Reginato, P.A.R.; Presotto, M.A.; Schneider, V.E. Estimativas de recarga aquífera com uso de filtros digitais em sub-bacias hidrográficas do Sistema Aquífero Serra Geral no estado do Rio Grande do Sul. Sci. Cum Ind. 2018, 6, 21–30. [Google Scholar] [CrossRef] [Green Version]
  2. Souza, E.; Pontes, L.M.; Fernandes Filho, E.I.; Schaefer, C.E.G.R.; Santos, E.E. Spatial and temporal potential groundwater recharge: The case of the doce river basin, Brazil. Rev. Bras. Cienc. Solo 2019, 43, 1–27. [Google Scholar] [CrossRef] [Green Version]
  3. Rahman, A.T.M.S.; Hosono, T.; Mazumnder, Q.H.; Jahan, C.S. Sustainable Groundwater Management in Context of Climate Change in Northwest Bangladesh. In Achievements and Challenges of Integrated River Basin Management; IntechOpen: London, UK, 2018; pp. 101–120. [Google Scholar]
  4. Alfranca, O. Groundwater Management Competitive Solutions: The Relevance of the Gisser-Sanchez Model. In Groundwater—Resource Characterisation and Management Aspects; IntechOpen: London, UK, 2019; pp. 1–13. [Google Scholar]
  5. Megdal, S.B. Invisible water: The importance of good groundwater governance and management. Npj Clean Water 2018, 1, 15. [Google Scholar] [CrossRef] [Green Version]
  6. Ministério do Meio Ambiente—MMA. Plano Nacional de Recursos Hídricos: Panorama e Estado dos Recursos Hídricos do Brasil; Ministério do Meio Ambiente, Secretaria de Recursos Hídricos: Brasília, Brazil, 2006; ISBN 8577380092.
  7. Gonçalves, V.F.M.; Manzione, R.L. Estimativa da recarga das águas subterrâneas no Sistema Aquífero Bauru (SAB). Geo UERJ 2019, 35, 2–19. [Google Scholar] [CrossRef]
  8. Cambraia-Neto, A.J.; Rodrigues, L.N. Evaluation of groundwater recharge estimation methods in a watershed in the Brazilian Savannah. Environ. Earth Sci. 2020, 79, 140. [Google Scholar] [CrossRef]
  9. Carrera-Hernández, J.J.; Smerdon, B.D.; Mendoza, C.A. Estimating groundwater recharge through unsaturated flow modelling: Sensitivity to boundary conditions and vertical discretization. J. Hydrol. 2012, 452, 90–101. [Google Scholar] [CrossRef]
  10. Mohan, C.; Western, A.W.; Wei, Y.; Saft, M. Predicting groundwater recharge for varying land cover and climate conditions—A global meta-study. Hydrol. Earth Syst. Sci. 2018, 22, 2689–2703. [Google Scholar] [CrossRef] [Green Version]
  11. Kubicz, J.; Kajewski, I.; Kajewska-Szkudlarek, J.; Dąbek, P.B. Groundwater recharge assessment in dry years. Environ. Earth Sci. 2019, 78, 555. [Google Scholar] [CrossRef] [Green Version]
  12. Scanlon, B.R.; Healy, R.W.; Cook, P.G. Choosing appropriate techniques for quantifying groundwater recharge. Hydrogeol. J. 2002, 10, 18–39. [Google Scholar] [CrossRef]
  13. Singh, A.; Panda, S.N.; Uzokwe, V.N.E.; Krause, P. An assessment of groundwater recharge estimation techniques for sustainable resource management. Groundw. Sustain. Dev. 2019, 9. [Google Scholar] [CrossRef]
  14. Mattos, T.S.; de Oliveira, P.T.S.; Lucas, M.C.; Wendland, E. Groundwater recharge decrease replacing pasture by Eucalyptus plantation. Water 2019, 11, 1213. [Google Scholar] [CrossRef] [Green Version]
  15. Walker, D.; Parkin, G.; Schmitter, P.; Gowing, J.; Tilahun, S.A.; Haile, A.T.; Yimam, A.Y. Insights from a multi-method recharge estimation comparison study. Groundwater 2019, 57, 245–258. [Google Scholar] [CrossRef] [Green Version]
  16. Carvalho, F.; Scopel, I. Escoamento superficial e recarga d’água subterrânea em diferentes usos do solo na microbacia do córrego do queixada. Rev. Caminhos Geogr. 2018, 19, 133–145. [Google Scholar] [CrossRef]
  17. Joshi, S.K.; Rai, S.P.; Sinha, R.; Gupta, S.; Densmore, A.L.; Rawat, Y.S.; Shekhar, S. Tracing groundwater recharge sources in the northwestern Indian alluvial aquifer using water isotopes (δ18O, δ2H and 3H). J. Hydrol. 2018, 559, 835–847. [Google Scholar] [CrossRef]
  18. Richards, L.A.; Magnone, D.; Boyce, A.J.; Casanueva-Marenco, M.J.; van Dongen, B.E.; Ballentine, C.J.; Polya, D.A. Delineating sources of groundwater recharge in an arsenic-affected Holocene aquifer in Cambodia using stable isotope-based mixing models. J. Hydrol. 2018, 557, 321–334. [Google Scholar] [CrossRef] [Green Version]
  19. Li, Z.; Jasechko, S.; Si, B. Uncertainties in tritium mass balance models for groundwater recharge estimation. J. Hydrol. 2019, 571, 150–158. [Google Scholar] [CrossRef]
  20. Parlov, J.; Kovač, Z.; Nakić, Z.; Barešić, J. Using water stable isotopes for identifying groundwater recharge sources of the unconfined alluvial Zagreb aquifer (Croatia). Water 2019, 11, 2177. [Google Scholar] [CrossRef] [Green Version]
  21. Ahmed, I.M.; Jalludin, M.; Razack, M. Hydrochemical and isotopic assessment of groundwater in the Goda mountains range system. Republic of Djibouti (horn of Africa). Water 2020, 12, 2004. [Google Scholar] [CrossRef]
  22. Tenenwurcel, M.A.; de Moura, M.S.; da Costa, A.M.; Mota, P.K.; Viana, J.H.M.; Fernandes, L.F.S.; Pacheco, F.A.L. An improved model for the evaluation of groundwater recharge based on the concept of conservative use potential: A study in the river Pandeiros Watershed, Minas Gerais, Brazil. Water 2020, 12, 1001. [Google Scholar] [CrossRef] [Green Version]
  23. Lacombe, G.; Douangsavanh, S.; Vongphachanh, S.; Pavelic, P. Regional assessment of groundwater recharge in the lower mekong basin. Hydrology 2017, 4, 60. [Google Scholar] [CrossRef] [Green Version]
  24. Yang, W.; Xiao, C.; Liang, X. Extraction method of baseflow recession segments based on second-order derivative of streamflow and comparison with four conventional methods. Water 2020, 12, 1953. [Google Scholar] [CrossRef]
  25. Rukundo, E.; Doğan, A. Dominant influencing factors of groundwater recharge spatial patterns in Ergene river catchment, Turkey. Water 2019, 11, 653. [Google Scholar] [CrossRef] [Green Version]
  26. Neukum, C.; Azzam, R. Impact of climate change on groundwater recharge in a small catchment in the Black Forest, Germany. Hydrogeol. J. 2012, 20, 547–560. [Google Scholar] [CrossRef]
  27. Wiebe, A.J.; Rudolph, D.L. On the sensitivity of modelled groundwater recharge estimates to rain gauge network scale. J. Hydrol. 2020, 585, 124741. [Google Scholar] [CrossRef]
  28. Hund, S.V.; Allen, D.M.; Morillas, L.; Johnson, M.S. Groundwater recharge indicator as tool for decision makers to increase socio-hydrological resilience to seasonal drought. J. Hydrol. 2018, 563, 1119–1134. [Google Scholar] [CrossRef]
  29. Ashaolu, E.D.; Olorunfemi, J.F.; PaulIfabiy, I.; Abdollahi, K.; Batelaan, O. Spatial and temporal recharge estimation of the basement complex in Nigeria, West Africa. J. Hydrol. Reg. Stud. 2020, 27, 100658. [Google Scholar] [CrossRef]
  30. Gonçalves, R.D.; Teramoto, E.H.; Engelbrecht, B.Z.; Alfaro Soto, M.A.; Chang, H.K.; van Genuchten, M.T. Quasi-Saturated Layer: Implications for Estimating Recharge and Groundwater Modeling. Groundwater 2019, 58, 432–440. [Google Scholar] [CrossRef]
  31. Ouyang, Y.; Jin, W.; Grace, J.M.; Obalum, S.E.; Zipperer, W.C.; Huang, X. Estimating impact of forest land on groundwater recharge in a humid subtropical watershed of the Lower Mississippi River Alluvial Valley. J. Hydrol. Reg. Stud. 2019, 26, 100631. [Google Scholar] [CrossRef]
  32. Petronici, F.; Pujades, E.; Jurado, A.; Marcaccio, M.; Borgatti, L. Numerical modelling of the Mulino Delle Vene Aquifer (Northern Italy) as a tool for predicting the hydrogeological system behavior under different recharge conditions. Water 2019, 11, 2505. [Google Scholar] [CrossRef] [Green Version]
  33. Bizhanimanzar, M.; Leconte, R.; Nuth, M. Catchment-Scale Integrated Surface Water-Groundwater Hydrologic Modelling Using. Water 2020, 12, 363. [Google Scholar] [CrossRef] [Green Version]
  34. Flinchum, B.; Banks, E.; Hatch, M.; Batelaan, O.; Peeters, L.; Pasquet, S. Identifying recharge under subtle ephemeral features in flat-lying semi-arid region using a combined geophysical approach. Hydrol. Earth Syst. Sci. Discuss. 2020, 24, 4353–4368. [Google Scholar] [CrossRef]
  35. Ali, M.; Mubarak, S. Approaches and Methods of Quantifying Natural Groundwater Recharge—A Review. Asian J. Environ. Ecol. 2017, 5, 1–27. [Google Scholar] [CrossRef]
  36. Chang, C.F.; Rubin, Y. Regionalization with hierarchical hydrologic similarity and ex situ data in the context of groundwater recharge estimation at ungauged watersheds. Hydrol. Earth Syst. Sci. 2019, 23, 2417–2438. [Google Scholar] [CrossRef] [Green Version]
  37. Monte-Mor, R.C.A.; Palmier, L.R.; Pinto, E.J.A.; Lima, J.E.S. Estabilidade temporal da distribuição espacial da umidade do solo em uma bacia intermitente no semiárido de Minas Gerais. Rev. Bras. Recur. Hidr. 2012, 17, 101–113. [Google Scholar] [CrossRef]
  38. Lima, J.E.F.W.; Montenegro, S.M.G.L.; de Montenegro, A.A.A.; Koide, S. Comparative hydrology: Relationships among physical characteristics, hydrological behavior, and results of the SWAT model in different regions of Brazil. Rev. Bras. Geogr. Física 2014, 7, 1187–1195. [Google Scholar]
  39. Creed, I.F.; Spargo, A.T.; Jones, J.A.; Buttle, J.M.; Adams, M.B.; Beall, F.D.; Booth, E.G.; Campbell, J.L.; Clow, D.; Elder, K.; et al. Changing forest water yields in response to climate warming: Results from long-term experimental watershed sites across North America. Glob. Chang. Biol. 2014, 20, 3191–3208. [Google Scholar] [CrossRef] [Green Version]
  40. Ibrahim, B.; Wisser, D.; Barry, B.; Fowe, T.; Aduna, A. Hydrological predictions for small ungauged watersheds in the Sudanian zone of the Volta basin in West Africa. J. Hydrol. Reg. Stud. 2015, 4, 386–397. [Google Scholar] [CrossRef] [Green Version]
  41. Santos, R.M.; Koide, S. Avaliação da recarga de águas subterrâneas em ambiente de cerrado com base em modelagem numérica do fluxo em meio poroso saturado. Rev. Bras. Recur. Hidr. 2016, 21, 451–465. [Google Scholar] [CrossRef]
  42. Araujo, D.L. Avaliação dos Impactos da Explotação de Águas Subterrâneas na Bacia do Ribeirão Rodeador por Meio de Simulação Integrada Entre os Modelos SWAT e MODFLOW. Master’s Thesis, Universidade de Brasília, Brasília, Brazil, 2018. Available online: https://repositorio.unb.br/handle/10482/34399 (accessed on 7 April 2020).
  43. Zomlot, Z.; Verbeiren, B.; Huysmans, M.; Batelaan, O. Spatial distribution of groundwater recharge and base flow: Assessment of controlling factors. J. Hydrol. Reg. Stud. 2015, 4, 349–368. [Google Scholar] [CrossRef] [Green Version]
  44. Ramires, T.; Manzione, R.L. Groundwater recharge estimation using water budget method for Bauru Aquifer system in a Cerrado environmental protection area. Appl. Res. Agrotechnol. 2019, 12, 25–36. [Google Scholar] [CrossRef]
  45. Silva, C.O.F. Modelagem Espacial da Recarga de Águas Subterrâneas Sob Diferentes Usos e Coberturas da Terra. Master’s Thesis, Universidade Estadual Paulista, São Paulo, Brazil, 2019. Available online: http://hdl.handle.net/11449/190710 (accessed on 19 April 2020).
  46. Charles, E.G.; Behroozi, C.; Schooley, J.; Hoffman, J.L. A Method for Evaluating Ground Water Recharge Areas in New Jersey; Trenton; 1993. Available online: https://www.nj.gov/dep/njgs/pricelst/gsreport/gsr32.pdf (accessed on 5 April 2020).
  47. Santos, R.M. Recarga de Águas Subterrâneas em Ambiente de Cerrado: Estudo com Base em Modelagem Numérica e Simulação Hidrológica em uma Bacia Experimental (Distrito Federal). Ph.D. Thesis, Universidade de Brasília, Brasília, Brazil, 2012. Available online: https://repositorio.unb.br/handle/10482/11253 (accessed on 9 April 2020).
  48. CPRM; Embrapa. Zoneamento Ecológico—Econômico da Região Integrada de Desenvolvimento do Distrito Federal e Entorno: Fase I; Ministério do Meio Ambiente: Rio de Janeiro, Brazil, 2003; Volume 1.
  49. Reatto, A.; Martins, É.D.S.; Cardoso, E.A.; Spera, S.T.; Carvalho, O.A., Jr.; De Silva, A.V.; Farias, M.F.R. Levantamento de Reconhecimento de Solos de Alta Intensidade do Alto Curso do Rio Descoberto; Distrito Federal, Brasília. 2003. Available online: http://ainfo.cnptia.embrapa.br/digital/bitstream/CPAC-2009/25006/1/bolpd_92.pdf (accessed on 11 April 2020).
  50. Liu, Y.B.; De Smedt, F. WetSpa Extension, a GIS-Based Hydrologic Model for Flood Prediction and Watershed Management Documentation and User Manual; Brussel, Belgium. 2004. Available online: https://www.vub.be/WetSpa/downloads/WetSpa_manual.pdf (accessed on 3 April 2020).
  51. Healy, R.W.; Cook, P.G. Using groundwater levels to estimate recharge. Hydrogeol. J. 2002, 10, 91–109. [Google Scholar] [CrossRef]
  52. Bailey, R.T. TUTORIAL SWAT-MODFLOW. Documentación para la Preparación del Acoplamiento SWAT-MODFLOW. 2015. Available online: https://swat.tamu.edu/media/115188/swat-modflow-tutorial-in-spanish.pdf (accessed on 24 April 2020).
  53. Brighenti, T.M. Modelagem hidrológica e avaliação de diferentes métodos de calibração para o modelo SWAT. Master’s Thesis, Universidade Federal de Santa Maria, São Paulo, Brazil, 2015. Available online: https://repositorio.ufsc.br/xmlui/handle/123456789/136347 (accessed on 17 April 2020).
  54. Arroio Junior, P.P. Aprimoramento das Rotinas E Parâmetros dos Processos Hidrológicos do Modelo Computacional Soil and Water Assessment Tool—Swat. Ph.D. Thesis, Universidade de São Paulo, São Paulo, Brazil, 2016. Available online: https://teses.usp.br/teses/disponiveis/18/18139/tde-25052017-084925/publico/TesePauloPonceArroioJunior.pdf (accessed on 29 April 2020).
  55. Yifru, B.A.; Chung, I.M.; Kim, M.G.; Chang, S.W. Assessment of groundwater recharge in agro-urban watersheds using integrated SWAT-MODFLOW model. Sustainbility 2020, 12, 6593. [Google Scholar] [CrossRef]
  56. Nossent, J.; Bauwens, W. Comparing SWAT and WetSpa on the River Grote Laak, Belgium. In Proceedings of the 4th International SWAT conference, Delft, The Netherlands, 2–6 July 2007. [Google Scholar]
  57. Polanco, E.I.; Fleifle, A.; Ludwig, R.; Disse, M. Improving SWAT model performance in the Upper Blue Nile River Basin using meteorological data integration and catchment scaling. Hydrol. Earth Syst. Sci. Discuss. 2017, 21, 4907–4926. [Google Scholar] [CrossRef] [Green Version]
  58. De Castro, K.B. Avaliação do Modelo SWAT na Simulaçaõ de Vazão em Bacia Agricola do Cerrado Intensamente Monitorada. Master’s Thesis, Universidade de Brasília, Brasília, Brazil, 2013. Available online: https://repositorio.unb.br/bitstream/10482/13863/1/2013_KassiaBatistaCastro.pdf (accessed on 25 April 2020).
  59. Gaspar, M.T.P. Sistema Aqüífero Urucuia: Caracterização Regional E Propostas De Gestão Sistema Aqüífero Urucuia. Ph.D. Thesis, Universidade de Brasília, Brasília, Brazil, 2006. Available online: https://repositorio.unb.br/handle/10482/6742 (accessed on 2 May 2020).
  60. Albuquerque, A.C.L.S.; Chaves, H.M.L. Estimativa de recarga da bacia do Rio das Fêmeas através de métodos manuais e automáticos Estimates of the recharge in the Femeas River basin by manual and automatic methods. Rev. Bras. Eng. Agrícola Ambient—Agriambi 2011, 15, 1123–1129. [Google Scholar] [CrossRef] [Green Version]
  61. Cunha, V.C.V. Avaliação da Interação Entre Águas Subterrâneas e Superficiais na Bacia do rio das Fêmeas, Sistema Aquífero Urucuia—Bahia, Centro de Desenvolvimento da Tecnologia Nuclear. Ph.D. Dissertations; 2017. Available online: http://rigeo.cprm.gov.br/xmlui/bitstream/handle/doc/18596/diss_viviane_cunha.pdf?sequence=1&isAllowed=y (accessed on 8 April 2020).
  62. Oliveira, P.T.S.; Wendland, E.; Nearing, M.A.; Scott, R.L.; Rosolem, R.; Da Rocha, H.R. The water balance components of undisturbed tropical woodlands in the Brazilian cerrado. Hydrol. Earth Syst. Sci. 2015, 19, 2899–2910. [Google Scholar] [CrossRef] [Green Version]
  63. Agência Nacional de Águas. Estudos hidrogeológicos e de vulnerabilidade do Sistema Aquífero Urucuia e proposição de modelo de gestão integrada compartilhada: Resumo executivo. 2017; 100p. Available online: https://metadados.ana.gov.br/geonetwork/srv/en/resources.get?id=499&fname=Resumo_Executivo_SAU.pdf&access=private (accessed on 16 April 2020).
  64. Agência Nacional de Águas. Estratégias de Manejo Sustentável dos Sistemas Aquíferos Urucuia e Areado e Conclusões. Relatório Final. Brasília; 2017; Volume 3. Available online: https://metadados.ana.gov.br/geonetwork/srv/en/resources.get?id=499&fname=Volume_3.pdf&access=private (accessed on 16 April 2020).
Figure 1. Location of the study area and the validation area.
Figure 1. Location of the study area and the validation area.
Water 12 03090 g001
Figure 2. Steps of the numerical modelling method (A: items of the conceptual hydrogeological model, B: steps of the numerical model implementation, C: database partition between modeling steps, D: steps of the calibration process, E: calibration process evaluation, F: recharge estimation by inverse modelling) [41,47].
Figure 2. Steps of the numerical modelling method (A: items of the conceptual hydrogeological model, B: steps of the numerical model implementation, C: database partition between modeling steps, D: steps of the calibration process, E: calibration process evaluation, F: recharge estimation by inverse modelling) [41,47].
Water 12 03090 g002
Figure 3. Hydrological processes and WETSPA (Water and Energy Transfer between Soils, Plants, and Atmosphere) model structure [50].
Figure 3. Hydrological processes and WETSPA (Water and Energy Transfer between Soils, Plants, and Atmosphere) model structure [50].
Water 12 03090 g003
Figure 4. (a) Spatial distribution of the simulated recharge; (b) Spatial distribution of the regionalized recharge.
Figure 4. (a) Spatial distribution of the simulated recharge; (b) Spatial distribution of the regionalized recharge.
Water 12 03090 g004
Figure 5. Behavior of the recharge in function of the factors “soil class”, “land use/cover”, and “slope”.
Figure 5. Behavior of the recharge in function of the factors “soil class”, “land use/cover”, and “slope”.
Water 12 03090 g005
Figure 6. (a) Relative error of the regionalized recharge and its relation with the factors (b) “slope”, (c) “soil class”, and (d) “land use/cover”.
Figure 6. (a) Relative error of the regionalized recharge and its relation with the factors (b) “slope”, (c) “soil class”, and (d) “land use/cover”.
Water 12 03090 g006
Table 1. Potential recharge rates for each combination of “soil class”, “land use/cover”, and “slope” for the Cerrado Biome in the Federal District.
Table 1. Potential recharge rates for each combination of “soil class”, “land use/cover”, and “slope” for the Cerrado Biome in the Federal District.
CodeSimulated Recharge
SWAT-MODFLOW
(% Annual Prec.)
Regionalized Recharge WETSPA
(% Annual Prec.)
Relative Error (%)
SWAT-MODFLOW “Versus” WETSPA
11156.9623.74−58.33
11256.6519.86−64.93
11356.8617.76−68.76
11455.8415.39−72.43
12133.5024.32−27.40
12230.3820.20−33.52
12332.8317.95−45.32
12432.1014.49−54.84
13130.9516.77−45.81
13230.0720.20−32.82
13333.2918.11−45.59
13432.8014.97−54.36
14159.8223.78−60.24
14260.5521.73−64.11
14360.9220.54−66.28
14459.9320.83−65.25
16159.5420.44−65.67
16259.8616.90−71.77
16360.4218.33−69.65
16459.8515.60−73.93
21153.6338.36−28.47
21253.0637.45−29.41
21351.0936.00−29.53
21449.6734.85−29.82
22153.5838.80−27.58
22253.1937.36−29.76
22350.6835.97−29.02
22451.1635.77−30.09
23146.2237.12−19.68
23245.1538.72−14.24
23345.6138.12−16.41
23443.8038.21−12.76
24150.9436.15−29.04
24251.9539.18−24.58
26149.3038.00−22.92
26250.0738.70−22.71
26349.7137.33−24.90
26448.1639.59−17.81
28129.2736.3023.99
28219.8638.9696.11
31148.8838.46−21.33
31232.1937.1915.51
33130.2637.7524.74
34151.6938.11−26.26
34253.4037.63−29.53
34356.5436.89−34.76
36148.2638.23−20.79
36247.8937.46−21.78
41138.9435.50−8.83
41238.6929.66−23.33
41337.8028.95−23.41
41437.1425.27−31.97
42139.1536.39−7.05
42238.9329.82−23.39
42337.6729.03−22.93
42421.5127.7428.99
43132.8435.869.20
43233.2228.21−15.09
43331.7027.94−11.85
43431.4024.92−20.63
44135.9834.97−2.79
44237.9726.93−29.08
44337.4725.59−31.70
44436.0522.63−37.24
46135.2536.112.45
46236.2229.31−19.08
46335.5228.76−19.03
46434.8523.97−31.20
47133.1235.627.54
47232.9929.18−11.55
48135.4934.84−1.84
51144.9338.12−15.15
51247.5837.69−20.79
51348.2936.86−23.67
51446.9631.10−33.76
53134.6538.039.74
53225.2937.0246.39
54149.5636.20−26.97
54252.6734.51−34.48
54352.6030.51−41.98
54450.1422.31−55.49
56154.1338.25−29.33
56249.9736.81−26.34
56353.1233.62−36.72
56453.1628.61−46.19
Table 2. Regionalized average effective recharge rates and average literature reference values obtained for similar environmental conditions (Cerrado Biome).
Table 2. Regionalized average effective recharge rates and average literature reference values obtained for similar environmental conditions (Cerrado Biome).
Abacus Average Effective Recharge
(% prec.)
MethodAverage Reference Values (% prec.)Method
15WTE15–25 [59]WTE
19NM20 [60]Baseflow
17 [61]Baseflow
26–30 [62]WTE
16 [63]Baseflow
17 [63]Baseflow
19 [42]NM
16 [7]WTE
24 [8]Baseflow
27 [8]WTE
Average = 17Average = 20
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Santos, R.M.d.; Koide, S.; Távora, B.E.; Araujo, D.L.d. Abacus to Predict Groundwater Recharge at Non-Instrumented Hydrographic Basins. Water 2020, 12, 3090. https://doi.org/10.3390/w12113090

AMA Style

Santos RMd, Koide S, Távora BE, Araujo DLd. Abacus to Predict Groundwater Recharge at Non-Instrumented Hydrographic Basins. Water. 2020; 12(11):3090. https://doi.org/10.3390/w12113090

Chicago/Turabian Style

Santos, Ronaldo Medeiros dos, Sérgio Koide, Bruno Esteves Távora, and Daiana Lira de Araujo. 2020. "Abacus to Predict Groundwater Recharge at Non-Instrumented Hydrographic Basins" Water 12, no. 11: 3090. https://doi.org/10.3390/w12113090

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