SEIS_FILO
SEIS_FILO copied to clipboard
SEISmological transdimensional inversion tools for Flat and Isotropic Layered structure in the Ocean
SEIS_FILO
SEISmological transdimensional inversion tools for Flat and Isotropic Layered structures in the Ocean
Copyright (C) 2019-2023 Takeshi Akuhara
The SEIS_FILO program package aims to carry out transdimensional joint inversion of surface waves and receiver functions for ocean-bottom observatories. The main features are:
- Transdimensional MCMC
- Parallel tempering
- Ocean layer
-
Supported input types
- Dispersion curves of the fundamental and higher mode Rayleigh waves
- Rayleigh wave ellipticity
- Rayleigh wave admittance
- P receiver functions (or Green's functions)
- S receiver functions (or Green's functions)
-
Model parameters
- Absolute Vp & Vs
- Vp and Vs anomalies relative to the reference
- Number of layers
- Layer depths
- Standard deviation of noise
Requirements
- FFTW library
- LAPACK library
- Open MPI
- Some Python modules listed in requirements.txt
Install
Type make
in the src
directory. Please edit the Makefile in accordance with your environment (i.e., compiler type and libarary paths).
Quick Guidance
Forward problem
- disper_fwd: Surface wave forward computation (phase and group velocities, ellipticity and admittance)
- recv_func_fwd: Receiver function forward computation
Inverse problem
- joint_inv: Surface wave and receiver function joint inversion by RJMCMC
See online documentation for more details.
Publications
Rayleigh wave dispersion curve and P receiver functions
- Akuhara, T., Yamashita, Y., Ohyanagi, S., Sawaki, Y., Yamada, T., & Shinohara, M. (2023). Shallow Low-Velocity Layer in the Hyuga-Nada Accretionary Prism and Its Hydrological Implications: Insights From a Passive Seismic Array. Journal of Geophysical Research: Solid Earth, 128(4). https://doi.org/10.1029/2022JB026298
Rayleigh wave dispersion curve, ellipticity, and S receiver functions
- Ai, S., Akuhara, T., Morishige, M., Yoshizawa, K., Shinohara, M., & Nakahigashi, K. (2023). Layered Evolution of the Oceanic Lithosphere Beneath the Japan Basin, the Sea of Japan. Journal of Geophysical Research: Solid Earth, 128(2). https://doi.org/10.1029/2022JB025581
Multimode Rayleigh wave dispersion curves
- Yamaya, L., Mochizuki, K., Akuhara, T., Nishida, K. (2021). Sedimentary structure derived from multi-mode ambient noise tomography with dense OBS network at the Japan Trench. Journal of Geophysical Research: Solid Earth, 126(6), e2021JB021789. https://doi.org/10.1029/2021JB021789
S receiver functions
- Akuhara, T., Nakahigashi, K., Shinohara, M., Yamada, T., Shiobara, H., Yamashita, Y., et al. (2021). Lithosphere–asthenosphere boundary beneath the Sea of Japan from transdimensional inversion of S-receiver functions. Earth, Planets and Space, 73(1), 171. https://doi.org/10.1186/s40623-021-01501-5