Abstract

Full-waveform inversion (FWI) is not yet a mature imaging technology for lithospheric imaging from teleseismic data. Therefore, its promise and pitfalls need to be assessed more accurately according to the specifications of teleseismic experiments. Three important issues are related to (1) the choice of the lithospheric parametrization for optimization and visualization, (2) the initial model and (3) the acquisition design, in particular in terms of receiver spread and sampling. These three issues are investigated with a realistic synthetic example inspired by the CIFALPS experiment in the Western Alps. Isotropic elastic FWI is implemented with an adjoint-state formalism and aims to update three parameter classes by minimization of a classical least-squares difference-based misfit function. Three different subsurface parametrizations, combining density (ρ) with P and S wave speeds (Vp and Vs) , P and S impedances (Ip and Is), or elastic moduli (λ and μ) are first discussed based on their radiation patterns before their assessment by FWI. We conclude that the (ρ, λ, μ) parametrization provides the FWI models that best correlate with the true ones after recombining a posteriori the (ρ, λ, μ) optimization parameters into Ip and Is. Owing to the low frequency content of teleseismic data, 1-D reference global models as PREM provide sufficiently accurate initial models for FWI after smoothing that is necessary to remove the imprint of the layering. Two kinds of station deployments are assessed: coarse areal geometry versus dense linear one. We unambiguously conclude that a coarse areal geometry should be favoured as it dramatically increases the penetration in depth of the imaging as well as the horizontal resolution. This results because the areal geometry significantly increases local wavenumber coverage, through a broader sampling of the scattering and dip angles, compared to a linear deployment.

Highlights

  • Building high-resolution and quantitative images of the lithosphere from body waves is a key challenge in earthquake seismology

  • This section reviews some basic principles of Full-waveform inversion (FWI) that will be useful to interpret the results of the numerical experiments carried out in this study. We review these principles with a matrix formalism suitable for the frequency-domain formulation of FWI (Pratt et al 1998), our teleseismic FWI is fully implemented in the time domain

  • We have analysed several theoretical and experimental factors that have a significant influence on teleseismic FWI for lithospheric imaging

Read more

Summary

Introduction

Building high-resolution and quantitative images of the lithosphere from body waves is a key challenge in earthquake seismology. From the geodynamical viewpoint, building images of the lithosphere with the best possible resolution according to the available frequency band (namely, of the order of the wavelength) is crucial to correlate tectonic deformation in the crust with deeper mantellic process arising in the asthenosphere. The teleseismic wavefield can fairly be approximated by an incident planar wavefield incoming from outside of the lithospheric target with different incidence angles and backazimuths. This planar configuration is conducive to a rather uniform illumination of the lithospheric target but might be a limiting resolution factor due to the limited angular illumination proved by plane-wave sources as opposed to point sources

Objectives
Methods
Conclusion
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call