Hostname: page-component-7c8c6479df-hgkh8 Total loading time: 0 Render date: 2024-03-29T02:28:55.730Z Has data issue: false hasContentIssue false

Rapid glacier sliding, reverse ice motion and subglacial water pressure during an autumn rainstorm

Published online by Cambridge University Press:  14 September 2017

T.J. Fudge
Affiliation:
Department of Earth and Space Sciences, University of Washington, Box 351310, Seattle, WA 98195-1310, USAE-mail: tjfudge@u.washington.edu
J.T. Harper
Affiliation:
Department of Geosciences, University of Montana, Missoula, MT 59812-1296, USA
N.F. Humphrey
Affiliation:
Department of Geology and Geophysics, University of Wyoming, Laramie, WY 82071-3006, USA
W.T. Pfeffer
Affiliation:
Institute of Arctic and Alpine Research, University of Colorado, Boulder, CO 80309-0450, USA
Rights & Permissions [Opens in a new window]

Abstract

Measurements of basal water pressure from 15 boreholes located at both local (tens of meters) and regional (kilometers) length scales were used to elucidate the pressure/sliding relationship during an autumn rapid motion event on Bench Glacier, Alaska, USA. The 8 day event had two distinct phases, each with a ten-fold speed-up with respect to winter velocity. The water pressure in all 15 boreholes varied synchronously during the speed-up. The first phase of rapid sliding began after a peak in basal water pressure and continued while the pressure was elevated and stable, or decreasing. The second phase of rapid sliding occurred when the basal water pressure was low but increasing, and terminated before the pressure peaked. Pressure and velocity do not appear unrelated, but the pressure/sliding relationship was not consistently linked to increasing, decreasing or a critical water pressure. The pressure variations and sliding accelerations are a response to a warm rainstorm, although equally large input events occurred in weeks prior with no apparent response. Drainage system evolution therefore appears to play a key role in both the acceleration and the pressure/velocity relationship. Basal cavity dynamics are likely responsible for three episodes of reverse (up-valley) motion observed after enhanced sliding.

Type
Research Article
Copyright
Copyright © The Author(s) [year] 2009

Introduction

Enhanced sliding due to seasonal water input is a long-established component of glacier motion (e.g. Reference WillisWillis, 1995) and has recently been documented in ice sheets (Reference Zwally, Abdalati, Herring, Larson, Saba and SteffenZwally and others, 2002; Reference Joughin, Das, King, Smith, Howat and MoonJoughin and others, 2008). However, the physical processes that govern the role of water in sliding remain poorly understood (Reference Fountain and WalderFountain and Walder, 1998; Reference HookeHooke, 1998). Despite the key role of water in basal processes, rapid sliding events do not necessarily coincide with times of maximum water input (Reference Truffer, Harrison and MarchTruffer and others, 2005). Instead, rapid sliding events often occur after large increases in the rate of water input, such as in the spring (e.g. Reference Mair, Nienow, Willis and SharpMair and others, 2001; Reference Harper, Humphrey, Pfeffer and LazarHarper and others, 2007), during rain events (e.g. Reference O’Neel, Echelmeyer and MotykaO’Neel and others, 2001) and during outburst floods (e.g. Reference Magnùsson, Rott, Björnsson and PálssonMagnùsson and others, 2007). The ‘antecedent conditions’ of the basal drainage therefore appear crucial to predicting sliding accelerations.

Sliding models often treat basal water pressure as the key water parameter dictating sliding (Reference LliboutryLliboutry, 1968; Reference RöthlisbergerRöthlisberger; 1972, Reference FowlerFowler, 1987; Reference SchoofSchoof, 2005). Indeed, some observational experiments have determined that basal water pressure correlates well with sliding (Reference JanssonJansson, 1995); other studies, however, have found poor correlations (e.g. Reference Iken and TrufferIken and Truffer, 1997; Reference Hanson, LeB Hooke and GraceHanson and others, 1998), and in some cases sliding has been related to water storage (Reference Kamb, Engelhardt, Fahnestock, Humphrey, Meier and StoneKamb and others, 1994; Reference MeierMeier and others, 1994). Discerning the relationship between sliding speed and water pressure is hampered by the difficulty of measuring basal pressure. Pressure is typically measured in a small number of spatially limited boreholes, or is sometimes inferred from records of surface uplift by assuming uplift results from a pressure increase (Reference Zwally, Abdalati, Herring, Larson, Saba and SteffenZwally and others, 2002).

Here we present detailed records of water pressure during a week-long autumn rapid-motion event on Bench Glacier, Alaska, USA. The subglacial water-pressure field is calculated from an array of 15 boreholes and considers both the local (tens of meters) and regional (kilometers) length scales. We find no consistent relationship between glacier motion and basal pressure at either length scale, suggesting that more complex basal processes govern the sliding accelerations.

Data Collection

Bench Glacier

