3D Parallel Elastodynamic Modeling of Large Subduction Earthquakes

Cabrera, Eduardo; Chávez, M.; Madariaga, Raúl; Perea, N.; Frisenda, M.

doi:10.1007/978-3-540-75416-9_50

Cited by 2 publications

(16 citation statements)

References 4 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…The 3D velocity-stress form of the elastic wave equation, consists of nine coupled, first order partial differential hyperbolic equations for the three particle velocity vector components and the six independent stress tensor components [1,5].…”

Section: D Wave Propagation Modeling and Its Algorithmmentioning

confidence: 99%

“…The domain was decomposed into small subdomains and distributed among a number of processors, using simple partitioning to give an equal number of grid points to each processor [1]. This approach is appropriate for the 3DFD wave propagation code, as large problems are too big to fit on a single processor [1].…”

Section: Parallel Implementation Of the 3dfd Algorithmmentioning

confidence: 99%

“…This approach is appropriate for the 3DFD wave propagation code, as large problems are too big to fit on a single processor [1]. The Message Passing Interface (MPI) was used to parallelize the 3DFD code [1]. In particular, MPI_Bcast, MPI_Cart_Shift and MPI_SendRecv instructions were used; the first two to communicate the geometry and physical properties of the problem, before starting the wave propagation loop, and the last to update the velocities and stresses calculated at each time step.…”

Section: Parallel Implementation Of the 3dfd Algorithmmentioning

confidence: 99%

“…Parallel I/O from MPI-2 was used in the code to read the earth model data by all processors and to write the velocity seismograms by the processors corresponding to the free surface of the physical domain [1], which is only a small percentage of the total number of processors. As this type of parallel I/O is machine independent, it fitted the benchmark experiment performed on the three platforms.…”

Section: Parallel Implementation Of the 3dfd Algorithmmentioning

confidence: 99%

“…1), poses both a numerical and a computational challenge, particularly because it requires enormous amounts of memory and storage, as well as an intensive use of computing resources. As the recurrence time estimated for this highly destructive type of event in Mexico is only a few decades, there is a seismological, engineering and socioeconomical interest in modeling them by using parallel computing [1]. In this paper, we present the results from benchmark studies performed on a recently optimized parallel 3D wave propagation staggered-grid finite difference code, using the Message Passing Interface (MPI).…”

Section: Introductionmentioning

confidence: 99%

See 4 more Smart Citations

Benchmark Study of a 3d Parallel Code for the Propagation of Large Subduction Earthquakes

Chávez

Cabrera

Madariaga

et al. 2008

Recent Advances in Parallel Virtual Machine and Message Passing Interface

Self Cite

View full text Add to dashboard Cite

Abstract. Benchmark studies were carried out on a recently optimized parallel 3D seismic wave propagation code that uses finite differences on a staggered grid with 2 nd order operators in time and 4 th order in space. Three dual-core supercomputer platforms were used to run the parallel program using MPI. Efficiencies of 0.91 and 0.48 with 1024 cores were obtained on HECToR (UK) and KanBalam (Mexico), and 0.66 with 8192 cores on HECToR. The 3D velocity field pattern from a simulation of the 1985 Mexico earthquake (that caused the loss of up to 30000 people and about 7 billion US dollars) which has reasonable agreement with the available observations, shows coherent, well developed surface waves propagating towards Mexico City.

show abstract

Section: D Wave Propagation Modeling and Its Algorithmmentioning

confidence: 99%

Section: Parallel Implementation Of the 3dfd Algorithmmentioning

confidence: 99%

Section: Parallel Implementation Of the 3dfd Algorithmmentioning

confidence: 99%

Section: Parallel Implementation Of the 3dfd Algorithmmentioning

confidence: 99%

Section: Introductionmentioning

confidence: 99%

See 3 more Smart Citations

Benchmark Study of a 3d Parallel Code for the Propagation of Large Subduction Earthquakes

Chávez

Cabrera

Madariaga

et al. 2008

Recent Advances in Parallel Virtual Machine and Message Passing Interface

Self Cite

View full text Add to dashboard Cite

show abstract

Low-Frequency 3D Wave Propagation Modeling of the 12 May 2008 Mw 7.9 Wenchuan Earthquake

Chávez¹,

Cabrera²,

Madariaga³

et al. 2010

Bulletin of the Seismological Society of America

View full text Add to dashboard Cite

The seismic potential of southern China is associated with the collision between the Indian and the Eurasian plates. This is manifested in the western Sichuan Plateau by several seismically active systems of faults, such as the Longmen Shan. The seismicity observed on the Longmen Shan fault includes recent events with magnitudes of up to 6.5, and the one of 12 May 2008 M w 7.9 Wenchuan earthquake. Herewith, as part of an ongoing research program, a recently optimized threedimensional (3D) seismic wave propagation parallel finite-difference code was used to obtain low-frequency (≤ 0:3 Hz) 3D synthetic seismograms for the Wenchuan earthquake. The code was run on KanBalam (Universidad Nacional Autónoma de México, Mexico) and HECToR (UK National Supercomputing Service) supercomputers. The modeling included the U.S. Geological Survey 40 × 315 km 2 kinematic description of the earthquake's rupture, embedded in a 2400 × 1600 × 300 km 3 physical domain, spatially discretized at 1 km in the three directions and a temporal discretization of 0.03 s. The compression and shear wave velocities and densities of the geologic structure used were obtained from recently published geophysical studies performed in the Sichuan region. The synthetic seismograms favorably compare with the observed ones for several station sites of the Seismological and Accelerographic Networks of China, such as MZQ, GYA, and TIY, located at about 90, 500, and 1200 km, respectively, from the epicenter of the Wenchuan event. Moreover, the comparisons of synthetic displacements with differential radar interferometry (DinSAR) ground deformation imagery, as well as of maximum velocity synthetic patterns with Mercalli modified intensity isoseist of the 2008 Wenchuan earthquake, are acceptable. 3D visualizations of the propagation of the event were also obtained; they show the source rupture directivity effects of the M w 7.9 Wenchuan event. Our results partially explain the extensive damage observed on the infrastructure and towns located in the neighborhood of the Wenchuan earthquake rupture zone.

show abstract

3D Parallel Elastodynamic Modeling of Large Subduction Earthquakes

Cited by 2 publications

References 4 publications

Benchmark Study of a 3d Parallel Code for the Propagation of Large Subduction Earthquakes

Benchmark Study of a 3d Parallel Code for the Propagation of Large Subduction Earthquakes

Low-Frequency 3D Wave Propagation Modeling of the 12 May 2008 Mw 7.9 Wenchuan Earthquake

Contact Info

Product

Resources

About