Instaseis Tutorial

Part 3: Plot a Record Section

Seismo-Live: http://seismo-live.org

Authors:

• Martin van Driel (@martinvandriel)
• Lion Krischer (@krischer)

---

Advanced Exercise 1: Plot Record Section

Use Instaseis to calculate a record section of your choice.

Notes

1. ObsPy Trace objects provide a function times() to conveniently get the time axis for plotting and the normalize() function to normalize seismograms to a common amplitude.
2. Use high pass filtering and deep sources if you want to enhance body waves
3. np.linspace() can help to generate equidistant stations
4. obspy.taup.TauPyModel can be used to compare with ray theoretical arrivals

---

Basic lines to set up the notebook and some paths.

%matplotlib inline
import matplotlib.pyplot as plt
import numpy as np
import os
import obspy
plt.style.use('ggplot')
plt.rcParams['figure.figsize'] = (10, 8)

Import Instaseis and open the database:

import instaseis
db = instaseis.open_db("data/database")

---

Solution