Millenial-Scale stable oxygen isotope record of ODP Site 172-1058

We use the oxygen isotopic composition of planktonic foraminifera Globigerinoides ruber (white) from Ocean Drilling Program Site 1058 in the subtropical northwestern Atlantic to construct a high-resolution (~800 year) climate record spanning the mid-Pleistocene climate transition (~410 ka to 1350 ka). We investigate whether or not millennial-scale instabilities in the proxy record are associated with the extent of continental glaciation. G. ruber d18O values display high-frequency fluctuations throughout the record, but the amplitude about mean glacial and interglacial d18O values increases at marine isotope stage (MIS) 22 (880 ka) and is highest during MIS 12. These observations support that millennial-scale climate instabilities are associated with ice sheet size. Time series analysis illustrates that these variations have significant concentration of spectral power centered on periods of ~10-12 ka and ~5 ka. The timing of these fluctuations agrees well, or coincides with, the periodicities of the second and fourth harmonics, respectively, of precessional forcing at the equator. An insolation-based origin of the millennial-scale instabilities would be independent of ice volume and explains the presence of these fluctuations before the mid-Pleistocene climate transition as well as during interglacial intervals (e.g., MIS 37 and 17). Because the amplitude of the millennial-scale variations increases during the mid-Pleistocene transition, feedback mechanisms associated with the growth of large, 100-ka-paced, polar ice sheets may be important amplifiers of regional surface water hydrographic changes.

Geochemistry on sediment core Co1014 from Rauer Group, Antarctica

The history of glacial advances and retreats of the East Antarctic ice sheet during the Holocene is not well-known, due to limited field evidence in both the marine and terrestrial realm. A 257-cm-long sediment core was recovered from a marine inlet in the Rauer Group, East Antarctica, 1.8 km in front of the present ice-sheet margin. Radiocarbon dating and lithological characteristics reveal that the core comprises a complete marine record since 4500 yr. A significant ice-sheet expansion beyond present ice margins therefore did not occur during this period.

Global compilation of last glacial maximum oxygen isotopic ratios for planktonic and benthic foraminifera

We use the fully coupled atmosphere-ocean three-dimensional model of intermediate complexity iLOVECLIM to simulate the climate and oxygen stable isotopic signal during the Last Glacial Maximum (LGM, 21 000 yr). By using a model that is able to explicitly simulate the sensor (d18O), results can be directly compared with data from climatic archives in the different realms. Our results indicate that iLOVECLIM reproduces well the main feature of the LGM climate in the atmospheric and oceanic components. The annual mean d18O in precipitation shows more depleted values in the northern and southern high latitudes during the LGM. The model reproduces very well the spatial gradient observed in ice core records over the Greenland ice-sheet. We observe a general pattern toward more enriched values for continental calcite d18O in the model at the LGM, in agreement with speleothem data. This can be explained by both a general atmospheric cooling in the tropical and subtropical regions and a reduction in precipitation as confirmed by reconstruction derived from pollens and plant macrofossils. Data-model comparison for sea surface temperature indicates that iLOVECLIM is capable to satisfyingly simulate the change in oceanic surface conditions between the LGM and present. Our data-model comparison for calcite d18O allows investigating the large discrepancies with respect to glacial temperatures recorded by different microfossil proxies in the North Atlantic region. The results argue for a trong mean annual cooling between the LGM and present (>6°C), supporting the foraminifera transfer function reconstruction but in disagreement with alkenones and dinocyst reconstructions. The data-model comparison also reveals that large positive calcite d18O anomaly in the Southern Ocean may be explained by an important cooling, although the driver of this pattern is unclear. We deduce a large positive d18Osw anomaly for the north Indian Ocean that contrasts with a large negative d18Osw anomaly in the China Sea between the LGM and present. This pattern may be linked to changes in the hydrological cycle over these regions. Our simulation of the deep ocean suggests that changes in d18Osw between the LGM and present are not spatially homogenous. This is supported by reconstructions derived from pore fluids in deep-sea sediments. The model underestimates the deep ocean cooling thus biasing the comparison with benthic calcite d18O data. Nonetheless, our data-model comparison support a heterogeneous cooling of few degrees (2-4°C) in the LGM Ocean.

