A hydrogeophysical model of the relationship between geoelectric and hydraulic parameters of anisotropic aquifers

An electrical resistivity method has been used to determine aquifer parameters in the Ganga-Yamuna interfluve in northern India. An existing relationship between the geoelectrical and hydraulic parameters has been modified for the case of an anisotropic aquifer. The hydrogeological framework in the upper part of the Ganga-Yamuna interfluve is evaluated by using existing relationships between hydraulic parameters and geoelectrical parameters for alluvial aquifers. On the basis of aquifer geometry, the area has been divided into two hydraulic units: the western Yamuna flood plain and the Ganga flood plain towards the east. The resistivity data collected in parts of the study area are first interpreted in terms of true resistivity and thicknesses of subsurface layers. The electrical parameters (resistivity and thicknesses) are subsequently correlated with the available pumping test data. Distinct correlations between transmissivity and modified transverse resistance are obtained for the two hydraulic units. A four-parameter model consisting of hydraulic conductivity, modified longitudinal resistivity, modified transverse resistance and hydraulic anisotropy is presented for the anisotropic aquifer underlain by conductive fine grained sediments. The model has been validated at a number of locations, where aquifer parameters are known from pumping test data.     

Evaluating the salinity distribution of a shallow coastal aquifer by vertical multielectrode profiling (Denmark)

A monitoring system, including five groups of piezometers and five vertical multielectrode profiling probes (VMEP), has been installed in an aquifer beneath a coastal dune in Denmark. In order to assess the salinity distribution within the aquifer, geoelectrical data were gathered in March, June and September 2008, by measuring a dipole-dipole and gradient array using multielectrode profiling. Interpretation of the processed resistivity data was performed by regularized inversion using a one-dimensional, horizontally layered model of formation resistivity. The standard deviation on estimated layer log-resistivity was 0.01–0.03. By estimating two parameters of a power function, observed fluid conductivities derived from samples of porewater were related to corresponding estimated formation resistivities. The conductivity profiles correlate with a winter situation in March with high sea level, active recharge and significant wave activity, causing increased hydraulic heads, a thicker freshwater lens and salt water overlying freshwater close to the sea. In June, the thickness of the freshwater lens is reduced due to less recharge and prevailing offshore winds, imposing density-stable conditions and a sharper transition between fresh and brackish water. During the autumn, aquifer recharge is enhanced and hydraulic heads increase, resulting in a thicker freshwater lens.     

Hydrogeological framework and estimation of aquifer hydraulic parameters using geoelectrical data: a case study from West Iran

Detailed local geological, geophysical, and hydrogeological investigations were carried out for the alluvial aquifer in the Kangavar basin, West Iran to delineate the architecture of different subsurface geological horizons using lithologs and generated vertical electrical sounding (VES) data. An attempt has also been made to estimate aquifer transmissivity from resistivity data. Forty VESs were recorded with the Schlumberger electrode configuration in the study area; 28 of these were selected for evaluation. The maximum current electrode spacing was 400–500 m. The data obtained were interpreted by computer iterative modeling with curve matching for calibration purposes. In order to ascertain the subsurface geological framework, the general distribution of resistivity responses of the geological formations was obtained and geoelectrical sections along a number of lines were prepared. Probable aquifer horizons from these sections were identified. The transmissivity of the unconfined aquifer was computed by determining the Dar-Zarrouk parameters (longitudinal unit conductance and transverse unit resistance) and were compared with the actual field transmissivity. The results showed a direct relation between aquifer transmissivity and modified transverse resistance.     

Real surface conductivity component as indicator for the hydraulic conductivity

Estimation of hydraulic conductivity from surface resistivity measurements is one of the most difficult and challenging hydrogeophysical targets. The promising side of this relation is the analogy between electric current flow and water flow, whereas the grand ambiguity is the non-dimensionality between both two quantities. Imaginary surface conductivity component is used recently to deduce the hydraulic conductivity via complex resistivity measurements. Since there are similar properties between imaginary (out-of-phase) and real (in-phase) surface conductivity components, the latter is used in this paper to predict the hydraulic conductivity. Two mathematical parameters were determined to express the electrical equivalent of hydraulic conductivity in sand and clay systems based on the mode of electrical double-layer formation in both systems. The reliability of the proposed method is tested through applying on two datasets representing sand and clay systems. The first dataset is a clean sand and gravel aquifer in the Keritis basin in Chania, Crete, Greece. The second is mostly clayey sand aquifer in Wadi El-Assuity, Egypt. Application of the present approach in these two cases resulted in promising nearly identical values with the measured hydraulic conductivity via pumping test or geometric hydraulic conductivity via grain size analysis.     

