Neogene aquifer properties specified through the interpretation of electrical sounding data,Salamiyeh region,central Syria

Twenty‐nine Schlumberger electrical soundings were carried out in the Salamiyeh region in Syria using a maximum current electrode separation of 1 km. Three soundings were made at existing boreholes for comparison. Aquifer parameters of hydraulic conductivity and transmissivity were obtained by analysing pumping test data from the existing boreholes. An empirical relationship between hydraulic conductivity determined from the pumping test and both resistivity and thickness of the Neogene aquifer has been established for these boreholes in order to calculate the geophysical hydraulic conductivity. A close agreement has been obtained between the computed hydraulic conductivity and that determined from the pumping test. The relationship established has, therefore, been generalized in the study area in order to evaluate hydraulic conductivity and transmissivity at all the points where geoelectrical measurements have been carried out. This generalization allows one to derive maps of the hydraulic conductivity and transmissivity in the study area based on geoelectrical measurements. These maps are important in future modelling processes oriented towards better exploitation of the aquifers. Copyright © 2006 John Wiley & Sons, Ltd.     

Geophysical Assessment of Groundwater Potential: A Case Study from Mian Channu Area,Pakistan

An integrated study using geophysical method in combination with pumping tests and geochemical method was carried out to delineate groundwater potential zones in Mian Channu area of Pakistan. Vertical electrical soundings (VES) using Schlumberger configuration with maximum current electrode spacing (AB/2 = 200 m) were conducted at 50 stations and 10 pumping tests at borehole sites were performed in close proximity to 10 of the VES stations. The aim of this study is to establish a correlation between the hydraulic parameters obtained from geophysical method and pumping tests so that the aquifer potential can be estimated from the geoelectrical surface measurements where no pumping tests exist. The aquifer parameters, namely, transmissivity and hydraulic conductivity were estimated from Dar Zarrouyk parameters by interpreting the layer parameters such as true resistivities and thicknesses. Geoelectrical succession of five‐layer strata (i.e., topsoil, clay, clay sand, sand, and sand gravel) with sand as a dominant lithology was found in the study area. Physicochemical parameters interpreted by World Health Organization and Food and Agriculture Organization were well correlated with the aquifer parameters obtained by geoelectrical method and pumping tests. The aquifer potential zones identified by modeled resistivity, Dar Zarrouk parameters, pumped aquifer parameters, and physicochemical parameters reveal that sand and gravel sand with high values of transmissivity and hydraulic conductivity are highly promising water bearing layers in northwest of the study area. Strong correlation between estimated and pumped aquifer parameters suggest that, in case of sparse well data, geophysical technique is useful to estimate the hydraulic potential of the aquifer with varying lithology.     

Using direct current resistivity sounding and geostatistics to aid in hydrogeological studies in the Choshuichi alluvial fan, Taiwan

Ground water reservoirs in the Choshuichi alluvial fan, central western Taiwan, were investigated using direct-current (DC) resistivity soundings at 190 locations, combined with hydrogeological measurements from 37 wells. In addition, attempts were made to calculate aquifer transmissivity from both surface DC resistivity measurements and geostatistically derived predictions of aquifer properties. DC resistivity sounding data are highly correlated to the hydraulic parameters in the Choshuichi alluvial fan. By estimating the spatial distribution of hydraulic conductivity from the kriged well data and the cokriged thickness of the correlative aquifer from both resistivity sounding data and well information, the transmissivity of the aquifer at each location can be obtained from the product of kriged hydraulic conductivity and computed thickness of the geoelectric layer. Thus, the spatial variation of the transmissivities in the study area is obtained. Our work is more comparable to Ahmed et al. (1988) than to the work of Niwas and Singhal (1981). The first "constraint" from Niwas and Singhal's work is a result of their use of linear regression. The geostatistical approach taken here (and by Ahmed et al. [1988]) is a natural improvement on the linear regression approach.     

Quaternary aquifer transmissivity derived from vertical electrical sounding measurements in the semi-arid Khanasser Valley Region,Syria

Vertical electrical sounding technique (VES) is used as an alternative approach to pumping test for computing the Quaternary aquifer transmissivity in the Khanasser Valley, Northern Syria. The method is inexpensive, easy and gives faster results with higher special resolution than the traditional pumping technique. The hydraulic conductivity values obtained using VES agree with the pumping test results, which in the Khanasser Valley vary between the order of 0.864 and 8.64 m/day (10⁻⁵ and 10⁻⁴ m/s). The probable location of the Quaternary aquifer in the Khanasser Valley is obtained through the transmissivity map derived from VES. The knowledge of transmissivity is fundamental for modeling and management processes in the Khanasser Valley. Other similar semiarid regions can benefit from the approach successfully applied in the study area.     

Effects of Improper Characterization of Aquifer Thickness on Estimates of Hydraulic Conductivity from Pumping Tests

Pumping test data for surficial aquifers are commonly analyzed under the assumption that the base of the aquifer corresponds to the bottom of the test wells (i.e., the aquifer is truncated). This practice can lead to inaccurate hydraulic conductivity estimates, resulting from the use of low saturated thickness values with transmissivity estimates, and not accounting for the effects of partially penetrating wells. Theoretical time-drawdown data were generated at an observation well in a hypothetical unconfined aquifer for various values of saturated thickness and were analyzed by standard curve-matching techniques. The base of the aquifer was assumed to be the bottom of the pumping and observation wells. The overestimation of horizontal hydraulic conductivity was found to be directly proportional to the error in assumed saturated thickness, and to the (actual) ratio of vertical to horizontal hydraulic conductivity (K_v/K_h). Inaccurately high estimates of hydraulic conductivity obtained by aquifer truncation can lead to overestimates of ground water velocity and contaminant plume spreading, narrow capture zone configuration estimates, and overestimates of available ground water resources.     

An inverse procedure to estimate transmissivity from heads and SP signals

Experimental hydraulic heads and electrical (self-potential) signals associated with a pumping test were used in an inverse model to estimate the transmissivity distribution of a real aquifer. Several works reported in the literature show that there is a relatively good linear relationship between the hydraulic heads in the aquifer and electrical signals measured at the ground surface. In this experimental test field, first, the current coupling coefficient was determined by the best fit between experimental and modeled self-potential signals at the end of the pumping phase. Soon afterward, with the hydraulic heads obtained from the self-potential signals, the transmissivity distribution of the aquifer was conditioned by means an inverse model based on the successive linear estimator (SLE). To further substantiate the estimated T field from the SLE analysis, we analyzed the drawdown rate, the derivative of the drawdown with respect to the ln(t), because the drawdown rate is highly sensitive to the variability in the transmissivity field. In our opinion, these results show that self-potential signals allow the monitoring of subsurface flow in the course of pumping experiments, and that electrical potentials serve as a good complement to piezometric observations to condition and characterize the transmissivity distribution of an aquifer.     

Estimation of anisotropic hydraulic conductivity using geophysical data in a coastal aquifer of Karnataka,India

Estimation of hydraulic parameters is essential to understand the interaction between groundwater flow and seawater intrusion. Though several studies have addressed hydraulic parameter estimation, based on pumping tests as well as geophysical methods, not many studies have addressed the problem with clayey formations being present. In this study, a methodology is proposed to estimate anisotropic hydraulic conductivity and porosity values for the coastal aquifer with unconsolidated formations. For this purpose, the one-dimensional resistivity of the aquifer and the groundwater conductivity data are used to estimate porosity at discrete points. The hydraulic conductivity values are estimated by its mutual dependence with porosity and petrophysical parameters. From these estimated values, the bilinear relationship between hydraulic conductivity and aquifer resistivity is established based on the clay content of the sampled formation. The methodology is applied on a coastal aquifer along with the coastal Karnataka, India, which has significant clayey formations embedded in unconsolidated rock. The estimation of hydraulic conductivity values from the established correlations has a correlation coefficient of 0.83 with pumping test data, indicating good reliability of the methodology. The established correlations also enable the estimation of horizontal hydraulic conductivity on two-dimensional resistivity sections, which was not addressed by earlier studies. The inventive approach of using the established bilinear correlations at one-dimensional to two-dimensional resistivity sections is verified by the comparison method. The horizontal hydraulic conductivity agrees with previous findings from inverse modelling. Additionally, this study provides critical insights into the estimation of vertical hydraulic conductivity and an equation is formulated which relates vertical hydraulic conductivity with horizontal. Based on the approach presented, the anisotropic hydraulic conductivity of any type aquifer with embedded clayey formations can be estimated. The anisotropic hydraulic conductivity has the potential to be used as an important input to the groundwater models.     

