Articles | Volume 16, issue 4
Research article
08 Apr 2022
Research article |  | 08 Apr 2022

A new Stefan equation to characterize the evolution of thermokarst lake and talik geometry

Noriaki Ohara, Benjamin M. Jones, Andrew D. Parsekian, Kenneth M. Hinkel, Katsu Yamatani, Mikhail Kanevskiy, Rodrigo C. Rangel, Amy L. Breen, and Helena Bergstedt

Data sets

Arctic lake transient electromagnetic (TEM) soundings 2016–2017 A. D. Parsekian, A. L. Creighton, B. M. Jones, A. Bondurant, and C. D. Arp

Floating and bedfast lake ice regimes across Arctic 25 Alaska using space-borne SAR imagery from 1992–2016 M. Engram

Short summary
New variational principle suggests that a semi-ellipsoid talik shape (3D Stefan equation) is optimum for incoming energy. However, the lake bathymetry tends to be less ellipsoidal due to the ice-rich layers near the surface. Wind wave erosion is likely responsible for the elongation of lakes, while thaw subsidence slows the wave effect and stabilizes the thermokarst lakes. The derived 3D Stefan equation was compared to the field-observed talik thickness data using geophysical methods.