Skip to content
/ pysep Public
forked from adjtomo/pysep

seismogram extraction and processing using obspy

License

Notifications You must be signed in to change notification settings

celsoa/pysep

 
 

Folders and files

NameName
Last commit message
Last commit date

Latest commit

 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 

Repository files navigation

SWEEP: a utility for seismic waveform discovery, extraction and processing

What it does

  • Search + fetch + process seismic waveforms from online repositories
  • Anywhere on Earth* for any time period and search radius
  • Prepare and QA waveforms for various analyses, in particular earthquake source estimation

(*NEXT: Earth's Moon and Mars)

How it works

  1. check for earthquakes automatically* or enter earthquake information manually
  2. fetch all possible waveform data across open networks (and some closed networks)
  3. Perform various checks on the data and prepare the waveforms for analysis

(*Using a built-in utility, or coupling with an external one)

Installation

git clone [email protected]:celsoa/sweep.git

Usage

cd sweep
python run_sweep.py <inputfile>
  • This command will extract the waveforms for the geographic location and time listed in inputfile
  • When the script finishes, it will save the waveforms in a folder with the same name as the time listed in inputfile
  • For example, if the time is 2000-12-24T23:59:00 the output folder will be named 20001224235900

Waveform corrections, QA checks, summary reports

  • correct for instrument response
  • demean, detrend (linear, quadratic, etc)
  • resample
  • pre- and post-filter waveforms (causal, acausal, lowpass, highpass, etc)
  • interpolate missing data
  • rescale waveform amplitudes (mm, cm, m, etc)
  • various status/summary/QA reports: waveform quality, SNR, spectrograms, instrument response, station map, station .xml file

Waveform processing for earthquake analysis

  • rotate to source-receiver frame (RTZ) using sensor orientations
  • rotate to UVW triaxial orthogonal frame of the T120PH Nanometrics sensor (useful for isolating spurious signals to sensor components)
  • add null traces so that there are always 3 components (required when rotating)
  • access embargoed data sets (requires user name and password)
  • write SAC files with headers
  • interface with other databases (see below)

Interface with various data centers and databases

Open data centers

  • IRIS
  • EIDA
  • NORSAR
  • Raspberry Shake

Restricted data centers or embargoed datasets (usually require a password)

  • IMS
  • ?

The LLNL database of nuclear explosion and earthquake waveforms

About

seismogram extraction and processing using obspy

Resources

License

Stars

Watchers

Forks

Releases

No releases published

Packages

No packages published

Languages

  • Python 98.2%
  • MATLAB 1.6%
  • Shell 0.2%