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
A lot of alternative Python packages for ambient noise cross-correlation exist. I refactored the receiver function functionality in a separate repository, see rf.
- 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