Earthquake Source and Seismic Wave

Textbooks:
Shearer, P. M. (2019). Introduction to seismology. Cambridge University Press.
GEOPHYS 210: Earthquake Seismology by Eric Dunham
Segall, P. (2010). Earthquake and volcano deformation. Princeton University Press.

Earthquake faults

Earthquakes may be idealized as movement across a planar fault of arbitrary orientation

  • strike: , the azimuth of the fault from north where it intersects a horizontal surface
  • dip: , the angle from the horizontal
  • rake: , the angle between the slip vector and the strike

Earthquake faults

Thrust faulting: reverse faulting on faults with dip angles less than 45
Overthrust faults: Nearly horizontal thrust faults
Strike-slip faulting: horizontal motion between the fault surfaces
Dip-slip faulting: vertical motion
Right-lateral strike–slip motion: standing on one side of a fault, sees the adjacent block move to the right
: left-lateral faulting
: right-lateral faulting
The San Andreas Fault: Right-lateral fault

Earthquake double couple

  • An earthquake is usually modeled as slip on a fault, a discontinuity in displacement across an internal surface in the elastic media.
  • Internal forces resulting from an explosion or stress release on a fault must act in opposing directions so as to conserve momentum.
  • A force couple is a pair of opposing point forces separated by a small distance
  • A double couple is a pair of complementary couples that produce no net torque

Moment tensor

We define the force couple as a pair of equal and opposite forces pointing in the direction and separated by a unit distance in the direction.

The magnitude of is the product of the force and the distance .

The condition that angular momentum be conserved requires that is symmetric (e.g., ).

Moment tensor

For example, right-lateral movement on a vertical fault oriented in the direction corresponds to the moment tensor representation

where is the scalar seismic moment:

where is the shear modulus, is the average fault displacement, and is the area of the fault.

The units for are Nm (or dynecm), the same as for force couples.

Global CMT catalog

Global Centroid Moment Tensor

Beach balls

Moment tensor

Because , there are two fault planes that correspond to a double-couple model.
In general, there are two fault planes that are consistent with distant seismic observations in the double-couple model.
The primary fault plane: The real fault plane
The auxiliary fault plane: The other plane

Eigenvectors

The moment tensor is a symmetric tensor, so it has three real eigenvalues and three orthogonal eigenvectors.

Tension axis
Pressure axis

Non-double couple sources

Isotropic part of :

where

Decomposing into isotropic and deviatoric parts:
where , free from isotropic sources by may contain non-double couple sources

Non-double couple sources

Diagonalize by rotating to coordinates of principal axes:

where .

Because , .
For a pure double couple source, and .

Non-double couple sources

We can decompose into a best-fitting double couple and a non-double couple part :

The complete decomposition of the original is:

Non-double couple sources

Example of the decomposition of a moment tensor into isotropic, best-fitting double couple, and compensated linear vector dipole terms:

The decomposition of into and is unique because we have defined as the best-fitting double couple, i.e., minimizing the CLVD part.

Non-double couple sources

A measure of the misfit between and a pure double-couple source is provided by the ratio to the remaining eigenvalue with the largest magnitude to the largest eigenvalue of

: pure double couple
: close to pure CLVD

Physically, non-double-couple components can arise from simultaneous faulting on faults of different orientations or on a curved fault surface.

Green's function

Review:

The momentum equation:

The stress-strain relation:

The strain:

Green's function

We consider a unit impulse source at point at time .

The unit force function is a useful concept because more realistic sources can be described as a sum of these force vectors.

For every , there is a unique that describes the Earth’s response, which could be computed if we knew the Earth’s structure to sufficient accuracy.

Green's function

We define the notation:

where is the Green's function, is the displacement, and is the force.

Assuming that can be computed, the displacement resulting from any body force distribution can be computed as the sum or superposition of the solutions for the individual point sources.

Green's function + Moment tensor

The moment tensor provides a general representation of the internally generated forces that can act at a point in an elastic medium. Although it is an idealization, it has proven to be a good approximation for modeling the distant seismic response for sources that are small compared to the observed seismic wavelengths.

So the displacement resulting from a force couple at is given by:

Radiation patterns

Solving for the Green's function is rather complicated. Here we consider the simple case of a spherical wavefront from an isotropic source.
Review: The solution for the P-wave potential:

where is the P-wave velocity, is the distance from the point source, and is the source-time function.

The displacement field = the gradient of the displacement potential:

Radiation patterns

Define as the delay time:

So the displacement field is:

The first term decays as , and is called the near-field term, which represents the permanent static displacement due to the source
The second term decays as , and is called the far-field term, which represents the dynamic response (the transient seismic waves radiated by the source that cause no permanent displacement)

