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.     

Monitoring the impact of de-icing salt on roadside soils with time-lapse resistivity measurements

Monitoring systems along roads are needed to facilitate decisions on improving protection of water resources and decreasing the impact of road-related pollutants on the roadside environment. This paper presents a monitoring system using permanently installed electrodes and monthly measurements of resistivity at a motorway in Sweden with heavy loads of de-icing salt. A significant increase in resistivity in the vadose zone with increasing distance from the road was shown in both sand and glacial till areas during the whole year. By measuring temporal variations in a less affected environment further from the road, a distinction could be made between more natural variations and variations due to de-icing salt and melting of roadside snowbanks. The highest resistivities occurred in October–November and the lowest in January–March, while the more natural resistivities showed an opposing temporal variation. The difference was up to 35% on a log-scale in the sand area during the latter period. Hence, the time-lapse resistivity measurements clearly showed a strong influence of de-icing salt on roadside soils and groundwater during winter and spring. The measurement system and the analysis methods proved useful for monitoring both the spatial and seasonal variation in resistivity.     

Diesel transport monitoring in simulated unconfined aquifers using miniature resistivity arrays

The principal aim of this study is to assess the scope of monitoring diesel plume migration in a scaled aquifer model with a miniaturised electrical resistivity array. Respectively 1000 and 500 ml of diesel were injected in both the unsaturated and water-saturated zones of a sand body overlying a clay aquitard, and diesel migration was monitored with a miniature electrode array and an off-the-shelf resistivity meter. Inverted time-lapse electrical resistivity tomography (ERT) data reflect downward and lateral spreading of the diesel plume away from the injection point in the unsaturated zone. Diesel was also imaged to spread upwards and laterally away from the injection point in the saturated zone, as controlled by capillary rise. In both cases later-time ERT images reflected preferential pooling of diesel on the water table, as well as vertical smearing of pooled diesel in response to simulated water-table fluctuations. Repeat fluid electrical conductivity (EC) and dissolved oxygen (DO) measurements validate the observed changes in bulk resistivity caused by both diesel injections. Artefacts introduced by 2D inversion of 3D contaminant transport were abound. Time-lapse ERT imaging of diesel transport is therefore inferred to be feasible and well-suited to complementing conventional techniques of intrusive site investigation, although time-lapse 3D or 4D ERT imaging is strongly advocated.     

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.     

Monitoring of wastewater percolation in unsaturated sandy soil using geoelectrical measurements at Gabal el Asfar farm,northeast Cairo,Egypt

The use of wastewater for irrigation in sandy soil increases the pollution risk of the soil and may infiltrate to the shallow groundwater aquifer. In such environment, some important parameters need to be obtained for monitoring the wastewater in the unsaturated zone over the aquifer. These parameters include clay content, heterogeneities of the upper soils, depth to the aquifer and the variations of groundwater quality. In the present work, the efficiency of DC resistivity method in forms of 1-D and 2-D measurements was studied for wastewater monitoring in the Gabal el Asfar farm, northeast of Cairo, Egypt. Forty-one Schlumberger soundings (VES) were performed then followed by three pole-dipole 2-D profiles along some considered regions within the area. The resistivity measurements were integrated with the boreholes, hydrogeological and hydrochemical (surface and groundwater samples) information to draw a clear picture for the subsurface conditions. The obtained results were presented as cross sections and 3-D visualization to trace the clay intercalations within the unsaturated zone. In addition, a vulnerability map was created using the obtained results from 1-D Schlumberger survey and confirmed with the 2-D resistivity profiling. The obtained results have shown that the 2-D resistivity imaging technique is a powerful tool for mapping the small-scale variability within the unsaturated zone and the wastewater infiltration. However, limitations of resistivity techniques were observed in the area with limited resistivity contrast such as thin clay layers with brackish water background. Under that condition, the measured pattern of resistivity distributions depends on the applied electrode array, electrode spacing and using the available geological information during the inversion process.     

时移电阻率反演模拟研究

