The assessment of age uncertainty in stratigraphically aligned records is a pressing need in paleoceanographic research. The alignment of ocean sediment cores is used to develop mutually consistent age models for climate proxies and is often based on the δ¹⁸O of calcite from benthic foraminifera, which records a global ice volume and deep water temperature signal. To date, δ¹⁸O alignment has been performed by manual, qualitative comparison or by deterministic algorithms. Here we present a hidden Markov model (HMM) probabilistic algorithm to find 95% confidence bands for δ¹⁸O alignment. This model considers the probability of every possible alignment based on its fit to the δ¹⁸O data and transition probabilities for sedimentation rate changes obtained from radiocarbon‐based estimates for 37 cores. Uncertainty is assessed using a stochastic back trace recursion to sample alignments in exact proportion to their probability. We applied the algorithm to align 35 late Pleistocene records to a global benthic δ¹⁸O stack and found that the mean width of 95% confidence intervals varies between 3 and 23 kyr depending on the resolution and noisiness of the record's δ¹⁸O signal. Confidence bands within individual cores also vary greatly, ranging from ~0 to >40 kyr. These alignment uncertainty estimates will allow researchers to examine the robustness of their conclusions, including the statistical evaluation of lead‐lag relationships between events observed in different cores.

We present a Python package geared toward the intuitive analysis and visualization of paleoclimate timeseries, Pyleoclim. The code is open‐source, object‐oriented, and built upon the standard scientific Python stack, allowing users to take advantage of a large collection of existing and emerging techniques. We describe the code's philosophy, structure, and base functionalities and apply it to three paleoclimate problems: (a) orbital‐scale climate variability in a deep‐sea core, illustrating spectral, wavelet, and coherency analysis in the presence of age uncertainties; (b) correlating a high‐resolution speleothem to a climate field, illustrating correlation analysis in the presence of various statistical pitfalls (including age uncertainties); (c) model‐data confrontations in the frequency domain, illustrating the characterization of scaling behavior. We show how the package may be used for transparent and reproducible analysis of paleoclimate and paleoceanographic datasets, supporting Findable, Accessible, Interoperable, and Reusable software and an open science ethos. The package is supported by an extensive documentation and a growing library of tutorials shared publicly as videos and cloud‐executable Jupyter notebooks, to encourage adoption by new users.

The progress of science is tied to the standardization of measurements, instruments, and data. This is especially true in the Big Data age, where analyzing large data volumes critically hinges on the data being standardized. Accordingly, the lack of community‐sanctioned data standards in paleoclimatology has largely precluded the benefits of Big Data advances in the field. Building upon recent efforts to standardize the format and terminology of paleoclimate data, this article describes the Paleoclimate Community reporTing Standard (PaCTS), a crowdsourced reporting standard for such data. PaCTS captures which information should be included when reporting paleoclimate data, with the goal of maximizing the reuse value of paleoclimate data sets, particularly for synthesis work and comparison to climate model simulations. Initiated by the LinkedEarth project, the process to elicit a reporting standard involved an international workshop in 2016, various forms of digital community engagement over the next few years, and grassroots working groups. Participants in this process identified important properties across paleoclimate archives, in addition to the reporting of uncertainties and chronologies; they also identified archive‐specific properties and distinguished reporting standards for new versus legacy data sets. This work shows that at least 135 respondents overwhelmingly support a drastic increase in the amount of metadata accompanying paleoclimate data sets. Since such goals are at odds with present practices, we discuss a transparent path toward implementing or revising these recommendations in the near future, using both bottom‐up and top‐down approaches.

We investigate the relationship between tropical Pacific and Southern Ocean variability during the Holocene using the stable oxygen isotope and magnesium/calcium records of cooccurring planktonic and benthic foraminifera from a marine sediment core collected in the western equatorial Pacific. The planktonic record exhibits millennial-scale sea surface temperature (SST) oscillations over the Holocene of ~0.5°C while the benthic δ¹⁸O_c document ~0.10‰ millennial-scale changes of Upper Circumpolar Deep Water (UCDW), a water mass which outcrops in the Southern Ocean. Solar forcing as an explanation for millennial-scale SST variability requires (1) a large climate sensitivity and (2) a long 400 year delayed response, suggesting that if solar forcing is the cause of the variability, it would need to be considerably amplified by processes within the climate system at least at the core location. We also explore the possibility that SST variability arose from volcanic forcing using a simple red noise model. Our best estimates of volcanic forcing falls short of reproducing the amplitude of observed SST variations although it produces power at low-frequency similar to that observed in the MD81 record. Although we cannot totally discount the volcanic and solar forcing hypotheses, we are left to consider that the most plausible source for Holocene millennial-scale variability lies within the climate system itself. In particular, UCDW variability coincided with deep North Atlantic changes, indicating a role for the deep ocean in Holocene millennial-scale variability.