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.

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.

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.

High-resolution late-glacial chronology for the Gerzensee lake record (Switzerland): δ18O correlation between a Gerzensee-stack and NGRIP

Oxygen-isotope variations were analyzed on bulk samples of shallow-water lake marl from Gerzensee, Switzerland, in order to evaluate major and minor climatic oscillations during the late-glacial. To highlight the overall signature of the Gerzensee δ18O record, δ18O records of four parallel sediment cores were first correlated by synchronizing major isotope shifts and pollen abundances. Then the records were stacked with a weighting depending on the differing sampling resolution. To develop a precise chronology, the δ18O-stack was then correlated with the NGRIP δ18O record applying a Monte Carlo simulation, relying on the assumption that the shifts in δ18O were climate-driven and synchronous in both archives. The established chronology on the GICC05 time scale is the basis for (1) comparing the δ18O changes recorded in Gerzensee with observed climatic and environmental fluctuations over the whole North Atlantic region, and (2) comparing sedimentological and biological changes during the rapid warming with smaller climatic variations during the Bølling/Allerød period. The δ18O record of Gerzensee is characterized by two major isotope shifts at the onset and at the termination of the Bølling/Allerød warm period, as well as four intervening negative shifts labeled GI-1e2, d, c2, and b, which show a shift of one third to one fourth of the major δ18O shifts at the beginning and end of the Bølling/Allerød. Despite some inconsistency in terminology, these oscillations can be observed in various climatic proxies over wide regions in the North Atlantic region, especially in reconstructed colder temperatures, and they seem to be caused by hemispheric climatic variations.

Response of chironomid assemblages to environmental change during the early Late-glacial at Gerzensee, Switzerland

Prior to ca. 14,660 yr BP, during the early Late-glacial (Oldest Dryas), larval assemblages of Chironomidae (Insecta: Diptera) in Gerzensee, Switzerland, were dominated by cold stenothermic taxa as well as by taxa typical of subalpine lakes today. This was the coldest period of the entire sequence. After ca. 14,660 yr BP, in the Late Glacial Interstadial (Bølling–Allerød), a temperature increase is recorded by a sharp rise in the oxygen-isotope ratio in lake marl and by an increase in the organic-matter content of the sediments. Changes in the chironomid fauna then are consistent with rising temperatures. This warming trend is interrupted between 14,070 and 13,940 yr BP, coinciding with the GI-1d cold oscillation, but the change in the chironomid assemblage is more consistent with a response to increasing lake depth and density of aquatic macrophytes than falling temperature. A rise in cold-adapted chironomid taxa between 13,840 and 13,710 yr BP suggests that summer air temperatures may have declined. Changes in the chironomid assemblage after 13,710 yr BP suggest a decline in submerged macrophytes coupled with a rise in lake productivity and summer temperature, although the latter is not reflected in the oxygen-isotope record. This suggests that there may have been increasing seasonality during this period when summer temperatures were rising, driven by rising summer insolation, and winters becoming cooler, which is largely reflected in the oxygen-isotope record. A decline in thermophilic chironomids and a rise in cold-adapted taxa after 13,180 yr BP suggest a response to cooling at the beginning of the Gerzensee Oscillation.

Vegetation responses to rapid warming and to minor climatic fluctuations during the Late-Glacial Interstadial (GI-1) at Gerzensee (Switzerland)

