Evaluation of a capture zone overlay method for designing groundwater remediation systems

 Simplified approaches are often used to model the removal of groundwater contamination. These approaches can yield poor remediation schemes because they incorrectly portray the effects of multiple pumping wells. In this study, a pumping configuration designed by graphically overlaying capture zones having an identical, quasi-elliptical shape was evaluated with a numerical mass transport model. After a 3-year period (within which the hypothetical aquifer was to be remediated) the contaminant mass had been reduced by 77%. Due to stagnation zones which developed between extraction wells, approximately 15 years of pumping was required to remediate the aquifer with the overlay configuration. An alternative design, consisting of an extraction well between two injection wells along the long axis of the plume, removed the contaminant within the 3-year design period. Received: 23 October 1995 · Accepted: 18 June 1996     

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.     

Temporal trends in concentrations of DBCP and nitrate in groundwater in the eastern San Joaquin Valley,California, USA

Temporal monitoring of the pesticide 1,2-dibromo-3-chloropropane (DBCP) and nitrate and indicators of mean groundwater age were used to evaluate the transport and fate of agricultural chemicals in groundwater and to predict the long-term effects in the regional aquifer system in the eastern San Joaquin Valley, California. Twenty monitoring wells were installed on a transect along an approximate groundwater flow path. Concentrations of DBCP and nitrate in the wells were compared to concentrations in regional areal monitoring networks. DBCP persists at concentrations above the US Environmental Protection Agency’s maximum contaminant level (MCL) at depths of nearly 40 m below the water table, more than 25 years after it was banned. Nitrate concentrations above the MCL reached depths of more than 20 m below the water table. Because of the intensive pumping and irrigation recharge, vertical flow paths are dominant. High concentrations (above MCLs) in the shallow part of the regional aquifer system will likely move deeper in the system, affecting both domestic and public-supply wells. The large fraction of old water (unaffected by agricultural chemicals) in deep monitoring wells suggests that it could take decades for concentrations to reach MCLs in deep, long-screened public-supply wells, however.     

Modeling approaches and strategies for data-scarce aquifers: example of the Dar es Salaam aquifer in Tanzania

Management of groundwater resources can be improved by using groundwater models to perform risk analyses and to improve development strategies, but a lack of extensive basic data often limits the implementation of sophisticated models. Dar es Salaam in Tanzania is an example of a city where increasing groundwater use in a Pleistocene aquifer is causing groundwater-related problems such as saline intrusion along the coastline, lowering of water-table levels, and contamination of pumping wells. The lack of a water-level monitoring network introduces a problem for basic data collection and model calibration and validation. As a replacement, local water-supply wells were used for measuring groundwater depth, and well-top heights were estimated from a regional digital elevation model to recalculate water depths to hydraulic heads. These were used to draw a regional piezometric map. Hydraulic parameters were estimated from short-time pumping tests in the local wells, but variation in hydraulic conductivity was attributed to uncertainty in well characteristics (information often unavailable) and not to aquifer heterogeneity. A MODFLOW model was calibrated with a homogeneous hydraulic conductivity field and a sensitivity analysis between the conductivity and aquifer recharge showed that average annual recharge will likely be in the range 80–100 mm/year.     

Delineation of groundwater protection zones based on tracer tests and transport modeling in alluvial sediments

 Regulations aiming to protect exploitable groundwater resources were edicted in Belgium a few years ago. Therefore, prevention and protection zones are defined by law and must be determined practically around each pumping well or spring, based on local hydrogeological conditions. The determination of hydrodynamic and hydrodispersive parameters, characterizing the local flow and transport properties of the aquifer, requires pumping and tracing tests. The interpretation of these field experiments, considering the heterogeneity of the geological layers, is performed through the use of numerical FEM simulations of the groundwater flow and pollutant transport conditions in a deterministic framework. After calibration of the model on experimental measurements, multiple simulations with contaminant injections at various points of the modeled domain allow the determination of the transfer time of the pollutant in the studied aquifer whilst taking the updated heterogeneity into account. On the basis of the computed transfer times in the saturated zone, the various prevention and protection areas can be assessed based on provisions of the law. Received: 27 June 1997 · Accepted: 29 July 1997     

Contamination,risk, and heterogeneity: on the effectiveness of aquifer remediation

