Cold-season temperatures in the European Alps during the past millennium: variability, seasonality and recent trends

This study presents a proxy-based, quantitative reconstruction of cold-season (mean October to May, TOct–May) air temperatures covering nearly the entire last millennium (AD 1060–2003, some hiatuses). The reconstruction was based on subfossil chrysophyte stomatocyst remains in the varved sediments of high-Alpine Lake Silvaplana, eastern Swiss Alps (46°27’N, 9°48′W, 1791 m a.s.l.). Previous studies have demonstrated the reliability of this proxy by comparison to meteorological data. Cold-season air temperatures could therefore be reconstructed quantitatively, at a high resolution (5-yr) and with high chronological accuracy. Spatial correlation analysis suggests that the reconstruction reflects cold season climate variability over the high- Alpine region and substantial parts of central and western Europe. Cold-season temperatures were characterized by a relatively stable first part of the millennium until AD 1440 (2σ of 5-yr mean values = 0.7 °C) and highly variable TOct–May after that (AD 1440–1900, 2σ of 5-yr mean values = 1.3 °C). Recent decades (AD, 1991-present) were unusually warm in the context of the last millennium (exceeding the 2σ-range of the mean decadal TOct–May) but this warmth was not unprecedented. The coolest decades occurred from AD 1510–1520 and AD 1880–1890. The timing of extremely warm and cold decades is generally in good agreement with documentary data representing Switzerland and central European lowlands. The transition from relatively stable to highly variable TOct–May coincided with large changes in atmospheric circulation patterns in the North Atlantic region. Comparison of reconstructed cold season temperatures to the North Atlantic Oscillation index (NAO) during the past 1000 years showed that the relatively stable and warm conditions at the study site until AD 1440 coincided with a persistent positive mode of the NAO. We propose that the transition to large TOct–May variability around AD 1440 was linked to the subsequent absence of this persistent zonal flow pattern, which would allow other climatic drivers to gain importance in the study area. From AD 1440–1900, the similarity of reconstructed TOct–May to reconstructed air pressure in the Siberian High suggests a relatively strong influence of continental anticyclonic systems on Alpine cold season climate parameters during periods when westerly airflow was subdued. A more continental type of atmospheric circulation thus seems to be characteristic for the Little Ice Age in Europe. Comparison of Toct–May to summer temperature reconstructions from the same study site shows that, as expected, summer and cold season temperature trends and variability differed completely throughout nearly the entire last 1000 years. Since AD 1980, however, summer and cold season temperatures show a simultaneous, strong increase, which is unprecedented in the context of the last millennium. We suggest that the most likely explanation for this recent trend is anthropogenic greenhouse gas (GHG) forcing.

A chrysophyte stomatocyst-based reconstruction of cold-season air temperature from Alpine Lake Silvaplana (AD 1500–2003); methods and concepts for quantitative inferences

Relatively little is known about past cold-season temperature variability in high-Alpine regions because of a lack of natural cold-season temperature proxies as well as under-representation of high-altitude sites in meteorological, early-instrumental and documentary data sources. Recent studies have shown that chrysophyte stomatocysts, or simply cysts (sub-fossil algal remains of Chrysophyceae and Synurophyceae), are among the very few natural proxies that can be used to reconstruct cold-season temperatures. This study presents a quantitative, high-resolution (5-year), cold-season (Oct–May) temperature reconstruction based on sub-fossil chrysophyte stomatocysts in the annually laminated (varved) sediments of high-Alpine Lake Silvaplana, SE Switzerland (1,789 m a.s.l.), since AD 1500. We first explore the method used to translate an ecologically meaningful variable based on a biological proxy into a simple climate variable. A transfer function was applied to reconstruct the ‘date of spring mixing’ from cyst assemblages. Next, statistical regression models were tested to convert the reconstructed ‘dates of spring mixing’ into cold-season surface air temperatures with associated errors. The strengths and weaknesses of this approach are thoroughly tested. One much-debated, basic assumption for reconstructions (‘stationarity’), which states that only the environmental variable of interest has influenced cyst assemblages and the influence of confounding variables is negligible over time, is addressed in detail. Our inferences show that past cold-season air-temperature fluctuations were substantial and larger than those of other temperature reconstructions for Europe and the Alpine region. Interestingly, in this study, recent cold-season temperatures only just exceed those of previous, multi-decadal warm phases since AD 1500. These findings highlight the importance of local studies to assess natural climate variability at high altitudes.

Identification of glacial meltwater runoff in a karstic environment and its implication for present and future water availability

Glaciers all over the world are expected to continue to retreat due to the global warming throughout the 21st century. Consequently, future seasonal water availability might become scarce once glacier areas have declined below a certain threshold affecting future water management strategies. Particular attention should be paid to glaciers located in a karstic environment, as parts of the meltwater can be drained by underlying karst systems, making it difficult to assess water availability. In this study tracer experiments, karst modeling and glacier melt modeling are combined in order to identify flow paths in a high alpine, glacierized, karstic environment (Glacier de la Plaine Morte, Switzerland) and to investigate current and predict future downstream water availability. Flow paths through the karst underground were determined with natural and fluorescent tracers. Subsequently, geologic information and the findings from tracer experiments were assembled in a karst model. Finally, glacier melt projections driven with a climate scenario were performed to discuss future water availability in the area surrounding the glacier. The results suggest that during late summer glacier meltwater is rapidly drained through well-developed channels at the glacier bottom to the north of the glacier, while during low flow season meltwater enters into the karst and is drained to the south. Climate change projections with the glacier melt model reveal that by the end of the century glacier melt will be significantly reduced in the summer, jeopardizing water availability in glacier-fed karst springs.

Quantitative high-resolution warm season rainfall recorded in varved sediments of Lake Oeschinen, northern Swiss Alps: calibration and validation AD 1901–2008