Estimation of regional-scale groundwater flow properties in the Bengal Basin of India and Bangladesh

Quantitative evaluation of management strategies for long-term supply of safe groundwater for drinking from the Bengal Basin aquifer (India and Bangladesh) requires estimation of the large-scale hydrogeologic properties that control flow. The Basin consists of a stratified, heterogeneous sequence of sediments with aquitards that may separate aquifers locally, but evidence does not support existence of regional confining units. Considered at a large scale, the Basin may be aptly described as a single aquifer with higher horizontal than vertical hydraulic conductivity. Though data are sparse, estimation of regional-scale aquifer properties is possible from three existing data types: hydraulic heads, ¹⁴C concentrations, and driller logs. Estimation is carried out with inverse groundwater modeling using measured heads, by model calibration using estimated water ages based on ¹⁴C, and by statistical analysis of driller logs. Similar estimates of hydraulic conductivities result from all three data types; a resulting typical value of vertical anisotropy (ratio of horizontal to vertical conductivity) is 10⁴. The vertical anisotropy estimate is supported by simulation of flow through geostatistical fields consistent with driller log data. The high estimated value of vertical anisotropy in hydraulic conductivity indicates that even disconnected aquitards, if numerous, can strongly control the equivalent hydraulic parameters of an aquifer system.     

Estimation of aquifer hydraulic parameters from surface geophysical measurements: a case study of the Upper Cretaceous aquifer, central Sinai, Egypt

The integration of geophysical data with direct hydrogeological measurements can provide a minimally invasive approach to characterize the subsurface at a variety of resolutions and over many spatial scales. The field of hydrogeophysics has attracted much attention during the last two decades. In this domain, the geophysical data inverted to geophysical models are interpreted in terms of the hydrogeology to serve as a basis for the definition of hydraulic models in the areas of interest. The hydraulic conductivity (K) value measured in a reference borehole has been combined with the electrical conductivity obtained from nearby geo-electromagnetic sounding data in the Cenomanian (Upper Cretaceous) aquifer, central Sinai, Egypt. The resulting relation was interpreted with Dar Zarrouk parameters to infer the transmissivity variations at other vertical electrical sounding locations, where K values are unknown. Coincident transient electromagnetic data have been adopted to increase accuracy while interpreting the aquifer geoelectrical properties. The results indicate that the transmissivity values in the aquifer of interest vary from 2,446 to 9,694 m²/day, and K varies from 12.9 to 57.0 m/day throughout the studied area.     

Geophysical characterization of the León-Chinandega aquifer,Nicaragua

Electrical geophysical surveys in the mode of vertical electrical soundings (VES) and continuous vertical electrical soundings (CVES) were conducted in the León-Chinandega plains, northwestern Nicaragua, in order to obtain detailed information about the geometry of the different hydrogeological layers in the aquifer and depth to the basement. A total of 51 VES were carried out within the plains. The results show a complex structure towards the north east of the area, and the southwestern part of the plains presents a smoother stratification. The geoelectrical measurements and borehole information indicate that the basement topography is characterized by hills and deep valleys with highly variable basement depths. Fifty CVES where done in a smaller area in the center of the plain. The resistivity data yielded considerable information revealing the existence of two main geo-electrical units. The combined interpretation of geological and geophysical data shows an environment typical of sedimentary volcanic coastal plains. The information collected during this investigation provides valuable data for estimating the fresh-water resources of the León-Chinandega aquifer system and for development of a groundwater management plan.     

Hydrogeophysical investigation and estimation of groundwater potentials of the Lower Palaeozoic to Precambrian crystalline basement rocks in Keffi area, north-central Nigeria, using resistivity methods