The effectiveness of aquifer remediation is typically expressed in terms of a reduction in contaminant concentrations relative to a regulated maximum contaminant level (MCL), and is usually confirmed by sparse monitoring data and/or simple model calculations. Here, the effectiveness of remediation is re-examined from a more thorough risk-based perspective that goes beyond the traditional MCL concept. A methodology is employed to evaluate the health risk to individuals exposed to contaminated household water that is produced from groundwater. This approach explicitly accounts for differences in risk arising from variability in individual physiology and water use, the uncertainty in estimating chemical carcinogenesis for different individuals, and the uncertainties and variability in contaminant concentrations within groundwater as affected by transport through heterogeneous geologic media. A hypothetical contamination scenario is developed as a case study in a saturated, alluvial aquifer underlying an actual Superfund site. A baseline (unremediated) human exposure and health risk scenario, as induced by contaminated groundwater pumped from this site, is predicted and compared with a similar estimate based upon pump-and-treat exposure intervention. The predicted reduction in risk in the remediation scenario is not an equitable one—that is, it is not uniform to all individuals within a population and varies according to the level of uncertainty in prediction. The importance of understanding the detailed hydrogeologic connections that are established in the heterogeneous geologic regime between the contaminated source, municipal receptors, and remediation wells, and its relationship to this uncertainty is demonstrated. Using two alternative pumping rates, we develop cost-benefit curves based upon reduced exposure and risk to different individuals within the population, under the presence of uncertainty.     

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.     

Aquifer contamination by chlorinated-VOCs: the case of an urban metropolis megasite overlying the Coastal Plain aquifer in Israel

The Tel Aviv metropolitan region (200 km²), located on the sandy and phreatic Mediterranean Coastal Plain aquifer in Israel, is an example of a contaminated megasite that is additionally complicated by being part of a densely populated urban environment. Despite being a major source of fresh potable water, no dedicated aquifer monitoring systems exist. The vast majority of tested water supply wells (82% of 90 wells) were found to be contaminated with one or more chlorinated volatile organic compounds (Cl-VOCs) such as trichloroethene and tetrachloroethene, originating from multiple sources, including former industrial sites. Temporal variations in Cl-VOC concentrations in water supply wells were large and unpredictable. Such aquifer contamination requires immediate action in terms of aquifer management and municipal water distribution. To cope with temporal and spatial variations in contaminant concentrations in water supply wells in an area where monitoring wells are absent, aquifer impact areas were defined based on the concentration of the Cl-VOC contaminant in pumping wells that was greatest in relation to its drinking water standard over a 3-year period. Such a map can be used to define water supply treatment and municipal well monitoring requirements, until an adequate monitoring system is established.     

某氨氮污染地下水体抽出 -处理系统优化模拟研究

抽出 -处理系统设计多侧重于考虑修复初期的效率,在修复后期通常效率低下,产生拖尾现象,其优化的关键在于布设的井群系统能否高效抽出受污染的地下水体。利用溶质运移数值模拟可为井群布设和抽水方案优化提供依据。本研究旨在优化我国北方某化肥厂高浓度氨氮污染的地下水体的抽出 -处理修复系统,节约时间和成本。在水文地质调查及氨氮浓度监测的基础上,综合考虑井数、抽水天数和总抽水量三个变量,采用中轴线法与三角形法结合的布井方法,利用GMS软件反复试算,筛选出三种较优抽水方案并进一步模拟优化,最终从中选出最优抽水方案。结果,相比最初方案（方案1）,最优方案（方案3）将修复周期缩短了23个月,抽水总量减少了约31.9×104 m3,而抽水井数量仅增加了1口。该模型进行了稳定流水位拟合验证和4期非稳定流实测溶质浓度验证,较符合实际。结果表明,针对抽水井数量不足引起的拖尾问题,关键因素在于合理的井位布设与分阶段的抽水模式。在修复过程中,及时对地下水中污染物进行监测,并随着污染羽变化过程及时调整抽水方案,保证高浓度区一直有抽水井进行较大流量抽水,可有效提高修复效率并缩短修复周期。     

Conceptual model of a coastal aquifer system in northern Greece and assessment of saline vulnerability due to seawater intrusion conditions

