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.

A sedimentological, faunal, and isotopic record of the middle-to-late Pliocene transition at DSDP Hole 80-548

Hydraulic piston coring at DSDP Site 548, on the upper continental slope southwest of Ireland, recovered a nearly complete Pliocene section spanning 103 m of sediment. The sediments are greenish gray carbonate-rich hemipelagites containing abundant nannofossils and foraminifers. Grain-size analysis demonstrates that the texture of the section is fairly constant, with most of the variation occurring in 63- to 32-µm and < 2-µm fractions. Previous research has shown that the middle-to-late Pliocene transition in the North Atlantic was marked by the appearance of the planktonic foraminiferal species Globorotalia inflata and by the first occurrence of significant quantities of ice-rafted sediment grains in deep-sea sediments. The latter is taken to represent the first important development of Northern Hemisphere glaciation. The first appearance of G. inflata is carefully documented for Site 548 and is demonstrated to be an evolutionary datum at this site, rather than an ecologically controlled first appearance. Surface ocean conditions represented in the sediment section spanning the appearance of G. inflata were strongly cyclic, resulting in large periodic changes in the abundances of Globorotalia puncticulata and N. acostaensis. The benthic foraminiferal population was studied in detail over the middle-to-upper Pliocene transition to establish the nature and behavior of the intermediate-depth water mass in the northeastern Atlantic at the time of ice-sheet growth in the Northern Hemisphere. This water mass is presently warm and saline, having its source in the Mediterranean Sea. The benthic data show that the intermediate-depth water mass was undergoing a series of progressive changes over the interval including the first appearance of G. inflata. These changes are particularly reflected in the relative abundances of Globocassidulina subglobosa (Brady), Uvigerina, and Ehrenbergina. Also, the mean size of individuals in the G. subglobosa populations shows systematic variation, indicating changing intermediate-depth water properties. Oxygen-isotope analyses show that the intermediate-depth water mass was cold during the middle-to-late Pliocene transition. This interpretation is supported by the relative abundances of benthic foraminiferal species. Hence, the intermediate-depth northeastern Atlantic water mass of the middle to late Pliocene was considerably different from the intermediate-depth water mass of the present.

Middle to late Miocene oxygen isotope stratigraphy of ODP site 1085 (SE Atlantic): new constrains on Miocene climate variability and sea-level fluctuations

The middle Miocene delta18O increase represents a fundamental change in earth's climate system due to a major expansion and permanent establishment of the East Antarctic Ice Sheet accompanied by some effect of deepwater cooling. The long-term cooling trend in the middle to late Miocene was superimposed by several punctuated periods of glaciations (Mi-Events) characterized by oxygen isotopic shifts that have been related to the waxing and waning of the Antarctic ice-sheet and bottom water cooling. Here, we present a high-resolution benthic stable oxygen isotope record from ODP Site 1085 located at the southwestern African continental margin that provides a detailed chronology for the middle to late Miocene (13.9-7.3 Ma) climate transition in the eastern South Atlantic. A composite Fe intensity record obtained by XRF core scanning ODP Sites 1085 and 1087 was used to construct an astronomically calibrated chronology based on orbital tuning. The oxygen isotope data exhibit four distinct delta18O excursions, which have astronomical ages of 13.8, 13.2, 11.7, and 10.4 Ma and correspond to the Mi3, Mi4, Mi5, and Mi6 events. A global climate record was extracted from the oxygen isotopic composition. Both long- and short-term variabilities in the climate record are discussed in terms of sea-level and deep-water temperature changes. The oxygen isotope data support a causal link between sequence boundaries traced from the shelf and glacioeustatic changes due to ice-sheet growth. Spectral analysis of the benthic delta18O record shows strong power in the 400-kyr and 100-kyr bands documenting a paleoceanographic response to eccentricity-modulated variations in precession. A spectral peak around 180-kyr might be related to the asymmetry of the obliquity cycle indicating that the response of the dominantly unipolar Antarctic ice-sheet to obliquityinduced variations probably controlled the middle to late Miocene climate system. Maxima in the delta18O record, interpreted as glacial periods, correspond to minima in 100-kyr eccentricity cycle and minima in the 174-kyr obliquity modulation. Strong middle to late Miocene glacial events are associated with 400-kyr eccentricity minima and obliquity modulation minima. Thus, fluctuations in the amplitude of obliquity and eccentricity seem to be the driving force for the middle to late Miocene climate variability.

Analysis of rocks from Shackleton Range

