Temperature reconstruction from 10 to 120 kyr b2k from the NGRIP ice core

In order to reconstruct the temperature of the North Greenland Ice Core Project (NGRIP) site, new measurements of δ15N have been performed covering the time period from the beginning of the Holocene to Dansgaard–Oeschger (DO) event 8. Together with previously measured and mostly published δ15N data, we present for the first time a NGRIP temperature reconstruction for the whole last glacial period from 10 to 120 kyr b2k (thousand years before 2000 AD) including every DO event based on δ15N isotope measurements combined with a firn densification and heat diffusion model. The detected temperature rises at the onset of DO events range from 5 °C (DO 25) up to 16.5 °C (DO 11) with an uncertainty of ±3 °C. To bring measured and modelled data into agreement, we had to reduce the accumulation rate given by the NGRIP ss09sea06bm timescale in some periods by 30 to 35%, especially during the last glacial maximum. A comparison between reconstructed temperature and δ18Oice data confirms that the isotopic composition of the stadial was strongly influenced by seasonality. We evidence an anticorrelation between the variations of the δ18Oice sensitivity to temperature (referred to as α) and obliquity in agreement with a simple Rayleigh distillation model. Finally, we suggest that α might be influenced by the Northern Hemisphere ice sheet volume.

The deglaciation history of the Simplon region (southern Swiss Alps) constrained by 10Be exposure dating of ice-molded bedrock surfaces

The deglaciation history of the Swiss Alps after the Last Glacial Maximum involved the decay of several ice domes and the subsequent disintegration of valley glaciers at high altitude. Here we use bedrock exposure dating to reconstruct the temporal and spatial pattern of ice retreat at the Simplon Pass (altitude: ∼2000 m) located 40 km southwest of the ‘Rhône ice dome’. Eleven 10Be exposure ages from glacially polished quartz veins and ice-molded bedrock surfaces cluster tightly between 13.5 ± 0.6 ka and 15.4 ± 0.6 ka (internal errors) indicating that the Simplon Pass depression became ice-free at 14.1 ± 0.4 ka (external error of mean age). This age constraint is interpreted to record the melting of the high valley glaciers in the Simplon Pass region during the warm Bølling–Allerød interstadial shortly after the Oldest Dryas stadial. Two bedrock samples collected a few hundred meters above the pass depression yield older 10Be ages of 17.8 ± 0.6 ka and 18.0 ± 0.6 ka. These ages likely reflect the initial downwasting of the Rhône ice dome and the termination of the ice transfluence from the ice dome across the Simplon Pass toward the southern foreland. There, the retreat of the piedmont glacier in Val d’Ossola was roughly synchronous with the decay of the Rhône ice dome in the interior of the mountain belt, as shown by 10Be ages of 17.7 ± 0.9 ka and 16.1 ± 0.6 ka for a whaleback at ∼500 m elevation near Montecrestese in northern Italy. In combination with well-dated paleoclimate records derived from lake sediments, our new age data suggest that during the deglaciation of the European Alps the decay of ice domes was approximately synchronous with the retreat of piedmont glaciers in the foreland and was followed by the melting of high-altitude valley glaciers after the transition from the Oldest Dryas to the Bølling–Allerød, when mean annual temperatures rose rapidly by ∼3 °C.

The first major incision of the Swiss Deckenschotter landscape

The Swiss Deckenschotter (“cover gravels”) is the oldest Quaternary units in the northern Swiss Alpine Foreland. They are a succession of glaciofluvial gravel layers intercalated with glacial and/or overbank deposits. This lithostratigraphic sequence is called Deckenschotter because it “covers” Molasse or Mesozoic bedrock and forms mesa-type hill-tops. Deckenschotter occurs both within and beyond the extent of the Last Glacial Maximum glaciers. The Swiss Deckenschotter consist of two sub-units: Höhere (Higher) and Tiefere (Lower) Deckenschotter. Although the Höhere Deckenschotter sub-unit (HDS) is topographically higher than the Tiefere Deckenschotter, it is older. The only available age for the Swiss Deckenschotter is 2.5–1.8 Ma based on mammal remains found in HDS at the Irchel site. In this study, we present an exposure age for the topographically lowest HDS, calculated from a cosmogenic 10Be depth-profile. Our results show that the first phase of the Deckenschotter glaciations in the Swiss Alps terminated at least 1,020+80−120 ka ago, which is indicated by a significant fluvial incision. This line of evidence seems to be close to synchronous with the beginning of the Mid-Pleistocene Revolution, when the frequency of the glacial-interglacial cyclicity changed from 41 to 100 ka and the amplitude from low to high, between marine isotope stages 23 and 22.

