Introduction
In China research on existing glaciers began in 1958. Early in the 1960’s the glaciers in West China were classified into continental type and maritime type according to their main properties and geographical environments. To the former belong those in most parts of Ch’ing-hai-Hsi-tsang kao-yuan (Qinghai-Xizang gao-yuan–Tibetan Plateau) from the middle Himalaya Shan to Ch’i lien Shan (Qilian Shan) and T’ien Shan (Tian Shan) with low precipitation and low temperature, while to the latter those in south-eastern parts of the Plateau with high monsoon precipitation (Reference Shi and XieShi and Xie, 1964; Reference ShiShi and others, 1980). Among them the continental-type glaciers have been researched rather more, and their temperature measurement is taken as a significant project in their investigation. Since the first measurement (Reference DolgushinDolgushin, 1961) taken in Crfi lien Shan (Qilian Shan) in 1959, temperature has so far been measured at 24 sites in 12 glaciers (Fig. 1). In many glaciers the temperature of the snow cover has also been measured. Furthermore, Chinese glaciologists have taken temperature measurements on Batura Glacier (lat. 36° 35’ N., long. 74°23’E.) (Batura Glacier Investigation Group, 1979) and Bait Bare Glacier (lat. 36° 20’ N., long. 74° 52’ E.) (Reference WangWang and others, in press) in Karakoram Mountains for a short time.
The ice temperature is usually measured using a copper resistance thermometer measured by a Wheatstone’s bridge. For laying the thermometer, a hand drill and a steam drill are used to make holes, the maximum depth of which reaches 18 m. Most drill holes are located in the ablation zone, while a few are in the accumulation zone.
The glaciers Ku-hsiang ping-ch’uan (Gu-xiang bing-chuan) and A-cha ping-ch’uan (A-zha bing chuan) are located in south-eastern Hsi-tsang (Xizang-Tibet) as Figure 1 shows. They have temperatures equal to 0°C or so and belong to the maritime type, with which we will not be dealing in this paper for the time being. The rest of the glaciers shown in Figure 1 are of continental type. Reference Shi and XieShi and Xie (1964), Division of Glaciology, Lanzhou Institute of Glaciology, Cryopedology, and Desert Research (1975), and Reference ShiShi and others (1980) have reviewed their temperature. Adding some of the latest data, we are going to undertake a comprehensive discussion of the temperature regime of the continental-type glaciers here.
Ice temperature regime in the active layer
Figure 2 shows the temperature stratification measured in No. 1 glacier, Urumchi river, weekly from June to September 1962, the first continuous measurements for a relatively long time in China. Figure 3 shows the stratification at the same place measured monthly from 1964 to 1966, which is characterized by having the longest continual measuring time to a greater depth in China.
Figure 4 shows the stratification measured in No. 5 glacier, Yang-lung he (Yang-long he), weekly in 1977 with the greatest depth and different levels for comparison.
On the basis of analysing the data of No. 5 glacier, Yang-lung he (Yang-long he), for describing the ice-layer temperature regime with a period of one year, i.e. ω/2π = la-1,(Figure 2 and 3) Reference Ren and HuangRen and Huang (1981) have derived a semi-empirical formula as follows,
where y is the coordinate of depth, with an origin on the glacial surface and a positive direction downwards; t the coordinate of time; T(y, t) is the temperature at a depth y and time t; Ts an apparent wave amplitude at the surface; ω/2π the frequency of temperature change; k the thermal diffusivity; β a factor for correcting the amplitude, being roughly 1.0; T0(y) the equilibrium temperature at depth y, namely the annual mean temperature. 30 April was chosen to be t= 0 so that the phase of the summer temperature wave gave a good fit to the sinusoid in Equation (1). By means of correlation analysis we obtain
(Fig. 4)where α is a constant and γ the vertical temperature gradient.
The above two formulae can be applied to the level near the equilibrium line. In addition the analysis cannot ensure that conduction is the dominant factor for heat exchange when y < 3 m, within which layer the heat exchange is complicated, and so there is some discrepancy between the Equation (1) (also Equation (2)) based on the heat-conduction equation and on the measured data.
Based on the picture described by Equation (1), we can divide the depth range within which we have temperature measurements into three layers: a surface layer, where an obvious temperature change occurs during all seasons of the year with a large amplitude and a variable sign of temperature gradient, so that the heat conducts sometimes upwards and sometimes downwards in this layer; a middle layer, where the temperature change becomes less and less during the year and the temperature gradient does not change sign, so the heat conduction always occurs in one direction in this layer; a deep layer, at the upper bound of which the yearly temperature change is less than 0.2 deg, so this layer can be considered to be lacking in yearly temperature change and the heat always conducts upwards. Sometimes the surface layer plus the middle layer is called the “active layer” (Reference AvsyukAvsyuk. 1956). The boundary between the surface layer and the middle layer is about 6 m deep. The lower bound of the active layer is at a depth of 15-20 m. Taking account of the fact that the maximum depth to which we have measured the temperature is merely 16 or 18 m, we take 16 m as the lower bound of the active layer for the moment. The temperature of the lower bound of the active layer, T0(16), can be considered an important characteristic of the glacial temperature regime, as we are going to mention later on.
Inference of temperature in the deep layer
Measured data in the deep layer are lacking. But we may make an inference from the active layer. One can assume roughly that all the heat generated by friction concentrates at the glacier bottom, so that the geothermal flux combining with the generated heat is considered to be the heat input at the bottom (Reference BuddBudd, 1969). For a cold glacier in thermal steady state, the entire heat input Q at the bottom per unit time, should transmit upwards in time, requiring a temperature gradient G, then,
where λ is the conductivity.
There are only three kinds of temperature boundary condition at the glacial bottom; (a) the temperature of the interface between ice and bedrock is below 0°C, ice and bedrock freeze together, namely are the ice temperature and temperature gradient at the interface respectively, and θ pressure-melting point); (b) the interface is at the pressure melting point, but the bottom ice layer is below pressure-melting point, namely T b=θ and 0 <(∂T/∂y)b <G. a part of the heat input transmits upwards and another part is expended by melting the bottom ice; (c) there is a temperate layer at the bottom; i.e. Tb=θ and (∂T/∂y)b=0 all the input is expended by melting without transmitting upwards.
We can approximately estimate the depth L at which the pressure-melting point occurs, as follows.
for the above three boundary conditions, assuming the γ within L to be the same as that within the active layer. There are the following relations between L and glacier thickness H, and between G and y, corresponding to the above three boundary conditions:
-
L >H and G = γ for condition (a);
-
L = H for condition (b);
-
L < H and G > γ for condition (c).
Equations (2) and (4) are only a very rough approximation to the actual temperature distribution because of upward or downward advection of ice except near the equilibrium line. Fortunately the annual surface accumulation/ablation rate of the continental type glaciers in China is quite low, the gradient of mass balance (activity index) ranges from 2.9 to 6.2 mm/m in "Tien Shan (Tian Shan) and Ch’i-lien Shan (Qilian Shan) for instance (Reference XieXie, 1980), so we can apply these equations approximately to some several hundred metres below equilibrium line.
Table I shows the values of the temperature T0(16) at the lower bound of the active layer, the vertical temperature gradient γ, and the depth L at which the pressure-melting point appears for several glaciers in China where the temperature has been measured to greater depth. The values in Table I are either calculated, if the observed depth is enough, or extrapolated, if the observed depth is less than 16 m, from the measured data by use of Equations (1), (2) and (4).
Now we are going to compare L with H as observed by the gravimeter method. In the glacier Hsi ch’iung-t’ai-lan ping-ch’uan (Xi qiong-tai-lan bing-chuan) H is observed to be 229 m at average altitude of 3 850 m and its maximum value to be 302 m. It will rise at the altitude shown in Table I and so be much greater than L as in condition (c) above mentioned. In Jung-pu ping-ch’uan (Rong-bu bing-chuan) the maximum of H is observed to be 146 m at the altitude of 5 350 m (Reference Xie and SuXie and Su, 1975), so H is greater than L as in condition (c) too. In No. 5 glacier, Yang-lung he (Yang-long he), there are no data for H. But in "Ch’i yi" ping-ch’uan ("Qi Yi" bing-chuan), an adjacent glacier a little bigger than No. 5 glacier, has been observed to be tens of metres with the maximum H of 93 m. Then No. 5 glacier, Yang-lung he (Yang long he), seems to belong to condition (c) or (b). Thus from the above comparison we can infer the temperature features in deep layers of continental-type glaciers in China as follows: (1) bottom ice and the glacier bed freeze together with permafrost below when the upper-bound temperature of the deep layer is quite low or the glacier is too thin; (2) the interface between ice and the bed is at the pressure melting point while the bottom ice layer still below the pressure-melting point when the upper-bound temperature is quite low but the glacier is very thick, or the glacier thickness not too large but the upper-bound temperature rather high; (3) there is a temperate ice layer in the lower part of the glacier on thawed ground (rock) when the upper-bound temperature is rather high and the glacier thickness rather great. Cases (2) and (3) usually predominate.
Discussion of vertical temperature gradient and basal sliding
From Table I one can see that γ values seem to reduce with decreasing elevation in the same glacier and to be 3.5-8.2 times the gradient needed for the geothermal heat flux to travel (Table I) upwards through the glacier. Recently Reference MūllerMuller (1976) has reported a vertical temperature gradient of 0.08 deg/m measured near the bottom of the glacier at site III (middle part) of White Glacier (lat. 80° N.) and considers that in this glacier there is basal sliding on most of the glacier bed. Reference Hooke and KociHooke and Koci (1978) report that the gradient between 19 and 23 m in hole S25 on Sukkertoppen Iskappe, Greenland, is 0.14 deg/ m, which corresponds to 7.5 times the gradient needed for the geothermal heat flux at the place in question. The data published by Reference VilesovVilesov (1961) show that in Lednik Tsentralniy Tuyuksu, Zailiyskiy Alatau, a gradient of 0.11 deg/m was measured between 10 and 20 m depth at No. 5 hole 600 m from the front, but the gradient decreased downards, and a temperature of -0.1 to -0.7 °C appeared at the bottom, 52 m in depth. It follows that in cold glaciers γ should be greater, several times more than the gradient needed for the geothermal flux. The problem is to extrapolate γ from the active layer to the deep layer. According to the circumstances in Lednik Tsentralniy Tuyuksu, which is similar to the continental-type glaciers in China, γ in the active layer is larger than that in the deep layer, so that the value of L calculated from Equation (4) probably is a little small. But the basic fact that basal sliding takes place extensively in glaciers in China cannot be changed in any case.
(The altitude of the equilibrium line of each glacier is shown in Table II)
The surface flow velocity observed in the middle and lower parts of glaciers in China shows that without exception the summer velocity is faster than the winter velocity, or than annual mean velocity. For instance, in Lao-hu-kou ping-ch’uan (Lao-hu-gou bing-chuan) the summer velocity is faster than winter velocity by a factor of 1.6 observed from 1960 to 1961; in No. 1 glacier, Urumchi river, on the average the summer velocity is faster than winter velocity by a factor of 1.05 to 1.90, according to Reference ZhangZhang (1965), who summarized the data observed from 1959 to 1962; in "Ch’i yi" ping-ch’uan ("Qi Yi" bing-chuan), the summer velocity is faster than the winter velocity by a factor of 1.06 on the average over seven sections observed from August 1976 to September 1977; in No. 5 glacier, Yang-lung he (Yang-long he), the summer velocity is faster than winter velocity by a factor of 1.22 on the average over eight sections observed from June 1977 to August 1978 ; in Hsi-ch’iung-t’ai-lan ping-ch’uan (Xi-qiong-tai-lan bing-chuan) the velocity from 10 July to 28 July 1977, was faster than that from July 1977 to July 1978 by a factor of 1.43 on the average over five sections. This phenomenon indicates that the pressuremelting point appears at the glacier bed, where basal sliding exists, while the melt water percolates to the bed during the ablation period, lubricating and promoting basal sliding.
The lower-bound temperature of the active layer at the elevation of the equilibrium line
Reference Shi and XieShi and Xie (1964), Reference Huang and YuanHuang and Yuan (1965), Division of Glaciology, Lanzhou Institute of Glaciology, Cryopedology and Desert Research (1975), Reference ShiShi and others (1980) and Reference Zhong-guo, Ke-xue-yuan, Lanzhou, Bing-chuan, Dong-tu, Yan-jiu-suo, Qilian Shan, Bing xue Liyong and Yan-jiu duiZhong-guo Ke-xue-yuan Lanzhou Bing-chuan Dong-tu Yan-jiu-suo (1980) have discussed the regional variation of glacier temperatures in China, giving an impression that the lowest ice temperatures arise in the western section of Ch’i-lien Shan (Qilian Shan), the place where glacier temperatures may be the lowest of all alpine glaciers located in the middle and low latitudes. From there eastwards to middle and eastern section of Ch’i-lien Shan (Qilian Shan), westwards to the middle and western sections of T’ien Shan (Tian Shan) and southwards to the middle Himalaya Shan, the glacial temperature will rise gradually, reflecting the increase of continentality towards the western section of Ch’i-lien Shan (Qilian Shan) from the east, west, and south.
In fact in Table I a comparison of T0(l6) is made. But it is better to take the lower-bound temperature, T0(16)E, of the active layer at the elevation of the equilibrium line as a standard for comparison. For this purpose the gradient dT0(16)/dz for T0(16) varying with altitude z must be sought.
The data from No. 5 glacier, Yang-lung he (Yang-long he), extending across the equilibrium line is what we need. We can get dT0(16)/dz = –1.26 deg/100 m between 4513m and 4648 m (Table II) for this glacier from Table I. The absolute value is greater than the annual mean vertical air-temperature gradient, estimated to be –0.74 deg/100 m, by 0.52 deg/100 m, which may result from upward advection and strain heating. Let us apply dT0(16)/dz = –1.25 deg/100 m to extrapolate T0(16)E from the elevations below the equilibrium line shown in Table I for the time being. The result is shown in Table II, the annual mean air temperature TE at the altitude of the equilibrium line is also shown for comparison. We can see from this Table that T0(16)E ranges from –5.5 to –7.7 °C. Mentioning the continental-type temperature regime of glaciers, Reference AvsyukAvsyuk (1956) had indicated that the annual mean temperature of the deep layer is always higher than the annual mean air temperature, and for Soviet Tyan’-Shan’ the former is 2 deg higher than the latter. Now we can exactly say that the lower-bound temperature of the active layer (usually lower than the mean temperature of the deep layer) at the altitude of the equilibrium line is 1.8 to 3.7 deg higher than the annual mean air temperature at the same height for continental-type glaciers in China.
Vertical zones of the temperature regime in glaciers
1. Ablation
Most temperature measurements we have taken are in the ablation zone, where a temperature below zero is always measured. The temperature regime of the active layer in the ablation zone can be described by Equation (1) except in the first 3 m, where a complicated heat exchange exists and a zero temperature casing arises in the warmest months. In many small-scale glaciers such as No. 1 glacier, Urumchi river, and No. 5 glacier, Yang-lung he (Yang-long he), a zero temperature casing of less than 0.5 m was found. The lower-bound temperature of the active layer is generally higher than the annual mean air temperature at the same height, for example, 3.7 deg higher at an altitude of 4513 m in No. 5 glacier, Yang lung he (Yang-long he). The pressure-melting point usually appears at the bottom layer of the glacier, at least at the interface between ice and bedrock. Basal sliding exists except at some termini of small glaciers and their margins or some extremely small glaciers like hanging glaciers, which may be entirely below zero and freeze onto bedrock.
The extremely large valley glaciers descend to a region where the climate is warm. The temperature in the active layer of its ablation zone will rise as elevation decreases. The rising rate will be greater than the vertical air-temperature gradient, so that the difference between the annual air temperature and the lower-bound temperature of the active layer gradually increases (of course the pressure melting point is a limit for the rise of the lower-bound temperature). Finally the ablation zone changes into a temperate glacier.
Table III shows the ice temperature of Batura Glacier measured in 1974 by the Batura Glacier Investigation Group, Lanzhou Institute of Glaciology and Cryopedology (1979). Its equilibrium line is located at an altitude between 4700 m and 5300 m and its terminus descends down to the bend of the Hunza River at 2 540 m a.s.l. One can see from this Table that its temperature rises as elevation decreases and the lower part of its ablation zone becomes a temperate glacier.
2. Infiltration-congelation zone (superimposed-iee zone)
Equation (1) can be applied to describe the temperature in the active layer except in the upper 3 m. The lower bound temperature of the active layer is lower than that in the ablation zone, but is higher than the annual mean air temperature at the same height by a difference smaller than that in the ablation zone, of 3.0 deg, for example, at an altitude of 4648 m in No. 5 glacier, Yang-lung he (Yang-long he). The pressure melting point generally occurs at the bottom except in some glaciers with extremely small thickness.
The heat exchange in the surface layer is complicated by accumulation, ablation, and the infiltration and recongelation of melt water. Table IV shows the temperature data for the snow layer taken at an altitude of 4 650 m in No. 5 glacier, Yang-lung he (Yang-long he), by means of a quartz thermometer with an accuracy up to ±0.01 deg. One can see from Table IV that the snow layer retains quite a low temperature before the melting season while its temperature rises rapidly once the melting season has come. During the main melting season at between 20 cm and 40 cm depth the snow temperature is near 0 °C all day long, with a minimum of –0.03 °C before dawn and a maximum over 0.00 °C at noon, resulting from the thermometer absorbing transmitted radiation.
3. Infiltration zone
Accumulation, ablation, and the infiltration and recongelation of melt water play a leading role in the surface layer, profoundly affecting the temperature regime in the active layer. The temperature of the entire active layer in this zone is higher than that in the other zones below it. For example, in No. 5 glacier, Yang-lung he (Yang-long he),T0(16) at the altitude of 4835 m is (Table V and Table VI) higher than that at the other two sites (Table I), hence the difference between T0(16) and the annual mean air temperature at the same height grows to 7.3 deg. Equation (1) can still be applied to describe the temperature in the active layer except in the upper 3 m. It will be at the pressure-melting point at the bottom in all probability. A temperature of 0°C will appear and remain in the upper snow layer in summer. In the snow layer nearest the surface the temperature will lower quickly in winter, with an extreme temperature lower than that in other zones below, but in inner snow layers, the temperature will not lower so much, because there is liquid water within it which releases latent heat as a result of recongelation. Thus a special temperature regime is formed: a quite low surface temperature followed by a quite high inner temperature; a large surface temperature amplitude accompanies a small inner temperature amplitude. Table V shows the temperature before the melting season measured in the infiltration zone of No. 1 glacier, Urumchi river. Table IV shows the temperature after the melting season measured in the same zone. One can see a rapid temperature rise before the melting season and a slow lowering after it, the latter, of course, results from the melt-water store and its recongelation. During the days of intense melting and infiltration there is a snow and ice layer of at least 4 m falling to the melting point.
(from Reference Huang and YuanHuang and Yuan, 1965)
(from Reference Huang and YuanHuang and Yuan, 1965)
It is to be expected that infiltration and recongelation would weaken and disappear with increasing height, and that the primary influence of climatic factors on ice temperature would then again be given full play.
Summing up the above ideas, we can draw a scheme to show how T 0(16) changes with glacial zones as Figure 5.
Acknowledgements
In this paper we cite many data from Lanzhou Institute of Glaciology and Cryopedology, Academia Sinica that have never previously been published, in particular the data of temperature in Hsi-ch’iung-t’ai lan ping-ch’uan (Xi-qiong tai-lan bing chuan) from Comrade Wang LilunFootnote *, of glacier thickness from Comrade Su ZhentFootnote †, of glacier flow velocity from Comrades Sun Zuozhe and Song GuopingFootnote ‡, and of air temperature gradient from Comrade Ding LiangfuFootnote §, we are grateful to them for their kind support.