Hostname: page-component-7c8c6479df-p566r Total loading time: 0 Render date: 2024-03-28T17:36:32.227Z Has data issue: false hasContentIssue false

The Basal Sliding of a Thinning Ice Sheet, Mizuho Plateau, East Antarctica

Published online by Cambridge University Press:  30 January 2017

Shinji Mae*
Affiliation:
National Institute of Polar Research, Kaga, Itabashi-ku, Tokyo, Japan
Rights & Permissions [Opens in a new window]

Abstract

The Japanese Antarctic Research Expedition observed the thinning of the ice sheet, about 70 cm/year, in Mizuho Plateau. The thinning observed is analysed using an equation of mass continuity. The result of the analysis indicates that the thinning is predominantly caused by the basal sliding and the basal sliding velocity is about 10 m/year. This sliding velocity is compared with the basal sliding velocity obtained by the calculation of the velocity due to internal deformation of ice sheet.

Résumé

Résumé

L′expédition antarctique japonaise a observé I’amincissement de la couverture de glace, qui est d′environ 70 cm par an, sur le Mizuho Plateau. Cet amincissement observé est analysé en utilisant une équation de continuité de masse. Les résultats de l′analyse indiquent que l′amincissement a pour cause prédominante le glissement au sol et que la vitesse du glissement au sol est d′environ 10 m/an. Cette vitesse de glissement est comparée avec celle obtenue par le calcul à partir de la déformation interne de la couverture de glace.

Zusammenfassung

Zusammenfassung

Die Japanese Antarctic Research Expedition stellte eine Ausdünnung des Eisschildes auf dem Mizuho-Plateau um ca 70cm pro Jahr fest. Die beobachtete Ausdünnung wird mit Hilfe einer Massenkontinuitäts-gleichung analysiert. Es zeigt sich, dass die Ausdünnung im wesentlichen durch das Gleiten am Untergrund verursacht wird, dessen Geschwindigkeit etwa 10 cm pro Jahr beträgt. Diese Gleitgeschwindigkeit wird mit der aus der inneren Deformation des Eisschildes berechneten verglichen.

Type
Research Article
Copyright
Copyright © International Glaciological Society 1979

1. Introduction

The Japanese Antarctic Research Expedition (JARE) carried out its glaciological research programme in Mizuho Plateau during the period from 1968 to 1975. The main objects of the programme were the measurements of mass balance, ice flow, and surface and sub-glacial topography in Mizuho Plateau Reference Shimizu and Ishida(Shimizu, 1978). The measurements of ice flow were made using the triangulation net set up along the traverse route shown in Figure 1 (Reference Mae and NaruseNaruse, 1978, Reference Naruse1979). The possible error of the survey resulting from refraction in the atmosphere and curvature of the Earth was counterbalanced because the vertical angles were measured twice from two positions of the triangulation net in opposite directions, and the angles were averaged. The triangulation net, of which the starting point was fixed on rock in the Yamato Mountains, was surveyed by the 10th JARE during the period between November and December 1969. It was resurveyed by the 14th JARE during the period between December 1973 and January 1974.

Fig. 1. Map of the Mizuho Plateau. A thick line represents the traverse route. Thin solid lines represent equal-elevation contours. Thin dashed lines represent the boundary of the Shirase Glacier drainage basin.

Reference Mae and NaruseNaruse (1978 Reference Naruse1979) analysed the survey result of the triangulation net and estimated the horizontal and vertical velocities of the ice surface and strain-rate at the surface. Taking into account the effect of the firn densification, he estimated the rate of change h of the ice thickness h. The mean-square error of the vertical movement is 18 cm/year at point A and 31 cm/year at point c. The mean-square error of h is approximately equal to that of the vertical movement. The results of this estimation revealed that the main part of the route in the Mizuho Plateau was thinning by 70 cm/year. Along the route, the accumulation-rate of snow was measured during this period Reference Naruse, Yokoyama and Naruse(Yokoyama, 1975) and the ice thickness and bedrock topography were observed using radio echo-sounding Reference Naruse, Yokoyama and Naruse(Naruse and Yokoyama, 1975; Reference Shimizu, Shimizu, Yoshimura, Naruse, Yokoyama and IshidaShimizu and others, 1978). Reinterpreting the radio echo observed, Reference MaeMae (1978) obtained improved bedrock topography.

