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.

Mg/Ca and stable isotope record of benthic foraminifera of the Paleocene-Eocene Thermal Maximum

A rapid increase in greenhouse gas levels is thought to have fueled global warming at the Paleocene-Eocene Thermal Maximum (PETM). Foraminiferal magnesium/calcium ratios indicate that bottom waters warmed by 4° to 5°C, similar to tropical and subtropical surface ocean waters, implying no amplification of warming in high-latitude regions of deep-water formation under ice-free conditions. Intermediate waters warmed before the carbon isotope excursion, in association with downwelling in the North Pacific and reduced Southern Ocean convection, supporting changing circulation as the trigger for methane hydrate release. A switch to deep convection in the North Pacific at the PETM onset could have amplified and sustained warming.

Faunal analysis and Mg/Ca records of planktic foraminifera in sediment core SO189/2_039KL

Climate phenomena like the monsoon system, El Niño Southern Oscillation (ENSO) and the Indian Ocean Dipole (IOD) are interconnected via various feedback mechanisms and control the climate of the Indian Ocean and its surrounding continents on various timescales. The eastern tropical Indian Ocean is a key area for the interplay of these phenomena and for reconstructing their past changes and forcing mechanisms. Here we present records of upper ocean thermal gradient, thermocline temperatures (TT) and relative abundances of planktic foraminifera in core SO 189-39KL taken off western Sumatra (0°47.400' S, 99°54.510' E) for the last 8 ka that we use as proxies for changes in upper ocean structure. The records suggest a deeper thermocline between 8 ka and ca 3 ka compared to the late Holocene. We find a shoaling of the thermocline after 3 ka, most likely indicating an increased occurrence of upwelling during the late Holocene compared to the mid-Holocene which might represent changes in the IOD-like mean state of the Indian Ocean with a more negative IOD-like mean state during the mid-Holocene and a more positive IOD-like mean state during the past 3 ka. This interpretation is supported by a transient Holocene climate model simulation in which an IOD-like mode is identified that involves an insolation-forced long-term trend of increasing anomalous surface easterlies over the equatorial eastern Indian Ocean.

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.

Mg/Ca, Sr/Ca and Ca isotope ratios in benthonic foraminifers of surface sediment samples

We analysed Mg/Ca, Sr/Ca and Ca isotope ratios of benthonic foraminifers from sediment core tops retrieved during several research cruises in the Atlantic Ocean, in order to improve the understanding of isotope fractionation and element partitioning resulting from biomineralisation processes and changes in ambient conditions. Species include foraminifers secreting tests composed of hyaline low magnesium calcite, porcelaneous high magnesium calcite as well as aragonite. Our results demonstrate systematic isotope fractionation and element partitioning patterns specific for these foraminiferal groups. Calcium isotope fractionation is similar in porcelaneous and hyaline calcite tests and both groups demonstrate the previously described anomaly with enrichment of heavy isotopes around 3 - 4 °C (Gussone and Filipsson, 2010). Calcium isotope ratios of the aragonitic species Hoeglundina elegans, on the other hand, are about 0.4 per mil lighter compared to the calcitic species, which is in general agreement with stronger fractionation in inorganic aragonite compared to calcite. However, the low and strongly variable Sr content suggests additional processes during test formation, and we propose that transmembrane ion transport or a precursor phase to aragonite may be involved. Porcelaneous tests, composed of high Mg calcite, incorporate higher amounts of Sr compared to hyaline low Mg calcite, in agreement with inorganic calcite systematics, but also porcelaneous tests with reduced Mg/Ca show high Sr/Ca. While calcium isotopes, Sr/Ca and Mg/Ca in benthonic foraminifers primarily appear to fractionate and partition with a dominant inorganic control, d44/40Ca temperature and growth rate dependencies of benthonic foraminifer tests favour a dominant contribution of light Ca by transmembrane transport relative to unfractionated seawater Ca to the calcifying fluid, thus controlling the formation of foraminiferal d44/40Ca and Sr/Ca proxy signals.

Mg/Ca and Sr/Ca ratios of Hoeglundina elegans from surface sediments

