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).

Glacial North Atlantic: Sea surface conditions reconstructed by GLAMAP 2000

The response of the tropical ocean to global climate change and the extent of sea ice in the glacial nordic seas belong to the great controversies in paleoclimatology. Our new reconstruction of peak glacial sea surface temperatures (SSTs) in the Atlantic is based on census counts of planktic foraminifera, using the Maximum Similarity Technique Version 28 (SIMMAX-28) modern analog technique with 947 modern analog samples and 119 well-dated sediment cores. Our study compares two slightly different scenarios of the Last Glacial Maximum (LGM), the Environmental Processes of the Ice Age: Land, Oceans, Glaciers (EPILOG), and Glacial Atlantic Ocean Mapping (GLAMAP 2000) time slices. The comparison shows that the maximum LGM cooling in the Southern Hemisphere slightly preceeded that in the north. In both time slices sea ice was restricted to the north western margin of the nordic seas during glacial northern summer, while the central and eastern parts were ice-free. During northern glacial winter, sea ice advanced to the south of Iceland and Faeroe. In the central northern North Atlantic an anticyclonic gyre formed between 45° and 60°N, with a cool water mass centered west of Ireland, where glacial cooling reached a maximum of >12°C. In the subtropical ocean gyres the new reconstruction supports the glacial-to-interglacial stability of SST as shown by CLIMAP Project Members (CLIMAP) [1981]. The zonal belt of minimum SST seasonality between 2° and 6°N suggests that the LGM caloric equator occupied the same latitude as today. In contrast to the CLIMAP reconstruction, the glacial cooling of the tropical east Atlantic upwelling belt reached up to 6°-8°C during Northern Hemisphere summer. Differences between these SIMMAX-based and published U37[k]- and Mg/Ca-based equatorial SST records are ascribed to strong SST seasonalities and SST signals that were produced by different planktic species groups during different seasons.

Benthic (C. wuellerstorfi) d18O, d13C and Cd/Ca in sediment core NEAP 4K across the last deglaciation

We present evidence that the characteristic chemical signature (based on coupled benthic foraminiferal Cd/Ca and d13C) of Antarctic Intermediate waters (AAIW) penetrated throughout the intermediate depths of the Atlantic basin to the high-latitude North Atlantic during the abrupt cooling events of the last deglaciation: Heinrich 1 and the Younger Dryas. AAIW may play the dynamic counterpart to the "bipolar seesaw" when near-freezing salty bottom waters from the Antarctic (AABW) sluggishly ventilate the deep ocean. Our data reinforce the concept that interglacial circulation is stabilized by salinity feedbacks between salty northern sourced deep waters (NADW) and fresh southern sourced waters (AABW and AAIW). Further, the glacial ocean may be susceptible to the more finely balanced relative densities of NADW and AAIW, due to either freshwater input or a reversal of the salinity gradient, such that the ocean is poised for NADW collapse via a negative salinity feedback. The unstable climate of the glacial period and its termination may arise from the closer competition for ubiquity at intermediate depths between northern and southern sourced intermediate waters.

Boron isotope and B/Ca ratios of planktonic foraminifera

Here a new analytical methodology is described for measuring the isotopic composition of boron in foraminifera using multicollector inductively coupled plasma mass spectrometry (MC-ICPMS). This new approach is fast (~10 samples analysed in duplicate per analytical session) and accurate (to better than 0.25 per mil at 95% confidence) with acceptable sample size requirements (1-3 mg of carbonate). A core top calibration of several common planktic and two benthic species from geographically widespread localities shows a very close agreement between the isotopic composition measured by MC-ICPMS and the isotopic composition of B(OH)-4 in seawater (as predicted using the recently measured isotopic equilibrium factor of 1.0272) at the depth of habitat. A down core and core top investigation of boron concentration (B/Ca ratio) shows that the partition coefficient is influenced by [CO2-3] complicating the application of this proxy. Nevertheless, it is demonstrated that these two proxies can be used to fully constrain the carbonate system of surface water in the Caribbean Sea (ODP Site 999A) over the last 130 kyr. This reconstruction shows that during much of the Holocene and the last interglacial period surface water at Site 999A was in equilibrium with the atmosphere with respect to CO2. During the intervening colder periods although the surface water pCO2 was lower than the Holocene, it was a minor to significant source of CO2 to the atmosphere possibly due to either an expansion of the eastern equatorial Atlantic upwelling zone, or a more local expansion of coastal upwelling in the southern Caribbean. Such reorganisation of the oceanic carbonate system in favour of a larger source of CO2 to the atmosphere from the equatorial ocean may require mechanisms responsible for lowering atmospheric CO2 during glacial periods to be more efficient than previously supposed.

Great Barrier Reef coral element/Ca and stable isotope data and U-Th ages from IODP Expedition 325

Tropical south-western Pacific temperatures are of vital importance to the Great Barrier Reef (GBR), but the role of sea surface temperatures (SSTs) in the growth of the GBR since the Last Glacial Maximum remains largely unknown. Here we present records of Sr/Ca and d18O for Last Glacial Maximum and deglacial corals that show a considerably steeper meridional SST gradient than the present day in the central GBR. We find a 1-2 °C larger temperature decrease between 17° and 20°S about 20,000 to 13,000 years ago. The result is best explained by the northward expansion of cooler subtropical waters due to a weakening of the South Pacific gyre and East Australian Current. Our findings indicate that the GBR experienced substantial meridional temperature change during the last deglaciation, and serve to explain anomalous deglacial drying of northeastern Australia. Overall, the GBR developed through significant SST change and may be more resilient than previously thought.

Dinoflagellate cyst assemblages, oxygen isotope records of foraminifera, and alkenone and Mg/Ca-based SST estimates for the time interval 3.40-3.18 Ma (Late Pliocene) from five North Atlantic and Caribbean sediment cores

The early Late Pliocene (3.6 to ~3.0 million years ago) is the last extended interval in Earth's history when atmospheric CO2 concentrations were comparable to today's and global climate was warmer. Yet a severe global glaciation during marine isotope stage (MIS) M2 interrupted this phase of global warmth ~3.30 million years ago, and is seen as a premature attempt of the climate system to establish an ice-age world. Our geochemical and palynological records from five marine sediment cores along a Caribbean to eastern North Atlantic transect show that increased Pacific-to-Atlantic flow via the Central American Seaway weakened the North Atlantic Current (NAC) and attendant northward heat transport prior to MIS M2. The consequent cooling of the northern high latitude oceans permitted expansion of the Greenland ice sheet during MIS M2, despite near-modern atmospheric CO2 concentrations. Before and after MIS M2, heat transport via the NAC was crucial in maintaining warm climates comparable to those predicted for the end of this century.

Stable oxygen isotope record and Sr/Ca ratios of corals from IODP Holes M0009D and M0023A

We present monthly resolved records of strontium/calcium (Sr/Ca) and oxygen isotope (d18O) ratios from well-preserved fossil corals drilled during the Integrated Ocean Drilling Program (IODP) Expedition 310 'Tahiti Sea Level' and reconstruct sea surface conditions in the central tropical South Pacific Ocean during two time windows of the last deglaciation. The two Tahiti corals examined here are uranium/thorium (U/Th)-dated at 12.4 and 14.2 ka, which correspond to the Younger Dryas (YD) cold reversal and the Bølling-Allerød (B-A) warming of the Northern Hemisphere, respectively. The coral Sr/Ca records indicate that annual average sea surface temperature (SST) was 2.6-3.1 °C lower at 12.4 ka and 1.0-1.6 °C lower at 14.2 ka relative to the present, with no significant changes in the amplitude of the seasonal SST cycle. These cooler conditions were accompanied by seawater d18O (d18Osw) values higher by ~0.8 per mill and ~0.6 per mill relative to the present at 12.4 and 14.2 ka, respectively, implying more saline conditions in the surface waters. Along with previously published coral Sr/Ca records from the island [Cohen and Hart (2004), Deglacial sea surface temperatures of the western tropical Pacific: A new look at old coral. Paleoceanography 19, PA4031, doi:10.1029/2004PA001084], our new Tahiti coral records suggest that a shift toward lower SST by ~1.5 °C occurred from 13.1 to 12.4 ka, which was probably associated with a shift toward higher d18Osw by ~0.2 per mill. Along with a previously published coral Sr/Ca record from Vanuatu [Corrège et al. (2004), Interdecadal variation in the extent of South Pacific tropical waters during the Younger Dyras event. Nature 428, 927-929], the Tahiti coral records provide new evidence for a pronounced cooling of the western to central tropical South Pacific during the Northern Hemisphere YD event.

Ground truthing the Cd/Ca-carbon isotope relationship in foraminifera of the Greenland-Iceland-Norwegian Seas

