ACE-A FORTRAN subroutine for analytical computation of effective grid parameters for finite-difference seismic waveform modeling with standard Earth models

Genti Toyokuni, Hiroshi Takenaka

Research output: Contribution to journalArticle

5 Citations (Scopus)

Abstract

Despite the broad use of the heterogeneous finite-difference (FD) method for seismic waveform modeling, accurate treatment of material discontinuities inside the grid cells has been a serious problem for many years. One possible way to solve this problem is to introduce effective grid elastic moduli and densities (effective parameters) calculated by the volume harmonic averaging of elastic moduli and volume arithmetic averaging of densities in grid cells. This scheme enables us to place a material discontinuity in an arbitrary position in the spatial grids. Standard Earth models have made a significant contribution to synthetic seismogram calculations with a variety of numerical procedures such as the FD method. For the FD computation of seismic waveform with these models, we must first ensure accurate treatment of material discontinuities in the radius (or depth). The present paper introduces a FORTRAN subroutine ACE which calculates effective parameters analytically for an arbitrary spatial region in either the radius or depth direction for four major standard Earth models, namely, the PREM, IASP91, SP6, and AK135. This program is intended for all FD users who are concerned with seismic wave simulation for these models.

Original languageEnglish
Pages (from-to)635-643
Number of pages9
JournalComputers and Geosciences
Volume35
Issue number3
DOIs
Publication statusPublished - Mar 2009
Externally publishedYes

Fingerprint

Subroutines
Earth (planet)
discontinuity
elastic modulus
finite difference method
Finite difference method
modeling
Elastic moduli
Seismic waves
synthetic seismogram
seismic wave
parameter
simulation
material

Keywords

  • Effective parameter
  • Finite-difference method (FDM)
  • Fortran
  • Seismic wave
  • Standard Earth model

ASJC Scopus subject areas

  • Information Systems
  • Computers in Earth Sciences

Cite this

@article{595fffaa3f7c490e8b0aa02c7afe9ed1,
title = "ACE-A FORTRAN subroutine for analytical computation of effective grid parameters for finite-difference seismic waveform modeling with standard Earth models",
abstract = "Despite the broad use of the heterogeneous finite-difference (FD) method for seismic waveform modeling, accurate treatment of material discontinuities inside the grid cells has been a serious problem for many years. One possible way to solve this problem is to introduce effective grid elastic moduli and densities (effective parameters) calculated by the volume harmonic averaging of elastic moduli and volume arithmetic averaging of densities in grid cells. This scheme enables us to place a material discontinuity in an arbitrary position in the spatial grids. Standard Earth models have made a significant contribution to synthetic seismogram calculations with a variety of numerical procedures such as the FD method. For the FD computation of seismic waveform with these models, we must first ensure accurate treatment of material discontinuities in the radius (or depth). The present paper introduces a FORTRAN subroutine ACE which calculates effective parameters analytically for an arbitrary spatial region in either the radius or depth direction for four major standard Earth models, namely, the PREM, IASP91, SP6, and AK135. This program is intended for all FD users who are concerned with seismic wave simulation for these models.",
keywords = "Effective parameter, Finite-difference method (FDM), Fortran, Seismic wave, Standard Earth model",
author = "Genti Toyokuni and Hiroshi Takenaka",
year = "2009",
month = "3",
doi = "10.1016/j.cageo.2008.05.005",
language = "English",
volume = "35",
pages = "635--643",
journal = "Computers and Geosciences",
issn = "0098-3004",
publisher = "Elsevier Limited",
number = "3",

}

TY - JOUR

T1 - ACE-A FORTRAN subroutine for analytical computation of effective grid parameters for finite-difference seismic waveform modeling with standard Earth models

AU - Toyokuni, Genti

AU - Takenaka, Hiroshi

PY - 2009/3

Y1 - 2009/3

N2 - Despite the broad use of the heterogeneous finite-difference (FD) method for seismic waveform modeling, accurate treatment of material discontinuities inside the grid cells has been a serious problem for many years. One possible way to solve this problem is to introduce effective grid elastic moduli and densities (effective parameters) calculated by the volume harmonic averaging of elastic moduli and volume arithmetic averaging of densities in grid cells. This scheme enables us to place a material discontinuity in an arbitrary position in the spatial grids. Standard Earth models have made a significant contribution to synthetic seismogram calculations with a variety of numerical procedures such as the FD method. For the FD computation of seismic waveform with these models, we must first ensure accurate treatment of material discontinuities in the radius (or depth). The present paper introduces a FORTRAN subroutine ACE which calculates effective parameters analytically for an arbitrary spatial region in either the radius or depth direction for four major standard Earth models, namely, the PREM, IASP91, SP6, and AK135. This program is intended for all FD users who are concerned with seismic wave simulation for these models.

AB - Despite the broad use of the heterogeneous finite-difference (FD) method for seismic waveform modeling, accurate treatment of material discontinuities inside the grid cells has been a serious problem for many years. One possible way to solve this problem is to introduce effective grid elastic moduli and densities (effective parameters) calculated by the volume harmonic averaging of elastic moduli and volume arithmetic averaging of densities in grid cells. This scheme enables us to place a material discontinuity in an arbitrary position in the spatial grids. Standard Earth models have made a significant contribution to synthetic seismogram calculations with a variety of numerical procedures such as the FD method. For the FD computation of seismic waveform with these models, we must first ensure accurate treatment of material discontinuities in the radius (or depth). The present paper introduces a FORTRAN subroutine ACE which calculates effective parameters analytically for an arbitrary spatial region in either the radius or depth direction for four major standard Earth models, namely, the PREM, IASP91, SP6, and AK135. This program is intended for all FD users who are concerned with seismic wave simulation for these models.

KW - Effective parameter

KW - Finite-difference method (FDM)

KW - Fortran

KW - Seismic wave

KW - Standard Earth model

UR - http://www.scopus.com/inward/record.url?scp=60249097467&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=60249097467&partnerID=8YFLogxK

U2 - 10.1016/j.cageo.2008.05.005

DO - 10.1016/j.cageo.2008.05.005

M3 - Article

AN - SCOPUS:60249097467

VL - 35

SP - 635

EP - 643

JO - Computers and Geosciences

JF - Computers and Geosciences

SN - 0098-3004

IS - 3

ER -