Chlorine stable isotopic composition of basement fluids from ODP Leg 168 sites

This paper presents chlorine stable isotope compositions (delta37Cl) of sediment pore waters collected by squeezing sediment cores from the sediment-basement interface along an East-West transect through the eastern flank of the Juan de Fuca Ridge (ODP Leg 168). These "near basement fluids" (NBF) are generally thought to be representative of low-temperature fluids circulating in the off-axis basaltic crust. The delta37Cl value of the fluid directly sampled from a flow at the base of Site 1026 (WSTP1026) is also reported. NBF display delta37Cl values between -2.09? and -0.12? relative to the Standard Mean Ocean Chloride (SMOC defined as 0?) and small variations in chlorinity (~4%). These data contrast with the homogeneity of delta37Cl values associated with highly variable chlorinities observed in high-temperature on-axis fluids [Bonifacie et al., 2005, doi:10.1016/j.chemgeo.2005.06.008]. The NBF delta37Cl values show a general decreasing trend with distance from the ridge-axis except for two fluids. When plotted against delta18O values, the delta37Cl of the NBF show two different trends. This paper discusses the possible contributions on NBF delta37Cl values of fluid-mixing, water-rock interactions and transport processes (diffusion, ion membrane filtration) that can occur in the igneous basement. However, as none of these processes can fully explain the observed delta37Cl variations, the potential effect of the sediment cover is also investigated. At site 1026, the interstitial pore fluid displays a delta37Cl signature significantly lower than that of the fluid discharge sample (-1.90? and -0.28?, respectively). This difference, demonstrated here cannot be an artifact of the sampling method, rather indicates the influence of the sediment cover on NBF delta37Cl values. The potential contributions of physical processes associated with transport/compaction (e.g., diffusion, ion membrane filtration, adsorption, ion exchange) on NBF delta37Cl values are qualitatively discussed here but require additional studies for further insights. However, this study indicates that "near basement fluids" (NBF) are not, at least for Cl isotopic compositions, necessarily as representative of fluids circulating in the basaltic crust as initially thought. These results add new constraints on Cl geodynamics and show that Cl-isotopes fractionate during low-temperature circulation of fluids in off-axis and off-margin flow contexts, but not to the extent observed for active margins. Fluids circulating at low-temperature in the magmatic and/or the sedimentary part of the oceanic crust might have played a major role on the delta37Cl evolution of seawater over geologic time.

Isotope ratios and trace element concentrations of basalts from DSDP Leg 15 sites (Appendix 1)

Formation of the Cretaceous Caribbean plateau, including the komatiites of Gorgona, has been linked to the currently active Galápagos hotspot. We use Hf-Nd isotopes and trace element data to characterise both the Caribbean plateau and the Galápagos hotspot, and to investigate the relationship between them. Four geochemical components are identified in the Galápagos mantle plume: two 'enriched' components with epsilon-Hf and epsilon-Nd similar to enriched components observed in other mantle plumes, one moderately enriched component with high Nb/Y, and a fourth component which most likely represents depleted MORB source mantle. The Caribbean plateau basalt data form a linear array in Hf-Nd isotope space, consistent with mixing between two mantle components. Combined Hf-Nd-Pb-Sr-He isotope and trace element data from this study and the literature suggest that the more enriched Caribbean end member corresponds to one or both of the enriched components identified on Galápagos. Likewise, the depleted end member of the array is geochemically indistinguishable from MORB and corresponds to the depleted component of the Galápagos system. Enriched basalts from Gorgona partially overlap with the Caribbean plateau array in epsilon-Hf vs. epsilon-Nd, whereas depleted basalts, picrites and komatiites from Gorgona have a high epsilon-Hf for a given epsilon-Nd, defining a high-epsilon-Hf depleted end member that is not observed elsewhere within the Caribbean plateau sequences. This component is similar, however, in terms of Hf-Nd-Pb-He isotopes and trace elements to the depleted plume component recognised in basalts from Iceland and along the Reykjanes Ridge. We suggest that the Caribbean plateau represents the initial outpourings of the ancestral Galápagos plume. Absence of a moderately enriched, high Nb/Y component in the older Caribbean plateau (but found today on the island of Floreana) is either due to changing source compositions of the plume over its 90 Ma history, or is an artifact of limited sampling. The high-epsilon-Hf depleted component sampled by the Gorgona komatiites and depleted basalts is unique to Gorgona and is not found in the Caribbean plateau. This may be an indication of the scale of heterogeneity of the Caribbean plateau system; alternatively Gorgona may represent a separate oceanic plateau derived from a completely different Pacific plume, such as the Sala y Gomez.

Plio-Pleistocene age-depth points and occurrence of hiatuses in ODP sites of the Southern Ocean

Changes in Atlantic deep water circulation were reconstructed by comparing the benthic foraminiferal delta13C record at ODP Site 1090 in the South Atlantic with similar records from the North Atlantic (Sites 982, 607, 925, 929) and deep Pacific (Site 849) oceans. Important deep water circulation changes occurred in the early Pleistocene at 1.55 Myr and during the Mid-Pleistocene Transition at 0.9 Myr. At 1.55 Myr, glacial delta13C values in the Southern Ocean became significantly lower than those in the deep Pacific, establishing a pattern that persisted throughout the late Pleistocene. We propose that the lowering of delta13C values of Southern Component Water (SCW) at this time resulted from expansion of sea ice and reduced ventilation of deep water during glacial periods after marine isotope stage 52. Accompanying this change in Southern Ocean deep water circulation was enhanced interhemispheric coupling between the North and South Atlantic after 1.55 Myr. At ~0.9 Myr, the magnitude of glacial-to-interglacial variabilityin delta13C increased and shifted to a longer frequency (100 kyr) along with oceanic delta18O (ice volume). Calculation of percent Northern Component Water (NCW) using Site 1090 as the SCW end member yielded 20-30% less reduction of NCW during glacial periods of the late Pleistocene. Also, a trend toward reduced glacial suppression of NCW during the past 400 kyr is not evident. The apparent decoupling of ice volume and deep water circulation reported previously maybe an artifact of using a Pacific, rather than a Southern Ocean, carbon isotopic record to calculate past mixing ratios of NCW and SCW.

Stable isotope record, Mg/Ca ratios and age model of St. Stephens Quarry, Alabama

In the largest global cooling event of the Cenozoic Era, between 33.8 and 33.5 Myr ago, warm, high-CO2 conditions gave way to the variable 'icehouse' climates that prevail today. Despite intense study, the history of cooling versus ice-sheet growth and sea-level fall reconstructed from oxygen isotope values in marine sediments at the transition has not been resolved. Here, we analyse oxygen isotopes and Mg/Ca ratios of benthic foraminifera, and integrate the results with the stratigraphic record of sea-level change across the Eocene-Oligocene transition from a continental-shelf site at Saint Stephens Quarry, Alabama. Comparisons with deep-sea (Sites 522 (South Atlantic) and 1218 (Pacific)) d18O and Mg/Ca records enable us to reconstruct temperature, ice-volume and sea-level changes across the climate transition. Our records show that the transition occurred in at least three distinct steps, with an increasing influence of ice volume on the oxygen isotope record as the transition progressed. By the early Oligocene, ice sheets were ~25% larger than present. This growth was associated with a relative sea-level decrease of approximately 105 m, which equates to a 67 m eustatic fall.

Stable carbon isotope in bottom water and Cibicidoides wuellerstorfi on the Ontong Java Plateau and Emperor Seamount

We have measured the carbon isotopic composition of dissolved inorganic carbon in bottom waters of the Ontong Java Plateau (western equatorial Pacific) and on the northern Emperor Seamounts (northwest Pacific). Each of these locations is several hundred miles from the nearest Geochemical Ocean Sections Study (GEOSECS) stations, and the observed delta13C values at each site differ substantially from regionally averaged GEOSECS delta13C profiles. We discuss the possible causes of these differences, including horizontal variability, near-bottom effects, and problems with the Pacific GEOSECS delta13C data. We also measured the isotopic composition (C and O) of core top C. wuellerstorfi from a depth transect of cores at each location. The delta18O data are used to verify that our samples are Holocene. Comparison of foraminiferal and bottom water delta13C values shows that this species faithfully records bottom water delta13C at both sites and demonstrates that there is no depth-related artifact in the dissolved inorganic carbon-C. wuellerstorfi delta13C relationship at these sites.

Magnetic susceptibility, XRF and color reflectance data of DSDP Hole 72-516F, cores 113 and 114

Deep marine successions of early Campanian age from DSDP site 516F drilled at low paleolatitudes in the South Atlantic reveal distinct sub-Milankovitch variability in addition to precession and eccentricity related variations. Elemental abundance ratios point to a similar 5 climatic origin for these variations and exclude a quadripartite structure - as observed in the Mediterranean Neogene - of the precession related cycles as an explanation for the inferred semi-precession cyclicity in MS. However, the semi-precession cycle itself is likely an artifact, reflecting the first harmonic of the precession signal. The sub-Milankovitch variability is best approximated by a ~ 7 kyr cycle as shown by 10 spectral analysis and bandpass filtering. The presence of sub-Milankovitch cycles with a period similar to that of Heinrich events of the last glacial cycle is consistent with linking the latter to low-latitude climate change caused by a non-linear response to precession induced variations in insolation between the tropics.

Planktic foraminifera from the Paleocene-Eocene Thermal Maximum at Shatsky Rise, Pacific Ocean

High-resolution biostratigraphic and quantitative studies of subtropical Pacific planktonic foraminiferal assemblages (Ocean Drilling Program, Leg 198 Shatsky Rise, Sites 1209 and 1210) are performed to analyse the faunal changes associated with the Paleocene-Eocene Thermal Maximum (PETM) at about 55.5 Ma. At Shatsky Rise, the onset of the PETM is marked by the abrupt onset of a negative carbon isotope excursion close to the contact between carbonate-rich ooze and overlying clay-rich ooze and corresponds to a level of poor foraminiferal preservation as a result of carbonate dissolution. Lithology, planktonic foraminiferal distribution and abundances, calcareous plankton and benthic events, and the negative carbon isotope excursion allow precise correlation of the two Shatsky Rise records. Results from quantitative analyses show that Morozovella dominates the assemblages and that its maximum relative abundance is coincident with the lowest delta 13C values, whereas subbotinids are absent in the interval of maximum abundance of Morozovella. The excursion taxa (Acarinina africana, Acarinina sibaiyaensis, and Morozovella allisonensis) first appear at the base of the event. Comparison between the absolute abundances of whole specimens and fragments of genera demonstrate that the increase in absolute abundance of Morozovella and the decrease of Subbotina are not an artifact of selective dissolution. Moreover, the shell fragmentation data reveal Subbotina to be the more dissolution-susceptible taxon. The upward decrease in abundance of Morozovella species and the concomitant increase in test size of Morozovella velascoensis are not controlled by dissolution. These changes could be attributed to the species' response to low nutrient supply in the surface waters and to concomitant changes in the physical and chemical properties of the seawater, including increased surface stratification and salinity. Comparison of the planktonic foraminiferal changes at Shatsky Rise to those from other PETM records (Sites 865 and 690) highlights significant similarities, such as the decline of Subbotina at the onset of the event, and discrepancies, including the difference in abundance of the excursion taxa. The observed planktonic foraminifera species response suggests a warm-oligotrophic scenario with a high degree of complexity in the ocean structure.

Major ion concentrations and chlorine isotope composition of water-soluble chloride and structurally bound chloride in DSDP/ODP serpentinized ultramafic rocks

Eight DSDP/ODP cores were analyzed for major ion concentrations and d37Cl values of water-soluble chloride (d37Clwsc) and structurally bound chloride (d37Clsbc) in serpentinized ultramafic rocks. This diverse set of cores spans a wide range in age, temperature of serpentinization, tectonic setting, and geographic location of drilled serpentinized oceanic crust. Three of the cores were sampled at closely spaced intervals to investigate downhole variation in Cl concentration and chlorine isotope composition. The average total Cl content of all 86 samples is 0.26±0.16 wt.% (0.19±0.10 wt.% as water-soluble Cl (Xwsc) and 0.09±0.09 wt.% as structurally bound Cl (Xsbc)). Structurally bound Cl concentration nearly doubles with depth in all cores; there is no consistent trend in water-soluble Cl content among the cores. Chlorine isotope fractionation between the structurally bound Cl**- site and the water-soluble Cl**- site varies from -1.08? to +1.16?, averaging to +0.21?. Samples with negative fractionations may be related to reequilibration of the water-soluble chloride with seawater post-serpentinite formation. Six of the cores have positive bulk d37Cl values (+0.05? to +0.36?); the other two cores (173-1068A (Leg-Hole) and 84-570) have negative bulk d37Cl values (-1.26? and -0.54?). The cores with negative d37Cl values also have variable Cl**-/SO4**2- ratios, in contrast to all other cores. The isotopically positive cores (153-920D and 147-895E) show no isotopic variation with depth; the isotopically negative core (173-1068A) decreases by ~1? with depth for both the water-soluble and structurally bound Cl fractions. Non-zero bulk d37Cl values indicate Cl in serpentinites was incorporated during original hydration and is not an artifact of seawater infiltration during drilling. Cores with positive d37Cl values are most likely explained by open system fractionation during hydrothermal alteration, with preferential incorporation of 37Cl from seawater into the serpentinite and loss of residual light Cl back to the ocean. Fluid / rock ratios were probably low as evidenced by the presence of water-soluble salts. The two isotopically negative cores are characterized by a thick overlying sedimentary package that was in place prior to serpentinization. We believe the low d37Cl values of these cores are a result of hydration of ultramafic rock by infiltrating aqueous pore fluids from the overlying sediments. The resulting serpentinites inherit the characteristic negative d37Cl values of the pore waters. Chlorine stable isotopes can be used to identify the source of the serpentinizing fluid and ultimately discern chemical and tectonic processes involved in serpentinization.

Isotopic record across the Aptian/Albian boundary of ODP Hole 171-1049C

Geochemical analyses of extraordinarily well preserved late Aptian-early Albian foraminifera from Blake Nose (Ocean Drilling Program Site 1049) reveal rapid shifts of d18O, d13C, and 87Sr/88Sr in the subtropical North Atlantic that may be linked to a major planktic foraminifer extinction event across the Aptian/Albian boundary. The abruptness of the observed geochemical shifts and their coincidence with a sharp lithologic contact is explained as an artifact of a previously undetected hiatus of 0.8-1.4 million years at the boundary contact, but the values before and after the hiatus indicate that major oceanographic changes occurred at this time. 87Sr/88Sr increase by ~0.000200, d13C values decrease by 1.5 per mil to 2.2 per mil, and d18O values decrease by ~1.0 per mil (planktics) to 0.5 per mil (benthics) across the hiatus. Further, both 87Sr/88Sr ratios and d18O values during the Albian are anomalously high. The 87Sr/88Sr values deviate from known patterns to such a degree that an explanation requires either the presence of inter-basin differences in seawater 87Sr/88Sr during the Albian or revision of the seawater curve. For d18O, planktic values in some Aptian samples likely reflect a diagenetic overprint, but preservation is excellent in the rest of the section. In well preserved material, benthic foraminiferal values are largely between 0.5 and 0.0 per mil and planktic samples are largely between 0.0 per mil to -1.0 per mil, with a brief excursion to -2.0 per mil during OAE 1b. Using standard assumptions for Cretaceous isotopic paleotemperature calculations, the d18O values suggest bottom water temperatures (at ~1000 -1500 m) of 8-10°C and surface temperatures of 10-14°C, which are 4-6°C and 10-16°C cooler, respectively, than present-day conditions at the same latitude. The cool subtropical sea surface temperature estimates are especially problematic because other paleoclimate proxy data for the mid-Cretaceous and climate model predictions suggest that subtropical sea surface temperatures should have been the same as or warmer than at present. Because of their exquisite preservation, whole scale alteration of the analyzed foraminifera is an untenable explanation. Our proposed solution is a high evaporative fractionation factor in the early Albian North Atlantic that resulted in surface waters with higher d18O values at elevated salinities than commonly cited in Cretaceous studies. A high fractionation factor is consistent with high rates of vapor export and a vigorous hydrological cycle and, like the Sr isotopes, implies limited connectivity among the individual basins of the Early Cretaceous proto-Atlantic ocean.

Gas content, composition of gas hydrate and volumetric gas/water ratios of ODP and DSDP sites

