Lab to Planet Geodynamics Jupyter Notebooks

• The Github repository includes teaching modules for simple driven harmonic oscillators, stick-slip friction, rate-state friction, the Lorentz system, heat conduction, Stokes flow (using an interface to FEniCX), as well as surface transport problems (with an interface to Landlab).

Megathrust Modeling Framework (MTMOD) Software pages

interact: boundary element fault interactions

An MPI parallel, PETSC and H-matrix (HMMVP, HACApK, HTOOLS) enabled boundary element, elastic half space fault system interaction and earthquake cycle code. Can compute displacements and stress given slip, or solve for slip given stress conditions, as well as simulate earthquake cycles based on rate state friction (five different evolution laws). Uses quad and triangular dislocations. Some say... it's the fastest RSF solver in town. The latest versions can be downloaded here

• interact code on GitHub
• README for the package
• README for the RSF cycle code
• releases

References/use cases:

• Becker, T. W. and Schott, B.: On boundary-element models of elastic fault interaction (abstract). Eos Trans. AGU, 83(47), Fall Meet. Suppl., Abstract NG62A-0925, 2002. (PDF)
• Becker, T. W., Hardebeck, J. L., and Anderson, G.: Constraints on fault slip rates of the southern California plate boundary from GPS velocity and stress inversions. Geophys. J. Int., 160, 634-650, 2005. (PDF)

Solid Earth Research and Teaching Environment (SEATREE)

We developed a python-based, modular, graphical and script driven interface to solid Earth geodynamics and seismology modeling tools. The goal is to provide easy access to research tools for teaching, and to facilitate interdisciplinary research. SEATREE includes an interface to the hc mantle flow tool and seismic mantle tomography tools. See the SEATREE project web page for the current status (updated 2025), and for access to the source code, see the

• download and installation from the SEATREE GitHub site
• CIG CIDER2016 Summer Program VirtualBox (use hc_user, no password, run SEATREE from the desktop) - newer version on GitHub

Reference:

• Milner, K., Becker, T. W., Boschi, L., Sain, J., Schorlemmer, D. and H. Waterhouse: The Solid Earth Research and Teaching Environment: a new software framework to share research tools in the classroom and across disciplines. Eos Trans. AGU, 90, 12, 2009. (PDF).

fstrack: Seismic anisotropy from flow predictor

Package with tools to analyze symmetry components of elastic tensors, predict synthetic waveforms and compute automated shear wave splitting along ray paths, and to track finite strain and predict LPO from mantle flow given on GMT/netcdf grds.

Includes Kaminski & Ribe (2001) type LPO texture formation modeling and Browaeys and Chevrot (2003) decomposition routines (as provided by DREX by Kaminski), as well as single and multiple layer SKS seismogram synthetics and splitting computation (codes from Schulte-Pelkum and Blackman). A version of fstrack was used for Becker et al. (2006a, 2006b, 2008).

Please keep in mind that this is a research product and so subject to change; installation and usage is described in README files, but may require some user intervention. Also note that a lot of the software in this package is copyrighted by other people, and you may wish to use the original codes (e.g. DREX) rather than our versions.

• fstrack code on GitHub
• complete tar balls:
  • latest fstrack release Allows compilation of anisotropy and LPO routines indendent from flow tracker (no reliance on GMT and HC).
  • fstrack-0.6.112725.tgz, stable release
  • fstrack-0.3.062614.tgz: old version, used for prior work up and including 2014.

The main references for fstrack are:

• Becker, T. W., Schulte-Pelkum, V., Blackman, D. K., Kellogg, J. B., and O'Connell, R. J.: Mantle flow under the western United States from shear wave splitting, Earth Planet. Sci. Lett., 247, 235-251, 2006. (PDF)
• Becker, T. W., Kustowski, B., Ekström, G.: Radial seismic anisotropy as a constraint for upper mantle rheology. Earth Planet. Sci. Lett., 267, 213-237, 2008. (PDF)

hc: Hager & O'Connell spherical flow computations

