The SCEC/USGS Dynamic Earthquake Rupture Code Verification Exercise

Harris, Ruth; Barall, Michael; Archuleta, Ralph J.; Dunham, Eric M.; Aagaard, Brad T.; Ampuero, Jean‐Paul; Bhat, Harsha S.; Cruz‐Atienza, V. M.; Dalguer, L. A.; Dawson, Phillip B.; Day, Steven M.; Duan, Benchun; Ely, G.; Kaneko, Yoshihiro; Kase, Yoshihiro; Lapusta, N.; Liu, Yajing; Ma, Shuo; Oglesby, D. D.; Olsen, Kim B.; Pitarka, Arben; Song, Seok Goo; Templeton, E. L.

doi:10.1785/gssrl.80.1.119

Cited by 240 publications

(204 citation statements)

References 11 publications

Supporting

Mentioning

201

Contrasting

Order By: Relevance

“…However, this approximation works quite well for strike-slip faults, which slip mostly in the along-strike direction x. In the code comparison exercise organized by the Southern California Earthquake Center (SCEC), simulations of dynamic rupture on a slip-weakening strike-slip fault in an elastic half-space have been compared for different numerical methods [Harris et al, 2004[Harris et al, , 2008. The comparison of our approach with other methods that can represent the true traction-free surface showed that this approximation captures most effects of the free surface and that the errors induced are relatively small and restricted to the region right next to the free surface.…”

Section: Model Of a Strike-slip Faultmentioning

confidence: 99%

“…Our approach builds on the model of Lapusta et al [2000], with a number of modifications required in three dimensions such as a different truncation procedure, and uses a boundary integral method. The dynamic part of our 3-D methodology has been validated through the Southern California Earthquake Center (SCEC) code comparison exercise and additional studies [Day et al, 2005;Harris et al, 2008].…”

Section: Introductionmentioning

confidence: 99%

See 1 more Smart Citation

Three‐dimensional boundary integral modeling of spontaneous earthquake sequences and aseismic slip

2009

View full text Add to dashboard Cite

[1] Fault processes involve complex patterns of seismic events and aseismic slip. This work develops a three-dimensional (3-D) methodology for simulating long-term history of spontaneous seismic and aseismic slip on a vertical planar strike-slip fault subjected to slow tectonic loading. Our approach reproduces all stages of earthquake cycles, from accelerating slip before dynamic instability, to rapid dynamic propagation of earthquake rupture, to postseismic slip, and to interseismic creep, including aseismic transients. We use the developed 3-D methodology to study interaction of fault slip with a small patch of higher normal stress over long-term slip history. For uniform initial prestress, dynamic rupture is significantly affected by the stronger patch in the first simulated event but not in subsequent ones. The change in behavior is due to redistribution of shear stress by prior slip, which demonstrates that distributions of fault strength and stress are related and illustrates the importance of simulating long slip histories even in studies of dynamic rupture. Despite no long-term effect on dynamic rupture, the small patch of higher normal stress influences nucleation processes and hence long-term slip patterns in the model. Comparison of the fully dynamic simulations and a widely used quasi-dynamic approach shows that the quasi-dynamic approach modifies long-term slip patterns in addition to resulting in much smaller slip velocities and rupture speeds during dynamic events. We show that the response of quasi-dynamic formulations with reduced radiation damping terms can be scaled to match the results of the standard quasi-dynamic formulation and hence cannot improve the comparison.Citation: Lapusta, N., and Y. Liu (2009), Three-dimensional boundary integral modeling of spontaneous earthquake sequences and aseismic slip,

show abstract

Section: Model Of a Strike-slip Faultmentioning

confidence: 99%

Section: Introductionmentioning

confidence: 99%

Three‐dimensional boundary integral modeling of spontaneous earthquake sequences and aseismic slip

2009

View full text Add to dashboard Cite

show abstract

“…The method implements Newmark damping (Hughes, 2000) and an optional thin viscous layer surrounding the fault zone (Day et al, 2005;Dalguer and Day, 2007) to suppress spurious highfrequency oscillations. For specific details about damping tuning strategies we refer to Rojas et al (2008) andBarall (2009). In all benchmarks presented here, FaultMod applies a grid-doubling technique to enable high resolution at the fault and its immediate surroundings, and coarser grid spacing away from the fault.…”

Section: Methodsmentioning

confidence: 99%

“…In contrast to the well verified and validated simulation of seismic wave propagation, the verification process of spontaneous dynamic rupture simulations suffers from the lack of analytical reference solutions. Therefore, we verify the performance of the ADER-DG method in the benchmark suite established by the Southern California Earthquake Center (SCEC) (Harris et al, 2009(Harris et al, , 2011. All simulation results presented here are available at http://scecdata.usc.edu/cvws/.…”

Section: Pelties Et Al: Dynamic Rupture Verification Of An Aderdgmentioning

confidence: 99%

Verification of an ADER-DG method for complex dynamic rupture problems

Pelties¹,

Gabriel²,

Ampuero

2014

Geosci. Model Dev.

View full text Add to dashboard Cite

Abstract. We present results of thorough benchmarking of an arbitrary high-order derivative discontinuous Galerkin (ADER-DG) method on unstructured meshes for advanced earthquake dynamic rupture problems. We verify the method by comparison to well-established numerical methods in a series of verification exercises, including dipping and branching fault geometries, heterogeneous initial conditions, bimaterial interfaces and several rate-and-state friction laws. We show that the combination of meshing flexibility and high-order accuracy of the ADER-DG method makes it a competitive tool to study earthquake dynamics in geometrically complicated setups.

show abstract

“…The overall requirement for a numerical code is to satisfy the three basic properties: i) the consistency of the discretized (algebraic) equations with respect to the original differential equations, ii) the stability and iii) the convergence of the numerical solution. The goodness of the obtained synthetic solution has to be validated through a systematic comparison against other numerical solutions, obtained independently and with different numerical algorithms (e.g., Bizzarri et al, 2001;Harris et al, 2009). Another essential feature of a numerical code is represented by the computation requests (or the computational efficiency), expressed in terms of memory requirements and CPU time.…”

Section: The Numerical Approachmentioning

confidence: 99%

Toward the Formulation of a Realistic Fault Governing Law in Dynamic Models of Earthquake Ruptures

Bizzarri¹

2010

Dynamic Modelling

View full text Add to dashboard Cite

The SCEC/USGS Dynamic Earthquake Rupture Code Verification Exercise

Cited by 240 publications

References 11 publications

Three‐dimensional boundary integral modeling of spontaneous earthquake sequences and aseismic slip

Three‐dimensional boundary integral modeling of spontaneous earthquake sequences and aseismic slip

Verification of an ADER-DG method for complex dynamic rupture problems

Toward the Formulation of a Realistic Fault Governing Law in Dynamic Models of Earthquake Ruptures

Contact Info

Product

Resources

About