Skip to content

Latest commit

 

History

History
35 lines (25 loc) · 3.01 KB

README.md

File metadata and controls

35 lines (25 loc) · 3.01 KB

MATnoise

MATLAB package to perform ambient noise tomography

This package consists of two parts: (1) Calculation of ambient noise cross-spectra and measuring interstation phase velocities and (2) inverting interstation velocities for 1D or 2D isotropic phase velocity maps and 1D or 2D azimuthal anisotropy.

1. Ambient noise & phase velocities (a1-a7)

  • a1 - Calculate ambient noise cross correlations (vertical, radial, transverse) in the frequency domain.

    Processing options:

    Assumed directory structure for data: {datadirectory}/{station}/{station}.{yyyy}.{jday}.{hh}.{mm}.{SS}.{COMP}.sac
    e.g.: mydata/CC05/CC05.2018.112.00.00.00.BDH.sac

  • a2 - Plot empirical Green's functions in the time domain

  • a3 - Apply group velocity window around the surface wave arrivals of interest. This was originally developed in order to remove the water column arrival from ocean bottom datasets but should also clean up land spectra.

  • a4 - Plot power spectral density of the cross spectra in order to determine the frequency content of the dominant signal.

  • a5 - Plot the real and imaginary parts of the cross spectra to evaluate where signal is best and degree of bias from inhomogeneous noise sources.

  • a6 - Fit J_0 Bessel functions to real part of cross-spectra to extract interstation phase velocity dispersion. Can use results from the linear radon transform (see ./mat-LRTdisp) as the starting dispersion model.

    Menke & Jin (2015) BSSA; DOI:10.1785/0120140245

  • a7 - Try fitting the phase velocities assuming 1D isotropic velocity variations and 1D azimuthal anisotropy. Azimuthal anisotropy includes options for both 2-theta and 4-theta variations.

2. ray_tomo: Ray tomography for 2D phase velocity maps (b0-b1)

  • b0 - Plot ray paths between station pairs colored by velocity.
  • b1 - Use simple travel-time ray tomography to invert for 2D maps of isotropic phase velocity and 1D or 2D azimuthal anisotropy.

Acknowledgements:

These codes have evolved over the years with the help of several people including Ge Jin, Natalie Accardo, and Pei-Ying (Patty) Lin.