The use of resistivity sounding and two-dimensional (2-D) resistivity imaging was investigated with the aim of delineating and estimating the groundwater potential in Keffi area. Rock types identified are mainly gneisses and granites. Twenty-five resistivity soundings employing the Schlumberger electrode array were conducted across the area. Resistivity sounding data obtained were interpreted using partial curve matching approach and 1-D inversion algorithm, RESIST version 1.0. The 2-D resistivity imaging was also carried out along two traverses using dipole–dipole array, and the data obtained were subjected to finite element method modeling using DIPRO inversion algorithm to produce a two-dimensional subsurface geological model. Interpretation of results showed three to four geoelectrical layers. Layer thickness values were generally less than 2 m for collapsed zone, and ranged from 5 to 30 m for weathered bedrock (saprolite). Two major aquifer units, namely weathered bedrock (saprolite) aquifer and fractured bedrock (saprock) aquifer, have been delineated with the latter usually occurring beneath the former in most areas. Aquifer potentials in the area were estimated using simple schemes that involved the use of three geoelectrical parameters, namely: depth to fresh bedrock, weathered bedrock (saprolite) resistivity and fractured bedrock (saprock) resistivity. The assessment delineated the area into prospective high, medium and low groundwater potential zones.     

Hydrogeophysical imaging of alluvial aquifers: electrostratigraphic units in the quaternary Po alluvial plain (Italy)

The integration of surface geological and geomorphological information with borehole point-data and geophysical (e.g., geoelectrical) images of the subsurface yields spatially consistent representations of alluvial aquifers heterogeneity at different scales, from depositional systems to basin fills. Such an approach requires a conceptual framework to match the stratigraphic units with their evidence from ground-based DC resistivity methods to effectively fill the gaps between sparse borehole data and to obtain valid representations of sedimentary heterogeneities. Such an approach is applied to characterize two sites of the Quaternary aquifers of the central Po Plain (Italy), which represent (1) the middle-upper Pleistocene braided to meandering river depositional systems sitting on Southalpine crust and (2) their down-current counterparts, where they are involved by the latest uplift and deformation due to the tectonic activity of the Apennine frontal thrusts. Electrical resistivity was considered as a proxy of the litho-textural properties of hydrofacies and their major hierarchical association at depth and was interpreted in accordance with the depth-decreasing resolution of ground-based resistivity methods. Thus, it was possible to identify the geophysical signature of hydrostratigraphic units through “Electrostratigraphic Units”, i.e., sedimentary volumes identified by resistivity contrasts that spatially preserve the vertical polarity. Hydrostratigraphy and electrostratigraphy were then joined together through a site-specific relationship between electrical resistivity and hydraulic conductivity, which takes into account the prevailing process of current conduction, the litho-textural properties of hydrofacies and the groundwater electrical conductivity. At the scales of aquifer systems and complexes, this approach permitted to establish the conceptual framework to match hydrostratigraphy, electrostratigraphy, average hydrodynamic properties and distribution of heterogeneities.     

Delineating subsurface heterogeneity at a loop of River Steinlach using geophysical and hydrogeological methods

Stratigraphic heterogeneity is a key controlling factor for flow and transport in groundwater systems. In this case study, we have combined two- and three-dimensional images of electric resistivity with hydrogeological data to delineate such stratigraphic heterogeneity at a loop of River Steinlach close to Tübingen, Germany, where we estimated spatially varying aquifer parameters using established evaluation techniques. We developed a three-layer stratigraphic model consisting of heterogeneous sandy gravel overlain by a fine alluvium including top soil, and underlain by sandstone weathered at its top. The sandy gravel makes up the aquifer unit while we interpret the sandstone as aquitard. We classified the area into two stratigraphic segments based on resistivity data. The Southern segment consists of a thicker alluvium where the aquifer material contains more clay than in the Northern segment. Two different stratigraphic structures were also delineated at the Southern and Northern segments, respectively. These structures are visible in inversion results as low-resistivity features trending in the SE–NW in the Southern segment and NE–NW in the Northern segment. Both structures may be previous flow paths of River Steinlach, and the low permeability material overlying them may impair flow exchange in the area. We performed slug and pumping tests to estimate the hydraulic conductivity (K) of the aquifer. The K distributions show slight variation with higher values obtained at wells within the Northern segment and close to the river bank. A qualitative comparison of measured hydraulic conductivities with the resistivity distribution shows good agreement of the spatial patterns. The stratigraphic and hydraulic heterogeneities delineated in this work are important for experimental and modeling studies of flow, transport, and hyporheic exchange at the site.     

