(Table 1) Two-way travel times and depths of picked internal reflection horizons (IRHs) of the NGRIP ice core

During the summer of 2003, a ground-penetrating radar survey around the North Greenland Icecore Project (NorthGRIP) deep ice-core drilling site (75°06' N, 42°20' W; 2957 m a.s.l.) was carried out using a shielded 250 MHz radar system. The drill site is located on an ice divide, roughly 300 km north-northwest of the summit of the Greenland ice sheet. More than 430 km of profiles were measured, covering a 10 km by 10 km area, with a grid centered on the drilling location, and eight profiles extending beyond this grid. Seven internal horizons within the upper 120 m of the ice sheet were continuously tracked, containing the last 400 years of accumulation history. Based on the age-depth and density-depth distribution of the deep core, the internal layers have been dated and the regional and temporal distribution of accumulation rate in the vicinity of NorthGRIP has been derived. The distribution of accumulation shows a relatively smoothly increasing trend from east to west from 145 kg/m**2/a to 200 kg/m**2/a over a distance of 50 km across the ice divide. The general trend is overlain by small-scale variations on the order of 2.5 kg/m**2/a/km, i.e. around 1.5% of the accumulation mean. The temporal variations of the seven periods defined by the seven tracked isochrones are on the order of +-4% of the mean of the last 400 years, i.e. at NorthGRIP ±7 kg/m**2/a. If the regional accumulation pattern has been stable for the last several thousand years during the Holocene, and ice flow has been comparable to today, advective effects along the particle trajectory upstream of NorthGRIP do not have a significant effect on the interpretation of climatically induced changes in accumulation rates derived from the deep ice core over the last 10 kyr.

Age determination of sediment cores from the northern Mid-Atlantic Ridge

We construct age models for a suite of cores from the northeast Atlantic Ocean by means of accelerator mass spectrometer dating of a key core, BOFS 5K, and correlation with the rest of the suite. The effects of bioturbation and foraminiferal species abundance gradients upon the age record are modeled using a simple equation. The degree of bioturbation is estimated by comparing modeled profiles with dispersal of the Vedde Ash layer in core 5K, and we find a mixing depth of roughly 8 cm for sand-sized material. Using this value, we estimate that age offsets between unbioturbated sediment and some foraminifera species after mixing may be up to 2500 years, with lesser effect on fine carbonate (< 10 µm) ages. The bioturbation model illustrates problems associated with the dating of 'instantaneous' events such as ash layers and the 'Heinrich' peaks of ice-rafted detritus. Correlations between core 5K and the other cores from the BOFS suite are made on the basis of similarities in the downcore profiles of oxygen and carbon isotopes, magnetic susceptibility, water and carbonate content, and via marker horizons in X radiographs and ash beds.

Geochemistry and age determination of zircons in Cretaceous to Quaternary sediments

In central Antarctica, drainage today and earlier back to the Paleozoic radiates from the Gamburtsev Subglacial Mountains (GSM). Proximal to the GSM past the Permian-Triassic fluvial sandstones in the Prince Charles Mountains (PCM) are Cretaceous, Eocene, and Pleistocene sediment in Prydz Bay (ODP741, 1166, and 1167) and pre-Holocene sediment in AM04 beneath the Amery Ice Shelf. We analysed detrital zircons for U-Pb ages, Hf-isotope compositions, and trace elements to determine the age, rock type, source of the host magma, and "crustal" model age (T(C)DM). These samples, together with others downslope from the GSM and the Vostok Subglacial Highlands (VSH), define major clusters of detrital zircons interpreted as coming from (1) 700 to 460 Ma mafic granitoids and alkaline rock, epsilon-Hf 9 to -28, signifying derivation 2.5 to 1.3 Ga from fertile and recycled crust, and (2) 1200-900 Ma mafic granitoids and alkaline rock, epsilon-Hf 11 to -28, signifying derivation 1.8 to 1.3 Ga from fertile and recycled crust. Minor clusters extend to 3350 Ma. Similar detrital zircons in Permian-Triassic, Ordovician, Cambrian, and Neoproterozoic sandstones located along the PaleoPacific margin of East Antarctica and southeast Australia further downslope from central Antarctica reflect the upslope GSM-VSH nucleus of the central Antarctic provenance as a complex of 1200-900 Ma (Grenville) mafic granitoids and alkaline rocks and older rocks embedded in 700-460 Ma (Pan-Gondwanaland) fold belts. The wider central Antarctic provenance (CAP) is tentatively divided into a central sector with negative ?Hf in its 1200-900 Ma rocks bounded on either side by positive epsilon-Hf. The high ground of the GSM-VSH in the Permian and later to the present day is attributed to crustal shortening by far-field stress during the 320 Ma mid-Carboniferous collision of Gondwanaland and Laurussia. Earlier uplifts in the ~500 Ma Cambrian possibly followed the 700-500 Ma assembly of Gondwanaland, and in the Neoproterozoic the 1000-900 Ma collisional events in the Eastern Ghats-Rayner Province at the end of the 1300-1000 Ma assembly of Rodinia.

A high-resolution record of the atmospheric CO2 concentration from 60-20 kyr BP from the Taylor Dome ice core, Antarctica

A high-resolution record of the atmospheric CO2 concentration from 60 to 20 thousand years before present (kyr BP) based on measurements on the ice core of Taylor Dome, Antarctica is presented. This record shows four distinct peaks of 20 parts per million by volume (ppmv) on a millennial time scale. Good correlation of the CO2 record with temperature reconstructions based on stable isotope measurements on the Vostok ice core (Antarctica) is found.

Age model and Plio-Pleistocene benthic stable oxygen isotope ratios of ODP Site 184-1143 in the Sourh China Sea

Based on benthic foraminiferal delta18O from ODP Site 1143, a 5-Myr astronomical timescale for the West Pacific Plio-Pleistocene was established using an automatic orbital tuning method. The tuned Brunhes/Matuyama paleomagnetic polarity reversal age agrees well with the previously published age of 0.78 Ma. The tuned ages for several planktonic foraminifer bio-events also agree well with published dates, and new ages for some other bio-events in the South China Sea were also estimated. The benthic delta18O from Site 1143 is highly coherent with the Earth's orbit (ETP) both at the obliquity and precession bands for the last 5 Myr, and at the eccentricity band for the last 2 Myr. In general, the 41-kyr cycle was dominant through the Plio-Pleistocene although the 23-kyr cycle was also very strong. The 100-kyr cycle became dominant only during the last 1 Myr. A comparison of the benthic delta18O between the Atlantic (ODP 659) and the East and West Pacific (846 and 1143) reveals that the Atlantic-Pacific benthic oxygen isotope difference ratio (Delta delta18OAtl-Pac) displays an increasing trend in three time intervals: 3.6-2.7 Ma, 2.7-2.1 Ma and 1.5-0.25 Ma. Each of the intervals begins with a rapid negative shift in Delta delta18OAtl-Pac, followed by a long period with an increasing trend, corresponding to the growth of the Northern Hemisphere ice sheet. This means that all three intervals of ice sheet growth in the Northern Hemisphere were accompanied at the beginning by a rapid relative warming of deep water in the Atlantic as compared to that of the Pacific, followed by its gradual relative cooling. This general trend, superimposed on the frequent fluctuations with glacial cycles, should yield insights into the processes leading to the boreal glaciation. Cross-spectral analyses of the Delta delta18OAtl-Pac with the Earth's orbit suggests that after the initiation of Northern Hemisphere glaciation at about 2.5 Ma, obliquity rather than precession had become the dominant force controlling the vertical structure or thermohaline circulation in the paleo-ocean.

Lagrangian back-tracing results for Kohnen station; output of time dependent numerical experiment with a 3D nested Antarctic ice sheet model

A nested ice flow model was developed for eastern Dronning Maud Land to assist with the dating and interpretation of the EDML deep ice core. The model consists of a high-resolution higher-order ice dynamic flow model that was nested into a comprehensive 3-D thermomechanical model of the whole Antarctic ice sheet. As the drill site is on a flank position the calculations specifically take into account the effects of horizontal advection as deeper ice in the core originated from higher inland. First the regional velocity field and ice sheet geometry is obtained from a forward experiment over the last 8 glacial cycles. The result is subsequently employed in a Lagrangian backtracing algorithm to provide particle paths back to their time and place of deposition. The procedure directly yields the depth-age distribution, surface conditions at particle origin, and a suite of relevant parameters such as initial annual layer thickness. This paper discusses the method and the main results of the experiment, including the ice core chronology, the non-climatic corrections needed to extract the climatic part of the signal, and the thinning function. The focus is on the upper 89% of the ice core (appr. 170 kyears) as the dating below that is increasingly less robust owing to the unknown value of the geothermal heat flux. It is found that the temperature biases resulting from variations of surface elevation are up to half of the magnitude of the climatic changes themselves.

