Parallel Earthquake Simulations on Large-Scale Multicore Supercomputers

Earthquakes are one of the most destructive natural hazards on our planet Earth. Hugh earthquakes striking offshore may cause devastating tsunamis, as evidenced by the 11 March 2011 Japan (moment magnitude Mw9.0) and the 26 December 2004 Sumatra (Mw9.1) earthquakes. Earthquake prediction (in terms of the precise time, place, and magnitude of a coming earthquake) is arguably unfeasible in the foreseeable future. To mitigate seismic hazards from future earthquakes in earthquake-prone areas, such as California and Japan, scientists have been using numerical simulations to study earthquake rupture propagation along faults and seismic wave propagation in the surrounding media on ever-advancing modern computers over past several decades. In particular, ground motion simulations for past and future (possible) significant earthquakes have been performed to understand factors that affect ground shaking in populated areas, and to provide ground shaking characteristics and synthetic seismograms for emergency preparation and design of earthquake-resistant structures. These simulation results can guide the development of more rational seismic provisions for leading to safer, more efficient, and economical50pt]Please provide V. Taylor author e-mail ID. structures in earthquake-prone regions.

[1]  Philip J. Maechling,et al.  TeraShake2: Spontaneous Rupture Simulations of Mw 7.7 Earthquakes on the Southern San Andreas Fault , 2008 .

[2]  Xingfu Wu,et al.  Parallel Finite Element Earthquake Rupture Simulations on Quad- and Hex-core Cray XT Systems , 2011 .

[3]  Dhabaleswar K. Panda,et al.  Scalable Earthquake Simulation on Petascale Supercomputers , 2010, 2010 ACM/IEEE International Conference for High Performance Computing, Networking, Storage and Analysis.

[4]  G. Mahinthakumar,et al.  A Hybrid Mpi-Openmp Implementation of an Implicit Finite-Element Code on Parallel Architectures , 2002, Int. J. High Perform. Comput. Appl..

[5]  J. Kristek,et al.  The finite-difference and finite-element modeling of seismic wave propagation and earthquake motion , 2007 .

[6]  Jing Zhu,et al.  Toward petascale earthquake simulations , 2009 .

[7]  Robert W. Graves,et al.  The SCEC Southern California Reference Three-Dimensional Seismic Velocity Model Version 2 , 2000 .

[8]  B. Duan.  Role of initial stress rotations in rupture dynamics and ground motion: A case study with implications for the Wenchuan earthquake , 2010 .

[9]  T. Hanks,et al.  Verifying a Computational Method for Predicting Extreme Ground Motion , 2011 .

[10]  Kim B. Olsen,et al.  On the implementation of perfectly matched layers in a three‐dimensional fourth‐order velocity‐stress finite difference scheme , 2003 .

[11]  G. G. Stokes.  "J." , 1890, The New Yale Book of Quotations.

[12]  Mario Cannataro,et al.  Parallel data intensive computing in scientific and commercial applications , 2002, Parallel Comput..

[13]  S. Day,et al.  Sensitivity Study of Physical Limits on Ground Motion at Yucca Mountain , 2010 .

