Biostratigraphy and Sr isotopic composition of sediments from DSDP Legs 12, 14, 15, 17, 20-22, and 30

We present the data used to construct the Cenozoic and Cretaceous portion of the Phanerozoic curve of seawater 87Sr/86Sr that had been given in summary form by W.H. Burke co-workers. All Cenozoic samples (128) and 22 Cretaceous samples are foram-nannofossil oozes and limestones from DSDP cores distributed among 13 sites in the Atlantic, Pacific and Indian Oceans, and the Caribbean Sea. Non-DSDP Cretaceous samples (126) include limestone, anhydrite and phosphate samples from North America, Europe and Asia. Determination of the 87Sr/86Sr value of seawater at particular times in the past is based on comparison of ratios derived from coeval marine samples from widely separated geographic areas. These samples are characterized by a wide variety of diagenetic and burial histories. The large size and cosmopolitan nature of the data set decreases the likelihood that, among coeval data, systematic error has been introduced by a similar pattern of diagenetic alteration of the ratios. There is good clustering of data points throughout the Cenozoic and Cretaceous curve. The consistency of data is illustrated by Cenozoic and Cretaceous data plots that include a separate symbol for each DSDP site and non-DSDP sample location. More than 98% of the data points are enclosed by upper and lower lines that define a narrow band. For any given time, the correct seawater ratio probably lies within this band. A line drawn within the band represents our estimate of the actual seawater ratio as a function of time. The general configuration of the Cenozoic and Cretaceous curve appears to be strongly influenced by the history of plate interactions and sea-floor spreading. Specific rises and falls in the 87Sr/86Sr of seawater, however, may be caused by a variety of factors such as variation in lithologic composition of the crust exposed to weathering, configuration and topographic relief of continents, volcanic activity, rate of sea-floor spreading, extent of continental inundation by epeiric seas, and variations in both climate and paleooceanographic conditions. Many or all of these factors are probably related to global tectonic processes, yet their combined effect on the temporal variation of seawater 87Sr/86Sr can complicate a direct platetectonic interpretation for portions of the seawater curve.

Pliocene stable oxygen and carbon isotope record of benthic foraminifera from ODP Site 138-846 in the eastern equatorial Pacific

Oxygen and carbon isotope ratios were measured in benthic foraminifers from the entire Pliocene and latest Miocene sections of Site 846, a 180-m section, at a sampling interval of 10 cm. This provides a temporal resolution of about 2500 yr. The documented continuity of the record is excellent. Using the time scale that was developed on the basis of orbital tuning of GRAPE density records, we observed a fairly constant phase relationship between delta18O and variations in the obliquity of Earth's rotational axis. A new numbering scheme for Pliocene isotope stages is proposed. This high-resolution delta18O record clarifies several interesting aspects of late Neogene climatic evolution, including a "glacial" event that may have caused the final Messinian desiccation of the Mediterranean Sea; one or more "interglacial" events that might have caused refilling of the Mediterranean; a well-resolved couplet of glacial events at about the age of the Sidujfall Subchron; interglacial extremes in the early part of the Gauss that could have resulted from either significant deglaciation on Antarctica or from warming of deep water; and a gradual ramp of increasingly extreme "glacial" events, starting at about the Kaena Subchron and culminating with delta18O stage 100 in the earliest Matuyama.

Chemical and Sr isotopic compositions and timing of carbonate mineral precipitation, ODP Leg 115

The strontium isotope ratios of authigenic carbonates from Indian Ocean sea-floor basalts have been used to determine the timing of carbonate mineral precipitation and fluid flow. The samples include calcites from 57.2 Ma crust from Ocean Drilling Project (ODP) Site 715, and calcites, aragonites, and siderites from 63.7 Ma crust from ODP Site 707. At Site 715, calcite precipitation may have begun at any time after the basalts cooled, and it continued until approximately 31 Ma, or 26 m.y. after basalt eruption. At Site 707, aragonite and siderite did not begin to precipitate until about 36 Ma, almost 30 m.y. after basalt eruption, and continued to precipitate until at least 30 and 28 Ma, respectively. Calcite precipitation began at approximately 32 Ma and continued until 22 Ma. These ages suggest that vein mineral deposition and low-temperature fluid circulation in the ocean crust may continue for much longer periods of time than previously observed.

Age determination and stable oxygen isotope record of ODP Hole 167-1017E

