4.02 Seismic Source Theory R Madariaga, Laboratoire de Ge´ologie Ecole Normale Supe´rieure, Paris, France
ã 2015 Elsevier B.V. All rights reserved.
4.02.1 Introduction 51 4.02.2 Seismic Wave Radiation from a Point Force: The Green Function 52 4.02.2.1 Seismic Radiation from a Point Source 52 4.02.2.2 Far-Field Body Waves Radiated by a Point Force 52 4.02.2.3 The Near Field of a Point Force 53 4.02.2.4 Energy Flow from Point Force Sources 53 4.02.2.5 The Green Tensor for a Point Force 54 4.02.3 Moment Tensor Sources 54 4.02.3.1 Radiation from a Point Moment Tensor Source 55 4.02.3.2 A More General View of Moment Tensors 55 4.02.3.3 Moment Tensor Equivalent of a Fault 56 4.02.3.4 Eigenvalues and Eigenvectors of the Moment Tensor 57 4.02.3.5 Seismic Radiation from Moment Tensor Sources in the Spectral Domain 57 4.02.3.6 Seismic Energy Radiated by Point Moment Tensor Sources 58 4.02.3.7 More Realistic Radiation Models 58 4.02.4 Finite Source Models 59 4.02.4.1 The Kinematic Dislocation Model 59 4.02.4.1.1 Haskell’s rectangular fault model 60 4.02.4.2 The Circular Fault Model 61 4.02.4.2.1 Kostrov’s self-similar circular crack 61 4.02.4.2.2 The kinematic circular source model of Sato and Hirasawa 62 4.02.4.3 Generalization of Kinematic Models and the Isochrone Method 62 4.02.5 Crack Models of Seismic Sources 63 4.02.5.1 Rupture Front Mechanics 64 4.02.5.2 Stress and Velocity Intensity 65 4.02.5.3 Energy Flow into the Rupture Front 65 4.02.5.4 The Circular Crack 65 4.02.5.5 The Dynamic Circular Fault in a Homogenous Medium 68 4.02.6 Conclusions 69 Acknowledgments 70 References 70
4.02.1 Introduction model was established in the early 1960s thanks to the obser- vational work of numerous seismologists and the crucial the- Earthquake source dynamics provides key elements for the oretical breakthrough of Maruyama (1963) and Burridge and prediction of ground motion and to understand the physics Knopoff (1964) who proved that a fault in an elastic model of earthquake initiation, propagation, and healing. The sim- was equivalent to a double-couple source. plest possible model of seismic source is that of a point source In this chapter, we shall review what we believe are the buried in an elastic half-space. The development of a proper essential results obtained in the field of kinematic earthquake model of the seismic source took more than 50 years since the rupture to date. In Section 4.02.2, we review the classical point first efforts by Nakano (1923) and colleagues in Japan. Earth- source model of elastic wave radiation and establish some quakes were initially modeled as simple explosions, then as the basic general properties of energy radiation by that source. In result of the displacement of conical surfaces, and finally as Section 4.02.3, we discuss the now classical seismic moment the result of fast transformational strains inside a sphere. In the tensor source. In Section 4.02.4, we discuss extended kine- early 1950s, it was recognized that P-waves radiated by earth- matic sources including the simple rectangular fault model quakes presented a spatial distribution similar to that pro- proposed by Haskell (1964, 1966) and a circular model that duced by single couples of forces, but it was very soon tries to capture some essential features of crack models. Sec- recognized that this type of source could not explain S-wave tion 4.02.5 introduces crack models without friction as models radiation (Honda, 1962). The next level of complexity was to of shear faulting in the earth. This will help establish some introduce a double-couple source, a source without resultant basic results that are useful in the study of dynamic models of force nor moment. The physical origin of the double-couple the earthquake source.
Treatise on Geophysics, Second Edition http://dx.doi.org/10.1016/B978-0-444-53802-4.00070-1 51 52 Seismic Source Theory
¼k k 4.02.2 Seismic Wave Radiation from a Point Force: where R x x0 is the distance from the source to the The Green Function observation point. Using the inverse Fourier transform, we can transform eqn [4] to the time domain to obtain the final There are many ways of solving the elastic wave equation for result different types of initial conditions, boundary conditions, sources, ð min ðÞt,R=b etc. Each of these methods requires a specific approach for every ðÞ¼1 ½ rrðÞ t ðÞ t t u R, t pr f 1R st d different problem that we would need to study. Ideally, we would 4 R=a like, however, to find a general solution method that would allow 1 1 + ½ ðÞrf rR R stðÞ R=a [5] us to solve any problem by a simple method. The basic building 4pra2 R block of such a general solution method is the Green function, the 1 1 solution of the following elementary problem: find the radiation + ½ f ðÞrf rR R stðÞ R=b 4prb2 R from a point source in a finite heterogeneous elastic medium. For simplicity, we consider first the particular case of a homogeneous This complicated-looking expression can be better under- elastic isotropic medium, for which we know how to calculate the stood considering each of its terms separately. The first line is Green function. This will let us establish a general framework for the near field, which comprises all the terms that decrease at studying more elaborate source models. long distance from the source faster than R 1. The last two lines are the far-field P and S spherical waves that decrease with 4.02.2.1 Seismic Radiation from a Point Source distance like R 1.
The simplest possible source of elastic waves is a point force of 4.02.2.2 Far-Field Body Waves Radiated by a Point Force arbitrary orientation located inside an infinite homogeneous, isotropic elasticpffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi body of density rp,ffiffiffiffiffiffiffiffi and elastic constants l and Much of the practical work of seismology is done in the far m. Let a ¼ ðÞl +2m =r and b ¼ m=r be the P- and S-wave field, at distances of several wavelengths from the source. speeds, respectively. Let us note u(x,t), the particle displace- When the distance R is large, only the last two terms in [5] ment vector. We have to find the solution of the elastodynamic are important. Under what conditions can we neglect the first wave equation: term of that expression with respect to the last two? For that @2 purpose, we notice that in [4], R appears always in the partic- r uxðÞ¼, t ðÞrr l + m ðÞuxðÞ, t + mr2uxðÞ, t + fxðÞ, t [1] ular combination oR/a or oR/b. Clearly, these two ratios @t2 determine the far-field conditions. Since a>b, we conclude under homogeneous initial conditions, that is, that the far field is defined by uxðÞ¼,0 ux_ ðÞ¼,0 0, and the appropriate radiation conditions o at infinity. In [1], f is a general distribution of force density as a R R a 1or l 1 function of position and time. For a point force of arbitrary orientation located at a point x0, the body force distribution is where l¼2pa/o is the wavelength of a P-wave of circular frequency o. The condition for the far field depends there- fxðÞ, t ¼ fstðÞdðÞx x0 [2] fore on the characteristic frequency or wavelength of the where s(t) is the source time function, the variation of the radiation. Thus, depending on the frequency content of the amplitude of the force as a function of time. And f is a unit signal SeðÞo , we will be in the far field for high-frequency vector in the direction of the unit force. waves, but we may be in the near field for the low- The solution of eqn [1] is easier to obtain in the Fourier frequency components. transformed domain. As is usual in seismology, we use the The far-field radiation from a point force is usually written following definition of the Fourier transform and its inverse: in the following, shorter form: ð1 eðÞ¼o ðÞ iot 1 1 ux, ux, t e dt uP ðÞ¼R, t ℜPstðÞ R=a 1 FF 4pra2 R ð [3] 1 1 [6] ðÞ¼ eðÞo iot o 1 1 ux, t ux, e d S ðÞ¼ ℜS ðÞ =b 2p 1 u R, t st R FF 4prb2 R Here and in the following, we will note Fourier transform ℜP ℜS with a tilde. where and are the radiation patterns of P- and S-waves, r ¼ After some lengthy work see, for example, Achenbach respectively. Noting that R eR, the unit vector in the radial (1975), we find the Green function in the Fourier domain: direction, we can write the radiation patterns in the following ℜP ¼ ℜS ¼ ¼ simplified form: fReR and fT f fReR where fR is the 1 1 esðÞo ioR ueðÞ¼R, o f rr 1+ e ioR=a radial component of the point force f and fT its transverse 4pr R o2 a component. o Thus, in the far field of a point force, P-waves propagate the i R ioR=b +1+b e radial component of the point force, while the S-waves prop- agate information about the transverse component of the 1 1 + ðÞrf rR ResðÞo e ioR=a S-wave. Expressing the amplitude of the radial and transverse pra2 4 R component of f in terms of the azimuth y of the ray with 1 1 + ½ f ðÞrf rR R esðÞo e ioR=b [4] respect to the applied force, we can rewrite the radiation pat- 4prb2 R terns in the simpler form Seismic Source Theory 53
P S ℜ ¼ cosyeR, ℜ ¼ sinyeT [7] n As we could expect from the natural symmetry of the prob- lem, the radiation patterns are axially symmetrical about the dS = R2 dW axis of the point force. P-waves from a point force have a typical dipolar radiation pattern, while S-waves have a toroidal R (doughnut-shaped) distribution of amplitudes. V
f 4.02.2.3 The Near Field of a Point Force S When oR/a is not large compared to one, all the terms in eqns [4] and [5] are of equal importance. In fact, both far- and near- field terms are of the same order of magnitude near the point Figure 1 Geometry for computing radiated energy from a point source. source. In order to find the small R behavior, it is preferable to The source is at the origin and the observer at a position defined go back to the frequency domain expression [4]. When R!0, by the spherical coordinates R,y,f, distance, polar angle, and azimuth. the term in brackets in the first line tends to zero. In order to calculate the near-field behavior, we have to expand the expo- ð 1 nentials to order R2. After some algebra, we find KtðÞ¼ ru_ 2dV [11] 2 V 1 1 ÂÃÀÁÀÁ uðÞR, t 5 ðÞrf rR R 1 b2=a2 + f 1+b2=a2 stðÞ prb2 while 8 R ð ÂÃ [8] 1 2 UtðÞ¼ lðÞr u +2me e dV [12] 2 ij ij This is the product of the source time function s(t) with the V static displacement produced by a point force of orientation f. is the strain energy change inside the same volume. The last term This is one of the most important results of static elasticity and in [10] is the rate of work of the force against elastic displace- is frequently referred to as the Kelvin or Somigliana solution ment. Equation [10] is the basic energy conservation statement (Aki and Richards, 2002). for elastic sources. It says that the rate of energy change inside The result [8] is quite interesting and somewhat unex- the body V is equal to the rate of work of the sources f plus the pected. The radiation from a point source decays like R 1 in inward energy flow across the boundary S. This energy balance the near field, exactly like the far-field terms. This result has equation can be extended to study energy flow for moment been remarked and extensively used in the formulation of tensor sources and for cracks (Kostrov and Das, 1988). regularized boundary integral equations for elastodynamics Let us note that in [10], energy flows into the body through (Fukuyama and Madariaga, 1995; Hirose and Achenbach, the surface S. In seismology, however, we are interested in the 1989). radiated seismic energy, which is the energy that flows out of the elastic body. The radiated seismic energy until a certain time t is ð ð t ðÞ¼ s _ 4.02.2.4 Energy Flow from Point Force Sources Es t dt ijuinjdS 0 S ð ð t _ A very important issue in seismology is the amount of energy ¼ KtðÞ DUtðÞ+ dt fiuidV [13] radiated by seismic sources. The flow of energy across any 0 V surface that encloses the point source must be the same, so where DU(t) is the strain energy change inside the elastic body that seismic energy is defined for any arbitrary surface. Let us from time t¼0tot.Ift is sufficiently long, so that all motion take the scalar product of eqn [1] with the particle velocity u_ inside the body has ceased, K(t)!0 and we get the simplest and integrate on a volume V enclosing the single point source possible expression located at x0: ð1 ð ð ð ð _ Es ¼ EsðÞ¼ 1 DU + dt fiuidV [14] _ € _ _ 0 V ruiuidV ¼ sij, juidV fiuidV [9] V V V Thus, total energy radiation is equal to the decrease in where we use dots to indicate time derivatives and the summa- internal energy plus the work of the sources against the elastic tion convention on repeated indices. In [9], we have rewritten deformation. Both terms in [14] contribute to seismic radia- the right-hand side of [1] in terms of the stresses tion as we will discuss in more detail for moment tensor s ¼le d me e ¼ @ @ ij ii ij +2 ij, whereÀÁ the strains are ij 1/2( jui + jui). sources and cracks. _ _ _ Using sij, jui ¼ sijui , j sijeij and Gauss’ theorem, we get Although we can use [14] to compute the seismic energy, it is the energy flow identity easier to evaluate it directly from the first term in [13] assuming ð ð that S is very far from the source, so that far-field approximation d ðÞ¼KtðÞ+ UtðÞ s u_ n dS + f u_ dV [10] [6] can be used in [13]. Consider, as shown in Figure 1, a cone dt ij i j i i S U of rays of cross section dO issued from the source around the where n is the outward normal to the surface S (see Figure 1). direction y,f. The energy crossing a section of this ray beam at In [10], K is the kinetic energy contained in volume V at time t: distance R from the source per unit time is given by the energy 54 Seismic Source Theory
_ _ 2 _ flow per unit solid angle es ¼ sijuinjR ,wheresij is the stress, ui Here, d(t) is Dirac’s delta, dij is Kronecker’s delta, and the the particle velocity, and n the normal to the surface dS¼R2dO. comma indicates derivative with respect to the component that _ We now use [6] in order to compute sij and u.Bystraightfor- follows it. ward differentiation and keeping only terms of order 1/R with Similarly, in the frequency domain, _ distance, we get sijnj rcui where rc is the wave impedance and h 1 1 1 c the appropriate wave speed. The energy flow rate per unit solid Ge ðÞ¼xjx ,o ðÞ1+ioRa e ioR=a ij 0 4pr R o2 angle for each type of wave is then , ij i o =b _ 2 _2 +1+ðÞioRb e i R esðÞ¼t rcR u ðÞR, t [15] ÀÁ Using the far-field approximation for the velocity derived 1 1 ioR=a + R,iR, j e from [6] and integrating around the source for the complete 4pra2 R Âà duration of the source, we get the total energy flow associated 1 1 ioR=b + dij R,iR, j e with P- and S-waves: 4prb2 R ð 1 1 EP ¼ RP s 2ðÞt dt for P waves s pra3 2 4 0 ð [16] 1 4.02.3 Moment Tensor Sources S ¼ 1 ℜS 2ðÞ Es 3 2 s t dt for S waves 4prb 0 Ð The Green function for a point force is the fundamental solu- hℜci ¼ p ℜc 2 O where 2 1/(4 ) O( ) d is the mean squared radiation tion of the equation of elastodynamics. However – except for a ¼ pattern for wave c {P,S}. Since the radiation patterns are the few rare exceptions – seismic sources are due to fast internal simple sinusoidal functions listed in [7], the mean square deformation in the Earth, for instance, faulting or fast phase radiation patterns are 1/3 for P-waves and 2/3 for the sum of changes on localized volumes inside the Earth. For a seismic the two components of S-waves. In [16], we assumed that source to be of internal origin, it has to have zero net force and _ðÞ¼ < st 0 for t 0. Finally, it is not difficult to verify that, since zero net moment. It is not difficult to imagine seismic sources _ s has units of force rate, Es andpEffiffiffiffiffiffip have units of energy. Noting that satisfy these two conditions: 3 3 that in the Earth, a is roughly 3b so that a ffi5b ; the amount X of energy carried by the S-waves emitted by a point force of f ¼ 0 source time function s(t) is close to ten times that carried by X [18]  ¼ P-waves. For double-couple sources, this ratio is much larger. f r 0 The simplest such sources are dipoles and quadrupoles. For instance, the so-called linear dipole is made of two identical 4.02.2.5 The Green Tensor for a Point Force point forces of strength f that act in opposite directions at two points separated by a very small distance h along the axes of the The Green function is a tensor formed by the waves radiated forces. The seismic moment of this linear dipole is M ¼fh. from a set of three point forces aligned in the direction of each 0 Experimental observation has shown that linear dipoles of this coordinate axis. For an arbitrary force of direction f, located at sort are not the most frequent models of seismic sources and, point x0 and source function s(t), we define the Green tensor furthermore, there does not seem to be any simple internal for elastic waves by deformation mechanism that corresponds to a pure linear ðÞ¼ðÞ j * ðÞ ux, t G x,t x0,0 f st dipole. It is possible to combine three orthogonal linear dipoles where the star indicates time-domain convolution. in order to form a general seismic source; any dipolar seismic We can also write this expression in the usual index source can be simulated by adjusting the strength of these three notation dipoles. It is obvious, as we will show later, that these three X dipoles represent the principal directions of a symmetrical ten- ðÞ¼ ðÞj * ðÞ ui x, t Gij x,t x0,0 fj st sor of rank two that we call the seismic moment tensor: j 2 3 Mxx Mxy Mxz in the time domain or 4 5 X M ¼ Mxy Myy Myz e e e uiðÞ¼x, o GijðÞxjx0,o fjsðÞo Mxz Myz Mzz j This moment tensor has a structure that is identical to that in the frequency domain. of a stress tensor, but it is not of elastic origin as we shall see The Green function itself can be easily obtained from the promptly. radiation from a point force [6] What do the off-diagonal elements of the moment tensor 1 1 represent? Let us consider a moment tensor such that all ele- GijðÞ¼x,tjx0,0 t½ HðÞ t R=a HtðÞ R=b pr ments are zero except Mxy. This moment tensor represents a 4 R , ij double couple, a pair of two couples of forces that turn in 1 1 ÀÁ + R R dðÞt R=a [17] opposite directions. The first of these couples consists in two 4pra2 R ,i , j forces of direction ex separated by a very small distance h in the Âà 1 1 d dðÞ =b direction y. The other couple consists in two forces of direction + 2 ij R,iR, j t R 4prb R ey with a small arm in the direction x. The moment of each of Seismic Source Theory 55 these couples is M the first pair has positive moment, and the xy, SV P second has a negative one. The conditions of conservation of total force and moment [18] are satisfied so that this source z model is fully acceptable from a mechanical point of view. In SH fact, as shown by Burridge and Knopoff (1964), the double R couple is the natural representation of a fault. One of the pair of forces is aligned with the fault, the forces indicate the direc- y tions of slip, and the arm is in the direction of the fault normal. q
4.02.3.1 Radiation from a Point Moment Tensor Source f Let us now use the Green functions obtained for a point force in order to calculate the radiation from a point moment tensor x source located at point x : 0 Figure 2 Radiation from a point double source. The source is at the
M0ðÞ¼r, t M0ðÞt dðÞx x0 [19] origin and the observer at a position defined by the spherical coordinates R,y,f, distance, polar angle, and azimuth. M0 is the moment tensor, a symmetrical tensor whose compo- nents are independent functions of time. We consider one of the components of the moment tensor, [22], we have introduced the standard notation OðÞ¼t st_ðÞfor for instance, M . This represents two point forces of direction i ij the source time function, the signal emitted by the source in separated by an infinitesimal distance h in the direction j. The j the far field. radiation of each of the point forces is given by the Green The term ℜc(y,f) is the radiation pattern, a function of the function G computed in [17]. The radiation from the M ij ij takeoff angle of the ray at the source. Let (R,y,f) be the radius, moment is then simply X colatitude, and azimuth of a system of spherical coordinates ukðÞ¼x, t Gki, jðÞx,tjx0,t *MijðÞt [20] centered at the source, respectively. It is not difficult to show ij that the radiation pattern is given by
c The complete expression of the radiation from a point ℜ ðÞ¼y, f eR M eR [23] moment tensor source can then be obtained from [17].We for P-waves, where e is the radial unit vector at the source. will be interested only on the far-field terms since the near field R Assuming that the z-axis at the source is vertical, so that y is is too complex to discuss here. We get for the far-field waves measured from that axis, S-waves are given by X 1 1 P _ uPðÞ¼R, t ℜ M ðÞt R=a ℜSVðÞ¼y f ℜSHðÞ¼y f i 4pra3 R ijk jk , ey M eR and , ef M eR [24] Xjk 1 1 [21] uSðÞ¼R, t ℜS M_ ðÞt R=b where ef and ey are unit vectors in spherical coordinates. Thus, i prb3 R ijk jk 4 jk the radiation patterns are the radial components of the moment tensor projected on spherical coordinates. ℜP ¼ ℜS ¼ d where ijk RiRjRk and ijk ( ij RiRj)Rk are the radiation With minor changes to take into account smooth variations patterns of P- and S-waves, respectively. We observe in [21] of elastic wave speeds in the Earth, these expressions are widely that the far-field signal carried by both P- and S-waves is the used to generate synthetic seismograms in the so-called far- time derivative of the seismic moment components, so that far- field approximation. The main changes that are needed are the field seismic waves are proportional to the moment rate of the use of travel time Tc(r,ro) instead of R/c in the waveform source. O(t Tc) and a more accurate geometric spreading factor Very often in seismology, it is assumed that the geometry of g(D,H)/a to replace 1/R, where a is the radius of the Earth the source can be separated from its time variation, so that the and g(D,H) is a tabulated function that depends on the angular moment tensor can be written in the simpler form distance D between the hypocenter and the observer and the source depth H. In most work with local earthquakes, the M0ðÞ¼t M stðÞ approximation [22] is frequently used with a simple correction where M is a time-invariant tensor that describes the geometry for free surface response. of the source and s(t) is the time variation of the moment, the source time function determined by seismologists. Referring to Figure 2, we can now write a simpler form of [21] 4.02.3.2 A More General View of Moment Tensors 1 ℜcðÞy, f ucðÞ¼x, t OðÞt R=c [22] What does a seismic moment represent? A number of mechan- 4prc3 R ical interpretations are possible. In the previous sections, we where R is the distance from the source to the observer. c stands introduced it as a simple mechanical model of double couples for either a for P-waves or b for shear waves (SH and SV). For and linear dipoles. Other authors Eshelby (1956) and Backus P-waves, uc is the radial component; for S-waves, it is the and Mulcahy (1976) have explained them in terms of the appropriate transverse component for SH or SV waves. In distribution of inelastic stresses (sometimes called stress ‘glut’). 56 Seismic Source Theory
Let us first notice that a very general distribution of force that satisfies the two conditions [18] necessarily derives from a symmetrical seismic moment density of the form e I fxðÞ¼r , t MxðÞ, t [25] where M(x,t) is the moment tensor density per unit volume. V Gauss’ theorem can be used to prove that such a force distri- bution has no net force nor moment. In many areas of applied mathematics, the seismic moment distribution is often termed a ‘double layer potential.’ We can now use [25] in order to rewrite the elastodynamic Δs eqn [1] as a system of first-order partial differential equations: M0 @ V r v ¼r s @t @ hi[26] Figure 3 Inelastic stresses or stress glut at the origin of the concept of T _ s ¼ lr vI + m ðÞrv + ðÞrv + M0 seismic moment tensor. We consider a small rectangular zone that @t undergoes a spontaneous internal deformation EI (top row). The elastic where v is the particle velocity and s is the corresponding stresses needed to bring it back to a rectangular shape are the moment tensor or stress glut (bottom row right). Once stresses are relaxed by elastic stress tensor. We observe that the moment tensor den- interaction with the surrounding elastic medium, the stress change is Ds sity source appears as an addition to the elastic stress rate s_. (bottom left). This is probably the reason that Backus and Mulcahy (1976) adopted the term ‘glut.’ In many other areas of mechanics, the 1960s). If the internal strain is produced in the sense of reduc- moment tensor is considered to represent the stresses produced ing applied stress – and reducing internal strain energy – then by inelastic processes. A full theory of these stresses was pro- stresses inside the source will decrease in a certain average posed by Eshelby (1956). Incidentally, the equation of motion sense. It must be understood, however, that a source of internal written in this form is the basis of some very successful numer- origin – like faulting – can only redistribute internal stresses. ical methods for the computation of seismic wave propagation During faulting, stresses reduce in the immediate vicinity of slip see, for example, Madariaga (1976), Virieux (1986), and zones but increase almost everywhere else. Madariaga et al. (1998). We can get an even clearer view of the origin of the moment tensor density by considering it as defining an inelastic strain tensor eI defined implicitly by 4.02.3.3 Moment Tensor Equivalent of a Fault
ðÞ¼ ld eI meI For a point moment tensor of type [27], we can write m0 ij ij kk +2 ij [27] ðÞ¼ ld eI meI dðÞ M0 ij ij kk +2 ij V x x0 [28] Many seismologists have tried to use eI in order to represent seismic sources. Sometimes termed ‘potency’ (Ben Menahem where V is the elementary source volume on which acts the and Singh, 1981), the inelastic strain has not been widely source. Let us now consider that the source is a very thin adopted even if it is a more natural way of introducing seismic cylinder of surface S, thickness h, volume V¼Sh, and unit source in bimaterial interfaces and other heterogeneous media. normal n. Letting the thickness of the cylinder tend to zero, For a recent discussion, see Ampuero and Dahlen (2005). the mean inelastic strain inside the volume V can be computed The meaning of eI can be clarified by reference to Figure 3. as follows: Let us make the following ‘gedanken’experiment. Let us cut an ÂÃ infinitesimal volume V from the source region. Next, we let it eI ¼ 1 D D lim ijh uinj + ujni [29] undergo some inelastic strain eI, for instance, a shear strain due h!0 2 to the development of internal dislocations as shown in the where Du is the displacement discontinuity (or simply the slip) figure. Let us now apply stresses on the borders of the internally across the fault volume. The seismic moment for the flat fault is deformed volume V so as to bring it back to its original shape. If then the elastic constants of the internally deformed volume V have ÂÃÀÁ ðÞ¼ ld D m D D not changed, the stresses needed to bring V back to its original M0 ij ij uknk + uinj + ujni S [30] shape are exactly given by the moment tensor components defined in [27]. This is the definition of seismic moment tensor: so that the seismic moment can be defined for a fault as the It is the stress produced by the inelastic deformation of a body product of an elastic constant with the displacement disconti- that is elastic everywhere. It should be clear that the moment nuity and the source area. Actually, this is the way the seismic tensor is not the same thing as the stress tensor acting on the moment was originally determined by Burridge and Knopoff fault zone. The latter includes the elastic response to the intro- (1964). If the slip discontinuity is written in terms of a direc- D ¼ duction of internal stresses as shown in the last row of Figure 3. tion of slip v and a scalar slip D, ui Dvi, we get ÀÁ The difference between the initial stresses before the internal ðÞ¼ d n l n n m M0 ij ij knk DS + inj + jni DS [31] deformation and those that prevail after the deformed body has been reinserted in the elastic medium is the stress change (or Most seismic sources do not produce normal displacement stress drop as originally introduced in seismology in the late discontinuities (fault opening) so that n n¼0 and the first Seismic Source Theory 57 term in [30] is equal to zero. In that case, the seismic moment better represented by their eigenvalues; separation into a fault tensor can be written as the product of a tensor with the scalar and a CLVD part is generally ambiguous. seismic moment M0 ¼mDS: Seismic moments are measured in units of Nm. Small ÀÁ earthquakes that produce no damage have seismic moments ðÞ¼ n n m M0 ij inj + jni DS [32] less than 1012 Nm, while the largest subduction events, like The first practical determination of the scalar seismic those of Chile in 1960, Alaska in 1964, and Sumatra in 2004, have moments of the order of 1022–1023 Nm. Large destructive moment is due to Aki (1966), who estimated M0 from seismic data recorded after the Niigata earthquake of 1966 in Japan. events like Izmit, Turkey, 1999; Chichi, Taiwan, 1999; and Determination of seismic moment has become the standard Landers, California, 1992 have moments of the order of 20 way in which earthquakes are measured. All sorts of seismo- 10 Nm. logical, geodetic, and geologic techniques have been used to Since the late 1930s, it became commonplace to measure earthquakes by their magnitude, a logarithmic measure of the determine M0. A worldwide catalog of seismic moment tensors was made available online by Harvard University (Dziewonski total energy radiated by the earthquake. Methods for measuring and Woodhouse, 1983). Initially, moments were determined radiated energy were developed by Gutenberg and Richter using by for the limited form [32], but nowadays, the full set of six well-calibrated seismic stations. At the time, the general proper- components of the moment tensor is regularly computed. ties of the radiated spectrum were not known and the concept of Let us remark that the restricted form of the moment moment tensor had not yet been developed. Since at present tensor [32] reduces the number of independent parameters time earthquakes are systematically measured using seismic of the moment tensor from 6 to 4. Very often, seismologists moments, it has become standard to use the following empirical use the simple fault model of the source moment tensor. In relation defined by Kanamori (1977) and Hanks and Kanamori practice, the fault is parameterized by the seismic moment (1979) to convert moment tensors into a magnitude scale: plus the three Euler angles of the fault plane. Following the log 10M0 ðÞ¼in Nm 1:5Mw + 9 [35] convention adopted by Aki and Richards (2002), these angles are defined as d the dip of the fault, f the strike of the fault Magnitudes are easier to retain and have a clearer meaning with respect to the north, and l the rake of the fault, which is for the general public than the more difficult concept of the angle of the slip vector with respect to the horizontal. moment tensor.
4.02.3.4 Eigenvalues and Eigenvectors of the Moment 4.02.3.5 Seismic Radiation from Moment Tensor Sources Tensor in the Spectral Domain Since the moment tensor is a symmetrical tensor of order 3, it In actual applications, the near-field signals radiated by earth- has three orthogonal eigenvectors with real eigenvalues, just quakes may become quite complex because of multipathing, like any stress tensor. These eigenvalues and eigenvectors are scattering, etc., so that the actually observed seismogram, say, the three solutions of u(t), resembles the source time function O(t) only at long periods. It is usually verified that complexities in the wave M v ¼ mv 0 propagation affect much less the spectral amplitudes in the Let the eigenvalues and eigenvector be Fourier transformed domain. Radiation from a simple point moment tensor source can be obtained from [22] by straight- mi, vi [33] forward Fourier transformation: Then, the moment tensor can be rewritten as ℜcðÞy f e 1 0, 0 e ioR=c X ucðÞ¼x, o OðÞo e [36] ¼ T 4prc3 R M0 mivi vi [34] e i where OðÞo is the Fourier transform of the source time func- Each eigenvalue–eigenvector pair represents a linear dipole. tion O(t). A well-known property of Fourier transform is that The eigenvalue represents the moment of the dipole. From the low-frequency limit of the transform is the integral of the extensive studies of moment tensor sources, it appears that source time function, that is, many seismic sources are very well represented by an almost ð1 OeðÞ¼o _ ðÞ ¼ pure double-couple model with m1 ¼ m3 and m2 ffi0. lim M0 t dt M0 o!0 A great effort for calculating moment tensors for deeper 0 sources has been made by several authors in recent years. It So that in fact, the low-frequency limit of the transform of appears that the non-double couple part is larger for these the displacement yields the total moment of the source. sources but that it does not dominate the radiation. For deep Unfortunately, the same notation is used to designate the sources, Knopoff and Randall (1970) proposed the so-called total moment release by an earthquake – M0 – and the time- compensated linear vector dipole (CLVD). This is a simple dependent moment M0(t). linear dipole from which we subtract the volumetric part so From the observation of many earthquake spectra, and that m1 +m2 +m3 ¼0. Thus, a CLVD is a source model where from the scaling of moment with earthquake size, Aki (1967) m2 ¼m3 ¼ 1/2m1. Radiation from a CLVD is very different and Brune (1970) concluded that the seismic spectra decayed from that from a double-couple model, and many seismolo- as o 2 at high frequencies. Although, in general, spectra are gists have tried to separate a double couple from CLVD com- more complex for individual earthquakes, a simple source ponents from the moment tensor. In fact, moment tensors are model can be written as follows: 58 Seismic Source Theory
OðÞ¼o M0 Seismic radiation is equal to the reduction in internal strain 2 [37] 1+ðÞo=o0 energy plus the work of the moment tensor against the elastic strain rate at the source. where o is the so-called corner frequency. In this simple 0 As we have already discussed for a point force, at any ‘omega-squared model,’ seismic sources are characterized by position sufficiently far from the source, energy flow per unit only two independent scalar parameters: the seismic moment solid angle is proportional to the square of local velocity [15]. M and the corner frequency o . Not all earthquakes have 0 0 Using the far-field approximation [36], we can follow the same displacement spectra as simple as [37], but the omega-squared steps as in [16] to express the radiated energy in terms of the model is a convenient starting point for understanding seismic seismic source time function as radiation. ð1 From [37], it is possible to compute the spectra of ground 1 2 E ¼ hiℜc O_ ðÞt dt O_ ðÞ¼o oOðÞ o c pr 5 2 velocity i . Ground velocity spectra have a peak 4 c 0 situated roughly at the corner frequency o . In actual earth- 0 where c stands again for P- or S-waves and, as in [18], quake spectra, this peak is usually broadened and contains Ð hℜii ¼1/(4p) (ℜi)2dO is the mean square radiation pattern. oscillations and secondary peaks; at higher frequencies, atten- 2 O We can now use Parseval’s theorem uation, propagation scattering, and source effects reduce the ð1 ð1 velocity spectrum. _ 2 1 e 2 OðÞt dt ¼ o2OðÞo do Finally, by an additional differentiation, we get the acceler- p € 0 0 ation spectra OðÞ¼ o o2OðÞo . This spectrum has an obvious problem; it predicts that ground acceleration is flat for arbi- in order to express the radiated energy in terms of the seismic trarily high frequencies. The acceleration spectrum usually spectrum as ð1 decays after a high-frequency corner identified as fmax. The 1 E ¼ hiℜc o2O2ðÞo do [41] origin of this high-frequency cutoff was a subject of discussion c p2r 5 2 4 c 0 in the 1990s, which was settled by the implicit agreement that f reflects the dissipation of high-frequency waves due to For Brune’s spectrum [37], the integral is max ð propagation in a strongly scattering medium, like the crust 1 p o2O2ðÞo o ¼ 2o3 and near-surface sediments. d M0 0 0 2 It is interesting to observe that [37] is the Fourier transform of so that radiated energy is proportional to the square of M0o0 jjo OðÞ¼t e 0t [38] moment. We can finally write 2 Ec ¼ 1 hiℜc 2 o3 This is a noncausal strictly positive function, is symmetrical 5 M0 0 [42] M0 16prc about the origin, and has an approximate width of 1/o0.By definition, the integral of the function is exactly equal to M0. This nondimensional relation makes no assumptions about o Even if this function is noncausal, it shows that 1/ 0 controls the rupture process at the source except that the spectrumpffiffiffi is of the width or duration of the seismic signal. At high frequencies, the form [37]. Noting that in the Earth, a is roughly 3b so that the function behaves like o 2. This is due to the slope discon- a5 ffi16b5, the amount of energy carried by the S-waves emitted tinuity of [38] at the origin. by a point moment tensor is close to 25 times that carried by We can also interpret [37] as the absolute spectral ampli- P-waves, if the source spectrum is the same for P- and S-waves. tude of a causal function. There are many such functions, one For finite sources, the partition of energy into P- and S-waves of them – proposed by Brune (1970) –is depends on the details of the rupture process.