Skip to main content
U.S. flag

An official website of the United States government

Perturbational and nonperturbational inversion of Rayleigh-wave velocities

June 7, 2017

The inversion of Rayleigh-wave dispersion curves is a classic geophysical inverse problem. We have developed a set of MATLAB codes that performs forward modeling and inversion of Rayleigh-wave phase or group velocity measurements. We describe two different methods of inversion: a perturbational method based on finite elements and a nonperturbational method based on the recently developed Dix-type relation for Rayleigh waves. In practice, the nonperturbational method can be used to provide a good starting model that can be iteratively improved with the perturbational method. Although the perturbational method is well-known, we solve the forward problem using an eigenvalue/eigenvector solver instead of the conventional approach of root finding. Features of the codes include the ability to handle any mix of phase or group velocity measurements, combinations of modes of any order, the presence of a surface water layer, computation of partial derivatives due to changes in material properties and layer boundaries, and the implementation of an automatic grid of layers that is optimally suited for the depth sensitivity of Rayleigh waves.

Publication Year 2017
Title Perturbational and nonperturbational inversion of Rayleigh-wave velocities
DOI 10.1190/geo2016-0397.1
Authors Matthew M. Haney, Victor C. Tsai
Publication Type Article
Publication Subtype Journal Article
Series Title Geophysics
Index ID 70188355
Record Source USGS Publications Warehouse
USGS Organization Volcano Science Center