Gas hydrate samples were recovered from four sites (Sites 994, 995, 996, and 997) along the crest of the Blake Ridge during Ocean Drilling Program (ODP) Leg 164. At Site 996, an area of active gas venting, pockmarks, and chemosynthetic communities, vein-like gas hydrate was recovered from less than 1 meter below seafloor (mbsf) and intermittently through the maximum cored depth of 63 mbsf. In contrast, massive gas hydrate, probably fault filling and/or stratigraphically controlled, was recovered from depths of 260 mbsf at Site 994, and from 331 mbsf at Site 997. Downhole-logging data, along with geochemical and core temperature profiles, indicate that gas hydrate at Sites 994, 995, and 997 occurs from about 180 to 450 mbsf and is dispersed in sediment as 5- to 30-m-thick zones of up to about 15% bulk volume gas hydrate. Selected gas hydrate samples were placed in a sealed chamber and allowed to dissociate. Evolved gas to water volumetric ratios measured on seven samples from Site 996 ranged from 20 to 143 mL gas/mL water to 154 mL gas/mL water in one sample from Site 994, and to 139 mL gas/mL water in one sample from Site 997, which can be compared to the theoretical maximum gas to water ratio of 216. These ratios are minimum gas/water ratios for gas hydrate because of partial dissociation during core recovery and potential contamination with pore waters. Nonetheless, the maximum measured volumetric ratio indicates that at least 71% of the cages in this gas hydrate were filled with gas molecules. When corrections for pore-water contamination are made, these volumetric ratios range from 29 to 204, suggesting that cages in some natural gas hydrate are nearly filled. Methane comprises the bulk of the evolved gas from all sites (98.4%-99.9% methane and 0%-1.5% CO2). Site 996 hydrate contained little CO2 (0%-0.56%). Ethane concentrations differed significantly from Site 996, where they ranged from 720 to 1010 parts per million by volume (ppmv), to Sites 994 and 997, which contained much less ethane (up to 86 ppmv). Up to 19 ppmv propane and other higher homologues were noted; however, these gases are likely contaminants derived from sediment in some hydrate samples. CO2 concentrations are less in gas hydrate than in the surrounding sediment, likely an artifact of core depressurization, which released CO2 derived from dissolved organic carbon (DIC) into sediment. The isotopic composition of methane from gas hydrate ranges from d13C of -62.5 per mil to -70.7 per mil and dD of -175 per mil to -200 per mil and is identical to the isotopic composition of methane from surrounding sediment. Methane of this isotopic composition is mainly microbial in origin and likely produced by bacterial reduction of bicarbonate. The hydrocarbon gases here are likely the products of early microbial diagenesis. The isotopic composition of CO2 from gas hydrate ranges from d13C of -5.7 per mil to -6.9 per mil, about 15 per mil lighter than CO2 derived from nearby sediment.

Simulated transition to agropastoralism in the Indus valley 7500-3000 BC

The Indus Valley Civilization (IVC) was one of the first great civilizations in prehistory. This bronze age civilization flourished from the end of the fourth millennium BC. It disintegrated during the second millennium BC; despite much research effort, this decline is not well understood. Less research has been devoted to the emergence of the IVC, which shows continuous cultural precursors since at least the seventh millennium BC. To understand the decline, we believe it is necessary to investigate the rise of the IVC, i.e., the establishment of agriculture and livestock, dense populations and technological developments 7000-3000 BC. Although much archaeologically typed information is available, our capability to investigate the system is hindered by poorly resolved chronology, and by a lack of field work in the intermediate areas between the Indus valley and Mesopotamia. We thus employ a complementary numerical simulation to develop a consistent picture of technology, agropastoralism and population developments in the IVC domain. Results from this Global Land Use and technological Evolution Simulator show that there is (1) fair agreement between the simulated timing of the agricultural transition and radiocarbon dates from early agricultural sites, but the transition is simulated first in India then Pakistan; (2) an independent agropas- toralism developing on the Indian subcontinent; and (3) a positive relationship between archeological artifact richness and simulated population density which remains to be quantified.

Chronology and petrology in MIS4-2 of DSDP Site 94-609

Ice-rafting evidence for a '1500-year cycle' sparked considerable debate on millennial-scale climate change and the role of solar variability. Here, we reinterpret the last 70,000 years of the subpolar North Atlantic record, focusing on classic DSDP Site 609, in the context of newly available raw data, the latest radiocarbon calibration (Marine09) and ice core chronology (GICC05), and a wider range of statistical methodologies. A ~1500-year oscillation is primarily limited to the short glacial Stage 4, the age of which is derived solely from an ice flow model (ss09sea), subject to uncertainty, and offset most from the original chronology. Results from the most well-dated, younger interval suggest that the original 1500 ± 500 year cycle may actually be an admixture of the ~1000 and ~2000 cycles that are observed within the Holocene at multiple locations. In Holocene sections these variations are coherent with 14C and 10Be estimates of solar variability. Our new results suggest that the '1500-year cycle' may be a transient phenomenon whose origin could be due, for example, to ice sheet boundary conditions for the interval in which it is observed. We therefore question whether it is necessary to invoke such exotic explanations as heterodyne frequencies or combination tones to explain a phenomenon of such fleeting occurrence that is potentially an artifact of arithmetic averaging.

Geochemical investigation and determination of hiatuses on Voering Plateau sediment records