Late Quaternary history of North Eurasian Norway spruce ( Picea abies ) and Siberian spruce ( Picea obovata ) inferred from macrofossils, pollen and cytoplasmic DNA variation

Aim We used combined palaeobotanical and genetic data to assess whether Norway spruce (Picea abies) and Siberian spruce (Picea obovata), two major components of the Eurasian boreal forests, occupied separate glacial refugia, and to test previous hypotheses on their distinction, geographical delimitation and introgression. Location The range of Norway spruce in northern Europe and Siberian spruce in northern Asia. Methods Pollen data and recently compiled macrofossil records were summarized for the Last Glacial Maximum (LGM), late glacial and Holocene. Genetic variation was assessed in 50 populations using one maternally (mitochondrial nad1) and one paternally (chloroplast trnT–trnL) inherited marker and analysed using spatial analyses of molecular variance (SAMOVA). Results Macrofossils showed that spruce was present in both northern Europe and Siberia at the LGM. Congruent macrofossil and pollen data from the late glacial suggested widespread expansions of spruce in the East European Plain, West Siberian Plain, southern Siberian mountains and the Baikal region. Colonization was largely completed during the early Holocene, except in the formerly glaciated area of northern Europe. Both DNA markers distinguished two highly differentiated groups that correspond to Norway spruce and Siberian spruce and coincide spatially with separate LGM spruce occurrences. The division of the mtDNA variation was geographically well defined and occurred to the east of the Ural Mountains along the Ob River, whereas the cpDNA variation showed widespread admixture. Genetic diversity of both DNA markers was higher in western than in eastern populations. Main conclusions North Eurasian Norway spruce and Siberian spruce are genetically distinct and occupied separate LGM refugia, Norway spruce on the East European Plain and Siberian spruce in southern Siberia, where they were already widespread during the late glacial. They came into contact in the basin of the Ob River and probably hybridized. The lower genetic diversity in the eastern populations may indicate that Siberian spruce suffered more from past climatic fluctuations than Norway spruce.

Ignoring heterozygous sites biases phylogenomic estimates of divergence times: implications for the evolutionary history of Microtus voles

Phylogenetic reconstruction of the evolutionary history of closely related organisms may be difficult because of the presence of unsorted lineages and of a relatively high proportion of heterozygous sites that are usually not handled well by phylogenetic programs. Genomic data may provide enough fixed polymorphisms to resolve phylogenetic trees, but the diploid nature of sequence data remains analytically challenging. Here, we performed a phylogenomic reconstruction of the evolutionary history of the common vole (Microtus arvalis) with a focus on the influence of heterozygosity on the estimation of intraspecific divergence times. We used genome-wide sequence information from 15 voles distributed across the European range. We provide a novel approach to integrate heterozygous information in existing phylogenetic programs by repeated random haplotype sampling from sequences with multiple unphased heterozygous sites. We evaluated the impact of the use of full, partial, or no heterozygous information for tree reconstructions on divergence time estimates. All results consistently showed four deep and strongly supported evolutionary lineages in the vole data. These lineages undergoing divergence processes split only at the end or after the last glacial maximum based on calibration with radiocarbon-dated paleontological material. However, the incorporation of information from heterozygous sites had a significant impact on absolute and relative branch length estimations. Ignoring heterozygous information led to an overestimation of divergence times between the evolutionary lineages of M. arvalis. We conclude that the exclusion of heterozygous sites from evolutionary analyses may cause biased and misleading divergence time estimates in closely related taxa.

Long distance dispersal shaped patterns of human genetic diversity in Eurasia

