Instaseis: instant global seismograms based on a broadband waveform database

We present a new method and implementation (Instaseis) to store global Green's functions in a database which allows for near-instantaneous (on the order of milliseconds) extraction of arbitrary seismograms. Using the axisymmetric spectral element method (AxiSEM), the generation of these databas...

Full description

Saved in:
Bibliographic Details
Published in:Solid earth (Göttingen) 2015-06, Vol.6 (2), p.701-717
Main Authors: van Driel, M, Krischer, L, Stähler, S. C, Hosseini, K, Nissen-Meyer, T
Format: Article
Language:English
Subjects:
Accuracy
Benchmarks
Broadband
Coefficients
Computation
Decomposition
Earth
Earth science
Earthquakes
Green's function
Green's functions
Interpolation
Mathematical models
Methods
Polynomials
Reciprocity
Resampling
Sampling
Seismograms
Source code
Spectral element method
Supercomputers
Waveforms
Citations: Items that this one cites
Items that cite this one
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:We present a new method and implementation (Instaseis) to store global Green's functions in a database which allows for near-instantaneous (on the order of milliseconds) extraction of arbitrary seismograms. Using the axisymmetric spectral element method (AxiSEM), the generation of these databases, based on reciprocity of the Green's functions, is very efficient and is approximately half as expensive as a single AxiSEM forward run. Thus, this enables the computation of full databases at half the cost of the computation of seismograms for a single source in the previous scheme and allows to compute databases at the highest frequencies globally observed. By storing the basis coefficients of the numerical scheme (Lagrange polynomials), the Green's functions are 4th order accurate in space and the spatial discretization respects discontinuities in the velocity model exactly. High-order temporal interpolation using Lanczos resampling allows to retrieve seismograms at any sampling rate. AxiSEM is easily adaptable to arbitrary spherically symmetric models of Earth as well as other planets. In this paper, we present the basic rationale and details of the method as well as benchmarks and illustrate a variety of applications. The code is open source and available with extensive documentation at www.instaseis.net.
ISSN:1869-9529
1869-9510
1869-9529
DOI:10.5194/se-6-701-2015