Documentation

rfsed is developed specifically to implement different techniques of analysing receiver functions from stations overlying sedimentary layers. The software is adaptable, efficient, and easy-to-use for different analysis of receiver functions obtained from stations overlying sedimentary layers.

Receiver functions techniques implemented in rfsed are:

• H-k stacking (one layer) of Zhu and Kanamori (2000)

• Sequential H-k stacking (two layers) of Yeck et al., (2013)

• Resonance filtering and modified H-k stacking of Yu et al., (2015)

• H-k stacking and Waveform Fitting (two-step method) of Akinremi et al., (2024)

• Analysis of the synthetic reciever functions with the above-mentioned methods.

Beside these receiver function methods, rfsed has the following features:

• Extracting earthquake data from local seismic record files.

• Multiprocessing options for waveform fitting and extracting earthquake data from local seismic record files.

• Creating publication quality figures for the results of the analysis.

Receiver functions streams in the rfsed are handled by the RFStream class of the rf open software, and it inherits a lot of useful methods from the Obspy class Stream. It is supported via the obspyh5 package. For more information on class RFStream, see documentation on rf. In the rfsed modules to extract earthquake data from local seismic records, read and write support for necessary metadata is provided for SAC, SeismicHanlder and HDF5 formats based on ObsPy.

This package is fully in python and common processing workflows are described in the Jupyter notebooks accompanying this package.

Quick Links

• GitHub Repository
• References

Getting Started

• License
• Installation
• Usage

Module

• Module hkZhu
• Module hkSeqYeck
• Module ReverbFilter
• Module hkYu
• Module synrf
• Module WaveformFitting
• Module WaveformFittingMultiproc
• Module ExtractEq
• Module ExtractEqMultiproc

Jupyter Notebooks

• Example 1: Downloading Earthquake Waveform Data from IRIS
• Example 2: Calculation of Receiver Functions
• Example 3: Download Catalog of earthquake events
• Example 4: Example of Extracting Earthquake Waveform Data from seismic data from local drive
• Example 5: Example of Extracting Earthquake Waveform Data from seismic data from local drive - Using Multiple Processors
• Example 6: H-K Stacking of Receiver Functions
• Example 7: Sequential H-K Stacking of Receiver Functions
• Example 8: Resonance Filtering and Modified H-K Stacking of Receiver Functions
• Example 9: Calculating Synthetic Reciever Function
• Example 10: Calculation and Analysis of Synthetic Receiver Functions
• Example 11: Waveform Fitting of Receiver Functions
• Example 12: Waveform Fitting of Receiver Functions - Using Multiple Processors