Mathematical Modeling of Column and Field Dense Nonaqueous Phase Liquid Tracer Tests

Mathematical models for the simulation of dense nonaqueous phase liquid tracer tests (DTTs) in laboratory columns and in the field are developed and examined. The DTT technique is a means of estimating the quantity of dense nonaqueous phase liquid (DNAPL) in a domain of interest in an aquifer. The two-dimensional field DTT model uses the Method of Principal Directions and an asymmetrical upwind algorithm for describing advective transport. Both models include diffusion transport of tracer into and from low-permeability porous structures such as clay lenses, as well as the mass transport kinetics of partitioning tracer to and from the DNAPL droplets. The dependence of the effluent tracer concentration curves on the parameters of the models is explored, and conclusions are drawn regarding the applicability of, and several possible problems with, the DTT technique. Model results indicate that the DTT performs well at locating distributed droplets of DNAPL, but is unlikely to be useful in the assessment of pooled DNAPL.     

RISK ASSESSMENT METHODOLOGY FOR KARST AQUIFERS: (1) ESTIMATING KARST CONDUIT-FLOW PARAMETERS

Quantitative ground-water tracing of conduit-dominated karst aquifers allows for reliable and practical interpretation of karst ground-water flow. Insights into the hydraulic geometry of the karst aquifer may be acquired that otherwise could not be obtained by such conventional methods as potentiometric-surface mapping and aquifer testing. Contamination of karst aquifers requires that a comprehensive tracer budget be performed so that karst conduit hydraulic-flow and geometric parameters be obtained. Acquisition of these parameters is necessary for estimating contaminant fate-and-transport. A FORTRAN computer program for estimating total tracer recovery from tracer-breakthrough curves is proposed as a standard method. Estimated hydraulic-flow parameters include mean residence time, mean flow velocity, longitudinal dispersivity, Peclet number, Reynolds number, and Froude number. Estimated geometric parameters include karst conduit sinuous distance, conduit volume, cross-sectional area, diameter, and hydraulic depth. These parameters may be used to (1) develop structural models of the aquifer, (2) improve aquifer resource management, (3) improve ground-water monitoring systems design, (4) improve aquifer remediation, and (5) assess contaminant fate-and-transport. A companion paper demonstrates the use of these hydraulic-flow and geometric parameters in a surface-water model for estimating contaminant fate-and-transport in a karst conduit. Two ground-water tracing studies demonstrate the utility of this program for reliable estimation of necessary karst conduit hydraulic-flow and geometric parameters.     

Cleanup of fractured rock aquifers: Implications of matrix diffusion

As contamination moves through a fractured rock aquifer, it tends to diffuse from the flowing fracture water into the rock's essentially stagnant pore water. This process tends both to retard a contamination plume's advance through a fractured rock aquifer and to substantially increase the difficulty of purging contamination from the aquifer. A mathematical model has been developed to evaluate the potential impact of this phenomenon upon water quality restoration in fractured rock aquifers. The numerical modeling reveals that cleanup of fractured rock aquifers will, in many cases, require many decades, even centuries, to achieve, particularly where substantial improvements in water quality are sought. The parameters which most strongly govern the degree to which matrix diffusion prolongs the aquifer restoration process are the rock's matrix porosity, fracture spacing, and matrix diffusivity, the chemical identity of the contaminant(s), and the length of time the aquifer has been contaminated.Since sedimentary rocks tend to have both relatively high matrix porosities and matrix diffusivities, it can be particularly difficult to purge contamination from sedimentary rock aquifers. Crystalline rocks, in contrast, typically have lower matrix porosities and matrix diffusivities, and therefore undergo more rapid cleanup. However, even in crystalline rocks, attainment of very high degrees of water quality improvement may be problematic. Numerical modeling also indicates that conventional groundwater pump and treat programs are not likely to be very effective in speeding up aquifer restoration if the rate limiting step in the process is diffusion of contaminants from the rock matrix.     

Migration of pollutants in groundwater. VI. Flushing of DNAPL droplets/ganglia

