Reply on EC1

Abstract

<strong class="journal-contentHeaderColor">Abstract.</strong> Previously developed software packages that generate probabilistic age models for ocean sediment cores are designed to use either age proxies (e.g., radiocarbon or tephra layers) or stratigraphic alignment (e.g., of benthic &delta;¹⁸O) and cannot combine age inferences from both techniques. Furthermore, many radiocarbon dating packages are not specifically designed for marine sediment cores and default settings may not accurately reflect the probability of sedimentation rate variability in the deep ocean, requiring subjective tuning of parameter settings. Here we present a new technique for generating Bayesian age models and stacks using ocean sediment core radiocarbon and benthic &delta;¹⁸O data, implemented in a software package named BIGMACS (Bayesian Inference Gaussian Process regression and Multiproxy Alignment of Continuous Signals). BIGMACS constructs multiproxy age models by combining age inferences from both radiocarbon ages and benthic &delta;¹⁸O stratigraphic alignment and constrains sedimentation rates using an empirically derived prior model based on 37 ¹⁴C-dated ocean sediment cores (Lin et al., 2014). BIGMACS also constructs continuous benthic &delta;¹⁸O stacks via a Gaussian process regression, which requires a smaller number of cores than previous stacking methods. This feature allows users to construct stacks for a region that shares a homogeneous deep water &delta;¹⁸O signal, while leveraging radiocarbon dates across multiple cores. Thus, BIGMACS efficiently generates local or regional stacks with smaller uncertainties in both age and &delta;¹⁸O than previously available techniques. We present two example regional benthic &delta;¹⁸O stacks and demonstrate that the multiproxy age models produced by BIGMACS are more precise than their single proxy counterparts.