Strontium and oxygen isotope data for DSDP Hole 44-390A and ODP Leg 171 holes

Differences in regional responses to climate fluctuations are well documented on short time scales (e.g., El Niño-Southern Oscillation), but with the exception of latitudinal temperature gradients, regional patterns are seldom considered in discussions of ancient greenhouse climates. Contrary to the expectation of global warming or global cooling implicit in most treatments of climate evolution over millions of years, this paper shows that the North Atlantic warmed by as much as 6°C (1.5% decrease in d18O values of planktic foraminifera) during the Maastrichtian global cooling interval. We suggest that warming was the result of the importation of heat from the South Atlantic. Decreasing North Atlantic d18O values are also associated with increasing gradients in planktic d13C values, suggesting increasing surface-water stratification and a correlated strengthening of the North Atlantic Polar Front. If correct, this conclusion predicts arctic cooling during the late Maastrichtian. Beyond implications for the Maastrichtian, these data demonstrate that climate does not behave as if there is a simple global thermostat, even on geologic time scales.

Stable isotope and Mg/Ca ratios of Paleogene benthic foraminifera in the western North Atlantic

Global cooling and the development of continental-scale Antarctic glaciation occurred in the late middle Eocene to early Oligocene (~38 to 28 million years ago), accompanied by deep-ocean reorganization attributed to gradual Antarctic Circumpolar Current (ACC) development. Our benthic foraminiferal stable isotope comparisons show that a large d13C offset developed between mid-depth (~600 meters) and deep (>1000 meters) western North Atlantic waters in the early Oligocene, indicating the development of intermediate-depth d13C and O2 minima closely linked in the modern ocean to northward incursion of Antarctic Intermediate Water. At the same time, the ocean's coldest waters became restricted to south of the ACC, probably forming a bottom-ocean layer, as in the modern ocean. We show that the modern four-layer ocean structure (surface, intermediate, deep, and bottom waters) developed during the early Oligocene as a consequence of the ACC.

Late Miocene carbon isotope and paleoproductivity record

We examine whether or not a relationship exists between the late Miocene carbon isotope shift (~7.6-6.6 Ma) and marine productivity at four sites from the Indian and Pacific Oceans (Ocean Drilling Program Sites 721, 1146, 1172, and 846). We use a multiproxy approach based on benthic foraminiferal accumulation rates, elemental ratios, and dissolution indices, and we compare these data to benthic foraminiferal d13C values measured on the same samples. Although some of these sites have been targeted previously in studies of either the late Miocene/early Pliocene "biogenic bloom" (Sites 721 and 846) or the late Miocene carbon isotope shift (Site 1172), our records are the first to establish paired proxy records of carbon isotopes and paleoproductivity allowing a direct assessment of a potential link. Our results indicate that at all sites, productivity increased sometime during the d13C shift; at three sites (721, 1146, and 846), productivity increased at the beginning of the shift. The correlation coefficients derived from linear regression between micropaleontologically derived productivity and foraminiferal d13C values are relatively high during the time interval containing the late Miocene d13C shift (and statistically significant at three of the sites). Carbon flux and isotope mass balance considerations illustrate that transfer of organic matter between the terrestrial and marine reservoirs together with enhanced oceanic upwelling best approximates observed changes in carbon isotope records and paleoproductivity. We note that long-term trend in the Site 846 paleoproductivity record can be correlated to the long-term trend in the Site 848 eolian flux reconstructions of Hovan (1995, doi:10.2973/odp.proc.sr.138.132.1995) hinting at a link between strengthened wind regime and productivity during the late Miocene.

Stable isotope record and age determination of ODP Site 121-758 in the northeastern Indian Ocean

Numerous studies have shown that delta18O records from benthic and planktonic foraminifera, primarily a proxy of global ice volume variations, reflect Milankovitch periodicities. To study climatic response to orbital forcing at Ocean Drilling Program site 758, we have generated continuous delta18O and delta13C records from a single benthic foraminiferal species Cibicides wuellerstorfi for the last 3.6 m.y. and extended the planktonic foraminiferal isotope records of Farrell and Janecek (1991, doi:10.2973/odp.proc.sr.121.124.1991) (0-2.5 Ma, based on Globigerinoides sacculifer) to 3.6 Ma (Chen, 1994). We then constructed an age model by matching, correlating and tuning the benthic delta18O record to a model simulation of ice volume (Imbrie and Imbrie, 1980, doi:10.1126/science.207.4434.943). The filtered 41- and 23-kyr signals based on the resultant astronomically tuned age model are highly correlated to obliquity (r=0.83) and precession (r=0.75), respectively. Although derived with methodology different from Shackleton et al. (1990) and Hilgen (1991, doi:10.1016/0012-821X(91)90206-W, 1991, doi:10.1016/0012-821X(91)90082-S), our results generally agree with their published astronomical timescales for the time interval from 0 to 3.0 Ma, providing additional support for the newly emerging chronology based on orbital tuning. Slight discrepancies exist in the time interval from 3.0 to 3.6 Ma, suggesting several possibilities, including differences in the approaches of orbital tuning and the relatively low amplitude of delta18O variations in our record. However, even if the discrepancies are due to the relatively low amplitude of the isotope signals in our record at 3.0-3.6 Ma, our resultant timescale as a whole does not adversely affect our evaluation of the paleoclimatology and paleoceanography of the Indian Ocean, such as the evolution of the 100-, 41- and 23-kyr cycles, and variation of global ice volume and deepwater temperature during the past 3.6 m.y.

Carbon and oxygen isotopic composition of benthic foraminifers from ODP Leg 167 sites

Three sites, drilled during Ocean Drilling Program (ODP) Leg 167, were chosen for detailed late Pleistocene paleoceanographic studies of intermediate water along the California margin. These sites are Site 1011 (Animal Basin, 31°17'N, 117°38'W, 2033 m water depth, 1600 m sill depth), Site 1012 (East Cortez Basin, 32°17'N, 118°23'W, 1783 m water depth, 1415 m sill depth), and Site 1018 (Guide Seamount, 36°59'N, 123°17'W, 2476 m water depth). Here we present carbon and oxygen isotopic measurements of benthic foraminifers from these three sites. We made 135 measurements from Site 1011, 387 measurements from Site 1012, and 231 measurements from Site 1018. This data report includes an explanation of the methods used to generate these isotopic records and the age models for each site. Detailed paleoceanographic interpretations of the isotopic records are currently under way.

Stable isotope record of foraminifera from the subarctic Pacific

Benthic (Uvigerina spp., Cibicidoides spp., Gyroidinoides spp.) and planktonic (N. pachyderma sinistral, G. bulloides) stable isotope records from three core sites in the central Gulf of Alaska are used to infer mixed-layer and deepwater properties of the late glacial Subarctic Pacific. Glacial-interglacial amplitudes of the planktonic delta18O records are 1.1-1.3 per mil, less than half the amplitude observed at core sites at similar latitudes in the North Atlantic; these data imply that a strong, negative deltaw anomaly existed in the glacial Subarctic mixed layer during the summer, which points to a much stronger low-salinity anomaly than exists today. If true, the upper water column in the North Pacific would have been statically more stable than today, thus suppressing convection even more efficiently. This scenario is further supported by vertical (i.e., planktic versus benthic) delta18O and delta13C gradients of >1 per mil, which suggest that a thermohaline link between Pacific deep waters and the Subarctic Pacific mixed layer did not exist during the late glacial. Epibenthic delta13C in the Subarctic Pacific is more negative than at tropical-subtropical Pacific sites but similar to that recorded at Southern Ocean sites, suggesting ventilation of the deep central Pacific from mid-latitude sources, e.g., from the Sea of Japan and Sea of Okhotsk. Still, convection to intermediate depths could have occurred in the Subarctic during the winter months when heat loss to the atmosphere, sea ice formation, and wind-driven upwelling of saline deep waters would have been most intense. This would be beyond the grasp of our planktonic records which only document mixed-layer temperature-salinity fields extant during the warmer seasons. Also we do not have benthic isotope records from true intermediate water depths of the Subarctic Pacific.

Oxygen and carbon isotopic composition in Eocene and Oligocene planktonic and benthic foraminifera from DSDP sediments

