The reflectivity method: a tutorial

Main Article Content

G. Muller


An extended reflectivity method is described by which complete seismograms for a point source in a layered half-space can be calculated. Starting with the differential equations and boundary conditions, the reflection and transmission of plane waves at layered media is treated first, followed by the synthesis of point-source wave fields. The frequency-domain displacements of the half-space surface are expressed as slowness integrals, and the most prominent parts of the integrands are the reflectivities of the layers below and above the point source and a function which is closely related to the transmissivity of the layers above the source. Reflectivities and transmissivities are calculated by recursive methods which are numerically stable for all frequencies and slownesses. Near- and far-field results are given for single-force and moment-tensor point sources. From the general results for the complete medium response, partial responses can easily be extracted, e.g. the original form of the reflectivity method which calculates only the response from the layers below the source. Thus, the extended reflectivity method has a flexibility which is not available if propagator methods are used for the calculation of the integrands. Various other aspects of seismogram calculation are addressed, such as extended sources, an earth-flattening transformation and the inclusion of absorption for constant and frequency-dependent Q. Theoretical seismograms are shown, first for body-wave propagation from explosions in a crustal model and in a model which came from seismic prospecting, and second for surface-wave propagation from a double-couple source. Due to the tutorial nature of this paper the derivations are mostly rather detailed. It is hoped that this will help interested newcomers to the field of theoretical seismograms to get started.

Google Scholar           ARK:



Article Details

How to Cite
Muller, G. (1985). The reflectivity method: a tutorial. Journal of Geophysics, 58(1), 153-174. Retrieved from


Aki, K., Richards, P.G. (1980) Quantitative seismology-theory and methods, vol. I and II. San Francisco: Freeman and Co.

Alekseev, A.S., Mikhailenko, B.G. (1980) The solution of dynamic problems of elastic wave propagation in inhomogeneous media by a combination of partial separation of variables and finite-difference methods. J. Geophys. 48:161-172

Apsel, R.J., Luco, J.E. (1983) On the Green's functions for a layered half-space, Part II. Bull. Seismol. Soc. Am. 73:931-951

Ben-Menahem, A., Vered, M. (1973) Extension and interpretation of the Cagniard-Pekeris method for dislocation sources. Bull. Seismol. Soc. Am. 63:1611-1636

Ben-Menahem, A., Singh, S.J. (1981) Seismic waves and sources. Springer, New York

Biswas, N.N., Knopoff, L. (1970) Exact earth-flattening calculation for Love waves. Bull. Seismol. Soc. Am. 60:1123-1137

Bouchon, M. (1979) Discrete wave number representation of elastic wave fields in three-space dimensions. J. Geophys. Res. 84:3609-3614

Bouchon, M. (1981) A simple method to calculate Green's functions for elastic layered media. Bull. Seismol. Soc. Am. 71:959-971

Brilstle, W., Muller, G. (1983) Moment and duration of shallow earthquakes from Love-wave modelling for regional distances. Phys. Earth Planet. Inter. 32:312-324

Campillo, M., Bouchon, M., Massinon, B. (1984) Theoretical study of the excitation, spectral characteristics, and geometrical attenuation of regional seismic phases. Bull. Seismol. Soc. Am. 74:79-90

Chapman, C.H. (1973) The earth. flattening transformation in body wave theory. Geophys. J. R. Astron. Soc. 35:55-70

Chapman, C.H. (1978) A new method for computing synthetic seismograms. Geophys. J. R. Astron. Soc. 54:481-518

Christensen, R.M. (1982) Theory of viscoelasticity - an introduction (2nd ed.). Academic Press, New York

Cormier, V.F. (1980) The synthesis of complete seismograms in an earth model specified by radially inhomogeneous layers. Bull. Seismol. Soc. Am. 70:691-716

Cormier, V.F., Richards, P.G. (1977) Full wave theory applied to a discontinuous velocity increase: the inner core boundary. J. Geophys. 43:3-31

Deichmann, N., Ansorge, J. (1983) Evidence for lamination in the lower continental crust beneath the Black Forest (Southwestern Germany). J. Geophys. 52:109-118

Dziewonski, A.M., Anderson, D.L. (1981) Preliminary reference earth model. Phys. Earth Planet. Inter. 25:297-356

Fuchs, K. (1968) The reflection of spherical waves from transition zones with arbitrary depth-dependent elastic moduli and density. J. Phys. Earth 16, Special Issue, 27-41

Fuchs, K. (1969) On the properties of deep crustal reflectors. J. Geophys. 35:133-149

Fuchs, K., Muller, G. (1971) Computation of synthetic seismograms with the reflectivity method and comparison with observations. Geophys. J. R. Astron. Soc. 23:417-433

Gerver, M., Markushevich, V. (1966) Determination of a seismic wave velocity from the travel-time curve. Geophys. J. R. Astron. Soc. 11:165-173

Ha, J. (1984) Recurrence relations for computing complete P and SV seismograms. Geophys. J. R. Astron. Soc. 79:863-873

