Lake Palcacocha, a moraine-dammed glacial lake above the city of Huaraz, Peru, threatens roughly 120,000 people with a glacial lake outburst flood (GLOF). Stuart-Smith, Roe, Li and Allen (2021) established, in a three-step analysis, that the retreat of Palcaraju glacier — which enlarged the lake and raised the hazard — is with very high confidence attributable to human-induced warming, making Palcacocha the paradigm case for attribution of a site-specific slow-onset hazard. That study itself left a precise open core: the attribution rests on a flowline ice model that omits lake–ice dynamics (calving, thermal conduction), debris cover, and subglacial processes; the pre-industrial (circa 1880) terminus position is an unobserved estimate; the hazard-ranking metrics exclude lake-volume growth and slope-stability drivers such as permafrost degradation and glacier de-buttressing; and the framework has not been demonstrated at any second site. Adjacent global-scale and glacier-retreat attribution advances exist in the 2021–2026 follow-up literature, but no work has closed any of these gaps, so the problem remains open.

Background

A glacial lake outburst flood (GLOF) is the sudden release of water impounded by a glacier or by a moraine — the ridge of debris deposited at a glacier's margin. Lake Palcacocha in the Cordillera Blanca, Peru, sits behind a moraine below Palcaraju glacier and directly above Huaraz, a city of about 120,000 people. In December 1941 a GLOF from the lake destroyed about a third of Huaraz and killed at least 1,800 people; the lake has since regrown and expanded as the glacier retreated. Whether that retreat — and the consequent increase in flood hazard — is attributable to human-induced climate change became the central scientific question of the climate-lawsuit Lliuya v RWE, in which a Huaraz farmer sought a proportional damage contribution from a German emitter; after ten years the Higher Regional Court of Hamm dismissed the claim on 28 May 2025, but the scientific question of how tightly emissions can be tied to a specific site's hazard remains live for future litigation and for adaptation planning.

Glacier attribution builds on a lineage of results. Marzeion, Cogley, Richter and Parkes (Science, 2014) showed from global mass-balance observations and models that the dominant share of modern glacier mass loss is anthropogenic. Roe, Christian and Marzeion (The Cryosphere, 2021) recast single-glacier attribution as a signal-to-noise problem: because a glacier's length integrates climate over its response time, its natural fluctuations amplify the climate signal, so an observed retreat can be compared against the distribution of retreats a glacier would show under purely natural climate variability (a null distribution built from synthetic stochastic climates that match the observed variance and year-to-year persistence of the local climate). Christian, Robel and Catania (The Cryosphere, 2022) extended probabilistic attribution to marine-terminating glaciers. Huggel, Carey, Emmer, Frey, Walker-Crawford and Wallimann-Helmer (NHESS, 2020) analysed the climatic and non-climatic drivers of Palcacocha risk in the litigation context.

Stuart-Smith, Roe, Li and Allen (Nature Geoscience 14, 85–90, 2021) supplied the first site-level attribution of GLOF hazard, in three steps. Step 1: using a global warming index scaled by control simulations of 41 CMIP5 models (CMIP5 being the fifth Coupled Model Intercomparison Project ensemble) and the HadCRUT4-CW regional temperature reconstruction, they found human-induced warming equal to 85–105% (5–95% confidence interval) of the observed warming of about 1 ∘C1\,^\circ\mathrm{C} since 1880. Step 2: they drove a flowline ice-dynamics model — a one-dimensional model along the glacier's centreline using the shallow-ice approximation (ice deforming under its own weight, with geometry averaged across the width) — with the observed and anthropogenic-only temperature series, where warming raises the equilibrium line altitude (the elevation separating net accumulation above from net ablation below). Comparing the modelled terminus evolution with observed 1940–2019 positions and with 10,000 years of synthetic natural variability, they concluded it is virtually certain (>99% probability) that the retreat cannot be explained by natural variability alone, with a central estimate that the retreat is entirely attributable to the warming trend. Step 3: under two established hazard-ranking methodologies, the GLOF hazard category rose from medium in the 19th century to very high today.

