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.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.