A tephrostratigraphic record for the last glacial-interglacial cycle from Lake Ohrid, Albania and Macedonia

ABSTRACT: Here we present a tephrostratigraphic record (core Co1202) recovered from the northeastern part of Lake Ohrid (Republics of Macedonia and Albania) reaching back to Marine Isotope Stage (MIS) 6. Overall ten horizons (OT0702-1 to OT0702-10) containing volcanic tephra have been recognised throughout the 14.94m long sediment succession. Four tephra layers were visible at macroscopic inspection (OT0702-4, OT0702-6, OT0702-8 and OT0702-9), while the remaining six are cryptotephras (OT0702-1, OT0702-2, OT0702-3, OT0702-5, OT0702-7 and OT0702-10) identified from peaks in K, Zr and Sr intensities, magnetic susceptibility measurements, and washing and sieving of the sediments. Glass shards of tephra layers and cryptotephras were analysed with respect to their major element composition, and correlated to explosive eruptions of Italian volcanoes. The stratigraphy and the major element composition of tephra layers and cryptotephras allowed the correlation of OT0702-1 to AD 472 or AD 512 eruptions of Somma-Vesuvius, OT0702-2 to the FL eruption of Mount Etna, OT0702-3 to the Mercato from Somma-Vesuvius, OT0702-4 to SMP1-e/Y-3 eruption from the Campi Flegrei caldera, OT0702-5 to the Codola eruption (Somma-Vesuvius or Campi Flegrei), OT0702-6 to the Campanian Ignimbrite/Y-5 from the Campi Flegrei caldera, OT0702- 7 to the Green Tuff/Y-6 eruption from Pantelleria Island, OT0702-8 to the X-5 eruption probably originating from the Campi Flegrei caldera, OT0702-9 to the X-6 eruption of generic Campanian origin, and OT0702-10 to the P-11 eruption from Pantelleria Island. The fairly well-known ages of these tephra layers and parent eruptions provide new data on the dispersal and deposition of these tephras and, furthermore, allow the establishment of a chronological framework for core Co1202 for a first interpretation of major sedimentological changes.

Lipid biomarkers in Holocene and glacial sediments from ancient Lake Ohrid (Macedonia, Albania)

Abstract. Organic matter preserved in Lake Ohrid sediments originates from aquatic and terrestrial sources. Its variable composition reflects climate-controlled changes in the lake basin’s hydrology and related organic matter export, i.e. changes in primary productivity, terrestrial plant matter input and soil erosion. Here, we present first results from lipid biomarker investigations of Lake Ohrid sediments from two near-shore settings: site Lz1120 near the southern shore, with low-lying lands nearby and probably influenced by river discharge, and site Co1202 which is close to the steep eastern slopes. Variable proportions of terrestrial n-alkanoic acids and n-alkanols as well as compositional changes of !- hydroxy acids document differences in soil organic matter supply between the sites and during different climate stages (glacial, Holocene, 8.2 ka cooling event). Changes in the vegetation cover are suggested by changes in the dominant chain length of terrestrial n-alkanols. Effective microbial degradation of labile organic matter and in situ contribution of organic matter derived from the microbes themselves are both evident in the sediments. We found evidence for anoxic conditions within the photic zone by detecting epicholestanol and tetrahymanol from sulphur-oxidising phototrophic bacteria and bacterivorous ciliates and for the influence of a settled human community from the occurrence of coprostanol, a biomarker for human and animal faeces (pigs, sheep, goats), in an early Holocene sample. This study illustrates the potential of lipid biomarkers for future environmental reconstructions using one of Europe’s oldest continental climate archives, Lake Ohrid.

Glacier response to the change in atmospheric circulation in the eastern Mediterranean during the Last Glacial Maximum

In this study, we document glacial deposits and reconstruct the glacial history in the Karagöl valley system in the eastern Uludağ in northwestern Turkey based on 42 cosmogenic 10Be exposure ages from boulders and bedrock. Our results suggest the Last Glacial Maximum (LGM) advance prior to 20.4 ± 1.2 ka and at least three re-advances until 18.6 ± 1.2 ka during the global LGM within Marine Isotope Stage-2. In addition, two older advances of unknown age are geomorphologically well constrained, but not dated due to the absence of suitable boulders. Glaciers advanced again two times during the Lateglacial. The older is exposure dated to not later than 15.9 ± 1.1 ka and the younger is attributed to the Younger Dryas (YD) based on field evidence. The timing of the glaciations in the Karagöl valley correlates well with documented archives in the Anatolian and Mediterranean mountains and the Alps. These glacier fluctuations may be explained by the change in the atmospheric circulation pattern during the different phases of North Atlantic Oscillation (NAO) winter indices.

