Home

Description

ObsPyck is a GUI application that is intended to cover the tasks in a standard analysis workflow for seismic events in seismological observatory practice. To achieve this ObsPyck heavily relies on several ObsPy modules.

Data Retrieval

ObsPyck can connect to different servers (FDSNWS, ArcLink, Seedlink, SeisHub, local SDS directory tree), automatically fetches waveforms and metadata (instrument responses, station coordinates) for given time span and network/station IDs and shows all fetched data in a fully interactive screen. Waveform data can also be visualized in form of spectrograms. There is also support for local files (all file formats recognized by ObsPy, metadata has to be supplied in either StationXML or Dataless SEED/XSEED).

Data Analysis

ObsPyck can then be used to set P and S wave picks (or custom defined phases) as well as set amplitude maxima/minima for magnitude estimation. These pick data can then be used to run various external programs (e.g. Hypo Inverse 2000, NLLoc, FocMec) that are addressed via system calls. Local magnitudes are calculated using ObsPy (including instrument corrections with fetched response information). The resulting event location and additional information is displayed together with station magnitudes and locations and in an interactive map.

Upload of Results

Finally the event results can be transformed to a XML document (structurally close to QuakeML) and uploaded if connected to a SeisHub or Jane server.

Installation

Recommended installation on Linux and OSX is via Anaconda. Install Anaconda first, following their instructions. Then add the conda-forge channel and install obspyck. Using a dedicated environment specifically for obspy is strongly recommended:

$ conda config --add channels conda-forge
$ conda create -n obspyck python=2.7 obspyck
Fetching package metadata .............
Solving package specifications: .

Package plan for installation in environment /home/megies/anaconda/envs/obspyck_tmp:

The following NEW packages will be INSTALLED:

    asn1crypto:         0.22.0-py27_0      conda-forge
    backports_abc:      0.5-py27_0         conda-forge
    ca-certificates:    2017.7.27.1-0      conda-forge
[...]
    numpy:              1.13.1-py27_0                 
    obspy:              1.1.0-py27_1       conda-forge
    obspyck:            0.5.0-py27_0       conda-forge
[...]

Proceed ([y]/n)? y

#
# To activate this environment, use:
# > source activate obspyck
#
# To deactivate an active environment, use:
# > source deactivate
#

Then to run obspyck, simply activate that environment and run obspyck:

$ source activate obspyck
(obspyck) $ obspyck -h

Alternatively, releases are also available on pypi and can be installed using pip.

Integrating with Hypoinverse (hypo2000)

Warning: Unfortunately, these instructions might be outdated, see issue tracker for various discussions on how to set up external programs: https://github.com/megies/obspyck/issues?q=is%3Aissue

Need to make a directory for location of Obspyck plugins (e.g. Hypoinverse, Nonlinloc, Focmec) - directory name is defined in the .cfg file – ‘pluginpath’
In the plug-in directory, you need the following: crustal velocity model (e.g. .crh); HypoInverse control file (e.g. .inp); hypoInverse binary executable file (see below)
Binary executable: source code available from ftp://ehzftp.wr.usgs.gov/klein/hyp1.40/; pre-compiled Linux version from http://faldersons.net/Software/Hypoinverse/Hypoinverse.html; Linux binary should be named “hyp2000__Linux__64bit” for Linux 64-bit or “hyp2000__Darwin__64bit” for Mac OSX
Hypoinverse working files are stored in temp directory determined by Python (normally /var/)

Screenshots

Warning: Unfortunately, screenshots might be outdated..

Clicking 'Overview' plots the Z trace of every stream.

After zooming in using the mouse wheel (shift + mouse wheel to zoom amplitude) picks can be set using the mouse buttons (or at the mouse position with user-defined keys).

After setting all picks the external location algorithm can be run e.g. by clicking 'do NLLoc'. The event location and information on residuals and magnitudes are displayed in an interactively zoomable map.