WaveDec: a tool for wavefield decomposition

WaveDec is a software for the analysis of seismic surface waves. It provides maximum likelihood estimates of wave parameters for Love and Rayleigh waves. All the sensor components are analysed jointly. The simultaneous presence of multiple modes of propagation and multiple wave types (i.e., Love and Rayleigh) is accounted for. Estimated wave parameters include: Love and Rayleigh wave wavenumber, Rayleigh wave ellipticity angle, wave amplitude, and direction of arrival.

An overview of the method is available in the Scientist’s Documentation, more details can be found in the referenced literature. See also the frequently asked questions.

The are are two main scripts:

• WaveDec.py can be used for the analysis of Love and Rayleigh waves and relies on the assumption of plane wavefronts (e.g. at large distance from the source). Translational and rotational motions can be jointly analysed. Applications include the analysis of ambient vibrations. See User’s guide: WaveDec.

• WaveDecActive.py can be used for the analysis of Rayleigh waves with circular wavefronts (e.g. at small distance from the source). Applications include the analysis of active surveys. See User’s guide: WaveDecActive.

Look at the Installation and at the A very quick start sections to run WaveDec as quickly as possible. A simple tool for picking, filtering and plotting is also provided, see User’s guide: wdPicker.

Some Usage Examples are also provided.

The algorithm was developed at the Swiss Seismological Service of ETH Zurich.

Download

The code can be cloned from the GitHub repository github.com/StefanoMarano/WaveDec

git clone https://github.com/StefanoMarano/WaveDec.git

or downloaded as ZIP archive.

Documentation and examples are available in these pages. For any questions related to the software, feel free to email me: .

Usage of the software is limited to non-commercial use (that is, educational and research purposes). If you use the software in your work, please acknowledge us, see How to cite us.

Last Changes

• 26/4/2020: GitHub repository and documentation are online

• 15/7/2017: Created a new page: Frequently asked questions (FAQ).

• 15/6/2017: Updated the examples concerning rotational motions, see Section Plane waves (translational and rotational motions).

• 29/5/2017: Following the acceptance for publication of our article “Analysis of Rayleigh Waves with Circular Wavefront: a Maximum Likelihood Approach”, the code WaveDecActive.py is released in the main distribution together with WaveDec.py.

• 12/5/2017: Added support for the option Tstart in WaveDec. See configuration file section.

Contents

• Scientist’s Documentation
  • Signal model
    • Love wave
    • Rayleigh wave
    • Relationship Between Classical H/V Ellipticity and Ellipticity Angle
  • Overview of the WaveDec method
    • Measurement model
    • Maximum likelihood estimation
    • Model selection
  • Example applications
    • SESAME M2.1
    • Davos (SDAK)
    • Zurich Airport (SKLW) - Passive
    • Zurich Airport (SKLW) - Active
    • Baar (SBAS)
    • Rayleigh wave ellipticity at different sites
• Frequently asked questions (FAQ)
  • What are the main advantages of WaveDec and WaveDecActive?
  • What are the limitations of WaveDec?
  • What are the limitations of WaveDecActive?
  • Why the analysed frequencies are not what I expected?
• Installation
  • Download
  • Requirements
    • GNU/Linux (Ubuntu >14.04)
    • GNU/Linux (Ubuntu <14.04)
    • Windows
    • Other platforms
  • Installation
    • GNU/Linux (Ubuntu)
    • Windows
  • Troubleshooting
• A very quick start
  • With no options
  • With few command line options
  • More options
• User’s guide: WaveDec
  • Input files
    • Configuration file
    • Data
  • Processing
    • Model selection: Love wave, Rayleigh wave, or nothing?
    • Speeding it up
  • Output files
    • Estimated wavefield parameters
    • Resolution limits
    • Array Layout
    • Log file
  • Analysis of the output files
• User’s guide: WaveDecActive
  • Input files
    • Configuration file
    • Data
  • Output files
    • Estimated wavefield parameters
    • Resolution limits
    • Array Layout
    • Log file
  • Analysis of the output files
• User’s guide: wdPicker
  • Analysis of the output files
    • Plotting
    • Filtering
    • Picking
    • Options
• Usage Examples
  • The M2.1 SESAME synthetic dataset (ambient vibrations)
  • The Brigerbad dataset (ambient vibrations)
  • Zurich Airport (active and passive)
  • Debugging examples
    • Plane waves (translational motions)
    • Plane waves (translational and rotational motions)
    • Circular wavefront (active seismic)
• Scientific publications related to WaveDec
  • How to cite us
  • Bibliography
• Code Documentation
  • Reading of input files
  • Synthetic wavefield generation
  • Estimation routines

Indices and tables

• Index

• Module Index

• Search Page