Age determinaton and age model of sediments from the Bay of Biscay

Continuous high-resolution mass accumulation rates (MAR) and X-ray fluorescence (XRF) measurements from marine sediment records in the Bay of Biscay (NE Atlantic) have allowed the determination of the timing and the amplitude of the 'Fleuve Manche' (Channel River) discharges during glacial stages MIS 10, MIS 8, MIS 6 and MIS 4-2. These results have yielded detailed insight into the Middle and Late Pleistocene glaciations in Europe and the drainage network of the western and central European rivers over the last 350 kyr. This study provides clear evidence that the 'Fleuve Manche' connected the southern North Sea basin with the Bay of Biscay during each glacial period and reveals that 'Fleuve Manche' activity during the glaciations MIS 10 and MIS 8 was significantly less than during MIS 6 and MIS 2. We correlate the significant 'Fleuve Manche' activity, detected during MIS 6 and MIS 2, with the extensive Saalian (Drenthe Substage) and the Weichselian glaciations, respectively, confirming that the major Elsterian glaciation precedes the glacial MIS 10. In detail, massive 'Fleuve Manche' discharges occurred at ca 155 ka (mid-MIS 6) and during Termination I, while no significant discharges are found during Termination II. It is assumed that a substantial retreat of the European ice sheet at ca 155 kyr, followed by the formation of ice-free conditions between the British Isles and Scandinavia until Termination II, allowed meltwater to flow northwards through the North Sea basin during the second part of the MIS 6. We assume that this glacial pattern corresponds to the Warthe Substage glacial maximum, therefore indicating that the data presented here equates to the Drenthe and the Warthe glacial advances at ca 175-160 ka and ca 150-140 ka, respectively. Finally, the correlation of our records with ODP site 980 reveals that massive 'Fleuve Manche' discharges, related to partial or complete melting of the European ice masses, were synchronous with strong decreases in both the rate of deep-water formation and the strength of the Atlantic thermohaline circulation. 'Fleuve Manche' discharges over the last 350 kyr probably participated, with other meltwater sources, in the collapse of the thermohaline circulation by freshening the northern Atlantic surface water.

Alkenone stable carbon isotopes, carbonate stable carbon and oxygen isotopes, planktic foraminifera boron isotopes and atmospheric CO2 calculations and box model results from Ras il-Pellegrin section, Malta

The development of a permanent, stable ice sheet in East Antarctica happened during the middle Miocene, about 14 million years (Myr) ago. The middle Miocene therefore represents one of the distinct phases of rapid change in the transition from the "greenhouse" of the early Eocene to the "icehouse" of the present day. Carbonate carbon isotope records of the period immediately following the main stage of ice sheet development reveal a major perturbation in the carbon system, represented by the positive d13C excursion known as carbon maximum 6 ("M6"), which has traditionally been interpreted as reflecting increased burial of organic matter and atmospheric pCO2 drawdown. More recently, it has been suggested that the d13C excursion records a negative feedback resulting from the reduction of silicate weathering and an increase in atmospheric pCO2. Here we present high-resolution multi-proxy (alkenone carbon and foraminiferal boron isotope) records of atmospheric carbon dioxide and sea surface temperature across CM6. Similar to previously published records spanning this interval, our records document a world of generally low (~300 ppm) atmospheric pCO2 at a time generally accepted to be much warmer than today. Crucially, they also reveal a pCO2 decrease with associated cooling, which demonstrates that the carbon burial hypothesis for CM6 is feasible and could have acted as a positive feedback on global cooling.

Holocene precipitation change in different monsoon sub-regions (time-slices and transient data) simulated by different global climate models, with links to model results

The recently proposed global monsoon hypothesis interprets monsoon systems as part of one global-scale atmospheric overturning circulation, implying a connection between the regional monsoon systems and an in-phase behaviour of all northern hemispheric monsoons on annual timescales (Trenberth et al., 2000). Whether this concept can be applied to past climates and variability on longer timescales is still under debate, because the monsoon systems exhibit different regional characteristics such as different seasonality (i.e. onset, peak, and withdrawal). To investigate the interconnection of different monsoon systems during the pre-industrial Holocene, five transient global climate model simulations have been analysed with respect to the rainfall trend and variability in different sub-domains of the Afro-Asian monsoon region. Our analysis suggests that on millennial timescales with varying orbital forcing, the monsoons do not behave as a tightly connected global system. According to the models, the Indian and North African monsoons are coupled, showing similar rainfall trend and moderate correlation in rainfall variability in all models. The East Asian monsoon changes independently during the Holocene. The dissimilarities in the seasonality of the monsoon sub-systems lead to a stronger response of the North African and Indian monsoon systems to the Holocene insolation forcing than of the East Asian monsoon and affect the seasonal distribution of Holocene rainfall variations. Within the Indian and North African monsoon domain, precipitation solely changes during the summer months, showing a decreasing Holocene precipitation trend. In the East Asian monsoon region, the precipitation signal is determined by an increasing precipitation trend during spring and a decreasing precipitation change during summer, partly balancing each other. A synthesis of reconstructions and the model results do not reveal an impact of the different seasonality on the timing of the Holocene rainfall optimum in the different sub-monsoon systems. They rather indicate locally inhomogeneous rainfall changes and show, that single palaeo-records should not be used to characterise the rainfall change and monsoon evolution for entire monsoon sub-systems.

