Observing lingering hyporheic storage using electrical resistivity: variations around stream restoration structures,Crabby Creek,PA

Time‐lapse geophysical surveys can map lingering hyporheic storage by detecting changes in response to saline tracer. Tracer tests were conducted in Crabby Creek, an urban stream outside Philadelphia, to examine the influence of stream restoration structures and variable sediment thickness. We compared electrical resistivity surveys with extensive well sampling (57 wells) in two 13.5‐m‐long reaches, each with a step drop created by a J‐hook. The two step drops varied in tracer behaviour, based on both the well data and the geophysical data. The well data showed more variation in arrival time where the streambed sediment was thick and was more uniform where sediment was thin. The resistivity in the reach with thin sediment showed lingering tracer in the hyporheic zone both upstream and downstream from the J‐hook. In the second reach where the sediment was thicker, the lingering tracer in the hyporheic zone was more extensive downstream from the J‐hook. The contrasting results between the two reaches from both methods suggested that sediments influenced hyporheic exchange more than the step at this location. Resistivity inversion differed from well data in both reaches in that it showed evidence for tracer after well samples had returned to background, mapping lingering tracer either upstream or downstream of a step. We conclude that resistivity surveys may become an important tool for hyporheic zone characterization because they provide information on the extent of slow moving fluids in the hyporheic zone, which have the potential to enhance chemical reactions. Copyright © 2012 John Wiley & Sons, Ltd.     

Delineating coastal groundwater discharge processes in a wetland area by means of electrical resistivity imaging, 224Ra and 222Rn

Coastal groundwater discharge (CGD) plays an important role in coastal hydrogeological systems as they are a water resource that needs to be managed, particularly in wetland areas. Despite its importance, identifying and monitoring CGD often presents physical and logistical constraints, restraining the application of more traditional submarine groundwater discharge surveying techniques. Here we investigate the capability of electrical resistivity imaging (ERI) in the Peníscola wetland (Mediterranean coast, Spain). ERI surveying made it possible to identify and delineate an ascending regional groundwater flow of thermal and Ra‐enriched groundwater converging with local flows and seawater intrusion. The continuous inputs of Ra‐rich groundwater have induced high activities of Ra isotopes and ²²²Rn into the marsh area, becoming among the highest previously reported in wetlands and coastal lagoons. Geoelectrical imaging enabled inferring focused upward discharging areas, leaking from the aquifer roof through a confining unit and culminating as spring pools nourishing the wetland system. Forward modelling over idealized subsurface configurations, borehole datasets, potentiometric records from standpipe piezometers, petrophysical analysis, and four natural and independent tracers (²²⁴Ra, ²²²Rn, temperature and salinity) permitted assessing the geoelectrical model and a derived hydrogeological pattern. The research highlights the potential of ERI to improve hydrogeological characterization of subsurface processes in complex contexts, with different converging flows. Additionally, a hydrogeological conceptual model for a groundwater‐fed coastal wetland was proposed, based on the integration of surveying datasets. Copyright © 2013 John Wiley & Sons, Ltd.     

Two-dimensional resistivity imaging and modelling in areas of complex geology

A system is described for the automatic measurement of electrical resistivity pseudo-sections. This comprises a linear array of up to 32 electrodes connected through a multicore cable to a computer controlled switching module and a resistivity meter. The processing of the measured sections to produce two-dimensional true resistivity images of the subsurface is briefly described. Some account is given of the capabilities and limitations of the technique. This is illustrated by a series of computed constant separation traverses for models of simple subsurface structures. Examples of processed images derived from sections measured in areas of relatively complex geology follow, a comparison being made of the interpretations obtained using an automatic imaging method and a manual iterative approach. It is concluded that with the equipment and software so far developed, in areas of modest subsurface geological complexity where some control is available and where the structures are essentially two-dimensional, then good approximations to the true geoelectric sections can be obtained down to depths of between 100 and 200 m.     

饱水和排水过程中岩石电阻率各向异性特征的电阻率成像法研究