Studies were made of the glacial geology and provenance of erratic in the Shackleton Range during the German geological expedition GEISHA in 1987/88, especially in the southern and northwestern parts of the range. Evidence that the entire Shackleton Range was once overrun by ice from a southerly to southeasterly direction was provided by subglacial erosional forms (e.g. striations, crescentic gouges, roches moutonnées) and erratics which probably orriginated in the region of the Whichaway Nunataks and the Pensacola Mountains in the southern part of the range. This probably happened during the last major expansion of the Anarctic polar ice sheet, which, on the basis of evidence from other parts of the continent, occurred towards the end of the Miocene. Till and an area of scattered erratics were mapped in the northwestern part of the range. These were deposited during a period of expansion of the Slessor Glacier in the Weichselian (Wisconsian) glacial stage earlier. This expansion was caused by blockage of the glacier by an expanded Filchner ice shelf which resulted from the sinking of the sea level during the Pleistocene, as demonstrated by geological studies in the Weddell Sea and along the coast of the Ross Sea. Studies of the erratics at the edges of glaciers provided information about rock concealed by the glacier.

Biomarker and x-ray diffraction data from Site U1313 (MIS 16 - 9)

A reconstruction of Milankovitch to millennial-scale variability of sea-surface temperature (SST) and sea-surface productivity in the Pleistocene mid-latitude North Atlantic Ocean (MIS 16-9) and its relationship to ice sheet instability was carried out on sediments from IODP Site U1313. This reconstruction is based on alkenone and n-alkane concentrations, Uk37' index, total organic carbon (TOC) and carbonate contents, X-Ray diffraction (XRD) data, magnetic susceptibility, and accumulation rates. Increased input of ice-rafted debris (IRD) occurred during MIS 16, 12, and 10, characterized by high concentrations of dolomite, quartz, and feldspars and elevated accumulation rates of terrigenous matter. Minimum input values of terrigenous matter, on the other hand, were determined for MIS 13 and 11. Peak values of dolomite, coinciding with quartz, plagioclase, and kalifeldspar peaks and maxima in long-chain n-alkanes indicative for land plants, are interpreted as Heinrich-like Events related to sudden instability of the Laurentide Ice Sheet during early and late (deglacial) phases of the glacials. The coincidence of increased TOC values with elevated absolute concentrations of alkenones suggest increased glacial productivity, probably due to a more southern position of the Polar Front. Alkenone-based SST reached absolute maxima of about 19°C during MIS 11.3 and absolute minima of <10°C during MIS 12 and 10. Within MIS 11, prominent cooling events (MIS 11.22 and 11.24) occurred. The absolute SST minima recorded directly before and after the glacial maxima MIS 10.2 and 12.2, are related to Heinrich-like Event meltwater pulses, as supported by the coincidence of SST minima and maxima in C37:4 alkenones and dolomite. These sudden meltwater pulses - especially during Terminations IV and V - probably caused a collapse of phytoplankton productivity as indicated by the distinct drop in alkenone concentrations. Ice-sheet disintegration and subsequent surges and outbursts of icebergs and meltwater discharge may have been triggered by increased insolation in the Northern High Latitudes.

Late Weichselian fluvial evolution of the southern Kara Sea

Glaciations had a profound impact on the global sea-level and particularly on the Arctic environments. One of the key questions related to this topic is, how did the discharge of the Siberian Ob and Yenisei rivers interact with a proximal ice sheet? In order to answer this question high-resolution (1-12 kHz), shallow-penetration seismic profiles were collected on the passive continental margin of the Kara Sea Shelf to study the paleo-drainage pattern of the Ob and Yenisei rivers. Both rivers incised into the recent shelf, leaving filled and unfilled river channels and river canyons/valleys connecting to a complex paleo-drainage network. These channels have been subaerially formed during a regressive phase of the global sea-level during the Last Glacial Maximum. Beyond recent shelf depths of 120 m particle transport is manifested in submarine channel-levee complexes acting as conveyor for fluvial-derived fines. In the NE area, uniform draping sediments are observed. Major morphology determining factors are (1) sea-level fluctuations and (2) LGM ice sheet influence. Most individual channels show geometries typical for meandering rivers and appear to be an order of magnitude larger than recent channel profiles of gauge stations on land. The Yenisei paleo-channels have larger dimensions than the Ob examples and could be originated by additional water release during the melt of LGM Putoran ice masses. Asymmetrical submarine channel-levee complexes with channel depths of 60 m and more developed, in some places bordered by glacially dominated morphology, implying deflection by the LGM ice masses. A total of more than 12,000 km of acoustic profiles reveal no evidence for an ice-dammed lake of greater areal extent postulated by several workers. Furthermore, the existence of the channel-levee complexes is indicative of unhindered sediment flow to the north. Channels situated on the shelf above 120-m water depth exhibit no phases of ponding and or infill during sea-level lowstand. These findings denote the non-existence of an ice sheet on large areas of the Kara Sea shelf.

