Hostname: page-component-7c8c6479df-5xszh Total loading time: 0 Render date: 2024-03-19T01:14:03.726Z Has data issue: false hasContentIssue false

Estimate of Glacier Ablation under a Debris Layer from Surface Temperature and Meteorological Variables

Published online by Cambridge University Press:  20 January 2017

M. Nakawo
Affiliation:
Geotechnical Section, Division of Building Research, National Research Council of Canada, Ottawa, Ontario K1A0R6, Canada
G.J. Young
Affiliation:
Snow and Ice Division, National Hydrology Research Institute, Inland Waters Directorate, Environment Canada, Ottawa, Ontario K1A0E7, Canada
Rights & Permissions [Opens in a new window]

Abstract

A simple model suggests that the ablation under a debris layer could be estimated from meteorological variables if the surface temperature data of the layer are available. This method was tested by analyzing the data obtained from experiments with artificial debris layers. Fairly good agreement was obtained between the estimated and the experimental data.

Résumé

Résumé

Evaluation de l’ablution d’un glacier sous une couche détritique à partir de la température de surface et de variables météorologiques. Il est possible, grâce à un modèle simple, d’évaluer l’ablation d’un glacier sous une couche de roches détritiques en fonction de variables météorologiques si l’on possède des données sur la température de surface de la couche détritique. Ce modèle a été vérifié lors de l’analyse des données provenant d’essais effectués au moyen de couches détritiques artificielles. Les données expérimentales sont venues corroborer les données obtenues à l’aide des modèles théoriques.

Zusammenfassung

Zusammenfassung

Abschätzung der Eisablation unter einer Schuttschicht mit Hilfe der Oberflächentemperatur und von meteorologischen Daten. Aus einem einfachen Modell geht hervor, dass die Ablation unter einer Schuttschicht aus meteorologischen Daten abgeschätzt werden kann, wenn die Oberflächentemperatur der Schicht bekannt ist. Diese Methode wurde mit Hilfe von Daten erprobt, die aus Versuchen mit kūnstlichen Schuttschichten hervorgingen. Zwischen den abgeschätzten und den experimentell gewonnenen Werten ergab sich eine recht gute Übereinstimmung.

Type
Research Article
Copyright
Copyright © International Glaciological Society 1982

Introduction

In order to evaluate glacier ablation under a debris layer. Reference Nakawo and YoungNakawo and Young (1981) proposed a simple model which was successfully employed in analyzing experimental data. With this model, ablation under a debris layer can be estimated from meteorological variables when the thermal resistance of the layer is known. Since it is difficult to determine directly the thermal resistance of a layer of unknown material in the field, it was suggested that the surface temperature of the debris layer may be used for estimating the thermal resistance and consequently the ablation under the layer.

This paper presents the results of testing the validity of the proposed method by comparing estimated data with field measurements. The symbols used are denned in Table 1.

Table I. Nomenclature

Model

The energy-balance equation at a debris surface, in which all the terms are taken to be positive downward, is given by

(1)

where

(2)
(3)
(4)

Assuming a steady temperature profile in the debris layer (i.e. a linear profile for a uniform layer), then

(5)

since the temperature at the ice-debris interface is 0°C and T s is in degrees Celsius. Neglecting the variation of the stored heat in the layer, and assuming no conduction of heat into the ice beneath, then

(6)

When condensation takes place, it is assumed that e s is equal to the saturation vapour pressure, which is a function of T s. As long as the debris surface is wet, this assumption is also made for periods when evaporation occurs. For a dry surface, on the other hand, e s is assumed to be equal to ea .

By combining Equations (1) through (6) and eliminating T s (and e s with the above assumptions), one can estimate r for a given R when F (or G and A if α is known), u a, T a, p, and e a are provided. This was demonstrated by Reference Nakawo and YoungNakawo and Young (1981). In most cases in the field, however, the value of R is unknown.

When T s is given instead, F, H, and E can be estimated (Equations (2) through (4)), allowing R to be determined by combining Equations (1) and (5). Once R is determined, r can be estimated for other periods using the procedure mentioned above. This is the method to be tested.

Experimental data

Experiments were carried out at Peyto Glacier (lat. 51° 41′N., long. 116° 33′W.) in the Rocky Mountains, Alberta, Canada from 20 to 22 August 1979. Meteorological variables during the measurement period are summarized in Table II. The data were collected using the procedures reported by Reference Nakawo and YoungNakawo and Young (1981).

