Numerical groundwater-flow modeling to evaluate potential effects of pumping and recharge: implications for sustainable groundwater management in the Mahanadi delta region,India

Process-based groundwater models are useful to understand complex aquifer systems and make predictions about their response to hydrological changes. A conceptual model for evaluating responses to environmental changes is presented, considering the hydrogeologic framework, flow processes, aquifer hydraulic properties, boundary conditions, and sources and sinks of the groundwater system. Based on this conceptual model, a quasi-three-dimensional transient groundwater flow model was designed using MODFLOW to simulate the groundwater system of Mahanadi River delta, eastern India. The model was constructed in the context of an upper unconfined aquifer and lower confined aquifer, separated by an aquitard. Hydraulic heads of 13 shallow wells and 11 deep wells were used to calibrate transient groundwater conditions during 1997–2006, followed by validation (2007–2011). The aquifer and aquitard hydraulic properties were obtained by pumping tests and were calibrated along with the rainfall recharge. The statistical and graphical performance indicators suggested a reasonably good simulation of groundwater flow over the study area. Sensitivity analysis revealed that groundwater level is most sensitive to the hydraulic conductivities of both the aquifers, followed by vertical hydraulic conductivity of the confining layer. The calibrated model was then employed to explore groundwater-flow dynamics in response to changes in pumping and recharge conditions. The simulation results indicate that pumping has a substantial effect on the confined aquifer flow regime as compared to the unconfined aquifer. The results and insights from this study have important implications for other regional groundwater modeling studies, especially in multi-layered aquifer systems.     

Hydrogeology and sustainable future groundwater abstraction from the Agua Verde aquifer in the Atacama Desert,northern Chile

The hyper-arid conditions prevailing in Agua Verde aquifer in northern Chile make this system the most important water source for nearby towns and mining industries. Due to the growing demand for water in this region, recharge is investigated along with the impact of intense pumping activity in this aquifer. A conceptual model of the hydrogeological system is developed and implemented into a two-dimensional groundwater-flow numerical model. To assess the impact of climate change and groundwater extraction, several scenarios are simulated considering variations in both aquifer recharge and withdrawals. The estimated average groundwater lateral recharge from Precordillera (pre-mountain range) is about 4,482 m³/day. The scenarios that consider an increase of water withdrawal show a non-sustainable groundwater consumption leading to an over-exploitation of the resource, because the outflows surpasses inflows, causing storage depletion. The greater the depletion, the larger the impact of recharge reduction caused by the considered future climate change. This result indicates that the combined effects of such factors may have a severe impact on groundwater availability as found in other groundwater-dependent regions located in arid environments. Furthermore, the scenarios that consider a reduction of the extraction flow rate show that it may be possible to partially alleviate the damage already caused to the aquifer by the continuous extractions since 1974, and it can partially counteract climate change impacts on future groundwater availability caused by a decrease in precipitation (and so in recharge), if the desalination plant in Taltal increases its capacity.     

利用环境同位素识别酒泉-张掖盆地地下水补给和水流系统

详细了解干旱区地下水的补给机制对地下水资源管理来说是非常重要的,天然环境同位素在过去的40年里广泛应用于解决有关地下水补给、流动等问题。笔者通过分析酒泉-张掖盆地水的环境同位素变化特征,识别地下水的补给和流动.研究结果表明,酒泉-张掖盆地地下水的补给源主要来自山区出山河流,补给发生在祁连山前戈壁带,在山前冲积扇带地下水主要为1963年以来的快速补给,而部分深部地下水为1952年以前补给。同时,环境同位素指示盆地存在深部区域流系统和浅部局域流系统,酒泉盆地和张掖盆地之间不存在明显的水力联系。该研究不仅对黑河流域地下水的开发管理有着重要意义,对我国西北类似的内陆盆地地下水的开发管理有着借鉴意义。     

Modelling hydrostratigraphy and groundwater flow of a fractured and karst aquifer in a Mediterranean basin (Salento peninsula, southeastern Italy)