High-resolution pollen analyses made on the same samples on which the ratios of oxygen isotopes were measured that provided the time scale and a temperature proxy after correlation to NorthGRIP. (1) A primary succession: The vegetation responded to the rapid rise of temperatures around 14,685 yr BP, with a primary succession on a decadal to centennial time scale. The succession between ca 15,600 and 13,000 yr BP included: (1.1.) The replacement of shrub-tundra by woodland of Juniperus and tree birch (around 14,665 yr BP) (1.2.) The response of Juniperus pollen to the shift in oxygen isotopes in less than 20 yr, (1.3.) A sequence of population increases of Hippophaë rhamnoides (ca 14,600 yr BP), Salix spp. (ca 14,600 yr BP), Betula trees (ca.14,480 yr BP), Populus cf. tremula (ca. 14,300 yr BP), and Pinus cf. sylvestris (ca. 13,830 yr BP). (2) Biological processes: Plants responded to the rapid increase of summer temperatures on all organisational levels: (2.1) Individuals may have produced more pollen (e.g. Juniperus); (2.2) Populations increased or decreased (e.g. Juniperus, Betula, later Pinus), and (2.3) Populations changed their biogeographical range and may show migrational lags. (2.4) Plant communities changed in their composition because the species pools changed through immigration and (local) extinction. Some plant communities may have been without modern analogue.These mechanisms require increasing amounts of time. (2.5) Processes on the level of ecosystems, with species interactions, may involve various time scales. Besides competition and facilitation, nitrogen fixation is discussed. (3) The minor fluctuations of temperature during the Late-Glacial Interstadial, which are recorded in δ18O, resulted in only very minor changes in pollen during the Aegelsee Oscillation (Older Dryas biozone, GI-1d) and the Gerzensee Oscillation (GI-1b). (4) Biodiversity: The afforestation at the onset of Bølling coincided with a gradual increase of taxonomic diversity up to the time of the major Pinus expansion.

Response of subfossil Cladocera in Gerzensee (Swiss Plateau) to early Late Glacial environmental change

Sub-fossil Cladocera were studied in a core from Gerzensee (Swiss Plateau) for the late-glacial periods of Oldest Dryas, Bølling, and Allerød. Cladocera assemblages were dominated by cold-tolerant littoral taxa Chydorus sphaericus, Acroperus harpae, Alonella nana, Alona affinis, and Alonella excisa. The rapid warming at the beginning of the Bølling (GI-1e) ca. 14,650 yr before present (BP: before AD 1950) was indicated by an abrupt 2‰ shift in carbonate δ18O and a clear change in pollen assemblages. Cladocera assemblages, in contrast, changed more gradually. C. sphaericus and A. harpae are the most cold-tolerant, and their abundance was highest in the earliest part of the record. Only 150–200 years after the beginning of the Bølling warming we observed an increase in less cold-tolerant A. excisa and A. affinis. The establishment of Alona guttata, A. guttata var. tuberculata, and Pleuroxus unicatus was delayed by ca. 350, 770, and 800 years respectively after the onset of the Bølling. The development of the Cladocera assemblages suggests increasing water temperatures during the Bølling/Allerød, which agrees with the interpretation by von Grafenstein et al. (2013-this issue) that decreasing δ18O values in carbonates in this period reflect increasing summer water temperatures at the sediment–water interface. Other processes also affected the Cladocera community, including the development and diversification of aquatic vegetation favourable for Cladocera. The record is clearly dominated by Chydoridae, as expected for a littoral core. Yet, the planktonic Eubosmina-group occurred throughout the core, with the exception of a period at ca. 13,760–13,420 yr BP. Lake levels reconstructed for this period are relatively low, indicating that the littoral location might have become too shallow for Eubosmina in that period.

The contribution of glacial archaeology to our knowledge of alpine passages - L'apport de l'archéologie des glaciers dans la connaissance des voies de passage alpins

Die europäischen Alpen sind in den Hochlagen trotz einem immensen Rückgang in den letzten Jahrzehnten immer noch stark vergletschert. Früher dachte man, dass die alpine Zone über rund 2000 m über Meer vom prähistorischen Menschen nicht begangen wurde. 1991 zeigte die Entdeckung des Mannes aus dem Eis beim Tisenjoch („Ötzi“) in über 3000 m Höhe, dass schon im Neolithikum Vorstösse bis in die vergletscherten Gebiete der Alpen stattgefunden haben. Die ältesten Spuren am Schnidejoch, einem Pass in den Berner Alpen der Schweiz, reichen bis in die Zeit zwischen 4800 und 4500 v.Chr. zurück. Der Pass wurde auch in der Frühen Bronzezeit benutzt, wie zahlreiche Objekte aus der Zeit zwischen etwa 2200 und 1600 v.Chr. belegen. Frühbronzezeitliche Funde liegen auch vom Lötschenpass, einem zweiten Passübergang in den Berner Alpen vor. Nördlich dieser Übergänge befinden sich die bekannten frühbronzezeitlichen Gräber des Berner Oberlandes (Region des unteren Thunersees), südlich davon stammen zahlreiche frühbronzezeitliche Funde aus Gräbern und Nekropolen im Rhonetal. Dank der Erhaltung von organischem Material bieten die Eisfundstellen wertvolle Einblicke zur Frequentierung der Hochalpen. Neben Bohlenwegen, Strassen und Brücken bilden Pässe wichtige Elemente des prähistorischen terrestrischen Transportsystems.