Late Weichselian local ice dome configuration and chronology in Northwestern Svalbard: early thinning, late retreat

The chronology and configuration of the Svalbard Barents Sea Ice Sheet (SBSIS) during the Late Weichselian (LW) are based on few and geographically scattered data. Thus, the timing and configuration of the SBSIS has been a subject of extensive debate. We present provenance data of erratic boulders and cosmogenic 10Be ages of bedrock and boulders from Northwest Spitsbergen (NWS), Svalbard to determine the thickness, configuration and chronology during the LW. We sampled bedrock and boulders of mountain summits and summit slopes, along with erratic boulders from coastal locations around NWS. We suggest that a local ice dome over central NWS during LW drained radially in all directions. Provenance data from erratic boulders from northern coastal lowland Reinsdyrflya suggest northeastward ice flow through Liefdefjorden. 10Be ages of high-elevation erratic boulders in central NWS (687–836 m above sea level) ranging from 18.3 ± 1.3 ka to 21.7 ± 1.4 ka, indicate that the centre of a local ice dome was at least 300 m thicker than at present. 10Be ages of all high-elevation erratics (>400 m above sea level, central and coastal locations) indicate the onset of ice dome thinning at 25–20 ka. 10Be ages from erratic boulders on Reinsdyrflya ranging from 11.1 ± 0.8 ka to 21.4 ± 1.7 ka, indicate an ice cover over the entire Reinsdyrflya during LW and a complete deglaciation prior to the Holocene, but apparently later than the thinning in the mountains. Lack of moraine deposits, but the preservation of beach terraces, suggest that the ice covering this peninsula possibly was cold-based and that Reinsdyrflya was part of an inter ice-stream area covered by slow-flowing ice, as opposed to the adjacent fjord, which possibly was filled by a fast-flowing ice stream. Despite the early thinning of the ice sheet (25–20 ka) we find a later timing of deglaciation of the fjords and the distal lowlands. Several bedrock samples (10Be) from vertical transects in the central mountains of NWS pre-date the LW, and suggest either ice free or pervasive cold-based ice conditions. Our reconstruction is aligned with the previously suggested hypothesis that a complex multi-dome ice-sheet-configuration occupied Svalbard and the Barents Sea during LW, with numerous drainage basins feeding fast ice streams, separated by slow flowing, possibly cold-based, inter ice-stream areas.

Water versus ice: The competing roles of modern climate and Pleistocene glacial erosion in the Central Alps of Switzerland

Recent studies have identified relationships between landscape form, erosion and climate in regions of landscape rejuvenation, associated with increased denudation. Most of these landscapes are located in non-glaciated mountain ranges and are characterized by transient geomorphic features. The landscapes of the Swiss Alps are likewise in a transient geomorphic state as seen by multiple knickzones. In this mountain belt, the transient state has been related to erosional effects during the Late Glacial Maximum (LGM). Here, we focus on the catchment scale and categorize hillslopes based on erosional mechanisms, landscape form and landcover. We then explore relationships of these variables to precipitation and extent of LGM glaciers to disentangle modern versus palaeo controls on the modern shape of the Alpine landscape. We find that in grasslands, the downslope flux of material mainly involves unconsolidated material through hillslope creep, testifying a transport-limited erosional regime. Alternatively, strength-limited hillslopes, where erosion is driven by bedrock failure, are covered by forests and/or expose bedrock, and they display oversteepened hillslopes and channels. There, hillslope gradients and relief are more closely correlated with LGM ice occurrence than with precipitation or the erodibility of the underlying bedrock. We relate the spatial occurrence of the transport- and strength-limited process domains to the erosive effects of LGM glaciers. In particular, strength-limited, rock dominated basins are situated above the equilibrium line altitude (ELA) of the LGM, reflecting the ability of glaciers to scour the landscape beyond threshold slope conditions. In contrast, transport-limited, soil-mantled landscapes are common below the ELA. Hillslopes covered by forests occupy the elevations around the ELA and are constrained by the tree line. We conclude that the current erosional forces at work in the Central Alps are still responding to LGM glaciation, and that the modern climate has not yet impacted on the modern landscape.

