Software:

1) Rayleigh wave polarization

Example of usage: Cascadia OBS

2) Power Spectral Density and Probability (dbpsdpdf)

1) Rayleigh wave polarization analysis software used to determine the orientation of ocean bottom seismometer horizontal components. Download software.

Stachnik et al. Determination of ocean bottom seismometer orientation via Rayleigh wave polarization. Accepted, SRL, 2012.

The scripts are written in a modular
fashion, processing one earthquake
for one station at a time. There are several accompanying wrapper
and plotting scripts that can be used to display the results and
calculate statistics from processing many earthquakes per station.

1) Extract waveforms containing surface wave arrival from db or other into ascii file t,z,n,e (css_wfextract.pl)

---------------------------------------------------------------------------------

2) Run CSpolar.pl to determine azimuth (between 0 and 360) for highest correlation between Hilbert transformed vertical component and radial component (computed from two horizontal components and assumed azimuth). Output of CSpolar is a file of time, azimuth, Szz, Szr, Srr, Czr, Czr2.

Time: Not used

Azimuth: 0-360

Szz: Zero lag amplitude of ZZ xcorr (see paper)

Szr:

Srr:

Czr: Value from Eq.2

Czr2: Value from Eq. 3

1) Extract waveforms containing surface wave arrival from db or other into ascii file t,z,n,e (css_wfextract.pl)

---------------------------------------------------------------------------------

2) Run CSpolar.pl to determine azimuth (between 0 and 360) for highest correlation between Hilbert transformed vertical component and radial component (computed from two horizontal components and assumed azimuth). Output of CSpolar is a file of time, azimuth, Szz, Szr, Srr, Czr, Czr2.

Time: Not used

Azimuth: 0-360

Szz: Zero lag amplitude of ZZ xcorr (see paper)

Szr:

Srr:

Czr: Value from Eq.2

Czr2: Value from Eq. 3

---------------------------------------------------------------------------------

3) Process many earthquakes.

4) Summarize results. Calculate statistics, determine any azimuthal bias.