从电阻率的角度研究岩石裂隙介质的各向异性特征是一种方便而有效的方法,但多限于空间单个点上的测量数据分析.通过在砂岩岩样上的饱水与排水实验以及同步进行的高密度电阻率成像监测,探讨了应用高密度电阻率成像法获得图像研究岩石各向异性特征的可能性,分析了饱水与排水过程中岩石电阻率在不同方向上的响应特性.结果表明,电阻率成像法在分析岩石裂隙介质的各向异性方面具有多方向成像和动态监测的优点,可以通过对不同方向上获得的电阻率图像的分析,提取出岩石沉积结构的空间分布模式,清晰地反映出岩石在饱水和排水过程中电阻率变化空间分布模式的各向异性特征.     

生活垃圾渗滤液污染砂土电阻率变化机制实验研究

为探讨含水多孔介质被垃圾渗滤液污染后的电阻率模型,实现基于土体电阻率的污染物量化分析,以不同粒径组成、不同初始含水量三种渗滤液污染砂土为样本进行了电阻率测试,基于实验结果对污染砂土电阻率变化机制进行了分析.结果显示污染砂土电阻率普遍随渗滤液侵入量增大呈指数降低;和相同含水程度未污染土相比,污染土电阻率值更低,降低幅度随初始含水量增大而减小.渗滤液污染导致的孔液电阻率变小、孔隙水含量增大是导电性增强的主导因素,与此同时污染过程也存在导致土体导电性降低的因素,降低因子受土的类型和初始含水量控制.涵盖这些因素的Archie型公式可作为渗滤液污染土的电性基本式,进而可利用这个基本式实现土体中污染物含量的量化分析.     

Electrical resistivity and seismic refraction studies for groundwater in the Jodhpura area,Rajasthan

Summary The results of geophysical investigations carried out for groundwater employing seismic refraction and electrical resistivity methods near Jodhpura (27^o4730N:75^o4240E) in the Jhunjhunu District, Rajasthan, are presented. Prominent resistivity low zones, revealed by the resistivity profiles, are indicative of deeper section of unconsolidated sediments overlying the bedrock and are usually favourable locations for groundwater. An equiresistivity and water-level contour map based on well data is furnished, indicating the fresh water zone boundary deduced geophysically. Longitudinal velocities for the various rock types occurring in and around the area are tabulated for general interest. A bedrock map prepared on the basis of depths estimated geophysically is incorporated. The depths to bedrock computed from electrical and/or seismic data compare well with those obtained from borehole data.     

Electrical resistivity imaging and hydrochemical analysis for groundwater investigation in Kuala Langat,Malaysia

ABSTRACT

Integrated two-dimensional electrical resistivity imaging (ERI) and hydrochemical surveys were used to investigate the groundwater alluvial aquifer in Kuala Langat, Malaysia. The study in the Langat basin considered the thickness of the aquifer, the depth of the bedrock, the regions influenced by seawater intrusion, and the monitoring of water levels. The resistivity imaging results show that the upper layer consists of clay, while the second layer is an aquifer whose thickness varies mostly in the range of 10–30 m, and in some cases extends to 40 m. The bedrock depth varies from 30 to 65 m. The chemical analyses were carried out on groundwater samples from nine boreholes collected between 2008 and 2012. The analyses indicate that the total dissolved solids (TDS) exceed 1000 mg L^-1 near the coastal area and are often less than 500 mg L^-1 further inland. The ERI and hydrochemical analyses reveal that groundwater in the study area, especially towards the coast, is a mixture of brackish and fresh waters.

EDITOR D. Koutsoyiannis; ASSOCIATE EDITOR M.D. Fidelibus     

Research Note: The mean sensitivity depth of the electrical resistivity method

