Plio-Pleistocene Mg/Ca and TEX86 sea surface temperature records of ODP Sites 184-1143 and 165-999,

The western warm pools of the Atlantic and Pacific oceans are a critical source of heat and moisture for the tropical climate system. Over the past five million years, global mean temperatures have cooled by 3-4 °C. Yet, current reconstructions of sea surface temperatures indicate that temperature in the warm pools has remained stable during this time. This stability has been used to suggest that tropical sea-surface temperatures are controlled by some sort of thermostat-like regulation. Here we reconstruct sea surface temperatures in the South China Sea, Caribbean Sea and western equatorial Pacific Ocean for the past five million years, using a combination of the Mg/Ca, TEXH86-and Uk'37 surface temperature proxies. Our data indicate that during the period of Pliocene warmth from about 5 to 2.6 million years ago, the western Pacific and western Atlantic warm pools were about 2 °C warmer than today. We suggest that the apparent lack of warming seen in the previous reconstructions was an artefact of low seawater Mg/Ca ratios in the Pliocene oceans. Taking this bias into account, our data indicate that tropical sea surface temperatures did change in conjunction with global mean temperatures. We therefore conclude that the temperature of the warm pools of the equatorial oceans during the Pliocene was not limited by a thermostat-like mechanism.

Benthic foraminiferal d18O and Mg/Ca studies of sediment cores off tropical NW Africa

Benthic foraminiferal d18O and Mg/Ca of sediment cores off tropical NW Africa are used to study the properties of Atlantic central waters during the Last Glacial Maximum (LGM) and Heinrich Stadial 1 (HS1). We combined our core top data with published results to develop a new Mg/Ca-temperature calibration for Planulina ariminensis, which shows a Mg/Ca-temperature sensitivity of 0.19 mmol/mol per °C. Estimates of the LGM and HS1 thermocline temperatures are comparable to the present-day values between 200 and 400 m water depth, but were 1.2-1.5°C warmer at 550-570 m depth. The HS1 thermocline waters (200-570 m depth) did not show any warming relative to the LGM. This is in contrast to previous climate model studies, which concluded that tropical Atlantic thermocline waters warmed significantly when Atlantic meridional overturning circulation was reduced. However, our results suggest that thermocline temperatures of the northeastern tropical Atlantic show no pronounced sensitivity to changes in the thermohaline circulation during glacial periods. In contrast, we find a significant increase in thermocline-water salinity during the LGM (200-550 m depth) and HS1 (200-400 m depth) with respect to the present-day, which we relate to changes in the wind-driven circulation. We infer that the LGM thermocline (200-550 m depth) and the HS1 upper thermocline (200-400 m depth) in the northeastern tropical Atlantic was ventilated by surface waters from the North Atlantic rather than the southern-sourced waters. This suggests that the frontal zone between the modern South Atlantic and North Atlantic Central Waters was probably shifted southward during the LGM and HS1.

B/Ca ratios of benthic foraminifera

We present modern B/Ca core-top calibrations for the epifaunal benthic foraminifer Nuttallides umbonifera and the infaunal Oridorsalis umbonatus to test whether B/Ca values in these species can be used for the reconstruction of paleo-D[[CO3]2-]. O. umbonatus originated in the Late Cretaceous and remains extant, whereas N. umbonifera originated in the Eocene and is the closest extant relative to Nuttallides truempyi, which ranges from the Late Cretaceous through the Eocene. We measured B/Ca in both species in 35 Holocene sediment samples from the Atlantic, Pacific and Southern Oceans. B/Ca values in epifaunal N. umbonifera (~ 85-175 µmol/mol) are consistently lower than values reported for epifaunal Cibicidoides (Cibicides) wuellerstorfi (130-250 µmol/mol), though the sensitivity of D[[CO3]2-] on B/Ca in N. umbonifera (1.23 ± 0.15) is similar to that in C. wuellerstorfi (1.14 ± 0.048). In addition, we show that B/Ca values of paired N. umbonifera and its extinct ancestor, N. truempyi, from Eocene cores are indistinguishable within error. In contrast, both the B/Ca (35-85 µmol/mol) and sensitivity to D[[CO3]2-] (0.29 ± 0.20) of core-top O. umbonatus are considerably lower (as in other infaunal species), and this offset extends into the Paleocene. Thus the B/Ca of N. umbonifera and its ancestor can be used to reconstruct bottom water D[[CO3]2?], whereas O. umbonatus B/Ca appears to be buffered by porewater [[CO3]2-] and suited for constraining long-term drift in seawater B/Ca.

