Geophysical Investigation of Seepage Beneath an Earthen Dam

A hydrogeophysical survey is performed at small earthen dam that overlies a confined aquifer. The structure of the dam has not shown evidence of anomalous seepage internally or through the foundation prior to the survey. However, the surface topography is mounded in a localized zone 150 m downstream, and groundwater discharges from this zone periodically when the reservoir storage is maximum. We use self‐potential and electrical resistivity tomography surveys with seismic refraction tomography to (1) determine what underlying hydrogeologic factors, if any, have contributed to the successful long‐term operation of the dam without apparent indicators of anomalous seepage through its core and foundation; and (2) investigate the hydraulic connection between the reservoir and the seepage zone to determine whether there exists a potential for this success to be undermined. Geophysical data are informed by hydraulic and geotechnical borehole data. Seismic refraction tomography is performed to determine the geometry of the phreatic surface. The hydro‐stratigraphy is mapped with the resistivity data and groundwater flow patterns are determined with self‐potential data. A self‐potential model is constructed to represent a perpendicular profile extending out from the maximum cross‐section of the dam, and self‐potential data are inverted to recover the groundwater velocity field. The groundwater flow pattern through the aquifer is controlled by the bedrock topography and a preferential flow pathway exists beneath the dam. It corresponds to a sandy‐gravel layer connecting the reservoir to the downstream seepage zone.     

Self-potential investigation of moraine dam seepage

Self-potential (SP) and electrical resistivity measurements are used to investigate seepage at a remote moraine dam in the Sierra Nevada of California. The site is a small terminal moraine impounding roughly 300,000 m³ of water at ~ 3400 m a.s.l. Suspicious fine sediment in a small lake at the dam's downstream toe prompted initial concerns that anomalous seepage may be eroding matrix material from the moraine. 235 individual SP measurements covering the surface of the dam were collected in order to investigate electrokinetic current sources resulting from seepage, while resistivity soundings probed moraine stratigraphy and suggest that the till contains interstitial ice. Contoured SP data reveal a non-uniform voltage distribution over the moraine dam and two distinct negative SP anomalies. The first, located in the central area of the moraine, shows a broad negative SP zone around the crest and increasingly positive SP moving downhill towards both the upstream and downstream toes. This anomaly can be explained by shallow gravitational groundwater flow in the near subsurface combined with upward groundwater flux through evapotranspiration; numerical simulation of the combined effect matches field data well. The second SP anomaly has a tightly localized distribution and can be explained by vertically descending flow into a bedrock fault conduit. Our conceptual seepage model suggests that flow travels from Dana Lake first at the boundary of ice-filled moraine and bedrock before converging on a concentrated channel in the subvertical fault zone. Positive SP near the dam abutments results from groundwater inflow from adjacent hillslopes. Combined analyses suggest that seepage erosion is not currently affecting the moraine dam, and that the sediment observed on the bed of the downstream toe lake is likely a remnant of past outflow events.     

SPONTANEOUS POTENTIALS AND RESISTIVITY SURVEYS APPLIED TO HYDROGEOLOGY IN A VOLCANIC AREA: CASE HISTORY OF THE CHAÎNE DES PUYS (PUY-DE-DÔME,FRANCE)1

Two geophysical methods [resistivity soundings and spontaneous potentials (SP)] are used to investigate aquifers in the Quaternary volcanic formations of the Chaîne des Puys (France). The interpretation of the resistivity soundings required a specific work of determination of the resistivity of the formations concerned. The basement topography, characterized by axial horst and graben structures and perpendicular palaeovalleys, was revealed. SP was developed experimentally in the field and theoretically. It is shown that the SP anomalies can be considered as double-layer potentials, with their source on the aquifers' interfaces. The most significant of them is the water table, which creates an SP anomaly of similar but inverted form. SP anomalies due to flows in unconfined aquifers, in palaeovalleys, and in volcanic ranges are computed and compared with the observed ones. These two methods have determined the extent and the boundary conditions of the different hydrogeological basins and have determined the major drainage axes and the groundwater flow pattern.     

Flow resistance equation for rills

In this paper, a new flow resistance equation for rill flow was deduced applying dimensional analysis and self‐similarity theory. At first, the incomplete self‐similarity hypothesis was used for establishing the flow velocity distribution whose integration gives the theoretical expression of the Darcy–Weisbach friction factor. Then the deduced theoretical resistance equation was tested by some measurements of flow velocity, water depth, cross section area, wetted perimeter, and bed slope carried out in 106 reaches of some rills shaped on an experimental plot. A relationship between the velocity profile, the channel slope, and the flow Froude number was also established. The analysis showed that the Darcy–Weisbach friction factor can be accurately estimated by the proposed theoretical approach based on a power–velocity profile.     

