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

Pelties, Christian; Gabriel, Alice‐Agnes; Ampuero, Jean‐Paul

doi:10.5194/gmd-7-847-2014

Cited by 61 publications

(49 citation statements)

References 73 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…Duan & Oglesby 2006;Ripperger et al 2007Ripperger et al , 2008Dalguer & Day 2009;Brietzke et al 2009;Pelties et al 2013) or velocity-and state-dependent friction laws with a static strength threshold (e.g. Ampuero & Ben-Zion 2008).…”

mentioning

confidence: 99%

On the initiation of sustained slip-weakening ruptures by localized stresses

Gális

Pelties

Kristek

et al. 2014

Geophysical Journal International

Self Cite

View full text Add to dashboard Cite

S U M M A R YNumerical simulations of dynamic earthquake rupture require an artificial initiation procedure, if they are not integrated in long-term earthquake cycle simulations. A widely applied procedure involves an 'overstressed asperity', a localized region stressed beyond the static frictional strength. The physical properties of the asperity (size, shape and overstress) may significantly impact rupture propagation. In particular, to induce a sustained rupture the asperity size needs to exceed a critical value. Although criteria for estimating the critical nucleation size under linear slip-weakening friction have been proposed for 2-D and 3-D problems based on simplifying assumptions, they do not provide general rules for designing 3-D numerical simulations. We conduct a parametric study to estimate parameters of the asperity that minimize numerical artefacts (e.g. changes of rupture shape and speed, artificial supershear transition, higher slip-rate amplitudes). We examine the critical size of square, circular and elliptical asperities as a function of asperity overstress and background (off-asperity) stress. For a given overstress, we find that asperity area controls rupture initiation while asperity shape is of lesser importance. The critical area obtained from our numerical results contrasts with published theoretical estimates when background stress is low. Therefore, we derive two new theoretical estimates of the critical size under low background stress while also accounting for overstress. Our numerical results suggest that setting the asperity overstress and area close to their critical values eliminates strong numerical artefacts even when the overstress is large. We also find that properly chosen asperity size or overstress may significantly shorten the duration of the initiation. Overall, our results provide guidelines for determining the size of the asperity and overstress to minimize the effects of the forced initiation on the subsequent spontaneous rupture propagation.

show abstract

mentioning

confidence: 99%

On the initiation of sustained slip-weakening ruptures by localized stresses

Gális

Pelties

Kristek

et al. 2014

Geophysical Journal International

Self Cite

View full text Add to dashboard Cite

show abstract

“…To discretize in time, implicit methods have already been employed on a few occasions (Pelties et al 2014) and in addition to adaptivity (Liu & Rice 2005), a physically motivated step size control has seen repeated use (Lapusta et al 2000;Lapusta & Liu 2009;Kaneko et al 2010;Barbot et al 2012).…”

Section: Introductionmentioning

confidence: 99%

“…Past developments involving rate-and-state friction include the use of boundary element methods which are well-suited for homogeneous problems (Perrin et al 1995;Lapusta et al 2000;Lapusta & Liu 2009;Kaneko et al 2010;Barbot et al 2012), spectral elements (Kaneko et al 2008) and discontinuous Galerkin for a regularized model (Pelties et al 2014) as spatial discretization schemes. To discretize in time, implicit methods have already been employed on a few occasions (Pelties et al 2014) and in addition to adaptivity (Liu & Rice 2005), a physically motivated step size control has seen repeated use (Lapusta et al 2000;Lapusta & Liu 2009;Kaneko et al 2010;Barbot et al 2012).…”

Section: Introductionmentioning

confidence: 99%

See 1 more Smart Citation

On the efficient and reliable numerical solution of rate-and-state friction problems

et al. 2016

View full text Add to dashboard Cite

“…Our dynamic source model incorporates new degree of realism by integrating a comprehensive set of geological and geophysical information such as high-resolution topography, rotating tectonic stresses, 3D velocity structure, depth-dependent bulk cohesion, and a complex intersecting fault geometry. Unifying aforementioned complexities is enabled by using SeisSol (www.seissol.org, Dumbser and ; Pelties et al [2014]), a software package specifically suited for handling complex geometries and for the efficient use on modern high-performance computing infrastructure [e.g., Heinecke et al, 2014;Uphoff et al, 2017]. This work extends recent models presented in Wollherr et al, 2018] which included complex fault geometries and off-fault plasticity but were restricted to 1D velocity structure, constantly oriented tectonic background stress and neglecting viscoelastic attenuation of the seismic wave field.…”

Section: Introductionmentioning

confidence: 99%

Landers 1992 "reloaded": Integrative dynamic earthquake rupture modeling

Gabriel

Wollherr

Schliwa

et al. 2018

Preprint

Self Cite

View full text Add to dashboard Cite

Key Points:• A new integrative physics-based simulation of the 1992 Landers earthquake reproduces a broad range of independent observation • Sustained dynamic rupture interconnecting the complex fault segments constraints pre-stress and fault strength • We discuss the interplay of rupture transfers, geometric fault complexity, initial stress, viscoelastic attenuation, and off-fault plasticity AbstractThe 1992 M w 7.3 Landers earthquake is perhaps one of the best studied seismic events. However, many aspects of the dynamics of the rupture process are still puzzling, e.g. how did rupture transfer between fault segments? We present 3D spontaneous dynamic rupture simulations of a new degree of realism, incorporating the interplay of fault geometry, topography, 3D rheology, off-fault plasticity and viscoelastic attenuation. The surprisingly unique scenario reproduces a broad range of observations, including final slip distribution, seismic moment-rate function, seismic waveform characteristics and peak ground velocities, as well as shallow slip deficits and mapped off-fault deformation patterns. Sustained dynamic rupture of all fault segments in general, and rupture transfers in particular, put strong constraints on amplitude and orientation of initial fault stresses and friction. Source dynamics include dynamic triggering over large distances and direct branching; rupture terminates spontaneously on most of the principal fault segments. We achieve good agreement between synthetic and observed waveform characteristics and associated peak ground velocities. Despite very complex rupture evolution, ground motion variability is close to what is commonly assumed in Ground Motion Prediction Equations. We examine the effects of variations in modeling parameterization, e.g. purely elastic setups or models neglecting viscoelastic attenuation, in comparison to our preferred model. Our integrative dynamic modeling approach demonstrates the potential of consistent in-scale earthquake rupture simulations for augmenting earthquake source observations and improving the understanding of earthquake source physics of complex, segmented fault systems.

show abstract

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

Cited by 61 publications

References 73 publications

On the initiation of sustained slip-weakening ruptures by localized stresses

On the initiation of sustained slip-weakening ruptures by localized stresses

On the efficient and reliable numerical solution of rate-and-state friction problems

Landers 1992 "reloaded": Integrative dynamic earthquake rupture modeling

Contact Info

Product

Resources

About