Patella.pdf 589 ANNALS OF GEOPHYSICS, VOL. 46, N. 3, June 2003 On the role of the J-E constitutive relationship in applied geoelectromagnetism Domenico Patella Dipartimento di Scienze Fisiche, Università «Federico II», Napoli, Italy Abstract In current applications of the Induced Polarization (IP) method, the Debye and Cole-Cole models are used to study relaxation and dispersion properties of rocks, though it is believed that this type of modelisation is confused and vague, because of the lack of a background physical description. In this paper, we show that the Debye model can physically be deduced as a consequence of the electrodynamic behaviour of a mixture of bound and unbound charged particles immersed in an external electric field. We also clarify that the Cole-Cole model is a synthetic model, which can physically be explained as a continuous distribution of Debye terms. Mailing address: Prof. Domenico Patella, Dipartimen- to di Scienze Fisiche, Università «Federico II», Complesso Universitario Monte Sant’Angelo, Via Cintia, 80126 Napoli, Italy; e-mail: patella@na.infn.it Key words rock electrodynamics – induced pola- rization – impedivity models 1. Introduction Electrical relaxation and dispersion in rocks are observed in the Time-Domain (TD) and Fre- quency-Domain (FD), respectively, using any standard Induced Polarization (IP) device. These effects are often investigated in mining and en- vironmental exploration. Figure 1a-c shows schematically the IP re- sponse of a polarizable medium, prospected, e.g., by a quadrupolar electrode device (fig. 1a). In the TD, a voltage transient is detected across the two receiving electrodes, after both the onset and the shutdown of a current step excitation through the two emitting electrodes, as depicted in fig. 1b. A remarkable IP feature is that at the onset of the excitation current, the voltage response suddenly jumps up from zero to a finite value and then gradually rises towards a stable value. Accordingly, at the shutdown of the current, the voltage suddenly falls down to a finite non-null value and then gradually vanishes. In the FD, one observes a voltage-to-current ratio which is a complex function of frequency, showing amplitude and phase spectra as in fig. 1c. Worth noting is the high-frequency asymptote of the amplitude curve, which is a horizontal straight- line generally placed at a finite non-vanishing level, lower than that corresponding to the low- frequency asymptote. Despite such different behaviour, relaxation in TD and dispersion in FD represent the same physical phenomenon and the relative responses are connected to each other via Fourier Transform (FT). Thus, IP observations are the evidence of a complex form of conduction of electricity in rocks, which cannot be synthesized by Ohm’s law. In fact, Ohm’s law states that the current density is linearly related to the electrical field by a factor s, known as the conductivity, which is assumed independent of time t in the TD, or frequency w in the FD. 590 Domenico Patella c Rigorously speaking, Ohm’s law can be assumed valid in rocks only when a Direct Current (DC) regime is utilized, as in the DC resistivity method. With variable currents, as in many natural or artificial source TD or FD electromagnetic (EM) methods, Ohm’s law can be adopted provided that relaxation and dis- persion effects are negligible, but in general this condition cannot be established beforehand. Therefore, in order to deal with the influence of IP in all EM methods, a generalized current density-to-electrical field relationship, including Ohm’s law as particular case, is needed. In earth materials, many mechanisms have been developed to explain IP effects (Marshall and Madden, 1959; Nilsson, 1971; Zonge and Wynn, 1975; Wong, 1979; Klein et al., 1984; Olhoeft, 1985; Wyller et al., 1992). Comprehensive treatments can be found in the review books by Wait (1959), Bertin and Loeb (1976), Sumner (1976) and Fink et al. (1990). Of the many empirical laws proposed so far, the most utilised in geophysics are the Debye (1928) TD and FD pair (see also Wait, 1959; Patella and Di Maio, 1989) (1.1a) (1.1b) where a n 0, tn 0, bn 0, Sa n = m with 0 £ m £ 1, b n = a n t n , and the Cole and Cole (1941) TD and FD pair (see also Pelton et al., 1978) (1.2a) (1.2b) with 0 £ c £1, t 0. In both pairs of formulas, t and w are the time and frequency, respectively, V(t) is the voltage during the TD relaxation phase, V0 is the steady Fig. 1a-c. The Induced Polarization phenomenology observed with a quadrupolar electrode array (a) in the time-domain (b) and in the frequency-domain (c). a b V t V tD n n n N ( ) = -( ) = Â0 1 a texp / r w r wb wt D n nn N i i ( ) = - + Ê ËÁ ˆ ¯̃= Â0 1 1 1 V t mV t nc CC n nc n ( ) = -( ) ( ) +( )= • Â0 0 1 1 / t G r w r w t w t CC c c m i i ( ) = - ( ) + ( ) È Î Í ˘ ˚ ˙0 1 1 591 On the role of the J-E constitutive relationship in applied geoelectromagnetism voltage at the end of the TD charging phase, m is the chargeability (Seigel, 1959), r(w) is the electric field-to-current density ratio in the FD, called impedivity (Patella, 1993), and r0 is the zero-frequency impedivity, i.e. the DC resistivity. Finally, G(.) is Euler’s function, i is the imaginary unit, and a n , t n , b n and c, t, m are assumed as heuristic parameters required to adapt the formulas to experimental data. It must be mentioned that the Cole-Cole model includes, as particular cases, Warburg (1899) and Madden and Cantwell (1967) models, for c = 0.5 and c = 0.25, respectively. Eqs. (1.1a) and (1.2a) can be derived from (1.1b) and (1.2b), respectively, using the relationship (1.3) For a background knowledge of the analytical properties of the Debye and Cole-Cole formulas, reference is made to Shuey and Johnson (1973), Patella and Ciminale (1979), Pelton et al. (1983, 1984) and Caputo (1993). However, Wait (1982) points out that this characterization of IP is rather confused and vague, because fundamentally non- physical descriptions are employed. To try to overcome this conceptual drawback, in this paper a mechanical approach is suggested to derive a generalised current density-to-electrical field relationship, allowing the physical properties of Debye and Cole-Cole descriptions to be fully investigated. 2. The J-E constitutive relationship in presence of polarization We consider Ampère-Maxwell equation in the TD written as (Stratton, 1941) (2.1) where h is the magnetic field, and d' and j' are the bound charges displacement field and the unbound charges current density, respectively. Postulating that j' and e, and d' and e are related to each other by a convolutional form typical of a linear, causal and time-invariant system, eq. (2.1) can be expanded as (2.2) where the causal functions σ̃ ' 't( ) and ε̃' 't( ) represent the electrical conductivity and dielectric permeability impulse responses, respectively. Applying the derivation rule of a convolution integral and putting (2.3) Ampère-Maxwell equation takes the compact form (2.4) The right-hand integral in eq. (2.4) defines the total current density vector j(t), say (2.5) As a consequence of eq. (2.3), the causal function σ̃ t( ) , which is called the admittivity impulse response, includes the electrical conductivity σ̃ ' t( ) and dielectric permeability ε̃' t( ) impulse responses. In the FD, putting with J(w), E(w), s(w), s'(w) and e'(w) the FT of j(t), e(t), σ̃ t( ) , σ̃ ' t( ) and ε̃' t( ), respectively, eq. (2.3) and eq. (2.5) are transformed, respectively, into (2.6) and . (2.7) Equation (2.7) can equivalently be put in the form (2.8) V t V e d i i t ( ) = − ( ) −∞ +∞ ∫ 0 0 0 1 2π ρ ρ ω ρ ω ω ω . ∇ × ( ) = ( ) + ∂ ( ) ∂h j dt t t t' ' / ∇ × ( ) = ( ) −( ) + ∂ ∂ ( ) −( ) ∞ ∞ ∫ ∫h e et t t t dt t t t t dt˜ ˜σ ε' ' ' ' ' ' ' ' 0 0 ˜ ˜ ˜ /σ σ εt t t t' ' ' ' ' '( ) = ( ) + ∂ ( ) ∂[ ] ∇ × ( ) = ( ) −( ) ∞ ∫h et t t t dt˜ .σ ' ' ' 0 j ' e ' 't t t t dt( ) = ( ) −( ) ∞ ∫ ˜ .σ 0 σ ω σ ω ω ε ω( ) = ( ) + ( )' 'i J Eω σ ω ω( ) = ( ) ( ) E Jω ρ ω ω( ) = ( ) ( ) 592 Domenico Patella where r(w) = 1/s(w). The frequency-dependent functions r(w) and s(w) are called here im- pedivity and admittivity, respectively, and represent a useful extension of the classical parameters of the resistivity r and conductivity s, respectively. Using now the identity (2.9) equation (2.2) can also be written as (2.10) Putting (2.11) Ampère-Maxwell eq. (2.1) finally takes also the compact form (2.12) whose right-hand integral defines now the total displacement field vector d(t), say (2.13) According to eq. (2.11), the causal function ε̃ t( ) , which is called permittivity impulse response, also includes the electrical conductivity and dielectric permeability impulse responses σ̃ ' t( ) and ε̃' t( ). In the FD, using the FT, eq. (2.11) and eq. (2.13) are transformed, respectively, into (2.14) and (2.15) where d(w) is the Dirac delta function and D(w) is the FT of d(t). Equation (2.15) can equivalently be put in the form (2.16) where h(w) = 1/e(w). The frequency-dependent parameters h(w) and e(w) are here called pre- ventivity and permittivity, respectively, and re- present an extension of the classical concept of dielectric permeability e. Equation (2.13) is fully equivalent to eq. (2.5): in fact, comparing eq. (2.3) with eq. (2.11), we readily derive the transformation formulas (2.17) (2.18) and finally, comparing eq. (2.4) with eq. (2.12) and eq. (2.5) with eq. (2.13), we readily get (2.19) Concluding, we can state that for linear, causal and time-invariant natural systems, electrical current density and displacement field are equivalent descriptive vectors of the electromagnetic prop- erties of matter. The bulk response can indif- ferently be synthesized either by the admittivity or the permittivity parameters, which are linked together with simple transformation formulas. By this formulation, a distinction between unbound and bound charge carriers, in relation to the long- or short-distance capacity of movement under the influence of an external force, becomes un- essential. In order to make Ampère-Maxwell equation applicable in practice we must define the physical model of the admittivity (permittivity) TD and FD functions, for which we need to study at first the electrodynamic behaviour of a charged particle. To this aim, it is worth explaining beforehand how we shall combine the bulk electrical properties of a geophysical ˜ ˜σ σ τ τ' ' ' ' ' t t d t ( ) = ∂ ∂ ( ) −∞ ∫ ∇ × ( ) = ∂ ∂ ( ) + ( )       −( ) −∞ ∞ ∫∫h ' ' ' e ' 't t d t t t dt t ˜ ˜ . ' σ τ τ ε 0 ˜ ˜ ˜ε σ τ τ εt d t t ' ' ' ' ' ( ) = ( ) + ( ) −∞ ∫ ∇ × ( ) = ∂ ∂ ( ) −( ) ∞ ∫h et t t t t dtε̃ ' ' ' 0 d et t t t dt( ) = ( ) −( ) ∞ ∫ ˜ .ε ' ' ' 0 ε ω πσ δ ω σ ω ω ε ω( ) = ( ) ( ) + ( ) + ( )' ' '0 i D Eω ε ω ω( ) = ( ) ( ) E Dω η ω ω( ) = ( ) ( ) ˜ ˜ /σ εt t t( ) = ∂ ( ) ∂ ˜ ˜ε σ τt t d t ( ) = ( ) −∞ ∫ j dt t t( ) = ∂ ( ) ∂/ . 593 On the role of the J-E constitutive relationship in applied geoelectromagnetism readily obtain the solution for R(w) as (3.3) Supposing that the system has K > 0 equal charged particles per unit of volume, by definition it is (Stratton, 1941; Parkhomenko, 1967) (3.4) which, using eq. (3.4), becomes (3.5) Comparing eq. (2.7) with eq. (3.5), we get (3.6) Though limited to only one species of charge car- rier, eq. (3.6) is a useful admittivity starting model. It describes a circuit-like cell with a resistance- capacitance-inductance series combination. In fact, where it is wm2 << m1 and m0 << wm1, eq. (3.6) reduces to (3.7) which describes a purely resistive cell, as as- sumed in DC geoelectrical methods. Where, instead, it is wm1 << m0 and w 2 m2 << m0 eq. (3.6) reduces to (3.8) describing a purely capacitive cell. Finally, if structure and the microscopic electrodynamic properties of a compound of charged particles. Indeed, in all equations so far developed, we have not indicated the spatial dependence of the constitutive parameters and field vectors, tac- itly admitting that the medium is a continuum. This is a useful, practical assumption at the scale of the geophysical field and laboratory experiments. However, as will be discussed later, the physical derivation of the analytical expression of the constitutive parameters requires a zooming at the smallest scale of the charge carriers, for which discontinuous properties must be introduced. Thus, the medium will be considered as an aggregate of elementary cells with unitary volume, each cell being sufficiently small to assimilate the medium to a continuum and sufficiently large to incorporate the full diversity of a charged mixture. 3. The electrodynamic equation of a charged particle Let us consider a charged particle confined in a medium immersed in an external electric field e(t). The trajectory r(t ) of the particle can be determined by solving the following general mechanical differential equation (3.1) where q and m2 are the charge and mass of the particle, m1 is a friction-like coefficient accounting for dissipative effects due to collisions and in- teractions with other charged particles, and m0 is an elastic-like coefficient accounting for recall effects on bound charges. Obviously, ml 0 for any l. In the FD eq. (3.1) transforms to (3.2) where R(w) is the FT of r (t). From eq. (3.2) we m d t dt q tλ λ λ λ = ∑ ( ) = ( ) 0 2 r e m i qλ λ λω ω ω = ∑ ( ) ( ) = ( ) 0 2 R E R Eω ω ω λ λ λ ( ) = ( ) ( ) = ∑ q m i 0 2 . J Rω ω ω( ) = ( )Kqi J Eω ω ω ω λ λ λ ( ) = ( ) ( ) = ∑ i Kq m i 2 0 2 . σ ω ω ωλ λ λ ( ) = ( ) = ∑ i Kq m i 2 0 2 . σ ω( ) = Kq m 2 1 σ ω ω( ) = i Kq m 2 0 594 Domenico Patella m1 << wm2 and m0 << w 2m2 eq. (3.6) reduces to (3.9) which represents the admittivity of a purely inductive cell. From eq. (3.6) more complex cases can also be discussed. If unbound charges (m0 = 0) are considered, the admittivity model results to be (3.10) which fits the behaviour of a circuit-like cell consisting of a resistance and inductance in series. This model is assumed to explain dispersion effects in metals (Stratton, 1941). For bound charges (m0 0) with negligible inertia (m2 ª 0) the admittivity reduces to (3.11) which characterises a circuit-like cell made of a resistance and capacitance in series. Finally, for bound charges (m0 0) with vanishing friction (m1 ª 0) the admittivity is (3.12) which characterises a circuit-like cell made of an inductance and capacitance in series. 4. The Debye impedivity model Let us consider now a system with two dif- ferent species of charge carriers and put with Kj, qj and mlj the number per unit of volume, the electrical charge and the passive coefficients, respectively, of the carriers of the j-th species ( j = 1,2). Such a pair of ionic species can, e.g., be the result of ionic dissociation of a salt dissolved in pore water. Since the total current density in the FD is now defined as (Parkhomenko, 1967) (4.1) we readily obtain the following analytical model of admittivity (4.2) For our purpose, the most interesting case to investigate is when one species ( j = 1) consists of unbound charges and the other species ( j = 2) of bound charges, both showing negligible inertia. This simplified physical set can already explain, at least qualitatively, non-resonant IP phenomena in rocks (Keller and Frischknecht, 1966). In fact, we have a two-branch parallel circuit-like cell, with a branch consisting of a single resistance and the other branch of a resistance-capacitance series combination, i.e. (4.3) Putting (4.4) (4.5) (4.6) with a,t, b 0 and t b, from eq. (4.3) we obtain the admittivity function as (4.7) and the impedivity function as (4.8) We have thus got admittivity and impedivity functions qualitatively reflecting the IP frequency spectral evolution from s s w w0 0 = ( ) = Æ lim a, or r0 = r w w 0 1= ( ) = Æ lim a , to s s w t t w• Æ• = ( ) = -( )lim a b a, or r r w t t w• Æ• = ( ) = -( ) £lim .b a a1 s w w w ( ) = + i Kq m i m 2 0 1 s w w w ( ) = - i Kq m m 2 0 2 2 J Rw w w( ) = ( ) = Âi K qj j j j 1 2 s w w w w ( ) = + -= Â i K q m i m m j j j j jj 2 0 1 2 21 2 s w w w ( ) = + + K q m i K q m i m 1 1 2 11 2 2 2 02 12 . a K q m = 1 1 2 11 t = +K q m K q m K q m 1 1 2 12 2 2 2 11 1 1 2 02 b K q m K q m = 2 2 2 11 1 1 2 02 s w wt w t ( ) = +( ) + -( ) a i i b 1 1 r w w t wt ( ) = + -( ) +( ) 1 1 i b a i . s w w ( ) = Kq i m 2 2 s w w ( ) = + Kq m i m 2 1 2 595 On the role of the J-E constitutive relationship in applied geoelectromagnetism Generalising, we can now consider within each elementary cell either a parallel or a series combination of N two-branch circuits, or a net- work of both combinations. For a parallel com- bination, the admittivity becomes (4.9) whereas for a series combination, the impedivity becomes (4.10) which takes the same form as in eq. (1.1b), putt- ing (4.11) (4.12) 5. The Cole-Cole impedivity model Shuey and Johnson (1973) showed that the analysis of the decay spectrum allows the pro- perties of an impedivity model to be more easily deduced and understood. The decay spectrum A(k) is defined as (5.1) where ρ γik m( ) is any given impedivity function, where w is replaced by ik m γ with g an arbitrarily small positive number. Using the Debye impedivity model given in eq. (1.1b), after some simple steps we get the Debye decay spectrum AD(k) as (Patella and Di Maio, 1989) (5.2) The limit in the right-hand side of eq. (5.2) defines the Dirac delta function d(k-1/t n ). Hence, we obtain at last the very important result (5.3) Now, considering the definition of the d-func- tion (Papoulis, 1962) applied to a generic decay spectrum A(k) (5.4) it follows that A(k) can always be considered as a continuous distribution of impulses with varying time constants, and thereof approximable to a sum as in eq. (5.3) with any desired accuracy. In fact, discretising eq. (5.4) using the trapezoidal rule for expressing the area under the curve and putting j = 1/t, we obtain the approximation (5.5) which conforms to eq. (5.3) putting An n n= ρ β τ0 2/ . This conclusion allows Cole-Cole, Warburg and Madden-Cantwell models to be considered as a continuous superposition of Debye terms with varying time constants distributed around k = l/t. Therefore, t in eq. (1.2b) takes the addition- al significance of main relaxation time constant, when the Cole-Cole model is considered as a sum of Debye terms (Pelton et al., 1983), i.e. as . (5.6) 6. Conclusions We have demonstrated that the Debye im- pedivity model for rocks can physically be de- duced as a consequence of the electrodynamic behaviour of a mixture of bound and unbound charged particles immersed in an external electrical field. We have also clarified that the Cole-Cole impedivity model is a synthetic σ ω ωτ ω τ ( ) = + + −( )=∑ a i i b n n n nn N 1 11 ρ ω ω τ ωτ ( ) = + −( ) +( )=∑ 1 11 i b a i n n n nn N ρ ρ ω ω0 0 1 1= ( ) = ( ) → = ∑lim an n N β ρn n nb a= 0 . A k ik ik i ( ) = −( ) − +( ) → lim γ ρ γ ρ γ π0 2 A k k D n nn N n ( ) = −( ) +[ ]= →∑ ρ β τ γ π τ γγ 0 2 1 0 2 21 lim . A k kD n nn N n( ) = −( ) = ∑ρ β τ δ τ0 2 1 1 . A k A k d( ) = ( ) −( ) ∞ ∫ ϕ δ ϕ ϕ 0 A k A kn n N n( ) ≈ −( ) = ∑ 1 1δ τ ρ ω ρ ωβ ωτ CC n nn N i i ( ) ≈ − +    = ∑0 1 1 1 596 Domenico Patella model, which can physically be explained as a continuous distribution of Debye terms. Though being the Debye and derived Cole- Cole impedivity models related to the physical parameters describing the motion of ionic spe- cies in rock capillaries under the influence of an external electric field, no possibility exists to distinguish each ionic species in a compound by means of only IP meaurements. Indeed, ad- mitting that one may be able to single out the N exponential decay terms in the TD voltage transient or the N complex fractions in the FD impedivity spectrum and obtain the set of a n ,t n or b n ,t n IP parameters, by no way the mechanical and electric parameters of bound and unbound charged particles can be determined singly. This conclusion contrasts the current belief that using IP it should be possible to distinguish ionic particles and to estimate their abundances, as, e.g., in environmental applications for mon- itoring contaminants in sediments. In this ap- plication field, only vague information on the possible existence of contaminants seems to be deducible from IP experiments. Weller and Börner (1996), e.g., have shown a reliable cor- relation between the imaginary part of ad- mittivity and porespace internal surface in sedi- ments, and concluded that the observed high sensitivity of the imaginary component to changes at the internal surface may be used as an indi- cator for contaminations. Acknowledgements Study performed in the frame of the INTAS project 97-2/612. REFERENCES BERTIN, J. and J. LOEB, (1976): Experimental and Theoretical Aspects of Induced polarization (Gebrüder-Bornträger, Berlin), 2 vols. CAPUTO, M. (1993): Free modes splitting and alterations of electrochemically polarizable media, Rend. Fis. Acc. Lincei, s.9, 4, 89-98. COLE, K.S. and R.H. COLE, (1941): Dispersion and absorption in dielectrics, J. Chem. Phys., 9, 341-351. DEBYE, P. (1928): Polar Molecules (Chemical Catalogue Company, New York). FINK, J.B., E.O. MCALISTER, B.K. STERNBERG, W.G. WIEDUWILT and S.H. WARD (Editors) (1990): In- duced polarization: applications and case histories, in Investigations in Geophysics, Society of Exploration Geophysicists, Tulsa, Okla, vol. 4, pp. 414. KELLER, G.V. and F.C. FRISCHKNECHT (1966): Electrical Methods in Geophysical Prospecting (Pergamon Press, Oxford), pp. 517. KLEIN, J.D., T. BIEGLER and M.D. HORNE (1984): Mineral interface processes in the method of induced polarization, Geophysics, 43, 1105-1114. LANDAU, L.D. and E.M. LIFŠITS (1976): Teoria dei Campi (Editori Riuniti, Roma) (translated from Russian Teorija Polja, MIR, Moscow). MADDEN, T.R. and T. CANTWELL (1967): Induced polar- ization, a review, in Mining Geophysics, SEG, Tulsa, vol. 2, 373-400. MARSHALL, D.J. and T.R. MADDEN (1959): Induced polarization: a study of its causes, Geophysics, 24, 790-816. NILSSON, B. (1971): A new combined resistivity and induced-polarization instrument and a new theory of the induced-polarization phenomenon, Geoexploration, 9, 35-54. OLHOEFT, G.R. (1985): Low frequency electrical properties, Geophysics, 50, 2492-2503. PAPOULIS, A. (1962): The Fourier Integral and its Ap- plications (McGraw-Hill, New York), pp. 318. PARKHOMENKO, E.I. (1967): Electrical Properties of Rocks (Plenum Press, New York), pp. 314. PATELLA, D. (1993): The methodological principles of magnetotellurics on generally dispersive media (in Italian), Ann. Geofi s., 36 (suppl. 5-6), 147-160. PATELLA, D. and M. CIMINALE (1979): The linear system theory applied to the electrical response of earth materials, Boll. Geofi s. Teor. Appl., 21, 227-236. PATELLA, D. and R. DI MAIO (1989): On the analysis of Cole- Cole relaxations transients in the induced polarization prospecting method, in Inverse Modeling in Exploration Geophysics, edited by A.VOGEL, R. GORENFLO, B. KUMMER and C.O. OFOEGBU (F. Vieweg & Sohn, Braunschweig/Wiesbaden), 205-219. PELTON, W.H., S.H. WARD, P.G. HALLOF, W.R. SILL and P.H. NELSON (1978): Mineral discrimination and removal of inductive coupling with multi-frequency IP, Geophysics, 43, 588-603. PE L T O N, W.H., W.R. SI L L and B.D. SM I T H (1983): Interpretation of complex resistivity and dielectric data. Part I, Geophys. Trans., 29, 297-330. PE L T O N, W.H., W.R. SI L L and B.D. SM I T H (1984): Interpretation of complex resistivity and dielectric data. Part II, Geophys. Trans., 30, 11-45. SEIGEL, H.O. (1959): Mathematical formulation and type curves for induced polarization, Geophysics, 24, 547-565. SHUEY, R.T. and M. JOHNSON (1973): On the phenomenology of electrical relaxation in rocks, Geophysics, 38, 37-48. STRATTON, J. (1941): Electromagnetic Theory (McGraw- Hill, New York), pp. 508. SUMNER, J.S. (1976): Principles of Induced Polarization for Geophysical Exploration (Elsevier, Amsterdam), pp. 277. 597 On the role of the J-E constitutive relationship in applied geoelectromagnetism WAIT, J.R. (Editor) (1959): Overvoltage Research and Geophysical Applications (Pergamon Press, Oxford), pp. 158. WAIT, J.R. (1982): Geo-electromagnetism (Academic Press, New York), pp. 268. WARBURG, E. (1899): Über das Verhalten sogennauter unpolarisierbarer Elektroden gegen Wechselstrom, Z. Chem., 67, 493. WELLER, A. and F.D. BÖRNER (1996): Measurements of spectral induced polarization for environmental purposes, Environ. Geol., 27, 329-334. WONG, J. (1979): An electrochemical model of the induced- polarization phenomenon in disseminated sulfide ores, Geophysics, 44, 1245-1265. WYLLER, J., N. WELLANDER, F. LARSON and D. PARASNIS (1992): Burger ’s equation as a model for the IP phenomenon, Geophys. Prosp., 40, 325-341. ZONGE, K.L. and J.C. WYNN (1975): Recent advances and applications in complex resistivity measurements, Geophysics, 40, 851-864.