We use cookies to ensure that we give you the best experience on our website. By continuing to browse this repository, you give consent for essential cookies to be used. You can read more about our Privacy and Cookie Policy.

Durham Research Online
You are in:

pyprop8: A lightweight code to simulate seismic observables in a layered half-space

Valentine, Andrew P. and Sambridge, Malcolm (2022) 'pyprop8: A lightweight code to simulate seismic observables in a layered half-space.', Journal of Open Source Software, 7 (76). p. 4217.


The package pyprop8 enables calculation of the response of a 1-D layered halfspace to a seismic source, and also derivatives (‘sensitivity kernels’) of the wavefield with respect to source parameters. Seismograms, seismic spectra, and measures of static displacement (e.g. GPS, InSAR and field observations) may all be simulated. The method is based on a ThompsonHaskell propagator matrix algorithm, described in O’Toole & Woodhouse (2011) and O’Toole et al. (2012). The package is entirely written in Python, dependent only on the mainstream libraries numpy (Harris et al., 2020) and scipy (Virtanen et al., 2020). As such, it is lightweight and easy to deploy across a variety of platforms, making it particularly suited to use for teaching and outreach purposes.

Item Type:Article
Additional Information:The paper describes an open-source software package which is available here: or archived as a snapshot with DOI here: Video/animation produced using the package can be found here:
Full text:(VoR) Version of Record
Available under License - Creative Commons Attribution 4.0.
Download PDF
Publisher Web site:
Publisher statement:Authors of papers retain copyright and release the work under a Creative Commons Attribution 4.0 International License (CC BY 4.0).
Date accepted:No date available
Date deposited:31 August 2022
Date of first online publication:26 August 2022
Date first made open access:31 August 2022

Save or Share this output

Look up in GoogleScholar