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

Mg/Ca and Sr/Ca ratios of planktonic foraminifera from North Atlantic surface sediments

Reliable temperature estimates from both surface and subsurface ocean waters are needed to reconstruct past upper water column temperature gradients and past oceanic heat content. This work examines the relationships between trace element ratios in fossil shells and seawater temperature for surface-dwelling foraminifera species, Globigerinoides ruber (white) and Globigerina bulloides, and deep-dwelling species, Globorotalia inflata, Globorotalia truncatulinoides (dextral and sinistral) and Pulleniatina obliquiloculata. Mg/Ca and Sr/Ca ratios in shells picked in 29 modern core tops from the North Atlantic Ocean are calibrated using calculated isotopic temperatures. Mg/Ca ratios on G. ruber and G. bulloides agree with published data and relationships. For deep-dwelling species, Mg/Ca calibration follows the equation Mg/Ca = 0.78 (±0.04) * exp (0.051 (±0.003) * T) with a significant correlation coefficient of R**2 = 0.74. Moreover, there is no significant difference between the different deep-dwellers analyzed. For the Sr/Ca ratio, the surface dwellers and P. obliquiloculata do not record any temperature dependence. For the Globorotalia species, the thermo dependence of Sr/Ca ratio can be described by a single linear relationship: Sr/Ca = (0.0182 (±0.001) * T) + 1.097 (±0.018), R**2 = 0.85. Temperature estimates with a 1 sigma error of ±2.0°C and ±1.3°C can be derived from the Mg/Ca and Sr/Ca ratios, respectively, as long as the Sr geochemistry in the ocean has been constant through time.

(Table 1) Temperature and foraminiferal Mg/Ca ratios of samples off Japan

