Introduction
Most isotopie studies dealing with high-altitude parts of Alpine glaciers have so far been based on shallow firn and ice cores which do not reach the bed (Reference Ambach, Eisner and PesslArnbach and others, 1972; Reference Jouzel, Merlivat and PourchetJouzel and others, 1977; Reference Oeschger, Schotterer, Stauffer, Haeberli and RöthlisbergerOeschger and others, 1977; Reference SchottererSchotterer and others, 1977, Reference Schotterer, Haeberli, Good, Oeschger and Röthlisberger1978; Stichler and others, 1982; Reference Haeberli, Schotterer, Wagenbach, Haeberli-Schwitter and BortenschlagerHaeberli and others, 1983; Reference Oerter, Baker, Stichler and RauertOerther and others, 1985) and, hence, can hardly give paleoclimatic information concerning time intervals reaching back more than a few years or decades. In this paper results of isotopie analyses in δD and δ18O are given for ice samples from a tunnel reaching the bed in the highest part of Titlisgletscher (Swiss Alps) at an altitude of 3020 m a.s.l. This exploratory sampling was carried out with particular emphasis on the very basal part. Additional values were obtained from a shallow core drilled from the surface above the tunnel. Questions are raised whether and how the results are interpretable from a paleoclimatic point of view.
Site Description
The investigated ice mass is the ice cover of the Chi i Titlis, a secondary summit of the Titlis near Engelberg, central Swiss Alps. The site is easily accessible by cable-car and the ice tunnel has been constructed by the cable-car company as a tourist attraction (Fig. 1). Glaciological studies, mainly carried out during consulting work for a governmental (PTT) radio tower on top of the mountain, included seismic refraction soundings on the geometry of the ice/rock-interface and the distribution of firn layers, thermal drilling, bore hole and tunnel temperature measurements, finite element modelling, and stake measurements near the bergschrund (Reference Haeberli, Iken and SiegenthalerHaeberli and others, 1979; Reference Haeberli and AleanHaeberli and Alean, 1985). The ice summit itself is a perennial comiche, consisting of well-layered ice roughly 25 m thick with a 15 m temperature of −0.7°C. Firn thickness increases downslope from 0 to about 15 m above the bergschrund. Subglacial permafrost temperatures, as measured in the rock tunnel connecting the cable-car station with the radio tower, increase steadily from about −1.5°C at the entrance on the north side to −0.5°C at the tower (rock wall exposed to the south, Fig. 2). Temperatures in the ice tunnel for tourists vary slightly around −1°C at the entrance, which is practically identical with the rock tunnel entrance, approaching melting point towards the bergschrund with increasing ice and firn overburden. Melt water percolation through crevasses in the bergschrund zone can be observed at the end of the ice tunnel. During construction of the rock tunnel, a rather gradual transition from debris-laden basal ice to ice-rich rock could be observed in places. Ice flow velocities at the surface are close to 1 m a×1. Sliding at the cold and dry bed above the bergschrund is assumed to be negligible.
Sampling
Samples were taken in the ice tunnel, near the entrance, where glacier thickness is a few metres only and a firn layer is absent. The ice contains about five air bubbles per cm3, some of them spherical with a diameter of 0.2-0.5 mm but most of them elongated, reaching lengths of about 15 mm. Fine mineral particles appear concentrated in sparse planes intersected by the tunnel walls but are more numerous and dispersed in the very basal layers.
Three types of samples were taken:
(1) A half cubic metre block of ice was cut from the very bottom of the tunnel wall about 15 m from the entrance at the only place where the bedrock was detected. From this block, different parts, a few cubic decimetres each, were cut in order to be stored at −20°C in portable freezers after being packed in plastic bags.
(2) On the opposite wall at a more easily accessible place, five stainless steel coring cylinders were hammered into the ice, regularly spaced from the bottom to the top of the wall, along the same vertical. Ice samples of 25 cm3 were transferred into fitted polyethylene bottles and also stored in a freezer. Before these two samplings, about 0.2 m of the tunnel wall was removed. Moreover, any contact with liquid water was carefully avoided.
(3) A 1 m-long ice core was taken with a CRREL auger from the glacier surface at a spot estimated to be above the second sampling site of the tunnel. The core was cut in three parts, put in plastic bags and stored in a freezer. All these original pieces of ice were resampled in a cold room to avoid contamination. At the present stage, 26 samples have been taken and analysed at the Centre d'Etudes Nucléaires de Saclay in France. The accuracy of the measurements is 0.5‰ for 8D and 0.17‰ for δ18O.
Analytical Results
The results of the isotopie analyses are presented on a graph in which δ18O is plotted against δD (Fig. 3). The 26 samples show a linear relationship with a slope S = 7.6 and a correlation coefficient of 0.98. This linear relationship is assimilable to the meteoric water line defined by Craig (Reference Craig1961). A slope value near 8 is considered in the literature as valid for ice which has not been submitted to major isotopie changes since its formation (Reference Dansgaard, Johnsen, Clausen and GundestrupDansgaard and others, 1973). A second graph (Fig. 4) presents the δD and δ18O values versus depth. This figure reveals a general trend for the δ-values to become less and less negative from the bottom towards the surface of the glacier with a shift of about 30‰ in δD and of about 3.5‰ in δ18O. From a closer examination of the graph, it can also be seen that a great part of the shift takes place within the deepest 20 cm, i.e. just above bedrock. In the lowest 10 cm of the ice, the δ-values indeed cluster around −115‰ in δD and −16‰ in δ18O whereas, just 10 cm above, they average −92.5‰ and −13‰ respectively; a particularly sharp gradient is evident.
Discussion
The increase of the δ-values from the botttom to the surface of the glacier as indicated above suggests a temperature rise since the formation of the deepest ice. Indeed, there is a direct relationship between temperature and δ-values. Theoretical considerations (Dansgaard, Reference Dansgaard1964) and data from the Mont Blanc Massif in the Alps (Reference Jouzel, Legrand, Pinglot, Pourchet and ReynaudJouzel and others, 1984) suggest that isotopie gradients of 4 to 6‰ per °C in δD and of 0.5 to 0.75‰ per °C in δ18O can be considered normal. Interpreting the isotopie shifts observed at Titlis according to these gradients implies a rise of 4.5 to 7°C in mean air temperature since the formation of the oldest ice sampled. Only rough estimates of minimum and maximum possible age of this transition can be derived from flow considerations. With a flowline length of 100–150 m perpendicular to the surface contours and, hence, along the marginal part of the ice, and an average surface velocity of 0.5 to 1 m a×1, perfect plasticity with n→∮ yields a lower bound value of 100 years; pure deformation with n = 1 and basal layer thickness = 1/20 to 1/50 the total ice thickness gives several thousand years. It therefore appears that the age of the oldest ice should be counted by centuries if not millennia, but is probably Holocene rather than (late-) Pleistocene. The climatic changes known for the period in question are comparable to 20th-century warming (0.5 to l°C, cf. Reference Gamper and SuterGamper and Suter, 1982; Haeberli, Reference Haeberli1990) and cannot account for temperature changes by several °C. Moreover, the fact that the measured δ-values correspond to the meteoric water line and, hence, show no evidence of basal melting and regelation along the flowpath further confirms that Neoglacial/Holocene variations of ice geometry and englacial temperatures must have been weak.
With increasing air temperature, however, percolation of melt water through snow and firn also increases and this may explain part of the observed isotopie shifts. There is indeed a change in the isotopie composition of melt water as it percolates through a snow-pack giving rise to a positive isotopie shift in the remaining snow. This process, which has been studied by Ȧrnason (1969) and Búason (1972) also explains the homogenization of the isotope content within the snow cover which is being transformed into glacier ice. Ȧrnason (Reference Ȧrnason1969, Reference Ȧrnason, Gat and Gonfiantini1981) compared two places in the firn area of Vatnajökull. In a place with abundant melt-water percolation, a positive shift of about 11.5‰ in δD was observed in the firn below the last summer surface whereas no shift was noticeable in a place where percolation was very limited. Ȧrnason concluded that the only way to explain the enrichment in the first case is to assume that the isotopie change occurs as the snow-pack recrystallizes. Moreover, he also proved by an experiment that the greater the proportion of the snow-pack which melts and produces percolating melt water, the greater is the positive shift mentioned above.
Other authors also have mentioned comparable isotopie shifts in the Alps. Arnbach and others (1972) studied firn cores in the temperate accumulation zone of Kesselwandferner. They observed the same type of isotopie shifts between the upper layers and the deeper ones, the latter being about 2‰ less negative in δ18O. The same authors also attributed a prominent role to melt water percolation in that situation, arguing that the value of the shift is lower in the highest part of the accumulation area than in the lower part. Moreover, Martinec and others (Reference Martinec, Moser, de Quervain, Rauert and Stichler1977) have studied a snow-pack with a lysimeter at Weissfluhjoch. They observed a shift of about 20‰ in δD within the lowest snow layer of the snow-pack between mid January when no percolation occurred and mid June when percolation was very effective. They noticed that this particular shift occurred without any change in the δD/δ18O-ratio, indicating that it cannot be due to evaporation effects.
Conclusions
In the ice of the high Alpine area of Titlisgletscher investigated here, isotopie shifts greater than those capable of explanation by direct temperature effect are observed. A substantial increase in melt water percolation since the formation of the deepest ice layer may help to explain the situation. The interpretation of isotopie shifts in terms of paleotemperature must thus be considered with caution in such environments. A quantitative treatment of the problem based on Búason's model (Reference Búason1972) would certainly be worthwhile but, unfortunately, cannot be carried out here because of the lack of information on the isotopie composition of present-day snow, rain water and melt water at the spot. Due to site modification by cableway construction, determination of representative flow and balance parameters is hardly possible today. The observed isotopie shift is nevertheless an interesting phenomenon and deserves further investigation at comparable sites.
Acknowledgements
Thanks are due to D. Gerber of the Engelberg-Titlis cable-car company for facilities given to the authors while sampling at the site, Dr J. Jouzel who made the isotopie analyses possible, Prof. R. Souchez who kindly made critical comments on the first draft of the paper, and Dr J.-L. Tison for his helpful participation to the ice sampling and for stimulating discussions.