Reply to ‘Comment on “Dependence of shear wave seismoelectrics on soil textures: a numerical study in the vadose zone by F.I. Zyserman, L.B. Monachesi and L. Jouniaux” by Revil, A.’

In this paper we reply to a the comment made by Revil (2017) on our paper (2017, Geophys. J. Int., 208), where we describe seismoelectric phenomena in the vadose zone based on the theory of Pride empirically extended for unsaturated conditions. We analyse and answer each one of the enumerated critics, and reaffirm the conclusions of our work. In particular, we prove that using the conductivity model suggested by Revil (2017) does not change our predictions significantly, contrary to what was argued in the comment. Further, in the light of previous and new theoretical and experimental results existing in the literature, we confirm the reasonability of having tested a non-monotonic saturation dependent streaming potential coefficient model besides the monotonic one, and discuss the suitability of assuming a linear relation between the permeability and the excess charge.Fil: Zyserman, Fabio Ivan. Consejo Nacional de Investigaciones Científicas y Técnicas; Argentina. Universidad Nacional de la Plata. Facultad de Ciencias Astronómicas y Geofísicas. Departamento de Geofísica Aplicada; ArgentinaFil: Monachesi, Leonardo Bruno. Consejo Nacional de Investigaciones Científicas y Técnicas; Argentina. Universidad Nacional de Río Negro. Sede Alto Valle. Instituto de Investigaciones en Paleobiología y Geología; ArgentinaFil: Jouniaux, L.. Centre National de la Recherche Scientifique; Franci

Seismo-electrics, electro-seismics, and seismo-magnetics for earth sciences

International audienceThe seismo-electromagnetic method (SEM) is used for non-invasive subsurface exploration. It shows interesting results for detecting fluids such as water, ice, oil, gas, CO_2 , and also to better characterise the subsurface in terms of porosity, permeability, and fractures. However, a limitation of this method is the low level of the induced 5 signals. We first describe SEM's theoretical background, and the role of some key parameters. We then detail recent studies on SEM, through theoretical and numerical developments, and through field and laboratory observations, to show that this method can bring advantages compared to classical geophysical methods

A review on electrokinetically induced seismo-electrics, electro-seismics, and seismo-magnetics for Earth sciences

The seismo-electromagnetic method (SEM) can be used for non-invasive subsurface exploration. It shows interesting results for detecting fluids such as water, oil, gas, CO2, or ice, and also help to better characterise the subsurface in terms of porosity, permeability, and fractures. However, the challenge of this method is the low level of the induced signals. We first describe SEM's theoretical background, and the role of some key parameters. We then detail recent studies on SEM, through theoretical and numerical developments, and through field and laboratory observations, to show that this method can bring advantages compared to classical geophysical methods.Facultad de Ciencias Astronómicas y Geofísica

SH Seismoelectric Response of a Glacier : An Analytic Study

In this work we derive the analytic solutions to the system of equations modeling, within the framework of Pride's theory, the seismic-to-electromagnetic conversions taking place in a glacial environment. Considering a one-dimensional approach, we set a pure shear horizontal wave seismic source on top of an elastic medium representing the glacier, which overlies a porous medium fully saturated with water, representing the glacier bed. The obtained solutions allow to separately represent and analyze the induced electromagnetic responses, the so called coseismic waves, for both the electric and magnetic fields along with the signals originated at the glacier bottom, the electric interface response, and magnetic interface response. We also propose approximate solutions, useful to be used in a fast inversion algorithm. We analyze the characteristics of the induced electromagnetic signals and their dependence on the type of glacier bed, considering an unconsolidated one and a consolidated one. The main results of the present paper are manifold, on the one hand, the mentioned analytic solutions, on the other hand, that the electric interface response originated at the glacier bottom is proportional to the electric current density at this depth, and depends on textural and electrical properties of the basement. We also showed that the amplitude of the electric interface response is three orders of magnitude higher than the amplitude of the electric coseismic field. This fact reinforces the idea proposed in our previous works that it would be interesting to test shear horizontal seismoelectrics as a possible geophysical prospecting and monitoring tool.Facultad de Ciencias Astronómicas y Geofísica

An analytical solution to assess the SH seismoelectric response of the vadose zone

We derive an analytical solution of the seismoelectric conversions generated in the vadose zone, when this region is crossed by a pure shear horizontal (SH) wave. Seismoelectric conversions are induced by electrokinetic effects linked to relative motions between fluid and porous media. The considered model assumes a 1D soil constituted by a single layer on top of a half-space in contact at the water table, and a shearing force located at the earth?s surface as the wave source. The water table is an interface expected to induce a seismoelectric interfacial response (IR). The top layer represents a porous rock in which porous space is partially saturated by water and air, while the half-space is completely saturated with water, representing the saturated zone. The analytical expressions for the coseismic fields and the interface responses, both electric and magnetic, are derived by solving Pride's equations with proper boundary conditions. An approximate analytical expression of the solution is also obtained, which is very simple and applicable in a fairly broad set of situations. Hypothetical scenarios are proposed to study and analyse the dependence of the electromagnetic fields on various parameters of the medium. An analysis of the approximate solution is also made together with a comparison to the exact solution. The main result of the present analysis is that the amplitude of the interface response generated at the water table is found to be proportional to the jump in the electric current density, which in turn depends on the saturation contrast, poro-mechanical and electrical properties of the medium and on the amplitude of the solid displacement produced by the source. This result is in agreement with the one numerically obtained by the authors, which has been published in a recent work. We also predict the existence of an interface response located at the surface, and that the electric interface response is several orders of magnitude bigger than the electric coseismic field, whereas it is the opposite using compressional waves as shown by theoretical and experimental results. This fact should encourage the performance of field and laboratory tests to check the viability of SHTE seismoelectrics as a near surface prospecting/monitoring tool.Facultad de Ciencias Astronómicas y Geofísica

Dependence of shear wave seismoelectrics on soil textures: a numerical study in the vadose zone

In this work, we study seismoelectric conversions generated in the vadose zone, when this region is traversed by a pure SH wave. We assume that the soil is a 1-D partially saturated lossy porous medium and we use the van Genuchten's constitutive model to describe the water saturation profile. Correspondingly, we extend Pride's formulation to deal with partially saturated media. In order to evaluate the influence of different soil textures we perform a numerical analysis considering, among other relevant properties, the electrokinetic coupling, coseismic responses and interface responses (IRs). We propose new analytical transfer functions for the electric and magnetic field as a function of the water saturation, modifying those of Bordes et al. and Garambois & Dietrich, respectively. Further, we introduce two substantially different saturation-dependent functions into the electrokinetic (EK) coupling linking the poroelastic and the electromagnetic wave equations. The numerical results show that the electric field IRs markedly depend on the soil texture and the chosen EK coupling model, and are several orders of magnitude stronger than the electric field coseismic ones. We also found that the IRs of the water table for the silty and clayey soils are stronger than those for the sandy soils, assuming a non-monotonous saturation dependence of the EK coupling, which takes into account the charged air-water interface. These IRs have been interpreted as the result of the jump in the viscous electric current density at the water table. The amplitude of the IR is obtained using a plane SH wave, neglecting both the spherical spreading and the restriction of its origin to the first Fresnel zone, effects that could lower the predicted values. However, we made an estimation of the expected electric field IR amplitudes detectable in the field by means of the analytical transfer functions, accounting for spherical spreading of the SH seismic waves. This prediction yields a value of 15 μV m-1, which is compatible with reported values.Facultad de Ciencias Astronómicas y GeofísicasConsejo Nacional de Investigaciones Científicas y Técnica

An analytical solution to assess the SH seismoelectric response of the vadose zone