Evaluation of methods for determination of hydraulic properties in an aquifer–aquitard system hydrologically connected to a river

Field and laboratory methods have been used to determine the hydraulic properties in a multiple-layer aquifer–aquitard system that is hydrologically connected to a river. First, hypothetical pumping tests in aquifer–aquitard systems were performed to evaluate the feasibility of MODFLOW-PEST in determining these parameters. Sensitivity analyses showed that: the horizontal hydraulic conductivity in the aquifer has the highest composite sensitivity; the vertical hydraulic conductivity has higher composite sensitivity than the horizontal hydraulic conductivity in the aquitard; and a partial penetration pumping well in an aquifer layer can improve the quality of the estimated parameters. This inverse approach was then used to analyze a pumping-recovery test conducted near the Platte River in southeastern Nebraska, USA. The hydraulic conductivities and specific yield were calculated for the aquitard and aquifer. The direct-push technique was used to generate sediment columns; permeameter tests on these columns produced the vertical hydraulic conductivities that are compatible with those obtained from the pumping-recovery test. Thus, the combination of the direct-push technique with permeameter tests provides a new method for estimation of vertical hydraulic conductivity. The hydraulic conductivity, determined from grain-size analysis, is smaller than the horizontal one but larger than the vertical one determined by the pumping-recovery test.     

Salinity mapping of coastal groundwater aquifers using hydrogeochemical and geophysical methods: a case study from north Kelantan,Malaysia

Integrated hydrogeochemical and geophysical methods were used to study the salinity of groundwater aquifers along the coastal area of north Kelantan. For the hydrogeochemical investigation, analysis of major ion contents of the groundwater was conducted, and other chemical parameters such as pH and total dissolved solids were also determined. For the geophysical study, both geoelectrical resistivity soundings and reflection seismic surveys were conducted to determine the characteristics of the subsurface and groundwater contained within the aquifers. The pH values range from 6.2 to 6.8, indicating that the groundwater in the study area is slightly acidic. Low content of chloride suggests that the groundwater in the first aquifer is fresh, with an average concentration of about 15.8 mg/l and high geoelectrical resistivity (>45 ohm m). On the other hand, the groundwater in the second aquifer is brackish, with chloride concentration ranging from 500 mg/l to 3,600 mg/l and very low geoelectrical resistivity (<45 ohm m) as well as high concentration of total dissolved solids (>1,000 mg/l). The groundwater in the third aquifer is fresh, with chloride concentrations generally ranging from 2 mg/l to 210 mg/l and geoelectrical resistivity of greater than 45 ohm m. Fresh and saltwater interface in the first aquifer is generally located directly in the area of the coast, but, for the second aquifer, both hydrogeochemical and geoelectrical resistivity results indicate that the fresh water and saltwater interface is located as far as 6 km from the beach. The considerable chloride ion content initially suggests that the salinity of the groundwater in the second aquifer is probably caused by the intrusion of seawater. However, continuous monitoring of the chloride content of the second aquifer indicated no significant changes with time, from which it can be inferred that the salinity of the groundwater is not affected by seasonal seawater intrusion. Schoeller diagrams illustrate that sulphate concentrations of the groundwater of the second aquifer are relatively low compared to those of the recent seawater. Therefore, this result suggests that the brackish water in the second aquifer is probably from ancient seawater that was trapped within the sediments for a long period of time, rather than due to direct seawater intrusion.     

Using Sequential Self-Calibration Method to Identify Conductivity Distribution: Conditioning on Tracer Test data