Grain size analyses, Ca/Ti and other selected elemental ratios from XRF core scanning of sediment cores from the Apulian Margin

An extensive, high-resolution, sedimentological-geochemical survey was done using geo-acoustics, XRF-core scans, ICP-AES, AMS 14C-dating and grain size analyses of sediments in 11 cores from the Gulf of Taranto, the southern Adriatic Sea, and the central Ionian Sea spanning the last 16 cal. ka BP. Comparable results were obtained for cores from the Gallipoli Shelf (eastern Gulf of Taranto), and the southern Adriatic Sea suggesting that the dominant provenance of Gallipoli Shelf sediments is from the western Adriatic mud belt. The 210Pb and 14C-dated high-accumulation-rate sediments permit a detailed reconstruction of climate variability over the last 16 cal. ka BP. Although, the Glacial-Interglacial transition is generally dry and stable these conditions are interrupted by two phases of increased detrital input during the Bølling-Allerød and the late Younger Dryas. The event during the Younger Dryas period is characterized by increased sediment inputs from southern Italian sources. This suggests that run-off was higher in southern- compared to northern Italy. At approximately ~ 7 cal. ka BP, increased detrital input from the Adriatic mud belt, related to sea level rise and the onset of deep water formation in the Adriatic Sea, is observed and is coincident with the end of sapropel S1 formation in the southern Adriatic Sea. During the mid-to-late Holocene we observed millennial-scale events of increased detrital input, e.g. during the Roman Humid Period, and of decreased detrital input, e.g., Medieval Warm Period. These dry/wet spells are consistent with variability in the North Atlantic Oscillation (NAO). A negative state of the NAO and thus a more advanced penetration of the westerlies into the central Mediterranean, that result in wet conditions in the research area concord with events of high detrital input e.g., during the Roman Humid Period. In contrast, a positive state of the NAO, resulting in dry conditions in the Mediterranean, dominated during events of rapid climate change such as the Medieval Warm Period and the Bronze Age.

Oxygen isotopes and Sr/Ca ratios of fossil corals from Bonaire

Explanations for the demise of the Classic Maya civilization on the Yucatán Peninsula during the Terminal Classic Period (TCP; CE 750-1050) are controversial. Multiyear droughts are one likely cause, but the role of the Caribbean Sea, the dominant moisture source for Mesoamerica, remains largely unknown. Here we present bimonthly resolved snapshots of reconstructed sea surface temperature (SST) and salinity (SSS) variability in the southern Caribbean from precisely dated fossil corals. Our fossil coral results from Bonaire indicate strong interannual to decadal SST and SSS variability in the southern Caribbean Sea during the TCP with multiyear extremes of high SSS and high SST that coincide with droughts on the Yucatán Peninsula. The results are best explained by changed Caribbean SST gradients affecting the Caribbean low-level atmospheric jet with consequences for Mesoamerican precipitation, possibly linked to changes in Atlantic Meridional Overturning Circulation strength. Our findings provide a new perspective on the anomalous hydrological changes on the Yucatán Peninsula during the TCP that complement the often-suggested southward displacement of the Intertropical Convergence Zone. We advocate for a strong role of ocean-atmosphere interactions in the Caribbean Sea related to the multiyear variations in Caribbean Sea surface conditions as an important driver of the spatially complex pattern of hydrological anomalies during the TCP.

Sable oxygen isotope and Mg/Ca ratios of planktonic foraminifera from DSDP Hole 31-292