电阻率法可用于地表浅层的探测,也可用于对动态地下目标进行监测.对于监测数据的反演,不同数据集的单独反演存在一定缺陷,为此本文在常规电阻率反演算法的基础上,推导了时移电阻率反演公式,实现了时移反演算法程序;为了论证时移反演算法对动态地下目标成像的优越性,建立一组多个正演模型,利用模拟数据进行单独反演和时移反演,并对比两种...     

Mapping lithological variations in a river basin of West Bengal,India using electrical resistivity survey: implications for artificial recharge

Groundwater is a treasured earth’s resource and plays an important role in addressing water and environmental sustainability. However, its overexploitation and wide spatial variability within a basin and/or across regions are posing a serious challenge for groundwater sustainability. Some parts of southern West Bengal of India are problematic for groundwater occurrence despite of high rainfall in this region. Characterization of an aquifer in this area is very important for sustainable development of water supply and artificial recharge. Electrical resistivity surveys using 1-D and 2-D arrays were performed at a regular interval from Subarnarekha River at Bhasraghat (south) to Kharagpur (north) to map the lithological variations in this area. Resistivity sounding surveys were carried out at an interval of 2–3 km. Subsurface resistivity variation has been interpreted using very fast simulated annealing (VFSA) global optimization technique. The analysis of the field data indicated that the resistivity variation with depth is suitable in the southern part of the area and corresponds to clayey sand. Interpreted resistivity in the northern part of the area is relatively high and reveals impervious laterite layer. In the southern part of the area resistivity varies between 15 and 40 Ωm at a depth below 30 m. A 2-D resistivity imaging conducted at the most important location in the area is correlated well with the 1-D results. Based on the interpreted resistivity variation with depth at different locations different types of geologic units (laterite, clay, sand, etc.) are classified, and the zone of interests for aquifer has been demarcated. Study reveals that southern part of the area is better for artificial recharge than the northern part. The presence of laterite cover in the northern part of the area restricts the percolation of rainwater to recharge the aquifer at depth. To recharge the aquifer at depth in the northern part of the area, rainwater must be sent artificially at depth by puncturing laterite layers on the top. Such studies in challenging areas will help in understanding the problems and finding its solution.     

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.     

Geoelectric resistivity sounding of riverside alluvial aquifer in an agricultural area at Buyeo,Geum River watershed,Korea: an application to groundwater contamination study

Twenty profiles of vertical electric soundings (VES) were obtained in a riverside alluvium at the Buyeo area, South Korea, to examine the variations of subsurface geology and associated groundwater chemistry. The combination of the VES data with the borehole data provided useful information on subsurface hydrogeologic conditions. The vestige of an ancient river channel (e.g. oxbow lake) was identified on the resistivity profiles by the lateral continuation of a near-surface perched aquifer parallel to the river. Such a perched aquifer is typically developed in the area with a clay-rich silty surface alluvium which prohibits the infiltration of oxygen. Therefore, groundwater below the oxbow lake shows a very low nitrate concentration and Eh values under the strong anoxic condition. The distribution of water resistivity is correlated with that of measured total dissolved solids concentration in groundwater, while the earth resistivity of the aquifer shows a significant spatial variation. It is interpreted that the earth resistivity of the aquifer is mainly controlled by the soil type rather than by the water chemistry in the study area.     

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.     

Hydrogeophysical investigations in a typical khondalitic terrain to delineate the kaolinised layer using resistivity imaging

