Strings and 1D Wave Equation: Important Concepts/Assumptions
Strings and 1D Wave Equation: Important Concepts/Assumptions
Important Concepts/Assumptions:
(1) Taylor expansion
a
a+h
x
(2) Continuum: continuous distribution of material in space.
(3) Lagrangian description: follows a particle and feel the
action, perfect for seismic instrument that senses particle 1
motion!
Shear Wave
Compressional Wave
2
Simple Proof of Wave Equation :
An idealized mathematical string extends in the x
direction, initially straight with tension force
along it.
Initial condition: u(x, t=0) =0 in y direction, where
u= displacement. Now plucking the string.
Problem: Describe the displacement u(x, t)
Solution:
y
Condition: u is vertical
displacement, dx is small
u(x+dx, t)
2
Vertical tension along y:
sin 2 sin 1 du
u(x, t)
Newton II:
1
dx
2u
F = ma = dx 2
t
2 u
x
sin 2 sin 1 = dx 2
t
Taylor expansion about 0:
sin()=-3/3!+5/5!-~ (coef = 0 1 0 -1 )
tan()= +3/3+25/15 ~
Hence sin=tan=
Another way to derive sin(x)=tan(x)
sin x cos x
limx 0 = limx 0 = 1 (used LHopitals rule)
x 1 2 2 3
tan x sec x 1/cos x
limx o = limx 0 = limx o =1
x 1 1
T = 2 /
= 2 / k
v = /k = /T
Body waves
Needle
Suppose there is 3 cm displacement in east-west component,
what does that really mean?
Answer: That means a particle positioned at the needle tip
moves by 3 cm, i.e., particles on the ground/rock moves
by 3 cm due to seismic wave. Seismometer is a
Lagrangian motion recorder (meaning observing motion
by sitting on a particle that moves).
Map view 2 (rotated)
Map view 1
Radial (R)
N-S
E-W
V
V
Transverse (T)
Projection of ray, say P
Projection of ray, say P
on Earths surface
on Earths surface
Cross-sectional view V
Extension of
after rotation
K
wave vector
Kv
P-SV plane
8
v1, 1
Kr
R
Typical Observations: E-W, N-S, Vertical ground motions
(not a natural description of particle motions)
(1)
Azimuth vs. Back-azimuth
N
station
N
Energy
from
earthquake
Energy to
----> earthquake azimuth
----> station azimuth
(Back-azimuth)
So does + = 360 ?? NO.
(2) Rotate by back-azimuth to natural particle
motion coordinates
N
L(radial)
E
Energy
from
S
9
South Pole
10
Seismic Observations
Rotated, a rough view
Body wave
Surface wave
Sep 26, 2000
(0:17:17.01 GMT time,
TUC tangential
lat=-17.5, lon = -178.0)
TUC radial
TUC vertical
11
Polarization/particle motion
What is it: Plotting the two recorded components of a given
seismic phase, say P, or SS, etc.
Why useful:
(1) V-R (vertical-radial) plot give first motion + arrival angle +
strength of vertical vs. radial components.
(2) R-T (radial vs. transverse) plot gives presence of anisotropy
(3) E-N (eastwest vs. northsouth) gives back-azimuth!
V-R R-T
12
(1) Notational confusion:
P wave, S wave --- related to particle
motions
PcP, ScS --- these are called Seismic
phases or seismic arrivals, they
belong to either P or S types of
particle motions
(2) Many S-type seismic phases come in preferred
polarizations,
S --- Strong on both Radial and
Transverse
SS --- Strong on Transverse
ScS --- Strong on Transverse
SP, PKP, or anything with P or K or I in
name --> SV type --> Radial
(3) P waves travel faster than S waves because
+ 2
Vp = Vs =
Where and are Lame parameters, is related to bulk
modulus of materials (coefficient on how easily a piece of
material can compressed) and is shear modulus that tells how
easily something can be torn or sheared.
> 0 and > 0 ----->
Vp > Vs for most materials
(4) Most Earth materials can be
approximated by what is called a
3
Poisson Solid where ~
-------> Vp ~
Vs
13
Distance Calculations in Global Earth Problems
Known:
Radius of Earth = 6371 km
Distance conversion to
degrees:
14
Advantage: Suppose we want to find the
perpendicular distance from the Smiley to the line,
by rotating to the Equatorial system, the latitude of
the rotated Smiley in the new Equatorial framework
IS the perpendicular distance we want.
= x 1 + x 2 + x 3
x1 x 2 x 3
Gradient of scalar field,
which makes sense to be
related to P.
P direction = wave
direction, compression
For comparison, divergence of a vector field ,
represents flux of vector field through a unit
volume:
1 2 3
= + + = i,i
x1 x 2 x 3