Abstract
The European Alps are an effective barrier for meridional moisture transport and are thus uniquely placed to record shifts in the North Atlantic storm track pattern associated with the waxing and waning of Late-Pleistocene Northern Hemisphere ice sheets. The lack of well-dated terrestrial proxy records spanning this time period, however, renders the reconstruction of past atmospheric patterns difficult. Here we present a precisely dated, continuous terrestrial record of meteoric precipitation in Europe between 30 and 14.7âka. In contrast to present-day conditions, our speleothem data provide strong evidence for preferential advection of moisture from the South across the Alps supporting a southward shift of the storm track during the local Last Glacial Maximum (that is, 26.5â23.5âka). Moreover, our age control indicates that this circulation pattern preceded the Northern Hemisphere precession maximum by ~3âka, suggesting that obliquity may have played a considerable role in the Alpine ice aggradation.
Similar content being viewed by others
Introduction
The Last Glacial Maximum (LGM) characterizes the climax of a pronounced sea-level minimum associated with the global ice extent during the Upper Pleistocene, between ca 26.5 and 19âka (refs 1, 2). While atmospheric circulation models support a southward shift of the North Atlantic storm tracks to ~40â°N (refs 3, 4) paleoclimate proxy records point to a global mean temperature decrease of between 4 and 7â°C when compared with preindustrial times5. In the European Alps, the extent of the regional ice sheet has been comprehensively documented by detailed geomorphological mapping6,7. Observations of trimlines reveal the presence of three major ice domes, all located south of the present climate divide, suggesting that precipitation was mostly associated with air advection from the South8. However, in the absence of continuous terrestrial palaeoclimate records from the northern Alpine mountain chain, it is yet unclear whether this precipitation pattern was associated with increased cyclogenesis in the western Mediterranean9 or if it relates to an overall shift in the storm trajectories from the eastern North Atlantic8. This uncertainty in regional atmospheric circulation patterns has major implications, primarily for the regional ice build-up and thus the timing of valley glacier advances but also for the distribution of temperature and precipitation responses and thus implicitly also for vegetation refugia.
Reconstructing LGM atmospheric circulation in the European Alps remains inherently difficult and requires accurately dated high-resolution climate proxy records. High-altitude ice cores certainly represent the most direct precipitation archive at mid-latitudes, but barely cover the Holocene10. Ice and permafrost conditions in the Alpine foreland further hinder formation/preservation of possible lake records and biogenic proxies during the LGM. Speleothems, secondary carbonate deposits found in caves, represent an alternative archive for Late Pleistocene climate history provided that the temperature in the subsurface stayed above freezing and local waterârock interactions gave rise to calcite-supersaturated water.
Here we present the first continuous continental European record of meteoric precipitation obtained from speleothems grown during Marine Isotope Stage 2. We show that our record reflects atmospheric circulation patterns in the Northern Hemisphere during the LGM, supporting a southward shift of the North Atlantic storm track about 3âka before the precession maximum. We conclude that the preferential advection of moisture from the South likely favoured regional differences in the ice build-up in the Alps, which suggests contrasting responses of valley glacier advances during the LGM.
Results
Study site
Two coeval speleothems (7H-2 and 7H-3) were sampled in the Sieben Hengste (7H) cave system (Bernese Alps, Switzerland; Fig. 1), a network comprising more than 150âkm of explored cave passages. The galleries formed within Urgonian carbonates (Schrattenkalk Formation), which are locally overlain by low-permeable Eocene sandstones (Supplementary Fig. 1). Several lines of evidence, including alignments of erratic boulders and glacial striations, point to an ice-covered hydrological catchment at 1,700â2,000âm a.s.l. during the LGM7. Based on geomorphological evidence from the nearby Napf palaeoglacier, the local Equilibrium Line Altitude (ELA) during the LGM is estimated at ca 1,250âm a.s.l. Despite the low temperatures expected during the LGM, summer rainfall and seasonal melting of winter-snow must have been sufficient to allow for recharge into the Sieben Hengste karst system and ensure speleothem deposition. Concentrated infiltration along major tectonic discontinuities likely maintained a locally constant temperature boundary condition in the epikarst close to the melting point (that is, 0â°C), whereas hydrological drainage ensured geothermal heat advection off the phreatic zone. Accordingly, the vertical temperature distribution in the Sieben Hengste karst system most likely followed the adiabatic lapse rate of humid air11, implying an isotherm of ca 1.0±0.5â°C at 1,500âm a.s.l.
Conceptual view of southerly moisture transport to the Alps during the Last Glacial Maximum (blue arrow) from the subtropical reservoir (blue shading). The advection of moisture is triggered by Rossby-wave breaking of the jet stream over Western Europe (red line), which is enhanced due to the presence of the Fennoscandian ice shield (light blue; 6) and the associated semi-permanent blocking high. Dust-mobilizing storms are induced further East by jet stream excursions propagating to Central and Eastern Asia (orange area and arrow). 7H: Sieben Hengste cave system; NGRIP: ice core; GIK-23415 and MD01-2443: deep sea sediment cores; Rhine and Tagliamento: palaeoglaciers.
Speleothem samples
The two stalagmites (Supplementary Fig. 2) were collected ca 215âm below the ground surface (that is, 1,540âm a.s.l.) in an alcove at the base of a vadose shaft. The candle-shaped stalagmites, 57- and 30-cm long show similar petrographies dominated by bundles of elongated columnar fascicular optic calcite crystals, commonly arranged in fans (Supplementary Fig. 3). This fabric is typical of elevated Mg/Ca ratios in the drip water and agrees with the chemical concentrations measured in the speleothem (Supplementary Table 1). Both, the calcite fabrics and the regular stalagmite morphologies suggest that the karst aquifer properties remained constant over decadal to millennial timescales. The regular recharge was mainly controlled by the glacier hydrology. Hence, the temperature in the cave remained largely constant, and we surmise that the 7H δ18O is primarily controlled by changes in the precipitation isotope composition.
The 7H-record
The 50 uranium-series ages obtained by multicollector inductively coupled plasma mass spectrometry are all in stratigraphic order and reveal continuous calcite deposition from 30.0 to 14.7âka b2k (Supplementary Table 2). Typical age uncertainties of less than 0.5% reflect the elevated 238U concentration of 1,200±300ângâgâ1 associated with low detrital thorium contents. The age model reveals an average growth rate of 40â70âμmâaâ1 (Supplementary Fig. 4). Sub-samples for δ18O and δ13C analyses were milled at 200âμm intervals along the speleothem growth axis providing an average resolution of 4.5 a.
The speleothem oxygen isotope composition (ca 4,200 analyses) ranges between δ18OVPDB â10.5 and â13.5â° and oscillates on decadal to millennial time scales (Fig. 2b). The stalagmite δ18O time series correlate in overlapping segments and, therefore, likely responded to the same hydrological forcing. This hydrological control is emphasized by a high local correlation (r2=0.54) between δ18O and δ13C, which suggests that any decrease in drip rate related to lower recharge forced a coherent C isotope change (Supplementary Fig. 5). Overall, δ13C values stayed remarkably constant at +4.4±0.4â°, indicating that the dissolved inorganic carbon was largely derived from the host rock. A composite δ18O record was obtained by using a Monte Carlo (MC) approach to find the best-possible correlation between adjacent time series (see Methods for details).
(a) Precession (brown) and obliquity (blue)23. (b) 7H-2 (red) and 7H-3 (brown) δ18O time series (21a low-pass filter) plotted with individual U/Th ages and measured errors (2Ï). (c) Alpine glacier advances reconstructed for Tagliamento and Rhine glaciers17,18. (d) North-eastern Atlantic sea-surface temperature15,16. (e) NGRIP dust record33. (f) NGRIP δ18O with a 100a low-pass filter12.
Discussion
The 7H δ18O record shows remarkable similarities to Northern Hemisphere palaeoclimate records, both at orbital and millennial timescales. Positive δ18O excursions of 7H (Fig. 2b) match the North Greenland Ice Core Project (NGRIP) δ18O record (Fig. 2f) for Interstadials (I) 3 and 4 within dating uncertainties12,13 (Table 1). The timing of the rapid δ18O increase at the onset of the Bølling (that is, I-1e at 14.70âka) is identical with the timing in NGRIP and led to a stalagmite growth stop due to the waning of the glacier cover in the hydrological catchment of the cave system. We argue that the high degree of similarity between the 7H-record and Greenland ice cores reflects the common variability of the North Atlantic climate system14. This interpretation is supported by coeval changes in the sea-surface temperature recorded off the Iberian Peninsula15 (core 2443 in Fig. 2d). This covariance is highest when the near-surface poleward heat transfer, between 26.3 and 23.5âka, was strongly reduced16(for example, core 23415 for comparison in Fig. 2d), suggesting internal feedback mechanisms associated with variations in the Atlantic Meridional Overturning Circulation.
In contrast to the NGRIP record, the 7H δ18O record is characterized by an absolute minimum centred at ca 25.3âka, consistent with the maximum Alpine glacier advance inferred from independently dated geomorphological field evidence17,18 (Fig. 2c). In addition, second-order changes in the 7H δ18O time-series correlate within dating uncertainty with minor glacier oscillations reconstructed from proglacial sediments17 and agree with a perialpine vegetation record from Les Echets, eastern France19, which has been interpreted to reflect primary changes in soil moisture availability. Based on these evidences, we surmise that the 7H δ18O record mirrors changes in the mass balance of LGM glaciers in the Alps and thus predominantly reflects changes in the Alpine precipitation pattern8.
The modern regime of moisture transport to Sieben Hengste is dominated by westerly advection (Supplementary Fig. 6), while about 10â30% of the annual precipitation is transported on average from the South across the Alpine mountain range. Intense southerly moisture advection is known to be related to pronounced Rossby-wave breaking west of the Alps, which induces enhanced meridional flow of moist and warm subtropical air towards the mountain range20. Such transport patterns may have been dominant during the LGM when both a larger (winter) sea-ice cover in the North Atlantic21 displaced the baroclinic zone further south and a semi-persistent blocking high over the Fennoscandian ice sheet deflected the North Atlantic jet stream substantially, leading to more frequent Rossby-wave breaking (Fig. 1). Since the tropics and subtropics experienced comparatively less temperature changes during glacial-interglacial cycles, the subtropical moisture reservoir could have provided sufficient moisture for precipitation to build up and maintain ice domes in the Alps with meridional moisture transport events. The strong correlation between the 7H-record and the Alpine ice mass balance (Fig. 2c) suggests that this atmospheric circulation pattern may have been a characteristic mechanism of moisture delivery to the Alps during the LGM.
Accordingly, we propose that, during the LGM, the isotopic composition of 7H primarily recorded changes in the position of the North Atlantic storm track and associated changes of moisture advection to the Alps. This hypothesis is generally consistent with global climate simulations suggesting an overall southward displacement of the storm track during the LGM3,4, accompanied by a significant precipitation increase in the mid-latitudes, in particular around the Iberian Peninsula21. Meanwhile, reconstructed meridional temperature gradients based on Alpine glacierâs ELA depression suggest cold and dry conditions in the northern Alps at the LGM climax22. The presence of liquid water in the 7H karst system at 1,750âm a.s.l., however, challenges an estimated mean annual air temperature of â10 to â15â°C at the ELA of glaciers22 and suggests that major precipitation events must have been associated with substantial advection of warm moist air, predominantly between spring and autumn. Owing to their transient nature, such events are hardly reflected in the mean annual air temperature record23 and the relatively small changes in the temperature of condensation only had a minor influence on the isotopic composition. Similarly, the extraordinarily cold and dry climate during the LGM makes it unlikely to explain the depletion in δ18Occ by a seasonal shift towards more winter precipitation. Instead, we propose that a higher fraction of precipitation was transported across the Alpine mountain range from the South because of a change in the predominant synoptic circulation pattern. The higher orographic barrier along the southerly transport path cause more moisture to condensate resulting in an orographically induced Rayleigh fractionation process (Supplementary Fig. 7). Accordingly, phases of more depleted δ18Occ are associated with the preferential advection of moisture from the South (see Supplementary Discussion), whereas less depleted δ18Occ values primarily reflect air masses reaching the northern Alps directly from the Northwest (Fig. 3). Based on the apparent correlation between the Remanzacco oscillations of the LGM South Alpine Tagliamento paleoglacier17 and the 7H-record, the 50% threshold is assessed at δ18Occ â12.4±0.1â°, which corresponds to ca â16.3â° on the Vienna Standard Mean Ocean Water (VSMOW) scale (Fig. 4). Because δ18O at the oceanic moisture source has not changed substantially during the LGM, a two end-member mixing model allows to derive that the contribution of moisture crossing the Alps may have varied between 25 and 65% during the LGM.
Rayleigh isotope fractionation along the northwestern (blue) and southern moisture transport trajectory (red) for the Last Glacial Maximum. Moisture evaporates from the source region in the Eastern North Atlantic and then follows a fractionation path with orographic rainout up to 1,700âm a.s.l. during northerly advection, or up to 3,000âm a.s.l. for the southerly advection across the Alpine main crest. δ18O values are expressed in VSMOW. Percentages reflect the proportion of precipitation reaching 7H along the northwestern and southern routes, respectively. The ranges reflect uncertainties on the interannual variability (see Supplementary Discussion for details). Based on the 7H δ18O record, we surmise that the proportion of North Atlantic moisture transported from the South to 7H during the LGM was two to three times higher than today.
The correlation with the Alpine Tagliamento glacier advances17 (Fig. 2c) suggests that low δ18O values are primarily associated with moisture transport from the South. The stippled line marks the approximate equilibrium between northwestern and southern moisture transport as identified from the Remanzacco oscillations. A dominant northwestern trajectory reflects a more northerly storm track position, and vice-versa. The secondary y axis shows the inferred isotopic composition of meteoric precipitation assuming calcite precipitation close to thermodynamic equilibrium34 in a cave environment just above freezing point.
Associating the (Alpine) LGM with a change in the synoptic circulation pattern rather than an insolation minimum enables alternative interpretations for the causes of major glacier advances in the Alps. We note that the 7H δ18O minimum predated the precession maximum at 22âka (ref. 24) by ca 3âka, suggesting a combined effect of precession and obliquity on the timing of the last Alpine glaciation, and thus possibly also on its global counterpart25. Evidence of de-phasing with respect to the precessional band has been reported from marine and terrestrial records26, but interpretations were mainly focused on the timing of deglaciation. Although changes in obliquity influence the latitudinal distribution of solar radiation and thus control glacier ablation27, its effect on ice accumulation is less obvious. Yet, the obliquity minimum at 29âka (ref. 24) implies reduced energy input at high latitudes. The associated cooling of the Northern Hemisphere likely favoured a drop in the ELA, as exemplified by the start of speleothem growth at 7H at 30âka. With tropical and subtropical areas being less affected by changes in obliquity, meridional temperature gradients and therefore also pressure gradients became steeper, strengthening the jet stream, which, together with an expansion of the Fennoscandian ice sheet, progressively migrated southwards3. This may have led to enhanced synoptic storm activity in North-Atlantic mid-latitudes, and possibly also in areas further downstream, because of jet stream excursions induced by the interaction between the jet and a quasi-stationary blocking over the Fennoscandian ice sheet (Fig. 1). The concomitant decline of the mid-latitude vegetation cover increased the potential global dust emissions by nearly 38% (ref. 28). An intriguing similarity is noted between the Ca record of NGRIP and the 7H δ18O record (Fig. 2e), which is particularly pronounced between 28 and 23âka (that is, râ¥0.60 between 25.5 and 23.0âka; Supplementary Fig. 8). Elevated dust concentrations in Greenland ice cores correlate with depleted δ18O values in 7H. This increased dust mobilization during periods when the storm track was dominantly south of the Alps supports a potential teleconnection between North Atlantic moisture transport and storminess in East Asian regions (Supplementary Discussion). Concurrently, the preferential transport of moisture towards the southern flank of the Alps favoured snow accumulation in the glaciological catchment of the large piedmont glaciers. These results are not only consistent with independent geomorphological evidence7,8 but also argue for a non-uniform ice build-up in the Alps during the LGM. Accordingly, we expect regional differences in the timing of valley glacier advances depending on the location of the glaciersâ accumulation areas.
We conclude that the Alpine glaciers responded in a sensitive way to Northern Hemisphere atmospheric circulation patterns during the LGM. Therefore, this mountain range represents a key location for identifying large-scale atmospheric reorganizations in the North Atlantic climate systems. The importance of moisture advection towards the Alps for glacier advances calls for high-resolution GCM simulations of the LGM climate, which could be tested extensively using 7H as a benchmark.
Methods
U/Th dating and age modelling
Chemical separation and MC-ICP-MS (Thermo-Finnigan Neptune) measurements of U and Th isotopic ratios were undertaken at the University of Minnesota using procedures similar to those described in ref. 29. The extent of detrital 230Th contamination was estimated and corrected for by measurement of the long-lived chemically equivalent 232Th and assuming a silicate bulk Earth initial 230Th/232Th atomic ratio of 4.4±2.2 à 10â6. Final ages are given as years before the year 2000 AD (a b2k). Age modelling was performed using StalAge7, an algorithm specifically designed for speleothems.
Speleothem geochemistry
Trace elements were measured by Laser ablation inductively coupled plasma mass spectrometry (LA-ICPMS) (RESOlution M-50 prototype at RHUL, 193ânm, ArF excimer) with a two-volume LA cell coupled to an Agilent 7500ce quadrupole ICPMS30 along a ca 8-cm-long profile at the base of the stalagmite (single track). Reported concentrations were calculated following ref. 31 using NIST610 as bracketing external standard and stoichiometric [Ca] of 40% m/m as internal standard.
Sub-samples for δ18O and δ13C analyses were milled at 200âμm intervals along the vertical growth axes of the stalagmites using a Merchantek micromill. Analyses of 0.05â0.35âmg calcite powders were performed at the University of Innsbruck on a ThermoFisher DeltaplusXL isotope ratio mass spectrometer with an analytical precision (1Ï) of 0.08â° for δ18O and 0.06â° for δ13C. Results are reported on the Vienna Pee Dee Belemnite (VPDB) scale and calibrated against NBS19.
The composite 7H-record
The composite 7H δ18O record (Fig. 4) was obtained using a MC approach applied on absolute age determinations to find the best correlation between adjacent time series (ISCAM32). We prescribed a point-wise linear interpolation between data points. The age model was determined from the highest correlation obtained over 100,000 MC simulations using a 50a smoothing. The significance estimation was performed using 2,000 pairs of artificially constructed first-order autoregressive time series (AR1), while each pair was scanned with 1,000 MC simulations to find the best correlation. Finally, the result was re-sampled at annual intervals, detrended and normalized to be compared with other proxy records.
Additional information
How to cite this article: Luetscher, M. et al. North Atlantic storm track changes during the Last Glacial Maximum recorded by Alpine speleothems. Nat. Commun. 6:6344 doi: 10.1038/ncomms7344 (2015).
References
Clark, P. U. et al. The Last Glacial Maximum. Science 325, 710â714 (2009).
Hughes, P. D., Gibbard, P. L. & Ehlers, J. Timing of glaciation during the last glacial cycle: evaluating the concept of a global âLast Glacial Maximumâ (LGM). Earth Sci. Rev. 125, 171â198 (2013).
Laîné, A. et al. Northern hemisphere storm tracks during the last glacial maximum in the PMIP2 ocean-atmosphere coupled models: energetic study, seasonal cycle, precipitation. Clim. Dynam. 32, 593â614 (2009).
Hofer, D., Raible, C. C., Merz, N., Dehnert, A. & Kuhlemann, J. Simulated winter circulation types in the North Atlantic and European region for preindustrial and glacial conditions. Geophys. Res. Lett. 39, L15805 (2012).
Schneider von Deimling, T., Ganopolski, A., Held, H. & Rahmstorf, S. How cold was the Last Glacial Maximum? Geophys. Res. Lett. 33, L14709 (2006).
Ehlers, J. & Gibbard, P. L. Quaternary Glaciations - Extent and Chronology: Part I: Europe Elsevier (2004).
Bini, A. et al. Switzerland during the Last Glacial Maximum (LGM) 1:500000 Federal Office of Topography Swisstopo (2009).
Florineth, D. & Schlüchter, C. Alpine evidence for atmospheric circulation patterns in Europe during the Last Glacial Maximum. Quat. Res. 54, 295â308 (2000).
Kuhlemann, J. et al. Regional synthesis of mediterranean atmospheric circulation during the Last Glacial Maximum. Science 321, 1338â1340 (2008).
Jenk, T. M. et al. A novel radiocarbon dating technique applied to an ice core from the Alps indicating late Pleistocene ages. J. Geophys. Res. 114, D14305 (2009).
Luetscher, M. & Jeannin, P.-Y. Temperature distribution in karst systems: the role of air and water fluxes. Terra Nova 16, 344â350 (2004).
Svensson, A. et al. The Greenland ice core chronology 2005, 15-42 ka. Part 2: Comparison to other records. Quat. Sci. Rev. 25, 3258â3267 (2006).
Rasmussen, S. O. et al. A stratigraphic framework for abrupt climatic changes during the last glacial period based on three synchronized Greenland ice core records: refining and extending the INTIMATE event stratigraphy. Quat. Sci. Rev. 106, 14â28 (2014).
Bond, G. et al. A pervasive millennial-scale cycle in North Atlantic Holocene and Glacial Climates. Science 278, 1257â1266 (1997).
Martrat, B. et al. Four climate cycles of recurring deep and surface water destabilizations on the Iberian Margin. Science 317, 502â507 (2007).
Weinelt, M. et al. Variability of North Atlantic heat transfer during MIS 2. Paleocenanography 18, 1071 (2003).
Monegato, G. et al. Evidence of a two-fold glacial advance during the last glacial maximum in the Tagliamento end moraine system (eastern Alps). Quat. Res. 68, 284â302 (2007).
Preusser, F., Graf, H. R., Keller, O., Krayss, E. & Schlüchter, C. Quaternary glaciation history of northern Switzerland. E&G. Quat. Sci. J. 60, 282â305 (2011).
Wohlfarth, B. et al. Rapid ecosystem response to abrupt climate changes during the last glacial period in western Europe, 40â16 ka. Geology 36, 407â410 (2008).
Martius, O., Zenklusen, E., Schwierz, C. & Davies, HC. Episodes of Alpine heavy precipitation with an overlying elongated stratospheric intrusion: a climatology. Int. J. Climatol. 26, 149â1164 (2006).
Strandberg, G., Brendefelt, J., Kjellström, E. & Smith, B. High-resolution regional simulation of the last glacial maximum climate in Europe. Tellus 63A, 107â125 (2011).
Haeberli, W. & Penz, U. An attempt to reconstruct glaciological and climatological characteristics of 18 ka BP ice age glaciers in and around the Swiss Alps. Z. Gletscherk. Glazialgeol. 21, 351â361 (1985).
Sodemann, H., Schwierz, C. & Wernli, H. Interannual variability of Greenland winter precipitation sources: Lagrangian moisture diagnostic and North Atlantic Oscillation influence. J. Geophys. Res. Atm 113, D3 (2008).
Berger, A. & Loutre, M. F. Insolation values for the climate of the last 10 million years. Quat. Sci. Rev. 10, 297â317 (1991).
Huybers, P. Combined obliquity and precession pacing of the late Pleistocene deglaciation. Nature 480, 229â232 (2011).
Drysdale, R. N. et al. Evidence for obliquity forcing of glacial Termination II. Science 325, 1527â1531 (2009).
Huybers, P. Early Pleistocene glacial cycles and the integrated summer insolation forcing. Science 313, 508â511 (2006).
Lunt, D. J. & Valdes, P. J. Dust deposition and provenance at the Last Glacial Maximum and present day. Geophys. Res. Lett. 29, 2085 (2002).
Shen, C. C. et al. High-precision and high-resolution carbonate 230Th dating by MC-ICP-MS with SEM protocols. Geochim. Cosmochim A 99, 71â86 (2012).
Müller, W., Shelley, M., Miller, P. & Broude, S. Initial performance metrics of a new custom-designed ArF excimer LA-ICPMS system coupled to a two-volume laser-ablation cell. J. Anal. Atom. Spectrom. 24, 209â214 (2009).
Longerich, H. P., Jackson, S. E. & Gunther, D. Laser ablation inductively coupled plasma mass spectrometric transient signal data acquisition and analyte concentration calculation. J. Anal. Atom. Spectrom. 11, 899â904 (1996).
Fohlmeister, J. A statistical approach to construct composite climate records of dated archives. Quatern. Geochron 14, 48â56 (2012).
Bigler, M. Hochauflösende Spurenstoffmessungen an Polaren Eisbohrkernen: Glaziochemische und Klimatische Prozessstudien PhD thesis, University of Bern 148 (2004).
Friedman, I. & O'Neil, J. R. in: Data of Geochemistry 6th edn ed. Fleischer M. 1â12US Geol. Survey Prof. Paper 440-KK (1977).
Acknowledgements
M.L. is recipient of an APART-fellowship of the Austrian Academy of Sciences (ÃAW). Partial support by the Austrian Science Fund (FWF, project P222780 to C.S.) is acknowledged. G. Moseley is thanked for providing an additional U/Th age of 7H-3 sample. C. Baroni, W. Haeberli, M. Sarnthein, A. Svensson and H. Wernli are acknowledged for inspiring discussions.
Author information
Authors and Affiliations
Contributions
M.L. designed the study, completed the field work together with F.H. and interpreted the results. R.B. carried out the U/Th analyses with support of H.C. and R.L.E. H.S. provided the moisture source analysis; C.S. generated the stable isotope data and contributed to the interpretation together with H.S., S.F. helped in the analysis and interpretation of the speleothem petrography and W.M. carried out the trace element analyses. M.L. wrote the manuscript and designed the figures with input from all co-authors.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Supplementary information
Supplementary Information
Supplementary Figures 1-8, Supplementary Tables 1-2, Supplementary Discussion and Supplementary References (PDF 952 kb)
Rights and permissions
This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the articleâs Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/
About this article
Cite this article
Luetscher, M., Boch, R., Sodemann, H. et al. North Atlantic storm track changes during the Last Glacial Maximum recorded by Alpine speleothems. Nat Commun 6, 6344 (2015). https://doi.org/10.1038/ncomms7344
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/ncomms7344
This article is cited by
-
Early expansion of mixed oak stands at 16,800â16,600 cal bp at a northern Italian glacial refugium in the Euganean Hills
Vegetation History and Archaeobotany (2024)
-
The role of Northern Hemisphere summer insolation for millennial-scale climate variability during the penultimate glacial
Communications Earth & Environment (2023)
-
Summer paleohydrology during the Late Glacial and Early Holocene based on δ2H and δ18O from Bichlersee, Bavaria
Scientific Reports (2023)
-
Geoheritage Assessment of a Hermit Cave System at Lakhamandal Valley, Upper Yamuna River Basin, NW Himalaya
Geoheritage (2023)
-
The response of the hydrological cycle to temperature changes in recent and distant climatic history
Progress in Earth and Planetary Science (2022)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.