Oxygen and carbon isotope ratios in Eocene and Oligocene planktonic and benthic foraminifera have been investigated from Atlantic, Indian, and Pacific Ocean locations. The major changes in Eocene-Oligocene benthic foraminiferal oxygen isotopes were enrichment of up to 1 per mil in 18O associated with the middle/late Eocene boundary and the Eocene/Oligocene boundary at locations which range from 1- to 4-km paleodepth. Although the synchronous Eocene-Oligocene 18O enrichment began in the latest Eocene, most of the change occurred in the earliest Oligocene. The earliest Oligocene enrichment in 18O is always larger in benthic foraminifera than in surface-dwelling planktonic foraminifera, a condition that indicates a combination of deep-water cooling and increased ice volume. Planktonic foraminiferal d18O does not increase across the middle/late Eocene boundary at our one site with the most complete record (Deep Sea Drilling Project Site 363, Walvis Ridge). This pattern suggests that benthic foraminiferal d18O increased 40 m.y. ago because of increased density of deep waters, probably as a result of cooling, although glaciation cannot be ruled out without more data. Stable isotope data are averaged for late Eocene and earliest Oligocene time intervals to evaluate paleoceanographic change. Average d18O of benthic foraminifera increased by 0.64 per mil from the late Eocene to the early Oligocene d18O maximum, whereas the average increase for planktonic foraminifera was 0.52 per mil. This similarity suggests that the Eocene/Oligocene boundary d18O increase was caused primarily by increased continental glaciation, coupled with deep sea cooling by as much as 2°C at some sites. Average d18O of surface-dwelling planktonic foraminifera from 14 upper Eocene and 17 lower Oligocene locations, when plotted versus paleo-latitude, reveals no change in the latitudinal d18O gradient. The Oligocene data are offset by ~0.45 per mil, also believed to reflect increased continental glaciation. At present, there are too few deep sea sequences from high latitude locations to resolve an increase in the oceanic temperature gradient from Eocene to Oligocene time using oxygen isotopes.

Geochemistry of middle Eocene-lower Oligocene pelagic sediments of ODP Hole 105-647A

Geochemical analyses of the middle Eocene through lower Oligocene lithologic Unit IIIC (260-518 meters below seafloor [mbsf]) indicate a relatively constant geochemical composition of the detrital fraction throughout this depositional interval at Ocean Drilling Program (ODP) Site 647 in the southern Labrador Sea. The main variability occurs in redox-sensitive elements (e.g., iron, manganese, and phosphorus), which may be related to early diagenetic mobility in anaerobic pore waters during bacterial decomposition of organic matter. Initial preservation of organic matter was mediated by high sedimentation rates (36 m/m.y.). High iron (Fe) and manganese (Mn) contents are associated with carbonate concretions of siderite, manganosiderite, and rhodochrosite. These concretions probably formed in response to elevated pore-water alkalinity and total dissolved carbon dioxide (CO2) concentrations resulting from bacterial sulfate reduction, as indicated by nodule stable-isotope compositions and pore-water geochemistry. These nodules differ from those found in upper Cenozoic hemipelagic sequences in that they are not associated with methanogenesis. Phosphate minerals (carbonate-fluorapatite) precipitated in some intervals, probably as the result of desorption of phosphorus from iron and manganese during reduction. The bulk chemical composition of the sediments differs little from that of North Atlantic Quaternary abyssal red clays, but may contain a minor hydrothermal component. The silicon/ aluminum (Si/Al) ratio, however, is high and variable and probably reflects original variations in biogenic opal, much of which is now altered to smectite and/or opal CT. An increase in the sodium/potassium (Na/K) ratio in the upper Eocene corresponds to the beginning of coarsergrained feldspar flux to the site, possibly marking the onset of more vigorous deep currents. Although the Site 647 cores provide a nearly complete high-resolution, high-latitude Eocene-Oligocene record, the high sedimentation rate and somewhat unusual diagenetic conditions have led to variable alteration of benthic foraminifers and fine-fraction carbonate and have overprinted the original stable-isotope records. Planktonic foraminifers are less altered, but on the whole, there is little chance of sorting out the nature and timing of environmental change on the basis of our stable-isotope analyses.

Benthic foraminiferal genera in DSDP Hole 95-612

On the basis of lithologic, foraminiferal, seismostratigraphic, and downhole logging characteristics, we identified seven distinctive erosional unconformities at the contacts of the principal depositional sequences at Site 612 on the New Jersey Continental Slope (water depth 1404 m). These unconformities are present at the Campanian/Maestrichtian, lower Eocene/middle Eocene, middle Eocene/upper Eocene, upper Eocene/lower Oligocene, lower Oligocene/upper Miocene, Tortonian/Messinian, and upper Pliocene/upper Pleistocene contacts. The presence of coarse sand or redeposited intraclasts above six of the unconformities suggests downslope transport from the adjacent shelf by means of sediment gravity flows, which contributed in part to the erosion. Changes in the benthic foraminiferal assemblages across all but the Campanian/Maestrichtian contact indicate that significant changes in the seafloor environment, such as temperature and dissolved oxygen content, took place during the hiatuses. Comparison with modern analogous assemblages and application of a paleoslope model where possible, indicate that deposition took place in bathyal depths throughout the Late Cretaceous and Cenozoic at Site 612. An analysis of two-dimensional geometry and seismic fades changes of depositional sequences along U.S.G.S. multichannel seismic Line 25 suggests that Site 612 was an outer continental shelf location from the Campanian until the middle Eocene, when the shelf edge retreated 130 km landward, and Site 612 became a continental slope site. Following this, a prograding prism of terrigenous debris moved the shelf edge to near its present position by the end of the Miocene. Each unconformity identified can be traced widely on seismic reflection profiles and most have been identified from wells and outcrops on the coastal plain and other offshore basins of the U.S. Atlantic margin. Furthermore, their stratigraphic positions and equivalence to similar unconformities on the Goban Spur, in West Africa, New Zealand, Australia, and the Western Interior of the U.S. suggest that most contacts are correlative with the global unconformities and sea-level falls of the Vail depositional model.

