Hostname: page-component-8448b6f56d-c4f8m Total loading time: 0 Render date: 2024-04-25T05:40:19.227Z Has data issue: false hasContentIssue false

Quantitative Determination of the Subglacial Hydrology of Two Alpine Glaciers

Published online by Cambridge University Press:  30 January 2017

David N. Collins*
Affiliation:
Department of Geography, University of Manchester, Manchester M13 9PL, England
Rights & Permissions [Opens in a new window]

Abstract

Two components of discharge through the internal hydrological systems of Alpine glaciers were separated on the basis of chemical composition of water. Some surface melt waters retain low solute contents after flowing without delay through conduits in which no chemical enrichment occurs, whereas those flowing slowly at the glacier bed have increased ionic concentrations. A simple mixing model was used to investigate temporal variations in the quantities of water routed through each of the two sub-systems. Electrical conductivity was taken as an indicator of melt-water composition and was monitored for periods during the summer ablation season of 1975 at Gornergletscher and of 1977 at Findelengletscher. At both glaciers, conductivity of melt waters varied diurnally inversely with discharge fluctuations, depending on the proportion of total discharge routed through the two sub-systems. Total discharge and the flow component routed rapidly through conduits within the glacier, a large proportion (50–80%) of total discharge, exhibited in-phase rhythmic diurnal hydrographs at the two glaciers. Distinctive subglacial hydrological regimes are contrasted. At Findelengletscher, the hydrographs of total discharge and of subglacial chemically enriched flow were in phase. At Gornergletscher, the subglacial hydrograph occurred with reverse asymmetry and out of phase. A possible interpretation is that water was temporarily stored in basal cavities during high total discharge. During the night, stored water was released, contributing much of the total discharge at times of low flow.

Résumé

Résumé

Les deux composantes du débit à travers le réseau hydrologique interne de glaciers alpins ont été distinguées sur la base de la composition chimique des eaux. Quelques eaux de fusion de surface n’ont qu’une faible teneur en matière dissoute après avoir coulé sans retard dans des cheneaux où il n’y a pas d’enrichissement chimique, tandis que celles qui ont coulé lentement le long du lit du glacier ont accru leur concentration en ions. Un simple modèle de mélange a servi pour rechercher les variations dans le temps des quantités d’eau empruntant l’un des deux sous-réseaux. La conductivité électrique fut choisie comme indicateur de la composition de l’eau de fusion et fut surveillée pour des périodes de la saison estivale d’ablation de 1975 au Gornergletscher et de 1977 au Findelengletscher. Sur les deux glaciers, la conductivité des eaux de fusion variait journellement en raison inverse des débits selon la proportion du débit total transitant dans les deux sous-réseaux. Le débit total et la composante de l’écoulement empruntant rapidement le réseau intra-glaciaire, une large part (50–80%) du débit total, montre pour les deux glaciers des hydrogrammes quotidiens rythmés en phase. Les régimes hydrologiques sous-glaciaires distincts sont opposés. Au Findelengletscher, les hydrogrammes du débit total et de l’écoulement sous-glaciaire, chimiquement enrichi étaient en phase. Au Gornergletscher, l’hydrogramme sous-glaciaire se produisit avec une asymétrie inverse et hors de phase. Une interprétation possible serait que l’eau fut temporairement stockée dans les cavités du fond pendant la période de fort débit total. Pendant la nuit, l’eau stockée était libérée contribuant au débit total au moment des basses eaux.

Zusammenfassung

Zusammenfassung

Auf der Grundlage der chemischen Zusammensetzung des Wassers liessen sich zwei Komponenten im Abfluss durch das innere hydrologische System von Alpengletschern unterscheiden. Ein Teil des oberflächlichen Schmelzwassers enthält geringe Lösungsanteile, wenn es ohne Verzögerung durch Kanäle geflossen ist, in denen keine chemische Anreicherung erfolgt; hingegen besitzt Wasser, das langsam über das Gletscherbett fliesst, eine erhöhte lonenkonzentration. Zur Untersuchung der zeitlichen Schwankungen in den Wassermengen, die durch jedes der beiden Teilsysteme strömen, wurde ein einfaches Mischmodell herangezogen. Als Indikator der Schmelzwasserzusammensetzung wurde die elektrische Leitfähigkeit benutzt und in Perioden während der sommerlichen Ablationszeit von 1975 am Gornergletscher und von 1977 am Findelengletscher aufgezeichnet. An beiden Gletschern schwankte die Leitfähigkeit der Schmelzwasser am Tag in umgekehrtem Rhythmus zu den Abflussschwankungen, jeweils in Abhängigkeit vom Anteil der Herkunft aus den beiden Teilsystemen am Gesamtabfluss. Der Gesamtabfluss und der Anteil, der schnell durch innere Kanäle des Gletschers strömt und mit 50–80% zum Gesamtabfluss beiträgt, zeigten phasengleiche rhythmische tägliche Pegelschwankungen an beiden Gletschern. Unterschiedliches Verhalten war jedoch für den subglazialen Abfluss zu beobachten. Am Findelengletscher waren die Aufzeichnungen für den Gesamtabfluss und für den subglazialen, chemisch angereicherten Abfluss phasengleich. Am Gornergletscher zeigte der subglaziale Abfluss umgekehrte Asymmetrie ohne Phasengleichheit. Dies lässt sich möglicherweise damit erklären, dass dort bei hohem Gesamtabfluss Wasser vorübergehend in Hohlräumen am Untergrund aufgestaut wurde. Während der Nacht wurde das aufgestaute Wasser freigegeben und bildete dann einen hohen Anteil am Gesamtabfluss, wenn dieser gering ist.

Type
Research Article
Copyright
Copyright © International Glaciological Society 1979

Introduction

