Hydrogeochemical modeling for groundwater management in arid and semiarid regions using MODFLOW and MT3DMS: A case study of the Jeffara of Medenine coastal aquifer,South-Eastern Tunisia

The study of water quality and the quantification of reserves and their variations according to natural and anthropogenic forcing is necessary to establish an adequate management plan for groundwater resources. For this purpose, a modeling approach is a useful tool that allows, after calibration phase and verification of simulation, and under different scenarios of forcing and operational changes, to estimate and control the groundwater quantity and quality. The main objective of this study is to collect all available data in a model that simulates the Jeffara of Medenine coastal aquifer system functioning. To achieve this goal, a conceptual model was constructed based on previous studies and hydrogeological investigations. The regional groundwater numerical flow model for the Jeffara aquifer was developed using MODFLOW working under steady-state and transient conditions. Groundwater elevations measured from the piezometric wells distributed throughout the study area in 1973 were selected as the target water levels for steady state (head) model calibration. A transient simulation was undertaken for the 42 years from 1973 to 2015. The historical transient model calibration was satisfactory, consistent with the continuous piezometric decline in response to the increase in groundwater abstraction. The developed numerical model was used to study the system's behavior over the next 35 years under various constraints. Two scenarios for potential groundwater extraction for the period 2015–2050 are presented. The predictive simulations show the effect of the increase of the exploitation on the piezometric levels. To study the phenomenon of salinization, which is one of the most severe and widespread groundwater contamination problems, especially in coastal regions, a solute transport model has been constructed by using MT3DMS software coupled with the groundwater flow model. The best calibration results are obtained when the connection with the overlying superficial aquifer is considered suggesting that groundwater contamination originates from this aquifer. Recommendations for water resource managers

• The results of this study show that Groundwater resources of Jeffara of Medenine coastal aquifer in Tunisia are under immense pressure from multiple stresses.
• The water resources manager must consider the impact of economic and demographic development in groundwater management to avoid the intrusion of saline water.
• The results obtained presented some reference information that can serve as a basis for water resources planning.
• The model runs to provide information that managers can use to regulate and adequately control the Jeffara of Medenine water resources.

     

Transient flow of slightly compressible fluids through double-porosity,double-permeability systems — A state-of-the-art review

The theory of transient flow of slightly compressible fluids through naturally fractured reservoirs based on the double porosity conceptualization is summarized. The main achievements in the theory of fluid flow in leaky aquifer systems which are closely related with the double-porosity, double-permeability problems are also addressed. The main emphasis of this review is the analytical treatment of these problems.     

Simulation of immiscible multiphase flow in porous media: A focus on the capillary fringe of oil-contaminated aquifers

This paper deals with the analysis of some aspects of the vertical and lateral migration of oil spills in the unsaturated and the capillary zone of a phreatic aquifer. Our motivation stems from the fact that such contamination represents a severe danger for ground-water resources all over the world and from the present acute problem of jet-fuel contamination in some location of Israel. In the present study, we shall focus our efforts on the analysis of the upper layers of the aquifer which are often subjected to the most significant oil contamination. Neglecting coupled processes effects such as dilution, adsorption and volatilization, also adopting Richard's assumption, a three-phase flow model is introduced with capillary heads of the water and the oil as variables. The resulting model which is coupled and strongly non-linear is solved using a vertical two-dimensional Finite-Element procedure together with a quasi-Newton optimization algorithm. Applying that scheme, various scenarios of oil migration in the unsaturated and the capillary zone were simulated. Some migration characteristics prediced by the numerical simulations are discussed. In particular, the dynamics of the water and oil phases during the migration process is discussed.     

Modelling of flow pattern in a shallow aquifer affected by reservoirs

A compartmental model is developed to estimate flow parameters of a shallow aquifer affected by water loads in surface reservoirs and to evaluate its nonsteady flow distribution. The method incorporates temporal piezometric head measurements and sampling of water for dissolved chemicals and isotopes analyses. Each compartment is governed by a set of equations describing the conservation of linear momentum and mass balance expressions for water, isotopes, and dissolved chemicals. The number of compartmental balance expressions always must be greater than that of the unknown flow parameters associated with each compartment. An optimization method is described to yield spatial distribution of aquifer storativity, transmissivity, porosity, leakage, and compliance coefficients and fluxes leaking into an aquifer's lower boundary. Future predictions of an aquifer's piezometric head distribution in a compartmental system is formulated on the basis of the estimated flow parameters and the leakage components.Compartmental modelling which incorporates concentrations of environmental tracers, may yield efficiency in computing resources and accuracy enhancement for predicting an aquifer's flow regime.     

Stochastic Analysis of Macrodispersion in a Semi‐Confined Aquifer

