AxiSEM: broadband 3-D seismic wavefields in axisymmetric media

Abstract. We present a methodology to compute 3-D global seismic wavefields for realistic earthquake sources in visco-elastic anisotropic media, covering applications across the observable seismic frequency band with moderate computational resources. This is accommodated by mandating axisymmetric background models that allow for a multipole expansion such that only a 2-D computational domain is needed, whereas the azimuthal third dimension is computed analytically on the fly. This dimensional collapse opens doors for storing space–time wavefields on disk that can be used to compute Frechet sensitivity kernels for waveform tomography. We use the corresponding publicly available AxiSEM ( www.axisem.info ) open-source spectral-element code, demonstrate its excellent scalability on supercomputers, a diverse range of applications ranging from normal modes to small-scale lowermost mantle structures, tomographic models, and comparison with observed data, and discuss further avenues to pursue with this methodology.

[1]  Peter Moczo,et al.  Time-frequency misfit and goodness-of-fit criteria for quantitative comparison of time signals , 2009 .

[2]  Andreas Fichtner,et al.  Theoretical background for continental‐ and global‐scale full‐waveform inversion in the time–frequency domain , 2008 .

[3]  J. Tromp,et al.  Three-dimensional structure of the African superplume from waveform modelling , 2005 .

[4]  Laurent Gizon,et al.  Time-Distance Helioseismology: The Forward Problem for Random Distributed Sources , 2002 .

[5]  PenShu Yeh.  Implementation of CCSDS Lossless Data Compression for Space and Data Archive Applications , 2002 .

[6]  Robert J. Geller,et al.  Complete synthetic seismograms for 3-D heterogeneous Earth models computed using modified DSM operators and their applicability to inversion for Earth structure , 2000 .

[7]  Guust Nolet.  A Breviary of Seismic Tomography: Preface , 2008 .

[8]  E. D. Mercerat,et al.  On the linearity of cross-correlation delay times in finite-frequency tomography , 2012 .

[9]  J. Tromp,et al.  Elastic imaging and time-lapse migration based on adjoint methods , 2009 .

[10]  Karin Sigloch,et al.  ObsPyLoad: A Tool for Fully Automated Retrieval of Seismological Waveform Data , 2013 .

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

[12]  B. Kennett,et al.  Traveltimes for global earthquake location and phase identification , 1991 .

[13]  R. Geller,et al.  Finite‐frequency structural sensitivities of short‐period compressional body waves , 2012 .

[14]  T. Nissen‐Meyer,et al.  Triplicated P-wave measurements for waveform tomography of the mantle transition zone , 2012 .

[15]  Nicholas Rawlinson,et al.  Seismic tomography: a window into deep Earth , 2010 .

[16]  A. Dziewoński,et al.  Anisotropic shear‐wave velocity structure of the Earth's mantle: A global model , 2008 .

[17]  Steven G. Johnson,et al.  The Design and Implementation of FFTW3 , 2005, Proceedings of the IEEE.

[18]  F. Dahlen,et al.  Traveltimes of waves in three-dimensional random media , 2003 .

[19]  Alexandre Fournier,et al.  A 2‐D spectral‐element method for computing spherical‐earth seismograms—II. Waves in solid–fluid media , 2008 .

[20]  Peter Messmer,et al.  Forward and adjoint simulations of seismic wave propagation on emerging large-scale GPU architectures , 2012, 2012 International Conference for High Performance Computing, Networking, Storage and Analysis.

[21]  Emmanuel Chaljub,et al.  Spectral element modelling of three-dimensional wave propagation in a self-gravitating Earth with an arbitrarily stratified outer core , 2003, physics/0308102.

[22]  D. Schorlemmer,et al.  Testing earthquake source inversion methodologies , 2011 .

[23]  Russ Rew,et al.  NetCDF: an interface for scientific data access , 1990, IEEE Computer Graphics and Applications.

[24]  Alexandre Fournier,et al.  Spherical‐earth Fréchet sensitivity kernels , 2007 .

[25]  Hiroshi Takenaka,et al.  Accurate and efficient modeling of global seismic wave propagation for an attenuative Earth model including the center , 2012 .

[26]  B. Romanowicz,et al.  On the numerical implementation of time-reversal mirrors for tomographic imaging , 2013, Geophysical Journal International.

[27]  J. Tromp,et al.  Theoretical Global Seismology , 1998 .

[28]  Alexandre Fournier,et al.  A two‐dimensional spectral‐element method for computing spherical‐earth seismograms – I. Moment‐tensor source , 2007 .

[29]  Paul Friberg,et al.  Near real‐time simulations of global CMT earthquakes , 2010 .

[30]  Heiner Igel,et al.  P‐SV wave propagation in the Earth's mantle using finite differences: Application to heterogeneous lowermost mantle structure , 1996 .

[31]  G. Nolet.  A Breviary of Seismic Tomography: Frontmatter , 2008 .

[32]  D. Komatitsch,et al.  Spectral-element simulations of global seismic wave propagation—I. Validation , 2002 .

[33]  Daoyuan Sun,et al.  Seismological support for the metastable superplume model, sharp features, and phase changes within the lower mantle , 2007, Proceedings of the National Academy of Sciences.

[34]  Jean-Pierre Vilotte,et al.  A Fourier-spectral element algorithm for thermal convection in rotating axisymmetric containers , 2005 .

[35]  Matthew B. Kennel.  KDTREE 2: Fortran 95 and C++ software to efficiently search for near neighbors in a multi-dimensional Euclidean space , 2004 .

[36]  B. Kennett,et al.  ELLIPTICITY CORRECTIONS FOR SEISMIC PHASES , 1996 .

