Conceptual and numerical modeling approach of the Guarani Aquifer System
In large aquifers, relevant for their considerable size, regional groundwater modeling remains challenging given geologic complexity and data scarcity in space and time. Yet, it may be conjectured that regional scale groundwater flow models can help in understanding the flow system functioning and the relative magnitude of water budget components, which are important for aquifer management. The Guaraní Aquifer System is the largest transboundary aquifer in South America. It contains an enormous volume of water; however, it is not well known, being difficult to assess the impact of exploitation currently used to supply over 25 million inhabitants. This is a sensitive issue because the aquifer is shared by four countries. Moreover, an integrated groundwater model, and therefore a global water balance, were not available. In this work, a transient regional scale model for the entire aquifer based upon five simplified, equally plausible conceptual models represented by different hydraulic conductivity parametrizations is used to analyze the flow system and water balance components. Combining an increasing number of hydraulic conductivity zones and an appropriate set of boundary conditions, the hypothesis of a continuous sedimentary unit yielded errors within the calibration target in a regional sense. The magnitude of the water budget terms resulted very similar for all parametrizations. Recharge and stream/aquifer fluxes were the dominant components representing, on average, 84.2% of total inflows and 61.4% of total outflows, respectively. However, leakage was small compared to stream discharges of main rivers. For instance, the simulated average leakage for the Uruguay River was 8 m 3 s −1 while the observed absolute minimum discharge was 382 m 3 s −1. Streams located in heavily pumped regions switched from a gaining condition in early years to a losing condition over time. Water is discharged through the aquifer boundaries, except at the eastern boundary. On average, pumping represented 16.2% of inflows while aquifer storage experienced a small overall increment. The model water balance indicates that the current rate of groundwater withdrawals does not exceed the rate of recharge in a regional sense.