Locationg Tectonic Tremors with Uncertainty Estimates
Copyright (C) 2023 Takeshi Akuhara
Welcome to the HypoTremorMCMC repository! This software, written in modern Fortran (F2008 style), is designed to determine the source locations of tremor signals with uncertainty estimates. For a detailed explanation of the method, please refer to our published paper:
(Akuhara et al. 2023, GJI).
Please note: This software is currently under development, and the documentation may not be complete at this time. Changes to the source code or documentation may occur without prior notification.
To install the software, ensure you have a Fortran compiler, OpenMPI, and FFTW installed on your system. To get started, clone the repository to your local machine by running:
git clone https://github.com/akuhara/HypoTremorMCMC.git
Next, compile the program by running the following command in the src directory of the repository:
cd HypoTremorMCMC/src
make
Please note that you may need to edit the Makefile to match your environment.
To estimate tremor locations, follow these six steps:
- Envelope Conversion: Read continuous seismic records of two horizontal components from SAC format files and convert them into a smoothed envelope.
mpirun -np [process number] hypo_tremor_convert [parameter file]
- Calculate Cross-correlation: Calculate cross-correlation functions of the envelope between all station pairs.
mpirun -np [process number] hypo_tremor_correlate [parameter file]
- Measure and Optimize Time- and Amplitude-difference: Measure the arrival time- and amplitude-difference between station pairs and then optimize them to obtain station-specific relative measurements.
mpirun -np [process number] hypo_tremor_measure [parameter file]
- Select Good-quality Events: Select events with good-quality based on rough estimates of propagation speed and attenuation strengths of seismic waves.
mpirun -np [process number] hypo_tremor_select [parameter file]
- Perform MCMC: Perform Bayesian inversion using Markov-chain Monte Carlo (MCMC) method to obtain the posterior probability of source locations.
mpirun -np [process number] hypo_tremor_mcmc [parameter file]
- Statistical Analysis: Extract statistical information from MCMC samples.
mpirun -np [process number] hypo_tremor_statistics [parameter file]
For detailed instructions on how to use this software, please visit the online documentation.
Contributions to the software are welcome. Please open an issue or submit a pull request on the GitHub repository.