Calcium and sodium concentrations of the Talos Dome ice core, Antarctica

In this study we report on new non-sea salt calcium (nssCa2+, mineral dust proxy) and sea salt sodium (ssNa+, sea ice proxy) records along the East Antarctic Talos Dome deep ice core in centennial resolution reaching back 150 thousand years (ka) before present. During glacial conditions nssCa2+ fluxes in Talos Dome are strongly related to temperature as has been observed before in other deep Antarctic ice core records, and has been associated with synchronous changes in the main source region (southern South America) during climate variations in the last glacial. However, during warmer climate conditions Talos Dome mineral dust input is clearly elevated compared to other records mainly due to the contribution of additional local dust sources in the Ross Sea area. Based on a simple transport model, we compare nssCa2+ fluxes of different East Antarctic ice cores. From this multi-site comparison we conclude that changes in transport efficiency or atmospheric lifetime of dust particles do have a minor effect compared to source strength changes on the large-scale concentration changes observed in Antarctic ice cores during climate variations of the past 150 ka. Our transport model applied on ice core data is further validated by climate model data. The availability of multiple East Antarctic nssCa2+ records also allows for a revision of a former estimate on the atmospheric CO2 sensitivity to reduced dust induced iron fertilisation in the Southern Ocean during the transition from the Last Glacial Maximum to the Holocene (T1). While a former estimate based on the EPICA Dome C (EDC) record only suggested 20 ppm, we find that reduced dust induced iron fertilisation in the Southern Ocean may be responsible for up to 40 ppm of the total atmospheric CO2 increase during T1. During the last interglacial, ssNa+ levels of EDC and EPICA Dronning Maud Land (EDML) are only half of the Holocene levels, in line with higher temperatures during that period, indicating much reduced sea ice extent in the Atlantic as well as the Indian Ocean sector of the Southern Ocean. In contrast, Holocene ssNa+ flux in Talos Dome is about the same as during the last interglacial, indicating that there was similar ice cover present in the Ross Sea area during MIS 5.5 as during the Holocene.

Ice-rafted debris of Late Cenozoic sediments from ODP Leg 104 holes

The abundance and composition of the upper Cenozoic terrigenous coarse-sand fraction (250 µm-2 mm) at ODP Sites 642, 643, and 644 were investigated to date the onset of significant ice-rafting in the Norwegian Sea, establish the regional chronology of ice-rafting, and determine the relative importance of global vs. regional controls on ice-rafting in this area. The first input of ice-rafted debris (IRD) occurs at approximately 2.9 Ma, with significant ice-rafting beginning at about 2.5 Ma. IRD abundances increase significantly in sediments younger than 0.9 Ma at all three holes, indicating climatic deterioration in the late Pleistocene. Differences in the timing of this IRD increase between holes result from regional patterns of IRD supply and surface circulation. Variations in IRD sources and dispersal patterns may also explain the slightly higher background level of IRD abundance at Hole 642B, a seaward site. Major peaks in the generalized IRD records from the Norwegian Sea are tentatively correlated to glacial stages or glacial-to-interglacial transitions in the globally defined oxygen isotope record. This correlation indicates the effect of global conditions on the regional climate of the Norwegian Sea, although the detailed IRD records at these sites are also affected by local/regional processes (e.g., circulation patterns and source area differences).

Ice-rafted debris in sites of ODP Leg 114