Uncertainty estimation of pathlines in ground water models

A method is proposed to estimate the uncertainty of the location of pathlines in two-dimensional, steady-state confined or unconfined flow in aquifers due to the uncertainty of the spatially variable unconditional hydraulic conductivity or transmissivity field. The method is based on concepts of the semianalytical first-order theory given in Stauffer et al. (2002, 2004), which allows estimates of the lateral second moment (variance) of the location of a moving particle. However, this method is reformulated in order to account for nonuniform recharge and nonuniform aquifer thickness. One prominent application is the uncertainty estimation of the catchment of a pumping well by considering the boundary pathlines starting at a stagnation point. In this method, the advective transport of particles is considered, based on the velocity field. In the case of a well catchment, backtracking is applied by using the reversed velocity field. Spatial variability of hydraulic conductivity or transmissivity is considered by taking into account an isotropic exponential covariance function of log-transformed values with parameters describing the variance and correlation length. The method allows postprocessing of results from ground water models with respect to uncertainty estimation. The code PPPath, which was developed for this purpose, provides a postprocessing of pathline computations under PMWIN, which is based on MODFLOW. In order to test the methodology, it was applied to results from Monte Carlo simulations for catchments of pumping wells. The results correspond well. Practical applications illustrate the use of the method in aquifers.     

Use of Groundwater Level Fluctuations near an Operating Water Supply Well to Estimate Aquifer Transmissivity

We developed a method to estimate aquifer transmissivity from the hydraulic-head data associated with the normal cyclic operation of a water supply well thus avoiding the need for interrupting the water supply associated with a traditional aquifer test. The method is based on an analytical solution that relates the aquifer's transmissivity to the standard deviation of the hydraulic-head fluctuations in one or more observation wells that are due to the periodic pumping of the production well. We analyzed the resulting analytical solution and demonstrated that when the observation wells are located near the pumping well, the solution has a simple, Dupuit like form. Numerical analysis demonstrates that the analytical solution can also be used for a quasi-periodic pumping of the supply well. Simulation of cyclic pumping in a statistically heterogeneous medium confirms that the method is suitable for analyzing the transmissivity of weakly or moderately heterogeneous aquifers. If only one observation well is available, and the shift in the phase of hydraulic-head oscillations between the pumping well and the observation well is not identifiable. Prior knowledge of aquifer's hydraulic diffusivity is required to obtain the value of the aquifer transmissivity.     

Conventional and combined pump-and-treat systems under nonuniform background flow

We investigate the performance of vertical hydraulic barriers in combination with extraction wells for the partial hydraulic isolation of contaminated aquifer areas. The potential advantage of such combinations compared to a conventional pump-and-treat system has already been demonstrated in a previous study. Here we extend the scope of the performance analysis to the impact of uncertainty in the regional flow direction as well as to highly heterogeneous aquifer transmissivity distributions. In addition, two new well-barrier scenarios are proposed and analyzed. The hydraulic efficiency of the scenarios is rated based on the expected (mean) reduction of the pumping rate that is required to achieve downgradient contaminant capture. The uncertain spatial distribution of aquifer transmissivity is considered by means of unconditioned Monte Carlo simulations. The significance of uncertain background flow conditions is incorporated by computing minimized pumping rates for deviations of the regional flow direction up to 30 degrees from a normative base case. The results give an answer on how pumping rates have to be changed for each barrier-well combination in order to achieve robust systems. It is exposed that in comparison to installing exclusively wells, the barrier-supported approach generally yields savings in the (average) pumping rate. The particular efficiency is shown to be highly dependent on the interaction of variance and integral scale of transmissivity distribution, well and barrier position, as well as direction of background flow.     

Estimating aquifer properties from the water level response to Earth tides

Water level fluctuations induced by tidal strains can be analyzed to estimate the elastic properties, porosity, and transmissivity of the surrounding aquifer material. We review underutilized methods for estimating aquifer properties from the confined response to earth tides. The earth tide analyses are applied to an open well penetrating a confined carbonate aquifer. The resulting range of elastic and hydraulic aquifer properties are in general agreement with that determined by other investigators for the area of the well. The analyses indicate that passive monitoring data from wells completed in sufficiently stiff, low porosity formations can provide useful information on the properties of the surrounding formation.     

Geometry and dynamics of the freshwater—seawater interface in a coastal aquifer in southeastern Spain

Abstract

