Past exam of the physics of the Earth as a Planet course of the University of Cambridge / 2016 / Paper 1 / 3

For uniform density $\rho >0$ and stiffness, linear momentum balance without body force is $\rho {\ddot{u}}_i={\partial }_j{\sigma }_{ij}$. The minor elastic symmetries make the two symmetric-strain terms equal, so ${\partial }_j{\sigma }_{ij}=c_{ijkl}{u}_{k,lj}$. Write $s=ct-{\hat{p}}_j{x}_j$ and $u_i=a_{i}f(s)$. Then ${\ddot{u}}_i=c^2{a}_i{f}^{\prime \prime }$ and $u_{k,lj}=a_k{\hat{p}}_l{\hat{p}}_j{f}^{\prime \prime }$. Thus a nontrivial plane-wave profile satisfiesThe density-normalized elastic acoustic tensor is consequentlyThe alternative dummy-index form $B_{ik}={\rho }^{-1}{c}_{ijkl}{\hat{p}}_j{\hat{p}}_l$ is identical. In ordinary conservative elasticity, major symmetry $c_{ikjl}=c_{jlik}$ gives $B_{ji}=B_{ij}$ after interchanging the dummy indices. This symmetry comes from a strain-energy density; it is an assumption of the ordinary elastic tensor, not a consequence of homogeneity alone.

For any such eigenpair the profile $f$ is arbitrary, since both sides carry the same $f^{\prime \prime }$. This nondispersive result uses linear, local, time-independent elasticity without attenuation or spatial dispersion. A profile with $f^{\prime \prime }=0$ is a trivial uniform or affine displacement field and cannot be used by itself to infer the wave-speed eigenvalue condition.

A real symmetric $3\times 3$ matrix has three real eigenvalues and a complete orthonormal basis of eigenvectors. For a stable elastic material, strong ellipticity makes ${\mathbf{a}}^{\mathsf{T}}B\mathbf{a}>0$ for every nonzero polarization, so its eigenvalues are positive squared wave speeds. There are therefore three orthogonal plane-wave polarizations in a given direction, with three potentially distinct speeds in an anisotropic medium. Degenerate eigenvalues allow arbitrary orthonormal combinations within their eigenspace; they do not remove polarization degrees of freedom. General anisotropic modes need not be exactly longitudinal or transverse.

For an isotropic medium let $\lambda =\kappa -2\mu /3$. Contract the isotropic stiffness with ${\hat{p}}_k{\hat{p}}_l$:UsinggivesFor $\mathbf{a}\parallel \hat{\mathbf{p}}$ the longitudinal projector has eigenvalue one, so $B\mathbf{a}={\alpha }^2\mathbf{a}$: this is the P wave. For either of two independent $\mathbf{a}\perp \hat{\mathbf{p}}$ the projector vanishes, so $B\mathbf{a}={\beta }^2\mathbf{a}$: these are two degenerate transverse S wave polarizations. An isotropic solid has one P speed and one S speed, but three polarization modes: one longitudinal and two transverse. A fluid limit with zero shear modulus has no propagating S wave, so the positive-speed conclusion requires a solid with positive shear stiffness.

Use geometric rays for body-wave arrivals and a traction-free upper surface. A compact isotropic explosion source has moment tensor proportional to ${\delta }_{ij}$ and produces direct P body radiation in a homogeneous isotropic medium. Thus the basic near-offset arrival is the direct upper-layer P wave, with time $X/{\alpha }_1$. A pure explosion is not a shear fault: it does not automatically emit a direct S wave of comparable strength. Interfaces can nevertheless convert its P energy into SV motion, and the free surface can support Rayleigh waves.

Downward P energy reflects from the interface to give a PP reflection returning to the receiver; for a horizontal interface its ray is symmetric and its time is $\sqrt{X^2+4H^2}/{\alpha }_1$. Other transmitted P energy enters the faster half-space. At $\sin i_c={\alpha }_1/{\alpha }_2$, a critical ray travels along the interface and feeds an upgoing seismic head wave, often labeled Pn in this simplified crustal model. Its ideal travel-time branch isIt requires enough offset for the two inclined upper-layer legs, $X\ge 2H\tan i_c$. Once $X$ exceeds the larger first-arrival crossover distance $X_c=2H\sqrt{({\alpha }_2+{\alpha }_1)/({\alpha }_2-{\alpha }_1)}$, this branch arrives before the direct upper-layer wave because it spends most of the path in faster material. A transmitted ray descending into an otherwise uniform infinite half-space does not return as a separate reflection without an additional reflector.

At oblique incidence, the bonded-interface conditions couple P and SV polarizations. A downgoing P ray can return as a converted SV reflection in the upper layer, with common horizontal slowness $p=\sin i_P/{\alpha }_1=\sin i_S/{\beta }_1$. There are transmitted P/SV branches and, after further interaction with the free surface or interface, mixed-mode reflections and layer reverberations. Conversion coefficients depend on density and impedance as well as wave speeds, and some vanish at normal incidence. The supplied inequalities do not order ${\beta }_2$ against ${\alpha }_1$, so one cannot universally draw a P-incident critical head ray travelling at the lower-layer S speed: that direct conversion construction requires ${\beta }_2>{\alpha }_1$. A shear head-wave phase can be present if the appropriate SV energy and ray conditions are available, but is not an inevitable direct explosion arrival from the stated data alone.

At later times a shallow explosion can excite dispersive Rayleigh-wave trains, with coupled P-SV motion decaying with depth. Short wavelengths mainly sample the slow upper layer, whereas long wavelengths penetrate the lower half-space; group and phase speeds need not be equal. Multiple reflected/converted body waves can overlap this coda. The shear-speed contrast also permits Love-wave modes structurally, but SH motion is decoupled from P-SV in this flat isotropic model. An ideal isotropic explosion does not excite Love waves here; appreciable SH/Love energy would require source asymmetry, scattering or a broken symmetry of the model.

Exact amplitudes, resolved arrival separation and attenuation cannot be calculated without densities, source spectrum, damping and instrument response. The qualitative robust features are a near-offset direct P arrival, reflected and converted phases, a faster-lower-layer P head wave that can become first at long range, reverberations and P-SV surface-wave energy. The ray diagram uses declared illustrative speed ratios solely to show the geometry; they are not additional supplied material constants.