Accurate and efficient modeling of global seismic wave propagation for an attenuative Earth model including the center

Genti Toyokuni, Hiroshi Takenaka

Research output: Contribution to journalArticle

10 Citations (Scopus)

Abstract

We propose a method for modeling global seismic wave propagation through an attenuative Earth model including the center. This method enables accurate and efficient computations since it is based on the 2.5-D approach, which solves wave equations only on a 2-D cross section of the whole Earth and can correctly model 3-D geometrical spreading. We extend a numerical scheme for the elastic waves in spherical coordinates using the finite-difference method (FDM), to solve the viscoelastodynamic equation. For computation of realistic seismic wave propagation, incorporation of anelastic attenuation is crucial. Since the nature of Earth material is both elastic solid and viscous fluid, we should solve stress-strain relations of viscoelastic material, including attenuative structures. These relations represent the stress as a convolution integral in time, which has had difficulty treating viscoelasticity in time-domain computation such as the FDM. However, we now have a method using so-called memory variables, invented in the 1980s, followed by improvements in Cartesian coordinates. Arbitrary values of the quality factor (Q) can be incorporated into the wave equation via an array of Zener bodies. We also introduce the multi-domain, an FD grid of several layers with different grid spacings, into our FDM scheme. This allows wider lateral grid spacings with depth, so as not to perturb the FD stability criterion around the Earth center. In addition, we propose a technique to avoid the singularity problem of the wave equation in spherical coordinates at the Earth center. We develop a scheme to calculate wavefield variables on this point, based on linear interpolation for the velocity-stress, staggered-grid FDM. This scheme is validated through a comparison of synthetic seismograms with those obtained by the Direct Solution Method for a spherically symmetric Earth model, showing excellent accuracy for our FDM scheme. As a numerical example, we apply the method to simulate seismic waves affected by hemispherical variations of P-wavespeed and attenuation in the top 300. km of the inner core.

Original languageEnglish
Pages (from-to)45-55
Number of pages11
JournalPhysics of the Earth and Planetary Interiors
Volume200-201
DOIs
Publication statusPublished - Jun 2012
Externally publishedYes

Fingerprint

seismic waves
finite difference method
seismic wave
wave propagation
wave equation
grids
wave equations
modeling
spherical coordinates
spacing
attenuation
viscoelasticity
Q factor
Cartesian coordinates
inner core
synthetic seismogram
elastic wave
seismograms
viscous fluids
convolution integrals

Keywords

  • Anelastic attenuation
  • Axisymmetric modeling
  • Earth's center
  • Finite-difference method (FDM)
  • Global seismology
  • Memory variables
  • Singularity
  • Waveform modeling

ASJC Scopus subject areas

  • Geophysics
  • Space and Planetary Science
  • Physics and Astronomy (miscellaneous)
  • Astronomy and Astrophysics

Cite this

@article{02df403124914ae39f9fb590f22e5144,
title = "Accurate and efficient modeling of global seismic wave propagation for an attenuative Earth model including the center",
abstract = "We propose a method for modeling global seismic wave propagation through an attenuative Earth model including the center. This method enables accurate and efficient computations since it is based on the 2.5-D approach, which solves wave equations only on a 2-D cross section of the whole Earth and can correctly model 3-D geometrical spreading. We extend a numerical scheme for the elastic waves in spherical coordinates using the finite-difference method (FDM), to solve the viscoelastodynamic equation. For computation of realistic seismic wave propagation, incorporation of anelastic attenuation is crucial. Since the nature of Earth material is both elastic solid and viscous fluid, we should solve stress-strain relations of viscoelastic material, including attenuative structures. These relations represent the stress as a convolution integral in time, which has had difficulty treating viscoelasticity in time-domain computation such as the FDM. However, we now have a method using so-called memory variables, invented in the 1980s, followed by improvements in Cartesian coordinates. Arbitrary values of the quality factor (Q) can be incorporated into the wave equation via an array of Zener bodies. We also introduce the multi-domain, an FD grid of several layers with different grid spacings, into our FDM scheme. This allows wider lateral grid spacings with depth, so as not to perturb the FD stability criterion around the Earth center. In addition, we propose a technique to avoid the singularity problem of the wave equation in spherical coordinates at the Earth center. We develop a scheme to calculate wavefield variables on this point, based on linear interpolation for the velocity-stress, staggered-grid FDM. This scheme is validated through a comparison of synthetic seismograms with those obtained by the Direct Solution Method for a spherically symmetric Earth model, showing excellent accuracy for our FDM scheme. As a numerical example, we apply the method to simulate seismic waves affected by hemispherical variations of P-wavespeed and attenuation in the top 300. km of the inner core.",
keywords = "Anelastic attenuation, Axisymmetric modeling, Earth's center, Finite-difference method (FDM), Global seismology, Memory variables, Singularity, Waveform modeling",
author = "Genti Toyokuni and Hiroshi Takenaka",
year = "2012",
month = "6",
doi = "10.1016/j.pepi.2012.03.010",
language = "English",
volume = "200-201",
pages = "45--55",
journal = "Physics of the Earth and Planetary Interiors",
issn = "0031-9201",
publisher = "Elsevier",

}

