pyBasinSubs

What is public today is BasinSubs (v0.1.0): R scripts Larry Syu-Heng Lai wrote for subsidence analysis of syn-orogenic strata in the Coastal Range of eastern Taiwan. The data and the results were produced with collaborators and used in the papers cited below.

R scripts

• Porosity-depth function.R calculates and plots a global average porosity–depth curve for sandy, muddy, and mean marine sediment.
• Age PB dSL-N.R and Age PB dSL-S.R fit a linear age–thickness model for the composite columns of the northern and southern Coastal Range. The columns are described in Lai et al. (2021) and Lai et al. (2022). The age models place foraminiferal paleobathymetry in time, including by interpolation and extrapolation. Around each age constraint (a unit boundary) the scripts average paleobathymetry and relative sea level, with uncertainties, for the subsidence calculation. Relative sea level is from Miller et al., Science Advances (doi:10.1126/sciadv.aaz1346).
• Subsidence Analysis N.R and Subsidence Analysis S.R decompact and backstrip the northern and southern sections, carry the errors through, and plot geohistory diagrams. Inputs are the porosity–depth functions (Pofun.csv), relative sea level at the age constraints (SL_Unit_N.csv or SL_Unit_S.csv), paleobathymetry at those constraints (PB_Unit_N.csv or PB_Unit_S.csv), unit thickness between the dated heights (Ages_N_new.csv or Ages_S_new.csv), and the lithology fractions in each unit (Litho_N.csv or Litho_S.csv).

Data and result folders

• Thickness. Unit thickness and its uncertainty between age constraints, paleobathymetry and the average at each constraint, and lithology fractions.
• Sealevel. Relative sea level from Miller et al., reorganized, plus the average and uncertainty at each unit boundary.
• Backstripping. Compiled porosity–depth functions and grain densities for sandy, muddy, and mean marine sediment, and the global average curves used in the analysis. Sources are in the supplement of Lai et al. (2022).
• Output results. The linear age models, the subsidence analysis, and long-term uplift and subsidence rates of Luzon arc crust in eastern Taiwan.

Citation

Lai, L.S.-H., Dorsey, R.J., Horng, C.-S., Chi, W.-R., Shea, K.-S., and Yen, J.-Y. (2022) Extremely rapid up-and-down motions of island arc crust during arc-continent collision. Communications Earth & Environment, 3, 100. https://doi.org/10.1038/s43247-022-00429-2

Lai, L.S.-H., Dorsey, R.J., Horng, C.-S., Chi, W.-R., Shea, K.-S., and Yen, J.-Y. (2021) Polygenetic mélange in the retrowedge foredeep of an active arc-continent collision, Coastal Range of eastern Taiwan. Sedimentary Geology, 418, 105901. https://doi.org/10.1016/j.sedgeo.2021.105901

Code: Zenodo, doi:10.5281/zenodo.5823613. Data: figshare, doi:10.6084/m9.figshare.19350530.