Investigation for high yielding water wells in the khondalitic terrain (graneti ferrous silliminite gnesiss) is mostly faced with the problem of identification of the extent of the depth of kaolinisation of the aquifer. The traditional Vertical Electrical Sounding survey, Seismic Refraction survey and Very Low Frequency Electromagnetic survey could not identify the kaolinisation of the aquifer in the present investigations. The Two Dimensional (2D) Resistivity and Induced Polarization (IP) Imaging surveys are attempted for the identification of kaolinised layer and depth of kaolinisation. Number of 2D Resistivity and IP Imaging profiles were conducted near Chipurupally in Vizianagaram district of Andhra Pradesh, India along successful and failed wells located within short distances. Resistivity and IP measurements were carried out using an ABEM SAS 4000 Terrameter. The resistivity and I.P. images have provided a clear view of the thickness of the highly weathered zone (kaolinised zone) at both successful and failed wells. The highly weathered zone is identified with the resistivity values below 25 ohm.m. The depth of highly weathered material at failed well is extended about 8–10 m more deeper than the successful wells at some places to as much as 20 m more deep at some other places. This extended deeper kaolinisation of the aquifer is responsible for failure of wells. Layers having resistivities between 25–65 Ohm.m are identified as aquifer layers which are composed of moderately weathered and fractured khondalitic suit of rocks (Garnti ferrous sillimanite/biotite gneiss). Layers with resistivities greater than 65 Ohm.m are interpreted to have basement characteristics belonging to the granite gneiss. Interestingly IP imaging has not provided any greater insights in delineating the kaolinistion of the aquifer when compared to resistivity Imaging, in fact resitivity imaging has shown greater depths of kaolinisation than IP Imaging.     

Integration of geoelectric and hydrochemical approaches for delineation of groundwater potential zones in alluvial aquifer

Geoelectric and hydrochemical approaches are employed to delineate the ground-water potential zones in District Okara, a part of Bari Doab, Punjab, Pakistan. Sixty-seven VES surveys are conducted with the Electrical Resistivity Meter. The resultant resistivity verses depth model for each site is estimated using computer-based software IX1D. Aquifer thickness maps and interpreted resistivity maps were generated from interpreted VES results. Dar-Zarrouk parameters, transverse resistance (TR), longitudinal conductance (SL) and anisotropy (λ) were also calculated from resistivity data to delineate the potential zones of aquifer. 70% of SL value is ≤3S, 30% of SL value is > 3S. According to SL and TR values, the whole area is divided into three potential zones, high, medium and low potential zones. The spatial distribution maps show that north, south and central parts of study area are marked as good potential aquifer zones. Longitudinal conductance values are further utilized to determine aquifer protective capacity of area. The whole area is characterized by moderate to good and up to some extent very good aquifer protective area on the basis of SL values. The groundwater samples from sixty-seven installed tube wells are collected for hydro-chemical analysis. The electrical conductivity values are determined. Correlation is then developed between the EC (μS/cm) of groundwater samples vs. interpreted aquifer resistivity showing R2 value 0.90.     

Comparison between 2D imaging and vertical electrical sounding in aquifer delineation: a case study of south and south west of Samawa City (IRAQ)

The vertical electrical sounding (VES) survey was carried out by using Schlumberger array in (12) points distributed along two profiles, except VES 6 was not located within any profile. The 2D imaging survey was carried out using Wenner array in four stations. The VES results revealed the presence of two groundwater aquifers. The first is Euphrates aquifer and the second is Dammam aquifer. While the results of 2D imaging survey distinguished two secondary aquifers in addition to previous aquifers located within Dammam formation. Also, 2D imaging gave a more accurate picture for the distribution of electrical horizons especially for the depths which ranges between 2.5 and 73 m. While the VES gives information for larger depths than that of 2D, which ranges between the earth’s surface and a depth 140 m. The geo-electrical sections of VES showed electrical horizons (layers) with sharp boundaries. But the 2D inverse models revealed the lateral and vertical variations of the resistivity within each horizon especially for shallow depths due to large volume of data available in 2D measurement. These variations were not shown in the geo-electrical sections of VES points. Comparison of VES and 2D imaging techniques revealed that the 2D imaging was the best for determined the shallow aquifers as observed from this study. It can be concluded that 2D resistivity imaging and VES can provide very useful guide for borehole drilling particularly where there is no existing well.     

Investigation of seawater intrusion in the Dibdibba Aquifer using 2D resistivity imaging in the area between Al-Zubair and Umm Qasr,southern Iraq

Electrical resistivity surveying for delineating seawater intrusion was performed in the Dibdibba aquifer in the area between the cities of Al-Zubair–Safwan and Al-Zubair–Umm Qasr in the vicinity of Khor AL-Zubair Channel, Basrah governorate, southern Iraq. Fourteen 2D resistivity profiles with a total length of 14 km were collected in the study area. The resistivity sections were compared with lithological data extracted from 11 boreholes. Thirty-nine groundwater samples were collected within the area and analyzed for chemical constituents; internal hydrogeological reports and unpublished studies were also evaluated. Results reveal the existence of three major resistivity layers, ranging from 0.1 to 130 Ωm at various depths and locations. The first layer has very low electrical resistivity (0.1–5 Ωm) representing a layer saturated with saltwater intruded from Khor AL-Zubair Channel. The second layer shows resistivity in the range of 5–130 Ωm, attributed to a transition zone and an unaffected zone saturated with brackish groundwater. The last resistivity layer (<?3 Ωm) represents coarse-grain sediments saturated with saline groundwater. Furthermore, a hard clay bed (Jojab) appears with a resistivity of 3–7 Ωm in all 2D imaging lines within a depth of 20–28 m. Electrical conductivity (EC) measurements from seven wells collected in 2014 and 2016 show a positive EC difference increasing landward with an average increase of 1927 µS/cm. In addition, six chemical relationships (Na/Cl, [Ca?+?Mg]/[HCO₃?+?SO₄], SO₄/HCO₃, SO₄/Cl, Mg/Ca and Cl/[HCO₃?+?CO₃]) are used to detect the source of salinity in groundwater. This study proves that extensive use of high-resolution 2D imaging sections, alongside lithological and hydrogeological data, can serve as a useful tool to delineate the boundaries between aquifers, identify hydraulic boundaries between groundwater with different salinities and allocate hard clay layers between the upper and lower Dibdibba aquifer. In general, the combination of 2D imaging and hydrochemistry enables conceptualization of the hydrogeological situation in the subsurface and characterization of the salinity source, here seawater intrusion, in the study area. There have been no studies published so far on the characteristics of saltwater intrusion in the study area, and this study is considered to be important for monitoring and studying the intrusion and regression of seawater.     

Estimation of hydraulic characteristics from electrical resistivity data in coastal aquifers of southern India

The applicability of geophysical techniques has been examined for evaluating aquifer properties like transmissivity and hydraulic conductivity of coastal aquifers, Tuticorin, Tamil Nadu. The pumping test data of 10 wells are interpreted by using forward modelling to obtain the aquifer characteristics in the study area. The available vertical electrical soundings (VES) data in the vicinity of the sites of pumping test have been interpreted; and true resistivity and thickness are determined at each site in the study area. Empirical relationships are established for estimating the hydraulic parameters from the electrical data.     

Application of electrical resistivity to map the stratigraphy and salinity of fluvio-deltaic aquifers: case studies from Bangladesh that reveal benefits and pitfalls

Fluvio-deltaic aquifers are the primary source of drinking water for the people of Bangladesh. Such aquifers, which comprise the Ganges-Brahmaputra-Meghna Delta, are hydrogeologically heterogeneous. Because of widespread groundwater quality issues in Bangladesh, it is crucial to know the hydrostratigraphic architecture and hydrochemistry, as some aquifer units are contaminated, whereas others are safe. Geophysical methods provide a potentially effective and noninvasive method for extensive characterization of these aquifers. This study applies and investigates the limitations of using electrical resistivity imaging (ERI) for mapping the hydrostratigraphy and salinity of an aquifer-aquitard system adjacent to the Meghna River. Some electrical resistivity (ER) sections showed excellent correlation between resistivity and grain size. These suggest that ERI is a powerful tool for mapping internal aquifer architecture and their boundaries with finer-grained aquitards which clearly appear as low-ER zones. However, in parts of some ER sections, variations in electrical properties were determined by porewater resistivity. In these cases, low ER was indicative of brine and did not indicate the presence of finer-grained materials such as silt or clay. Accordingly, the following hydrostratigraphic zones with different resistivities were detected: (1) aquifers saturated with fresh groundwater, (2) a regional silt/clay aquitard, and (3) a deeper brine-saturated formation. In addition, shallow silt/clay pockets were detected close to the river and below the vadose zone. ERI is thus a promising technique for mapping aquifers versus aquitards; however, the observations are easily confounded by porewater salinity. In such cases, borehole information and groundwater salinity measurements are necessary for ground-truthing.

     