Dynamics of the last four glacial terminations recorded in Lake Van, Turkey

A well-dated suite of Lake Van climate-proxy data covering the last 360 ka documents environmental changes over 4 glacial/interglacial cycles in Eastern Anatolia, Turkey. The picture of cold and dry glacials and warm and wet interglacials emerging from pollen, organic carbon, authigenic carbonate content, elemental profiling by XRF and lithological analyses is inconsistent with classical interpretation of ox- ygen isotopic composition of carbonates pointing to a more complex pattern in Lake Van region. Detailed analysis of glacial terminations allows for the constraining of a depositional model explaining different patterns observed in all the proxies. We hypothesize that variations in relative contribution of rainfall, snowmelt and glacier meltwater recharging the basin have a very important role for all sedimentary processes in Lake Van. Lake level of glacial Lake Van, predominantly fed by snowmelt, was low, the water column was oxic, and carbonates precipitating in the epilimnion recorded the light isotopic signature of inflow. During terminations, increasing rainfall and significant supply of mountain glaciers' meltwater contributed to lake level rise. Increased rainfall enhanced density gradients in the water column, and hindered mixing leading to development of bottom-water anoxia. Carbonates precipitating during terminations show large fluctuations in their isotopic composition. Full interglacial conditions in Lake Van are characterized by high or slowly falling lake level. Rainfall and snowmelt feed the lake but due to re-established mixing, the isotopic composition of authigenic carbonates is heavier and closer to that of evaporation-influenced lake water than that of runoff representing snowmelt and atmospheric precipitation.

Validation of climate model-inferred regional temperature change for late-glacial Europe

Comparisons of climate model hindcasts with independent proxy data are essential for assessing model performance in non-analogue situations. However, standardized palaeoclimate data sets for assessing the spatial pattern of past climatic change across continents are lacking for some of the most dynamic episodes of Earth’s recent past. Here we present a new chironomid-based palaeotemperature dataset designed to assess climate model hindcasts of regional summer temperature change in Europe during the late-glacial and early Holocene. Latitudinal and longitudinal patterns of inferred temperature change are in excellent agreement with simulations by the ECHAM-4 model, implying that atmospheric general circulation models like ECHAM-4 can successfully predict regionally diverging temperature trends in Europe, even when conditions differ significantly from present. However, ECHAM-4 infers larger amplitudes of change and higher temperatures during warm phases than our palaeotemperature estimates, suggesting that this and similar models may overestimate past and potentially also future summer temperature changes in Europe.

Terrestrial and submarine evidence for the extent and timing of the Last Glacial Maximum and the onset of deglaciation on the maritime-Antarctic and sub-Antarctic islands

This paper is the maritime and sub–Antarctic contribution to the Scientific Committee for Antarctic Research (SCAR) Past Antarctic Ice Sheet Dynamics (PAIS) community Antarctic Ice Sheet reconstruction. The overarching aim for all sectors of Antarctica was to reconstruct the Last Glacial Maximum (LGM) ice sheet extent and thickness, and map the subsequent deglaciation in a series of 5000 year time slices. However, our review of the literature found surprisingly few high quality chronological constraints on changing glacier extents on these timescales in the maritime and sub–Antarctic sector. Therefore, in this paper we focus on an assessment of the terrestrial and offshore evidence for the LGM ice extent, establishing minimum ages for the onset of deglaciation, and separating evidence of deglaciation from LGM limits from those associated with later Holocene glacier fluctuations. Evidence included geomorphological descriptions of glacial landscapes, radiocarbon dated basal peat and lake sediment deposits, cosmogenic isotope ages of glacial features and molecular biological data. We propose a classification of the glacial history of the maritime and sub–Antarctic islands based on this assembled evidence. These include: (Type I) islands which accumulated little or no LGM ice; (Type II) islands with a limited LGM ice extent but evidence of extensive earlier continental shelf glaciations; (Type III) seamounts and volcanoes unlikely to have accumulated significant LGM ice cover; (Type IV) islands on shallow shelves with both terrestrial and submarine evidence of LGM (and/or earlier) ice expansion; (Type V) Islands north of the Antarctic Polar Front with terrestrial evidence of LGM ice expansion; and (Type VI) islands with no data. Finally, we review the climatological and geomorphological settings that separate the glaciological history of the islands within this classification scheme.

