Earthquakes: The Quake Project

Jacobo Bielak, Omar Ghattas,
Volkan Akcelik, Hesheng Bao, Ioannis Epanomeritakis, Loukas Kallivokas,
Eui Joong Kim, David O'Hallaron, Jonathan Shewchuk, Tiankai Tu, Jifeng Xu,
Carnegie Mellon University
CMU scientists used the PSC's T3E and now up to 3000 processors on our terascale system, Lemieux, for problem sizes up to 100 million elements, to investigate how soil composition affects ground motion during earthquakes.

Simulation of 1994 Northridge Earthquake Aftershock

These first images and animations show their simulation of a 1994 Northridge earthquake aftershock. This area is an alluvial basin- soil and soft rock contained by a bowl shaped space within denser rock.
The following animation shows a surface view of displacement amplitude and is composited over landsat data (provided by Bill Harbert, U of Pittsburgh) to show the local terrain.
Still frame from earthquake simulation
		animation.
animation by G.Foss, PSC
play the animation (mpeg - 1.24meg)
This volume rendering shows how the shock wave radiates upwards until it reaches less dense rock, then becomes amplified as it enters the basin.
Still frame from earthquake simulation animation.
animation by G.Foss, PSC
play the animation(mpeg - 1.1meg)
These next two animations show visualizations of simulated ground motion of an earthquake generated along a vertical (strike-slip) fault. Oblique view from above and below ground, showing horizontal velocity amplitude pattern in the vicinity of the fault, at increasing times after onset of earthquake. Earthquake source originates near the far edge of the computational domain and travels along the fault toward the near edge. Seismic waves travel both perpendicularly to the fault and downstream from the epicenter. At first, ground motion is strongest in the direction parallel to the fault (initial red crescent moving away from fault), and later in the perependicular direction (initial green cresecent becoming red as it moves downstream). Peak fault-normal velocity is twice as large as peak fault-parallel velocity.
Still frame from earthquake simulation animation.
animation by G.Foss, PSC
play the animation (mpeg - 5.5meg)
Still frame from earthquake simulation animation.
animation by G.Foss, PSC
play the animation (mpeg - 4.3meg)
The two views below show the results of inversion of synthetic surface seismograms for a portion of the San Fernando Valley, California. The inverse problem is solved on a 1293 grid using an acoustic wave propagation forward model. Each animation compares the Inverted (left) with Target (right) compressional wave velocity. (wave velocity km/s, domain 32x32x16 [depth] km)
Still frame from earthquake simulation animation.
animation by G.Foss, PSC
play the animation (mpeg - 0.74meg)
Still frame from earthquake simulation animation.
animation by G.Foss, PSC
play the animation (mpeg - 2.6meg)
This image shows inversion of synthetic surface seismograms for soil stiffness of San Fernando basin model. (read left->right, top->bottom) First three images show progression of inverted models. Final image shows target basin.
image from earthquake simulation.
image by G.Foss, PSC
select image for larger size
Next is animation of surface response from a simulation of the 1994 Northridge earthquake in an 80x80x30 km3 Los Angeles Basin model. Color indicates magnitude of horizontal component of velocity.
Still frame from earthquake simulation animation.
image by Eui Joong Kim, CMU / animation by G.Foss, PSC
play the animation (mpeg - 4meg)

Lastly a simulation of the 1995 Kobe, Japan main shock

Still frame from earthquake simulation animation.
animation by G.Foss, PSC
play the animation (mpeg - 7.4meg)
More information: