Petascale Computing for Future Breakthroughs in Global Seismology

Will the advent of “petascale” computers be relevant to research in global seismic tomography? We illustrate here in detail two possible consequences of the expected leap in computing capability. First, being able to identify larger sets of differently regularized/parameterized solutions in shorter times will allow to evaluate their relative quality by more accurate statistical criteria than in the past. Second, it will become possible to compile large databases of sensitivity kernels, and update them efficiently in a non-linear inversion while iterating towards an optimal solution. We quantify the expected computational cost of the above endeavors, as a function of model resolution, and of the highest considered seismic-wave frequency.

[1]  R. Rosenfeld.  Nature , 2009, Otolaryngology--head and neck surgery : official journal of American Academy of Otolaryngology-Head and Neck Surgery.

[2]  C. Tape,et al.  Surface wave tomography: global membrane waves and adjoint methods , 2007 .

[3]  Carl Tape,et al.  Finite‐frequency tomography using adjoint methods—Methodology and examples using membrane surface waves , 2007 .

[4]  G. Soldati,et al.  Global seismic tomography and modern parallel computers , 2006 .

[5]  T. Jordan,et al.  Full 3D Waveform Tomography: A Comparison Between the Scattering-Integral and Adjoint- Wavefield Methods , 2006 .

[6]  J. Montagner,et al.  A SPICE Blind Test to Benchmark Global Tomographic Methods , 2006 .

[7]  Guust Nolet,et al.  A catalogue of deep mantle plumes: New results from finite‐frequency tomography , 2006 .

[8]  M. Sambridge,et al.  Trans-dimensional inverse problems, model comparison and the evidence , 2006 .

[9]  L. Boschi.  Global multiresolution models of surface wave propagation: comparing equivalently regularized Born and ray theoretical solutions , 2006 .

[10]  N. Sleep.  Mantle plumes from top to bottom , 2006 .

[11]  Lapo Boschi,et al.  On the relevance of Born theory in global seismic tomography , 2006 .

[12]  C. Tape,et al.  Towards Seismic Tomography Based Upon Adjoint Methods , 2005 .

[13]  Barbara Romanowicz,et al.  Towards global earth tomography using the spectral element method: a technique based on source stacking , 2005 .

[14]  E. Boschi,et al.  Linear or nonlinear rheology in the mantle: a 3D finite-element approach to postglacial rebound modeling , 2005 .

[15]  Qinya Liu,et al.  Tomography, Adjoint Methods, Time-Reversal, and Banana-Doughnut Kernels , 2004 .

[16]  E. Engdahl,et al.  Finite-Frequency Tomography Reveals a Variety of Plumes in the Mantle , 2004, Science.

[17]  C. Tape,et al.  Waves on a Spherical Membrane , 2003 .

[18]  A. Dziewoński,et al.  Distinct seismic anisotropy at the centre of the Earth , 2003 .

[19]  Chen Ji,et al.  A 14.6 billion degrees of freedom, 5 teraflops, 2.5 terabyte earthquake simulation on the Earth Simulator , 2003, ACM/IEEE SC 2003 Conference (SC'03).

[20]  David R. Anderson,et al.  Model selection and multimodel inference : a practical information-theoretic approach , 2003 .

[21]  D. Komatitsch,et al.  The Spectral-Element Method, Beowulf Computing, and Global Seismology , 2002, Science.

[22]  Lapo Boschi,et al.  A comparison of tomographic and geodynamic mantle models , 2002 .

[23]  A. Dziewoński,et al.  Global Event Location with Full and Sparse Data Sets Using Three-dimensional Models of Mantle P-wave Velocity , 2001 .

[24]  John S. J. Hsu,et al.  Bayesian Methods: An Analysis for Statisticians and Interdisciplinary Researchers , 1999 .

[25]  Harmen Bijwaard,et al.  Closing the gap between regional and global travel time tomography , 1998 .

[26]  S. Grand.  Mantle shear structure beneath the Americas and surrounding oceans , 1994 .

[27]  T. Tanimoto.  Modelling curved surface wave paths: membrane surface wave synthetics , 1990 .

[28]  Clifford M. Hurvich,et al.  Regression and time series model selection in small samples , 1989 .

[29]  Adam M. Dziewonski,et al.  Mapping the lower mantle: Determination of lateral heterogeneity in P velocity up to degree and order 6 , 1984 .

[30]  John H. Woodhouse,et al.  Mapping the upper mantle: Three‐dimensional modeling of earth structure by inversion of seismic waveforms , 1984 .

[31]  H. Akaike.  A new look at the statistical model identification , 1974 .