Table II. Meteorological variables during the experimental period

The ablation rate under debris layers was observed at six plots prepared artificially with debris materials collected from the supraglacial debris of the glacier. Each plot was 0.3 m square with a layer thickness h given in Table III. The ablation at the plots during a given period was determined by measuring the increase in the relative distance between the debris surface and a taut string installed over the plots. The results are also compiled in Table III.

Table III. Meteorological variables during the experimental period

Surface temperature was measured by a thermistor inserted within a few millimetres of the surface of the debris. This measurement was made only twice in the daytime, but this is considered to be satisfactory for a test of the method as the weather was very stable during the experiments (e.g. atmospheric pressure was almost constant at 803.1 ± 1 mbar). The observed values for T s of the measurements are shown in Table IV.

Table IV. Values for the meteorological variables when the surface temperature was measured

Test results and discussion

The albedo of the layers was not determined. The debris material had a relatively dark colour and its albedo, when dry, was considered to be about 0.1 to 0.2 (Reference PenndorfPenndorf, 1956; Reference GeigerGeiger, 1961). For a wet surface, the albedo decreases by about 20% (Reference GeigerGeiger, 1961), but it was still considered to be in the range of 0.1 to 0.2.

By substituting the values of T s (Table III) and meteorological variables (Table IV) into Equations(2) through (5), and combining with Equation (1), R was estimated for each plot. The estimated values are shown in Table III. Uncertainty in R was caused both by the uncertainty in the albedo and by the difference in the two estimates for 21 and 22 August. It should be noted that E in Equations (1) and (4) was assumed to be zero in the calculation for plot F because its surface was dry. Thermal conductivity Km estimated from these R values was in the range 1.4 to 2.6 W m −1 deg −1 m (Table III) which are comparable with the values obtained for various soils (e.g. Reference KerstenKersten, 1949; Reference Penner, Penner, Johnston and GoodrichPenner and others, 1975; Reference JumikisJumikis, 1977).

The ablation rate was calculated using the data on meteorological variables given in Table II, and assuming α = 0.1 and α = 0.2. The results are plotted against R in Figure 1 (a) for the first two periods and (b) for the latter two periods. The solid and dashed lines are for wet and dry surfaces respectively in the daytime. The short dashed lines correspond to the estimates for night time. The observed ablation rates in Table III are also plotted in Figure 1 using the R values given in Table III.

Fig.1. Ablation rate versus thermal resistance during 20–21 August(a)and 21–22 August(b).Solid and open circles are observed data during night time and day time respectively. Short dashed lines show the estimation from meteoro logical variables for night–time. Solid and dashed lines represent the calculation for wet and dry surface respectively during daytime.

The agreement between calculation and observation is fairly good, although there are some discrepancies. It is considered that the disagreement could be attributed to uncertainty in the estimates of R because there were few measurements of T s and there was the uncertainty as to whether the temperature was in a steady state at the time of the observations. The errors involved in the measurement of T s could also cause an error in the determination of R, particularly when R is large. For plots E and F, for example, a 0.5 deg difference in T s would result in 10 × 10 3 and 5 × 10−3 deg W−1 difference respectively in the value of R when the modified Ts is applied through Equations (1) to (5).

Another source of disagreement between the calculated and the observed data is the uncertainty in the value of /β. The value of 4.89 J m−3 deg−1 is an average compiled by Reference NaruseNaruse and others (1970) from the data for β obtained at various surfaces of glaciers, snow fields, and artificial basins. The original data for β scattered in a range of ± 1.16 J m−3 deg−1 around the mean value. Owing to the presence of the experimental plots, the surface roughness of these plots on Peyto Glacier was greater than that of a natural glacier surface. This would result in a larger value of β than for a natural surface. The value of β at the plot could therefore have been larger than 4.89 J m−3 deg −1. Using a larger value of β would result in a larger ablation rate for a given R; if R is large, however, an increase in the value of β has little effect. The value of β is also dependent on wind stratification. Log–linear profiles of wind speed and temperature were found applicable at the glacier (Reference DerikxDerikx, [1975]; Reference Munro and DaviesMunro and Davies, 1977, Reference Munro and Davies1978). In the present experiments advection could have played an important role in heat exchange at the surfaces of the plots, since the area of the plots was small. However, the determination of the value of β taking the advection term into consideration is a very complex problem.

