Modeling the sensitivity of precipitation oxygen isotopes to the Last Glacial Maximum boundary conditions: A study using Community Atmosphere Model (CAM3.0)

To understand the validity of d18O proxy records as indicators of past temperature change, a series of experiments was conducted using an atmospheric general circulation model fitted with water isotope tracers (Community Atmosphere Model version 3.0, IsoCAM). A pre-industrial simulation was performed as the control experiment, as well as a simulation with all the boundary conditions set to Last Glacial Maximum (LGM) values. Results from the pre-industrial and LGM simulations were compared to experiments in which the influence of individual boundary conditions (greenhouse gases, ice sheet albedo and topography, sea surface temperature (SST), and orbital parameters) were changed each at a time to assess their individual impact. The experiments were designed in order to analyze the spatial variations of the oxygen isotopic composition of precipitation (d18Oprecip) in response to individual climate factors. The change in topography (due to the change in land ice cover) played a significant role in reducing the surface temperature and d18Oprecip over North America. Exposed shelf areas and the ice sheet albedo reduced the Northern Hemisphere surface temperature and d18Oprecip further. A global mean cooling of 4.1 °C was simulated with combined LGM boundary conditions compared to the control simulation, which was in agreement with previous experiments using the fully coupled Community Climate System Model (CCSM3). Large reductions in d18Oprecip over the LGM ice sheets were strongly linked to the temperature decrease over them. The SST and ice sheet topography changes were responsible for most of the changes in the climate and hence the d18Oprecip distribution among the simulations.

Analysis of the iron oxide assemblages in the ANDRILL AND-1B drill core, Ross Sea, Antarctica

The AND-1B drill core recovered a 13.57 million year Miocene through Pleistocene record from beneath the McMurdo Ice Shelf in Antarctica (77.9°S, 167.1°E). Varying sedimentary facies in the 1285 m core indicate glacial-interglacial cyclicity with the proximity of ice at the site ranging from grounding of ice in 917 m of water to ice free marine conditions. Broader interpretation of climatic conditions of the wider Ross Sea Embayment is deduced from provenance studies. Here we present an analysis of the iron oxide assemblages in the AND-1B core and interpret their variability with respect to wider paleoclimatic conditions. The core is naturally divided into an upper and lower succession by an expanded 170 m thick volcanic interval between 590 and 760 m. Above 590 m the Plio-Pleistocene glacial cycles are diatom rich and below 760 m late Miocene glacial cycles are terrigenous. Electron microscopy and rock magnetic parameters confirm the subdivision with biogenic silica diluting the terrigenous input (fine pseudo-single domain and stable single domain titanomagnetite from the McMurdo Volcanic Group with a variety of textures and compositions) above 590 m. Below 760 m, the Miocene section consists of coarse-grained ilmenite and multidomain magnetite derived from Transantarctic Mountain lithologies. This may reflect ice flow patterns and the absence of McMurdo Volcanic Group volcanic centers or indicate that volcanic centers had not yet grown to a significant size. The combined rock magnetic and electron microscopy signatures of magnetic minerals serve as provenance tracers in both ice proximal and distal sedimentary units, aiding in the study of ice sheet extent and dynamics, and the identification of ice rafted debris sources and dispersal patterns in the Ross Sea sector of Antarctica.

Radiocarbon ages and clay minerals in sediment cores obtained during a 2009 Jan Mayen cruise to the Norwegian Sea

The origin of two acoustic sediment units has been studied based on lithological facies, chronology and benthic stable isotope values as well as on foraminifera and clay mineral assemblages in six marine sediment cores from Kveithola, a small trough west of Spitsbergenbanken on the western Barents Sea margin. We have identified four time slices with characteristic sedimentary environments. Before c. 14.2 cal. ka, rhythmically laminated muds indicate extensive sea ice cover in the area. From c. 13.9 to 14.2 cal. ka, muds rich in ice-rafted debris were deposited during the disintegration of grounded ice on Spitsbergenbanken. From c. 10.3 to 13.1 cal. ka, sediments with heterogeneous lithologies suggest a shifting influence of suspension settling and iceberg rafting, probably derived from a decaying Barents Sea Ice Sheet in the inner-fjord and land areas to the north of Kveithola. Holocene deposition was episodic and characterized by the deposition of calcareous sands and shell debris, indicative of strong bottom currents. We speculate that a marked erosional boundary at c. 8.2 cal. ka may have been caused by the Storegga tsunami. Whilst deposition was sparse during the Holocene, Kveithola acted as a sediment trap during the preceding deglaciation. Investigation of the deglacial sediments provides unprecedented details on the dynamics and timing of glacial retreat from Spitsbergenbanken.

Age determination of sediment cores from the Weddell Sea

We report here the results of a study aimed at providing radiometric age control on glacial events in the Weddell Sea during the late Quaternary. Sediment cores from the eastern continental shelf, where the East Antarctic ice sheet was grounded, have recovered glacial-marine sediments resting on tills and the latter deposits predate the isotope stage 2 last glacial maximum. Sediment cores from the continental slope and rise sampled a prominent ice-rafted debris layer, and radiocarbon ages indicate that this ice-rafting event took place prior to 26 000 yr B.P. Thus, the combined data indicate that significant deglaciation of the Weddell Sea continental shelf took place prior to the last glacial maximum. Our data also suggest that the ice masses that border the Weddell Sea are more extensive than they were during the previous glacial minimum.

Proglacial river dataset from the Akuliarusiarsuup Kuua River northern tributary, Southwest Greenland, 2008 - 2013, version 2.0

The dataset contains measurements of river stage and discharge for one sites along the Akuliarusiarsuup Kuua River's northern tributary, with 30 minute temporal resolution between June 2008 and August 2013 This river is a tributary to the Watson River discharging into Kangerlussuaq Fjord by the town of Kangerlussuaq, Southwest Greenland. Additional data of water temperature, air pressure are also provided. Compared to version 1.0 of the dataset, this dataset used a total of 36 in situ discharge observations collected between 2008 and 2012 to construct the rating curve. Furthermore, data of Station AK-004-001 between 2010-09-06T11:30 to 2010-09-07T13:30 have been removed from version 2.0 because these values were likely caused by backflow when a jokulhlaup from a large glacier dammed lake caused increased water levels in the downstreams lake. Thus, data measured at AK-004-001 between 2010-09-06T11:30 to 2010-09-07T13:30 are not representative for the AK-004 catchment.

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.