Core top samples from Atlantic (Little Bahama Banks (LBB)) and Pacific (Hawaii and Indonesia) depth transects have been analyzed in order to assess the influence of bottom water temperature (BWT) and aragonite saturation levels on Mg/Ca and Sr/Ca ratios in the aragonitic benthic foraminifer Hoeglundina elegans. Both the Mg/Ca and Sr/Ca ratios in H. elegans tests show a general decrease with increasing water depth. Although at each site the decreasing trends are consistent with the in situ temperature profile, Mg/Ca and Sr/Ca ratios in LBB are substantially higher than in Indonesia and Hawaii at comparable water depths with a greater difference observed with increasing water depth. Because we find no significant difference between results obtained on "live" and "dead" specimens, we propose that these differences are due to primary effects on the metal uptake during test formation. Evaluation of the water column properties at each site suggests that in situ CO3 ion concentrations play an important role in determining the H. elegans Mg/Ca and Sr/Ca ratios. The CO3 ion effect is limited, however, only to aragonite saturation levels ([DeltaCO3]aragonite) below 15 µmol/kg. Above this level, temperature exerts a dominant effect. Accordingly, we propose that Mg/Ca and Sr/Ca in H. elegans tests can be used to reconstruct thermocline temperatures only in waters oversaturated with respect to the mineral aragonite using the following relationships: Mg/Ca = (0.034 ± 0.002)BWT + (0.96 ± 0.03) and Sr/Ca = (0.060 ± 0.002)BWT + (1.53 ± 0.03) (for [DeltaCO3]aragonite > 15 µmol/kg). The standard error associated with these equations is about ±1.1°C. Reconstruction of deeper water temperatures is complicated because in undersaturated waters, changes in Mg/Ca and Sr/Ca ratios reflect a combination of changes in [CO3] and BWT. Overall, we find that Sr/Ca, rather than Mg/Ca, in H. elegans may be a more accurate proxy for reconstructing paleotemperatures.

Supplementary material: XRF core scanning data, CaCO3 values, foraminiferal carbon and oxygen isotope data, Mg/Ca data, clay mineralogy and helium isotope results

The Paleocene-Eocene Thermal Maximum (PETM; ~56 Ma) is associated with abrupt climate change, carbon cycle perturbation, ocean acidification, as well as biogeographic shifts in marine and terrestrial biota that were largely reversed as the climatic transient waned. We report a clear exception to the behavior of the PETM as a reversing climatic transient in the eastern North Atlantic (Deep-Sea Drilling Project Site 401, Bay of Biscay) where the PETM initiates a greatly prolonged environmental change compared to other places on Earth where records exist. The observed environmental perturbation extended well past the d13C recovery phase and up to 650 kyr after the PETM onset according to our extraterrestrial 3He-based age-model. We observe a strong decoupling of planktic foraminiferal d18O and Mg/Ca values during the PETM d13C recovery phase, which in combination with results from helium isotopes and clay mineralogy, suggests that the PETM triggered a hydrologic change in western Europe that increased freshwater flux and the delivery of weathering products to the eastern North Atlantic. This state change persisted long after the carbon-cycle perturbation had stopped. We hypothesize that either long-lived continental drainage patterns were altered by enhanced hydrological cycling induced by the PETM, or alternatively that the climate system in the hinterland area of Site 401 was forced into a new climate state that was not easily reversed in the aftermath of the PETM.

Calcite saturation and yield and Mg/Ca, Mn/Ca, Al/Ca, Fe/Ca for four species of planktic foraminifera cleaned using Mg-cleaning and Cd-cleaning methods

Four species of planktic foraminifera from core-tops spanning a depth transect on the Ontong Java Plateau were prepared for Mg/Ca analysis both with (Cd-cleaning) and without (Mg-cleaning) a reductive cleaning step. Reductive cleaning caused etching of foraminiferal calcite, focused on Mg-rich inner calcite, even on tests which had already been partially dissolved at the seafloor. Despite corrosion, there was no difference in Mg/Ca of Pulleniatina obliquiloculata between cleaning methods. Reductive cleaning decreased Mg/Ca by an average (all depths) of ~ 4% for Globigerinoides ruber white and ~ 10% for Neogloboquadrina dutertrei. Mg/Ca of Globigerinoides sacculifer (above the calcite saturation horizon only) was 5% lower after reductive cleaning. The decrease in Mg/Ca due to reductive cleaning appeared insensitive to preservation state for G. ruber, N. dutertrei and P. obliquiloculata. Mg/Ca of Cd-cleaned G. sacculifer appeared less sensitive to dissolution than that of Mg-cleaned. Mg-cleaning is adequate, but SEM and contaminants (Al/Ca, Fe/Ca and Mn/Ca) show that Cd-cleaning is more effective for porous species. A second aspect of the study addressed sample loss during cleaning. Lower yield after Cd-cleaning for G. ruber, G. sacculifer and N. dutertrei confirmed this to be the more aggressive method. Strongest correlations between yield and Delta[CO3^2-] in core-top samples were for Cd-cleaned G. ruber (r = 0.88, p = 0.020) and Cd-cleaned P. obliquiloculata (r = 0.68, p = 0.030). In a down-core record (WIND28K) correlation, r, between yield values > 30% and dissolution index, XDX, was -0.61 (p = 0.002). Where cleaning yield < 30% most Mg-cleaned Mg/Ca values were biased by dissolution.

