From Rayleigh to Park: A Century of Surface Wave Dispersion Mapping

Surface wave dispersion mapping represents a fundamental methodology in geophysics used to determine the mechanical properties of the shallow subsurface. By analyzing how the velocity of seismic waves varies with frequency, researchers can infer the shear-wave velocity profiles of soil and rock without invasive drilling. This process relies on the inherent property of surface waves—primarily Rayleigh and Love waves—to propagate through heterogeneous media at depths proportional to their wavelengths.

The evolution of this discipline began in the late 19th century with purely theoretical mathematical derivations and culminated in the late 20th century with the development of sophisticated digital signal processing techniques. Modern applications now bridge the gap between theoretical seismology and civil engineering, facilitating non-destructive testing for infrastructure integrity and the detection of buried hazards. The field is currently characterized by the use of multichannel arrays and automated inversion algorithms that transform raw ground motion data into high-resolution subsurface models.

Timeline

• 1885:Lord Rayleigh publishesOn Waves Propagated along the Plane Surface of an Elastic Solid, mathematically predicting the existence of surface waves that bear his name.
• 1911:A.E.H. Love publishesSome Problems of Geodynamics, describing transverse surface waves (Love waves) that occur in layered media.
• 1950s:The development of the Haskell-Thomson transfer matrix method provides a mathematical framework for calculating dispersion in multilayered media.
• 1980s:Introduction of Spectral Analysis of Surface Waves (SASW) using two-receiver configurations for geotechnical site characterization.
• 1999:Choon Park, Richard Miller, and Jianghai Xia introduce Multichannel Analysis of Surface Waves (MASW), significantly improving signal-to-noise ratios and data reliability.
• 2010s–Present:Integration of Full Waveform Inversion (FWI) and ambient noise tomography for urban seismic hazard mapping.

Background

Surface waves are a category of seismic waves that travel along the interface between two media, such as the Earth's surface and the atmosphere. Unlike body waves (P-waves and S-waves) which travel through the interior of a medium, surface waves are confined to the near-surface zone. Their energy decays exponentially with depth, meaning their influence is most pronounced within one wavelength of the boundary.

The primary characteristic that makes surface waves useful for subsurface imaging is dispersion. In a vertically heterogeneous medium, where physical properties like density and elasticity change with depth, surface waves of different frequencies travel at different phase velocities. High-frequency waves have shorter wavelengths and are confined to shallower layers, while low-frequency waves have longer wavelengths and penetrate deeper. By measuring the velocity of these different frequencies, geophysicists can construct a dispersion curve, which serves as the basis for calculating the depth-dependent properties of the ground.

The Rayleigh and Love Wave Distinction

Surface waves are generally categorized into two types based on their particle motion. Rayleigh waves, discovered by Lord Rayleigh in 1885, involve an elliptical motion in a vertical plane parallel to the direction of wave propagation. These waves can exist in a simple elastic half-space. Love waves, identified by A.E.H. Love in 1911, involve horizontal motion perpendicular to the direction of propagation. Love waves require a layered structure—specifically a lower-velocity layer overlying a higher-velocity layer—to exist. In practical site characterization, Rayleigh waves are more frequently utilized because they are easily generated by vertical impacts, such as a sledgehammer strike or a weight drop.

The Mathematical Foundations: 1885–1911

The mathematical process of surface wave analysis began when John William Strutt, better known as Lord Rayleigh, examined the equations of motion for an isotropic elastic solid. He sought to understand how waves would behave at a free surface. His 1885 paper demonstrated that a specific type of wave could propagate along the surface with a velocity slightly lower than that of the shear wave. Rayleigh noted that these waves would be the most prominent feature on a seismogram due to their two-dimensional spreading, which results in less geometric attenuation than three-dimensional body waves.

