Elastic waves propagate through solids by deforming the material elastically and transferring energy via strain-stress coupling. The wave equation derives from Newton's second law applied to continuous media, yielding plane wave solutions with velocity depending on elastic moduli and density. Seismic P and S waves are the two fundamental modes of elastic wave propagation in 3D solids.
From your study of the 1D wave equation and seismic waves, you know that disturbances can propagate through materials, and that seismic P- and S-waves have different velocities and particle motions. Elastic wave propagation in solids gives you the mathematical framework to understand *why* those differences exist — deriving wave speeds and wave modes from the fundamental mechanical properties of materials.
The key idea is treating the solid as a continuous elastic medium. Real solids are made of atoms separated by angstroms, but seismic wavelengths span meters to kilometers — ten orders of magnitude larger. At these scales, the discrete atomic structure is invisible, and the solid can be modeled as a continuous field of stress and strain. Newton's second law applied to an infinitesimal volume element gives: ρ ∂²u/∂t² = ∇·σ, where u is the displacement field, ρ is density, and σ is the stress tensor. The constitutive relation (generalized Hooke's law) then connects stress to strain: σ = C : ε, where C is the elastic stiffness tensor and ε is the strain tensor. Combining these two equations yields the elastic wave equation — a PDE governing how displacement disturbances evolve in space and time.
For an isotropic solid (one whose properties are the same in all directions), the stiffness tensor simplifies to just two independent parameters: the bulk modulus K (resistance to volumetric compression) and the shear modulus G (resistance to shear deformation). The wave equation then splits into two independent modes. Compressional (P-wave) motion involves volume changes — particles move back and forth along the direction of propagation — and travels at v_P = sqrt((K + 4G/3)/ρ). Shear (S-wave) motion involves no volume change — particles move perpendicular to the propagation direction — and travels at v_S = sqrt(G/ρ). Because K and G are both positive for any solid and G appears with a positive coefficient in v_P, P-waves are always faster than S-waves in isotropic media.
The dependence on G explains the S-wave behavior you encountered in seismology. For a liquid, G = 0 — liquids cannot resist sustained shear deformation because they flow. Substituting G = 0 gives v_S = 0: shear waves cannot exist in liquids. P-waves still propagate because liquids do resist compression (K > 0). This is the rigorous foundation for the seismological observation that S-waves disappear at the boundary with Earth's liquid outer core.
Plane wave solutions — displacement fields of the form u = A exp(i(k·x − ωt)) — are the natural solutions to the elastic wave equation. The dispersion relation (the relationship between wavenumber k and frequency ω) is non-dispersive for these bulk modes in a homogeneous solid: all frequencies travel at the same speed, which is why seismic body waves arrive as sharp pulses rather than smeared-out signals. Surface waves (Rayleigh and Love) behave differently — they are dispersive, with different frequencies traveling at different speeds — but that requires boundary conditions at a free surface and is the subject of the next topics in this sequence.