Transient simulations of the carbon and nitrogen dynamics in northern peatlands: from the Last Glacial Maximum to the 21st century

The development of northern high-latitude peatlands played an important role in the carbon (C) balance of the land biosphere since the Last Glacial Maximum (LGM). At present, carbon storage in northern peatlands is substantial and estimated to be 500 ± 100 Pg C (1 Pg C = 1015 g C). Here, we develop and apply a peatland module embedded in a dynamic global vegetation and land surface process model (LPX-Bern 1.0). The peatland module features a dynamic nitrogen cycle, a dynamic C transfer between peatland acrotelm (upper oxic layer) and catotelm (deep anoxic layer), hydrology- and temperature-dependent respiration rates, and peatland specific plant functional types. Nitrogen limitation down-regulates average modern net primary productivity over peatlands by about half. Decadal acrotelm-to-catotelm C fluxes vary between −20 and +50 g C m−2 yr−1 over the Holocene. Key model parameters are calibrated with reconstructed peat accumulation rates from peat-core data. The model reproduces the major features of the peat core data and of the observation-based modern circumpolar soil carbon distribution. Results from a set of simulations for possible evolutions of northern peat development and areal extent show that soil C stocks in modern peatlands increased by 365–550 Pg C since the LGM, of which 175–272 Pg C accumulated between 11 and 5 kyr BP. Furthermore, our simulations suggest a persistent C sequestration rate of 35–50 Pg C per 1000 yr in present-day peatlands under current climate conditions, and that this C sink could either sustain or turn towards a source by 2100 AD depending on climate trajectories as projected for different representative greenhouse gas concentration pathways.

The response of atmospheric nitrous oxide to climate variations during the last glacial period

Detailed insight into natural variations of the greenhouse gas nitrous oxide (N2O) in response to changes in the Earth's climate system is provided by new measurements along the ice core of the North Greenland Ice Core Project (NGRIP). The presented record reaches from the early Holocene back into the previous interglacial with a mean time resolution of about 75 years. Between 11 and 120 kyr BP, atmospheric N2O concentrations react substantially to the last glacial-interglacial transition (Termination 1) and millennial time scale climate variations of the last glacial period. For long-lasting Dansgaard/Oeschger (DO) events, the N2O increase precedes Greenland temperature change by several hundred years with an increase rate of about 0.8-1.3 ppbv/century, which accelerates to about 3.8-10.7 ppbv/century at the time of the rapid warming in Greenland. Within each bundle of DO events, the new record further reveals particularly low N2O concentrations at the approximate time of Heinrich events. This suggests that the response of marine and/or terrestrial N2O emissions on a global scale are different for stadials with and without Heinrich events.

A reconstruction of atmospheric carbon dioxide and its stable carbon isotopic composition from the penultimate glacial maximum to the last glacial inception

The reconstruction of the stable carbon isotope evolution in atmospheric CO2 (δ13Catm), as archived in Antarctic ice cores, bears the potential to disentangle the contributions of the different carbon cycle fluxes causing past CO2 variations. Here we present a new record of δ13Catm before, during and after the Marine Isotope Stage 5.5 (155 000 to 105 000 yr BP). The dataset is archived on the data repository PANGEA® (www.pangea.de) under 10.1594/PANGAEA.817041. The record was derived with a well established sublimation method using ice from the EPICA Dome C (EDC) and the Talos Dome ice cores in East Antarctica. We find a 0.4‰ shift to heavier values between the mean δ13Catm level in the Penultimate (~ 140 000 yr BP) and Last Glacial Maximum (~ 22 000 yr BP), which can be explained by either (i) changes in the isotopic composition or (ii) intensity of the carbon input fluxes to the combined ocean/atmosphere carbon reservoir or (iii) by long-term peat buildup. Our isotopic data suggest that the carbon cycle evolution along Termination II and the subsequent interglacial was controlled by essentially the same processes as during the last 24 000 yr, but with different phasing and magnitudes. Furthermore, a 5000 yr lag in the CO2 decline relative to EDC temperatures is confirmed during the glacial inception at the end of MIS5.5 (120 000 yr BP). Based on our isotopic data this lag can be explained by terrestrial carbon release and carbonate compensation.

