FD3D_TSN: A fast and simple code for dynamic rupture simulations with GPU acceleration

We introduce FD3D_TSN, an open-source Fortran code for 3D dynamic earthquake rupture modeling based on the staggered grid 4th order finite difference method employing a regular cubical spatial discretization. Slip weakening and fast velocity weakening rate-and-state fault friction laws are combined with vertical planar fault geometry, orthogonal to a planar free surface. FD3D_TSN demonstrates good agreement with other methods in a range of benchmark exercises of the Southern California Earthquake Center‐U.S. Geological Survey (SCEC/USGS) dynamic rupture code verification project. Efficient GPU acceleration using the OpenACC framework yields a factor of ten speed-up in terms of time to solution compared to a single-core solution for presented hardware (CPU Intel i9-9900K and GPU Nidia RTX 2070). The software is fast and easy-to-use and suited explicitly for data-driven applications requiring a high number of forward simulations, such as dynamic source inversion or probabilistic ground motion modeling. The code is freely available for use by the scientific community and may be incorporated in physics-based earthquake source imaging and seismic hazard assessment, or used for teaching

@article{id2522,
  author = {Premus, Jan and Gallovi\v{c}, Franti\v{s}ek and Hanyk, Ladislav and Gabriel, Alice-Agnes},
  doi = {10.1785/0220190374},
  journal = {Seismological Research Letters},
  language = {en},
  note = {open access: http://geo.mff.cuni.cz/{\textasciitilde}gallovic/abst/Premus.etal.SRL.2020.pdf},
  number = {5},
  pages = {2881{\textendash}2889},
  title = {FD3D\_TSN: A fast and simple code for dynamic rupture simulations with GPU acceleration},
  url = {https://pubs.geoscienceworld.org/ssa/srl/article/doi/10.1785/0220190374/587424/FD3D-TSN-A-Fast-and-Simple-Code-for-Dynamic},
  volume = {91},
  year = {2020},
}

%O Journal Article
%A Premus, Jan
%A Gallovič, František
%A Hanyk, Ladislav
%A Gabriel, Alice-Agnes
%R 10.1785/0220190374
%J Seismological Research Letters
%G en
%O open access: http://geo.mff.cuni.cz/~gallovic/abst/Premus.etal.SRL.2020.pdf
%N 5
%P 2881–2889
%T FD3D_TSN: A fast and simple code for dynamic rupture simulations with GPU acceleration
%U https://pubs.geoscienceworld.org/ssa/srl/article/doi/10.1785/0220190374/587424/FD3D-TSN-A-Fast-and-Simple-Code-for-Dynamic
%V 91
%D 2020