This paper refers to the development of a conceptual model for the management of a coastal aquifer in northern Greece. The research presents the interpretation and analysis of the quantitative (groundwater level recordings and design of piezometric maps) regime and the formation of the upcone within the area of investigation. Additionally it provides the elaboration of the results of chemical analyses of groundwater samples (physicochemical parameters, major chemical constituents and heavy metals and trace elements) of the area which were taken in three successive irrigation periods (July–August 2003, July–August 2004 and July 2005), in order to identify areas of aquifer vulnerability. The study identifies the areas where ion exchange phenomena occur, as well as the parts of the aquifer where the qualitative degradation of the aquifer system is enhanced. The paper, finally, assesses the lack of any scientific groundwater resources management of the area by the local water authorities, as well as the current practices of the existing pumping conditions scheme as applied by groundwater users.An erratum to this article can be found at     

Hypothesis testing of buoyant plume migration using a highly parameterized variable-density groundwater model at a site in Florida, USA

A highly parameterized variable-density groundwater flow and solute transport model was developed to test multiple hypotheses for upward movement of treated wastewater (effluent) injected into a saline coastal aquifer in southeastern Florida, USA. The model was designed to assess risk to a drinking-water aquifer above the zone of injection, where monitoring wells have detected effluent. The model-based analysis accommodated geological and data complexity, including the observed presence of effluent in upper monitoring wells, but not in lower monitoring wells, thereby giving the appearance of the effluent having bypassed geological layers. The modeling approach included the application of multiple methodologies to reduce model run times during parameter estimation while providing detailed calibrated model(s) that can be used to assess the potential capacity for different mechanisms of effluent migration. The methods included use of a semi-analytical equation to quickly calculate initial concentrations, parallelization of model runs over multiple processors when calibrating, and utilization of the concepts of singular value decomposition and Tikhonov regularization to accommodate a high level of parameterization complexity. The results reveal that vertical effluent migration could occur as diffuse flow through heterogeneous confining units; however, flow through a channelized pathway caused by well construction appears to be more likely.     

Sensitivity Analysis of a Combined Groundwater Flow and Solute Transport Model Using Local-Grid Refinement: A Case Study

Combining groundwater flow models with solute transport models represents a common challenge in groundwater resources assessments and contaminant transport modeling. Groundwater flow models are usually constructed at somewhat larger scales (involving a coarser discretization) to include natural boundary conditions. They are commonly calibrated using observed groundwater levels and flows (if available). The groundwater solute transport models may be constructed at a smaller scale with finer discretization than the flow models in order to accurately delineate the solute source and the modeled target, to capture any heterogeneity that may affect contaminant migration, and to minimize numerical dispersion while still maintaining a reasonable computing time. The solution that is explored here is based on defining a finer grid subdomain within a larger coarser domain. The local-grid refinement (LGR) implemented in the Modular 3D finite-difference ground-water flow model (MODFLOW) code has such a provision to simulate groundwater flow in two nested grids: a higher-resolution sub-grid within a coarse grid. Under the premise that the interface between both models was well defined, a comprehensive sensitivity and uncertainty analysis was performed whereby the effect of a parameter perturbation in a coarser-grid model on transport predictions using a higher-resolution grid was quantified. This approach was tested for a groundwater flow and solute transport analysis in support of a safety evaluation of the future Belgian near-surface radioactive waste disposal facility. Our reference coarse-grid groundwater flow model was coupled with a smaller fine sub-grid model in two different ways. While the reference flow model was calibrated using observed groundwater levels at a scale commensurate with that of the coarse-grid model, the fine sub-grid model was used to run a solute transport simulation quantifying concentrations in a hypothetical well nearby the disposal facility. When LGR coupling was compared to a one-way coupling, LGR was found to provide a smoother flow solution resulting in a more CPU-efficient transport solution. Parameter sensitivities performed with the groundwater flow model resulted in sensitivities at the head observation locations. These sensitivities identified the recharge as the most sensitive parameter, with the hydraulic conductivity of the upper aquifer as the second most sensitive parameter in regard to calculated groundwater heads. Based on one-percent sensitivity maps, the spatial distribution of the observations with the highest sensitivities is slightly different for the upper aquifer hydraulic conductivity than for recharge. Sensitivity analyses were further performed to assess the prediction scaled sensitivities for hypothetical contaminant concentrations using the combined groundwater flow and solute transport models. Including all pertinent parameters into the sensitivity analysis identified the hydraulic conductivity of the upper aquifer as the most sensitive parameter with regard to the prediction of contaminant concentrations.     

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.     