Field observations were made at Bench Glacier, a temperate glacier located in the Chugach mountain range of south– central Alaska. The 7 km glacier terminates at 1100m elevation, and has an equilibrium-line altitude of about 1400m (Fig. 1). Bench Glacier has a relatively simple geometry with no major tributaries. The average surface slope is about 10˚, with one steeper section just above the equilibrium line sloping 20–30˚. Radar profiles of the bed show a parabolic cross-section without major riegels and that the bed has an average slope in the ablation area of 3.5˚ (Reference Bradford and HarperBradford and Harper, 2005). The ice depth increases up-glacier throughout the ablation area to a maximum of about 200 m near the equilibrium line. Borehole video inspection, penetrometer tests and careful monitoring of the drill tip as it encountered the bed showed only rare patches of till, so we classify the bed as ‘hard’. One major stream exits the glacier at the terminus, with average streamflow during the summer of about 5m3 s–1 (Reference AndersonAnderson and others, 2004) and at least an order of magnitude less in the winter.

Fig. 1. Map of Bench Glacier and measurement locations. Circles show boreholes; the star shows the GPS ‘rover’ receiver; the triangle indicates the position of the meteorological station and GPS ‘base’ receiver; the solid line indicates the approximate equilibrium line; dashed lines indicate elevation contours of 100 ft (∽30.5 m). Borehole sites are named by the distance, in meters, from the terminus. Regional water pressure is averaged from measurements at the five sites; local water pressure is averaged from measurements in the nine boreholes at site 2880.

Borehole water levels

A total of 43 boreholes were drilled to the bed (140–190 m) and instrumented during spring of 2002 and 2003. The boreholes were drilled using hot-water methods during May and early June when the glacier was snow-covered. We assume the boreholes to be nearly vertical based on detailed inclinometry measurements of boreholes drilled using the same equipment and procedures (Reference Harper, Humphrey, Pfeffer, Huzurbazar, Bahr and WelchHarper and others, 2001). On intersection with the bed, the drill was reversed and readvanced repeatedly in an effort to penetrate possible englacial debris near the bed.

Boreholes were spaced as pairs 20m apart at different sites spanning the entire length of the glacier, and one 16-borehole grid in the center of the ablation area. Sites are named as the distance (in meters) up-glacier from the terminus. In the autumn of 2003, pressure sensors in 15 boreholes remained operational, four from 2002 and 11 from 2003 (Fig. 1). The boreholes from 2003 were open at the surface, whereas the tops of the boreholes from 2002 had been frozen shut during the previous winter. No sensors were operational in boreholes located above site 2880.

A pressure transducer was lowered to within 1 m of the bed in each borehole and connected to a data logger at the surface. The resolution of water-level differences between boreholes is ∽2m; this results mostly from errors derived from surface lowering by ablation, although we fixed cables at the surface to minimize this effect. The data loggers recorded measurements every 15 min during the autumn season.

Ice motion

Glacier motion was recorded with two Trimble Geomatic R5700 global positioning system (GPS) receivers logging at 30 s intervals and installed as a ‘rover’ and ‘base’. The rover was installed in the center of site 2880 with piers drilled into the ice, and the base station was placed on rock on a ridge adjacent to the glacier. The receivers were operational throughout the summer and were reset between days 246 and 248. The base station malfunctioned on day 250, whereas the rover operated until day 328. Two processing techniques (described below) were used: (1) Trimble Geomatics Office for the full record (Fig. 2); and (2) Precise Point Positioning (PPP) for more precise locations during the autumn speed-up (Figs 3 and 4).

Fig. 2. Data spanning the summer to winter transition. (a) Horizontal surface motion recorded by GPS at site 2880; (b) vertical surface motion recorded by GPS at site 2880; (c) local (thick curve) and regional (thin curve) average water level (see Fig. 1 caption); (d) temperature and rainfall. GPS positions processed with Trimble Geomatics Office; horizontal (∽10 cm) and vertical (∽20 cm) errors are not shown. The speed-up event (days 270–278) clearly stands out in all records. Prior to the speed-up (days 248–270), the horizontal velocity slowed and the surface elevation had downward motion. High temperatures and significant rain earlier in the year (days 217–227) failed to produce a similar speed-up. The average water levels are only considered representative of the average basal pressure from days 270–276 (see Fig. 5).

Fig. 3. Autumn speed-up event. (a) Horizontal surface motion recorded by GPS at site 2880; error bars represent 95% confidence interval; (b) local (solid curve) and regional (dashed curve) water pressure; horizontal dotted lines are pre-event pressure; (c) temperature and rainfall. GPS positions processed using PPP. The speed-up is divided into two phases. During phase I, rapid sliding occurred for ∽36 hours (first shaded box) during elevated pressure and after a pressure spike. During phase II, rapid sliding occurred for ∽16 hours (second shaded box) while the pressure was increasing. Velocity slowed during phase II before pressure reached a maximum.

Fig. 4. Reverse-motion events. Horizontal motion (a) and vertical motion (b) recorded by GPS at site 2880 processed with PPP. Error bars represent 95% confidence interval (CI). Shading shows times of reverse motion. The amount of reverse motion was 15 cm (4–26 cm, when considering 95% CI) at the end of phase I, 21 cm (8–34 cm) during phase II, and 21 cm (11–31 cm) at the end of phase II. Vertical motion cannot be determined for individual episodes because of the large errors. The total downward surface motion from days 272–277 was 62 cm (26–98 cm).

