Outline
Comptes Rendus

Internal geophysics (Applied geophysics)
Review of self-potential methods in hydrogeophysics
Comptes Rendus. Géoscience, Volume 341 (2009) no. 10-11, pp. 928-936.

Abstracts

The self-potential (SP) method is a passive geophysical method based on the natural occurrence of electrical fields on the Earth's surface. Combined with other geophysical methods, SP surveys are especially useful for localizing and quantifying groundwater flows and pollutant plume spreading, and estimating pertinent hydraulic properties of aquifers (water table, hydraulic conductivity). Laboratory experiments have shown that the involved coupling coefficients mainly depend on the fluid chemistry, conductivity and pH, and on the soil or rock properties. The interpretation of SP observations can be done qualitatively, for instance, by correlation of SP gradients with water fluxes (through electrokinetics) or salt fluxes (through electro-diffusion). In recent years, the interpretation has been improved with the help of modelling or/and inversion of the Poisson equation and endeavours to estimate hydraulic parameters by means of the intensity of electric current sources caused by underground flows.

Les potentiels spontanés (PS) sont les potentiels électriques mesurés de façon passive à la surface de la terre. En association à d’autres méthodes géophysiques, les relevés PS sont particulièrement utiles pour localiser et quantifier les flux d’eau, les panaches de polluant et pour estimer les propriétés hydrauliques des aquifères. Les études de laboratoire ont montré que les coefficients de couplage en jeu dépendent principalement de la composition de l’eau, de sa conductivité, de son pH et de la surface du sol. L’interprétation des observations PS peut se faire qualitativement, en corrélant les gradients PS avec les flux d’eau (électrocinétique) ou avec les flux de sel (électrodiffusion). Récemment, cette interprétation a été améliorée par la modélisation ou/et l’inversion de l’équation de Poisson et vise à estimer les paramètres hydrauliques en terme d’intensité de source de courant électrique induite par les flux d’eau.

Metadata
Received:
Accepted:
Published online:
DOI: 10.1016/j.crte.2009.08.008
Keywords: Self-potential, Electrokinetics, Electrochemical, Water table, Aquifer, Pollution, Inversion
Mot clés : Potentiel spontané, Électrocinétique, Électrochimie, Nappe phréatique, Aquifère, Pollution, Inversion

Laurence Jouniaux 1; Alexis Maineult 2; Véronique Naudet 3; Marc Pessel 4; Pascal Sailhac 1