Water stable oxygen isotope records from six Antarctic ice core sites for the present and last interglacial periods

We compare the present and last interglacial periods as recorded in Antarctic water stable isotope records now available at various temporal resolutions from six East Antarctic ice cores: Vostok, Taylor Dome, EPICA Dome C (EDC), EPICA Dronning Maud Land (EDML), Dome Fuji and the recent TALDICE ice core from Talos Dome. We first review the different modern site characteristics in terms of ice flow, meteorological conditions, precipitation intermittency and moisture origin, as depicted by meteorological data, atmospheric reanalyses and Lagrangian moisture source diagnostics. These different factors can indeed alter the relationships between temperature and water stable isotopes. Using five records with sufficient resolution on the EDC3 age scale, common features are quantified through principal component analyses. Consistent with instrumental records and atmospheric model results, the ice core data depict rather coherent and homogenous patterns in East Antarctica during the last two interglacials. Across the East Antarctic plateau, regional differences, with respect to the common East Antarctic signal, appear to have similar patterns during the current and last interglacials. We identify two abrupt shifts in isotopic records during the glacial inception at TALDICE and EDML, likely caused by regional sea ice expansion. These regional differences are discussed in terms of moisture origin and in terms of past changes in local elevation histories, which are compared to ice sheet model results. Our results suggest that elevation changes may contribute significantly to inter-site differences. These elevation changes may be underestimated by current ice sheet models

Specific conductivity and hydrogen ions measured on two ice core (NEEM and NGRIP)

A stratigraphy-based chronology for the North Greenland Eemian Ice Drilling (NEEM) ice core has been derived by transferring the annual layer counted Greenland Ice Core Chronology 2005 (GICC05) and its model extension (GICC05modelext) from the NGRIP core to the NEEM core using 787 match points of mainly volcanic origin identified in the electrical conductivity measurement (ECM) and dielectrical profiling (DEP) records. Tephra horizons found in both the NEEM and NGRIP ice cores are used to test the matching based on ECM and DEP and provide five additional horizons used for the timescale transfer. A thinning function reflecting the accumulated strain along the core has been determined using a Dansgaard-Johnsen flow model and an isotope-dependent accumulation rate parameterization. Flow parameters are determined from Monte Carlo analysis constrained by the observed depth-age horizons. In order to construct a chronology for the gas phase, the ice age-gas age difference (Delta age) has been reconstructed using a coupled firn densification-heat diffusion model. Temperature and accumulation inputs to the Delta age model, initially derived from the water isotope proxies, have been adjusted to optimize the fit to timing constraints from d15N of nitrogen and high-resolution methane data during the abrupt onset of Greenland interstadials. The ice and gas chronologies and the corresponding thinning function represent the first chronology for the NEEM core, named GICC05modelext-NEEM-1. Based on both the flow and firn modelling results, the accumulation history for the NEEM site has been reconstructed. Together, the timescale and accumulation reconstruction provide the necessary basis for further analysis of the records from NEEM.

Constraints on the magnitude and patterns of ocean cooling at the Last Glacial Maximum - Supplementary material

Observation-based reconstructions of sea surface temperature from relatively stable periods in the past, such as the Last Glacial Maximum, represent an important means of constraining climate sensitivity and evaluating model simulations. The first quantitative global reconstruction of sea surface temperatures during the Last Glacial Maximum was developed by the Climate Long-Range Investigation, Mapping and Prediction (CLIMAP) project in the 1970s and 1980s. Since that time, several shortcomings of that earlier effort have become apparent. Here we present an updated synthesis of sea surface temperatures during the Last Glacial Maximum, rigorously defined as the period between 23 and 19 thousand years before present, from the Multiproxy Approach for the Reconstruction of the Glacial Ocean Surface (MARGO) project. We integrate microfossil and geochemical reconstructions of surface temperatures and include assessments of the reliability of individual records. Our reconstruction reveals the presence of large longitudinal gradients in sea surface temperature in all of the ocean basins, in contrast to the simulations of the Last Glacial Maximum climate available at present.

