back

Anselme Borgeaud

Data scientist with a PhD in computational seismology and a strong interest in financial markets. Lausanne, Switzerland

Probabilistic seismic imaging using a Neighbourhood algorithm

pydsm can compute several seismic models in parallel, which is used in pytomo to run probabilistic waveform inversion with a neighbourhood algorithm (Sambridge, 1999). A few examples of using pytomo to probe the Earth's seismic structure are shown below.

D" discontinuity

The video below shows an example of synthetic inversion with a neighbourhood algorithm for the shear-velocity increase in the D'' layer and the depth of the D'' discontinuity at the bottom of the Earth's mantle (left: voronoi cells colored by waveform misfit; middle: target, and best inverted model for the current iteration; right: waveforms for the target, and current best inverted model).

D" anisotropy

The same neighbourhood algorithm can be used to infer shear-wave anisotropy (VSH and VSV) in the D" layer. A synthetic inversion is shown in the video below (left: voronoi cells colored by waveform misfit; middle: target, and best inverted model for the current iteration; right: radial and transverse component waveforms for the target, and current best inverted model).