Age determination and geochemical composition of sediment core DA04-31P

There is limited knowledge pertaining to the history of the Greenland Ice Sheet (GIS) during the last glacial-interglacial transition as it retreated from the continental margins to an inland position. Here we use multiproxy data, including ice-rafted debris (IRD); planktonic isotopes; alkenone temperatures; and tephra geochemistry from the northern Labrador Sea, off southwest Greenland, to investigate the deglacial response of the GIS and evaluate its implications for the North Atlantic deglacial development. The results imply that the southern GIS retreated in three successive stages: (1) early deglaciation of the East Greenland margins, by tephra-rich IRD that embrace Heinrich Event 1; (2) progressive retreat during Allerød culminating in major meltwater releases (d18O depletion of 1.2 per mil) at the Allerød-Younger Dryas transition (12.8-13.0 kyr B.P.); and (3) a final stage of glacial recession during the early Holocene (~9-11 kyr B.P.). Rather than indicating local temperatures of ambient surface water, the alkenones likely were transported to the core site by the Irminger Current. We attribute the timing of GIS retreat to the incursion of warm intermediate waters along the base of grounded glaciers and below floating ice shelves on the continental margin. The results lend support to the view that GIS meltwater presented a forcing factor for the Younger Dryas cooling.

Age determination and model of sediment core RAPiD-15-4P

Greenland ice core records indicate that the last deglaciation (~7-21 ka) was punctuated by numerous abrupt climate reversals involving temperature changes of up to 5°C-10°C within decades. However, the cause behind many of these events is uncertain. A likely candidate may have been the input of deglacial meltwater, from the Laurentide ice sheet (LIS), to the high-latitude North Atlantic, which disrupted ocean circulation and triggered cooling. Yet the direct evidence of meltwater input for many of these events has so far remained undetected. In this study, we use the geochemistry (paired Mg/Ca-d18O) of planktonic foraminifera from a sediment core south of Iceland to reconstruct the input of freshwater to the northern North Atlantic during abrupt deglacial climate change. Our record can be placed on the same timescale as ice cores and therefore provides a direct comparison between the timing of freshwater input and climate variability. Meltwater events coincide with the onset of numerous cold intervals, including the Older Dryas (14.0 ka), two events during the Allerød (at ~13.1 and 13.6 ka), the Younger Dryas (12.9 ka), and the 8.2 ka event, supporting a causal link between these abrupt climate changes and meltwater input. During the Bølling-Allerød warm interval, we find that periods of warming are associated with an increased meltwater flux to the northern North Atlantic, which in turn induces abrupt cooling, a cessation in meltwater input, and eventual climate recovery. This implies that feedback between climate and meltwater input produced a highly variable climate. A comparison to published data sets suggests that this feedback likely included fluctuations in the southern margin of the LIS causing rerouting of LIS meltwater between southern and eastern drainage outlets, as proposed by Clark et al. (2001, doi:10.1126/science.1062517).

Topographic model of Greenland (elevation lines in PanMap and ASCII format)

In this paper, a new high-resolution elevation model of Greenland, including the ice sheet as well as the ice free regions, is presented. It is the first published full coverage model, computed with an average resolution of 2 km and providing an unprecedented degree of detail. The topography is modeled from a wide selection of data sources, including satellite radar altimetry from Geosat and ERS 1, airborne radar altimetry and airborne laser altimetry over the ice sheet, and photogrammetric and manual map scannings in the ice free region. The ice sheet model accuracy is evaluated by omitting airborne laser data from the analysis and treating them as ground truth observations. The mean accuracy of the ice sheet elevations is estimated to be 12-13 m, and it is found that on surfaces of a slope between 0.2° and 0.8°, corresponding to approximately 50% of the ice sheet, the model presents a 40% improvement over models based on satellite altimetry alone. On coastal bedrock, the model is compared with stereo triangulated reference points, and it is found that the model accuracy is of the order of 25-35 m in areas covered by stereo photogrammetry scannings and between 200 and 250 m elsewhere.