SHaxi animations

Movie Credits:

This page contains animations of seismic wave propagation. The global animations in this section are animations of SH-wave propagation. These animations were created with the axi-symmetric finite difference code SHaxi. The code is presented in: Jahnke, G., Thorne, M.S., Cochard, A., and Igel, H. Global SH-wave propagation using a parallel axi-symmetric finite difference scheme, in review GJI, 2008.

All computations were performed on the GCC linux cluster at ASU, the Hitachi SR800 at the Leibniz Rechenzentrum, or ICEBERG at the Arctic Region Supercomputing Center.

The SHaxi source code can be downloaded at the SPICE webpage here.

Movie formats:

All animations on this page are available as Quicktime Movies. To play the movies, the Quicktime Player is required. For best playback saving the movies to your harddrive is recommended. To download the movies right click on the image | Save Target As...

In the Quicktime player, selecting 'Movie | Loop' will continuously loop the movies.

Movies are now available in .m4v format for use with Itunes/Ipod.

01 Homogeneous Velocities; 0 km source depth
homogeneous_earth_0km.mov (25 Mb) Quicktime format
homogeneous_earth_0km.m4v (566 Kb) Ipod format
homogeneous_earth_0km.ogg Ogg format(best for Linux)
      • Vs: 2.0 km/sec
      • Density: 4.0 g/cm3
      • Source Depth: 0 km
      • Dominant Period: 30 sec
      • Plot Range: -60° < θ < 60°
      • Non-linear scale: 0.5
      • Z-grid points: 2,000
      • Theta-grid points: 10,000

02 Homogeneous Velocities; 1400 km source depth
homogeneous_earth_1400km.mov (28.5 Mb) Quicktime format
homogeneous_earth_1400km.m4v (1.1 Mb) Ipod format
homogeneous_earth_1400km.ogg Ogg format(best for Linux)
      • Vs: 2.0 km/sec
      • Density: 4.0 g/cm3
      • Source Depth: 1400 km
      • Dominant Period: 30 sec
      • Plot Range: -60° < θ < 60°
      • Non-linear scale: 0.5
      • Z-grid points: 2,000
      • Theta-grid points: 10,000

03 Homogeneous Velocities; Low velocity Layer
LowV_Layer.mov (34.6 Mb) Quicktime format
LowV_Layer.m4v (1.8 Mb) Ipod format
LowV_Layer.ogg Ogg format(best for Linux)
      • Homogeneous velocity model with a low velocity layer. Model has density of 4.0 g/cm3 and Vs is 3.5 km/sec except in low velocity layer. The low velocity layer is contained between radii 5200 and 5800 km and has a velocity of 2.5 km/sec.
      • Source Depth: 0 km
      • Dominant Period: 30 sec
      • Plot Range: -70° < θ < 70°
      • Non-linear scale: 0.5
      • Z-grid points: 2,000
      • Theta-grid points: 10,000

04 Homogeneous Velocities; High velocity Layer
HiV_Layer.mov (39.8 Mb) Quicktime format
HiV_Layer.m4v (1.8 Mb) Ipod format
HiV_Layer.ogg Ogg format(best for Linux)
      • Homogeneous velocity model with a high velocity layer. Model has density of 4.0 g/cm3 and Vs is 3.5 km/sec except in high velocity layer. The high velocity layer is contained between radii 5200 and 5800 km and has a velocity of 4.5 km/sec.
      • Source Depth: 0 km
      • Dominant Period: 30 sec
      • Plot Range: -85° < θ < 85°
      • Non-linear scale: 0.5
      • Z-grid points: 2,000
      • Theta-grid points: 10,000

05 Homogeneous Velocities; Low velocity body (50 km radius)
Low_Vbody_50km.mov (30.6 Mb) Quicktime format
Low_Vbody_50km.m4v (842 Kb) Ipod format
Low_Vbody_50km.ogg Ogg format(best for Linux)
      • Vs 3.5 km/sec except in the low velocity body which has Vs reduction of 5% (Velocity = 3.325 km/sec). Radius of low velociy region is 50.0 km, and is centered at a depth of 1371 km, and epicentral distance of 30°. 30 sec Dominant Period corresponds to a wavelength of 105 km in the homogeneous velocity model. Thus, the size of the low velocity body is on the order of 1-dominant wavelength.
      • Source Depth: 0 km
      • Dominant Period: 30 sec
      • Plot Range: -10° < θ < 75°
      • Non-linear scale: 0.5
      • Z-grid points: 2,000
      • Theta-grid points: 10,000

06 Homogeneous Velocities; Low velocity body (150 km radius)
Low_Vbody_150km.mov (29.5 Mb) Quicktime format
Low_Vbody_150km.m4v (937 Kb) Ipod format
Low_Vbody_150km.ogg Ogg format(best for Linux)
      • Vs 3.5 km/sec except in the low velocity body which has Vs reduction of 5% (Velocity = 3.325 km/sec). Radius of low velociy region is 50.0 km, and is centered at a depth of 1371 km, and epicentral distance of 30°. 30 sec Dominant Period corresponds to a wavelength of 105 km in the homogeneous velocity model. Thus, the size of the low velocity body is on the order of 3-dominant wavelengths.
      • Source Depth: 0 km
      • Dominant Period: 30 sec
      • Plot Range: -10° < θ < 75°
      • Non-linear scale: 0.5
      • Z-grid points: 2,000
      • Theta-grid points: 10,000