An optimized multi-proxy, multi-site Antarctic ice and gas orbital chronology (AICC2012): 120--800 ka

An accurate and coherent chronological framework is essential for the interpretation of climatic and environmental records obtained from deep polar ice cores. Until now, one common ice core age scale had been developed based on an inverse dating method (Datice), combining glaciological modelling with absolute and stratigraphic markers between 4 ice cores covering the last 50 ka (thousands of years before present) (Lemieux-Dudon et al., 2010). In this paper, together with the companion paper of Veres et al. (2013), we present an extension of this work back to 800 ka for the NGRIP, TALDICE, EDML, Vostok and EDC ice cores using an improved version of the Datice tool. The AICC2012 (Antarctic Ice Core Chronology 2012) chronology includes numerous new gas and ice stratigraphic links as well as improved evaluation of background and associated variance scenarios. This paper concentrates on the long timescales between 120–800 ka. In this framework, new measurements of δ18Oatm over Marine Isotope Stage (MIS) 11–12 on EDC and a complete δ18Oatm record of the TALDICE ice cores permit us to derive additional orbital gas age constraints. The coherency of the different orbitally deduced ages (from δ18Oatm, δO2/N2 and air content) has been verified before implementation in AICC2012. The new chronology is now independent of other archives and shows only small differences, most of the time within the original uncertainty range calculated by Datice, when compared with the previous ice core reference age scale EDC3, the Dome F chronology, or using a comparison between speleothems and methane. For instance, the largest deviation between AICC2012 and EDC3 (5.4 ka) is obtained around MIS 12. Despite significant modifications of the chronological constraints around MIS 5, now independent of speleothem records in AICC2012, the date of Termination II is very close to the EDC3 one.

High-resolution glacial and deglacial record of atmospheric methane by continuous-flow and laser spectrometer analysis along the NEEM ice core

The Greenland NEEM (North Greenland Eemian Ice Drilling) operation in 2010 provided the first opportunity to combine trace-gas measurements by laser spectroscopic instruments and continuous-flow analysis along a freshly drilled ice core in a field-based setting. We present the resulting atmospheric methane (CH4) record covering the time period from 107.7 to 9.5 ka b2k (thousand years before 2000 AD). Companion discrete CH4 measurements are required to transfer the laser spectroscopic data from a relative to an absolute scale. However, even on a relative scale, the high-resolution CH4 data set significantly improves our knowledge of past atmospheric methane concentration changes. New significant sub-millennial-scale features appear during interstadials and stadials, generally associated with similar changes in water isotopic ratios of the ice, a proxy for local temperature. In addition to the midpoint of Dansgaard–Oeschger (D/O) CH4 transitions usually used for cross-dating, sharp definition of the start and end of these events brings precise depth markers (with ±20 cm uncertainty) for further cross-dating with other palaeo- or ice core records, e.g. speleothems. The method also provides an estimate of CH4 rates of change. The onsets of D/O events in the methane signal show a more rapid rate of change than their endings. The rate of CH4 increase associated with the onsets of D/O events progressively declines from 1.7 to 0.6 ppbv yr−1 in the course of marine isotope stage 3. The largest observed rate of increase takes place at the onset of D/O event #21 and reaches 2.5 ppbv yr−1.

Coupled W–Os–Pt isotope systematics in IVB iron meteorites: In situ neutron dosimetry for W isotope chronology