This study presents new evidence of when and how the Western Pacific Warm Pool (WPWP) was established in its present form. We analyzed planktic foraminifera, oxygen isotopes, and Mg/Ca ratios in upper Miocene through Pleistocene sediments collected at Deep Sea Drilling Program (DSDP) Site 292. These data were then compared with those reported from Ocean Drilling Program (ODP) Site 806. Both drilling sites are located in the western Pacific Ocean. DSDP Site 292 is located in the northern margin of the modern WPWP and ODP Site 806 near the center of the WPWP. Three stages of development in surface-water conditions are identified in the region using planktic foraminferal data. During the initial stage, from 8.5 to 4.4 Ma, Site 806 was overlain by warm surface water but Site 292 was not, as indicated by the differences in faunal compositions and sea-surface temperature (SST) between the two sites. In addition, the vertical thermal gradient at Site 292 was weak during this period, as indicated by the small differences in the delta18O values between Globigerinoides sacculifer and Pulleniatina spp. During stage two, from 4.4 to 3.6 Ma, the SST at Site 292 rapidly increased to 27 °C, but the vertical thermal gradient had not yet be strengthened, as shown by Mg/Ca ratios and the presence of both mixed-layer dwellers and thermocline dwellers. Finally, a warm mixed layer with a high SST ca. 28 °C and a strong vertical thermal gradient were established at Site 292 by 3.6 Ma. This event is marked by the dominance of mixed-layer dwellers, a high and stable SST, and a larger differences in the delta18O values between G. sacculifer and Pulleniatina spp. Thus, evidence of surface-water evolution in the western Pacific suggests that Site 292 came under the influence of the WPWP at 3.6 Ma. The northward expansion of the WPWP from 4.4 to 3.6 Ma and the establishment of the modern WPWP by 3.6 Ma appear to be closely related to the closure of the Indonesian and Central American seaways.

Mg/Ca ratios of Cenozoic benthic foraminifera from the Indian and South Atlantic Ocean

We investigate the evolution of Cenozoic climate and ice volume as evidenced by the oxygen isotopic composition of seawater (delta18Osw) derived from benthic foraminiferal Mg/Ca ratios to constrain the temperature effect contained in foraminiferal delta18O values. We have constructed two benthic foraminiferal Mg/Ca records from intermediate water depth sites (Ocean Drilling Program sites 757 and 689 from the subtropical Indian Ocean and the Weddell Sea, respectively). Together with the previously published composite record of Lear et al. (2002, doi:10.1126/science.287.5451.269) and the Neogene record from the Southern Ocean of Billups and Schrag (2002, doi:10.1029/2000PA000567), we obtain three, almost complete representations of the delta18Osw for the past 52 Myr. We discuss the sensitivity of early Cenozoic Mg/Ca-derived paleotemperatures (and hence the delta18Osw) to assumptions about seawater Mg/Ca ratios. We find that during the middle Eocene (~ 49-40 Ma), modern seawater ratios yield Mg/Ca-derived temperatures that are in good agreement with the oxygen isotope paleothermometer assuming ice-free conditions. Intermediate waters cooled during the middle Eocene reaching minimum temperatures by 40 Ma. The corresponding delta18Osw reconstructions support ice growth on Antarctica beginning by at least 40 Ma. At the Eocene/Oligocene boundary, Mg/Ca ratios (and hence temperatures) from Weddell Sea site 689 display a well-defined maximum. We caution against a paleoclimatic significance of this result and put forth that the partitioning coefficient of Mg in benthic foraminifera may be sensitive to factors other than temperature. Throughout the remainder of the Cenozoic, the temporal variability among delta18Osw records is similar and similar to longer-term trends in the benthic foraminiferal delta18O record. An exception occurs during the Pliocene when delta18Osw minima in two of the three records suggest reductions in global ice volume that are not apparent in foraminiferal delta18O records, which provides a new perspective to the ongoing debate about the stability of the Antarctic ice sheet. Maximum delta18Osw values recorded during the Pleistocene at Southern Ocean site 747 agree well with values derived from the geochemistry of pore waters (Schrag et al., 1996, doi:10.1126/science.272.5270.1930) further highlighting the value of the new Mg/Ca calibrations of Martin et al. (2002, doi:10.1016/S0012-821X(02)00472-7) and Lear et al. (2002, doi:10.1016/S0016-7037(02)00941-9) applied in this study. We conclude that the application of foraminiferal Mg/Ca ratios allows a refined view of Cenozoic ice volume history despite uncertainties related to the geochemical cycling of Mg and Ca on long time scales.

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

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.