Radiation patterns

The first term decays as , and is called the near-field term, which represents the permanent static displacement due to the source

The second term decays as , and is called the far-field term, which represents the dynamic response - the transient seismic waves radiated by the source.
These waves cause no permanent displacement, and their displacements are given by the first time derivative of the source-time function.

Radiation patterns

Most seismic observations are made at sufficient distance from faults that only the far-field terms are important.
Consider the far-field P-wave displacement for a double-couple source, assuming the fault is in the plane with motion in the direction:

The far-field radiation pattern for P-waves

  • The fault plane and the auxiliary fault plane form nodal lines of zero motion.
  • The outward pointing vectors in the compressional quadrant.
  • The inward pointing vectors in the dilatational quadrant.
  • The tension (T axis) is in the middle of the compressional quadrant;
  • The pressure (P axis) is in the middle of the dilatational quadrant.

The far-field radiation pattern for S-waves

The far-field S displacements:

where is the S-wave velocity.

There are no nodal planes, but there are nodal points.
S-wave polarities generally point toward the T axis and away from the P axis.

Beach balls

Plotting beach balls

Projection of the focal sphere onto an equal-area lower-hemisphere map.
The numbers around the outside circle show the fault strike in degrees.
The circles show fault dip angles with 0° dip (a horizontal fault) to 90° dip (a vertical fault).
The curved line ABC shows the intersection of the fault with the focal sphere.

Review: Basic types of faulting

Basic types of faulting

20250330233905

Body force

First-motion polarity

20250330234223

Radiation pattern

20250330234626

Radiation pattern

20250330234649

Earthquake Focal Mechanism

Earthquake Focal Mechanism

Earthquake Focal Mechanism

Earthquake Focal Mechanism

Earthquake Focal Mechanism

Earthquake Focal Mechanism

P/T axes v.s. Compressional/Dilatational quadrants

The Pressure/Tension (P/T) axes are defined inside the beach ball;
The Compressional/Dilatational quadrants are defined outside the beach ball;

The tension (T axis) is in the middle of the compressional quadrant;
The pressure (P axis) is in the middle of the dilatational quadrant.

Moment tensor and Beach ball

20250331000705

Moment tensor and Radiation pattern

Moment tensor decomposition

Magnitude

The magnitude of the equivalent body forces is
The scalar seismic moment of the earthquake; units of dyn-cm, or N-m

Time-dependent seismic moment

The Haskell source model

20250331234705

The rise time: describes the slip duration at any point on a fault.

Rupture propagation

For a long, narrow fault, we assume that the rupture propagates along the fault of length from left to right at a rupture velocity
The total duration of the rupture is

20250331235148

Far field the apparent rupture time

20250331235242

The apparent rupture duration time for P-waves:
Rupturing directly toward us:
Rupturing directly away from us:

Doppler effect

General apparent rupture duration

The apparent rupture duration for a seismic phase with local horizontal phase velocity at the observing station as

where is the station azimuth relative to the rupture direction.

The changes in as a function of receiver location are termed directivity effects.

Rupture Length

Since and
The rupture length is

The true rupture duration is

The average rupture velocity is

Example: 2004 Sumatra Earthquake Directivity

20250401002059
High-frequency (2–4 Hz) envelopes from teleseismic P-wave observations of the 2004 Sumatra earthquake.

Example: 2004 Sumatra Earthquake Directivity

20250401002059
We have s, ; and assume , so:

  • km
  • s
  • km/s

Rupture velocity

The rupture velocity is generally observed to be somewhat less than the shear-wave velocity for most earthquakes.
Anomalously fast ruptures sometimes exceed the local S-wave velocity and are termed supershear ruptures.

Earthquake rupture simulation


Eric Dunham

Earthquake faults

Earthquakes may be idealized as movement across a planar fault of arbitrary orientation

  • strike: , the azimuth of the fault from north where it intersects a horizontal surface
  • dip: , the angle from the horizontal
  • rake: , the angle between the slip vector and the strike

Earthquake Focal Mechanism

Moment tensor decomposition


Earthquake rupture

The Haskell source model

: the rise time, describes the slip duration at any point on a fault.
: the duration of the rupture ()

Earthquake Magnitude

How to quantify the size of an earthquake?

  • For historical reasons the most well-known measure of earthquake size is the earthquake magnitude.
  • Derived from the largest amplitude that is recorded on seismograms.
  • There are now many different types of magnitude scales, but all are connected in some way to the earliest definitions of magnitude.
20250402232519

