Introduction

Surge-type glaciers oscillate between periods of rapid movement and periods of relative inactivity or quiescence (Reference PostPost, 1969). The exact mechanisms driving this flow instability are not fully understood, though recent studies suggest that surging is likely a result of: changing conditions in the sediments at the ice–bed interface (e.g. Reference Clarke, Collins and ThompsonClarke and others, 1984); changing thermal conditions (e.g. Reference RobinRobin, 1955; Reference ClarkeClarke,1976; Reference Fowler, Murray and NgFowler and others, 2001); and/or changing hydrological conditions (e.g. Reference KambKamb, 1987). What is increasingly apparent is that the complexities of surge initiation, the mechanisms for sustaining fast flow and surge termination cannot all be explained by any single model. To understand better the dynamics of surge-type glaciers, more direct observations are needed of the surge cycle for a range of surging glaciers. This is particularly difficult in remote areas, where surges are difficult to monitor, and in Arctic areas, where the period of quiescence is often >150 years (Reference Dowdeswell, Hamilton and HagenDowdeswell and others, 1991). Alternative methods must be developed to investigate the development and termination of surges in such areas. A number of studies have attempted to identify surge-type glaciers using diagnostic sedimentary structure in the glacier body and the glacier forefield (e.g. Reference SharpSharp, 1985; Reference Evans and ReaEvans and Rea, 1999). If certain sedimentary structures can be uniquely linked to different phases of the surge then identification of these structures may allow us to recognize different modes of surging. The problem is that sedimentary structures are only available for study at the glacier surface and at the glacier front. Ground-penetrating radar (GPR) is one method for extending the study area to permit suites of structure to be identified within the ice matrix.

GPR is a geophysical technique developed for the non-invasive investigation of subsurface features. Electromagnetic wave energy from a transmitter unit on the ground surface is partially reflected from subsurface features exhibiting electrical properties contrasting with those of the surrounding matrix (Daniels and others, 1988). In glacier ice, marked contrasts in electrical properties are a result of changes in the temperature of the ice (e.g. where a warm ice layer exists below cold surface ice (Reference BjörnssonBjörnsson and others, 1996)), changes in water content (Reference Murray, Stuart, Fry, Gamble and CrabtreeMurray and others, 2000a, Reference Murrayb), changes in bubble content of the ice (Reference Macheret and ZhuravlevMacheret and Zhuravlev, 1982), changes in the chemical composition of the ice (Reference GudmandsenGudmandsen,1975), or varying debris concentration in the ice (e.g. Reference Arcone, Lawson and DelaneyArcone and others,1995). Radar can also be used to map features formed as a result of brittle fracture in the ice, such as crevasses (e.g. Reference Jezek, Bentley and CloughJezek and others,1979; Reference Glover and ReesGlover and Rees, 1992), faults (Reference Murray, Gooch and StuartMurray and others, 1997) and shear zones (Reference Nobes, Leary, Hochstein and HenrysNobes and others,1994).

Few GPR surveys have been carried out to investigate structures produced by the incorporation of sediments in surge-type valley glaciers. One glacier that has been studied is Bakaninbreen, Svalbard. Basal and englacial sedimentary features were mapped, many of which were interpreted to be sediment-filled thrust features 0.6–1.1 m thick, although displacement could only be demonstrated across one feature (Reference Murray, Gooch and StuartMurray and others, 1997, Reference Murray2000b). Some of these englacial features correspond to surface outcrops mapped by Reference Hambrey, Dowdeswell, Murray and PorterHambrey and others (1996), and they appear to transport considerable volumes of sediment from the basal zone (Reference MurrayMurray and others, 2000b).

