Younger Dryas and Allerød summer temperatures at Gerzensee (Switzerland) inferred from fossil pollen and cladoceran assemblages

Linear- and unimodal-based inference models for mean summer temperatures (partial least squares, weighted averaging, and weighted averaging partial least squares models) were applied to a high-resolution pollen and cladoceran stratigraphy from Gerzensee, Switzerland. The time-window of investigation included the Allerød, the Younger Dryas, and the Preboreal. Characteristic major and minor oscillations in the oxygen-isotope stratigraphy, such as the Gerzensee oscillation, the onset and end of the Younger Dryas stadial, and the Preboreal oscillation, were identified by isotope analysis of bulk-sediment carbonates of the same core and were used as independent indicators for hemispheric or global scale climatic change. In general, the pollen-inferred mean summer temperature reconstruction using all three inference models follows the oxygen-isotope curve more closely than the cladoceran curve. The cladoceran-inferred reconstruction suggests generally warmer summers than the pollen-based reconstructions, which may be an effect of terrestrial vegetation not being in equilibrium with climate due to migrational lags during the Late Glacial and early Holocene. Allerød summer temperatures range between 11 and 12°C based on pollen, whereas the cladoceran-inferred temperatures lie between 11 and 13°C. Pollen and cladocera-inferred reconstructions both suggest a drop to 9–10°C at the beginning of the Younger Dryas. Although the Allerød–Younger Dryas transition lasted 150–160 years in the oxygen-isotope stratigraphy, the pollen-inferred cooling took 180–190 years and the cladoceran-inferred cooling lasted 250–260 years. The pollen-inferred summer temperature rise to 11.5–12°C at the transition from the Younger Dryas to the Preboreal preceded the oxygen-isotope signal by several decades, whereas the cladoceran-inferred warming lagged. Major discrepancies between the pollen- and cladoceran-inference models are observed for the Preboreal, where the cladoceran-inference model suggests mean summer temperatures of up to 14–15°C. Both pollen- and cladoceran-inferred reconstructions suggest a cooling that may be related to the Gerzensee oscillation, but there is no evidence for a cooling synchronous with the Preboreal oscillation as recorded in the oxygen-isotope record. For the Gerzensee oscillation the inferred cooling was ca. 1 and 0.5°C based on pollen and cladocera, respectively, which lies well within the inherent prediction errors of the inference models.

Multiple causes of the Younger Dryas cold period

The Younger Dryas cooling event disrupted the overall warming trend in the North Atlantic region during the last deglaciation. Climate change during the Younger Dryas was abrupt and thus provides insights into the sensitivity of the climate system to perturbations. The sudden Younger Dryas cooling has traditionally been attributed to a shutdown of the Atlantic Meridional Overturning Circulation by meltwater discharges. However, alternative explanations such as strong negative radiative forcing14 and a shift in atmospheric circulation have also been offered. Here we investigate the importance of these different forcings in coupled climate model experiments constrained by data assimilation. We find that the Younger Dryas climate signal as registered in proxy evidence is best simulated using a combination of processes: a weakened Atlantic Meridional Overturning Circulation, moderate negative radiative forcing and an altered atmospheric circulation. We conclude that none of the individual mechanisms alone provide a plausible explanation for the Younger Dryas cold period. We suggest that the triggers for abrupt climate changes such as the Younger Dryas are more complex than suggested so far, and that studies on the response of the climate system to perturbations should account for this complexity.

Quantification of biotic responses to rapid climatic changes around the Younger Dryas — a synthesis

To assess the presence or absence of lags in biotic responses to rapid climatic changes, we: (1) assume that the δ18O in biogenically precipitated carbonates record global or hemispheric climatic change at the beginning and at the end of the Younger Dryas without any lag at our two study sites of Gerzensee and Leysin, Switzerland; (2) derive a time scale by correlating the δ18O record from these two sites with the δ18O record of the GRIP ice core; (3) measure δ18O records in ostracods and molluscs to check the record in the bulk samples and to detect possible hydrological changes; (4) analyse at Gerzensee and Leysin as well as at two additional sites (that lack carbonates and hence a δ18O record) pollen, plant macrofossils, chironomids, beetles and other insects, and Cladocera; (5) estimate our sampling resolution using the GRIP time scale for the isotope stratigraphies and the biostratigraphies; and (6) summarise the major patterns of compositional change in the biostratigraphies by principal component analysis or correspondence analysis. We conclude that, at the major climatic shifts at the beginning and end of the Younger Dryas, hardly any biotic lags occur (within the sampling resolution of 8–30 years) and that upland vegetation responded as fast as aquatic invertebrates. We suggest that the minor climatic changes associated with the Gerzensee and Preboreal oscillations were weakly recorded in the biostratigraphies at the lowland site, but were more distinct at higher altitudes. Individualistic responses of plant and animal species to climatic change may reflect processes in individuals (e.g. productivity and phenology), in populations (e.g. population dynamics), in spatial distributions (e.g. migrations), and in ecosystems (e.g. trophic state). We suggest that biotic responses may be telescoped together into relatively short periods (50 to 150 years), perhaps disrupting functional interactions among species and thus destabilising ecosystems.