Groundwater occurrence and flow pattern in the Enugu coal-mine area,Anambra State,Nigeria

Abstract

The distribution of major geological units, static water level data, water chemistry data, and observations of surface features influenced by groundwater seepage were used to ascertain the nature of groundwater occurrence and flow pattern in the Enugu coal field, Nigeria. Considerations of the geological units, the static groundwater levels and groundwater seepages in the mines indicate that the coal sequence is a multiaquifer system in which sandstone and coal aquifers alternate with shale aquitards. Based on the hydraulic head data, the groundwater flow is predominantly downwards. Groundwater velocity calculation across the multiaquifer system using the Darcy equation gave a flow velocity of about 1 m day^?1. For groundwater systems, such a calculated velocity is considered high. The high velocity is most probably due to the high fracture porosity as well as the presence of other stratigraphic and structural features such as alluvial fills that provide high hydraulic conductivity pathways across the aquifer system. The pattern of groundwater inflow into the mines is also influenced by these stratigraphie and structural features.     

Investigation of Seepage Meter Measurements in Steady Flow and Wave Conditions

Water exchange between surface water and groundwater can modulate or generate ecologically important fluxes of solutes across the sediment‐water interface. Seepage meters can directly measure fluid flux, but mechanical resistance and surface water dynamics may lead to inaccurate measurements. Tank experiments were conducted to determine effects of mechanical resistance on measurement efficiency and occurrence of directional asymmetry that could lead to erroneous net flux measurements. Seepage meter efficiency was high (average of 93%) and consistent for inflow and outflow under steady flow conditions. Wave effects on seepage meter measurements were investigated in a wave flume. Seepage meter net flux measurements averaged 0.08 cm/h—greater than the expected net‐zero flux, but significantly less than theoretical wave‐driven unidirectional discharge or recharge. Calculations of unidirectional flux from pressure measurements (Darcy flux) and theory matched well for a ratio of wave length to water depth less than 5, but not when this ratio was greater. Both were higher than seepage meter measurements of unidirectional flux made with one‐way valves. Discharge averaged 23% greater than recharge in both seepage meter measurements and Darcy calculations of unidirectional flux. Removal of the collection bag reduced this net discharge. The presence of a seepage meter reduced the amplitude of pressure signals at the bed and resulted in a nearly uniform pressure distribution beneath the seepage meter. These results show that seepage meters may provide accurate measurements of both discharge and recharge under steady flow conditions and illustrate the potential measurement errors associated with dynamic wave environments.     

A Tracer Test to Characterize Treatment of TCE in a Permeable Reactive Barrier

A tracer test was conducted to characterize the flow of groundwater across a permeable reactive barrier constructed with plant mulch (a biowall) at the OU‐1 site on Altus Air Force Base, Oklahoma. This biowall is intended to intercept and treat groundwater contaminated by trichloroethylene (TCE) in a shallow aquifer. The biowall is 139‐m long, 7.3‐m deep, and 0.5‐m wide. Bromide was injected from an upgradient well into the groundwater as a conservative tracer, and was subsequently observed breaking through in monitoring wells within and downgradient of the biowall. The bromide breakthrough data demonstrate that groundwater entering the biowall migrated across it, following the slope of the local groundwater surface. The average seepage velocity of groundwater was approximately 0.06 m/d. On the basis of the Darcy velocity of groundwater and geometry of the biowall, the average residence time of groundwater in the biowall was estimated at 10 d. Assuming all TCE removal occurred in the biowall, the reduction in TCE concentrations in groundwater across the biowall corresponds to a first‐order attenuation rate constant in the range of 0.38 to 0.15 per d. As an independent estimate of the degradation rate constant, STANMOD software was used to fit curves through data on the breakthrough of bromide and TCE in selected wells downgradient of the injection wells. Best fits to the data required a first‐order degradation rate constant for TCE removal in the range of 0.13 to 0.17 per d. The approach used in this study provides an objective evaluation of the remedial performance of the biowall that can provide a basis for design of other biowalls that are intended to remediate TCE‐contaminated groundwater.     

SNMR联合SP用于滑动带水体赋存状态的探测研究

