Introduction
Surges of polythermal glaciers are thought to be initiated when cold regions of the glacier bed are raised to the pressure-melting point, and basal meltwater is produced more rapidly than it can be evacuated (Reference MurrayMurray and others, 2000; Reference A.C., Murray and NgFowler and others, 2001). Following surge initiation, evolution of the glacier flow field and patterns of surge propagation can follow a variety of pathways. Surges of land-terminating glaciers usually begin in the upper reservoir zone, then propagate down-glacier as a kinematic wave (e.g. Reference Meier and PostMeier and Post, 1969), whereas some tidewater glacier surges have been observed to begin on the lower tongue, then propagate upstream (e.g. Reference Luckman, Murray and StrozziLuckman and others, 2002; Reference Dowdeswell and BenhamDowdeswell and Benham, 2003). Peak ice velocities and rates of surge propagation also vary widely between glaciers. Modelling work by Reference A.C., Murray and NgFowler and others (2001) has identified some of the factors that are likely to control patterns and rates of surge propagation, although observational data to support further model testing and development are still very sparse.
In this paper, we focus on patterns of surge propagation on Bakaninbreen and Paulabreen, two confluent glaciers in Svalbard. Both glaciers have surged within the last two decades, but in each case the surge affected only one branch of the glacier. We present geomorphological, historical and speleological evidence for a persistent subglacial drainage channel between the two glaciers, and argue that efficient drainage along this corridor served to limit lateral surge propagation.
Bakaninbreen and Paulabreen
Bakaninbreen and Paulabreen are confluent glaciers that drain into the fjord of Rindersbukta, an arm of Van Mijenfjorden in central Spitsbergen (Fig. 1). In their lower reaches, the glaciers flow on either side of a prominent medial moraine extending northwest from their confluence below the nunatak Siggerudfjella. At the time of writing, most of the front of Paulabreen calves into Rindersbukta, whereas Bakaninbreen is entirely land-based. In 1898, the glaciers (which were then collectively known as Paula’s Glacier) terminated close to the mouth of Rindersbukta (Reference KjellströmKjellström, 1901; Reference OttesenOttesen and others, 2008). For the following century, the glaciers progressively calved back into the fjord, retreating ∽10 km between 1899 and 2003. Bakaninbreen and Paulabreen surged in 1985–95 and 2003–05, respectively. The surge of Bakaninbreen terminated before the surge front reached the terminus, which remained inactive. In contrast, the surge of Paulabreen resulted in a terminus advance of 1.9 km, temporarily reversing the retreat of ice in the fjord.
At the onset of the surge of Bakaninbreen, a steep ramp developed at the down-glacier boundary of fast-flowing ice, ∽8 km from the terminus (Reference Dowdeswell, Hamilton and HagenDowdeswell and others, 1991). The ramp was 25 m high during the early stages of the surge, but grew to 60 m high as the surge propagated down-glacier during the next few years. Between 1985 and 1989, the rate of surge-front propagation varied between ∽1 and ∽1.8 kma–1, but after 1989 the rate progressively declined and by 1994/95 was only 1.8–3.0ma–1 (Reference Murray, Dowdeswell, Drewry and FrearsonMurray and others, 1998; Fig. 2). Although the surge front propagated down-glacier for over 6 km and affected the whole width of the Bakaninbreen flow unit, it did not propagate laterally across the medial moraine onto Paulabreen. In addition, the medial moraine between the two glaciers underwent only very small lateral displacement during the surge, despite significant longitudinal strain of Bakaninbreen ice (Fig. 3a and b).
The mechanisms of surge propagation at Bakaninbreen have been studied in detail by Reference Murray and PorterMurray and others (1997, Reference Murray, Dowdeswell, Drewry and Frearson1998, Reference Murray2000) and Reference Murray and PorterMurray and Porter (2001). Radio-echo soundings in 1986 and 1996 and temperature measurements in 1994 and 1995 showed that the ice ahead of the surge front was frozen to the bed, whereas the rapidly flowing ice behind the front was warm-based, with the lowermost ∽40m of the glacier at or very close to the pressure-melting point. The glacier bed consisted of a mixture of glaciomarine muds and glacial till, at least 0.2 m thick. Reference MurrayMurray and others (2000) argued that these sediments did not completely thaw during the surge, and that basal water was trapped between the glacier sole and permafrost. The resulting elevated water pressures reduced the frictional strength of the bed, encouraging basal motion. Down-glacier propagation of the surge bulge was also associated with longitudinal compression and vertical thickening of the ice, part of which was accommodated by brittle failure. Reference Murray, Gooch and StuartMurray and others (1997, Reference Murray2000) documented a series of coherent englacial structures dipping up-glacier at around 21–35˚, both ahead of and behind the surge front, which were interpreted as thrust faults. Most of the structures terminated below the surface, but a few cropped out on the glacier surface, where they contained basal sediment and intermittently discharged pressurized water. Further evidence for brittle failure at the surge front was provided by analysis of seismic emission recordings made in 1987 (Reference Stuart, Murray, Brisbourne, Styles and ToonStuart and others, 2005). Three distinct types of emissions were identified, interpreted as the result of shallow faulting up-glacier of the surge front, brittle failure of cold ice, and resonance in water-filled cracks ahead of the surge front. Propagation of the surge was therefore associated with thawing of the basal ice and the uppermost part of the bed, and down-glacier extension of the region affected by high basal water pressures. Water intermittently escaped from the bed to the surface along thrust faults, encouraging lower basal water pressures and heat loss from the bed. As a result, this ‘leaky surge front’ appears to have limited the rate of surge propagation and contributed to surge termination (Reference MurrayMurray and others, 2000). Escape of water through gaps in the subglacial permafrost also appears to have played a role in surge termination (Reference Smith, Murray, Davison, Clough, Woodward and JiskootSmith and others, 2002).
The recent surge of Paulabreen originated in Skobreen, a tributary on the southwestern flank of the glacier (Fig. 1). The time of surge onset is not known, but the appearance of new crevasses in the upper part of Skobreen in 1990 suggests that the glacier may have undergone a long initiation phase prior to full onset (Reference SundSund, 2006). By summer 2003, most of Skobreen was affected by the surge, with a prominent surge bulge close behind the moraine loop where the glacier enters Paulabreen. At that time, ice in the upper part of Skobreen was drawn down and covered with transverse crevasses, and longitudinal crevasses were present on the surge bulge. Between 8 August 2003 and 3 April 2005 the surge propagated down-glacier and activated Paulabreen as far as its terminus. The mean rate of surge propagation was ∽7.8md–1, or ∽2.8kma–1, approximately double the maximum rate of surge propagation on Bakaninbreen. Measured on Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) images from July 2003 and July 2005, the moraine loop of the small glacier Sokkbreen (Fig. 1) was displaced ∽1.4 km, corresponding to a mean ice velocity of ∽1.6md–1, or ∽0.7km a–1. The surge front arrived at the terminus in early April 2005, and between then and August 2005 the glacier front advanced at a mean rate of 4.8 md–1. Surge termination occurred sometime between December 2005 and February 2006. The total advance was 1.9 km. Saturated till was widely exposed around the glacier margins at the end of the surge, and, given the common geological, geographical and climatic environment of Paulabreen and Bakaninbreen, it is likely that Paulabreen surged by the thermally regulated, soft-bed surging mechanism proposed by Reference A.C., Murray and NgFowler and others (2001) and Reference Murray, Strozzi, Luckman, Jiskoot and ChristakosMurray and others (2003).
Similarly to the 1985–95 surge of Bakaninbreen, the lateral margin of the Paulabreen surge closely followed the medial moraine. At the maximum extent of the surge, the boundary between surging and non-surging ice (delineated by the region of heavy crevassing) trended obliquely down-glacier from the confluence of Skobreen and Paulabreen to the medial moraine, whereupon it turned sharply down-glacier (Fig. 3d). For the first 2 km, the surge did not extend across the medial moraine, which remained in its pre-surge position. Farther down-glacier, however, the moraine was progressively displaced to the right (relative to glacier flow), with a maximum lateral displacement of 350m near the terminus. Crevasses extend across the displaced moraine, affecting a wedge-shaped area of Bakaninbreen with a maximum width of 150m at the terminus.
The two glaciers therefore behaved independently during their recent surge history. The 1985–95 surge was confined entirely to Bakaninbreen, and Paulabreen was not affected. Conversely, the 2003–05 surge affected only Skobreen and Paulabreen. The fact that Paulabreen was not activated during the Bakaninbreen surge is particularly notable, since the former glacier was dramatically activated by the surge of a tributary glacier only 10 years later. The independent surge behaviour of Bakaninbreen and Paulabreen thus strongly suggests that the medial moraine delineates the position of an effective barrier between them, which acts to constrain the lateral propagation of surging ice.
Evidence for a Subglacial Conduit Between Bakaninbreen and Paulabreen
A 250 kHz swath bathymetric survey of Rindersbukta (Reference OttesenOttesen and others, 2008) revealed a major meandering ridge on the fjord floor, extending from close to the present glacier terminus to close to the 1898 glacier limit (Fig. 4). The ridge is sharp-crested, 5–10m high and is overprinted by recessional moraine ridges. A second, confluent ridge trends northward from the southwest shore of Rindersbukta, just south of the point where the former tributary glacier Scheelebreen flowed into the fjord. Historical maps and photographs of Paulabreen show that the positions of both ridges coincide with persistent deep embayments in the ice front, varying from narrow notches up to 900 m long to more open V-shapes. In our experience, embayments in the fronts of tidewater glaciers in Svalbard are typically associated with upwellings of turbid meltwater, marking the location of effluxes from subglacial conduits. Upwelling water undercuts the ice margin, locally enhancing calving and margin retreat. The association between the embayments and upwellings is supported by observations made by Reference CösterCöster (1925), who reported water emerging from the back of two embayments in the margin of Paulabreen in 1924. The correspondence between the location of the fjord-floor ridges, ice-front embayments and observed water efflux points leads us to conclude that the ridges are eskers, and mark the position of persistent subglacial conduits (see also Reference OttesenOttesen and others, 2008). The efflux points of the main conduit (i.e. the deepest point of the embayments) were consistently located at the terminus of the medial moraine between Paulabreen and Bakaninbreen.
To determine whether a subglacial conduit exists at the upper end of the Bakaninbreen–Paulabreen medial moraine, the glacier confluence was visited in September 2007. A conduit entrance was discovered in the floor of a large hollow on the glacier surface at the up-glacier end of the moraine (marked ‘C’ in Fig. 1). The conduit was entered in April 2008. A steep drop of 10m led directly to the upper end of a Nye channel cut in bedrock, which consisted of thin-bedded sandstones and shales (Fig. 5a). The lowermost 1–2m of the glacier ice had high debris content (including clasts and bands of silt ‘clots’), and was frozen to the bed. The Nye channel was box-shaped, 3–5m wide and typically 4 m deep. The height of the ice roof above the top of the Nye channel varied from ∽1 to ∽10m (Fig. 5b). The conduit was followed for 100 m downstream until progress was abruptly halted by the collapse of thin ice overlying deep water. The accessed part of the channel had a meandering planform and a stepped long profile, with gently sloping sections separated by steep steps up to a few metres high.
Access was also gained to two feeder passages on the Paulabreen side, one near the upper end of the Nye channel and the other ∽90m down-glacier of the entrance. Both passages consisted of narrow, incised canyons graded to the Nye channel. The upper passage was surveyed for 200m and took the form of a meandering canyon typically 10m deep. For much of its length, the passage floor was incised ∽1–2m into the glacier bed. The roof consisted of masses of snow or icicles plugging the upper levels of the canyon, or ‘canyon sutures’ formed where opposing walls have been brought into contact by ice creep. This passage continued, but was not surveyed due to time constraints. The lower feeder passage averaged 0.5m wide and 7m high. Our survey was terminated ∽15m from the main channel at a vertical nick point, but the passage continued above. The two feeder passages are interpreted as ‘cut-and-closure’ conduits formed by the incision of supraglacial streams followed by closure of the upper levels by a combination of blown snow, meltwater refreezing and ice creep. The characteristics of cut-and-closure conduits have been described in detail by Reference Gulley, Benn, Müller and LuckmanGulley and others (2009). In September 2007, water was also seen to enter the head of the Nye channel from a supraglacial channel draining the Bakaninbreen side. The catchment of this channel is small, however, and most water entering the Nye channel is routed via the upper feeder passage on the Paulabreen side.
The presence of water in the conduit in April 2008, in very cold weather, leads us to conclude that water may persist year-round, preventing complete conduit closure during the winter. This means that the conduit can be reactivated each summer when surface meltwater enters its upper end, even though at least part of the adjacent bed is below the pressure-melting point. We have frequently encountered englacial and subglacial water pockets in cold ice in wintertime in Svalbard, and attribute its presence to either solute enrichment during partial freezing (cf. Reference Hodson, Mumford, Kohler and WynnHodson and others, 2005) or inefficient heat loss by conduction into the surrounding ice and bed (cf. Reference Murray and PorterMurray and Porter, 2001). It is notable in this respect that the subglacial conduit persisted during a long period during which both Paulabreen and Bakaninbreen were quiescent, when much of the bed was presumably below the pressure-melting point.
Discussion
According to Reference MurrayMurray and others (2000), surge-front propagation at Bakaninbreen occurred by (1) strain thickening of the ice ahead of the surge front in response to compressive stresses, by a combination of ice creep and overthrusting;
(2) warming of the ice and bed by strain heating; (3) injection of pressurized water into the bed ahead of the surge front, increasing basal lubrication; and (4) warming of the bed by refreezing of water injected into thrust faults. Retention of pressurized meltwater at the glacier sole is a necessary condition for maintaining fast flow (e.g. Reference KambKamb and others, 1985; Reference A.C., Murray and NgFowler and others, 2001; Reference Harrison and PostHarrison and Post, 2003), and slowing of surge-front propagation and surge termination appear to have resulted from the release of water along thrust faults and through the bed.
We hypothesize that the conduit below the Bakaninbreen–Paulabreen medial moraine played a key role in constraining the lateral extent of the surges of the two glaciers. When pressurized water beneath surging ice reached the conduit beneath the medial moraine, basal water could have been discharged efficiently to the glacier terminus. This could have occurred during all seasons of the year, if we are correct in concluding that the conduit is a perennial feature, retaining stored water throughout the winter. By allowing water to escape along the bed below the medial moraine, the conduit would have acted as a low-pressure subglacial corridor between the two glaciers, switching off surge propagation mechanisms (3) and (4). It is notable that the only place where the surge propagated from one flow unit to the other and the moraine displaced laterally was on the lower part of Bakaninbreen during the 2003–05 Paulbreen surge. At this locality, the surface of Bakaninbreen is ∽20–30m lower than that of Paulabreen, so that transverse compressive stress gradients across the moraine would have been much larger than farther up-glacier where the two glacier surfaces are closer in elevation. Limited propagation of the surge across the medial moraine could therefore occur in this zone by mechanisms (1) and (2). Alternatively, the conduit may have been less efficient at evacuating water from the lower glacier, decreasing its effectiveness as a pressure barrier.
If we are correct, it may be asked why discharge of water along the conduit did not lead to surge termination, as appears to have been the case for the ‘leaky surge front’ described by Reference MurrayMurray and others (2000). One possible explanation is that significant water losses to the conduit were confined to a relatively narrow zone, and that elsewhere below the surging part of the glacier production of meltwater by strain heating was able to maintain high basal water pressures.
We believe that the presence of a subglacial conduit below the medial moraine is not fortuitous, because confluences in glacier ablation zones can be expected to be favourable sites for surface-to-bed drainage. Reference GudmundssonGudmundsson (1997) has shown that ice at confluences is subject to transverse compression and longitudinal extension, the former because of flow convergence, and the latter because ice is advected away from the no-slip lateral boundary into a mid-stream position where the restraining effect of lateral drag is reduced. In addition, water supply tends to be relatively high at glacier confluences, where lateral supraglacial drainage systems converge. The combination of stressed ice and a plentiful water supply is conducive to hydrologically assisted fracturing, allowing surface crevasses to penetrate to the glacier bed (Reference Alley, Dupont, Parizek and AnandakrishnanAlley and others, 2005; Reference Van der VeenVan der Veen, 2007; Reference Catania, Neumann and PriceCatania and others, 2008; Reference Benn, Gulley, Luckman, Adamek and GlowackiBenn and others, inn press). Alternatively, water can leak along the ice–bed interface below the confluence zone (Reference Huss, Bauder, Werder, Funk and HockHuss and others, 2007). It is unclear whether surface-to-bed drainage at the Bakaninbreen–Paulabreen confluence originally formed by hydrologically assisted fracturing or leakage along the bed. Extensive hoar frost on the passage walls obscured much of the structural detail, making interpretation difficult. Examination of surface structures on aerial photographs, however, supports the idea that water initially entered the glacier down an overdeepened fracture, then incised headward to create the observed feeder passages.
Because of the association between glacier confluences and surface-to-bed drainage, we anticipate that other examples may be found of independently surging flow units within confluent glacier systems. Interestingly, during the 1992/93–95 surge of Sortebræ, East Greenland, surging activity was confined to one flow unit in the early part of the event but spread across the whole glacier at a later stage (Reference Pritchard, Murray, Luckman, Strozzi and BarrPritchard and others, 2005). The reasons for this pattern are not known, although it suggests that barriers to surge propagation may be effective in some circumstances but break down in others.
Acknowledgements
We are grateful for funding and logistical support from The University Centre in Svalbard (UNIS) and Store Norske Spitsbergen Kulkompani (SNSK). Additional funding was provided by Sigma Xi, the Evolving Earth Foundation, the Cleveland, OH, Grotto Science Fund and the US National Speleological Society. We also thank A. Bergström for assistance during exploration of the subglacial conduit.