Peak Last Glacial weathering intensity on the North American continent recorded by the authigenic Hf isotope composition of North Atlantic deep-sea sediments

We have retrieved radiogenic hafnium (Hf) isotope compositions (ɛHf) from authigenic Fe–Mn oxyhydroxides of deep northwest Atlantic sediments deposited over the past 26 ka to investigate the oceanic evidence of changes in dissolved weathering inputs from NE America during the last deglaciation. The extraction of seawater-derived Hf isotopic compositions from Fe–Mn oxyhydroxides is not a standard procedure. Comparisons between the Al/Hf ratios and Hf isotopic compositions of the chemically extracted authigenic phase on the one hand, and those of the corresponding detrital fractions on the other, provide evidence that the composition of past seawater has been reliably obtained for most sampled depths with our leaching procedures. This is endorsed most strongly by data for a sediment core from 4250 m water depth at the deeper Blake Ridge, for which consistent replicates were produced throughout. The Hf isotopic composition of the most recent sample in this core also closely matches that of nearby present day central North Atlantic seawater. Comparison with previously published seawater Nd and Pb isotope compositions obtained on the same cores shows that both Hf and Pb were released incongruently during incipient chemical weathering, but responded differently to the deglacial retreat of the Laurentide Ice Sheet. Hafnium was released more congruently during peak glacial conditions of the Last Glacial Maximum (LGM) and changed to typical incongruent interglacial ɛHf signatures either during or shortly after the LGM. This indicates that some zircon-derived Hf was released to seawater during the LGM. Conversely, there is no clear evidence for an increase in the influence of weathering of Lu-rich mineral phases during deglaciation, possibly since relatively unradiogenic Hf contributions from feldspar weathering were superimposed. While the authigenic Pb isotope signal in the same marine sediment samples traced peak chemical weathering rates on continental North America during the transition to the Holocene a similar incongruent excursion is notably absent in the Hf isotope record. The early change towards more radiogenic ɛHf in relation to the LGM may provide direct evidence for the transition from a cold-based to a warm-based Laurentide Ice Sheet on the Atlantic sector of North America.

Tracer applications of noble gas radionuclides in the geosciences

Noble gas radionuclides, including 81Kr (t1/2 = 229,000 years), 85Kr (t1/2 = 10.8 years), and 39Ar (t1/2 = 269 years), possess nearly ideal chemical and physical properties for studies of earth and environmental processes. Recent advances in Atom Trap Trace Analysis (ATTA), a laser-based atom counting method, have enabled routine measurements of the radiokrypton isotopes, as well as the demonstration of the ability to measure 39Ar in environmental samples. Here we provide an overview of the ATTA technique, and a survey of recent progress made in several laboratories worldwide.We review the application of noble gas radionuclides in the geosciences and discuss how ATTA can help advance these fields, specifically: determination of groundwater residence times using 81Kr, 85Kr, and 39Ar; dating old glacial ice using 81Kr; and an 39Ar survey of the main water masses of the oceans, to study circulation pathways and estimate mean residence times. Other scientific questions involving a deeper circulation of fluids in the Earth's crust and mantle are also within the scope of future applications. We conclude that the geoscience community would greatly benefit from an ATTA facility dedicated to this field, with instrumentation for routine measurements, as well as for research on further development of ATTA methods.