As global temperatures continue their upward trajectory due to anthropogenic climate change, one of the most conspicuous impacts is the rapid melting of glaciers, which historically have been key repositories of fresh water. However, beyond the visible retreat of ice, the terrain that lies beneath these glaciers remains largely uncharted in fine detail. The team’s study ingeniously uses high-resolution digital elevation models combined with glacier outlines to reconstruct the exposed landscape beneath current ice masses, thereby revealing a topography primed for lake formation once the ice recedes.

The researchers’ methodology involved compiling glacier-free digital elevation models (DTMs) for all ice-covered terrains globally. These DTMs negate the presence of ice, effectively depicting the “bare Earth” beneath. This approach is innovative because it circumvents the traditional limitations encountered in glaciated regions where surface ice distorts or obscures bedrock features essential for predicting post-glacial hydrological patterns. By overlaying the glacier-free topography with existing glacier extents, the team could simulate the potential basins where water might concentrate.

Their results are astonishing in scale. The study identifies thousands of previously unrecognized depressions and basins capable of becoming proglacial lakes in the near future. These potential lake sites are not randomly distributed but cluster distinctly in certain glaciated regions, particularly in the Himalayas, the Andes, Alaska, and parts of the Arctic. This highlights the varying regional vulnerabilities to newly emerging water bodies, which could profoundly alter local environments and human activities.

Crucially, the study underscores how the formation of such lakes could modify downstream hydrology, impacting river discharge regimes and sediment transport. As glaciers withdraw, the newly formed lakes could act as buffers, storing meltwater temporarily but also posing risks of sudden outburst floods should their natural dams fail. Predicting the locations, sizes, and volumes of these lakes is therefore essential for risk assessment and adaptation planning in vulnerable mountainous communities.

Moreover, the work highlights the ecological ramifications of these new aquatic environments. Emerging lakes will create novel habitats, potentially fostering biodiversity and altering existing ecosystems abruptly. Species adapted to cold glacier-fed streams may face habitat fragmentation or loss, while new aquatic niches could emerge, challenging prevailing conservation strategies. Scientists and policymakers must anticipate these ecological cascades to mitigate negative impacts on biodiversity and ecosystem services.

The researchers emphasize the dynamic nature of glacier retreat and lake formation processes; these phenomena are not instantaneous but evolve over decades or even centuries. However, with current acceleration rates in glacier mass loss, the transition from ice-covered terrain to lake-dotted landscapes might proceed rapidly in several hotspots around the globe. This temporal insight gives urgency to mapping and monitoring efforts as hydrological infrastructures may face profound transformations within a generation.

In addition to the environmental and ecological consequences, the study provides a crucial tool for water resource management. In many regions, glacier meltwater supplements river flow during dry seasons, supporting agriculture and human consumption. The formation of stable lakes could modulate seasonal water availability, potentially buffering against drought or conversely complicating water distribution networks depending on lake longevity and outflow patterns.

The availability of detailed glacier-free topography also serves as a valuable baseline to improve climate model projections and glacio-hydrological simulations. Prior models often estimated lake potential based on incomplete or ice-covered terrain data, limiting their predictive accuracy. The new dataset, openly shared by the authors, will facilitate enhanced integrated assessments, fostering interdisciplinary collaboration among climate scientists, hydrologists, ecologists, and hazard modelers.

Interestingly, the paper also delves into the geological implications of deglaciation-induced lake formation. Newly exposed terrain beneath ice can possess complex bedrock structures, potentially influencing sedimentation processes within nascent lakes. Moreover, post-glacial rebound, the gradual uplift of land previously depressed by ice weight, may alter drainage paths and lake morphology over time, introducing another layer of complexity in predicting lake evolution.

The global scope of the study sets it apart from previous regional analyses, providing a holistic view of glacial landscapes transitioning to proglacial lake systems worldwide. This comprehensive vantage point is invaluable as it exposes global patterns and regional disparities in lake formation potential driven by diverse climatic, topographic, and geological conditions.

While the study represents a significant leap forward, the authors acknowledge limitations, such as the challenges in precisely modeling small-scale terrain features beneath extraordinarily thick ice and uncertainties associated with future climate scenarios. They advocate for continuous refinement of remote sensing technologies and ground-based observations to enhance resolution and validation of glacier-free topographies.

The implications of this research extend beyond pure science. For communities residing downstream of shrinking glaciers, early warning systems for glacier lake outburst floods will become increasingly important. Likewise, hydroelectric projects relying on glacier-fed reservoirs may need to revise risk assessments considering possible sudden lake formations or changes in meltwater dynamics.