[14]  Kwan-Liu Ma,et al.  Visualizing Very Large-Scale Earthquake Simulations , 2003, ACM/IEEE SC 2003 Conference (SC'03).

[15]  Xingfu Wu,et al.  An OpenMP Approach to Modeling Dynamic Earthquake Rupture Along Geometrically Complex Faults on CMP Systems , 2009, 2009 International Conference on Parallel Processing Workshops.

[16]  Kengo Nakajima.  OpenMP / MPI Hybrid vs. Flat MPI on the Earth Simulator: Parallel Iterative Solvers for Finite Element Method , 2003, ISHPC.

[17]  Xingfu Wu,et al.  Performance Analysis and Optimization of Parallel Scientific Applications on CMP Cluster Systems , 2008, 2008 International Conference on Parallel Processing - Workshops.

[18]  Christopher R. Bradley,et al.  Memory-Efficient Simulation of Anelastic Wave Propagation , 2001 .

[19]  Robert B. Ross,et al.  Using MPI-2: Advanced Features of the Message Passing Interface , 2003, CLUSTER.

[20]  Ralph J. Archuleta,et al.  The Three-Dimensional Dynamics of Dipping Faults , 2000 .

[21]  D. Oglesby,et al.  Erratum: ``Dynamics of dip-slip faulting: Explorations in two dimensions'' , 2000 .

[22]  Yoshiaki Ida,et al.  Cohesive force across the tip of a longitudinal‐shear crack and Griffith's specific surface energy , 1972 .

[23]  S. Day,et al.  Comparison of finite difference and boundary integral solutions to three‐dimensional spontaneous rupture , 2005 .

[24]  T. Tu,et al.  From Mesh Generation to Scientific Visualization: An End-to-End Approach to Parallel Supercomputing , 2006, ACM/IEEE SC 2006 Conference (SC'06).

[25]  A. Pitarka,et al.  The SCEC/USGS Dynamic Earthquake Rupture Code Verification Exercise , 2012 .

[26]  Kenneth W. Hudnut,et al.  The Shakeout scenario , 2008 .

[27]  Jonathan P. Stewart,et al.  Broadband simulations for Mw 7.8 southern San Andreas earthquakes: Ground motion sensitivity to rupture speed , 2008 .

[28]  Michelle R. Hribar,et al.  Balancing Load versus Decreasing Communication: Parameterizing the Tradeoff , 2001, J. Parallel Distributed Comput..

[29]  Xingfu Wu,et al.  Performance characteristics of hybrid MPI/OpenMP implementations of NAS parallel benchmarks SP and BT on large-scale multicore supercomputers , 2011, PERV.

[30]  Monica D. Kohler,et al.  Mantle Heterogeneities and the SCEC Reference Three-Dimensional Seismic Velocity Model Version 3 , 2003 .

[31]  S. Day,et al.  Inelastic strain distribution and seismic radiation from rupture of a fault kink , 2008 .

[32]  S. Day.  Three-dimensional simulation of spontaneous rupture: The effect of nonuniform prestress , 1982, Bulletin of the Seismological Society of America.

[33]  Philip J. Maechling,et al.  ShakeOut‐D: Ground motion estimates using an ensemble of large earthquakes on the southern San Andreas fault with spontaneous rupture propagation , 2009 .

[34]  Chris H. Q. Ding,et al.  An element-based concurrent partitioner for unstructured finite element meshes , 1996, Proceedings of International Conference on Parallel Processing.

[35]  Xingfu Wu,et al.  Performance Analysis and Optimization of Parallel Scientific Applications on CMP Clusters , 2009, Scalable Comput. Pract. Exp..

[36]  Oglesby,et al.  Earthquakes on dipping faults: the effects of broken symmetry , 1998, Science.

[37]  B. Duan,et al.  Heterogeneous fault stresses from previous earthquakes and the effect on dynamics of parallel strike‐slip faults , 2006 .

[38]  W. Marsden.  I and J , 2012 .

[39]  David R. O'Hallaron,et al.  Scalable Parallel Octree Meshing for TeraScale Applications , 2005, ACM/IEEE SC 2005 Conference (SC'05).

[40]  R. G. Stockwell,et al.  An unusual mesospheric bore event observed at high latitudes over Antarctica , 2006 .

[41]  Phillip A. Laplante.  Performance Analysis and Optimization , 2004 .

[42]  A. Ruina.  Slip instability and state variable friction laws , 1983 .

[43]  D. J. Andrews,et al.  Rupture velocity of plane strain shear cracks , 1976 .

[44]  P. Maechling,et al.  Strong shaking in Los Angeles expected from southern San Andreas earthquake , 2006 .

[45]  B. Duan,et al.  Nonuniform prestress from prior earthquakes and the effect on dynamics of branched fault systems , 2007 .

[46]  Nobuki Kame,et al.  Effects of prestress state and rupture velocity on dynamic fault branching , 2002 .

[47]  Hideo Aochi,et al.  Three‐dimensional nonplanar simulation of the 1992 Landers earthquake , 2002 .

[48]  James H. Dieterich,et al.  Modeling of rock friction: 2. Simulation of preseismic slip , 1979 .

[49]  Xingfu Wu,et al.  Parallel Simulations of Dynamic Earthquake Rupture along Geometrically Complex Faults on CMP Systems , 2011 .

[50]  J. Dieterich.  Modeling of rock friction: 1. Experimental results and constitutive equations , 1979 .

[51]  Gideon Juve,et al.  The ShakeOut earthquake scenario: Verification of three simulation sets , 2010 .