A note on the equivalence of three major propagator algorithms for computational stability and efficiency
Yanlu Ma1, 2, Rongjiang Wang3, Huilan Zhou4
1 Institute for Geology, Mineralogy and Geophysics, Ruhr-University, Bochum 44801, Germany 2 China Earthquake Networks Center, Beijing 100045, China 3 GeoForschungsZentrum Potsdam (GFZ), Telegrafenberg, D-14473 Potsdam, Germany 4 College of Earth Sciences, Graduate University of Chinese Academy of Sciences, Beijing 100049, China
A note on the equivalence of three major propagator algorithms for computational stability and efficiency
Yanlu Ma1, 2, Rongjiang Wang3, Huilan Zhou4
1 Institute for Geology, Mineralogy and Geophysics, Ruhr-University, Bochum 44801, Germany 2 China Earthquake Networks Center, Beijing 100045, China 3 GeoForschungsZentrum Potsdam (GFZ), Telegrafenberg, D-14473 Potsdam, Germany 4 College of Earth Sciences, Graduate University of Chinese Academy of Sciences, Beijing 100049, China
摘要It is shown in this note that the three methods, the orthonormalization method, the minor matrix method and the recursive reflection-transmission matrix method are closely related and solve the numerical instability in the original Thomson-Haskell propagator matrix method equally well. Another stable and efficient method based on the orthonormalization and the Langer block-diagonal decomposition is presented to calculate the response of a horizontal stratified model to a plane, spectral wave. It is a numerically robust Thomson-Haskell matrix method for high frequencies, large layer thicknesses and horizontal slownesses. The technique is applied to calculate reflection-transmission coefficients, body wave receiver functions and Rayleigh wave dispersion.
Abstract:It is shown in this note that the three methods, the orthonormalization method, the minor matrix method and the recursive reflection-transmission matrix method are closely related and solve the numerical instability in the original Thomson-Haskell propagator matrix method equally well. Another stable and efficient method based on the orthonormalization and the Langer block-diagonal decomposition is presented to calculate the response of a horizontal stratified model to a plane, spectral wave. It is a numerically robust Thomson-Haskell matrix method for high frequencies, large layer thicknesses and horizontal slownesses. The technique is applied to calculate reflection-transmission coefficients, body wave receiver functions and Rayleigh wave dispersion.
基金资助:National Natural Science Foundation of China (Nos. 40374009 and 40574024)
通讯作者:
Yanlu Ma
E-mail: yanlu.ma@rub.de
引用本文:
Yanlu Ma, Rongjiang Wang, Huilan Zhou. A note on the equivalence of three major propagator algorithms for computational stability and efficiency[J]. 《地震学报》英文版, 2012, 25(1): 55-64.
Yanlu Ma, Rongjiang Wang, Huilan Zhou. A note on the equivalence of three major propagator algorithms for computational stability and efficiency. Earthquake Science, 2012, 25(1): 55-64.