Challenges in predicting Greenland supraglacial lake drainages at the regional scale

A leading hypothesis for the mechanism of fast supraglacial lake drainages is that transient extensional stresses briefly allow crevassing in otherwise compressional ice flow regimes. Lake water can then hydrofracture a crevasse to the base of the ice sheet, and river inputs can
maintain this connection as a moulin. If future ice sheet models are to accurately represent moulins, we must understand their formation processes, timescales, and locations. Here,
we use remote-sensing velocity products to constrain the relationship between strain rates and lake drainages across ∼ 1600 ksq.m in Pâkitsoq, western Greenland, between 2002–2019. We find significantly more extensional background strain rates at moulins associated with fast-draining lakes than at slow-draining or non-draining lake moulins. We test whether moulins in more extensional background settings drain their lakes earlier, but we find insignificant correlation. To investigate the frequency at which strain-rate transients are associated with fast lake drainage, we examined Landsat-derived strain rates over 16 and 32 d periods at moulins associated with 240 fast-lake-drainage events over 18 years. A low signal-to-noise ratio, the presence of water, and the multi-week repeat cycle obscured any resolution of the hypothesized transient strain rates. Our results support the hypothesis that transient strain rates drive fast lake drainages. However, the current generation of ice sheet velocity products, even when stacked across hundreds of fast lake drainages, cannot resolve these transients. Thus, observational progress in understanding lake drainage initiation will rely on field-based tools such as GPS networks and photogrammetry.