(Table S3) Benthic foraminiferal Mg/Ca, Li/Ca and stable isotope ratios of ODP Site 122-761

The sensitivities of benthic foraminiferal Mg/Ca and Li/Ca to bottom water temperature and carbonate saturation state have recently been assessed. Here we present a new approach that uses paired Mg/Ca and Li/Ca records to calculate simultaneous changes in temperature and saturation state. Using previously published records, we first use this approach to document a cooling of deep ocean waters associated with the establishment of the Antarctic ice sheet at the Eocene-Oligocene climate transition. We then apply this approach to new records of the Middle Miocene Climate Transition from ODP Site 761 to estimate variations in bottom water temperature and the oxygen isotopic composition of seawater. We estimate that the oxygen isotopic composition of seawater varied by ~1 per mil between the deglacial extreme of the Miocene Climatic Optimum and the glacial maximum following the Middle Miocene Climate Transition, indicating large amplitude variations in ice volume. However, the longer-term change between 15.3 and 12.5 Ma is marked by a ~1°C cooling of deep waters, and an increase in the oxygen isotopic composition of seawater of ~0.6 per mil. We find that bottom water saturation state increased in the lead up to the Middle Miocene Climate Transition and decreased shortly after. This supports decreasing pCO2 as a driver for global cooling and ice sheet expansion, in agreement with existing boron isotope and leaf stomatal index CO2 records but in contrast to the published alkenone CO2 records.

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

Dinoflagellate cyst assemblages, Mg/Ca based reconstructed temperatures, stable carbon isotopes, and geochemical parameters of sediment core GeoB9508-5

In order to investigate a possible connection between tropical northeast (NE) Atlantic primary productivity, Atlantic meridional overturning circulation (AMOC), and drought in the Sahel region during Heinrich Stadial 1 (HS1), we used dinoflagellate cyst (dinocyst) assemblages, Mg/Ca based reconstructed temperatures, stable carbon isotopes (d13C) and geochemical parameters of a marine sediment core (GeoB 9508-5) from the continental slope offshore Senegal. Our results show a two-phase productivity pattern within HS1 that progressed from an interval of low marine productivity between ~ 19 and 16 kyr BP to a phase with an abrupt and large productivity increase from ~ 16 to 15 kyr BP. The second phase is characterized by distinct heavy planktonic d13C values and high concentrations of heterotrophic dinocysts in addition to a significant cooling signal based on reconstructions of past sea surface temperatures (SST). We conclude that productivity variations within HS1 can be attributed to a substantial shift of West African atmospheric processes. Taken together our results indicate a significant intensification of the North East (NE) trade winds over West Africa leading to more intense upwelling during the last millennium of HS1 between ~ 16 and 15 kyr BP, thus leaving a strong imprint on the dinocyst assemblages and sea surface conditions. Therefore, the two-phase productivity pattern indicates a complex hydrographic setting suggesting that HS1 cannot be regarded as uniform as previously thought.

Mg/Ca-temperature proxy in benthic foraminifera from the Florida Strait

Over the past decade, the ratio of Mg to Ca in foraminiferal tests has emerged as a valuable paleotemperature proxy. However, large uncertainties remain in the relationships between benthic foraminiferal Mg/Ca and temperature. Mg/Ca was measured in benthic foraminifera from 31 high-quality multicore tops collected in the Florida Straits, spanning a temperature range of 5.8° to 18.6°C. New calibrations are presented for Uvigerina peregrina, Planulina ariminensis, Planulina foveolata, and Hoeglundina elegans. The Mg/Ca values and temperature sensitivities vary among species, but all species exhibit a positive correlation that decreases in slope at higher temperatures. The decrease in the sensitivity of Mg/Ca to temperature may potentially be explained by Mg/Ca suppression at high carbonate ion concentrations. It is suggested that a carbonate ion influence on Mg/Ca may be adjusted for by dividing Mg/Ca by Li/Ca. The Mg/Li ratio displays stronger correlations to temperature, with up to 90% of variance explained, than Mg/Ca alone. These new calibrations are tested on several Last Glacial Maximum (LGM) samples from the Florida Straits. LGM temperatures reconstructed from Mg/Ca and Mg/Li are generally more scattered than core top measurements and may be contaminated by high-Mg overgrowths. The potential for Mg/Ca and Mg/Li as temperature proxies warrants further testing.

