590 likes | 1.14k Views
Differential wave equation and seismic events . Sean Ford & Holly Brown Berkeley Seismological Laboratory. Outline. Holly and Sean 101 Quick intro: Wave equation Monitoring nuclear tests Predicting ground motion for a future event on the Hayward Fault Hayward Fault tour. Introduction.
E N D
Differential wave equation and seismic events Sean Ford & Holly Brown Berkeley Seismological Laboratory
Outline • Holly and Sean 101 • Quick intro: Wave equation • Monitoring nuclear tests • Predicting ground motion for a future event on the Hayward Fault • Hayward Fault tour
Introduction • Seismic sources • Earthquake • Explosion Slip on a plane Pressure pulse on a sphere
Seismic wave equation • Start with Newton • Add ‘constitutive equation’ to relate stress to strain to displacement
Seismic wave equation • Can decompose to P and S wave solutions of the form • Wave equation can be solved by plane-wave
Finite differences • The wave equation is solved on a computer by using a discrete representation of the differential equation
Possible project: Sumatra earthquake Time (sec)
Possible project Time (sec)
Possible project T Time (sec)
Seismic sources • Earthquake - slip on a plane • Explosion - pressure pulse on a sphere Compressional (P-wave) radiation pattern Shear (S-wave) radiation pattern No volume change Compressional (P-wave) radiation constant No Shear (S-wave) radiation Volume change
North Korea Nuclear Test From Walter et al. (2008)
Seismic moment tensors M represents all possible force couple components due to a seismic source in a cartesian coordinate system Necessary to have two force couples (double couple, DC), so that angular momentum is conserved in the source sphere, which leads to Mij = Mji and the moment tensor is symmetric
Seismic moment tensors Model Source M Couples Focal Mechanism Double-couple (DC) x Strike-slip y z Compensated linear vector dipole (CLVD) x y Ring Fault z x Isotropic Explosion y z
Moment tensor inversion in matrix form d = Gm m = vector of 6 independent moment tensor elements m = (GTG)-1GTd
Source-type plot We calculate source-type plot parameters (Hudson et al., 1989)
Source-type plot We calculate source-type plot parameters (Hudson et al., 1989) Explosion Implosion
Source-type plot We calculate source-type plot parameters (Hudson et al., 1989) Explosion +CLVD DC -CLVD Implosion
Source-type plot We calculate source-type plot parameters (Hudson et al., 1989) Explosion HOYA +CLVD DC -CLVD Little Skull Implosion
Source-type plot We calculate source-type plot parameters (Hudson et al., 1989) Explosion HOYA HOYA +CLVD DC -CLVD Little Skull Little Skull Implosion
Source-type plot We calculate source-type plot parameters (Hudson et al., 1989) Explosion HOYA HOYA +CLVD DC -CLVD Little Skull Little Skull Implosion
Crandall Canyon, Utah Site of 6 Aug 08 Mine Collapse
Use method for nuclear test explosions to find source at Crandall Canyon
Crandall Canyon event plotted near where a closing crack or collapse source should plot Dominantly implosional Crandall Canyon
Crandall Canyon event plotted near where a closing crack or collapse source should plot Dominantly implosional But also some shear portion
Hayward fault • GoogleEarth tour • 1868 Hayward event M6.8 • Scenario events from USGS http://earthquake.usgs.gov/regional/nca/simulations/hayward/ • San Pablo epicenter • Fremont epicenter • View from Concord