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.