High-resolution, well-calibrated records of lake sediments are critically important for quantitative climate reconstructions, but they remain a methodological and analytical challenge. While several comprehensive paleotemperature reconstructions have been developed across Europe, only a few quantitative high-resolution studies exist for precipitation. Here we present a calibration and verification study of lithoclastic sediment proxies from proglacial Lake Oeschinen (46°30′N, 7°44′E, 1,580 m a.s.l., north–west Swiss Alps) that are sensitive to rainfall for the period AD 1901–2008. We collected two sediment cores, one in 2007 and another in 2011. The sediments are characterized by two facies: (A) mm-laminated clastic varves and (B) turbidites. The annual character of the laminae couplets was confirmed by radiometric dating (210Pb, 137Cs) and independent flood-layer chronomarkers. Individual varves consist of a dark sand-size spring-summer layer enriched in siliciclastic minerals and a lighter clay-size calcite-rich winter layer. Three subtypes of varves are distinguished: Type I with a 1–1.5 mm fining upward sequence; Type II with a distinct fine-sand base up to 3 mm thick; and Type III containing multiple internal microlaminae caused by individual summer rainstorm deposits. Delta-fan surface samples and sediment trap data fingerprint different sediment source areas and transport processes from the watershed and confirm the instant response of sediment flux to rainfall and erosion. Based on a highly accurate, precise and reproducible chronology, we demonstrate that sediment accumulation (varve thickness) is a quantitative predictor for cumulative boreal alpine spring (May–June) and spring/summer (May–August) rainfall (rMJ = 0.71, rMJJA = 0.60, p < 0.01). Bootstrap-based verification of the calibration model reveals a root mean squared error of prediction (RMSEPMJ = 32.7 mm, RMSEPMJJA = 57.8 mm) which is on the order of 10–13 % of mean MJ and MJJA cumulative precipitation, respectively. These results highlight the potential of the Lake Oeschinen sediments for high-resolution reconstructions of past rainfall conditions in the northern Swiss Alps, central and eastern France and south-west Germany.

Atmospheric CO₂ response to volcanic eruptions: The role of ENSO, season, and variability

Tropical explosive volcanism is one of the most important natural factors that significantly impact the climate system and the carbon cycle on annual to multi-decadal time scales. The three largest explosive eruptions in the last 50�years�Agung, El Chichón, and Pinatubo�occurred in spring/summer in conjunction with El Niño events and left distinct negative signals in the observational temperature and CO2 records. However, confounding factors such as seasonal variability and El Niño-Southern Oscillation (ENSO) may obscure the forcing-response relationship. We determine for the first time the extent to which initial conditions, i.e., season and phase of the ENSO, and internal variability influence the coupled climate and carbon cycle response to volcanic forcing and how this affects estimates of the terrestrial and oceanic carbon sinks. Ensemble simulations with the Earth System Model (Climate System Model 1.4-carbon) predict that the atmospheric CO2 response is �60 larger when a volcanic eruption occurs during El Niño and in winter than during La Niña conditions. Our simulations suggest that the Pinatubo eruption contributed 11�±�6 to the 25�Pg terrestrial carbon sink inferred over the decade 1990�1999 and �2�±�1 to the 22�Pg oceanic carbon sink. In contrast to recent claims, trends in the airborne fraction of anthropogenic carbon cannot be detected when accounting for the decadal-scale influence of explosive volcanism and related uncertainties. Our results highlight the importance of considering the role of natural variability in the carbon cycle for interpretation of observations and for data-model intercomparison.

Variations of clinopyroxene/melt element partitioning during assimilation of olivine/peridotite by low-Mg diorite magma

The effects of crystal chemistry and melt composition on the control of clinopyroxene/melt element partitioning (D) during the assimilation of olivine/peridotite by felsic magma have been investigated in Mesozoic high-Mg diorites from North China. The assimilation resulted in significant increase of Mg, Cr and Ni and only slight (< 30%) decrease of incompatible elements of the magma, and the compositional variations have been mirrored by the normally and reversely zoned clinopyroxene microphenocrysts formed at the early stage of the magma evolution. The Mg# [100 × Mg / (Mg + Fe)] values of the reversely zoned clinopyroxenes increase from 65 to 75 in the core to 85–90 in the high-Mg midsection, and reduce back to 73–79 at the rim. Trace element profiles across all these clinopyroxene domains have been measured by LA-ICP-MS. The melt trace element composition has been constrained from bulk rock analyses of the fine-grained low- and high-Mg diorites. Clinopyroxene/melt partition coefficients for rare earth elements (REE) and Y in the high-Mg group zonings (Mg# > 73–79, DDy < 1.2) are positively correlated with tetrahedral IVAl and increase by a factor of 3–4 as tetrahedral IVAl increases from 0.01 to 0.1 per formula unit (pfu). These systematic variations are interpreted to be controlled by the clinopyroxene composition. In contrast, partition coefficients for low-Mg group zonings (Mg# < 75–79, DDy > 1.2) are elevated by up to an order of magnitude (for REE and Y) or more (for Zr and Hf) at similar IVAl, indicating dominant control of melt composition/structure. DZr and DHf show a larger sensitivity to the compositional change of crystal and melt than DREE. DTi values for the low- and high-Mg zonings show a uniform dependence on IVAl. DSr and DLi are insensitive to the compositional change of clinopyroxene and melt, resulting in Sr depletions in the clinopyroxene zonings with elevated REE without crystallization of plagioclase. Our observations show that crystal chemistry and melt composition/structure may alternatively control clinopyroxene/melt partitioning during the assimilation of peridotite by felsic magma, and may be useful for deciphering clinopyroxene compositions and related crust–mantle processes.