Hostname: page-component-848d4c4894-nr4z6 Total loading time: 0 Render date: 2024-05-07T01:55:13.327Z Has data issue: false hasContentIssue false
  • English
  • Français

On the temperature gradient in the upper part of cold ice sheets

Published online by Cambridge University Press:  30 January 2017

K. Philberth  and
B. Fédérer
K. Philberth
Affiliation:
Eidg. Institut fur Schnee- und Lawinenforschung, Davos, Switzerland
B. Fédérer
Affiliation:
Eidg. Institut fur Schnee- und Lawinenforschung, Davos, Switzerland
Rights & Permissions [Opens in a new window]

Abstract

The influence of the surface slope on the temperature profile in the upper part of a cold ice sheet can be described by a function which is independent of the geothermal heat and the heat of friction. This function is calculated for the two-dimensional and the axisymmetric cases. In the two-dimensional case its simplest form is proportional to the horizontal velocity and to the height above bedrock reduced by a constant; another form of this function is approximately proportional to the square of velocity and height.

L'influence de l'inclinaison de la surface sur le profil de température dans la partie supérieure d'une calotte glaciaire froide peut être décrite par une fonction qui est indépendante des chaleurs géothermique et de friction. Cette fonction est calculée pour les cas bidimensionnel et axisymmetrique. Dans le cas bidimensionnel la forme la plus simple de cette fonction est proportionnelle à la vitesse horizontale et à l'altitude sur le lit rocheux diminuée d'une constante; une autre forme est donnée qui est approximativement proportionnelle au carré de ces deux facteurs.

Zusammenfassung

Zusammenfassung

Der Einfluss der Oberflächenneigung auf das Temperaturprofil im oberen Teil eines kalten Eisschildes lässt sich als Funktion besehreiben, die unabhängig ist von der Erd- und der Reibungswärme. Diese Funktion ist für den zweidimensionalen und für den rotationssymetrischen Fall angegeben. Im zweidimensionalen Fail ist ihre einfachste Form proportional der Horizontalgeschwindigkeit und der um eine Konstante reduzierten Höhe üher dem Felsboden; eine andere Form ist angenähert proportional dem Quadrat dieser beiden Grössen.

Type
Short Notes
Information
Journal of Glaciology , Volume 13 , Issue 67 , 1974 , pp. 148 - 151
Copyright
Copyright © International Glaciological Society 1974

I. Introduction

The negative temperature gradient in the upper part of large ice sheets has two sources: The movement of the ice in the x-direction and climatic changes. Both effects result in a temperature decrease with depth of the order of 1 deg. It is possible to draw conclusions about climatic changes if the influence of the movement is determined precisely.

Theoretical and field work on this subject has been published by Reference RobinRobin (1955, Reference Robin1970), Bogoslovski (1958), Reference RadokRadok (1959), Reference WeertmanWeertman (1968), Reference BuddBudd (1969), Reference Budd and RadokBudd and others ([1970]), Reference Radok, Radok, Jenssen and BuddRadok and others {1970), Reference Budd, Budd, Jenssen and RadokBudd and Radok (1971). It seems, however, that the following useful solution for the upper part of ice sheets has been overlooked so far.

II. General solution for the two-dimensional case

The homogeneous form of the heat transport equation under steady-state conditions is

(1)

where κ is the diffusivity of ice, y the height above bedrock, x the distance from the ice divide, vu and v/z the vertical and horizontal velocities respectively.

In the upper part of the ice sheet, vz is considered to be independent of y (Reference PhilberthB. Philberth, 1956; Haefeli, 1961; Reference WeertmanWeertman, 1968; Reference Dansgaard and JohnsenDansgaard and Johnsen, 1969[b]; Reference PhilberthPhilberth, I972[b]), and

can be written

Using the continuity equation we can write