Nonetheless, the general agreement between the calculated and observed values suggests that glacier ablation under a debris layer can be predicted from meteorological and surface temperature measurements. To obtain a good prediction, it is recommended, as pointed out by Reference KrausKraus ([1975 ]), that special attention be paid to surface roughness which is sometimes very large at stagnant areas near termini of glaciers (see e.g. Reference Iwata, Iwata, Watanabe and FushimiIwata and others, 1980). A continuous record of surface temperature as well as observations on temperature profile in the debris layer would also improve the prediction.

Acknowledgements

The authors would like to thank N. Maeno and H. Narita of the Institute of Low Temperature Science, Hokkaido University, for their assistance in the field, and L. W. Gold of the Division of Building Research, National Research Council of Canada, for encouraging them to write this paper.

References

Derikx, A.L. [1975]. The heat balance and associated runoff from an experimental site on a glacier tongue. [Union Géodésique et Géophysique Internationale. Association Internationale des Sciences Hydrologiques. Commission des Neiges et Glaces.] Symposium. Neiges et glaces. Actes du colloque de Moscow, août 1971,p. 5969.(IAHS–AISH Publication No. 104.)Google Scholar
Geiger, R. 1961. Das Klima der hodennahen Luftschicht. Vierte Auflage. Braunschweig, Friedrich Viewcg und Sohn. [English translation: The climate near the ground. Translated by Scripta Technica, Inc. Cambridge, Mass., Harvard University Press, 1965.]Google Scholar
Iwata, S. and others. 1980. Surface morphology in the ablation area of the Khumbu Glacier. by Iwata, S. Watanabe, O. and Fushimi, H. Seppyō, Vol. 41, Special Issue, p. 917.Google Scholar
Jumikis, A.r. 1977. Thermal geotechnics. New Brunswick, N.J., Rutgers University Press. Google Scholar
Kersten, M, S. 1949. Thermal properties of soils. Bulletin of the University of Minnesota Engineering Experinment Station. No. 28. Google Scholar
Kondo, J. 1967. Analysis of solar radiation and downward long–wave radiation data in Japan. Scientific Report of Tohoku University. Ser. 5, Vol. 18, No. 3, p. 91124.Google Scholar
Kraus, H. [1975.] An energy balance model for ablation in mountainous areas. [Union Géodésique el Géophysique Internationale. Association Internationale des Sciences Hydrologiques. Commission de Neiges el Glaces. Symposium. Neiges el glaces. Actes du colloque de Moscow, août 1971.p. 7482.(IAHS-AISH Publication No.104).Google Scholar
Munro, D.S. and Davies, J.A. 1977. An experimental study of the glacier boundary layer over melting ice. Journal Of Glaciology, Vol. 18, No. 80. p. 425–36.Google Scholar
Munro, D.S. and Davies, J.A. 1978. On fitting the log-linear model to wind speed and temperature profiles over a melting glacier. Boundary–Layer Meteorology, Vol. 15, No. 4, p. 423–37.CrossRefGoogle Scholar
Nakawo, M. and Young, G.J. 1981. Field experiments to determine the effect of a debris layer on ablation of glacier ice. Annals of Glaciology, Vol. 2. p. 8591.CrossRefGoogle Scholar
Naruse, R. and others. 1970. Kion-yūsetsu no yagai kenkyū [Field studies on snow melt due to sensible heal transfer from the atmosphere]. [By] R. Naruse, H. Ōura. [and] K. Kojima. Teion-kagaku: Low Temperature Science. Ser. A.[No.] 28, p. 191202.Google Scholar
Penndorf, r. 1956. Luminous reflectance (visual albedo) of natural objects. Bulletin of the American Meteorological Society, Vol. 37, No. 4, p. 142–44.Google Scholar
Penner, E. and others. 1975. Thermal conductivity laboratory studies of some Mackenzie Highway soils, by Penner, E. Johnston, G, H. and Goodrich, L. E. Canadian Geotechnical Journal. Vol. 12, No. 3, p. 271–88.Google Scholar
Figure 0

Table I. Nomenclature

Figure 1

Table II. Meteorological variables during the experimental period

Figure 2

Table III. Meteorological variables during the experimental period

Figure 3

Table IV. Values for the meteorological variables when the surface temperature was measured

Figure 4

Fig.1. Ablation rate versus thermal resistance during 20–21 August(a)and 21–22 August(b).Solid and open circles are observed data during night time and day time respectively. Short dashed lines show the estimation from meteoro logical variables for night–time. Solid and dashed lines represent the calculation for wet and dry surface respectively during daytime.