OrientPy is a package containing Python tools to estimate seismometer
orientation using automated (and manual) processing of earthquake data. These
methods are particularly useful for broadband ocean-bottom seismic stations,
but are also applicable to broadband land stations or shorter period instruments
(depending on the method selected). Currently the toolbox includes the following methods:
DL (Doran and Laske, 2017): Based on Rayleigh-wave polarization at a range of frequencies and for the two fundamental mode Rayleigh wave orbits.
BNG (Braunmiller, Nabelek and Ghods, 2020): Based on P-wave polarization of earthquakes at regional teleseismic distances.
LKSS (Lim et al., 2018): Based on a harmonic decomposition of radial and transverse receiver functions near zero lag times.
Each method can be used independently to produce an estimate of station orientation, in
terms of the azimuth of seismic component ‘1’ (or ‘N’). The code uses
StDb package for querying and building a station database
used in command-line scripts. Tutorials are provided in the sections below.
DL is heavily based on the code
Doran and Laske (2017). The original code can be found
LKSS is part of the package but has not been fully tested
yet and has no accompanying documentation. We don’t recommend its usage at this time.