Receiver functions#

Synthetic receiver functions (teleseismic body wave)#

1D modelling#

  • Raysum: Ray-theoretical modelling of teleseismic waves in dipping, anisotropic structures

    • PyRaysum: Teleseismic body wave modeling through stacks of (dipping/anisotropic) layers

    • pyfwrd: A forward modelling code for surface wave, receiver functions and shear wave splitting, given tilted hexagonal symmetric media

  • Telewavesim: Teleseismic body wave modeling through stacks of (submarine/anisotropic) layers

    • RMATRIX is invoked: Finds the R/T coefficient matrices for a stack of anisotropic layers.

  • trftn96 in CPS330: Make a P-wave receiver function for a transverse isotropic model

  • RFTN C.J. Ammon’s Code: An Overview of Receiver-Function Analysis

2D modelling#

  • PSV Hybrid RF: Calculating synthetic RF in two-dimensional localized hetergeneous structures based on PSV Hybrid method (GRT-FD)

3D modolling#

  • SPECFEM3D_Cartesian: Simulates seismic wave propagation at the local or regional scale and performs full waveform imaging (FWI) or adjoint tomography based upon the spectral-element method (SEM). The injection of teleseismic wavefield is avaliable with SEM-FK and AxisSEM method. See COUPLE_WITH_INJECTION_TECHNIQUE in Par_file

    • SEM-FK: Plane wavefield injection coupling with Cartesian SEM domain. Refer to Specfem3D and Full Waveform Inversion for usage.


      Tong, P., Komatitsch, D., Tseng, T.-L., Hung, S.-H., Chen, C.-W., Basini, P., and Liu, Q. (2014), A 3-D spectral-element and frequency-wave number hybrid method for high-resolution seismic array imaging, Geophys. Res. Lett., 41, 7025– 7034, doi:10.1002/2014GL061644.

    • Specfem3d-AxiSEM: 3-D seismic wavefields propagation in axisymmetric media and inject into SEM domain.


      T. Nissen-Meyer, M. van Driel, S. C. Staehler, K. Hosseini, S. Hempel, L. Auer, A. Colombi and A. Fournier: “AxiSEM: broadband 3-D seismic wavefields in axisymmetric media”, Solid Earth, 5, 425-445, 2014 doi:10.5194/se-5-425-2014

      Monteiller V, Beller S, Plazolles B, et al. On the validity of the planar wave approximation to compute synthetic seismograms of teleseismic body waves in a 3-D regional model[J]. Geophysical Journal International, 2021, 224(3): 2060-2076.

Receiver function Process#


  • Seispy: Python Module for seismology and receiver functions

    • Batch calculation

    • GUI picking

    • H-k stacking

    • CCP stacking

    • Azimuth anisotropic estimation

    • Slant stacking

  • rf: Python framework for receiver function analysis.


    Tom Eulenfeld (2020), rf: Receiver function calculation in seismology, Journal of Open Source Software, 5(48), 1808, doi:10.21105/joss.01808

    • Batch calculation

    • CCP stacking

  • RfPy

    • Gaussian weighted CCP stacking

    • H-k stacking

    • Back-azimuth harmonics


  • Funclab: Matlab based GUI for handling receiver functions


    Porritt, R. W. and Miller, M. S., (2018), Updates to FuncLab, a Matlab based GUI for handling receiver functions. Computers and Geoscience, 111, 260-271, doi:10.1016/j.cageo.2017.11.022

  • SplitRFLab: A Matlab toolbox of processing receiver functions and shear wave splitting


    Xu M, Huang H, Huang Z, et al. SplitRFLab: A MATLAB GUI toolbox for receiver function analysis based on SplitLab[J]. Earthquake Science, 2016, 29(1): 17-26.

    • RF codes invoked from processRFmatlab: Matlab functions and scripts for working with receiver functions

  • Crazyseismic: A MATLAB GUI‐based software package for passive seismic data preprocessing


    Yu, C.Q., Zheng, Y.C., Shang, X.F. (2017), Crazyseismic: A MATLAB GUI-based software package for passive seismic data preprocessing, Seismol. Res. Lett., 88(2A), 410-415, doi:10.1785/0220160207


  • CPS330: Collection of programs for calculating theorectical seismogram, receiver function, surface wave dispersion curve et al.


    Herrmann, R. B. (2013) Computer programs in seismology: An evolving tool for instruction and research, Seism. Res. Lettr. 84, 1081-1088, doi:10.1785/0220110096

  • hk & ccp: Prof. Zhu Lupei’s code

    • hk1.3: Receiver function package (deconvolution and H-k stacking)

    • ccp1.0: Common-Conversion-Point (CCP) stacking of receiver functions

  • H-k-c: Generalized H-k after harmonic correction on receiver functions


    Li, J., Song, X., Wang, P., & Zhu, L. (2019). A generalized H-k method with harmonic corrections on Ps and its crustal multiples in receiver functions. J. Geophys. Res. Solid Earth, 124(4), 3782-3801.

Receiver function (joint) inversion#

Bayesian inversion#

  • MC3deconv: Bayeisan inversion to recover Green’s functions of receiver-side structures from teleseismic waveforms


    T. Akuhara, M. Bostock, A. Plourde, M. Shinohara (2019) Beyond Receiver Functions: Green’s Function Estimation by Trans-Dimensional Inversion and Its Application to OBS Data, Journal of Geophysical Research: Solid Earth, doi:10.1029/2018JB016499

  • BayHunter: McMC transdimensional Bayesian inversion of surface wave dispersion and receiver functions


    Dreiling, Jennifer; Tilmann, Frederik (2019): BayHunter - McMC transdimensional Bayesian inversion of receiver functions and surface wave dispersion. GFZ Data Services. doi:10.5880/GFZ.2.4.2019.001

  • RfSurfHmc: Joint inversion of Receiver Function and Surface Wave Disperion by Hamilton Monte Carlo Method

Gradient inversion#

  • CPS330 Collection of programs for calculating theorectical seismogram, receiver function, surface wave dispersion curve et al.

    • rftn96 for receiver function inversion

    • joint96 for joint inversion with surface wave dispersion