The control exerted by the hydrostratigraphic structure on aquifer recharge, groundwater flow and discharge along the coastal areas of a Mediterranean basin (Salento peninsula, about 5,000?km² wide, southern Italy) is assessed through the development and application of a groundwater flow model based on the reconstruction of the hydrostratigraphic architecture at the regional scale. The hydrostratigraphic model, obtained by processing surface and subsurface data, is applied to map the top of the main aquifer, which is hosted in the deep hydrostratigraphic unit corresponding to Cretaceous and Oligocene limestones with complex geometrical relationships with the sea. It is also used to estimate the aquifer recharge, which occurs by percolation through overlying younger sediments with low permeability. These data are completed with information about the soil use to estimate water abstraction for irrigation and with literature data to estimate the water abstraction for drinking and industrial purposes. The above-sketched conceptual model is the basis for a finite difference groundwater 2D pseudo-stationary flow model, which assumes the following fundamental approximations: the fractured and karst limestone hydrostratigraphic unit can be approximated, at the model scale, as a continuous medium for which the discrete Darcy??s law is valid; the transition zone between salt and fresh water is so small with respect to the grid spacing that the Ghyben?CHerzberg??s approximation for a sharp interface can be applied. Along the coastline different boundary conditions are assigned if the top of the limestone hydrostratigraphic unit lies either above the sea level (the aquifer has a free surface and fresh water is drained), or below the sea level (the aquifer is under pressure and the contact with sea occurs off-shore). The groundwater flow model correctly predicts the areas where the aquifer is fully saturated with salt water.     

Groundwater resource sustainability in the Wadi Watir delta, Gulf of Aqaba, Sinai, Egypt

The Wadi Watir delta, in the arid Sinai Peninsula, Egypt, contains an alluvial aquifer underlain by impermeable Precambrian basement rock. The scarcity of rainfall during the last decade, combined with high pumping rates, resulted in degradation of water quality in the main supply wells along the mountain front, which has resulted in reduced groundwater pumping. Additionally, seawater intrusion along the coast has increased salinity in some wells. A three-dimensional (3D) groundwater flow model (MODFLOW) was calibrated using groundwater-level changes and pumping rates from 1982 to 2009; the groundwater recharge rate was estimated to be 1.58?×?10⁶ m³/year. A variable-density flow model (SEAWAT) was used to evaluate seawater intrusion for different pumping rates and well-field locations. Water chemistry and stable isotope data were used to calculate seawater mixing with groundwater along the coast. Geochemical modeling (NETPATH) determined the sources and mixing of different groundwaters from the mountainous recharge areas and within the delta aquifers; results showed that the groundwater salinity is controlled by dissolution of minerals and salts in the aquifers along flow paths and mixing of chemically different waters, including upwelling of saline groundwater and seawater intrusion. Future groundwater pumping must be closely monitored to limit these effects.     

Regional groundwater flow model for Abu Dhabi Emirate: scenario-based investigation

Despite the continuous increase in water supply from desalination plants in the Emirate of Abu Dhabi, groundwater remains the major source of fresh water satisfying domestic and agricultural demands. Groundwater has always been considered as a strategic water source towards groundwater security in the Emirate. Understanding the groundwater flow system, including identification of recharge and discharge areas, is a crucial step towards proper management of this precious source. One main tool to achieve such goal is a groundwater model development. As such, the main aim of this paper is to develop a regional groundwater flow model for the surficial aquifer in Abu Dhabi Emirate using MODFLOW. Up to our knowledge, this is the first regional numerical groundwater flow model for Abu Dhabi Emirate. After steady state and transient model calibration, several future scenarios of recharge and pumping are simulated. Results indicate that groundwater pumping remains several times higher than aquifer recharge from rainfall, which provides between 2 and 5% of total aquifer recharge. The largest contribution of recharge is due to subsurface inflow from the eastern Oman Mountains. While rainfall induced groundwater level fluctuation is absent in the western coastal region, it reaches a maximum of 0.5 m in the eastern part of the Emirate. In contrast, over the past decades, groundwater levels have declined annually by 0.5 m on average with local extremes spanning from 93 m of decline to 60 m of increase. Results also indicate that a further decrease in groundwater levels is expected in most of Emirate. At other few locations, upwelling of groundwater is expected due to a combination of reduced pumping and increased infiltration of water from nonconventional sources. Beyond results presented here, this regional groundwater model is expected to provide an effective tool to water resources managers in Abu Dhabi. It will help to accurately estimate sustainable extraction rates, assess groundwater availability, and identify pathways and velocity of groundwater flow as crucial information for identifying the best locations for artificial recharge.     

Groundwater flow modeling of Quaternary aquifer Ras Sudr,Egypt

