Seispy: Python module for batch calculation and postprocessing of receiver functions

Abstract

Seispy is a graphical interface Python module for receiver function (RF) calculation and postprocessing in seismological research. Automated workflows of RF calculations facilitate processing large volume of different types of seismic data. The graphical user interface enables an intuitive and straightforward evaluation of RF quality. All parameters about the preprocessing for RF estimation can be adjusted based on user preference. Water-level frequency-domain deconvolution and iterative time-domain deconvolution for RF estimation are available in Seispy. The current version of Seispy contains five main modules for the postprocessing of RF, such as H-κ stacking, crustal anisotropic estimation, harmonic decomposition, and 2D and 3D common conversion point (CCP) stacking. The CCP stacking in the different application scenarios can be handled by a rich collection of modules, such as time-to-depth conversion, 2D or 3D CCP stacking, and adaptive station or bin selection for CCP stacking profiles in a dense seismic array or a linear seismic array. As a Python module, functions in the Seispy can be called easily in Python scripts for other purposes. The modular design allows new functionality to be added in a collaborative development environment. Licensed under GPLv3, Seispy allow users and developers to freely use, change, share, and distribute copies of the package.