Reference MaeMae (1977) and Reference Mae and NaruseMae and Naruse (1978) investigated a possible cause of the thinning of the ice sheet in Mizuho Plateau and reported their preliminary conclusion that the thinning was not caused by a decrease of the accumulation-rate due to a change in climate but by basal sliding of which the velocity was about 10 m/year. In this paper the preliminary work is developed and the stability of the ice sheet in Mizuho Plateau is discussed.

2. Basal Sliding Velocity Estimated using the Equation of Mass Continuity

The values of h calculated by Reference Mae and NaruseNaruse (1978) are shown in Figure 2. The error in the vertical movement caused by refraction in the atmosphere and curvature of the Earth was eliminated and the mean-square error of the vertical movement is 18 cm/year at point A (Fig. 1) and 31 cm/year at point c (Fig. 1)Reference Mae and Naruse (Naruse, 1978). The error in h is approximately equal to that of the vertical movement. Thus it is seen from Figure 2 that h is about zero between points A and B and about –70 cm/year between points A and C.

Fig. 2. The rate of change h of the ice thickness along the traverse route.

Reference Naruse, Yokoyama and NaruseYokoyama (1975) measured the accumulation-rate of snow along the route. The accumulation-rate d of ice calculated from Yokoyama′s measurement is about 7 cm/year on average along the route and even the maximum value of ȧ is about 30 cm/year. In the coastal area near Syowa Station (Fig. 1) where cyclones can occur and the amount of snow precipitation is large, ȧ rarely exceeds about 50 cm/year (Reference Yamada, Yamada, Okuhira, Yokoyama, Watanabe and IshidaYamada and others, 1978). This suggests that h is too large to be compensated by the present value of snow accumulation-rate observed in the Mizuho Plateau and temporal variation of the accumulation-rate cannot account for the large value of – h.

For the study of the thickness change of the ice sheet, a very useful equation is that of mass continuity. The equation of mass continuity for an ice sheet is given by Reference Shumskiy and BauerShumskiy and Bauer (1965). Taking the x-axis along the horizontal component of flow velocity, they obtain the equation

(1)

where m is the mass of a vertical column with unit area of horizontal cross-section at a given point (x, y), ρ(z) is the density of ice and firn, ε̇ is strain-rate, is the accumulation-rate at bottom, u is the velocity in the x direction, and α is the slope. Roman subscripts s (for surface) and b (for bottom) denote the position where the value is to be taken. The ice thickness is given by h = za—Zb.

In the present case, since the firn densification is eliminated in the calculation of , p is approximately constant and equal to the density of ice and ∂p/∂x = 0. αs and αb < 1, so tan αs ~ αs and tan ~ a∞. If we assume that the vertical strain-rate is independent of depth and the freezing or melting rate is neglected, Equations (1) can be rewritten

(2)

Substituting the observed values of into Equation (2), we obtain ubαb which is shown in Figure 3.

Fig. 3. ub ∞b (thick line) and (thin line) along the traverse route.

The flow lines on the traverse route generally converge slightly Reference Mae and Naruse(Naruse, 1978). However, when the convergence or divergence of the flow lines is large, the strain-rate, , is extraordinarily large or small, as shown in Figure 3; we should not use Equations (2) at such points. Excluding the large strain-rates, Figure 3 indicates that ubαb ≈ 0 between points A and B and ubαb = 55 cm/year in average between points A and c. This large value of ubαb between points A and c suggests that ub is significantly large and the thinning of the ice sheet may be predominantly caused by the basal sliding because The result of the radio echo-sounding shows that the slope of the bedrock is large below the route while it is small down-glacier Reference Shimizu, Shimizu, Yoshimura, Naruse, Yokoyama and Ishida(Shimizu and others, 1978; Reference MaeMae, 1978). Reference MaeMae (1977) estimated αb at point c to be about 4X 10−2. Since ubαb = 45 cm/year at point c, ub is about 10 m/year.