The failure of the local base station required the use of two more distant but continuously operating base stations, Potato Point (55km east) and Cape Hinchinbrook (103 km south-southeast), for processing with Trimble Geomatics Office. The GPS records were processed into 12 hour positions for the full length of the record. Errors were large because of the long baselines between the rover at site 2880 and the base stations; the horizontal and vertical errors were typically 10 and 20 cm, respectively.

The PPP system available through the Canadian Natural Resources Department was used to achieve more precise positions for data covering the autumn speed-up (Reference Heroux and KoubaHeroux and Kouba, 2001), days 269–278. It is these more precise data that will be compared in detail with the water-level records during the speed-up. PPP processing uses corrected satellite orbits and does not require a fixed reference station. Static 4 hour positions were calculated, as this was the minimum time resolution with acceptable positional errors. The mean, two-sigma errors for latitude, longitude and elevation are 3.9 cm, 7.3 cm and 14.5 cm, respectively. To check for biasing errors that might be introduced by processing a moving glacier in static mode (Reference KingKing, 2004), we also processed using kinematic methods with the PPP software and did not find any consistent biasing in the static positions.

Meteorological data

A meteorological station was installed on a ridge adjacent to Bench Glacier (Fig. 1); it recorded temperature, rainfall, solar radiation, wind speed and wind direction in 4 hour averages from May to December, 2003. The ridge site is higher than the immediately adjacent glacier surface but roughly equal in elevation to site 2880. The main difference between the ridge and the glacier is the lack of katabatic winds at the ridge during warm, clear-sky days. A pressure transducer was placed in the outlet stream at the end of summer but was not recovered the following spring; therefore we do not have streamflow estimates during the autumn.

Results

Records from the period mid-July to mid-November reveal a distinct seasonal evolution of borehole water-level patterns and glacier motion in addition to several significant meteorological events (Fig. 2). The week-long autumn speed-up (days 270–278) clearly stands out within the records as an unusual event. The horizontal motion in Figures 2a, 3a and 4a is plotted as the distance along the vector of the main direction of flow. This vector is the best-fit normal line to the positions. A velocity time series was not computed due to errors of the individual positions and the noise associated with differentiating digital data. Two length scales of the average borehole water level are shown in Figures 2b and 3b. The local water level is the average of the nine boreholes at site 2880 (Fig. 1). The regional water level is the average of the mean water level at each of the five sites.

To understand this event within the framework of seasonal changes, we first investigate when the average borehole water level can be considered representative of a basal water pressure. We then examine the glacier motion, basal pressure and meteorological data in more detail during the pre-speed-up, speed-up and post-speed-up periods.

Basal water pressure

Determining basal water pressure from borehole water-level measurements is not straightforward. First, the boreholes act as manometers only if the volume of water at the bed is large in comparison with the borehole and any water flow through the borehole does not significantly affect the subglacial hydrology system (e.g. Reference Harper, Humphrey and GreenwoodHarper and others, 2002). Second, the water level in an individual borehole might not be representative of average basal conditions (e.g. Reference FountainFountain, 1994; Reference Gordon, Sharp, Hubbard, Smart, Ketterling and WillisGordon and others, 1998; Reference Fudge, Humphrey, Harper and PfefferFudge and others, 2008). For studies examining the relationship of basal water pressure and sliding, the relevant length scale for basal pressure is not known but has been suggested to be one to a few ice thicknesses based on the longitudinal stress-coupling lengths (e.g. Reference AndersonAnderson and others, 2004; Reference Howat, Tulaczyk, Waddington and BjornssonHowat and others, 2008). Therefore, a reliable estimate of average basal water pressure can be achieved only when many borehole water-level records over a large area show similar fluctuations.

The water levels of four different boreholes are shown in Figure 5. Two of the records are from site 2880 (16-borehole grid), whereas the other two are from sites 1730 and 1060. Any of the 15 boreholes could have been chosen, as all show a similar pattern (in fact, these four were chosen because they did not overlap extensively and could be discerned in the figure). Before the autumn speed-up, boreholes reached a stable water level, with diurnal fluctuations in water level lasting as late as day 260 in a couple of boreholes. During the speed-up, the water levels in all 15 boreholes follow the same basic pattern: a sudden increase and sustained high levels (270–272), a broad decline before a double-humped peak (272–275), and a final decline (275–276). On day 276, the borehole water levels no longer vary in unison as the water levels return irregularly to the high, stable levels typical of winter (Reference Fudge, Harper, Humphrey and PfefferFudge and others, 2005) over a period of days to weeks.

Fig. 5. Water levels from four boreholes at sites 2880 (thin solid and dotted curves), 1730 (thick curve) and 1030 (dashed curve). (a) Summer to winter transition; (b) autumn speed-up. All 15 boreholes show the same pattern of variations in water level during the autumn speed-up from day 270 to day 276. Only during this period are the average water levels considered representative of the average basal water pressure relevant to basal sliding.

All of the borehole water levels acted in concert during the speed-up (days 270–276) and are therefore representative of the basal water pressure on the length scales of both tens of meters and kilometers. Such widespread connectivity between boreholes is rare and is observed at only one other time: during the spring speed-up (Reference Harper, Humphrey, Pfeffer, Fudge and O’NeelHarper and others, 2005, Reference Harper, Humphrey, Pfeffer and Lazar2007). Note that the values of the water level and the magnitude of the variations in water level vary between boreholes. These differences are likely determined by the specifics of how each borehole is connected to the subglacial hydrology system, such as whether the borehole intersects the stoss or lee side of a bedrock bump. Therefore, the water-level averages might not yield the exact magnitude of the basal pressure field, but do accurately record how the basal pressure field varies through time.