07 PREM model; 500 km source depth
prem_500km_12s.mov (19.2 Mb) Quicktime format
prem_500km_12s.m4v (1.6 Mb) Ipod format
prem_500km_12s.ogg Ogg format(best for Linux)
      • Model Type: PREM
      • Source Depth: 500 km
      • Dominant Period: 12 sec
      • Plot Range: -10° < θ < 140°
      • Non-linear scale: 0.3
      • Z-grid points: not recorded
      • Theta-grid points: not recorded

08 PREM model; 500 km source depth
shinjuku.mov (21.3 Mb) Quicktime format
shinjuku.m4v (1.1 Mb) Ipod format
shinjuku.ogg Ogg format(best for Linux)
      • Model Type: PREM
      • Source Depth: 500 km
      • Dominant Period: 15 sec
      • Plot Range: -180° < θ < 180°
      • Non-linear scale: not recorded
      • Z-grid points: not recorded
      • Theta-grid points: not recorded

09 PREM model with topography on the D" discontinuity
Ddoubleprime_interaction.mov (27.0 Mb) Quicktime format
Ddoubleprime_interaction.m4v (3.7 Mb) Ipod format
Ddoubleprime_interaction.ogg Ogg format(best for Linux)
      • Model Type: Model Type: PREM w/ D"; discontinuity based on Thomas et al. (2004)
      • Source Depth: 500 km
      • Dominant Period: 10 sec
      • Plot Range Inset: -5° < θ < 90°
      • Plot Range Zoom (theta): 25° < θ < 55°
      • Plot Range Zoom (radius): 3480 < r < 4500 km
      • Non-linear scale: 0.3
      • Z-grid points: 2,750
      • Theta-grid points: 20,000
      • Reference: Thorne, M.S., Lay, T., Garnero, E.J., Jahnke, G., Igel, H., Seismic imaging of the laterally varying D" region beneath the Cocos Plate, Geophysical Journal International (170), pp. 635-648, doi: 10.1111/j.1365-246X.2006.03279.x, 2007. (pdf - 9.8 Mb)

10 PREM model with lens of post-perovskite in lower mantle
pPv_lens_animation.mov (7.0 Mb) Quicktime format
pPv_lens_animation.m4v (2.2 Mb) Ipod format
pPv_lens_animation.ogg Ogg format(best for Linux)
      • Model Type: Model Type: PREM w/ lens of pPv.
      • Source Depth: 600 km
      • Dominant Period: 10 sec
      • Reference: Lay, T., Hernlund, J., Garnero, E.J., Thorne, M.S., A lens of post-perovskite and CMB heat flux in an iron-rich pile in D" beneath the Central Pacific, Science (314), 1272-1276, 2006. (pdf - 323 kb) (online supplemental material - 1.6 Mb)

11 Clip from Naked Science: Earth's Core
naked_science.mov (21.4 Mb) Quicktime format
naked_science.m4v (984 Kb) Ipod format
naked_science.ogg Ogg format(best for Linux)
      • Clip of Ed Garnero from National Geographic Episode - Naked Science: The Earth's Core. Background animation is SH-wave propagation in PREM model at dominant period of 15 sec. Animation file is shinjuku.mov

PSVaxi animations

PSVaxi animations coming soon.

Regional animations

Regional animations were computed with the E3D finite difference code.

01 Regional animation to demonstrate pure wavefield without scattering.
noscatter.mov (4 Mb) Quicktime format
noscatter.m4v (1.5 Mb) Ipod format
noscatter.ogg Ogg format(best for Linux)
      • Regional model with a 1 sec dominant period Ricker source-time function and explosion source. Amplitudes are color-coded red and green for vertical and horizontal component energy respectively.
      • Background model is homogeneous with Vs = 4.0 km/sec, Vp = 5.0 km/sec, and density = 3.0 g/cm3.
      • Source is centered at 40 km depth and 40 km away from the left grid boundary.
      • Total grid dimensions are 80x120 km.

02 Regional animation to demonstrate wavefield with scattering.
scattering.mov (6.6 Mb) Quicktime format
scattering.m4v (1.7 Mb) Ipod format
scattering.ogg Ogg format(best for Linux)
      • Regional model with a 1 sec dominant period Ricker source-time function and explosion source. Amplitudes are color-coded red and green for vertical and horizontal component energy respectively.
      • Background model is homogeneous with Vs = 4.0 km/sec, Vp = 5.0 km/sec, and density = 3.0 g/cm3.
      • A von Karman autocorrelation function with 1 km correlation length is used to construct the random media. 3% RMS Vp perturbations are applied to the model.
      • Source is centered at 40 km depth and 40 km away from the left grid boundary.
      • Total grid dimensions are 80x120 km.

03 Seismic Refraction - Two layer model
Two Layer.mov (16.5 Mb) Quicktime format
Two_Layer.m4v (1.5 Mb) Ipod format
noscatter.ogg Ogg format(best for Linux)
      • Regional model with a 5 Hz dominant period Ricker source-time function and explosion source. Shown is the divergence of the wavefield so only P-wave energy is highlighted.
      • Background model is a simple 2-layer model. Both wavefield and raypaths for direct, reflected, and critically refracted arrivals are shown.
      • Travel Time curves for the animation setup can be downloaded HERE.

Other movie downloads

Other movie files will be placed here in the future.

Contact

  • email: michael.thorne at utah.edu