In order to examine whether the paleoceanographic nutrient proxies, d13C and cadmium/calcium in foraminiferal calcite, are well coupled to nutrients in the region of North Atlantic Deep Water formation, we present da ta from two transects of the Greenland-Iceland-Norwegian Seas. Along Transect A (74.3°N, 18.3°E to 75.0°N, 12.5°W, 15 stations), we measured phosphate and Cd concentrations of modern surface sea water. Along Transect B (64.5°N, 0.7°W to 70.4°N, 18.2°W, 14 stations) we measured Cd/Ca ratios and d13C of the planktonic foraminifera Neogloboquadrina pachyderma sinistral in core top sediments. Our results indicate that Cd and phosphate both vary with surface water mass and are well correlated along Transect A. Our planktonic foraminiferal d13C data indicate similar nutrient variation with water mass along Transect B. Our Cd/Ca data hint at the same type of nutrient variability, but interpretations are hampered by low values close to the detection limit of this technique and therefore relatively large error bars. We also measured Cd and phosphate concentrations in water depth profiles at three sites along Transect A and the d13C of the benthic foraminifera Cibicidoides wuellerstorfi along Transect B. Modern sea water depth profiles along Transect A have nutrient depletions at the surface and then constant values at depths greater than 100 meters. The d13C of planktonic and benthic foraminifera from Transect B plotted versus depth also reflect surface nutrient depletion and deep nutrient enrichment as seen at Transect A, with a small difference between intermediate and deep waters. Overall we see no evidence for decoupling of Cd/Ca ratio and d13C in foraminiferal calcite from water column nutrient concentrations along these transects in a region of North Atlantic Deep Water formation.

Sr/Ca ratios of coccoliths from the Paleocene-Eocene Thermal Maximum

During the Paleocene-Eocene Thermal Maximum (PETM), rapid release of isotopically light C to the ocean-atmosphere system elevated the greenhouse effect and warmed temperatures by 5-7 °C for 105 yr. The response of the planktic ecosystems and productivity to the dramatic climate changes of the PETM may represent a significant feedback to the carbon cycle changes, but has been difficult to document. We examine Sr/Ca ratios in calcareous nannofossils in sediments spanning the PETM in three open ocean sites as a new approach to examine productivity and ecological shifts in calcifying plankton. The large heterogeneity in Sr/Ca among different nannofossil genera indicates that nannofossil Sr/Ca reflects primary productivity-driven geochemical signals and not diagenetic overprinting. Elevated Sr/Ca ratios in several genera and constant ratios in other genera suggest increased overall productivity in the Atlantic sector of the Southern Ocean during the PETM. Dominant nannofossil genera in tropical Atlantic and Pacific sites show Sr/Ca variations during the PETM which are comparable to background variability prior to the PETM. Despite acidification of the ocean there was not a productivity crisis among calcifying phytoplankton. We use the Pandora ocean box model to explore possible mechanisms for PETM productivity change. If independent proxy evidence for more stratified conditions in the Southern Ocean during the PETM is robust, then maintenance of stable or increased productivity there likely reflects increased nutrient inventories of the ocean. Increased nutrient inventories could have resulted from climatically enhanced weathering and would have important implications for burial rates of organic carbon and stabilization of climate and the carbon cycle.

(Table 1) Calcium isotope composition, Mg/Ca, Sr/Ca and U/Ca ratios of calcareous cysts of Thoracosphaera heimii

We investigated the influences of temperature, salinity and pH on the calcium isotope as well as trace and minor element (uranium, strontium, magnesium) to Ca ratios on calcium carbonate cysts of the calcareous dinoflagellate species Thoracosphaera heimii grown in laboratory cultures. The natural habitat of this species is the photic zone (preferentially at the chlorophyll maximum depth) of temperate to tropical oceans, and it is abundant in deep-sea sediments over the entire Cenozoic. In our experiments, temperatures ranged from 12 to 30 °C, salinity from 36.5 to 38.8 and pH from 7.9 to 8.4. The delta44/40Ca of T. heimii cysts resembles that of other marine calcifiers, including coccolithophores, foraminifers and corals. However, its temperature sensitivity is considerably smaller and statistically insignificant, and T. heimii might serve as a recorder of changes in seawater delta44/40Ca over geologic time. The Sr/Ca ratios of T. heimii cysts show a pronounced temperature sensitivity (0.016 mmol/mol °C**-1) and have the potential to serve as a palaeo-sea surface temperature proxy. No clear temperature- and pH-dependences were observed for Mg/Ca. U/Ca seems to be influenced by temperature and pH, but the correlations change sign at 23 °C and pH 8.2, respectively.

Stable oxygen isotope record and Mg/Ca ratios at the Mid-Pleistocene Climate Transition, ODP Site 181-1123

Earth's climate underwent a fundamental change between 1250 and 700 thousand years ago, the Mid-Pleistocene Transition (MPT), when the dominant periodicity of climate cycles changed from 41,000 to 100,000 years in the absence of significant change in orbital forcing. Over this time, an increase occurred in the amplitude of change of deep ocean foraminiferal oxygen isotopic ratios, traditionally interpreted as defining the main rhythm of ice ages although containing large effects of changes in deep-ocean temperature. We have separated the effects of decreasing temperature and increasing global ice volume on oxygen isotope ratios. Our results suggest that the MPT was initiated by an abrupt increase in Antarctic ice volume at 900 ka. We see no evidence of a pattern of gradual cooling but near-freezing temperatures occur at every glacial maximum.