Pre-speed-up

Based on glacier motion, the period prior to the speed-up event had two distinctly different modes (Fig. 2). Early in the record (days ∽220–245), ice velocity was 2.8 cm d–1. The diurnal variations in the average borehole water levels result from strong diurnal swings in a few boreholes while the majority showed a nearly constant water level. Large head-gradients often existed between closely spaced boreholes, and water-level fluctuations were not closely correlated between boreholes (Reference Fudge, Humphrey, Harper and PfefferFudge and others, 2008). Several significant warm-weather and rain events occurred during this period, but none had a discernible influence on either the pressure or velocity records.

Beginning on day 248, the ice displacement changed to a new mode of significantly slower horizontal motion of only 1.6 cmd–1. Coincident with the slow horizontal motion was downward vertical motion of roughly 30 cm (Fig. 2b). Later, during winter conditions (days 279–328), the horizontal velocity was a steady 2 cmd–1, and we assume, supported by calculations with Glen’s flow law, that this is the deformational velocity. The unusually slow velocity during the pre-speed-up period amounts to nearly 9 cm of motion due to internal deformation that was offset either by a change in the stress field or by reverse motion at the bed. No significant change in the borehole water-level records occurred on day 248, although the last two boreholes experiencing diurnal water-level fluctuations stopped doing so during this period. The accumulated rainfall in the 10 days after day 248 was 0.4 cm, only one-tenth of the amount of precipitation in the preceding 10 days. The temperature dropped significantly on day 255 and conditions remained cold and dry until the speed-up.

Speed-up (days 270–278)

A sliding acceleration began gradually on day 271. Figure 3 shows details of the horizontal positions, basal pressure and meteorological conditions during the speed-up event. The error bars on the horizontal positions represent the 95% confidence interval (CI). Based on the positional record, we partition the speed-up event into phase I (270.7–273.9) and phase II (273.9–278). During the first phase, rapid sliding lasted for ∽36 hours (271.1–272.6) and had a gradually increasing rate of motion. During the second phase, rapid sliding occurred over ∽16 hours (273.9–274.6) before leveling off. The second phase also included a short-duration reversal in the direction of motion at day 275.6. In fact, both phases ended with motion directed up-valley. Reverse motion has been observed infrequently on a number of glaciers (Reference WillisWillis, 1995) and recently studied in considerable detail at Gornergletscher, Switzerland (Reference Sugiyama, Bauder, Weiss and FunkSugiyama and others, 2007). We address its implications in the discussion.

We present the magnitude of reverse motion as the most likely value and include the range of values with a 95% CI in parentheses. The reverse motion at the end of phase I was 15cm (4–22 cm), during phase II it was 21cm (8–34 cm) and at the end of phase II it was 21 cm (11–31 cm); the durations were 16, 8 and 28 hours, respectively. We feel that the reverse motion at the end of both phases of the speed-up is robust because there are multiple measurements recording the progression of reverse motion. The reverse motion from 275.3 to 275.6 has a short duration (two measurements), but we find no reason to believe these positions less than the others.

The measurements of surface raising and lowering are limited by the large errors associated with the vertical positions. Uplift of 38cm (9–74 cm) occurred at the start of phase I (270.9–272.3), which roughly corresponds to the period of rapid sliding. Whether uplift occurred during the rapid sliding of phase II cannot be determined, but total accumulated downward displacement, after the phase I uplift until the end of phase II (272.3–277), was 62cm (26–98 cm). Specifics of the vertical motion during the individual episodes of reverse horizontal motion cannot be determined, as the errors allow anywhere from ∽20cm of upward to ∽50 cm of downward displacement (Fig. 4b).

The local and regional water-pressure records behaved similarly during the entire speed-up period. The higher pressures during both phases of the speed-up occurred when temperatures were high and/or rainfall was occurring, and the pressure troughs coincided with reduced temperature and rainfall. Both regional and local water pressures dropped at the end of phase I, with a greater drop in regional pressure to just under 70% of flotation pressure. During phase II, the pressure showed irregular peaks and then a large drop that reduced water pressure along the glacier to less than 65% of flotation pressure.

Post-speed-up (days 278–287)

After day 277, the glacier returned to steady motion and no other rapid-motion events occurred throughout the remainder of autumn. The temperature returned to near freezing and no rainfall was recorded, implying that the precipitation had turned to snow. The boreholes no longer experienced synchronous variations in water level.

Discussion and Conclusion

Water pressure and sliding speed

Our dataset provides an opportunity to compare surface velocity before, during and after a significant speed-up event with both locally averaged (tens of meters) and regionally averaged (several kilometers) measurements of borehole water levels. The first significant finding is that rapid sliding occurred only when the borehole water levels had similar fluctuations, implying a widespread connectivity of the subglacial hydrology system. During this period (days 270–276), the borehole water levels are representative of an average basal water pressure. At first inspection, the relationship between sliding speed and basal water pressure (Fig. 3a and b) shows close association. The complex multiphase speed-up with interludes of slow or reverse motion has associated pressure increases or decreases. During this event, the pressure is substantially more related to velocity than during the reoccurring ten-fold increase in speed we observe at Bench Glacier during the onset of spring melt conditions (Reference Harper, Humphrey, Pfeffer and LazarHarper and others, 2007). In spring, the rapid sliding occurs while the basal pressure either shows no change or decreases slightly.