Geochemistry of reduction haloes in ODP Holes 135-834A and 135-837A

The sediment column overlying basement in the Lau Basin consists of a sequence of volcaniclastic turbidites interbedded with hemipelagic clayey nannofossil mixed sediments, overlain in turn by a sequence of hemipelagic clayey nannofossil oozes containing sporadic calcareous turbidites. The clayey nannofossil oozes and mixed sediments are pervasively stained by hydrothermally derived iron and manganese oxyhydroxides. Sharply defined, lighter colored bands occur in the hemipelagic sediments, immediately beneath some (but by no means all) volcaniclastic and calcareous turbidites. These are identified as reduction haloes, of a type previously identified in quite different turbidite/pelagic sequences. The haloes are attributed to the burial of labile surficial Corg by turbidites, followed by the remineralization of this Corg with Mn and Fe oxyhydroxides as electron acceptors. The resultant characteristic Mn and Fe concentration/depth profiles are described, and a model is proposed for their development. The color alteration of the halo is ascribed to the removal of Mn oxyhydroxides, because, although the Fe content fluctuates through the haloes, this does not appear to affect their color. Other elements (Co, Cu, and Ni) are also at low concentration levels in the haloes like Mn, consistent with remobilization and migration out of the halo section, although the profile shapes are not identical with those of Mn. The behavior of V is distinctive in that it appears to have migrated into the haloes to be enriched there. Haloes are unlikely to form if turbidite emplacement is erosive and removes the near-surface layer, which generally is the most fluid part of the sediment and contains the highest levels of reactive Corg to drive the reduction process. Conversely, the presence of a halo implies that emplacement of the overlying turbidite did not significantly erode the pre-existing sediment/water interface.

Digital elevation model (DEM) of the ice sheet in western Dronning Maud Land, Antarctica

In this paper, a new digital elevation model (DEM) is derived for the ice sheet in western Dronning Maud Land, Antarctica. It is based on differential interferometric synthetic aperture radar (SAR) from the European Remote Sensing 1/2 (ERS-1/2) satellites, in combination with ICESat's Geoscience Laser Altimeter System (GLAS). A DEM mosaic is compiled out of 116 scenes from the ERS-1 ice phase in 1994 and the ERS-1/2 tandem mission between 1996 and 1997 with the GLAS data acquired in 2003 that served as ground control. Using three different SAR processors, uncertainties in phase stability and baseline model, resulting in height errors of up to 20 m, are exemplified. Atmospheric influences at the same order of magnitude are demonstrated, and corresponding scenes are excluded. For validation of the DEM mosaic, covering an area of about 130,000 km**2 on a 50-m grid, independent ICESat heights (2004-2007), ground-based kinematic GPS (2005), and airborne laser scanner data (ALS, 2007) are used. Excluding small areas with low phase coherence, the DEM differs in mean and standard deviation by 0.5 +/- 10.1, 1.1 +/- 6.4, and 3.1 +/- 4.0 m from ICESat, GPS, and ALS, respectively. The excluded data points may deviate by more than 50 m. In order to suppress the spatially variable noise below a 5-m threshold, 18% of the DEM area is selectively averaged to a final product at varying horizontal spatial resolution. Apart from mountainous areas, the new DEM outperforms other currently available DEMs and may serve as a benchmark for future elevation models such as from the TanDEM-X mission to spatially monitor ice sheet elevation.

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.

Lithogeochemical and isotope (Rb-Sr, Sm-Nd, Pb-Pb whole rock and Lu-Hf zircon) composition of samples from the Shackleton Range, East Antarctica