Sedimentology and Cd/Ca ratios of benthic foraminifera from ODP Site 172-1060 in the western Nort Atlantic

The Western Boundary Undercurrent (WBUC), off eastern America, is an important component of the Atlantic Meridional Overturning circulation and is the principal route for southward transport of North Atlantic waters and southward return of Southern Source Water (SSW). Here a direct flow speed proxy (mean grain size of the sortable silt) is used to infer the vigour of flow of the palaeo-WBUC at Blake Outer Ridge, (ODP Site 1060, depth 3481 m) during Marine Isotope Stage (MIS) 3. The overall shape of the flow speed proxy record shows a complex pattern of variability, with generally more vigorous flow and larger-scale flow variations between 35 and 60 ka than in the younger part of MIS 3 and MIS 2 (b35 ka). Six events of reduced bottom flow vigour (Slow Events, SEs) occur. These appear uncorrelated with Heinrich events, but are instead synchronous with the warming phases of Antarctic Warm Events A-1 to A-4 (with one new one, A-1a and one poorly defined, 'A-0'). This indicates that Antarctic climate exerts a stronger control on deep flow vigour in the North Atlantic during MIS 3 than Northern Hemisphere climate. The correspondence of SEs with Antarctic warming suggests a weaker WBUC flow due to reduced volume flux at SSW source or reduced SSW density. Because the variability of the lower limb of the WBUC was not connected to sharp North Atlantic changes in temperature, it is unlikely that the Dansgaard/Oeschger cycles were associated with a mode of MOC variation involving wholeocean overturn, but more likely with perturbations of only the shallow Glacial Gulf Stream-Glacial Northern Source Intermediate Water cell. Nutrient proxies (benthic carbon isotopes and Cd/Ca of Uvigerina peregrina) at this site show similar trends to the GRIP delta18O record. This correlation has previously been attributed mainly to hydrographic and flow changes but is here shown to be better explained by variations in surface ocean productivity and subsequent decomposition of 12C rich organic material on the sea floor.

Stable oxygen isotope. Mg/Ca and Temperature record of Caribbean and Pacific sediments

The final phase of the closure of the Panamanian Gateway and the intensification of Northern Hemisphere Glaciation (NHG) both occurred during the Late Pliocene. Glacial-interglacial (G-IG) variations in sea level might, therefore, have had a significant impact on the remaining connections between the East Pacific and the Caribbean. Here, we present combined foraminiferal Mg/Ca and d18O measurements from Ocean Drilling Program (ODP) Site 1241 from the East Pacific and ODP Site 999 from the Caribbean. The studied time interval covers the first three major G-IG Marine Isotope Stages (MIS 95-100, ~2.5 Ma) after the intensification of NHG. Analyses were performed on the planktonic foraminifera Neogloboquadrina dutertrei and Globigerinoides sacculifer, representing water mass properties in the thermocline and the mixed-layer, respectively. Changes in sea water temperature, relative salinity, and water column stratification strongly suggest that the Panamanian Gateway temporarily closed during glacial MIS 98 and 100, as a result of changes in ice volume equivalent to a drop in sea level of 60-90 m. Reconstructed sea surface temperatures (SST) from G. sacculifer show a glacial decrease of 2.5°C at Site 1241, but increases of up to 3°C at Site 999 during glacial MIS 98 and 100 suggesting that the Panamanian Gateway closed during these glacial periods. The Mg/Ca-temperatures of N. dutertrei remain relatively stable in the East Pacific, but do show a 3°C warming in the Caribbean at the onset of these glacial periods suggesting that the closing of the gateway also changed the water column stratification. We infer that the glacial closure of the gateway allowed the Western Atlantic Warm Pool to extend into the southern Caribbean, increasing SST (G. sacculifer) and deepening the thermocline (N. dutertrei). Additionally, ice volume appears to have become large enough during MIS 100 to survive the relatively short lasting interglacial MIS 99 so that the gateway remained closed. Towards the end of MIS 98, during MIS 97 and into MIS 96 temperatures on both sides are mostly similar suggesting water masses exchanged again. Additionally, Caribbean variations in SST and d18Owater follow a precession-like cyclicity rather than the obliquity-controlled variations characteristic of the East-Pacific and many other tropical areas, suggesting that regional atmospheric processes related to the trade winds and the Intertropical Convergence Zone (ITCZ) had a dominant impact in the Caribbean.