Closer inspection of the data, however, reveals that the relationship between pressure and velocity is complex, with pressure events sometimes leading and sometimes lagging velocity events. Several complications in particular stand out.

1. Onset of rapid sliding:

Phase I: A pressure increase preceded the velocity increase by 4–8 hours.

Phase II: The start of rapid motion coincided with increasing pressure.

2. Sustained rapid sliding:

Phase I: Occurred while the pressure was steady and elevated (above pre-speed-up levels).

Phase II: Occurred while the pressure was increasing but low (below pre-speed-up levels).

3. End of rapid sliding:

Phase I: Occurred 4–8 hours after the pressure began dropping.

Phase II: Occurred as the pressure continued to rise for 8–12 hours.

4. Reverse motion:

End of phase I: Was concurrent with dropping pressure.

During phase II (275.6): Began about 8 hours after a dip in pressure.

End of phase II: Occurred while water levels were recovering from low values.

During the first phase of the speed-up, the water-pressure/sliding relationship appears sensitive to a critical pressure predicted by theory (Reference IkenIken, 1981; Reference Truffer and IkenTruffer and Iken, 1998). A preliminary pressure spike initiated sliding, which continued while the pressure remained elevated. However, the rapid sliding of phase II did not correspond to the pressure attaining a critical value because the pressure continued to increase after the sliding had stopped. During the two phases of the speed-up, the sliding speed was not consistently linked to increasing, decreasing or a critical water pressure. This is not a unique finding, as Reference Hanson, LeB Hooke and GraceHanson and others (1998) observed – at Storglaciären, Sweden – that velocity both preceded and lagged pressure changes. Hence, high or increasing water pressure may have been necessary for the rapid sliding, but was not sufficient.

Basal water storage is another key variable that influences glacier sliding (e.g. Reference Iken, Röthlisberger, Flotron and HaeberliIken and others, 1983; Reference MeierMeier and others, 1994; Reference Vieli, Jania, Blatter and FunkVieli and others, 2004). Water storage can be estimated by analyzing bed separation from measurements of surface uplift (e.g. Reference Howat, Tulaczyk, Waddington and BjornssonHowat and others, 2008) and by calculating meltwater input and outflow (e.g. Reference Bartholomaus, Anderson and AndersonBartholomaus and others, 2008). The surface uplift during the rapid sliding in phase I (270.9–272.3) suggests that bed separation and water storage are important during the autumn speed-up. However, the large errors in the vertical positions and the inability to constrain vertical strain with only a single measurement of velocity prevent a detailed analysis of bed separation here. With no record of the outlet stream discharge, we cannot estimate the water storage through the difference of inputs and outputs. Further, the borehole water-level records cannot straightforwardly be used to infer water storage, as Reference Harper, Humphrey, Pfeffer and LazarHarper and others (2007) showed that bed separation occurred without a corresponding increase in borehole water levels.

Synchronous variations of the basal pressure field occurred at only one other time during the year: the annual spring speed-up (Reference Harper, Humphrey, Pfeffer, Fudge and O’NeelHarper and others, 2005, Reference Harper, Humphrey, Pfeffer and Lazar2007). The synchronous behavior between 15 boreholes spaced kilometers apart implies that a large area of the basal drainage system was interconnected and had high transmissivity. A large amount of water at the bed promotes basal water connectivity, but connectivity is also promoted by enhanced sliding, which opens cavities at the bed. Connectivity in the spring commenced after an initial displacement of the ice (Reference Harper, Humphrey, Pfeffer and LazarHarper and others, 2007). Due to the compressed timing of the autumn event and fewer borehole measurements, the exact timing of the onset of connectivity relative to ice displacement is less clear during the autumn speed-up. Nevertheless, rapid sliding at Bench Glacier occurred only during periods of widespread connectivity of the basal pressure field. This suggests that basal connectivity is a key component of rapid sliding for Bench Glacier, and might also be critical for other glaciers.

Basal cavity dynamics

High melt and rainfall conditions similar to those observed during the autumn speed-up occurred earlier in the summer (days 215–225), yet no increase in motion was observed. A likely explanation is that the capacity of the drainage system was greater during summer and surface water input was routed through efficient basal conduits. The lack of water input during the 2 week period of cold and dry conditions (with a lowering solar angle) that preceded the autumn speed-up likely reduced the capacity of the drainage system. Indeed, the downward surface motion (and associated reduced horizontal motion) may be the result of an efficient subglacial drainage network deteriorating into an inefficient one. The large increase in water input caused by the warm, wet weather beginning on day 270 might then have overwhelmed the inefficient drainage network. The synchronous behavior of all 15 boreholes along nearly 2 km of the ablation area indicates that water was forced into a widespread and highly transmissive configuration, perhaps with a feedback from ice displacement. This was likely a system of very well-linked cavities.