All holes drilled during Leg 114 contained ice-rafted debris. Analysis of samples from Hole 699A, Site 701, and Hole 704A yielded a nearly complete history of ice-rafting episodes. The first influx of ice-rafted debris at Site 699, on the northeastern slope of the Northeast Georgia Rise, occurred at a depth of 69.94 m below seafloor (mbsf) in sediments of early Miocene age (23.54 Ma). This material is of the same type as later ice-rafted debris, but represents only a small percentage of the coarse fraction. Significant ice-rafting episodes occurred during Chron 5. Minor amounts of ice-rafted debris first reached Site 701, on the western flank of the Mid-Atlantic Ridge (8.78 Ma at 200.92 mbsf), and more arrived in the late Miocene (5.88 Ma). The first significant quantity of sand and gravel appeared at a depth of 107.76 mbsf (4.42 Ma). Site 704, on the southern part of the Meteor Rise, received very little or no ice-rafted debris prior to 2.46 Ma. At this time, however, the greatest influx of ice-rafted debris occurred at this site. This time of maximum ice rafting correlates reasonably well with influxes of ice-rafted debris at Sites 701 (2.24 Ma) and 699 (2.38 Ma), in consideration of sample spacing at these two sites. These peaks of ice rafting may be Sirius till equivalents, if the proposed Pliocene age of Sirius tills can be confirmed. After about 1.67 Ma, the apparent mass-accumulation rate of the sediments at Site 704 declined, but with major fluctuations. This decline may be the result of a decrease in the rate of delivery of detritus from Antarctica due to reduced erosive power of the glaciers or a northward shift in the Polar Front Zone, a change in the path taken by the icebergs, or any combination of these factors.

First results from stable-istope measurements on ice cores from Dronning Maud Land

The European Project for Ice Coring in Antarctica (EPICA) focuses on the drilling of two deep ice cores, the first at Dome C and the second at Kohnen station (75°00' S, 0°04' E) in Dronning Maud Land (DML). This paper deals with stable-isotope records from ice cores drilled in DML. In the first season, the deep EPICA DML core reached a depth of 450 m, recovering ice approximately 7000 years old. Generally, the d18O record indicates a stable Holocene climate and shows low variability. However, during the last 4000 years (based on a preliminary time-scale) the d18O values decrease continuously by about 0.6%, and the deuterium excess values increase by about 0.5%. The correlation between d18O and the deuterium excess d is investigated for a 50m long core section and the near-surface snow. High-pass filtered profiles are positively correlated, whereas the correlation between low-pass filtered profiles is negative. A post-depositional effect due to diffusion processes can be seen in a sub-annually resolved profile from snow-pit samples. Changes in the seasonality of the evolution of the snow cover and the consequences for stable-isotope content are demonstrated with data from ice core B31.

Data from the EPICA Dome C ice core EDC

The Antarctic Vostok ice core provided compelling evidence of the nature of climate, and of climate feedbacks, over the past 420,000 years. Marine records suggest that the amplitude of climate variability was smaller before that time, but such records are often poorly resolved. Moreover, it is not possible to infer the abundance of greenhouse gases in the atmosphere from marine records. Here we report the recovery of a deep ice core from Dome C, Antarctica, that provides a climate record for the past 740,000 years. For the four most recent glacial cycles, the data agree well with the record from Vostok. The earlier period, between 740,000 and 430,000 years ago, was characterized by less pronounced warmth in interglacial periods in Antarctica, but a higher proportion of each cycle was spent in the warm mode. The transition from glacial to interglacial conditions about 430,000 years ago (Termination V) resembles the transition into the present interglacial period in terms of the magnitude of change in temperatures and greenhouse gases, but there are significant differences in the patterns of change. The interglacial stage following Termination V was exceptionally long - 28,000 years compared to, for example, the 12,000 years recorded so far in the present interglacial period. Given the similarities between this earlier warm period and today, our results may imply that without human intervention, a climate similar to the present one would extend well into the future.

Age model of ODP Hole 107-654A

It has long been speculated that glacio-eustatic sea level drop may have caused or contributed to the isolation and consequent desiccation of the Mediterranean in the late Miocene (the 'Messinian salinity crisis'). Ocean Drilling Program site 654 on the Sardinia margin of the Tyrrhenian Sea is the first deep-sea drill site to penetrate through upper Messinian evaporites into lower Messinian/upper Tortonian open marine sediments, and thus offers a unique opportunity to date the onset of the salinity crisis. A reexamination of the magnetostratigraphic, biostratigraphic, and stable isotope-stratigraphic constraints on the preevaporite sediments of site 654 has yielded two possible ages for the contact between salinity crisis sediments and the underlying normal marine sediments. One magnetostratigraphic interpretation plus the biostratigraphically determined position of the Tortonian/Messinian boundary imply a date of about 6.2 Ma for the youngest presalinity crisis sediments. An alternative magnetostratigraphic interpretation plus the carbon isotope stratigraphy imply a date of about 5.2 Ma. The younger of these dates coincides with a delta18O spike in open ocean sediments [Keigwin, 1987 doi:10.1029/PA002i006p00639], which is attributed to increased ice volume.