seismic tomography: art or science? frederik j simons princeton university
TRANSCRIPT
![Page 1: Seismic tomography: Art or science? Frederik J Simons Princeton University](https://reader037.vdocument.in/reader037/viewer/2022110322/56649d215503460f949f6fa4/html5/thumbnails/1.jpg)
Seismic tomography: Art or science?
Frederik J SimonsPrinceton University
![Page 2: Seismic tomography: Art or science? Frederik J Simons Princeton University](https://reader037.vdocument.in/reader037/viewer/2022110322/56649d215503460f949f6fa4/html5/thumbnails/2.jpg)
What’s inside the Earth?
Dalton, Nature 2003
Only seismic waves have actually been there, done that
www.glencoe.com
www.usgs.gov
![Page 3: Seismic tomography: Art or science? Frederik J Simons Princeton University](https://reader037.vdocument.in/reader037/viewer/2022110322/56649d215503460f949f6fa4/html5/thumbnails/3.jpg)
This looks more complicated than it is;and that’s my point.
![Page 4: Seismic tomography: Art or science? Frederik J Simons Princeton University](https://reader037.vdocument.in/reader037/viewer/2022110322/56649d215503460f949f6fa4/html5/thumbnails/4.jpg)
This looks simpler than it is;and that’s my point.
![Page 5: Seismic tomography: Art or science? Frederik J Simons Princeton University](https://reader037.vdocument.in/reader037/viewer/2022110322/56649d215503460f949f6fa4/html5/thumbnails/5.jpg)
X-Ray attenuation tomography
Projections from all angles:X-ray intensity
Reconstructed image:X-ray attenuation constants
![Page 6: Seismic tomography: Art or science? Frederik J Simons Princeton University](https://reader037.vdocument.in/reader037/viewer/2022110322/56649d215503460f949f6fa4/html5/thumbnails/6.jpg)
![Page 7: Seismic tomography: Art or science? Frederik J Simons Princeton University](https://reader037.vdocument.in/reader037/viewer/2022110322/56649d215503460f949f6fa4/html5/thumbnails/7.jpg)
Seismic wavespeed tomography
Projections from all angles:Waveforms and arrival times
Reconstructed image:Wavespeed variations
![Page 8: Seismic tomography: Art or science? Frederik J Simons Princeton University](https://reader037.vdocument.in/reader037/viewer/2022110322/56649d215503460f949f6fa4/html5/thumbnails/8.jpg)
Forward modeling of the wave field, Part I:
Ray tracing, most 1-D
Before After
Buland, BSSA, 1983Kennett, GJI, 1995 Bullen & Bolt, 1985
![Page 9: Seismic tomography: Art or science? Frederik J Simons Princeton University](https://reader037.vdocument.in/reader037/viewer/2022110322/56649d215503460f949f6fa4/html5/thumbnails/9.jpg)
Forward modeling of the wave field, Part II:
Normal-mode summation, 1-D
AfterBefore
Dahlen & Tromp, 1998 Simons, Lithos, 1999
![Page 10: Seismic tomography: Art or science? Frederik J Simons Princeton University](https://reader037.vdocument.in/reader037/viewer/2022110322/56649d215503460f949f6fa4/html5/thumbnails/10.jpg)
Forward modeling of the wave field, Part III:
Spectral-element methods, 3-DBefore
After
Komatitsch, GJI, 2002
![Page 11: Seismic tomography: Art or science? Frederik J Simons Princeton University](https://reader037.vdocument.in/reader037/viewer/2022110322/56649d215503460f949f6fa4/html5/thumbnails/11.jpg)
![Page 12: Seismic tomography: Art or science? Frederik J Simons Princeton University](https://reader037.vdocument.in/reader037/viewer/2022110322/56649d215503460f949f6fa4/html5/thumbnails/12.jpg)
Non-continuous source coverage
The CMT catalog of large events
![Page 13: Seismic tomography: Art or science? Frederik J Simons Princeton University](https://reader037.vdocument.in/reader037/viewer/2022110322/56649d215503460f949f6fa4/html5/thumbnails/13.jpg)
Source location – (in)extricably linked
After
Before
Schaff, JGR, 2002Source relocation is big business.
![Page 14: Seismic tomography: Art or science? Frederik J Simons Princeton University](https://reader037.vdocument.in/reader037/viewer/2022110322/56649d215503460f949f6fa4/html5/thumbnails/14.jpg)
![Page 15: Seismic tomography: Art or science? Frederik J Simons Princeton University](https://reader037.vdocument.in/reader037/viewer/2022110322/56649d215503460f949f6fa4/html5/thumbnails/15.jpg)
Fermat’s Principle at Work for you
Zhao, PEPI, 2004
![Page 16: Seismic tomography: Art or science? Frederik J Simons Princeton University](https://reader037.vdocument.in/reader037/viewer/2022110322/56649d215503460f949f6fa4/html5/thumbnails/16.jpg)
The reference Earth: Radial models
![Page 17: Seismic tomography: Art or science? Frederik J Simons Princeton University](https://reader037.vdocument.in/reader037/viewer/2022110322/56649d215503460f949f6fa4/html5/thumbnails/17.jpg)
… and at least some of it is true…
Karki et al., Rev. Geoph., 2001 Jackson, 1998
![Page 18: Seismic tomography: Art or science? Frederik J Simons Princeton University](https://reader037.vdocument.in/reader037/viewer/2022110322/56649d215503460f949f6fa4/html5/thumbnails/18.jpg)
![Page 19: Seismic tomography: Art or science? Frederik J Simons Princeton University](https://reader037.vdocument.in/reader037/viewer/2022110322/56649d215503460f949f6fa4/html5/thumbnails/19.jpg)
Menke, 1989
![Page 20: Seismic tomography: Art or science? Frederik J Simons Princeton University](https://reader037.vdocument.in/reader037/viewer/2022110322/56649d215503460f949f6fa4/html5/thumbnails/20.jpg)
Picking the right continent
A dense path coverage minimizes the amount of a priori information needed
Receiver coverage
Simons, GJI 2002
![Page 21: Seismic tomography: Art or science? Frederik J Simons Princeton University](https://reader037.vdocument.in/reader037/viewer/2022110322/56649d215503460f949f6fa4/html5/thumbnails/21.jpg)
![Page 22: Seismic tomography: Art or science? Frederik J Simons Princeton University](https://reader037.vdocument.in/reader037/viewer/2022110322/56649d215503460f949f6fa4/html5/thumbnails/22.jpg)
Regularization: the Mathematics
![Page 23: Seismic tomography: Art or science? Frederik J Simons Princeton University](https://reader037.vdocument.in/reader037/viewer/2022110322/56649d215503460f949f6fa4/html5/thumbnails/23.jpg)
Regularization: the Physics
Such “fat” rays sample more of the Earth and thus we need fewer of them to have a well-constrained tomographic problem.
Dahlen, GJI, 2002
![Page 24: Seismic tomography: Art or science? Frederik J Simons Princeton University](https://reader037.vdocument.in/reader037/viewer/2022110322/56649d215503460f949f6fa4/html5/thumbnails/24.jpg)
Regularization: the Art
Too much?Too smooth?
Too little?Too rough?
Simons, Lithos, 1999
![Page 25: Seismic tomography: Art or science? Frederik J Simons Princeton University](https://reader037.vdocument.in/reader037/viewer/2022110322/56649d215503460f949f6fa4/html5/thumbnails/25.jpg)
How to interpret seismic models
Pillet, PEPI, 1999
![Page 26: Seismic tomography: Art or science? Frederik J Simons Princeton University](https://reader037.vdocument.in/reader037/viewer/2022110322/56649d215503460f949f6fa4/html5/thumbnails/26.jpg)
Demand to see the ray paths
Pillet, PEPI, 1999
![Page 27: Seismic tomography: Art or science? Frederik J Simons Princeton University](https://reader037.vdocument.in/reader037/viewer/2022110322/56649d215503460f949f6fa4/html5/thumbnails/27.jpg)
Nature isn’t always kind
Wolfe, Nature, 1997Shen, Nature, 1998
![Page 28: Seismic tomography: Art or science? Frederik J Simons Princeton University](https://reader037.vdocument.in/reader037/viewer/2022110322/56649d215503460f949f6fa4/html5/thumbnails/28.jpg)
Seismic anisotropy Wave speeds depend on
propagation direction and polarization:
No surprise: elasticity maps stress and strain, and both depend on three directions
![Page 29: Seismic tomography: Art or science? Frederik J Simons Princeton University](https://reader037.vdocument.in/reader037/viewer/2022110322/56649d215503460f949f6fa4/html5/thumbnails/29.jpg)
Polarization anisotropy• The particles of Love and Raleigh surface
waves move in orthogonal directions
• SH and SV body waves sometimes exhibit clear splitting
Azimuthal anisotropy• It’s usually very hard to separate whether
the time difference arises from an anisotropic direction or an isotropic wave speed difference (aka heterogeneity)
![Page 30: Seismic tomography: Art or science? Frederik J Simons Princeton University](https://reader037.vdocument.in/reader037/viewer/2022110322/56649d215503460f949f6fa4/html5/thumbnails/30.jpg)
Why is this so hard?
It’s very hard to tell whether a phase comes in early because it went through a fast patch or because it came in a fast direction – heterogeneity and anisotropy “trade off.”
![Page 31: Seismic tomography: Art or science? Frederik J Simons Princeton University](https://reader037.vdocument.in/reader037/viewer/2022110322/56649d215503460f949f6fa4/html5/thumbnails/31.jpg)
Questions to ask of the tomographer
• How is the forward model computed?
• What is the ray coverage?
• What (sort of) damping did you use?
• What does velocity estimation trade off with?
• What is the grid size / the correlation length?
• How are different data sets weighted?
• How far is the final from the starting model?
• Does the starting model have discontinuities?
• How is the surface/depth parameterization
• Is your sensitivity 1-D, 2-D,or 3-D?
![Page 32: Seismic tomography: Art or science? Frederik J Simons Princeton University](https://reader037.vdocument.in/reader037/viewer/2022110322/56649d215503460f949f6fa4/html5/thumbnails/32.jpg)
Journey to Middle Earth, Part I:
The continental lithosphere
Gung, Nature, 2003
Simons, GRL, 2002
![Page 33: Seismic tomography: Art or science? Frederik J Simons Princeton University](https://reader037.vdocument.in/reader037/viewer/2022110322/56649d215503460f949f6fa4/html5/thumbnails/33.jpg)
Journey to Middle Earth, Part II:
Subduction zones
Replumaz, EPSL, 2004
![Page 34: Seismic tomography: Art or science? Frederik J Simons Princeton University](https://reader037.vdocument.in/reader037/viewer/2022110322/56649d215503460f949f6fa4/html5/thumbnails/34.jpg)
Journey to Middle Earth, Part III:
Deep mantle plumes
Montelli, Science, 2004
![Page 35: Seismic tomography: Art or science? Frederik J Simons Princeton University](https://reader037.vdocument.in/reader037/viewer/2022110322/56649d215503460f949f6fa4/html5/thumbnails/35.jpg)
What does it all mean? Part I:
Temperature anomalies
Goes, JGR, 2002110 km
![Page 36: Seismic tomography: Art or science? Frederik J Simons Princeton University](https://reader037.vdocument.in/reader037/viewer/2022110322/56649d215503460f949f6fa4/html5/thumbnails/36.jpg)
What does it all mean? Part II:
Compositional anomalies
Perry, GJI, 2003150 km
Fe/(Fe + Mg) T
![Page 37: Seismic tomography: Art or science? Frederik J Simons Princeton University](https://reader037.vdocument.in/reader037/viewer/2022110322/56649d215503460f949f6fa4/html5/thumbnails/37.jpg)
What does it all mean? Part III:
Deformation in the mantle
Simons, EPSL, 2003
Fossil Contemporaneous
![Page 38: Seismic tomography: Art or science? Frederik J Simons Princeton University](https://reader037.vdocument.in/reader037/viewer/2022110322/56649d215503460f949f6fa4/html5/thumbnails/38.jpg)
Conclusions• Ultimately, seismology can only tell us where,
or in which direction, wave propagation is faster or slower than a reference model
• The non-seismologist has to know the basics of inverse problem modeling, understand the sometimes poor constraints, and be critical
• Improvements are being made: better data, better forward models, better inversions
• As much as with the a posteriori interpretation, the community needs to help defining a priori acceptable starting models
![Page 39: Seismic tomography: Art or science? Frederik J Simons Princeton University](https://reader037.vdocument.in/reader037/viewer/2022110322/56649d215503460f949f6fa4/html5/thumbnails/39.jpg)
![Page 40: Seismic tomography: Art or science? Frederik J Simons Princeton University](https://reader037.vdocument.in/reader037/viewer/2022110322/56649d215503460f949f6fa4/html5/thumbnails/40.jpg)
![Page 41: Seismic tomography: Art or science? Frederik J Simons Princeton University](https://reader037.vdocument.in/reader037/viewer/2022110322/56649d215503460f949f6fa4/html5/thumbnails/41.jpg)
![Page 42: Seismic tomography: Art or science? Frederik J Simons Princeton University](https://reader037.vdocument.in/reader037/viewer/2022110322/56649d215503460f949f6fa4/html5/thumbnails/42.jpg)
More equations, for completeness
![Page 43: Seismic tomography: Art or science? Frederik J Simons Princeton University](https://reader037.vdocument.in/reader037/viewer/2022110322/56649d215503460f949f6fa4/html5/thumbnails/43.jpg)
![Page 44: Seismic tomography: Art or science? Frederik J Simons Princeton University](https://reader037.vdocument.in/reader037/viewer/2022110322/56649d215503460f949f6fa4/html5/thumbnails/44.jpg)
![Page 45: Seismic tomography: Art or science? Frederik J Simons Princeton University](https://reader037.vdocument.in/reader037/viewer/2022110322/56649d215503460f949f6fa4/html5/thumbnails/45.jpg)
![Page 46: Seismic tomography: Art or science? Frederik J Simons Princeton University](https://reader037.vdocument.in/reader037/viewer/2022110322/56649d215503460f949f6fa4/html5/thumbnails/46.jpg)
![Page 47: Seismic tomography: Art or science? Frederik J Simons Princeton University](https://reader037.vdocument.in/reader037/viewer/2022110322/56649d215503460f949f6fa4/html5/thumbnails/47.jpg)