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.

Stable oxygen and carbon isotope record of foraminifera of ODP SIte 184-1146

Ocean Drilling Program Site 1146 was drilled within a small rift basin on the midcontinental slope of the northern South China Sea. It is located at 19°27.4'N, 116°16.37'E, in 2092 m water depth. This site was drilled to recover records of Asian monsoon variability into the middle Miocene with temporal resolution sufficient for orbital-scale analyses. Here we present oxygen and carbon isotopic measurements of planktonic foraminifers (Globigerinoides ruber) and benthic foraminifers (Uvigerina peregrina and Cibicides wuellerstorfi) as well as a preliminary age model for the top 185 meters composite depth (mcd).

Age model and Plio-Pleistocene benthic stable oxygen isotope ratios of ODP Site 184-1143 in the Sourh China Sea

Based on benthic foraminiferal delta18O from ODP Site 1143, a 5-Myr astronomical timescale for the West Pacific Plio-Pleistocene was established using an automatic orbital tuning method. The tuned Brunhes/Matuyama paleomagnetic polarity reversal age agrees well with the previously published age of 0.78 Ma. The tuned ages for several planktonic foraminifer bio-events also agree well with published dates, and new ages for some other bio-events in the South China Sea were also estimated. The benthic delta18O from Site 1143 is highly coherent with the Earth's orbit (ETP) both at the obliquity and precession bands for the last 5 Myr, and at the eccentricity band for the last 2 Myr. In general, the 41-kyr cycle was dominant through the Plio-Pleistocene although the 23-kyr cycle was also very strong. The 100-kyr cycle became dominant only during the last 1 Myr. A comparison of the benthic delta18O between the Atlantic (ODP 659) and the East and West Pacific (846 and 1143) reveals that the Atlantic-Pacific benthic oxygen isotope difference ratio (Delta delta18OAtl-Pac) displays an increasing trend in three time intervals: 3.6-2.7 Ma, 2.7-2.1 Ma and 1.5-0.25 Ma. Each of the intervals begins with a rapid negative shift in Delta delta18OAtl-Pac, followed by a long period with an increasing trend, corresponding to the growth of the Northern Hemisphere ice sheet. This means that all three intervals of ice sheet growth in the Northern Hemisphere were accompanied at the beginning by a rapid relative warming of deep water in the Atlantic as compared to that of the Pacific, followed by its gradual relative cooling. This general trend, superimposed on the frequent fluctuations with glacial cycles, should yield insights into the processes leading to the boreal glaciation. Cross-spectral analyses of the Delta delta18OAtl-Pac with the Earth's orbit suggests that after the initiation of Northern Hemisphere glaciation at about 2.5 Ma, obliquity rather than precession had become the dominant force controlling the vertical structure or thermohaline circulation in the paleo-ocean.

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.

Chemical and isotopic composition of interstitial waters and hydrocarbon composition at ODP Site 184-1146

Bottom-simulating reflectors were observed beneath the southeastern slope of the Dongsha Islands in the South China Sea, raising the potential for the presence of gas hydrate in the area. We have analyzed the chemical and isotopic compositions of interstitial water, headspace gas, and authigenic siderite concretions from Site 1146. Geochemical anomalies, including a slight decrease of chlorine concentration in interstitial water, substantial increase of methane concentration in headspace gas, and 18O enrichment in the authigenic siderite concretion below 400 meters below seafloor are probably caused by the decomposition of gas hydrate. The low-chlorine pore fluids contain higher molecular-weight hydrocarbons and probably migrate to Site 1146 along faults or bedded planes.

Benthic foraminiferal faunas of Miocene to Holocene sediments of ODP Site 184-1143 from the southern South China Sea