. Hence
, where yr , the “reduced height” is equal to
. The term h0(x) can be visualized as the height below which the ice has zero velocity and above which it is governed by the block flow law. Usually h0(x) is taken to be zero (Reference RobinRobin, 1955; Reference WeertmanWeertman, 1968; Reference BuddBudd, 1969; Reference Radok, Radok, Jenssen and BuddRadok and others, 1970). In this paper we shall use the relation:
(2)

where h = h(x) is the total thickness of the ice sheet, vzm is the mean value of the horizontal velocity over the total ice sheet and vx is its value in the upper part of the ice sheet. Corresponding to yr we shall use the “reduced ice-thickness”

. The ratio hr/h is equal to νxm/vz . For the station Jarl Joset, central Greenland, this ratio is 0.9 (Reference PhilberthPhilberth, 1972[b]).

Using

and neglecting the heat diffusivity in the x-direction, Equations (1) can now be written
(3)

The simplest non-trivial solution of this equation is

(4)

where C 1, and C 0 and C 2 used below, are constants.

If

is small either with respect to
or to 2 (for station Jarl Joset the first term is 0.15; the second term is > 10, down to a depth of 1 km), the following approximate solution can also be used
(5)

Solutions with higher powers of uz need not be considered for practical use.

The combination of Equations (4) and (5) leads to the general solution

(6)

III. General solution for divergent movement

A generalization of Equations (6) for cases of divergent movement can be obtained by introducing a factor ε characterizing the divergence of the streamlines

(7)

In the two dimensional case ε = 1, If the accumulation a and the total depth of the ice h are independent of x, ε is 2 in the axisymmetric case and between 1 and 2 in intermediate cases. If a = a(x) and h = h(x), however, in the axisymmetric case ε is equal to

where z is the horizontal direction orthogonal to the streamline. If these last expressions are not constant a mean value along the streamline has to be used.

IV. Application of the two-dimensional solution

Differentiating Equations (6) with respect to x, using T = Ts (surface temperature) and y = h(x), hence yr = hr , yields

(8)

where we introduced the substitution

, with a = a(x) denoting the ice value of the accumulation rate. In Equations (8) the term
is neglected for the reasons explained in connection with Equations (5). Equations (8) allows us to determine the constants C 1, and C 2.

Under normal conditions dT s/dx, the horizontal gradient of the surface temperature, and a(x) increase with increasing x. If they increase at an equal rate, C- becomes zero. Using Equations (8) with C 2 = 0, Equations (6) becomes

(9)

Near the surface, yrhr approximates y/h.

The solutions (4) and (9) do not take into account the heat of friction and the geothermal heat. On the other hand, Reference RobinRobin (1955), Reference Dansgaard, Johnsen and WeertmanDansgaard and Johnsen (1969[a]) and Reference Philberth and FedererPhilberth and Federer {1971) calculated vertical temperature profiles taking Ts as independent of x but considering the geothermal and frictional heats. Let Tg be such a profile, which is normalized by the addition of a constant so that it becomes zero at the surface.

The real temperature profile taking into account all three influences, is obtained simply by adding Tg to Equations (6), (7) and (9) respectively. This can be explained in the following way: Tg is the solution of a differential equation, which differs from Equations (1) only by a function of x and y on the right-hand side, expressing the heat of friction. In the case of a linear differential equation, the sum of the solutions of the homogeneous plus the inhomogeneous forms is also a solution of the inhomogeneous form. At the surface T g is zero (T = T g), but for y r = 0 the function T and its vertical gradient are negligibly small with respect to T g and its vertical gradient.

Of practical significance is the fact that in the upper part (region with negative temperature gradient) not only T but also T g, can be calculated in a simple way. In the upper part T g can be taken as independent of x (e.g. according to Reference RobinRobin, 1955, Equations (8) or Reference Philberth and FedererPhilberth and Federer, 1971, table I). This can be verified as follows: The upper part of the ice sheet is influenced by a heat of friction which is much smaller than the heat of friction produced below point x, because it originated at a considerably smaller x and vs. Therefore its influence is normally much smaller than the (x-independent) geothermal heat and can be neglected or approximated by a function which does not depend on x. The horizontal gradient

