Scalable Earthquake Simulation on Petascale Supercomputers

Yifeng Cui, Kim B. Olsen, Thomas H. Jordan, Kwangyoon Lee, Jun Zhou, Patrick Small, Daniel Roten, Geoffrey P. Ely, Dhabaleswar K. Panda, Amit Chourasia, John Levesque, Steven M. Day, & Philip J. Maechling

Published 2010, SCEC Contribution #1443

Petascale simulations are needed to understand the rupture and wave dynamics of the largest earthquakes at shaking frequencies required to engineer safe structures (> 1 Hz). Toward this goal, we have developed a highly scalable, parallel application (AWP-ODC) that has achieved “M8”: a full dynamical simulation of a magnitude-8 earthquake on the southern San Andreas fault up to 2 Hz. M8 was calculated using a uniform mesh of 436 billion 40-m3 cubes to represent the three-dimensional crustal structure of Southern California, in a 800 km by 400 km area, home to over 20 mil-lion people. This production run producing 360 sec of wave propagation sustained 220 Tflop/s for 24 hours on NCCS Jaguar using 223,074 cores. As the largest-ever earthquake simulation, M8 opens new territory for earthquake science and engineering—the physics-based modeling of the largest seismic hazards with the goal of reducing their potential for loss of life and property.

Citation
Cui, Y., Olsen, K. B., Jordan, T. H., Lee, K., Zhou, J., Small, P., Roten, D., Ely, G. P., Panda, D. K., Chourasia, A., Levesque, J., Day, S. M., & Maechling, P. J. (2010). Scalable Earthquake Simulation on Petascale Supercomputers. Oral Presentation at SC10 International Conference for HPC, Networking and Analysis. doi: 10.1109/SC.2010.45.