Previous investigations at the surge-type glacier Kongsvegen, Svalbard, have identified suites of englacial sedimentary structures. How this sediment is incorporated into the ice is currently unresolved. It could be by thrusting (Reference Bennett, Hambrey, Huddart and GhienneBennett and others, 1996; Reference Hambrey, Bennett, Dowdeswell, Glasser and HuddartHambrey and others, 1999), folding sub-parallel to glacier flow (Reference Glasser, Hambrey, Crawford, Bennett and HuddartGlasser and others, 1998; Reference Hambrey, Bennett, Dowdeswell, Glasser and HuddartHambrey and others, 1999) or crevasse filling (Reference Woodward, Murray and McCaigWoodward and others, 2002). This paper details GPR investigations of the sediment structures near the terminus of Kongsvegen. GPR results will be used to develop models of sediment emplacement and models for the surge of Kongsvegen. First, the ability of GPR to image structure in the ice of Kongsvegen will be investigated. Second, mechanisms for the development of the structure identified using GPR will be developed. These mechanisms will then be used to construct a model for surge initiation, the period of fast flow and then surge termination.

Field Site

Kongsvegen (78°48′ N, 12°59′ E) is a 27 km long, 189 km² surge-type glacier in northwest Spitsbergen (Reference Hagen, Liest l, Roland and JørgensenHagen and others, 1993) (Fig. 1), which currently terminates in a partially grounded and partially calving front. The grounded section of the front forms a cliff face 5–20 m high running northwest–southeast for >1 km along the fjord margin (at 160–185° to the ice-flow direction). The glacier surface profile parallel to the cliff face dips at ∼2° along the fjord, and in its lower reaches is almost horizontal perpendicular to the ice cliff. The calving margin of Kongsvegen is confluent with the principal flow unit of Kronebreen and Infantfonna at the head of Kongsfjorden. During periods of advance, the complex is confluent with another stream of Kronebreen to the north of Colleth gda. The Kongsvegen flow stream is currently quiescent, with a flow rate of ∼3 m a⁻¹ (Reference Melvold and HagenMelvold and Hagen, 1998). By contrast, Kronebreen is currently one of the fastest-flowing tidewater glaciers in Svalbard, with a flow speed of ∼750 m a⁻¹ (Reference Lefauconnier, Hagen and RudantLefauconnier and others, 1994; Reference Melvold and HagenMelvold and Hagen,1998).

Fig. 1. (a) The tidewater margin of Kongsvegen, Svalbard, showing the position of the margin from 1948 and 1990 aerial photographs and the current dominance of Kronebreen at the tidewater margin. The box marks the location of the photograph shown in (b). (b) Oblique photograph, with the location of the GPR grid marked.The grounded cliff face running along the margins of Kongsfjorden and the proglacial sediment structures in the marginal area of Kongsvegen are visible. The inset box shows details of the GPR grid, arrows marking the direction of survey for each of the lines.

Interpretation of aerial photographs suggests that the last surge of Kongsvegen occurred between 1936 and 1956 (Reference Woodward, Murray and McCaigWoodward and others, 2002). Oblique aerial photographs from 1936 show the catchment in a pre-surge state, with the medial moraine complex separating Kronebreen from Infantfonna and Kongsvegen close to the Botnfjellbreen margin of Kongsvegen, as it is today. By comparison, the 1948 photograph shows a wider and more heavily crevassed Kongsvegen flow unit. The 1948 photograph also shows the maximum recorded advance of the tidewater margin, thought to be a result of the push effects of the advance of the surge front propagating down the Kongsvegen flow unit (Reference Woodward, Murray and McCaigWoodward and others, 2002). By 1956 the Kongsvegen flow unit covered two-thirds of the joint tidewater margin. Subsequently the tidewater complex has retreated, Kongsvegen has entered quiescence and Kronebreen has gradually widened to almost completely squeeze out the Kongsvegen tongue (Reference Melvold and HagenMelvold and Hagen,1998).

Methodology

Data collection

A grid of 200 MHz GPR lines was collected along the margin of Kongsvegen during July 1996 (Fig. 1) using a Pulse EKKO 100 GPR system. Two types of survey were undertaken, common-offset (CO) and common-midpoint (CMP). During CO surveys the antennae, fixed perpendicularly to the survey line on a plastic sledge, were moved progressively along the profile in step mode, collecting traces at 0.25 m intervals along the profile. During CMP surveys, the antennae were moved progressively further apart at an equal offset from the central image point. Reflections from CMP surveys were used to calculate the velocity of propagation of electromagnetic energy in the subsurface. This allows a depth scale to be added to processed survey lines. Three 300 m long parallel lines, approximately 11 m apart, were collected parallel to the grounded cliff face at the glacier terminus. These lines are denoted low line (LL), mid-line (ML) and upper line (UL) (Fig. 1b). Five cross-lines, each approximately 65 m long, were also surveyed, spaced every 55 m. These lines are denoted XLA–XLE. During data collection, the locations of surface water and streams, and the angle at which any visible trace or sediment-filled fracture crossed the GPR lines, were noted, along with the location of the structure, and where possible an estimate of its dip.