Oxygen isotopic (d18O) climatic stratigraphy and radiocarbon chronology, at high resolution, have been used to establish an age model for Ocean Drilling Program Hole 1017E, a continuous 25-m sequence of hemipelagic sediments from the continental slope (956 m water depth), east of Point Arguella, Southern California. The upper part of Hole 1017E from ~33 ka (7.445 mbsf) was dated using 13 calendar-corrected radiocarbon ages of mixed planktonic foraminiferal assemblages. Benthic oxygen isotopic stratigraphy records a continuous 130-k.y. sequence ranging from marine isotope Stage 6 to the present day. The benthic d18O curve, representing the last two interglacial and glacial cycles, closely resembles the well-dated, deep-sea reference sequence, providing a detailed chronologic framework. Sedimentation rates remained relatively constant throughout the sequence at ~18 cm/k.y. and were sufficiently rapid to provide considerable potential for high-resolution paleoceanographic/paleoclimatic investigations. Planktonic foraminiferal oxygen isotopic stratigraphy based on the surface-dwelling form Globigerina bulloides defines an almost complete sequence of interstadial/stadial oscillations (Dansgaard/Oeschger cycles [D/O]). Combined use of radiocarbon chronology, deep-sea oxygen isotopic datums, and visual pattern matching has enabled us to identify the sequence of D/O cycles as described for the Greenland (GRIP2) ice core. This has strengthened the stratigraphic framework for the last 60 k.y. in the sequence as a basis for further paleoenvironmental investigations.

Chemical and Sr isotopic compositions of pore fluids and marine carbonates from ODP Hole 130-807A

The 87Sr/86Sr ratios and Sr concentrations in sediment and pore fluids are used to evaluate the rates of calcite recrystallization at ODP Site 807A on the Ontong Java Plateau, an 800-meter thick section of carbonate ooze and chalk. A numerical model is used to evaluate the pore fluid chemistry and Sr isotopes in an accumulating section. The deduced calcite recrystallization rate is 2% per million years (%/Myr) near the top of the section and decreases systematically in older parts of the section such that the rate is close to 0.1/age (in years). The deduced recrystallization rates have important implications for the interpretation of Ca and Mg concentration profiles in the pore fluids. The effect of calcite recrystallization on pore fluid chemistry is described by the reaction length, L, which varies by element, and depends on the concentration in pore fluid and solid. When L is small compared to the thickness of the sedimentary section, the pore fluid concentration is controlled by equilibrium or steady-state exchange with the solid phase, except within a distance L of the sediment-water interface. When L is large relative to the thickness of sediment, the pore fluid concentration is mostly controlled by the boundary conditions and diffusion. The values of L for Ca, Sr, and Mg are of order 15, 150, and 1500 meters, respectively. L_Sr is derived from isotopic data and modeling, and allows us to infer the values of L_Ca and L_Mg. The small value for L_Ca indicates that pore fluid Ca concentrations, which gradually increase down section, must be equilibrium values that are maintained by solution-precipitation exchange with calcite and do not reflect Ca sources within or below the sediment column. The pore fluid Ca measurements and measured alkalinity allow us to calculate the in situ pH in the pore fluids, which decreases from 7.6 near the sediment-water interface to 7.1+/-0.1 at 400-800 mbsf. While the calculated pH values are in agreement with some of the values measured during ODP Leg 130, most of the measurements are artifacts. The large value for L_Mg indicates that the pore fluid Mg concentrations at 807A are not controlled by calcite-fluid equilibrium but instead are determined by the changing Mg concentration of seawater during deposition, modified by aqueous diffusion in the pore fluids. We use the pore fluid Mg concentration profile at Site 807A to retrieve a global record for seawater Mg over the past 35 Myr, which shows that seawater Mg has increased rapidly over the past 10 Myr, rather than gradually over the past 60 Myr. This observation suggests that the Cenozoic rise in seawater Mg is controlled by continental weathering inputs rather than by exchange with oceanic crust. The relationship determined between reaction rate and age in silicates and carbonates is strikingly similar, which suggests that reaction affinity is not the primary determinant of silicate dissolution rates in nature.

Os and Re concentrations and Os isotope ratios for DSDP and ODP metalliferous sediments from the Pacific Ocean (Table 1)