(Table S1) Original element/Ca data for the OGA and Tahiti corals

[1] Previous studies have demonstrated the potential for the Li content of coral aragonite to record information about environmental conditions, but no detailed study of tropical corals exists. Here we present the Li and Mg to Ca ratios at a bimonthly to monthly resolution over 25 years in two modern Porites corals, the genus most often used for paleoclimate reconstructions in the tropical Indo-Pacific. A strong relationship exists between coral Li/Ca and locally measured SST, indicating that coral Li/Ca can be used to reconstruct tropical SST variations. However, Li/Ca ratios of the skeleton deposited during 1979-1980 do not track local SST well and are anomalously high in places. The Mg/Ca ratios of this interval are also anomalously high, and we suggest Li/Ca can be used to reconstruct tropical SST only when Mg/Ca data are used to carefully screen for relatively rare biological effects. Mg/Li or Li/Mg ratios provide little advantage over Li/Ca ratios, except that the slope of the Li/Mg temperature relationship is more similar between the two corals. The Mg/Li temperature relationship for the coral that experienced a large temperature range is similar to that found for cold water corals and aragonitic benthic foraminifera in previous studies. The comparison with data from other biogenic aragonites suggests the relationship between Li/Mg and water temperature can be described by a single exponential relationship. Despite this hint at an overarching control, it is clear that biological processes strongly influence coral Li/Ca, and more calibration work is required before widely applying the proxy.

(Appendix) Stable carbon and oxygen isotopes, and Mg/Ca ratios of foraminifera from ODP Hole 154-925A

To better understand the links between the carbon cycle and changes in past climate over tectonic timescales we need new geochemical proxy records of secular change in silicate weathering rates. A number of proxies are under development, but some of the most promising (e.g. palaeoseawater records of Li and Nd isotope change) can only be employed on such large samples of mono-specific foraminifera that application to the deep sea sediment core archive becomes highly problematic. "Dentoglobigerina" venezuelana presents a potentially attractive target for circumventing this problem because it is a typically large (> 355 ?m diameter), abundant and cosmopolitan planktic foraminifer that ranges from the early Oligocene to early Pliocene. Yet considerable taxonomic and ecological uncertainties associated with this taxon must first be addressed. Here, we assess the taxonomy, palaeoecology, and ontogeny of "D." venezuelana using stable isotope (oxygen and carbon) and Mg/Ca data measured in tests of late Oligocene to early Miocene age from Ocean Drilling Program (ODP) Site 925, on Ceara Rise, in the western equatorial Atlantic. To help constrain the depth habitat of "D." venezuelana relative to other species we report the stable isotope composition of selected planktic foraminifera species within Globigerina, Globigerinoides, Paragloborotalia and Catapsydrax. We define three morphotypes of "D." venezuelana based on the morphology of the final chamber and aperture architecture. We determine the trace element and stable isotope composition of each morphotype for different size fractions, to test the validity of pooling these morphotypes for the purposes of generating geochemical proxy datasets and to assess any ontogenetic variations in depth habitat. Our data indicate that "D." venezuelana maintains a lower thermocline depth habitat at Ceara Rise between 24 and 21 Ma. Comparing our results to published datasets we conclude that this lower thermocline depth ecology for the Oligo-Miocene is part of an Eocene-to-Pliocene evolution of depth habitat from surface to sub-thermocline for "D." venezuelana. Our size fraction data advocate the absence of photosymbionts in "D." venezuelana and suggest that juveniles calcify higher in the water column, descending into slightly deeper water during the later stages of its life cycle. Our morphotype data show that d18O and d13C variation between morphotypes is no greater than within-morphotype variability. This finding will permit future pooling of morphotypes in the generation of the "sample hungry" palaeoceanographic records.