akuhara / SEIS_FILO

SEISmological transdimensional inversion tools for Flat and Isotropic Layered structure in the Ocean

Geek Repo:Geek Repo

Github PK Tool:Github PK Tool

SEIS_FILO

SEISmological transdimensional inversion tools for Flat and Isotropic Layered structures in the Ocean

LOGO

Build Status codecov GitHub Docker Cloud Build Status Documentation Status DOI

Copyright (C) 2019-2023 Takeshi AkuharaORCID


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

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

About

SEISmological transdimensional inversion tools for Flat and Isotropic Layered structure in the Ocean

License:GNU General Public License v3.0


Languages

Language:Fortran 89.6%Language:Python 8.5%Language:Makefile 1.3%Language:Shell 0.2%Language:Dockerfile 0.2%Language:AMPL 0.2%