Inversion of receiver function by wavelet transformation

Introduction Receiver function has been extensively applied in studying S wave velocity of crust and up-per mantle for about 20 years (Owens, et al, 1987; LIU, et al, 1996), which is a time series ob-tained by the deconvolution of vertical component from horizontal component for teleseismic P waveform. Receiver function represents the teleseismic P plane wave response of crust and upper mantle beneath seismic station, from which the source and propagation effects are removed. Receiver funct…

Inversion of receiver function by wavelet transformation

A new method for receiver function inversion by wavelet transformation is presented in this paper. Receiver func-tion is expanded to different scales with different resolution by wavelet transformation. After an initial model be-ing taken, a generalized least-squares inversion procedure is gradually carried out for receiver function from low tohigh scale, with the inversion result for low order receiver function as the initial model for high order. Aneighborhood containing the global minimum is firstly searched from low scale receiver function, and will gradu-ally focus at the global minimum by introducing high scale information of receiver function. With the gradual ad-dition of high wave-number to smooth background velocity structure, wavelet transformation can keep the inver-sion result converge to the global minimum, reduce to certain extent the dependence of inversion result on theinitial model, overcome the nonuniqueness of generalized least-squares inversion, and obtain reliable crustal andupper mantle velocity with high resolution.