GPR profiles were processed using GRADIX 1.10 software (Interpex Ltd). The processing sequence was optimized using CO profile LL (Fig. 2), then applied to all profiles. Drift correction and time-zero correction were applied, followed by dewow (4 MHz highpass) and bandpass filtering (trapezoidal filter set at 30–60–250–500 MHz) to remove noise at the high and low ends of the spectrum. Background removal in a window of 1000 traces was applied to remove coherent horizontal banding from the radar profile. Gain (either a constant gain or a time-varying automatic gain adjusted within 40 ns blocks) was applied to enhance returns from deeper reflectors whose signal strength is diminished by losses from attenuation, geometric spreading, scattering and energy partitioning at shallow interfaces. Depth and elevation correction allow the finished output to match surface topography with a true depth scale.

Fig. 2. Processing flow for GPR data. A description of processing parameters is given in the text. Migration has not been applied to all GPR profiles.

Reflections from near-horizontal layers were identified from CO profiles, and these reflectors were then picked manually from CMP surveys. This allowed calculation of radar wave velocities of 0.168 ± 0.002 m ns⁻¹. CMP velocities were used to apply depth and elevation corrections.

All profiles were also migrated using an F-K migration algorithm with a constant velocity of 0.168 m ns⁻¹, though not all profiles in this paper are displayed with migration applied. As the GPR grid was not dense enough for the reliable application of three-dimensional migration algorithms, GRADIX’s two-dimensional algorithm was used, which assumes all energy is returned from in-line reflectors. However, the profiles did not cross all features at 90° to their trend, thereby violating the two-dimensional assumption. Furthermore, the velocity of 0.168 ± 0.002 m ns⁻¹ was derived from propagation in cold ice, and is therefore inappropriate for migrating parts of the profile containing sediments or warm ice. GRADIX supports migration using only one-dimensional velocity models, but the dominant reflectors are not horizontal, and lateral variations in velocity are evidently present. Thus migration has not always successfully focused the reflection images where there are apparently large volumes of included sediment, or in the deeper (warmer) parts of the section.

Results

The GPR profiles can be divided into three distinct sequences by two high-amplitude, sub-horizontal reflectors denoted SHR1 and SHR2 (Fig. 3). On LL, ML and UL, SHR1 dips up-glacier at ∼3° and has opposite phase to the direct wave, indicating a higher dielectric permittivity (lower velocity) for the layers below the reflector than for those above it. The reflector is continuous except for occasional small steps in the profile and disturbances where steeply dipping structures rise from the reflector (e.g. at 93.3–96.5 m on LL). On the cross-lines XLA and XLB the reflector dips towards the fjord at 3–5°, less steeply than the trend of the glacier surface. On XLC–XLE the reflector dips towards the fjord at a slightly steeper angle, > 5°, and is approximately parallel to the surface.

Fig. 3. GPR profiles collected from the grid shown in Figure 1. The low line (LL) and the three cross-lines XLA, XLC and XLE are displayed. Processing parameters are as recorded in Figure 2, though profiles are unmigrated and have constant gain. Two prominent sub-horizontal reflectors, SHR1 and SHR2, are identified on LL, as are surface scatter (S), dipping structures above SHR1 (DS) and shallow reflectors (SR) below SHR1. Flow direction is from left to right for LL, and out of the page for XLA, XLC and XLE.The depth-scale bar (in m) on LL applies to all parts of the figure.

The second sub-horizontal reflector, SHR2, also dips down-glacier at 0.5–3.5° and has opposite phase to the direct wave. From 0 to 80 m on the up-glacier orientated lines, the reflector is discontinuous and appears to be made up of a number of sub-parallel reflectors. From 80 to 240 m the reflector is more continuous and more steeply dipping, before again becoming somewhat discontinuous after 240 m as the layer rises to meet and merge with SHR1. On XLA–XLD the reflector dips towards the fjord parallel to the general trend of the glacier surface. On XLE the reflector is barely visible but seems to be slightly below SHR1.