Until now, a simple formula to estimate the depth of investigation of the electrical resistivity method that takes into account the positions of all of the electrodes for a general four‐electrode array has not been available. While the depth sensitivity function of the method for a homogeneous infinite half‐space is well known, previous attempts to use it to characterize the depth of investigation have involved calculating its peak and median, both of which must be determined numerically for a general four‐electrode array. I will show that the mean of the sensitivity function, which has not been considered previously, does admit a very simple mathematical formula. I compare the mean depth with the median and peak sensitivity depths for some common arrays. The mean is always greater than or equal to the median that is always greater than the peak. All three measures give reasonable estimates to the depths of actual structures for most circumstances. I will further show that, for 1D soundings, the use of the mean sensitivity depth as the pseudo‐depth assigns an apparent resistivity to a given pseudo‐depth that is consistent between different arrays. One consequence of this is that smoother depth soundings are obtained as “clutches,” caused by a change in the depth sensitivity due to moving the potential electrodes, are effectively removed. I expect that the mean depth formula will be a useful “rule of thumb” for estimating the depth of investigation before the resistivity structure of the ground is known.     

Effect of high-pressures on the electrical resistivity of natural zeolites from Deccan Trap, Maharashtra, India

We report here the electrical resistivity measurements on two natural zeolites–natrolite and scolecite (from the Killari borehole, Maharashtra, India) as a function of pressure up to 8 GPa at room temperature. High-pressure electrical resistivity studies on hydrous alumino-silicate minerals are very helpful in understanding the role of water in deep crustal conductivities obtained from geophysical models. The results obtained by magneto-telluric (MT) soundings and direct current resistivity surveys, along with the laboratory data on the electrical resistivity of minerals and rocks at high-pressure–temperature are used to determine the electrical conductivity distribution in continental lithosphere. The electrical resistivity of natural natrolite decreases continuously from 2.9 × 10⁹ Ω cm at ambient condition to 7.64 × 10² Ω cm at 8 GPa, at room temperature. There is no pressure-induced first order structural phase transitions in natrolite, when it is compressed in non-penetrating pressure transmitting medium up to 8 GPa. On the other hand scolecite exhibits a pressure-induced transition, with a discontinuous decrease of the electrical resistivity from 2.6 × 10⁶ to 4.79 × 10⁵ Ω cm at 4.2 to 4.3 GPa. The observed phase transition in scolecite is found to be irreversible. Vibrational spectroscopic and X-ray diffraction studies confirm the amorphous nature of the high-pressure phase. The results of the present high-pressure studies on scolecite are in good agreement with the high-pressure Raman spectroscopic data on scolecite. The thermo gravimetric studies on the pressure-quenched samples show that the samples underwent a pressure-induced partial dehydration. Such a pressure-induced partial dehydration, which has been observed in natural scolecite could explain the presence of high conductive layers in the earth's deep-crust.     

Some considerations on electrical resistivity imaging for characterization of waterbed sediments

The paper focuses on defining the performance and limits of ERI in the detection and sedimentary characterization of near-bottom thin layers. The analysis of the resolution of floating and submerged cables, and the effect of the accuracy of a priori information (resistivity and thickness) in the data inversion, is based on theory, models and actual data. Theoretical models show that the actual reconstruction of the near water-bottom sediments, in terms of geometry and resistivity, can be obtained only with the submerged cable, however, the data, unlike that acquired with the floating cable, require a priori information on water resistivity and thickness for the data inversion. Theoretical forward models based on wrong a priori water thickness and resistivity information influence the inverted model in different ways, depending on the under- and over-estimation of water resistivity and thickness, and the resistivity contrast of the water–solid layer; however a water–solid resistivity contrast of less than 2 and within 10% of error in water resistivity has no effect. Overestimating water resistivity depicts a ground similar to the actual ground in terms of resistivity, more so than the underestimation of water resistivity. Moreover, the data inversion is less influenced by water parameter error in the case of low resistivity contrast in the water–solid layer, than it is for high resistivity contrast. Wenner and Schlumberger arrays give comparable results, while a dipole–dipole array seems to be more sensitive to the accuracy of apparent resistivity measurements and a priori information on water.The theoretical considerations were validated by actual data acquired with a submerged cable on the Tiber River. The study has shown that if highly accurate measurements are made of water thickness and resistivity, then electrical resistivity imaging from the submerged cable can be used in addition to, or even to substitute, seismic data for the reconstruction of the features and sedimentary characterization of near-bed sediments where seismic data fail to give a suitable resolution.     

Determination of electrical resistivity,its anisotropy and heterogeneity on drill cores: a new method1

