skip to main content
DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Accelerating full waveform inversion via source stacking and cross-correlations

Abstract

SUMMARY Accurate synthetic seismic wavefields can now be computed in 3-D earth models using the spectral element method (SEM), which helps improve resolution in full waveform global tomography. However, computational costs are still a challenge. These costs can be reduced by implementing a source stacking method, in which multiple earthquake sources are simultaneously triggered in only one teleseismic SEM simulation. One drawback of this approach is the perceived loss of resolution at depth, in particular because high-amplitude fundamental mode surface waves dominate the summed waveforms, without the possibility of windowing and weighting as in conventional waveform tomography. This can be addressed by redefining the cost-function and computing the cross-correlation wavefield between pairs of stations before each inversion iteration. While the Green’s function between the two stations is not reconstructed as well as in the case of ambient noise tomography, where sources are distributed more uniformly around the globe, this is not a drawback, since the same processing is applied to the 3-D synthetics and to the data, and the source parameters are known to a good approximation. By doing so, we can separate time windows with large energy arrivals corresponding to fundamental mode surface waves. This opens the possibility ofmore » designing a weighting scheme to bring out the contribution of overtones and body waves. It also makes it possible to balance the contributions of frequently sampled paths versus rarely sampled ones, as in more conventional tomography. Here we present the results of proof of concept testing of such an approach for a synthetic 3-component long period waveform data set (periods longer than 60 s), computed for 273 globally distributed events in a simple toy 3-D radially anisotropic upper mantle model which contains shear wave anomalies at different scales. We compare the results of inversion of 10 000 s long stacked time-series, starting from a 1-D model, using source stacked waveforms and station-pair cross-correlations of these stacked waveforms in the definition of the cost function. We compute the gradient and the Hessian using normal mode perturbation theory, which avoids the problem of cross-talk encountered when forming the gradient using an adjoint approach. We perform inversions with and without realistic noise added and show that the model can be recovered equally well using one or the other cost function. The proposed approach is computationally very efficient. While application to more realistic synthetic data sets is beyond the scope of this paper, as well as to real data, since that requires additional steps to account for such issues as missing data, we illustrate how this methodology can help inform first order questions such as model resolution in the presence of noise, and trade-offs between different physical parameters (anisotropy, attenuation, crustal structure, etc.) that would be computationally very costly to address adequately, when using conventional full waveform tomography based on single-event wavefield computations.« less

Authors:
 [1];  [2];  [3]
  1. University of California, Department of Earth and Planetary Science, Berkeley, CA, USA, Collège de France, Paris, France, Institut de Physique du Globe de Paris, France
  2. University of California, Department of Earth and Planetary Science, Berkeley, CA, USA
  3. University of California, Department of Earth and Planetary Science, Berkeley, CA, USA, Now at Google Inc., Mountainview, CA, USA
Publication Date:
Sponsoring Org.:
USDOE
OSTI Identifier:
1573084
Resource Type:
Published Article
Journal Name:
Geophysical Journal International
Additional Journal Information:
Journal Name: Geophysical Journal International Journal Volume: 220 Journal Issue: 1; Journal ID: ISSN 0956-540X
Publisher:
Oxford University Press
Country of Publication:
United Kingdom
Language:
English

Citation Formats

