Glacial sedimentation in ODP Leg 119 samples

Although scientific evidence prior to that from ODP Leg 119 indicates the presence of an ice sheet on East Antarctica by at least the earliest Oligocene, the question as to the size and stability of that initial ice sheet is still contested. Current hypotheses include (1) the presence of a small ice sheet in the earliest Oligocene with stepwise growth during the Neogene, (2) the presence of a continental-sized ice sheet in the late middle Eocene with no major evidence of subsequent deglaciation, and (3) the presence of glacial ice in the earliest Oligocene with a major ice sheet during the mid-Oligocene, followed by growth and decay of several ice sheets with characteristics similar to the temperate ice sheets of the Pleistocene of North America but with changes over a longer time scale (millions of years vs. 100,000 yr). Principal results from Leg 119 suggest the presence of significant late middle and late Eocene glaciation in East Antarctica and the presence of a continental-size ice sheet in East Antarctica during the earliest Oligocene. Although the Leg 119 results provide only glimpses of the Neogene glacial history of East Antarctica, they do provide evidence of fluctuations in the extent of the ice sheet and the waxing and waning of glaciers across the Prydz Bay shelf during the later part of the late Miocene and Pliocene.

Geotechnical properties of shelf sdiments from the Prydz Bay

This paper presents a geotechnical characterization of the glacigenic sediments in Prydz Bay, East Antarctica, based on the shipboard physical properties data obtained during Leg 119, combined with results of land-based analyses of 24 whole-round core samples. Main emphasis is placed on the land-based studies, which included oedometer consolidation tests, triaxial and simple shear tests for undrained shear strength, permeability tests in oedometer and triaxial cell, Atterberg limits, and grain-size analyses. The bulk of the tested sediments comprise overconsolidated diamictites of a relatively uniform lithology. The overconsolidation results from a combination of glacial loading and sediment overburden subsequently removed by extensive glacial erosion of the shelf. This leads to downhole profiles of physical properties that have been observed not to change as a function of the thickness of present overburden. A number of fluctuations in the parameters shows a relatively systematic trend and most likely results from changes in the proximity to the ice sheet grounding line in response to variations in the glacial regime. Very low permeabilities mainly result from high preconsolidation stresses (Pc'). Pc' values up to 10,000 kPa were estimated from the oedometer tests, and empirical estimates based on undrained shear strengths (up to 2500 kPa) indicate that the oedometer results are conservative. The diamictites generally classify as inactive, of low to medium plasticity, and they consolidate with little deformation, even when subjected to great stresses. This is the first report of geotechnical data from deep boreholes on the Antarctic continental shelf, but material of similar character can also be expected in other areas around the Antarctic.

Mg/Ca ratios of Cenozoic benthic foraminifera from the Indian and South Atlantic Ocean

We investigate the evolution of Cenozoic climate and ice volume as evidenced by the oxygen isotopic composition of seawater (delta18Osw) derived from benthic foraminiferal Mg/Ca ratios to constrain the temperature effect contained in foraminiferal delta18O values. We have constructed two benthic foraminiferal Mg/Ca records from intermediate water depth sites (Ocean Drilling Program sites 757 and 689 from the subtropical Indian Ocean and the Weddell Sea, respectively). Together with the previously published composite record of Lear et al. (2002, doi:10.1126/science.287.5451.269) and the Neogene record from the Southern Ocean of Billups and Schrag (2002, doi:10.1029/2000PA000567), we obtain three, almost complete representations of the delta18Osw for the past 52 Myr. We discuss the sensitivity of early Cenozoic Mg/Ca-derived paleotemperatures (and hence the delta18Osw) to assumptions about seawater Mg/Ca ratios. We find that during the middle Eocene (~ 49-40 Ma), modern seawater ratios yield Mg/Ca-derived temperatures that are in good agreement with the oxygen isotope paleothermometer assuming ice-free conditions. Intermediate waters cooled during the middle Eocene reaching minimum temperatures by 40 Ma. The corresponding delta18Osw reconstructions support ice growth on Antarctica beginning by at least 40 Ma. At the Eocene/Oligocene boundary, Mg/Ca ratios (and hence temperatures) from Weddell Sea site 689 display a well-defined maximum. We caution against a paleoclimatic significance of this result and put forth that the partitioning coefficient of Mg in benthic foraminifera may be sensitive to factors other than temperature. Throughout the remainder of the Cenozoic, the temporal variability among delta18Osw records is similar and similar to longer-term trends in the benthic foraminiferal delta18O record. An exception occurs during the Pliocene when delta18Osw minima in two of the three records suggest reductions in global ice volume that are not apparent in foraminiferal delta18O records, which provides a new perspective to the ongoing debate about the stability of the Antarctic ice sheet. Maximum delta18Osw values recorded during the Pleistocene at Southern Ocean site 747 agree well with values derived from the geochemistry of pore waters (Schrag et al., 1996, doi:10.1126/science.272.5270.1930) further highlighting the value of the new Mg/Ca calibrations of Martin et al. (2002, doi:10.1016/S0012-821X(02)00472-7) and Lear et al. (2002, doi:10.1016/S0016-7037(02)00941-9) applied in this study. We conclude that the application of foraminiferal Mg/Ca ratios allows a refined view of Cenozoic ice volume history despite uncertainties related to the geochemical cycling of Mg and Ca on long time scales.

Benthic foraminifers, isotopes and varimax factors at DSDP Holes 72-517 and 72-518

Pliocene changes in the vertical water mass structure of the western South Atlantic are inferred from changes in benthic foraminiferal assemblages and stable isotopes from DSDP Holes 516A, 517, and 518. Factor analysis of 34 samples from Site 518 reveals three distinct benthic foraminiferal assemblages that have been associated with specific subsurface water masses in the modern ocean. These include a Nuttalides umbonifera assemblage (Factor 1) associated with Antarctic Bottom Water (AABW), a Globocassidulina subglobosa-Uvigerina peregrina assemblage (Factor 2) associated with Circumpolar Deep Water (CPDW), and an Oridorsalis umbonatus-Epistominella exigua assemblage associated with North Atlantic Deep Water (NADW). Bathymetric gradients in d13C between Holes 516A (1313 m), 517 (2963 m), and 518 (3944 m) are calculated whenever possible to monitor the degree of similarity and/or difference in the apparent oxygen utilization (AOU) of water masses located at these depths during the Pliocene. Changes in bathymetric d13C gradients coupled with benthic foraminiferal assemblages record fundamental changes in the vertical water mass structure of the Vema Channel during the Pliocene from 4.1 to 2.7 Ma. At Site 518, the interval from 4.1 to 3.6 Ma is dominated by the N. umbonifera (Factor 1) and O. umbonatus-E. exigua (Factor 3) assemblages. The d13C gradient between Holes 518 (3944 m) and 516A (1313 m) undergoes rapid oscillations during this interval though no permanent increase in the gradient is observed. However, d13C values at Site 518 are clearly lighter during this interval. These conditions may be related to increased bottom water activity associated with the re-establishment of the West Antarctic Ice Sheet in the late Gilbert Chron (-4.2 to 3.6 Ma) (Osborn et al., 1982). The interval from 3.6 to 3.2 Ma is marked by a dominance of the G. subglobosa-U. peregrina (Factor 2) assemblage and lack of a strong d13C gradient between Holes 518 (3944 m) and 516A (1313 m). We suggest that shallow circumpolar waters expanded to depths of a least 3944 m (Site 518) during this time. The most profound faunal and isotopic change occurs at 3.2 Ma, and is marked by dominance of the N. umbonifera (Factor 1) and O. umbonatus-E. exigua (Factor 3) assemblages, a 1.1 per mil enrichment in d18O, and a large negative increase in the d13C gradient between Holes 518 and 516A. These changes at Site 518 record the vertical displacement of circumpolar waters by AABW and NADW. This change in vertical water mass structure at 3.2 Ma was probably related to a global cooling event and/or final closure of the Central American seaway. A comparison of the present-day d13C structure of the Vema Channel with a reconstruction between 3.2 and 2.7 Ma indicates that circulation patterns during this late Pliocene interval were similar to those of the modern western South Atlantic.

(Table S3) Benthic foraminiferal Mg/Ca, Li/Ca and stable isotope ratios of ODP Site 122-761

The sensitivities of benthic foraminiferal Mg/Ca and Li/Ca to bottom water temperature and carbonate saturation state have recently been assessed. Here we present a new approach that uses paired Mg/Ca and Li/Ca records to calculate simultaneous changes in temperature and saturation state. Using previously published records, we first use this approach to document a cooling of deep ocean waters associated with the establishment of the Antarctic ice sheet at the Eocene-Oligocene climate transition. We then apply this approach to new records of the Middle Miocene Climate Transition from ODP Site 761 to estimate variations in bottom water temperature and the oxygen isotopic composition of seawater. We estimate that the oxygen isotopic composition of seawater varied by ~1 per mil between the deglacial extreme of the Miocene Climatic Optimum and the glacial maximum following the Middle Miocene Climate Transition, indicating large amplitude variations in ice volume. However, the longer-term change between 15.3 and 12.5 Ma is marked by a ~1°C cooling of deep waters, and an increase in the oxygen isotopic composition of seawater of ~0.6 per mil. We find that bottom water saturation state increased in the lead up to the Middle Miocene Climate Transition and decreased shortly after. This supports decreasing pCO2 as a driver for global cooling and ice sheet expansion, in agreement with existing boron isotope and leaf stomatal index CO2 records but in contrast to the published alkenone CO2 records.

Ba/Ca ratios in Neogloboquadrina pachyderma of sedimet core 94B17

Down-core samples of the planktonic foraminifer Neogloboquadrina pachyderma sinistral from the Mendeleyev Ridge in the western Arctic Ocean have been analyzed for Ba/Ca and d18O. The apparent distribution coefficient for N. pachyderma sin. is estimated at DBa = 0.22 ± 0.02. A meltwater event is identified at around 11.8 14C kyr BP and is coincident with elevated Ba/Ca ratios. The barium enrichment is believed to be the result of enhanced weathering and erosion following glaciation. Additionally, barium may have desorbed from shelf sediments as sea level rose. Changes in Ba/Ca correlate with the retreat of the Laurentide Ice Sheet and the evolution of the Mackenzie River drainage basin. Therefore maximum Ba/Ca in Arctic surface waters at 11.8 ka may be indicative of an increase in the export of freshwater from the Arctic to the North Atlantic, potentially contributing to the onset of the Younger Dryas. This work suggests that Ba/Ca in planktonic foraminifera may be a useful indicator of the timing and processes associated with deglaciation.

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.

Clay mineralogy of ODP Site 188-1165

The Ocean Drilling Program Leg 188 Site 1165 was drilled on the Wild Drift on the Continental Rise off Prydz Bay, East Antarctica to a total depth of 999.1 meters below seafloor (mbsf). It recovered an extensive suite of terrigenous and hemipelagic sediments of early Miocene to Pleistocene age. Of special interest in this study is the sediment column between 0 and 50 mbsf, which consists of a well-preserved section of Pliocene-Pleistocene-age sediments that was sampled at 10-cm intervals. Multiproxy study of this interval could show possible intervals of expansions of the ice-sheet across the continental shelves and express the climatic evolution in Antarctica, particularly during the 'middle' Pliocene warm period (3.15 to 2.85 Ma) which may provide an indication of how the Earth responds to a rise of its surface temperature. According to the existing age model, the upper 50 mbsf stratigraphic sequence of Site 1165 reaches back to ~4.9 Ma. Throughout this interval, the clay-mineral content is characterized by fluctuations of individual clay minerals, particularly smectite and chlorite. The smectite concentration varies mainly between 0% and 30%. Illite fluctuates less between 50% and 80%, and kaolinite varies mainly between 10% and 20%. Chlorite concentrations are mainly 0% to 10%. There is also a noticeable change in magnetic susceptibility at ~34 mbsf that is clearly indicated in the composition of the clay-mineral suite. At this level, smectite decreases and illite, kaolinite and chlorite show some variability. In particular, there is a slight but persistent increase in chlorite. The results from the Plio-Pleistocene transition, with evenly fluctuating smectite and illite contents and the gradually increasing chlorite content, may indicate cooler conditions compared to the mid-Pliocene conditions. Slight increase in illite content and decrease in smectite content towards Pleistocene supports the previous assumption. The results from the mid-Pliocene with the increasing smectite content and decreasing illite content may indicate warmer and possibly interglacial conditions.

(Table S1) Flux gate discharge and centerline discharge of Helheim, Kangerdlugssuaq and Jakobshavn glacier, Greenland

Acceleration of Greenland's three largest outlet glaciers, Helheim, Kangerdlugssuaq and Jakobshavn Isbræ, accounted for a substantial portion of the ice sheet's mass loss over the past decade. Rapid changes in their discharge, however, make their cumulative mass-change uncertain. We derive monthly mass balance rates and cumulative balance from discharge and surface mass balance (SMB) rates for these glaciers from 2000 through 2010. Despite the dramatic changes observed at Helheim, the glacier gained mass over the period, due primarily to the short duration of acceleration and a likely longer-term positive balance. In contrast, Jakobshavn Isbræ lost an equivalent of over 11 times the average annual SMB and loss continues to accelerate. Kangerdlugssuaq lost over 7 times its annual average SMB, but loss has returned to the 2000 rate. These differences point to contrasts in the long-term evolution of these glaciers and the danger in basing predictions on extrapolations of recent changes.

Antarctic dataset in NetCDF format

The dataset described in this document has been put together for the purposes of numerical ice sheet modelling of the Antarctic Ice Sheet (AIS), containing data on the ice sheet configuration (e.g. ice surface and ice thickness) and boundary conditions, such as the surface air temperature and accumulation. It is now possible to download a community ice sheet model (e.g. Glimmer-CISM, Rutt et al., 2009 doi:10.1029/2008JF001015), but without adequate data it is difficult to utilise such models. More specifically, ice sheet models that are initialised and run forward from the present day ice sheet configuration, need input data to represent the present-day ice sheet configuration as closely as possible (unlike those spun-up from ice free conditions, which only require the bed/bathymetry). Whilst the BEDMAP dataset (Lythe et al., 2001) was a step forward when it was made, there are a number of inconsistencies within the dataset (see Section 3), and since its release, more data has become available. The dataset described here incorporates some major new datasets (e.g. AGASEA/BBAS ice thickness, Nitsche et al. (2006) bathymetry doi:10.1029/2007GC001694), but by no means incorporates all the new data available. This considerable task is left for a 'BEDMAP2', (an updated version of BEDMAP), however, the processing carried out in this document illustrates the requirements of a dataset for the purpose of high resolution ice sheet modelling, and bridges the gap until a BEDMAP2 is published. It is envisaged, however, that updated versions of the data set will be made available periodically when new regional data sets become available and can be readily incorporated.

Table 5 Comparison of surface elevations by Seasat and optical survey along EGIG line

Surface elevation maps of the southern half of the Greenland subcontinent are produced from radar altimeter data acquired by the Seasat satellite. A summary of the processing procedure and examples of return waveform data are given. The elevation data are used to generate a regular grid which is then computer contoured to provide an elevation contour map. Ancillary maps show the statistical quality of the elevation data and various characteristics of the surface. The elevation map is used to define ice flow directions and delineate the major drainage basins. Regular maps of the Jakobshavns Glacier drainage basin and the ice divide in the vicinity of Crete Station are presented. Altimeter derived elevations are compared with elevations measured both by satellite geoceivers and optical surveying.

Magnetic hysteresis parameters of ODP Site 178-1096 (Table T1)

During Ocean Drilling Program (ODP) Leg 178, we drilled three sites on sediment drifts deposited on the continental rise on the western margin of the Antarctic Peninsula. These hemipelagic drifts were targeted for their potential to preserve a continuous record of the behavior of the West Antarctic Ice Sheet over the last 10 m.y. It has been proposed that drift development is linked to advances and retreats of the Antarctic continental ice sheet (Pudsey and Camerlenghi, 1998, doi:10.1017/S0954102098000376, and references therein; Barker, Camerlenghi, Acton, et al., 1999, doi:10.2973/odp.proc.ir.178.1999). However, the sediment is characterized by a very low carbonate content, with foraminifers restricted to very narrow intervals. This lack of carbonate precludes the construction of a delta18O or CaCO3 stratigraphy, depriving these sites of an important chronologic tool and global ice volume proxy.

Detriral flux of sediment cores in the southeast Indian Ocean

Because of a close relationship between detrital flux variations and magnetic susceptibility (MS) flux (MS cm**3 of bulk sediment multiplied by the linear sedimentation rate) variations in the southeast Indian basin of the southern ocean, MS flux profiles have been used to examine the spatial and temporal detrital flux changes in this basin during the last climatic cycle. Results indicate a general increase in detrital material input during the coldest periods, suggesting a widespread phenomenon, at least on the basin scale. Mineralogical data, geochemical data, and 87Sr/86Sr isotopic ratios have been used to determine the origin and transport mechanisms responsible for increased detrital flux during glacial periods. Mineralogical and geochemical data show that these glacial 'highs' are due to increases in both Kerguelen-Crozet volcanic and Antarctic detrital inputs. The 87Sr/86Sr isotopic composition of the >45-µm fraction indicates that the Kerguelen-Crozet province contributes to at least 50% of the coarse particule input to the west. This contribution decreases eastward to reach less than 10%. These tracers clearly indicate that the Crozet-Kerguelen province was a major source region of detrital in the western part of the basin during glacial times. In contrast, material of Antarctic origin is well represented in the whole basin (fine and coarse fractions). Because of the minor amount of coarse particles in the sediments, volcanic particles from Kerguelen and crustal particles from Antarctica have most probably been transported by the Antarctic bottom water current and/or the Circumpolar deepwater current during glacial periods as is the case today. Nevertheless, the presence of coarse particles even in low amount suggests also a transport by ice rafting (sea-ice and icebergs), originated from both Kerguelen and Antarctic sources. However, the relative importance of both hydrographic and ice-rafting modes of transport cannot be identified accurately with our data. During low sea level stands (glacial maximum periods), increasing instability and erosion of the continental platform and shallow plateaus could have resulted in a more efficient transfer of crustal and volcano-detrital material to the Southeast Indian basin. At the same time, extension of the grounded ice shelves over the continental margins and increase in the erosion rate of the Antarctic ice sheet could have induced a greater input of ice rafted detritus (IRD) to southern ocean basins. Enhancement of the circumpolar deepwater current strength might have also carried a more important flux of detrital material from Kerguelen. However, an increase in the bottom water flow is not necessarily required.

Neogene benthic foraminiferal faunas of a southwest Pacific bathyal depth transect

Temporal changes in benthic foraminiferal assemblages were quantitatively examined (> 63 µm fraction) in four southwest Pacific deep-sea Neogene sequences in a depth transect between approximately 1300 and 3200 m to assist in evaluating paleoeeanographic history. The most conspicuous changes in benthic foraminiferal assemblages occurred in association with paleoclimatic changes defined at least in part by oxygen isotopic changes. The largest, centered at ~15 Ma (early Middle Miocene), is represented by an increase in the relative frequencies of Epistominella exigua, which underwent a major upward depth migration at that time. This was contemporaneous with the well-known positive oxygen isotopic shift in the early Middle Miocene. In Sites 588 and 590, most of the increase in relative abundances of E. exigua occurred during the middle to later part of the ~80 shift, following major growth of the east Antarctic ice sheet. Later assemblage changes occurred at 8.5 and 6.5 Ma. These associations indicate that the benthic foraminiferal assemblages in this depth transect largely adjusted to changes in deep waters related to Antarctic cryospheric evolution. In general, the Neogene benthic foraminiferal assemblages in this region underwent little change during the last 23 million years. This faunal conservatism suggests that deep-sea environments underwent relatively little change in the southwest Pacific during much of the Neogene. Although paleoceanographic changes did occur, partly in response to highlatitude cryospheric evolution, these were not of sufficient magnitude to create major deep-sea faunal changes in this part of the ocean. The benthic foraminiferal assemblages are dominated by individuals smaller than 150 µm. Most taxonomic turnover occurred in the larger (> 150 µm) size fractions.