Subglacial water is of importance in both the hydrology and dynamics of glaciers. A hydrological network at the bed of a glacier must allow movement of water produced by melting at places with high pressure to places of low pressure where refreezing occurs, and permit the transfer of large quantities of melt water down-glacier. Variations in the rate of sliding of a temperate glacier are probably affected markedly by the availability of subglacial water. The nature of a subglacial hydrological system remains problematical, and water films (Reference WeertmanWeertman, 1964), conduits (Reference NyeNye, 1973; Reference Behrens, Behrens, Bergmann, Moser, Ambach and JochumBehrens and others, 1975), and water-filled cavities (Reference LliboutryLliboutry, 1968) have been postulated. It is generally accepted that water flows in networks of englacial and subglacial conduits (Reference RöthlisbergerRöthlisberger, 1972; Reference ShreveShreve, 1972; Reference Krimmel, Krimmel, Tangborn and MeierKrimmel and others, 1973) though the location, structure, and functioning of such networks have not been settled.

Because the stream emerging from a glacier portal transports large quantities of sediment, flow within a glacier is inferred to occur at the bed in a network consisting of permanent channels incised into bedrock and temporary conduits incised upwards into the ice, which may suffer closure from meeting bedrock protuberances as the glacier moves (Reference NyeNye, 1973). Tunnels in ice tend to close under the ice overburden pressure but this tendency is opposed by melting of the walls by frictional heat produced by the flow of water, by melting if the water is initially warmer than ice, and because of water pressure in the tunnel if the cross-section is fully occupied. From calculations of theoretical water pressures, Reference RöthlisbergerRöthlisberger (1972) showed that water must flow in main arteries and that basal conduits can co-exist with lateral streams at the hydraulic grade line. The size of these channels probably adapts slowly to the quantity of melt water entering the system, such that at times of high discharge, there is insufficient time for capacity to adjust to supply and water pressure within the conduit network is increased.

The few direct observations of subglacial water channels have been made from adduction galleries designed to capture subglacial torrents for hydro-electric developments (Reference Vivian and ZumsteinVivian and Zumstein, 1973; Reference VivianVivian, 1977). These observations show that subglacial drainage is unstable, with stream courses continually changing positions. Indirect observations of water flow, using dye and salt tracers from individual moulins at the surface, to the glacier portal indicate rapid passage of melt waters across the ablation area (Reference StenborgStenborg, 1969; Reference Krimmel, Krimmel, Tangborn and MeierKrimmel and others, 1973), probably under open-channel flow conditions (Reference Behrens, Behrens, Bergmann, Moser, Ambach and JochumBehrens and others, 1975), and Reference Ambach, Ambach, Behrens, Bergmann and MoserAmbach and others (1972) considered that no large water-filled cavities existed at the bed of Hintereisferner, Ötztal Alps. Basal water-pressure measurements are difficult to interpret, since bore holes may or may not penetrate the subglacial conduit system (Reference HodgeHodge, 1976). From analysis of hydrographs of the river draining from Gornergletscher, Reference EllistonElliston (1973) considered that flow resulted from diurnal variations in water pressure of englacial reservoirs caused by variations in rates of water supply. Studies of the distinctive diurnal hydrographs of melt streams draining from glacier portals are inconclusive with respect to measurement of water flow at the bed.

This study aims to apply a method for the separation of the proportion of total discharge routed through the subglacial hydrological system to two Alpine valley glaciers. Water-quality characteristics are used as a means of determining the hydrographs of flow at the beds of the glaciers and understanding the functioning of their internal hydrological systems. Temporal observations were undertaken in order to evaluate changing subglacial flow capacities, and any interaction between conduit flow and temporary storage of water which may occur in cavities existing at the ice–rock interface (Reference LliboutryLliboutry, 1968). In non-glacierized catchments, portions of total discharge from delayed flow through the ground-water system and quick flow arising from precipitation have been determined from measurements of temporal changes of stream hydrochemistry (Reference Pinder and JonesPinder and Jones, 1969). In glacierized catchments, the chemical composition of ice and snow results from atmospheric influences, whereas run-off quality is determined also by terrestrial factors. This problem is complex because of the interaction in run-off of ice and snow melt, precipitation on the glacier surface and over the non-glacierized area, and a possible ground-water component.

Hydrochemical separation of components of flow within glaciers

Components of water flow through glaciers

During the summer ablation season, flow is principally derived from the melting of snow and ice at the glacier surface, since internal and subglacial melting together produce negligible quantities of water. Streams on the ice-free slopes of the catchment provide a significant contribution to total discharge only until all the snow has melted, and at times of heavy rainfall. At times with no precipitation in July and August, variations in the chemical quality of melt waters leaving the glacier portal reflect the interaction of water with sources of solutes within or beneath the ice, and any ground-water contribution.

Diurnal inverse variations of solute concentration of portal melt streams with discharge, during summer, result from dilution of solute-rich melt waters by relatively pure waters flowing rapidly through the glacier from the surface (Reference Rainwater and GuyRainwater and Guy, 1961; Reference CollinsCollins, [1978]). At least two routes for water flow through a glacier are suggested. Water derived from the surface can flow through the glacier without delay and without chemical enrichment. Since ice away from glacier soles is relatively solute-free (Reference RenaudRenaud, [1952]; Reference Hallet, Hallet, Lorrain and SouchezHallet and others, 1978), flow through an englacial, sediment-free, ice-walled conduit system would provide a route in which enrichment was minimized. Water routed through the subglacial zone, probably at overall slower rates, would become increasingly concentrated with ions on account of contact with the solute-rich basal hydrochemical environment (Reference Souchez, Souchez, Lorrain and LemmensSouchez and others, 1973; Reference Souchez and LorrainSouchez and Lorrain, 1975). Determinations of chemical composition have shown that subglacial waters in cavities beneath Alpine glaciers are relatively solute-rich (Reference Vivian and ZumsteinVivian and Zumstein, 1973). Observations in immediate pro-glacial areas from which temperate glaciers have recently retreated have demonstrated the occurrence of both calcite and silicate precipitation at glacier soles (Reference Ford, Ford, Fuller and DrakeFord and others, 1970; Reference HalletHallet, 1975), which reflect high solute concentrations in subglacial waters.

The subglacial hydrological system has channels which may be moraine-walled or cut in bedrock, and in which sediment transport occurs.

