skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


This content will become publicly available on March 29, 2026

Title: TOWARDS A COMPREHENSIVE ANALYSIS OF LATE-GLACIAL LAKE LEVELS WITHIN THE LAKE SUPERIOR BASIN
The increasing availability of LiDAR digital elevation models (DEMs) for Canada and the United States is improving our ability to recognize and detail the geomorphological record of former glacial lakes revealed by strandlines. More accurate strandline maps refine our understanding of glacial-isostatic adjustment (GIA), reveal information about deglacial ice margins, help identify rapid drawdown events and catastrophic flooding by ice dam failure, and advance models that connect waterplane curvature to age. Here we present new strandline maps in the Lake Superior Basin based on LiDAR DEMs from Wisconsin, Michigan, and Ontario. When combined with published records from Minnesota, they comprise the most extensive set of strandlines detailing former lake levels from any lake basin. Multiple lakes and lake phases are associated with these strandlines, including glacial lakes Wrenshall, Ontonagon, Duluth, Algonquin and Minong, as well as Holocene Lake Nipissing. The relative degree of strandline development varies widely, likely due to factors such as the availability of gravel and sand, length of time of formation, fetch, and local shoreline geometry (e.g. bay or headland). Strandlines for at least one level of glacial Lake Duluth are much more extensive on the southern shore, suggesting formation from only a very few numbers of storms with winds from the north. New LiDAR DEMs from the Sault Ste. Marie region are particularly helpful because they allow us to determine isobases (lines of equal GIA) across the primary basin outlet with high confidence. In general, our work supports the work of prior researchers, but adds details that are easily missed by field surveys or when using coarser resolution DEMs.  more » « less
Award ID(s):
2218464
PAR ID:
10618294
Author(s) / Creator(s):
; ;
Publisher / Repository:
2025 Joint Northeastern and North-Central Section Meeting
Date Published:
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract As ice sheets load Earth's surface, they produce ice‐marginal depressions which, when filled with meltwater, become proglacial lakes. We include self‐consistently evolving proglacial lakes in a glacial isostatic adjustment (GIA) model and apply it to the Laurentide ice sheet over the last glacial cycle. We find that the locations of modeled lakes and the timing of their disappearance is consistent with the geological record. Lake loads can deflect topography by >10 m, and volumes collectively approach 30–45 cm global mean sea‐level equivalent. GIA increases deglaciation‐phase lake volume up to five‐fold and average along‐ice‐margin depth ≤90 m compared to glaciation‐phase ice volume analogs—differences driven by changes in the position and size of the peripheral bulge. Since ice‐marginal lake depth affects grounding‐line outflow, GIA‐modulated proglacial lake depths could affect ice‐sheet mass loss. Indeed, we find that Laurentide ice‐margin retreat rate sometimes correlates with proglacial lake presence, indicating that proglacial lakes aid glacial collapse. 
    more » « less
  2. Oneida Lake, New York, is the remnant of Glacial Lake Iroquois, a large proglacial lake that delivered fresh water to the Atlantic Ocean during the last deglaciation. The formation of Glacial Lake Iroquois and its subsequent drainage into the Atlantic Ocean via the Mohawk Valley was a significant shift in the routing of Laurentide Ice Sheet meltwater to the east instead of south via the Allegheny or Susquehanna Rivers. Catastrophic drainage of Glacial Lake Iroquois into the Atlantic Ocean via the Champlain Valley is interpreted as the meltwater pulse responsible for the Intra-Allerod cold stadial. Therefore, understanding the evolution of Glacial Lake Iroquois has significant implications for understanding late Pleistocene paleoclimate. High-resolution CHIRP seismic reflection data provides insight into the evolution of Glacial Lake Iroquois and Oneida Lake. Three seismic units image distinct stages of the Oneida Basin. Unit 1 is interpreted as proglacial lake deposits that overlie glacial till. Unit 2 is interpreted as sediments deposited when the Oneida Basin became isolated from Glacial Lake Iroquois and Unit 3 is interpreted as lacustrine sediments of the modern lake. Distally sourced turbidites possibly triggered by seismic activity or ice sheet meltwater pulses are represented as reflection- free acoustic facies that infill topographic lows and range in thickness from ~1–5m within otherwise conformable proglacial lake deposits. Local slump deposits imaged at the boundary between Unit 1 and 2 were likely triggered by the drainage of Glacial Lake Iroquois. Wave cut terraces indicative of a low stand on the upper bounding surface of Unit 2 are likely the result of drier conditions during the Holocene Hypsithermal. Furthermore, preservation of this low stand suggests a rapid rise in lake level, possibly the result of the same transition to a wetter climate responsible for the Nipissing transgression observed in the Laurentian Great lakes. 
    more » « less
  3. Abstract Proglacial lakes, whose numbers have been growing around the world, may drive accelerated glacier retreat and provide valuable records of past glacier and climatic changes. Despite their importance, few studies have investigated the sedimentary properties and processes acting within large proglacial lakes. Lago Argentino (LArg) is a 1,500 km2ice‐contact lake on the eastern flank of the Southern Patagonian Icefield. Here, we describe the results from a detailed analysis of 47 sediment cores obtained throughout this lake basin, supplemented with remotely sensed data. We show that: (a) LArg exhibits a seasonal variation in sediment properties (varves); (b) varve formation results from three distinct processes, driven by seasonal changes in glacial sediment input, seasonal changes in fluvial sediment input, and seasonal variations in lake mixing; and (c) distance from glacier calving fronts provides the first‐order control on sediment grain size and accumulation rate. Our findings highlight the exceptional preservation of annual laminations within proglacial lakes, their potential for reconstructing past glacier changes, and their relevance for forecasting future glacier–lake interactions. 
    more » « less
  4. The glaciers of the Cordillera Blanca, Peru, are rapidly retreating and thinning as a result of climate change, altering the timing, quantity and quality of water available to downstream users. Furthermore, increases in the number and size of proglacial lakes associated with these melting glaciers is increasing potential exposure to glacier lake outburst floods (GLOFs). Understanding how these glaciers are changing and their connection to proglacial lake systems is thus of critical importance. Most satellite data are too coarse for studying small mountain glaciers and are often affected by cloud cover, while traditional airborne photogrammetry and lidar are costly. Recent developments have made unmanned aerial vehicles (UAVs) a viable and potentially transformative method for studying glacier change at high spatial resolution, on demand and at relatively low cost.Using a custom designed hexacopter built for high-altitude (4000–6000 m a. s. l. ) operation, we completed repeat aerial surveys (2014 and 2015) of the debris-covered Llaca Glacier tongue and proglacial lake system. High-resolution orthomosaics (5 cm) and digital elevation models (DEMs) (10 cm) were produced and their accuracy assessed. Analysis of these datasets reveals highly heterogeneous patterns of glacier change. The most rapid areas of ice loss were associated with exposed ice cliffs and meltwater ponds on the glacier surface. Considerable subsidence and low surface velocities were also measured on the sediments within the pro-glacial lake, indicating the presence of extensive regions of buried ice and continued connection to the glacier tongue. Only limited horizontal retreat of the glacier tongue was observed, indicating that measurements of changes in aerial extent alone are inadequate for monitoring changes in glacier ice quantity. 
    more » « less
  5. Nearly 25% of all lakes on earth are located at high latitudes. These lakes are formed by a combination of thermokarst, glacial, and geological processes. Evidence suggests that the origin of periglacial lake formation may be an important factor controlling the likelihood of lakes to drain. However, geospatial data regarding the spatial distribution of these dominant Arctic and subarctic lakes are limited or do not exist. Here, we use lake-specific morphological properties using the Arctic Digital Elevation Model (DEM) and Landsat imagery to develop a Thermokarst lake Settlement Index (TSI), which was used in combination with available geospatial datasets of glacier history and yedoma permafrost extent to classify Arctic and subarctic lakes into Thermokarst (non-yedoma), Yedoma, Glacial, and Maar lakes, respectively. This lake origin dataset was used to evaluate the influence of lake origin on drainage between 1985 and 2019 in northern Alaska. The lake origin map and lake drainage datasets were synthesized using five-year seamless Landsat ETM+ and OLI image composites. Nearly 35,000 lakes and their properties were characterized from Landsat mosaics using an object-based image analysis. Results indicate that the pattern of lake drainage varied by lake origin, and the proportion of lakes that completely drained (i.e., >60% area loss) between 1985 and 2019 in Thermokarst (non-yedoma), Yedoma, Glacial, and Maar lakes were 12.1, 9.5, 8.7, and 0.0%, respectively. The lakes most vulnerable to draining were small thermokarst (non-yedoma) lakes (12.7%) and large yedoma lakes (12.5%), while the most resilient were large and medium-sized glacial lakes (4.9 and 4.1%) and Maar lakes (0.0%). This analysis provides a simple remote sensing approach to estimate the spatial distribution of dominant lake origins across variable physiography and surficial geology, useful for discriminating between vulnerable versus resilient Arctic and subarctic lakes that are likely to change in warmer and wetter climates. 
    more » « less