sito - Scripts for receiver function analysis and ambient noise cross-correlation - based on ObsPy
license: MIT license
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