Tungsten isotope compositions of magmatic iron meteorites yield ages of differentiation that are within ±2 Ma of the formation of CAIs, with the exception of IVB irons that plot to systematically less radiogenic compositions yielding erroneously old ages. Secondary neutron capture due to galactic cosmic ray (GCR) irradiation is known to lower the ε182W of iron meteorites, adequate correction of which requires a measure of neutron dosage which has not been available, thus far. The W, Os and Pt isotope systematics of 12 of the 13 known IVB iron meteorites were determined by MC-ICP-MS (W, Os, Pt) and TIMS (Os). On the same dissolutions that yield precise ε182W, stable Os and Pt isotopes were determined as in situ neutron dosimeters for empirical correction of the ubiquitous cosmic-ray induced burn-out of 182W in iron meteorites. The W isotope data reveal a main cluster with ε182W of ∼−3.6, but a much larger range than observed in previous studies including irons (Weaver Mountains and Warburton Range) that show essentially no cosmogenic effect on their ε182W. The IVB data exhibits resolvable negative anomalies in ε189Os (−0.6ε) and complementary ε190Os anomalies (+0.4ε) in Tlacotepec due to neutron capture on 189Os which has approximately the same neutron capture cross section as 182W, and captures neutrons to produce 190Os. The least irradiated IVB iron, Warburton Range, has ε189Os and ε190Os identical to terrestrial values. Similarly, Pt isotopes, which are presented as ε192Pt, ε194Pt and ε196Pt range from +4.4ε to +53ε, +1.54ε to −0.32ε and +0.73ε to −0.20ε, respectively, also identify Tlacotepec and Dumont as the most GCR-damaged samples. In W–Os and W–Pt isotope space, the correlated isotope data back-project toward a 0-epsilon value of ε192Pt, ε189Os and ε190Os from which a pre-GCR irradiation ε182W of −3.42±0.09 (2σ) is derived. This pre-GCR irradiation ε182W is within uncertainty of the currently accepted CAI initial ε182W. The Pt and Os isotope correlations in the IVB irons are in good agreement with a nuclear model for spherical irons undergoing GCR spallation, although this model over-predicts the change of ε182W by ∼2×, indicating a need for better W neutron capture cross section determinations. A nucleosynthetic effect in ε184W in these irons of −0.14±0.08 is confirmed, consistent with the presence of Mo and Ru isotope anomalies in IVB irons. The lack of a non-GCR Os isotope anomaly in these irons requires more complex explanations for the production of W, Ru and Mo anomalies than nebular heterogeneity in the distribution of s-process to r-process nuclides.

Neutron capture on Pt isotopes in iron meteorites and the Hf–W chronology of core formation in planetesimals

The short-lived 182Hf–182W isotope system can provide powerful constraints on the timescales of planetary core formation, but its application to iron meteorites is hampered by neutron capture reactions on W isotopes resulting from exposure to galactic cosmic rays. Here we show that Pt isotopes in magmatic iron meteorites are also affected by capture of (epi)thermal neutrons and that the Pt isotope variations are correlated with variations in 182W/184W. This makes Pt isotopes a sensitive neutron dosimeter for correcting cosmic ray-induced W isotope shifts. The pre-exposure 182W/184W derived from the Pt–W isotope correlations of the IID, IVA and IVB iron meteorites are higher than most previous estimates and are more radiogenic than the initial 182W/184W of Ca–Al-rich inclusions (CAI). The Hf–W model ages for core formation range from +1.6±1.0 million years (Ma; for the IVA irons) to +2.7±1.3 Ma after CAI formation (for the IID irons), indicating that there was a time gap of at least ∼1 Ma between CAI formation and metal segregation in the parent bodies of some iron meteorites. From the Hf–W ages a time limit of <1.5–2 Ma after CAI formation can be inferred for the accretion of the IID, IVA and IVB iron meteorite parent bodies, consistent with earlier conclusions that the accretion of differentiated planetesimals predated that of most chondrite parent bodies.

Glacial forcing of central Indonesian hydroclimate since 60,000 y B.P.

The Indo-Pacific warm pool houses the largest zone of deep atmospheric convection on Earth and plays a critical role in global climate variations. Despite the region’s importance, changes in Indo-Pacific hydroclimate on orbital timescales remain poorly constrained. Here we present high-resolution geochemical records of surface runoff and vegetation from sediment cores fromLake Towuti, on the island of Sulawesi in central Indonesia, that continuously span the past 60,000 y.We show that wet conditions and rainforest ecosystems on Sulawesi present during marine isotope stage 3 (MIS3) and the Holocene were interrupted by severe drying between ∼33,000 and 16,000 y B.P. when Northern Hemisphere ice sheets expanded and global temperatures cooled. Our record reveals little direct influence of precessional orbital forcing on regional climate, and the similarity between MIS3 and Holocene climates observed in Lake Towuti suggests that exposure of the Sunda Shelf has a weaker influence on regional hydroclimate and terrestrial ecosystems than suggested previously. We infer that hydrological variability in this part of Indonesia varies strongly in response to high-latitude climate forcing, likely through reorganizations of the monsoons and the position of the intertropical convergence zone. These findings suggest an important role for the tropical western Pacific in amplifying glacial–interglacial climate variability.