Stable isotope record of Cibicidoides wuellerstorfi from Pleistocene sediments of DSDP Site 86-502 in the Caribbean Sea

Variations in the stable isotopic composition of benthic foraminifera from Deep-Sea Drilling Project (DSDP) site 502B in the Caribbean Sea are used to reconstruct Atlantic intermediate water circulation variability over the last 1.2 m.y. Comparison of this record with other North Atlantic benthic isotope records indicates that Atlantic intermediate water circulation was relatively enhanced during glacial maxima when North Atlantic deep water (NADW) production was reduced. However, a simple, compensatory relationship between intermediate and deepwater circulation is not apparent. Geochemical models have shown that such changes in ocean circulation can affect atmospheric CO2 levels by changing vertical nutrient and alkalinity profiles. The Delta delta13C difference between Caribbean site 502B and deep equatorial Pacific site 677 is highly coherent and in phase with ice volume. Like the delta18O record, there is an increase in amplitude (40%) and a large increase in 100 kyr power after 0.7 Ma. The 1.2? Delta delta13C amplitude scales to 70 ppm V in atmospheric CO2 using Boyle's (1986) box model result. The implied increase in CO2 amplitude after 0.7 Ma may suggest a positive feedback role in effecting the higher-amplitude climatic fluctuations which characterize the last 0.7 m.y.

Sedimentation rates calculated on surface sediment samples from different site of the Atlantic and Pacific Oceans (Table 1)

Radiocarbon ages on CaCO3 from deep-sea cores offer constraints on the nature of the CaCO3 dissolution process. The idea is that the toll taken by dissolution on grains within the core top bioturbation zone should be in proportion to their time of residence in this zone. If so, dissolution would shift the mass distribution in favor of younger grains, thereby reducing the mean radiocarbon age for the grain ensemble. We have searched in vain for evidence supporting the existence of such an age reduction. Instead, we find that for water depths of more than 4 km in the tropical Pacific the radiocarbon age increases with the extent of dissolution. We can find no satisfactory steady state explanation and are forced to conclude that this increase must be the result of chemical erosion. The idea is that during the Holocene the rate of dissolution of CaCO3 has exceeded the rain rate of CaCO3. In this circumstance, bioturbation exhumes CaCO3 from the underlying glacial sediment and mixes it with CaCO3 raining from the sea surface.

Water temperature reconstructions based on paired Uk'37 and Tex86 records, and paired surface and subsurface Mg/Ca records

This dataset contains the collection of available published paired Uk'37 and Tex86 records spanning multi-millennial to multi-million year time scales, as well as a collection of Mg/Ca-derived temperatures measured in parallel on surface and subsurface dwelling foraminifera, both used in the analyses of Ho and Laepple, Nature Geoscience 2016. As the signal-to-noise ratios of proxy-derived Holocene temperatures are relatively low, we selected records that contain at least the last deglaciation (oldest sample >18kyr BP).

Compilation of planktonic foraminiferal census counts from the Southern Hemisphere Oceans

We present an improved database of planktonic foraminiferal census counts from the Southern Hemisphere Oceans (SHO) from 15°S to 64°S. The SHO database combines 3 existing databases. Using this SHO database, we investigated dissolution biases that might affect faunal census counts. We suggest a depth/[DCO3]2- threshold of ~3800 m/[DCO3]2- = ~-10 to -5 µmol/kg for the Pacific and Indian Oceans, and ~4000 m/[DCO3]2- = ~0 to 10 µmol/kg for the Atlantic Ocean, under which core-top assemblages can be affected by dissolution and are less reliable for paleo-sea surface temperature (SST) reconstructions. We removed all core-tops beyond these thresholds from the SHO database. This database has 598 core-tops and is able to reconstruct past SST variations from 2° to 25.5°C, with a root mean square error of 1.00°C, for annual temperatures. To inspect dissolution affects SST reconstruction quality, we tested the data base with two "leave-one-out" tests, with and without the deep core-tops. We used this database to reconstruct Summer SST (SSST) over the last 20 ka, using the Modern Analog Technique method, on the Southeast Pacific core MD07-3100. This was compared to the SSST reconstructed using the 3 databases used to compile the SHO database. Thus showing that the reconstruction using the SHO database is more reliable, as its dissimilarity values are the lowest. The most important aspect here is the importance of a bias-free, geographic-rich, database. We leave this dataset open-ended to future additions; the new core-tops must be carefully selected, with their chronological frameworks, and evidence of dissolution assessed.