The reverse motion at the surface may also be the result of subglacial cavity dynamics. The repeated association of the reverse motion with shutdown of rapid sliding may imply a loss of water storage in cavities if the magnitude of surface elevation loss can explain the 15–21 cm of horizontal reverse motion observed. With a bed slope at site 2880 of approximately 0.15, and assuming the water exerted a force normal to the bed, 1 m of elevation loss would be required to produce 15 cm of horizontal (backwards) motion. However, the water–ice contact in basal cavities is not necessarily parallel to the mean bed slope. Reference IkenIken (1981) showed that cavities on an idealized sinusoidal bed might occupy the space between the next peak of the sine wave and the next inflection point, which effectively steepens the water–ice contact. As the wavelength for a given amplitude of the sine wave decreases, or the cavity becomes more step-like, the slope of the water–ice contact steepens. A two- to three-fold increase in effective slope from the mean slope is plausible. At slopes of 0.3 and 0.45, we find that the downward surface motion required to achieve 15cm of horizontal (backwards) motion is 50 and 33 cm, respectively.

Although large errors in the vertical GPS measurements preclude a detailed analysis of elevation changes associated with the reverse motion, they do show downward surface motion of 62 cm (26–98 cm) from day 272 to day 277 (Fig. 4b). Approximately 15 cm of lowering is due to normal down-valley motion, so 47 cm (12–83 cm) can be attributed to a loss of bed separation if no other factors affected the vertical position. This is similar to the 33–50 cm, computed above, that would be necessary to explain 15 cm of reverse motion. This analysis shows only that the rapid drainage of water from cavities is a plausible mechanism for a significant amount of the reverse motion; it does not consider individual events (because of high measurement errors) or exclude other mechanisms. In fact, it is unlikely that all of the surface lowering is the result of changes at the bed because vertical strain is an important component of surface uplift measurements (e.g. Reference GudmundssonGudmundsson, 2002). Further, if the reverse motion occurred at the large end of our measured ranges, or the downward motion at the small end, then drainage of basal water can explain only a fraction of the reverse motion.

Elastic deformation and subglacial sediment deformation were proposed as other mechanisms for reverse motion by Reference Sugiyama, Bauder, Weiss and FunkSugiyama and others (2007). As Bench Glacier is underlain by little till, we can eliminate subglacial sediment deformation as a possibility. A detailed estimate of the elastic deformation potentially contributing to the horizontal reverse motion is difficult because we cannot calculate the longitudinal stress from a point-velocity measurement. However, assuming standard flow-law values from Reference PatersonPaterson (1994), and that the highest velocity observed (∽1md–1) occurred when another part of the glacier 1 km away did not move, gives a longitudinal stress of ∽105 Pa. Using 8.3×109 Pa for Young’s modulus for ice (Reference HookeHooke, 1998), we calculate that a stress of this magnitude could produce ∽0.01m of instantaneous elastic strain. The observed relaxation times of 8–28 hours imply that no more than 15% (∽0.002 m) of the strain can be recovered elastically (Reference Sugiyama, Bauder, Weiss and FunkSugiyama and others, 2007, equation 1). Therefore, we find it difficult to attribute any significant amount of reverse motion to elastic deformation.

The above calculations suggest that reverse motion smaller than 15cm can be plausibly explained by rapidly draining cavities and a loss of bed separation, implying that either the observed reverse motion was not at the high end of our possible range or other unaccounted-for processes were at play. If a loss of bed separation was the dominant mechanism of reverse motion, we note that only the reverse motion at the end of phase I occurred during declining water pressure. The reverse motion at the end of phase II coincided with increasing water levels, although variations were no longer synchronous between boreholes, and cavities might have become isolated. These results illustrate that reverse-motion events provide a new perspective on subglacial cavity dynamics during rapid sliding and highlight the need for more observations.

Acknowledgements

We thank the scientific editor, T. Murray, and two reviewers, D. Benn and T. Bartholomaus, for helpful comments that improved the manuscript. We also thank the many people who conducted fieldwork on Bench Glacier. The input of the glaciology group at University of Washington, particularly D. Power and A. Rasmussen, was of great help.

References