We derive an analytical solution of the seismoelectric conversions generated in the vadose zone, when this region is crossed by a pure shear horizontal (SH) wave. Seismoelectric conversions are induced by electrokinetic effects linked to relative motions between fluid and porous media. The considered model assumes a 1D soil constituted by a single layer on top of a half-space in contact at the water table, and a shearing force located at the earth?s surface as the wave source. The water table is an interface expected to induce a seismoelectric interfacial response (IR). The top layer represents a porous rock in which porous space is partially saturated by water and air, while the half-space is completely saturated with water, representing the saturated zone. The analytical expressions for the coseismic fields and the interface responses, both electric and magnetic, are derived by solving Pride's equations with proper boundary conditions. An approximate analytical expression of the solution is also obtained, which is very simple and applicable in a fairly broad set of situations. Hypothetical scenarios are proposed to study and analyse the dependence of the electromagnetic fields on various parameters of the medium. An analysis of the approximate solution is also made together with a comparison to the exact solution. The main result of the present analysis is that the amplitude of the interface response generated at the water table is found to be proportional to the jump in the electric current density, which in turn depends on the saturation contrast, poro-mechanical and electrical properties of the medium and on the amplitude of the solid displacement produced by the source. This result is in agreement with the one numerically obtained by the authors, which has been published in a recent work. We also predict the existence of an interface response located at the surface, and that the electric interface response is several orders of magnitude bigger than the electric coseismic field, whereas it is the opposite using compressional waves as shown by theoretical and experimental results. This fact should encourage the performance of field and laboratory tests to check the viability of SHTE seismoelectrics as a near surface prospecting/monitoring tool.Facultad de Ciencias Astronómicas y Geofísica

Reply to ‘Comment on “Dependence of shear wave seismoelectrics on soil textures: a numerical study in the vadose zone by F.I. Zyserman, L.B. Monachesi and L. Jouniaux” by Revil, A.’

In this paper we reply to a the comment made by Revil (2017) on our paper (2017, Geophys. J. Int., 208), where we describe seismoelectric phenomena in the vadose zone based on the theory of Pride empirically extended for unsaturated conditions. We analyse and answer each one of the enumerated critics, and reaffirm the conclusions of our work. In particular, we prove that using the conductivity model suggested by Revil (2017) does not change our predictions significantly, contrary to what was argued in the comment. Further, in the light of previous and new theoretical and experimental results existing in the literature, we confirm the reasonability of having tested a non-monotonic saturation dependent streaming potential coefficient model besides the monotonic one, and discuss the suitability of assuming a linear relation between the permeability and the excess charge.Facultad de Ciencias Astronómicas y GeofísicasConsejo Nacional de Investigaciones Científicas y Técnica

A 1-D modelling of streaming potential dependence on water content during drainage experiment in sand

The understanding of electrokinetics for unsaturated conditions is crucial for numerous of geophysical data interpretation. Nevertheless, the behaviour of the streaming potential coefficient C as a function of the water saturation Sw is still discussed. We propose here to model both the Richards' equation for hydrodynamics and the Poisson's equation for electrical potential for unsaturated conditions using 1-D finite element method. The equations are first presented and the numerical scheme is then detailed for the Poisson's equation. Then, computed streaming potentials (SPs) are compared to recently published SP measurements carried out during drainage experiment in a sand column. We show that the apparent measurement of DV / DP for the dipoles can provide the SP coefficient in these conditions. Two tests have been performed using existing models for the SP coefficient and a third one using a new relation. The results show that existing models of unsaturated SP coefficients C(Sw) provide poor results in terms of SP magnitude and behaviour. We demonstrate that the unsaturated SP coefficient can be until one order of magnitude larger than Csat, its value at saturation. We finally prove that the SP coefficient follows a non-monotonous behaviour with respect to water saturation. Key words: Electrical properties; Electromagnetic theory; Hydrogeophysics; Hydrology; Permeability and porosity; electrokinetic; streaming potential; self-potential; water content; water saturation; unsaturated condition; finite element modelin

A review on electrokinetically induced seismo-electrics, electro-seismics, and seismo-magnetics for Earth sciences

The seismo-electromagnetic method (SEM) can be used for non-invasive subsurface exploration. It shows interesting results for detecting fluids such as water, oil, gas, CO2, or ice, and also help to better characterise the subsurface in terms of porosity, permeability, and fractures. However, the challenge of this method is the low level of the induced signals. We first describe SEM's theoretical background, and the role of some key parameters. We then detail recent studies on SEM, through theoretical and numerical developments, and through field and laboratory observations, to show that this method can bring advantages compared to classical geophysical methods.Facultad de Ciencias Astronómicas y Geofísica