TY - JOUR

T1 - Accurate and efficient modeling of global seismic wave propagation for an attenuative Earth model including the center

AU - Toyokuni, Genti

AU - Takenaka, Hiroshi

PY - 2012/6

Y1 - 2012/6

N2 - We propose a method for modeling global seismic wave propagation through an attenuative Earth model including the center. This method enables accurate and efficient computations since it is based on the 2.5-D approach, which solves wave equations only on a 2-D cross section of the whole Earth and can correctly model 3-D geometrical spreading. We extend a numerical scheme for the elastic waves in spherical coordinates using the finite-difference method (FDM), to solve the viscoelastodynamic equation. For computation of realistic seismic wave propagation, incorporation of anelastic attenuation is crucial. Since the nature of Earth material is both elastic solid and viscous fluid, we should solve stress-strain relations of viscoelastic material, including attenuative structures. These relations represent the stress as a convolution integral in time, which has had difficulty treating viscoelasticity in time-domain computation such as the FDM. However, we now have a method using so-called memory variables, invented in the 1980s, followed by improvements in Cartesian coordinates. Arbitrary values of the quality factor (Q) can be incorporated into the wave equation via an array of Zener bodies. We also introduce the multi-domain, an FD grid of several layers with different grid spacings, into our FDM scheme. This allows wider lateral grid spacings with depth, so as not to perturb the FD stability criterion around the Earth center. In addition, we propose a technique to avoid the singularity problem of the wave equation in spherical coordinates at the Earth center. We develop a scheme to calculate wavefield variables on this point, based on linear interpolation for the velocity-stress, staggered-grid FDM. This scheme is validated through a comparison of synthetic seismograms with those obtained by the Direct Solution Method for a spherically symmetric Earth model, showing excellent accuracy for our FDM scheme. As a numerical example, we apply the method to simulate seismic waves affected by hemispherical variations of P-wavespeed and attenuation in the top 300. km of the inner core.

AB - We propose a method for modeling global seismic wave propagation through an attenuative Earth model including the center. This method enables accurate and efficient computations since it is based on the 2.5-D approach, which solves wave equations only on a 2-D cross section of the whole Earth and can correctly model 3-D geometrical spreading. We extend a numerical scheme for the elastic waves in spherical coordinates using the finite-difference method (FDM), to solve the viscoelastodynamic equation. For computation of realistic seismic wave propagation, incorporation of anelastic attenuation is crucial. Since the nature of Earth material is both elastic solid and viscous fluid, we should solve stress-strain relations of viscoelastic material, including attenuative structures. These relations represent the stress as a convolution integral in time, which has had difficulty treating viscoelasticity in time-domain computation such as the FDM. However, we now have a method using so-called memory variables, invented in the 1980s, followed by improvements in Cartesian coordinates. Arbitrary values of the quality factor (Q) can be incorporated into the wave equation via an array of Zener bodies. We also introduce the multi-domain, an FD grid of several layers with different grid spacings, into our FDM scheme. This allows wider lateral grid spacings with depth, so as not to perturb the FD stability criterion around the Earth center. In addition, we propose a technique to avoid the singularity problem of the wave equation in spherical coordinates at the Earth center. We develop a scheme to calculate wavefield variables on this point, based on linear interpolation for the velocity-stress, staggered-grid FDM. This scheme is validated through a comparison of synthetic seismograms with those obtained by the Direct Solution Method for a spherically symmetric Earth model, showing excellent accuracy for our FDM scheme. As a numerical example, we apply the method to simulate seismic waves affected by hemispherical variations of P-wavespeed and attenuation in the top 300. km of the inner core.

KW - Anelastic attenuation

KW - Axisymmetric modeling

KW - Earth's center

KW - Finite-difference method (FDM)

KW - Global seismology

KW - Memory variables

KW - Singularity

KW - Waveform modeling

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

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

U2 - 10.1016/j.pepi.2012.03.010

DO - 10.1016/j.pepi.2012.03.010

M3 - Article

AN - SCOPUS:84860676053

VL - 200-201

SP - 45

EP - 55

JO - Physics of the Earth and Planetary Interiors

JF - Physics of the Earth and Planetary Interiors

SN - 0031-9201

ER -