Groundwater-flow modeling in the Yucatan karstic aquifer, Mexico

The current conceptual model of the unconfined karstic aquifer in the Yucatan Peninsula, Mexico, is that a fresh-water lens floats above denser saline water that penetrates more than 40 km inland. The transmissivity of the aquifer is very high so the hydraulic gradient is very low, ranging from 7–10 mm/km through most of the northern part of the peninsula. The computer modeling program AQUIFER was used to investigate the regional groundwater flow in the aquifer. The karstified zone was modeled using the assumption that it acts hydraulically similar to a granular, porous medium. As part of the calibration, the following hypotheses were tested: (1) karstic features play an important role in the groundwater-flow system; (2) a ring or belt of sinkholes in the area is a manifestation of a zone of high transmissivity that facilitates the channeling of groundwater toward the Gulf of Mexico; and (3) the geologic features in the southern part of Yucatan influence the groundwater-flow system. The model shows that the Sierrita de Ticul fault, in the southwestern part of the study area, acts as a flow barrier and head values decline toward the northeast. The modeling also shows that the regional flow-system dynamics have not been altered despite the large number of pumping wells because the volume of water pumped is small compared with the volume of recharge, and the well-developed karst system of the region has a very high hydraulic conductivity. Electronic Publication     

Impact of model complexity and multi-scale data integration on the estimation of hydrogeological parameters in a dual-porosity aquifer

This study investigates the impact of model complexity and multi-scale prior hydrogeological data on the interpretation of pumping test data in a dual-porosity aquifer (the Chalk aquifer in England, UK). In order to characterize the hydrogeological properties, different approaches ranging from a traditional analytical solution (Theis approach) to more sophisticated numerical models with automatically calibrated input parameters are applied. Comparisons of results from the different approaches show that neither traditional analytical solutions nor a numerical model assuming a homogenous and isotropic aquifer can adequately explain the observed drawdowns. A better reproduction of the observed drawdowns in all seven monitoring locations is instead achieved when medium and local-scale prior information about the vertical hydraulic conductivity (K) distribution is used to constrain the model calibration process. In particular, the integration of medium-scale vertical K variations based on flowmeter measurements lead to an improvement in the goodness-of-fit of the simulated drawdowns of about 30%. Further improvements (up to 70%) were observed when a simple upscaling approach was used to integrate small-scale K data to constrain the automatic calibration process of the numerical model. Although the analysis focuses on a specific case study, these results provide insights about the representativeness of the estimates of hydrogeological properties based on different interpretations of pumping test data, and promote the integration of multi-scale data for the characterization of heterogeneous aquifers in complex hydrogeological settings.     

Evaluation of alternative conceptual models for groundwater modelling

This study evaluates the alternative conceptual models for groundwater modelling. A true model was created with a synthetic alluvial fan-plain hydrogeological framework. Various alternative conceptual models were evaluated for groundwater flow simulations. The first alternative model is a single aquifer layer model; the second alternative model is a 3-layer aquifer model; and the third model is a 5-layer model consisting of 3 aquifers separated by 2 aquitards. All models could fit very well to the observations with optimized values of hydraulic conductivities. However, the single aquifer layer model can only compute water balance components with good accuracy. The 3-layer aquifer model can be used for water balance computation and groundwater head simulation with small errors. The 5-layer model is capable of simulating water budget, groundwater head distribution and travel times with high accuracy. Multi-model analysis found only the 3rd alternative model superior.     

Groundwater quality of porous aquifers in Greece: a synoptic review

Greece is dependent on groundwater resources for its water supply. The main aquifers are within carbonate rocks (karstic aquifers) and coarse grained Neogene and Quaternary deposits (porous aquifers). The use of groundwater resources has become particularly intensive in coastal areas during the last decades with the intense urbanization, tourist development and irrigated land expansion. Sources of groundwater pollution are the seawater intrusion due to over-exploitation of coastal aquifers, the fertilizers from agricultural activities and the disposal of untreated wastewater in torrents or in old pumping wells. In the last decades the total abstractions from coastal aquifers exceed the natural recharge; so the aquifer systems are not used safely. Over-exploitation causes a negative water balance, triggering seawater intrusion. Seawater intrusion phenomena are recorded in coastal aquifer systems. Nitrate pollution is the second major source of groundwater degradation in many areas in Greece. The high levels of nitrate are probably the result of over-fertilization and the lack of sewage systems in some urban areas.     

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.     

A method for the stochastic modeling of karstic systems accounting for geophysical data: an example of application in the region of Tulum, Yucatan Peninsula (Mexico)