Delaunay方法的改进及其在地下水污染监测网设计中的应用

地下水污染监测网的设计包括取样点在空间上的采样位置和时间上的取样频率这两方面的确定,其目的是为了准确刻画污染羽在含水层中随时间的变化状况。介绍了一种新的基于Delaunay三角形剖分方法用于地下水污染监测网在空间上的去冗余设计。同时,针对Delaunay方法没有考虑含水层参数变化和不能新增必要监测点的缺陷,提出了改进措施。改进的Delaunay方法可以直接与污染物的运移模型耦合,通过获取的含水层参数能够比较好地检验监测网在空间上的去冗余设计效果。实例研究表明,该方法在保证监测精度的基础上,能够找到最优的空间取样位置点,同时与改进前相比,利用该方法得到的地下水污染监测网能更加客观地反映污染物在含水层中的空间分布特征。     

水力层析法与克立金法估算非均质含水层渗透系数场比较

含水层非均质性空间分布特征的识别,是准确模拟地下水流和污染物运移的前提.基于室内非均质含水层砂箱实验,分别利用水力层析法和克立金插值法刻画了非均质含水层渗透系数场.研究结果表明:(1) 水力层析法与克立金法相比,不仅可以更好地刻画非均质含水层渗透系数场,还可以更高精度地预测地下水流过程;(2) 水力层析抽水实验中,通过增加抽水实验组数可以有效地提高水力层析参数反演的精度,但是抽水实验组数增加到一定程度以后,再增加抽水组数不会显著提升参数反演的效果.后续需要进一步研究水力层析抽水实验合适的组数,进一步对抽水井进行优化布设.      

Conceptual and numerical models for sustainable groundwater management in the Thaphra area, Chi River Basin, Thailand

Sustainable management of groundwater resources is vital for development of areas at risk from water-resource over-exploitation. In northeast Thailand, the Phu Thok aquifer is an important water source, particularly in the Thaphra area, where increased groundwater withdrawals may result in water-level decline and saline-water upconing. Three-dimensional finite-difference flow models were developed with MODFLOW to predict the impacts of future pumping on hydraulic heads. Four scenarios of pumping and recharge were defined to evaluate the system response to future usage and climate conditions. Primary model simulations show that groundwater heads will continue to decrease by 4–12?m by the year 2040 at the center of the highly exploited area, under conditions of both increasing pumping and drought. To quantify predictive uncertainty in these estimates, in addition to the primary conceptual model, three alternative conceptual models were used in the simulation of sustainable yields. These alternative models show that, for this case study, a reasonable degree of uncertainty in hydrostratigraphic interpretation is more impactful than uncertainty in recharge distribution or boundary conditions. The uncertainty-analysis results strongly support addressing conceptual-model uncertainty in the practice of groundwater-management modeling. Doing so will better assist decision makers in selecting and implementing robust sustainable strategies.     

Modeling of seawater intrusion in a coastal aquifer of Karaburun Peninsula,western Turkey

Seawater intrusion is a major problem to freshwater resources especially in coastal areas where fresh groundwater is surrounded and could be easily influenced by seawater. This study presents the development of a conceptual and numerical model for the coastal aquifer of Karareis region (Karaburun Peninsula) in the western part of Turkey. The study also presents the interpretation and the analysis of the time series data of groundwater levels recorded by data loggers. The SEAWAT model is used in this study to solve the density-dependent flow field and seawater intrusion in the coastal aquifer that is under excessive pumping particularly during summer months. The model was calibrated using the average values of a 1-year dataset and further verified by the average values of another year. Five potential scenarios were analyzed to understand the effects of pumping and climate change on groundwater levels and the extent of seawater intrusion in the next 10 years. The result of the analysis demonstrated high levels of electrical conductivity and chloride along the coastal part of the study area. As a result of the numerical model, seawater intrusion is simulated to move about 420 m toward the land in the next 10 years under “increased pumping” scenario, while a slight change in water level and TDS concentrations was observed in “climate change” scenario. Results also revealed that a reduction in the pumping rate from Karareis wells will be necessary to protect fresh groundwater from contamination by seawater.     

The dependence of stream depletion by seasonal pumping on various hydraulic characteristics and engineering factors