Three distinct, spatially separated crustal terranes have been recognised in the Shackleton Range, East Antarctica: the Southern, Eastern and Northern Terranes. Mafic gneisses from the Southern Terrane provide geochemical evidence for a within-plate, probably back-arc origin of their protoliths. A plume-distal ridge origin in an incipient ocean basin is the favoured interpretation for the emplacement site of these rocks at c. 1850 Ma, which, together with a few ocean island basalts, were subsequently incorporated into an accretionary continental arc/supra-subduction zone tectonic setting. Magmatic underplating resulted in partial melting of the lower crust, which caused high-temperature granulite-facies metamorphism in the Southern Terrane at c. 1710-1680 Ma. Mafic and felsic gneisses there are characterised by isotopically depleted, positive Nd and Hf initials and model ages between 2100 and 2000 Ma. They may be explained as juvenile additions to the crust towards the end of the Palaeoproterozoic. These juvenile rocks occur in a narrow, c. 150 km long E-W trending belt, inferred to trace a suture that is associated with a large Palaeoproterozoic accretionary orogenic system. The Southern Terrane contains many features that are similar to the Australo-Antarctic Mawson Continent and may be its furthermost extension into East Antarctica. The Eastern Terrane is characterised by metagranitoids that formed in a continental volcanic arc setting during a late Mesoproterozoic orogeny at c. 1060 Ma. Subsequently, the rocks experienced high-temperature metamorphism during Pan-African collisional tectonics at 600 Ma. Isotopically depleted zircon grains yielded Hf model ages of 1600-1400 Ma, which are identical to Nd model ages obtained from juvenile metagranitoids. Most likely, these rocks trace the suture related to the amalgamation of the Indo-Antarctic and West Gondwana continental blocks at ~600 Ma. The Eastern Terrane is interpreted as the southernmost extension of the Pan-African Mozambique/Maud Belt in East Antarctica and, based on Hf isotope data, may also represent a link to the Ellsworth-Whitmore Mountains block in West Antarctica and the Namaqua-Natal Province of southern Africa. Geochemical evidence indicates that the majority of the protoliths of the mafic gneisses in the Northern Terrane formed as oceanic island basalts in a within-plate setting. Subsequently the rocks were incorporated into a subduction zone environment and, finally, accreted to a continental margin during Pan-African collisional tectonics. Felsic gneisses there provide evidence for a within-plate and volcanic arc/collisional origin. Emplacement of granitoids occurred at c. 530 Ma and high-temperature, high-pressure metamorphism took place at 510-500 Ma. Enriched Hf and Nd initials and Palaeoproterozoic model ages for most samples indicate that no juvenile material was added to the crust of the Northern Terrane during the Pan-African Orogeny but recycling of older crust or mixing of crustal components of different age must have occurred. Isotopically depleted mafic gneisses, which are spatially associated with eclogite-facies pyroxenites, yielded late Mesoproterozoic Nd model ages. These rocks occur in a narrow, at least 100 km long, E-W trending belt that separates alkaline ocean island metabasalts and within-plate metagranitoids from volcanic arc metabasalts and volcanic arc/syn-collisional metagranitoids in the Northern Terrane. This belt is interpreted to trace the late Neoproterozoic/early Cambrian Pan-African collisional suture between the Australo-Antarctic and the combined Indo-Antarctic/West Gondwana continental blocks that formed during the final amalgamation of Gondwana.

Helium isotope ratios and sedimentation rate models of ODP holes

A continuous age model for the brief climate excursion at the Paleocene-Eocene boundary has been constructed by assuming a constant flux of extraterrestrial 3He (3He[ET]) to the seafloor. 3He[ET] measurements from ODP Site 690 provide quantitative evidence for the rapid onset (and short duration (<120 kyr) of global warming and of the associated disturbance to the Earth's surficial carbon budget at this time. These observations support astronomically calibrated age models indicating extremely rapid release of isotopically light carbon, possibly from seafloor methane hydrate, as the proximal cause of the event. However, the 3He[ET] technique indicates a previously unrecognized and extreme increase in sedimentation rate coincident with the return of climate proxies to pre-event values. The 3He[ET]-based age model thus suggests a far more rapid recovery from the climatic perturbation than previously proposed or predicted on the basis of the modern carbon cycle, and so may indicate additional or accelerated mechanisms of carbon removal from the ocean-atmosphere system during this period. 3He[ET] was also measured at ODP Site 1051 to test the validity of the Site 690 chronology. Comparison of these data sets seems to require removal of several tens of kyr of sediment within the climatic excursion at Site 1051, an observation consistent with sediment structures and previous age modeling efforts. The Site 1051 age model shows a ~30 kyr period in which climate proxies return toward pre-event values, after which they remain invariant for ~80 kyr. If this rise represents the recovery interval identified at Site 690, then the 3HeET-based age models of the two sites are in good agreement. However, alternative interpretations are possible, and work on less disrupted sites is required to evaluate the reliability of the proposed new chronology of the climate excursion. Regardless of these details, this work shows that the 3HeET technique can provide useful independent evidence for the development and testing of astronomically calibrated age models.