Introduction
Storglaciären (Fig. 1) is a small, predominantly temperate valley glacier in northern Sweden that has been studied since Reference SchyttSchytt (1959, Reference Schytt1966, Reference Schytt1968) began the still-continuing mass-balance measurements on it in the spring of 1946. Recent work has included velocity measurements over periods as short as one day (Reference Hooke, LeB, Hooke, Brzozowski and BrongeHooke and others, 1983; unpublished data), bore-hole deformation surveys (Reference Hooke, Hooke, Holmlund and IversonHooke and others, 1987; paper in preparation by V. Pohjola and R. LeB. Hooke), sensitive strain and tilt measurements (paper in preparation by P. Jansson and R. LeB. Hooke), water-balance studies (Reference Östling and HookeÖstling and Hooke, 1986), and dye-trace experiments (Hooke and others, in press; Seaberg and others, in press).
The purpose of the present study was to determine the nature of the substrate beneath the glacier near the middle of the ablation area, using geoelectrical techniques. To circumvent problems posed by the highly resistive ice cover, the measurements were made in a series of bore holes that penetrated to the base of the glacier (Fig. 1). This technique was pioneered by Haeberli and Reference Haeberli and FischFisch(1984). Radio-echo soundings (Reference BjörnssonBjörnsson, 1981) indicated that the ice was between 120 and 140 m thick in the area chosen.
The geoelectric method
Standard procedures for resistivity-sounding experiments are described in any geophysical prospecting textbook (e.g. Reference DobrinDobrin, 1976). Only a short overview will be given here.
In the resistivity method, either a direct current or an alternating current of very low frequency is applied to the ground by two electrodes, commonly referred to as the current electrodes. The potential conditions of the surface depend on the conductivity of the material beneath it. By measuring the potential difference between two additional electrodes, commonly referred to as the potential electrodes, a quantity known as the apparent resistivity, ρa, can be determined from ρa = KV/I where K is a geometric factor, defined later, V is the voltage difference, and I is the total current. ρa has the units of resistance times length and is most commonly given in ohm meters. By moving the current electrodes further apart. deeper layers of material gradually have a greater influence on the value of ρa. The actual resistivity of a layer, ρ, can be calculated from the apparent resistivity.
Procedure
Measurements conducted on the surface of a glacier would not reveal useful information on the subglacial bed characteristics because temperate or near-temperate ice is extremely resistant. Reference RöthlisbergerRöthlisberger (1967, p. 604), for example, reported values of ρ > 107 ohm m from several temperate glaciers, and King and others (in press) obtained values of 3–4 x 107 ohm m from measurements on Storglaciären. The ice thus shields the deeper material. To overcome this problem, a number of holes were drilled to the base of the glacier using hot-water drilling techniques (Reference Iken, Iken, Röthlisberger and HutterIken and others, 1977). By lowering electrodes down these holes, direct contact with the ice-bed interface was possible. Because the resistivity of glacial ice is so high, it can be assumed that the current entering it is negligible. The analysis can thus be conducted as if the geometry were an electrical half-space, similar to an air-rock interface.
To determine whether a particular bore hole reached the bottom, a simple electrical measurement was taken. The equipment used for this purpose, as well as in all other electrical measurements, was an ABEM Terrameter SAS 300B. A current was applied between the electrode being lowered down the hole and an electrode in the moraine ~300 m north of the study area (N’ in Fig. 1). Because glacier ice is an almost perfect insulator, the circuit was completed only when (and if) the electrode made contact with the subglacial material. (The fact that electrical contact was achieved only when the electrodes reached the bottoms of the holes, and not before, supports the assumption above that little current traveled through the ice itself.) Of 14 holes, ten were thus found to have reached bottom. Two of these refroze near the surface at times when electrodes were not in place. Electrodes were left permanently in place in the remaining eight, thus avoiding further problems with refreezing. Contact with the bottom was sporadic and unreliable in one of these (86–16; Fig. 1), however, and detailed inclinometry measurements in another (86–5) indicated that its bottom was displaced so far from the line of profile that it could not be used. Thus, of the holes drilled in 1986, only the six shown as solid lines in Figure 1 (inset B) were used in the experiment. A seventh hole (83–1; Fig, 1), drilled in 1983 and cased with aluminum tubing, was also used in one series of measurements (Fig. 2). The leads from the electrodes were gathered at a single point near the middle of the array so measurements could be made easily in a short period of time, even under adverse weather conditions.
Because the drilling was labor– and energy-intensive, the number of possible electrode placements was limited. This necessitated certain modifications of standard electrode arrangements. The half-Schlumberger arrangement (Reference Orellano and MooneyOrellana and Mooney, 1966) was found to be the most appropriate for the situation. In this arrangement, as applied herein, the potential electrodes were set in two of the bore holes, and one of the current electrodes was set in the moraine at N (Fig. 1), far enough away to be considered at infinity. The other current electrode was then moved to bore holes farther away from the potential electrodes. With this geometry, the geometrical factor K is given by
where MN is the spacing between the potential electrodes and L is the distance from the moving current electrode to the midpoint between the two potential electrodes.
The thickness of any conducting layer can be determined if the minimum value of L is of the same order as the layer thickness. The electrode spacing can then increase in a geometric progression (Fig. 1). As we expected that any till layer present would be thin, the potential electrodes were placed in holes that were 1.0 m apart at the surface, and the closest current electrode was in a hole 2 m from the hole containing the nearest potential electrode. Thus, L ∝ 2.5 m for this measurement. For measurements in the most distant hole (83–1)MN was increased to ~13m. Through rough inclinometry measurements and study of the internal consistency of the resistivity measurements, we established that the holes used were sufficiently straight that errors due to uncertainty in MN and L were negligible.
Each electrical measurement was averaged 16 times, first with the positive electrode at infinity and then with the negative electrode at infinity. These two values were then averaged. By means of this procedure, the influence of any natural telluric currents is cancelled. In each measurement we recorded a slight telluric current that was consistent in direction.
Before turning to the results, we point out that, although ice thicknesses determined by drilling were generally consistent with Reference BjörnssonBjörnsson’s (1981) radio-echo measurements, holes 86‒3 and 86‒4 suggested the presence of a 25 m deep gorge about 10 m south of hole 86–5. In view of the inclinometry results from hole 86–5, however, we suspect that holes 86–3 and 86‒4 may also have deviated significantly from vertical, and that this “gorge” is not real. Drilling techniques were modified after drilling hole 86–5, so subsequent holes were straighter
Discussion
Resistivity measurements were taken on 24 and 26 July 1986, and are plotted in Figure 2. The steep almost linear slopes of the lines in the plots indicate the presence of a thin layer of low-resistivity material, henceforth referred to as till, overlying a thicker layer with much higher resistivity. We assume that the latter is the amphibolite bedrock in the area. The resistivity of the amphibolite was measured in the pro-glacial environment and found to be ≈14 000 ohm m. As this is effectively infinite compared with the resistivity of the till, we can use the S-rule (Reference Orellano and MooneyOrellano and Mooney, 1966) to estimate the thickness of the latter: ~0.7 m. The slopes of the lines in Figure 2 are, however, somewhat steeper than is theoretically possible in a homogeneously layered medium. This suggests the presence of a lateral inhomogeneity of some kind. Southward thinning of the till or the presence of a topographic feature such as the “gorge” could have this effect. Analysis taking this into consideration (Reference FominaFomina, 1958; Reference MooneyMooney, 1958; Reference kumar and ChowdaryKumar and Chowdary, 1977) results in an estimated thickness of only ~0.35 m. There is considerable uncertainty involved in making such estimates based on electrical resistivity measurements; among other things, the interpretations involve a subjective assessment of what is realistically possible. Thus, while we feel that these estimates of the layer thickness are of the right order of magnitude, we do not claim significantly greater precision.
There is a marked difference in resistivity between the two series of measurements. On 24 July, ρ was 20 ohm m and on 26 July it was 3.5 ohm m. Between these 2 days a rainfall of 23.3 mm occurred. Because the electrodes were permanently emplaced, changing conditions at the bed are the only explanation that can account for the change in ρ. This indicates that the physical character of the low-resistivity layer changes rapidly in response to changes in water input. Increases in water input apparently increase the water content of the layer, thus decreasing the resistivity.
At the bottoms of two of the bore holes, Lindstrand (unpublished) collected samples of water containing suspended sediment. The material was predominantly silt with as much as 20% clay and less than 5‒10% sand. In bore hole 86‒12, a glob of compact clayey sediment was recovered on the bottom of the sampler. We infer that this was part of the bed itself. The composition of this material was found to be similar to that of the tills in front of the glacier. ρ was also measured on these pro-glacial tills and found to be ~500 ohm m. The much lower resistivity of the subglacial till layer suggests that its water content was substantially higher.
Conclusions
Knowledge of bed characteristics is important in any glaciological research, whether it be focused on movement, sliding, erosion, or water flow. Resistivity-sounding experiments can provide valuable information on the nature of this material. Our measurements on Storglaciären strongly suggest that there is a layer of till, ~0.4–0.7 m thick, beneath the experiment site. This layer probably extends to the moraine (N or N′ ) to the north and may extend as far as hole 83–1 to the south (Figs 1 and 2) but the resistivity data do not demonstrate this conclusively. The water content of this layer is high relative to the tills in front of the glacier, as indicated by the lower resistivities. The character of the layer adjusts rapidly to changes in water input.
Bore-hole deformation studies (Reference Hooke, Hooke, Holmlund and IversonHooke and others, 1987; paper in preparation by V. Pohjola and R. LeB. Hooke) suggest that the glacier is decoupled from the bed in late July and August at the site of the present experiment, and also at two additional sites nearer the glacier center line. We infer that the till layer is deformable and that this is responsible for the decoupling at the site of the present experiment. By extension, we further infer that a till layer is probably also present at the other two sites. Dye-trace studies (Hooke and others, in press) and bore-hole water-level measurements (unpublished data) indicate that water-flow velocities through this part of the glacier, an area where the bed is markedly overdeepened (Fig. 1), are unusually low. Continuous conduits of substantial size probably do not exist. This leads us to believe that a sub-glacial layer of deformable till is probably present beneath a large fraction of the upper half of the ablation area.
Acknowledgements
We owe a warm debt of thanks to the late H. Mooney, teacher and colleague, without whose invaluable advice during the planning stages, this project could never have been completed in its present form. V. Chandler reviewed the paper carefully; his critical comments and those of W.Haeberli resulted in significant improvements. A. Fleetwood and W. Horn assisted with instrumentation. O. Lindstrand provided unpublished data. Our co-workers at Tarfala provided encouragement and assistance during the field work.
Financial support was provided by the University of Minnesota Undergraduate Research Opportunities Program, the U.S. National Science Foundation (grant DPP–8414190), the Swedish Natural Sciences Research Council, the Swedish Society of Anthropology and Geography, and the University of Stockholm.