skip to main content

Title: Recent warming reverses forty-year decline in catastrophic lake drainage and hastens gradual lake drainage across northern Alaska

Lakes represent as much as ∼25% of the total land surface area in lowland permafrost regions. Though decreasing lake area has become a widespread phenomenon in permafrost regions, our ability to forecast future patterns of lake drainage spanning gradients of space and time remain limited. Here, we modeled the drivers of gradual (steady declining lake area) and catastrophic (temporally abrupt decrease in lake area) lake drainage using 45 years of Landsat observations (i.e. 1975–2019) across 32 690 lakes spanning climate and environmental gradients across northern Alaska. We mapped lake area using supervised support vector machine classifiers and object based image analyses using five-year Landsat image composites spanning 388 968 km2. Drivers of lake drainage were determined with boosted regression tree models, using both static (e.g. lake morphology, proximity to drainage gradient) and dynamic predictor variables (e.g. temperature, precipitation, wildfire). Over the past 45 years, gradual drainage decreased lake area between 10% and 16%, but rates varied over time as the 1990s recorded the highest rates of gradual lake area losses associated with warm periods. Interestingly, the number of catastrophically drained lakes progressively decreased at a rate of ∼37% decade−1from 1975–1979 (102–273 lakes draining year−1) to 2010–2014 (3–8 lakes draining more » year−1). However this 40 year negative trend was reversed during the most recent time-period (2015–2019), with observations of catastrophic drainage among the highest on record (i.e. 100–250 lakes draining year−1), the majority of which occurred in northwestern Alaska. Gradual drainage processes were driven by lake morphology, summer air and lake temperature, snow cover, active layer depth, and the thermokarst lake settlement index (R2adj= 0.42, CV = 0.35,p< 0.0001), whereas, catastrophic drainage was driven by the thawing season length, total precipitation, permafrost thickness, and lake temperature (R2adj= 0.75, CV = 0.67,p< 0.0001). Models forecast a continued decline in lake area across northern Alaska by 15%–21% by 2050. However these estimates are conservative, as the anticipated amplitude of future climate change were well-beyond historical variability and thus insufficient to forecast abrupt ‘catastrophic’ drainage processes. Results highlight the urgency to understand the potential ecological responses and feedbacks linked with ongoing Arctic landscape reorganization.