Geochemistry, lithology, and bulk density of Cricket Flat paleosol samples

The high-resolution marine isotope climate record indicates pronounced global cooling during the Langhian (16-13.8 Ma), beginning with the warm middle Miocene climatic optimum and ending with significant Antarctic ice sheet expansion and the transition to "icehouse" conditions. Terrestrial paleoclimate data from this interval is sparse and sometimes conflicting. In particular, there are gaps in the terrestrial record in the Pacific Northwest during the late Langhian and early Serravallian between about 14.5 and 12.5 Ma. New terrestrial paleoclimate data from this time and region could reconcile these conflicting records. Paleosols are particularly useful for reconstructing paleoenvironment because the rate and style of pedogenesis is primarily a function of surface environmental conditions; however, complete and well-preserved paleosols are uncommon. Most soils form in erosive environments that are not preserved, or in environments such as floodplains that accumulate in small increments; the resulting cumulic soils are usually thin, weakly developed, and subject to diagenetic overprinting from subsequent soils. The paleosol at Cricket Flat in northeastern Oregon is an unusually complete and well-preserved paleosol from a middle Miocene volcanic sequence in the Powder River Volcanic Field. An olivine basalt flow buried the paleosol at approximately 13.8 ± 0.6 Ma, based on three 40Ar/39Ar dates on the basalt. We described the Cricket Flat paleosol and used its physical and chemical profile and micromorphology to assess pedogenesis. The Cricket Flat paleosol is an Ultisol-like paleosol, chemically consistent with a high degree of weathering. Temperature and rainfall proxies suggest that Cricket Flat received 1120 ± 180 mm precipitation y-1 and experienced a mean annual temperature of 14.5 ± 2.1 °C during the formation of the paleosol, significantly warmer and wetter than today. This suggests slower cooling after the middle Miocene climatic optimum than is seen in the existing paleosol record.

Age determination and stable isotope record of planktonic foraminifera from sediments of the Baffin Shelf

Stable isotopic values on planktonic foraminifera in a suite of cores from basins across the SE Baffin Shelf are used to extract a record of meltwater events during Termination I deglaciation. Resolution and Hatton basins lie on the SE Baffin Shelf at water depths > 500 m, seaward of major conduits for ice drainage from the eastern sector of the Laurentide Ice Sheet (LIS). Accelerator mass spectrometry 14C dates are used to constrain our chronology of events in ten cores. In Resolution Basin, three cores have 14C AMS dates on foraminifera of > 20 ka at their bases; whereas Hatton Basin cores terminate in sediments < 13 kyr. Sedimentation rates varied between 0.1 to 4.5 m/ka. Stable oxygen and carbon isotopic ratios were obtained on 146 samples of the planktonic foraminifera Neogloboquadrina pachyderma (Ehrenberg) sinistral, from seven of the ten cores. No evidence was found to indicate that test morphology or size affected delta18O. Between 7 and 13.5 ka the surface water on the shelf was on average 1 per mil lower than the open ocean signal. Significant temporal variations were found in both delta18O and delta13C. Evidence for significant low delta18O events occurred between 13 and 8 ka. The delta13C record from the planktonic foraminifera suggests a threefold division of events between 13 and 7 ka, with positive values between 10.8 and 13.0 ka, negative values between 9 and 10.8 ka, and positive values from 7 to 9 ka. The delta18O data suggest the presence of meltwater on the shelf some 3,000 years prior to the first late glacial dates on terrestrial deglaciation (at circa 10.4 ka). "Hudson Strait must be the real key to the importance of the calving process during deglaciation, because it is potentially the largest marine outlet for the Laurentide Ice Sheet and because it leads into the very center of the ice sheet.....the rates of calving through Hudson Strait during the period of initial ?18O rise unfortunately are unknown." W. F. Ruddiman (1987, p. 151)

Oxygen and Carbon stable isotope ratios for Site U1334