Models for describing the flushing of DNAPL from contaminated aquifers are developed, and the dependence of the calculated cleanup times on the model parameters is explored. Diffusion transport from isolated DNAPL droplets, from low-permeability porous spherical domains containing distributed DNAPL droplets, and from low-permeability porous planar lamellae containing distributed DNAPL is analyzed, and the resulting expressions then coupled with the equations for advective transport of dissolved VOC by means of natural uniform flow and a system of injection and recovery wells generating a two-dimensional flow field. The models are readily run on currently available microcomputers. The results of computations with the models are consistent with the severe tailing and slow rates of remediation which are generally observed when DNAPLs are removed by flushing.     

Spatial and temporal properties of environmental exposure assessments related to ground water contamination

This paper describes a new methodology for modeling contaminant transport in ground water for (1) better quantifying the magnitude of exposure in a contaminated aquifer, (2) characterizing the spatial and temporal variation of exposure in a heterogeneous aquifer, and (3) providing more tools for characterizing population exposure. This methodology is also applied to a semi-hypothetical case study.     

Hydrostratigraphy and hydrogeology of the western part of Maira area, Khyber Pakhtunkhwa, Pakistan: a case study by using electrical resistivity

Hydrostratigraphy and hydrogeology of the Maira vicinity is important for the characterization of aquifer system and developing numerical groundwater flow models to predict the future availability of the water resource. Conventionally, the aquifer parameters are obtained by the analysis of pumping tests data which provide limited spatial information and turn out to be costly and time consuming. Vertical electrical soundings and pump testing of boreholes were conducted to delineate the aquifer system at the western part of the Maira area, Khyber Pakhtun Khwa, Pakistan. Aquifer lithology in the eastern part of the study area is dominated by coarse sand and gravel whereas the western part is characterized by fine sand. An attempt has been made to estimate the hydraulic conductivity of the aquifer system by establishing a relationship between the pumping test results and vertical electrical soundings by using regression technique. The relationship is applied to the area along the resistivity profiles where boreholes are not drilled. Our findings show a good match between pumped hydraulic conductivity and estimated hydraulic conductivity. In case of sparse borehole data, regression technique is useful in estimating hydraulic properties for aquifers with varying lithology.     

Comparison of Groundwater Level Estimation Using Neuro-fuzzy and Ordinary Kriging

Water level in aquifer plays the main role in groundwater modeling as one of the input data. In practice, due to aspects of time and cost, data monitoring of water levels is conducted at a limited number of sites, and interpolation technique such as kriging is widely used for estimation of this variable in unsampled sites. In this study, the efficiency of the ordinary kriging (OK) and adaptive network-based fuzzy inference system (ANFIS) was investigated in interpolation of groundwater level in an unconfined aquifer in the north of Iran. The results showed that ANFIS model is more efficient in estimating the groundwater level than OK.     

Groundwater modeling of Saq Aquifer Buraydah Al Qassim for better water management strategies

Saudi Arabia is an arid country. It has limited water supplies. About 80?C90% of water supplies come from groundwater, which is depleting day by day. It needs appropriate management. This paper has investigated groundwater modeling of Saq Aquifer in Buraydah Al Qassim to estimate the impact of its excessive use on depletion of Saq Aquifer. MODFLOW model has been used in this study. Data regarding the aquifer parameters was measured by pumping tests. Groundwater levels and discharge of wells in the area for the year 2008 and previous record of year 1999 have been collected from Municipal Authority of Buraydah. Location of wells was determined by Garmin. The model has been run for different sets of pumping rates to recommend an optimal use of groundwater resources and get prolonged life of aquifer. Simulations have been made for a long future period of 27?years (2008?C2035). Model results concluded that pumping from the Saq Aquifer in Buraydah area will result into significant cones of depression if the existing excessive pumping rates prevail. A drawdown up to 28?m was encountered for model run for 27?years for existing rates of pumping. Aquifer withdrawals and drawdowns will be optimal with the conservation alternative. The management scheme has been recommended to be adopted for the future protection of groundwater resources in Kingdom of Saudi Arabia.     

Estimating contaminant attenuation half-lives in alluvial groundwater systems