We report new 187Os/186Os data and Re and Os concentrations in metalliferous sediments from the Pacific to construct a composite Os isotope seawater evolution curve over the past 80 m.y. Analyses of four samples of upper Cretaceous age yield 187Os/186Os values of between 3 and 6.5 and 187Re/186Os values below 55. Mass balance calculations indicate that the pronounced minimum of about 2 in the Os isotope ratio of seawater at the K-T boundary probably reflects the enormous input of cosmogenic material into the oceans by the K-T impactor(s). Following a rapid recovery to 187Os/186Os of 3.5 at 63 Ma, data for the early and middle part of the Cenozoic show an increase in 187Os/186Os to about 6 at 15 Ma. Variations in the isotopic composition of leachable Os from slowly accumulating metalliferous sediments show large fluctuations over short time spans. In contrast, analyses of rapidly accumulating metalliferous carbonates do not exhibit the large oscillations observed in the pelagic clay leach data. These results together with sediment leaching experiments indicate that dissolution of non-hydrogenous Os can occur during the hydrogen peroxide leach and demonstrate that Os data from pelagic clay leachates do not always reflect the Os isotopic composition of seawater. New data for the late Cenozoic further substantiate the rapid increase in the 187Os/186Os of seawater during the past 15 Ma. We interpret the correlation between the marine Sr and Os isotope records during this time period as evidence that weathering within the drainage basin of the Ganges-Brahmaputra river system is responsible for driving seawater Sr and Os toward more radiogenic isotopic compositions. The positive correlation between 87Sr/86Sr and U concentration, the covariation of U and Re concentrations, and the high dissolved Re, U and Sr concentrations found in the Ganges-Brahmaputra river waters supports this interpretation. Accelerating uplift of many orogens worldwide over the past 15 Ma, especially during the last 5 Ma, could have contributed to the rapid increase in 187Os/186Os from 6 to 8.5 over the past 15 Ma. Prior to 15 Ma the marine Sr and Os record are not tightly coupled. The heterogeneous distribution of different lithologies within eroding terrains may play an important role in decoupling the supplies of radiogenic Os and Sr to the oceans and account for the periods of decoupling of the marine Sr and Os isotope records.

Thallium isotope composition in pyrites

This paper presents the first study of Tl isotopes in early diagenetic pyrite. Measurements from two sections deposited during the Toarcian Ocean Anoxic Event (T-OAE, ~183 Ma) are compared with data from Late Neogene (<10 Ma) pyrite samples from ODP legs 165 and 167 that were deposited in relatively oxic marine environments. The Tl isotope compositions of Late Neogene pyrites are all significantly heavier than seawater, which most likely indicates that Tl in diagenetic pyrite is partially sourced from ferromanganese oxy-hydroxides that are known to display relatively heavy Tl isotope signatures. One of the T-OAE sections from Peniche in Portugal displays pyrite thallium isotope compositions indistinguishable from Late Neogene samples, whereas samples from Yorkshire in the UK are depleted in the heavy isotope of Tl. These lighter compositions are best explained by the lack of ferromanganese precipitation at the sediment-water interface due to the sulfidic (euxinic) conditions thought to be prevalent in the Cleveland Basin where the Yorkshire section was deposited. The heavier signatures in the Peniche samples appear to result from an oxic water column that enabled precipitation of ferromanganese oxy-hydroxides at the sediment-water interface. The Tl isotope profile from Yorkshire is also compared with previously published molybdenum isotope ratios determined on the same sedimentary succession. There is a suggestion of an anti-correlation between these two isotope systems, which is consistent with the expected isotope shifts that occur in seawater when marine oxic (ferromanganese minerals) fluxes fluctuate. The results outlined here represent the first evidence that Tl isotopes in early diagenetic pyrite have potential to reveal variations in past ocean oxygenation on a local scale and potentially also for global oceans. However, much more information about Tl isotopes in different marine environments, especially in anoxic/euxinic basins, is needed before Tl isotopes can be confidently utilized as a paleo-redox tracer.

1.2 Ma stacked benthic foraminiferal carbon isotope record

Pleistocene stable carbon isotope (d13C) records from surface and deep dwelling foraminifera in all major ocean basins show two distinct long-term carbon isotope fluctuations since 1.00 Ma. The first started around 1.00 Ma and was characterised by a 0.35 per mil decrease in d13C values until 0.90 Ma, followed by an increase of 0.60 per mil lasting until 0.50 Ma. The subsequent fluctuation started with a 0.40 per mil decrease between 0.50 and 0.25 Ma, followed by an increase of 0.30 per mil between 0.25 and 0.10 Ma. Here, we evaluate existing evidence and various hypotheses for these global Pleistocene d13C fluctuations and present an interpretation, where the fluctuations most likely resulted from concomitant changes in the burial fluxes of organic and inorganic carbon due to ventilation changes and/or changes in the production and export ratio. Our model indicates that to satisfy the long-term 'stability' of the Pleistocene lysocline, the ratio between the amounts of change in the organic and inorganic carbon burial fluxes would have to be close to a 1:1 ratio, as deviations from this ratio would lead to sizable variations in the depth of the lysocline. It is then apparent that the mid-Pleistocene climate transition, which, apart from the glacial cycles, represents the most fundamental change in the Pleistocene climate, was likely not associated with a fundamental change in atmospheric pCO2. While recognising that high frequency glacial/interglacial cycles are associated with relatively large (100 ppmv) changes in pCO2, our model scenario (with burial changes close to a 1:1 ratio) produces a maximum long-term variability of only 20 ppmv over the fluctuation between 1.00 and 0.50 Ma.