Hage, H. (1983) Velocity constraints for the inner core inferred from long-period PKP amplitudes. Phys. Earth Planet. Inter. 31:171-185

Harvey, D.J. (1981) Seismogram synthesis using normal mode superposition: the locked mode approximation. Geophys. J. R. Astron. Soc. 66:37-69

Helmberger, D.V. (1968) The crust-mantle transition in the Bering Sea. Bull. Seismal. Soc. Am. 58:179-214

Hron, F., Mikhailenko, B.G. (1981) Numerical modeling of nongeometrical effects by the Alekseev-Mikhailenko method. Bull. Seismal. Soc. Am. 71:1011-1029

Ingate, S.F., Bock, G., Kind, R. (1983) Synthesis of complete SH seismograms. Geophys. J. R. Astron. Soc. 75:261-274

Kennett, B.L.N. (1974) Reflections, rays and reverberations. Bull. Seismal. Soc. Am. 65:1643-1651

Kennett, B.L.N. (1979) Theoretical reflection seismograms for an elastic medium. Geophys. Prospecting 27:301-321

Kennett, B.L.N. (1980) Seismic waves in a stratified medium II. theoretical seismograms. Geophys. J. R. Astron. Soc. 61:1-10

Kennett, B.L.N. (1983) Seismic wave propagation in stratified media. Cambridge University Press

Kennett, B.L.N., Kerry, N.J. (1979) Seismic waves in a stratified half-space. Geophys. J. R. Astron. Soc. 57:557-583

Kind, R. (1978) The reflectivity method for a buried source. J. Geophys. 44:603-612

Kind, R. (1985) The reflectivity method for different source and receiver structures and comparison with GRF data. J. Geophys. 58:146-152

Kind, R., Odom, R.I. (1983) Improvements to layer matrix methods. J. Geophys. 53:127-130

Korn, M. (1985) Kombination von Integraltransformationen und Finiten-Differenzen-Methoden zur Berechnung der Wellenausbreitung in geschichteten Medien. Ph.D. thesis, University of Frankfurt

Korn, M., Muller, G. (1983) Comparison of the Alekseev-Mikhailenko method and the reflectivity method. Geophys. J. R. Astron. Soc. 72:541-556

Luco, J.E., Apsel, R.J. (1983) On the Green's functions for a layered half-space, Part I. Bull. Seismal. Soc. Am. 73:909-929

Minster, J.B. (1980) Anelasticity and attenuation. Proceedings of Enrico Fermi International School of Physics, A. Dziewonski and E. Boschi, eds.: pp. 152-212. North Holland, Amsterdam

Muller, G. (1969) Theoretical seismograms for some types of point sources in layered media. Part III: Single force and dipole sources of arbitrary orientation. J. Geophys. 35:347-371

Muller, G. (1973) Seismic moment and long-period radiation of underground nuclear explosions. Bull. Seismal. Soc. Am. 63:847-857

Muller, G. (1977a) Earth-flattening approximation for body waves derived from geometric ray theory-improvements, corrections and range of applicability. J. Geophys. 42:429-436

Muller, G. (1977b) Correction. J. Geophys. Res. 82:2541-2542

Muller, G. (1983) Rheological properties and velocity dispersion of a medium with power-law dependence of Q on frequency. J. Geophys. 54:20-29

Olson, A.H., Orcutt, J.A., Frazier, G.A. (1984) The discrete wavenumber/finite element method for synthetic seismograms. Geophys. J. R. Astron. Soc. 77:421-460

O'Neil, M.E., Hill, D.P. (1979) Causal absorption: its effect on synthetic seismograms computed by the reflectivity method. Bull. Seismal. Soc. Am. 69:17-26

Panza, G.F. (1985) Synthetic seismograms: the Rayleigh waves modal summation. J. Geophys. 58:125-145

Pilant, W.L. (1979) Elastic waves in the earth. Elsevier, Amsterdam

Rial, J.A., Cormier, V.F. (1980) Seismic waves at the epicenter's antipode. J. Geophys. Res. 85:2661-2668

Sandmeier, K.-J., Wenzel, F. (1985) Synthetic seismograms for a complex crustal model. Geophys. Res. Lett. 12 (in press)

Spudich, P., Ascher, U. (1983) Calculation of complete theoretical seismograms in vertically varying media using collocation methods. Geophys. J. R. Astron. Soc. 75:101-124

Temme, P., Muller, G. (1982) Numerical simulation of vertical seismic profiling. J. Geophys. 50:177-188

Wang, C.Y., Herrmann, R.B. (1980) A numerical study of P-, SVand SH-wave generation in a plane layered medium. Bull. Seismol. Soc. Am. 70:1015-1036

Woodhouse, J.H. (1980) Efficient and stable methods for performing seismic calculations in stratified media. Proceedings of Enrico Fermi International School of Physics, A. Dziewonski and E. Boschi, eds.: pp. 127-151. North Holland, Amsterdam