/SiteAmplification

Calculate site amplification given a 1D model including body and surface waves

Primary LanguagePython

SiteAmplification

Yan Yang (yanyang@caltech.edu), Spring 2020, Caltech

Citation:

Yang, Yan, James W. Atterholt, Zhichao Shen, Jack B. Muir, Ethan F. Williams, and Zhongwen Zhan. "Sub‐kilometer correlation between near‐surface structure and ground motion measured with distributed acoustic sensing." Geophysical Research Letters 49, no. 1 (2022): e2021GL096503.

Bowden, Daniel C., and Victor C. Tsai. "Earthquake ground motion amplification for surface waves." Geophysical Research Letters 44, no. 1 (2017): 121-127.

This is to calculate site amplification given a 1D model.

Vertically incident body wave site amplification is computed using propagation matrix method (Aki & Richards, 2002).

Horizontally propagating surface wave site amplification is including body and surface waves (Bowden & Tsai, 2017).

The surface wave integrals are calcualated using mat_disperse function from the MATLAB package forked from Yiran Ma

benchmark.m should produce a figure like this:

Screenshot 2023-02-02 at 14 39 30

which is the same as Figure 1C in Bowden and Tsai (2017):

Screenshot 2023-02-02 at 14 39 00