An iterative inverse method, the sequential self-calibration method, is developed for mapping spatial distribution of a hydraulic conductivity field by conditioning on nonreactive tracer breakthrough curves. A streamline-based, semi-analytical simulator is adopted to simulate solute transport in a heterogeneous aquifer. The simulation is used as the forward modeling step. In this study, the hydraulic conductivity is assumed to be a deterministic or random variable. Within the framework of the streamline-based simulator, the efficient semi-analytical method is used to calculate sensitivity coefficients of the solute concentration with respect to the hydraulic conductivity variation. The calculated sensitivities account for spatial correlations between the solute concentration and parameters. The performance of the inverse method is assessed by two synthetic tracer tests conducted in an aquifer with a distinct spatial pattern of heterogeneity. The study results indicate that the developed iterative inverse method is able to identify and reproduce the large-scale heterogeneity pattern of the aquifer given appropriate observation wells in these synthetic cases.     

The use of vertical electrical sounding resistivity method for the location of low salinity groundwater for irrigation in Chaj and Rachna Doabs

A geoelectrical resistivity survey using vertical electrical sounding (VES) was conducted at Chaj Doab (land between rivers Jhelum and Chenab, Pakistan) and Rachna Doab (land between rivers Chenab and Ravi, Pakistan), with the objective of investigating groundwater conditions. A total of 90 sites were selected with 43 sites in Chaj and 47 sites in Rachna Doabs. The resistivity meter (ABEM Terrameter SAS 4000, Sweden) was used to collect the VES data by employing a Schlumberger electrode configuration, with half current electrode spacings (AB/2) ranging from 2 to 180 m and the potential electrode (MN) from 1 to 40 m. The field data were interpreted using the Interpex IX1D computer software and the resistivity versus depth models for each location was estimated. The outputs of subsurface layers with resistivities and thickness presented in contour maps and 3-D views by using SURFER software were created. A total of 102 groundwater samples from nearby hydrowells at different depths were collected to develop a correlation between the aquifer resistivity of VES and the electrical conductivity (EC) of the groundwater and to confirm the resulted geophysical resistivity models. From the correlation developed, it was observed that the groundwater salinity in the aquifer may be considered low and so safe for irrigation if resistivity >45 Ω m, and marginally fit for irrigation having resistivity between 25 and 45 Ω m. The study area has resistivities from 3.9 to 2,222 Ω m at the top of the unsaturated layer, between 1.21 and 171 Ω m, in the shallow aquifers, and 0.14–152 Ω m in the deep aquifers of the study area. The results indicate that the quality of groundwater is better near the rivers and in the shallow layers compared to the deep layers.     

Geo-electrical delineations within Kubanni basin, north-central Nigeria

The research site is the whole landmass of the Federal College of Education, Zaria, seated on basement complex of north-central Nigeria. Direct current resistivity geophysical method was employed to characterise parameters such as the basement depth and topography, aquifer depth and thickness, weathered basement distribution as well as mapping of orientations of fractures and faults present in the premises using radial sounding technique. The conventional vertical electrical sounding (VES) Schlumberger array was carried out at 40 stations, and eight of which were radial stations. Radial sounding was used to establish resistivity anisotropy which gives clue for the choice of consistent VES profile direction used throughout the fieldwork. Results from the resistivity interpretation suggest three layers in most parts of the premises with some minor occurrence of two and four layers. The first layer (topsoil) has its thickness ranging between 3.5 and 14.0 m; second layer (weathered basement) thickness ranges between 9.0 and 36.5 m, while the third layer (fresh basement) is deepest (40.1 m) towards the eastern corner of the area. The aquifer depth ranges from 1.5 to 4.0 m with a thickness range of 5.0 to 14.0 m. The thickest aquifer occurs around the centre to the west in the area. Results from radial sounding show presence of resistivity anisotropy, an insight to fracturing and faulting; this is more pronounced around the west-central part of the premises.     

Estimating regional-scale fractured bedrock hydraulic conductivity using discrete fracture network (DFN) modeling

