Seismo-Live is an interactive online learning platform for seismology. It utilizes Jupyter Notebooks. Most links on the front page lead to tools and projects that are used on Seismo-Live. The glowing green button will launch a light-weight remote virtual machine with ready to use interactive notebooks. Click on it to get started.

Once you clicked the green button, an overview of all available notebooks will be shown in a new tab. Choose any. If you have never seen this before, the Python Courses ones also contains a quick introduction on how to use the notebooks. There is also some official documentation on how to use the notebooks and lots of tutorials can be found in the internet.

go to seismo-live


Measurement of Rayleigh wave ellipticity from polarisation of ambient noise and inversion for the S-wave velocity profile.

get DOP-E





Originally, MSNoise was a “Python Package for Monitoring Seismic Velocity Changes using Ambient Seismic Noise”. With the release of MSNoise 1.4, and because of the Plugin Support, we could call MSNoise: “Measuring with Seismic Noise”.

The standard MSNoise workflow is designed to go from seismic data archives to dv/v curves. The monitoring is achieved by computing the cross-correlation of continuous seismic records for each pair of a network and by studying the changes in the crosscorrelation function relative to a reference.

The goal of the “suite” is to provide researchers with an efficient processing tool, while keeping the need for coding to a minimum and avoiding being a black box. Moreover, as long as the in- and outputs of each step are respected, they can easily be replaced with one’s own codes ! (See Workflow).

Plugins can be added and extend the standard workflow from any steps, e.g. using MSNoise as a cross-correlation toolbox until the stack step, and then branching to the workflow provided by one’s plugin.

get MSNoise



ObsPy is an open-source project dedicated to provide a Python framework for processing seismological data. It provides parsers for common file formats, clients to access data centers and seismological signal processing routines which allow the manipulation of seismological time series.

The goal of the ObsPy project is to facilitate rapid application development for seismology.

get ObsPy






Instaseis: Instant Global Seismograms Based on a Broadband Waveform Database

Instaseis calculates broadband seismograms from Green’s function databases generated with AxiSEM and allows for near instantaneous (on the order of milliseconds) extraction of seismograms. Using the 2.5D axisymmetric spectral element method, the generation of these databases, based on reciprocity of the Green’s functions, is very efficient and is approximately half as expensive as a single AxiSEM forward run. Thus this enables the computation of full databases at half the cost of the computation of seismograms for a single source in the previous scheme and hence allows to compute databases at the highest frequencies globally observed. By storing the basis coefficients of the numerical scheme (Lagrange polynomials), the Green’s functions are 4th order accurate in space and the spatial discretization respects discontinuities in the velocity model exactly. On top, AxiSEM allows to include 2D structure in the source receiver plane and readily includes other planets such as Mars.





AxiSEM is a parallel spectral-element method for 3D (an-)elastic, anisotropic and acoustic wave propagation in spherical domains. It requires axisymmetric background models and runs within a 2D computational domain, thereby reaching all desired highest observable frequencies (up to 2Hz) in global seismology. The Fortran2003 (OpenMP, MPI) code scales very well on supercomputers and is available here under the GPL license.

The code is also hosted and maintained at the Computational Infrastructure for Geodynamics (CIG), and code development takes place at CIG's github repository.

get AxiSEM


A Matlab package for phase-velocity dispersion measurement from ambient-noise correlations

GSpecDisp is an interactive package to measure phase-velocity dispersion of surface waves extracted by noise correlations. The input data are time-domain cross-correlations in SAC format.

GspecDisp can be run on any operation system with Matlab and also its Mapping toolbox. It is written based on the new Matlab Graphics which are introduced in Release 2014b.

get GSpecDisp




SPECFEM3D Cartesian simulates acoustic (fluid), elastic (solid), coupled acoustic/elastic, poroelastic or seismic wave propagation in any type of conforming mesh of hexahedra (structured or not.) It can, for instance, model seismic waves propagating in sedimentary basins or any other regional geological model following earthquakes. It can also be used for non-destructive testing or for ocean acoustics.

SPECFEM3D_Cartesian version 2.0 uses the continuous Galerkin spectral-element method, which can be seen as a particular case of the discontinuous Galerkin technique with optimized efficiency owing to its tensorized basis functions, to simulate forward and adjoint coupled acoustic-(an)elastic seismic wave propagation on arbitrary unstructured hexahedral meshes.






CUBIT is a full-featured software toolkit for robust generation of two- and three-dimensional finite element meshes (grids) and geometry preparation. Its main goal is to reduce the time to generate meshes, particularly large hex meshes of complicated, interlocking assemblies. It is a solid-modeler-based preprocessor that meshes volumes and surfaces for finite element analysis. Mesh generation algorithms include:

  • Quadrilateral and triangular paving
  • 2D and 3D mapping
  • Hex sweeping and multi-sweeping
  • Tet meshing
  • Many special purpose primitives.

CUBIT also contains many algorithms for controlling and automating much of the meshing process, such as

  • Automatic scheme selection
  • Interval matching
  • Sweep grouping
  • Sweep verification

And, of course, CUBIT also includes state-of-the-art smoothing algorithms.




Mineos computes synthetic seismograms in a spherically symmetric non-rotating Earth by summing normal modes.

Attenuation, gravity and transversal anisotropy effects may be optionally taken into account. The package computes mode eigenfrequencies and eigenfunctions, Green's functions and synthetic seismograms for sources with a moment tensor defined in time as a step-function.





The IRIS Synthetics Engine (Syngine) is a webservice that quickly returns synthetic seismograms custom requested by the user. The IRIS DMC stores precalculated TB scale databases of Green’s Functions for several different 1D reference models. The Green’s Functions were generated by Tarje Nissen-Meyer & colleagues using AxiSEM (Nissen-Meyer et al., 2014), an axi-symmetric spectral element method for 3D (an-)elastic, anisotropic and acoustic wave propagation in spherical domains. It requires axisymmetric background models and runs within a 2D computational domain, thereby reaching all desired highest observable frequencies (up to 2Hz) in global seismology. The requested synthetic seismograms are calculated using Instaseis (van Driel et al., 2015), a python library built for storing wavefield databases with a fast extraction algorithm. By using Instaseis and obspy, Syngine provides users with fully customizable synthetic seismograms for any source-receiver combination with a variety of source options including point source or user supplied moment tensors (or GCMT).

get syngine



The Adaptable Seismic Data Format (ASDF) is a modern file format intended for researchers and analysts. It combines the capability to create comprehensive data sets including all necessary meta information with high-performance parallel I/O for the most demanding use cases. Implementations for C/Fortran as well as Python are available.

get ASDF