Data compilation on the biological response to ocean acidification: environmental and experimental context of data sets and related literature

The exponential growth of studies on the biological response to ocean acidification over the last few decades has generated a large amount of data. To facilitate data comparison, a data compilation hosted at the data publisher PANGAEA was initiated in 2008 and is updated on a regular basis (doi:10.1594/PANGAEA.149999). By January 2015, a total of 581 data sets (over 4 000 000 data points) from 539 papers had been archived. Here we present the developments of this data compilation five years since its first description by Nisumaa et al. (2010). Most of study sites from which data archived are still in the Northern Hemisphere and the number of archived data from studies from the Southern Hemisphere and polar oceans are still relatively low. Data from 60 studies that investigated the response of a mix of organisms or natural communities were all added after 2010, indicating a welcomed shift from the study of individual organisms to communities and ecosystems. The initial imbalance of considerably more data archived on calcification and primary production than on other processes has improved. There is also a clear tendency towards more data archived from multifactorial studies after 2010. For easier and more effective access to ocean acidification data, the ocean acidification community is strongly encouraged to contribute to the data archiving effort, and help develop standard vocabularies describing the variables and define best practices for archiving ocean acidification data.

(Table 2) Radiocarbon dates from lakes on Fildes Peninsula, King George Island

Precise relative sea level (RSL) data are important for inferring regional ice sheet histories, as well as helping to validate numerical models of ice sheet evolution and glacial isostatic adjustment. Here we develop a new RSL curve for Fildes Peninsula, South Shetland Islands (SSIs), a sub-Antarctic archipelago peripheral to the northern Antarctic Peninsula ice sheet, by integrating sedimentary evidence from isolation basins with geomorphological evidence from raised beaches. This combined approach yields not only a Holocene RSL curve, but also the spatial pattern of how RSL change varied across the archipelago. The curve shows a mid-Holocene RSL highstand on Fildes Peninsula at 15.5 m above mean sea level between 8000 and 7000 cal a BP. Subsequently RSL gradually fell as a consequence of isostatic uplift in response to regional deglaciation. We propose that isostatic uplift occurred at a non-steady rate, with a temporary pause in ice retreat ca. 7200 cal a BP, leading to a short-lived RSL rise of ~1 m and forming a second peak to the mid-Holocene highstand. Two independent approaches were taken to constrain the long-term tectonic uplift rate of the SSIs at 0.22-0.48 m/ka, placing the tectonic contribution to the reconstructed RSL highstand between 1.4 and 2.9 m. Finally, we make comparisons to predictions from three global sea level models.

Compilation of 220 sediment core d13C data from the glacial Atlantic Ocean

We compare a compilation of 220 sediment core d13C data from the glacial Atlantic Ocean with three-dimensional ocean circulation simulations including a marine carbon cycle model. The carbon cycle model employs circulation fields which were derived from previous climate simulations. All sediment data have been thoroughly quality controlled, focusing on epibenthic foraminiferal species (such as Cibicidoides wuellerstorfi or Planulina ariminensis) to improve the comparability of model and sediment core carbon isotopes. The model captures the general d13C pattern indicated by present-day water column data and Late Holocene sediment cores but underestimates intermediate and deep water values in the South Atlantic. The best agreement with glacial reconstructions is obtained for a model scenario with an altered freshwater balance in the Southern Ocean that mimics enhanced northward sea ice export and melting away from the zone of sea ice production. This results in a shoaled and weakened North Atlantic Deep Water flow and intensified Antarctic Bottom Water export, hence confirming previous reconstructions from paleoproxy records. Moreover, the modeled abyssal ocean is very cold and very saline, which is in line with other proxy data evidence.

(Table S1) Age determination of Atlantic Ocean sediment cores