地下水的空间分布及运移规律是建立地下水流系统和分析渗流场规律的重要因素.具有直接找水特性的地面核磁共振方法(SNMR)可以用于确定地下水的空间分布,而自然电场法(SP)是一种可以探测地下水流动方向的物探方法,这两种方法的联合使用为地下水研究提供了一种较好的工作模式.本文以滑动带水体探测为例,阐述该工作模式的应用效果.首...     

Improved imaging of a karst aquifer using focused source electromagnetic and differentially normalized method: a qualitative analysis

Geophysical investigations using conventional techniques applied to groundwater exploration can often present strong limitations involving high financial costs, complex acquisition logistics and high ambiguity in results. Dispersion of the electric current flow, induced polarization) effects, cultural noises and shallow lateral heterogeneities represent the main problems faced by geoelectric methods in these types of surveys. Moreover, elements such as intrusions and mineralization at different depths may be responsible for signal attenuation as well as high resistivity in unsaturated zones and complex three-dimensional formations or clayey zones cause variations in the electric current. The focused source electromagnetic and differentially normalized method approaches can help to solve some these issues. Aiming at a higher signal-to-noise ratio, the focused source electromagnetic method and approaches of the differentially normalized method, first applied to petroleum exploration, are tested on a groundwater target, in a karst environment sectioned by a diabase dyke. We performed the processing and analysis on real IP resistivity profiling data acquired with two-way dipole-dipole array, guided by magnetic data acquired on the same profile, mapping a diabase dyke. The inversion of focused source electromagnetic method/differentially normalized method was not performed, instead that we converted the induced polarization–resistivity data to a differential signal to qualitatively prove the presence of aquifer. Joint interpretation of focused source electromagnetic method curves and inverted two-dimensional induced polarization–resistivity sections allowed for precise delineation of a conductive zone associated with the karst aquifer, le magnetics allowed for the definition of a neighbour dyke. The techniques have great potential in the aid of groundwater exploration, contributing substantially to the reduction of interpretation ambiguity. Focused source electromagnetic method/differentially normalized method/ approaches show that a simple linear combination of the conventional geoelectric data is able to remove the geological noise and provide the vertical focusing of the electric current.     

Geophysical and Hydrologic Studies of Lake Seepage Variability

Variations in lake seepage were studied along a 130 m shoreline of Mirror Lake NH. Seepage was downward from the lake to groundwater; rates measured from 28 seepage meters varied from 0 to ?282 cm/d. Causes of this variation were investigated using electrical resistivity surveys and lakebed sediment characterization. Two‐dimensional (2D) resistivity surveys showed a transition in lakebed sediments from outwash to till that correlated with high‐ and low‐seepage zones, respectively. However, the 2D survey was not able to predict smaller scale variations within these facies. In the outwash, fast seepage was associated with permeability variations in a thin (2 cm) layer of sediments at the top of the lakebed. In the till, where seepage was slower than that in the outwash, a three‐dimensional resistivity survey mapped a point of high seepage associated with heterogeneity (lower resistivity and likely higher permeability). Points of focused flow across the sediment–water interface are difficult to detect and can transmit a large percentage of total exchange. Using a series of electrical resistivity geophysical methods in combination with hydrologic data to locate heterogeneities that affect seepage rates can help guide seepage meter placement. Improving our understanding of the causes and types of heterogeneity in lake seepage will provide better data for lake budgets and prediction of mass transfer of solutes or contaminants between lakes and groundwater.     

Out‐of‐plane effects in 2D borehole‐to‐surface resistivity tomography and applications in mineral exploration

In this paper, we discuss the effects of anomalous out‐of‐plane bodies in two‐dimensional (2D) borehole‐to‐surface electrical resistivity tomography with numerical resistivity modelling and synthetic inversion tests. The results of the two groups of synthetic resistivity model tests illustrate that anomalous bodies out of the plane of interest have an effect on two‐dimensional inversion and that the degree of influence of out‐of‐plane body on inverted images varies. The different influences are derived from two cases. One case is different resistivity models with the same electrode array, and the other case is the same resistivity model with different electrode arrays. Qualitative interpretation based on the inversion tests shows that we cannot find a reasonable electrode array to determine the best inverse solution and reveal the subsurface resistivity distribution for all types of geoelectrical models. Because of the three‐dimensional effect arising from neighbouring anomalous bodies, the qualitative interpretation of inverted images from the two‐dimensional inversion of electrical resistivity tomography data without prior information can be misleading. Two‐dimensional inversion with drilling data can decrease the three‐dimensional effect. We employed two‐ and three‐dimensional borehole‐to‐surface electrical resistivity tomography methods with a pole–pole array and a bipole–bipole array for mineral exploration at Abag Banner and Hexigten Banner in Inner Mongolia, China. Different inverse schemes were carried out for different cases. The subsurface resistivity distribution obtained from the two‐dimensional inversion of the field electrical resistivity tomography data with sufficient prior information, such as drilling data and other non‐electrical data, can better describe the actual geological situation. When there is not enough prior information to carry out constrained two‐dimensional inversion, the three‐dimensional electrical resistivity tomography survey is the better choice.     

