Description: Palaeoclimate data hold the unique promise of providing a long-term perspective on climate change and as such can serve as an important benchmark for climate models. However, palaeoclimate data have generally been archived with insufficient standardisation and metadata to allow for transparent and consistent uncertainty assessment in an automated way. Thanks to improved computation capacity, transient palaeoclimate simulations are now possible, calling for data products containing multi-parameter time series rather than information on a single parameter for a single time slice. Efforts are underway to simulate a complete glacial–interglacial cycle using general circulation models (https://www.palmod.de/, last access: 6 May 2020), and to confront these simulations with palaeoclimate data, we have compiled a multi-parameter marine palaeoclimate data synthesis that contains time series spanning 0 to 130 000 years ago. We present the first version of the data product that focuses exclusively on time series for which a robust chronology based on benthic foraminifera δ18O and radiocarbon dating is available. The product contains 896 time series of eight palaeoclimate parameters from 143 individual sites, each associated with rich metadata, age–depth model ensembles, and information to refine and update the chronologies. This version contains 205 time series of benthic foraminifera δ18O; 169 of benthic foraminifera δ13C; 131 of seawater temperature; 174 and 119 of planktonic foraminifera δ18O and δ13C; and 44, 38 and 16 of carbonate, organic carbon and biogenic silica content, respectively. The data product is available in three formats (R, LiPD and netCDF) facilitating use across different software and operating systems and can be downloaded at https://doi.org/10.1594/PANGAEA.908831 (Jonkers et al., 2019). This data descriptor presents our data synthesis strategy and describes the contents and format of the data product in detail. It ends with a set of recommendations for data archiving.

Description: Deepwater circulation significantly changed during the last deglaciation from a shallow to a deep-reaching overturning cell. This change went along with a drawdown of isotopically light waters into the abyss and a deep ocean warming that changed deep ocean stratification from a salinity-to a temperature-controlled mode. Yet, the exact mechanisms causing these changes are still unknown. Furthermore, the long-standing idea of a complete shutdown of North Atlantic deepwater formation during Heinrich Stadial 1 (HS1) (17.5–14.6 kyr BP) remains prevalent. Here, we present a new compilation of benthic δ13C and δ18O data from the North Atlantic at high temporal resolution with consistent age models, established as part of the international PAGES working group OC3, to investigate deepwater properties in the North Atlantic. The extensive compilation, which includes 105 sediment cores, reveals different water masses during HS1. A water mass with heavy δ13C and δ18O signature occupies the Iceland Basin, whereas between 20 and 50°N, a distinct tongue of 18O depleted, 13C enriched water reaches down to 4000 m water depths. The heavy δ13C signature indicates active deepwater formation in the North Atlantic during HS1. Differences in its δ18O signature indicate either different sources or an alteration of the deepwater on its southward pathway. Based on these results, we discuss concepts of deepwater formation in the North Atlantic that help to explain the deglacial change from a salinity-driven to a temperature-driven circulation mode. Highlights • Spatial analyses of benthic δ13C and δ18O data from OC3 Atlantic compilation for HS1. • Heavy δ13C, light δ18O waters migrated into deep western North Atlantic basin during HS1. • Active deepwater formation between 30 and 60°N in the North Atlantic during HS1. • New concepts for transport of isotopically light δ18O into deep ocean. • Major contribution of North Atlantic waters to deglacial deep ocean stratification changes.

Description: Despite its great ecological importance, the main factors governing tree cover in tropical savannas as well as savanna-forest boundaries are still largely unknown. Here we address this issue by investigating marine sediment records of long-chain n-alkane stable carbon (δ13Cwax) and hydrogen (δDwax) isotopes from a core collected off eastern tropical South America spanning the last ca. 45 thousand years. While δ13Cwax is a proxy for the main photosynthetic pathway of terrestrial vegetation, tracking the relative proportion of C3 (mainly trees) versus C4 (mainly grasses) plants, δDwax is a proxy for continental precipitation, tracking the intensity of rainfall. The investigated core was collected off the mouth of the São Francisco River drainage basin, a tropical savanna-dominated region with dry austral autumn, winter and spring. On top of millennial-scale changes, driven by anomalies in the amount of precipitation associated with Heinrich Stadials, we identify a marked obliquity control over the expansion and contraction of tree and grass cover. During periods of maximum (minimum) obliquity, trees (grasses) reached maximum coverage. We suggest that maximum (minimum) obliquity decreased (increased) the length of the dry season allowing (hampering) the expansion of tree-dominated vegetation. Periods of maximum obliquity induced an anomalous heating (cooling) of the summer (winter) hemisphere that in combination with a delayed response of the climate system slightly increased autumn precipitation over the São Francisco River drainage basin, through a shift of the Intertropical Convergence Zone towards or further into the anomalously heated hemisphere. We found that atmospheric CO2 concentration has only a secondary effect on tree cover. Our results underline the importance of the dry season length as a governing factor in the long-term control of tree cover in tropical landscapes.

Description: Although the global environmental impact of Laurentide Ice-Sheet destabilizations on glacial climate during Heinrich Events is well-documented, the mechanism driving these ice-sheet instabilities remains elusive. Here we report foraminifera-based subsurface (~150 m water depth) ocean temperature and salinity reconstructions from a sediment core collected in the western subpolar North Atlantic, showing a consistent pattern of rapid subsurface ocean warming preceding the transition into each Heinrich Event identified in the same core of the last 27,000 years. These results provide the first solid evidence for the massive accumulation of ocean heat near the critical depth to trigger melting of marine-terminating portions of the Laurentide Ice Sheet around Labrador Sea followed by Heinrich Events. The repeated build-up of a subsurface heat reservoir in the subpolar Atlantic closely corresponds to times of weakened Atlantic Meridional Overturning Circulation, indicating a precursor role of ocean circulation changes for initiating abrupt ice-sheet instabilities during Heinrich Events. We infer that a weaker ocean circulation in future may result in accelerated interior-ocean warming of the subpolar Atlantic, which could be critical for the stability of modern, marine-terminating Arctic glaciers and the freshwater budget of the North Atlantic.

Description: We present a climatology of the near-sea-surface temperature (NSST) anomaly and the sea-ice extent during the Last Glacial Maximum (LGM, 23 000–19 000 years before present) mapped on a global regular 1∘×1∘ grid. It is an extension of the Glacial Atlantic Ocean Mapping (GLAMAP) reconstruction of the Atlantic NSST based on the faunal and floral assemblage data of the Multiproxy Approach for the Reconstruction of the Glacial Ocean Surface (MARGO) project and several recent estimates of the LGM sea-ice extent. Such a gridded climatology is highly useful for the visualization of the LGM climate, calculation of global and regional NSST averages, and estimation of the equilibrium climate sensitivity, as well as a boundary condition for atmospheric general circulation models. The gridding of the sparse NSST reconstruction was done in an optimal way using the Data-Interpolating Variational Analysis (DIVA) software, which takes into account the uncertainty in the reconstruction and includes the calculation of an error field. The resulting Glacial Ocean Map (GLOMAP) confirms the previous findings by the MARGO project regarding longitudinal and meridional NSST differences that were greater than today in all oceans. Taken at face value, the estimated global and tropical cooling would imply an equilibrium climate sensitivity at the lower end of the currently accepted range. However, because of anticipated changes in the seasonality and thermal structure of the upper ocean during the LGM as well as uneven spatial sampling, the estimated cooling and implied climate sensitivity are likely to be biased towards lower values.

Description: We present a global atlas of downcore foraminiferal oxygen and carbon isotope ratios available at https://doi.org/10.1594/PANGAEA.936747 (Mulitza et al., 2021a). The database contains 2106 published and previously unpublished stable isotope downcore records with 361 949 stable isotope values of various planktic and benthic species of Foraminifera from 1265 sediment cores. Age constraints are provided by 6153 uncalibrated radiocarbon ages from 598 (47 %) of the cores. Each stable isotope and radiocarbon series is provided in a separate netCDF file containing fundamental metadata as attributes. The data set can be managed and explored with the free software tool PaleoDataView. The atlas will provide important data for paleoceanographic analyses and compilations, site surveys, or for teaching marine stratigraphy. The database can be updated with new records as they are generated, providing a live ongoing resource into the future.

Description: We present the first version of the Ocean Circulation and Carbon Cycling (OC3) working group database, of oxygen and carbon stable isotope ratios from benthic foraminifera in deep ocean sediment cores from the Last Glacial Maximum (LGM, 23-19 ky) to the Holocene (〈10 ky) with a particular focus on the early last deglaciation (19-15 ky BP). It includes 287 globally distributed coring sites, with metadata, isotopic and chronostratigraphic information, and age models. A quality check was performed for all data and age models, and sites with at least millennial resolution were preferred. Deep water mass structure as well as differences between the early deglaciation and LGM are captured by the data, even though its coverage is still sparse in many regions. We find high correlations among time series calculated with different age models at sites that allow such analysis. The database provides a useful dynamical approach to map physical and biogeochemical changes of the ocean throughout the last deglaciation.

Description: Laboratory experiments showed that the isotopic fractionation of δ13C and of δ18O during calcite formation of planktic foraminifera are species-specific functions of ambient CO32- concentration. This effect became known as the carbonate ion effect (CIE), whose role for the interpretation of marine sediment data will be investigated here in an in-depth analysis of the 13C cycle. For this investigation, we constructed new 160 kyr long mono-specific stacks of changes in both δ13C and δ18O from either the planktic foraminifera Globigerinoides ruber (rub) or Trilobatus sacculifer (sac) from 112 and 40 marine records, respectively, from the wider tropics (latitudes below 38°). Both mono-specific time series Δ(δ13Crub) and Δ(δ13Csac) are very similar to each other, and a linear regression through a scatter plot of both data sets has a slope of ∼ 0.99 – although the laboratory-based CIE for both species differs by a factor of nearly 2, implying that they should record distinctly different changes in δ13C, if we accept that the carbonate ion concentration changes on glacial–interglacial timescales. For a deeper understanding of the 13C cycle, we use the Solid Earth version of the Box model of the Isotopic Carbon cYCLE (BICYLE-SE) to calculate how surface-ocean CO32- should have varied over time in order to be able to calculate the potential offsets which would by caused by the CIE quantified in culture experiments. Our simulations are forced with atmospheric reconstructions of CO2 and δ13CO2 derived from ice cores to obtain a carbon cycle which should at least at the surface ocean be as close as possible to expected conditions and which in the deep ocean largely agrees with the carbon isotope ratio of dissolved inorganic carbon (DIC), δ13CDIC, as reconstructed from benthic foraminifera. We find that both Δ(δ13Crub) and Δ(δ13Csac) agree better with changes in simulated δ13CDIC when ignoring the CIE than those time series which were corrected for the CIE. The combination of data- and model-based evidence for the lack of a role for the CIE in Δ(δ13Crub) and Δ(δ13Csac) suggests that the CIE as measured in laboratory experiments is not directly transferable to the interpretation of marine sediment records. The much smaller CIE-to-glacial–interglacial-signal ratio in foraminifera δ18O, when compared to δ13C, prevents us from drawing robust conclusions on the role of the CIE in δ18O as recorded in the hard shells of both species. However, theories propose that the CIE in both δ13C and δ18O depends on the pH in the surrounding water, suggesting that the CIE should be detectable in neither or both of the isotopes. Whether this lack of role of the CIE in the interpretation of planktic paleo-data is a general feature or is restricted to the two species investigated here needs to be checked with further data from other planktic foraminiferal species.

Description: Slowdowns of the Atlantic meridional overturning circulation during Heinrich stadials (HS) caused reductions in cross-equatorial heat transport, southward shifts of the Intertropical Convergence Zone and intensification of precipitation over eastern tropical South America. While these changes are well described, the associated spatial sea surface temperature (SST) patterns are still unclear. Here, we analyze proxy data to assess changes in ocean atmosphere interactions during HS over the western South Atlantic. Our SST proxy records show contrasting patterns between the tropical (warming) and the subtropical (no change) western South Atlantic during HS. We propose that the distinct SST behavior in the subtropics was associated with the cloud cover effect of the South Atlantic Convergence Zone over the upper ocean.

Description: Sediment core GeoB18530-1 was collected during MSM39 cruise in 2014 in the subpolar western North Atlantic. X-ray fluorescence (XRF) core-scanning measurements of GeoB18530-1 were performed with an Avaatech XRF Core Scanner at the MARUM - University of Bremen. For Mg/Ca analyses, ~100 individuals of the foraminiferal species N. pachyderma sin. were picked from sieved sediment samples of the size fraction 〉250 μm. Mg/Ca analyses were performed with an axial-viewing ICP-OES Varian 720 ES (Inductively Coupled Plasma Optical Emission Spectrometry) at GEOMAR Helmholtz-Centre for Ocean Research in Kiel. Stable isotopes (δ18O) of N. pachyderma sin. were measured with Thermo Scientific MAT 253 mass spectrometers equipped with an automated Kiel IV Carbonate Preparation Device at MARUM, University of Bremen and GEOMAR Helmholtz-Centre for Ocean Research in Kiel. We report 20 accelerator mass spectrometry radiocarbon dates (AMS14C ages) from samples of N. pachyderma sin. (19 samples) and one mixed sample (Globigerina bulloides and N. pachyderma sin.) from core GeoB18530-1 spanning the last ~ 35,000 years.