[37]  Jean-Pierre Vilotte,et al.  Solving elastodynamics in a fluid-solid heterogeneous sphere: a parallel spectral element approximation on non-conforming grids , 2003 .

[38]  Simon Stähler,et al.  Fully probabilistic seismic source inversion - Part 1: Efficient parameterisation , 2013 .

[39]  A. Fichtner,et al.  Efficient numerical surface wave propagation through the optimization of discrete crustal models—a technique based on non-linear dispersion curve matching (DCM) , 2008 .

[40]  Guust Nolet,et al.  Fréchet kernels for finite‐frequency traveltimes—II. Examples , 2000 .

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

[42]  T. Nissen‐Meyer,et al.  Seismic wave propagation in fully anisotropic axisymmetric media , 2014 .

[43]  J. Tromp,et al.  Misfit functions for full waveform inversion based on instantaneous phase and envelope measurements , 2011 .

[44]  Barbara Romanowicz,et al.  The three‐dimensional shear velocity structure of the mantle from the inversion of body, surface and higher‐mode waveforms , 2000 .

[45]  S. Rost.  Deep Earth: Core-mantle boundary landscapes , 2013 .

[46]  D. Komatitsch,et al.  Spectral-element simulations of global seismic wave propagation: II. Three-dimensional models, oceans, rotation and self-gravitation , 2002 .

[47]  T. Nissen‐Meyer,et al.  Wave propagation in 3D spherical sections: effects of subduction zones , 2002 .

[48]  Heiner Igel,et al.  Global SH-wave propagation using a parallel axisymmetric spherical finite-difference scheme: Application to whole mantle scattering , 2008 .

[49]  Martin van Driel,et al.  Optimized viscoelastic wave propagation for weakly dissipative media , 2014 .

[50]  C. Bassin,et al.  The Current Limits of resolution for surface wave tomography in North America , 2000 .

[51]  T. Nissen‐Meyer,et al.  Savani: A variable resolution whole‐mantle model of anisotropic shear velocity variations based on multiple data sets , 2014 .

[52]  Dapeng Zhao,et al.  Seismic structure and origin of hotspots and mantle plumes , 2001 .

[53]  F. Dahlen.  The Normal Modes of a Rotating, Elliptical Earth , 1968 .

[54]  Heiner Igel,et al.  Seismic imaging of the laterally varying D″ region beneath the Cocos Plate , 2007 .

[55]  T. Nissen‐Meyer,et al.  Seismic waveform inversion for core–mantle boundary topography , 2014 .

[56]  J. Montagner,et al.  Residual homogenization for seismic forward and inverse problems in layered media , 2013 .

[57]  Dimitri Komatitsch,et al.  A hybrid method to compute short-period synthetic seismograms of teleseismic body waves in a 3-D regional model , 2013 .

[58]  E. R. Engdahl,et al.  Constraints on seismic velocities in the Earth from traveltimes , 1995 .

[59]  Thomas J. Owens,et al.  The TauP Toolkit: Flexible Seismic Travel-Time and Raypath Utilities , 1999 .

[60]  Wolfgang Friederich,et al.  COMPLETE SYNTHETIC SEISMOGRAMS FOR A SPHERICALLY SYMMETRIC EARTH BY A NUMERICAL COMPUTATION OF THE GREEN'S FUNCTION IN THE FREQUENCY DOMAIN , 1995 .

[61]  John H. Woodhouse,et al.  S40RTS: A degree-40 shear-velocity model for the mantle from new Rayleigh wave dispersion, teleseismic traveltime and normal-mode splitting function measurements , 2011 .

[62]  John H. Day,et al.  Implementation of CCSDS Lossless Data Compression in HDF , 2002 .

[63]  F. Scherbaum,et al.  Acoustic simulation of P‐wave propagation in a heterogeneous spherical earth: numerical method and application to precursor waves to PKPdf , 2000 .

[64]  Tom Chen,et al.  Design and implementation , 2006, IEEE Commun. Mag..

[65]  Emanuele Casarotti,et al.  Forward and adjoint simulations of seismic wave propagation on fully unstructured hexahedral meshes , 2011 .

[66]  Jean-Pierre Vilotte,et al.  Spectral Element Analysis in Seismology , 2007 .

[67]  Philippe Roux,et al.  Teleseismic correlations of ambient seismic noise for deep global imaging of the Earth , 2013 .

[68]  Kenji Kawai,et al.  Complete synthetic seismograms up to 2 Hz for transversely isotropic spherically symmetric media , 2006 .

[69]  T. Nissen‐Meyer,et al.  Seismic waveform sensitivity to global boundary topography , 2012 .

[70]  Lion Krischer,et al.  ObsPy: A Python Toolbox for Seismology , 2010 .

[71]  Dimitri Komatitsch,et al.  High-resolution seismic array imaging based on an SEM-FK hybrid method , 2014 .

[72]  H. Philip Crotwell,et al.  An Educational Resource for Visualizing the Global Seismic Wave Field , 2013 .

[73]  B. Romanowicz,et al.  Observations of changing anisotropy across the southern margin of the African LLSVP , 2013 .

[74]  Hiroshi Takenaka,et al.  FDM computation of seismic wavefield for an axisymmetric earth with a moment tensor point source , 2006 .

[75]  D. L. Anderson,et al.  Preliminary reference earth model , 1981 .

[76]  J. Woodhouse,et al.  Calculation of seismic displacement fields in self-gravitating earth models—applications of minors vectors and symplectic structure , 2008 .

[77]  B. Romanowicz,et al.  Superplumes from the Core-Mantle Boundary to the Lithosphere: Implications for Heat Flux , 2002, Science.