Computational Infrastructure for GeoDynamics
The Computational Infrastructure for (CIG) develops, supports, and disseminates community-accessible software for the geodynamics research community. CIG software supports a variety of geodynamic research from mantle and core dynamics, to crustal and earthquake dynamics, to magma migration and seismology. CIG is a community-governed organization that is committed to develop and maintain the geodynamics community through community participation across this research spectrum.

SPECFEM3D_GLOBE simulates three-dimensional global and regional seismic wave propagation based upon the spectral-element method (SEM)
Princeton University's Global ShakeMovie uses SPECFEM3D_GLOBE to make public near real time global seismic waveform propagation visualizations.

Global ShakeMovie, Virginia Earthquake, August 23, 2011

CitcomS is a finite element code written in C that solves for thermo-chemical convection within a spherical shell. It can solve for problems within either a full or a restricted spherical domain. Although the code is capable of solving many different kinds of convection problems using the flexibility of finite elements, there are aspects of CitcomS which make it well-suited for solving problems in which the plate tectonic history is incorporated.

CitcomS Animations

Serpentine Wave Propagation - LLNL
WPP is a computer program for simulating seismic wave propagation on parallel machines. WPP implements substantial capabilities for 3-D seismic modeling, with a free surface condition on the top boundary, non-reflecting far-field boundary conditions on the other boundaries, point force and point moment tensor source terms with many predefined time dependencies, fully 3-D heterogeneous material model specification, elastic or visco-elastic materials, output of synthetic seismograms in the SAC format, output of GMT scripts for laying out simulation information on a map, and output of 2-D slices of (derived quantities of) the solution field as well as the material model.

WPP and its supporting libraries are built on Intel based desktops and laptops running LINUX and OSX. It has also been built on various supercomputers such as the large Linux clusters at LLNL (currently zeus, atlas and sierra), as well as on IBM BG/L and BG/P systems. WPP is likely to build if the underlying third party libraries can be built. Gnu: g++/gcc/gfortran versions 4.3 to 4.5 Intel: icpc/icc/ifort versions 9.1 to 11.1

User’s guide to WPP version 2.1.5

High Performance Computing at CNSI UC Santa Barbara

E3D, 3-D Elastic Seismic Wave Propagation Code - Nuclear Energy Agency
E3D is capable of simulating seismic wave propagation in a 3D heterogeneous earth. Seismic waves are initiated by earthquake, explosive, and/or other sources. These waves propagate through a 3D geologic model, and are simulated as synthetic seismograms or other graphical output.


JRG - Seismic Processing on Any Computer, Open-Source Software - John N. Louie
The Resource Geology Seismic Processing System for Java (JRG) is a basic reflection processing package with superb graphics, 3-d and crooked-line capabilities, SEG-Y and sound file I/O, and a friendly GUI that runs on any machine.

Converting E3D Images to Movies - John N. Louie
This tutorial will help you convert sequences of images output by Shawn Larsen’s E3D seismic-modeling code to movies that you can play on a computer, like in a Powerpoint presentation.

ModelAssembler Community Seismic Modeling Environment - John N. Louie
ModelAssembler is an open-source, Java-based velocity-model gridding code that can integrate scattered and heterogeneous geophysical data sets. The ModelAssembler has a tutorial graphical interface that sets up the input file for earthquake-scenario modeling runs with LLNL's E3D or WPP codes.

John N. Louie Youtube Earthquake Scenarios

Pahrump M7.2 Training Scenario
Animation of 0.1-Hz wave propagation from a 45-km-long Pahrump Valley fault rupture and shaking in Las Vegas, Nevada. Timelapse compression of 100 seconds of wave propagation into a 10-sec movie at 30 frames/sec.

Visualization Toolkit (VTK)
The Visualization Toolkit (VTK) is an open-source, freely available software system for 3D computer graphics, image processing and visualization. VTK consists of a C++ class library and several interpreted interface layers including Tcl/Tk, Java, and Python. VTK supports a wide variety of visualization algorithms including: scalar, vector, tensor, texture, and volumetric methods; and advanced modeling techniques such as: implicit modeling, polygon reduction, mesh smoothing, cutting, contouring, and Delaunay triangulation.

ParaView - Open Source Scientific Visualization
ParaView is an open source, freely available program for parallel, interactive, scientific visualization. It has a client–server architecture to facilitate remote visualization of datasets, and generates level of detail (LOD) models to maintain interactive framerates for large datasets. It is an application built on top of the Visualization Tool Kit (VTK) libraries. Where VTK is a set of libraries that provide visualization services for data, task, and pipeline parallelism, ParaView is an application designed for data parallelism on shared-memory or distributed-memory multicomputers and clusters. It can also be run as a single-computer application.

Global Seismic Wave Propagation Simulation
Contributors from ICES, The University of Texas at Austin:
Carsten Burstedde, Omar Ghattas, James R. Martin, Georg Stadler, Lucas C. Wilcox
Visualization by Gregory D. Abram, Texas Advanced Computing Center, The University of Texas at Austin.
Abstract: Modeling the propagation of seismic waves through the earth is an essential first step to inferring the structure of earth's interior. The propagation of waves through the oceans and outer core is modeled with the acoustic wave equation, while the elastic wave equation is employed for the crust, mantle, and inner core. We create a computational mesh in parallel by partitioning the earth into warped hexagonal elements whose size is adapted to the local seismic wavelengths. We employ a highly parallel and scalable discontinuous Galerkin method to discretize the governing elastodynamics equations on the adapted mesh. The mesh is partitioned for parallel execution using a space-filling curve. The visualization corresponds to a simulation with a central source frequency of 1/85 Hz and contains about 100 million degrees of freedom.