Stable isotope ratios of foraminifera from last glacial-interglacial sediments of the southern high latitudes

A set of numerical equations is developed to estimate past sea surface temperatures (SST) from fossil Antarctic diatoms. These equations take into account both the biogeographic distribution and experimentally derived silica dissolution. The data represent a revision and expansion of a floral data base used previously and includes samples resulting from progressive opal dissolution experiments. Factor analysis of 166 samples (124 Holocene core top and 42 artificial samples) resolved four factors. Three of these factors depend on the water mass distribution (one Subantarctic and two Antarctic assemblages); factor 4 corresponds to a 'dissolution assemblage'. Inclusion of this factor in the data analysis minimizes the effect of opal dissolution on the assemblages and gives accurate estimates of SST over a wide range of biosiliceous dissolution. A transfer function (DTF 166/34/4) is derived from the distribution of these factors versus summer SST. Its standard error is +/- 1°C in the -1 to +10 °C summer temperature range. This transfer function is used to estimate SST changes in two southern ocean cores (43°S and 55°S) which cover the last climatic cycle. The time scale is derived from the changes in foraminiferal oxygen and carbon isotopic ratios. The reconstructed SST records present strong analogies with the air temperature record over Antarctica at the Vostok site, derived from changes in the isotopic ratio of the ice. This similarity may be used to compare the oceanic isotope stratigraphy and the Vostok time scale derived from ice flow model. The oceanic time scale, if taken at face value, would indicate that large changes in ice accumulation rates occurred between warm and cold periods.

Pb isotopes in individual sand-sized (>150 µm) ice-rafted feldspars deposited at subpolar North Atlantic Ocean Deep Sea Drilling Project Site 611 during the late Pliocene (Marine Isotope Stages, MIS, G6 through 100)

The onset of abundant ice-rafted debris (IRD) deposition in the Nordic Seas and subpolar North Atlantic Ocean 2.72 millions of years ago (Ma) is thought to record the Pliocene onset of major northern hemisphere glaciation (NHG) due to a synchronous advance of North American Laurentide, Scandinavian and Greenland ice-sheets to their marine calving margins during marine isotope stage (MIS) G6. Numerous marine and terrestrial records from the Nordic Seas region indicate that extensive ice sheets on Greenland and Scandinavia increased IRD inputs to these seas from 2.72 Ma. The timing of ice-sheet expansion on North America as tracked by IRD deposition in the subpolar North Atlantic Ocean, however, is less clear because both Europe and North America are potential sources for icebergs in this region. Moreover, cosmogenic-dating of terrestrial tills on North America indicate that the Laurentide Ice Sheet did not extend to ~39°N until 2.4 ±0.14 Ma, at least 180 ka after the onset of major IRD deposition at 2.72 Ma. To address this problem,we present the first detailed analysis of the geochemical provenance of individual sand-sized IRD deposited in the subpolar North Atlantic Ocean between MIS G6 and 100 (~2.72-2.52 Ma). IRD provenance is assessed using laser ablation lead (Pb) isotope analyses of single ice-rafted (>150 mm) feldspar grains. To track when an ice-rafting setting consistent with major NHG first occurred in the North Atlantic Ocean during the Pliocene intensification of NHG (iNHG), we investigate when the Pb-isotope composition (206Pb/204Pb, 207Pb/204Pb, 208Pb/204Pb) of feldspars deposited at DSDP Site 611 first resembles that determined for IRD deposited at this site during MIS 100, the oldest glacial for which there exists convincing evidence for widespread glaciation of North America. Whilst Quaternary-magnitude IRD fluxes exist at Site 611 during glacials from 2.72 Ma, we find that the provenance of this IRD is not constant. Instead, we find that the Pb isotope composition of IRD at our study site is not consistent with major NHG until MIS G2 (2.64 Ma). We hypothesise that IRD deposition in the North Atlantic Ocean prior to MIS G2 was dominated by iceberg calving from Greenland and Scandinavia. We further suggest that the grounding line of continental ice on Northeast America may not have extended onto the continental shelf and calved significant numbers of icebergs to the North Atlantic Ocean during glacials until 2.64 Ma.

Reconstructed ice-sheet limit in the Kara and Barents Sea area at about 50,000 years ago (Markhida Line)

Glacial landforms in northern Russia, from the Timan Ridge in the west to the east of the Urals, have been mapped by aerial photographs and satellite images supported by field observations. An east-west trending belt of fresh hummock-and-lake glaciokarst landscapes has been traced to the north of 67°N. The southern boundary of these landscapes is called the Markhida Line, which is interpreted as a nearly synchronous limit of the last ice sheet that affected this region. The hummocky landscapes are subdivided into three types according to the stage of postglacial modification: Markhida, Harbei and Halmer. The Halmer landscape on the Uralian piedmont in the east is the freshest, whereas the westernmost Markhida landscape is more eroded. The west- east gradient in morphology is considered to be a result of the time-transgressive melting of stagnant glacier ice and of the underlying permafrost. The pattern of ice-pushed ridges and other directional features reflects a dominant ice flow direction from the Kara Sea shelf. Traces of ice movement from the central Barents Sea are only discernible in the Pechora River left bank area west of 50°E. In the Polar Urals the horseshoe-shaped end moraines at altitudes of up to 560 m a.s.l. reflect ice movement up-valley from the Kara Ice Sheet, indicating the absence of a contemporaneous ice dome in the mountains. The Markhida moraines, superimposed onto the Eemian strata, represent the maximum ice sheet extent in the western part of the Pechora Basin during the Weichselian. The Markhida Line truncates the huge arcs of the Laya-Adzva and Rogovaya ice-pushed ridges protruding to the south. The latter moraines therefore reflect an older ice advance, probably also of Weichselian age. Still farther south, fluvially dissected morainic plateaus without lakes are of pre-Eemian age, because they plunge northwards under marine Eemian sediments. Shorelines of the large ice-dammed Lake Komi, identified between 90 and 110 m a.s.l. in the areas south of the Markhida Line, are radiocarbon dated to be older than 45 ka. The shorelines, incised into the Laya-Adzva moraines, morphologically interfinger with the Markhida moraines, indicating that the last ice advance onto the Russian mainland reached the Markhida Line during the Middle or Early Weichselian, before 45 ka ago.

Planktic and benthic 14C reservoir ages, calibrated by a suite of 14C plateaus in the glacial-to-deglacial Suigetsu atmospheric 14C record

We here present a compilation of planktic and benthic 14C reservoir ages for the Last Glacial Maximum (LGM) and early deglacial from 11 key sites of global ocean circulation in the Atlantic and Indo-Pacific Ocean. The ages were obtained by 14C plateau tuning, a robust technique to derive both an absolute chronology for marine sediment records and a high-resolution record of changing reservoir/ventilation ages (Delta14C values) for surface and deep waters by comparing the suite of planktic 14C plateaus of a sediment record with that of the atmospheric 14C record (Sarnthein et al., 2007, doi:10.1029/173GM13). Results published thus far used as atmospheric 14C reference U/Th-dated corals, the Cariaco planktic record, and speleothems (Fairbanks et al., 2005, doi:10.1016/j.quascirev.2005.04.007; Hughen et al., 2006, doi:10.1016/j.quascirev.2006.03.014; Beck et al., 2001, doi:10.1023/A:1008175728826). We have now used the varve-counted atmospheric 14C record of Lake Suigetsu terrestrial macrofossils (Ramsey et al., 2012, doi:10.1126/science.1226660) to recalibrate the boundary ages and reservoir ages of the seven published records directly to an atmospheric 14C record. In addition, the results for four new cores and further planktic results for four published records are given. Main conclusions from the new compilation are: (1) The Suigetsu atmospheric 14C record on its varve counted time scale reflects all 14C plateaus, their internal structures and relative length previously identified, but implies a rise in the average 14C plateau age by 200-700 14C yr during LGM and early deglacial times. (2) Based on different 14C ages of coeval atmospheric and planktic 14C plateaus, marine surface water Delta14C may have temporarily dropped to an equivalent of ~0 yr in low-latitude lagoon waters, but reached >2500 14C yr both in stratified subpolar waters and in upwelled waters such as in the South China Sea. These values differ significantly from a widely assumed constant global planktic Delta14C value of 400 yr. (3) Suites of deglacial planktic Delta14C values are closely reproducible in 14C records measured at neighboring core sites. (4) Apparent deep-water 14C ventilation ages (equivalents of benthic Delta14C), deduced from the sum of planktic Delta14C and coeval benthic-planktic 14C differences, vary from 500 up to >5000 yr in LGM and deglacial ocean basins.

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.

