Open AccessComputation of complete synthetic seismograms for laterally heterogeneous models using the Direct Solution Method
Author(s) -
Cummins Phil R.,
Takeuchi Nozomu,
Geller Robert J.
Publication year - 1997
Publication title -
geophysical journal international
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.302
H-Index - 168
eISSN - 1365-246X
pISSN - 0956-540X
DOI - 10.1111/j.1365-246x.1997.tb00983.x
Subject(s) - seismogram , computation , inversion (geology) , algorithm , waveform , perturbation (astronomy) , computer science , synthetic data , mathematics , geology , seismology , physics , telecommunications , radar , quantum mechanics , tectonics
SUMMARY We use the Direct Solution Method (DSM) together with the modified operators derived by Geller & Takeuchi (1995) and Takeuchi, Geller & Cummins (1996) to compute complete synthetic seismograms and their partial derivatives for laterally heterogeneous models in spherical coordinates. The methods presented in this paper are well suited to conducting waveform inversion for 3‐D Earth structure. No assumptions of weak perturbation are necessary, although such approximations greatly improve computational efficiency when their use is appropriate. An example calculation is presented in which the toroidal wavefield is calculated for an axisymmetric model for which velocity is dependent on depth and latitude but not longitude. The wavefield calculated using the DSM agrees well with wavefronts calculated by tracing rays. To demonstrate that our algorithm is not limited to weak, aspherical perturbations to a spherically symmetric structure, we consider a model for which the latitude‐dependent part of the velocity structure is very strong.