Introduction
Thermal aspects are fundamental to understanding the formation, existence, preservation, and stability of ice in the ground. Rich information and correspondingly specialized literature is available concerning permafrost in polar lowlands (e.g.Reference Gold and Lachenbruch Gold and Lachenbruch, 1973; Reference Lachenbruch, Cladouhos, Saltus and SennesetLachenbruch and others, 1988), where increasing attention is now being paid to effects of ongoing secular warming trends (Reference OsterkampOsterkamp, 1984;Reference Smith Smith, 1988). In contrast, much less is known about the thermal conditions of mountain permafrost at lower latitudes, despite the fact that ice-rich permafrost at a high elevation is common. Creep of such permafrost on slopes is widespread and easily recognizable from the large number of striking, lava-stream-like rock glaciers (e.g.Reference Wahrhaftig and Cox Wahrhaftig and Cox, 1959;Reference Barsch Barsch, 1978; Reference HaeberliHaeberli, 1985; Reference Giardino and VitekGiardino and Vitek, 1988).
Within the framework of a recently started core-drilling project, two bore holes were drilled into the permafrost of the active rock glacier Murtèl 1 (Fig. 1) near Piz Corvatsch, Engadin/Grisons, eastern Swiss Alps (Reference Haeberli, Huder, Keusen, Pika and RöthlisbergerHaeberli and others, 1988). The present paper describes the bore-hole temperatures measured during the first year of observation, discusses the thermal properties of the frozen material, and presents an analysis of the heat-flow conditions within the creeping rock-glacier permafrost. A complete description of the measurements has been given by Reference Vonder MöhllVonder Muhll (unpublished). Long-term monitoring of bore-hole temperatures is planned in order to improve the thermal analysis and to investigate possible effects of anticipated warming trends with respect to the stability of mountain permafrost on steep slopes.
Investigated Site
The active rock glacier Murtèl I is situated close to the lower boundary of the local permafrost distribution and, in fact, presently advances over permafrost-free terrain (Reference HÖlzleHölzle, unpublished). A number of measurements and soundings had previously been carried out at this easily accessible site (Reference BarschBarsch, 1973, 1977; Reference Barsch and HellBarsch and Hell, 1975; Reference King, Fisch, Haeberli and WaechterKing and others, 1987). Surface altitude at the drill site is 2670 m a.s.l. and exposure is to the north-west. The active-layer thickness on the rock glacier is 2–4 m and perennially frozen, ice-rich sediments reach the slightly overdeepened bedrock. Surface movements are on the order of a few centimeters per year, reflecting compressive flow within the lower part of the rock-glacier lobe, where striking ogive-like transverse ridges exist.
In April, May, and June 1987, an experimental hole (bore hole I, 21 m) and a deep hole (bore hole II, 62.5 m) — separated by a horizontal distance of about 2 m — were drilled into the rock-glacier permafrost. To minimize thermal disturbance of the cores and the bore hole during drilling operations, a triple-tube coring system was used in combination with cold-air cooling (cf.Reference Lange Lange, 1973). The last piece of ice was found at 51 m depth, heavily fissured and highly permeable, non-frozen bedrock was reached at about 52 m, impermeable rocks were encountered near 54 m, and the drilling stopped at 62.5 m. A set of bore-hole logs was obtained and a plastic tube with magnetic rings was installed down to a depth of 58 m, where loose material had obstructed the bore hole, in order to monitor horizontal as well as vertical bore-hole deformation. Both bore holes remained air-filled (no bore-hole fluid) and were protected with an impermeable lid against atmospheric influences and snow avalanches. Figure 2 illustrates preliminary results of core and bore-hole observations which are helpful in the geothermal analysis presented below (cf. Reference Haeberli, Gamper, Zimmermann and KienholzHaeberli and others, 1988, 1989): an upper layer extending from the permafrost table down to 28 m consists of massive ice and supersaturated sands and silts with an extremely high ice content (>90% ice by volume); a lower, saturated layer from 32 to 50 m mainly contains coarse granodiorite blocks with ice-filled pores; 75% (4 cm year−1) of the total horizontal movement at the surface (6 cm year−1) takes place in super-saturated frozen fine material of the transition zone between these two layers, and the lower (structured) layer does not measurably deform. Vertical compression of about 1 cm year−1 in a layer immediately above the main shear horizon enables roughly surface-parallel flow to take place in the upper part of the rock-glacier permafrost.
Bore-Hole Temperatures
Bore-hole temperatures were measured with two types of negative-temperature-coefficient (NTC) thermistors: Fenwall UUA 41J1 with a resistance of approximately 3.27 χ 104 Ω at 0°C, and Yellow Springs Instrument YSI 44006 with a resistance of about 2.95 χ 104 Ω at 0°C; both thermistors have a temperature coefficient of about 5% per 0C Thermistors were placed in a metallic cylinder for protection against water and mechanical influences. High-precision calibration of the instruments was done by the Swiss Federal Office of Metrology at Wabern, Bern. The accuracy of thermistor readings in the laboratory is better than ±0.01°C. Readings under field conditions using a Fluke 8060A digital multimeter gave absolute values assumed to be correct within the nearest tenth and relative differences measured with the same thermistor within the nearest hundredth of a degree Centigrade.
A wire equipped with thermistors at intervals of I m was placed in the open bore hole I down to 21 m depth. In the deep bore hole II, thermistors, again at 1 m intervals, were mounted on the outside of the plastic tube used for strain measurements, which could be lowered into the bore hole down to 58 m depth. In the lower part of this hole, from the bottom up to 30 m below the surface, the space between the plastic tube and the bore hole wall was filled with concrete, whereas in the upper part between 30 m and the permafrost table water was used, which quickly turned into ice. Within the concrete filling all thermistors developed an instability and were rendered unreliable. The cause is not exactly known, but it may be related to mechanical damage to the cable insulation by a magnetic ring during installation of the plastic tube or to subsequent chemical interaction with the concrete. To obtain temperatures periodically at depths below 30 m, a chain with thermistors at intervals of 2 m was lowered into the plastic tube down to its bottom. After a few hours of thermal adjustment to the surrounding temperature as checked with thermistors which are firmly installed above 30 m, temperatures were measured and the chain was pulled up 1 m five times at appropriate time intervals enabling two to three thermistor readings to be taken at every meter. Convection within the air-filled bore holes is assumed to be negligible (cf.Reference Haeberli and Funk Haeberli and Funk, in press). After instrumentation (1 July 1987), measurements were made daily during 1 week. Measurement intervals were extended continuously until August 1987 and readings were taken every 3 weeks thereafter. Temperatures below 30 m were measured four times a year.
In bore hole I, temperatures were determined between 13 July 1987 and 30 July 1988. During this period, the two bore holes showed similar temperatures. Besides transient effects from drilling disturbance — which lasted less than 6d in the case of the open bore hole I as verified by bottom-hole temperature measurements on 14–15 May 1987 — temperature differences between the two bore holes were less than 0.1°C. The following considerations apply to the data from bore hole II. Measurement procedures in the this deep bore hole began 2d after completion of the borehole instrumentation. The main thermal disturbance from drilling activity resulted from the infilling of concrete. Comparison with the development of temperatures in the much less disturbed ice-filled section above 30 m and below the depth of seasonal fluctuations at about 20 m indicates that the dissipation of warming due to the concrete infilling agrees well with the theoretical concepts proposed byReference Lachenbruch and Brewer Lachenbruch and Brewer (1959): the thermal disturbance was reduced to a few hundredths of a degree after about half a year (Fig. 3).
Figure 4 illustrates the temperature of the rock glacier as measured during the first year after thermal stabilization. Temperature varies seasonally down to a depth of about 20 m (= depth of Zero Annual Amplitude, ZAA); it remains negative throughout the year below the permafrost table at about 3 m and steadily increases with depth from the ZAA down to the permafrost base at 52 m. The following features of the profile deserve special attention:
1. Within the uppermost 15m, mean annual temperatures are not aligned with the relatively straight temperature profile at greater depth, but curve towards warmer temperatures. Above about 10 m the temperature gradient is even inverted, probably reflecting a marked but most recent development towards warmer surface tmperatures. “Equilibrium surface temperature” as estimated by extrapolating the linear otemperature profile at depth to the surface is around −3°C, while mean surface temperature measured in 1987–88 is close to −1°. Similar thermal offsets have been described in soil temperature profiles from northern Canada and the U.S.S.R. (Reference Burn and SmithBurn and Smith, 1987).
2. Between the ZAA and a depth of 30 m the temperature gradient is 69 ± 11 °C km−1, reducing to 47±5°CkirT1 between 30 and 50 m. This break in the temperature profile at around 30 m depth will be discussed below together with the results from heat-flow calculations.
3. The deepest fragment of ice was detected at 51 m depth. The temperature measured there is −0.08 °C, i.e. only about twice the value (-0.035 °C) expected from pressure-induced depression of the freezing point. Hence, no evidence is found in the rock-glacier permafrost for strong chemical or soil-particle effects on freezing-point depression as sometimes observed in polar permafrost (cf. Reference Osterkamp and PayneOsterkamp and Payne, 1981; Reference Collett and BirdCollet and Bird, 1988).
4. Between 52 and 55 m, i.e. at the permafrost base and within the fissured top zone of bedrock, seasonal fluctuations at slightly positive temperatures can be observed. This phenomenon indeed hints at the possible presence of circulating ground water fed from the surface by seasonal snow and ice melt, possibly also from circulating air.
5. At the bottom of the hole, temperature is again slightly negative (-0.01 °C). This inversion of the temperature gradient between 55 and 58 m depth remains fairly stable in time, leaving the possibility open that permafrost conditions exist again at greater depth, i.e. beneath the permeable bedrock zone, and that the bore hole penetrates a non-frozen zone (talik) between 52 and 55 m.
Figure 5 shows measured temperature variations above the ZAA. At the permafrost table — where the mean temperature during the first measurement year was about −1.3°C high-frequency changes show up and the seasonal temperature variation is strongly asymmetrical. The “zero curtain” due to effects of latent-heat exchange during melting (spring) and freezing (autumn) of the active layer in 1988 has approximately the same duration (cf. Reference Lachenbruch, Brewer, Greene and MarshallLachenbruch and others, 1962). Because of the insulating effect of the thick snow cover, winter temperatures remain nearly constant at −4°C as expected from earlier snow-temperature measurements for rapid permafrost mapping (cf.Reference King, Fisch, Haeberli and Waechter King and others, 1987;Reference HÖlzle Hölzle, unpublished), whereas summer temperatures fluctuate more rapidly and with higher amplitudes. At greater depths, high-frequency variations are filtered out, temperature variations become progressively more symmetrical, and below about 9–10 m closely resemble sine curves. Due probably to transient effects from special weather conditions during the first measurement year, however, the period of these sine curves was about 10 months and the temperature minimum in 1988 was slightly warmer than that in 1987, although the latter was still influenced by thermal disturbance from the drilling activity. The seasonal temperature maximum propagating downwards was shifted by about half a year at 10 m depth with respect to the surface. The amplitude A of the oscillation was attenuated with depth; following the relation
where m represents the slope of the regression line in Figure 5 (correlation coefficient = 0.9978).
Thermal Properties
Thermal conductivity K and thermal diffusivity к are
the two main parameters characterizing the thermal properties of the investigated permafrost. Due to the influence of unfrozen water near the melting temperature and corresponding problems with sample disturbance from core recovery and storage, exact and representative values cannot easily be determined. In order to get realistic estimates for the investigated bore hole, the problem was approached in two ways: by laboratory measurements on frozen material and by theoretical analysis of the seasonal temperature variations observed in the bore hole as described above.
Thermal conductivity K was measured in a cold laboratory using a Quick Thermal conductivity Meter (QTM) composed of a flat-iron-like thermal probe coupled with an electronic processor. The measuring principle is based on a line source on a semi-infinite medium. Only the probe had to be in the cool-room. During measurement, the heater wire’s temperature increased by about 10–15°C. To prevent the sample from reaching the melting point and to avoid corresponding effects from latent-heat exchange, the cold laboratory was set at a temperature of −18°C. The thermal conductivity of two rocks (granodiorite and amphibolite) from the rock-glacier surface, an ice sample, and two ice-saturated samples containing silt, sand, and gravel was determined. The ice and ice/rock samples were prepared in the laboratory using ordinary fresh water and fine material collected at the rock-glacier front. The number of randomly oriented measurements on each sample varies between 28 and 54 (Table I).
Assuming that the rock glacier’s permafrost is composed of about 50% amphibolite and 50% granodiorite rocks, the mean thermal conductivity for the rock mixtures in question can be estimated at about 2.5 W m−1 °C−1, a value which does not significantly differ from the value of ice. Reference Lachenbruch, Sass, Marshall and Moses jrLachenbruch and others (1982) proposed an equation which expresses the thermal conductivity K of an aggregate in terms of the thermal conductivities of its (n) constituents k1 K2 …. Kn and their respective volume fractions Φ1,, Φ2, &, Φn
Since thermal conductivities of rocks and ice are similar, K =1 Kl = K2. However, measured values of thermal conductivity for supersaturated debris mixtures (3.0 ± 0.3 W m 1 C−1) are higher than both, ice as well as rock values. The reason for this discrepancy is not yet clear, but it may be related to increased quartz concentrations in the weathered fine material from the rock-glacier front.
In the absence of sources or sinks of latent heat and of convective heat transport through percolating melt water, temperature T at depth : as a function of time t in a supersaturated permafrost with a presumably low permeability is controlled by heat conduction following the equation (Reference Carslaw and JaegerCarslaw and Jaeger, 1959;Reference Tautz Tautz, 1971):
To calculate thermal diffusivity к from the observed seasonal variations of bore-hole temperatures by analyzing the attenuation of temperature amplitudes and the phase lag of the corresponding sine functions, Fourier series were fitted to three measured temperature-time curves at depths of 6.6, 9.6, and 12.6 m. The following equations were used:
Attenuation of amplitude: Phase lag:
Phase lag:
where ω = 2π/р, p is the period, h represents the vertical separation between the two thermistors corresponding to temperature curves with amplitudes A1 (upper) and A2 (lower), and δ means the phase lag between the two sine functions.
For the uppermost 20 m, the thermal diffusivity к can be determined from the slope m of the straight line in Figure 5, assuming that the period p is 1 year:
Thermal conductivity K, thermal diffusivity к, density p, and mass specific heat c are related by (Reference Gold and LachenbruchGold and Lachenbruch, 1973):
Considering that the product pc for ice (1.10−3°C−1) and granitic rocks (1.30Jm−3°C−1; cf.Reference Carslaw and Jaeger Carslaw and Jaeger, 1959) does not differ much and based on the fact that the ice content is 90–100% by volume within the considered depth range, thermal conductivity K can be calculated from thermal diffusivity к:
Table II is a summary of calculated thermal conductivities in Wm−1° C−1 for several layers in a range, where ice content amounts to about 95% (measured £ for pure ice-2.5 ± 0.2 W m−1 0C−1.
The scatter of the calculated values is considerable (roughly ±50%). The best agreement with the laboratory measurements is reached with the values for the uppermost 20 m as a whole and for 6.6–9.6 m depth. All theoretically calculated values may be influenced by transient effects due to ongoing trends in the development of surface temperatures.
Heat Flow
Besides surface temperature and thermal properties — and disregarding the complex processes (phase transitions, convection, etc.) within the non-saturated active layer — a number of factors influence the temperature distribution and heat flow within the creeping permafrost of an active rock glacier (Fig. 6): (1) geothermal heat from the underlying bedrock, (2) advective heat from sub-permafrost ground water or air, (3) viscous dissipation from deformation and sliding, (4) advective heat from horizontal displacements and vertical compression, (5) “sedimentation effects” due to the build-up of the rock-glacier body in time (rising surface), and (6) paleoclimatic fluctuations. In addition, three-dimensional effects of the extreme surface roughness with its thermal implications (wind-blown crests and snow-filled troughs) could influence the thermal conditions observed in the bore hole. In view of this complexity, the following one-dimensional consideration of heat flow along the bore hole cannot be more than a first attempt to describe the relative importance of some of the processes involved.
Conductive heat flow q along the bore hole is the product of the temperature gradient and the thermal conductivity. Thermal conductivity must be assumed to be constant with depth, since it is identical for ice and encountered rocks within the above-described uncertainties of its determination. Hence, variations in calculated values of q are essentially controlled by variations in the measured temperature gradient. In the zones with annual temperature variations above the ZAA and beneath the permafrost base, temperature gradients strongly change during the year. In both cases, however, temperature values averaged over the entire observation period indicate negative q values, i.e. heat flow is predominantly towards deeper layers. As mentioned before, these heat-flow inversions point to a most recent warming at the surface and a possible continuation of permafrost beneath the bottom of the hole, i.e. within the underlying bedrock. The latter feature could involve transient effects from secular warming and interference of ground water or air.
The calculated heat flow between the ZAA and the permafrost base at 52 m varies between about 100 and 220 mWm−2 with an overall average of around 150 mWm−2. Mean geothermal heat flow in Switzerland is 85 mW пГ with peak values occurring in relation to ascending thermal water (Schinznach Bad and Oberuzwil: 150mWnT2) and, hence, being due to convective-flow systems (Reference Bodmer and RybachBodmer and Rybach, 1985). In comparison, the value observed in the Murtèl bore hole is remarkably high. The geothermal map of Switzerland (Reference Bodmer and RybachBodmer and Rybach, 1984) gives only two measurements for the eastern Grisons, suggesting a regional geothermal heat flow at Murtèl I of about 60–110 mWm2 or one- to two-thirds of the measured average at Murtèl. An important local heat source must therefore exist at the base of the permafrost. The seasonal temperature variations at this depth indeed point to the possible influence of advective heat from relatively rapidly circulating sub-permafrost ground water or air.
At around 30 m depth, temperature gradient changes quite sharply, indicating an apparent increase in heat flow from a value of 120 ± 30 mW −2> between 30 and 50 m to one of 175 ± 15 mW ra−2 between 30 m and the ZAA. Interpretation of this feature involves several possibilities but remains highly uncertain. Undetected variations in thermal conductivity are the most probable reason. However, assuming a simplified two-layer model with pure ice (K = 2.5 Wm−1 ଁC−1, q = 170 mW m−2) on top and q = constant with depth, thermal conductivity of the deeper layer would have to be about 3.6 W m−1° °C−1. This value is considerably higher than the described laboratory values. It would involve a marked difference between the K values of ice and rocks, respectively, possibly corresponding to a change in rock-quartz content with depth of about 35%. Other explanations may even be less plausible. Viscous dissipation from permafrost deformation within the main shear zone represents a heat source which can be calculated from the shear stress and the measured deformational velocity at the top of the shear zone (Reference ZotikovZotikov, 1986) but remains by three orders of magnitude below the required value. The influence of advection from up-stream parts of the creeping permafrost can be estimated using the Peclet number P representing the ratio between characteristic time constants for vertical heat diffusion (T1) and horizontal heat advection (TJL cf.Reference Barsch Blatter, 1987); T1 is given by L2 1/K and T2 by L2/u, where L1 is the characteristic distance for vertical heat diffusion (30 m in the present case), k is the thermal diffusivity, L2 is the characteristic distance for horizontal heat advection, and u is the characteristic horizontal flow velocity (5–10 cm year−1 in the present case). P is close to 1 if L2 is less than a few meters — an order of magnitude which renders strong advective effects rather implausible. Effects of three-dimensional surface topography seem improbable in view of the quite sharp change in dT/dz, and the same is true for paleoclimatic effects which would have to involve accelerated cooling during the past few years. The observed vertical compression immediately above the main shear zone would also be able to enhance the vertical heat flow in the upper layer, but it takes place at a considerable distance from the main break in dt/dz. These brief considerations illustrate the relative importance of some processes influencing the heat-flow conditions within an active rock glacier.
Discussion and Conclusions
The most detailed observations reported so far with regard to thermal conditions within creeping mountain permafrost are those by Reference Johnson and NicklingJohnson and Nickling (1979) from a partially vegetation-covered active rock glacier near Kluane Lake, Yukon Territory, Canada, and the shallow-depth measurements byReference Haeberli Haeberli (1985) from the active Gruben rock glacier in the Swiss Alps. The bore-hole measurements down to 17 m depth near Kluane Lake indicated rapid warming and degradation of permafrost, coupled with thermokarstic development, whereas the average surface temperature at Gruben remained close to −1 °C over several years. Drilling through the active rock glacier Murtèl I now promotes a more complete discussion of temperature, thickness, and stability of rock-glacier permafrost.
Average temperature for recent years at the permafrost table of the rock glacier is estimated at about −3°C. Considering the fact that the twentieth century is among the warmest periods within post-glacial time, permafrost conditions have probably existed during the whole of the Holocene and, hence, during the entire time of rock-glacier formation in the Alps (cf. more detailed discussion by Reference HaeberliHaeberli (1985)).
The thermally defined base of the permafrost at 52 m depth is very close to the base of the ice-rich permafrost as given by the debris/bedrock interface near 50 m depth. It can easily be estimated in this case (cf.Reference Osterkamp Osterkamp, 1984;Reference Haeberli Haeberli, 1985) that melting at the permafrost base as a consequence of surface temperature suddenly rising to 0°C would take place after about a decade only. Complete melting of the permafrost by geothermal heat alone would then take several thousand years. With surface temperatures rising even further, i.e. above 0°C, melting from the surface downwards could be faster, but permafrost degradation would even then take many decades if not centuries to be completed. These orders of magnitude clearly show that rock-glacier permafrost has a long “memory” and that its thermal state must be explained with due consideration to past developments.
In a number of Arctic permafrost bore holes, thermal anomalies have developed as a consequence of recent temperature changes, with reduced or even inverted heat flow above about 150 m depth due to warming in the early twentieth century and with normal to increased heat flow within the uppermost 50 m due to cooling since around 1950 (Reference Lachenbruch, Cladouhos, Saltus and SennesetLachenbruch and others, 1988). Permafrost in the Alps is thought to have warmed up as well following the end of the Little Ice Age around 1850; its thickness and temperature profile may therefore not be in equilibrium with present-day conditions at the surface (Reference HaeberliHaeberli, 1985). Permafrost thickness at Murtèl I could easily have been IOOm or more during the Little Ice Age, if one assumes average geothermal heat flow, no thermal influence from ground water, and^ a mean surface temperature colder than today by about 1°C. The possibility that relict permafrost from the Little Ice Age indeed still exists below the depth reached by the bore hole — as indicated by the inversion of the heat flow near the bottom of the hole — is therefore not unreasonable. Slight cooling of surface temperatures during the past decades may have increased the heat flow above 50 m depth. At some time in the past, ground water or air circulating within the heavily fissured top layer of the bedrock may have started to form a talik with a sharp break in the temperature profile and, hence, to increase even further the heat flow in the permafrost above (cf.Reference Echelmeyer Echelmeyer, 1987). Long-term monitoring of bore-hole temperatures is continuing to permit time-dependent modelling of such scenarios and documentation of future trends.
With regard to the fact that the front of rock glacier Murtèl I is advancing over permafrost-free terrain, the mean annual surface temperature of about −3°C at the drill site, about 200 m behind the front, is remarkably low. Moreover, it is very close to the mean annual air temperature as calculated from weather stations, despite the thick winter snow cover usually existing at this site. Rock-glacier permafrost obviously reflects an especially cold micro-climatic environment and extreme gradients of surface-temperature conditions must exist at the rock-glacier front. As neither air temperature nor incoming solar radiation changes dramatically between the drill site and the terrain in front of the rock glacier, surface characteristics must be assumed to play an important role. The prevailing absence of systematic studies about active-layer processes on rock glaciers, however, renders all explanations speculative. It is possible that the permafrost of rock glaciers creates its own micro-climate by various processes such as heat loss through water run-off and, hence, that critical threshold values could exist with regard to its stability. Such important aspects remain to be investigated.
Acknowledgements
The study was funded by a special grant from ETH Zurich (permafrost core drilling, project 0230.70). The help of P. Gnos (Vaw) and a number of colleagues during the preparation and installation of the thermistors is gratefully acknowledged. Thanks are also due to Professor S. Müller and Professor L. Rybach from the Geophysical Institute οΓ ETH Zürich, as well as to Dr A. Iken and Dr M. Funk (VAW) for critical comments on the manuscript and helpful suggestions. B. Nedela assisted with the drawings and S. Braun edited the English of the present paper.