However, Rayleigh’s model assumed a homogeneous half-space, which does not exhibit dispersion; in his model, all frequencies traveled at the same speed. It was not until A.E.H. Love addressed the problem of a thin layer atop a half-space that the dispersive nature of surface waves was mathematically formalized. Love’s work showed that the interaction between layers of different velocities caused the wave speed to depend on frequency. This realization laid the groundwork for modern lithological characterization, as it implied that the velocity-frequency relationship (the dispersion curve) contained the structural information of the subsurface.

From Manual Analysis to Digital Processing

For much of the early 20th century, the analysis of surface waves was a labor-intensive manual process confined to global seismology and the study of the Earth's crustal thickness. Researchers used paper seismograms and hand-calculated arrival times to estimate phase velocities. The Society of Exploration Geophysicists (SEG) archives contain numerous records from this era, showing how analysts would manually pick peaks and troughs of waveforms to determine frequency content. This manual approach was limited by the accuracy of the timing instruments and the complexity of interfering wave modes.

The Rise of Active Source Methods

In the 1980s, the development of the Spectral Analysis of Surface Waves (SASW) method transitioned these concepts to the engineering scale. SASW used a source and two receivers to measure the phase difference between signals. While effective for simple sites, SASW struggled with ambient noise and the presence of "higher modes"—harmonics of surface waves that travel at higher velocities. Because the two-receiver method could not easily distinguish between these modes, the resulting subsurface profiles often contained significant errors.

The 1999 Milestone: Multichannel Analysis

A major shift occurred in 1999 with the publication of "Multichannel analysis of surface waves" by Park, Miller, and Xia in the journalGeophysics. This paper introduced the Multichannel Analysis of Surface Waves (MASW) technique, which utilized 12, 24, or more receivers simultaneously. By using a multichannel array, geophysicists could apply wavefield transformation techniques (such as the f-k transform or the phase-shift method) to separate the signal from noise in the frequency-wave-number domain.

"The multichannel approach allows for the recognition and isolation of various types of seismic events... This makes the dispersion imaging more strong against noise and non-fundamental mode interference."

The introduction of MASW allowed for the rapid collection of data in the field. Instead of moving two receivers multiple times, a linear array of geophones could be deployed once, and a single seismic shot could capture many frequencies. This efficiency led to the widespread adoption of surface wave mapping in geotechnical engineering for measuring $V_s30$—the average shear-wave velocity in the upper 30 meters—which is a critical parameter for seismic building codes.

Contemporary Inversion and Automated Software

Today, the transition from raw data to a subsurface model is largely handled by automated inversion software. The inversion process is a mathematical optimization where the computer generates a theoretical model of the Earth, calculates its predicted dispersion curve, and compares it to the observed data. The software then iteratively adjusts the model’s elastic moduli, density, and thickness of layers until the predicted and observed curves match.

Applications in Infrastructure and Anomaly Detection

The practical application of surface wave characteristics extends beyond soil layering. In civil engineering, the discipline is used for the non-destructive testing of concrete and asphalt. By analyzing the dispersion of induced high-frequency waves, engineers can detect delamination in bridge decks or voids under pavement. Because surface waves are sensitive to the stiffness of the material, they are uniquely suited to identifying the softening of foundations or the presence of buried utilities and tunnels.

Furthermore, microtremor analysis—using ambient city noise rather than a controlled strike—has allowed for subsurface imaging in urban environments where traditional seismic sources are prohibited. This reliance on the "background hum" of the Earth to map lithology represents the most recent evolution of the century-old theories established by Rayleigh and Love.

Practical Challenges in Heterogeneous Media

Despite the advancements in software, the study of surface waves in highly heterogeneous solid-state media remains complex. Lateral variations in geology can cause scattering and diffraction of the wavefield, complicating the 1D assumption used in most inversion algorithms. Research at facilities like the Surface Wave Hub continues to focus on the precise calibration of accelerometers and the development of 2D and 3D inversion schemes to account for these complexities. These efforts ensure that the interpretation of ground-motion signatures remains accurate, even in the presence of complex stratigraphy and engineered material interfaces.