The Oligocene-Miocene transition (OMT) (~23 Ma) is interpreted as a transient global cooling event, associated with a large-scale Antarctic ice sheet expansion. Here we present a 2.23 Myr long high-resolution (~3 kyr) benthic foraminiferal oxygen and carbon isotope (d18O and d13C) record from Integrated Ocean Drilling Program Site U1334 (eastern equatorial Pacific Ocean), covering the interval from 21.91 to 24.14 Ma. To date, five other high-resolution benthic foraminiferal stable isotope stratigraphies across this time interval have been published, showing a ~1 per mil increase in benthic foraminiferal d18O across the OMT. However, these records are still few and spatially limited and no clear understanding exists of the global versus local imprints. We show that trends and the amplitudes of change are similar at Site U1334 as in other high-resolution stable isotope records, suggesting that these represent global deep water signals. We create a benthic foraminiferal stable isotope stack across the OMT by combining Site U1334 with records from ODP Sites 926, 929, 1090, 1264, and 1218 to best approximate the global signal. We find that isotopic gradients between sites indicate interbasinal and intrabasinal variabilities in deep water masses and, in particular, note an offset between the equatorial Atlantic and the equatorial Pacific, suggesting that a distinct temperature gradient was present during the OMT between these deep water masses at low latitudes. A convergence in the d18O values between infaunal and epifaunal species occurs between 22.8 and 23.2 Ma, associated with the maximum d18O excursion at the OMT, suggesting climatic changes associated with the OMT had an effect on interspecies offsets of benthic foraminifera. Our data indicate a maximum glacioeustatic sea level change of ~50 m across the OMT.

Benthic foraminifers in Mesozoic and Cenozoic sediments of the southwestern Atlantic

Benthic foraminiferal assemblages in Mesozoic and Cenozoic sediments were studied at Sites 511, 512, 513, and 514 drilled during Leg 71 in the southwestern Atlantic on the Maurice Ewing Bank and in the Argentine Basin. Benthic foraminifers in almost all stratigraphic subdivisions of Sites 511 and 512 reflect the gradual subsidence of the Falkland Plateau from shelf depths in the Barremian-Albian, when a semiclosed basin with restricted circulation of water masses and anaerobic conditions existed, to lower bathyal depths in the Late Cretaceous and Cenozoic, with an abrupt acceleration at the boundary of Lower and Upper Cretaceous. The composition, distribution, and preservation of Late Cretaceous assemblages of benthic foraminifers suggest considerable fluctuations of the foraminiferal lysocline and the CCD. This is evidenced by dissolution facies and foraminiferal assemblages in which agglutinated and resistant calcareous forms predominated during high stands of the CCD and by calcareous facies in which rich assemblages of calcareous species predominated during low stands. The highest position of the CCD on the Plateau (less than 1500-2000 m) was in the late Cenomanian, Turonian, and Coniacian. In the Santonian and Campanian the CCD was at depths below 1500-2000 meters. At the end of the Campanian the CCD shifted again to depths comparable with those of Cenomanian and Turonian time. In the latest Campanian and the Maestrichtian the CCD was low and nanno-foraminiferal oozes with a rich assemblage of benthic foraminifers accumulated. Foraminiferal assemblages at Sites 513 and 514 in the Argentine Basin also testify to oceanic subsidence from lower bathyal depths in the Oligocene to abyssal ones at present. This process was complicated by the influence of geographical migrations of the Polar Front caused by extensions of the ice sheet in the Antarctic after the opening of the Drake Passage during the Oligocene. In Mesozoic and Cenozoic deposits of the Falkland Plateau and the Argentine Basin seven assemblages of benthic foraminifers were distinguished by age: early-middle Albian, middle-late Albian, Late Cretaceous (including four groups), middle Eocene, late Eocene-early Miocene, middle-late Miocene, and Pliocene-Quaternary. The Albian assemblages contain many species common to the foraminiferal fauna of the Austral Biogeographical Province. The Late Cretaceous assemblage contains, along with Austral species, species common to foraminifers of North America, Western Europe, the Russian platform, and the south of the U.S.S.R. Deep-sea cosmopolitan species prevail in Cenozoic assemblages.

Distribution of foraminifera and other components in ODP Leg 188 sites

During Leg 188 of the Ocean Drilling Program (ODP), employing JOIDES Resolution, we drilled holes at three sites in the southern Indian Ocean in and near Prydz Bay, East Antarctica, between 28 January and 29 February 2000. The objectives of the voyage were to: - Core through sediments deposited when Antarctica underwent the transition from "greenhouse" to the modern "icehouse" state late in the Eocene or early in the Oligocene, at sites obtaining their sediment from the currently subglacial Gamburtsev Mountains that probably were the site of nucleation of the ice sheet (principally Site 1166); - Obtain a sediment record from times at which major changes in the ice sheet volume and characteristics took place as judged from oxygen isotope records, especially at ~23.7 Ma (Oligocene/Miocene boundary), 12-16 Ma (middle Miocene), and 2.7 Ma (late Pliocene) (mainly Site 1165); and - Sample through the upper Pliocene and Quaternary in an attempt to document fluctuations in the extent of the ice sheet over the continental shelf during the Quaternary (especially Site 1167). Paleogene foraminifer-bearing marine sections were not intersected, and thus discussion of marine sections is restricted to the Neogene. Foraminifers are not major contributors to Leg 188 chronostratigraphy but contribute to paleoenvironmental interpretation, to issues such as carbonate compensation depth (CCD) effects and source and history of sediment, and provide a basis for Sr and d18O studies. Chronostratigraphy for the various sections was compiled from diatoms, radiolarians, and paleomagnetism (Shipboard Scientific Party, 2001, doi:10.2973/odp.proc.ir.188.101.2001). Foraminifers were sporadic rather than continuous except in short intervals; however, the Neogene foraminifers from the region are very poorly known and the new records proved to be of significant value in paleoenvironmental interpretation. Only at Site 1167 did drilling intersect a section that yielded foraminifers virtually throughout. Other than for the very young section at each site, there is virtually no continuity of assemblages between sites and thus each section is treated here as separate and unrelated.

Sedimentology and Nd isotope ratios of sediment core PS72/410-1

The freshwater budget of the Arctic Ocean is a key component governing the deep water formation in the North Atlantic and the global climate system. We analyzed the isotopic composition of neodymium (epsilon-Nd) in authigenic phases of marine sediments on the Mendeleev Ridge in the western Arctic Ocean spanning an estimated time interval from present to about 75 ka BP. This continuous record was used to reconstruct the epsilon-Nd of the polar deep water (PDW) and changes in freshwater sources to the PDW through time. Three deviations in epsilon-Nd from a long term average of -10.2 were identified at estimated 46-51, 35-39 and 13-21 ka BP. The estimated 46-51 ka BP event can be traced to bursting of ice-dammed lakes accompanying the collapse of the Barents-Kara Ice Sheet, which would have released radiogenic Nd to the eastern Arctic Ocean. The cyclonic surface circulation in the eastern Arctic Ocean must have been stronger than at present for the event to be recorded on the Mendeleev Ridge. For the 35-39 and 13-21 ka BP events, it is likely that the Laurentide Ice Sheet (LIS) supplied the unradiogenic freshwater. The configuration of the anticyclonic circulation in the western Arctic was probably similar to today or expanded eastward. Our simple mass balance calculations suggest that large amounts of freshwater were released but due to significant deep water formation within the Arctic Ocean, the effect on the formation of NADW was probably minor.

Orbital-resolution Pliocene-Pleistocene simulations of CO2, sea level, global temperature and benthic d18O, and d11B-based proxy-CO2 data

