Lamb's problem

Lamb's problem. In the field of theoretical seismology and wave mechanics, it is a foundational boundary value problem concerning the generation of elastic waves in a half-space by a point source. First formulated by the British applied mathematician Horace Lamb in a seminal 1904 paper, it seeks the complete wavefield—encompassing both P-waves and S-waves—resulting from an impulsive force applied at or beneath the surface of a homogeneous, isotropic, elastic medium. Its solution provides critical insight into the nature of seismic waves, Rayleigh waves, and head waves, serving as a benchmark for validating numerical methods in geophysics and non-destructive testing.

Definition and historical context

The problem was precisely articulated by Horace Lamb in his publication in the Proceedings of the London Mathematical Society. Lamb, renowned for his work in fluid dynamics and author of the influential textbook Hydrodynamics, turned his analytical prowess to the propagation of disturbances in an elastic solid. The historical context lies in the late 19th and early 20th-century development of continuum mechanics, alongside contemporaneous work by figures like Lord Rayleigh on surface waves and Augustin-Jean Fresnel on wave theory. Its formulation was driven by fundamental questions in theoretical physics about how energy radiates from a localized source, predating the instrumental needs of modern exploration geophysics. The problem's elegance and complexity made it a classic touchstone, linking the mathematical theories of James Clerk Maxwell and George Gabriel Stokes to practical observations of earthquake-generated ground motion.

Mathematical formulation

The medium is modeled as a homogeneous, isotropic, linearly elastic half-space occupying the region z ≥ 0, with the free surface at z = 0. The governing equations are the Navier-Cauchy equations of motion, derived from Hooke's law and the strain-displacement relations. The source is typically an idealized point force, often vertical or horizontal, applied at the surface or at some depth, represented mathematically using the Dirac delta function in space and time. The Lamé parameters (λ and μ) define the material's elastic properties, from which the velocities of the P-wave and S-wave are derived. The boundary conditions require that the stress tensor components vanish at the free surface, while Sommerfeld radiation conditions are imposed to ensure outgoing waves at infinity, preventing energy from reflecting back from the mathematical boundaries of the model.

Analytical solutions

Lamb's own solution employed integral transform methods, specifically the application of the Laplace transform in time and a double Fourier transform in the horizontal spatial coordinates. The solution in the transform domain reveals distinct contributions corresponding to different wave types: the direct body waves, waves reflected from the free surface, and the critically refracted head wave (or P-S wave). The inverse transformation to obtain the physical wavefield in the time-space domain is highly non-trivial, involving the evaluation of complex integrals and leading to expressions containing Bessel functions and Heaviside step functions. Notable extensions were provided by later researchers like Leon Knopoff and Keiiti Aki, who derived explicit formulas for displacements at the free surface. The solution famously shows the arrival sequence of the P-wave, S-wave, and the dispersive Rayleigh wave, which decays slowly with distance.

Applications in geophysics

The problem serves as the fundamental Green's function for a half-space, making it indispensable in theoretical seismology. It is used to model and understand the complete seismogram from a near-surface source, such as an explosion or a small earthquake. Insights from its solution underpin the interpretation of seismic data in oil exploration, particularly for understanding wave amplitudes and phases in vertical seismic profiling and reflection seismology. It provides the theoretical basis for identifying phases like the Stoneley wave in borehole acoustics and for calibrating instruments used by organizations like the United States Geological Survey. Furthermore, it is critical in seismic hazard assessment for predicting ground motion from crustal earthquakes and in the design of structures to withstand vibrations.

Modern computational approaches

While the analytical solution provides profound insight, its complexity limits its use for heterogeneous or anisotropic media. Modern approaches rely heavily on high-performance computational methods. Techniques such as the finite difference method, spectral element method, and boundary element method are now routinely employed to solve generalized versions of Lamb's problem for complex subsurface models. These numerical solutions are validated against the classical analytical result. Institutions like the California Institute of Technology and Stanford University have developed advanced codes for seismic wave propagation simulation. Furthermore, the problem is a standard benchmark in software packages like SPECFEM3D and is used in research on inverse problems and full-waveform inversion conducted by companies such as Schlumberger and CGG (company).

Category:Elasticity (physics) Category:Seismology Category:Partial differential equations Category:Wave mechanics