Solutes may be added to subglacial melt waters from bedrock, subglacial morainic sediments, sediment-laden basal ice and from suspended sediments in transit in streams. Rapid desorption, ion-exchange, and solution within the subglacial hydrochemical environment is assumed to bring subglacial waters to a uniform equilibrium concentration. The supply of ionic material is maintained by fracture of mineral lattices of basal particles and the bed by glacial erosion. Crushing and frictional wear of subglacial morainic layers mobilize ionic particles as the sediments are continually deformed. Most effective enrichment probably occurs where slowly moving dilute melt waters first enter sedimentary environments, especially in circulation through pores and capillaries in basal moraine. Subglacial spring water, considerably enriched chemically during slow seepage through bedrock cracks and pores, may add to the solute load of water at the bed.

Large quantities of englacially routed waters flowing from moulins and in arterial conduits located at the bed, under the ablation area, dilute waters circulating from the subglacial system. No enrichment occurs in arterial conduits because of large volumes and high rates of flow, and limited areas of contact with materials supplying ions. In this paper, “englacial” is used to describe the conduit system in which chemical enrichment does not occur. Englacial does not necessarily define conduit location, since arterial tunnels located at the bed form part of this system. The “subglacial" hydrological system is that in which enrichment takes place, and of necessity is located at the ice–bedrock interface but excludes basal arterial conduits.

Quantitative model

Concentration of solutes in total run-off results from the contributions of water and dissolved material from the two components of flow:

(1)

where Q represents the run-off proportions, C the total dissolved solids content or the concentration of an individual ion in solution, the subscripts refer to the total discharge from the glacier (t), flow routed through the subglacial system (s), and flow through the englacial system (e), and Q t = Q s+Q e.

Equation (1) can be solved for the portion of total discharge routed through the subglacial system:

(2)

Mass-balance Equation (2) can be used to determine subglacial flow (which separated on a chemical basis includes any ground-water flow) provided that both components of total run-off have a uniform density of 1 Mg m–3, no chemical reaction occurs when the components mix, and there is no reaction between suspended sediments from the subglacial system and englacially routed water when they meet. Sedimentary particles in suspension retain significant quantities of ions adsorbed to their surfaces, even in transit away from the portal (Reference Lorrain and SouchezLorrain and Souchez, 1972), suggesting that solution of particles (Reference SlattSlatt, 1972) and cation exchange between suspended load and melt waters (Reference CollinsCollins, [1978]) play only a minor role modifying the chemical composition of melt waters once flow is in large channels.

parameter estimation

In order to determine temporal variations in the subglacial flow of water, continuous records of both discharge and solute concentrations were required. Electrical conductivity was adopted as a measure of the total dissolved load of melt waters, since it can be monitored continuously with portable field instruments. Measured electrical conductivity is the sum of the conductances of the individual ions. Although conductivity may be used only as a measure of overall chemical composition, because the proportions of individual ions are not constant through time, inaccuracy resulting from this variability is considered to be negligible in the operation of this model. The conductivity (C e) of englacially routed water was taken as an average value of measurements of the input: ice and snow-melt waters on the glacier surface. Measurement of the electrical conductivity characterizing the hydrochemical environment beneath the glacier (C s) was not possible, and this parameter was obtained indirectly. During recession flow, when surface ablation ceases or is greatly reduced (for example, during the lying of snow from summer precipitation), the reduction of dilute melt-water input and the steady draining of englacially and subglacially stored water allows the solute concentration of total discharge from the portal to increase towards that pertaining in the subglacial system. The maximum conductivity repeatedly recorded during recession flows provides a minimum estimate of the true subglacial value, though some dilution probably occurs at all times during the ablation season.

Field measurements

Study area

The observations were made in the adjacent catchments of Findelengletscher and Gornergletscher, Canton Wallis, Switzerland (Fig. 1). Each glacier has only one pro-glacial meltwater stream emerging from its snout. The catchment of the Findelenbach extends over 24.9 km2, 76.7% of which is glacierized, and the stream is gauged 0.4 km from the snout of Findelengletscher. The Gornergletscher catchment has an area of 82 km2, which is currently 83.7% glacierized and the Gornera is gauged about 1 km from the snout. Characteristics of the catchment areas and glaciers are given in Table 1. Both catchments have considerable altitudinal range, from heights over 4 000 m a.s.l. to about 2 000 m a.s.l. Gornergletscher receives several tributary glaciers from southern snow-fields, including Unterer Theodul-gletscher, Zwillingsgletscher, and Grenzgletscher. Approximately 50% of the glacier surface of each of Findelengletscher and Gornergletscher lies in the ablation zone in summer (Reference Müller, Müller, Caflisch and MüllerMüller and others, 1976).

Fig. 1. Map showing the catchments of the Findelenbach and Gornera, and the sites of gauging stations and sampling sites.

Table I. Characteristics of the Findelengletscher and Gornergletscher catchments

Both catchments are underlain by rocks of igneous and metamorphic origins, granite, gneiss, schists, serpentine, amphibolite, and gabbro (Reference Lütschg, Lütschg, Huber, Huber and de QuervainLütschg and others, 1950), with a very small outcrop of calcareous rocks around Gornergrat (Reference BearthBearth, 1953). The bedrock topography of Gornergletscher is known in detail, but discrepancies between ice thicknesses determined by drilling and seismic methods suggest that a subglacial moraine layer may be up to 50 m thick at the centre of the glacier (Reference Bezinge, Bezinge, Perreten and SchaferBezinge and others, 1973).

Electrical conductivity

Records of conductivity of melt waters were obtained using immersed Sproule electrolytic cells of carbon electrodes in resin probes (cell constant = 1.0) attached to modified Walden Precision Apparatus CM 25 conductivity meters. The probes were positioned in the turbulent main flow lines of the Gornera and Findelenbach, at depths suitable to prevent the probe being left out of water during lowest discharges and to minimize the possibility of damage by bed-load movement. Probes were removed frequently to ensure that electrodes were intact, and that no sediment became lodged. Cells were calibrated against 0.01 M KCl at 25°C before and after use, and showed no temporal changes. Conductivity was monitored at a site 250 m from the glacier portal at Gornergletscher, and within 20 m of the ice margin at Findelengletscher.