We studied Mg/Ca in high-Mg, shallow-water benthic foraminifera in culture and in samples from natural environments, in order to evaluate the expression of latitudinal and seasonal temperature variability in Mg/Ca in their tests. We cultured Planoglabratella opercularis (d'Orbigny) and Quinqueloculina yabei Asano under controlled temperature (10°-25°C) and salinity (30-38) conditions. Both species show a linear correlation between Mg/Ca and temperature, but they differ in temperature sensitivity. Salinity does not significantly influence Mg/Ca. In the samples collected in nature, Mg/Ca and seawater temperatures are positively correlated, but there are more complexities than in the records for cultured specimens due to such factors as seasonal fluctuations in temperature. We conclude that Mg/Ca ratios in monospecific benthic foraminiferal samples may be used as a reliable temperature proxy, if the lifetime of the species is taken into account.

Mg/Ca and stable oxygen isotope ratios of Uvigerina spp. from the Southern Ocean

The sensitivity to temperature of Mg/Ca ratios in the shallow-infaunal benthic foraminifera Uvigerina spp. has been assessed. Core-top calibrations over ~1-20 °C show a range in sensitivity of 0.065-0.084 mmol/mol/°C but few data are available spanning the temperature range anticipated in deep-sea records over glacial-interglacial cycles. In contrast to epibenthic foraminiferal species, carbonate ion saturation appears not to affect Mg/Ca significantly. A method based on estimating the ratio of the temperature sensitivity of foraminiferal Mg/Ca to that of d18Ocalcite shows that sensitivity for Mg/Ca at the high end of the observed core-top range (~0.1 mmol/mol/°C) is required for consistency with LGM-Holocene differences in each property as constrained by independent proxy data. This is supported by a Mg/Ca record for Uvigerina spp. generated for the Southern Ocean over the past 440,000 years from Ocean Drilling Program Site 1123 (Chatham Rise, New Zealand). The record shows variability that correlates with climate oscillations. The LGM deep ocean temperature derived from the Mg/Ca record is -1.1 ± 0.3 °C. Transformation to temperature allows estimates to be made of changes in bottom water temperature and seawater d18O and comparison made with literature records. Analysis reveals a ~2.5-kyr lead in the record of temperature over calcite d18O and a longer lead over seawater d18O. This is a reflection of larger phase offsets at eccentricity periods; phase offsets at tilt and precession are within error zero.

Mg/Ca ratio of Globorotalia tumida and subsurface temperature reconstruction for ODP Site 130-806

The tropical Pacific thermocline strength, depth, and tilt are critical to tropical mean state and variability. During the early Pliocene (~3.5 to 4.5 Ma), the Eastern Equatorial Pacific (EEP) thermocline was deeper and the cold tongue was warmer than today, which resulted in an mean state with a reduced zonal sea surface temperature gradient, or El Padre. However, it is unclear whether the deep thermocline was a local feature of the EEP or a basin-wide condition with global implications. Our measurements of Mg/Ca of Globorotalia tumida in a western equatorial Pacific site indicate Pliocene subsurface temperatures warmer than today; thus, El Padre included a basin-wide thermocline that was relatively warm, deep, and weakly tilted. At ~4 Ma, thermocline steepening was coupled to cooling of the cold tongue. Since ~4 Ma, the basin-wide thermocline cooled/shoaled gradually, with implications for thermocline feedbacks in tropical dynamics and the interpretation of TEX86-derived temperatures.

Sr/Ca and Cd/Ca ratios of benthic foraminifera in surface deep-sea sediments

Measurements of Sr/Ca of benthic foraminifera show a linear decrease with water depth which is superimposed upon significant variability identified by analyses of individual foraminifera. New data for Cd/Ca support previous work in defining a contrast between waters shallower and deeper than ~2500 m. Measured element partition coefficients in foraminiferal calcium carbonate relative to sea water (D) have been described by means of a one-box model in which elements are extracted by Rayleigh distillation from a biomineralization reservoir that serves for calcification with a constant fractionation factor (alpha), such that D = (1 - f**alpha)/(l - f), where f is the proportion of Ca remaining after precipitation. A modification to the model recognises differences in element speciation. The model is consistent with differences between D[Sr], D[Ba], and D[Cd] in benthic but not planktonic foraminifera. Depth variations in D for Sr and Ba are consistent with the model, as are differences in depth variation of D[Cd] in calcitic and aragonitic benthic foraminifera. The shallower depth variations may reflect increasing calcification rates with increasing water depth to an optimum of about 2500 m. Observations of unusually lower DCd for some deep waters, not accompanied by similar [Sr], or D[Ba] may be because of dissolution or a calcification response to a lower carbonate saturation state.

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.

Potassium and argon concentrations and age determination of ODP Site 125-782 and 125-786

K-Ar whole-rock ages have been obtained for 30 samples from Sites 782 and 786, Ocean Drilling Program Leg 125 in the Izu-Bonin (Ogasawara) forearc region. They form a trimodal spread of ages between 9 Ma and 44 Ma and are, with a few exceptions, consistent with the inferred lithostratigraphy. The ages have been interpreted in terms of at least two distinct episodes of magmatic and/or hydrothermal activity. A group of ten samples, including the lava flows, gave an isochron age of 41.3 ± 0.5 Ma (middle-late Eocene). This is thought to represent the age of the principal magmatic development of the volcanic forearc basement, and is comparable to published ages on equivalent rocks from other parts of the forearc basement high (e.g., the Ogasawara Islands). It may be significant that this age is slightly younger than the timing of major plate reorganization in the Western Pacific at about 43 Ma. This was followed by a minor episode of intrusive magmatism at 34.6 ± 0.7 Ma (early Oligocene) which appears to have reset the ages of some of the earlier units. This event probably corresponds to the initiation of rifting of the "proto-arc" to form the Parece Vela Basin. Boninitic samples were erupted during both episodes of magmatism, the earlier being of low-Ca boninite type and the later being of medium- and high-Ca types. It is also possible that a third episode of intrusive magmatism affected the Izu-Bonin forearc region at both Sites 782 and 786 at about 17 Ma. This would be consistent with magmatic activity elsewhere in the region during the Miocene, associated with the end of active spreading in the Parece Vela Basin and the start of arc activity in the West Mariana Ridge.

Planktonic foraminiferal Cd/Ca ratios

Records of Cd/Ca in planktonic foraminiferal calcite of Globigerinoides bulloides in cores from the Subantarctic region of the Southern Ocean show large glacial-interglacial variations with lower Cd/Ca (by 0.06-0.10 µmol/mol) at glacial times. Interpretation of these records in terms of lower dissolved phosphate and inferred higher glacial nutrient utilization has significant implications for glacial atmospheric carbon dioxide (pCO2) draw-down. However, box core-top data for G. bulloides in the North Atlantic suggest that the incorporation of Cd into planktonic foraminifera relative to seawater (DCd) is temperature sensitive (DCd=0.637 exp 0.15T). When the Subantarctic planktonic Cd/Ca records are corrected for this temperature dependence, they show little or no glacial-interglacial diferences. If, as seems likely, this observation can be interpreted to indicate a minimal change (< 0.5 µmol/kg) in surface water phosphate concentrations, then the explanation for lowered glacial pCO2 must be looked for elsewhere.

Magnetic properties of ODP Site 175-1075

We present sediment magnetic and chemical analysis of cyclic ocean sediments of the upwelling region of the Lower Congo Basin (equatorial Atlantic). We investigated two >100-k.y. intervals from Ocean Drilling Program Site 1075 to analyze the hysteresis properties, sources of magnetic susceptibility, anhysteretic remanent magnetizations, thermomagnetic behavior, and element concentrations of Fe, Ca, Ti, Mn, and K using an X-ray fluorescence (XRF) core scanner. The upper interval was sampled between 14 and 32 meters composite depth (mcd; 0.09-0.21 Ma) and the lower between 141 and 163 mcd (1.31-1.54 Ma) at a resolution of 20 cm, which represents a temporal resolution of 2.0 and 1.3 k.y., respectively. XRF core-scanner data were acquired at 5-cm intervals. The measurements show that ferri(o)magnetic minerals have no significant influence on the cyclicity of the magnetic susceptibility, which is dominated by paramagnetic and diamagnetic minerals and reflects changes of sediment input from the Congo River. The Fe, Ti, K, and Mn concentrations covary with the magnetic susceptibility where high concentrations of these elements correlate with intervals of high susceptibility and low concentrations with intervals of low susceptibility. The Ca counts correlate well with the calcium carbonate concentration but do not show the same cyclicity as the other elements or the susceptibility. With the exception of the Ca concentration, which is significantly higher in the upper interval, and the magnetic grain size, which indicates that less fine grained magnetite is present in the lower interval, no significant differences in the properties of the upper and the lower intervals were detected.

Primary composition, alteration, and origin of Cretaceous volcaniclastic rocks at DSDP Site 89-585

An upper Aptian to middle Albian series of volcaniclastic rocks more than 300 m thick was drilled at Site 585 in the East Mariana Basin. On the basis of textural and compositional (bulk-rock chemistry, primary and secondary mineral phases) evidence, the volcaniclastic unit is subdivided into a lower (below 830 m sub-bottom) and an upper (about 670-760 m) sequence; the boundary in the interval between is uncertain owing to lack of samples. The rocks are dominantly former vitric basaltic tuffs and minor lapillistones with lesser amounts of crystals and basaltic lithic clasts. They are mixed with shallow-water carbonate debris (ooids, skeletal debris), and were transported by mass flows to their site of deposition. The lower sequence is mostly plagioclase- and olivine-phyric with lesser amounts of Ti-poor clinopyroxene. Mineralogical and bulk-rock chemical data indicate a tholeiitic composition slightly more enriched than N-MORB (normal mid-ocean ridge basalt). Transport was by debris flows from shallow-water sites, as indicated by admixed ooids. Volcanogenic particles are chiefly moderately vesicular to nonvesicular blocky shards (former sideromelane) and less angular tachylite with quench plagioclase and pyroxene, indicating generation of volcanic clasts predominantly by spalling and breakage of submarine pillow and/or sheet-flow lavas. The upper sequence is mainly clinopyroxene- and olivine-phyric with minor plagioclase. The more Ti-rich clinopyroxene and the bulk-rock analyses show that the moderately alkali basaltic composition throughout is more mafic than the basal tholeiitic sequence. Transport was by turbidity currents. Rounded epiclasts of crystalline basalts are more common than in the lower sequence, and, together with the occurrence of oxidized olivine pseudomorphs and vesicular tachylite, are taken as evidence of derivation from eroded subaerially exposed volcanics. Former sideromelane shards are more vesicular than in the lower sequence; vesicularity exceeds 60 vol.% in some clasts. The dominant clastic process is interpreted to be by shallow-water explosive eruptions. All rocks have undergone low-temperature alteration; the dominant secondary phases are "palagonite," chlorite/smectite mixed minerals, analcite, and chabazite. Smectite, chlorite, and natrolite occur in minor amounts. Phillipsite is recognized as an early alteration product, now replaced by other zeolites. During alteration, the rocks have lost up to 50% of their Ca, compared with a fresh shard and fresh glass inclusions in primary minerals, but have gained much less K, Rb, and Ba than expected, indicating rapid deposition prior to significant seafloor weathering.

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.

Carbon chemistry of oceanic crust

Carbonate mineral precipitation in the upper oceanic crust during low-temperature, off-axis, hydrothermal circulation is investigated using new estimates of the bulk CO2 content of seven DSDP/ODP drill cores. In combination with previously published data these new data show: (i) the CO2 content of the upper ~ 300 m of the crust is substantially higher in Cretaceous than in Cenozoic crust and (ii) for any age of crust, there is substantially more CO2 in Atlantic (slow-spreading) than Pacific (intermediate- to fast-spreading) crust. Modelling the Sr-isotopic composition of the carbonates suggests that > 80% of carbonate mineral formation occurs within < 20 Myr of crust formation. This means that the higher CO2 content of Cretaceous crust reflects a secular change in the rate of CO2 uptake by the crust. Oxygen isotope derived estimates of carbonate mineral precipitation temperatures show that the average and minimum temperature of carbonate precipitation was ~10 °C higher temperatures in the Cretaceous than in the Cenozoic. This difference is consistent with previous estimates of secular change in bottom seawater temperature. Higher fluid temperature within the crust will have increased reaction rates potentially liberating more basaltic Ca and hence enhancing carbonate mineral precipitation. Additionally, if crustal fluid pH is controlled by fluid-rock reaction, the higher Ca content of the Cretaceous ocean will also have enhanced carbonate mineral precipitation. New estimates of the rate of CO2 uptake by the upper ocean crust during the Cenozoic are much lower than previous estimates.