A Physical Model of Moulin Evolution on the Greenland Ice Sheet

Nearly all proglacial water discharge from the Greenland Ice Sheet is routed englacially via moulins. Identification of these moulins in high-resolution imagery is a frequent topic of study, but the processes controlling how and where moulins form remain poorly understood. Because moulins may reasonably compose approximately 10-15% of the englacial-subglacial hydrologic system, the evolution and shape of moulins can alter both the timing and variability of meltwater inputs to the bed. This evolution can impact both the form of the subglacial hydrologic system and associated response of ice motion. Here, we develop a physical model of moulin formation and evolution to constrain the role of englacial processes in shaping the form and structure of the subglacial hydrologic system. Within this model, moulin geometry is controlled by a balance of viscous and elastic deformation and is dependent on that deformation, refreezing, and the dissipation of turbulent and sensible heat energy. All of which are dependent on the characteristics of the available supraglacial meltwater and the surrounding ice. We find moulin geometry is responsive to changes in these parameters over the course of hours to days, indicating that diurnal and multi-day variations in melt can substantially alter the geometry of a moulin and, consequently, the pressure-discharge relationship at the bed of the ice sheet. Therefore, there is no single moulin shape that can appropriately represent englacial storage across the Greenland Ice Sheet.