The authors themselves enumerated what this attribution does not include. The flowline model omits lake–ice dynamics (ice calving into the lake and thermal conduction between water and ice), debris cover (surface rock debris that insulates ice and modifies melt rates), subglacial till (deformable sediment under the ice) and the possibility that past ice flow altered the bed topography, so the model is not expected to reproduce the transient details of terminus evolution. The pre-industrial (circa 1880) terminus position L0L_{0} was set conservatively at 4,250 m from moraine evidence rather than direct observation. The hazard rankings do not account for lake-volume growth or for dynamic slope-stability processes — permafrost degradation (thawing of permanently frozen ground that weakens rock slopes) and glacier de-buttressing (loss of the lateral support a glacier provides to adjacent rock walls) — both of which raise the probability of avalanches into the lake. Local climate persistence, associated with El Niño variability, and its representation in CMIP5 models carry uncertainty. And the framework, though proposed as transferable, has been demonstrated only at Palcacocha.

Problem Statement

The problem left open by Stuart-Smith et al. (2021) and still open in the literature is to make site-level attribution of anthropogenic climate change influence on GLOF hazard process-complete, uncertainty-quantified, and transferable. Its single overarching objective is a robust quantitative statement of how much of the current outburst-flood hazard at Lake Palcacocha — and, through the same method, at other deglaciating alpine sites — is attributable to human-induced warming, with the physical processes and boundary conditions the original analysis omitted explicitly included. Concretely, the objective requires: (1) incorporate lake–ice dynamics (glacier calving and thermal conduction), debris cover, and subglacial processes (till deformation, past-ice-flow-altered bed topography) into the ice-flow modeling of Palcaraju glacier's retreat, and determine how the transient terminus evolution, the attributable fraction of retreat, and its confidence bounds change relative to the flowline-only analysis of Stuart-Smith et al. (2021); (2) constrain the pre-industrial (circa 1880) terminus position L0L_{0} and its influence on the attribution — including the early-impact claim for the 1941 flood — together with the uncertainty in local climate variability and persistence (for example El Niño-related interannual autocorrelation and its representation in climate models) that sets the natural-variability null distribution; (3) extend the hazard quantification so that changes in triggering probability and flood magnitude account for lake-volume growth and for slope-stability processes — permafrost degradation and glacier de-buttressing — rather than categorical rankings that omit them; and (4) demonstrate the same three-step attribution (regional warming attribution, glacier-retreat attribution against a natural-variability null, hazard quantification) at least one other deglaciating alpine site, handling that site's specific data gaps. Components (1)–(3) jointly deliver the robustness of the Palcacocha attribution statement; component (4) tests whether that statement generalizes as a method. All four are linked by the common objective and are jointly what the source work identifies as necessary to tighten quantitative attributions and hazard projections.