Oxygen isotopes of lake marl at Gerzensee and Leysin (Switzerland), covering the Younger Dryas and two minor oscillations, and their correlation to the GRIP ice core

The ratio of oxygen isotopes is a temperature proxy both in precipitation and in the calcite of lacustrine sediments. The very similar oxygen-isotope records from Greenland ice cores and European lake sediments during the Last Glacial Termination suggest that the drastic climatic changes occurred quasi-simultaneously on an extra-regional, probably hemispheric scale. In order to study temporal relations of the different parameters recorded in lake sediments, for example biotic response times to rapid climatic changes, a precise chronology is required. In unlaminated lake sediments there is not yet available a method to provide a high-resolution chronology, especially for periods with radiocarbon plateaux. Alternatively, an indirect time scale can be constructed by linking the lake stratigraphy with other well-dated climate records. New oxygen-isotope records from Gerzensee and Leysin, with an estimated sampling resolution of between 15 and 40 years, match the Greenlandic isotope record in many details. Under the assumption that the main variations in temperature and thus in oxygen isotopes occurred about simultaneously in Greenland and Switzerland, we have assigned a time scale to the lake sediments of Gerzensee and Leysin by wiggle-matching their stable-isotope records with those of Greenland ice cores, which are among the best dated climatic archives. We estimate a precision of 20 to 100 years during the Last Glacial Termination.

Isotope signature of the Younger Dryas and two minor oscillations at Gerzensee (Switzerland): palaeoclimatic and palaeolimnologic interpretation based on bulk and biogenic carbonates

Oxygen- and carbon-isotope ratios in the carbonate of benthic ostracodes (Pseudocandona marchica) and molluscs (Pisidium ssp.) were measured across the transitions bordering the Younger Dryas chronozone in littoral lacustrine cores from Gerzensee (Switzerland). The specific biogenic carbonate records confirm the major shifts already visible in the continuous bulk-carbonate oxygen-isotope record (δ18OCc). If corrected for their vital offsets, oxygen-isotope ratios of Pisidium and juvenile P. marchica, both formed in summer, are almost identical to δ18OCc. This bulk carbonate is mainly composed of encrustations of benthic macrophythes (Chara ssp.), also mainly produced during summer. Adult P. marchica, which calcify in winter, show consistently higher δ18O, larger shifts across both transitions, and short positive excursions compared with the summer forms, especially during early Preboreal. Despite such complexity, the δ18O of adult P. marchica probably reflects more accurately the variations of the δ18O of former lake water because, during winter, calcification temperatures are less variable and the water column isotopically uniform. The difference between normalised δ18O of calcite precipitated in winter to that formed in summer can be used to estimate the minimum difference between summer and winter water temperatures. In general, the results indicate warmer summers during the late Allerød and early Preboreal compared with the Younger Dryas. Altogether, the isotopic composition of lake water (δ18OL) and of the dissolved inorganic carbonate (δ13CDIC) reconstructed from adult Pseudocandona marchica, as well as the seasonal water temperature contrasts, indicate that the major shifts in the δ18O of local precipitation at Gerzensee were augmented by changes of the lake's water balance, with relatively higher evaporative loss occurring during the Allerød compared with the Younger Dryas. It is possible that during the early Preboreal the lake might even have been hydrologically closed for a short period. We speculate that such hydrologic changes reflect a combination of varying evapotranspiration and a rearrangement of groundwater recharge during those climatic shifts.

NGRIP CH4 concentration from 120 to 10 kyr before present and its relation to a δ15N temperature reconstruction from the same ice core

During the last glacial cycle, Greenland temperature showed many rapid temperature variations, the so-called Dansgaard–Oeschger (DO) events. The past atmospheric methane concentration closely followed these temperature variations, which implies that the warmings recorded in Greenland were probably hemispheric in extent. Here we substantially extend and complete the North Greenland Ice Core Project (NGRIP) methane record from the Preboreal Holocene (PB) back to the end of the last interglacial period with a mean time resolution of 54 yr. We relate the amplitudes of the methane increases associated with DO events to the amplitudes of the local Greenland NGRIP temperature increases derived from stable nitrogen isotope (δ15N) measurements, which have been performed along the same ice core (Kindler et al., 2014). We find the ratio to oscillate between 5 parts per billion (ppb) per °C and 18 ppb °C−1 with the approximate frequency of the precessional cycle. A remarkably high ratio of 25.5 ppb °C−1 is reached during the transition from the Younger Dryas (YD) to the PB. Analysis of the timing of the fast methane and temperature increases reveals significant lags of the methane increases relative to NGRIP temperature for DO events 5, 9, 10, 11, 13, 15, 19, and 20. These events generally have small methane increase rates and we hypothesize that the lag is caused by pronounced northward displacement of the source regions from stadial to interstadial. We further show that the relative interpolar concentration difference (rIPD) of methane is about 4.5% for the stadials between DO events 18 and 20, which is in the same order as in the stadials before and after DO event 2 around the Last Glacial Maximum. The rIPD of methane remains relatively stable throughout the full last glacial, with a tendency for elevated values during interstadial compared to stadial periods.