Well-dated benthic foraminifer oxygen isotopic records (d18O) from different water depths and locations within the Atlantic Ocean exhibit distinct patterns and significant differences in timing over the last deglaciation. This has two implications: on the one hand, it confirms that benthic d18O cannot be used as a global correlation tool with millennial-scale precision, but on the other hand, the combination of benthic isotopic records with independent dating provides a wealth of information on past circulation changes. Comparing new South Atlantic benthic isotopic data with published benthic isotopic records, we show that (1) circulation changes first affected benthic d18O in the 1000-2200 m range, with marked decreases in benthic d18O taking place at ~17.5 cal. kyr B.P. (ka) due to the southward propagation of brine waters generated in the Nordic Seas during Heinrich Stadial 1 (HS1) cold period; (2) the arrival of d18O-depleted deglacial meltwater took place later at deeper North Atlantic sites; (3) hydrographic changes recorded in North Atlantic cores below 3000 m during HS1 do not correspond to simple alternations between northern- and southern-sourced water but likely reflect instead the incursion of brine-generated deep water of northern as well as southern origin; and (4) South Atlantic waters at ~44°S and ~3800 m depth remained isolated from better-ventilated northern-sourced water masses until after the resumption of North Atlantic Deep Water (NADW) formation at the onset of the Bølling-Allerod, which led to the propagation of NADW into the South Atlantic.

(Table 2) Age determination of Aegean Sea sediments

The modern Aegean Sea is an important source of deep water for the eastern Mediterranean. Its contribution to deep water ventilation is known to fluctuate in response to climatic variation on a decadal timescale. This study uses marine micropalaeontological and stable isotope data to investigate longer-term variability during the late glacial and Holocene, in particular that associated with the deposition of the early Holocene dysoxic/anoxic sapropel S1. Concentrating on the onset of sapropel-forming conditions, we identify the start of 'seasonal' stratification and highlight a lag in d18O response of the planktonic foraminifer N. pachyderma to termination T1b as identified in the d18O record of G. ruber. By use of a simple model we determine that this offset cannot be a function of bioturbation effects. The lag is of the order of 1 kyr and suggests that isolation of intermediate/deep water preceded the start of sapropel formation by up to 1.5 kyr. Using this discovery, we propose an explanation for the major unresolved problem in sapropel studies, namely, the source of nutrient supply required for export productivity to reach levels needed for sustained sapropel deposition. We suggest that nutrients had been accumulating in a stagnant basin for 1-1.5 kyr and that these accumulated resources were utilized during the deposition of S1. In addition, we provide a first quantitative estimate of the diffusive (1/e) mixing timescale for the eastern Mediterranean in its "stratified" sapropel mode, which is of the order of 450 years.

(Table S1) Diatom abundance according to ice cycles in ANDRILL AND1-1B drill core

The influence of Antarctica and the Southern Ocean on Late Pliocene global climate reconstructions has remained ambiguous due to a lack of well-dated Antarctic-proximal, paleoenvironmental records. Here we present ice sheet, sea-surface temperature, and sea ice reconstructions from the ANDRILL AND-1B sediment core recovered from beneath the Ross Ice Shelf. We provide evidence for a major expansion of an ice sheet in the Ross Sea that began at ~3.3 Ma, followed by a coastal sea surface temperature cooling of ~2.5°C, a stepwise expansion of sea ice, and polynya-style deep mixing in the Ross Sea between 3.3 and 2.5 Ma. The intensification of Antarctic cooling resulted in strengthened westerly winds and invigorated ocean circulation. The associated northward migration of Southern Ocean fronts has been linked with reduced Atlantic Meridional Overturning Circulation by restricting surface water connectivity between the ocean basins, with implications for heat transport to the high latitudes of the North Atlantic. While our results do not exclude low-latitude mechanisms as drivers for Pliocene cooling, they indicate an additional role played by southern high-latitude cooling during development of the bipolar world.

Age determination of sediment core W8709A-13

Northeast Pacific benthic foraminiferal d18O and d13 reveal repeated millennial-scale events of strong deep-sea ventilation (associated with nutrient depletion and/or high gas exchange) during stadial (cool, high ice volume) episodes from 10 to 60 ka, opposite the pattern in the deep North Atlantic. Two climate mechanisms may explain this pattern. North Pacific surface waters, chilled by atmospheric transmission from a cold North Atlantic and made saltier by reduced freshwater vapor transports, could have ventilated the deep Pacific from above. Alternatively, faster turnover of Pacific bottom and mid-depth waters, driven by Southern Ocean winds, may have compensated for suppressed North Atlantic Deep Water production during stadial intervals. During the Younger Dryas event (~11.6-13.0 cal ka), ventilation of the deep NE Pacific (~2700 m) lagged that in the Santa Barbara Basin (~450 m) by >500 years, suggesting that the NE Pacific was first ventilated at intermediate depth from above and then at greater depth from below. This apparent lag may reflect the adjustment time of global thermohaline circulation.