The University of Southampton
University of Southampton Institutional Repository

Fast calculation of spatial sensitivity kernels for scattered waves in arbitrary heterogeneous media using graph theory

Fast calculation of spatial sensitivity kernels for scattered waves in arbitrary heterogeneous media using graph theory
Fast calculation of spatial sensitivity kernels for scattered waves in arbitrary heterogeneous media using graph theory
P-to-S and S-to-P receiver functions are widely used to constrain the seismic discontinuity structures of the Earth. Typically, receiver functions are projected to the depth and location of conversion assuming a 1-D layered Earth structure. Receiver function finite frequency sensitivity kernels have the potential to increase resolution. Here we present a method for rapidly calculating the P- and S-wave receiver function sensitivity kernels, based upon the shortest path method and Dijkstra's algorithm to calculate the traveltime fields, and accounting for geometrical spreading in heterogeneous media. The validity of the approach is evaluated by comparing with amplitudes derived from a finite difference elastic full waveform simulation in a complex subduction zone geometry. We show P-to-S and S-to-P kernels calculated using our method for three examples cases: a half space, a layer with topography, and a sinusoidal discontinuity. We also demonstrate the kernel recovery of discontinuities with these topographies by inverting synthetic data from SPECFEM2D. We find that P-to-S kernels recover the structure of strong topography better than S-to-P kernels, although S-to-P kernels may be useful in some situations. P-to-S kernels also show better recovery of the amplitude of the discontinuities in comparison to S-to-P, although both typically achieve values within a few percent of the input model. The computational cost of our approach for improved kernel calculation in heterogeneous media is up to a few tens of seconds per station for typical regional scale models on the scale of several 100s of kilometres.
Body waves, Numerical approximations and analysis, Wave scattering and diffraction
0956-540X
654-672
Bogiatzis, Petros
8fc5767f-51a2-4d3f-aab9-1ee9cfa9272d
Rychert, Catherine
70cf1e3a-58ea-455a-918a-1d570c5e53c5
Harmon, Nicholas
10d11a16-b8b0-4132-9354-652e72d8e830
Xie, Yujiang
77c46c7b-1aa6-4534-bca1-8c6a3dd40705
Bogiatzis, Petros
8fc5767f-51a2-4d3f-aab9-1ee9cfa9272d
Rychert, Catherine
70cf1e3a-58ea-455a-918a-1d570c5e53c5
Harmon, Nicholas
10d11a16-b8b0-4132-9354-652e72d8e830
Xie, Yujiang
77c46c7b-1aa6-4534-bca1-8c6a3dd40705

Bogiatzis, Petros, Rychert, Catherine, Harmon, Nicholas and Xie, Yujiang (2022) Fast calculation of spatial sensitivity kernels for scattered waves in arbitrary heterogeneous media using graph theory. Geophysical Journal International, 230 (1), 654-672. (doi:10.1093/gji/ggac078).

Record type: Article

Abstract

P-to-S and S-to-P receiver functions are widely used to constrain the seismic discontinuity structures of the Earth. Typically, receiver functions are projected to the depth and location of conversion assuming a 1-D layered Earth structure. Receiver function finite frequency sensitivity kernels have the potential to increase resolution. Here we present a method for rapidly calculating the P- and S-wave receiver function sensitivity kernels, based upon the shortest path method and Dijkstra's algorithm to calculate the traveltime fields, and accounting for geometrical spreading in heterogeneous media. The validity of the approach is evaluated by comparing with amplitudes derived from a finite difference elastic full waveform simulation in a complex subduction zone geometry. We show P-to-S and S-to-P kernels calculated using our method for three examples cases: a half space, a layer with topography, and a sinusoidal discontinuity. We also demonstrate the kernel recovery of discontinuities with these topographies by inverting synthetic data from SPECFEM2D. We find that P-to-S kernels recover the structure of strong topography better than S-to-P kernels, although S-to-P kernels may be useful in some situations. P-to-S kernels also show better recovery of the amplitude of the discontinuities in comparison to S-to-P, although both typically achieve values within a few percent of the input model. The computational cost of our approach for improved kernel calculation in heterogeneous media is up to a few tens of seconds per station for typical regional scale models on the scale of several 100s of kilometres.

Text
ggac078 - Version of Record
Available under License Creative Commons Attribution.
Download (9MB)

More information

Accepted/In Press date: 22 February 2022
Published date: 24 February 2022
Keywords: Body waves, Numerical approximations and analysis, Wave scattering and diffraction

Identifiers

Local EPrints ID: 456611
URI: http://eprints.soton.ac.uk/id/eprint/456611
ISSN: 0956-540X
PURE UUID: dea721fb-6827-4357-8875-041035d2c9cf
ORCID for Petros Bogiatzis: ORCID iD orcid.org/0000-0003-1902-7476
ORCID for Nicholas Harmon: ORCID iD orcid.org/0000-0002-0731-768X

Catalogue record

Date deposited: 05 May 2022 16:52
Last modified: 06 May 2022 01:50

Export record

Altmetrics

Contributors

Author: Petros Bogiatzis ORCID iD
Author: Nicholas Harmon ORCID iD
Author: Yujiang Xie

Download statistics

Downloads from ePrints over the past year. Other digital versions may also be available to download e.g. from the publisher's website.

View more statistics

Atom RSS 1.0 RSS 2.0

Contact ePrints Soton: eprints@soton.ac.uk

ePrints Soton supports OAI 2.0 with a base URL of http://eprints.soton.ac.uk/cgi/oai2

This repository has been built using EPrints software, developed at the University of Southampton, but available to everyone to use.

We use cookies to ensure that we give you the best experience on our website. If you continue without changing your settings, we will assume that you are happy to receive cookies on the University of Southampton website.

×