# Seispy Documentation¶

Seispy is a Python module for processing seismological data and calculating Receiver Functions. The advanced functions are available to improve the Obspy.

I have been writing Seispy when I was a master student. At first, I wanted to calculate Receiver Functions in Python, but there is no suitable toolkit. Fortunately, The Obspy provided mounts of APIs for processing seismic data, so I ported codes for calculating Receiver Functions from Matlab to Python. Today increased functions have been added to Seispy to further process seismic data over than Obspy.

## Libraries¶

• `seispy.distaz`: Calculate distance and azimuth (by the lithospheric seismology program at USC).

• `seispy.geo`: Tiny codes of geophysics.

• `seispy.bootstrap`: Bootstrap confidence interval estimation (by scikits-bootstrap)

• `seispy.decon`: Functions of deconvolution transferred from iwbailey/processRFmatlab including

• Iterative time domain deconvolution method (Ligorría and Ammon 1999 BSSA).

• Water level frequency domain deconvolution method (CJ. Ammon 1991 BSSA)

• `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.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.

Configurations & Commands:

# References¶

Ammon C J. The isolation of receiver effects from teleseismic P waveforms[J]. Bulletin-Seismological Society of America, 1991, 81(6): 2504-2510.

Krischer L, Megies T, Barsch R, et al. ObsPy: A bridge for seismology into the scientific Python ecosystem[J]. Computational Science & Discovery, 2015, 8(1): 014003.

Ligorría J P, Ammon C J. Iterative deconvolution and receiver-function estimation[J]. Bulletin of the seismological Society of America, 1999, 89(5): 1395-1400.

Xu M, Huang H, Huang Z, et al. Insight into the subducted Indian slab and origin of the Tengchong volcano in SE Tibet from receiver function analysis[J]. Earth and Planetary Science Letters, 2018, 482: 567-579.