1 UMR 7516, CNRS, université de Strasbourg (EOST), institut de physique du globe de Strasbourg, 5, rue René-Descartes, 67000 Strasbourg, France
2 UMR 7154, CNRS, université Paris-Diderot, institut de physique du globe de Paris, 4, place Jussieu, case 89, 75252 Paris cedex 05, France
3 EA 4134 (géosciences hydrosciences matériaux constructions), laboratoire Ghymac, université Bordeaux 1, avenue des Facultés, 33405 Talence cedex, France
4 UMR 8148, CNRS, laboratoire IDES, université de Paris-Sud XI, 91405 Orsay, France
@article{CRGEOS_2009__341_10-11_928_0,
     author = {Laurence Jouniaux and Alexis Maineult and V\'eronique Naudet and Marc Pessel and Pascal Sailhac},
     title = {Review of self-potential methods in hydrogeophysics},
     journal = {Comptes Rendus. G\'eoscience},
     pages = {928--936},
     publisher = {Elsevier},
     volume = {341},
     number = {10-11},
     year = {2009},
     doi = {10.1016/j.crte.2009.08.008},
     language = {en},
}
TY  - JOUR
AU  - Laurence Jouniaux
AU  - Alexis Maineult
AU  - Véronique Naudet
AU  - Marc Pessel
AU  - Pascal Sailhac
TI  - Review of self-potential methods in hydrogeophysics
JO  - Comptes Rendus. Géoscience
PY  - 2009
SP  - 928
EP  - 936
VL  - 341
IS  - 10-11
PB  - Elsevier
DO  - 10.1016/j.crte.2009.08.008
LA  - en
ID  - CRGEOS_2009__341_10-11_928_0
ER  - 
%0 Journal Article
%A Laurence Jouniaux
%A Alexis Maineult
%A Véronique Naudet
%A Marc Pessel
%A Pascal Sailhac
%T Review of self-potential methods in hydrogeophysics
%J Comptes Rendus. Géoscience
%D 2009
%P 928-936
%V 341
%N 10-11
%I Elsevier
%R 10.1016/j.crte.2009.08.008
%G en
%F CRGEOS_2009__341_10-11_928_0
Laurence Jouniaux; Alexis Maineult; Véronique Naudet; Marc Pessel; Pascal Sailhac. Review of self-potential methods in hydrogeophysics. Comptes Rendus. Géoscience, Volume 341 (2009) no. 10-11, pp. 928-936. doi : 10.1016/j.crte.2009.08.008. https://comptes-rendus.academie-sciences.fr/geoscience/articles/10.1016/j.crte.2009.08.008/

Version originale du texte intégral

1 Introduction

Measurements of self-potential (SP) are performed with non-polarizable electrodes, in various geophysical contexts and at different scales, from centimetric to kilometric, at the Earth surface or in boreholes [30,31]. The SP method has been used for a variety of geophysical applications [51]: to characterize active volcanic areas [3,14], to detect and monitor groundwater flow [13,45], to image karstic conduits [62] or contaminant plumes [36] and to study landslides [39,46]. The charge polarisation mechanisms proposed to explain SP anomalies are associated with electrokinetic, electrochemical, thermoelectric, redox, and piezoelectric effects. The general equation for coupled flows can be written [43]

Ji=jLijXj(1)
where the fluxes Ji (of charges, matter, heat, etc) are related to the various forces Xj (gradients of electrical potential, pressure, temperature, etc) through the coupling coefficients Lij (“phenomenological coefficients” [12] or “conductivities” [59]). In this article, we detail the electrokinetic and electrochemical origins of the SP signals, followed by a review of SP inversion techniques.

2 Self-potentials with electrokinetic origin

In steady-state, and under isothermal conditions, the electric flow (Je, in ampere per meter cube) is related to the pressure gradient (grad P) and to the electric potential gradient (grad V) by the following coupling equation:

Je=L22gradV+L21gradP.(2)

The first term on the right hand side in Eq. (2) represents Ohm's law. The coefficient L22 is the rock conductivity σr. Although the coupling coefficient L21 could be measured, the electrokinetic coefficient Cs is more often measured in laboratory, by applying a fluid flow (ΔP) and by measuring the induced electric potential (ΔV) (Fig. 1) [21,47].

Cs=ΔVΔP=εζησf,(3)
where σf and ɛ are the fluid conductivity and the dielectric constant of the fluid, ζ is the zeta potential (within the double-layer at the interface between the rock and the fluid, see [11]), and η the dynamic viscosity of the fluid. Eq. (2) can be written:
Je=σrgradVCsσrgradP.(4)

Fig. 1

Example of measured streaming potential ΔV when the applied driving pore pressure is 0.1 MPa in one way (positive signal) and in the reverse way (negative signal) within a sandstone sample. Water resistivity is 720 Ω m, and permeability is 6.5 × 10−14 m2.

Exemple de mesure de potentiel d’électrofiltration ΔV, quand on applique une différence de pression de 0,1 MPa dans un sens (signal positif) ou en sens inverse (signal négatif) sur un échantillon de grès de Fontainebleau. La résistivité de l’eau est de 720 Ω m et la perméabilité du grès de 6,5 × 10−14 m2.

Due to this electrokinetic origin, the corresponding SP signal is sometimes called streaming potential. The electrokinetic coefficient is also directly involved in seismoelectrics effects [6,18]. Note that the electrokinetic coefficient is not expected to depend on permeability or porosity, although this can be observed in some cases. Indeed, the electrokinetic coefficient can be affected through an additional surface conductivity (which is neglected in Eq. (3)) when water conductivity is very low. Based on laboratory studies, electrokinetic coefficient values are roughly ∼ −10−6 to −10−7 V/Pa (or −10−2 to −10−3 V/m) for water conductivity values from 0,01 to 0,1 S/m (and pH around 7), for quartz sands, sandstones, granites, volcanic ashes. The electrokinetic coefficient depends on the pH value: the isoelectric point (IEP) is the pH value at which the electrokinetic coefficient is zero. The IEP value is about pH 2∼3 for a quartz-water interface, and about pH 8∼11 for a calcite-water interface [22]. Below the IEP, the electrokinetic coefficient is positive and, above the IEP, it is negative. For most of the earth sciences applications, values of the electrokinetic coefficients are negative (meaning that positive ions are carried in the flow direction), because involved pH values are usually above the IEP. However, with calcite or volcanic rocks, the electrokinetic coefficient can be positive [23]. The presence of clay can also modify the electrokinetic coefficient, depending on the type of clay and its amount. The zeta potential value (Eq. (3)) is not directly measured, but is deduced from streaming potential or electro-osmotic measurements, using some hypotheses. Therefore, the zeta potential values available in the literature for rock samples have to be taken with caution, and it is recommended to use electrokinetic coefficient values rather than zeta potential values for electrokinetic modelling.

Based on theoretical considerations, it has been proposed [8,54] that the electrokinetic coefficient is inversely proportional to the effective saturation. However, the few observations published up to now [21] show that the electrokinetic coefficient is either non-dependent or proportional to the effective saturation as:

Cs=ΔVΔP=εζησfSen(5)

Based on this electrokinetic effect, SP observations can yield an estimate of aquifer hydraulic properties [9,53]. The time variations of SP can be identified and associated with time-varying fluid flow from metric to kilometric scale [13,45]. The positive anomalies observed in active areas are used to define hydrothermal zones [14,23]. Moreover, SP are monitored during hydraulic tests in boreholes [10,30,31,52] and can provide an estimate of the fracture aperture [24] or permeability [34].

3 Self-potentials with electrochemical origin

Electrochemical sources can be attributed to several phenomena. The common one is the diffusion of ions due to a concentration gradient between two regions. This source current is balanced by a current flowing throughout the Earth conductivity structure, so that the total current density (Eq. (1)) is divergence-free. Another electrochemical source mechanism is associated with redox processes occurring in ore bodies and contaminant plumes. This source current must also conform to the requirement that the total electric current density is divergence-free. In this section, we present the theoretical background, laboratory and field SP investigations related to these two electrochemical sources.

3.1 Junction potential

The junction potential results from the separation of the ionic charges occurring across a concentration gradient provided that the anions and cations have different ionic mobility. In the case of a non-moving fluid, the generalized ionic diffusion equation for the ith ionic species, i = 1 N, can be written:

Cit=divDigradCi+siuidivCigradV,(6)
where Ci is the concentration, ui the ionic mobility, si the valence times the sign of the charge, Di = RTui/(Ae) is the so-called ionic diffusion coefficient, T the absolute temperature, e the absolute elementary charge, R the molar gas constant and A the Avogadro constant. The Eq. (6), combined with the electroneutrality condition, provides a non-linear system governing the evolution of the Ci, i = 1 N, and V. In the simple case of a monovalent, symmetric and binary salt such as Na+Cl, the re-arrangement of these equations leads to the Nernst-Einstein relation for Fickian diffusion:
Ct=divDmgradC=div2RTAeuCluNauCl+uNagradC,(7)
where Dm is the molecular diffusion coefficient, which then leads to the Planck-Henderson equation:
gradV=αmgradCC=RTAeuCluNauCl+uNagradCC,(8)
where αm is the fluid junction coupling coefficient in fluid. In a porous medium where surface conduction can be neglected, the effective diffusion coefficient D is usually taken equal to Dm divided by the tortuosity τ and the effective coupling coefficient α to αm multiplied by the porosityφ. Note that Eq. (8) is a local equation to be integrated over the considered domain in order to find the macroscopic potential difference.

The SP response of diffusion through clayey materials is quite different due to their high intrinsic charge, large specific area, high cationic exchange capacity and small pore throats. To briefly describe the process: since the counter-ions are attracted to the mineral surface and the pore throats are small, the co-ions are electrostatically blocked at the entrance of the throats. This triggers a charge separation and thus a so-called “membrane potential”, and decreases the diffusion process. The expressions of α and D must then be modified in consequence [50].

The junction potential is generally neglected when the surface observations are interpreted [7]. Indeed, since α is around a few millivolts (5.24 for NaCl), the expected SP amplitude in non-clayey materials is very low (i.e., < 20 mV) for “usual” concentration gradients. Within sand, maximum SP differences of around 13 mV were observed for an initial NaCl concentration step of 2250/1 [26]. In contrast, log analysts interpret the borehole SP signals as a combination of membrane (associated with clay layers) and junction (associated with sand or sandstones layer) potentials, and neglect the other SP sources.

When the transport process is not restricted to pure diffusion but also includes a movement of the fluid, Eq. (7) has to be replaced by the advection–dispersion equation and D by the hydrodynamic dispersion coefficient. For sand, Maineult et al. [27] reported that, in this case, the SP variations with respect to the initial state (i.e., prior to the establishment of the gradient concentration) are the sum of two contributions: (a) the junction potential variation (Eq. (8)) and (b) the variation of the electrokinetic response, resulting from the local change in the salinity σf and in the ζ–potential that modifies the electrokinetic coefficient Cs (Eq. (5)). Sandberg et al. [55] measured surface signals up to 100 mV over a NaCl injection 3–6 m deep in unconsolidated sediments with an initial concentration ratio of 150/1, and proposed a good estimate of the fluid velocity, but without modelling. They, therefore, suggested that SP could be an alternative method to the salt fronts at the field scale. Recently, it has been shown that SP measurements could be used to estimate not only the fluid velocity but also the dispersion coefficient of granular media [27]. Moreover, SP measurements can be used to reconstruct mixing and reaction zones between fluids [29], or to study the effect of chemical reactions between the fluid and a reactive matrix [28].

3.2 “Electro-redox” potential

When the chemical potential is a redox potential, the transfer of electrons through an electronic conductor can generate large SP anomalies in the surrounding conductive medium. Strong negative SP anomalies with magnitudes usually reaching a few hundred millivolts have been reported over ore deposits for more than 50 years, and explained thanks to a geobattery model [4]. In this model, the ore body connecting two regions of different redox potential participates directly in oxidizing and reducing reactions occurring at the bottom and top of the ore body, respectively. Then, the ore body serves as an electronic conductor to transfer electrons from the anode (bottom part) to the cathode (upper part). In contaminated environments, redox processes also occur since the contaminant is biodegraded through microbially mediated redox reactions. Despite the long history of the SP method, its application to the study of contaminant plumes is very recent [33,36,37,40,42].

In the South of France, the Entressen landfill, which is responsible for an organic-rich contaminant plume spreading in a shallow unconfined aquifer, has been studied through SP measurements. A strong negative SP anomaly (several hundreds of millivolts) has been identified downstream of the landfill [36]. After removing the electrokinetic component associated with the groundwater flow, the residual SP anomaly, called the “electro-redox” potential by the authors, has been shown to be directly proportional to the redox potential [37]. More recently, new SP and electrical conductivity measurements have been performed in the first kilometre downstream of this landfill (Fig. 2). The electrical conductivity map was drawn up using a Geonics EM31 in vertical magnetic dipole mode, an electromagnetic two-coil system with an investigation depth of about 5.5 m; it highlights the spreading of the contaminant plume and is in good agreement with water conductivity values measured in boreholes. As for previous measurements, the SP map clearly shows negative anomalies with respect to the reference station located at the eastern border of the landfill. Based on these two geophysical maps, the contaminant plume can be located.

Fig. 2

Observations performed downstream the Entressen landfill: (a): electrical conductivity map (in millisiemens per metre) and water conductivity measured in boreholes (black points); (b): SP map (in millivolts) with the redox potential (Eh) values measured in boreholes (black points). Water conductivity and redox potential measured upstream the landfill are respectively 70 mS/m and 150 mV.

Mesures obtenues en aval du centre d’enfouissement technique d’Entressen : (a) : carte de conductivité électrique (en millisiemens par mètre), et valeurs de conductivité électrique de l’eau mesurée en puits (points noirs) ; (b) : carte PS (en millivolts) avec les valeurs de potentiel redox (Eh) mesuré en puits (points noirs). La conductivité de l’eau et le potentiel redox mesurés en amont du terrain sont 70 mS/m et 150 mV, respectivement.

To further evaluate the potential role of bacteria in the generation of an electrical current, sandbox experiments were made with sulphate-reducing bacteria and organic nutrients [35]. The results showed a linear relationship between SP signals and redox potential variations, apparently related to bacterial activity. Therefore, the source current density Jer associated with this so-called “electro-redox” effect can be expressed as follows [2]:

Jer=σEH,(9)

where σ is the electrical conductivity of the volume characterized by a large variation in the redox potential EH between the contaminated aquifer and the vadose zone. Based on these results and on the SP theory developed for massive ore deposits [56], a biogeobattery model is proposed to explain this “electro-redox” effect observed on contaminant plumes [2,35,41]. In this model (Fig. 3), electron donating (organic matter oxidation) and accepting (O2, sulphate, nitrate, ferric iron reduction) processes can be respectively associated with the anode and cathode of the geobattery. Electron transfer through this system is proposed to be driven by two possible mechanisms: (i) through extracellular appendages called “nanowires” and/or (ii) through the biominerals resulting from precipitation of metallic particles induced by bacterial activity. Indeed, under electron acceptor limitation, in order to reach other electron acceptors, specific bacteria have been shown to produce structured nanowires known to serve as electron transfer pathways and, thus, to produce an electrical field [41,48]. These two processes can therefore serve as electrical connections that link the oxidizing and reducing zones, thus permitting electron transfer, and producing an electrical field. Such a biogeobattery was numerically modelled on the Entressen landfill [2]. Field data as well as numerical modelling indicate that the redox contribution can be explained by dipoles distributed throughout the water table with a strength proportional to the difference in the redox potentials between the aquifer and the vadose zone. This model was successfully applied to invert SP signals measured over the contaminant plume of the Entressen landfill [36] in terms of redox potential distribution [25].

Fig. 3

Conceptual biogeobattery model describing the “electro-redox” effect that generates electrical field driven by microbial activity and therefore induces SP anomalies. Oriented microbial nanowires and/or biominerals are conductors for electron produced in the reduced zone (anode) and consumed in the oxidized zone (cathode). An electric current is then flowing from the anode to the cathode of this system.

Modèle conceptuel de biogéobatterie schématisant l’effet « électro-redox » qui génère un champ électrique lié à l’activité microbienne et qui, par conséquent, induit des anomalies PS. Les filaments protéiques extracellulaires et/ou les biominéraux orientés sont des conducteurs pour les électrons produits dans la région réductrice (anode) et captés dans la région oxydée (cathode). Un courant électrique circule alors de l’anode vers la cathode de ce système.

All these results have encouraged research in the emerging discipline called biogeophysics, which refers to the application of geophysical methods to observe microbial processes in the subsurface. From recent field and laboratory studies, bacteria have been shown to play an important role in the modification of geoelectrical parameters, in particular, through mineral weathering, growth of microbial cells, biomineralization, and microbial nanowires [1,41].

4 Self-potentials inversion

Let us consider the interpretation of SP data of electrokinetic origin. The fundamental equation connecting the electric streaming potential V to the hydraulic pressure P follows the conservation of the total electric flux Je defined in Eq. (2); in steady-state and without direct electric current source, it is:

Je=σrgradVCsσrgradP=0.(10)

This is a Poisson equation for the electric potential V where the divergence of the hydraulic potential gradient is a primary electrokinetic source term EKS [15,59]:

σrgradV=EKS,(11a)
where
EKS=CsσrgradP.(11b)

In specific cases, such as electrokinetics within rock samples used in the laboratory, a constant hydraulic pressure gradient is applied between the upper and lower boundaries of a vertical homogeneous porous cylinder. In these conditions, the water is uniformly distributed within the sample. Therefore, the water diffusion is in steady-state regime and there is no internal source. Thus, the gradient of P is constant within the sample and the conservation of the total electric flux Je involves the proportionality of the SP gradients to the Darcy velocity. It is possible in such a case to simply obtain proportionality between potential differences measured at upper and lower boundaries of the sample and the Darcy velocity. This proportionality also occurs in the field in the case of a 1D infiltration experiment [61].

However, in the general underground case with heterogeneity and preferential flow paths, there are differences between parameters of the hydraulic and electric modelling with different sources, conductivity and boundary conditions, so there is a non unique conversion of SP gradients into Darcy velocities. This conversion is an ill-posed inverse problem that needs appropriate methods. Let us recall the theory. Fitterman [15] and later Fournier [16] were the first to consider solutions to the Poisson equations ((10) or (11a), (11b)) using the total electric potential ψ = V+CsP and potential integrals in case of heterogeneous media with contrasts in the electrokinetic coefficient Cs or the electric conductivity σr. Especially, Claude Fournier considered the first integral formulae for the streaming potential caused by aquifers in tabular media. The SP solution for an unconfined aquifer in a medium of constant electric conductivity provides a classical convolution integral [17]:

V(r)=CSρgSnG(rr)H(r)dr2,(12)
where G(r−r)=1/(r–r) is the Green function for the homogeneous half-space (r at the ground surface and r′ on the piezometric surface S), ∂/∂n is the partial derivative in the normal direction to S, H is the piezometric head, that can be related to pressure (=P/ρg). Fournier's equation shows that SP data near the ground of a homogeneous electrical medium are similar to the potential field produced by dipoles distributed throughout the water table with a strength that is proportional to the piezometric depth. Similar equations have been introduced to propose images of SP sources in the underground via correlation integrals [44,49]. Other theories, also based on classical convolution integrals of potential field theory, have been developed to inverse SP data in the wavelet domain, in order to identify location, intensity and type of singularities of causative underground hydraulic flows [19,53]. Gibert and Sailhac [20] have recently commented Patella's correlation approach to demonstrate that the so-called probability of tomography defines images of SP data in the wavelet domain that must not be considered as underground images of SP sources; appropriate inversion is necessary to achieve underground images. By inversion, we mean the trial-and-error procedure that minimizes differences between observed SP data and modelled SP data obtained for specific model parameters [60]. For instance, SP data from the Vulcarolo fissure in the Etna volcano were interpreted by inversion of the wavelet transform SP data, to estimate the water flux within the Vulcarolo fissure [53]. More recently, one finds several applications of actual inversion of the electric Poisson application (Eq. (11a)), with tomography of SP sources intensity [32,58], allowing the inversion code to consider a heterogeneous electric conductivity distribution, as observed by means of electrical resistivity tomography (ERT). Another method consists in using the Particle swarm optimisation (PSO), which is a global optimisation strategy that simulates the social behaviour observed in a flock (swarm) of birds searching for food [57]. Recently, Naudet et al. [38] used the Particle Swarm algorithm to estimate the water table elevation from SP.

In hydrogeology and soil physics, images of electrokinetic sources (EKS in Eq. (11a)) can be useful as a first step, but it is probably not the best way forward. Indeed, the actual water flux and hydraulic parameters are related to the pressure field obtained by solving equation (11b) using electrokinetic sources of Eq. (11a). Currently, there are just a very few examples of applications to SP data involving tomography of actual hydrology parameters. A typical example of an interesting hydrogeological application has been considered [9] to determine hydrogeological parameters from SP data recorded at a producing well in an unconfined aquifer (using classical data from [5]). Darnet et al. [9] used a genetic algorithm for the inversion and instead of Fournier's equation of the SP solution for an unconfined aquifer in a medium of constant electric conductivity, they used the Dupuit hypothesis and considered a 3D integral formula:

V(r,z=0)=CS2πQπKS2Vh02+QπKSlnrr03/2rrcosθr2+rsinθ2+z2dθdrdz(13)
where the integration in cylindrical coordinates includes the whole volume between the top of the unconfined aquifer and the ground surface, Q is the hydraulic flux at the well head, KS is the hydraulic conductivity of the saturated rock and h0 is the piezometric head at a given radial distance r0 (related to the Dupuit hypothesis). A clear sensitivity of SP data to Cs and to the ratio Q/KS has been demonstrated [9]. In the vadose zone a clear sensitivity of SP data to Cs and to the ratio Q/KS was also shown by numerical simulations of 2D infiltration, as well as the sensitivity to the specific length α and to the effective saturation at the beginning of the infiltration experiment [54].

5 Conclusions

There is an increasing interest in the use of the SP method for non-invasive characterizing of subsurface groundwater flow and contamination. The interpretation of SP signals is underlined by the understanding of the main physical process involved in their origin. When the electrokinetic origin is involved, a better quantification of the electrokinetic coefficient as a function of the water content is needed. We also hope that, ultimately, the SP method evolves into a non-invasive tool for characterization and long-term monitoring of microbial processes. Regarding inversion, moving tomography of the electrokinetic source term into tomography of some hydraulic parameters such as the ratio of the fluid flow to the permeability is the real challenge of the next decade.

Acknowledgements

This research was supported by the CNRS. The CNRS-INSU is thanked for its support regarding the EC2CO-BGPhy program and ECCO-PNRH-Waterscan program. Laurence Jouniaux and Alexis Maineur do not acknowledge their supervisory ministry for its action against the French public system for research and education. This is IPGP contribution no 2466.


References

[1] E. Atekwana; D.D. Werkema; E.A. Atekwana Biogeophysics: the effects of microbial processes on geophysical properties of the shallow subsurface, Hydrogeophysics, Volume 71 (2006), pp. 161-193

[2] T. Arora; N. Linde; A. Revil; J. Castermant Non-intrusive characterization of the redox potential of landfill leachate plumes from self-potential data, J. Contam. Hydrol. (2007), pp. 274-292

[3] M. Aubert; Q.Y. Atangana Self-potential method in hydrogeological exploration of volcanic areas, Groundwater, Volume 34 (1996), pp. 1010-1016

[4] J. Bigalke; E.W. Grabner The geobattery model: a contribution to large-scale electrochemistry, Electrochim. Acta, Volume 42 (1997), pp. 3443-3452

[5] V.A. Bogoslovsky; A.A. Ogilvy Deformations of natural electric fields near drainage structures, Geophys. Prospect., Volume 21 (1973), pp. 716-723

[6] C. Bordes; L. Jouniaux; M. Dietrich; J.-P. Pozzi; S. Garambois First laboratory measurements of seismo-magnetic conversions in fluid-filled Fontainebleau sand, Geophys. Res. Lett., Volume 33 (2006), p. L01302

[7] R.F. Corwin; D.B. Hoover The self-potential method in geothermal exploration, Geophysics, Volume 44 (1979), pp. 226-245

[8] M. Darnet; G. Marquis Modelling streaming potential (SP) signals induced by water movement in the vadose zone, J. Hydrol., Volume 285 (2004), pp. 114-124

[9] M. Darnet; G. Marquis; P. Sailhac Estimating aquifer hydraulic properties from the inversion of surface Streaming Potential (SP) anomalies, Geophys. Res. Lett., Volume 30 (2003), pp. 12-21

[10] M. Darnet; G. Marquis; P. Sailhac Hydraulic stimulation of geothermal reservoirs: fluid flow, electric potential and microseismicity relationships, Geophys. J. Int., Volume 166 (2006), pp. 438-444

[11] J.A. Davis; R.O. James; J. Leckie Surface ionization and complexation at the Oxide/water interface, J. Colloid Interface Sci., Volume 63 (1978), pp. 480-499

[12] S.R. de Groot, P. Mazur, Non-linear thermodynamics. Dover Publications. 1983. 510 p.

[13] C. Doussan; L. Jouniaux; J.-L. Thony Temporal variations of SP and unsaturated water flow in loam and clay soils: a seasonal field study, J. Hydrol., Volume 267 (2002), pp. 173-185

[14] A. Finizola; J.-F. Lénat; O. Macedo; D. Ramos; J.-C. Thouret; F. Sortino Fluid circulation and structural discontinuities inside Misti volcano (Peru) inferred from self-potential measurements, J. Volc. Geoth. Res., Volume 135 (2004), pp. 343-360

[15] D.V. Fitterman Electrokinetic and magnetic anomalies associated with dilatant regions in a layered Earth, J. Geophys. Res., Volume 83 (1978), pp. 5923-5928

[16] C. Fournier, Méthodes géoélectriques appliquées à l’hydrogéologie en région volcanique (Chaîne-des-Puys, Massif central français), Développement de la méthode des potentiels spontanés en hydrogéologie, Thèse de 3e cycle (1983) 157 p.

[17] C. Fournier 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), Geophys. Prospect., Volume 37 (1988), pp. 647-668

[18] S. Garambois; M. Ditrich Seismoelectric wave conversions in porous media: field measurements and transfer function analysis, Geophysics, Volume 66 (2001), pp. 1417-1430

[19] D. Gibert; M. Pessel Identification of sources of potential fields with the continuous wavelet transform: application to SP profiles, Geophys. Res. Lett., Volume 28 (2001), pp. 1863-1866

[20] D. Gibert; P. Sailhac Comment on the paper by A. Jardani, J.P. Dupont, and A. Revil (2006), Self-potential signals associated with preferential groundwater flow pathways in sinkholes, J. Geophys. Res., Volume 111 (2008), p. B09204

[21] X. Guichet; L. Jouniaux; J.-P. Pozzi Streaming potential of a sand column in partial saturations conditions, J. Geophys. Res., Volume 108 (2003), p. 2141

[22] X. Guichet; L. Jouniaux; N. Catel Modification of streaming potential by precipitation of calcite in a sand-water system: laboratory measurements in the pH range from 4 to 12, Geophys. J. Int., Volume 166 (2006), pp. 445-460

[23] H. Hase; T. Hashimoto; S. Sakanaka; W. Kanda; Y. Tanaka Hydrothermal system beneath Aso volcano as inferred from self-potential mapping and resistivity structure, J. Volc. Geoth. Res., Volume 143 (2005), pp. 259-277

[24] C.W. Hunt; M.H. Worthington Borehole electrokinetic responses in fracture dominated hydraulically conductive zones, Geophys. Res. Lett., Volume 27 (2000), pp. 1315-1318

[25] N. Linde; A. Revil Inverting residual self-potential data for redox potentials of contaminant plumes, Geophys. Res. Lett., Volume 34 (2007), p. L14302

[26] A. Maineult; Y. Bernabé; P. Ackerer Electrical response of flow, diffusion and advection in a laboratory sand-box, Vadose Zone J., Volume 3 (2004), pp. 1180-1192

[27] A. Maineult; Y. Bernabé; P. Ackerer Detection of advected concentration and pH fronts from self-potential measurements, J. Geophys. Res., Volume 110 (2005), p. B11205

[28] A. Maineult; L. Jouniaux; Y. Bernabé Influence of the mineralogical composition on the self-potential response to advection of KCl concentration fronts through sand, Geophys. Res. Lett., Volume 33 (2006), p. L24311

[29] A. Maineult; Y. Bernabé; P. Ackerer Detection of advected, reacting redox fronts from self-potential measurements, J. Contam. Hydrol., Volume 86 (2006), pp. 32-52

[30] A. Maineult; E. Strobach; J. Renner Self-potential signals induced by periodic pumping tests, J. Geophys. Res., Volume 113 (2008), p. B01203

[31] G. Marquis; M. Darnet; P. Sailhac; A.K. Singh; A. Gérard Surface electric variations induced by deep hydraulic stimulation: an example from the Soultz HDR site, Geophys. Res. Lett., Volume 29 (2002) | DOI

[32] B.J. Minsley, Modeling and inversion of self-potential data, PhD Thesis, Massachusetts Institute of Technology (2007) 251 p.

[33] B.J. Minsley; J. Sogade; F.D. Morgan Three-dimensional self-potential inversion for subsurface DNAPL contaminant detection at the Savannah River Site, South Carolina, Water Resour. Res., Volume 43 (2007), p. W04429

[34] H. Murakami; T. Hashimoto; N. Oshiman; S. Yamaguchi; Y. Honkura; N. Sumitomo Electrokinetic phenomena associated with a water injection experiment at the Nojima fault on Awaji Island, Japan, Isl. Arc, Volume 10 (2001), pp. 244-251

[35] V. Naudet; A. Revil A sandbox experiment to investigate bacteria-mediated redox processes on self-potential signals, Geophys. Res. Lett., Volume 32 (2005), p. L11405

[36] V. Naudet; A. Revil; J.-Y. Bottero; P. Bégassat Relationship between self-potential (SP) signals and redox conditions in contaminated groundwater, Geophys. Res. Lett., Volume 30 (2003)

[37] V. Naudet; A. Revil; E. Rizzo; J.-Y. Bottero; P. Bégassat Groundwater redox conditions and conductivity in a contaminant plume from geoelectrical investigations, Hydrol. Earth Syst. Sci., Volume 8 (2004), pp. 8-22

[38] V. Naudet; J.L. Fernández-Martínez; E. García-Gonzalo; J.P. Fernández-Álvarez Estimation of water table from self-potential data using Particle Swarm Optimization (PSO), SEG, Volume 27 (2008), pp. 1203-1207 ([expanded abstract])

[39] V. Naudet; M. Lazzari; A. Perrone; A. Loperte; S. Piscitelli; V. Lapenna Integrated geophysical and geomorphological approach to investigate the snowmelt-triggered landslide of Bosco Piccolo village (Basilicata, southern Italy), Eng. Geol., Volume 98 (2008), pp. 156-167

[40] R.E. Nimmer; J.L. Osiensky Direct current and self-potential monitoring of an evolving plume in partially saturated fractured rock, J. Hydrol., Volume 267 (2002), pp. 258-272

[41] D. Ntarlagiannis; E.A. Atekwana; E.A. Hill; Y. Gorby Microbial nanowires: is the subsurface “hardwired”?, Geophys. Res. Lett., Volume 34 (2007), p. L17305

[42] J.E. Nyquist; C.E. Corry Self-potential: the ugly duckling of environmental geophysics Lead, Edge, Volume 21 (2002), pp. 446-451

[43] J. Overbeek Electrochemistry of the double layer (H.R. Kruyt, ed.), Colloid Science, 1, Irreversible Systems, Elsevier, New York, 1952 (pp. 115)

[44] D. Patella Introduction to ground surface self-potential tomography, Geophys. Prospect., Volume 45 (1997), pp. 653-681

[45] F. Perrier; M. Trique; B. Lorne; J.-P. Avouac; S. Hautot; P. Tarits Electric potential variations associated with lake level variations, Geophys. Res. Lett., Volume 25 (1998), pp. 1955-1958

[46] A. Perrone; A. Iannuzzi; V. Lapenna; P. Lorenzo; S. Piscitelli; E. Rizzo; F. Sdao High-resolution electrical imaging of the Varco d’Izzo earthflow (southern Italy), J. Appl. Geophys., Volume 56 (2004), pp. 17-29

[47] J.-P. Pozzi; L. Jouniaux Electrical effects of fluid circulation in sediments and seismic prediction, C. R. Acad. Sci. Paris, Ser. II, Volume 318 (1994), pp. 73-77

[48] G. Reguera; K.P. Nevin; J.S. Nicoll; S.F. Cocalla; T.L. Woodard; D.R. Lovley Biofilm and nanowire production leads to increased current in geobacter sulfurreducens fuel cells, Appl. Environ. Microbiol., Volume 72 (2006), pp. 7345-7348

[49] A. Revil; V. Naudet; J. Nouzaret; M. Pessel Principle of electrography applied to self-potential electrokinetic sources and hydrogeological applications, Water Resour. Res., Volume 39 (2003), p. 1114

[50] A. Revil; P. Leroy; K. Titov Characterization of transport properties of argillaceous sediments: application to the Callovo-Oxfordian argillite, J. Geophys. Res., Volume 110 (2005), p. B06202

[51] A. Revil; K. Titov; C. Doussan; V. Lapenna Applications of the self-potential method to hydrological problems In: Applied hydrogeophysics, Springer Netherlands, NATO Ser. IV, Volume 71 (2006), pp. 255-292

[52] E. Rizzo; B. Suski; A. Revil; S. Straface; S. Troisi Self-potential signals associated with pumping tests experiments, J. Geophys. Res., Volume 109 (2004), p. B10203

[53] P. Sailhac; G. Marquis Forward and inverse modelling of SP anomalies caused by subsurface flow, Geophys. Res. Lett., Volume 28 (2001), pp. 1851-1854

[54] P. Sailhac; M. Darnet; G. Marquis Electrical streaming potential measured at the ground surface: forward modeling and inversion issues for monitoring infiltration and characterizing the vadose zone, Vadose Zone J., Volume 3 (2004), pp. 1200-1206

[55] S.K. Sandberg; L.D. Slater; R. Versteeg An integrated geophysical investigation of the hydrology of an anisotropic unconfined aquifer, J. Hydrol., Volume 267 (2003), pp. 227-243

[56] M. Sato; H.M. Mooney The electrochemical mechanism of sulfide self-potentials, Geophysics, Volume 25 (1960), pp. 226-249

[57] R. Shaw; S. Srivastava Particle swarm optimization: a new tool to invert geophysical data, Geophysics, Volume 72 (2007), p. F75-F83

[58] M.R. Sheffer, Forward modelling and inversion of streaming potential for the interpretation of hydraulic conditions from self-potential data, PhD Thesis, Univ. of British Columbia (2007) 207 p.

[59] W.R. Sill Self-potential modeling from primary flows, Geophysics, Volume 48 (1983), pp. 76-86

[60] A. Tarantola; B. Valette Inverse problems = quest for information, J. Geophys., Volume 50 (1982), pp. 159-170

[61] J.-L. Thony; P. Morat; G. Vachaud; J.-L. Le Mouël Field characterization of the relationship between electrical potential gradients and soil water flux, C. R. Acad. Sci. Paris, Ser. IIa, Volume 325 (1997), pp. 317-321

[62] Z. Wanfang; B.F. Beck; J.B. Spephenson Investigation of groundwater flow in karst area using component separation of natural potential measurements, Environ. Geol., Volume 37 (1999), pp. 19-25


Comments - Policy