PyBacktrack 1.4 released

PyBacktrack 1.4 is now available, and adds support for generating paleobathymetry grids from submerged present-day crust.

PyBacktrack is available as a Python package and a Docker image. Installation instructions can be found online in the pyBacktrack documentation.

Changes since version 1.3:

  • Can now generate paleobathymetry grids:
    • Submerged oceanic and continental present-day crust is backtracked and reconstructed.
      • Similar to drill sites, but 2D instead of 1D.
    • Continental subsidence model obtains rifting period from builtin rift grids:
    • All submerged crust is assigned a single lithology, whereas drill sites have multiple stratigraphic layers.
    • Water depths are negative (below sea level) in paleobathymetry grids, but positive in backtracked drill site output.
    • Supports multiple CPUs to reduce running time.
  • Supplementary scripts can be installed:
    • Previously was `misc/` folder in Github repository, but is now installable.
    • Added script to preferentially merge paleobathymetry grids produced by pybacktrack with externally produced paleobathymetry grids (after adding dynamic topography to external grids).
    • Added script to generate rift start/end time grids from a global deforming model (used to generate builtin rift grids).
    • Added script to extract present day trenches from a global deforming model (used internally to avoid deep bathymetry near trenches).
    • Added script to convert sediment thickness into water thickness based on porosity assuming a single “average ocean floor sediment” lithology.
  • Backtracking:
  • Updated example notebooks:
  • Updated functions and classes (API reference):
    • New paleobathymetry functions to:
      • Generate paleobathymetry by reconstructing and backtracking sediment-covered crust through time,
      • write reconstructed paleobathymetry as NetCDF grids.
    • Can decompact at any age (not just ages at stratigraphic boundaries).
      • Assumes a constant sediment deposition rate within each stratigraphic layer.
    • Dynamic topography:
      • Now supports multiple point locations:
        • Previously supported just a single drill site location.
      • Slight change in how grids are interpolated (see notes under DynamicTopography.sample).
        • Note: there is only a difference when interpolating between grid times (no difference at grid times).
      • New class InterpolateDynamicTopography:
        • For just interpolating time-dependent dynamic topography mantle frame grids.
          • Used in above-mentioned script for merging paleobathymetry grids.
        • Existing class DynamicTopography both reconstructs and interpolates.

Python package: