Bayesian Dynamic Finite-Fault Inversion: 1. Method and Synthetic Test

Dynamic earthquake source inversions aim to determine the spatial distribution of initial stress and friction parameters leading to dynamic rupture models that reproduce observed ground motion data. Such inversions are challenging, particularly due to their high computational burden, thus so far only few attempts have been made. Here we introduce a novel method providing a representative sample of acceptable dynamic models at an affordable computational cost, from which dynamic source parameters and their uncertainties can be assessed. The method assumes a linear slip-weakening friction law and spatially variable prestress, strength and critical slip weakening distance along the fault. The inverse problem is formulated in a Bayesian framework and the posterior probability density function is sampled using the Parallel Tempering Monte Carlo algorithm. We use a highly efficient forward solver combining 3D finite difference code for dynamic rupture simulation on a simplified geometry to compute slip rates, and pre-calculated Green’s functions to compute ground motions. We demonstrate the performance of the proposed method on a community benchmark test for source inversion. We find that the dynamic parameters are resolved well within the uncertainty, especially in large slip areas. The overall relative uncertainty of the dynamic parameters is rather large, reaching ~50% of the averaged values. In contrast, the kinematic rupture parameters (rupture times, rise times, slip values), also well-resolved, have relatively lower uncertainties of ~10%. We conclude that incorporating physics-based constraints, such as an adequate friction law, may serve also as an effective constraint on the rupture kinematics in finite-fault inversions.

@article{id2447,
  author = {Gallovic, Frantisek and Valentova, Lubica and Ampuero, Jean-Paul and Gabriel, Alice-Agnes},
  doi = {10.1029/2019JB01751},
  journal = {Journal of Geophysical Research: Solid Earth},
  language = {en},
  title = {Bayesian Dynamic Finite-Fault Inversion: 1. Method and Synthetic Test},
  url = {https://eartharxiv.org/tmjv4/},
  year = {2019},
}

%O Journal Article
%A Gallovic, Frantisek
%A Valentova, Lubica
%A Ampuero, Jean-Paul
%A Gabriel, Alice-Agnes
%R 10.1029/2019JB01751
%J Journal of Geophysical Research: Solid Earth
%G en
%T Bayesian Dynamic Finite-Fault Inversion: 1. Method and Synthetic Test
%U https://eartharxiv.org/tmjv4/
%D 2019