adg vol5 n02 Teis 393_399.pdf ANNALS OF GEOPHYSICS, VOL. 45, N. 2, April 2002 393 Electromagnetic radiation related to dislocation dynamics in a seismic preparation zone Roman Teisseyre and Tomasz Ernst Institute of Geophysics, Polish Academy of Sciences, Warszawa, Poland Abstract Electromagnetic emission is observed frequently before earthquakes as high noise level in VLF and ULF bands. We present theoretical considerations on electromagnetic radiation caused by dislocation dynamics in the preseismic micro-sources (micro-crackings) located in an earthquake preparation zone. Some of these micro-sources could be located near the ground surface and their electromagnetic signals could be accessible in some recording stations. The examples of the numerically simulated induction and radiation fields are given and one example of the observed radio-noise recording is shown. 1. Introduction The preseismic micro-source activity may be followed by a rebound motion resulting in the main break: the earthquake. The dynamics of micro-sources is related to the accelerated motion of dislocations and to coalescences of the opposite dislocation arrays (micro-cracks). The theory of stress and dislocation evolution is based on the equation of motion with the self-stress field (continuum with defect distribution) on the slip- fault plane; this equation is supplemented with the source/sink function entering in a similar way as the body forces in the equation of motion for displacements. This theory of stress evolution in a pre-seismic zone combined with the physics of the charge dislocations leads to precursory electromagnetic emission. 2. Evolution of stresses Most breaks start along the pre-existing faults or micro-faults. Considering a plane problem containing a fault, we may assume that the rigidity on a fault µ* is smaller than the bulk value µ of the surrounding rocks: µ* < µ. For the boundary between the fault zone and the intact rocks y = 0, we can take the stress- dislocations relation (see further on) which leads us to the condition *µ* = µ. This is due to the continuity of tangent stress derivatives. We obtain a jump of dislocation density on a fault: áá * > . (2.1) The preseismic micro-source activity may be Mailing address: Prof. Roman Teisseyre, Polish Acade- my of Sciences, Institute of Geophysics, ul. Ks. Janusza 64, 01-402 Warszawa, Poland; e-mail: rt@igf.edu.pl Key words stress evolution – instabilities – charge dislocations – electromagnetic emission µ = 1 S x 394 Roman Teisseyre and Tomasz Ernst followed by a rebound motion resulting in a main break: earthquake. The dynamic processes on a fault, or on micro-faults, can be described by the equation of motion for stresses including the interaction with the self-stress fields (dislocation stresses). We follow the results obtained in the papers by Teisseyre (1996, 2001a,b) and Teisseyre and Yamashita (1999). According to Kossecka and DeWitt (1977), the motion equation is expressed by the balance of divergence of elastic stresses and elastic acceleration. Elastic acceleration is introduced as a difference between the total acceleration and self (plastic) acceleration (the total acceleration can be expressed by time derivatives of displacements). We obtain in this way the motion equation for the elastic continuum with defects (2.2) where S are the elastic stresses, is the elastic velocity, u is the total displacement, S is the self/ plastic velocity (total velocity: u . = + S), is the source/sink function ( . corresponds to a body force). In 1D, a dislocation density is related to stresses S on a fault and also to a slip as follows: (2.3) where S F is the resistance stress (friction stress - in the fracture phase) related to slip progress on fault, is here a certain reference thickness: when computing dislocation density related to slip process we count the dislocations comprised in this thickness. The commonly used stick-slip model, together with the elastodynamic equation, eliminates the self-stress fields from the objects on a fault plane. In our approach the self stresses in interaction with elastic stresses govern a slip evolution process. On a fault (x, t) the dynamic processes cause a rapid increase in the self-stress fields expressed by formation of dislocations and dislocation arrays (cracks – in the fracture phase). The total stresses can be split into elastic stress field and self-stress field, which for this case correspond to wave field and dislocation stress field: S T = S + S S . For dislocation density , dislocation velocity V and also for the self- stresses, we may assume, according to condition (2.1), a rapid decrease away from the fault plane (2.4) where a is a constant. The dislocation velocity can be also related to the local stress field (Mataga et al., 1987). Finally, for the near fault zone we come to the following equations for dislocation density (the 1D case {x, t}; Teisseyre and Yamashita, 1999; Teisseyre, 2001a): (2.5) where c is shear wave velocity. The dynamic evolution is related to the dislocation fields and it depends mainly on the source/sink function . This function describes the nucleation of the new dislocations and the coalescence processes (mutual annihilation of dislocations of the opposite signs – the process equivalent to coalescence of two neighboring dislocated elements or formation of a crack). The source/sink function shall first of all de- scribe the coalescence processes between the dislocation arrays of the opposite signs (the coalescence of the arrays well describes a nucleation process of a crack material frac- turing). We may relate this function to a stress = + t t div d d d d S µ = 1 S, µ = ( )S S F x y t ay x t V x y t ay V x t S x y t ay S x ys , , , , , , , , , ( ) = [ ] ( ) ( ) = [ ] ( ) ( ) = [ ] ( ) exp exp exp 4 4 4 r r r = + 2 2t t t sdiv d d d S u d = 1 , 1 c t x = + 395 Electromagnetic radiation related to dislocation dynamics in a seismic preparation zone surface curvature, or a gradient of dislocation density (inverse of a mean distance between the opposite dislocation groups). Attraction forces between the opposite dislocations can effectively produce a weakening effect, and in this sense the function Π acts equivalently to the friction weakening laws. For a review related to source/sink function Π, see the paper by Teisseyre and Nagahama (1998). Here, we use the modified form of this function (Teisseyre, 2001b) in terms of dis- location density (2.6) The first term describes the formation process of new small elements of dislocations (stresses decrease as 1/r2 ) and proportional to dislocation density α and to square of stress curvature, while the second term describes coalescence process of the opposite dislocation lines (stresses de- crease as 1/r) and proportional to square of dis- location density α 2 and to stress curvature; ε is a constant. 3. Numerical simulations of dislocation evolution Following the paper by Teisseyre (2001b) we briefly present the results of the numerical simulations of dislocation evolution, assuming that a micro-source is represented by a microdomain {x}, {t} of the dimensions {0,2π .10−3[m]}, {0,1.6 . 10−5[s]}, and the initial conditions chosen as follows: According to eq. (2.5) and definition (2.6), we obtain the equation The numerical simulations show some instabilities in time which might be identified with the microseismic events acting in a premonitory domain. The computed fields remain much smoother along the x-axis. 4. Electric current and electromagnetic field at micro-source The electric polarization induced by dislocation processes depends on the square of dislocation density, because both the electric charge and the coalescence break distance between the dislocations are proportional to dislocation density α. For polarization density and its derivatives we obtain while, according to Ogawa et al. (1985), for electromagnetic exitations we have the following formulae expressing the near or electrostatic field, as well as the induction and radiation electric fields related to the dipoles activated in rocks (4.2) Π Π= ∂ ∂ − ∂ ∂         0 2α α εα α( ) x x . ( ) = ( )0 2 0α α πt, , . p = 1 2 2κα , p = καα˙ ˙ , p = +( )2κ αα α˙̇ ˙̇ ˙ (4.1) E p r near = 1 4 0 3πε , E p c r ind = 1 4 0 0 2πε ˙ , E p c r rad = 1 4 0 0 2πε ˙̇ α πx x, .0 0 05 10 16( ) = − −( )[ ] ⋅     0.2 sin m ∂α ∂ x t t ,( ) + α ∂ α ∂ x t x t x . . , , + − ⋅ ( )( ) ( ) +− −10 0 41 0 1 102 6 α αx t x t. , ,+ ⋅ ( ) ( ) =−1 5 10 11 α ∂α ∂ x t x t x . . , , − ⋅ ( ) ( ) ⋅−4 63951 5 10 13 α ∂α ∂ x t x t x , , .⋅ ( ) − ( ) −10 3( )2 396 Roman Teisseyre and Tomasz Ernst where 0 is the dielectric constant; the time derivatives are related to the source {x i 0} and related activity time t0; c0 is light velocity. Of course, we are aware of the fact that the different processes related to excitation of the electromagnetic field take place in a precursory time domain and during an earthquake event (Shevtsova, 1984; Varotsos and Alexopoulos, 1984, 1986; Gokhberg et al., 1985; Varotsos et al., 1992; Park, 1993; Teisseyre, 1992, 1995). For review on the dislocation processes and formation of electric charges, see Whitworth (1975). Here, we try to include also the electro- magnetic emission from the pre-seismic micro- sources and to explain this emission by dislo- cation dynamic processes. 5. Numerical estimations for electric current and electromagnetic radiation We refer to estimations presented by Teisseyre (2001b). An order of magnitude of dislocation density at the time of micro-fracturing may be compared to density at the tip of dislocation array. For the dislocation density treated as a product of line density and the Burgers vector (the definition used in theoretical papers) and referred here to the band width and the n dislocations with the Burgers vector 0 (lattice constant), densely packed in an element d0 = m 0 , we obtain (5.1) while for the density of dislocation lines (the definition of line density used in papers on laboratory experiments, as a reference value we may cite its order for the elastic regime of a crystal body: 1011[m 2]) (5.2) The maximum of this density at local micro- fraturing (n m) may theoretically approach max (corresponding to the line density at a crack tip). The dislocation dynamics, as discussed in our model, can cause the appearance and activation of the electromagnetic dipoles in preseismic sources located in an earthquake preparation zone, but not exactly at the future main seismic spot. Based on numerical sim- ulations for the evolution dislocation fields, we can present here the numerical simulation of the dipole intensities in the pre-seismic sources; the dipole source intensities as a function of time for the induction and radiation terms (4.1) are presented in figs. 1 and 2 in the following scales: in the scale ; for x i = 1,2,3,4 . 10 3[m] and {t} in the scale [10 5 s], and in the scale for x i = 1,2,3,4 . 10 3[m] and {t} in the scale[10 5 s]. These electromagnetic phenomena evidently precede a seismic event; the sharp increases of the electric current source intensities both for the induction and radiation fields are observed. For density electric current, or induction field in a microfracture source we obtain the estimate (Teisseyre, 2001b): while for the source intensity of radiation field x t x ti i, ˙ ,( ) ( ) 10 117 2m s x t x t x ti i i˙ , , ˙̇ ,( ) + ( ) ( )2 10 122 2 2 ; m s j ṗ ˙= = [ ] < [ ]10 100 9 2 6 2A m A m = n d 0 0 , = n m 0 0 1 , max [ ]106 1m = 0 = 0 0 1n m max = [ ]1016 2m . << << 397 Electromagnetic radiation related to dislocation dynamics in a seismic preparation zone Fig. 1. Source intensities of electromagnetic induction field (x,t) . (x,t): the plots of (x i ,t) . (x i ,t) in the scale 1017 [1/m2s]; for x i =1,2,3,4 .10 3[m] and {t} in the scale [10 5s]. Fig. 2. Source intensities of electromagnetic radiation field . 2(x,t)+ (x,t) .. (x,t) : the plots of . 2(x i ,t)+ (x i ,t) .. (x i ,t) in the scale 1022 [1/m2s2]; for x i =1,2,3,4 .10 3[m] and {t} in the scale [10 5s]. 398 Roman Teisseyre and Tomasz Ernst we obtain Such intensity peaks may appear at the isolated small areas of the precursory microfracturings; activation of such local peaks of the and . fields causes the electromagnetic emission. 6. Example of the observed electromagnetic excitations There exist numerous examples of radio-noise excitations before earthquake events; there, we present one example. At L’Aquila Geophysical Observatory (Central Apennines) a radio-noise recording system operating at frequencies of 10, 20 and 40 kHz has been installed (cooperation between the Institute of Geophysics, Polish Academy of Sciences, Warszawa, and the Istituto Nazionale di Geofisica e Vulcanologia, Roma). Several seismic events with magnitudes over 4 and in some cases more than 5 were felt at the end of 1997 at distances of 60-80 km from the Observatory. In fig. 3 (Meloni et al., 2001; Ernst et al., 2001) we present the radio-noise recording (here for 20 kHz; similar records have been obtained for other frequencies) for the time when these seismic events occurred (marked on the plots). We shall note that for several years since the station was installed no such signals have appeared in the radio-noise observations; throughout this time, the seismicity was very low, with earthquakes reaching magnitudes of only M � 3. Fig. 3. VLF/recording (20 kHz) of radio-noise at L’Aquila Observatory (after Meloni et al., 2001, and Ernst et al., 2001): September and October, 1997; two greater seismic events with M 5 are marked by the black circles, straight intervals correspond to the gaps of recording. t j =0 d d ˙̇̇ ˙ ˙̇ .p = +( ) [ ]< [ ]10 109 2 2 11 2A m A sm 399 Electromagnetic radiation related to dislocation dynamics in a seismic preparation zone REFERENCES ERNST, T., J. JANKOWSKI, R. TEISSEYRE, A. MELONI and P. PALANGIO (2001): Observations of electromagnetic radiation in the Central Apennines seismic region, Acta Geophys. Pol., 49, 43-53. GOKHBERG, M.B., I.L. GUFELD, N.I. GERSHENZON and V.A. PILIPENKO (1985): Electromagnetic effects during rupture of the Earth’s crust, Izv. Akad. Nauk SSSR, Fiz. Zemli, 1, 72-87 (in Russian); Izv. Earth Phys., 21 (1), 52-63. KOSSECKA, E. and R. DEWITT (1977): Disclination kinematic, Arch. Mech., 29, 633-651. MATAGA, P.A., L.B. FREUND and J.W. HUTCHINSON (1987): Crack tip plasticity in dynamic fracture, J. Phys. Chem. Solid, 48, 985-1005. MELONI, A., D. DI MAURO, G. MELE, P. PALANGIO, T. ERNST, and R. TEISSEYRE (2001): Evolution of magnetotelluric, total magnetic field, and VLF field parameters in Central Italy: relations to local seismic activity, Ann. Geofis., 44 (2), 383-394. OGAWA, T., K. OIKE and T.MIURA (1985): Electromagnetic radiations from rocks, J. Geophys. Res., 98 (D4), 6245-6249. PARK, S.K., M.S. JOHNSTON, T.R. MADDEN, F.D. MORGAN and F.F. MORRISON (1993): Electromagnetic precursors to earthquakes in the ULF band; a review of observations and mechanisms, Rev. Geophys., 31 (2), 117-132. SHEVTSOVA, I.N. (1984): Charging of dislocations during deformation of crystals with ionic-type bonds, Izv. Earth Phys., 20 (8), 643-648. TEISSEYRE, R. (1992): Earthquake premonitory processes: Evolution of stresses and electric current generation, Terra Nova, 4, 509-513. TEISSEYRE, R. (1995): Electric field generation in earthquake premonitory process, in Theory of Earthquake Premonitory and Fracture Processes, edited by R. TEISSEYRE, 282-292, PWN, Warszawa. TEISSEYRE, R. (1996): Motion and flow equation for stresses, Acta Geophys. Pol., 44, 19-29. TEISSEYRE, R. (2001a): Evolution, propagation and diffusion of dislocation fields, in Earthquake Thermodynamics and Phase Transformations in the Earth’s Interior, edited by R. TEISSEYRE and E. M A J E W S K I (Academic Press, San Diego, San Francisco, NewYork, Boston, London, Sydney, Tokyo), International Geophysical Series, vol. 76, 167-198. TEISSEYRE, R. (2001b): Dislocation dynamics and related electromagnetic excitation, Acta Geophys. Pol., 49, 55-73. TEISSEYRE, R. and H. NAGAHAMA (1998): Dislocation field evolution and dislocation source/sink function, Acta Geophys. Pol., 46, 13-33. TEISSEYRE, R. and T. YAMASHITA (1999): Splitting stress motion equation into seismic wave and fault-related fields, Acta Geophys. Pol., 47, 135-147. VAROTSOS, P. and A. ALEXOPOULOS (1984): Physical properties of the variations of the electric field of the earth preceding earthquakes. II - Determination of epicenter and magnitude, Tectonophysics, 110, 99-125. VAROTSOS, P. and A. ALEXOPOULOS (1986): Thermo- dynamics of Point Defects and their Relation with Bulk Properties (North-Holland Publ. Comp., Amsterdam - New York), pp. 474. VAROTSOS, P., N.G. BOGRIS and A. KYRITSIS (1992): Comments on the depolarization currents stimulated by variations of temperature or pressure, J. Phys. Chem. Solids, 53 (8), 1007-1011. WHITWORTH, R.W. (1975): Charged dislocations in ionic crystals, Adv. Phys., 24 (2), 203-304.