During the past five million yrs, benthic d18O records indicate a large range of climates, from warmer than today during the Pliocene Warm Period to considerably colder during glacials. Antarctic ice cores have revealed Pleistocene glacial-interglacial CO2 variability of 60-100 ppm, while sea level fluctuations of typically 125 m are documented by proxy data. However, in the pre-ice core period, CO2 and sea level proxy data are scarce and there is disagreement between different proxies and different records of the same proxy. This hampers comprehensive understanding of the long-term relations between CO2, sea level and climate. Here, we drive a coupled climate-ice sheet model over the past five million years, inversely forced by a stacked benthic d18O record. We obtain continuous simulations of benthic d18O, sea level and CO2 that are mutually consistent. Our model shows CO2 concentrations of 300 to 470 ppm during the Early Pliocene. Furthermore, we simulate strong CO2 variability during the Pliocene and Early Pleistocene. These features are broadly supported by existing and new d11B-based proxy CO2 data, but less by alkenone-based records. The simulated concentrations and variations therein are larger than expected from global mean temperature changes. Our findings thus suggest a smaller Earth System Sensitivity than previously thought. This is explained by a more restricted role of land ice variability in the Pliocene. The largest uncertainty in our simulation arises from the mass balance formulation of East Antarctica, which governs the variability in sea level, but only modestly affects the modeled CO2 concentrations.

North Atlantic IODP Site 303-U1308 Middle Pleistocene Heinrich Events record

Heinrich events are well documented for the last glaciation, but little is known about their occurrence in older glacial periods of the Pleistocene. Here we report scanning XRF and bulk carbonate d18O results from Integrated Ocean Drilling Program Site U1308 (reoccupation of Deep Sea Drilling Project Site 609) that are used to develop proxy records of ice-rafted detritus (IRD) for the last ~1.4 Ma. Ca/Sr is used as an indicator of IRD layers that are rich in detrital carbonate (i.e., Heinrich layers), whereas Si/Sr reflects layers that are poor in biogenic carbonate and relatively rich in detrital silicate minerals. A pronounced change occurred in the composition and frequency of IRD at ~640 ka during marine isotope stage (MIS) 16, coinciding with the end of the middle Pleistocene transition. At this time, "Hudson Strait" Heinrich layers suddenly appeared in the sedimentary record of Site U1308, and the dominant period of the Si/Sr proxy shifted from 41 ka prior to 640 ka to 100 ka afterward. The onset of Heinrich layers during MIS 16 represents either the initiation of surging of the Laurentide Ice Sheet (LIS) off Hudson Strait or the first time icebergs produced by this process survived the transport to Site U1308. We speculate that ice volume (i.e., thickness) and duration surpassed a critical threshold during MIS 16 and activated the dynamical processes responsible for LIS instability in the region of Hudson Strait. We also observe a strong coupling between IRD proxies and benthic d13C variation at Site U1308 throughout the Pleistocene, supporting a link between iceberg discharge and weakening of thermohaline circulation in the North Atlantic.

Mg/Ca values for planktonic foraminifera N. pachyderma and G. bulloides from ocean sediment core MD02-2496 on the Vancouver margin, Northeastern Pacific Ocean

A new, high-resolution planktonic foraminiferal Mg/Ca-based ocean temperature record has been generated for deep sea core MD02-2496, sited offshore of Vancouver Island, Western Canada during the last deglaciation (21-12 ka). The relationship between Cordilleran Ice Sheet (CIS) retreat and changing regional ocean temperatures has been reconstructed through glaciomarine sediments in MD02-2496 that capture tidewater glacier response to surface ocean thermal forcing. At CIS maximum extent, the marine margin of the ice sheet advanced onto the continental shelf. During this interval, ocean temperatures recorded by surface ocean dwelling Globigerina bulloides remained a relatively constant ~7.5°C while subsurface dwelling Neogloboquadrina pachyderma (s.) recorded temperatures of ~5°C. These ocean temperatures were sufficiently warm to induce significant melt along the tidewater ice terminus similar to modern Alaskan tidewater glacial systems. During the deglacial retreat of the CIS, the N. pachyderma temperature record shows two distinct warming steps of ~2 and 2.5°C between 17.2-16 and 15.5-14 ka respectively, coincident with ice rafting events from the CIS, while G. bulloides records an ~3°C warming from 15 to14 ka. We hypothesize that submarine melting resulting from relatively warm ocean temperatures was an important process driving ice removal from CIS tidewater glaciers during the initial stages of deglaciation.