Dust fluxes and iron fertilization in Holocene and Last Glacial Maximum climates

Mineral dust aerosols play a major role in present and past climates. To date, we rely on climate models for estimates of dust fluxes to calculate the impact of airborne micronutrients on biogeochemical cycles. Here we provide a new global dust flux data set for Holocene and Last Glacial Maximum (LGM) conditions based on observational data. A comparison with dust flux simulations highlights regional differences between observations and models. By forcing a biogeochemical model with our new data set and using this model's results to guide a millennial-scale Earth System Model simulation, we calculate the impact of enhanced glacial oceanic iron deposition on the LGM-Holocene carbon cycle. On centennial timescales, the higher LGM dust deposition results in a weak reduction of <10?ppm in atmospheric CO2 due to enhanced efficiency of the biological pump. This is followed by a further ~10?ppm reduction over millennial timescales due to greater carbon burial and carbonate compensation.

Gridded climate data from 5 Global Climate Models (GCM) of the Last Glacial Maximum (LGM) downscaled to 30 arc seconds for Europe, with links to NetCDF files

Studies on the impact of historical, current and future global change require very high-resolution climate data (less or equal 1km) as a basis for modelled responses, meaning that data from digital climate models generally require substantial rescaling. Another shortcoming of available datasets on past climate is that the effects of sea level rise and fall are not considered. Without such information, the study of glacial refugia or early Holocene plant and animal migration are incomplete if not impossible. Sea level at the last glacial maximum (LGM) was approximately 125m lower, creating substantial additional terrestrial area for which no current baseline data exist. Here, we introduce the development of a novel, gridded climate dataset for LGM that is both very high resolution (1km) and extends to the LGM sea and land mask. We developed two methods to extend current terrestrial precipitation and temperature data to areas between the current and LGM coastlines. The absolute interpolation error is less than 1°C and 0.5 °C for 98.9% and 87.8% of all pixels for the first two 1 arc degree distance zones. We use the change factor method with these newly assembled baseline data to downscale five global circulation models of LGM climate to a resolution of 1km for Europe. As additional variables we calculate 19 'bioclimatic' variables, which are often used in climate change impact studies on biological diversity. The new LGM climate maps are well suited for analysing refugia and migration during Holocene warming following the LGM.

Mineralogy of ODP Site 105-645