Most previous attempts at reconstructing the past history of human populations did not explicitly take geography into account, or considered very simple scenarios of migration and ignored environmental information. However, it is likely that the Last Glacial Maximum (LGM) affected the demography and the range of many species, including our own. Moreover, long-distance dispersal (LDD) may have been an important component of human migrations, allowing fast colonization of new territories and preserving high levels of genetic diversity. Here, we use a high-quality microsatellite dataset genotyped in 22 populations to estimate the posterior probabilities of several scenarios for the settlement of the Old World by modern humans. We considered models ranging from a simple spatial expansion to others including LDD and a LGM-induced range contraction, as well as Neolithic demographic expansions. We find that scenarios with LDD are much better supported by data than models without LDD. Nevertheless, we show evidence that LDD events to empty habitats were strongly prevented during the settlement of Eurasia. This unexpected absence of LDD ahead of the colonization wave front could have been caused by an Allee effect, either due to intrinsic causes such as an inbreeding depression built during the expansion, or to extrinsic causes such as direct competition with archaic humans. Overall, our results suggest only a relatively limited effect of the LGM-contraction on current patterns of human diversity. This is in clear contrast with the major role of LDD migrations, which have potentially contributed to the intermingled genetic structure of Eurasian populations.

A 22,000 14C year BP sediment and pollen record of climate change from Laguna Miscanti (23°S), northern Chile

Lake sediments and pollen, spores and algae from the high-elevation endorheic Laguna Miscanti (22°45′S, 67°45′W, 4140 m a.s.l., 13.5 km2 water surface, 10 m deep) in the Atacama Desert of northern Chile provide information about abrupt and high amplitude changes in effective moisture. Although the lack of terrestrial organic macrofossils and the presence of a significant 14C reservoir effect make radiocarbon dating of lake sediments very difficult, we propose the following palaeoenvironmental history. An initial shallow freshwater lake (ca. 22,000 14C years BP) disappeared during the extremely dry conditions of the Last Glacial Maximum (LGM; 18,000 14C years BP). That section is devoid of pollen. The late-glacial lake transgression started around 12,000 14C years BP, peaked in two phases between ca. 11,000 and <9000 14C years BP, and terminated around 8000 14C years BP. Effective moisture increased more than three times compared to modern conditions (∼200 mm precipitation), and a relatively dense terrestrial vegetation was established. Very shallow hypersaline lacustrine conditions prevailed during the mid-Holocene until ca. 3600 14C years BP. However, numerous drying and wetting cycles suggest frequent changes in moisture, maybe even individual storms during the mid-Holocene. After several humid spells, modern conditions were reached at ca. 3000 14C years BP. Comparison between limnogeological data and pollen of terrestrial plants suggest century-scale response lags. Relatively constant concentrations of long-distance transported pollen from lowlands east of the Andes suggest similar atmospheric circulation patterns (mainly tropical summer rainfall) throughout the entire period of time. These findings compare favorably with other regional paleoenvironmental data.

Great Barrier Reef coral element/Ca and stable isotope data and U-Th ages from IODP Expedition 325

Tropical south-western Pacific temperatures are of vital importance to the Great Barrier Reef (GBR), but the role of sea surface temperatures (SSTs) in the growth of the GBR since the Last Glacial Maximum remains largely unknown. Here we present records of Sr/Ca and d18O for Last Glacial Maximum and deglacial corals that show a considerably steeper meridional SST gradient than the present day in the central GBR. We find a 1-2 °C larger temperature decrease between 17° and 20°S about 20,000 to 13,000 years ago. The result is best explained by the northward expansion of cooler subtropical waters due to a weakening of the South Pacific gyre and East Australian Current. Our findings indicate that the GBR experienced substantial meridional temperature change during the last deglaciation, and serve to explain anomalous deglacial drying of northeastern Australia. Overall, the GBR developed through significant SST change and may be more resilient than previously thought.

Reconstructed glaciation limits in Eurasia (digitized from map Fig.1)

Based on field investigations in northern Russia and interpretation of offshore seismic data, we have made a preliminary reconstruction of the maximum ice-sheet extent in the Barents and Kara Sea region during the Early/Middle Weichselian and the Late Weichselian. Our investigations indicate that the Barents and Kara ice sheets attained their maximum Weichselian positions in northern Russia prior to 50 000 yr BP, whereas the northeastern flank of the Scandinavian Ice Sheet advanced to a maximum position shortly after 17000 calendar years ago. During the Late Weichselian (25 000-10000 yr BP), much of the Russian Arctic remained ice-free. According to our reconstruction, the extent of the ice sheets in the Barents and Kara Sea region during the Late Weichselian glacial maximum was less than half that of the maximum model which, up to now, has been widely used as a boundary condition for testing and refining General Circulation Models (GCMs). Preliminary numerical-modelling experiments predict Late Weichselian ice sheets which are larger than the ice extent implied for the Kara Sea region from dated geological evidence, suggesting very low precipitation.