It is impossible to estimate αb below the traverse route but it appears that αb is approximately equal to αb at point c. Therefore, between points A and c the velocity of basal sliding may be 10 m/year and between points A and B it may be zero.

We assume that are independent of depth. Reference Dansgaard and JohnsenDansgaard and Johnsen (1969) calculated the temperature profile in the ice sheet and compared it with that observed at Camp Century, Greenland. From the comparison they showed that the profile was considerably improved if they used a vertical strain-rate that was uniform down to some 400 m or so above the bottom of the ice, and then decreased to zero at bedrock, instead of using a uniform strain-rate (the ice thickness at Camp Century is about 1400 m). If the model with a non-uniform strain-rate is right, the term of Equations (1) decreases by about 20% while ubαb increases by about 6 cm/year between points A and B. If we take 0.7 Mg/m3 for ρs , the summation ρsȧ + ρsus tan αs of Equation (1) decreases by about 5 cm/year but the other terms of Equations (1) hardly change, so that Ubαb decreases by about 5 cm/year between points A and c about 3 cm/year between points A and B. Thus, the variation of ubαb due to the assumptions is within the mean-square error of h and can be ignored in the estimation of ub .

3. Basal Sliding Velocity Estimated Using the Flow Law of Ice

The surface velocity us is given by

(3)

where u p is the velocity due to the plastic deformation of the ice mass. The calculation of us is complicated because of the dependence of strain-rate upon temperature. The temperature profile in the Mizuho Plateau has not yet been observed. However, Nishio and Mae (in press) calculated the temperature profile using a steady-state model used by Reference RobinRobin (1955) under an assumption that the icesheet was in equilibrium, because Reference MaeMae (1977) and Reference Mae and NaruseMae and Naruse (1978) considered that the thinning commenced within 103 years B.P. and F. Nishio (private communication) proposed, on the basis of the analysis of the temperature observations of a core hole at Mizuho Station, that the thinning commenced about a few hundred years ago. The profile calculated shows that the base temperature is at the melting point between points A and c and the base is frozen near the Yamato Mountains. This may be evidence that basal sliding occurs between points A and c.

H. Shōji (private communication; an abstract is presented in this symposium, Reference Shōji and HigashiShōji and Higashi (1979)) carried out deformation tests under high hydrostatic pressures with core-ice specimens obtained from various depths at Byrd Station in order to examine the deformation behaviour of ice masses inside the ice sheet. The results of the test revealed that above a depth of about 1200 m the ice was soft because of air bubbles, while below this depth the ice became hard because of the closure of air bubbles, and the relation between the shear strain-rate on the slip plane γ̇xz and the shear stress τxz is as follows: depth d ⩽ 1200 m

where Q = 66 kJ/mol and R is the gas constant; 1200 m ⩽ dD (the depth where the temperature is —8°C)

where Q = 66 KJ/mol; D ⩽ d

where Q = 126 kJ/mol.

Using Shōji′s result and the temperature profile obtained by Nishio and Mae (in press) and assuming that the ice-sheet flow is laminar, up is calculated as shown in Figure 4. The sliding velocity ub is given by us—up; it is about 11 to 15m/year, which is approximately equal to the value for uh obtained in the previous section. If we assume that the stress-strain-rate relation is expressed by

throughout the whole depth and ub is about 10 m/year, we obtain B = 3 X 103 bar−3 s−I and Q = 66 kJ/mol. These values are approximately equal to those used by Reference HughesHughes (1973) for his calculation of the internal deformation velocity of the West Antarctic ice sheet.

Reference PatersonPaterson (1977) obtained the stress-strain-rate relation using the measured bore-hole closure rates. The data of the rates cover temperatures between — 16°C and —28°C anddepths down to about 300 m. The obtained relation is given by

whereQ = 54 kJ/mol. For a given stress the strain-rates measured by Shōji are 2 to 8 times larger than Paterson′s result. If we use Paterson′s result, the calculated up is 3 to 4 times smaller than the up shown in Figure 4.