Anderson, R.S. and 6 others. 2004. Strong feedback between hydrology and sliding of a small alpine glacier. J. Geophys. Res., 109(F3), FO3005. (10.1029/2004JF000120.)Google Scholar
Bartholomaus, T.C., Anderson, R.S. and Anderson, S.P.. 2008. Response of glacier basal motion to transient water storage. Nature Geosci., 1(1), 33–37.CrossRefGoogle Scholar
Bradford, J.H. and Harper, J.T.. 2005. Wave field migration as a tool for estimating spatially continuous radar velocity and water content in glaciers. Geophys. Res. Lett., 32(8), L08502. (10.1029/2004GL021770.)CrossRefGoogle Scholar
Fountain, A.G. 1994. Borehole water-level variations and implications for the subglacial hydraulics of South Cascade Glacier, Washington State, U.S.A. J. Glaciol., 40(135), 293–304.CrossRefGoogle Scholar
Fountain, A.G. and Walder, J.S.. 1998. Water flow through temperate glaciers. Rev. Geophys., 36(3), 299–328.CrossRefGoogle Scholar
Fowler, A.C. 1987. Sliding with cavity formation. J. Glaciol., 33(115), 255–267.Google Scholar
Fudge, T.J., Harper, J.T., Humphrey, N.F. and Pfeffer, W.T.. 2005. Diurnal water-pressure fluctuations: timing and pattern of termination below Bench Glacier, Alaska, USA. Ann. Glaciol., 40, 102–106.CrossRefGoogle Scholar
Fudge, T.J., Humphrey, N.F., Harper, J.T. and Pfeffer, W.T.. 2008. Diurnal fluctuations in borehole water levels: configuration of the drainage system beneath Bench Glacier, Alaska, USA. J. Glaciol., 54(185), 297–306.CrossRefGoogle Scholar
Gordon, S., Sharp, M., Hubbard, B., Smart, C., Ketterling, B. and Willis, I.. 1998. Seasonal reorganization of subglacial drainage inferred from measurements in boreholes. Hydrol. Process., 12(1), 105–133.3.0.CO;2-#>CrossRefGoogle Scholar
Gudmundsson, G.H. 2002. Observations of a reversal in vertical and horizontal strain-rate regime during a motion event on Unteraargletscher, Bernese Alps, Switzerland. J. Glaciol., 48(163), 566–574.Google Scholar
Hanson, B., LeB Hooke, R. and Grace, E.M. Jr., 1998. Short-term velocity and water-pressure variations down-glacier from a riegel, Storglaciären, Sweden. J. Glaciol., 44(147), 359–367.CrossRefGoogle Scholar
Harper, J.T., Humphrey, N.F., Pfeffer, W.T., Huzurbazar, S., Bahr, D.B. and Welch, B.C.. 2001. Spatial variability in the flow of a valley glacier: deformation of a large array of boreholes. J. Geophys. Res., 106(B5), 8547–8562.Google Scholar
Harper, J.T., Humphrey, N.F. and Greenwood, M.C.. 2002. Basal conditions and glacier motion during the winter/spring transition, Worthington Glacier, Alaska, U.S.A. J. Glaciol., 48(160), 42–50.CrossRefGoogle Scholar
Harper, J.T., Humphrey, N.F., Pfeffer, W.T., Fudge, T. and O’Neel, S.. 2005. Evolution of subglacial water pressure along a glacier’s length. Ann. Glaciol., 40, 31–36.CrossRefGoogle Scholar
Harper, J.T., Humphrey, N.F., Pfeffer, W.T. and Lazar, B.. 2007. Two modes of accelerated glacier sliding related to water. Geophys. Res. Lett., 34(12), L12503. (10.1029/2007/GL030233.)CrossRefGoogle Scholar
Heroux, P. and Kouba, J.. 2001. GPS precise point positioning using IGS orbit products. Phys. Chem. Earth, 26(6–8), 573–578.Google Scholar
Hooke, R.LeB. 1998. Principles of glacier mechanics. Upper Saddle River, NJ, Prentice Hall.Google Scholar
Howat, I.M., Tulaczyk, S., Waddington, E. and Bjornsson, H.. 2008. Dynamic controls on glacier basal motion inferred from surface ice motion. J. Geophys Res., 113(F3).Google Scholar
Iken, A. 1981. The effect of the subglacial water pressure on the sliding velocity of a glacier in an idealized numerical model. J. Glaciol., 27(97), 407–421.CrossRefGoogle Scholar
Iken, A. and Truffer, M.. 1997. The relationship between subglacial water pressure and velocity of Findelengletscher, Switzerland, during its advance and retreat. J. Glaciol., 43(144), 328–338.Google Scholar
Iken, A., Röthlisberger, H., Flotron, A. and Haeberli, W.. 1983. The uplift of Unteraargletscher at the beginning of the melt season – a consequence of water storage at the bed? J. Glaciol., 29(101), 28–47.CrossRefGoogle Scholar
Jansson, P. 1995. Water pressure and basal sliding on Storglaciären, northern Sweden. J. Glaciol., 41(138), 232–240.Google Scholar
Joughin, I., Das, S.B., King, M.A., Smith, B.E., Howat, I.M. and Moon, T.. 2008. Seasonal speedup along the western flank of the Greenland Ice Sheet. Science, 320(5877), 781–783.Google Scholar
Kamb, B., Engelhardt, H., Fahnestock, M.A., Humphrey, N., Meier, M. and Stone, D.. 1994. Mechanical and hydrologic basis for the rapid motion of a large tidewater glacier. 2. Interpretation. J. Geophys. Res., 99(B8), 15,231–15,244.Google Scholar
King, M. 2004. Rigorous GPS data-processing strategies for glaciological applications. J. Glaciol., 50(171), 601–607.CrossRefGoogle Scholar
Lliboutry, L. 1968. General theory of subglacial cavitation and sliding of temperate glaciers. J. Glaciol., 7(49), 21–58.CrossRefGoogle Scholar
Magnùsson, E., Rott, H., Björnsson, H. and Pálsson, F.. 2007. The impact of jökulhlaups on basal sliding observed by SAR interferometry on Vatnajökull, Iceland. J. Glaciol., 53(181), 232–240.CrossRefGoogle Scholar
Mair, D., Nienow, P., Willis, I. and Sharp, M.. 2001. Spatial patterns of glacier motion during a high-velocity event: Haut Glacier d’Arolla, Switzerland. J. Glaciol., 47(156), 9–20.CrossRefGoogle Scholar
Meier, M. and 9 others. 1994. Mechanical and hydrologic basis for the rapid motion of a large tidewater glacier. 1. Observations. J. Geophys. Res., 99(B8), 15,219–15,229.Google Scholar
O’Neel, S., Echelmeyer, K.A. and Motyka, R.J.. 2001. Short-term flow dynamics of a retreating tidewater glacier: LeConte Glacier, Alaska, U.S.A. J. Glaciol., 47(159), 567–578.Google Scholar
Paterson, W.S.B. 1994. The physics of glaciers. Third edition. Oxford, Elsevier.Google Scholar
Röthlisberger, H. 1972. Water pressure in intra- and subglacial channels. J. Glaciol., 11(62), 177–203.CrossRefGoogle Scholar
Schoof, C. 2005. The effect of cavitation on glacier sliding. Proc. R. Soc. London, Ser. A, 461(2055), 609–627.Google Scholar
Sugiyama, S., Bauder, A., Weiss, P. and Funk, M.. 2007. Reversal of ice motion during the outburst of a glacier-dammed lake on Gornergletscher, Switzerland. J. Glaciol., 53(181), 172–180.CrossRefGoogle Scholar
Truffer, M. and Iken, A.. 1998. The sliding velocity over a sinusoidal bed at high water pressure. J. Glaciol., 44(147), 379–382.CrossRefGoogle Scholar
Truffer, M., Harrison, W.D. and March, R.S.. 2005. Correspondence. Record negative glacier balances and low velocities during the 2004 heatwave in Alaska, USA: implications for the interpretation of observations by Zwally and others in Greenland. J. Glaciol., 51(175), 663–664.CrossRefGoogle Scholar
Vieli, A., Jania, J., Blatter, H. and Funk, M.. 2004. Short-term velocity variations on Hansbreen, a tidewater glacier in Spitsbergen. J. Glaciol., 50(170), 389–398.CrossRefGoogle Scholar
Willis, I.C. 1995. Intra-annual variations in glacier motion: a review. Progr. Phys. Geogr., 19(1), 61–106.Google Scholar
Zwally, H.J., Abdalati, W., Herring, T., Larson, K., Saba, J. and Steffen, K.. 2002. Surface melt-induced acceleration of Greenland ice-sheet flow. Science, 297(5579), 218–222.CrossRefGoogle ScholarPubMed
Figure 0