Cores from the upper 70 meters below seafloor (mbsf) (upper Pleistocene) at Ocean Drilling Program (ODP) Site 645 in Baffin Bay show dramatic meter-scale changes in color and mineralogy. Below this interval, mineralogical changes are more gradual to the top of the Miocene at about 550 mbsf. The Pliocene-Pleistocene section can be divided into five facies: Facies 1 - massive, poorly sorted, gravel-bearing muds; Facies 2 - gray silty clays and silty muds; Facies 3 - laminated detricarbonate silty muds; Facies 4 - silty sand and sandy silt; and Facies 5 - poorly sorted muddy sands and silty muds. Facies 4 and 5 are restricted to the Pliocene section below depths of about 275 mbsf. The mineralogical/color cycles in the upper 70 mbsf are the result of alternations between Facies 2 and three lithotypes of Facies 1: lithotype A - tan-colored, carbonate-rich, gravel-bearing mud; lithotype B - weak, red-colored, gravel-bearing mud rich in sedimentary rock fragments; and lithotype C - gray, gravel-bearing mud. A fourth lithotype, D, is restricted to depths of 168-275 mbsf and is dark gray, carbonate-poor, gravel-bearing mud. We believe that all lithotypes of Facies 1 and the sand and gravel fractions of Facies 2 and 3 were deposited by ice rafting. Depositional processes for Facies 4 and 5 probably include ice rafting and bottom- and turbidity-current transport. Data from petrographic analyses of light and heavy sand-sized grains and X-ray analyses of silt- and clay-size fractions suggest that tan-colored sediments (lithotype A of Facies 1; Facies 3) were derived mainly from Paleozoic carbonates of Ellesmere, Devon, and northern Baffin islands. Weak red sediments (lithotype B) contain significant red sedimentary clasts, reworked quartzarenite grains and clasts, and rounded colorless garnets, all derived from Proterozoic sequences of the Borden and Thule basins, and from minor Mesozoic red beds. Other sediments in the upper 335 mbsf at Site 645 contain detritus from a heterogeneous mixture of sources, including Precambrian shield terranes around Baffin Bay. Sediments from 335 to 550 mbsf (Facies 5) are rich in friable sedimentary clasts and detrital micas and contain glauconite and, in a few samples, reworked diatoms. These components suggest derivation from poorly consolidated Mesozoic-Tertiary sediments in coastal outcrops and beneath the modern shelves of northeastern Baffin Island and western Greenland. For the upper Pleistocene section (about 0-100 mbsf), marked mineralogical cyclicity is attributed to fluctuating glacial margins, calving rates, and iceberg melting rates, particularly around the northern end of Baffin Bay. Tan-colored, carbonate-rich units were derived at times of maximum advance of glaciers on Ellesmere and Devon islands, during relatively warm intervals induced by incursion of warm Atlantic surface water into the bay. At the beginning of these warmer episodes, most icebergs were contributed by glaciers near sea level around the Arctic channels, which resulted in deposition of weak red, ice-rafted units rich in Proterozoic sedimentary clasts.

Biogenic opal and carbonate concentration in late Neogene sediments of ODP Site 114-704

A preliminary composite depth section was generated for Site 704 by splicing Holes 704A and 704B together over the interval 0-350 mbsf (0-9 m.y.). High-resolution carbonate and opal data from the cores were correlated with the calcium and silicon signals from the GST logging run in Hole 704B to identify missing and disturbed intervals in the cores. Paleomagnetic and biostratigraphic age boundaries were then transferred to the composite depth records to obtain an age model, and sedimentation rates were calculated by linear interpolation between datums. Algorithms relating measured dry-bulk density to carbonate content and depth were generated to produce predicted values of density for every sample. Accumulation rates of bulk, carbonate, opal, and terrigenous sediment components were then computed to generate a record of sediment deposition on the Meteor Rise that has a resolution of better than 200,000 yr for the period from 8.6 to 1.0 m.y. From 8.6 to 2.5 m.y., bulk-accumulation rates on the Meteor Rise averaged less than 2 g/cm**2/1000 yr and were dominated by carbonate deposition. The first significant opal deposition (6.0 m.y.) punctuated a brief (less than 0.6 Ma) approach of the Polar Front Zone (PFZ) northward that heralded a period of increasing severity of periodic carbonate dissolution events (terrigenous maxima) that abruptly terminated at 4.8 m.y. (base of the Thvera Subchron), synchronous with the reflooding of the Mediterranean after the Messinian salinity crisis. From 4.8 to 2.5 m.y., carbonate again dominated deposition, and the PFZ was far south except during brief northward excursions bracketing 4.2-3.9, 3.3-2.9, and 2.8-2.7 m.y. At 2.5 m.y., all components of bulk-accumulation rates increased dramatically (up to 15 g/cm2/1000 yr), and by 2.4 m.y., a pattern of alternating, high-amplitude carbonate and opal cyclicity marked the initiation of rapid glacial to interglaci·l swings in the position of the PFZ, synchronous with the "onset" of major Northern Hemisphere glaciation. Both mass-accumulation rates and the amplitude of the cycles decreased by about 2 m.y., but opal accumulation rates remained high up through the base of the Jaramillo (0.98 m.y.). From 1.9 to 1 m.y., the record is characterized by moderate amplitude fluctuations in carbonate and opal. This record of opal accumulation rates is interpreted as a long-term "Polar Front Indicator" that monitors the advance and retreat of the opal-rich PFZ northward (southward) toward (away from) the Meteor Rise in the subantarctic sector of the South Atlantic Ocean. The timing of PFZ migrations in the subantarctic South Atlantic Ocean is remarkably similar to Pliocene-Pleistocene climate records deduced from benthic oxygen isotope records in the North Atlantic Ocean (Raymo et al., 1989, doi:10.1029/PA004i004p00413; Ruddiman et al., 1989, doi:10.1029/PA004i004p00353). These include northward migrations during "cold" intervals containing strong glacial isotope stages (2.4-2.3, 2.1-2.0, 1.95-1.55, 1.45-1.30 m.y. and at about 1.13 and 1.09 m.y.) and southward migrations during "warm" intervals containing weak glacial and/or strong interglacial stages (2.45-2.40, 2.30-2.10, 2.00-1.95, 1.52-1.45, 1.30-1.18, 1.11, and 1.06-0.93 m.y.). Although our preliminary composite record is not continuous (some stages are obviously missing), there is hope that future work will identify these missing intervals in the as yet incomplete Hole 704B and will extend this high-resolution Southern Hemisphere climate record back to 8.6 m.y.