Fig. 4. The horizontal surface velocity, Us Reference Mae and Naruse(from Naruse, 1978), and the velocity up due to plastic deformation along the traverse route are represented by and ×. The thick line represents the smoothed curve of up.

4. Basal Shear Stress

The thinning caused by basal sliding suggests that the ice sheet of the Mizuho Plateau may be in an unstable state at present. The basal shear stress τb is an important parameter to examine the state of ice sheet. It is given by

(4)

Substituting the observed values of h and αs along the traverse route into Equations (4), we obtain τb, as shown in Figure 5. It is seen from Figure 5 that τb is approximately 1 bar. The deviation of τb from 1 bar is mainly caused by the fluctuation of αs due to the micromorphology of the ice surface.

Fig. 5. The basal shear stress τb, along the traverse route.

The main part of the traverse route flows into the Shirase Glacier, as shown in Figure 1. Then, τb along the central flow line of the Shirase Glacier drainage basin and τb along the ridge of the basin, the Syowa Station flow line, are calculated based on the surface and bedrock morphology shown in Figure 6 and JARE observations. τb along the central flow line increases with increasing distance from the coast, reaches a maximum value of 1.8 bar around 100 km from the coast, and then gradually decreases. On the other hand, τb along the Syowa Station flow line decreases drastically near the coast and shows a slight decrease in the inland area. The increase in τb along the central flow line for 100 km from the coast arises from the increase in h, while the decrease in Tb beyond 100 km from the coast is caused by the decrease in αs. The curve of τb along the central flow line is very similar to those of glaciers in a pre-surge state Reference Robin and Weertman(Robin and Weertman, 1973). This suggests that the ice sheet in the Mizuho Plateau may not be in a stable state.

Fig. 6. The surface and bedrock topography along the central flow line.

5. Discussion

The estimation of the basal sliding velocity indicates that basal sliding is occurring over a main part of the traverse route; that is at least over that section of the route in the northern part of the Shirase Glacier drainage basin, basal sliding occurs and the base is wet. This may be a reason for the rapid motion of the Shirase Glacier, us = 2.2 km/year Reference Fujiwara and Yoshida(Fujiwara and Yoshida, 1972). From Figure 1 it is seen that the surface elevation is lower in the central area of the Shirase Glacier drainage basin than near the boundary of the basin. If the lowering of the surface is caused by thinning of the ice sheet due to basal sliding, the thinning commenced about 2 X 102 years ago because the lowering is about 100 m at the traverse route. This result for the age of the thinning is supported by analysis of the temperatures of a core hole at Mizuho Station. Therefore, we can propose that the area of basal sliding and wet base has been expanding from the northern part of the basin and reached the traverse route about 200 years ago provided that the thinning-rate was constant.

If the basal sliding is caused by a rise of ice temperature at the base, the following two causes are conceivable:

  • (1) a decrease in the accumulation-rate,

  • (2) a rise of air temperature.

Reference RobinRobin (1955) calculated the temperature profile of an ice sheet based on a simple steady-state model of heat conduction, and showed that the temperature difference between the surface and base increases as the accumulation decreases. A propagation time of temperature change from the surface to the base is 104 to 105 years. However, we have no information on accumulation-rate in the Mizuho Plateau since 104 to 105 years B.P., and it is impossible at present to determine which is the major cause of the rise of the base temperature.

It is conceivable that the Shirase Glacier should be classified as a fast-moving outlet glacier of the ice sheet and that the ice sheet of the Shirase Glacier drainage basin has the basal shear-stress curve shown in Figure 7 in a steady way. If this hypothesis were right, we could not explain the thinning of the ice sheet observed along the traverse route, nor our conclusion that the thinning commenced about 2 X 102 years agoand the wet-base area has been expanding from the down-glacier part of the basin.

Fig. 7. The basal shear stress τb τalong the central flow line and the Syowa Station flow line.