Chronology of Lateglacial ice flow reorganization and deglaciation in the Gotthard Pass area, Central Swiss Alps, based on cosmogenic 10Be and in situ 14C

We reconstruct the timing of ice flow reconfiguration and deglaciation of the Central Alpine Gotthard Pass, Switzerland, using cosmogenic 10Be and in situ14C surface exposure dating. Combined with mapping of glacial erosional markers, exposure ages of bedrock surfaces reveal progressive glacier downwasting from the maximum LGM ice volume and a gradual reorganization of the paleoflow pattern with a southward migration of the ice divide. Exposure ages of ∼16–14 ka (snow corrected) give evidence for continuous early Lateglacial ice cover and indicate that the first deglaciation was contemporaneous with the decay of the large Gschnitz glacier system. In agreement with published ages from other Alpine passes, these data support the concept of large transection glaciers that persisted in the high Alps after the breakdown of the LGM ice masses in the foreland and possibly decayed as late as the onset of the Bølling warming. A younger group of ages around ∼12–13 ka records the timing of deglaciation following local glacier readvance during the Egesen stadial. Glacial erosional features and the distribution of exposure ages consistently imply that Egesen glaciers were of comparatively small volume and were following a topographically controlled paleoflow pattern. Dating of a boulder close to the pass elevation gives a minimum age of 11.1 ± 0.4 ka for final deglaciation by the end of the Younger Dryas. In situ14C data are overall in good agreement with the 10Be ages and confirm continuous exposure throughout the Holocene. However, in situ14C demonstrates that partial surface shielding, e.g. by snow, has to be incorporated in the exposure age calculations and the model of deglaciation.

Lateglacial and early Holocene glaciation in the tropical Andes caused by La Niña-like conditions

The response of the tropics to North Atlantic cold events, such as Heinrich Event I (H-I, ∼ 17–15 ka) and the Younger Dryas (YD, 12.7–11.5 ka), is still one of the most tantalizing, yet unresolved issues in paleoclimatology. The advent of surface exposure dating has therefore instigated the establishment of glacial chronologies in the tropical Andes to investigate potential climate teleconnections. Here, we present new exposure ages from the Cordillera Cochabamba (17°17′S), Bolivia, that reveal glacial advances during H-I and YD, as well as during the Early Holocene. Our chronology correlates well with cold sea surface temperatures in the eastern tropical Pacific, which indicates that La Niña-like conditions, i.e. forcings intrinsic to the tropics, played a key role for moisture advection and glaciation in the tropical Andes.

Fingerprints of changes in the terrestrial carbon cycle in response to large reorganizations in ocean circulation

CO2 and carbon cycle changes in the land, ocean and atmosphere are investigated using the comprehensive carbon cycle-climate model NCAR CSM1.4-carbon. Ensemble simulations are forced with freshwater perturbations applied at the North Atlantic and Southern Ocean deep water formation sites under pre-industrial climate conditions. As a result, the Atlantic Meridional Overturning Circulation reduces in each experiment to varying degrees. The physical climate fields show changes qualitatively in agreement with results documented in the literature, but there is a clear distinction between northern and southern perturbations. Changes in the physical variables, in turn, affect the land and ocean biogeochemical cycles and cause a reduction, or an increase, in the atmospheric CO2 concentration by up to 20 ppmv, depending on the location of the perturbation. In the case of a North Atlantic perturbation, the land biosphere reacts with a strong reduction in carbon stocks in some tropical locations and in high northern latitudes. In contrast, land carbon stocks tend to increase in response to a southern perturbation. The ocean is generally a sink of carbon although large reorganizations occur throughout various basins. The response of the land biosphere is strongest in the tropical regions due to a shift of the Intertropical Convergence Zone. The carbon fingerprints of this shift, either to the south or to the north depending on where the freshwater is applied, can be found most clearly in South America. For this reason, a compilation of various paleoclimate proxy records of Younger Dryas precipitation changes are compared with our model results. The proxy records, in general, show good agreement with the model's response to a North Atlantic freshwater perturbation.