Arctic lakes located in permafrost regions are susceptible to catastrophic drainage. In this study, we reconstructed historical lake drainage events on the western Arctic Coastal Plain of Alaska between 1955 and 2017 using USGS topographic maps, historical aerial photography (1955), and Landsat Imagery (ca. 1975, ca. 2000, and annually since 2000). We identified 98 lakes larger than 10 ha that partially (>25% of area) or completely drained during the 62‐year period. Decadal‐scale lake drainage rates progressively declined from 2.0 lakes/yr (1955–1975), to 1.6 lakes/yr (1975–2000), and to 1.2 lakes/yr (2000–2017) in the ~30,000‐km² study area. Detailed Landsat trend analysis between 2000 and 2017 identified two years, 2004 and 2006, with a cluster (five or more) of lake drainages probably associated with bank overtopping or headward erosion. To identify future potential lake drainages, we combined the historical lake drainage observations with a geospatial dataset describing lake elevation, hydrologic connectivity, and adjacent lake margin topographic gradients developed with a 5‐m‐resolution digital surface model. We identified ~1900 lakes likely to be prone to drainage in the future. Of the 20 lakes that drained in the most recent study period, 85% were identified in this future lake drainage potential dataset. Our assessment of historical lake drainage magnitude, mechanisms and pathways, and identification of potential future lake drainages provides insights into how arctic lowland landscapes may change and evolve in the coming decades to centuries.

Lakes and drained lake basins (DLBs) together cover up to ∼80% of the western Arctic Coastal Plain of Alaska. The formation and drainage of lakes in this continuous permafrost region drive spatial and temporal landscape dynamics. Postdrainage processes including vegetation succession and permafrost aggradation have implications for hydrology, carbon cycling, and landscape evolution. Here, we used surface nuclear magnetic resonance (NMR) and transient electromagnetic (TEM) measurements in conjunction with thermal modeling to investigate permafrost aggradation beneath eight DLBs on the western Arctic Coastal Plain of Alaska. We also surveyed two primary surface sites that served as nonlake affected control sites. Approximate timing of lake drainage was estimated based on historical aerial imagery. We interpreted the presence of taliks based on either unfrozen water estimated with surface NMR and/or TEM resistivities in DLBs compared to measurements on primary surface sites and borehole resistivity logs. Our results show evidence of taliks below several DLBs that drained before and after 1949 (oldest imagery). We observed depths to the top of taliks between 9 and 45 m. Thermal modeling and geophysical observations agree about the presence and extent of taliks at sites that drained after 1949. Lake drainage events will likely become more frequent in the future due to climate change and our modeling results suggest that warmer and wetter conditions will limit permafrost aggradation in DLBs. Our observations provide useful information to predict future evolution of permafrost in DLBs and its implications for the water and carbon cycles in the Arctic.

Lake formation and drainage are pervasive phenomena in permafrost regions. Drained lake basins (DLBs) are often the most common landforms in lowland permafrost regions in the Arctic (50% to 75% of the landscape). However, detailed assessments of DLB distribution and abundance are limited. In this study, we present a novel and scalable remote sensing-based approach to identifying DLBs in lowland permafrost regions, using the North Slope of Alaska as a case study. We validated this first North Slope-wide DLB data product against several previously published sub-regional scale datasets and manually classified points. The study area covered >71,000 km², including a >39,000 km² area not previously covered in existing DLB datasets. Our approach used Landsat-8 multispectral imagery and ArcticDEM data to derive a pixel-by-pixel statistical assessment of likelihood of DLB occurrence in sub-regions with different permafrost and periglacial landscape conditions, as well as to quantify aerial coverage of DLBs on the North Slope of Alaska. The results were consistent with previously published regional DLB datasets (up to 87% agreement) and showed high agreement with manually classified random points (64.4–95.5% for DLB and 83.2–95.4% for non-DLB areas). Validation of the remote sensing-based statistical approach on the North Slope of Alaska indicated that it may be possible to extend this methodology to conduct a comprehensive assessment of DLBs in pan-Arctic lowland permafrost regions. Better resolution of the spatial distribution of DLBs in lowland permafrost regions is important for quantitative studies on landscape diversity, wildlife habitat, permafrost, hydrology, geotechnical conditions, and high-latitude carbon cycling.

Thermokarst lake dynamics, which play an essential role in carbon release due to permafrost thaw, are affected by various geomorphological processes. In this study, we derive a three-dimensional (3D) Stefan equation to characterize talik geometry under a hypothetical thermokarst lake in the continuous permafrost region. Using the Euler equation in the calculus of variations, the lower bounds of the talik were determined as an extremum of the functional describing the phase boundary area with a fixed total talik volume. We demonstrate that the semi-ellipsoid geometry of the talik is optimal for minimizing the total permafrost thaw under the lake for a given annual heat supply. The model predicting ellipsoidal talik geometry was compared to talik thickness observations using transient electromagnetic (TEM) soundings in Peatball Lake on the Arctic Coastal Plain (ACP) of northern Alaska. The depth : width ratio of the elliptical sub-lake talik can characterize the energy flux anisotropy in the permafrost, although the lake bathymetry cross section may not be elliptic due to the presence of near-surface ice-rich permafrost. This theory suggests that talik development deepens lakes and results in more uniform horizontal lake expansion around the perimeter of the lakes, while wind-induced waves and currents are likely responsible for the elongation and orientation of shallow thermokarst lakes without taliks in certain regions such as the ACP of northern Alaska.

Permafrost warming and degradation is well documented across the Arctic. However, observation‐ and model‐based studies typically consider thaw to occur at 0°C, neglecting the widespread occurrence of saline permafrost in coastal plain regions. In this study, we document rapid saline permafrost thaw below a shallow arctic lake. Over the 15‐year period, the lakebed subsided by 0.6 m as ice‐rich, saline permafrost thawed. Repeat transient electromagnetic measurements show that near‐surface bulk sediment electrical conductivity increased by 198% between 2016 and 2022. Analysis of wintertime Synthetic Aperture Radar satellite imagery indicates a transition from a bedfast to a floating ice lake with brackish water due to saline permafrost thaw. The regime shift likely contributed to the 65% increase in thermokarst lake lateral expansion rates. Our results indicate that thawing saline permafrost may be contributing to an increase in landscape change rates in the Arctic faster than anticipated.