A Caltech Library Service

Community-Driven Code Comparisons for Three-Dimensional Dynamic Modeling of Sequences of Earthquakes and Aseismic Slip

Jiang, Junle and Erickson, Brittany A. and Lambert, Valère R. and Ampuero, Jean-Paul and Ando, Ryosuke and Barbot, Sylvain D. and Cattania, Camilla and Dal Zilio, Luca and Duan, Benchun and Dunham, Eric M. and Gabriel, Alice-Agnes and Lapusta, Nadia and Li, Duo and Li, Meng and Liu, Dunyu and Liu, Yajing and Ozawa, So and Pranger, Casper and van Dinther, Ylona (2022) Community-Driven Code Comparisons for Three-Dimensional Dynamic Modeling of Sequences of Earthquakes and Aseismic Slip. Journal of Geophysical Research. Solid Earth, 127 (3). Art. No. e2021JB023519. ISSN 2169-9313. doi:10.1029/2021JB023519.

[img] PDF - Published Version
Creative Commons Attribution Non-commercial.

[img] PDF (3 Feb 2022) - Submitted Version
Creative Commons Attribution Non-commercial No Derivatives.

[img] PDF (Data Set S1) - Supplemental Material
Creative Commons Attribution Non-commercial.

[img] PDF (Data Set S2) - Supplemental Material
Creative Commons Attribution Non-commercial.


Use this Persistent URL to link to this item:


Dynamic modeling of sequences of earthquakes and aseismic slip (SEAS) provides a self-consistent, physics-based framework to connect, interpret, and predict diverse geophysical observations across spatial and temporal scales. Amid growing applications of SEAS models, numerical code verification is essential to ensure reliable simulation results but is often infeasible due to the lack of analytical solutions. Here, we develop two benchmarks for three-dimensional (3D) SEAS problems to compare and verify numerical codes based on boundary-element, finite-element, and finite-difference methods, in a community initiative. Our benchmarks consider a planar vertical strike-slip fault obeying a rate- and state-dependent friction law, in a 3D homogeneous, linear elastic whole-space or half-space, where spontaneous earthquakes and slow slip arise due to tectonic-like loading. We use a suite of quasi-dynamic simulations from 10 modeling groups to assess the agreement during all phases of multiple seismic cycles. We find excellent quantitative agreement among simulated outputs for sufficiently large model domains and small grid spacings. However, discrepancies in rupture fronts of the initial event are influenced by the free surface and various computational factors. The recurrence intervals and nucleation phase of later earthquakes are particularly sensitive to numerical resolution and domain-size-dependent loading. Despite such variability, key properties of individual earthquakes, including rupture style, duration, total slip, peak slip rate, and stress drop, are comparable among even marginally resolved simulations. Our benchmark efforts offer a community-based example to improve numerical simulations and reveal sensitivities of model observables, which are important for advancing SEAS models to better understand earthquake system dynamics.

