Geochemistry and Mg-isotopes of carbonate sediments and porefluids of ODP Site 130-807

Deep-sea pore fluids are potential archives of ancient seawater chemistry. However, the primary signal recorded in pore fluids is often overprinted by diagenetic processes. Recent studies have suggested that depth profiles of Mg concentration in deep-sea carbonate pore fluids are best explained by a rapid rise in seawater Mg over the last 10-20 Myr. To explore this possibility we measured the Mg isotopic composition of pore fluids and carbonate sediments from Ocean Drilling Program (ODP) site 807. Whereas the concentration of Mg in the pore fluid declines with depth, the isotopic composition of Mg in the pore fluid increases from -0.78 per mil near the sediment-water interface to -0.15 per mil at 778 mbsf. The Mg isotopic composition of the sediment, with few important exceptions, does not change with depth and has an average d26Mg value of -4.72 per mil. We reproduce the observed changes in sediment and pore-fluid Mg isotope values using a numerical model that incorporates Mg, Ca and Sr cycling and satisfies existing pore-fluid Ca isotope and Sr data. Our model shows that the observed trends in magnesium concentrations and isotopes are best explained as a combination of two processes: a secular rise in the seawater Mg over the Neogene and the recrystallization of low-Mg biogenic carbonate to a higher-Mg diagenetic calcite. These results indicate that burial recrystallization will add Mg to pelagic carbonate sediments, leading to an overestimation of paleo-temperatures from measured Mg/Ca ratios. The Mg isotopic composition of foraminiferal calcite appears to be only slightly altered by recrystallization making it possible to reconstruct the Mg isotopic composition of seawater through time.

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.

Sample locations, CTD data, mean spring and annual temperatures and salinities from World Ocean Atlas (WOA, 2001), and ICP-OES Mg/Ca data

Due to its strong gradient in salinity and small temperature gradient the Mediterranean provides an ideal setting to study the impact of salinity on the incorporation of Mg into foraminiferal tests. We have investigated tests of Globorotalia inflata and Globigerina bulloides in plankton tow and core top samples from the Western Mediterranean using ICP-OES for bulk analyses and LA-ICP-MS for analyses of individual chambers in single specimens. Mg/Ca observed in G. inflata are consistent with existing calibrations, whereas G. bulloides had significantly higher Mg/Ca than predicted, particularly in core top samples from the easterly stations. Scanning Electron Microscopy and Laser Ablation ICP-MS revealed secondary overgrowths on some tests, which could explain the observed high core top Mg/Ca. We suggest that the Mediterranean intermediate and deep water supersaturated with respect to calcite cause these overgrowths and therefore increased bulk Mg/Ca. However, the different species are influenced by diagenesis to different degrees probably due to different test morphologies. Our results provide new perspectives on reported anomalously high Mg/Ca in sedimentary foraminifera and the applicability of the Mg/Ca paleothermometry in high salinity settings, by showing that (1) part of the signal is generated by precipitation of inorganic calcite on the foraminifer test due to increased calcite saturation state of the water and (2) species with high surface-to-volume shell surfaces are potentially more affected by secondary Mg-rich calcite encrustation.

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.

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.

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