This paper presents a new method for the measurement on core samples of their electrical resistivity, its anisotropy and heterogeneity. The equipment used has been developed in the field laboratory of the German Continental Deep Drilling program KTB in the north-east of Bavaria on the western rim of the Bohemian Massif. The apparatus measures the resistivity at a fixed frequency as a function of the drill core azimuth and along the core by moving point electrode configurations. From these azimuth and depth dependences, mean values of resistivity and additional information about its anisotropy and heterogeneity are determined. Geometrical averaging is used, because the resistivity data follows a log normal distribution. The quantitative parameters ‘azimuth factor’, corresponding to horizontal anisotropy, and ‘heterogeneity factor’ are introduced. The depth logs of resistivity, azimuth factor and heterogeneity factor, measured on cores obtained from the KTB main drill hole (gneisses and amphibolites) at depths between 4150 m and 8080 m are presented. The geometrically averaged mean values of resistivity of gneisses and amphibolites are in the same range (? 10³Ωm). The resistivities tend to decrease with depth. The stress release of the drill cores during recovery produces microcracks which may partially account for this effect. Reduced resistivities (down to 150 Ωm) within an amphibolite core correlate with an alteration zone. One sample of this core displays alteration from fresh to completely altered. This sample is also electrically heterogeneous (heterogeneity factor ? 2). Other samples with uniform low alteration are more homogeneous heterogeneity factor ? 1.4). In general, higher anisotropies are observed in gneisses (mean azimuth factor 2.8), lower anisotropies in amphibolites (mean azimuth factor 1.3). Examples of isotropic and homogeneous samples, as well as anisotropic and heterogeneous samples are also presented.     

Agricultural and forested land use impacts on floodplain shallow groundwater temperature regime

Given the importance of groundwater temperature to the biogeochemical health of aquatic ecosystems, a floodplain study was implemented to improve understanding of rural land use impacts on shallow groundwater (SGW) temperature. Study sites included a historic agricultural field (Ag) and bottomland hardwood forest (BHF), each with nine piezometers in an 80 × 80 m grid. Piezometers were equipped with pressure transducers to monitor SGW temperature and level at 30 min intervals during the 2011, 2012, 2013, and 2014 water years. The study is one of the first to utilize long‐term, continuous, automated, in situ monitoring to investigate rural land use impacts on shallow groundwater temperatures. Average SGW temperature during the study period was 11.1 and 11.2 °C at the Ag and BHF sites, respectively. However, temperature range at the Ag site was 72% greater than at the BHF site. Results indicate a greater responsiveness to seasonal climate fluctuations in Ag site SGW temperature related to absence of forest canopy. Patterns of intra‐site groundwater temperature differences at both study sites illustrate the influence of stream–aquifer thermal conduction and occasional baseflow reversals. Considering similar surface soil temperature amplitudes and low average groundwater flow values at both sites, results suggest that contrasting rates of plant water use, groundwater recharge, and subsurface hydraulic conductivity are likely mechanistic causes for the observed SGW temperature differences. Results highlight the long‐term impact of forest removal on subsurface hydrology and groundwater temperature regime. Copyright © 2015 John Wiley & Sons, Ltd.     

The inversion of data from very large three-dimensional electrical resistivity tomography mobile surveys

New developments in mobile resistivity meter instrumentation have made it possible to survey large areas with dense data coverage. The mobile system usually has a limited number of electrodes attached to a cable that is pulled along behind an operator so that a large area can be covered within a short time. Such surveys can produce three-dimensional datasets with hundreds of thousands of electrodes positions and data points. Similarly, the inverse model used to interpret the data can have several hundred thousand cells. It is impractical to model such large datasets within a reasonable time on microcomputers used by many small companies employing standard inversion techniques. We describe a model segmentation technique that subdivides the finite-element mesh used to calculate the apparent resistivity and Jacobian matrix values into a number of smaller meshes. A fast technique that optimizes the calculation of the Jacobian matrix values for multi-channel systems was also developed. A one-dimensional wavelet transform method was then used to compress the storage of the Jacobian matrix, in turn reducing the computer time and memory required to solve the least-squares optimization equation to determine the inverse model resistivity values. The new techniques reduce the calculation time and memory required by more than 80% while producing models that differ by less than 1% from that obtained using the standard inversion technique with a single mesh. We present results using a synthetic model and a field dataset that illustrates the effectiveness of the proposed techniques.     

