Scripts for receiver function analysis and ambient noise cross-correlation
Python FORTRAN Other
Latest commit afca090 Jun 19, 2015 @trichter change readme
Failed to load latest commit information.
bin correct for wrong sensitivities Jun 3, 2015
farmtool initial commit Jul 26, 2012
ppsdtool initial commit Jul 26, 2012
rftool initial commit Jul 26, 2012
src initial commit Jul 26, 2012
tests move/delete scripts Nov 16, 2012
tests2 initial commit Jul 26, 2012
util zerotest Apr 30, 2014
LICENSE update license Jun 3, 2015
README.rst change readme Jun 19, 2015
__init__.py initial commit Jul 26, 2012
colormap.py update Nov 15, 2012
data.py zerotest Apr 30, 2014
debug.py initial commit Jul 26, 2012
events.py major update Jan 16, 2014
imaging.py major update Jan 16, 2014
map.py correct for wrong sensitivities Jun 3, 2015
noise_migration.py initial commit Jul 26, 2012
noisexcorr.py zerotest Apr 30, 2014
psd.py initial commit Jul 26, 2012
rf.py update Feb 14, 2013
seismometer.py update Nov 15, 2012
stations.py major update Jan 16, 2014
stream.py update+relicense Aug 26, 2013
trace.py update Feb 14, 2013
xcorr.py recent changes Sep 11, 2012

README.rst

sito - Scripts for receiver function analysis and ambient noise cross-correlation - based on ObsPy

author: Tom Eulenfeld
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