Global mantle circulation solver following Hager & O'Connell (1981). This software is a modular, C-language rewrite of Bernhard Steinberger's FORTRAN code and developed by Becker, Craig O'Neill, and others. The code has been tested for computation of velocities, tractions, and the geoid given incompressible, Newtonian flow in the mantle for only radially varying viscosity.

The spherical harmonic format from Becker & Boschi (2002) can be read by HC to compute velocities for spherical harmonic models, and a graphical user interface to HC is provided by SEATREE.

Availability and documentation:

• HC is one of the modules of SEATREE which also provides a graphical user interface, editing tools, and example input data. I recommend to download HC with the SEATREE package, both also included fully compiled and installed in
  • HC code on GitHub
  • CIG CIDER2016 Summer Program VirtualBox (use hc_user, no password, run SEATREE from the desktop). contains old release
• Tar balls:
  • releases
  • hc-1.1.2.050824.tgz

Reference:

• Becker, T. W., O'Neill, C., and Steinberger, B. (2009): HC, a global mantle circulation solver. Available online at https://github.com/geodynamics/hc, accessed 03/2022.

shansyn: spherical harmonics packages for GMT/Netcdf grids

My spherical harmonics analysis/synthesis programs shana and shsyn which use GMT (versions 4.X.X. and 6+) respectively Netcdf grd files for I/O, and interoperate with the seismic tomography models from Becker & Boschi (2002). The most recent code is here:

• shansyn code repository on GitHub
• tar balls of the code
  • latest release
  • shansyn.1.6-070923.tgz (new version, all GMT versions)
  • shansyn.1.2.7-012818.tgz (old version, only for GMT4)

Reference:

• Becker, T. W. and Boschi, L.: A comparison of tomographic and geodynamic mantle models, Geochem., Geophys., Geosys., 3(1), 1003, doi:10.1029/2001GC000168, 2002. (PDF)

Unified Geodynamics Earth Science Computing Environment (UGESCE)

A VirtualBox Linux installation (i.e. a repacking of existing software) that contains fstrack, SEATREE (2016 version, new one available), iGMT, GMT, and a range of Earth science software and datasets.

• UGESCE web site (outdated, 2016)

iGMT: Graphical user interface and GMT script generator

A graphical user interface for the Generic Mapping Tools (GMT) written in Tcl/Tk and used around the world to teach GMT. Useful to learn and teach GMT, and to generate simple example scripts for general Earth science map making.

• iGMT on GitHub
• current release
• igmt_v1.2-20080724.tgz

Reference:

• Becker, T. W. and Braun, A.: New program maps geoscientific datasets interactively. EOS Transactions AGU, 79 (42), 505-506, 1998. (PDF)

• bin_catalog source code on GitHub
• the code contains C language rewrites of software by
  • Michael, A.J., 1984. Determination of stress from slip data: Faults and folds, J. Geophys. Res. 89, 11.517-11.526.
  • Michael, A.J., 1987. Use of focal mechanisms to determine stress: A control study, J. Geophys. Res. 92, 357-368.
  • Vavryčuk, V., 2014. Iterative joint inversion for stress and fault orientations from focal mechanisms, Geophysical Journal International, 199, 69-77.

Also contains a range of AWK scripts to convert polygons, earthquake hypocenters, and focal mechanisms/beach balls into VTK.

• recent tar file with grd2vtk scripts
• GitHub repository for grd2vtk tools

Minor scripts, tools, patches, etc.

• plot_ternary: A GMT script to generate ternary diagrams
  
• ssrf_surface: spherical GMT surface gridder
  A GMT surface-like wrapper routine for SSRFPACK by Renka which allows to interpolate irregular data on a sphere. (Now included, differently, in GMT itself.) You can read more on my iGMT project page, or download the latest version:
  • ssrf.010506.tar.gz
• ndk2meca.awk: Gawk script to convert ndk to psmeca
  The script will convert Global CMT ndk format output to GMT psmeca format and select mechanism types.
• CitcomS and CitcomCU extensions (e.g. GMT grd support, VTK formatting, anisotropic viscosity, strain tracking) are all updated in the Citcom CIG GitHub repository.