Estimating bedrock hydraulic conductivity of regional fractured aquifers is challenging due to a lack of aquifer testing data and the presence of small and large-scale heterogeneity. This study provides a novel approach for estimating the bedrock hydraulic conductivity of a regional-scale fractured bedrock aquifer using discrete fracture network (DFN) modeling. The methodology is tested in the mountainous Okanagan Basin, British Columbia, Canada. Discrete fractures were mapped in outcrops, and larger-scale fracture zones (corresponding to lineaments) were mapped from orthophotos and LANDSAT imagery. Outcrop fracture data were used to generate DFN models for estimating hydraulic conductivity for the fractured matrix (K _m). The mountain block hydraulic conductivity (K _mb) was estimated using larger-scale DFN models. Lineament properties were estimated by best fit parameters for a simulated pumping test influenced by a fracture zone. Unknown dip angles and directions for lineaments were estimated from the small-scale fracture sets. Simulated K _m and K _mb values range from 10^–8 to 10^–7?m/s and are greatest in a N–S direction, coinciding with the main strike direction of Okanagan Valley Fault Zone. K _mb values also decrease away from the fault, consistent with the decrease in lineament density. Simulated hydraulic conductivity values compare well with those estimated from pumping tests.     

Assessment of integrated electrical resistivity data on complex aquifer structures in NE Nuba Mountains – Sudan

Two geophysical techniques were integrated to map the groundwater aquifers on complex geological settings, in the crystalline basement terrain in northeast Nuba Mountains. The water flow is structurally controlled by the northwest–southeast extensional faults as one of several in-situ deformational patterns that are attributed to the collision of the Pan-African oceanic assemblage of the Nubian shield against the pre-Pan African continental crust to the west. The structural lineaments and drainage systems have been enhanced by the remote sensing technique. The geophysical techniques used are: vertical electrical soundings (VES) and electrical resistivity tomography (ERT), in addition to hydraulic conductivity measurements. These measurements were designed to overlap in order to improve the producibility of the geophysical data and to provide a better interpretation of the hydrogeological setting of the aquifer complex structure. Smooth and Block inversion schemes were attempted for the observed ERT data to study their reliability in mapping the different geometries in the complex subsurface. The VES data was conducted where ERT survey was not accessible, and inverted smoothly and merged with the ERT in the 3D resistivity grid. The hydraulic conductivity was measured for 42 water samples collected from the distributed dug wells in the study area; where extremely high saline zones were recorded and have been compared to the resistivity values in the 3D model.     

Integrated geoelectrical resistivity and hydrogeochemical methods for delineating and mapping heavy metal zone in aquifer system

A novel study on using geoelectrical resistivity, soil property, and hydrogeochemical analysis methods for delineating and mapping of heavy metal in aquifer system is presented in this paper. A total of 47 surveys of geoelectrical resistivity with Wenner configuration were conducted to determine the subsurface and the groundwater characteristics. The groundwater sample from 53 existing wells and 2 new wells has been analyzed to derive their water chemical content. The chemical analysis was done on the soil sample obtained from new two wells and from selected locations. The water and soil chemical analysis results from the new two wells were used as calibration in resistivity interpretation. The occurrence of heavy metal in aquifer system was expected to detect using the geoelectrical resistivity survey for the whole study area. The result of groundwater analysis shows that the groundwater sample contains a relatively low concentration of Fe (<?0.3 mg/L) elongating from the south up to the middle region. While in the middle and the northwestern, Fe concentration is relatively high (around 12 mg/L). Chemical analysis of soil sample shows that in the lower resistivity zone (<?18 Ωm), Al and Fe concentrations are comparatively high with an average of 68,000 and 40,000 mg/kg, respectively. Starting from the middle to the northwestern zone, the resistivity value appears to be low. It is definitely caused by higher Al and Fe concentration within the soil, and it is supported also by lower total anion content in the groundwater. While the resistivity value of more than 40 Ωm in aquifers is obtained in the zone which Fe concentration is relatively lower in the soil but not present in the groundwater. Correlation Fe concentration in the soil and Fe concentration in the groundwater sample shows the trend of positively linear; however, the Al concentration in soil has no correlation with Al content in groundwater. Finally, the probability of high heavy metal zone in the aquifer system is easily delineated by the distribution of geoelectrical resistivity presented in depth slice shapes which extend from the Boundary Range Composite Batholith in the north to the northwest.