The eastern coast of the Yucatan Peninsula, Mexico, contains one of the most developed karst systems in the world. This natural wonder is undergoing increasing pollution threat due to rapid economic development in the region of Tulum, together with a lack of wastewater treatment facilities. A preliminary numerical model has been developed to assess the vulnerability of the resource. Maps of explored caves have been completed using data from two airborne geophysical campaigns. These electromagnetic measurements allow for the mapping of unexplored karstic conduits. The completion of the network map is achieved through a stochastic pseudo-genetic karst simulator, previously developed but adapted as part of this study to account for the geophysical data. Together with the cave mapping by speleologists, the simulated networks are integrated into the finite-element flow-model mesh as pipe networks where turbulent flow is modeled. The calibration of the karstic network parameters (density, radius of the conduits) is conducted through a comparison with measured piezometric levels. Although the proposed model shows great uncertainty, it reproduces realistically the heterogeneous flow of the aquifer. Simulated velocities in conduits are greater than 1 cm?s^?1, suggesting that the reinjection of Tulum wastewater constitutes a pollution risk for the nearby ecosystems.     

Impact of groundwater withdrawals on the interaction of multi-layered aquifers in the Viterbo geothermal area (central Italy)

The impact of groundwater withdrawals on the interaction between multi-layered aquifers with different water qualities in the Viterbo geothermal area (central Italy) was studied. In this area, deep thermal waters are used to supply thermal spas and public pools. A shallow overlying aquifer carries cold and fresh water, used for irrigation and the local drinking-water supply. Starting with a conceptual hydrogeological model, two simplified numerical models were implemented: a steady-state flow model of the entire groundwater system, and a steady-state flow and heat transport model of a representative area, which included complex interactions between the aquifers. The impact of increased withdrawals associated with potential future development of the thermal aquifer must be considered in terms of the water temperature of the existing thermal sources. However, withdrawals from the shallow aquifer might also influence the discharge of thermal sources and quality of the water withdrawn from the shallow wells. The exploitation of the two aquifers is dependent on the hydraulic conductivity and thickness of the intervening aquitard, which maintains the delicate hydrogeological equilibrium. Effective methods to control this equilibrium include monitoring the vertical gradient between the two aquifers and the residual discharge of natural thermal springs.     

Modeling seawater intrusion in overexploited aquifers in the absence of sufficient data: application to the aquifer of Nea Moudania,northern Greece

In many coastal areas, overexploitation of groundwater resources has led both to the quantitative degradation of local aquifers and the deterioration of groundwater quality due to seawater intrusion. To investigate the behavior of coastal aquifers under these conditions, numerical modeling is usually implemented; however, the proper implementation of numerical models requires a large amount of data, which are often not available due to the time-consuming and costly process of obtaining them. In the present study, the investigation of the behavior of coastal aquifers under the lack of adequate data is attempted by developing a methodological framework consisting of a series of numerical simulations: a steady-state, a false-transient and a transient simulation. The sequence and the connection between these simulations constitute the backbone of the whole procedure aimed at adjusting the various model parameters, as well as obtaining the initial conditions for the transient simulation. The validity of the proposed methodology is tested through evaluation of the model calibration procedure and the estimation of the simulation errors (mean error, mean absolute error, root mean square error, mean relative error) using the case of Nea Moudania basin, northern Greece. Furthermore, a sensitivity analysis is performed in order to minimize the error estimates and thus to maximize the reliability of the models. The results of the whole procedure affirm the proper implementation of the developed methodology under specific conditions and assumptions due to the lack of sufficient data, while they give a clear picture of the aquifer’s quantitative and qualitative status.     

High‐resolution conceptual hydrogeological model of complex basaltic volcanic islands: a Mayotte,Comoros, case study

We present a high‐resolution conceptual hydrogeological model for complex basaltic volcanic islands based on Mayotte Island in the Comoros. Its geological structure and hydrogeological functioning are deduced from a large dataset: geological mapping, geophysics, some forty new boreholes, piezometric data, hydraulic conductivity, hydrochemical data, etc. We describe previously unknown deep cut‐and‐fill palaeovalleys. The resulting conceptual geological and hydrogeological model of the island is very different from the Hawaiian model, in that it lacks a low‐elevation basal aquifer and dyke‐impounded high‐level aquifers. It is closer to the Canary Islands model, which has, however, not yet been described at a high‐resolution scale. It does not have a continuous aquifer, but rather a discontinuous succession of perched aquifers separated by aquicludes and aquitards. This results more from the complex geological structure of the island, which has experienced several phases of volcanism, erosion and weathering, than from its age, but is also a result of the high‐resolution scale of the model. High‐resolution conceptual modelling is now necessary to solve problems of applied geology and hydrogeology.     

Effects of different boundary conditions on the simulation of groundwater flow in a multi-layered coastal aquifer system (Taranto Gulf,southern Italy)

The evaluation of the accuracy or reasonableness of numerical models of groundwater flow is a complex task, due to the uncertainties in hydrodynamic properties and boundary conditions and the scarcity of good-quality field data. To assess model reliability, different calibration techniques are joined to evaluate the effects of different kinds of boundary conditions on the groundwater flow in a coastal multi-layered aquifer in southern Italy. In particular, both direct and indirect approaches for inverse modeling were joined through the calibration of one of the most uncertain parameters, namely the hydraulic conductivity of the karst deep hydrostratigraphic unit. The methodology proposed here, and applied to a real case study, confirmed that the selection of boundary conditions is among the most critical and difficult aspects of the characterization of a groundwater system for conceptual analysis or numerical simulation. The practical tests conducted in this study show that incorrect specification of boundary conditions prevents an acceptable match between the model response to the hydraulic stresses and the behavior of the natural system. Such effects have a negative impact on the applicability of numerical modeling to simulate groundwater dynamics in complex hydrogeological situations. This is particularly important for management of the aquifer system investigated in this work, which represents the only available freshwater resource of the study area, and is threatened by overexploitation and saltwater intrusion.     

State of knowledge of coastal aquifer management in South America

A comparative analysis of the existing hydrogeological and management information from 15 coastal aquifers in South America was performed in order to obtain insight into common features of the sub-continent coastal zones. Some knowledge from other areas has been incorporated. There is a very variable degree of knowledge and management practice, ranging from almost no data and no action (the most common case), to sound conceptual models about aquifer behaviour and comprehensive management actions such as relocation of abstractions, pumping brackish groundwater, and aquifer vulnerability mapping. Some common features are: intensive groundwater exploitation; lack of characterization studies to support resource planning and management; lack of monitoring networks; and the need for raising awareness within society and its involvement in resource planning and management action programmes. Quality and quantity problems arising in heavily populated areas associated with coastal aquifers in South America point to unsustainable groundwater development. The sustainable use of those aquifers must rely on adequate evaluation of aquifer characteristics and monitoring.     

Assessing the saltwater remediation potential of a three-dimensional,heterogeneous, coastal aquifer system

This paper evaluates the remediation potential of a salinized coastal aquifer by utilizing a scenario simulation. Therefore, the numerical model OpenGeoSys is first validated against analytical and experimental data to represent transient groundwater level development and variable density saline intrusion. Afterwards, a regional scale model with a three-dimensional, heterogeneous hydrogeology is calibrated for a transient state and used to simulate a best-case scenario. Water balances are evaluated in both the transient calibration and scenario run. Visualization techniques help to assess the complex model output providing valuable insight in the occurring density-driven flow processes. Furthermore, modeling and visualization results give information on the time scale for remediation activities and, due to limitations in data quality and quantity reveal potential for model improvement.     

Review: Simulation-optimization models for the management and monitoring of coastal aquifers

The literature on the application of simulation-optimization approaches for management and monitoring of coastal aquifers is reviewed. Both sharp- and dispersive-interface modeling approaches have been applied in conjunction with optimization algorithms in the past to develop management solutions for saltwater intrusion. Simulation-optimization models based on sharp-interface approximation are often based on the Ghyben-Herzberg relationship and provide an efficient framework for preliminary designs of saltwater-intrusion management schemes. Models based on dispersive-interface numerical models have wider applicability but are challenged by the computational burden involved when applied in the simulation-optimization framework. The use of surrogate models to substitute the physically based model during optimization has been found to be successful in many cases. Scalability is still a challenge for the surrogate modeling approach as the computational advantage accrued is traded-off with the training time required for the surrogate models as the problem size increases. Few studies have attempted to solve stochastic coastal-aquifer management problems considering model prediction uncertainty. Approaches that have been reported in the wider groundwater management literature need to be extended and adapted to address the challenges posed by the stochastic coastal-aquifer management problem. Similarly, while abundant literature is available on simulation-optimization methods for the optimal design of groundwater monitoring networks, applications targeting coastal aquifer systems are rare. Methods to optimize compliance monitoring strategies for coastal aquifers need to be developed considering the importance of monitoring feedback information in improving the management strategies.     

