There are several publically available software modules that can be useful for confirming or correcting data and metadata. If you calculate better clock drift or sensor orientation values than those provided in the data center, inform the data center of the new values and the method you used to calculate them. It would be best if all of these codes worked on the same data/metadata formats (for example, SeisComp Data Structure and StationXML)
Few (none?) of these tools are OBS-specific, but they can be used for OBS-specific tasks
obsinfo is a tool for creating StationXML files that include OBS-specific information
The instrument clock drift (relative to a reference instrument) can be calculated using noise correlation (REF_SITE) or using the observed drift in hypocenter traveltime residuals.
Determining the horizontal orientation of the instrument using teleseisms, local earthquakes and/or active seismic sources. A common method used is the Rayleigh wave polarization method, outlined in the Stachnik et al., 2012 paper1 with code made available on the OBSIP website.
Removal of noise caused by infragravity waves (based on pressure measurements) and dynamic tilt (based on correlation with horizontal channel noise).
A program to input continuous data-center-level data and a shot file (in some standardized format) and output SEG-Y files.
| Subject | Software | Comments |
|---|---|---|
| Modify miniseed | msmod | Application of leapseconds, header modification |
| " | qedit | Linear and piecewise linear time correction, application of leap seconds, header modification |
| " | GIPP tools | toolbox, e.g. modify header, cut in pieces, export header to ASCII, etc. |
| Clock drift | sara hable's code? | |
| " | justCorrel (Hanneman 2014) | |
| " | ObsPy (Bayreuther 2010, Krischer 2015) | |
| " | Integrated Seismic Program | |
| " | Fortran code after Weemstra (2020) | |
| Sensor orientation | ppol | Based Scholz et al. (2017), plus event location uncertainty |
| " | OBSIP orientation | Stachnik et al. code |
| " | OrientPy | |
| Noise removal | ATACR | matlab version |
| " | ATACR | python version |
| " | tiskitpy | Wayne Crawford |
| " | BRUIT-FM Toolbox | Improved windowing selection and data stacking for more accurate transfer function calculation. |
| " | [Compy] | |
| Station localization | OBSrange | matlab and python |
| Conversion to active seismic data format | Does not exist, should be fairly simple using obspy |
Footnotes
-
Stachnik, J.C., A.F. Sheehan, D.W. Zietlow, Z. Yang, J. Collins, and A. Ferris (2012), Determination of New Zealand Ocean Bottom Seismometer Orientation via Rayleigh- Wave Polarization: Seismol. Res. Lett., 83, 704–712, doi:10.1785/0220110128. ↩