« less
; ;
Award ID(s):
1806213 1929170 1928048
Publication Date:
Journal Name:
Environmental Research Letters
Page Range or eLocation-ID:
Article No. 124019
IOP Publishing
Sponsoring Org:
National Science Foundation
More Like this
  1. 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%more »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.« less
  2. Abstract

    Lakes are an important ecosystem component and geomorphological agent in northern high latitudes and it is important to understand how lake initiation, expansion and drainage may change as high latitudes continue to warm. In this study, we utilized Landsat Multispectral Scanner System images from the 1970s (1972, 1974, and 1975) and Operational Land Imager images from the 2010s (2013, 2014, and 2015) to assess broad-scale distribution and changes of lakes larger than 1 ha across the four permafrost zones (continuous, discontinuous, sporadic, and isolated extent) in western Alaska. Across our 68 000 km2study area, we saw a decline in overall lake coverage across all permafrost zones with the exception of the sporadic permafrost zone. In the continuous permafrost zone lake area declined by −6.7% (−65.3 km2), in the discontinuous permafrost zone by −1.6% (−55.0 km2), in the isolated permafrost zone by −6.9% (−31.5 km2) while lake cover increased by 2.7% (117.2 km2) in the sporadic permafrost zone. Overall, we observed a net drainage of lakes larger than 10 ha in the study region. Partial drainage of these medium to large lakes created an increase in the area covered by small water bodies <10 ha, in the form of remnant lakesmore »and ponds by 7.1% (12.6 km2) in continuous permafrost, 2.5% (15.5 km2) in discontinuous permafrost, 14.4% (74.6 km2) in sporadic permafrost, and 10.4% (17.2 km2) in isolated permafrost. In general, our observations indicate that lake expansion and drainage in western Alaska are occurring in parallel. As the climate continues to warm and permafrost continues to thaw, we expect an increase in the number of drainage events in this region leading to the formation of higher numbers of small remnant lakes.

    « less
  3. Abstract
    Forty-five years (i.e. 1975-2019) of Landsat observations were used to map the spatiotemporal patterns of lake drainage in northern Alaska. All Landsat data was pre-processed by the United States Geological Survey and downloaded by google earth engine in a radiometrically, atmospherically, and geometrically terrain-corrected state. We used Landsat surface reflectance products acquired from the Multispectral Scanner (MSS), Terrestrial Mapper (TM), Enhanced Terrestrial Mapper Plus (ETM&#43;), and Operational Land Imager (OLI) sensors to compute eight image mosaics for the ice-free period (June 15 to September 1) at five-year time-periods. This data product represents the change in lake area between time-periods or epochs. All data used to spatially identify patterns of lake change are presented in a map (LakeDrainChg.tif), where the associated morphometric controls on drainage are summarized for each of the ~33,000 lake boundaries (Lake_Drainage.shp). All data can also be viewed within Google Earth Engine (;users/mjlara71/LakeDrainage; accessed on 13 September 2021) or clone Git repository (git clone ; accessed on 13 September 2021).
  4. Abstract

    Thermokarst lakes accelerate deep permafrost thaw and the mobilization of previously frozen soil organic carbon. This leads to microbial decomposition and large releases of carbon dioxide (CO2) and methane (CH4) that enhance climate warming. However, the time scale of permafrost-carbon emissions following thaw is not well known but is important for understanding how abrupt permafrost thaw impacts climate feedback. We combined field measurements and radiocarbon dating of CH4ebullition with (a) an assessment of lake area changes delineated from high-resolution (1–2.5 m) optical imagery and (b) geophysical measurements of thaw bulbs (taliks) to determine the spatiotemporal dynamics of hotspot-seep CH4ebullition in interior Alaska thermokarst lakes. Hotspot seeps are characterized as point-sources of high ebullition that release14C-depleted CH4from deep (up to tens of meters) within lake thaw bulbs year-round. Thermokarst lakes, initiated by a variety of factors, doubled in number and increased 37.5% in area from 1949 to 2009 as climate warmed. Approximately 80% of contemporary CH4hotspot seeps were associated with this recent thermokarst activity, occurring where 60 years of abrupt thaw took place as a result of new and expanded lake areas. Hotspot occurrence diminished with distance from thermokarst lake margins. We attribute older14C ages of CH4released from hotspot seepsmore »in older, expanding thermokarst lakes (14CCH420 079 ± 1227 years BP, mean ± standard error (s.e.m.) years) to deeper taliks (thaw bulbs) compared to younger14CCH4in new lakes (14CCH48526 ± 741 years BP) with shallower taliks. We find that smaller, non-hotspot ebullition seeps have younger14C ages (expanding lakes 7473 ± 1762 years; new lakes 4742 ± 803 years) and that their emissions span a larger historic range. These observations provide a first-order constraint on the magnitude and decadal-scale duration of CH4-hotspot seep emissions following formation of thermokarst lakes as climate warms.

    « less
  5. Abstract
    Lakes are abundant features on coastal plains of the Arctic and most are termed &#34;thermokarst&#34; because they form in ice-rich permafrost and gradually expand over time. The dynamic nature of thermokarst lakes also makes them prone to catastrophic drainage and abrupt conversion to wetlands, called drained thermokarst lake basins (DTLBs). Together, thermokarst lakes and DTLBs cover up to 80% of arctic lowland regions, making understanding their response to ongoing climate change essential for coastal plain environmental assessment. Datasets presented here document water level and temperature (surface and ground) regimes for a large (38 sites) array of lake with high drainage potential and lake basin (DTLBS), which have already drained, located on differing terrain units of Alaska&#39;s Arctic Coastal Plain. Lake data was measured along deep protected shorelines using pressure transducers to record hourly water level and bed temperature. Wetland (DTLB) data was also measured with pressure transducers and ground thermistors at 25 and 100 centimeters (cm) depth. Of special interest at some DTLB sites was the potential occurrence of snow-dam outburst events during the early summer snowmelt periods. In these cases, pressure transducers were set to log at 10 minute intervals for this period. All data archived here areMore>>