depends on dh/dx and da/dx; but Reference WeertmanWeertman's (1968) Equation (6b) for dT B/dx = 0 (his dT u/dx) shows
to be very small in the range where T — T, is small.

V. An example: station jarl joset (Lat. 33° 30' W., Long. 71° 20' N., a 865 m A.S.L.

General values:

K = 40 m2 a−1; temperature lapse rate λ = 9.5 deg km−1.

Local values :

h = 2.500 km; ht = 2.250 km; x = 125 km (Reference PhilberthPhilberth, 1972[b]).

Values between ice divide and Jarl Joset:

a = 0.30 m a−1 ice value (independent of x; Reference Federer, Federer and SuryFederer and others, 1970), ε - 1 (two dimensional case; Reference PhilberthPhilberth, in press), height of surface = (10) (Mäker, 1964; Reference LliboutryLliboutry, 1968; Reference Philberth and FedererPhilberth and Federer, 1970).

Derived relations :

Multiplication of Equations (10) by λ yields

according to Equations (6) we have:

The comparison of Tfor y r = h r with T g yields:

For x = 125 km (Jarl Joset) the result is:

(11)

Comparison with measured values:

At the depth of 615 m (yr = 1 635 m) -29-30° C and at 1 005 m (y r, = 1 245 m) -30.00° C have been measured by the thermal probe method (Reference PhilberthPhilberth, 1962, Reference Philberth and Federer1970); that is a difference of 0.70 deg. For these two depths the Equations (11) yields a difference of 0.52 deg and the function T g, yields a difference of -0.32 deg. Hence the total amount of the theoretical difference for steady-state conditions is 0.20 deg.

The measured value (0.70 deg) and theoretical value (0.20 deg) differ by 0.50 deg, which can be explained by palaeoclimatic changes. If a temperature jump θ (end of ice age) is assumed to be at 10 000 years B.P. (Reference Dansgaard, Dansgaard, Johnsen and LangwayDansgaard and others, 1969), the 0.50 deg difference corresponds to θ 5 deg, if the jump is assumed to be at 12000 years B.P., it corresponds to θ = 6 deg (Reference PhilberthPhilberth, 1972[a]).

References

Bogoslovskiy, V. N. 1958. The temperature conditions (regime) and movement of the Antarctic glacial shield. Union Géodésique et Géophysique Internationale. Association Internationale d'Hydrologie Scientifique. Symposium de Chamonix, 16—24 sept. 1958, P- 287–305.Google Scholar
Budd, W. F. 1969. The dynamics of ice masses. ANARE Scientific Reports. Ser. A(IV). Glaciology. Publication No. 108.Google Scholar
Budd, W. F., and Radok, U. 1971. Glaciers and other large ice masses. Reports on Progress in Physics, Vol. 34, No. 1, p. 1–70.Google Scholar
Budd, W. F., and others. [1970.] The extent of basal melting in Antarctica, by Budd, W. [F.], Jenssen, D. and Radok, U.. Polarforschung, Bd. 6, Jahrg. 39, Nr. 1, 1969, p. 293–306.Google Scholar
Dansgaard, W., and Johnsen, S. J. 1969[a]. Comment on paper by Weertman, J., “Comparison between measured and theoretical temperature profiles of the Camp Century, Greenland, borehole”. Journal of Geophysical Research, Vol. 74, No. 4, p. 1109–10.Google Scholar
Dansgaard, W., and Johnsen, S. J. 1969[b]. 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
Dansgaard, W., and others. 1969. One thousand centuries of climatic record from Camp Century on the Greenland ice sheet, by Dansgaard, W., Johnsen, S. J., J. MØler and Langway, C. C., Jr. Science, Vol. 166, No. 3903, p. 377–81.Google Scholar
Federer, B., and others. 1970. Outflow and accumulation of ice in Jarl-Joset station, Greenland, by Federer, B., Sury, H. V., K. Philberth and M. [R.] de Quervain. Journal of Geophysical Research, Vol. 75, No. 24, p. 4567–69. Haefeli, R. 1961. Contribution to the movement and the form of ice sheets in the Arctic and Antarctic. Journal of Glaciology, Vol. 3, No. 30, p. 1133–51. Google Scholar
Lliboutry, L. A. 1968. Steady-state temperatures at the bottom of ice sheets and computation of the bottom ice flow law from the surface profile. Journal of Glaciology, Vol. 7, No. 51, p. 363–76.Google Scholar
Mälzer, H. 1964. Das Nivellement über das grönländische Inlandeis der Internationalen Glaziologischen Grönland-Expedition 1959. Meddelelser om GrØnland, Bd. 173, Nr. 7.Google Scholar
Philberth, B. 1956. Beseitigung radioaktiver Abfallsubstanzen. Atomkern-Energie (München), 1. jahrg., Ht. 11–12, p. 396–400.Google Scholar
Philberth, K. 1962. Une méthode pour mesurer les températures a l'intérieur d'un inlandsis. Comptes Rendus Hebdomadaires des Séances de l’Académie des Sciences (Paris), Tom. 254, No. 22, p. 3881—83.Google Scholar
Philberth, K. 1970. Thermische Tiefbohrung in Zentralgrönland. Umschat in Wissenschafl und Technik, 1970, Ht. 16, p. 515–16.Google Scholar
Philberth, K. 1972[a]. Factors affecting deep ice temperatures. Nature, Physical Science, Vol. 237, No. 72, p. 44–45.Google Scholar
Philberth, K. 1972[b]. Über den inneren Warmehaushalt in mächtigen Eisschilden. Polarforschung, Bd. 7, 42. Jahrg., Nr. 1, p. n-17.Google Scholar
Philberth, K. In press. Die thermische Tiefbohrung in Station Jarl-Joset und ihre theoretische Auswertung. Meddelelser om GrØnland. Google Scholar
Philberth, K., and Federer, B. 1970. A note on the surface profile of the Greenland ice sheet. Journal of Glaciology, Vol. 9, No. 55, p. 150–53.Google Scholar
Philberth, K., and Federer, B. 1971. On the temperature profile and the age profile in the central part of cold ice sheets. Journal of Glaciology, Vol. 10, No. 58, p. 3–14.Google Scholar
Radok, U. 1959. Temperatures in polar ice caps. Nature, Vol. 184, No. 4692, p. 1056–57.Google Scholar
Radok, U-, and others, 1970. Steady-state temperature profiles in ice sheets, by Radok, U., Jenssen, D. and Budd, W. [F.|. [Union Géodésique et Géophysique Internationale. Association Internationale d'Hydrologie Scientifique.] [International Council of Scientific Unions. Scientific Committee on Antarctic Research. International Association of Scientific Hydrology. Commission of Snow and Ice.] International Symposium on Antarctic Glaciological Exploration (ISAGE), Hanover, New Hampshire, U.S.A., 3–7 September 1963, p. 151–65.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. 1970. Stability of ice sheets as deduced from deep temperature gradients. [Union Géodésique et Géophysique Internationale. Association Internationale d'Hydrologie Scientifique.] [International Council of Scientific Unions. Scientific Committee on Antarctic Research. International Association of Scientific Hydrology. Commission of Snow and Ice.] International Symposium on Antarctic Glaciological Exploration (ISAGE), Hanover, New Hampshire, U.S.A., 3–7 September 1963, p. 141–51.Google Scholar
Weertman, J. 1968. Comparison between measured and theoretical temperature profiles of the Camp Century, Greenland, borehole. Journal of Geophysical Research, Vol. 73, No. 8, p. 2691–700.Google Scholar