The near-surface region of all the profiles is dominated by shallow diffraction patterns, which originate at the glacier surface and produce scatter (noise) on the GPR profile (Fig. 3). Some of these features appear to be symmetric and migrate successfully, while others are asymmetric and cannot be migrated successfully, leaving migration “smiles” and residual noise as the most obvious signal of the shallow structures in migrated sections.

The section of the lines below the near-surface and above SHR1 is dominated by asymmetric diffraction patterns suggesting steeply dipping structures, occurring every 10–15 m, and similar to those shown in Figure 4. Before migration, these reflectors seem to stem from a diffraction pattern with a straight limb on the up-glacier side. Below the diffraction pattern are a series of hyperbolae running towards the bed. When migrated, the upper hyperbola becomes a point reflector with a steeply dipping linear structure running from the point reflector towards the bed (Fig. 4a). The linear part of the reflector has moved up-dip and steepened considerably. The hyperbolae collapse to point reflectors along the plane of the structure. After migration, features do not appear to cross SHR1, while migration smiles below the reflector are a common feature, presumably due to the difficulties of the constant velocity model as described above, or because the reflections from the structures are not truly hyperbolic.

Fig. 4. Dipping reflectors above SHR1 seen in migrated sections (a, c) and in the field (b, d). (a) After migration, the reflector on LL shows a strong, steeply dipping structure, with migration smiles on the up-glacier side of the feature.The reflector shadows any reflection from SHR1 and has no expression below SHR1. From field observations, the feature can be identified as the sedimentary inclusion on the cliff face pictured in (b). The sedimentary structure runs parallel to sedimentary mounds melting out in the forefield. (c, d) Another steeply dipping structure on LL which outcrops on the cliff face.The feature is composed of a series of steeply dipping, sediment-rich fractures, with shallow-angled structures diverging from the main fractures.

Below SHR1 a series of shallow, 20–30°, up-glacier dipping structures are apparent, which continue below SHR2. A number of structures, such as that dipping up-glacier from 104 m on LL, are visible before migration and remain prominent post-migration (Fig. 5).

Fig. 5. Dipping reflectors below SHR1: (a) migrated data from LL; (b) interpretation of the migrated GPR plot. The reflector crosses SHR2, where it appears to steepen slightly towards the surface, suggesting that the speed of propagation of electromagnetic waves is slower below SHR2 than above it. These structures are much shallower than the steeply dipping structure found above SHR1, are weaker reflectors and are much less frequent on the GPR traces.

Interpretation of Structure in the GPR Grid

The descriptions of the GPR lines indicate that there are a number of sets of structure with very different radar signatures (Fig. 3). These features are: near-surface diffraction patterns (S); two major sub-horizontal reflectors (SHR1 and SHR2); steeply dipping features (DS) above SHR1; and localized shallow-dipping features (SR) below SHR1. The near-surface diffraction patterns (S in Fig. 3) are produced by surface features, such as streams, sediment-filled deformation structures and cracks. Many of the diffraction patterns are asymmetric, as GPR lines do not cross the structures perpendicularly and because the structures are dipping.

Glacier Surge Mechanisms

The GPR results (summarized in Fig. 6) allow a model for structural development during the surge of Kongsvegen to be proposed:

Fig. 6. Schematic summary of GPR results. The subsurface has been divided and separated by the dominant sub-horizontal reflectors SHR1 and SHR2. In the glacier ice above SHR1, schematic sedimentary structures are shown. SHR1 represents the glacier bed, and none of the sedimentary structures appear to cross this interface. Below the bed a frozen wedge of sediment is found, the base of which is defined by SHR2.The sediment is frozen to the base of the glacier by the winter cold wave which penetrates to a depth of ∼14 m into glacier ice.The second sub-horizontal interface, SHR2 represents the interface between ice and sediment frozen by the winter cold wave, and underlying, unfrozen, possibly saturated sediment. Some structures appear to cross SHR2. These are interpreted as thrust planes in the subglacial sediments.