Records were obtained for the Gornera during the ablation season of 1975, for a large surface melt stream on Gornergletscher during August 1976, and for the Findelenbach during August 1977. Several breaks resulted from recorder instrumental failures. Conductivity was measured in the presence of suspended sediment and was not corrected to a standard temperature, since individual melt-water samples were found to increase in conductance with temperature at different rates (Reference CollinsCollins, [1978]). Water temperatures ranged from 0.1 to 1.2°C in the Gornera, were consistently 0.1 °C in the supraglacial stream, and ranged from 0.1 to 4.0°C in the Findelenbach. The maximum errors introduced by the above-mentioned non-standardization are about 8% for the Gornera and 15% for the Findelenbach.

Samples of snow and ice melt waters were collected at several locations on the glacier surfaces. Conductivity was measured immediately on collection using platinum electrodes in glass dip cells (cell constants accurately calibrated in the range 0.10–1.65) attached to a WPA E1 conductance module.

Discharge

Both the Gornera and Findelenbach are gauged at the intakes of hydro-electric schemes. Between June and September, Grande Dixence, S.A. record continuous hydrographs of total discharges in flumes. Hourly mean discharges have been used in this paper, to remove short-term irregularities in flow.

Variations of electrical conductivity of melt waters

Ranges of variation

Ranges of variations of electrical conductivity of snow and ice melt waters and of outwash streams are shown in Table II for periods during summer ablation seasons. The results show very low conductivity values for all supraglacial streams irrespective of glacier and source. These dilute waters represent the atmospheric ionic input at the time of precipitation, as subsequently modified during firnification and regelation. The minimum conductivity recorded in the Gornera is almost within the range of the surface melt-water determinations, whereas for the Findelenbach the minimum is considerably higher. The maximum measured conductivity of the Findelenbach is considerably higher than that of the Gornera. Reference Lütschg, Lütschg, Huber, Huber and de QuervainLütschg and others (1950) considered that the ground-water component of winter discharge of the Findelenbach resulted in its having several individual ionic concentrations about ten times higher than those of melt streams draining neighbouring glaciers. The nature of the igneous and meta-morphic rocks would suggest that ground-water flow is unlikely in both catchments. Reference SlattSlatt (1972) found that bedrock type was not an important variable determining melt-stream hydrochemistry. In this study, only relative concentrations of melt waters within each stream are of importance.

Table II. Summary of measured electrical conductivity of melt waters

Diurnal variations

Discharge hydrographs and continuous records of electrical conductivity of the run-off in the Findelenbach are shown in Figure 2 for the period 1–24 August 1977. Discharge exhibits the usual distinctive diurnal rhythm of flow with the daily peaks superimposed on a slowly fluctuating background flow. The hydrographs show little short-term fluctuation and the only precipitation event to influence markedly the discharge occurred between 00.00 h and 01.00 h on 17 August, when 2.5 mm of rain fell over the catchment. Diurnal hydrographs are asymmetrical, showing a rapid rise to peak in response to the daily increase in energy supply for ablation followed by a slower decline. Maximum discharge occurred between 17.00 h and 19.00 h and minimum flow between 08.00 h and 10.00 h. Daily variations of conductivity of melt water occur in inverse phase compared with discharge. As flow increases, conductivity decreases rapidly as surface ablation waters pass quickly to the portal. During the night, solute concentration builds up slowly as the proportion of chemically-enriched subglacial water increases.

Fig. 2. Discharge hydrographs and chemographs of recorded electrical conductivity of melt water in the Findelenbach, draining from Findelengletscher from 1 to 24 August 1977, and precipitation at a gauge close to the glacier snout.

A similar inverse relationship between daily variations of discharge and electrical conductivity was observed in the Gornera (Fig. 3). The proportions of daily flows accounted for by the peaks were smaller than for the Findelenbach. Daily conductance minima often fell beneath 10.0 μS cm–1 (Fig. 3), suggesting that during the peak discharge periods much of the flow of the Gornera was derived from water passing through the englacial network without undergoing chemical change. Both the proportions and amounts of subglacial flow are reduced at times of high surface input to the glacier hydrological system.

Fig. 3. Discharge hydrographs and chemographs of recorded electrical conductivity of melt water in the Gornera, Gornergletscher, 20 July–2 September 1975.

Estimates of conductivity subglacial waters

The conductivities of subglacial water were taken as the highest observed conductivities during recession flows of the melt streams draining from the snouts. Records for the Gornera show that 44.0 μS cm–1 was a repeated maximum value of conductivity during several periods of snow-fall-induced recession in 1974 (Reference CollinsCollins, [1978]) and 1975. During observations at Findelengletscher, such recession occurred for only one period (21–24 August 1977). The maximum observed conductivity was 68.2 μS cm–1. A longer period of observation, extending into the winter recession flow, is needed to check the accuracy of these estimates as indicators of chemical conditions beneath the two glaciers.

Hydrographs of the components of discharge

Findelengletscher

The proportions of total discharge flowing through the englacial and subglacial hydro-logical networks, calculated from Equation (2), are shown in Figure 4. During the build-up of flow from 1 to 5 August 1977, the subglacial component of flow remained constant with a discharge of about 1.0 m3 s–1. The englacial conduits enabled rapid transport of ablation melt waters to the snout during the day, and the hydrographs of the flow through the englacial system are perfectly in phase with those of total flow, and are asymmetrical. Discharge from the englacial system continues throughout the night, suggesting that some water is delayed in flow, probably water arising from melting in the lower parts of the accumulation zone.

Fig. 4. Proportions of the total discharge of the Findelenbach routed through englacial and subglacial conduits of Findelengletscher calculated from measurements of electrical conductivity of melt waters and discharge.

Increased total discharge from 6 to 13 August was accompanied by a slight increase of mean daily flow in the subglacial conduits to 1.4 m3 s–1, with low-amplitude diurnal variation of subglacial flow, having maxima and minima in phase with those of the total discharge. A high subglacial flow on 11 August appears to have suddenly altered the subglacial conduit capacity. Subsequent incomplete data, while not showing evidence of diurnal rhythm, illustrate the increase of capacity of flow of subglacial conduits to about 2.2 m3s–1 as ablation was maintained. Reduction of energy supply during overcast days at the end of the observation period was marked by reduction of the flow through the subglacial system. The englacial system produced hydrographs of increased peak flows with considerable temporal variation, some flow at night and maintaining asymmetry. The asymmetry of flow probably results from the time–concentration curve resulting from the shape and extent of the glacier.

In early August, the proportion of englacial flow increased steadily to a maximum of about 60% of the total discharge. The subglacial system became increasingly important as the season progressed, and also carried most of the discharge when surface inputs were low. Following the high subglacial flow on 11 August, the proportion of water routed subglacially increased from 40 to 60%. Diurnal variations of proportions of water routed through the subglacial system are great, reflecting the importance of englacial conduits as routes for the transmission of much of the diurnal melt component through the ablation area.

Gornergletscher

The diurnal rhythmic component of discharge of the Gornera occurs in phase with variations of flow through the englacial conduits of Gornergletscher as at Findelengletscher (Fig. 5). Both total flow and the englacial component have asymmetrical hydrographs, with steep rising limbs concurrent with increasing melt at the surface in the morning. Rising-limb gradients are steeper for the hydrographs of the englacial component and recession of peak flow decreases more rapidly than the curve of total discharge of the Gornera. The proportion of flow in the englacial system occurring during diurnal peaks is greater than for total discharge. Some flow continues through the englacial system after ablation has ceased, probably supplied by ablation in the higher remote parts of the tributary glaciers, traversing long distances to the portal. During times of sustained ablation, 60–80% of the total flow was routed through the englacial conduits.

Fig. 5. Components of flow in englacial and subglacial conduits of Gornergletscher derived from measurements of discharge and electrical conductivity in the Gornera.

In contrast with Findelengletscher, as total flow through Gornergletscher varies diurnally, calculated discharge through the subglacial hydrological system shows an inverse hydrograph with reverse asymmetry. During afternoon and evening, when discharge of the Gornera is highest, both the quantity and proportion of flow through the subglacial conduits are much reduced. As total and englacial flows decline overnight, the discharge routed through the subglacial system increases approximately exponentially but fairly irregularly. The peak of subglacial discharge is rapidly terminated as flow through the englacial system builds up in the morning.

The diurnal cycle of subglacial flow was superimposed on a steady background basal discharge of 1.5–3.5 m3 s–1, which occurred throughout the period when ablation produced rhythmic diurnal variations of total flow in the Gornera. Maximum daily discharge in basal conduits varied from 3.2 to 8.3 m3 s–1 during sustained ablation. Flow in subglacial conduits may increase to two or three times its daily minimum during mid-morning, and subsequently decrease rapidly by about 80% in 7 h. The minimum percentage of total flow passing subglacially occurred at the same time as both minimum subglacial discharge and maximum total discharge, and vice versa, suggesting a delicate interdependence between englacial and subglacial components. The daily ranges of proportions of flow routed subglacially were greater than for Findelengletscher. During continuing ablation, between 10 and 25% of total flow always flows subglacially, but for a few hours each day 30–70% may be so routed. Total subglacial flow is variable from day to day, and this variability over-rides any trend which might result from the evolution of subglacial conduits during the ablation season. Day-by-day variations of total discharge and subglacial discharge are not released, implying glacial self-regulation of subglacial flow rather than hydrological or hydraulic pressure control.

The usual diurnal rhythm of flow of the Gornera was curtailed for three separate periods. On 15–16 August 1975, discharge remained at the early evening level throughout the night. The usual peak flow of the subglacial system was suppressed at this time, while increased flow through the englacial system produced a subsidiary peak between 00.00 h and 03.00 h, maintaining high flow until the onset of ablation in the morning. This event shows the control of subglacial flow by englacial flow within Gornergletscher. Since no precipitation occurred on these days, the increased englacial discharge probably resulted from the sudden draining of an englacial water body. Low temperatures and cloud cover from 8 to 13 August reduced the englacial contribution of total flow, but a break in the recording sequence has prevented the determination of the impact of precipitation. Diurnal flow regimes were maintained in both the en- and subglacial systems.

During the period 22–27 August, marked depletion occurred because of raised albedo and delayed run-off resulting from summer snow-fall. Reference EllistonElliston (1973) attributed the continuation of depletion flow for several days after snow-fall to the draining of englacial reservoirs. The flow of the englacial system effectively ceased after one day on 23 August, showing that water passes readily without storage through the conduit network. When some snow melt occurred on 24 August, the englacial conduits allowed a small volume of melt water quick access to the glacier portal. As englacial flow decreased, flow from the subglacial network increased for 16 h to a maximum of 6.4 m3 s–1 before subglacial storage became subject to depletion. During recession flow from 22 to 27 August, the subglacial component increased to 90% of the total.

From 28 August to the end of the observation period, slow recovery of flow through the englacial system resulted in low-amplitude daily peaks with an increasing background flow, suggesting flow retention in englacial conduits, either resulting from decreased capacity due to closure by ice-overburden pressure, or because low flows are associated with low velocities in open channel flow. Diurnal fluctuation of subglacial system discharge only re-commenced on 1 September, probably as depleted subglacial reservoirs re-filled.

Discussion

The simple two-component mixing model appears to give a useful impression of the contrasting behaviour of water beneath two Alpine glaciers, on the basis of measurements of chemical composition of melt waters. The method has the advantage of integrating the effects of water flow throughout the bed area, and provides an indication of average conditions at the bed, in contrast to isolated measurements of pressure in bore holes or moulins or dye tracing from individual moulins. Further, temporal variations of subglacial hydrological conditions can be determined, although no secular trend in flow was found at Gornergletscher. The separation of subglacial flow does not distinguish between groundwater and other components of basal discharge. Unfortunately, the model does not allow the estimation of total water storage at the base of a glacier but enables determination of rates of addition to and depletion from reservoirs.

