Abstract
Terrestrial evidence of ice-sheet meltwater events during the last deglaciation are exceptionally rare. However, they are a valuable complement to marine records that suggest immense ice-sheet meltwater pulses discharging into the Northeast Pacific, termed Siku events. Here we describe and date a terrestrial Cordilleran Ice Sheet meltwater archive in a karstified glacial outwash plain during the early last deglaciation, showing semi-continuous meltwater events between about 20,000 and 17,000 years ago.
This is a preview of subscription content, access via your institution
Access options
Access Nature and 54 other Nature Portfolio journals
Get Nature+, our best-value online-access subscription
$29.99 / 30 days
cancel any time
Subscribe to this journal
Receive 12 print issues and online access
$259.00 per year
only $21.58 per issue
Buy this article
- Purchase on SpringerLink
- Instant access to full article PDF
Prices may be subject to local taxes which are calculated during checkout


Similar content being viewed by others
Data availability
Devil’s Sequence OSL and grain-size data are available from the NOAA National Centers for Environmental Information website at https://doi.org/10.25921/y45g-ps62 (ref. 29).
References
Denton, G. H. et al. The last glacial termination. Science 328, 1652–1656 (2010).
Carlson, A. E. & Winsor, K. Northern Hemisphere ice-sheet responses to past climate warming. Nat. Geosci. 5, 607–613 (2012).
Abe-Ouchi, A. et al. Insolation-driven 100,000-year glacial cycles and hysteresis of ice-sheet volume. Nature 500, 190–193 (2013).
Robinson, A. & Goelzer, H. The importance of insolation changes for paleo ice sheet modelling. Cryosphere 8, 1419–1428 (2014).
Broecker, W. S. Massive iceberg discharges as triggers for global climate change. Nature 372, 421–424 (1994).
Walczak, M. H. et al. Phasing of millennial-scale climate variability in the Pacific and Atlantic Oceans. Science 370, 716–720 (2020).
Lesnek, A. J., Briner, J. P., Baichtal, J. F. & Lyles, A. S. New constraints on the last deglaciation of the Cordilleran Ice Sheet in coastal Southeast Alaska. Quat. Res. 96, 140–160 (2020).
Dalton, A. S. et al. Deglaciation of the North American ice sheet complex in calendar years based on a comprehensive database of chronological data: NADI-1. Quat. Sci. Rev. 321, 108345 (2023).
Eyles, N. The role of meltwater in glacial processes. Sediment. Geol. 190, 257–268 (2006).
Wilson, F. H., Hults, C P., Mull, C. G. & Karl, S. M. Geologic Map of Alaska Scientific Investigations Map 3340 (US Geological Survey, 2015).
Hansen, B. C. S. & Engstrom, D. R. Vegetation history of Pleasant Island, southeastern Alaska, since 13,000 yr B.P. Quat. Res. 46, 161–175 (1996).
Wilcox, P. S., Fowell, S. J. & Baichtal, J. F. A mild Younger Dryas recorded in southeastern Alaska. Arct. Antarct. Alp. Res. 52, 236–247 (2020).
Lacourse, T., Mathewes, R. W. & Fedje, D. W. Late-glacial vegetation dynamics of the Queen Charlotte Islands and adjacent continental shelf, British Columbia, Canada. Palaeogeogr. Palaeoclimatol. Palaeoecol. 226, 36–57 (2005).
Praetorius, S. K. et al. The role of Northeast Pacific meltwater events in deglacial climate change. Sci. Adv. 6, eaay2915 (2020).
Cowan, E. A. et al. Sediment controls dynamic behavior of a Cordilleran Ice Stream at the Last Glacial Maximum. Nat. Commun. 11, 1826 (2020).
Lesnek, A. J., Briner, J. P., Lindqvist, C., Baichtal, J. F. & Heaton, T. H. Deglaciation of the Pacific coastal corridor directly preceded the human colonization of the Americas. Sci. Adv. 4, eaar5040 (2018).
Wilcox, P. S., Spötl, C., Honkonen, J. & Edwards, R. L. A Walker switch mechanism driving millennial-scale climate variability. Innov. Geosci. 1, 100026 (2023).
Zhang, S. et al. Thermal coupling of the Indo-Pacific warm pool and Southern Ocean over the past 30,000 years. Nat. Commun. 13, 5457 (2022).
MacAyeal, D. R. Binge/purge oscillations of the Laurentide ice sheet as a cause of the North Atlantic’s Heinrich events. Paleoceanogr. Paleoclimatol. 8, 775–784 (1993).
Berger, A. L. Long-term variations of daily insolation and Quaternary climatic changes. J. Atmos. Sci. 35, 2362–2367 (1978).
Bøtter-Jensen, L., Andersen, C. E., Duller, G. A. T. & Murray, A. S. Developments in radiation, stimulation and observation facilities in luminescence measurements. Radiat. Meas. 37, 535–541 (2003).
Murray, A. et al. Optically stimulated luminescence dating using quartz. Nat. Rev. Methods Primers 1, 72 (2021).
Murray, A. S. & Wintle, A. G. Luminescence dating of quartz using an improved single-aliquot regenerative-dose protocol. Radiat. Meas. 32, 57–73 (2000).
Cunningham, A. C. & Wallinga, J. Selection of integration time intervals for quartz OSL decay curves. Quat. Geochronol. 5, 657–666 (2010).
Galbraith, R. F., Roberts, R. G., Laslett, G. M., Yoshida, H. & Olley, J. M. Optical dating of single and multiple grains of quartz from Jinmium rock shelter, northern Australia. I. Experimental design and statistical models. Archaeometry 41, 339–364 (1999).
Faegri, K. & Iversen, J. Textbook of Pollen Analysis (Wiley, 1989)
Stockmarr, J. Tables with spores used in absolute pollen analysis. Pollen Spores 13, 615–621 (1971).
Juggins, S. C2 Data Analysis v.1.4.2. Software for ecological and palaeoecological data analyses and visualization (Newcastle Univ., 1991).
Wilcox, P. S. & Meyer, M. NOAA/WDS Paleoclimatology – Devil’s Canopy Cave (Alaska) sediment analysis from 17–19 ka. NOAA National Centers for Environmental Information https://doi.org/10.25921/y45g-ps62 (2025).
Acknowledgements
This work was funded by Austrian Science Fund (FWF) grant FP338960. We are grateful to A. Harris, J. Baichtal, C. DeCelle and the Tongass National Forest geology program for their continued support for this work. In addition, we thank J. Honkonen for invaluable field-work assistance and T. Ivanova for precious support in preparation of the pollen samples. Finally, we thank C. Spötl and J. Monroe for providing valuable feedback from an earlier version of this paper.
Author information
Authors and Affiliations
Contributions
Conceptualization: P.S.W.; funding acquisition of (FWF) grant FP338960: P.S.W.; project administration/supervision: P.S.W.; methodology design: P.S.W.; conducted OSL dating: M.C.M.; conducted pollen analysis: D.F.; prepared visual components: P.S.W. and M.C.M.; writing–original draft: P.S.W.; writing–reviewing and editing: P.S.W., M.C.M. and D.F.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature Geoscience thanks Maureen Walczak and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Primary Handling Editor: Alison Hunt, in collaboration with the Nature Geoscience team.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Extended data
Extended Data Fig. 1 Profile map of Devils Canopy Cave.
Blue arrows denote modern-day streamflow. Water from modern-day stream exits Devils Canopy into a large network of other caves, some being unexplored. Red star denotes location of Devils sequence in Devils Canopy Room. Also shown is the boulder-sealed passage where sediments from Devils Sequence are hypothesized to have originated before being blocked. Green dashed arrow shows possible pathway of former sediment-laden ice-sheet meltwater in DC Room. Map modified from original map by Kevin and Carlene Allred.
Extended Data Fig. 2 Illustration of Devils Sequence from Devils Canopy Cave.
Unit 1 is composed of diamicton while unit 2 is composed of normally graded bedding of fine/medium sand. Red circles denote OSL age sampling locations with corresponding sample names and ages.
Extended Data Fig. 3 Images of DC Room and Devils Sequence.
(a) Devils Sequence exposed with human for scale. (b) The possible former ice-sheet meltwater pathway in DC Room (red circle). (c) Selected segment of unit 1 composed of diamicton. (d) Selected segment of unit 2 composed of normally graded bedding of find/medium sand.
Extended Data Fig. 4 Images of Devils Canopy Cave surroundings.
(a) Thickly vegetated temperate rainforest in the vicinity of Devils Canopy Cave with Entrance 1 (see Extended Data Fig. 1) located ~2 m to the left of people (red X). (b) Modern-day stream into Devils Canopy Cave (red X).
Extended Data Fig. 5 Pollen diagram of Devils Sequence.
Absolute values of pollen, spores and algae are shown. Total pollen and fern spores concentration is shown in black.
Extended Data Fig. 6 Examples of well-preserved microfossils found in the Devils Sequence.
(a) Picea cf. sitchensis; (b) Pinus cf. contorta; (c) Betula; (d) Alnus; (e) Chenopodiaceae; (f) Fern monolete spore; (g) Algae (cf. snow algae).
Extended Data Fig. 7 Sediments from other caves in the region showing similarities with Devils Sequence from Devils Canopy Cave.
(a) Raven Cave, ~20 km northeast Devils Canopy Cave. (b) El Captain Cave, ~10 km northwest of Devils Canopy Cave (Fig. 1).
Supplementary information
Supplementary Information
Supplementary information on OSL dating methods/results of Devil’s Canopy Cave sediment.
Rights and permissions
Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
About this article
Cite this article
Wilcox, P.S., Meyer, M.C. & Festi, D. Semi-continuous release of Cordilleran Ice Sheet meltwater between 20,000 and 17,000 years ago. Nat. Geosci. (2025). https://doi.org/10.1038/s41561-025-01694-4
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41561-025-01694-4