Examples

The Pyrocko framework covers many aspects of seismological data handling and management. This section offers tutorials and real-world examples on how to utilize the framework.

Jupyter Notebooks

Example notebooks are available in the Pyrocko Notebooks Repository - Contributions to the collection are highly welcome!

Scripts collection

Annotated scripts and real-world use cases using Pyrocko, grouped roughly by how deep into the framework they go.

Getting started

Data access and handling basics - start here if you are new to Pyrocko.

• Squirrel-powered data access
  • Tutorial: writing a Squirrel based tool to calculate hourly RMS values
• Seismic traces
  • Load, filter and save
  • Visual inspection of traces
  • Creating a trace object from scratch
  • Extracting part of a trace (trimming)
  • Time shifting a trace
  • Resampling a trace
  • Convert SAC to MiniSEED
  • Convert MiniSEED to ASCII
  • Finding the comparative misfits of multiple traces
  • Restitute to displacement using poles and zeros
  • Restitute to displacement using SEED RESP response
• Metadata read & write
  • StationXML import
  • Pyrocko stations to StationXML
  • Create a StationXML file with flat displacement responses
• Downloading seismic data (FDSN)
  • Downloading and restituting waveform data
  • StationXML data manipulation
• Earthquake catalog
  • QuakeML import
  • Creating QuakeML from scratch
  • Searching an online earthquake catalog
• Snuffler markers
  • Read markers from file exported by Snuffler

Core toolbox

Datasets, coordinates, plotting, and other tools used across most Pyrocko-based projects.

• Traveltime calculation and raytracing
  • Calculate traveltimes in layered media
  • Calculate traveltimes in heterogeneous media
  • Traveltime table interpolation
• Accessing crustal velocity databases
  • Crust 2.0 Database
  • Global Crustal Database
• Moment tensor conversions
  • Convert moment tensor components to strike, dip and rake
  • Strike, dip and rake to moment tensor
• GPS (GNSS) data handling
  • Loading a GNSS campaign from CSV files
  • Loading and mapping of GNSS campaign data from UNR
• Geographical datasets
  • Topography
  • GSHHG coastal database
  • Tectonic plates and boundaries (PB2003)
  • Global strain rate (GSRM1)
• Plotting functions
  • Generating topographic maps with automap
  • Plotting beachballs (focal mechanisms)
  • Hudson’s source type plot
  • Source radiation plot
  • Pseudo dynamic rupture - slip map, slip movie, source plots
• GmtPy - A Python interface to GMT 4/5
  • Introduction
  • Tutorial
  • Random remarks
• Dataset management - Squirrel
  • Reorganizing a dataset into hour-files
  • Downsampling a whole dataset
  • Converting a dataset from Mini-SEED to SAC format
• Geodesic functions
  • Distance between points on earth
  • Azimuth and backazimuth
  • Conversion to carthesian coordinates
  • Relative carthesian coordinates to latitude and longitude
• ObsPy compatibility
  • Using the Snuffler on obspy.Stream
• Working with the guts package
  • What is Guts?
  • Defining classes and composition
  • Creating objects
  • Validation
  • Regularization
  • Serialization
  • Serialization of attributes with defaults
  • Deserialization
  • Inheritance
  • Default value factories
  • Computed properties
  • Customized serialization
  • Numerical data with guts_array
  • Cloning objects
  • Reactivity
  • How Pyrocko makes use of Guts

Advanced functionality

Green’s function based forward modeling and source models for seismic and geodetic problems.

• Forward modeling synthetic seismograms and displacements
  • Calculate synthetic seismograms from a local GF store
  • Calculate synthetic seismograms using the Pseudo Dynamic Rupture
  • Calculate Coulomb Failure Stress (CFS) changes using the Pseudo Dynamic Rupture
  • Calculate spatial surface displacement from a local GF store
  • Calculate spatial surface displacement using subfault dislocations
  • Calculate spatial surface displacement using the Pseudo Dynamic Rupture
  • Calculate spatial surface displacement and export Kite scenes
  • Calculate forward model of thrust faulting and display wrapped phase
  • Combining dislocation sources
  • Modelling viscoelastic static displacement
  • Creating a custom Source Time Function (STF)
• Kinematic/Dynamic source parameter modeling/inversion
  • Calculate subfault dislocations from tractions with Okada half-space equation