Some inaccuracy in the calculated proportion of total flow routed through the subglacial system of a glacier will follow from errors resulting from the choice of conductivity values to represent subglacial hydrochemical environment conditions (C s) and englacial or surface melt-water characteristics (C e). If different estimates were used for these parameters, the percentage of total flow routed subglacially would be changed, with consequent effects on the actual quantities of melt waters calculated as passing through both the subglacial and englacial systems. In the calculations, the value of C s was taken as the maximum conductivity determined in the melt stream draining from each glacier during summer, and is a minimum estimate of C s. For Gornergletscher, should for example 88.0 μS cm–1 be the true value of conductivity associated with the subglacial environment (rather than 44.0 μS cm–1 used in the calculations), both the proportion of total flow routed subglacially (Q s/Q t) and the actual quantity Q s would be reduced by about 51% throughout the range of total discharges in comparison with the results presented here. The choice of value of C e also has an effect on the proportion of discharge passing subglacially and on Q s, reduction of which varies directly with Q t and inversely with C t. Taking C e as 4.0 μS cm–1, instead of 2.0, with C s = 44.0, the values of Q s given in this paper would be reduced by approximately 21% at highest total discharges, the percentage reduction in Q s declining towards zero as discharge decreases. Q e would be subject to corresponding numerical increases in respect of both parameter estimation errors since Q e = Q tQ s. Although there may be errors in estimation of parameters, the true shapes of the hydrographs of subglacial flow will be similar to those presented here (albeit for reduced absolute discharges and with reduced or accentuated amplitudes).

The results of separation of the subglacial component of flow at Findelengletscher provide evidence of seasonal changes in the nature of the internal drainage system of the glacier. The calculated subglacial discharges for the beginning of the observation period support the observations by Reference Behrens, Behrens, Bergmann, Moser, Rauert, Stichler, Ambach, Eisner and PesslBehrens and others (1971) that a steady basal component of flow is maintained beneath Alpine glaciers, with only slight diurnal variations. However, when total discharge from the glacier is high, the subglacial system acquires an in-phase diurnal rhythm of flow, alongside that of the englacial system, which responds rapidly to changes of input caused by surface hydrometeorological conditions. Interconnection between subglacial and englacial systems is probably widespread, allowing increased water pressures in the englacial hydro-logical system occurring at the times of highest discharges to cause in-phase variation of subglacial discharge. During the observation period, the subglacial drainage network was well developed, and water was able to pass through without storage in either cavities or subglacial moraine. Some cavities may occur at the bed but they are integrated with conduit water flow. Similar conditions were inferred from dye-tracing experiments at Hintereisferner (11.2 km2, length 6.9 km) (Reference Ambach, Ambach, Behrens, Bergmann and MoserAmbach and others, 1972; Reference Behrens, Behrens, Bergmann, Moser, Ambach and JochumBehrens and others, 1975).

The functioning of the internal hydrological system of Gornergletscher appears to depend on variations of total water supply from the surface. As water supply increases with ablation in the morning, flow through the englacial system rises rapidly. Not all supplied melt water can run off through the englacial conduits, which are probably adjusted to a mean discharge over several days (Reference RöthlisbergerRöthlisberger, 1972). Consequently, water pressure builds up in the conduits, when rates of channel widening fail to keep pace with increasing supply. Some increase in flow will occur on account of the pressure head. Subglacial water pressure also rises, assuming a network interconnecting englacial and subglacial systems. Subglacial water pressure will be relatively high at all times but it becomes increased at times of high water supply. Hydraulic pressure in the basal conduits may exceed the hydrostatic pressure of ice around the channel margins. The ice-pressure distribution over an irregular bed will not be uniform, and areas of reduced overburden pressure may occur down-stream of upstanding protuberances. Ice can separate from the bed down-stream of irregularities to form cavities (Reference LliboutryLliboutry, 1968).

It is suggested that water is diverted under pressure into cavities, which continue to expand provided water supply is maintained. Channel-margin storage in cavities can account for decreased conduit flow when water pressure is relatively high. The thick morainic layers beneath Gornergletscher may also store water which is forced at high pressure into inter-granular spaces. Although flow is reduced in the subglacial network, the actual quantity of water held at the bed is greatest during high water pressures. Evidence for diurnal fluctuations of water pressures at the bed has been obtained at Gornergletscher (Reference Bezinge, Bezinge, Perreten and SchaferBezinge and others, 1973) and at other glaciers (Reference Vivian and ZumsteinVivian and Zumstein, 1973; Reference HodgeHodge, 1976). When ablation decreases in late afternoon, the supply of surface melt waters is reduced, and falling subglacial water pressure is exceeded by the ice-overburden pressure. Water from the cavities and moraine, now chemically enriched, is returned slowly to the conduits.

The daily fluctuation of water storage beneath Gornergletscher suggests that about 30% of the total water flowing within the glacier spends at least 12 h in traversing through the basal conduits. The development of water-filled cavities at the bed may be related to high rates of glacier sliding (Reference LliboutryLliboutry, 1968), and the difference between water pressure in conduits and ice-overburden pressure, which also determines the closure rate of conduits (Reference RöthlisbergerRöthlisberger, 1972). High water pressures may result from the head of water within the thick glacier, or relate to the occurrence of closed depressions found in the long profile of the glacier (Reference BearthBearth, 1953). Water pressure appears to be transmitted from the englacial conduits to the glacier bed.

The importance of the non-enriched englacial components, as separated by the mixing model, supports Röthlisberger’s hypothesis of the existence of major arterial conduits. While these are probably located at the glacier bed, significant englacial cavities have been proved to exist within South Cascade Glacier (Reference HodgeHodge, 1976). The possible existence of drainage conduits away from the bed in the lower areas of ablation zones of Alpine glaciers would be of considerable practical importance in the subglacial collection of waters for hydro-electric development. Irregularity of the identified subglacial component suggests that the role of ground-water springs beneath glaciers is insignificant in summer, in contrast with late winter conditions previously reported (Reference Lütschg, Lütschg, Huber, Huber and de QuervainLütschg and others, 1950; Reference StenborgStenborg, 1965).

Conclusion

Separation of components of melt-water flow through Alpine glaciers on the basis of chemical composition has allowed some insight into the possible functioning of their internal hydrological systems. Although the parameters of the simple mass-balance equation may not be accurate, the two-component model of flow has enabled the identification of two contrasting basal hydrological regimes within Alpine valley glaciers. The model provides an overall approach to integrating the hydrological characteristics over wide areas of the glacier beds.