Item Type:Article
Related URLs:
URLURL TypeDescription Paper ItemBP4 and BP5 ItemSEAS ItemGARNET
Jiang, Junle0000-0002-8796-5846
Erickson, Brittany A.0000-0001-9457-8572
Lambert, Valère R.0000-0002-6174-9651
Ampuero, Jean-Paul0000-0002-4827-7987
Ando, Ryosuke0000-0002-6205-3699
Barbot, Sylvain D.0000-0003-4257-7409
Cattania, Camilla0000-0003-0031-1696
Dal Zilio, Luca0000-0002-5642-0894
Duan, Benchun0000-0003-2622-3811
Dunham, Eric M.0000-0003-0804-7746
Gabriel, Alice-Agnes0000-0003-0112-8412
Lapusta, Nadia0000-0001-6558-0323
Li, Duo0000-0002-8641-337X
Li, Meng0000-0001-9459-7167
Liu, Dunyu0000-0001-6740-3080
Liu, Yajing0000-0002-5323-8077
Ozawa, So0000-0002-6507-6925
Pranger, Casper0000-0002-8943-896X
van Dinther, Ylona0000-0002-4273-8287
Alternate Title:Community-Driven Code Comparisons for Three-Dimensional Dynamic Modeling of Sequences of Earthquakes and Aseismic Slip (SEAS)
Additional Information:© 2022 The Authors. This is an open access article under the terms of the Creative Commons Attribution-NonCommercial License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited and is not used for commercial purposes. Issue Online: 10 March 2022; Version of Record online: 10 March 2022; Accepted manuscript online: 03 March 2022; Manuscript accepted: 28 February 2022; Manuscript revised: 04 February 2022; Manuscript received: 01 November 2021. We thank Michael Barall for maintaining the SEAS online platform and Ruth Harris for providing experience from the SCEC/USGS code verification exercises for dynamic rupture simulations. We thank Steve Day, Ruth Harris, Associate Editor, and Editor Rachel Abercrombie for reviewing this manuscript. J. Jiang and B. A. Erickson designed the benchmark problems and organized the workshops for code verification exercises. J. Jiang analyzed all simulation results and led the writing of the manuscript. All remaining authors provided feedback on benchmark design, participated in the benchmark exercises (listed in Table 2), and/or helped with revising the manuscript. V. Lambert additionally helped with the early tests of benchmark problems; the other authors are listed alphabetically. The work of J. Jiang and B. A. Erickson on this project and SEAS-themed workshops have been supported by the Southern California Earthquake Center (SCEC). This research is SCEC Contribution No. 11680. SCEC is funded by National Science Foundation (NSF) Cooperative Agreement EAR-1600087 and United States Geological Survey (USGS) Cooperative Agreement G17AC00047. The simulations with BICyclE (J. Jiang) were conducted with the award EAR170014 from the Extreme Science and Engineering Discovery Environment (XSEDE), which is supported by NSF Grant No. ACI-1548562. The simulations with BICyclE (V. Lambert and N. Lapusta) were carried out on the High Performance Computing Center cluster of the California Institute of Technology with the grant support from SCEC. The simulations with Unicycle and Motorcycle are supported in part by the NSF Grant No. EAR-1848192. The simulations with ESAM are supported by Natural Sciences and Engineering Research Council of Canada (NSERC) Discovery Grant RGPIN-2018-05389. The simulations with HBI are supported in part by Japan Society for the Promotion of Science (JSPS) KAKENHI Grant No. 19K04031.The simulations with GARNET are supported by the Dutch Research Council (NWO) grant DEEP.NL.2018.037, Swiss National Science Foundation (SNSF) Grant 200021-169880, EU-MC ITN Grant 604713, and EU ERC StG 852992. The simulations with TriBIE are supported by the European Union's Horizon 2020 research and innovation program (TEAR ERC Starting Grant No. 852992) and NSF Grant No. EAR-2121666. The simulations with EQsimu were conducted with the advanced computing resources provided by Texas A&M High Performance Research Computing and by the Texas Advanced Computing Center (TACC) at the University of Texas at Austin. Data Availability Statement: Descriptions of benchmarks BP4 and BP5 are available at (SEAS_BP4_QD.pdf and SEAS_BP5.pdf) and included as supplementary materials. The simulation data analyzed in this paper are available at The SEAS online platform ( hosts the time series and rupture time data for the full suite of simulations. See publications in Table 2 for the availability and repositories of numerical codes. GARNET is available at
Group:Center for Geomechanics and Mitigation of Geohazards (GMG), Division of Geological and Planetary Sciences, Seismological Laboratory
Funding AgencyGrant Number
Southern California Earthquake Center (SCEC)UNSPECIFIED
Natural Sciences and Engineering Research Council of Canada (NSERC)RGPIN-2018-05389
Japan Society for the Promotion of Science (JSPS)19K04031
Nederlandse Organisatie voor Wetenschappelijk Onderzoek (NWO)DEEP.NL.2018.037
Swiss National Science Foundation (SNSF)200021-169880
European Research Council (ERC)604713
European Research Council (ERC)852992
Subject Keywords:earthquake source processes; aseismic slip; crustal faulting; numerical modeling; code verification; community benchmarks
Other Numbering System:
Other Numbering System NameOther Numbering System ID
Southern California Earthquake Center11680
Issue or Number:3
Record Number:CaltechAUTHORS:20211206-184347093
Persistent URL:
Official Citation:Jiang, J., Erickson, B. A., Lambert, V. R., Ampuero, J.-P., Ando, R., Barbot, S. D., et al. (2022). Community-driven code comparisons for three-dimensional dynamic modeling of sequences of earthquakes and aseismic slip. Journal of Geophysical Research: Solid Earth, 127, e2021JB023519.
Usage Policy:No commercial reproduction, distribution, display or performance rights in this work are provided.
ID Code:112228
Deposited By: Tony Diaz
Deposited On:07 Dec 2021 17:49
Last Modified:15 Nov 2022 19:02

Repository Staff Only: item control page