Joining statistics and geophysics for assessment and uncertainty quantification of three-dimensional seismic Earth models

Carène Larmat, Monica Maceira, David M. Higdon, Dale N. Anderson

Research output: Contribution to journalArticlepeer-review

1 Scopus citations

Abstract

Seismic inversions produce seismic models, which are 3-dimensional (3D) images of wave velocity of the entire planet retrieved by fitting seismic measurements made on records of past earthquakes or other seismic events. Computing power of the TeraFlop era, along with the dataflow from new, very dense, seismic arrays, has led to a new generation of 3D seismic Earth models with an unprecedented level of resolution. Here we compare two recent models of western United States from the Dynamic North America (DNA) seismic imaging effort. The two models only differ in the wave propagation that was used for their inversion: one is based on ray theory (RT), and the other on finite frequency (FF). We evaluate the two models using an independent numerical method and statistical tests. We show that they differ in how they produce seismic signals from a subset of earthquakes that were used in the original inversion and were recorded on the US array. This is especially true for measurements done in the Yellowstone area which has a large negative seismic anomaly. This result is of importance for seismologists who have been debating on the practical benefit of using FF in ill-posed Earth inversions. Model evaluation, such as the one reported here, represents an opportunity for collaboration between geophysical and statistical communities. More opportunities should arise with the upcoming Exascale era, which will provide enough computational power to explore together several sources of errors in models with thousands of parameters, opening the way of uncertainty quantification of seismic models.

Original languageEnglish
Pages (from-to)277-289
Number of pages13
JournalStatistical Analysis and Data Mining
Volume10
Issue number5
DOIs
StatePublished - Oct 2017

Funding

We thank R.M. Allen, M. Obrebski, and R. Porritt for kindly providing us the DNA09 model as well as for fruitful discussions. We also thank IRIS for making the data so easily accessible to us (data from II, IU, TA, and US networks were used in this study), the developers of the Generic Mapping Tools software (Wessel and Smith, 1991) and Basemap. We especially thank the organizers of the CoDA conference for the opportunity of cross-disciplinary handshake. We also acknowledge the CIG program and the developers who make available SPECFEM3D. This research was supported by U.S. Department of Energy by Los Alamos National Laboratory under contract DE-AC52-06NA25396/LA12-SignalPropagation-NDD2Ab.

Keywords

  • permutation test
  • seismic models
  • seismology
  • spectral element method
  • t-test
  • uncertainty quantification

Fingerprint

Dive into the research topics of 'Joining statistics and geophysics for assessment and uncertainty quantification of three-dimensional seismic Earth models'. Together they form a unique fingerprint.

Cite this