Physical Origins of Anelasticity and Attenuation in Rock I
Total Page:16
File Type:pdf, Size:1020Kb
2.17 Properties of Rocks and Minerals – Physical Origins of Anelasticity and Attenuation in Rock I. Jackson, Australian National University, Canberra, ACT, Australia ª 2007 Elsevier B.V. All rights reserved. 2.17.1 Introduction 496 2.17.2 Theoretical Background 496 2.17.2.1 Phenomenological Description of Viscoelasticity 496 2.17.2.2 Intragranular Processes of Viscoelastic Relaxation 499 2.17.2.2.1 Stress-induced rearrangement of point defects 499 2.17.2.2.2 Stress-induced motion of dislocations 500 2.17.2.3 Intergranular Relaxation Processes 505 2.17.2.3.1 Effects of elastic and thermoelastic heterogeneity in polycrystals and composites 505 2.17.2.3.2 Grain-boundary sliding 506 2.17.2.4 Relaxation Mechanisms Associated with Phase Transformations 508 2.17.2.4.1 Stress-induced variation of the proportions of coexisting phases 508 2.17.2.4.2 Stress-induced migration of transformational twin boundaries 509 2.17.2.5 Anelastic Relaxation Associated with Stress-Induced Fluid Flow 509 2.17.3 Insights from Laboratory Studies of Geological Materials 512 2.17.3.1 Dislocation Relaxation 512 2.17.3.1.1 Linearity and recoverability 512 2.17.3.1.2 Laboratory measurements on single crystals and coarse-grained rocks 512 2.17.3.2 Stress-Induced Migration of Transformational Twin Boundaries in Ferroelastic Perovskites 513 2.17.3.3 Grain-Boundary Relaxation Processes 513 2.17.3.3.1 Grain-boundary migration in b.c.c. and f.c.c. Fe? 513 2.17.3.3.2 Grain-boundary sliding 514 2.17.3.4 Viscoelastic Relaxation in Cracked and Water-Saturated Crystalline Rocks 516 2.17.4 Geophysical Implications 517 2.17.4.1 Dislocation Relaxation 517 2.17.4.1.1 Vibrating string model 517 2.17.4.1.2 Kink model 518 2.17.4.2 Grain-Boundary Processes 519 2.17.4.3 The Role of Water in Seismic Wave Dispersion and Attenuation 519 2.17.4.4 Bulk Attenuation in the Earth’s Mantle 520 2.17.4.5 Migration of Transformational Twin Boundaries as a Relaxation Mechanism in the Lower Mantle? 520 2.17.4.6 Solid-State Viscoelastic Relaxation in the Earth’s Inner Core 521 2.17.5 Summary and Outlook 521 References 522 Nomenclature d grain size; characteristic separation of a spacing between adjacent Peierls vacancy sources valleys f frequency; the reciprocal of du/dx ¼ ad period of kink potential measured tan (Figure 4(c)) parallel to dislocation line g dimensionless multiplier of Tm/T in b Burgers vector exponent for thermal activation of relaxation time 493 494 Physical Origins of Anelasticity and Attenuation in Rock h separation of closely spaced disloca- EP activation energy describing position tions in subgrain wall (period) of dissipation peak hj jth coefficient in Fourier series repre- Eup binding energy of dislocation to pin- senting grain-boundary topography ning point k Boltzmann constant G shear modulus ml effective mass per unit length of Gm free energy of kink migration dislocation GHS mean of the Hashin–Shtrikmann n exponent in term tn representing bounds on the shear modulus of a transient creep in the Andrade model polycrystalline aggregate p internal variable providing link between GR relaxed value of shear modulus (ane- stress and anelastic strain lastic relaxation only) pL(L) distribution of dislocation segment GRS Reuss average of single-crystal elastic lengths constants for shear ps small fluctuation in pressure P GU unrelaxed shear modulus p0 stress-dependent equilibrium value Hint kink interaction energy s Laplace transform variable; distance Hk formation energy for a single geome- measured along dislocation line; trical kink exponent describing grain-size sensi- Hkp kink-pair formation energy À1 tivity of Q Hm enthalpy of kink migration t time Hmc activation enthalpy for multicomponent u displacement of dislocation segment, diffusion or position of grain boundary, in J(t) creep function à y-direction J1 real part of J à u0 x-independent amplitude of u; position J2 negative imaginary part of J of straight dislocation segment under JU unrelaxed compliance ("/) static stress Jà dynamic compliance um maximum displacement for newly Kf bulk modulus of fluid formed kink pair Km bulk modulus of melt wk width of a single geometrical kink Ks bulk modulus of solid matrix x coordinate measuring distance parallel K bulk modulus of composite of coexist- to dislocation line or grain boundary ing phases x1 unstable equilibrium spacing of kinks L length of dislocation segment between of opposite sign adjacent pinning points y coordinate measuring distance per- Lc critical dislocation segment length for pendicular to dislocation line dislocation multiplication A measure of anisotropy or inter-phase L subgrain diameter variability in elastic moduli or thermal N number of terms retained in Fourier expansivity involved in thermoelastic series representing grain-boundary relaxation strength topography B dislocation drag coefficient Nv number of dislocation segments of COH concentration of hydroxyl ions given length L per unit volume D diffusivity (either of matter or heat) P of internal variable p D() normalised distribution of anelastic QÀ1 inverse quality factor (¼ tan), a mea- relaxation times sure of strain energy dissipation À1 Db grain-boundary diffusivity QD height of the Debye dissipation peak of Dk kink diffusivity the standard anelastic solid Dmc multicomponent diffusivity R radius of curvature of dislocation E activation energy for relaxation time bowed out under applied stress; Ed elastic strain energy per unit length of gas constant; radius of grain-edge dislocation or line tension tubule Physical Origins of Anelasticity and Attenuation in Rock 495 Sm entropy of kink migration Poisson’s ratio T absolute temperature k0 attempt frequency of dislocation Tm melting temperature vibration To oscillation period (¼ 2/!) crystal density À þ U equilibrium distance for grain-bound- k , k densities of kinks of opposite sign ary sliding (number per unit length of dislocation U(u) Peierls potential line) À þ X chemical composition k0 , k0 as above, under conditions of zero Y Youngs modulus stress eq negative of exponent describing fre- k equilibrium kink density: number of quency dependence of QÀ1 kinks per unit length of dislocation line b aspect ratio of grain-boundary region stress, usually shear stress (¼ /d) n normal stress e aspect ratio (minimum/maximum P Peierls stress dimension) of ellipsoidal inclusion t tectonic shear stress t aspect ratio of grain-edge tubule up unpinning shear stress (¼ 2R/d) 0 amplitude of sinusoidally time-varying coefficient of term tn representing stress transient creep in the Andrade model relaxation time multiplicative constant of order unity A anelastic relaxation time s multiplicative constant d duration of transient diffusional creep phase lag between stress and strain; dr timescale for draining of cylindrical width of grain-boundary region (<<d) laboratory specimen G anelastic relaxation of the shear e relaxation time for elastically accom- modulus modated grain-boundary sliding J anelastic relaxation of the compliance M Maxwell (viscous) relaxation time P width in pressure of binary loop s,t relaxation time for melt squirt between T0 width of melting interval grain-edge tubules Ts pressure-induced perturbation to soli- s,e relaxation time for melt squirt between dus and liquidus temperatures ellipsoidal inclusions pressure-induced perturbation in melt melt fraction; the angle betweeen a fraction dislocation segment and its Peierls " strain valley "a anelastic strain ! angular frequency "e elastic strain !0 (angular) resonance frequency of dis- "m maximum anelastic strain due to location segment migration of geometrical kinks on a À Gamma function dislocation segment Á anelastic relaxation strength "0 amplitude of sinusoidally time-varying à dislocation density, that is, dislocation strain length per unit volume ¼ NvL viscosity or effective viscosity Ãm density of mobile dislocations 1 viscosity of dashpot in parallel with molecular volume of diffusing species spring in the anelastic element of the (/)0 fractional density contrast between Burgers model coexisting polymorphs b grain-boundary viscosity (/)a fractional change in density resulting f fluid viscosity from small change in proportions of m melt viscosity coexisting phases; anelastic volu- wavelength of periodic grain-boundary metric strain topography or of seismic wave (/)e fractional change in density caused by k kink mobility (¼ Dk/kT) pressure ps 496 Physical Origins of Anelasticity and Attenuation in Rock 2.17.1 Introduction Departures from elastic behavior have long been studied in metals. Notable early successes came from The Earth and its constituent materials are subjected the application of resonance (pendulum) methods at to naturally applied stresses that vary widely in inten- fixed frequencies near 1 Hz that revealed low- sity and timescale. At sufficiently high temperatures, temperature dissipation peaks attributable to various nonelastic behavior will be encountered – even at the types of point-defect relaxation (e.g., Zener, 1952; low stress amplitudes of seismic wave propagation, Nowick and Berry, 1972). Until recently most of the tidal forcing, and glacial loading (Figure 1). For har- available information concerning elastic wave speeds in monic loading, nonelastic behavior is manifest in a geological materials and their variations with pressure phase lag between stress and strain giving rise to strain and temperature came from ultrasonic and opto-acous- energy dissipation and associated frequency depen- tic techniques at frequencies at least six orders of dence (dispersion) of the relevant modulus or wave magnitude higher than those of teleseismic waves speed. In the case of glacial loading, the transient creep (Figure 1). However, there is growing recognition of of mantle material of sufficiently low viscosity allows the need for a thorough understanding of viscoelastic time-dependent relaxation of the stresses imposed by behavior in minerals and rocks at low strains and seis- the growth and decay of ice sheets with consequent mic frequencies.