Romanowicz, Barbara, Chen, Li-Wei, and French, Scott W. Accelerating full waveform inversion via source stacking and cross-correlations. United Kingdom: N. p., 2019. Web. doi:10.1093/gji/ggz437.
Romanowicz, Barbara, Chen, Li-Wei, & French, Scott W. Accelerating full waveform inversion via source stacking and cross-correlations. United Kingdom. doi:10.1093/gji/ggz437.
Romanowicz, Barbara, Chen, Li-Wei, and French, Scott W. Mon . "Accelerating full waveform inversion via source stacking and cross-correlations". United Kingdom. doi:10.1093/gji/ggz437.
@article{osti_1573084,
title = {Accelerating full waveform inversion via source stacking and cross-correlations},
author = {Romanowicz, Barbara and Chen, Li-Wei and French, Scott W.},
abstractNote = {SUMMARY Accurate synthetic seismic wavefields can now be computed in 3-D earth models using the spectral element method (SEM), which helps improve resolution in full waveform global tomography. However, computational costs are still a challenge. These costs can be reduced by implementing a source stacking method, in which multiple earthquake sources are simultaneously triggered in only one teleseismic SEM simulation. One drawback of this approach is the perceived loss of resolution at depth, in particular because high-amplitude fundamental mode surface waves dominate the summed waveforms, without the possibility of windowing and weighting as in conventional waveform tomography. This can be addressed by redefining the cost-function and computing the cross-correlation wavefield between pairs of stations before each inversion iteration. While the Green’s function between the two stations is not reconstructed as well as in the case of ambient noise tomography, where sources are distributed more uniformly around the globe, this is not a drawback, since the same processing is applied to the 3-D synthetics and to the data, and the source parameters are known to a good approximation. By doing so, we can separate time windows with large energy arrivals corresponding to fundamental mode surface waves. This opens the possibility of designing a weighting scheme to bring out the contribution of overtones and body waves. It also makes it possible to balance the contributions of frequently sampled paths versus rarely sampled ones, as in more conventional tomography. Here we present the results of proof of concept testing of such an approach for a synthetic 3-component long period waveform data set (periods longer than 60 s), computed for 273 globally distributed events in a simple toy 3-D radially anisotropic upper mantle model which contains shear wave anomalies at different scales. We compare the results of inversion of 10 000 s long stacked time-series, starting from a 1-D model, using source stacked waveforms and station-pair cross-correlations of these stacked waveforms in the definition of the cost function. We compute the gradient and the Hessian using normal mode perturbation theory, which avoids the problem of cross-talk encountered when forming the gradient using an adjoint approach. We perform inversions with and without realistic noise added and show that the model can be recovered equally well using one or the other cost function. The proposed approach is computationally very efficient. While application to more realistic synthetic data sets is beyond the scope of this paper, as well as to real data, since that requires additional steps to account for such issues as missing data, we illustrate how this methodology can help inform first order questions such as model resolution in the presence of noise, and trade-offs between different physical parameters (anisotropy, attenuation, crustal structure, etc.) that would be computationally very costly to address adequately, when using conventional full waveform tomography based on single-event wavefield computations.},
doi = {10.1093/gji/ggz437},
journal = {Geophysical Journal International},
number = 1,
volume = 220,
place = {United Kingdom},
year = {2019},
month = {10}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record
DOI: 10.1093/gji/ggz437

Citation Metrics:
Cited by: 1 work
Citation information provided by
Web of Science

Save / Share:

Works referenced in this record:

Emergence of broadband Rayleigh waves from correlations of the ambient seismic noise: CORRELATIONS OF THE SEISMIC NOISE
journal, April 2004

  • Shapiro, N. M.; Campillo, M.
  • Geophysical Research Letters, Vol. 31, Issue 7
  • DOI: 10.1029/2004GL019491

Waveforms of Long-Period Body Waves In A Slightly Aspherical Earth Model
journal, January 1993


Earth's Background Free Oscillations
journal, March 1998


Simultaneous inversion for 3-D variations in shear and bulk velocity in the mantle
journal, March 1997


Imaging the North American continent using waveform inversion of global and USArray data
journal, September 2014


Earth's Correlation Wavefield: Late Coda Correlation
journal, April 2018

  • Phạm, Thanh‐Son; Tkalčić, Hrvoje; Sambridge, Malcolm
  • Geophysical Research Letters, Vol. 45, Issue 7
  • DOI: 10.1002/2018GL077244

A three-dimensional radially anisotropic model of shear velocity in the whole mantle
journal, October 2006


Seismic structure of the European upper mantle based on adjoint tomography
journal, February 2015

  • Zhu, Hejun; Bozdağ, Ebru; Tromp, Jeroen
  • Geophysical Journal International, Vol. 201, Issue 1
  • DOI: 10.1093/gji/ggu492

Three-dimensional full waveform inversion of short-period teleseismic wavefields based upon the SEM–DSM hybrid method
journal, May 2015

  • Monteiller, Vadim; Chevrot, Sébastien; Komatitsch, Dimitri
  • Geophysical Journal International, Vol. 202, Issue 2
  • DOI: 10.1093/gji/ggv189

Comparison of global waveform inversions with and without considering cross-branch modal coupling
journal, June 1995


An anisotropic shear velocity model of the Earth's mantle using normal modes, body waves, surface waves and long-period waveforms
journal, October 2014

  • Moulik, P.; Ekström, G.
  • Geophysical Journal International, Vol. 199, Issue 3
  • DOI: 10.1093/gji/ggu356

Theory of multisource crosstalk reduction by phase-encoded statics: Theory of multisource crosstalk reduction
journal, January 2011


Anisotropic shear-wave velocity structure of the Earth's mantle: A global model
journal, January 2008

  • Kustowski, B.; Ekström, G.; Dziewoński, A. M.
  • Journal of Geophysical Research, Vol. 113, Issue B6
  • DOI: 10.1029/2007JB005169

Partitioned waveform inversion and two-dimensional structure under the network of autonomously recording seismographs
journal, January 1990


On establishing the accuracy of noise tomography travel-time measurements in a realistic medium
journal, September 2009


The global CMT project 2004–2010: Centroid-moment tensors for 13,017 earthquakes
journal, June 2012

  • Ekström, G.; Nettles, M.; Dziewoński, A. M.
  • Physics of the Earth and Planetary Interiors, Vol. 200-201
  • DOI: 10.1016/j.pepi.2012.04.002

Fast computation of synthetic seismograms within a medium containing remote localized perturbations: a numerical solution to the scattering problem
journal, November 2016

  • Masson, Yder; Romanowicz, Barbara
  • Geophysical Journal International, Vol. 208, Issue 2
  • DOI: 10.1093/gji/ggw412

Spectral-element simulations of global seismic wave propagation-I. Validation
journal, May 2002


Radial anisotropy of the North American upper mantle based on adjoint tomography with USArray
journal, July 2017

  • Zhu, Hejun; Komatitsch, Dimitri; Tromp, Jeroen
  • Geophysical Journal International, Vol. 211, Issue 1
  • DOI: 10.1093/gji/ggx305

Earth's Background Free Oscillations
journal, May 2013


The three-dimensional shear velocity structure of the mantle from the inversion of body, surface and higher-mode waveforms
journal, December 2000


The effects of 3-D anelasticity (Q) structure on surface wave amplitudes: 3-D anelastic effects on surface wave amplitudes
journal, February 2012


Extracting seismic core phases with array interferometry: EXTRACT CORE PHASES WITH INTERFEROMETRY
journal, March 2013

  • Lin, Fan-Chi; Tsai, Victor C.; Schmandt, Brandon
  • Geophysical Research Letters, Vol. 40, Issue 6
  • DOI: 10.1002/grl.50237

Hessian kernels of seismic data functionals based upon adjoint techniques: Hessian kernels
journal, February 2011


Waveform Tomography Reveals Channeled Flow at the Base of the Oceanic Asthenosphere
journal, September 2013


Hybrid-domain simultaneous-source full waveform inversion without crosstalk noise
journal, September 2018

  • Zhang, Qingchen; Mao, Weijian; Zhou, Hui
  • Geophysical Journal International, Vol. 215, Issue 3
  • DOI: 10.1093/gji/ggy366

AxiSEM: broadband 3-D seismic wavefields in axisymmetric media
journal, January 2014

  • Nissen-Meyer, T.; van Driel, M.; Stähler, S. C.
  • Solid Earth, Vol. 5, Issue 1
  • DOI: 10.5194/se-5-425-2014

A global shear velocity model of the upper mantle from fundamental and higher Rayleigh mode measurements: GLOBAL SHEAR WAVE VELOCITY DISTRIBUTION
journal, October 2012

  • Debayle, Eric; Ricard, Yanick
  • Journal of Geophysical Research: Solid Earth, Vol. 117, Issue B10
  • DOI: 10.1029/2012JB009288

Global tomography using seismic hum
journal, December 2015

  • Haned, A.; Stutzmann, E.; Schimmel, M.
  • Geophysical Journal International, Vol. 204, Issue 2
  • DOI: 10.1093/gji/ggv516

The Iceland–Jan Mayen plume system and its impact on mantle dynamics in the North Atlantic region: Evidence from full-waveform inversion
journal, April 2013

  • Rickers, Florian; Fichtner, Andreas; Trampert, Jeannot
  • Earth and Planetary Science Letters, Vol. 367
  • DOI: 10.1016/j.epsl.2013.02.022

Determination of earthquake source parameters from waveform data for studies of global and regional seismicity
journal, April 1981

  • Dziewonski, A. M.; Chou, T. -A.; Woodhouse, J. H.
  • Journal of Geophysical Research: Solid Earth, Vol. 86, Issue B4
  • DOI: 10.1029/JB086iB04p02825

On the numerical implementation of time-reversal mirrors for tomographic imaging
journal, December 2013

  • Masson, Yder; Cupillard, Paul; Capdeville, Yann
  • Geophysical Journal International, Vol. 196, Issue 3
  • DOI: 10.1093/gji/ggt459

Lithospheric architecture of the South-Western Alps revealed by multiparameter teleseismic full-waveform inversion
journal, November 2017

  • Beller, S.; Monteiller, V.; Operto, S.
  • Geophysical Journal International, Vol. 212, Issue 2
  • DOI: 10.1093/gji/ggx216

Inversion of seismic reflection data in the acoustic approximation
journal, August 1984


An automated time-window selection algorithm for seismic tomography
journal, July 2009


Broad plumes rooted at the base of the Earth's mantle beneath major hotspots
journal, September 2015


Towards global earth tomography using the spectral element method: a technique based on source stacking
journal, August 2005


Introduction to the spectral element method for three-dimensional seismic wave propagation
journal, December 1999


Fast full-wavefield seismic inversion using encoded sources
journal, November 2009

  • Krebs, Jerome R.; Anderson, John E.; Hinkley, David
  • GEOPHYSICS, Vol. 74, Issue 6
  • DOI: 10.1190/1.3230502

Box tomography: localized imaging of remote targets buried in an unknown medium, a step forward for understanding key structures in the deep Earth
journal, April 2017

  • Masson, Yder; Romanowicz, Barbara
  • Geophysical Journal International, Vol. 211, Issue 1
  • DOI: 10.1093/gji/ggx141

Box Tomography: first application to the imaging of upper-mantle shear velocity and radial anisotropy structure beneath the North American continent
journal, February 2018

  • Clouzet, P.; Masson, Y.; Romanowicz, B.
  • Geophysical Journal International, Vol. 213, Issue 3
  • DOI: 10.1093/gji/ggy078

Confirmation of a change in the global shear velocity pattern at around 1000 km depth
journal, September 2017

  • Durand, S.; Debayle, E.; Ricard, Y.
  • Geophysical Journal International, Vol. 211, Issue 3
  • DOI: 10.1093/gji/ggx405

Methods for inversion of body-wave waveforms for localized three-dimensional seismic structure and an application to D′′ structure beneath Central America
journal, January 2014

  • Kawai, Kenji; Konishi, Kensuke; Geller, Robert J.
  • Geophysical Journal International, Vol. 197, Issue 1
  • DOI: 10.1093/gji/ggt520

On the computation of long period seismograms in a 3-D earth using normal mode based approximations
journal, November 2008


Whole-mantle radially anisotropic shear velocity structure from spectral-element waveform tomography
journal, September 2014

  • French, S. W.; Romanowicz, B. A.
  • Geophysical Journal International, Vol. 199, Issue 3
  • DOI: 10.1093/gji/ggu334

Global adjoint tomography: first-generation model
journal, September 2016

  • Bozdağ, Ebru; Peter, Daniel; Lefebvre, Matthieu
  • Geophysical Journal International, Vol. 207, Issue 3
  • DOI: 10.1093/gji/ggw356

On the robustness of global radially anisotropic surface wave tomography
journal, January 2010

  • Ferreira, A. M. G.; Woodhouse, J. H.; Visser, K.
  • Journal of Geophysical Research, Vol. 115, Issue B4
  • DOI: 10.1029/2009JB006716

Mapping the upper mantle: Three-dimensional modeling of earth structure by inversion of seismic waveforms
journal, July 1984

  • Woodhouse, John H.; Dziewonski, Adam M.
  • Journal of Geophysical Research: Solid Earth, Vol. 89, Issue B7
  • DOI: 10.1029/JB089iB07p05953

Coupling the spectral element method with a modal solution for elastic wave propagation in global earth models
journal, January 2003


Preliminary reference Earth model
journal, June 1981


Inferring upper-mantle structure by full waveform tomography with the spectral element method: SEM based global tomography
journal, March 2011


Global shear speed structure of the upper mantle and transition zone
journal, April 2013

  • Schaeffer, A. J.; Lebedev, S.
  • Geophysical Journal International, Vol. 194, Issue 1
  • DOI: 10.1093/gji/ggt095

Free oscillations and surface waves of an aspherical Earth
journal, December 1986


The deep roots of the western Pyrenees revealed by full waveform inversion of teleseismic P waves
journal, May 2016

  • Wang, Yi; Chevrot, Sébastien; Monteiller, Vadim
  • Geology, Vol. 44, Issue 6
  • DOI: 10.1130/G37812.1