Python module of seismology and receiver functions
Seispy is a Python module for processing seismological data and calculating Receiver Functions. The advanced functions are available to improve the Obspy.
For the use of the Seispy package, please cite as:
For 3D time-difference correction, please also consider citing:
Xu, M., Huang, H., Huang, Z., Wang, P., Wang, L., Xu, M., ... & Yuan, X. (2018). Insight into the subducted Indian slab and origin of the Tengchong volcano in SE Tibet from receiver function analysis. Earth and Planetary Science Letters, 482, 567-579.
Xu, M., Huang, Z., Wang, L., Xu, M., Mi, N., & Yu, D. (2020). Lateral variation of the mantle transition zone beneath the Tibetan plateau: Insight into thermal processes during Indian–Asian collision. Physics of the Earth and Planetary Interiors, 301, 106452.
For 2D and 3D CCP stacking, please also consider citing:
See Seispy documentation in detail.
seispy.distaz
: Calculate distance and azimuth credited by the lithospheric seismology program at USC, but numpy.ndarray
operations are supported.
seispy.geo
: Tiny codes of geophysics.
seispy.decon
: Functions of deconvolution transferred from iwbailey/processRFmatlab including
seispy.rf
: Procedure for RF calculation. The functions of match_eq
, search_eq
invoked obspy.core.UTCDateTime
and obspy.clients
from the Obspy.
seispy.eq
: RF processing for each event, which invoked obspy.io.sac
, obspy.signal
, obspy.taup
and obspy.core.Stream
from the Obspy.
seispy.hk
: H-k stacking for single station (Zhu and Kanamori 2000 JGR).
seispy.rfani
: A joint method for crustal anisotropic calculation (Liu and Niu 2011 GJI).
seispy.slantstack
: Slant stacking for single station (Tauzin et al., 2008)
seispy.rfcorrect
: Subsequent process of RFs including moveout correction and time to depth conversion (1D and 3D) (see Xu et al., 2018 EPSL)
seispy.ccpprofile
: CCP stacking along a profile.
seispy.ccp3d
: 3-D CCP stacking with extracting depth D410 and D660.
prf
: Calculate PRFs for a station.pickrf
: Pick PRFs with virtual quality control after the calculation.plotrt
: Plot PRFs with R and T components order by back-azimuth.plotr
: Plot PRFs with R component order by back-azimuth.hk
: H-Kappa stacking for estimating Moho depth and crustal Vp/Vs.rf2depth
: Convert PRFs to depth axis.ccp_profile
: Stack PRFs along a profile with a CCP stacking method.ccp3d
: Stack PRFs with spaced bins.rfani
: Estimating crustal anisotropy with a joint method.rfharmo
: Harmonic decomposition to extract constant component of RF and plot dip/anisotropic components.pickdepth
: Pick depth of stacked PRFs following ccp3d
veltxt2mod
: Create 3D velocity model with numpy.lib.npyio.NpzFile
format from a ASCII table file.downlod_catalog
: Download catalogs from FDSN web-service.gen_rayp_lib
: Generate a rap-parameter library with depth of source and epicentral distance.setpar
: Set up the values in configure files.