pyBasinSubs

pyBasinSubs reconstructs how a sedimentary basin sank and filled through time by backstripping its layers. It's the Python version, and an extension, of my original R code, BasinSubs, which I wrote for our work on the Coastal Range of eastern Taiwan (Lai et al., 2022, Communications Earth & Environment).

pyBasinSubs brings that backstripping workflow into Python and goes further: Monte Carlo uncertainty, thermal history and maturity, 2-D flexural backstripping and section restoration, CSV templates and apps for your own data.

BasinSubs is the original set of R scripts from Larry Syu-Heng Lai’s University of Oregon research on the Coastal Range of eastern Taiwan. It is tied to that research: the data and results were made with collaborators for the 2022 Communications Earth & Environment paper. It stays archived as published, on Zenodo, under the BSD-3-Clause license.

• Porosity and age models. A global average porosity–depth curve, and linear age–thickness models that place paleo-water depth and relative sea level at each dated unit boundary.
• Subsidence analysis. Decompaction and backstripping of the northern and southern sections, with errors carried through, and geohistory plots.
• Data and results. Unit thickness, paleo-water depth, sea level, porosity–depth data, and the long-term uplift and subsidence rates of the arc crust.

pyBasinSubs is a separate, general Python package.