Recently, Ras Sudr (the delta of Wadi Sudr) area received a great amount of attention due to different development expansion activities directed towards this area. Although Quaternary aquifer is the most prospective aquifer in Ras Sudr area, it has not yet completely evaluated. The present work deals with the simulation of the Quaternary groundwater system using a three-dimensional groundwater flow model. MODFLOW code was applied for designing the model of the Ras Sudr area. This is to recognize the groundwater potential as well as exploitation plan of the most prospective aquifer in the area. The objectives were to determine the hydraulic parameters of the Quaternary aquifer, to estimate the recharge amount to the aquifer, and to determine the hydrochemistry of groundwater in the aquifer. During this work, available data has been collected and some field investigation has been carried out. Groundwater flow model has been simulated using pilot points’ method. SEAWAT has been also applied to simulate the variable-density flow and sea water intrusion from the west. It can be concluded that: (1) the direction of groundwater flow is from the east to the west, (2) the aquifer system attains a small range of log-transformed hydraulic conductivity. It ranges between 3.05 and 3.35 m/day, (3) groundwater would be exploited by about 6.4 × 10⁶ m³/year, (4) the estimated recharge accounts for 3 × 10⁶ m³/year, (5) an estimated subsurface flow from the east accounts for 2.7 × 10⁶ m³/year, (6) the increase of total dissolved solids (TDS) most likely due to dilution of salts along the movement way of groundwater from recharge area to discharge area in addition to a contribution of sea water intrusion from the west. Moreover, it is worth to note that a part of TDS increase might be through up coning from underlying more saline Miocene sediments. It is recommended that: (1) any plan for increasing groundwater abstraction is unaffordable, (2) reliable estimates of groundwater abstraction should be done and (3) automatic well control system should be made.     

Groundwater vulnerability map of the Chrzanów karst-fissured Triassic aquifer (Poland)

A map shows intrinsic vulnerability to pollution of the Chrzanów karst-fissured aquifer (273 km²) in the southern part of Poland. This aquifer is intensively drained by numerous intakes and Zn-Pb ore mines. A DRASTIC-type parametric system was applied for groundwater vulnerability evaluation. Vulnerability assessment is based on six factors (depth to groundwater table, lithology of the unsaturated zone, net recharge, hydraulic conductivity of the aquifer, groundwater flow velocity, aquifer thickness). For the final vulnerability map construction at the scale of 1:50,000, a combination of the aquifer simulation model (using MODFLOW) and a geographical information system was applied. Maps of the net recharge, hydraulic conductivity of the aquifer and groundwater flow velocity were derived by aquifer modelling. Based on the vulnerability index (21-182), six relative vulnerability classes were selected. Reliability of the map has been verified.     

Evaluation of the importance of clay confining units on groundwater flow in alluvial basins using solute and isotope tracers: the case of Middle San Pedro Basin in southeastern Arizona (USA)

As groundwater becomes an increasingly important water resource worldwide, it is essential to understand how local geology affects groundwater quality, flowpaths and residence times. This study utilized multiple tracers to improve conceptual and numerical models of groundwater flow in the Middle San Pedro Basin in southeastern Arizona (USA) by determining recharge areas, compartmentalization of water sources, flowpaths and residence times. Ninety-five groundwater and surface-water samples were analyzed for major ion chemistry (water type and Ca/Sr ratios) and stable (¹⁸O, ²H, ¹³C) and radiogenic (³H, ¹⁴C) isotopes, and resulting data were used in conjunction with hydrogeologic information (e.g. hydraulic head and hydrostratigraphy). Results show that recent recharge (<60 years) has occurred within mountain systems along the basin margins and in shallow floodplain aquifers adjacent to the San Pedro River. Groundwater in the lower basin fill aquifer (semi confined) was recharged at high elevation in the fractured bedrock and has been extensively modified by water-rock reactions (increasing F and Sr, decreasing ¹⁴C) over long timescales (up to 35,000 years BP). Distinct solute and isotope geochemistries between the lower and upper basin fill aquifers show the importance of a clay confining unit on groundwater flow in the basin, which minimizes vertical groundwater movement.     

Surface and subsurface conceptual model of an arid environment with respect to mid- and late Holocene climate changes

The water demand in arid regions is commonly covered by groundwater resources that date back to more humid periods of the Pleistocene and Holocene. Within the investigated arid part of SE Saudi-Arabia information about climate, groundwater levels, and pumping rates are only available for regions where groundwater extractions occur at present-day. For the prediction of the impact of long-term climate changes on groundwater resources an understanding of the hydrogeological and hydrological past and the development of the aquifers is necessary. Therefore, all available information about hydrology and hydrogeology for the past 10,000 years BP were collected and compiled to a conceptual model of the aquifer development on the Arabian Peninsula since the last Ice-Age. The climatic history was displayed by changes in precipitation, temperature and recharge during the mid-S and late Holocene. The hydrogeological development is described by groundwater ages, sea level fluctuations, movement of the coastline, and the development of sabkhas. The most sensitive parameter to describe the development of aquifer system is recharge. Present-day recharge was calculated with the hydrological model system HEC-HMS accounting for current precipitation, temperature, wind, soil types, and geomorphology. With respect to changes in precipitation and temperature over the past 10,000 years the temporal and spatial variability of groundwater recharge was calculated using empirical equations valid for semi-arid and arid settings. Further inflow into the groundwater system results from surface water infiltration in wadi beds, while natural outflow from the groundwater system occurs by discharge to the Gulf, evaporation from sabkhas, and spring discharge. Backward predictions can be verified by sedimentological observations of palaeo-river systems and lakes indicating that groundwater levels reached temporarily the surface under wetter climate conditions and ¹⁴C groundwater ages displaying groundwater residence times.     

Assessing the vulnerability of over-exploited volcanic aquifer systems using multi-parameter analysis,Toluca Basin,Mexico

The potential for contamination of groundwater in an over-exploited system led to this investigation on aspects of aquifer vulnerability in the Toluca Basin, Mexico. This study involved the use of various hydrogeological, geochemical and isotopic tools to better understand the susceptibility to contamination where heavy groundwater extraction along with industrial extensive development is concentrated. Geochemical and isotopic analyses of groundwater from production wells along the mountain boundary show little evidence of contamination at present time. Radiocarbon and tritium data collected in piezometers and wells penetrating the Lower basalt aquifer indicate that the groundwater is of varying age and quality. These data, in addition to stable isotope and chemical data also indicate evidence of contamination in the Lower aquifer; which could be associated with fast recharge in the mountains or injection of waste water by illegal wells. Evidence found in this study of low downward gradients in the valley fill sediments and the presence of low permeability aquitard deposits suggest that the aquifer system is currently not highly threatened by surface sources of contamination. However, where urban and industrial development occurs, in the centre of the valley and along the flanks of the Basin, the vulnerability of the aquifer system increases significantly.     

Desalination of a sedimentary rock aquifer system invaded by Pleistocene Champlain Sea water and processes controlling groundwater geochemistry

The objective of this study was to identify geochemical processes and Quaternary geological events responsible for the variations in groundwater geochemistry observed in a sedimentary rock aquifer system, including brackish to saline groundwater. Inorganic constituents and environmental isotopes were analyzed for 146 groundwater samples. Dissolution of carbonates dominates in recharge areas, resulting in Ca-, Mg-HCO₃ groundwater. Further along flow paths, under confined conditions, Ca²⁺–Na⁺ ion exchange causes groundwater evolution to Na-HCO₃ type. Na-Cl groundwater is also found and it falls on a seawater mixing line. Using conservative tracers, Cl⁻ and Br⁻, the original Champlain Sea water is shown to have been, in the region, a mixture of about 34% seawater and 66% freshwater, a composition still retained by some groundwater. Na-Cl groundwater thus results from mixing with former Champlain Sea water and also from solute diffusion from overlying marine clay. The system is thus found to be at different stages of desalinization, from the original Champlain Sea water still present in hydraulically stagnant areas of the aquifer to fully flushed conditions in parts, where more flow occurs, especially in recharge zones. The geochemical processes are integrated within the hydrogeological context to produce a conceptual geochemical evolution model for groundwater of the aquifer system.     

Groundwater flow modelling of Hindon-Yamuna interfluve region, Western Uttar Pradesh

The study area Hindon -Yamuna interfluve region is underlain by a thick pile of unconsolidated Quaternary alluvial deposits and host multiple aquifer system. Excessive pumping in the last few decades, mainly for irrigation, has resulted in a significant depletion of the aquifer. Therefore, proper groundwater management of Hindon-Yamuna interfluve region is necessary. For effective groundwater management of a basin it is essential that careful zone budget study should be carried out. Keeping this in view, groundwater flow modelling was attempted to simulate the behavior of flow system and evaluate zone budget. Visual MODFLOW, pro 4.1 is used in this study to simulate groundwater flow. The model simulates groundwater flow over an area of about 1345 km² with a uniform grid size of 1000 m by 1000 m and contains three layers, 58 rows and 37 columns. The horizontal flows, seepage losses from unlined canals, recharge from rainfall and irrigation return flows were applied using different boundary packages available in Visual MODFLOW, pro 4.1. The river — aquifer interaction was simulated using the river boundary package. Simulated pumping rates of 500 m³/day, 1000 m³/day and 1500 m³/day were used in the pumping well package.The zone budget for the steady state condition of study area indicated that the total annual direct recharge is 416.10 MCM and the total annual groundwater draft through pumping is of the order of 416.63 MCM. Two scenarios were considered to predict aquifer system response under different conditions. Sensitivity analysis on model parameters was conducted to quantitatively evaluate the impact of varying model inputs. Based on the results obtained from the sensitivity analysis, it was found that the model is more sensitive to hydraulic conductivity and recharge parameter. Present study deals with importance of groundwater modelling for planning, design, implementation and management of groundwater resources.     