The contact between freshwater and seawater in coastal aquifers is studied using a relatively simple model for homogeneous aquifers. However, for real aquifers it is not so simple. The desalination plant built to supply water to the city of Almería is situated over the aquifer in the southern part of the River Andarax Delta. Its design capacity is 1100 L s^?1, and it is supplied from boreholes pumping water from beneath the freshwater—seawater contact in this aquifer. Well logs kept over a period of two years have allowed us to accurately define the interface geometry of the freshwater—seawater contact. Lithological data collected from 31 boreholes have also indicated the existence of strata with low hydraulic conductivity, within others of high conductivity. During a simultaneous pumping test of six wells with 690 L s^?1 total discharge, electrical conductivity measurements showed the influx of seawater 6–10 m below sea level and a drawdown of the interface in the piezometers closest to the pumping wells.     

Random errors in the hydrogeological parameters of a two-layer stratum evaluated by multi-well pumping tests

Random errors in estimates of flow and storage parameters of a two-layer stratum, derived from multi-well pumping test data are examined. It was found that the relative random error with a confidence probability of 0.954 is appreciably less than 100% only for the estimates of aquifer transmissivity, while the estimates of the hydraulic conductivity and the specific yield of the overlying deposits show relative random errors with the same confidence probability in excess of 100%. Therefore, the random errors of parameters are to be given as multiples of the established values.     

GIS integration of hydrogeological and geoelectrical data for groundwater potential modeling in the western part of greater Kushtia district of Bangladesh

Geographical Information System (GIS) has been used in this paper to delineate groundwater resources potential in the western part of greater Kushtia district of Bangladesh, where urgent attention for augmentation of irrigation water supply is required. Thematic maps of transmissivity, hydraulic conductivity, specific yield, net recharge, aquifer thickness, surface water bodies, aquifer resistivity, overburden aquitard thickness and its resistivity have been prepared and assigned weight according to their relative importance using Analytical Hierarchical Process for the preparation of groundwater potential model. Since the values within each thematic map vary significantly, they are classified into various ranges or types and assigned ratings. Finally, the thematic maps are integrated using GIS to prepare the groundwater potential map for the study area in terms of Ground Water Potential Index (GWPI). The evolved map indicates that 22.51% of the study area have GWPI more than 0.70 and therefore, have excellent prospective for exploitation. About 69.12% of the area with GWPI ranging from 0.50 to 0.70 is also quite promising for groundwater abstraction, while the rest 8.37% area having GWPI below 0.50 indicates moderate potential. The obtained map of groundwater potential is found in good agreement with the yields of available pumping test data.     

A Potential-Based Inversion of Unconfined Steady-State Hydraulic Tomography

The importance of estimating spatially variable aquifer parameters such as transmissivity is widely recognized for studies in resource evaluation and contaminant transport. A useful approach for mapping such parameters is inverse modeling of data from series of pumping tests, that is, via hydraulic tomography. This inversion of field hydraulic tomographic data requires development of numerical forward models that can accurately represent test conditions while maintaining computational efficiency. One issue this presents is specification of boundary and initial conditions, whose location, type, and value may be poorly constrained. To circumvent this issue when modeling unconfined steady-state pumping tests, we present a strategy that analyzes field data using a potential difference method and that uses dipole pumping tests as the aquifer stimulation. By using our potential difference approach, which is similar to modeling drawdown in confined settings, we remove the need for specifying poorly known boundary condition values and natural source/sink terms within the problem domain. Dipole pumping tests are complementary to this strategy in that they can be more realistically modeled than single-well tests due to their conservative nature, quick achievement of steady state, and the insensitivity of near-field response to far-field boundary conditions. After developing the mathematical theory, our approach is first validated through a synthetic example. We then apply our method to the inversion of data from a field campaign at the Boise Hydrogeophysical Research Site. Results from inversion of nine pumping tests show expected geologic features, and uncertainty bounds indicate that hydraulic conductivity is well constrained within the central site area.     

The value of subsidence data in ground water model calibration