One aspect of describing contamination in an alluvial aquifer is estimating changes in concentrations over time. A variety of statistical methods are available for assessing trends in contaminant concentrations. We present a method that extends trend analysis to include estimating the coefficients for the exponential decay equation and calculating contaminant attenuation half-lives. The conceptual model for this approach assumes that the rate of decline is proportional to the contaminant concentration in an aquifer. Consequently, the amount of time to remove a unit quantity of the contaminant inventory from an aquifer lengthens as the concentration decreases. Support for this conceptual model is demonstrated empirically with log-transformed time series of contaminant data. Equations are provided for calculating system attenuation half-lives for non-radioactive contaminants. For radioactive contaminants, the system attenuation half-life is partitioned into the intrinsic radioactive decay and the concentration reduction caused by aquifer processes. Examples are presented that provide the details of this approach. In addition to gaining an understanding of aquifer characteristics and changes in constituent concentrations, this method can be used to assess compliance with regulatory standards and to estimate the time to compliance when natural attenuation is being considered as a remediation strategy. A special application of this method is also provided that estimates the half-life of the residence time for groundwater in the aquifer by estimating the half life for a conservative contaminant that is no longer being released into the aquifer. Finally, the ratio of the half-life for groundwater residence time to the attenuation half-life for a contaminant is discussed as a system-scale retardation factor which can be used in analytical and numerical modeling.     

Implications and concerns of deep-seated disposal of hydrocarbon exploration produced water using three-dimensional contaminant transport model in Bhit Area, Dadu District of Southern Pakistan

A three-dimensional contaminant transport model has been developed to simulate and monitor the migration of disposal of hydrocarbon exploration produced water in Injection well at 2,100 m depth in the Upper Cretaceous Pab sandstone, Bhit area in Dadu district of Southern Pakistan. The regional stratigraphic and structural geological framework of the area, landform characteristics, meteorological parameters, and hydrogeological milieu have been used in the model to generate the initial simulation of steady-state flow condition in the underlying aquifer’s layers. The geometry of the shallow and deep-seated characteristics of the geological formations was obtained from the drilling data, electrical resistivity sounding surveys, and geophysical well-logging information. The modeling process comprised of steady-state simulation and transient simulation of the prolific groundwater system of contamination transport after 1, 10, 30 years of injection. The contaminant transport was evaluated from the bottom of the injection well, and its short- and long-term effects were determined on aquifer system lying in varying hydrogeological and geological conditions.     

Arsenic-safe alternate aquifers and their hydraulic characteristics in contaminated areas of Middle Ganga Plain, Eastern India

Arsenic groundwater contamination exceeding 0.05 mg/l affecting the Newer Alluvial tracts of Patna and Bhojpur, the two worst affected districts located in the Middle Ganga Plain in the Bihar state, has been studied The area is underlain by two-tier Quaternary aquifer systems within a depth of 300 m below ground level, separated by a 15?C32-m-thick clay/sandy clay aquitard. The upper part (<50 m depth) of the shallow aquifer system is arsenic-contaminated. The deeper aquifer system (lying below 120?C130 m depth) exhibits low arsenic load (max 0.0035 mg/l), having hydraulic conductivity between 64.88 and 82.04 m/day. Groundwater in the deeper aquifer occurs under semi-confined to confined condition due to poor hydraulic conductivity of the middle clay (4.7 × 10^???2???7.2 × 10^???3 m/day). Hydraulic head of the deeper aquifer remains close to the surface than the shallow aquifer. The two aquifer systems in the Newer Alluvium are replaced by a thick single aquifer system in the adjoining Older Alluvium, within a depth of 330 m below ground. In the arsenic-contaminated area, deeper aquifer is protected by a middle clay, which may be developed for community drinking water supply by deep tube wells having a yield capacity of 150 m³/h.     

Modelling Tracer Dispersion from Landfills

Several wind tunnel experiments of tracer dispersion from reduced-scale landfill models are presented in this paper. Different experimental set-ups, hot-wire anemometry, particle image velocimetry and tracer concentration measurements were used for the characterisation of flow and dispersion phenomena nearby the models. The main aim of these experiments is to build an extensive experimental data set useful for model validation purposes. To demonstrate the potentiality of the experimental data set, a validation exercise on several mathematical models was performed by means of a statistical technique. The experiments highlighted an increase in pollutant ground level concentrations immediately downwind from the landfill because of induced turbulence and mean flow deflection. This phenomenon turns out to be predominant for the dispersion process. Tests with a different set-up showed an important dependence of the dispersion phenomena from the landfill height and highlighted how complex orographic conditions downwind of the landfill do not affect significantly the dispersion behaviour. Validation exercises were useful for model calibration, improving code reliability, as well as evaluating performances. The Van Ulden model proved to give the most encouraging results.     