A geoelectrical and hydrogeological study for the assessment of groundwater resources in Wadi Al Bih, UAE

Drilling information, historical water table levels, groundwater salinity records of the existing water wells in Wadi Al Bih area, United Arab Emirates, were stored in a geodatabase and used to characterize the geological and hydrogeological settings of this area. A 2D earth resistivity imaging survey was conducted for the first time in the Northern UAE to determine the potential of the Quaternary aquifer and its groundwater quality in the areas where there are no monitoring wells. The results of the chemical analyses of the collected groundwater samples together with the inversion results of the resistivity data were used to draw a total salinity map and determine the spatial variations in groundwater quality. The inversion results of the 2D earth resistivity imaging data indicated that the Quaternary aquifer in the study area is in a good connection with the underlying carbonate aquifer. It also indicated that the carbonate aquifer is of major regional and vertical extension and it contains the fresh water in this area. The data stored in the developed database were used to produce different types of geopotential maps.     

Using DC resistivity method to characterize the geometry and the salinity of the Plioquaternary consolidated coastal aquifer of the Mamora plain,Morocco

The Direct Current resistivity method was applied to the consolidated coastal Plioquaternary aquifer of Mamora plain, located on the Atlantic coast of Morocco. The aim is to determine the depth of the base of the saturated zone in the aquifer and to help in imaging lateral and vertical distribution of groundwater salinity. The geoelectric survey showed four geoelectric formations with the following electrical resistivities from top to bottom: 20–80, 200–2000, 200–300 and 5–70 Ohm m. The latter designates the basement of the aquifer constituted of marls and sandy marls. The mean resistivity of 250 Ohm m designates the aquifer formation. It decreases to less than 25% of its initial values for the soundings near the shore, reflecting the oceanic impact on the aquifer formation resistivity. The contour map shows that the basement of the saturated zone in the aquifer is deeper in the Northwest near the Sebou River estuary with values up to 70 m below sea level. This results in a larger thickness of the saturated zone of the aquifer leading to a consequent hydraulic potential. On the other hand, it has been deduced that the extent of marine intrusion inside the continent can be governed by human activities, natural properties and substratum geometry of the aquifer as well as by ecological factors. An optimal network of electrical soundings has been proposed for the monitoring of saltwater intrusion.     

基于二维NLCG反演的水文地质结构辨识研究

地下水调查是草原水资源合理开发利用的前提,而探明水文地质结构对地下水调查尤为重要.笔者充分利用大地电磁法(MT)和地面核磁共振(SNMR)各自的优点,尝试采用MT和SNMR对水文地质结构进行探测,并通过二维非线性共轭梯度(NLCG)反演解译了电性与地层两者的对应关系;实现了典型层状含水层模型和局部含水层模型正反演试验;...     

Hydrogeological modeling constraints provided by geophysical and geochemical mapping of a chlorinated ethenes plume in northern France

A methodological approach is described which combines geophysical and geochemical data to delineate the extent of a chlorinated ethenes plume in northern France; the methodology was used to calibrate a hydrogeological model of the contaminants’ migration and degradation. The existence of strong reducing conditions in some parts of the aquifer is first determined by measuring in situ the redox potential and dissolved oxygen, dissolved ferrous iron and chloride concentrations. Electrical resistivity imaging and electromagnetic mapping, using the Slingram method, are then used to determine the shape of the pollutant plume. A decreasing empirical exponential relation between measured chloride concentrations in the water and aquifer electrical resistivity is observed; the resistivity formation factor calculated at a few points also shows a major contribution of chloride concentration in the resistivity of the saturated porous medium. MODFLOW software and MT3D99 first-order parent–daughter chain reaction and the RT3D aerobic–anaerobic model for tetrachloroethene (PCE)/trichloroethene (TCE) dechlorination are finally used for a first attempt at modeling the degradation of the chlorinated ethenes. After calibration, the distribution of the chlorinated ethenes and their degradation products simulated with the model approximately reflects the mean measured values in the observation wells, confirming the data-derived image of the plume.