Fig. 1. Map of Bench Glacier and measurement locations. Circles show boreholes; the star shows the GPS ‘rover’ receiver; the triangle indicates the position of the meteorological station and GPS ‘base’ receiver; the solid line indicates the approximate equilibrium line; dashed lines indicate elevation contours of 100 ft (∽30.5 m). Borehole sites are named by the distance, in meters, from the terminus. Regional water pressure is averaged from measurements at the five sites; local water pressure is averaged from measurements in the nine boreholes at site 2880.

Figure 1

Fig. 2. Data spanning the summer to winter transition. (a) Horizontal surface motion recorded by GPS at site 2880; (b) vertical surface motion recorded by GPS at site 2880; (c) local (thick curve) and regional (thin curve) average water level (see Fig. 1 caption); (d) temperature and rainfall. GPS positions processed with Trimble Geomatics Office; horizontal (∽10 cm) and vertical (∽20 cm) errors are not shown. The speed-up event (days 270–278) clearly stands out in all records. Prior to the speed-up (days 248–270), the horizontal velocity slowed and the surface elevation had downward motion. High temperatures and significant rain earlier in the year (days 217–227) failed to produce a similar speed-up. The average water levels are only considered representative of the average basal pressure from days 270–276 (see Fig. 5).

Figure 2

Fig. 3. Autumn speed-up event. (a) Horizontal surface motion recorded by GPS at site 2880; error bars represent 95% confidence interval; (b) local (solid curve) and regional (dashed curve) water pressure; horizontal dotted lines are pre-event pressure; (c) temperature and rainfall. GPS positions processed using PPP. The speed-up is divided into two phases. During phase I, rapid sliding occurred for ∽36 hours (first shaded box) during elevated pressure and after a pressure spike. During phase II, rapid sliding occurred for ∽16 hours (second shaded box) while the pressure was increasing. Velocity slowed during phase II before pressure reached a maximum.

Figure 3

Fig. 4. Reverse-motion events. Horizontal motion (a) and vertical motion (b) recorded by GPS at site 2880 processed with PPP. Error bars represent 95% confidence interval (CI). Shading shows times of reverse motion. The amount of reverse motion was 15 cm (4–26 cm, when considering 95% CI) at the end of phase I, 21 cm (8–34 cm) during phase II, and 21 cm (11–31 cm) at the end of phase II. Vertical motion cannot be determined for individual episodes because of the large errors. The total downward surface motion from days 272–277 was 62 cm (26–98 cm).

Figure 4

Fig. 5. Water levels from four boreholes at sites 2880 (thin solid and dotted curves), 1730 (thick curve) and 1030 (dashed curve). (a) Summer to winter transition; (b) autumn speed-up. All 15 boreholes show the same pattern of variations in water level during the autumn speed-up from day 270 to day 276. Only during this period are the average water levels considered representative of the average basal water pressure relevant to basal sliding.