The impact of salt diapirs on the quality of carbonate karst waters,Bastak, Iran

The Gavbast karstic aquifer located in southern Iran is in direct contact with an exposed salt diapir. To assess the influence of the diapir on the quality of groundwater in the karstic aquifer, electrical conductivity, total dissolved solids, flow rate, temperature and major ion concentrations were measured at 57 sampling sites, including springs, surface waters and wells. A conceptual model of groundwater flow is proposed for the Gavbast karstic aquifer based on the geological setting, water budget, local base of erosion, and hydrochemistry of the sampling sites. The model suggests two subbasins in the Gavbast Anticline draining into two distinct discharging alluvial sections. Unexpectedly, groundwater discharging from the carbonate Gavbast aquifer is saline or brackish and water is of chloride type. The study indicates that the source of salinity of the Gavbast aquifers is infiltration of surface diapir-derived brine into the aquifer. The contribution of the diapir brine in the Gavbast karst aquifer is calculated about 4 L/s, using chloride mass balance. Construction of salt basins to evaporate brine discharging from the diapir springs is proposed to reduce the salinity of karst water. A row of strategically placed wells in the Gavbast karst aquifer would potentially exploit large volumes of fresh groundwater before it is contaminated by the salt. Such low-cost remediation should allow the agricultural exploitation of 40 km² of currently barren land.     

Characterization of the groundwater flow system in the hillside and coastal aquifers of the Mersin-Tarsus region (Turkey)

In the region between Mersin and Tarsus cities, located along the Mediterranean Sea coast in southern Turkey, the demand for groundwater has increased dramatically as the available surface water supplies have already been developed. Fundamental information is required to characterize the existing groundwater system in this area in order to establish a sustainable groundwater-use policy. For this purpose, hydrochemical and environmental isotopic data were collected and integrated with available geological and hydrogeological information to develop a conceptual model of the system. Results, backed up mainly by depleted stable isotope composition and infinitesimal tritium content, suggest that most of the groundwater along the coastal zone is supplied by the neighboring mountain belt while local precipitation has also contributes to aquifer recharge. The validation of the conceptual perspective by a steady-state numerical groundwater flow model reveals that about 90% of the recharge to the aquifer system is supplied by the deep flow of karstic groundwater fed from the Taurus Mountains. Monitoring of changes in the recharge regime of the mountain sector seems to be critical in establishing future groundwater use policies.     

Development and parameterisation of a complex hydrogeological model based on high-resolution direct-push data

Ongoing developments in geological and hydrogeological investigation techniques, especially direct-push methods, have led to an increase in the quality, density and spatial resolution of data available from such investigations. This has created new challenges in the development of numerical models in terms of accurately and efficiently translating detailed and complex conceptual models into effective numerical models. Suitable geometrical and numerical modelling tools are essential in order to meet these challenges. This paper describes the development of a three-dimensional hydrogeological flow model for a contaminated site near Berlin, Germany, based on high-resolution geological data obtained principally using direct-push methods. The available data were first interpreted to construct a detailed GIS-based geological model, which formed the basis of the conceptual site model. The conceptual model was then translated into a geometrical model, which was used to create a finite element numerical model. An innovative geometry object-based approach enabled the complex structural details of the conceptual model to be accurately reproduced in the numerical model domain. The resulting three-dimensional steady-state unconfined flow model was successfully calibrated using external automated calibration software, whereby parameter values for groundwater recharge and hydraulic conductivity were determined.     

Tidal controls on coastal groundwater conditions: field investigation of a macrotidal system

Previous studies of coastal groundwater resources have not properly explored the hydraulic conditions at coastal and estuarine boundaries, despite recognised influences of oceanic and estuarine water-level fluctuations (e.g. tides and waves) on groundwater near these boundaries. Such influences may have important implications for determining submarine groundwater discharge and seawater intrusion. In this paper, oceanic and estuarine controls on the hydrology of coastal aquifers are characterised for a macrotidal system—the Pioneer Valley coastal plain, northeastern Australia. The tidal water-table over-height (tide-induced increase in average water-table height at the coastline) is quantified at three locations and compared with theoretical estimates, which assume simplified physiographical conditions compared with those encountered at the field sites. The results indicate that local geological conditions, beach morphology and characteristics of tidal forcing control the behaviour of nearshore groundwater within the system. Existing analytical and numerical solutions that are commonly applied as first-pass estimates are found to be insufficient for predicting observed tidal water-table over-height in the Pioneer Valley, due to the sediment heterogeneities, non-uniform beach slopes and large tidal ranges of the system. The study reveals the spatial and temporal variability in tidally influenced hydraulic heads at the estuarine and coastal boundaries of the aquifer, and provides estimates of tidal water-table over-height up to 2.41 m during spring tides. These findings highlight the complexity of coastal groundwater systems, and the need to incorporate appropriate nearshore and near-estuary boundary conditions in models of regional groundwater flow in coastal aquifers.     

The importance of alternative conceptual models for simulation of concentrations in a multi-aquifer system

Four different conceptual models based on alternative geological interpretations were formulated for a shallow 600 km² aquifer system in Denmark comprising Quaternary deposits. Each of the four models was calibrated against groundwater heads and discharge measurements through inverse modeling. Subsequently, the transport capabilities of the four models were compared to 32 concentration measurements of environmental tracers (tritium ³H, helium-3 ³He, chlorofluorocarbons CFC11, CFC12 and CFC113). The flow simulations showed only minor differences in spatial head distribution associated with alternative conceptualizations despite the complexity of the aquifer system and the significant differences in geological interpretations. The models, however, showed major differences in predictions of the age of the groundwater and environmental tracer concentrations, differences that are seen as an effect of model structure uncertainty, because no additional calibrations to these data were performed. A single conceptualization may be adequate in characterizing the natural behavior of a field system after calibration, because the calibration procedure is able to compensate for errors in the data or in the conceptual model through biased parameter values. However, once extrapolation beyond the calibration base is attempted, different conceptual model formulations result in significantly different results. Consequently, it is crucial to take model conceptual uncertainty into account when making predictions beyond the calibration base.     

Compilation simulation of surface water and groundwater resources using the SWAT-MODFLOW model for a karstic basin in Iran

Modeling of karstic basins can provide a better understanding of the interactions between surface water and groundwater, a more accurate estimation of infiltrated water amount, and a more reliable water balance calculation. In this study, the hydrological simulation of a karstic basin in a semiarid region in Iran was performed in three different stages. In the first stage, the original SWAT model was used to simulate surface-water flow. Then, the SWAT-MODFLOW conjunctive model was implemented according to the groundwater characteristics of the study area. Finally, due to the karstic characteristics of the region and using the CrackFlow (CF) package, the SWAT-MODFLOW-CF conjunctive model was developed to improve the simulation results. The coefficient of determination (R²) and the Nash-Sutcliffe efficiency coefficient (NSE) as error evaluation criteria were calculated for the models, and their average values were 0.63 and 0.57 for SWAT, 0.68 and 0.61 for SWAT-MODFLOW, 0.73 and 0.7 for SWAT-MODFLOW-CF, respectively. Moreover, the mean absolute error (MAE) and root mean squared error (RMSE) of the calibration for groundwater simulation using the SWAT-MODFLOW model were 1.23 and 1.77 m, respectively. These values were 1.01 and 1.33 m after the calibration of the SWAT-MODFLOW-CF model. After modifying the CF code and keeping the seams and cracks open in both dry and wet conditions, the amount of infiltrated water increased and the aquifer water level rose. Therefore, the SWAT-MODFLOW-CF conjunctive model can be proposed for use in karstic areas containing a considerable amount of both surface water and groundwater resources.

     

基于广义径向流模型的非均质孔隙含水层井流试验分析

Theis模型、Dupuit模型等经典井流模型分析非均质含水层井流试验数据有一定的局限性,获取的参数不适合表征非均质含水层特性。而GRF模型可以获取含水层流动特性的数据。相比于Theis模型、Dupuit模型,GRF模型更能表征非均质含水层特性。以黄石东湖新村棋盘洲长江大桥的抽水试验数据为例,采用Theis模型和GRF模型计算含水层渗透系数,结合实际水文地质条件,对比分析不同方法计算的水文地质参数,并计算含水层水流维数和表观压力传导系数（K_f /S_sf）。结果表明:研究区含水层为细砂夹条带状黏土透镜体的非均质含水层,采用GRF模型计算结果更符合实际情况,渗透系数为 4.09×10?3cm/s;含水层水流维数为1.61,地下水为双线性流动状态,含水层对抽水试验的响应主要受黏土条带控制;观测井和抽水主井的K_f /S_sf呈非线性相关,进一步验证了含水层的非均质性。在非均质孔隙含水层中,应用多孔联合非稳定GRF井流试验方法不仅能确定水文地质参数,并且能丰富对含水层特性的认知。