Global seismic attenuation imaging using full-waveform inversion: a comparative assessment of different choices of misfit functionals

Journal Article · Sat Oct 14 00:00:00 EDT 2017 · Geophysical Journal International

DOI:https://doi.org/10.1093/gji/ggx442· OSTI ID:1523659

Karaoğlu, Haydar ^[1]; Romanowicz, Barbara ^[2]

---

1. Inst. de Physique du Globe de Paris, Paris (France)
2. Inst. de Physique du Globe de Paris, Paris (France); Univ. of California, Berkeley, CA (United States). Berkeley Seismological Lab.; College de France, Paris (France)

We present the results of synthetic tests that aim at evaluating the relative performance of three different definitions of misfit functionals in the context of 3-D imaging of shear wave attenuation in the earth’s upper mantle at the global scale, using long-period full-waveform data. The synthetic tests are conducted with simple hypothetical upper-mantle models that contain Q_μ anomalies centred at different depths and locations, with or without additional seismic velocity anomalies. To build synthetic waveform data sets, we performed simulations of 50 events in the hypothetical (target) models, using the spectral element method, filtered in the period range 60–400 s. The selected events are chosen among 273 events used in the development of radially anisotropic model SEMUCB-WM1 and recorded at 495 stations worldwide. The synthetic Z-component waveforms correspond to paths and time intervals (fundamental mode and overtone Rayleigh waves) that exist in the real waveform data set. The inversions for shear attenuation structure are carried out using a Gauss–Newton optimization scheme in which the gradient and Hessian are computed using normal mode perturbation theory. The three different misfit functionals considered are based on time domain waveform (WF) and waveform envelope (E-WF) differences, as well as spectral amplitude ratios (SA), between observed and predicted waveforms. We evaluate the performance of the three misfit functional definitions in the presence of seismic noise and unresolved S-wave velocity heterogeneity and discuss the relative importance of physical dispersion effects due to 3-D Q_μ structure. We observed that the performance of WF is poorer than the other two misfit functionals in recovering attenuation structure, unless anelastic dispersion effects are taken into account in the calculation of partial derivatives. WF also turns out to be more sensitive to seismic noise than E-WF and SA. Overall, SA performs best for attenuation imaging. Our tests show that it is important to account for 3-D elastic effects (focusing) before inverting for Q_μ. Additionally, we show that including high signal-to-noise ratio overtone wave packets is necessary to resolve Q_μ structure at depths greater than 250 km.

• Research Organization: Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States). National Energy Research Scientific Computing Center (NERSC)
• Sponsoring Organization: USDOE
• OSTI ID: 1523659
• Journal Information: Geophysical Journal International, Vol. 212, Issue 2; ISSN 0956-540X
• Publisher: Oxford University Press
• Country of Publication: United States
• Language: English

Related Subjects

58 GEOSCIENCES
Seismic attenuation
Seismic tomography