Compensation pumping is used to alleviate deficiencies in streamflow discharge during dry seasons. Short-term groundwater pumping can use aquifer storage instead of catchment-zone water until the drawdown reaches the edge of the stream. The capacitance is a complex, dimensionless parameter of an aquifer system that defines the delayed effect on streamflow when there is groundwater pumping. This parameter is a function of aquifer hydraulic characteristics, pumping time, and distance between the well and stream edge; the latter can involve stream leakance and vertical leakance of an associated aquitard. Three typical hydraulic cases of combined water systems (major catchment-zone wells close to the stream and compensation pumping wells) were classified depending on their capacitance structure (i.e. the relationship between surface water and groundwater): (1) perfect hydraulic connection between the stream and aquifer; (2) imperfect hydraulic connection between the stream and aquifer; and (3) essentially imperfect hydraulic connection between the stream and the underlying confined aquifer. The impact of various hydraulic characteristics and engineering factors on stream depletion was examined by conceptual and numerical modeling. To predict the suitability and efficiency of a combined water system application, regression tests were undertaken on unit stream depletion and capacitance, and power dependencies were defined.     

Case history of one of the few successful Superfund remediation sites: A site at Salinas, California, USA

In 1985, a former tire manufacturing plant surrounded by agricultural fields in the Salinas Valley was designated a Superfund Site by the US Environmental Protection Agency. The plant had been operating for seventeen years, from 1963 until 1980. When dismantling of the plant was started, it was determined that toxic hydrocarbon solvents and oils from the plant had contaminated soil and groundwater in alluvial deposits alongside the plant. It was determined later that the groundwater contamination also lay beneath the agricultural fields in a narrow groundwater plume that extends about 4.3 km downgradient from the plant. Because of the complex architecture of the aquifer system, the gradient, and extensive pumping of agricultural wells, the contaminants migrated northwestward and downward to deeper levels away from the plant.

The agricultural fields are underlain by an unconfined shallow aquifer and by a system of confined aquifers that extend to more than 180 m below surface. Aquifers are discontinuous beds of channel sand and gravel; confining beds are overbank clay and silt, and estuarine clay.

Geophysical data, logs of existing agricultural and other wells, and careful consideration of the stratigraphic architecture of the depositional environment provided the basis for a conceptual hydrogeologic model and for locating characterization wells for detailed visual and geophysical logging and hydrologic testing. Successive refinements of the characterization by sequential installation of wells indicated optimal locations for installation of extraction and monitoring wells. Validity of the concept of the hydrogeologic regime was verified by close match of predictions made by modeling with the later results of pumping from the extraction wells in a pump-and-treat system.

Successful remediation was accomplished by analyzing data from 110 agricultural wells, the few domestic water wells, nearly 200 sequentially installed stratigraphic-characterization and monitoring wells, 25 extraction wells, and by close cooperation among federal, state and local agencies, and the ranchers and growers.

Total contaminants recovered from activated-charcoal strippers of the treatment system totalled < 230 kg. Large quantities were harmlessly volatilized and dispersed into the atmosphere by air strippers and by agricultural sprinkling systems spraying water onto the fields. Crop testing showed no contamination of food crops. The activity has taken seven years and has cost more than US$22 million.     

Evaluation of a pumping test of the Snake River Plain aquifer using axial-flow numerical modeling

The Snake River Plain aquifer in southeast Idaho is hosted in a thick sequence of layered basalts and interbedded sediments. The degree to which the layering impedes vertical flow has not been well understood, yet is a feature that may exert a substantial control on the movement of contaminants. An axial-flow numerical model, RADFLOW, was calibrated to pumping test data collected by a straddle-packer system deployed at 23 depth intervals in four observation wells to evaluate conceptual models and estimate properties of the Snake River Plain aquifer at the Idaho National Engineering and Environmental Laboratory. A delayed water-table response observed in intervals beneath a sediment interbed was best reproduced with a three-layer simulation. The results demonstrate the hydraulic significance of this interbed as a semi-confining layer. Vertical hydraulic conductivity of the sediment interbed was estimated to be about three orders of magnitude less than vertical hydraulic conductivity of the lower basalt and upper basalt units. The numerical model was capable of representing aquifer conceptual models that could not be represented with any single analytical technique. The model proved to be a useful tool for evaluating alternative conceptual models and estimating aquifer properties in this application. Electronic Publication