A2.16

Theoretical Geophysics
Part II, 2002

(i) State the equations that relate strain to displacement and stress to strain in a linear, isotropic elastic solid.

In the absence of body forces, the Euler equation for infinitesimal deformations of a solid of density ρ\rho is

ρ2uit2=σijxj\rho \frac{\partial^{2} u_{i}}{\partial t^{2}}=\frac{\partial \sigma_{i j}}{\partial x_{j}}

Derive an equation for u(x,t)\mathbf{u}(\mathbf{x}, t) in a linear, isotropic, homogeneous elastic solid. Hence show that both the dilatation θ=u\theta=\boldsymbol{\nabla} \cdot \mathbf{u} and the rotation ω=u\boldsymbol{\omega}=\nabla \wedge \mathbf{u} satisfy wave equations and find the corresponding wave speeds α\alpha and β\beta.

(ii) The ray parameter p=rsini/vp=r \sin i / v is constant along seismic rays in a spherically symmetric Earth, where v(r)v(r) is the relevant wave speed (α(\alpha or β)\beta) and i(r)i(r) is the angle between the ray and the local radial direction.

Express tani\tan i and sec ii in terms of pp and the variable η(r)=r/v\eta(r)=r / v. Hence show that the angular distance and travel time between a surface source and receiver, both at radius RR, are given by

Δ(p)=2rmRprdr(η2p2)1/2,T(p)=2rmRη2rdr(η2p2)1/2\Delta(p)=2 \int_{r_{m}}^{R} \frac{p}{r} \frac{d r}{\left(\eta^{2}-p^{2}\right)^{1 / 2}} \quad, \quad T(p)=2 \int_{r_{m}}^{R} \frac{\eta^{2}}{r} \frac{d r}{\left(\eta^{2}-p^{2}\right)^{1 / 2}}

where rmr_{m} is the minimum radius attained by the ray. What is η(rm)\eta\left(r_{m}\right) ?

A simple Earth model has a solid mantle in R/2<r<RR / 2<r<R and a liquid core in r<R/2r<R / 2. If α(r)=A/r\alpha(r)=A / r in the mantle, where AA is a constant, find Δ(p)\Delta(p) and T(p)T(p) for P\mathrm{P}-arrivals (direct paths lying entirely in the mantle), and show that

T=R2sinΔAT=\frac{R^{2} \sin \Delta}{A}

[You may assume that duuu1=2cos1(1u)\int \frac{d u}{u \sqrt{u-1}}=2 \cos ^{-1}\left(\frac{1}{\sqrt{u}}\right).]

Sketch the TΔT-\Delta curves for P\mathrm{P} and PcP\mathrm{PcP} arrivals on the same diagram and explain briefly why they terminate at Δ=cos114\Delta=\cos ^{-1} \frac{1}{4}.