FASTDASH: An Implementation of 3D Earthquake Cycle Simulation on Complex Fault Systems Using the Boundary Element Method Accelerated by H-matrices

Fuente: arXiv
Salvato in:
Dettagli Bibliografici
Autori principali: Cheng, Jinhui, Bhat, Harsha S., Almakari, Michelle, Lecampion, Brice, Peruzzo, Carlo
Natura: Preprint
Pubblicazione: 2025
Soggetti:
Accesso online:
Tags: Aggiungi Tag
Nessun Tag, puoi essere il primo ad aggiungerne!!
_version_ 1866909622953050112
author Cheng, Jinhui
Bhat, Harsha S.
Almakari, Michelle
Lecampion, Brice
Peruzzo, Carlo
author_facet Cheng, Jinhui
Bhat, Harsha S.
Almakari, Michelle
Lecampion, Brice
Peruzzo, Carlo
contents Fault systems have geometrically complex structures in nature, such as stepovers, bends, branches, and roughness. Many geological and geophysical studies have shown that the geometrical complexity of fault systems in nature decisively influences the initiation, arrest, and recurrence of seismic and aseismic events. However, a vast majority of models of slip dynamics are conducted on planar faults due to algorithmic limitations. We develop a 3D quasi-dynamic slip dynamics model to overcome this restriction. The calculation of the elastic response due to slip is a matrix-vector multiplication in boundary element method, which can be accelerated by using hierarchical matrices. The computational complexity is reduced from the order of O(N^2) to O(N log N), where N is the number of degree of freedoms used. We validate our code with a static crack analytical solution and the SEAS benchmark/validation exercise from Southern California Earthquake Center. We further employ this method on a realistic fault system with complex geometry that was reactivated during the 2023 Kahramanmaras - Turkiye doublet earthquakes, generating slip sequences that closely match real observations.
format Preprint
id arxiv_https___arxiv_org_abs_2505_02398
institution arXiv
publishDate 2025
record_format arxiv
spellingShingle FASTDASH: An Implementation of 3D Earthquake Cycle Simulation on Complex Fault Systems Using the Boundary Element Method Accelerated by H-matrices
Cheng, Jinhui
Bhat, Harsha S.
Almakari, Michelle
Lecampion, Brice
Peruzzo, Carlo
Geophysics
Fault systems have geometrically complex structures in nature, such as stepovers, bends, branches, and roughness. Many geological and geophysical studies have shown that the geometrical complexity of fault systems in nature decisively influences the initiation, arrest, and recurrence of seismic and aseismic events. However, a vast majority of models of slip dynamics are conducted on planar faults due to algorithmic limitations. We develop a 3D quasi-dynamic slip dynamics model to overcome this restriction. The calculation of the elastic response due to slip is a matrix-vector multiplication in boundary element method, which can be accelerated by using hierarchical matrices. The computational complexity is reduced from the order of O(N^2) to O(N log N), where N is the number of degree of freedoms used. We validate our code with a static crack analytical solution and the SEAS benchmark/validation exercise from Southern California Earthquake Center. We further employ this method on a realistic fault system with complex geometry that was reactivated during the 2023 Kahramanmaras - Turkiye doublet earthquakes, generating slip sequences that closely match real observations.
title FASTDASH: An Implementation of 3D Earthquake Cycle Simulation on Complex Fault Systems Using the Boundary Element Method Accelerated by H-matrices
topic Geophysics
url https://arxiv.org/abs/2505.02398