An efficient Matlab script to calculate heterogeneous anisotropically elastic wave propagation in three dimensions

Computers & Geosciences
By:

Links

Abstract

We have created a second-order finite-difference solution to the anisotropic elastic wave equation in three dimensions and implemented the solution as an efficient Matlab script. This program allows the user to generate synthetic seismograms for three-dimensional anisotropic earth structure. The code was written for teleseismic wave propagation in the 1-0.1 Hz frequency range but is of general utility and can be used at all scales of space and time. This program was created to help distinguish among various types of lithospheric structure given the uneven distribution of sources and receivers commonly utilized in passive source seismology. Several successful implementations have resulted in a better appreciation for subduction zone structure, the fate of a transform fault with depth, lithospheric delamination, and the effects of wavefield focusing and defocusing on attenuation. Companion scripts are provided which help the user prepare input to the finite-difference solution. Boundary conditions including specification of the initial wavefield, absorption and two types of reflection are available. ?? 2005 Elsevier Ltd. All rights reserved.
Publication type Article
Publication Subtype Journal Article
Title An efficient Matlab script to calculate heterogeneous anisotropically elastic wave propagation in three dimensions
Series title Computers & Geosciences
DOI 10.1016/j.cageo.2005.06.019
Volume 32
Issue 2
Year Published 2006
Language English
Larger Work Type Article
Larger Work Subtype Journal Article
Larger Work Title Computers and Geosciences
First page 259
Last page 264
Google Analytic Metrics Metrics page
Additional publication details