If a large amount of thinning occurs in West Antarctica, the mechanism of that thinning may be different from that in the Mizuho Plateau. The West Antarctic ice sheet forms one continuous body of ice with the Ross Ice Shelf. The base of the West Antarctic ice sheet over much of its area is at the startling depth of 0.5 to 1 km below sea-level. Moreover, the bedrock would remain below the sea-level over a large area even if all the ice were removed and iso-static rebound occurred. Contrary to this, the base of the ice sheet in the Mizuho Plateau is above or at sea-level. For example, the base below the traverse route is about 500 m above sea-level and the base at the maximum τb is at sea-level.

Discussion

W. F. BUDD: It is difficult to calculate the internal deformation because one needs to know the temperature distribution and the ice crystallography. What is the error on your estimate for internal deformation?

S. Mae: I think that the error of the calculated internal deformation is mainly caused by the calculation ofthe vertical temperature profile because I used a simple model. I believe that the stress–strain-rate relation obtained by Reference Shōji and HigashiShōji and Higashi (1979) is more reliable than the temperature profile, and the ice crystallography causes a rather smaller error than the temperature profile.

R. H. Thomas: Accumulation measurements in the catchment area of Shirase Glacier indicate that total accumulation is greater than outflow. How do you reconcile this with your estimates of thinning?

Mae: If the estimation of the total ice accumulation is correct, the positive mass balance in the Shirase Glacier drainage basin indicates that a thickening may take place down-glacier. We have not measured the thickness change in the other part and it is not clear whether thickening is taking place or not.

References

