JpGU-AGU Joint Meeting 2020

Presentation information

[E] Oral

S (Solid Earth Sciences ) » S-SS Seismology

[S-SS03] Seismological advances in the ocean

convener:Takashi Tonegawa(Research and Development center for Earthquake and Tsunami, Japan Agency for Marine-Earth Science and Technology), Takeshi Akuhara(Earthquake Research Institute, University of Tokyo), Pascal Audet(University of Ottawa), Taka'aki Taira(Berkeley Seismological Laboratory, University of California, Berkeley)

[SSS03-05] Development of seismic inversion tool for flat and isotropic layered structures in the ocean (SEIS-FILO)

*Takeshi Akuhara1 (1.Earthquake Research Institute, The University of Tokyo)

Keywords:Transdimensional inversion, Receiver function, Surface wave , Ocean-bottom seismometer, Software

Seismological studies have approximated the Earth’s structure using the stratification of flat and isotropic layers. Such one-dimensional models are suitable for single-station measurements, including receiver functions, surface wave dispersion curves, and H/V ratios, and have been widely used to illuminate a wide range of structures from shallow basin to deep mantle structures. Recently a growing number of studies come to use more sophisticated methods (e.g., full-waveform inversion) for modeling laterally heterogeneous structures of the Earth. Nevertheless, the modeling by one-dimensional structure will not lose its position because of the less computational cost and no need for dense seismic instrument deployments. Especially for the ocean area, insufficient coverage by station networks will force to use such single-station methods.

Here, I present a computer program package termed SEIS-FILO (seismic inversion tool for flat and isotropic layered structures in the ocean). The software performs transdimensional inversion of receiver functions and surface wave dispersion curves, with the effects of the seawater column fully considered. The state-of-the-art Bayesian probabilistic sampling method known as reversible-jump Markov-chain Monte Carlo is used, where the number of model parameters (i.e., the number of layers) is automatically selected based on the data complexity. The parallel tempering technique is also implemented to enhance the conversion rate. The source codes are written in a modern Fortran programming style (i.e., objective-oriented programming) for a low-maintenance and easy expansion in the future.

Currently, the software undergoes alpha testing: independent/joint inversions of receiver functions and fundamental Rayleigh wave dispersion curves are possible to perform. This alpha version is available at the GitHub repository (