Modeling oscillations in connected glacial lakes

Abstract Mountain glaciers and ice sheets often host marginal and subglacial lakes that are hydraulically connected through subglacial drainage systems. These lakes exhibit complex dynamics that have been the subject of models for decades. Here we introduce and analyze a model for the evolution of g...

Full description

Bibliographic Details
Published in:Journal of Glaciology
Main Authors: Stubblefield, Aaron G., Creyts, Timothy T., Kingslake, Jonathan, Spiegelman, Marc
Format: Article in Journal/Newspaper
Language:English
Published: Cambridge University Press (CUP) 2019
Subjects:
Online Access:http://dx.doi.org/10.1017/jog.2019.46
https://www.cambridge.org/core/services/aop-cambridge-core/content/view/S0022143019000467
Description
Summary:Abstract Mountain glaciers and ice sheets often host marginal and subglacial lakes that are hydraulically connected through subglacial drainage systems. These lakes exhibit complex dynamics that have been the subject of models for decades. Here we introduce and analyze a model for the evolution of glacial lakes connected by subglacial channels. Subglacial channel equations are supplied with effective pressure boundary conditions that are determined by a simple lake model. While the model can describe an arbitrary number of lakes, we solve it numerically with a finite element method for the case of two connected lakes. We examine the effect of relative lake size and spacing on the oscillations. Complex oscillations in the downstream lake are driven by discharge out of the upstream lake. These include multi-peaked and anti-phase filling–draining events. Similar filling–draining cycles have been observed on the Kennicott Glacier in Alaska and at the confluence of the Whillans and Mercer ice streams in West Antarctica. We further construct a simplified ordinary differential equation model that displays the same qualitative behavior as the full, spatially-dependent model. We analyze this model using dynamical systems theory to explain the appearance of filling–draining cycles as the meltwater supply varies.