Search USGS⌕ Search
Home/Archive/USGS

USGS · 70188355

Perturbational and nonperturbational inversion of Rayleigh-wave velocities

Matthew M. Haney·Victor C. Tsai
View original source ↗Export metadata ↓

Abstract

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.

Explore related subjects

GeophysicsArticle

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Matthew M. Haney, Victor C. Tsai. 2017. Perturbational and nonperturbational inversion of Rayleigh-wave velocities. https://doi.org/10.1190/geo2016-0397.1

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

USGS2023

Cooperative joint inversion of magnetotelluric and microseismic data for imaging the Geysers geothermal field, California, USA

The Geysers geothermal field located in northern California, USA, is the world’s largest electricity-generating geothermal facility. To delineate the spatio-temporal distribution of reservoir steam and recharge water, we have collected microseismic and magnetotelluric (MT) data using a dense array of stations in 2021. The microseismic and MT data have been inverted together using a 3D cooperative joint inversion workflow. The joint inversion exploits a cross-gradient structural constraint because electrical conductivity structures observed in the geothermal field are strongly correlated with

California↗
USGS2022

Inversion of induced polarization-affected towed-transient electromagnetic data in a lateritic regolith geology: A case study from western Tanzania

For several decades, induced polarization (IP) effects on transient electromagnetic (TEM) responses have been observed. These effects can manifest as late-time negative transients or as rapidly decaying curves and are usually associated with highly polarizable bodies. If neglected, IP effects can lead to erroneous resistivity models. Recent work allows IP effects to be incorporated into the inversion of TEM data on a more routine basis. In a recent field survey in western Tanzania, strongly IP-affected TEM signals are observed using a towed-transient electromagnetic (tTEM) system. The survey have been carried out to locate drinking water resources in a weathered regolith setting. In these settings, an inversion of tTEM data using a resistivity-only forward model (i.e., IP neglected) cannot fit the data and severely limits the value of the TEM data for hydrogeologic interpretation. To account for IP effects, we have applied a modified version of the Cole-Cole model called the maximum phase angle (MPA) model to invert IP-affected tTEM data. The MPA model incorporates four inversion model parameters: resistivity ( ⁠ ρ "> ρ ), MPA ( ⁠ ϕ max "> ϕ max ⁠ ), relaxation time ( ⁠ τ "> τ ), and frequency exponent ( ⁠ c "> c ). The MPA model fits the data well and improves the reliability of the resistivity model. In much of the surveyed region, the inverted models using MPA display a three-layer system consisting of an upper resistive laterite layer of varying thickness and an intermediate polarizable conductive unit overlying more resistive weathered basement rocks. The conductive polarizable layer is interpreted as a chemically weathered saprolite separating the surficial and deeper aquifers. Overall, tTEM inversion results provide a local understanding of groundwater systems, especially in such regions with very limited subsurface knowledge.

Kaguruka, Kitagata↗
USGS2022

Evidence-based guidelines for protective actions and earthquake early warning systems

Earthquake early warning systems (EEW) are becoming increasingly available or in development throughout the world. With public alerting in Mexico, Japan, Taiwan, and parts of the United States, it is important to provide evidence-based recommendations for protective action so people can protect themselves when they receive an alert. Best-practice warning communication research suggests that providing a protective action will increase the efficacy of the message. However, given the diversity of earthquakes and building types, as well as social and cultural contexts where these systems exist, the question is: what is the best protective action to recommend? The answer lies in maximizing life-saving protective actions during an earthquake event requires both contextually relevant messaging and widespread public education about appropriate protective actions under a range of conditions. By researching previous earthquake injury literature, examining current best practices and public education campaigns, key protective actions may be determined and used to increase the life-saving potential of earthquake early warning systems.

Geophysics↗