skip to main content


Title: A new Stefan equation to characterize the evolution of thermokarst lake and talik geometry
Abstract. Thermokarst lake dynamics, which play an essential role in carbon releasedue to permafrost thaw, are affected by various geomorphological processes.In this study, we derive a three-dimensional (3D) Stefan equation tocharacterize talik geometry under a hypothetical thermokarst lake in thecontinuous permafrost region. Using the Euler equation in the calculus ofvariations, the lower bounds of the talik were determined as an extremum ofthe functional describing the phase boundary area with a fixed total talikvolume. We demonstrate that the semi-ellipsoid geometry of the talik isoptimal for minimizing the total permafrost thaw under the lake for a givenannual heat supply. The model predicting ellipsoidal talik geometry wascompared to talik thickness observations using transient electromagnetic(TEM) soundings in Peatball Lake on the Arctic Coastal Plain (ACP) ofnorthern Alaska. The depth : width ratio of the elliptical sub-lake talik cancharacterize the energy flux anisotropy in the permafrost, although the lakebathymetry cross section may not be elliptic due to the presence ofnear-surface ice-rich permafrost. This theory suggests that talikdevelopment deepens lakes and results in more uniform horizontal lakeexpansion around the perimeter of the lakes, while wind-induced waves andcurrents are likely responsible for the elongation and orientation ofshallow thermokarst lakes without taliks in certain regions such as the ACPof northern Alaska.  more » « less
Award ID(s):
1806213 1806202 1806287 1820883
NSF-PAR ID:
10321230
Author(s) / Creator(s):
; ; ; ; ; ; ; ;
Date Published:
Journal Name:
The Cryosphere
Volume:
16
Issue:
4
ISSN:
1994-0424
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Methane (CH4) release to the atmosphere from thawing permafrost contributes significantly to global CH4emissions. However, constraining the effects of thaw that control the production and emission of CH4is needed to anticipate future Arctic emissions. Here are presented robust rate measurements of CH4production and cycling in a region of rapidly degrading permafrost. Big Trail Lake, located in central Alaska, is a young, actively expanding thermokarst lake. The lake was investigated by taking two 1 m cores of sediment from different regions. Two independent methods of measuring microbial CH4 production, long term (CH4accumulation) and short term (14C tracer), produced similar average rates of 11 ± 3.5 and 9 ± 3.6 nmol cm−3 d−1, respectively. The rates had small variations between the different lithological units, indicating homogeneous CH4production despite heterogeneous lithology in the surface ~1 m of sediment. To estimate the total CH4production, the CH4production rates were multiplied through the 10–15 m deep talik (thaw bulb). This estimate suggests that CH4 production is higher than emission by a maximum factor of ~2, which is less than previous estimates. Stable and radioactive carbon isotope measurements showed that 50% of dissolved CH4in the first meter was produced further below. Interestingly, labeled14C incubations with 2‐14C acetate and14C CO2indicate that variations in the pathway used by microbes to produce CH4depends on the age and type of organic matter in the sediment, but did not appear to influence the rates at which CH4 was produced. This study demonstrates that at least half of the CH4produced by microbial breakdown of organic matter in actively expanding thermokarst is emitted to the atmosphere, and that the majority of this CH4is produced in the deep sediment.

     
    more » « less
  2. 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 seeps 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.

     
    more » « less
  3. Abstract

    The impact of permafrost thaw on hydrologic, thermal, and biotic processes remains uncertain, in part due to limitations in subsurface measurement capabilities. To better understand subsurface processes in thermokarst environments, we collocated geophysical and biogeochemical instruments along a thaw gradient between forested permafrost and collapse‐scar bogs at the Alaska Peatland Experiment site near Fairbanks, Alaska. Ambient seismic noise monitoring provided continuous high‐temporal resolution measurements of water and ice saturation changes. Maps of seismic velocity change identified areas of large summertime velocity reductions nearest the youngest bog, indicating potential thaw and expansion at the bog margin. These results corresponded well with complementary borehole nuclear magnetic resonance measurements of unfrozen water content with depth, which showed permafrost soils nearest the bog edges contained the largest amount of unfrozen water along the study transect, up to 25% by volume. In situ measurements of methane within permafrost soils revealed high concentrations at these bog‐edge locations, up to 30% soil gas. Supra‐permafrost talik zones were observed at the bog margins, indicating talik formation and perennial liquid water may drive lateral bog expansion and enhanced permafrost carbon losses preceding thaw. Comparison of seismic monitoring with wintertime surface carbon dioxide fluxes revealed differential responses depending on time and proximity to the bogs, capturing the controlling influence of subsurface water and ice on microbial activity and surficial emissions. This study demonstrates a multidisciplinary approach for gaining new understanding of how subsurface physical properties influence greenhouse gas production, emissions, and thermokarst development.

     
    more » « less
  4. Abstract

    The ongoing global temperature rise enhances permafrost thaw in the Arctic, allowing Pleistocene‐aged frozen soil organic matter to become available for microbial degradation and production of greenhouse gases, particularly methane. Here, we examined the extent and mechanism of anaerobic oxidation of methane (AOM) in the sediments of four interior Alaska thermokarst lakes, which formed and continue to expand as a result of ice‐rich permafrost thaw. In cores of surface (~ 1 m) lake sediments we quantified methane production (methanogenesis) and AOM rates using anaerobic incubation experiments in low (4°C) and high (16°C) temperatures. Methanogenesis rates were measured by the accumulation of methane over ~ 90 d, whereas AOM rates were measured by adding labeled‐13CH4and measuring the produced dissolved inorganic13C. Our results demonstrate that while methanogenesis was vigorous in these anoxic sediments, AOM was lower by two orders of magnitude. In almost all sediment depths and temperatures, AOM rates remained less than 2% of the methanogenesis rates. Experimental evidence indicates that the AOM is strongly related to methanogens, as the addition of a methanogens' inhibitor prevented AOM. Variety of electron acceptor additions did not stimulate AOM, and methanotrophs were scarcely detected. These observations suggest that the AOM signals in the incubation experiments might be a result of enzymatic reversibility (“back‐flux”) during CH4production, rather than thermodynamically favorable AOM. Regardless of the mechanism, the quantitative results show that near surface (< 1‐m) thermokarst sediments in interior Alaska have little to no buffer mechanisms capable of attenuating methane production in a warming climate.

     
    more » « less
  5. Abstract

    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 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.

     
    more » « less