Updated Cariaco Basin 14C Calibration Dataset from 0–60 cal kyr BP

ABSTRACT We present new updates to the calendar and radiocarbon (14C) chronologies for the Cariaco Basin, Venezuela. Calendar ages were generated by tuning abrupt climate shifts in Cariaco Basin sediments to those in speleothems from Hulu Cave. After the original Cariaco-Hulu calendar age model was published, Hulu Cave δ18O records have been augmented with increased temporal resolution and a greater number of U/Th dates. These updated Hulu Cave records provide increased accuracy as well as precision in the final Cariaco calendar age model. The depth scale for the Ocean Drilling Program Site 1002D sediment core, the primary source of samples for 14C dating, has been corrected to account for missing sediment from a core break, eliminating age-depth anomalies that afflicted the earlier calendar age models. Individual 14C dates for the Cariaco Basin remain unchanged from previous papers, although detailed comparisons of the Cariaco calibration dataset to those from Hulu Cave and Lake Suigetsu suggest that the Cariaco marine reservoir age may have shifted systematically during the past. We describe these recent changes to the Cariaco datasets and provide the data in a comprehensive format that will facilitate use by the community.


INTRODUCTION
Calibration of the 14 C timescale over its full length has advanced remarkably in recent decades. Several independent datasets now extend beyond~50 cal kyr BP (Hughen et al. 2006;Bronk Ramsey et al. 2012;Cheng et al. 2018), displaying detailed shared structure and improving confidence in the general trends of the 14 C calibration curve, though differences still remain (Reimer et al. 2013 in this issue). Ocean Drilling Program (ODP) Site 1002, drilled in the Cariaco Basin off the coast of Venezuela, possesses high deposition rate sediments that are intermittently laminated, with annual varves during the last deglaciation and Holocene, but predominantly bioturbated sediments throughout the last Glacial period. Cariaco sediments contain detailed paleoclimate records showing abrupt millennial-scale climate shifts throughout the last Glacial period and beyond (Peterson et al. 2000). These millennial climate changes display abrupt shifts between quasi-stable states that occurred on timescales of decades or less (Hughen et al. 1996). Similar changes have been reported from other high-resolution archives, including Greenland ice cores (Johnsen et al. 1992) and speleothems from China (Wang et al. 2001), that possess calendar-age chronologies from layer-counting or U/Th dating. General circulation model (GCM) results show that reducing North Atlantic Meridional Overturning Circulation (AMOC), believed to have caused abrupt changes in high-latitude temperature, also result in "downstream" shifts in tropical Atlantic trade wind strength and monsoon rainfall over Asia (e.g., Rind et al. 1986;Manabe and Stouffer 1997). The Cariaco high-resolution paleoclimate records have therefore been "tuned" to the Greenland ice core (Meese et al. 1997) and Hulu Cave speleothem (Wang et al. 2001) records to provide calendar age models. Such paleoclimate tuning is based upon the identification of a set of tie-points, relating to shared abrupt palaeoclimatic events and assumed to be approximately contemporaneous, in both the Cariaco record and the Greenland/Hulu records from which chronological information can be transferred. The resulting Cariaco Basin 14 C calibration datasets (Hughen et al. 2004a(Hughen et al. , 2006 were used to help build the IntCal04, IntCal09 and IntCal13 calibration curves (Reimer et al. 2004(Reimer et al. , 2009(Reimer et al. , 2013Hughen et al. 2004b).
Here, we present the results of tuning Cariaco Basin sediment reflectance data to revised, highresolution Hulu Cave δ 18 O records (Cheng et al. 2016) with updated U/Th chronologies (Cheng et al. 2018). This improvement in Hulu Cave resolution enables us to re-evaluate and improve our ties. Given these updated ties, we then use the same Gaussian Process methodology  as used for the IntCal13 update, to propagate uncertainties in both the precise synchronicity of the selected ties and the Hulu U/Th dating through to the new Cariaco timescale. The updated Cariaco Basin 14 C calibration dataset is used to help build the IntCal20, SHCal20, and Marine20 14 C calibration curves (Reimer et al. 2020 in this issue;Hogg et al. 2020 in this issue;Heaton et al. 2020b in this issue), and the data are provided here with replicate 14 C measurements reported as averages as well as individual dates.
Use of marine records in reconstructing atmospheric radiocarbon is complicated by the systematic 14 C depletion within surface-ocean environments known as a marine reservoir age (MRA), which may vary over time. To incorporate the Cariaco 14 C record into IntCal20, changes in the MRA for the Cariaco Basin through time were modeled as an additive spline and estimated during curve construction by comparison with all of the other constituent IntCal20 datasets (Heaton et al. 2020a in this issue). This additive approach enables the Cariaco Basin to contribute strength to the IntCal20 curve in the case where 14 C features are shared with other IntCal20 datasets; while reducing the impact of potentially spurious variation seen only in a single set. Further, it provides insight of potential independent interest into the magnitude and duration of MRA variability within the Cariaco Basin.
Notation: In this paper, both 14 C ages and calendar ages are reported relative to AD 1950 (= 0 BP, before present). To distinguish between them, 14 C ages are given in units " 14 C yrs BP" while calendar/calibrated ages are given as "cal yr BP" or cal kyr BP (thousands of calibrated years before present). All quoted uncertainties on values, e.g. 420 ± 50 14 C yrs, refer to the 1σ level.

Calendar Chronology
The Cariaco Basin ODP sediment record has no calendar age scale of its own, and rapid climate shifts in Cariaco and other, absolute-age, paleoclimate archives have been used to transfer calendar age models to the Cariaco record. Originally, the calendar-age chronology for the Site 1002 sediment record was constructed by aligning climate changes with the annual layer counted GISP2 Greenland ice core δ 18 O record (Meese et al. 1997) (Figure 1), based on evidence for synchronous timing of rapid climate shifts within the timescale of one or two decades (Hughen et al. 1996(Hughen et al. , 2004c. Tie-points were based on rapid stadialinterstadial transitions ("warmings") in Greenland, which have counterparts (reduced trade winds and increased precipitation) in the tropical North Atantic region (Hughen et al. 2004a). In Hughen et al. (2006), the Cariaco record, including transitions both into and out of millennial events, was tuned to analogous hydrological changes in U/Th-dated speleothem δ 18 O records from Hulu Cave (Wang et al. 2001). The Hulu Cave δ 18 O records possessed relatively low temporal resolution, and there was some uncertainty in identifying precisely analogous signals in the two records, especially in the late Glacial period~15-25 cal kyr BP (Hughen et al. 2006). The Hulu Cave speleothem records have since been updated with higher temporal resolution δ 18 O measurements (Cheng et al. 2016) and a denser sequence of U/Th dates (Cheng et al. 2018), and the revised records show even stronger agreement with Cariaco Basin paleoclimate ( Figure 2).

Sediment Depth Scale
The depth scale previously used for Cariaco Basin paleoclimate and 14 C records used in IntCal was based on ODP Site 1002 Hole D, from which the majority of samples for 14 C dating were obtained. A core-break occurred in Hole D near 14 mbsf, which was bridged using continuous sediments from adjacent Hole E (Hughen et al. 2004a). For consistency with sample identifications at ODP sediment core repositories, the Hole D depth scale was maintained, despite creating an anomaly in the age-depth curve around 14 mbsf (30-32 cal kyr BP) ( Figure 3). Here, we have revised the depth scale by restoring 56 cm of missing sediment to the Hole D depth scale between 14.2 and 14.76 mbsf. Both original and adjusted depth scales are reported in Table 1.

Tie-Points
In this study, the tie points used to align the Cariaco Basin and Hulu Cave paleoclimate records were selected primarily at the rapid transitions between the tropical analogs of "stadial" and "interstadial" millennial climate events (i.e, "warmings" at high latitudes and increased precipitation in the tropics) ( Figure 2). These are generally the same tie points used in Hughen et al. (2006), except that most of the less-rapid transitions between "interstadial" and "stadial" events (i.e, high-latitude "coolings" or tropical "dryings") have been eliminated due to the lower precision of the matches. We note that during the later Glacial period, from 15-25 cal kyr BP, the tie points are less distinct, with greater variability in the Cariaco record than in Hulu. The remaining individual tie points were evaluated to identify anomalous  Close-up view of Cariaco-Hulu calendar age plotted versus Cariaco Basin sediment depth. 56 cm of sediment missing from a core-break below 14 mbsf have been restored to the original depth scale (gray), and the revised age-depth model (blue) showed an anomalous tie point around 32 cal kyr BP (red). This tie-point has been removed from the current age model. changes in sedimentation rate resulting from the revised depth scale. All identified tie-points were assumed to represent approximately synchronous events in the Hulu Cave and Cariaco Basin records, with an uncertainty in this contemporaneity of ±100 cal yrs (1σ) for each tie incorporated into chronology construction.
One tie point, related to a relatively minor climate event correlated with Greenland Interstadial (GIS) 5.1 immediately preceding Heinrich Event (HE) 3 (Rasmussen et al. 2014) (Figure 2), created a sharp change in sedimentation rate at 30-32 cal kyr BP ( Figure 3) and likely contributed to anomalously old calendar ages within this section in the Cariaco record and IntCal13 (Reimer et al. 2013). This spurious link has therefore been removed, eliminating the anomalous sedimentation rates and resulting in a smoother age-depth curve overall ( Figure 4). The tie points were used as input into a Gaussian process model ) to create a continuous calendar age-depth model for the Cariaco Basin with associated age uncertainty estimates that propagate both our uncertainty in the Hulu timescale and the precise contemporaneity of the chosen ties.

C Chronology
The same 14 C data set used previously to help constrain 14 C calibration in IntCal09 and IntCal13 (Reimer et al. 2009(Reimer et al. , 2013) is plotted here versus the updated Cariaco Basin calendar chronology ( Figure 5). In this plot, we have adjusted the Cariaco 14 C data sets based on a constant marine reservoir age (MRA) of 420 ± 50 14 C yrs . In addition to 14 C age, the absolute concentration of 14 C (Δ 14 C), corrected for isotopic fractionation and calendar age and reported relative to 1890 atmospheric values (Stuiver and Polach 1977), is plotted to amplify potential variability through time. The Cariaco data are plotted together with Hulu Cave (Cheng et al. 2018) and Lake Suigetsu (Bronk Ramsey et al. 2012) 14 C calibration data for comparison. The Hulu data set has been corrected for dead carbon fraction (DCF) under an assumption of a constant DCF of 450 ± 70 14 C yrs to enable comparison.
All three datasets reveal close agreement, including prominent structures previously reported for the Glacial and Deglacial periods (e.g., Hughen et al. 2006;Bronk Ramsey et al. 2012;Cheng et al. 2018). These include the Laschamps geomagnetic minimum and period of rapid 14 C production from~42-41 cal kyr BP, large oscillations from 40-30 cal kyr BP also attributed to geomagnetic variability (Laj et al. 2002;Hughen et al. 2004a), and a gradual overall decline after~25 cal kyr BP due to increasing geomagnetic field strength and the end of Glacial climate conditions ( Figure 5). The Lake Suigetsu data show a relatively higher degree of scatter than the other two records, and for clarity we show the Cariaco and Hulu data sets alone (Figure 6), although each plot is replicated including the Suigetsu data in the Supplementary Information (Figures S1-S7).

Marine Reservoir Age
Comparing the newly updated Cariaco Basin 14 C dataset to Hulu Cave 14 C in Figures 5 and 6, both plotted under the assumption of a constant level of MRA/DCF 14 C depletion, does however reveal several periods where apparent offsets exist, mostly consisting of younger 14 C ages (higher Δ 14 C values) for the Cariaco Basin data. These are most evident around 18, 23, and 30 cal kyr BP. Potential explanatory factors for these periodic offsets include periodic misalignment of the Cariaco tuning; uncertainties in the calendar age chronologies which are not indicated in Figure 5 but are incorporated in construction of the IntCal20 curve; and/or that an assumption of constant MRA/DCF depletion over time is not appropriate.
Both systems sample atmospheric 14 C indirectly. The world's ocean sequesters 14 C away from replenishment at the surface long enough for significant decay to occur, and thus surface waters appear older than the contemporary atmosphere with a global average, present day, MRA of 400 years, although with substantial spatial variability (Reimer and Reimer 2001). The Cariaco Basin MRA has been measured in recent, pre-bomb sediments of known calendar age (Hughen et al. 1996), and by comparison with overlapping tree rings (Kromer and Spurk 1998;Kromer et al. 2004) further back in time during Deglaciation, to yield a value of 420 ± 50 14 C yrs . Speleothems also exhibit 14 C depletion as groundwater feeding Hulu Cave filters through overlying soil and bedrock, absorbing " 14 Cdead" CO 2 and giving the growing speleothems a DCF equivalent to~450 ± 70 14 C yrs (Southon et al. 2012).
Over the period from 11-14 cal kyr BP, it is possible to obtain an estimate of MRA within the Cariaco Basin that is independent of tuning using the varved section of the record . This section of the Cariaco record, with its varve-counting chronology, can be compared with radiocarbon measurements from Northern Hemisphere (NH) tree-rings taken from the IntCal20 database  in this issue) that provide a direct record of atmospheric radiocarbon. Consequently, this estimate is also independent of the Hulu Cave's DCF. In Figure 7 we plot the resultant MRA estimate, obtained via a Bayesian spline. As we see, the MRA is relatively constant despite past abrupt climatic changes, such as the end of the Younger Dryas (YD) cold event around 11.5 cal kyr BP, but shows a very clear reduction in MRA during the onset of the YD around 12.8 cal kyr BP where it drops by several hundred years to near zero. For the Hulu cave speleothem, an analogous comparison to NH atmospheric trees indicates that the DCF with Hulu has remained highly constant (see Figure 6 in Heaton et al. 2020b in this issue) over the entire 11-14 cal kyr BP period. Figure 4 Cariaco-Hulu age-depth models shown over the entire time scale. The revised agedepth curve (blue) is smoother overall and shows fewer abrupt changes than previously (gray).
Figure 5 14 C calibration dataset from the Cariaco Basin (blue) plotted versus Hulu Cave (red) and Lake Suigetsu (green) datasets. (upper) Datasets shown as age-age plots. (lower) Data plotted as original 14 C concentration (Δ 14 C), following corrections for calendar age decay and fractionation. The datasets show fairly close agreement from~32-10 cal kyr BP, but increased scatter prior to that time. All 14 C uncertainties plotted as 1σ and include estimated uncertainties in MRA and DCF for Cariaco Basin and Hulu Cave, respectively.

K A Hughen & T J Heaton
The above independent comparisons of the varved Cariaco and Hulu records against NH treerings that provide direct sampling of atmospheric 14 C suggest that the periodic (and short lived) offsets between Hulu and the tuned Cariaco are likely due to MRA changes as opposed to mistuning.

Incorporating Cariaco Basin into the IntCal20 Curve
To incorporate the Cariaco Basin 14 C dataset into the IntCal20 curve, we need to model its MRA changes over time. For the other marine datasets used for the estimation of IntCal20, such MRA estimates were obtained using the three-dimensional Large Scale Geostrophic Ocean General Circulation Model (LSG OGCM) (Butzin et al. 2017 in this issue). However, this model does not have the resolution to simulate the unique topography of the Cariaco Basin, as evidenced by its inability to recreate the reduction in MRA seen around 12.8 cal kyr BP in Figure 7.
Consequently, we adopt a different approach to MRA modeling for the Cariaco Basin. We consider the Cariaco MRA as an unknown 14 C offset, via an additional Bayesian spline, and seek to adaptively estimate its value during curve construction itself based upon comparison to the other datasets 1 . At each iterative step of our Markov Chain Monte Carlo curve construction approach (Heaton et al. 2020a Figure 7 Cariaco MRA estimate obtained via a Bayesian spline when comparing the varved section of the record from 11-14 cal kyr BP to atmospheric tree-ring data from the IntCal20 database. The MRA estimate we obtain over this period is independent from both tuning and Hulu Cave's DCF. Each plotted point corresponds to the observed offset between an individual 14 C determination from the Cariaco varved record and the IntCal20 curve, which in this period is based only on tree-ring samples. The observed offsets are shown with 1σ uncertainties representing the calendar age uncertainty of the varved Cariaco 14 C determinations, and 1σ uncertainties in 14 C which incorporate both the uncertainty in the Cariaco 14 C sample and the IntCal20 curve at that calendar age. The MRA is relatively constant across the abrupt climate change at the YD termination (~11.5 cal kyr BP) but shows a clear reduction to near zero during the YD onset (~12.8 cal kyr BP). 1 For completeness, as for all the marine datasets, when incorporating the Cariaco samples into IntCal20 we further add an independent 14 C uncertainty to each 14 C determination to account for any additional short-term uncertainty not represented by our chosen MRA model. The size of this additional MRA variability is dataset dependent and based upon the increase needed to make that specific marine dataset consistent with the independent atmospheric NH treering determinations in the 11-14 cal kyr BP time period. See  this issue) for more detail.
Cariaco Basin 14 C Calibration Dataset 1029 estimate for the atmospheric 14 C calibration curve and our estimate for the Cariaco MRA based upon its observed offset from that curve.
This adaptive approach does introduce some questions regarding identifiability of 14 C changes within the Cariaco Basin record-to what extent is a feature seen within the Cariaco record due to a change in MRA within the basin, and to what extent an atmospheric event? Typically, through our chosen approach, we would expect features only seen within the Cariaco Basin to be considered a basin-specific MRA effect; while those additionally seen in other records would be considered atmospheric and reproduced in our atmospheric IntCal20 estimate.
As such, the Cariaco Basin plays a somewhat unique role in the construction of the IntCal20 curve. The additional flexibility provided by our adaptive modeling of the MRA will hopefully resolve the periods of systematic offsets between the Cariaco Basin and Hulu datasets (as seen in Figures 5 and 6) that may be due to short term MRA changes. However, the inclusion of Cariaco Basin still provides us with a better ability to recognize and retain shared features seen across multiple records and hence provide a more robust atmospheric 14 C curve-the intended aim of the IntCal project.

Inferred MRA Changes within Cariaco Basin over Time
As a consequence of this adaptive approach to modeling Cariaco Basin's MRA, alongside construction of the IntCal20 curve, we simultaneously obtain a posterior estimate of the MRA within the Cariaco Basin over time. This is shown in Figure 8a and is of potential independent interest. We discuss this posterior estimate below, illustrating the changing offset in the Cariaco Basin compared to Hulu cave. In building the IntCal20  in this issue) and Marine20 (Heaton et al. 2020b in this issue) curves, marine datasets were corrected for changes in MRA through time. Global MRAs were calculated using an ocean GCM (OGCM)  in this issue) and applied to individual datasets according to region. However, the modeled MRA variability for the Cariaco Basin appears much too high (Figure 8b), possibly due the OGCM's inability to capture the unique bathymetry and hydrographic setting of the basin. We therefore use the posterior estimate of the MRA to correct the Cariaco 14 C calibration curve back through time  in this issue). Due to the density of Hulu Cave measurements in the combined IntCal20 database, and the assumption of an approximately constant DCF over time within each Hulu speleothem in IntCal20 construction, this 14 C offset between the Cariaco and Hulu data is the primary driver for the obtained MRA estimate.
Looking closely at both our posterior MRA estimate and the apparent deviations between Cariaco Basin and Hulu Cave 14 C through time, we can evaluate the timing of 14 C change relative to rapid climate shifts in the paleoclimate record. As discussed above in Figure 7 comparing against trees that directly sample the atmosphere, the YD onset is clearly seen to coincide with an apparent reduction of Cariaco MRA by several hundred years. This is also seen when the Cariaco data is compared against the Hulu data alone (Figure 9). Similarly, Cariaco MRA appears to drop by several hundred years during the time of Heinrich Event 1 (Figure 10). The MRA decrease follows event H1b and persists throughout the transition into event H1a. However, there is no apparent change in MRA during the climatically more severe H1a event (Figure 10). The calendar age of this MRA decrease is bracketed by tie points at the end of H1b and the onset of H1a (17,148 and 16,150 cal yr BP, respectively), and falls within the older half of an anomalous "grey layer" that bears evidence of increased freshwater input and fluvial runoff (Yurco 2010).  Figure 8 (a) Posterior mean estimate (shown in black with the 95% credible interval) of Cariaco Basin's MRA through time, obtained simultaneously to IntCal20 curve construction. This MRA estimate is obtained adaptively as an additive Bayesian spline at the same time as the IntCal20 curve is constructed based upon the observed offset between Cariaco and other individual 14 C calibration datasets used in IntCal20. The procedure captures the same timing, direction and relative magnitude of Cariaco MRA changes as seen in direct comparisons with Hulu Cave data (following figures), likely due to the density of Hulu Cave measurements in the combined IntCal20 datasets and the assumption of constant DCF depletion. (b) Cariaco posterior mean MRA estimate shown together with the estimates for the nearest open-ocean site to the Cariaco Basin obtained by the LSG OGCM under three different climate scenarios when driven by a preliminary estimate of atmospheric Δ 14 C obtained from the Hulu cave speleothems alone  in this issue). Regional first-order Hulu-based LSG OGCM estimates were used to incorporate the other marine datasets into IntCal20  in this issue) but do not appear appropriate for the Cariaco Basin, as they overestimate the variability and magnitude of MRA values through time.      However, in order to explain the Cariaco 14 C age plateau as resulting from rapidly increased sedimentation (from terrestrial runoff),~63 cm of sediment would have had to be deposited nearly instantaneously soon after event H1b. In the absence of evidence for such an extreme sedimentation event, this MRA decrease around the time of H1 appears real. Although the precise calendar age of this feature has relatively high uncertainties from the tuning process, the relative timing of the MRA decrease well before the climatic H1a stadial suggests that increased atmospheric Δ 14 C associated with reduced AMOC (e.g., Hughen et al. 1998) may not have played a role here. Rather, the apparent reduction in Cariaco MRA is likely associated with decreased mixing of deep waters due to increased stratification in surface waters of the tropical North Atlantic or Cariaco Basin region.

K A Hughen & T J Heaton
The general timing of periods of reduced Cariaco MRAs reveals that an episode occurred approximately every 6-7 cal kyrs, suggestive of the intervals between Heinrich Events (Rasmussen et al. 2014). To test for a potential link between climate and MRA changes, we examine the detailed evidence for reduced Cariaco MRA during earlier Heinrich events (Figures 11-14). Although we do in fact observe reduced Cariaco MRA near the occurrence of Heinrich Events, the timing in general is off, with reduced MRA beginning up to 2000 cal yrs after the onset of the climatic events (Figures 11-14). In addition to reduced MRA, we also see evidence for a period of increased MRA between 28 and 25 cal kyr BP (Figure 15). Although there was likely an overall increase in MRA during the Glacial period, due to reduced carbon reservoirs in the atmosphere, biosphere and shallow ocean sediments (e.g., Hughen et al. 2004a), this would not explain a short interval of increased MRA from 25-28 cal kyr BP. This interval instead appears to be part of the cyclical pattern of Cariaco MRA changes ( Figure 8a) and may represent a period of relatively weak stratification and strong deep mixing (upwelling).
We note that, due to potential variability in both Hulu DCF and Cariaco MRA through time, the magnitude and even exact timing of changes in the Cariaco Basin MRA cannot be determined with precision. In particular, our current estimates are heavily based upon the assumption that the DCF within Hulu Cave remained approximately constant from 0-55 cal kyr BP. If global MRA were adjusted upward due to Glacial boundary conditions, it would increase the apparent magnitude of periodic Cariaco MRA reductions, and even slightly broaden the durations of the events. Regardless, the general pattern of quasi-periodic cycles of MRA variability appears to be robust (Figures 8, 16). Overall, the disagreement in relative timing between reduced MRA in the Cariaco Basin and Heinrich Events (Figure 16) precludes reduced AMOC during H events as a potential cause. The data suggest cyclic decreases in Cariaco MRA to values below 100 14 C years, despite Glacial boundary conditions that should have somewhat increased surface ocean MRA around the globe (Hughen et al. 2004a;Butzin et al. 2020 in this issue). These low reservoir ages are an extraordinary finding without analog in the modern ocean and are difficult to reconcile with physical or geochemical models of ocean circulation. At a minimum, these results suggest periodic swings in surface stratification in the Cariaco Basin region with times of greatly reduced mixing of 14 C-depleted "older" water from below. However, a clearer vision of true atmospheric 14 C throughout the Glacial period may be necessary in order to draw a firmer conclusion. The Lake Suigetsu record provides such an atmospheric record, but the magnitude of scatter in the current dataset precludes the identification of MRA shifts, and extension of dendrochronologies such as from floating Kauri sequences (Turney et al. 2010) may be the best solution.

CONCLUSIONS
The 14 C calibration data set from intermittently varved sediments of the Cariaco Basin has been substantially updated, with changes to calendar and 14 C chronologies. Revising the depth scale to reflect sediment missing from a core break provided a more accurate view of age-depth changes over the length of the record. A re-evaluation of tie points between existing Cariaco and newly updated Hulu Cave paleoclimate records resulted in the removal of one anomalous tie point, and better agreement between 14 C datasets. Close comparison of Cariaco and Hulu 14 C data shows periodic episodes of reduced Cariaco MRA, but coeval paleoclimate records show no direct relationship with the occurrence of Heinrich Events in the high-latitude North Atlantic region. The determination of the exact magnitude and duration of MRA variability is sensitive to long-term background changes in both Cariaco MRA and Hulu Cave DCF. However, an apparently cyclic pattern of Cariaco MRA decreases, possibly driven by increased surface ocean stratification, appears to be a robust feature of the updated record. An atmospheric 14 C dataset with less scatter than the Lake Suigetsu record, potentially from floating Kauri sequences, is needed to determine the precise history of MRA variability in the Cariaco Basin and elsewhere. Nevertheless, the Cariaco 14 C calibration dataset presented here incorporates several improvements over previous versions and more accurately reflects changes in both marine and atmospheric 14 C through time.