Modeling of nitrate concentration in groundwater using artificial intelligence approach—a case study of Gaza coastal aquifer

Nitrate concentration in groundwater is influenced by complex and interrelated variables, leading to great difficulty during the modeling process. The objectives of this study are (1) to evaluate the performance of two artificial intelligence (AI) techniques, namely artificial neural networks and support vector machine, in modeling groundwater nitrate concentration using scant input data, as well as (2) to assess the effect of data clustering as a pre-modeling technique on the developed models' performance. The AI models were developed using data from 22 municipal wells of the Gaza coastal aquifer in Palestine from 2000 to 2010. Results indicated high simulation performance, with the correlation coefficient and the mean average percentage error of the best model reaching 0.996 and 7 %, respectively. The variables that strongly influenced groundwater nitrate concentration were previous nitrate concentration, groundwater recharge, and on-ground nitrogen load of each land use land cover category in the well's vicinity. The results also demonstrated the merit of performing clustering of input data prior to the application of AI models. With their high performance and simplicity, the developed AI models can be effectively utilized to assess the effects of future management scenarios on groundwater nitrate concentration, leading to more reasonable groundwater resources management and decision-making     

Remediating Over-Produced and Contaminated Aquifers by Artificial Recharge from Surface Waters

Recharging from open-surface water resources is investigated as a method for remediation of over-produced and contaminated aquifers. The hydraulic and contaminant modeling using the Modflow-2000 and MT3DMS simulation software are resorted for evaluation of the available options. Methodology is developed and illustrated for reviewing the alternative remediation approaches for optimal selection of the best remediation approach to aquifer recharging from surface waters. The effectiveness of the present systematic approach is demonstrated by identifying the best choice among the multifunctional artificial recharge options available for the Great Forest Park in Debrecen, Hungary.     

A Compartmental Model for Stable Time-Dependent Modeling of Surface Water and Groundwater

Water resource system planners make decisions that guide water management policy. The fundamental tools for assessing management and infrastructure strategies are combined hydro-economic models of river basins (RBHE models). These models have improved the economic efficiency of water use in situations of competition for scarce water resources. In RBHE models, a groundwater model is coupled with surface water models of the various water resources. Today, the groundwater models used in an RBHE model can be of two types: cell models or numerical models. Cell models are easy to use, but they are too simple to realistically describe the geology and hydrology of the area under investigation. Numerical models, in contrast, are closer to the physical behavior of the aquifer. However, the vast quantity of data to be analyzed makes them impractical for many management scenarios. Moreover, the calibrations of these high-resolution models are generally difficult and sensitive to the variation of parameters, especially when boundary conditions are dynamic. This is the case when dynamic river data or dynamic surface lake data are present. In this work, a compartmental cell model is built on the hydrogeology of the aquifer. In this approach, the hydrogeology of the aquifer and the dynamic boundary conditions are treated with separate models. A general mathematical formulation is presented where the calibration stage, the validation stage, and the prediction stage are formulated as a series of sub-model calibrations and solved using a general least squares routine. With this approach, it becomes possible to treat both the water level and the pumping rate in each cell as variables to be predicted. In most of the models, the pumping rates are known and the goals of the computation are to estimate the groundwater level. However, when for political or technical reasons access to some of the wells is difficult, the pumping rates are only partially known. Then, both groundwater levels and pumping rates are variables to be predicted by the groundwater model. A computer program was developed using MATLAB, with a Visual Basic graphical user interface using COM technology to access the advanced mathematical libraries. The approach was implemented with a real case study of the Yarkon–Taninim aquifer in Israel. The results indicate that the method is more stable than the classical approach.     

Development of a Model for Simulation of Solute Transport in a Stream–Aquifer System

A model for simulation of solute transport in a dynamic stream–aquifer system is developed by integrating four existing sub-models. Interface packages are created to link the sub-models. The developed model is successfully used to simulate chloride transport in the stream–aquifer system of the Arkansas River and the Equus Beds Aquifer in Kansas and demonstrates that chloride concentration in the aquifer decreases in the vicinity of the simulated channel over time.     