References

  1. Stuart-Smith, R. F., Roe, G. H., Li, S. & Allen, M. R. Increased outburst flood hazard from Lake Palcacocha due to human-induced glacier retreat. Nature Geoscience 14, 85–90 (2021). https://doi.org/10.1038/s41561-021-00686-4
  2. Marzeion, B., Cogley, J. G., Richter, K. & Parkes, D. Attribution of global glacier mass loss to anthropogenic and natural causes. Science 345, 919–921 (2014). https://doi.org/10.1126/science.1254702
  3. Roe, G. H., Christian, J. E. & Marzeion, B. On the attribution of industrial-era glacier mass loss to anthropogenic climate change. The Cryosphere 15, 1889–1905 (2021). https://doi.org/10.5194/tc-15-1889-2021
  4. Christian, J. E., Robel, A. A. & Catania, G. A probabilistic framework for quantifying the role of anthropogenic climate change in marine-terminating glacier retreats. The Cryosphere 16, 2725–2743 (2022). https://doi.org/10.5194/tc-16-2725-2022
  5. Huggel, C., Carey, M., Emmer, A., Frey, H., Walker-Crawford, N. & Wallimann-Helmer, I. Anthropogenic climate change and glacier lake outburst flood risk: local and global drivers and responsibilities for the case of Lake Palcacocha, Peru. Natural Hazards and Earth System Sciences 20, 2175–2193 (2020). https://doi.org/10.5194/nhess-20-2175-2020
  6. Emmer, A., Allen, S., Carey, M., Frey, H., Huggel, C., Korup, O., Mergili, M., Sattar, A., Veh, G., Chen, T., Cook, S., Correas-Gonzalez, M., Das, S., Diaz Moreno, A., Drenkhan, F., Fischer, M., Immerzeel, W., Izagirre, E., Joshi, R., Kougkoulos, I., Kuyakanon Knapp, R., Li, D., Majeed, U., Matti, S., Moulton, H., Nick, F., Piroton, V., Rashid, I., Reza, M., Ribeiro de Figueiredo, A., Riveros, C., Shrestha, F., Shrestha, M., Steiner, J., Walker-Crawford, N., Wood, J. & Yde, J. Progress and challenges in glacial lake outburst flood research (2017–2021): a research community perspective. Natural Hazards and Earth System Sciences 22, 3041–3061 (2022). https://doi.org/10.5194/nhess-22-3041-2022
  7. Lüning, S., Gałka, M., Bamonte, F., García-Rodríguez, F. & Vahrenholt, F. Attribution of modern Andean glacier mass loss requires successful hindcast of pre-industrial glacier changes. Journal of South American Earth Sciences 119, 104024 (2022). https://doi.org/10.1016/j.jsames.2022.104024
  8. Stuart-Smith, R. F., Roe, G. H., Li, S. & Allen, M. R. Comment on 'Attribution of modern Andean glacier mass loss requires successful hindcast of pre-industrial glacier changes' by Sebastian Lüning et al. Journal of South American Earth Sciences 133, 104692 (2024). https://doi.org/10.1016/j.jsames.2023.104692
  9. Lüning, S., Gałka, M., Bamonte, F., García-Rodríguez, F. & Vahrenholt, F. Reply to 'Comment on attribution of modern Andean glacier mass loss requires successful hindcast of pre-industrial glacier changes'. Journal of South American Earth Sciences 133, 104693 (2024). https://doi.org/10.1016/j.jsames.2023.104693
  10. Stansell, N. D., Abbott, M. B., Bezada Díaz, M., Licciardi, J. M., Mark, B. G., Polissar, P. J., Rodbell, D. T. & Shutkin, T. Y. Pre-industrial Holocene glacier variability in the tropical Andes as context for anthropogenically driven ice retreat. Global and Planetary Change 229, 104242 (2023). https://doi.org/10.1016/j.gloplacha.2023.104242
  11. Zhang, T., Wang, W., Kougkoulos, I., Cook, S., Li, S., Iribarren-Anacona, P., Watson, C., An, B. & Yao, T. High frequency of moraine-dammed lake outburst floods driven by global warming. Nature Communications 16, 11173 (2025). https://doi.org/10.1038/s41467-025-67650-3
  12. Veh, G., Wang, B., Zirzow, A., Schmidt, C., Lützow, N., Steppat, F., Zhang, G., Vogel, K., Geertsema, M., Clague, J. J. & Korup, O. Progressively smaller glacier lake outburst floods despite worldwide growth in lake area. Nature Water 3, 271–283 (2025). https://doi.org/10.1038/s44221-025-00388-w
  13. Ochoa-Sánchez, A., Stone, D., Drenkhan, F., Mendoza, D., Gualán, R. & Huggel, C. Detection and attribution of climate change impacts in coupled natural-human systems in the Andes. Communications Earth & Environment 6, 314 (2025). https://doi.org/10.1038/s43247-025-02092-9
  14. Undorf, S., Dimitrova, A., Harrington, L. J., Hegerl, G. C., Jézéquel, A., Kimutai, J., Lo, Y. T. E., Mengel, M., Perkins-Kirkpatrick, S., Pietroiusti, R., Pinto, I., Schurer, A., Scussolini, P., Sippel, S., Stone, D., Stuart-Smith, R., Gudmundsson, L., Huber, V., Colón-González, F. J., Dellicour, S., Erazo, D., Fischer, E., Frieler, K., Gornott, C., Hawkins, E., Holden, P., Koren, G., Lampe, S., Mitchell, D., Murken, L., New, M., Nkwasa, A., Seneviratne, S. I., Wehner, M. & Thiery, W. Approaches, challenges and applications of climate change impact attribution. Nature Reviews Earth & Environment 7, 469–488 (2026). https://doi.org/10.1038/s43017-026-00798-8
  15. Bradley, A. T., Bett, D. T., Williams, C. R., Arthern, R. J., Holland, P. R., Byrne, J., Edwards, T. & Adhikari, M. Detection and attribution of the role of anthropogenic climate change in industrial-era retreat of Pine Island Glacier. The Cryosphere 20, 3443–3465 (2026). https://doi.org/10.5194/tc-20-3443-2026
  16. Wentz, J., Merner, D., Franta, B., Lehmen, A. & Frumhoff, P. Research Priorities for Climate Litigation. Earth's Future 11, e2022EF002928 (2022). https://doi.org/10.1029/2022EF002928
  17. Nogrady, B. Securing climate justice in the courtroom. Nature 644, S1–S3 (2025). https://doi.org/10.1038/d41586-025-02623-6