Mg/Ca ratios in planktonic foraminifera

The primary Mg/Ca ratio of foraminiferal shells is a potentially valuable paleoproxy for sea surface temperature (SST) reconstructions. However, the reliable extraction of this ratio from sedimentary calcite assumes that we can overcome artifacts related to foraminiferal ecology and partial dissolution, as well as contamination by secondary calcite and clay. The standard batch method for Mg/Ca analysis involves cracking, sonicating, and rinsing the tests to remove clay, followed by chemical cleaning, and finally acid-digestion and single-point measurement. This laborious procedure often results in substantial loss of sample (typically 30-60%). We find that even the earliest steps of this procedure can fractionate Mg from Ca, thus biasing the result toward a more variable and often anomalously low Mg/Ca ratio. Moreover, the more rigorous the cleaning, the more calcite is lost, and the more likely it becomes that any residual clay that has not been removed by physical cleaning will increase the ratio. These potentially significant sources of error can be overcome with a flow-through (FT) sequential leaching method that makes time- and labor-intensive pretreatments unnecessary. When combined with time-resolved analysis (FT-TRA) flow-through, performed with a gradually increasing and highly regulated acid strength, produces continuous records of Mg, Sr, Al, and Ca concentrations in the leachate sorted by dissolution susceptibility of the reacting material. Flow-through separates secondary calcite from less susceptible biogenic calcite and clay, and further resolves the biogenic component into primary and more resistant fractions. FT-TRA reliably separates secondary calcite (which is not representative of original life habitats) from the more resistant biogenic calcite (the desired signal) and clay (a contaminant of high Mg/Ca, which also contains Al), and further resolves the biogenic component into primary and more resistant fractions that may reflect habitat or other changes during ontogeny. We find that the most susceptible fraction of biogenic calcite in surface dwelling foraminifera gives the most accurate value for SST and therefore best represents primary calcite. Sequential dissolution curves can be used to correct the primary Mg/Ca ratio for clay, if necessary. However, the temporal separation of calcite from clay in FT-TRA is so complete that this correction is typically <=2%, even in clay-rich sediments. Unlike hands-on batch methods, that are difficult to reproduce exactly, flow-through lends itself to automation, providing precise replication of treatment for every sample. Our automated flow-through system can process 22 samples, two system blanks, and 48 mixed standards in <12 hours of unattended operation. FT-TRA thus represents a faster, cheaper, and better way to determine Mg/Ca ratios in foraminiferal calcite.

Benthic foraminiferal Cd/Ca and isotope record from the South Atlantic

A depth transect of cores from 1268 to 3909 m water depth in the western South Atlantic are ideally situated to monitor the interocean exchange of deep water and variations in the relative strength of northern and southern sources of deep water production. Benthic foraminiferal Cd/Ca and d13C data suggest that Glacial North Atlantic Intermediate Water (GNAIW) extended at least as far south as 28°S in the western South Atlantic. The core of nutrient-depleted water was situated at ~1500 m, above and below water masses with higher nutrient concentrations. When examined in conjunction with published paired Cd/Ca and d13C from intermediate depth cores from other basins, it appears that the extent of GNAIW influence on the intermediate waters of the world's oceans was less than suggested previously. Differentiating among possible pathways for the glacial deep ocean (>3 km) requires a better understanding of the controls on Cd/Ca and d13C values of benthic foraminifera.

Mg/Ca and stable oxygen measurement of the past 27 Ma of ODP Hole 120-747A