Beryllium-10 and aluminium-26 concentration and production rates in glacial deposits from the Quartermain Mountains

Many glacial deposits in the Quartermain Mountains, Antarctica present two apparent contradictions regarding the degradation of unconsolidated deposits. The glacial deposits are up to millions of years old, yet they have maintained their meter-scale morphology despite the fact that bedrock and regolith erosion rates in the Quartermain Mountains have been measured at 0.1-4.0 m/Ma. Additionally, ground ice persists in some Miocene-aged soils in the Quartermain Mountains even though modeled and measured sublimation rates of ice in Antarctic soils suggest that without any recharge mechanisms ground ice should sublimate in the upper few meters of soil on the order of 10**3 to 10**5 years. This paper presents results from using the concentration of cosmogenic nuclides beryllium-10 (10Be) and aluminum-26 (26Al) in bulk sediment samples from depth profiles of three glacial deposits in the Quartermain Mountains. The measured nuclide concentrations are lower than expected for the known ages of the deposits, erosion alone does not always explain these concentrations, and deflation of the tills by the sublimation of ice coupled with erosion of the overlying till can explain some of the nuclide concentration profiles. The degradation rates that best match the data range 0.7-12 m/Ma for sublimation of ice with initial debris concentrations ranging 12-45% and erosion of the overlying till at rates of 0.4-1.2 m/Ma. Overturning of the tills by cryoturbation, vertical mixing, or soil creep is not indicated by the cosmogenic nuclide profiles, and degradation appears to be limited to within a few centimeters of the surface. Erosion of these tills without vertical mixing may partially explain how some glacial deposits in the Quartermain Mountains maintain their morphology and contain ground ice close to the surface for millions of years.

Paleoenvironmental changes in the western Mediterranean since the last glacial maximum: High resolution multiproxy record from the Algero-Balearic basin

The present study uses a multiproxy approach in order to further understand the evolution of climate responses in the western Mediterranean as of the Last Glacial Maximum. Sediments from ODP Site 975 in the Algero-Balearic basin have been analysed at high resolution, both geochemically andmineralogicallly. The resulting data have been used as proxies to establish a sedimentary regime, primary marine productivity, the preservation of the proxies and oxygen conditions. Fluctuations in detrital element concentrations were mainly the consequence of wet/arid oscillations. Productivity has been established using Ba excess, according to which marine productivity appears to have been greatest during cold events Heinrich 1 and Younger Dryas. The S1 time interval was not as marked by increases in productivity as was the eastern Mediterranean. In contrast, the S1 interval was first characterized by a decreasing trend and then by a fall in productivity after the 8.2 ky BP dry-cold event. Since then productivity has remained low. Here we report that there was an important redox event in this basin, probably a consequence of the major oceanographic circulation change occurring in the western Mediterranean at 7.7 ky BP. This circulation change led to reventilation as well as to diagenetic remobilization of redox-sensitive elements and organic matter oxidation. Comparisons between our paleoceanographic reconstruction for this basin and those regarding other Mediterranean basins support the hypothesis that across the Mediterranean there were different types of responses to climate forcing mechanism. The Algero-Balearic basin is likely to be a key area for further understanding of the relationships between the North Atlantic and the eastern Mediterranean basins.