Petascale High Order Dynamic Rupture Earthquake Simulations on Heterogeneous Supercomputers

Heinecke, Alexander; Breuer, Alexander; Rettenberger, Sebastian; Bäder, Michael; Gabriel, Alice‐Agnes; Pelties, Christian; Bode, Arndt; Barth, William L.; Liao, Xiangke; Vaidyanathan, Karthikeyan; Smelyanskiy, Mikhail; Dubey, Pradeep

doi:10.1109/sc.2014.6

Cited by 111 publications

(95 citation statements)

References 35 publications

Supporting

Mentioning

Contrasting

Unclassified

Order By: Relevance

“…• , consistent with previous static and dynamic modeling studies of the full or southern-central fault system [Madden et al, 2013;Heinecke et al, 2014;Wollherr et al, 2018].…”

Section: On-fault Initial Stressessupporting

confidence: 88%

“…Harris et al, 2018]. End-to-end optimization Heinecke et al, 2014;Breuer et al, 2015Breuer et al, , 2016Rettenberger and Bader , 2015;Rettenberger et al, 2016] targeting high efficiency on high-performance computing infrastructure includes a ten-fold speedup by an efficient local time-stepping algorithm [Uphoff et al, 2017]. Viscoelastic rheologies are incorporated using an offline code-generator to compute matrix products in a computationally highly efficient way.…”

Section: Methodsmentioning

confidence: 99%

“…Peyrat et al, 2001;Aochi and Fukuyama, 2002;Fliss et al, 2005;Heinecke et al, 2014;Wollherr et al, 2018]. Together with detailed analysis of the recorded strong ground motions [e.g., Campbell and Bozorgnia, 1994], rupture transfer mechanisms [e.g., Wesnousky, 2006;Madden and Pollard , 2012;Madden et al, 2013] and potential energy release [e.g., Dreger , 1994;Wald and Heaton, 1994] a comprehensive picture of the source kinematics and macroscopic earthquake properties has been developed.…”

Section: Introductionmentioning

confidence: 99%

“…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%

See 3 more Smart Citations

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

“…• , consistent with previous static and dynamic modeling studies of the full or southern-central fault system [Madden et al, 2013;Heinecke et al, 2014;Wollherr et al, 2018].…”

Section: On-fault Initial Stressessupporting

confidence: 88%

Section: Methodsmentioning

confidence: 99%

Section: Introductionmentioning

confidence: 99%

Section: Introductionmentioning

confidence: 99%

See 2 more Smart Citations

Landers 1992 "reloaded": Integrative dynamic earthquake rupture modeling

Gabriel

Wollherr

Schliwa

et al. 2018

Preprint

Self Cite

View full text Add to dashboard Cite

show abstract

“…Some state-of-the-art works have been presented at the International Conference on High Performance Computing, Networking, Storage and Analysis (SC) conference series (cf. Chhugani et al, 2012;BermejoMoreno et al, 2013;Heinecke et al, 2014;Rudi et al, 2015) and special extreme scaling workshops (cf. Brömmel et al, 2015.…”

Section: Introductionmentioning

confidence: 99%

Asynchronous communication in spectral-element and discontinuous Galerkin methods for atmospheric dynamics – a case study using the High-Order Methods Modeling Environment (HOMME-homme_dg_branch)

Jamroz¹,

Klöfkorn²

2016

Geosci. Model Dev.

View full text Add to dashboard Cite

Abstract. The scalability of computational applications on current and next-generation supercomputers is increasingly limited by the cost of inter-process communication. We implement non-blocking asynchronous communication in the High-Order Methods Modeling Environment for the time integration of the hydrostatic fluid equations using both the spectral-element and discontinuous Galerkin methods. This allows the overlap of computation with communication, effectively hiding some of the costs of communication. A novel detail about our approach is that it provides some data movement to be performed during the asynchronous communication even in the absence of other computations. This method produces significant performance and scalability gains in large-scale simulations.

show abstract

Off‐fault deformations and shallow slip deficit from dynamic rupture simulations with fault zone plasticity

Roten

Olsen

Day

2017

Geophysical Research Letters

View full text Add to dashboard Cite

Kinematic source inversions of major (M≥7) strike‐slip earthquakes show that the slip at depth exceeds surface displacements measured in the field, and it has been suggested that this shallow slip deficit (SSD) is caused by distributed plastic deformation near the surface. We perform dynamic rupture simulations of M 7.2–7.4 earthquakes in elastoplastic media and analyze the sensitivity of SSD and off‐fault deformation (OFD) to rock quality parameters. While linear simulations clearly underpredict observed SSD and OFDs, nonlinear simulations for a moderately fractured fault damage zone predict a SSD of 44–53% and OFDs of 39–48%, consistent with the 30–60% SSD and 46 ± 10% (1σ) OFD reported for the 1992 M 7.3 Landers earthquake. Both SSD and OFDs are sensitive to the quality of the fractured rock mass inside the fault damage zone, and surface rupture is almost entirely suppressed in poor quality material.

show abstract

Petascale High Order Dynamic Rupture Earthquake Simulations on Heterogeneous Supercomputers

Cited by 111 publications

References 35 publications

Landers 1992 "reloaded": Integrative dynamic earthquake rupture modeling

Landers 1992 "reloaded": Integrative dynamic earthquake rupture modeling

Asynchronous communication in spectral-element and discontinuous Galerkin methods for atmospheric dynamics – a case study using the High-Order Methods Modeling Environment (HOMME-homme_dg_branch)

Off‐fault deformations and shallow slip deficit from dynamic rupture simulations with fault zone plasticity

Contact Info

Product

Resources

About