We present a high-resolution simulation of the 2004 Sumatra-Andaman earthquake, including non-linear frictional failure on a megathrustsplay fault system. Our method exploits unstructured meshes capturing the complicated geometries in subduction zones that are crucial to understand large earthquakes and tsunami generation. These up-to-date largest and longest dynamic rupture simulations enable analysis of dynamic source effects on the seafloor displacements.

SeisSol simulates earthquake dynamics by coupling seismic wave propagation and dynamic rupture simulations with high order accuracy on fully adaptive, unstructured meshes. In this paper we present an optimization of SeisSol's I/O implementations to establish a workflow that supports petascale simulations on large unstructured datasets. Our implementations can handle meshes with more than 1 billion cells and 660 billion degrees of reedom. The results show that SeisSol can initialize the mesh structure within 35 seconds on 2048 SuperMUC nodes from our new optimized mesh format. For the wave field output we implemented carefully tuned I/O routines based on HDF5 and MPI-IO. With an aggregation strategy we are able to increase the write bandwidth from 832 MiB/s to 6.7 GiB/s on 2048 SuperMUC nodes.

We present an end-to-end optimization of the innovative Arbitrary high-order DERivative Discontinuous Galerkin (ADER-DG) software SeisSol targeting Intel® Xeon Phi coprocessor platforms, achieving unprecedented earthquake model complexity through coupled simulation of full frictional sliding and seismic wave propagation. SeisSol exploits unstructured meshes to flexibly adapt for complicated geometries in realistic geological models. Seismic wave propagation is solved simultaneously with earthquake faulting in a multiphysical manner leading to a heterogeneous solver structure. Our architecture aware optimizations deliver up to 50% of peak performance, and introduce an efficient compute-communication overlapping scheme shadowing the multiphysics computations. SeisSol delivers near-optimal weak scaling, reaching 8.6 DP-PFLOPS on 8,192 nodes of the Tianhe-2 supercomputer. Our performance model projects reaching 18 -- 20 DP-PFLOPS on the full Tianhe-2 machine. Of special relevance to modern civil engineering needs, our pioneering simulation of the 1992 Landers earthquake shows highly detailed rupture evolution and ground motion at frequencies up to 10 Hz.