In conclusion, Frank et al.’s innovative mapping of global glacier-free topography offers an unprecedented glimpse into the landscapes soon to emerge as glaciers fade. Their findings illuminate an extensive potential for future lakes that will reshape hydrology, ecosystems, and hazards in mountain and polar environments. This study not only advances scientific frontiers but also equips society with crucial knowledge to navigate the complex transition ahead in a warming world. Vigilant monitoring, integrated modeling, and proactive adaptation strategies will be essential to harness opportunities and mitigate risks posed by the new lakes born from disappearing ice.

Subject of Research: Global glacier retreat and potential lake formation in deglaciated terrain

Article Title: Global glacier-free topography reveals a large potential for future lakes in presently ice-covered terrain

Article References:
Frank, T., van Pelt, W.J.J., Rounce, D.R. et al. Global glacier-free topography reveals a large potential for future lakes in presently ice-covered terrain. Nat Commun 17, 3985 (2026). https://doi.org/10.1038/s41467-026-72548-9

Image Credits: AI Generated

DOI: https://doi.org/10.1038/s41467-026-72548-9

Summary of Synthesis of Glacial Erosion Rates and Related Data:

1. Data Collection and Scope:

• Glacial erosion rates were synthesized from published data primarily measuring the removal of material from glacier beds (164 contemporary rates + 17 Holocene rates).
• Only data from topographically constrained extant glaciers and ice caps with dominantly warm basal thermal regimes were included.
• Measurements come from four methods:
  • Sediment discharge from meltwater streams.
  • Bulk sediment in proglacial areas.
  • Bedrock erosion via terrestrial cosmogenic nuclide dating.
  • Instrumental measurements from extant glaciers.
• Multiple measurements from the same glacier were included to reflect the range of observed erosion rates.

2. Glaciological Variables:

• Variables like ice mass length, area, median elevation, surface slope, glacier termination type, form, surge behaviour were sourced primarily from the Randolph Glacier Inventory (RGI25).
• Missing variables were supplemented by digital elevation models (SRTM and ASTER GDEM).
• Surface and sliding velocities and ice thicknesses were included where available, preferring multi-year averages temporally matching erosion measurements.
• Where velocity/ice thickness not available from original sources, they were estimated from published velocity datasets or ice thickness models.
• Sliding velocity (U_sliding) was estimated from surface speed (U_surface) and local surface slope (β) using a constrained ratio and linear interpolation.

3. Geological Variables:

• Glacier bed lithology was extracted from the Global Lithological Map (GLiM) using glacier centroids.
• Lithologies categorized into:
  • Crystalline (igneous and metamorphic)
  • Sedimentary (carbonates, sandstones, shales)
  • Volcanics (volcanics, basalts)
• Seismicity approximated using PGA (Peak Ground Acceleration) from the Global Seismic Hazard Assessment dataset, filling gaps for polar regions where needed.

4. Climatic Variables:

• Mean Annual Air Temperature (MAAT) and Mean Annual Precipitation (MAP) were extracted from the CHELSA climate dataset for modern erosion rates and from PaleoClim for Holocene rates, using glacier centroids.

5. Dataset and Analysis Limitations:

• Methods differ in the processes they capture, potentially biasing erosion rates (e.g., instrumental measurements focus on abrasion only).
• Bulk sediment and sediment discharge methods might underestimate or misattribute erosion due to sediment storage effects or partial sediment load measurement.
• Timescale dependence is mitigated by focusing on short-term (~350 years or less) erosion rates; no significant correlation between rate and timespan was observed.
• Temporal mismatches between velocity and erosion rate measurements remain a limitation, though methods to mitigate this were employed.
• Velocity data were selected preferentially to match erosion measurement periods, with sensitivity analyses confirming this approach improved model performance.
• Surge-type glaciers were underrepresented (<15% of data), which may limit extrapolation for those glacier types.

Additional Notes:

• Velocity estimation: sliding velocity is a proportion of surface speed scaled by slope β constrained between 0.001 and 0.1, interpolated linearly.
• Velocity variability and surge behaviour introduce further complexities.
• Data synthesis incorporated multiple data quality and spatial resolution sources to achieve as comprehensive and consistent a dataset as possible.

Application of This Dataset:

• Enables near-global estimation of glacial erosion rates.
• Can be used to investigate relationships between glacial erosion and glaciological, geological, climatic, and seismic variables.
• Provides a basis for “velocity–erosion” rules and erosion scaling models relevant to Earth surface processes, landscape evolution, and glaciology.

If you have any specific questions regarding this synthesis, such as the methods for velocity estimation, geological categorizations, or limitations on data interpretation, please let me know!