We explore the applicability of paired Mg/Ca and 18O/16O measurements on benthic foraminifera from Southern Ocean site 747 to paleoceanographic reconstructions on pre-Pleistocene timescales. We focus on the late Oligocene through Pleistocene (27-0 Ma) history of paleotemperatures and the evolution of the d18O values of seawater (d18Osw) at a temporal resolution of ~100-200 kyr. Absolute paleotemperature estimates depend on assumptions of how Mg/Ca ratios of seawater have changed over the past 27 Myr, but relative changes that occur on geologically brief timescales are robust. Results indicate that at the Oligocene to Miocene boundary (23.8 Ma), temperatures lag the increase in global ice-volume deduced from benthic foraminiferal d18O values, but the smaller-scale Miocene glaciations are accompanied by ocean cooling of -1°C. During the mid-Miocene phase of Antarctic ice sheet growth (~15-13 Ma), water temperatures cool by ~3°C. Unlike the benthic foraminiferal d18O values, which remain relatively constant thereafter, temperatures vary (by 3°C) and reach maxima at ~12 and ~8.5 Ma. The onset of significant Northern Hemisphere glaciation during the late Pliocene is synchronous with an ~4°C cooling at site 747. A comparison of our d18Osw curve to the Haq et al. (1987, doi:10.1126/science.235.4793.1156 ) sea level curve yields excellent agreement between sequence boundaries and times of increasing seawater 18O/16O ratios. At ~12-11 Ma in particular, when benthic foraminiferal d18O values do not support a further increase in ice volume, the d18Osw curve comes to a maximum that corresponds to a major mid-Miocene sea level regression. The agreement between the character of our Mg/Ca-based d18Osw curve and sequence stratigraphy demonstrates that benthic foramaniferal Mg/Ca ratios can be used to trace the d18Osw on pre-Pleistocene timescales despite a number of uncertainties related to poorly constrained temperature calibrations and paleoseawater Mg/Ca ratios. The Mg/Ca record also highlights that deep ocean temperatures can vary independently and unexpectedly from ice volume changes, which can lead to misinterpretations of the d18O record.

Sr/Ca ratios of Cenozoic benthic foraminifera

A Cenozoic multi-species record of benthic foraminiferal calcite Sr/Ca has been produced and is corrected for interspecific offsets (typically less than 0.3 mmol/mol) and for the linear relationship between decreasing benthic foraminiferal Sr/Ca and increasing water depth. The water depth correction, determined from Holocene, Late Glacial Maximum and Eocene paleowater-depth transects, is ~0.1 mmol/mol/km. The corrected Cenozoic benthic foraminiferal Sr/Ca record ranges from 1.2 to 2.0 mmol/mol, and has been interpreted in terms of long-term changes in seawater Sr/Ca, enabling issues related to higher-resolution variability in Sr/Ca to be ignored. We estimate that seawater Sr/Ca was ~1.5 times modern values in the late Cretaceous, but declined rapidly into the Paleogene. Following a minimum in the Eocene, seawater Sr/Ca increased gradually through to the present day with a minimum superimposed on this trend centered in the late Miocene. By assuming scenarios for changing seawater calcium concentration, and using published carbonate accumulation rate data combined with suitable values for Sr partition coefficients into carbonates, the seawater Sr/Ca record is used to estimate global average river Sr fluxes. These fluxes are used in conjunction with the seawater strontium isotope curve and estimates of hydrothermal activity/tectonic outgassing to calculate changes in global average river 87Sr/86Sr through the Cenozoic. The absolute magnitude of Sr fluxes and isotopic compositions calculated in this way are subject to relatively large uncertainties. Nevertheless, our results suggest that river Sr flux increased from 35 Ma to the present day (roughly two-fold) accompanied by an overall increase in 87Sr/86Sr (by ~0 to 0.001). Between 75 and 35 Ma, river 87Sr/86Sr also increased (by ~0.001 to 0.002) but was accompanied by a decrease (two- to three-fold) in river Sr flux.