We show that the finite-differencing technique based on the consecutive application of the central difference operator to spatial derivatives, a standard well-known technique that has been commonly used in the seismological literature for solving the elastic equation of motion, can also be used to obtain a stable time-domain, finite-difference scheme for solving the anelastic equation of motion. We compare the results of the scheme for a heterogeneous medium with those of the time-domain finite-difference scheme previously developed by Emmerich and Korn and find that they agree very closely. We show, analytically, that in the case of a homogeneous medium, the two schemes give identical numerical results for certain zero initial conditions. The scheme based on the standard technique uses more computer time and memory than the scheme of Emmerich and Korn. However, from a theoretical viewpoint, it is easier to analyze, as it is developed solely with a familiar standard method.

This content is PDF only. Please click on the PDF icon to access.

First Page Preview

First page PDF preview
You do not currently have access to this article.