-
Notifications
You must be signed in to change notification settings - Fork 1
Home
Welcome to the tsuquakes wiki!
This set of scripts has been used in my tsunami earthquake research. I am evaluating the rupture parameters of tsunami earthquakes for use in early tsunami warning (ETW). Near-field geophysical data currently only exist for the 2010 M7.8 Mentawai tsunami earthquake, so I am generating synthetic data patterned after this event. I am varying the rupture parameters stress drop, rise time, and rupture velocity, to determine a mean and std deviation for these parameters that best reproduces a tsunami earthquake scenario. The rupture parameters are evaluated by comparing various intensity measures (IMs) and spectra of the synthetic data to the observed.
The main steps are:
Generate suites of synthetic rupture scenarios and waveforms using code from this repository https://github.com/dmelgarm/MudPy.git (commonly referred to as FakeQuakes). The basic script that runs all of this code is mentawai.fq.py but it can be varied to have different values for rupture parameters. More detailed instructions on how to fill out the script are found on the next wiki page: FakeQuakes Script.
Compute the following IMs: PGD, PGA, PGV, tPGD and tPGA. Bin the Fourier spectra and obtain avg values for the bins. The script used to do all of this is synthetic_calc_mpi.py. This is the parallelized version of the script; however, there is an unparallelized version called synthetic_calc.py. IMs and spectra for the observed data are also needed for Step 3. The script mentawai.py computes these.
Calculate and plot the residuals between the observed and synthetic data. The script used to do this is residuals.py