Some subglacial build-up of water in cavities occurs beneath some glaciers, but this phenomenon may not be general. Englacial and arterial conduits transmit significant proportions of total discharge irrespective of subglacial hydrological conditions. A better knowledge of subglacial hydrology may result from measurements of basal water pressures alongside observations of chemical composition of melt waters of several glaciers with differing sizes, thicknesses, and bedrock characteristics in order to determine the relationships between conduit flow and cavity storage.

Acknowledgements

The author wishes to thank A. Bezinge and J. P. Perreten of Grande Dixence, S.A., who have generously provided field support and made available discharge limnigraphs for the Findelenbach and Gornera.

Discussion

A. Iken: On Findelengletscher diurnal variations of horizontal surface velocity were observed in 1973. This seems to indicate that the glacier possessed a well interconnected subglacial system in the following sense: during the time of largest daily melt-water production, a high pressure builds up in the subglacial system and simultaneously in a dense network of cavities, thus affecting the sliding velocity. Some hours later the system drains quickly.

D. N. Collins: This observation is in agreement with the results I have described. A well-developed subglacial conduit system is indicated for Findelengletscher carrying a greater proportion of total discharge than that of Gornergletscher. In-phase discharge variations through both en- and subglacial systems suggests that cavity development is limited, but that there is an interconnected conduit network in which water pressure will increase at times of high discharge from the glacier portal.

C. C. Smart: It is possible to calculate the volume of water stored subglacially during the diurnal cycle you observe?

Collins: No. It is possible to make estimates only of the amount of water derived from storage in subglacial conduits and cavities; some water may remain in the cavities at all times and not be returned into conduit flow when conduit water pressure lowers.

T. Stenborg: Just a comment on your use of the terms “englacial" and “subglacial”. What you call englacial flow in fact comprises englacial and arterial subglacial flow, and your subglacial flow means the subglacial water outside at least the major subglacial arteries. Is that correct?

Collins: Englacial channels are envisaged as being those in which there is no chemical enrichment of melt water, whereas in subglacial conduits suspended sediments, morainic channel margins, and bedrock provide readily available solute supply sources. In the main arterial channels, I have assumed that no further chemical enrichment takes place when dilute waters from the ice-walled englacial channel system join. This is justified by the observations of melt waters appearing at the portal with solute characteristics similar to those of surface melt water, and that significant quantities of ions remain surface-adsorbed on sediment particles in suspension in pro-glacial streams. The main arterial channel can therefore be thought of as forming part of the englacial system, although it is located at the ice–rock interface. The englacial system effectively transmits water rapidly without solute enrichment, whereas flow in the non-arterial subglacial system is slower and involves chemical change.

References

