Package: rfsed
A python package for receiver functions analysis and dealing with sediment effectsrfsed is developed specifically to implement different techniques ofanalysing receiver functions
from stations overlying sedimentary layers.The package isadaptable, efficient, and easy-to-usefor different analysis receiver functions
obtained from stations overlying sedimentary layers.Receiver functions techniques implemented in rfsed are:H-k stacking of Zhu and Kanamori (2000)Sequential H-k stacking 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:Creating publication quality figures for the results of the analysis.Extracting earthquake data from local seismic record files.Multiprocessing options for waveform fitting and extracting earthquake data from local seismic record files.Receiver functions streams in the rfsed are handled by the 'RFStream' class of therfopen 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 onrf.
The rf framework is included in the rfsed package for completness under The MIT License (MIT) Copyright (c) 2013-2019 Tom Eulenfeld. 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 onObsPy.The receiver function forward calculation in rfsed is done using SynR, which is modified after theSynSeis module - seispyproject (under the GNU GENERAL PUBLIC LICENSE).Installation of this packageInstallation from PyPiThe easiest way to install rfsed is viapip_::pipinstallrfsedGetting startedAccess rfsed's documentationhere.rfsed comes with tutorials that demonstrates all its methods. You can find those in theexamples/directory.Reporting Bugs / Contact the developersThis version is an early release of rfsed. If you encounter any issues or unexpected behaviour, pleaseopen an issueon GitHub.Questions?If you have any questions about the package, please use thediscussions featureContributingAll contributions are welcome ... e.g. report bugs, discuss or add new features.Citing rfsedIf you found this package useful, please consider citing it.Akinremi S., van der Meijde, M., Thomas, C., Afonso, J. C., Ruigrok, E., & Fadel, I. (2024). Waveform fitting of receiver functions for enhanced retrieval of crustal structure in the presence of sediments. Journal of Geophysical Research: Solid Earth (under review)ReferencesYeck, W. L., Sheehan, A. F., & Schulte-Pelkum, V. (2013). Sequential h-k stacking to obtain accurate crustal thicknesses beneath sedimentary basins. Bulletin of the Seismological Society of America, 103 , 2142-2150. doi: 10.1785/0120120290Yu, Y., Song, J., Liu, K. H., & Gao, S. S. (2015). Determining crustal structure beneath seismic stations overlying a low-velocity sedimentary layer using receiver functions. Journal of Geophysical Research: Solid Earth, 120 , 3208-3218. doi: 10.1002/2014JB011610Zhu, L., & Kanamori, H. (2000). Moho depth variation in southern California from teleseismic receiver functions. Journal of Geophysical Research: Solid Earth, 105, 2969-2980. doi: 10.1029/1999jb900322Related receiver function projectsrfincluding calculation of receiver functionsseispyincluding hk-stackingRFPyincluding hk-stacking, harmonic decompositionBayHunterinversion of receiver functions and surface wave dispersiontelewavesimsyntheticsPyGLImERincluding common conversion point imaging