Stable carbon and oxygen isotope record of planktonic and benthic foraminifera at ODP Site 111-677 in the Panama Basin

Oxygen and carbon isotope ratio measurements are presented for Globigerinoides ruber and for benthic species (mainly Uvigerina spp.) in the Pleistocene and uppermost Pliocene section of ODP Hole 677A in the Panama Basin. This provides the best available continuous Pleistocene stable isotope records from any location, fully justifying the recoring of DSDP Site 504. Oxygen isotope stage 22 (age about 0.85 Ma) was of similar magnitude to the most extensive glacials of the Brunhes and constitutes a logical base for the middle Pleistocene. Oxygen isotope stages as defined by Ruddiman et al. (1986, doi:10.1016/0012-821X(86)90024-5) and by Raymo et al. (1989, doi:10.1029/PA004i004p00413) back to stage 104 are recognized. Although the internationally agreed base of the Quaternary at or near stage 62 (about 1.6 Ma) is not marked by a major isotopic event, it does approximate the base of a regime characterized by highly regular 41,000-yr climate cycles. The records at Site 677 are ideal for time-series analyses and will permit a new attempt to develop a chronology for the early Pleistocene based on tuning to the orbital frequencies. The carbon isotope records also appear to contain considerable variance at orbital frequencies throughout the sequence analyzed.

Benthic foraminiferal stable isotope stratigraphy of ODP Site 138-846 in the tropical Pacific Ocean

A stable-isotope stratigraphy at Site 846 (tropical Pacific, 3°06'S, 90°49'W, 3307 m water depth), based on the benthic foraminifers Cibicides wuellerstorfi and Uvigerina peregrina, yields a high-resolution record of deep-sea delta18O and delta13C over the past 1.8 Ma, with an average sampling interval of 3 k.y. Variance in the delta18O and delta13C records is concentrated in the well-known orbital periods of 100, 41, and 23 k.y. In the 100-k.y. band, both isotopic signals grow from relatively low amplitudes prior to 1.2 Ma, to high amplitudes in the late Quaternary since 0.7 Ma. The amplitude of delta18O and especially of delta13C decreases in the 41-k.y. band as it grows in the 100-k.y. band, consistent with a transfer of energy into an orbitally-paced internal oscillation. A weak 30-k.y. rhythm, present in both delta18O and delta13C, may reflect nonlinear interaction between the 41-k.y. and 100-k.y. bands in the evolving climate system. In the 23-k.y. and 19-k.y. bands associated with orbital precession, delta18O and delta13C are not coherent with each other on long time scales, and do not evolve like the 100-k.y. and 41-k.y. bands. This suggests that the source of the growing 100-k.y. oscillation is not a nonlinear response to precession, in contrast to predictions of some climate models. Sedimentation rates at this site also vary with a strong 100-k.y. cycle. Unlike the isotope records, the amplitude of 100-k.y. variations in sedimentation rate is relatively constant over the past 1.8 Ma, ranging from about 15 to 70 m/m.y. Prior to 0.9 Ma, sedimentation rates co-vary with orbital eccentricity, rather than with global climate as reflected by delta18O or delta13C. A source of this 100-k.y. cycle of sedimentation rate in the absence of similar ice volume fluctuations may be precessional heating of equatorial land masses, which in an energy balance climate model drives variations of monsoonal climates with a 100-k.y. rhythm. For the interval younger than 0.9 Ma, high sedimentation rates in the 100-k.y. band are consistently associated with glacial stages. This change of pattern suggests that when the amplitude of glacial cycles become large enough, their global effects overpower a local monsoon-driven variation in sedimentation rate at Site 846.