954 Journal of Glaciology, Vol. 57, No. 205, 2011
Assessing the summer water budget of a moulin basin in the Sermeq Avannarleq ablation region, Greenland ice sheet
Daniel McGRATH,1 William COLGAN,1 Konrad STEFFEN,1 Phillip LAUFFENBURGER,2 James BALOG3 1Cooperative Institute for Research in Environmental Sciences, University of Colorado at Boulder, Boulder, Colorado 80309-0216, USA E-mail: [email protected] 2Aerospace Engineering Sciences Department, University of Colorado at Boulder, Boulder, Colorado, 80309-0429, USA 3Extreme Ice Survey, 1435 Yarmouth Avenue, Boulder, Colorado 80304-4338, USA
ABSTRACT. We provide an assessment of the supraglacial water budget of a moulin basin on the western margin of the Greenland ice sheet for 15 days in August 2009. Meltwater production, the dominant input term to the 1.14 0.06 km2 basin, was determined from in situ ablation measurements. The dominant water-output terms from the basin, accounting for 52% and 48% of output, respectively, were moulin discharge and drainage into crevasses. Moulin discharge exhibits large diurnal variability (0.017–0.54 m3 s–1) with a distinct late-afternoon peak at 16:45 local time. This lags peak meltwater production by 2.8 4.2 hours. An Extreme Ice Survey time-lapse photography sequence complements the observations of moulin discharge. We infer, from in situ observations of moulin geometry, previously published borehole water heights and estimates of the temporal lag between meltwater production and observed local ice surface uplift (‘jacking’), that the transfer of surface meltwater to the englacial water table via moulins is nearly instantaneous (<30 min). We employ a simple crevasse mass-balance model to demonstrate that crevasse drainage could significantly dampen the surface meltwater fluctuations reaching the englacial system in comparison to moulin discharge. Thus, unlike crevasses, moulins propagate meltwater pulses to the englacial system that are capable of overwhelming subglacial transmission capacity, resulting in enhanced basal sliding.
1. INTRODUCTION glacier, resulting in surface uplift interpreted as bed separ- The Greenland ice sheet (GrIS) is currently losing 200– ation. This process reduces basal friction and increases basal 250 Gt a–1 of ice through a combination of an increasingly sliding velocities (Iken and others, 1983; Willis and others, negative surface mass balance and enhanced ice discharge 1996; Anderson and others, 2004; Bartholomaus and others, from major outlet glaciers (Hanna and others, 2008; Thomas 2008). Similar to alpine glaciers, enhanced basal sliding in and others, 2009; Van den Broeke and others, 2009). Both in the GrIS has been shown to continue as long as meltwater situ GPS and interferometric synthetic aperture radar (InSAR) production exceeds subglacial transmission capacity, creat- satellite observations in the western ablation zone of the ing conditions of positive net water storage in the subglacial GrIS demonstrate a distinct annual ice velocity cycle, in environment (Colgan and others, 2011). During the melt which peak velocities are observed during the summer melt season, conduits enlarge by melting from the frictional season (Zwally and others, 2002; Joughin and others, 2008; heating of the flowing water, which allows the subglacial Bartholomew and others, 2010; Colgan and others, 2011). hydrologic system to evolve to accommodate larger melt- On daily and seasonal timescales, higher surface velocities water fluxes later in the melt season (Ro¨thlisberger, 1972; have been attributed to enhanced basal sliding, which Hock and Hooke, 1993). As a result, a transition occurs mid- occurs when meltwater input exceeds the transmission melt season as the subglacial transmission capacity exceeds capacity of the subglacial hydrologic network (Iken and meltwater input, and water is efficiently drained via low- others, 1983; Anderson and others, 2004; Shepherd and pressure channels (Cutler, 1998; Bartholomew and others, others, 2009; Bartholomew and others, 2010). 2010). Overlaid on this seasonal meltwater pattern is a diurnal meltwater input cycle, resulting in a daily cycle in which input exceeds transmission capacity (Schoof, 2010). 1.1. Glacier hydrology This cycle drives localized uplift and enhanced basal sliding Strong diurnal and seasonal variations in meltwater produc- in the GrIS several hours after peak surface meltwater pro- tion imply that the englacial and subglacial hydrologic duction in the late afternoon (Shepherd and others, 2009). networks are seldom in steady state, but rather constantly This is followed by a surface lowering and subsequent adjusting to changing input volumes (Hock and Hooke, decrease in ice velocity, presumably as the meltwater input 1993; Cutler, 1998; Bartholomaus and others, 2008). At the volume falls below the efficiency of the subglacial hydro- beginning of the melt season, meltwater is delivered to a logic network (Shepherd and others, 2009). quiescent subglacial hydrologic network, which has largely The cumulative effect of enhanced diurnal and seasonal closed during the winter through creep closure (Nye, 1953). velocities is an increase in total annual ice displacement, the As a result of initial meltwater input exceeding subglacial magnitude of which is positively correlated with modeled transmission capacity, subglacial water pressure increases, meltwater production (Zwally and others, 2002). The sea- driving water into a distributed cavity network beneath the sonal acceleration is most significant ( 50% increase above McGrath and others: Summer water budget of a moulin basin 955 mean annual ice velocity) in land-terminating regions of the thus meltwater drainage per crevasse is substantially smaller ice sheet, and less significant ( 10–15% increase) in outlet than for moulins, which typically drain a well-developed glaciers, where velocity is more directly related to changes catchment basin. Under-appreciated components of the in back-stress at the tidewater terminus (Howat and others, supraglacial hydrologic system are the widespread smaller 2005; Joughin and others, 2008). Observations show that ice (5–90 cm) surface and englacial cracks, which we refer to as velocities in the ablation zone respond quickly to short-term ‘fractures’ to differentiate them from typical crevasses. These changes in meltwater production, although significant fractures are observed to persist to depths of 70 m and uncertainty exists in predicting how the GrIS hydrology typically have near-vertical orientations, with surface ex- system evolves over decadal and longer timescales to pressions that are not preferentially aligned with flow changing volumes of meltwater. One study finds no long- direction. Borehole observations in temperate glaciers term (decadal) increase in mean annual ice velocity despite suggest that surface fractures penetrate to significant depth significant seasonal melt-driven accelerations occurring ( 130 m), consisting of near-vertical ( 708) features 0.3– during a 17 year period (Van de Wal and others, 2008). 20 cm wide (Fountain and others, 2005). Water flow has The apparently conflicting nature of these observations been observed to be slow (1–2 cm s–1) in these fractures, but highlights the uncertainty in this aspect of ice-sheet they frequently intersect other fractures, suggesting glaciers evolution and certainly warrants long-term observations in may be analogous to ‘fractured rock-type aquifers’ with high order to develop a unifying explanation. hydraulic conductivity (Fountain and others, 2005; Colgan The supraglacial hydrologic cycle is most pronounced and others, 2011). Fountain and Walder (1998) suggest that along the marginal zone of the GrIS where the relatively in temperate glaciers the majority of water storage occurs in high surface slope (2–58), complex surface topography and the englacial, rather than subglacial, hydrologic system. The relatively high ablation rates differ substantially from the vast englacial system likely consists of a combination of well- interior accumulation zone. Melting of snow and glacier ice connected and discrete voids and fractures, which together in the ablation zone produces water that flows across the ice create a bulk macroporosity of 0.1–10% (Pohjola, 1994; surface, developing supraglacial streams, which can drain Harper and Humphrey, 1995; Huss and others, 2007). into surface depressions to form supraglacial lakes or drain Observations of bubble-free blue-ice bands, both at the directly into moulins (Box and Ski, 2007). Meltwater can surface (exposed as the surface ablates) and at depth with also drain into crevasses or smaller fractures, where it may borehole cameras, suggest that meltwater stored in the seasonally refreeze (Catania and others, 2008). Supraglacial englacial system frequently refreezes (Pohjola, 1994; Harper streams evolve from interconnected runnels into an arbor- and Humphrey, 1995). escent network on the ice surface, incising through thermal erosion at a rate that exceeds the surface ablation rate (Knighton, 1981; Marston, 1983). On account of the latent 2. FIELD SITE energy contained in liquid water, even modest water The site of this investigation is a supraglacial catchment basin temperatures (0.005–0.018C) can result in channel incision in the Sermeq Avannarleq ablation zone on the western flank rates of 3.8–5.8 cm d–1 (Pinchak, 1972; Marston, 1983). of the GrIS, 6 km from the ice margin (Fig. 1). Mean surface Meltwater lakes are a common feature of the GrIS ablation rate is 1.65 m w.e. a–1 (Fausto and others, 2009). The ablation zone. They typically range in size from a few basin is 360 m below the regional equilibrium-line altitude hundred meters to >2 km in diameter, with mean water (ELA). The regional ELA is 1140 m a.s.l., although inter- depths of 2–5 m (Box and Ski, 2007). The presence of annual variations over the past decade range from 1000 to meltwater on the ice surface, and the subsequent increase in 1400 m a.s.l. (Mernild and others, 2010). The ice in this water depth as melt lakes develop, reduces surface albedo, region is believed to be temperate, based on nearby (<10 km) and increases shortwave radiation absorption, amplifying borehole temperatures (Thomsen and Thorning, 1992). The melt rates (Perovich and others, 2002). Further, meltwater supraglacial hydrologic network of this basin is arborescent, lakes are observed to frequently drain over short timescales with multiple tributaries collecting into a single large stream (hours to days) at discharge rates exceeding 300 m3 s–1 (Box that discharges into a moulin at 69.5548 N, 49.8998 W and and Ski, 2007; Das and others, 2008). 776 m elevation. Ice thickness at this site is 530 m. Field Conversely, other supraglacial basins have established observations at our study site in August 2007, 2008 and 2009 stream networks that terminate in a moulin, which indicate that the main trunk stream has maintained an commonly consists of an initial vertical shaft and subsequent approximate geometry of 1–4 m width and 1–6 m depth, plunge pools (Holmlund, 1988; Gulley, 2009). Moulins are although water has never been observed to fill the entire presumed to form from hydrofracturing of water-filled stream geometry. Panchromatic WorldView-1 satellite im- crevasses (Boon and Sharp, 2003; Alley and others, 2005; agery (50 cm pixel resolution) of the study site was acquired Van der Veen, 2007; Das and others, 2008) and are believed on 15 July 2009. Manual identification of connecting to persist for multiple years in locations fixed by bedrock supraglacial tributaries was used to delineate the extent of geometry (Catania and Neumann, 2010). Unlike crevasse the basin. We estimate a basin area, A, of 1.14 0.06 km2. drainage, moulins concentrate the surface meltwater pro- Comparison of this recent satellite imagery with a historical duced over a large area and deliver it to the englacial system topographic map (Thomsen, 1986; Thomsen and others, at a single point. 1988) suggests that since 1985 the supraglacial tributaries Surface crevasses are ubiquitous features of the ablation have occupied the same general positions and the moulin zone that form when the longitudinal strain rate exceeds the location has changed by <250 m. From this we speculate that critical fracture toughness of the ice or in response to thermal this supraglacial channel network is a recurring annual stress in the spring/early summer (Sanderson, 1978). These feature of the local hydrologic system. The moulin is located features, which are typically linear in nature and transect in the basin bottom (at the lowest elevation), although a slopes, prevent large catchment areas from developing, and series of abandoned moulins exist as a ‘string of pearls’, 956 McGrath and others: Summer water budget of a moulin basin
Fig. 1. Supraglacial stream network in West Greenland (inset) overlaid on panchromatic WorldView-1satellite imagery (acquired 15 July 2009) with elevation shading from the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) global digital elevation model. The locations of the stream gauging station (G), moulin (M) and automatic weather station (A) are identified. advected towards the ice margin at an annual ice velocity of fractures similar to those observed on alpine glaciers 100ma–1. These abandoned moulins are not regularly (Fountain and others, 2005). spaced and are thus not necessarily annual features. Field observations in 2008 found englacial conduits transporting water into an abandoned moulin shaft at 70 m depth, 3. METHODS suggesting while surface entry points may change, the 3.1. Data acquisition englacial drainage system is relatively persistent (Catania A stream gauging station was deployed on the main and Neumann, 2010). supraglacial stream 30 m upstream of the moulin between Large crevasses ( 1 m wide) are observed in the World- 3 and 17 August 2009 (day of year (DOY) 215–229; Fig. 2a). View-1 imagery at the periphery of the basin, but are not It recorded water surface height with a Campbell Scientific present in the basin bottom. Refrozen fractures of bubble- SR-50 sonic sensor and stream velocity with a Geopacks free ice (5–90 cm wide) are observed at the surface across MFP51 Flowmeter. The sonic instrument measured instant- the basin bottom and at depth ( 70 m) within a large, aneous water surface height at 15 min intervals, while abandoned moulin in the basin (Fig. 2a and b), suggesting stream velocity was measured as the mean over 15 min the widespread existence of refrozen englacial voids and intervals. The sonic sensor was installed 1 m above the water surface, while the flowmeter was installed at the center width of the stream channel at a fixed height. While the absolute position of the impeller was fixed, its relative position fluctuated within the water column as both the stage height varied and the channel incised downwards (Fig. 2a). An automatic weather station (AWS) deployed in 2007 was located 160 m from the gauging station and measured ice surface height (Campbell Scientific SR-50), air temperature (Vaisala HMP50) and wind speed and direction (RM Young 050103). An Extreme Ice Survey digital camera (Nikon D200, 20 mm lens) was installed 15 m downstream of the gauging station and recorded photographs every 15 min (Animation 1). These in situ measurements allow us to assess the water budget of this supraglacial stream basin. 3.2. Water mass-balance model Using conservation of water mass, we may formulate our Animation 1. Animation of Extreme Ice Survey time-lapse photo- stream basin-scale water budget as a simple balance of graphs (every 15 min) of supraglacial stream and gauging station. inputs and outputs: Full animation available at www.igsoc.org/hyperlink/10J209_animation1.mp4. 0 ¼ ðÞ IMELT þ IRAIN ðÞþQMOULIN þ QCREVASSE þ QE ÁS, ð1Þ McGrath and others: Summer water budget of a moulin basin 957
Fig. 2. (a) Oblique aerial photo of the supraglacial stream and gauging station ( 30 m upstream from moulin). Note numerous refrozen fractures crossing stream channel. The solar panel is 42 cm 27 cm and water flow is from right to left in photo. (b) Englacial fractures filled with dark-blue bubble-free refrozen meltwater on the side-wall of a moulin shaft at 40 m depth.
–3 where IMELT and IRAIN are water inputs to the basin due to water (1000 kg m ), respectively. Although the AWS was not ablation and rainfall, QMOULIN, QCREVASSE and QE are water equipped to record liquid precipitation, we estimate water fluxes from the basin due to moulin discharge, crevasse input due to rainfall from the in situ Extreme Ice Survey drainage and evaporation, and ÁS is the rate of change in photographic record. The photographic record indicates that surface storage. From the absence of supraglacial meltwater only two brief rainfall events occurred over the 15 day study ponds within the study area we assume there is no period. The first event, 1.25 0.25 hours in duration, significant change in surface water storage over the study occurred on 3 August 2009 while the second event, period (i.e. ÁS 0). Field data may be used to constrain all 4.25 0.25 hours in duration, occurred on 11 August water-budget terms except QCREVASSE, which includes drain- 2009. Both these rain events fit the description of being age via both crevasses and smaller moulins (<0.25 m ‘light’ in intensity (<2.5 mm h–1; American Meteorological diameter) within the basin. We treat this as a free term Society, http://amsglossary.allenpress.com/glossary/search? when solving the water budget, so QCREVASSE incorporates id=rain1, so we estimate the intensity, R, of the events as the uncertainty in all other budget terms. We reformulate the 0.5 0.25 and 1.0 0.5 mm h–1 respectively. Even assuming basin-scale water budget as extreme intensities (10 mm h–1), the brevity of these events precludes rainfall from being a major term in the cumulative Q ¼ ðÞ I þ I ðÞQ þ Q : ð2Þ CREVASSE MELT RAIN MOULIN E water budget. Combining rainfall intensity and basin area Meltwater production due to ablation was calculated allows the rate of water input due to rainfall to be estimated Á Á from observed ice surface height change, ZS/ t, during the according to study period (Fig. 3). Ice surface height measurements were ¼ : ð Þ made every hour. No significant snow or ice deposition IRAIN RA 4 occurred during the study period (Animation 1). Surface Water output due to moulin discharge was determined height change at the AWS location was assumed to be from gauge station data, consisting of stage height and water representative of ablation across the basin due to minimal velocity measurements (Fig. 4a and b). Unlike terrestrial variations in slope, aspect and surface albedo (by August, all river gauging where the bed may be assumed to be constant snow has melted and the surface is bare ice). Thus, we over short timescales, in supraglacial stream gauging the bed calculate the hourly meltwater production as incises into the ice at an appreciable rate. Thus, both the water surface and bed elevations vary relative to the fixed ÁZS i I ¼ A , ð3Þ impeller elevation through time. The cross-sectional profile MELT Át w of the supraglacial channel was measured upon installation –3 where i and w represent the density of ice (900 kg m ) and (day 215). The time-lapse photographic record (Animation 1) 958 McGrath and others: Summer water budget of a moulin basin
Fig. 3. Observed ice surface height, ZS (black line), and corresponding instantaneous basin-wide meltwater production from in situ surface height measurements (cf. Equation (3); blue line) over the 15 day study period. Open circles are values interpolated by a high-order polynomial fit during times of instrument error.
and measurements upon gauging-station removal suggest and changing water surface height, we develop a site- that the assumption of a constant geometry that incised specific rating curve using velocity values when the impeller downward is justified. The stream incised into the ice was located at 0.37 0.03% of water column height. This is surface at a rate of 3.3 0.47 cm d–1. This was calculated as consistent with the assumption that the depth-averaged the slope of a linear best fit through daily minimum stream velocity in turbulent flow occurs at 37% of water column surface elevation values over the 15 day study period height (Reynolds number 2000; Anderson and Anderson, 2 (Fig. 4a). By constraining the stream bed elevation and 2010). We apply the rating curve (QMOULIN = 1.234H – recording the water surface elevation and cross-sectional 0.303H; n = 88; r2 = 0.99; Fig. 5b) to calculate moulin profile (Fig. 5a), the cross-sectional area of the stream, discharge at times when the impeller is outside this range ACROSS, can be determined at any time-step via numerical of water column heights and after the impeller failed on day integration. 220. Depth-averaged velocity values, u, were multiplied by A theoretical rate of water output due to evaporation the cross-sectional stream area (Fig. 4c) to determine moulin was calculated based on latent heat flux: discharge: ¼ E : ð Þ ¼ : ð Þ QE 6 QMOULIN uACROSS 5 LS w –2 However, due to the fixed elevation of the impeller in space We assume a latent heat flux, E, of 6.8 W m , which was
Fig. 4. (a) Water surface height over the 15 day period relative to the impeller elevation (0 m). Black line is stream bottom, which incised at 3.3 0.47 cm d–1. Incision rate is taken as the slope of the minimum stream surface height over the 15 day study period (red lines). (b) Observed stream velocity, u, at all times (cyan) and when located at 0.37 0.03% of the water column height, H, and used to develop the rating curve (blue). Impeller failed on day 220. (c) Calculated instantaneous discharge, QMOULIN (bold line), error (thin line) of the supraglacial stream as it enters the moulin. McGrath and others: Summer water budget of a moulin basin 959
Fig. 5. (a) Cross-sectional area of the supraglacial stream at the gauging station (measured DOY 215). (b) Rating curve error used to 2 2 calculate stream discharge (QMOULIN = 1.234H – 0.303H; r = 0.99). Open circles are measured values used to construct rating curve (n = 88).
observed in August 2000 at the nearby JAR2 AWS ( 17 km at the elevation of the moulin ( 1.65 m w.e.; Fausto and south-southwest at 570 m elevation; Box and Steffen, others, 2009) over the duration of the melt season (Colgan 2001). The latent heat of sublimation, LS, is taken as and others, 2011). Within each day, ablation is distributed 2.834 106 Jkg–1. over a 12 hour period with a secondary sine function, while A comprehensive error analysis for each term and the the remaining 12 hours of the day experience no ablation. complete water budget is included in the Appendix. The model is solved at a 1 hour time-step, to resolve diurnal meltwater drainage into the idealized crevasse, using 3.3. Crevasse mass-balance model MATLAB’s semi-implicit ‘stiff’ solver (ode15s). We developed a simple one-dimensional model to describe We do not impose a channel or conduit geometry to the mass balance of both water and ice within an idealized calculate water discharge from the crevasse to the englacial crevasse over multiple melt seasons. This model is meant to hydrology system. Instead, we employ a linear reservoir serve as a proof of concept that crevasses can significantly model (cf. Flowers and Clarke, 2002) and simply assign a dampen the diurnal cycle of meltwater entering the mean residence time for water within the crevasse network, englacial system in comparison to moulins. However, this and thus neglect a possible influence of the englacial simplified model is neither meant to close the water budget drainage system on crevasse discharge. Modeled crevasse nor explicitly determine the percentage of refrozen versus pseudo-discharge is therefore equal to the mass of water discharged meltwater. Further, we do not examine melt- within the crevasse, Mw, divided by this mean residence water-driven crevasse hydrofracture (Van der Veen, 2007; Das and others, 2008). Our idealized crevasse has a simple triangular geometry with a depth of 10 m and a width of 0.5 m (Animation 2). The mass of liquid-phase water within a crevasse, Mw, may be described by three terms: (1) the rate at which surface meltwater enters the crevasse, i w, (2) the rate at which water drains from the crevasse, Mw/ res, and (3) the rate at which liquid water freezes into solid ice within the crevasse, q /Lf: dMw ¼ Mw q ð Þ i w , 7 dt res Lf 3 –1 where i is the rate of external meltwater input (m h ), res is the mean residence time of water within the crevasse network (hours), q is the heat flux into the ice per unit area –1 (J h ; Equation (9)) and Lf is the latent heat of fusion (333 550 J kg–1). The rate of external meltwater input can be calculated according to
i ¼ asdcdy, ð8Þ where a is the modeled surface ablation rate (m w.e. h–1), d s c Animation 2. Top: Masses of ice, M , and water, M , within the is the mean crevasse spacing (taken as 25 m; Phillips and w i idealized crevasse through time. Bottom: Crevasse water mass- others, 2010) and dy is a unit width (1 m). We assume that all balance terms: surface meltwater input, i w, crevasse discharge, meltwater produced on both sides of the crevasse within a Mw/ res, and refreezing into solid ice, q /L. Right: Schematic distance of dc/2 enters the crevasse. We model daily ablation showing transient refrozen ice (cyan) and water (blue) levels within rate by approximating the annual ablation cycle with a sine the crevasse. Full animation available at function that integrates to the observed annual ablation value www.igsoc.org/hyperlink/10J209_animation2.mov. 960 McGrath and others: Summer water budget of a moulin basin
Fig. 6. Time series of daily water budget components for the moulin basin over the 15 day study period. The total input (blue) is the sum of IMELT and IRAIN (cyan) while the total loss (red) is the sum of evaporation (green), QCREVASSE (black) and QMOULIN.