visualizing seismic wave propagation peter van kekenuniversity of michigan jeroen tromp dimitri...

21
Visualizing seismic wave propagation Peter van Keken University of Michigan Jeroen Tromp Dimitri Komatitsch Caltech Shalini Venkataraman Nicholas Schwarz Luc Renambot Jason Leigh University of Illinois at Chicago American Geophysical Union, Fall 2003, ED31E-01

Post on 23-Jan-2016

216 views

Category:

Documents


0 download

TRANSCRIPT

Page 1: Visualizing seismic wave propagation Peter van KekenUniversity of Michigan Jeroen Tromp Dimitri Komatitsch Caltech Shalini Venkataraman Nicholas Schwarz

Visualizing seismic wave propagation

Peter van Keken University of MichiganJeroen TrompDimitri Komatitsch CaltechShalini VenkataramanNicholas SchwarzLuc RenambotJason Leigh University of Illinois at Chicago

American Geophysical Union, Fall 2003, ED31E-01

Page 2: Visualizing seismic wave propagation Peter van KekenUniversity of Michigan Jeroen Tromp Dimitri Komatitsch Caltech Shalini Venkataraman Nicholas Schwarz

Seismology in intro geophysicsUniversity of Michigan GS420

• Elastic deformation/wave equation• Refraction, reflection• 1D earth structure• Earthquakes and seismograms• Tomography• Attenuation, anisotropy

• www.geo.lsa.umich.edu/~keken/420.html

Page 3: Visualizing seismic wave propagation Peter van KekenUniversity of Michigan Jeroen Tromp Dimitri Komatitsch Caltech Shalini Venkataraman Nicholas Schwarz
Page 4: Visualizing seismic wave propagation Peter van KekenUniversity of Michigan Jeroen Tromp Dimitri Komatitsch Caltech Shalini Venkataraman Nicholas Schwarz

From Shearer (1999)

Page 5: Visualizing seismic wave propagation Peter van KekenUniversity of Michigan Jeroen Tromp Dimitri Komatitsch Caltech Shalini Venkataraman Nicholas Schwarz

PcP: average P velocity in the mantle

PcP

P w

ave

shad

ow z

one

shadow zone

PKIKP

PKIKP: averageP velocity of the Earth

From Shearer (1999)

Page 6: Visualizing seismic wave propagation Peter van KekenUniversity of Michigan Jeroen Tromp Dimitri Komatitsch Caltech Shalini Venkataraman Nicholas Schwarz

Ideal: Develop a good understanding of seismic wave propagation

Greatly aided by improvements in:- data acquisition

- computational power- computer visualization

Page 7: Visualizing seismic wave propagation Peter van KekenUniversity of Michigan Jeroen Tromp Dimitri Komatitsch Caltech Shalini Venkataraman Nicholas Schwarz

Courtesy ExxonMobil

Example: active source seismics

Page 8: Visualizing seismic wave propagation Peter van KekenUniversity of Michigan Jeroen Tromp Dimitri Komatitsch Caltech Shalini Venkataraman Nicholas Schwarz

Examples: Global seismology

• Alan Jones: Seismic Waves– Seismic phases, seismograms– www.geol.binghamton.edu/faculty/jones

• Saadia Baker and Michael Wysession– Synthetic mantle shear waves, seismograms– Normal mode summation– http://epsc.wustl.edu/~saadia

Page 9: Visualizing seismic wave propagation Peter van KekenUniversity of Michigan Jeroen Tromp Dimitri Komatitsch Caltech Shalini Venkataraman Nicholas Schwarz
Page 10: Visualizing seismic wave propagation Peter van KekenUniversity of Michigan Jeroen Tromp Dimitri Komatitsch Caltech Shalini Venkataraman Nicholas Schwarz

01347121723Jo

nes:

sei

smic

wav

es

Page 11: Visualizing seismic wave propagation Peter van KekenUniversity of Michigan Jeroen Tromp Dimitri Komatitsch Caltech Shalini Venkataraman Nicholas Schwarz

Saa

dia

and

Wys

essi

on

Page 12: Visualizing seismic wave propagation Peter van KekenUniversity of Michigan Jeroen Tromp Dimitri Komatitsch Caltech Shalini Venkataraman Nicholas Schwarz
Page 13: Visualizing seismic wave propagation Peter van KekenUniversity of Michigan Jeroen Tromp Dimitri Komatitsch Caltech Shalini Venkataraman Nicholas Schwarz

New approach: SPECFEM3D

• Spectral finite element code (Komatitsch and Tromp, GJI, 2002;

Komatitsch, Tromp and Ritsema, Science, 298, 2002

www.geoframework.org)

• ‘Complete’ model: allows introduction of

– Ellipticity, anisotropy

– 3D mantle and crustal structure

• Massively parallel computation – Linux PC clusters at Caltech, Michigan

• Challenging code for visualization– Irregular grid

– Massive amount of data: O(Terabyte)

Page 14: Visualizing seismic wave propagation Peter van KekenUniversity of Michigan Jeroen Tromp Dimitri Komatitsch Caltech Shalini Venkataraman Nicholas Schwarz
Page 15: Visualizing seismic wave propagation Peter van KekenUniversity of Michigan Jeroen Tromp Dimitri Komatitsch Caltech Shalini Venkataraman Nicholas Schwarz

50s lowpass filter 14/8/95 New Britain (Mw=6.7)

Page 16: Visualizing seismic wave propagation Peter van KekenUniversity of Michigan Jeroen Tromp Dimitri Komatitsch Caltech Shalini Venkataraman Nicholas Schwarz

Quiz: Name the phases labeled A through F (white/red indicate positive amplitude only)

P S

Page 17: Visualizing seismic wave propagation Peter van KekenUniversity of Michigan Jeroen Tromp Dimitri Komatitsch Caltech Shalini Venkataraman Nicholas Schwarz

Vanuatu, 26/11/99

Page 18: Visualizing seismic wave propagation Peter van KekenUniversity of Michigan Jeroen Tromp Dimitri Komatitsch Caltech Shalini Venkataraman Nicholas Schwarz

Bolivia, 9/6/94

Page 19: Visualizing seismic wave propagation Peter van KekenUniversity of Michigan Jeroen Tromp Dimitri Komatitsch Caltech Shalini Venkataraman Nicholas Schwarz

3D = GeowallGS420 Fall 2002

Page 20: Visualizing seismic wave propagation Peter van KekenUniversity of Michigan Jeroen Tromp Dimitri Komatitsch Caltech Shalini Venkataraman Nicholas Schwarz

EVL Volume Rendering

• Use color and transparancy

(Alternative: PC based volume renderer HVR www.lcse.umn.edu)

5123 8 bit 5123 16 bit

Page 21: Visualizing seismic wave propagation Peter van KekenUniversity of Michigan Jeroen Tromp Dimitri Komatitsch Caltech Shalini Venkataraman Nicholas Schwarz

Quiz: Name the phases labeled A through F

SK

PP

PKP

Love

PPS

ScPP

ScSP ScSS