Journal of Earth and Environmental Sciences (ISSN: 2577-0640)

Article / review article

"Brief Historical Tour of Glacier Ice on Earth and its Role in Climate Dynamics1 "

Kolumban Hutter1*, Dietmar Gross1

Department of Mechanics, Darmstadt University of Technology, Darmstadt, Germany

*Corresponding author: Kolumban Hutter, Department of Mechanics, Darmstadt University of Technology, Darmstadt, Germany. Email:

Received Date: 18 June, 2019; Accepted Date: 17 July, 2019; Published Date: 25 July, 2019


The science of the physics of ice crystals started almost 500 years ago with Cardano, followed by Kepler, Hooke, Dalton, the two Braggs (father and son) and Pauling. Here, we are concerned with naturally formed ice of glaciers, ice sheets and ice shelves. We do not discuss lake, river, sea and atmospheric ices, even though they equally play a role in today’s General Circulation Models (GCMs) for the Earth.

The study of glacier ice started with the observation of the behavior of the dynamics of Alpine glaciers, specifically that they are not rigid, but moving bodies that deform. A ladder, left in 1788 at the icefall of the Col de Géant by de Saussure was 44 years later found at the three-glacier merge of the Mer de Glace, corresponding to a mean velocity of 375 feet/year. Hugi and Agassiz measured the motion of rocks on the middle moraine of the Unteraar-gletscher and found similar values. However, at early times their location of origin was an enigma; so they were named erratic. Forbes and Tyndall found by position measurements with theodolites that velocities of glacier surface objects are (i) considerably larger in summer than in winter; (ii) larger on the glacier surface than at depth and (iii) larger in the middle of the glacier-width than close to the boundaries, as reported by Helmholtz (1865), [77]. On this basis Rendu and Forbes were the first to identify similarities of glacier flows with streams of very viscous fluids. It lasted until the 1950s when physicists had postulated the constitutive equation for isotropic ice as a non-Newtonian power law fluid (Nye 1952), [134], that was experimentally verified by Glen (1952, 1953), [55,56], and Steinemann (1954, 1956, 1958), [157-159]. Orowan (1949), [136], likely prompted the British contributions. This law is now known as Glen’s flow law. Various second grade fluid alternatives (e.g. Man and others (1985, 1987, 1992 2010), [116-118], provide extensions to capture primary and secondary creep. Comparing results of creep tests from distinct experimentalists disclose unexplained disparities of stress relations.

All large-scale mass Initial Boundary Value Problems (IBVP)s are viscous, heat conducting incompressible Stokes problems, using the Glen-Steinemann flow law, which may involve cold ice regions with the ice temperature below melting, separated (close to the basal boundary) by a Class-I mixture of temperate ice, i.e. ice with water inclusions. The dynamics of the cold-temperate transition surface between the two regions is governed by a (simplified) Clausius-Clapeyron equation (Fowler & Larson (1978, 1980), [49-50] Hutter (1980, 1983), [84-86], Hutter, et al. (1981, 1988), [82, 88], [Blatter & Hutter (1991), [15]).

As for boundary conditions at the free surface, snow accumulation and radiative heat must be parameterized as functions of position and time or these quantities must be taken over from concurrent large-scale flow models of the atmosphere. One must proceed in a similar way with the ice ocean interface on floating portions of ice shelves. Here, general ocean-circulation models are operative. Moreover, at the ice-shelf front, mass loss of ice by calving must be parameterized, which is significant for marine glaciers, ice shelves, and their dynamic response in such scenarios.


To begin, let us briefly present those articles that treat some aspects of the historical development of the physics of glacier flow. We collect articles in which the focus is devoted to the history of glacier and ice sheet flows approximately during the past 3 centuries but also including the past 30 years. Our focus is the Earth bound slow movement of ice masses and how this creeping flow was, in earlier centuries, intellectually developed and principally rationalized, eventually as a nonlinear viscous heat-conducting fluid. Of similar significance as the early material science properties of ice was the recognition of the concept of ice age(s). For the 19th century scientists, the enigmatic appearance of erratic boulders led to the recognition of the concept of ice age(s). This concept, first spelled out around 1830 by the German poet Goethe, see Cameron, D. (1965), [24], was only formally rationalized as the planetological cause of the ice ages that is coupled with Earth internal processes affecting the time evolution of the Earth’s climate.

The most significant memoir, in which the physics of glaciers is reviewed, is Garry Clarke’s (1987) ‘A short history of scientific investigations on glaciers’, [27]. The author gives a tour of ‘all’ glacier related processes, starting in the first quarter of the 19th century, but primarily focusing on the years ~1940-1987 and describing mostly research performed in the past 50 years related to the physics of glaciers. This work is much geared to accumulate the varied knowledge separated by subjects; it is excellently summarized. Its 19th century coverage is similar to ours, but with a broader coverage of subtopics and has perhaps less of a natural-philosophical emphasis than ours.

The paper ‘Decent of glaciers: some early speculations on glacier flow and ice physics by John Walker and E.D. Waddington (1988), [171], concentrates on the early work of the motion of glaciers as observed in the 19th century. They cover the literature of the period from ~1820 to 1870, notably listing the memoirs of Agassiz (1837, 1840, 1842), [3-5] de Charpentier (1841), [26] Forbes (1841, 1842, 1843), [ 45-47], covering essentially what was stated by Clarke (1987), [27].

The memoir ‘Life, death and afterlife of the extrusion flow theory’ by E.D. Waddington (2010), [170], is devoted to a specialized subject, namely extrusion, which was a fallacious model of ice flow. In our manuscript, we do not deal with such phenomena, i.e. processes that failed to be recognized as realistic.

Two additional memoirs in the Journal of Glaciology are devoted to the hydraulic component of Alpine glacier flow. Gwenn Flowers (2010), [43], devotes her paper to such related works by Almut Iken (1933-2018), [93-95], Similarly, Christina Hulbe, et al. (2010), [81], devoted their paper to Women in glaciology. Both these articles are at best tangibly related to our present concern of glacial flows.

Thus, it transpires that our focus is, first, the birth of the description of glacier ice flows as the motion of a heat-conducting fluid and, second, the presentation of the necessity of extending the flow law of ice from its nascent postulation beyond the Glen-Steinemann flow law as a general nonlinear viscous fluid. We stress circumstances where we believe that present research activities do not strictly enough draw the inferences that should be drawn on the basis of the findings.

To understand how the large ice masses and glaciers on Earth contribute to the climate scenarios, one must first securely know to what extent the climate driving of the solar system interacted, and will interact, with the Earth-interior processes, and how these affect the various components: the atmosphere, the solid Earth, the oceans and hydrosphere, as snow, ice, water and vapor, not to forget the biosphere, which also has a significant anthropogenic component. Unquestionably, there is also an anthropogenic component, as demonstrated by industrial activities since approximately the year 1780 and expressed by the equivalent  concentration in the atmosphere, about 280 ppm (parts per million) in the year 1780 and 415 ppm on May 19. 2019 as measured by the Mauna Loa Observatory in Hawaii.

This brief account of the most significant components of the Earth’s climate system only touches those parts of the entire picture, as they chiefly developed in the last 90 years or so, beginning with Milankowitch’s thesis (1930), [124]. It consists of three distinct but concurrent processes: eccentricity, obliquity, and precession. Eccentricity, the elliptical cycle of variation in shape of the Earth’s orbit around the Sun is about a 100.000-year cycle. Obliquity is the cycle of axial tilt of the Earth’s rotation axis toward or away from the Sun and varies from about 22 to 24.5 degrees, usually taken as 23.5 degrees with about a 26.000-year cycle. As these parameters change, so does the amount of sunlight that hits different latitudes on the Earth; for a NASA illustration of these, see Gustovich (2018), [74].

As a scientific focus, climate dynamics on Earth and its anthropogenic coupling came last. In the period from the mid-17th to the mid-19th century, glaciers were first thought to be rigid objects. That they somehow move was mainly recognized in the 18th century, in parts by conjecturing that the volume expansion of water in the freezing process would push the ice downward in their valleys. Later in the 18th century this belief was replaced by the hypothesis that glaciers would slide over their beds. Moreover, edgy rocks and dirt on the glacier surface, which were seen to move, brought the impetus for the assumption that glaciers move and deform. Careful geodetic measurement with theodolites in the 18th and 19th century, finally, disclosed the creeping character of this motion as deformable like a dough. A further enigma was still the inexplicable presence of erratic boulders within an environment of unmatchable petrography. The solution of this mystery led finally to the hypothesis of earlier ice ages, a fundamental topic of climate dynamics, postulated in the 18th and 19th century and rationally explained in the 20th century.

To be able to describe the deforming motion of glaciers calls for material science of the ice. The determination of the material properties of polycrystalline ice as a nonlinear viscous, heat-conducting body - at early times often denoted as plastic - is vital for the quantification of its melting processes. This concerns the substance and its appearances as ice, water and vapor. These are described by the thermodynamic principles (first and second law). It so happens that the Earth’s climate conditions through the past few million years have driven the phase changes of  so sensibly that these processes reign not only the various forms of life on Earth, but also its ice coverage, the evolution of the sea level due to natural external variations as well as anthropogenic Earth internal climate changes. These latter will be the cause of roughly 80 cm sea level rise until the end of the 21st century, according to the IPCC (Intergovernmental Panel on Climate Change) Report.

This brief overview of how the phenomena ‘glaciers and ice sheets’ triggered the curiosity of mankind to eventually become a respected part of the natural sciences of even high significance for the humans’ living on this planet. This picture must be complemented by the material sciences. These ought to describe the  phases as objects of mathematical physics. Only with the addition of the description of the material behavior of ice, the ice processes on Earth can rationally be described; only with these, quantitative predictions can be made. This phase begins with the recognition that isotropic polycrystalline ice under slow creep can be modeled as a power law heat-conducting fluid. It goes back to experiments conducted in the 1950s by John Glen (1952, 1953), [55, 56], and S. Steinemann (1954, 1956, 1958), [157-159]. This law reads for simple shear

where is the rate of strain (rate of the shear angle), T is the temperature and τ is the shear stress. The scalar function f is generally written as the product

A(T) is called rate factor, and n (generally chosen as n = 3) determines the exponent of the power law. In the

glaciological community, the law (1) is called Glen’s flow law (we call it the Glen-Steinemann flow law, because the law was approximately simultaneously and independently proposed in the 50s of the last century by both researchers in the UK and Switzerland, respectively [1]). However, in material science, depending on the field of science, it is attributed to different distinguished scientists, in rheology to Graham (1850), Guoy (1910). Ostwald-de Waele [2] brought it forward in (1929) in the context of colloidal fluids and Reiner already before (1929) in rheological problems. In metallurgy, it is called Norton’s flow law (1929), and in plasticity theory Orowan had used it. In Glaciology John Nye (1953), [134], had demonstrated to the applied glaciologists, how it could be ‘derived’ from simple material postulates of viscous fluid behavior. It is likely that Orowan, who had also worked on constitutive behavior of plastic materials at that time, had directed Glen and Nye into such a nonlinear law; all three were then working at the Cavendish Laboratory at Cambridge University, UK, where Glen and Nye were junior collaborators.

So, the Glen-Steinemann flow law has precursors in not too distant related fields. We scrutinize its performance in comparison to laboratory and in-situ experiments and will also report on a number of attempts to the parameterization of the proposed laws. The likely best form of this viscous thermo-mechanical material description will then at last be incorporated in full-scale computations of glaciers and large ice masses in regional and completely Earth-embracing physical-mathematical-numerical models. This is in order to try to understand the climate relevant processes on our Globe for computational reproduction of the past climate variation due to the thermo-mechanical input data from extra-Earth processes, as well as to predict of future reaction.

Early Deformation and Sliding Models of Alpine Glaciers

The proposals for the enigmatic large ice masses by Johann Jakob Scheuchzer [3] and others in the 18th century can only be understood, if we accept that prior to the mid-19th century no quantification of constitutive relations was available to the material scientists interested in the physics of large ice masses. Nevertheless, glaciers as natural phenomena at high altitudes of mountainous regions were attractive, if simply for their mystic appearance. Prior to the 17th century, not even the deformability of glacier ice was recognized, certainly not admitted. Glaciers were postulated to be rigid objects, see e.g., Moraltus, Johann (1669), [105]. Only valley inhabitants close to them accepted some deformability, because e.g. they observed their snouts to advance or retreat, but could hardly rationalize a cause, and if they had been able to, they were too isolated in their valleys to get in contact with musing intellectual hikers. These modest people often knew it better.

So, with this background in the year 1705, Johann Jakob Scheuchzer visited Swiss Glaciers and proposed a theory on their motion. He knew from physics that water is expanding in the freezing process to ice, and that the ‘Force’ of expansion is so large that cartridge rounds which are filled with water that freezes, are blasted into pieces. Scheuchzer assumed that the water in glacier fissures and crevasses that freezes would extend with such excessive power that its force will unquestionably push the glacier downward. This concept, later often called ‘Dilatation Theory’, was also adopted by Jean de Charpentier [4], his brother Toussaint de Charpentier [5] and Louis Agassiz [6] and others. Basic thought of this concept was the belief that glaciers are permanent storehouses of coldness, capable to freeze all water that percolates through them. Interesting and strange in this model is that the volume expansion in the freezing process of ice is the only cause of the downward motion of the ice of glaciers. As a question at hindsight one might ask whether these scientists had ever looked at the water outlets of Alpine glacier termini with their sub-glacier rivers and generally substantial water discharge that is often rather large, certainly not zero during summer, as they required by their early deformation postulate.

About in the year 1760, Altmann [7] & Gruner [8] brought forward their opinion that glaciers would move by means of sliding along their beds. Both published their contributions independently, Altmann in 1751 and Gruner in 1760. Almost 40 years later, Horace Bénédict de Saussure revived this concept, which became the ‘de Saussure Theory’, even though, according to Tyndall (1878, p 185), he was never chiefly involved in it.

Cautious Birth of the ‘Plasticity Theory’

The 18th century also brought the first ideas about the deformability of ice sitting on solid ground. The concept of sliding should simply have brought the science mountaineers to the postulation of the deformability of glacier ice. Indeed, glaciers generally mostly move on non-planar beds; this fact, together with the existing forward motion requires deformability to maintain the sliding hypothesis. Yet, none of the above mentioned scientists attributes the terms ‘viscosity’ and ‘kneadability’ to characterize the deformability of the ice, even though, according to Tyndall (1878, p. 185), [164], the appearances of many glaciers suggest these terminologies, if it were not so contradictory to any daily experience with ice. Needless to say that the concept of viscosity was known since Newton had introduced it in the Principia, (first edition 1687, third edition 1726), [132].

In spite of this, these kinds of plastic concepts found their defenders. In the year 1773, André César Bordier from Geneva published a small booklet ‘Voyage pittoresque aux Glacieres de Savoyes’, Genève, 1773 (‘Pituresque journey into the Glaciers of Savoy’). He advocates for an over-all view of the ice mass which moves as a whole from high to lower altitudes in a manner as seen with other fluids. ‘Do not let us view the ice as an immobile and stiff material, similar to mollified wax, which, to a certain degree, is flexible and extensible’. This is likely the first occasion where the kneadability of glacier ice is mentioned.

However, according to Tyndall, André César Bordier’s [9] concepts were unheard in the natural scientific community in the 1770s. They were reborn more than 60 years later by a successor of better scientific training. This person was a catholic priest and later bishop of Annecy, Louis Rendu [10]. In the year 1840 he submitted ‘Théorie des Glaciers de la Savoie’, [146], to the Royal Academy of Savoy. Tyndall, was fascinated by Rendu’s writings and states that it is not known, whether he ever saw Bordier’s works, probably not, because he never refers to it. Tyndall quotes a few of Rendu’s statements to illustrate the preciseness of his expressions:

· Between Mer de Glace and a river, there exists such a perfect similarity that it is impossible to find a circumstance in a glacier, which would not equally occur in a river.

·       In flows of water, the movement is not uniform, neither with regard to the width nor to the depth.

·     The friction at the bed and at the sides, paired with the effect of obstructions, makes the motion to differ from position to position; only toward the middle of the free surface, one reaches the full motion.

·    There exists a large set of facts, which seem to force us to believe that glacier ice possesses some sort of extensibility, which allows it to adjust to the local circumstances, to thin, to swell and to contract as if it were a soft dough.

To corroborate these inferences, Rendu requested careful measurements of the glacier motion; he did not perform these himself, but asked his mountain guide to observe the motion of a particular boulder at the glacier side, which, during 5 years moved 40 feet/year. Other boulder position measurements corroborated the motion. This is a first indication to search for motion patterns of glaciers, which we shall subject to a detailed analysis later on.

To summarize

First descriptions of the motion and deformation of glaciers are due to Bordier in 1773, but the flexibility and extensibility with a characterization of viscous or plastic behavior was recognized by Rendu as late as 1840, [143], and afterwards, but could not yet be phrased in terms of a constitutive relation such as the Glen-Steinemann law.

Erratic Boulders

Since the 19th century, the edged boulders, apparently incoherently distributed in mountainous territory are named glacial erratic, which explicitly suggests for us that they were once moved by glaciers and deposited at their present positions when the carrier ice retreated. These rocks suggested to the mountaineers misleading or inexplicable or false behaving, briefly err-behaving. The causal connection had first to be recognized, and it influenced at least subconsciously the early understanding of the existence of the Ice Ages. According to Krüger (2009), [105], the oldest written statement on erratic boulders dates from 1301 and is reported by ‘Henricus dictus von dem Steine’ [11]. The first assured written statement was found in Johannes Guler’s ‘Raetia’ [12] In his script Krüger (2009), [105] tells that Guler von Wyneck mentions a colossal rock in Veltlin, in the Italian Alps, for which he could not see, where it may have broken off. With the growing development of modern geology in the mid-18th century the general interest in the erratic boulders grew. In 1727/28 Moritz-Anton Capeller [13], physician of the town Lucerne, recognized the alpine origin of the erratic boulders, deposited in the mountainous foreland. However, he only published his recognition in 1767 in his ‘Pilatii montis historia’, shortly before his death.

Initially, the contemporary scientists as e.g. the notary Abraham Schellhammer [14] (1675-1735) and Albrecht von Haller [15] both from Berne assigned the ultimate cause of the positional spreading of the erratic boulders to the biblical deluge or (later) to a number of exceptional floods. Tobias Krüger (2009), (2013), [106, 107], reports a number of exceptional and extreme interpretations, even by scientists of the status of Horace Bénédict de Saussure (1740-1799) and many others.

One interpretation by the Prussian geologist Christian Leopold von Buch [16], why an erratic boulder found in the Jura Eastern slopes, but of likely origin in the upper Valais is as follows: He states that originally far in the past, today’s Rhone Valley was blocked at Saint Maurice by a rocky Rigel (bedrock bar) between the Dents du Midi and the Dents des Morcles. Behind, a gigantic lake was formed up to the mountain peaks. At the breakdown of this natural Rigel, the water masses were set free with hardly imaginable power; they threw huge rocks as far as to the Jura (Krüger (2009), [106], notes that this information was only published after von Buch’s death in 1867). We will refrain here from reporting further erroneous propositions of the displaced positions of erratic boulders. Interesting, however, is that in the mid-18th century, more than 100 years earlier, in 1742, the engineer and geographer Pierre Martel (1701-1767) from Geneva, reports of large displaced rocks found during his journey to the Valley of Chamonix. Its inhabitants had apparently explained to him that the ‘Glacier du Dois’ once carried these rock pieces down the valley. Krüger (2009), [106], concludes: For this reason, at the present state of knowledge [meaning 2009], the inhabitants of the Savoy Alps were the first, who established a causal relation between glaciers and position-alien rocks.

In the last decade of the 18th century, James Hutton [17] devoted his working time to geology and authored the treatise ‘Theory of the Earth’, originally planned as four volumes. He adored Saussure’s imagination in his theory of erratic boulders, but rejected Saussure’s concept that erratic boulders reached their places of final position before the erosion had formed today’s valleys. He writes ‘there would then have been immense valleys of ice sliding down in all directions towards the lower country and carrying large blocks of granite to a great distance where they would be an object of admiration after ages, conjecturing from whence or how they came.’ As primary cause (be aware, existence of earlier ice ages was not known at this time) for the extended glaciers of the Alps, Hutton assumed that the latter had been substantially larger at earlier times. Therefore, it then must have been colder there, which favored an intensified glaciation [because atmospheric temperature decreases with height; this was known at that time]. Damaging of Hutton’s concept was that it required a tremendous erosion rate to excavate the suspected high-level plane and to form the valleys from the initially much higher mountains. This would last an excessively long time. Positive is that he is one of the first, who suggested a connection between the phenomenon of erratic boulders and the spacious earlier glaciation.

After Hutton’s death, the theologian and mathematician John Playfair continued his work [18]. He summarized this work in his book of the year 1802, [142]. His message is the same as that of Hutton. Krüger, [106], quotes him as follows:

‘For the moving of large masses of rock, the most powerful engines which nature employs are without doubt the glaciers, those lakes or rivers of ice which are formed in the highest valleys of the Alps, and other mountains of first order. These great masses are in perpetual motion, undermined by the influx of heat from the Earth, and impelled down the declivities on which they rest by their own enormous weight, together with that of the innumerable fragments of rock, which are loaded. These fragments they gradually transport to their utmost boundaries, where a formidable wall ascertains the magnitude, and attests the force, of the great engine by which it was erected’, from ‘John Play fair, Illustration of the Huttonian theory of the Earth’, Edinburgh 1802, p 388ff §348, [142].

Similar observations and interpretations were also made by natural scientists from other countries. In Bavaria the geologist, mineralogist and physicist Mathias von Flurl [19] recognized the position-alien character of the erratic boulders in Bavaria. In the year 1809, the Bavarian physician and astronomer Franz von Paula Gruithuisen [20] devoted some of his research to erratic boulders in the Alps. He knew of the significance of such boulders, and rejected the hypothesis of large floods to be the cause, but kept the transport hypothesis of the apparently unmotivated spreading by water; so, the heritage of the deluge as the principle cause was still behind his thinking. Krüger (2009), [104], cites von Paula Gruithuisen with arguments which today must be interpreted as abstruse geological and hydrological thinking. He, further muses that von Paula Gruithuisen did not seem to know that glaciers do slowly flow down their valleys. They must be immobile like frozen lakes, and a geological force must exist, which pushes the glaciers, including their edgy rocks, upwards to transport the latter large distances. Evidently, such a concept seems today completely unrealistic. The biblical deluge can best serve in this situation as a ‘deus ex machina’ to resolve this intellectual dilemma.

About at the same time, during the early 19th century, Jean-Pierre Perraudin [21] and Siméon Gilliéron [22] came to the conclusion that “The glaciers of our mountain chains […] had at earlier times a considerably larger extent than today [sic: at Perraudin’s and Gilliéron’s time]. Our entire valley was up to a high altitude above the Dranse, a river in the Valley de Bagnes, occupied by a single huge glacier that extended far down to Martigny, as the rock boulder evidenced, which were found in the vicinity of this town and which are too big to ever have been transported by water” (translated form Krüger (2009), [106], by KH). Apparently, Perraudin’s thesis of 1815 was only based upon positions where erratic boulders were found. This was an indirect recognition that positions of erratic boulders were determined by the earlier glaciation of the valley area, in which they were found.

To summarize:

In the 18th century the contemporary scientists had abandoned the biblical deluge hypothesis in favor of a series of flood events as the causes of the enigmatic positioning of the erratic boulders. Only after the 1740s, first suppositions of the transport of erratic boulders by glaciers appeared, but the principal postulated causes were largely based upon rather unrealistic theoretical concepts, which even were contradictory. The question, how glaciers would move - even the problem whether they would move at all - remained controversial. This remained so in the 19th century. Reflected on from the 21st century, this cannot be a surprise as the climate variations were not yet understood and a continuum mechanical formulation for material behavior had not yet been established.

Systematic Measurement of the Deformation and Motion of Glaciers

Erratic boulders were of help in searching for the deformation and motion of glaciers as a whole. Direct observation of objects on the glacial surface and at certain depths below it, in crevasses and fissures, on the other hand, brought a breakthrough. The detailed observations and measurements were done by James David Forbes [23,44-46], and John Tyndall [24] and are expertly described by Tyndall (1860, 1878) [163] and -- in lecture notes ‘Eis und Gletscher’ (‘Ice and Glaciers’), held by Herrmann von Helmholt [77], in February 1865 in Frankfurt am Main and Heidelberg. On pp.108, 109 of these notes he writes. ‘We have so far compared glaciers according to their manifestation with streams; this similarity is, however, not just an external one; on the contrary, the ice of glaciers moves in fact forward, similar to water in a stream, only slower […]. Since, namely, the ice at its lower end is incessantly lessened by melting, it would soon have disappeared, if new mass would not continuously move forward from above, which would by snowfalls in the firn fields repeatedly be renewed’ (translation KH). Von Helmholtz reports on a number of such convincing situations which we itemize below.

Isolated rock boulders on the surface of a glacier can be used as identifiers of material surface points of glaciers. This observational procedure can be used for any object positioned on the ice and left there, but repeatedly observed for its position.

· In the year 1827 Hugi [25] erected a hut on the middle moraine of the Unteraargletscher, to perform there some measurements. The position of this hut was determined by him and by Agassiz 14 years later, in the year 1841; it stood then 4884 feet (~1602m) down-glacier. A somewhat smaller velocity was found by Agassiz for his own hut of the same glacier, [3-5].

· Measuring positions of the surface boulders by theodolites, made by Forbes and Tyndall, from day to day disclosed for the Mer de Glace that during summer the velocities were of the order of 20-35 inch/day (54.7-96.7 cm/day); during winter, they were about half as large; at the ice surface they were larger than at depth, and at the sides they were considerably smaller than in-between’. (von Helmholtz 1865, [77], pp 110-11, translation by (KH))

· The observation of the kind on motion of glaciers allows also an interpretation, in which orientation crevasses must be formed; since, namely, the forward motion of the ice at different positions is not the same, the relative distances of these mass points change with time. Now, because the ice between any two such points cannot be arbitrarily stretched, it will be fractured and form fissures or crevasses’ (Helmholtz 1865, [77], p 112, translation KH). Von Helmholtz, quoting Tyndall’s conviction that the glacier is hardly able to resist tension but rather ruptures under the influence of such tensions, illustrates this in (Figure 2).

· Both, Tyndall and von Helmholtz, see (Figure 3), beautifully explain how moraines and dust stripes illustrate the deforming motion of the glacier system ‘Mer de Glace’.

The three contributing ice-arms are ‘Glacier du Géant’, ‘Glacier de Léchaut’ and ‘Glacier du Talèfre’, which unite as ‘Mer de Glace’. Side moraines transporting boulders, fallen from the side flanks of the mountains onto the moving ice, in the three upper glacier arms become in the ‘Mer de Glace’ interior moraines, except for the outermost two side moraines to the far orientation crevasses must be formed; since, namely, the forward motion of the ice at the far left and far right. The dust stripes in the ‘Glacier du Géant’ and further down in the ‘Mer de Glace’ all come from the icefall fed by the ‘Col de Géant’ (g in panel (b) of (Figure 3)). The material falling down this icefall -- varying from winter to summer -- forms the ice-dust mingling at the base of the icefall, from which the stripes emerge as displayed in (Figure 3), panels (a) and (b). These stripes have in the downstream stretch of the ‘Glacier de Géant’ a curved appearance across the glacier width that amplifies further down in the ‘Mer de Glace’. These bands of dust and the spread of rocks and stones appear as alternating grey and whitish stripes, likely indicating annual repetitions of the ice, as first noted by Forbes. All this is impressive manifestation that velocities at the sides are smaller when compared with the free surface motion in the interior, as so described by Tyndall in the year 1842, [163].

The dust stripes, indicating the displacement distribution of the free surface are in the ‘Mer de Glace’, limited to the Western part of the glacier; the Eastern part is covered by interior moraines but no curved dusty stripes are seen - at most an onset of spreading that grows in the down-glacier direction [26]. Corroboration of this interpretation is given in detail by Tyndall (1873, pp 93-123). He and his collaborators demonstrated this for the ‘Mer de Glace’, the Grindelwald-, Aletsch- and Morteratsch-Gletscher. The method was to insert wooden sticks at several points along an initially straight line into the snowy surface across the glacier (generally 6 to 10 per transect and to determine by theodolite their positions at consecutive times (generally one to a few days apart), see [Figures 4,5].

This allowed evaluating the travelled distance for each wooden stick at consecutive times. This meticulous procedure was the verification of the bold assertion that ‘the ice of a glacier flows slowly, and similarly to a stream of a very viscous substance, as e.g. honey, coal-tar, or a thick tonic pulp’. The ice does not simply slide over the bed, just as a rigid body does, which slides down a slope; it rather bends and displaces itself in itself and, even though it glides over the base of the valley, the parts, which are in contact with the bottom surface and the side walls of the valley, are effectively slowed down, caused by the significant friction. By contrast, the middle of the free surface of the glacier, which is farthest from the bottom and the walls of the valley, move faster. Louis Rendu and James David Forbes were amongst the first to emphasize the similarity of glacier flows with the flow of a fluid substance’ (Helmholtz, 1865, [77], p 116, translation KH).

Thus, by the 1860s the scientists focusing on the physics of the glaciers had acquired convincing knowledge that glaciers were slowly moving and deforming similarly as water in rivers, creeping soil down mountain slopes, hot lava as rivers in volcanic eruptions, honey or polymeric fluid substances. However, a mathematical formulation of such motions as continuous media was not available at that time. Neither was there any mention that inertial terms in Newton’s law in river motions have a sizeable influence, whilst they can safely be ignored in the determination of glacier motions.

To summarize

Systematic measurements of the motion of the ice on the glacier surface at selected points on the surface, at its middle line, close to its side boundaries in fissures and crevasses below the surface were consistently done by Forbes and Tyndall using theodolitic positioning at consecutive days. The evaluation of these measurements led to the interpretation that glacier motion fields are analogous to the motion of terrestrial surface river flows, a fact explicitly spelled out by Rendu in the first half of the 19th century, see [143].

Birth of the Concept of an Ice Age

According to Rowlinson [27] (1971), [145] the years from 1840 to 1870 were amongst the most fruitful ones in the history of physics. As for the solution of the problem of th position of the erratic boulders the systematic measurements of the motion and deformation of glaciers primarily done by Agassiz, Forbes and Tyndall were important as was Rendu’s ‘Théeorie des Gaciers de la Savoie’ (1840), [143]. Forbes had formulated a model of viscous deformation, which was criticized, because it conflicted with the sliding model, but Forbes successfully showed that the sliding model was in conflict with the (experimental) facts, which he had discovered (Rowlinson 1971, [130]). By 1851 Forbes’ measurements were summarized in his article of 1855 for the eighth edition of the encyclopedia Britannica. Abbreviated as evident, this statement reads

“Each portion of a glacier moves…in a continuous manner…The ice in the middle part of the glacier moves much faster than that near the sides or banks; also the surface moves faster than the bottom…. The variation of velocity (as in a river) is most rapid near the sides…. The glacier, like a stream, has its pools and rapids. Where it is embayed by rocks it accumulates, its declivity decreases, and its velocity at the same time…and the increased temperature of the air favor the motion of the ice…. The velocity does not, however, descend to nothing even in the depth of winter…. [These] circumstances of motion … appear to be reconcilable with the assumption of what may be called the Viscous or Plastic Theory of glacier motion, and with that alone … the veined or ribbon structure of the ice is the result of internal forces, by which one portion of the ice is dragged past another…. The veined structure is unquestionably the result of the struggle between the rigidity of the ice and the quasi-fluid character of the motion impressed upon it. That it is so evident not only from the direction of the laminae, but from their becoming distinct exactly in proportion to their nearness to the point where the bruise is necessarily strongest”.

Rowlinson’s paper is in large parts a detailed historical account of the dispute and disagreement over how the motion of glaciers had to be physically interpreted. “Tyndall’s criticism centered on the word ‘viscous’, which Forbes never defined clearly. He apparently understood it to mean resistance to shear stress, but used the terms viscous, plastic and semi-fluid indifferently…Tyndall interpreted it differently.” By viscosity he understands that property of a semi-fluid body, which permits of its being drawn out when subjected to a force of tension. He denied that ice could be permanently stretched without braking and so dismissed Forbes’ viscosity as apparent, not real. Instead he proposed his theory of ‘fracture and regelation’…” (Rowlinson, 1971, [145]). He describes in detail this somewhat bitter fight that lasted for years; it also bore elements of jealousy, because Forbes feared that Tyndall (the younger of the two) might get the Copley Medal of the Royal Society in 1859. William Thomson 1824-1907 (later (1892) Lord Kelvin), apparently was also involved in this maneuver.

Thermodynamicists of the 19th century played equally a role in the process of creating a realistic physical basis for the motion of glaciers. At that time, Julius Robert Mayer (1814 - 1878) [119] (a physician in Heilbronn, Germany) had formulated the First Law of Thermodynamics - then called the mechanical law of heat – and Rudolf Clausius had introduced the Second Law of Thermodynamics [28,29], as did Lord Kelvin [161], who introduced in 1848, at the age of 24, the absolute temperature. This was known already earlier, since ice floats on water. So one knew of the expansion of water in the freezing process to ice, and it is, thus, not so surprising that Johann Jakob Scheuchzer and later Jean de Charpentier, [26] and Louis Agassiz, [4,5] at least for a while assumed that only freezing of intra-glacial water would cause the downward motion of glacier ice.

Nothing of this provides a hint as to the concept of ice ages; neither Agassiz nor de Charpentier nor any other of the glacier scientists of that time proposed the concept of very cold climates in the past. Today the creation of the idea ‘ice age’ is largely still attributed to Louis Agassiz, but this is erroneous: Surprisingly, it is Johann Wolfgang von Goethe (1749-1832), one of the world’s greatest poets, who was equally a respected scientist. Dorothy Cameron (1965), [24], in Vol. 5 of the Journal of Glaciology tells the story. Goethe as a scientist was the first to attribute the transport of erratic blocks to glaciers; he spoke out that an ice sheet covered Northern Germany; furthermore, he is given the credit to first having believed in an ice age. Cameron (1965), [24], states that Louis Agassiz (1837) and Jean de Charpentier (1841) were correct when they assigned to Goethe the credit for being the first to have conceived the concept [28] of the ice sheet. There are very few written attributes to this fact. We have only found Reinhard Hederich (1898), [76], being mentioned in a footnote of Alfred Krehbiel’s dissertation of (1902), [104]. The comment is: ‘…neither should Goethe be forgotten, who built himself noteworthy correct views of glaciers and erratic bolders…’ Similarly, also Forbes (who was not directly addressing erratic boulders or the concept of ice ages) was well aware of Goethe’s scientific significance. On the title page of his ‘Travels through the Alpes’ he cites Geothe’s sentence: ‘Sage mir was du an diesen kalten und starren Liebhabereyen gefunden hast’ (tell me what you believe to have found in these cold, rigid fondnesses).

The reference to Goethe is also missing in the reviews by Clarke (1987), [27], Walker & Waddington (1988), [171] and Waddington (2010), [170].

Juvenile Struggle with a Creep Law for Ice

As explained in the last section, it was in the 19th century that natural scientists recognized that glaciers are slowly moving ice masses. Helmholtz (1865), [77] and Tyndall (1860, 1878), [163,164], summarized the knowledge as of ~1870; they beautifully describe the slow birth of this new understanding of glacier flows as movements of a very viscous material. For the recognition of this property of deformability a large number of scientists were involved, besides Tyndall, at least Forbes and Rendu ought to be mentioned. At the mid 19-hundreds, further progress in developing a physically based theory for the motion of large ice masses was hampered by two facts,

   ·  The nonexistence of a formal observer-invariant material theory, coupled with the basic physical laws,

          ·  The primitive state of experimental techniques for the determination of the constitutive relations (here primarily for the Cauchy stress tensor).

The second half of the 19th and the early 20th century were needed to develop these two physical, mathematical and engineering-type specialties. As a prelude to the glaciological activities starting shortly before 1950, it is important to note that the material sciences of continuous media -- fluids and solids -- as they developed in the second quarter of the 20th century were active in a nascent field, devoted to the description of the mainly continuous deformations of bodies to external driving elements such as forces, temperature, etc. This new science was coined rheology and became quickly fashionable among chemists, material scientists, applied mathematicians and physicists. The important features in the field of material creep to be solved were how e.g. the stress tensor can be expressed by deformation measures, such as strain and/or rate of strain tensors, temperature, etc. in a materially objective manner, i.e. observer invariant form. Markus Reiner (1886 – 1976), an Austrian-Israeli (civil) engineer, was an early protagonist of the description of the creep behavior of fluids and developed prior to 1929 together with his scientific associate Ms. R. Riwlin the Reiner-Riwlin constitutive relation for fluid creep. Expressed in a nutshell: Let t be the stress tensor and t´ its deviator and the strain rate tensor [29]. These are symmetric 3x3-tensors, and in a volume-preserving (i.e. incompressible) material, D is a deviator (its trace or first invariant vanishes). So, a possible constitutive relation must be expressible as, where the function f, when evaluated in terms of and T, must equally be a deviator. The explicit form of f, satisfying these requirements can be constructed with relatively simple methods of 3x3-matrix algebra and is given by

in which both sides of this relation are deviators,  are scalars given by

and   are the second and third invariants of  (while the first invariant  vanishes). The above expressions are in conformity with isotropy. We ask the reader to accept the above polynomial expression as the most general expression for  under the stated conditions. In ensuing developments, the term involving  will be called the collinear (or affine) term to D for obvious reasons. The second term, involving , may then be called the quadratic term.

  of this relation are deviators,  are scalars given by

The utmost majority of scientific works in glaciology has been done for the case that  and that  does not depend on . Investigations have been done in which , but it will be shown that for polycrystalline isotropic ice an extended dependence also involving  can consistently only be introduced when the quadratic term involving  is included.

The remainder of this historical review will be devoted to the test of adequacy of the Reiner-Riwlin fluid as a constitutive model for the creep of polycrystalline isotropic ice. Naturally, this historical episode only dates back to the last 60-70 years with important contributions done since ~1980. We do this in this historical article against the recommendation of two reviewers and the handling Scientific Editor of the J. Glaciology. Referee 1 of that version of this paper stated

“There is an important borderline between and ,

 and I [sic: the referee] believe that the paper tends to cross the border from history toward the side of  expressing a scientific opinion. I view this to be a point where the paper is no longer about history but is about new science being advocated…”

We as authors disagree with this request that the two parts ought to be separated. Quite contrary: Historical facts ought to be used for the consequences, which they imply. Trivially, likely all facts of our manuscript have happened in the past (with the credits given to others) and are therefore part of the history.

Applied glaciologists have almost exclusively employed the Reiner-Riwlin fluid model in its reduced form


in which  is a power law expression of , multiplied with a temperature dependent rate factor. For temperatures distant from the melting point this dependence is expressed by an Arrhenius-type relation, and close to the melting point generally constructed by optimizing a functional ansatz to observed temperature data. Polynomial fits to temperature data have also been done and show generally a better agreement with data points. However, when comparing optimally determined expressions from one experimental site to another (or from one author group to another) agreement is, more often than not, disastrous. Smith and Morland (1981), [130] have documented this. Despite this demonstration, practically all glaciologists not involved in the (experimental or theoretical) determination of the ‘flow law for ice’, apply the power law unquestioned as if it were the commandment of a glaciological deity. Nye’s paper of (1952), [135], may well serve as the trigger of such attitudes. In that paper, two basic assumptions for the strain-rate--stress relation were that (i) this relation did not depend upon the third invariant and (ii) the strain rate deviator is collinear to the stress deviator. Given the prerequisite that Orowan was a leading scientist at the Cavendish Laboratory in Cambridge at that time (see his paper of (1949), [136]) involved with plastic modeling of materials, and Nye and Glen were junior members there, it is but natural that Glen came up with a power law. Let us now illustrate, how the last 30-40 years were vital for our thesis that polycrystalline ice ought to be modeled as a full Reiner-Riwlin fluid or even as a more general fluid.

The scene began in the 1950s with the recognition that creep of ice in large ice masses must be treated as an isotropic nonlinear viscous or plastic body. Creep tests on polycrystalline ice were done first by Perutz (1950), [141]; Glen (1952, ‘53, ‘54, ’55, ‘58, ‘74), [55-61] and Steinemann (1954, ‘56, ‘58), [157-159], followed by Gold (1958), [61-63], Mellor & Smith (1967), [120], Mellor & Testa (1969a,b), [121,122], Jacka (1984), [97], with reviews given by Kuo (1972), [108], Shumskiy (1974), [154] and Mellor (1979), [123]. The commonly used stress systems are tension, compression and simple shear, although combined states of stress had also been looked at. In addition, the attitude of material scientists were, that experimental results obtained with such simple stress systems could be extrapolated to 2- and 3-dimensional stress systems.

Typical creep curves for polycrystalline ice under simple stress states (shear, tension/compression) are as shown in (Figure 6) with primary, secondary and tertiary creep, where

·  Primary creep is decelerating,

·  Secondary creep is steady

·   and tertiary creep is accelerating, merging into another steady state or leading to rupture.

Experiments by Steinemann (1958), [153], see (Figure 6a,b) show how creep curves look like in these ranges and (Figure 7a,b) shows analogous curves in doubly logarithmic representation. Since in these latter plots plotted

against ln( ), are not straight lines, these results provide corroboration for the fact that the power law behavior with constant exponent cannot describe the appropriate strain-rate relation [30]. Steinemann writes


Note, we write for the one-dimensional strain rate now  (in lieu of |D|) and  (in lieu of ) as this is the notation employed by Steinemann and Glen and used by engineers. Steinemann and Glen approximate these expressions as


requesting that this law is restricted to a finite stress range   , with different constant values of n in different stress ranges of lower  and upper .

The above representations (2) and (3) are restricted to uniaxial normal stress and simple shear, but the structure of these formulas was wished to be also applicable to 3-dimensional stress states. For density preserving materials one may write an objective strain-rate stress relation for polycrystalline isotropic ice as a Reiner-Riwlin fluid; i.e., if t is the Cauchy stress and its deviator, the stretching-stress relation reads


in which both sides of this relation are deviators and


 are second and third invariants of . This relation was made plausible above when we gave a wordy explanation of the Reiner-Riwlin [31] constitutive relation. Glacial theorists must in principle have known the relation in the 1950s. John Nye (1952, 1953), [121,122], took this to postulate that,

                 (i) the third invariant, , has no effect on the material response and,

                 (ii) the strain rate tensor and stress deviator are co-axial.

These postulates require that  (see, equation (4)) and


valid as the general constitutive relation for density preserving polycrystalline ice.

As his example, Steinemann chose a combined unilateral compression plus shear test as follows


which, under ideal experimental performance of a volume preserving isotropic material, yields


Steinemann did perform such creep tests and measured but not ; so, his experiments are insufficient to determine the functions  and  in (4) by experiment. However, one can test whether Nye’s postulates, leading to (6) are reasonable. His co-axiality postulate between  and  requires that

for which case, according to , we obtain


in which (11) follows by use of (6). Moreover, at , we have assumed isothermal conditions. If one now plots the values of  against  one should obtain a single curve, independent of whether compression is acting together with shear or not. Glen (1958) [59] performed this comparison, which is illustrated in (Figure 8); the figure caption reports on details. It is clearly seen that the experimental points cannot be regarded to approximate a single curve. It follows, co-axiality of the constitutive part of the stress,  and the strain-rate tensor D cannot be maintained as a postulate of a creep law in combined stress states of isotropic polycrystalline ice. We shall below have the occasion to see further arguments, which support this evidence. For multi-axial stress--strain-rate relations there is no experimental support for Nye’s co-axiality assumption between  and D. It is interesting that this more than 60 year old finding has not yet seriously been taken up in applied glaciology.

The combined compression plus shear test with the expressions (7) and (8) provides a possibility to study the relevance of the third invariant(s) as independent constitutive variables. Equations  and  imply

 • in simple shear with and , we have  and ,

 • in uni-axial compression/tension without shear,  and , we have  and

 It follows that independence of the constitutive relation for the Cauchy stress deviator from the third invariant implies that creep laws obtained for compression/extension and simple shear, respectively, must be brought to coincidence. If this coincidence cannot be verified by experiments to within acceptable errors, then a dependence of the stress strain-rate relation on the third invariant is likely.

To summarize

Multi-axial deformation tests (shear plus compression) show that isotropic ice does not satisfy Nye’s assumption that the creep law for the Cauchy stress is independent of the third invariant, see (Figure 8).

Questionable Consolidation of the Co-axiality Hypothesis in the 1960s to 1980s

·    The period after the formulation of the power law by Glen and Steinemann was not characterized by searching for, and extension of, the flow law to a general constitutive relation of a nonlinear fluid. Experimentalists rather seemed to be trapped in the simplicity and beauty of Nye’s two postulates, [135]. The power law was replaced by functions better matching the experiments (e.g. sinh (•) to the  power, see Barnes, et al. (1971), [10]). Another disadvantage of these ‘generalized flow laws’ was the problem of infinite viscosity (or vanishing fluidity) at zero stress. More precisely, the constitutive part of the stress would infinitely quickly change at stress initiation. This is unphysical and can be removed by replacing the power law for viscosity and fluidity, respectively, by e.g. polynomials, which include a constant term that mimics the finite viscosity at zero stress. These polynomials also allow a better fit of the coefficients of the creep law with experimental data, as demonstrated by Smith and Morland (1981), [155]. These authors work with effective strain-rates and stresses defined by


                               and write the constitutive law (6) for isothermal processes as  This can then be written as


where  accounts for a finite viscosity at zero stress, whereas  are reference scales and


The coefficients  for data of a number of experiments studied by Smith & Morland (1981) are determined by using an objective function minimizing the sum of the distances squared of all data points from the mathematical expressions (10) and (11), respectively.

(Figure 9) collects measured -data points from experiments performed by Glen (1953), [56], Mellor & Testa (1969a,b), [121,122], and Steinemann (1956), [158]. Panel (a) suggests a tremendous spread of the data; those obtained by different experimenters or from different machines hint at a separate coherent monotonic growth of the -functions, but obviously also point at data incoherence from one experimental site to another one. This interpretation is supported by the two -laws presented by Butkovich & Landauer (1960), [19] from their own experiments. They corroborate the accelerated monotonic growth but global inconsistence with the data of Glen, Steinemann and Mellor & Testa. The coherence of the 3-term polynomial representations for the data of separate experimentalists is demonstrated in (Figure 10).

To summarize:

·  Comparison of creep data of polycrystalline isotropic ice with co-axial strain-rate stress relations and independence of these relations of the third invariants, showed that the accelerated monotonic increase of  with  is a common feature of all these creep laws.

·  The power law description of the -relation is generally less accurate than polynomials, provided the number of terms in the polynomials, is small (J = 2 in our case).

·  The rate factor a(T) of the constitutive relation varies between different experimental sites.


The above findings suggest that in addition to the temperature dependence at least a further variable, of which the cause may be physical, instrumental or experimental, must also influence the strain-rate stress relation, viz.,


Many experiments are generally made with grains of pure ice mingled with pure water and then frozen, the rate factor E(?) in the last formula remains presently a mystery, because it does not even tell the variable. Incidentally, eq. (15), also states that the assumption of the thermo-rheological simplicity assumption of the constitutive relation for the stress tensor is likely invalid, see Morland & Lee (1960), [127]. We conclude:

  For simple shear and uniaxial compression Nye’s postulates give exactly the same class of constitutive functions, but if the ψ-functions determined by the experiments differ from one another; this difference is an indication that  must play a role.

  Steinemann performed this kind of combined experiment, measured  but not . However, this situation still allows at least to check whether the co-axiality postulate between  and  is satisfied. It follows from the definition of the second invariant of D and the co-axiality assumption of D,


On the other hand, equation (11) can for isothermal processes be written as


Glen constructed (Figure 8) more than 60 years ago and demonstrated that Nye’s postulate ought to be abandoned.

It is time to do this!

Recent climate relevant research that employs constitutive relations of polycrystalline isotropic ice as a viscous material generally assumes co-axiality of the strain-rate tensor and the Cauchy stress deviator that hardly deviates from the simple power law of the Glen-Steinemann type. This essentially remained so, even though facts to the contrary are well known to the experimental specialists. There is a necessity for a more general form of the constitutive relation with non-coaxiality of the strain rate tensor and the Cauchy stress deviator (.) The simplest such constitutive relations are the Reiner-Riwlin fluids, which for density preserving fluids have one of the forms


In these formulae, different from eq. (4), the stretching tensor  has been written as , where  is a temperature dependent rate factor. Moreover,

and and p is the pressure.


Morland and Staroszczyk (2019) [127] – whom we follow here - have performed a Gedanken experiment; it is a combination of compression and shear deformation, and it serves as a suggestion to catch the interest of experimentalists. So, let X and x be position vectors of ice particles in the reference and present configurations, respectively; then, compressing or extending and shear deformations can be described as


Here,  are normal stretchings and κ is the shear strain rate in the xz-direction. With these relations it is straightforward to derive for the velocity v


from which




Li & Jacka (1996), [121], Warner et. al. (1999), [172], Treverrow, et al. (2012), [162] and Budd, et al. (2013), [17], conducted experiments that fit the above conditions and performed these with the assumption


The requirement  can be tested for appropriateness, provided the lateral deformation is measured. This apparently has not been possible with the apparatus used by Budd, et al. (2013). [17] 

With the above expressions (19)-(21), restricted by relation (22), and the stress expressions (we write ) we obtain


and, upon using , we may deduce


in which  and  are functions of  and . With  given by , equations (24) & (25) allow evaluation of as follows: [32]


 are the applied stresses. In the Melbourne experiments (Budd, et al. (2013), [17]), the lateral extension is unconfined; this corresponds to


Moreover, eliminating between these two expressions  and then substituting the result into (26) with allows evaluation of  and then , (for details see Morland & Staroszczyk (2019), [130], viz.,

Finally, substituting these results into  leads to the equation


 This equation delivers  as a function of the measured quantities.

To summarize: The above results, due to Morland & Staroszczyk (2019), [130], are important for the following reasons:

1.     In the experiments the strain-rates  &  and the stresses  &  are measured or monitored. These are the quantities arising on the right-hand sides of (28); they can be calculated from the data obtained during the experimental performances. On the other hand,  are functions of  &  which, on the basis of (21) and (22), are given by


    which are equally point-wise known for these experiments.

2.  So, provided the experiment can be conducted as intended and described above, we have pointwise knowledge of


 The parameter j=1,2… J counts the number of experimental points. Kriging is popular method to  determine the mathematical expressions of  as functions of


3.    Equation (29) is independent of the response functions ; since ,  and  are measured or monitored, this relation determines  directly in terms of the measured data. On the other hand, once the –functions are determined by optimization from eq. (30),  can be evaluated from , using the optimized -functions. If this expression does not ‘reasonably’ reproduce measured -values, then the Reiner-Riwlin constitutive relation is too restrictive as a constitutive equation for glacier ice.

4.   Prerequisites of the results of the solution are that  requesting that the normal strain-rates  in the -direction are prevented, whereas those in the -direction balance those in the -direction (because of volume preserving). Since the normal stress component in the -direction is compressive, the corresponding - and -components of the strain-rates are extensional (but that in the -direction is prevented ( ), so  is a pressure.


It is evident that in principle combined compression-shear tests can be conducted, in which the scalar functions  of a Reiner-Riwlin fluid for ice can be determined. Morland & Staroszczyk (2019), [130], used data from Budd, et al. (2013), [13], and plot the data-related quantities  against  and obtain the plot in (Figure 11). If the ice specimens in these experiments were insensitive to variations of  these experimental points should lie on the -axis. Some experimental points indeed lie on this axis, but not all. Morland (2007), [113], shows that for independence of the -functions of  a unique relation  ) must exist, but co-axiality of the strain-rate and stress deviators, viz.,


cannot be maintained. So, (Figure 11) supports the Reiner-Riwlin structure of the constitutive relation         .

The result also confirms for the Budd, et al. (2013) experiments, [17], that the deviation from co-axiality is not small. In almost all cases ice is treated as volume preserving ( ) and the momentum equation is applied in the Stokes approximation, i.e., the acceleration terms are ignored.

Applications to Climate Dynamics 

Field equations for cold and temperate ice regions

Glaciers, ice sheets and ice shelves in climate relevant analyses are mostly treated as polycrystalline isotropic power law fluids, but restricted to constitutive relations of the form ‘ice plus inclusion of water’. This is so, even though careful experiments have shown that a complete Reiner-Riwlin structure with non-coaxiality of the stress-strain rate relation is likely better modeling the creep behavior of isotropic ice. The concentration of the water in the ice varies according to the thermal regime that prevails as a result of input of geothermal heat from the interior of the Earth, the heat flow exchange at the ice-atmosphere interface and the water production per unit mixture volume due to viscous heat of the mixture. Cold ice is ice with vanishing water content and temperature below the freezing point, while temperate ice is at the pressure melting point, in which the temperature is related to the pressure.

Furthermore, the contact surface of the cold and temperate ice – this surface is called the Cold Temperate Transition surface (CTS) at which the Clausius-Clapeyron equation must hold, is an internal singular surface. The thermomechanical conditions that are described by the jump conditions of the field equations in the cold and temperate subsets of the polythermal ice mass determine locally, how the melting and freezing processes evolve along the CTS. On the other hand, within the temperate ice the production of moisture equals the production of water mass per unit volume.

The governing field equations of ice in the cold ice region are the balances of mass, momentum (in the Stokes approximation) and internal energy. The constitutive relations describing the creeping flow are here given by the Reiner-Riwlin strain-rate--stress relation, eqs. (4) & (5) (i.e. not the usual Glen-Steinemann law), and the heat flux vector q is given by Fourier’s law of heat conduction, whereas the stress evolution is governed by the Stokes equations, and the temperature evolution follows from the internal energy balance. These equations are for a density preserving fluid of the form


and  are given by eqs. (5), while κ is generally assumed to be at most a function of T. In eqs. (32) v, p and T are the unknown fields. 

Correspondingly, in temperate ice regions, the physical laws are a class-I binary mixture; i.e. the equations comprise of the balances of mass, momentum and energy for the mixture as a whole plus a balance for the water mass. Interpreting the equations on the left-hand side of equations (32) as the mass, momentum and energy balances of the mixture as a whole, then these equations must be complemented by the balance of the water mass,


in which w is the moisture content per unit mass, j is the moisture flux and C the moisture production rate per unit volume of the mixture. For this mixture model constitutive relations are needed for j, C and the -relation, apart from the Fourier law for q. For j we postulate Fick’s law


where  is the diffusivity, and  may equally be a possible option. Moreover, considering that all internal energy production of the mixture is instantly used up by melting, energy and mass production can be interrelated with the aid of the latent heat of fusion L per unit volume of the mixture.  can, thus, be identified with L, the latent heat of fusion, while in the energy equation the heat flux contribution  is ignored, as  is small, since , where  is the Clausius-Clapeyron constant ( ) that is small. So, C =  and


in which eq. (4) has been substituted, here adequate for a Reiner-Riwlin fluid. Moreover, employing the Caley-Hamilton theorem for the stress deviator , it can be shown that


Substituting this into the above formula yields

It follows that


If Nye’s co-linearity assumption of D and  is imposed, then , and


as stated e.g. by Hutter (1983), [86]. Furthermore, according to Lliboutry (1979), [113], the presence of melt-water can significantly affect the constitutive relation of the stress deviator. He concludes that the constitutive equation for the creep law remains formally valid but coefficients depend now on the moisture content w and not on T, which is related to the pressure via the Clausius-Clapeyron equation. For the Reiner-Riwlin, fluid this means that


Analyses of the constitutive functions for the creep laws for temperate ice have also been given with more details than here by Blatter & Hutter (1991), [15], Fowler & Larson (1978), (1980), [49, 50],, Greve (1995, ’97, 2000), [67-70], Hutter (1982, ‘83), 2019), [83,86,90], Hutter, et al. (1988), [88], however only for the situation that the strain-rate is collinear to the stress deviator. For the Reiner-Riwlin fluid, the equation for the moisture content takes the form


Note that C depends also on the third stress deviator invariant, if stress and strain-rate are not co-axial. In any concrete situation, the union of these statements defines the initial boundary value problems (IBVP) that must mathematically be solved to arrive at a set of quantifications of the variables as functions of space and time that will provide information on the climate relevant questions for which the IBVPs were formulated.

Phase change properties of a viscous-heat conducting fluid are well known, see e.g. Hutter (1983 or 2019), [86,90], but need carefully be introduced for materials, which are kinematically handled as density preserving. The difficulty is that there cannot be a pressure melting formula and neither a Clausius-Clapeyron formula in a density preserving material. We, thus, will treat here compressible viscous fluids and present phase change properties of such materials and then will come back to the case of a density preserving material.

A phase change surface of a continuous material is defined as a special singular surface at which the temperature and the velocity component that is tangential to this surface are continuous. Thus,


where  is the speed tangential to , and  denotes a unit vector in the direction of . Apart from (40), also the ordinary jump conditions of the balance laws of mass, momentum and energy hold, as does a jump condition of entropy. These laws are


in which  is the velocity vector of a particle, and  is the velocity vector of the corresponding geometric point on the singular surface . These equations can be found in any book on continuum mechanics, e.g. Hutter (1983), [86], or Hutter & Wang (2016, 2018), [89]. Incorporating (40) into (41) yields


Hutter (2019), [90] gives proof of these relations. In this proof, the fact that the components of v and u that are tangential to  are continuous, have been employed.

In thermostatic equilibrium one has  It follows, the term on the right-hand side of  must vanish. This can be achieved by setting either one of the two factors equal to zero. We set , because the other alternative would lead to complicated expressions, not well interpretable; in words: the surface of phase change in thermodynamic equilibrium is material.

Next, eliminating  between  and then using  leads to the entropy jump condition

When both phases are in thermostatic equilibrium, then

The pressures on both sides of  equilibrate themselves, so eq. (40) yields

an equation that is usually written as


in which  is called equilibrium chemical potential.[33] Hence, in thermodynamic equilibrium, pressure and chemical potential are continuous across the phase change surface . Another familiar form of (43) is


 is called the equilibrium latent heat of fusion, , implying


It consists of the jump of the internal energy across  plus the power of working of the stresses due to the specific volume change across . A different way of writing it is


The interpretation of this equation is facilitated, if the relations


are applied. Both,  and  are functions of, i.e. the densities of the two phases and the temperature T which is continuous across. Therefore, at a given temperature, eqs. (47) may serve as equations to determine , respectively, as functions of T [34]. Once these densities are known, eq. (46) may be used to determine  as a function of the temperature alone. Hence, eq. (46) may be written as


Differentiating this expression with respect to the temperature yields (see Hutter 1983), [79],


The Clausius-Clapeyron equation is the inverse of this expression,


Since the temperature variations in problems involving phase changes are usually small, c in (50) is generally regarded as a constant.

Notice that in a density preserving material  and with