The Role of Groundwater for Lake‐Water Quality and Quantification of N Seepage

The heterogeneous nature of both groundwater discharge to a lake (inflow) and nitrate concentrations in groundwater can lead to significant errors in calculations of nutrient loading. Therefore, an integrated approach, combining groundwater flow and transport modelling with observed nitrate and ammonium groundwater concentrations, was used to estimate nitrate loading from a catchment via groundwater to an oligotrophic flow‐through lake (Lake Hampen, Denmark). The transport model was calibrated against three vertical nitrate profiles from multi‐level wells and 17 shallow wells bordering a crop field near the lake. Nitrate concentrations in groundwater discharging to the lake from the crop field were on average 70 times higher than in groundwater from forested areas. The crop field was responsible for 96% of the total nitrate loading (16.2 t NO₃/year) to the lake even though the field only covered 4.5% of the catchment area. Consequently, a small change in land use in the catchment will have a large effect on the lake nutrient balance and possible lake restoration. The study is the first known attempt to estimate the decrease of nitrate loading via groundwater to a seepage lake when an identified catchment source (a crop field) is removed.     

Real-Time Geoelectric Monitoring of Seepage into Sand and Clay Layer

The migration of groundwater in rock and soil can appear as abnormalities in geoelectric fields. It is therefore important to study the characteristics of seepage in porous media by measuring the geoelectric field signatures. In this study, a physical model with layers of sand and clay was constructed and an electrical resistivity meter was used to examine the changes in the geoelectric field parameters during the infiltration process. The results show that the infiltration could be detected based on the geoelectric signatures including temporal changes through the spontaneous potential, excitation currents, and apparent resistivity. Specifically, the spontaneous potential was reduced by 100 to 200 mV when the water reaches an electrode. During the second water injection in the experiment, the measured spontaneous potential of all the electrodes recovered to the previous extreme values that range from −200 to −550 mV, thus indicating a “memory” effect. With stepwise changes in the excitation current, it was possible to determine the seepage velocity in sand and clay layer. The apparent resistivity is reduced to less than 400 Ωm when the infiltration reaches the electrodes. These results indicate the potential for real-time monitoring of water flow.     

Missing role of groundwater in water and nutrient cycles in the shallow eutrophic lake Kasumigaura,Japan

A process‐based model was developed, NICE‐LAKE (NIES (National Institute for Environmental Studies) Integrated Catchment‐based Ecohydrology), which includes interactions between surface water, canopy, unsaturated water, aquifer, lake and rivers, and used it to model the shallow eutrophic Lake Kasumigaura in Japan. By estimating the spatial distribution of the hydrological cycle, the model shows that groundwater withdrawal greatly affects groundwater distribution and seepage and indirectly influences lake water level. The simulated seepage agrees excellently with the budget‐derived value calculated from the observed groundwater level, lake level and isotope analyses. The model showed that groundwater seepage and groundwater concentrations are important contributors to the eutrophication of Lake Kasumigaura, an important contribution not recognized in previous studies of the lake. Groundwater entering the lake from the north side is contaminated with high concentrations of nitrate and ammonia from intensive pig and cattle raising and cultivated fields. The simulation showed that this high nitrogen load plays an important role in the eutrophication of the lake (the nitrogen load in inflowing groundwater is 30% of river inflow and 4 times that from wastewater treatment plants) in spite of government policies to prevent overland flow of nutrients into the lake. Our results show that NICE‐LAKE is a powerful tool for forecasting how the water quality of the lake will be affected by the (illegal) disposal of animal excreta in the surrounding open fields. Copyright © 2007 John Wiley & Sons, Ltd.     

Hydrodynamic and isotopic investigations for evaluating the mechanisms and amount of groundwater seepage through a rockslide dam