Fig. 7. Schematic showing the development of structure during the surge of Kongsvegen. (a) As the surge front arrives, ice and sediment are thrust due to compression between fast-moving surgingice up-glacierand stagnant icedown-glacier. A water and sediment pulse moves down-glacier withthesurge, producing décollement between the ice and sediment. (b) Deforming basal sediments allow fast flow during the surge as the flow regime becomes extensional. (c) Sediments injected into fractures in the ice during the extensional flow regime produce crevasse-fill structures and purge the bed of saturated sediments duringsurge termination.

While it is recognized that the GPR data used to construct this model cover a small spatial area of the glacier surface, the results from the GPR are consistent with other observations, suggesting that results can be extrapolated across the catchment. The model has two important implications for the development of the structures in the ice of Kongsvegen. First, it suggests that thrust faulting only occurred during the initial compressive regime. Thrust faults are likely to be difficult to recognize in the current ice, as they will have been moved down-glacier from their expression in the subglacial sediments and also will have been reorientated by velocity gradients within the ice during and after the surge (Reference Woodward, Murray and McCaigWoodward and others, 2002). Furthermore, the likely inability of compressive thrust planes to entrain large volumes of sediment suggests they may be difficult to identify as isolated structures in the GPR profiles for the most part. Second, we suggest large amounts of sediment were injected into the ice of Kongsvegen into voids opened by extensional crevasses. Thus most of the structure visible in the ice in the marginal areas of Kongsvegen is thought to originate from crevasse filling.

Studies of surge-type glaciers, which oscillate between short periods of extremely rapid movement and long periods of relative inactivity or quiescence, suggest that fast flow is associated with basal conditions (Reference Clarke, Collins and ThompsonClarke and others, 1984; Reference KambKamb, 1987; Reference RaymondRaymond, 1987). In the case of surge-type glaciers overlying a sedimentary bed, a rapidly deforming, saturated, sediment substrate may induce and/or support surging (Reference Clarke, Collins and ThompsonClarke and others, 1984). Surge-type glaciers overlying sedimentary beds often entrain large amounts of sediment (Reference ClappertonClapperton, 1975). Therefore, the mechanisms by which sediment is incorporated into the ice during surge may play a significant role in the speed of termination of fast flow by controlling the time taken to reduce the availability of basal sediment and water. The model proposed above suggests that some of the mechanisms that drove the surge of Kongsvegen can be elucidated from studies of the englacial and subglacial sedimentary structures. The termination of the surge may have been assisted by the change in basal conditions when saturated sediments, which had supported fast flow, were injected into open crevasses.

Conclusions

GPR investigations near a grounded cliff face at the margins of Kongsvegen have enabled the identification of a number of internal reflectors, which represent sedimentary structure within the glacier and the underlying sedimentary bed. Reflectors include: the glacier bed; an underlying thermal interface between frozen and unfrozen sediments; steeply dipping structures imaged above the bed, which are interpreted as crevasse fills and correlate with sediment-rich features seen on the cliff face; and shallow structures visible in the sediment below the bed, which are interpreted as thrust planes formed during the compressional phase of the surge.

The GPR results have allowed the development of a model for structural development during the surge of Kongsvegen. During the 1940s, surge ice moved from a reservoir, propagated down-glacier into slow-flowing ice, creating a zone of compression which resulted in thrusting. Remnant thrust faults are recorded in the subglacial sediments. With the arrival of fast-flowing ice during the surge, the sedimentological and hydrological regime at the basal interface changed rapidly, and the ice decoupled from the basal sediments. Rapid flow during the surge resulted in the extensive fracture of the ice. Sediment is thought to have been injected into these extensional crevasses from the glacier bed, exhausting saturated sediments from the basal interface, possibly accelerating surge termination.

As such, this model suggests that thrusting was an important occurrence during the surge of Kongsvegen, but the thrusting was of little consequence in developing sedimentary landforms in the glacier ice or forefield environment. This suggests that the majority of sedimentary structures visible in the glacier ice and forefield were produced by crevasse filling of extensional crevasses opened during the surge, which have been extensively reorientated by flow during and after the 1940s surge of Kongsvegen.