The accurate estimation of aquifer parameters such as transmissivity and specific storage is often an important objective during a ground water modeling investigation or aquifer resource evaluation. Parameter estimation is often accomplished with changes in hydraulic head data as the key and most abundant type of observation. The availability and accessibility of global positioning system and interferometric synthetic aperture radar data in heavily pumped alluvial basins can provide important subsidence observations that can greatly aid parameter estimation. The aim of this investigation is to evaluate the value of spatial and temporal subsidence data for automatically estimating parameters with and without observation error using UCODE-2005 and MODFLOW-2000. A synthetic conceptual model (24 separate cases) containing seven transmissivity zones and three zones each for elastic and inelastic skeletal specific storage was used to simulate subsidence and drawdown in an aquifer with variably thick interbeds with delayed drainage. Five pumping wells of variable rates were used to stress the system for up to 15 years. Calibration results indicate that (1) the inverse of the square of the observation values is a reasonable way to weight the observations, (2) spatially abundant subsidence data typically produce superior parameter estimates under constant pumping even with observation error, (3) only a small number of subsidence observations are required to achieve accurate parameter estimates, and (4) for seasonal pumping, accurate parameter estimates for elastic skeletal specific storage values are largely dependent on the quantity of temporal observational data and less on the quantity of available spatial data.     

A brief review of kriging and its application to optimal interpolation and observation well selection

ABSTRACT

The current state of kriging in subsurface hydrology is critically reviewed. In an application to a region where boreholes already exist, methods of optimal location of additional observation wells for geophysical parameter investigation and optimal interpolation for the purpose of solving the inverse problem are investigated. The particular case of the location of wells for the measurements of transmissivity and hydraulic head in the Kennet Valley Chalk aquifer, UK, is examined. Results of interpolation of measured hydraulic conductivity values by kriging are compared with results from a standard graphical package for interpolation. Reference is also made to the distribution obtained by the inverse method (in which the conductivity distribution is obtained from the head distribution). On the basis of the application, the conditional simulation (in which the generated data are both consistent with field values and the field statistical structure) is deemed to be the best. It is also found that different methods of interpolation give widely different distributions in the case of hydraulic conductivity. It is suggested that the kriged map or conditional map of the transmissivity should serve as the basis for regional discretization to which corrections via the inverse model may be made.     

Use of Improved Hydrologic Testing and Borehole Geophysical Logging Methods for Aquifer Characterization

Depth-discrete aquifer in formal ion was obtained using recently developed adaptations and improvements to conventional characterization techniques. These improvements included running neutron porosity and hulk density geophysical logging tools through a cased hole, performing an enhanced point-dilution tracer test for monitoring tracer concentration as a function of Lime and depth, and using pressure derivatives for diagnostic and quantitative analysis of constant rate discharge lest data. Data results from the use of these techniques were used to develop a conceptual model of a heterogeneous aquifer. Depth-discrete aquifer information was required to effectively design field-scale deployment and monitoring of an in situ bioremediation technology.
Geophysical logging and point-dilution tracer test results provided the relative distribution of porosity and horizontal hydraulic conductivity, respectively, with depth and correlated well. Hydraulic pumping tests were conducted to estimate mean values for transmissivity and effective hydraulic conductivity, Tracer lest and geophysical logging results indicated that ground water flow was predominant in the upper approximate 10 feet of the aquifer investigated. These results were used to delineate a more representative interval thickness for estimating effective hydraulic conductivity. Hydraulic conductivity, calculated using this representative interval, was estimated lo be 73 ft/d, approximately three limes higher than that calculated using the full length of the screened test interval.     

Identification of Groundwater Parameters Using an Adaptative Multiscale Method

The identification of groundwater parameters in heterogeneous systems is a major challenge in groundwater modeling. Flexible parameterization methods are needed to assess the complexity of the spatial distributions of these parameters in real aquifers. In this article, we introduce an adaptative parameterization to identify the distribution of hydraulic conductivity within the large‐scale (4400 km²) Upper Rhine aquifer. The method is based on adaptative multiscale triangulation (AMT) coupled with an inverse problem procedure that identifies the parameters' distributions by reducing the error between measured and simulated heads. The AMT method has the advantage of combining both zonation and interpolation approaches. The AMT method uses area‐based interpolation rather than an interpolation based on stochastic features. The method is applied to a standard 2D groundwater model that takes into account the interactions between the aquifer and surface water bodies, groundwater recharge, and pumping wells. The simulation period covers 204 months, from January 1986 to December 2002. Recordings at 109 piezometers are used for model calibration. The simulated heads are globally quite accurate and reproduce the main dynamics of the system. The local hydraulic conductivities resulting from the AMT method agree qualitatively with existing local experimental observations across the Rhine aquifer.