Ambach, W., and others. 1972. Markierungsversuche am innern Abfluss-system des Hintereisferners (Ötztaler Alpen), Ambach, von W., Behrens, H., Bergmann, H., und Moser, H.. Zeitschrift für Gletscherkunde und Glazialgeologie, Bd. 8, Ht. 1–2, p. 13745.Google Scholar
Bearth, P. 1953. Erläuterungen. In Schwcizerische Geologische Kommission. Geologischer Atlas der Schweiz. 1: 25000. Bern, Kommissionsverlag Kümmerly und Frey AG, Blatt 535 Zermatt (Atlasblatt 29).)Google Scholar
Behrens, H., and others. 1971. Study of the discharge of Alpine glaciers by means of environmental isotopes and dye tracers, by Behrens, H., Bergmann, H., Moser, H., Rauert, W., Stichler, W. and Ambach, W., Eisner, H., Pessl, K.. zeitschrift für Gletscherkunde und Glazialgeologie, Bd. 7, Ht. 1–2, p. 79102.Google Scholar
Behrens, H. and others. 1975. On the water channels of the internal drainage system of the Hintereisferner, Ötztal Alps, Austria, by Behrens, H., Bergmann, H., Moser, H., Ambach, W., and Jochum, O.. Journal of Glaciology, Vol. 14, No. 72, p. 37582.CrossRefGoogle Scholar
Bezinge, A., and others. 1973. Phénomènes du lac glariaire du Gorner, Bezinge, par A., Perreten, J. P., et Schafer, F.. Union Géodésique et Géophysique Internationale. Association Internationale d’Hydrologie Scientifique. Commission de Neiges et Glaces. Symposium on the Hydrology of Glaciers, Cambridge, 7–13 September 1969, p. 6578. (Publication No. 95 de l’Association Internationale d’Hydrologie Scientifique.)Google Scholar
Collins, D. N. [1978.] Hydrology of an Alpine glacier as indicated by the chemical composition of meltwater. Zeitschrift für Gletscherkunde und Glazialgeologie, Bd. 13, Ht. 1–2, 1977, p. 21938.Google Scholar
Elliston, G. R. 1973. Water movement through the Gornergletscher. Union Géodésique et Géophysique Internationale. Association Internationale d’Hydrologie Scientifique. Commission de Neiges et Glaces. Symposium on the Hydrology of Glaciers, Cambridge, 7–13 September 1969, p. 7984. (Publication No. 95 de l’Association Internationale d’Hydrologie Scientifique.)Google Scholar
Ford, D. C, and others. 1970. Calcite precipitates at the soles of temperate glaciers, [by] Ford, D. C., Fuller, P. G., Drake, J. J., Nature, Vol. 226, No. 5244, p. 44142.Google Scholar
Hallet, B. 1975. Subglacial silica deposits. Nature, Vol. 254, No. 5502, p. 68283.CrossRefGoogle Scholar
Hallet, B., and others. 1978. The composition of basal ice from a glacier sliding over limestones, [by] Hallet, B., Lorrain, R. [D.], and Souchez, R. [A.]. Geological Society of America. Bulletin, Vol. 89, No. 2, p. 31420.Google Scholar
Hodge, S. M. 1976. Direct measurement of basal water pressures: a pilot study. Journal of Glaciology, Vol. 16, No. 74, p. 20518.CrossRefGoogle Scholar
Krimmel, R.M., and others. 1973. Water flow through a temperate glacier, [by] Krimmel, R. M., Tangborn, W. V., and Meier, M. F.. (In [International Hydrological Decade.] The role of snow and ice in hydrology. Proceedings of the Banff symposia, September 1975. Paris, UNESCO; Geneva, WMO; Budapest, IAHS, Vol. 1, p. 40116.)Google Scholar
Lliboutry, L. A. 1968. General theory of subglacial cavitation and sliding of temperate glaciers. Journal of Glaciology, Vol. 7, No. 49, p. 2158.CrossRefGoogle Scholar
Lorrain, R. D., and Souchez, R. A, 1972. Sorption as a factor in the transport of major cations by meltwaters from an Alpine glacier. Quaternary Research, Vol. 2, No. 2, p. 25356.CrossRefGoogle Scholar
Lütschg, O., and others. 1950. Zum Wasserhaushalt des schweizer Hochgcbirges. 9. Kapitel. Zur Hydrologie, Chemie, und Geologie der winterlichen Gletscherabflüsse der schweizer Alpen, Lütschg, von O., Huber, P., Huber, H., und de Quervain, M. [R.]. Beiträdge zur Geologie der Schweiz. Geotechnische Serie. Hydrologie, Lief. 4.Google Scholar
Müller, F., and others. 1976. Firn und Eis der schweizer Alpen. [By] Müller, F., Caflisch, T., Müller, G.. Eidg. Technische Hochschule Zürich. Geographisches Institut. Publikation Nr. 57.Google Scholar
Nye, J. F. 1973. Water at the bed of a glacier. Union Géodésique et Géophysique Internationale. Association Internationale d’Hydrologie Scientifique. Commission de Neiges et Glaces. Symposium on the Hydrology of Glaciers, Cambridge, 7–13 September 1969, p. 18994. (Publication No. 95 de l’Association Internationale d’Hydrologie Scientifique.)Google Scholar
Pinder, G. F., and Jones, J. F. 1969. Determination of the groundwater component of peak discharge from the chemistry of total runoff. Water Resources Research, Vol. 5, No. 2, p. 43845.Google Scholar
Rainwater, F. H., and Guy, H. P. 1961. Some observations on the hydrochemistry and sedimentation of the Chamberlin Glacier area, Alaska. U.S. Geological Survey. Professional Paper 414-C.CrossRefGoogle Scholar
Renaud, A. [1952.] Nouvelle contribution à l’étude du grain de glacier. Union Géodésique et Géophysique Internationale. Association Internationale d’Hydrologie Scientifique. Assemblee générale de Bruxelles, 1951, Tom 1, p. 20611. (Publication No. 32 de l’Association Internationale d’Hydrologie.)Google Scholar
Röthlisberger, H. 1972. Water pressure in intra- and subglacial channels. Journal of Glaciology, Vol. 11, No. 62, p. 177203.CrossRefGoogle Scholar
Shreve, R. L. 1972. Movement of waters in glaciers. Journal of Glaciology, Vol. 11, No. 62, p. 20514.Google Scholar
Slatt, R. M. 1972. Geochemistry of meltwater streams from nine Alaskan glaciers. Geological Society of America. Bulletin, Vol. 83, No. 4, p. 112531.Google Scholar
Souchez, R. A., and Lorrain, R. D, 1975. Chemical sorting effect at the base of an Alpine glacier. Journal of Glaciology, Vol. 14, No. 71, p. 26165.CrossRefGoogle Scholar
Souchez, R. A., and others. 1973. Refreezing of interstitial water in a subglacial cavity of an Alpine glacier as indicated by the chemical composition of ice, by Souchez, R. A. Lorrain, R. D. and Lemmens, M. M. Journal of Glaciology, Vol. 12, No. 66, p. 45359.CrossRefGoogle Scholar
Stenborg, T. 1965. Problems concerning winter run-off from glaciers. Geografiska Annaler, Vol. 47A, No. 3, p. 14184.CrossRefGoogle Scholar
Stenborg, T. 1969. Studies of the internal drainage of glaciers. Geografiska Annaler, Vol. 51A, Nos 1–2, p. 1341.Google Scholar
Vivian, R. A. 1977. Tourism, summer ski-ing, hydroelectricity, and protection of the public in the French Alpine glacial area: the development of an applied glaciology. Journal of Glaciology, Vol. 19, No. 81, p. 63942.CrossRefGoogle Scholar
Vivian, R. A. and Zumstein, J. 1973. Hydrologie sous-glaciaire au glacier d’Argentière (Mont-Blanc, France). Union Géodésique et Géophysique Internationale. Association Internationale d’Hydrologie Scientifique. Commission de Neiges et Glaces. Symposium on the Hydrology of Glaciers, Cambridge, 7–13 September 1969, p. 5364. (Publication No. 95 de l’Association Internationale d’Hydrologie Scientifique.)Google Scholar
Weertman, J. 1964. The theory of glacier sliding. Journal of Glaciology, Vol. 5, No. 39, p. 287303.Google Scholar
Figure 0

Fig. 1. Map showing the catchments of the Findelenbach and Gornera, and the sites of gauging stations and sampling sites.

Figure 1

Table I. Characteristics of the Findelengletscher and Gornergletscher catchments

Figure 2

Table II. Summary of measured electrical conductivity of melt waters

Figure 3

Fig. 2. Discharge hydrographs and chemographs of recorded electrical conductivity of melt water in the Findelenbach, draining from Findelengletscher from 1 to 24 August 1977, and precipitation at a gauge close to the glacier snout.

Figure 4

Fig. 3. Discharge hydrographs and chemographs of recorded electrical conductivity of melt water in the Gornera, Gornergletscher, 20 July–2 September 1975.

Figure 5

Fig. 4. Proportions of the total discharge of the Findelenbach routed through englacial and subglacial conduits of Findelengletscher calculated from measurements of electrical conductivity of melt waters and discharge.

Figure 6

Fig. 5. Components of flow in englacial and subglacial conduits of Gornergletscher derived from measurements of discharge and electrical conductivity in the Gornera.