A Generalized Analytical Model for Crosswind-Integrated Concentrations with Ground-Level Deposition in the Atmospheric Boundary Layer

A generalized mathematical model describing the crosswind-integrated concentrations is presented for dispersion of pollutants emitted from a continuous source in the atmospheric boundary layer with deposition to the ground surface. The model is based on a solution of the resulting two-dimensional steady state advection-diffusion equation with deposition to the ground surface. It considers the horizontal wind speed as a generalized function of vertical height above the ground surface and eddy diffusivity as a function of both downwind distance from the source and vertical height. Various special cases of the model are deduced. A sensitivity analysis of the model prediction of the ground-level crosswind-integrated concentrations with deposition velocity is performed. Various issues and limitations associated with this work are discussed. The model is evaluated with the observations of a depositing tracer obtained from Hanford dual diffusion experiment in stable conditions. The statistical measures show that the present model, by considering deposition, is performing well with the observations. The model is giving an over-predicting trend with the observations and predicts 100 % cases within a factor of two. On the other hand, the consideration of a non-depositing condition for a depositing tracer yields the severe over-prediction and thus, it introduces the significant errors in the model prediction. The selection of the lower boundary condition for a depositing tracer at the height of deposition surface gives better prediction than those at a height of surface roughness length.     

Hydrochemical and geoelectrical investigation of the coastal shallow aquifers in El-Omayed area, Egypt

Monitoring and assessment of the coastal aquifers are becoming a worldwide concern for the need of additional and sustainable water resources to satisfy demographic growth and economic development. A hydrochemical and geoelectrical investigation was conducted in the El-Omayed area in the northwestern coast of Egypt. The aim of the study was to delineate different water-bearing formations, provide a general evaluation of groundwater quality, and identify the recharge sources in aquifers. Thirty-seven water samples were collected and chemically analyzed from the sand dune accumulations and oolitic limestone aquifers. Fifteen profiles of vertical electrical soundings (VESs) were obtained in the oolitic limestone aquifer to examine the variations of subsurface geology and associated groundwater chemistry. The groundwater reserves in the El-Omayed area are mainly contained in sand dune accumulations and oolitic limestone aquifers. The aquifer of sand dune accumulations contains freshwater of low salinity (average total dissolved solids (TDS)?=?974 mg/l). Groundwater of oolitic limestone aquifer is slightly brackish (average TDS?=?1,486 mg/l). Groundwater of these aquifers can be used for irrigation under special management for salinity control, and regular leaching as indicated by electrical conductivity and sodium adsorption ratio. Results of VES interpretation classified the subsurface sequence of oolitic limestone aquifer into four geoelectric zones, with increasing depth, calcareous loam, gypsum, oolitic limestone, and sandy limestone. Oolitic limestone constitutes the main aquifer and has a thickness of 12–32 m.     

我国污染场地含水层监测现状与技术研发趋势

通过分析我国污染场地含水层监测现状和污染场地环境管理需求,指出我国污染场地含水层监测的不足之处并探讨了我国污染场地含水层监测技术研发趋势.当前,我国污染场地含水层监测主要以地下水监测为主,且监测仪器设备简陋陈旧、技术手段落后,不能满足数量日益增长的污染场地环境管理需求,急需发展一批新型、高效、经济、实用的监测技术和监测...     

Arsenic mobilization and attenuation by mineral-water interactions: implications for managed aquifer recharge

Managed aquifer recharge (MAR) has potential for addressing deficits in water supplies worldwide. It is also widely used for preventing saltwater intrusion, maintaining the groundwater table, and augmenting ecological stream flows, among many other beneficial environmental applications. However, field MAR sites have experienced arsenic mobilization from aquifer formation minerals due to induced changes in groundwater chemistry. To address this environmental concern, it is crucial to understand the potential sources and sinks impacting arsenic mobilization. This paper outlines important mineral-water interactions that can occur at MAR sites. Detailed information on minerals of concern, physiochemical processes for arsenic mobilization or attenuation, and the potential impact of microbial activity and hydrology on these processes is provided. Based on these mineral-water interactions, guidelines for predicting arsenic mobility are presented, and recommendations are made concerning MAR site monitoring. The review emphasizes important aspects in correlating interfacial reactions to reactive transport modeling and elucidating future challenges, a first step toward developing safer and more sustainable MAR operations.