Marine sediments are typically subject to bioturbation from benthic macrofauna, which disrupt sediments and can have a substantial impact on the amplitude and timing of climate proxy signals recorded in marine sediment cores. This study explores a computational approach to quantify and remove these post-depositional effects from marine climate proxy records. The bioturbation process is modeled as a linear time-invariant filter building upon prior work (Guinasso and Schinck, 1975). Using forward modeling we find that, given modeled mixing layer depths ≥ 5 cm, simulated centennial scale climate variability is generally not preserved even when the sedimentation rate in our model is above 15 cm/kyr. On the scale of ice ages (10⁴-10⁵ yr), the observed effects of bioturbation have the smallest impact, considering the event scale and typical sedimentation rates and mixed layer depth in the deep sea. For millennial scale events, the signal attenuation strongly depends on the event scale and specific bioturbation and sedimentation parameters. To account for the bioturbation effect on a given climate proxy series, a deconvolution method is proposed. We apply the approach to three individual benthic foraminifera oxygen isotope records to reconstruct ${\delta }^{18}\text{O}$ across the last interglacial (Marine Isotope Stage 5e; MIS 5e), and to assess the magnitude of the deglacial shift between MIS 6 and MIS 5e. In the highest resolution records (Site MD952042, Site 1012), the deconvolution results reveal an increase in the amplitude of the recovered deglacial ${\delta }^{18}\text{O}$ shift by approximately 0.1‰, with implications for reconstructed global average sea level variations and temperature estimates from benthic foraminifera ${\delta }^{18}\text{O}$. The recovered signal at Site 677 suggests high amplitude ${\delta }^{18}\text{O}$ variability during MIS 5e, but exhibits larger uncertainty due to a lower sampling frequency, lower sedimentation rate, and an inferred higher noise level. While uncertainties on the deconvolution results can be large, the results provide a new view of the potential impact of bioturbation on paleoclimate reconstruction. Future work can reduce the uncertainty in the deconvolution estimates, through quantitative integration of additional constraints, provided by our knowledge of the climate and depositional systems.

海洋沉积物通常受到底栖大型动物的生物扰动，这会破坏沉积物，并可能对海洋沉积物核心中记录的气候替代信号的振幅和时间产生重大影响。这项研究探索了一种计算方法，可以从海洋气候代用记录中量化和消除这些沉积后的影响。生物扰动过程被建模为基于先前工作的线性时不变滤波器（Guinasso和Schinck，1975）。使用正演模型，我们发现，在模型化混合层深度≥5 cm的情况下，即使我们模型中的沉积速率高于15 cm / kyr，模拟百年尺度的气候变异性通常也不会保留。在冰河时代的规模上（10 ⁴ -10 ⁵yr），考虑到事件规模和典型的沉积速率以及深海的混合层深度，观察到的生物扰动影响最小。对于千禧年规模的事件，信号衰减在很大程度上取决于事件的规模以及特定的生物扰动和沉降参数。考虑到生物扰动对给定气候代理序列的影响，提出了一种反卷积方法。我们将该方法应用于三个单独的底栖有孔虫氧同位素记录以进行重建${\delta }^{\mathrm{18岁}}\text{Ø}$跨最后一个冰间期（海洋同位素第5e阶段； MIS 5e），并评估MIS 6和MIS 5e之间的冰期转变幅度。在最高分辨率的记录中（Site MD952042，Site 1012），反卷积结果表明，恢复的冰期振幅增加了。${\delta }^{\mathrm{18岁}}\text{Ø}$ 偏移约0.1‰，对重构的全球平均海平面变化和底栖有孔虫的温度估计有影响 ${\delta }^{\mathrm{18岁}}\text{Ø}$。站点677处恢复的信号表明幅度较高${\delta }^{\mathrm{18岁}}\text{Ø}$MIS 5e期间的可变性，但由于较低的采样频率，较低的沉降速率和推断的较高噪声水平而表现出较大的不确定性。尽管反卷积结果的不确定性可能很大，但这些结果为生物扰动对古气候重建的潜在影响提供了新的观点。通过我们对气候和沉积系统知识的了解，通过对其他约束条件进行定量整合，未来的工作可以减少反褶积估算中的不确定性。