In this paper, the macroscopic dispersion resulting from one and twodimensional flows through a semiconfined aquifer with spatially variable hydraulic conductivity K which is represented by a stationary (statistically homogeneous) random process is analyzed using the spectral representation technique. Stochastic fluctuation equations of the steady flow and solute transport are solved to construct the macroscopic dispersive flux and evaluate the resulting macrodispersivity tensor in terms of the leakage factor and input covariances describing the hydraulic conductivity in a semiconfined aquifer bounded by a leaky layer above and an impervious stratum below. The macrodispersivity tensor is studied using some convenient forms of the log hydraulic conductivity process. The sensitivity of the resulting macrodispersivity to the input covariances is discussed along with the influence of the leakage factor for both one and twodimensional flows. It is found that the longitudinal macrodispersivities are increased due to the presence of leakage, while the transverse macrodispersivities are reduced due to leakage.     

The Kinematic Flow Structure for the Gvirtzman–Gorelick In Situ VOC Remediation System

Potential theory and Stokes' stream function techniques are used to investigate the flow structure around the recirculation system developed by Gvirtzman and Gorelick (1992, 1993), which consists of an extraction well and a gallery (trench) for the recharge of treated water to the aquifer. Analytical formulas are derived for the drawdown, velocity, and stream function for a model in which the extraction well is modeled as a uniformly distributed line sink and the gallery is modeled as a uniformly distributed ring source. Travel times are reported for water particles traveling along the streamlines containing 50 and 90% of the flow for various degrees of well penetration and various radii of the ring source. The travel times along the streamline resulting in the shortest travel time (not necessarily the shortest path) are also reported for various degrees of well penetration and various radii of the ring source. The method completely eliminates the use of numerical finite-difference or finite-element methods and can be used for optimization of technological parameters of this remediation system.     

Estimating the spatial distribution of artificial groundwater recharge using multiple tracers

Stable isotopes of water, organic micropollutants and hydrochemistry data are powerful tools for identifying different water types in areas where knowledge of the spatial distribution of different groundwater is critical for water resource management. An important question is how the assessments change if only one or a subset of these tracers is used. In this study, we estimate spatial artificial infiltration along an infiltration system with stage–discharge relationships and classify different water types based on the mentioned hydrochemistry data for a drinking water production area in Switzerland. Managed aquifer recharge via surface water that feeds into the aquifer creates a hydraulic barrier between contaminated groundwater and drinking water wells. We systematically compare the information from the aforementioned tracers and illustrate differences in distribution and mixing ratios. Despite uncertainties in the mixing ratios, we found that the overall spatial distribution of artificial infiltration is very similar for all the tracers. The highest infiltration occurred in the eastern part of the infiltration system, whereas infiltration in the western part was the lowest. More balanced infiltration within the infiltration system could cause the elevated groundwater mound to be distributed more evenly, preventing the natural inflow of contaminated groundwater.

Dedicated to Professor Peter Fritz on the occasion of his 80th birthday     

Genetic Algorithms for Optimizing the Remediation of Contaminated Aquifer

Produced water constitutes a large amount of waste fluids during the production operation of an oil field. Underground injection for disposing the wastewater from hydrocarbon production is an engineering problem due to the possibility of leakage of injected pollutant material from receiving medium to a drinking water source. This paper describes a method for optimization of polluted aquifer remediation design using one of the artificial intelligence optimization methods, namely Genetic Algorithms (GAs). As a case study, the contaminated area was created by using a groundwater transport simulator, which is based on Method of Characteristics (MOC). Then, the developed computer program was run to find the optimum solution for remediation, and the solution yielded from the program was verified by using a groundwater simulator. The plume was captured and the concentration level of chloride ion within the aquifer was diminished by using extraction wells. The analytical model approach provided different alternatives for appropriate isolation of plume. GAs were used as an optimization technique for making a decision among the alternatives, by considering operation time, number of wells, pumping rate and drawdown as decision variables and constraints.     

垂直与水平渗透作用下潜水非稳定渗流运动规律

对河渠边界控制的半无限含水层，建立垂向入渗与河渠水平渗透共同作用下的潜水非稳定渗流模型；利用Boussinesq第一线性化方法，通过Laplace变换，给出模型的解析解．证明相关经典公式与模型特定解之间的转换关系，分析经典公式适用范围．根据模型解，逐一定量研究下述变量，如垂向入渗强度、河渠水位变动幅度、含水层结构参数如给水度和导压系数、计算点与边界之间的距离，对渗流过程的影响．这些变量的变化，对潜水位获最大上升速度的时间产生延迟效应；论证一些变量间产生等效延迟效应的条件．根据解的数学特征，讨论其对应的物理意义和潜水位变动规律．     

时变垂向入渗影响下河渠-潜水非稳定渗流模型的解及应用

在河渠边界控制下的半无限含水层中,对时变垂向入渗影响下的潜水非稳定渗流模型,利用Boussinesq第一线性化方法,通过Laplace变换并结合卷积原理,导出模型的解析解.根据不同水文地质条件下解的数学特征,建立相应的含水层参数求解方法;在此基础上,建立计算河渠与含水层之间水量交换的公式,以及计算潜水蒸发强度的递推公式.以安徽淮北平原某河流-潜水含水层系统为例,阐述上述方法的计算过程与步骤.