sito - Scripts for receiver function analysis and ambient noise cross-correlation - based on ObsPy
These scripts were written in the scope of my PhD at GFZ Postdam/FU Berlin. I do not maintain them any longer, but of course you are free to copy&paste code of interest.
The following publicatons are based upon results of these scripts:
Richter, T., C. Sens-Schönfelder, R. Kind, and G. Asch (2014), Comprehensive observation and modeling of earthquake and temperature related seismic velocity changes in northern Chile with passive image interferometry, J. Geophys. Res. Solid Earth, 119, 4747–4765, doi:10.1002/2013JB010695. pdf1
Richter, T. (2014), Temporal Variations of Crustal Properties in Northern Chile Analyzed with Receiver Functions and Passive Image Interferometry, dissertation, FU Berlin. pdf2
I refactored the receiver function functionality in a separate repository, see rf. Some of the ambient noise correlation code is re-used in the yam package. Alternative Python packages for ambient noise correlation exist.
- Required:
- python 2.7
- numpy
- scipy
- matplotlib
- basemap
- ipython
- obspy
- seispy
- progressbar
sito consist of the following major modules:
data: Class for data handling (You can set up your own data calss here) events: Class for handling event data (better: use new Catalog class in ObsPy) imaging: Customized plotting map: Customized station map plotting noisexcorr: correlation of noise rf: Receiver function calculation station: Class for handling station data stream: Class derived from obspy.trace.Stream with custom methods trace: Class derived from obspy.trace.Trace with custom methods util: Utilities like deconvolution, polarisation, rotation, pspier (mainly imported from other projects) xcorr: Cross correlation functions
minor modules:
noise_migration: migrate noise phases in correlations backwards (experimental) debug: ipython debugging functionality seismometer: PAZ of STS-2 and Wood-Anderson