Correlation of vertical electrical sounding and electrical borehole log data for groundwater exploration

We have correlated the longitudinal unit conductance CL obtained from interpreted vertical electrical sounding data with the formation resistivity R_t and the formation resistivity factor F, obtained by carrying out electrical borehole logging. Interpreted geophysical data of eleven soundings and two electrical borehole log records are used for the analysis. The geophysical data used were acquired in a sedimentary basin. The study area is called Lower Maner Basin located in the province of Andhra Pradesh, India. Vertical electrical soundings were carried out using a Schlumberger configuration with half current electrode separation varying from 600–1000 m. For logging the two boreholes, a Widco logger‐model 3200 PLS was used. True formation resistivity R_t was calculated from a resistivity log. Formation resistivity factor F was also calculated at various depths using R_t values. An appreciable inverse relation exists between the correlated parameters. The borehole resistivity R_t and the formation resistivity factor F decrease with the increase in the longitudinal unit conductance CL. We have shown the use of such a relation in computing borehole resistivity R_t and formation resistivity factor F at sites that posses only vertical electrical sounding data, with a fair degree of accuracy. Validation of the correlation is satisfactory. Scope for updating the correlation is discussed. Significance and applications of the relation for exploration of groundwater, namely to update the vertical electrical sounding data interpretation by translating the vertical electrical sounding data into electrical borehole log parameters, to facilitate correlations studies and to estimate the porosity (φ), permeability (K) and water saturation S_w of water bearing zones are discussed.     

基于三维高密度电法的地质BIM模型应用研究

地质体三维可视化BIM(Building Information Model)技术具有直观呈现地质体的三维赋存情况、实现地质体的空间计算与分析等众多功能,为工程设计和施工提供可视化的基础操作平台.目前,地质BIM模型的构建主要依赖大量钻孔资料和地质测绘成果.而由这些点状信息构建的三维地质模型易出现孔间盲区,且大量钻孔费时费力、破坏三维地质体结构,难以满足精细化、快捷、无损的地质建模需求.针对此问题,本文以豫北某拟建水库大坝工程为依托,利用三维高密度电法对坝址区左岸杂填土进行探测,物探解译成果结合少量钻孔资料快速精确构建工区地质模型.在此基础上,将大坝输水洞设计参数引入模型,直观展示沿洞线地质体与输水洞的空间关系.通过该方法计算施工开挖土方量比传统地质断面法多6000 m3,系孔间盲区基岩面下降形成基岩凹槽所致.研究表明,三维高密度电法成果结合少量钻孔资料可快速实现三维地质BIM精细化建模,为工程后期预算、施工地质预报以及构筑物优化设计提供可靠地质依据,也为精细三维地质BIM建模提供了一种新思路.     

地震前后祁连山中段深部电性变化的跟踪研究

利用在青藏高原北缘第一个长约４００ｋｍ的ＭＴ重复监测剖面上所获得的地震前后最新的地表视电阻率资料，用一维ＭＴ反演方法得到了剖面古方电性结构随时间的变化。着重分析相当于大陆多震层深度范围内介质电阻率参数在地震前后的变化。结果表明，地下电阻率的变化呈现分区特征，并与剖面附近地震的强度及位置有关。     

The use of non-conventional CPTe data in determination of 3-D electrical resistivity distribution

The spatial distribution of the electrical resistivity data provides useful information for investigating and modeling the fluid transport processes. 3D electrical resistivity distribution provides information about water flow and changes in electrical resistivity of the pore fluid.Therefore, to assist in understanding and modeling of the fluid transport process, 3D spatial distribution of the electrical resistivity data with the corresponded 3D geological section were mapped and interpreted in the test site located in western Germany. A process of deriving electrical resistivity values from the mechanical and radioactive parameters of cone penetration tests (CPT) and geological information of boreholes was presented. A reliable method which gives accurate resistivity values in cases of near surface sediments was introduced. Then a field test was executed where the calculated resistivity values were compared with the measured CPTe resistivity data. The CPTe (cone penetration test with electrical extension) data were also used in correlating to the ERT (electrical resistivity tomography) data. Consequently, obtained dense CPT surveys give us the possibility to determine a high resolution resistivity distribution of the investigated area.     