大清河流域平原区地下水人工补给潜力与补给方式分析

为了有效提升大清河流域平原区地下水水位,亟需在此区域开展地下水人工补给工程,并确定合理的建设位置及有效的补给方式。首先基于研究区可利用补给水源、地下水位、地表高程、地表坡度及与河道距离5个指标的分布特征,构建地下水补给潜力评价体系,采用ArcGIS空间分析功能对研究区进行了地下水人工补给潜力区划;然后在此评价体系基础上,在典型人工补给高潜力区进一步开展系列野外现场试验,探讨适宜可行的地下水人工补给方式。结果表明：研究区西北部及南部河道附近区域开展人工补给工程潜力较高,而中部、北部及西南部远离河道的区域潜力较低。高潜力区——白沟引河地段包气带及含水层渗透性良好,整体渗透系数均在5 m/d左右或更高,适宜地表补给,但河床渗透性较差,渗透系数基本在0.01~0.09 m/d间,若通过河道补给需配合清淤等措施。其中,在上游及中游沿岸适宜将河道水通过生态水渠引至修建的地表入渗池或借助天然渗坑内入渗补给,在中下游沿岸区域适宜将补给水进行严格的水处理后采用井灌方式补给,在白沟引河中下游河道适宜修建拦水坝,利用河道进行入渗补给。     

Why conceptual groundwater flow models matter: a trans-boundary example from the arid Great Basin, western USA

Spring and Snake valleys, western USA, are scheduled for development and groundwater export to Las Vegas, Nevada (USA). New work, compared to published studies, illustrates the critical role of conceptual models to underpin water withdrawals in arid regions. Interbasin flow studies suggest that 30–55?% of recharge to Snake Valley arrives from adjacent Spring Valley. This study, however, suggest little or no interbasin flow; rather, Spring and Snake valleys comprise separate systems. Contrary to expectation, δD and δ¹⁸O contours are perpendicular to proposed interbasin flow paths. ¹⁴C age gradients up to 10?ka along interbasin flow paths indicate that old waters are not displaced by such fluxes. ¹⁴C and ³H patterns indicate local recharge occurs in adjacent mountain ranges and is transferred to basin-fill by losing streams, mountain front recharge, and upward leakage from carbonate bedrock beneath basins. The choice of conceptual models is critical for groundwater development. Simple analyses of water withdrawals indicate that monitoring discharges at desert springs is an inadequate protective measure. Once flows decline, recovery is lengthy even if pumping is stopped. The conceptual framework behind quantitative evaluations of sustainable yield is critical to determine the ability of a groundwater system to deliver sustained withdrawals.     

Hydrogeological behaviour of some fault zones in a carbonate aquifer of Southern Italy: an experimentally based model

Complex flow circulation patterns are likely to be present in fault‐controlled groundwater flow systems, such as carbonate aquifers. Nevertheless, not much information is available for faults in carbonates, and their hydrogeological behaviour is often neglected in conceptual and numerical models. The understanding of this aspect of subsurface fluid flow has been improved in a carbonate aquifer, where hydrogeological investigations at site scale demonstrated the existence of fault zones that act as barriers. The hydraulic conductivity of the fault core is as low as that of siliciclastic rocks that represent the regional aquitard of the carbonate aquifer. Despite the lower permeability, the fault zones allow a significant groundwater flowthrough and a good interdependence of piezometric heads upgradient and downgradient of the faults. Because of this discontinuous heterogeneity, the aquifer looks like a basins‐in‐series system, where seasonal springs can be detected along some fault zones, as a function of groundwater level fluctuations.     

Hydrostratigraphic framework and hydrogeological behaviour of the Mancha Oriental System (SE Spain)