The sediments recovered on ODP Leg 104 have been reported to be characterized by hiatuses. The hiatuses were defined by biostratigraphy and were believed to be caused by erosion related to temporary changes of bottom current composition and velocity. They have been associated with major paleoenvironmental changes, reorganization of global deep water production, and increased bottom water flows. Because of the importance of hiatuses for ongoing research, we decided to closely investigate the sedimentation history for the most significant Pliocene and Miocene biostratigraphic hiatuses by sedimentologic and geochemical means. The sedimentologic studies include clay mineral distributions, grain size data, and organic carbon concentrations. The geochemical studies include determination of 87/86Sr ratios, 10Be and Ir concentrations. The results of the sedimentologic studies suggest either that paleoenvironmental changes associated with hiatuses are not represented in the preserved sediments, or that the hiatuses are an artifact of interpretation of the biostratigraphic data. Strontium isotopes indicate continuous sedimentation for the interval investigated at Site 642, an interpretation confirmed by the steady decline in 10Be. 87/86Sr ratios in the interval from above and below proposed hiatuses H 2.2/2.3 and H2.1/2.2 at Site 643 display stronger changes with depth than expected by steady sedimentation. Ir data for this same interval indicate reduced sedimentation rates. Combining both, sedimentologic and geochemical evidence, the proposed hiatuses could not be confirmed and may represent preservation artifacts.

(Table 1) 87Sr/86Sr ratios and ages of DSDP Site 144A samples

Glassy Turonian foraminifera preserved in clay-rich sediments from the western tropical Atlantic yield the warmest equivalent d18O sea-surface temperatures (SSTs) yet reported for the entire Cretaceous-Cenozoic. We estimate Turonian SSTs that were at least as warm as (conservative mean ~30 °C) to significantly warmer (warm mean ~33 °C) than those in the region today. However, if independent evidence for high middle Cretaceous pCO2 is reliable and resulted in greater isotopic fractionation between seawater and calcite because of lower sea-surface pH, our conservative and warm SST estimates would be even higher (32 and 36°C, respectively). Our new tropical SSTs help reconcile geologic data with the predictions of general circulation models that incorporate high Cretaceous pCO2 and lend support to the hypothesis of a Cretaceous greenhouse. Our data also strengthen the case for a Turonian age for the Cretaceous thermal maximum and highlight a 20-40 m.y. mismatch between peak Cretaceous-Cenozoic global warmth and peak inferred tectonic CO2 production. We infer that this mismatch is either an artifact of a hidden Turonian pulse in global ocean-crust cycling or real evidence of the influence of some other factor on atmospheric CO2 and/or SSTs. A hidden pulse in crust cycling would explain the timing of peak Cretaceous-Cenozoic sea level (also Turonian), but other factors are needed to explain high-frequency (~10-100 k.y.) instability in middle Cretaceous SSTs reported elsewhere.

Early Cenozoic Atlantic and Pacific benthic foraminiferal isotopes from 12 sites

Oxygen and carbon isotope records are important tools used to reconstruct past ocean and climate conditions, with those of benthic foraminifera providing information on the deep oceans. Reconstructions are complicated by interspecies isotopic offsets that result from microhabitat preferences (carbonate precipitation in isotopically distinct environments) and vital effects (species-specific metabolic variation in isotopic fractionation). We provide correction factors for early Cenozoic benthic foraminifera commonly used for isotopic measurements (Cibicidoides spp., Nuttallides truempyi, Oridorsalis spp., Stensioina beccariiformis, Hanzawaia ammophila, and Bulimina spp.), showing that most yield reliable isotopic proxies of environmental change. The statistical methods and larger data sets used in this study provide more robust correction factors than do previous studies. Interspecies isotopic offsets appear to have changed through the Cenozoic, either (1) as a result of evolutionary changes or (2) as an artifact of different statistical methods and data set sizes used to determine the offsets in different studies. Regardless of the reason, the assumption that isotopic offsets have remained constant through the Cenozoic has introduced an 1-2°C uncertainty into deep sea paleotemperature calculations. In addition, we compare multiple species isotopic data from a western North Atlantic section that includes the Paleocene-Eocene thermal maximum to determine the most reliable isotopic indicator for this event. We propose that Oridorsalis spp. was the most reliable deepwater isotopic recorder at this location because it was best able to withstand the harsh water conditions that existed at this time; it may be the best recorder at other locations and for other extreme events also.