基于不等式约束的三维电阻率探测混合反演方法

三维电阻率探测的线性反演和非线性反演中均存在着多解性的固有难题.电阻率线性反演方法的效率较高,但反演结果对初始模型的依赖性较强,易陷入局部极小;而非线性反演方法不依赖初始模型,但搜索效率极低,尚未见到关于三维电阻率非线性反演的文献.针对上述问题,融合线性与非线性反演方法的互补优势,提出了最小二乘法（线性方法）与改进遗传算法（非线性方法）相结合的混合反演方法的概念和思想.首先,提出了将介质电阻率变化范围作为不等式约束引入反演方程的思路,以实现压制多解性、提高可靠性的目标.提出了宽松不等式约束和基于钻孔推断的局部严格不等式约束的获取及定义方法.在此基础上,分别提出了基于不等式约束的最小二乘线性反演方法和遗传算法非线性反演方法.其次,对于遗传算法在变异搜索方向控制、初始群体产生等方面进行了改进,优化了其搜索方向和初始群体多样性.然后,提出了混合反演方法及其实现方案,利用改进遗传算法进行第一阶段反演,发挥其对初始模型的依赖程度低的优势,搜索到最优解附近的空间,输出当前最优个体;利用最小二乘法进行第二阶段反演,将遗传算法得到的当前最优个体作为初始模型,在最优解附近空间执行高效率的局部线性搜索,最终实现地电结构的三维成像.最后,开展了合成数据与实际工程算例验证,与传统最小二乘方法进行了对比,发现混合反演方法在压制多解性、摆脱初始模型依赖和提高反演效果方面有较好效果.

     

A watershed-scale study of climate change impacts on groundwater recharge (Annapolis Valley,Nova Scotia,Canada)

Abstract

The potential impacts of future climate change on the evolution of groundwater recharge are examined at a local scale for a 546-km² watershed in eastern Canada. Recharge is estimated using the infiltration model Hydrologic Evaluation of Landfill Performance (HELP), with inputs derived from five climate runs generated by a regional climate model in combination with the A2 greenhouse gas emissions scenario. The model runs project an increase in annual recharge over the 2041–2070 period. On a seasonal basis, however, a marked decrease in recharge during the summer and a marked increase during the winter are observed. The results suggest that increased evapotranspiration resulting from higher temperatures does not offset the large increase in winter infiltration. In terms of individual water budget components, clear differences are obtained for the different climate change scenarios. Monthly recharge values are also found to be quite variable, even for a given climate scenario. These findings are compared with results from two regional-scale studies.

Editor D. Koutsoyiannis; Associate editor M. Besbes     

A comparison of peak and plate electrodes in electrical resistivity tomography: application to the chalky groundwater of the Beauvais aquifer (northern part of the Paris basin,France)

Electrical resistivity tomography was used in order to explore an experimental site of the LaSalle Beauvais Polytechnic Institute (France). The test was conducted along a profile line of 315 m length, using 64 electrodes deployed at an inter‐electrode spacing of 5 m, and the data were recorded using gradient, Wenner and pole–dipole arrays. The performance of plate electrodes (non‐conventional flat‐based) is compared with the performance of peak electrodes (conventional spike). The hydrogeophysical investigation of the chalk aquifer system of Beauvais shows that the performance of plate electrodes is satisfactory, leading to inversions of small root‐mean‐square errors. Peak and the plate electrodes were tested before and after injection of a salt tracer in the piezometer of the experimental site. The study demonstrates the usefulness of plate electrodes (efficient, less time consuming) and the possibility of aquifer characterization by a salt tracer. Copyright © 2010 John Wiley & Sons, Ltd.