Richter Magnitude (Local magnitude )

The original magnitude scale is based on the maximum amplitude recorded on a standard Wood-Anderson torsion seismograph.

: the amplitude of the reference event
: the epicentral distance

Richter Magnitude (Local magnitude )

An approximate empirical formula has been derived for at different ranges.
The local magnitude can be calculated by

where is the displacement amplitude in microns (10 m) and X is in kilometers.

  • Events below about are generally not felt
  • Significant damage to structures in California begins to occur at about
  • A earthquake implies amplitude 100 times greater than a event.

Global earthquakes: body wave magnitude

where A is the ground displacement in microns, T is the dominant period of the measured waves, is the epicentral distance in degrees, and Q is an empirical function of range and event depth h.

  • Why ?
  • h?

Global earthquakes: surface wave magnitude

For Rayleigh waves on vertical instruments:

Since the strongest Rayleigh wave arrivals are generally at a period of 20 s, this expression is often written as

  • Note that this equation is applicable only to shallow events
  • surface wave amplitudes are greatly reduced for deep events.

Magnitude saturation

The Haskell fault model

20250331235242

20250401002831

Source spectra

A boxcar pulse in the time domain produces a sinc function in the frequency domain:

Source spectra

The far-field amplitude spectrum for the Haskell fault model may be expressed as

where is a scaling term that includes geometrical spreading, etc

where .

Source spectra

We can approximate for and for

The amplitude spectrum for the Haskell fault model.

Magnitude saturation

Moment magnitude

The saturation of the and scales for large events helped motivate development of the moment magnitude

where is the moment measured in N-m.

  • The advantage of the scale is that it is clearly related to a physical property of the source and it does not saturate for even the largest earthquakes.
  • One unit increase in corresponds to a times increase in the moment.
  • A earthquake releases about 1000 times more energy than a event.

Magnitude as a function of moment

20250402234333

USGS Magnitude Types; Latest earthquake

20250402235133
20250402235213

The intensity scale

The local strength of ground shaking as determined by damage to structures and the perceptions of people who experienced the earthquake.

  • One earthquake can have different intensities at different locations.

USGS Latest Earthquakes

The Mercalli intensity scale (MMI)

20250403001052

JMA intensity scale

The Earthquake Cycle

  • Elastic rebound

Spring-block model (notebook)

20250407155433

spring slider

Spring-block model

Earthquake recurrence model

Stress drop

Moment

Stress drop

large D small A v.s. small D large A?

20250409232116

Circular crack model

20250409232316

Moment

Stress drop

Far-field seismic pulse

Area under displacement pulse is related to seismic moment

Radiation patterns, e.g.:

Far-field seismic pulse shape

20250409232944

Self-Similar Earthquake Scaling

Assuming dimensions are scaled proportionally, displacement D will increase by b

Self-Similar Earthquake Scaling

20250409234921

Self-Similar Earthquake Scaling

20250403002414

Self-similar Earthquake Scaling

Are different sized events self-similar?

  • Self-similarity predicts smaller events are “shifted” versions of larger events.

  • Following for self-similarity

Self-similar Earthquake Scaling

  • Stacking spectra following , spectral shapes are similar within uncertainties.
  • Self-similarity implies that apparent stress is size independent

Brune model (1970)

From model fit obtained:

  • Corner frequency ()
  • Radiated energy ()

Corner frequency

20250410000325

Far-field seismic pulse shape

20250409232623

Earthquake Energy Partitioning

20250409234107

Earthquake Energy Partitioning

20250403002307

Earthquake Scaling: static and dynamic

20250410000604

Stess and Radiated Energy

Is the stress drop close to the absolute stress level (weak faults) or just a small portion (strong faults)?

Radiation efficiency

For simple slip-weakening friction

where is called the scaled energy and is dimensionless.

Seismic efficiency

The radiation efficiency should not be confused with the seismic efficiency

The seismic efficiency is more difficult to estimate than the radiation efficiency because it depends upon the poorly constrained absolute stress level on the fault.

Apparent Stress (Abercrombie, 1995)

20250410001145

Apparent Stress (Ide and Beroza, 2001)

20250410000915

Seismic wave propagation (notebook)

P-wave S-wave

![height:600px](./assets/ElasticRebound.jpeg)

--- ![height:500px](./assets/Screenshot%202023-08-24%20at%2010.46.37.png) --- ![height:600px](./assets/Screenshot%202023-08-24%20at%2010.48.17.png)

---

## Earthquake recurrence model ![width:900px](./assets/Screenshot%202023-08-24%20at%2016.07.47.png) ( Shimazaki and Nakata, 1980)