Barnes, P., and others. 1971. Friction and creep of polycrystalline ice, by Barnes, P., Tabor, D., and Walker, J. C. F.. Proceedings of the Royal Society of London, Ser. A, Vol. 324, No. 1557, p. 127–55.Google Scholar
Dansgaard, W., and Johnsen, S. J. 1969. A flow model and a time scale for the ice core from Camp Century, Greenland. Journal of Glaciology, Vol. 8, No. 53, p. 215–23.Google Scholar
Fujiwara, K., and Yoshida, Y. 1972. Ryutso-Horumu-wan tōgan no hyōga ryūdō sokutei [Ice flow measurements on the east coast of Lützow-Holm Bay [Antarctica]]. Nankyoku Shiryō: Antarctic Record, No. 44, p. 7992.Google Scholar
Hughes, T. J. 1973. Is the West Antarctic ice sheet disintegrating? Journal of Geophysical Research, Vol. 78, No. 33, p. 7884910.Google Scholar
Mae, S. 1977. Hyōshō oyobi hyōga no hyōkōhenka no gen-in ni tsuite—nankyoku Mizuho-kōgen hyōshō to Nepāru Himaraya Kunbu hyōga no baai [The variation of the thickness of the Mizuho Plateau ice sheet, East Antarctica, and the Khumbu glacier, Nepal Himalaya]. Seppyō, Vol. 39, No. 3, p. 117–24.Google Scholar
Mae, S. 1978. The bedrock topography deduced from multiple radar echoes observed in the Mizuho Plateau, East Antarctica. Nankyoku Shiryō: Antarctic Record, No. 61, p. 2331.Google Scholar
Mae, S., and Naruse, R. 1978. Possible causes of ice sheet thinning in the Mizuho Plateau. Nature, Vol. 273, No. 5660, p. 291–92.Google Scholar
Naruse, R. 1978. Surface flow and strain of the ice sheet measured by a triangulation chain in Mizuho Plateau. (In Ishida, T., ed. Glaciological studies in Mizuho Plateau—west Enderby Land, East Antarctica, 1969–1975. Tokyo, National Institute of Polar Research, p. 198226. (Memoirs of National Institute of Polar Research. Special Issue No. 7. ))Google Scholar
Naruse, R. 1979. Thinning of the ice sheet in Mizuho Plateau, East Antarctica. Journal of Glaciology, Vol. 24, No. 90, p. 4552.Google Scholar
Naruse, R., and Yokoyama, K. 1975. Position, elevation, and ice thickness of stations. (In Naruse, R., ed. Glaciological research program in Mizuho Plateau—west Enderby Land, East Antarctica. Pt. 3, 1973–1974. Tokyo, National Institute of Polar Research, p. 747. (Japanese Antarctic Research Expedition. JARE Data Reports, No. 28 (Glaciology).))Google Scholar
Nishio, F., and Mae, S. In press. Temperature profile in the base ice area of the Yamato Mountains. [Memoirs of National Institute of Polar Research.]Google Scholar
Paterson, W. S. B. 1977. Secondary and tertiary creep of glacier ice as measured by borehole closure rates. Reviews of Geophysics and Space Physics, Vol. 15, No. 1, p. 4755.Google Scholar
Robin, G. de Q. 1955. Ice movement and temperature distribution in glaciers and ice sheets, Journal of Glaciology, Vol. 2, No. 18, p. 523–32.Google Scholar
Robin, G. de Q., and Weertman, J. 1973. Cyclic surging of glaciers. Journal of Glaciology, Vol. 12, No. 64, p.318.Google Scholar
Shimizu, H. 1978. Outline of the studies of the glaciological research program in Mizuho Plateau, East Antarctica, 1969–1975. (In Ishida, T., ed. Glaciological studies in Mizuho Plateau—west Enderby Land, East Antarctica, 1969–1975. Tokyo, National Institute of Polar Research, p. 113. (Memoirs of National Institute of Polar Research. Special Issue No.7.))Google Scholar
Shimizu, H., and others. 1978. Morphological feature of the ice sheet in Mizuho Plateau, [by] Shimizu, H., Yoshimura, A., Naruse, R., and Yokoyama, K. (In Ishida, T., ed. Glaciological studies in Mizuho Plateau—west Enderby Land, East Antarctica, 1969–1975. Tokyo, National Institute of Polar Research, p. 1425. (Memoirs of National Institute of Polar Research. Special Issue No. 7.))Google Scholar
Shōji, H., and Higashi, A. 1979. Mechanical properties of Antarctic deep-core ice. Journal of Glaciology, Vol. 24, No. 90, p. 487–89,Google Scholar
Shumskiy, P. A., and Bauer, A. 1965. Issledovaniye sovremennykh izmeneniy tsentral′noy chasti lednikovogo pokrova vostochnoy Antarktidy v 1964 g. [Study of recent changes in the central part of the ice cap in East Antarctica in 1964]. Informatsionnyy Byulleten′ Sovetskoy Antarkticheskoy Ekspeditsii, No. 51, p. 311–16.Google Scholar
Yamada, T., and others. 1978. Distribution of accumulation measured by the snow stake method in Mizuho Plateau [by] Yamada, T., Okuhira, F., Yokoyama, K., and Watanabe, O.. (In Ishida, T., ed. Glaciological studies in Mizuho Plateau—west Enderby Land, East Antarctica, 1969–1975. Tokyo, National Institute of Polar Research, p. 125–39. (Memoirs of National Institute of Polar Research. Special Issue No. 7.))Google Scholar
Yokoyama, K. 1975. Net accumulation by stake measurements. (In Naruse, R., ed. Glaciological research program in Mizuho Plateau—west Enderby Land, East Antarctica. Pt. 3, 1973–1974. Tokyo, National Institute of Polar Research, p. 6282. (Japanese Antarctic Research Expedition. JARE Data Reports, No. 28 (Glaciology).))Google Scholar
Figure 0

Fig. 1. Map of the Mizuho Plateau. A thick line represents the traverse route. Thin solid lines represent equal-elevation contours. Thin dashed lines represent the boundary of the Shirase Glacier drainage basin.

Figure 1

Fig. 2. The rate of change h of the ice thickness along the traverse route.

Figure 2

Fig. 3. ub ∞b (thick line) and (thin line) along the traverse route.

Figure 3

Fig. 4. The horizontal surface velocity, Us (from Naruse, 1978), and the velocity up due to plastic deformation along the traverse route are represented by and ×. The thick line represents the smoothed curve of up.

Figure 4

Fig. 5. The basal shear stress τb, along the traverse route.

Figure 5

Fig. 6. The surface and bedrock topography along the central flow line.

Figure 6

Fig. 7. The basal shear stress τb τalong the central flow line and the Syowa Station flow line.