This study addresses the influence of landslide dams on surface water drainage and groundwater flow. In the study area of Scanno Lake and Sagittario River (Central Italy), a limestone rockslide‐avalanche formed a lake, which has an outlet that is occasionally active, showing infiltration into the rockslide dam. Several springs are present at the lake's base and are partly fed by seepage through the rockslide debris. Piezometric surveys, discharge measurements, pumping tests and chemical analyses are tools used to build a conceptual model of the groundwater flow and to evaluate the flow through the rockslide debris. Seasonal water isotopic signatures validate the assumed model, showing a mixing of infiltration recharge and groundwater seepage throughout the rockslide debris. Various recharge areas have been found for springs, pointing out those directly fed by the rockslide debris aquifer. Hypotheses about seasonal groundwater mixing between the regional carbonate aquifer and the rockslide debris aquifer are supported by isotope results. Seasonal changes in groundwater table level due to recharge and surface losses from seasonal outlet have been correlated with isotopic groundwater composition from the rockslide debris aquifer and the downstream springs; this relationship highlights the role of the rockslide dam body on the hydrodynamics of the studied area. Relationships between surface waters and groundwater in the area have been completely understood on the basis of water isotopic fingerprinting, finally obtaining a complete evaluation of groundwater renewable resources and its regimen. Copyright © 2010 John Wiley & Sons, Ltd.     

Upscaling Point Velocity Measurements to Characterize a Glacial Outwash Aquifer

Small‐scale point velocity probe (PVP)‐derived velocities were compared to conventional large‐scale velocity estimates from Darcy calculations and tracer tests, and the possibility of upscaling PVP data to match the other velocity estimates was evaluated. Hydraulic conductivity was estimated from grain‐size data derived from cores, and single‐well response testing or slug tests of onsite wells. Horizontal hydraulic gradients were calculated using 3‐point estimators from all of the wells within an extensive monitoring network, as well as by representing the water table as a single best fit plane through the entire network. Velocities determined from PVP testing were generally consistent in magnitude with those from depth specific data collected from multilevel monitoring locations in the tracer test, and similar in horizontal flow direction to the average hydraulic gradient. However, scaling up velocity estimates based on PVP measurements for comparison with site‐wide Darcy‐based velocities revealed issues that challenge the use of Darcy calculations as a generally applicable standard for comparison. The Darcy calculations were shown to underestimate the groundwater velocities determined both by the PVPs and large‐scale tracer testing, in a depth‐specific sense and as a site‐wide average. Some of this discrepancy is attributable to the selective placement of the PVPs in the aquifer. Nevertheless, this result has important implications for the design of in situ treatment systems. It is concluded that Darcy estimations of velocity should be supplemented with independent assessments for these kinds of applications.     

Imaging Hydrological Processes in Headwater Riparian Seeps with Time‐Lapse Electrical Resistivity

Delineating hydrologic and pedogenic factors influencing groundwater flow in riparian zones is central in understanding pathways of water and nutrient transport. In this study, we combined two‐dimensional time‐lapse electrical resistivity imaging (ERI) (depth of investigation approximately 2 m) with hydrometric monitoring to examine hydrological processes in the riparian area of FD‐36, a small (0.4 km²) agricultural headwater basin in the Valley and Ridge region of east‐central Pennsylvania. We selected two contrasting study sites, including a seep with groundwater discharge and an adjacent area lacking such seepage. Both sites were underlain by a fragipan at 0.6 m. We then monitored changes in electrical resistivity, shallow groundwater, and nitrate‐N concentrations as a series of storms transitioned the landscape from dry to wet conditions. Time‐lapse ERI revealed different resistivity patterns between seep and non‐seep areas during the study period. Notably, the seep displayed strong resistivity reductions (～60%) along a vertically aligned region of the soil profile, which coincided with strong upward hydraulic gradients recorded in a grid of nested piezometers (0.2‐ and 0.6‐m depth). These patterns suggested a hydraulic connection between the seep and the nitrate‐rich shallow groundwater system below the fragipan, which enabled groundwater and associated nitrate‐N to discharge through the fragipan to the surface. In contrast, time‐lapse ERI indicated no such connections in the non‐seep area, with infiltrated rainwater presumably perched above the fragipan. Results highlight the value of pairing time‐lapse ERI with hydrometric and water quality monitoring to illuminate possible groundwater and nutrient flow pathways to seeps in headwater riparian areas.