Particle concentration and size distribution in the NGRIP ice core

A novel laser microparticle detector used in conjunction with continuous sample melting has provided a more than 1500 m long record of particle concentration and size distribution of the NGRIP ice core, covering continuously the period approximately from 9.5-100 kyr before present; measurements were at 1.65 m depth resolution, corresponding to approximately 35-200 yr. Particle concentration increased by a factor of 100 in the Last Glacial Maximum (LGM) compared to the Preboreal, and sharp variations of concentration occurred synchronously with rapid changes in the delta18O temperature proxy. The lognormal mode µ of the volume distribution shows clear systematic variations with smaller modes during warmer climates and coarser modes during colder periods. We find µ ~ 1.7 µm diameter during LGM and µ ~ 1.3 µm during the Preboreal. On timescales below several 100 years µ and the particle concentration exhibit a certain degree of independence present especially during warm periods, when µ generally is more variable. Using highly simplifying considerations for atmospheric transport and deposition of particles we infer that (1) the observed changes of µ in the ice largely reflect changes in the size of airborne particles above the ice sheet and (2) changes of µ are indicative of changes in long range atmospheric transport time. From the observed size changes we estimate shorter transit times by roughly 25% during LGM compared to the Preboreal. The associated particle concentration increase from more efficient long range transport is estimated to less than one order of magnitude.

Model results of ocean physical field of the reference state and the targets in NetCDF format

In a feasibility study, the potential of proxy data for the temperature and salinity during the Last Glacial Maximum (LGM, about 19 000 to 23 000 years before present) in constraining the strength of the Atlantic meridional overturning circulation (AMOC) with a general ocean circulation model was explored. The proxy data were simulated by drawing data from four different model simulations at the ocean sediment core locations of the Multiproxy Approach for the Reconstruction of the Glacial Ocean surface (MARGO) project, and perturbing these data with realistic noise estimates. The results suggest that our method has the potential to provide estimates of the past strength of the AMOC even from sparse data, but in general, paleo-sea-surface temperature data without additional prior knowledge about the ocean state during the LGM is not adequate to constrain the model. On the one hand, additional data in the deep-ocean and salinity data are shown to be highly important in estimating the LGM circulation. On the other hand, increasing the amount of surface data alone does not appear to be enough for better estimates. Finally, better initial guesses to start the state estimation procedure would greatly improve the performance of the method. Indeed, with a sufficiently good first guess, just the sea-surface temperature data from the MARGO project promise to be sufficient for reliable estimates of the strength of the AMOC.

Geological investigations on sediment core MD07-3128

Glacial millennial-scale paleoceanographic changes in the Southeast Pacific and the adjacent Southern Ocean are poorly known due to the scarcity of well-dated and high resolution sediment records. Here we present new surface water records from sediment core MD07-3128 recovered at 53°S off the Pacific entrance of the Strait of Magellan. The alkenone-derived sea surface temperature (SST) record reveals a very strong warming of ca. 8°C over the last Termination and substantial millennial-scale variability in the glacial section largely consistent with our planktonic foraminifera oxygen isotope (d18O) record of Neogloboquadrina pachyderma (sin.). The timing and structure of the Termination and some of the millennial-scale SST fluctuations are very similar to those observed in the well-dated SST record from ODP Site 1233 (41°S) and the temperature record from Drowning Maud Land Antarctic ice core supporting the hemispheric-wide Antarctic timing of SST changes. However, differences in our new SST record are also found including a long-term warming trend over Marine Isotope Stage (MIS) 3 followed by a cooling toward the Last Glacial Maximum (LGM). We suggest that these differences reflect regional cooling related to the proximal location of the southern Patagonian Ice Sheet and related meltwater supply at least during the LGM consistent with the fact that no longer SST cooling trend is observed in ODP Site 1233 or any SST Chilean record. This proximal ice sheet location is documented by generally higher contents of ice rafted debris (IRD) and tetra-unsaturated alkenones, and a slight trend toward lighter planktonic d18O during late MIS 3 and MIS 2.