Deep-sea benthic foraminiferal assemblages from Ocean Drilling Program (ODP) Site 1143 located in the southern South China Sea (SCS) were investigated to evaluate the relationship between faunal composition patterns and paleoceanographic changes during the last 6 million years (late Miocene to Holocene). We used multivariate statistics (correspondence analysis) to analyze carbon-flux-related changes in assemblage composition of benthic foraminifers. Additional proxies for carbon flux and deep-water ventilation include delta13C records of epifaunal Cibicidoides wuellerstorfi and infaunal Uvigerina peregrina var. dirupta and Melonis pompilioides, benthic foraminiferal accumulation rates (BFARs), diversity indices, and relative abundances of indicator species. We observe three significant benthic faunal changes in the southern South China Sea during the last 6 million years. Strong fluctuations in BFAR and relative abundance of productivity indicator species between glacial and interglacial stages after the mid-Pleistocene revolution (MPR) at approximately 0.9 Ma, indicating stronger seasonal carbon flux fluctuations, are accompanied by the extinction of such species as Stilostomella spp. Increases in carbon flux indicator species are coupled with an overall decrease in benthic foraminifer diversity around 3.0 Ma in the late Pliocene. This may indicate increasing carbon flux in a period of productivity maximum caused by enhanced offshore upwelling from intensified winter monsoon wind strength.

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.

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.

Benthic foraminiferal Mg/Ca ratios

Core-top samples from different ocean basins have been analyzed to refine our current understanding of the sensitivity of benthic foraminiferal calcite magnesium/calcium (Mg/Ca) to bottom water temperatures (BWT). Benthic foraminifera collected from Hawaii, Little Bahama Bank, Sea of Okhotsk, Gulf of California, NE Atlantic, Ceara Rise, Sierra Leone Rise, the Ontong Java Plateau, and the Southern Ocean covering a temperature range of 0.8 to 18°C were used to revise the Cibicidoides Mg/Ca-temperature calibration. The Mg/Ca-BWT relationship of three common Cibicidoides species is described by an exponential equation: Mg/Ca = 0.867 ± 0.049 exp (0.109 ± 0.007 * BWT) (stated errors are 95% CI). The temperature sensitivity is very similar to a previously published calibration. However, the revised calibration has a significantly different preexponential constant, resulting in different predicted absolute temperatures. We attribute this difference in the preexponential constant to an analytical issue of accuracy. Some genera, notably Uvigerina, show apparently lower temperature sensitivity than others, suggesting that the use of constant offsets to account for vital effects in Mg/Ca may not be appropriate. Downcore Mg/Ca reproducibility, as determined on replicate foraminiferal samples, is typically better than 0.1 mmol/mol (2 S.E.). Thus, considering the errors associated with the Cibicidoides calibration and the downcore reproducibility, BWT may be estimated to within ±1°C. Application of the revised core-top Mg/Ca-BWT data to Cenozoic foraminiferal Mg/Ca suggests that seawater Mg/Ca was not more than 35% lower than today in the ice-free ocean at 50 Ma.

Pollen ages and depths in ODP Site 184-1144

A high-resolution pollen record (sampling interval averages 820 years) has been obtained from ODP Site 1144 (water depth 2037 m), northern South China Sea. The 504-m sequence (in composition length) covers the last 1.03 million years according to micropaleontological and isotopic stratigraphy. The pollen assemblages are characterized by high proportions of Pinus and herb pollen, and by their frequent alternations. Based on these alternations, 29 pollen zones have been recognized that are closely correlated to the Marine Oxygen Isotope Stages (MIS) 1-29. Pinus- dominant pollen zones correspond to interglacial periods with lighter delta18O values, while herb-marked ones relate to the heavier delta18O stages assigned to glacials. Judging from the pollen data, the exposed northern continental shelf of the South China Sea during the glacials was covered by grassland, and the extensive northern shelf has formed only since MIS 6 (ca. 150 ka), probably as a result of tectonic subsidence. Tree pollen influx values are indicative of winter monsoon which began to intensify 600 ka ago. The summer monsoon variations can be approximated by the fern percentage within the total pollen and spore abundance, and the result shows high values in general occurring at interglacials, with the maxima at MIS 15, 5e and 1. The relatively high fern percentage with smaller amplitude in variations before 600 ka might suggest more stable humid conditions before the intensification of winter monsoon.