Introduction
The origin of glacier basal ice layers (B.I.L.) has held the attention of many authors. Several studies were made on a theoretical basis (Weertman,. 1961; Lliboutry 1964-65, p. 689). Direct observations were also carried out. For instance, at the base of the Blue Glacier (Washington, U.S.A.), crystallographic investigations supported by laboratory experiments led to the recognition of regelation phenomena (Reference Kamb and LaChapelleKamb and LaChapelle, 1964) Other field observations have confirmed this (Reference BoultonBoulton, 1970; Reference ClappertonClapperton, 1975; Reference VivianVivian, 1975) and have permitted a better understanding of ice formation at the glacier-bedrock interface, although in these studies the chemistry of basal ice was not investigated The present authors (Reference Souchez and LorrainSouchez and Lorrain, 1975 suggested that a knowledge of the chemical composition of basal ice would help in understanding its formation and clarifying the processes operating at this interface Approximately at the same time, the effect of subglacial chemical Presses on regelation phenomena and glacier sliding was studied theoretically by Hallet (1976[a], [b]). This author has put forward the role played by the phase relations of the subglacial Water-ice system.
In a recent paper (Hallet and others, in press), we and Hallet have shown a very good agreement between the chemical composition of basal ice from an Alpine glacier sliding over limestone and that predicted from his model.
In this paper, we present the results of chemical analyses of ice sampled at the base of several Alpine glaciers. In many of the cases studied, these results clearly indicate the nature of the processes by which basal ice forms.
Distribution and Characteristics of the B.I.L.
We have examined five examples of B.I.L. located under four Alpine glaciers all situated in Switzerland (Fig. 1). In each case, these B.I.L. are visibly different from the “normal” glacier ice which overlays them. In the one observed at the base of Glacier de Trient which is located on the Swiss side of the Mont Blanc massif, a homogeneous firn layer about I m thick appears at the base of the glacier front. It is very bubbly and whiter than glacier ice free of foliation, and composed of compressed snow partially transformed into ice. This feature is characteristic of very active advancing glaciers. Glacier de 1 rient has been advancing from 1958 until now (except during the 1960/61 period) at an average rate of 17.4 m/year. Avalanche snow is overridden by the glacier and progressively transformed into a frontal layer of basal firn or ice (Kasser and Aellen, 1974).
The other examples show that the B.I.L. is generally stratified. It is composed of a sequence of different types of layers with variable thickness. Three elementary layers constitute the basis for such sequences but are not always present together:
-
(a) a bubble-poor ice layer about 0.2 to 3 cm thick, nearly particle-free,
-
(b) a very dirty layer of about the same thickness containing more mineral particles than
-
(c) a very bubbly but particle-free ice layer, the thickness of which is very variable from one sequence to another.
It was possible to observe two different types of such sequences at the base of Glacier dc Tsanfleuron which is located, on limestones, north-west of the town of Sion in Can on Valais. At the base of the northern margin of the small tongue which terminates this plateau glacier in its north-eastern part, we studied a B.I.L. about 20 cm thick, composed of thin layers of bubble-poor ice (1 or 2 cm thick). It contains dispersed rock fragments up to 15 cm lone Along the southern margin of the same tongue, a completely different situation occurs In contrast with the northern margin, which consists of a sub-vertical ice ramp the southern margin is a relatively gentle slope and, from place to place, the glacier loses contact will, the bedrock. Some natural subglacial tunnels are thus accessible and allow the B.I.L. to be seen very clearly. It consists of a wedge of ice; its thickness reaching 1 or 2 m at the glacier margin diminishes gradually towards the interior, disappearing completely after about 8 to 10m The three elementary layers described above are present with a striking dominance of thick bubble-rich layers (up to 20 cm). These layers are separated by dirty layers locally associated with thin bubble-free ice (Fig. 2).
A somewhat different sequence has been encountered at the base of the right side of Glacier dc Tsijiore Nouve which is flowing over gneissic rocks, at the head of the Arolla branch of the Val d'Hérens (Valais). The site is located just below a large ice fall. At this point access to the B.I.L. is possible between the valley slope and the glacier margin because of à bend in the valley. The sequence of ice layering may be clearly seen in Figure 3. The bubble-rich layers are thick (up to 20 cm) and are separated from each other by a series of bubble-poor layers, each of them being overlain by a thin particle layer. The bubble-poor layers have a mean thickness of about 1 cm (maximum 2 cm) and the particle layers, of about 0.5 cm (maximum 1 cm). Each series of these two types of layers reaches a thickness of 10 to 15 cm. The whole sequence was observable over a height of about 2 m. Small rock fragments (up to 5 cm long) are dispersed in the layers.
Another sequence is to be found in a tunnel dug at the base of the Grubengletscher which is flowing over gneissic rocks in the Saastal. The tunnel is excavated partly in the ground moraine of the glacier and partly in the ice itself to channel the water of an ice-dammed lake. The tongue of the glacier is flowing in a permafrost zone as pointed out by Haeberli (1976). Temperature measurements made by this author indicate that the ice at this point is clearly below the pressure-melting point. In this last example, the sequence has a total thickness of about 1 m and is composed of bubble-poor layers about 2 cm thick interbedded with thin dirty layers and, from place to place, with an isolated thin bubbly layer. The dirty layers have a mean thickness of about 0.5 cm (maximum 1 cm). They often contain rock fragments which can reach twice the normal layer thickness or more so that they interrupt one or more bubble-poor ice layers (Fig. 4). The isolated bubbly layers appear in a proportion of about one per twenty bubble-poor layers and have a mean thickness of about 1 cm.
This description of different types of stratified B.I.L. shows that particle layers are associated with bubble-poor layers and that bubble-rich layers represent a variable fraction of the sequence, from about a half to a negligible proportion.
Methods of Sampling and Analysis
From the five B.I.L. described above, we sampled the four which are stratified and determined the composition of major cations in the ice by analysing the samples, in liquid form, by atomic absorption spectrophotometry.
The methods used to collect and treat the samples were devised to reduce contamination to a level comparable with the background noise level of the instrument and to control chemical changes caused by melting the samples. All samples were taken with the aid of a "Teflon"-coated USA-CRREL 3 inch (7.5 cm) ice corer. The cores were sliced along layer boundaries with a non-contaminant saw.
To determine the main cationic composition, two ice cores were taken at each sampling site. One was transferred in the frozen state to a press-filtered air-supplied laboratory; the other was treated in a field laboratory using a procedure we have described elsewhere (Hallet and others, in press). Samples having a particle content of more than 0.5 g/kg ice were systematically eliminated because of the potential importance of ion exchange and dissolution phenomena (Souchez and others, 1973). No significant differences were noticed between these two groups of samples.
210 Pb activity measurements were, for logistic reasons, limited to the south side of the tongue of the Glacier de Tsanfleuron. For these measurements, the chemical and low-level counting techniques were the same as those previously described by Crozaz and others (1964) and Crozaz (unpublished), and indeed the same apparatus was used. About 5 kg of ice were taken at each sampling site. Before melting, a Pb-carrier in concentrated nitric acid was added to the samples. Because of potential influence of melt water in a temperate glacier, no attempt was made to date the ice. Only a general comparison between the 210 Pb activity of the B.I.L. and of the ice immediately above was attempted. The validity of such an approach is discussed with the results below.
Cationic Composition of the B.I.L.
In our study with B. Hallet (Hallet and others, in press), we were able to explain quantitatively the calcium content of the B.I.L. from the north side of Glacier de Tsanfleuron from consideration of the phase relations in regelating subglacial solutions in approximate equilibrium with subglacial calcite deposits appearing on the lee side of rock protuberances. The calculated eutectic Ca concentration determined by the equation
(where the carbon dioxide partial pressure, pCQ2, is expressed in bars and the Ca concentration in eqI1), used in conjunction with an experimental distribution coefficient (0.01) which characterizes the solute partitioning during growth, can indeed yield estimates of the Ca content of the basal ice.
The minimum calcium content of the B.I.L. which would result if regelation is important can be easily estimated. If the subglacial water is at atmospheric pressure and in equilibrium with the carbon dioxide partial pressure of the atmosphere we obtain, using Equation ( 1) with 0.75 bar for the pressure (the altitude of the region is more than 2 500 m a.s.l.) and 0.03% for the carbon dioxide volume content, o.000 225 bar for pC02 and thus 0.001 3 eq l_I for the Ca content of the freezing subglacial water. Taking into account a distribution coefficient of 0.01, the Ca content of basal ice would be 13 μeq l"1. The value 0.01 for the distribution coefficient is rather low. In Hallet's experiments (Hallet, 1976[a]) and in our own unpublished ones, values of 0.02 were found, or even up to 0.05 if the rate of freezing is increased, and these will probably be more representative. This means that 13 $$$ μeq l"* Ca is the minimum value for basal ice if the B.I.L. is to be explained by regelation in equilibrium with calcite.
Several B.I.L. studied in the Alps do not reach such a Ca concentration. The B.I.L. of Glacier de Tsijiore Nouve described above is worth considering in this respect. As a result of weathering of silicate rocks by carbon dioxide and of limited occurrence of calcite in this environment, hydrocarbonate and calcium ions are present in the subglacial water, and therefore it is theoretically possible to reach by freezing the concentration level at which calcium carbonate will precipitate. Table I gives data for the cationic composition of three selected types of ice: glacier ice, and the bubbly and bubble-poor ice layers of the B.I.L. Values higher than 13 μeq I−1 Ca are not reached and thus the preceding model cannot apply. Table I also shows that the bubble-poor ice of the B.I.L. is distinctly different from the other types by a higher total cation content and a higher (Na+K)/{Ca+Mg) ratio.
Members of our laboratory (Lemmens and Roger, in press) have investigated the relationships between chemical concentration and discharge in the melt stream of the glacier frontal zone.
Good negative correlations exist if each cationic content is plotted against discharge on logarithmic scales: correlation coefficients are respectively —0.82 for Na, —0.69 for K, —0.95 for Ca, —0.95 for Mg, and are significant at a confidence level of 0.999. The ratio (Na+K)/(Ca+Mg) slightly decreases as concentration increases; its maximum value is about 0.10, which is not far from the values obtained on supraglacier streams (0.12). On the other hand, experiments were carried out at o°C in a cold laboratory leaving in contact for various periods the fine particles (<5θ μπι) of the morainic deposits from Glacier de Tsijiore Nouve and melted ice from the same glacier. The results indicate that the (Na+K)/(Ca + Mg) ratio decreases in the course of time, a value as low as 0.18 being recorded after only 1 min of contact. The two facts developed above indicate that water in contact with ice and with the bedrock in the area of Glacier de Tsijiore Nouve never reach a (Na + K)/(Ca+Mg) ratio as high as 0.58, which is the mean value for bubble-poor ice of the B.I.L. A possible reduction in the (Ca+Mg) content of the water by sub-glacial carbonate deposition is precluded because the calcium concentration in the ice is much lower than would be obtained by freezing water saturated with CaC03. In addition, there is no calcite deposition on the lee side of rock obstacles. Selective incorporation of cations during freezing will decrease the (Na+K)/( Ca + Mg) ratio in ice as compared to water since the divalent ions are preferentially incorporated (Malo and Baker, 1968). Simple refreczing of melt water circulating at the ice-bedrock interface cannot thus explain the higher (Na + K)/(Ca+Mg) ratio of bubble-poor ice of the B.I.L.
Robin (1976) discusses the pressure-melting mechanism within the basal ice mass and indicates that excess water is formed in zones of high-pressure ice, up-stream of obstacles. If this water is squeezed out of the ice by the pressure, a simple heat pump exists that will tend to cool the basal ice. This mechanism can produce cold patches at the ice-rock interface. Robin suggests that air bubbles would tend to be expelled along with the water being squeezed out of the capillaries, which partially explains the bubble-poor nature of the B.I.L. Robin's process probably occurs at the base of Glacier de Tsijiore Nouve but is not in itself sufficient to explain the peculiar composition of the bubble-poor ice layers of the B.I.L. Indeed, if water and air are expelled, the chemical content of the ice is likely to change because of selective flushing-out of ions. Thus a lower total ion content and a lower (Na+K)/ (Ca+Mg) ratio may be expected due to a more rapid migration of alkalis (which have larger diffusion coefficients than Ca or Mg), which is the reverse of the situation observed here. As a corollary, the squeezed water may be selectively enriched in alkalis by the same mechanism. On the other hand, forcing of squeezed water across a mud layer at the base may also produce an increase in (Na + K)/(Ca+Mg) ratio (Souciiez and others, 1973; Souchez and Lorrain, 1975). Thus, the peculiar chemical characteristics of bubble-poor ice from the B.I.L. may be explained by the refreezing of squeezed water circulating at the ice-bedrock interface. This process involves basal accretion of ice. The squeezed water does not have to be considerably enriched before freezing. Indeed, location of this B.I.L. along the margin of Glacier de Tsijiore Nouve indicates that cold air flowing beneath the glacier can explain the refreezing and thus the rate of freezing can be quite rapid so that a 5 to 10% distribution coefficient can apply (Gross, 1968).
The higher total ion content and the higher (Na+K)/(Ca-|-Mg) ratio of bubble-poor ice of the B.I.L. as compared with interbedded bubbly ice is not a situation unique to Glacier de Tsijiore Nouve. The same characteristics apply to Grubengletscher (table II), where refreezing is mainly due to the existence of permafrost beneath the glacier.
The bubbly ice of the B.I.L. at Glacier de Tsijiore Nouve and at the Grubengletscher has a chemical composition in major cations intermediate between glacier ice and bubble-poor ice layers. More information is clearly needed to explain this. In view of what will be developed in the next section, a wind-drift origin is possible. Such an origin would be in accordance with the greater thickness of bubbly ice in the B.I.L. at Glacier de Tsijiore Nouve which is developed along the marginal part of the glacier.
210 Pb Activity of the B.I.L.
This investigation was carried out on ice samples from a natural subglacial tunnel extending into the wedge-shaped B.I.L. and glacier ice of the south side of Glacier de Tsanfleuron. The :2l0Pb activity in dpm/kg of ice from the B.I.L. is compared in Table III with the activity of glacier ice situated immediately above. Percolation of melt water through small crevices in the ice can modify the ^"Pb activity by migration of this radionuclide along with the water. Therefore, considering temperate glacier ice as a closed system may be a rough approximation. However the internal consistency of the results in Table III probably signifies that this approximation is reasonable for the southern side of the Glacier de Tsanfleuron tongue. The filtration role played by the particles in the ice which tend to absorb migrating radionuclides (Prantl and others, 1973) must also be considered. This effect can be discounted in the present case since particles were carefully eliminated before sample melting had taken place, and since dirt extracted from a known volume of ice only shows negligible activity due to 210 Pb.
Table III shows that values of activity for glacier ice are well grouped around the mean, whereas for the wedge of basal ice below, two groups can be distinguished: low values of activity are encountered near the margin, probably because of more influence of melt water, and high values—higher than for the glacier ice above—are obtained from samples near the end of the subglacial tunnel. Melt water may have a highly variable 210 Pb activity depending on its origin (old glacier ice, recent snowbank) and on the relative mobility of this radio-nuclide during partial melting.
The fact that basal ice at the extremity of the wedge has substantially higher 2l0Pb activity values than glacier ice immediately above, substantiates an origin for that ice by basal accretion. Blowing snow deposited in a marginal crevasse under glacier ice is a normal feature of this margin of Glacier de Tsanfleuron. By compression induced by the overlying glacier ice, blowing snow, rock fragments, and regelating melt water may form several basal layers in the course of time. Such a B.I.L. is of course much more recent than the glacier ice above and this situation explains the contrast in 210 Pb activity. Some values encountered in the B.I.L. reach the 210 Pb activity of a fresh snowfall (mean value: 4 dpm/kg) indicating that the accretion process described is operating at the present time.
Conclusions
On the basis of location, stratigraphy, and chemistry, two origins can be suggested for the B.I.L. of Alpine glaciers:
-
(1) A basal accretion process occurs by refreezing of melt water squeezed out of the ice because of pressure fluctuations which change the pressure-melting point. Accretion by refreezing requires an open thermal system and is possible either because of the contribution of cold air even in the central part of the glacier {Vivian, 1975) or because of the presence of cold patches (Robin, 1976) or permafrost (Haerberli, 1976) in the down-glacier direction.
-
(2) A basal lateral accretion process occurs at the border of the glacier due to a supply of blowing snow and melt water from this snow at the interface where a décollement takes place.
In certain cases, the two processes may together contribute to the constitution of a single B.I.L.
Acknowledgements
A grant received from the Belgian Science Foundation for this study is gratefully acknowledged. Easy access to Glacier de Tsijiore Nouve was provided by Grande Dixence S.A. Thanks for this are due to Ing. Bezinge and Ing. Arlettaz. We are also indebted to Dr Barsch, Dr Haeberli, and Dr Rothlisberger for the possibility of working at the Grubengletscher sub-glacial tunnel and for fruitful discussions.