The Mancha Oriental System (MOS) has a surface area of 7,260 km², making it one of the largest carbonate aquifer systems in Spain. The system sustains about 1,000 km² of irrigated crops and supplies groundwater to 275,000 inhabitants. The economic transformation brought about by the development of extensive irrigated cropland has led to a water-balance disequilibrium of about 75 million m³/year. This input–output deficit has negative consequences in the quantity and quality of the available resources, in the river–aquifer relationship, and in the associated ecosystems as well. To understand the hydrogeological behaviour of the system, it is necessary to design a conceptual model. Further, the conceptualisation of a groundwater flow system is a requirement of the European Water Framework Directive for the characterisation of groundwater bodies. The robustness of the conceptual model depends heavily on the user capability of representing the real system. In this work, a multidisciplinary approach has been used to represent the three-dimensional geological framework and the groundwater flow conceptualisation of the MOS. Data management and three-dimensional visualisation have been carried out by means of geographical information system (GIS) tools and software for contouring and three-dimensional surface mapping.     

A conceptual and numerical model for groundwater management: a case study on a coastal aquifer in southern Tuscany, Italy

Ongoing hydrogeological research aims to develop a correct management model for the Plio-Pleistocene multi-aquifer system of the Albegna River coastal plain (southern Tuscany, Italy); overexploitation of this aquifer for irrigation and tourism has caused seawater intrusion. The conceptual model is based on field and laboratory data collected during the 1995–2003 period. Meteoric infiltration and flows from the adjoining carbonate aquifer recharge the aquifer. Natural outflow occurs through a diffuse flow into the sea and river; artificial outflow occurs through intensive extraction of groundwater from wells. Water exchanges in the aquifer occur naturally (leakage, closing of aquitard) and artificially (multiscreened wells). The aquifer was represented by a three-dimensional finite element model using the FEFLOW numerical code. The model was calibrated for steady-state and transient conditions by matching computed and measured piezometric levels (February 1995–February 1996). The model helped establish that seawater intrusion is essentially due to withdrawals near the coast during the irrigation season and that it occurs above all in the Osa-Albegna sector, as well as along the river that at times feeds the aquifer. The effects of hypothetical aquifer exploitation were assessed in terms of water budget and hydraulic head evolution.     

地下水机器学习方法研究

随着我国地下水监测工作的高速发展,高频率高密度水位监测数据的出现催生了对其进行深入信息挖掘的需求。在传统地下水模型研究中,地下水水位监测值常位于模型构建过程的下游,当水位监测的时空密度逐渐增大时,新增信息无法有效传导至模型的规划阶段并指导概念模型的修订。文章提出了一种地下水系统补排边界的识别方法,在不建立地下水数值模型的前提下,以监测井空间位置为节点,按照德劳内原则建立三角网格。在此网格系统中,首先定义一个水力梯度变换函数gradF,以求取网格中任意位置的水力梯度;借鉴机器学习领域的优化算法,使用水力梯度场驱动含水层中随机分布质点的运行轨迹,并以此推断和识别区域内地下水补给和排泄边界。在环境地学计算平台EnviFusion-CGS中实现,并构建了详细工作流程。以山东省青岛市大沽河中下游含水层为示范区,对含水系统的补给区和排泄区的空间分布及其动态变化进行了分析,取得了良好效果。本研究为构建和修订已有含水层概念模型提供了新思路。     

Assessment of groundwater recharge and water fluxes of the Guarani Aquifer System, Brazil

The groundwater recharge and water fluxes of the Guarani Aquifer System in the state of Sao Paulo in Brazil were assessed through a numeric model. The study area (6,748 km²) comprises Jacaré-Guaçú and Jacaré-Pepira River watersheds, tributaries of the Tietê River in the central region of the state. GIS based tools were used in the storage, processing and analysis of data. Main hydrologic phenomena were selected, leading to a groundwater conceptual model, taking into account the significant outcrops occurring in the study area. Six recharge zones were related to the geologic formation and structures of the semi-confined and phreatic aquifer. The model was calibrated against the baseflows and static water levels of the wells. The results emphasize the strong interaction of groundwater flows between watersheds and the groundwater inflow into the rivers. It has been concluded that lateral groundwater exchanges between basins, the deep discharges to the regional system, and well exploitation were not significant aquifer outflows when compared to the aquifer recharge. The results have shown that the inflows from the river into the aquifer are significant and have the utmost importance since the aquifer is potentially more vulnerable in these places.