Hostname: page-component-848d4c4894-cjp7w Total loading time: 0 Render date: 2024-07-01T04:24:36.885Z Has data issue: false hasContentIssue false

The Dynamics of the Antarctic Ice Sheet

Published online by Cambridge University Press:  20 January 2017

W.F. Budd
Affiliation:
Meteorology Department, University of Melbourne, Parkville, Victoria 3052, Australia
D. Jenssen
Affiliation:
Meteorology Department, University of Melbourne, Parkville, Victoria 3052, Australia
Rights & Permissions [Opens in a new window]

Abstract

A three-dimensional dynamic, thermodynamic ice-sheet model has been developed to simulate the past, present, and future behaviour of the Antarctic ice sheet. The present ice velocities depend on the deep ice temperatures which in turn depend on the past changes of the ice sheet, including surface temperature, accumulation rate, and ice thickness. The basal temperatures are also strongly dependent on the geothermal heat flux. The model has therefore been used to study the effect on the basal temperatures, of changes to the geothermal heat flux, as well as the past changes of surface temperature and accumulation rate based on results obtained from the Vostok deep ice core. The model is also used to compute the distribution of surface velocity required to balance the present accumulation rate and the dynamics velocity based on the stress, temperature, and flow properties of ice, for the internal deformation, plus a component due to ice sliding. These velocities are compared to observed surface velocities in East Antarctica to assess the state of balance and the performance of the dynamics formulation.

Type
Research Article
Copyright
Copyright © International Glaciological Society 1989

Introduction

Prognostic dynamic models of the whole Antarctic ice sheet are needed for comprehensive global atmosphere-ice-ocean models to study long-term changes of climate and sea-level. A large part of the coastal ice flux from the Antarctic is channelled into narrow deep channels forming relatively fast-flowing glaciers or ice streams. Earlier dynamic modelling of the Antarctic ice sheet by Reference Budd and SmithBudd and Smith (1982) used a resolution of 100 km which was not adequate to simulate the differentiation into the high-speed flow of the outlet glaciers. It was therefore decided to increase the resolution to 20 km based on the digitization of the surface and bedrock from the maps of the SPRI Map folio by Reference DrewryDrewry (1983). Although these maps do not have accurate data over the whole of the continent, they provide a useful basis to start from, and the regions where the data are reliable are indicated by the coverage of flight lines and traverses on the maps. A first use of this digitized data set for computing the ice velocities which would be required to balance the present surface snow-accumulation regime (the “balance velocities”) was described by Reference Budd and SmithBudd and Smith (1985). Further uses of the high-resolution data for thermodynamic and dynamic modelling have been given by Reference Budd, Jenssen and SmithBudd and others (1984, 1987) and Reference Radok, Barry, Jenssen, Keen, Kiladis and MclnnesRadok and others (1986). This modelling computes the “dynamics velocities” from the stress and temperature distributions and the flow properties of ice. To allow more rapid computations of the full three-dimensional representation of the temperature and velocity distribution, the model has been reformulated by Jenssen and has been used to simulate the present ice sheet and its past and future changes.

A useful test for the performance of the model is the degree to which it can simulate the large number of surface velocities of the ice sheet now becoming available, e.g. from Reference Hamley, Smith and YoungHamley and others (1985) and Hamley (unpublished). The model computations of surface-dynamics velocity require deep ice temperatures which depend on the past history of the ice sheet. The model has therefore been run forward from 500 000 B.P. with a variety of past histories to study the impact on the basal temperatures. The results from the deep Vostok ice core have been used as a guide to the past changes in the interior of the ice sheet over the last 160 000 years for such variables as the surface temperature or accumulation rate, cf. Reference LoriusLorius and others (1985), Reference JouzelJouzel and others (1987), and Yiou and others (Reference Yiou, Raisbeck, Bourles, Lorius and Barkov1985).

The purpose of the present paper is to examine the computed surface velocities in relation to the velocities required for balance (with the present accumulation) and also to the observed velocities, in order to provide an indication of the present state of balance and the degree to which the observed velocities can be simulated by the dynamics. The sensitivity of the dynamics to the major unknowns such as the geothermal flux, the past accumulation rates, and the ice-flow relation is also examined.

Ice sliding also becomes important, particularly near the coast, and the observed flow of the outlet glaciers provides a check on the model ice-sliding formulation as given by Reference Mclnnes and BuddMclnnes and Budd (1984) and Reference Budd, Mclnnes, Jenssen, Smith, Veen van der and OeriemansBudd and others (1987) based on the empirical studies of Reference Budd, Jenssen, Radok, Budd, Keage and BlundyBudd and others (1979).

The details of earlier versions of the model have been described by Reference Radok, Barry, Jenssen, Keen, Kiladis and MclnnesRadok and others (1982, 1986). Some of the essential components are described in the Appendices. Only a brief outline is given here to indicate the basic principles of the model and the major changes from previous versions.

Outline of the Model

The model covers a domain of 281 = 281 grid points with 20 km spacing over the Antarctic region with 31 points in the vertical through the ice and the possibility for an additional 20 points in the bedrock reaching several kilometres depth. The equations for heat conduction and ice velocity are solved together, as given in Appendix I, allowing for internal heating varying with depth in the ice sheet. The horizontal shear stress (t xz) is taken as the only significant stress for ice deformation and heating except for basal sliding (b) in which case the effective normal stress is important and the frictional heating occurs at the basal interface. Vertical velocities are computed from the horizontal velocity profile (as given in Appendix I), and both vertical and horizontal advection are computed explicitly via finite differences from surrounding points. Space is not adequate here to present a complete coverage of the model but further details are given in the Appendices. The work is an advance from the earlier three-dimensional modelling of Reference Jenssen, Budd, Smith and RadokJenssen and others (1985). Reference Budd and JackaBudd and Jenssen (1987), and Reference Radok, Jenssen and MclnnesRadok and others (1987), in that flow lines are no longer explicitly required (but can be derived) since the model reaches a solution by iteration over the whole domain. A steady-state solution may be independently obtained by running the model for a long time (500 ka) with constant boundary conditions. Non-steady-state conditions can be derived from an initial state by computing forward in time with some prescribed varying boundary conditions.

Examples of the use of the model have been to run steady-state solutions for the present regime with different geothermal gradients, different past accumulation regimes, different ice-flow parameters, and different sea-level or bedrock conditions. Only a few examples are given here to examine the model's simulation of the present regime. By comparing both the model-computed balance velocity and dynamics velocity with observations, it should be possible to assess the degree to which steady state currently holds and the degree to which the computed velocities match reality. To compute dynamics velocities for the present regime, we first compute temperature distributions from prescribed computed balance velocities. Similarly, the temperature distribution has also been computed from lower accumulation rates in the past, as indicated from the Vostok ice-core results.

Assessment of Output

To begin with. Figure 1 shows the computed balance velocities based on the present regime (involving the surface elevations, accumulation rate and ice thickness). It should be noted that a change of a constant factor to the accumulation rate would simply result in the same constant-factor scale change in the balance velocities. In this regard, lower accumulation during the ice age would correspond to proportionally lower velocities. The Vostok ice-core data indicate that past accumulation rates have ranged between the present levels and about half those levels. The model results are therefore initially obtained for the present configuration with prescribed velocities first corresponding to the present values and then to half the present values. This allows us to examine the effects of the velocity distribution on the temperature. Similarly, the model can be computed with different prescribed accumulations. The model can also be run with a time series of the input variables changing as prescribed over the entire surface which is an improvement on the flow-line approach of Reference Budd and YoungBudd and Young (1983).

Fig. 1. Balance velocities. The distribution of balance velocities is shown as required for steady state with the present accumulation, ice thickness, and elevation contours.

In regard to the basal ice temperatures which have such a strong impact on the computed dynamics velocities, the prescribed geothermal heat flux (ϕ G) is of great importance. By examining the model results for a range of different values of this gradient, it appears that a reasonable central value to use is about ϕ G = 5.17 × 10−2 W m−2 [1.23 μcal cm−2 s−1]. This value converts to a basal temperature gradient (depending on the temperature) of roughly γG = 2.25°C/102 m (see Table I). A geothermal gradient of this amount has also been found to provide a close fit to the observed deep temperature profile at Vostok (Jenssen and Budd, in press). Figure 2 shows the pattern of

Table I Geothermal heat flux and temperature gradients in ice

Ice conductivity, at 0°C = 2.24 W m−1 °C−1.

Geothermal heat flux unit 1 μcal cm−2 s−1 = 4.187 × 10−2 W m−2.

computed basal temperatures derived from the prescribed geothermal fluxes ϕ G = 4.62, 5.17, 5.75 × 10−2W m−2, based on the steady state under the present regime. Similar calculations have been carried out with the accumulation rate half its present value or the velocities half the balance velocities or both of these changes together, as shown in Figure 2(d).

The dynamics velocity d is computed from the internal deformation (V i) plus the sliding velocity (V s)

The internal velocity is computed from the integration of the shear strain-rate (

) through the ice related to the shear stress (τxy) and temperature (θ). The flow relation used is of the form

where the temperature-dependent factor (θ) is given by the table determined for tertiary flow as reported by Budd and Jacka (1987, in press). The numerical values are similar to those given by Reference Paterson and BuddPaterson and Budd (1982) but increased by about a factor of 3 for tertiary flow rates.

The sliding velocity (s) is taken as the general form derived empirically by Reference Budd and YoungBudd and others (1979), with an additional factor for the effect of temperature:

where is the effective normal stress (above buoyancy) relative to pressure-melting point, and are simple monotonie functions, θ is the temperature, and ν ≈ 0.1°C−1. The particular form used here is given in Appendix II, Equation (A2-4).

The total dynamics velocity is shown in Figure 3 and the sliding component in Figure 4. The sliding velocity is generally small for temperatures below pressure melting. The sliding velocity also increases strongly towards the coast as the normal stress decreases and thus allows high-speed flow in the ice streams even though the shear stress may decrease.

A comparison of the total dynamics velocity with the balance velocity shows general broad similarity but with some notable regions of difference. The largest such region is east of the Filchner Ice Shelf where the dynamics velocities are much below the balance velocities. The discrepancy is more than can be reasonably expected from the uncertainty in the ice-flow parameters or the temperature calculations. It is suggested that either the input data (from the SPRI maps where coverage is sparse) are grossly in error or the ice sheet there is grossly out of balance and may be still building up in response to the post-ice-age accumulation increases.

Fig. 2. Basal temperature. The computed basal temperatures are shown derived for steady state based on the present regime for the following geothermal flux in mW m−2: (a) ϕ G = 46.2, (b) ϕ G = 51.7, (c) ϕ G = 57.5. In (d) the corresponding results are for ϕ G = 51.7 but for the accumulation rate and balance velocity distributions at one-half of their present-day values.

Fig. 3. Dynamics velocity. The velocity distribution computed from the temperature distribution obtained using ϕ G = 50.4 mW m−2 and the flow and sliding relations given in Appendix I. There are broad similarities with the balance velocities but notable regions of difference such as up-stream of the Shackleton Range in Dronning Maud Land.

Fig. 4. Sliding velocity. The distribution of sliding velocity computed from the shear and effective normal stresses shows low values inland but high values near the coast particularly in the region of the deep ice streams and outlet glaciers.

Fig. 5. Comparison with observations in East Antarctica. For the region of the IAGP surveys, the observed velocities from Hamley and others are shown in (b). The computed balance velocities are shown in (c) and the dynamics velocities in (d). The model elevations are shown in (a).

Fig. 6. Balance and dynamics velocities versus observed. The balance velocities are plotted against the observed velocities from Figure 5 in (a). The corresponding dynamics velocities are plotted against the observed in (b). Although considerable scatter exists, both show reasonable trends which can be used to assess balance and the effective flow parameters.

Finally, in Figure 5 a comparison is made with the East Antarctic region covered by the IAGP (International Antarctic Glaciology Project). Here, the observed velocities from Reference Budd and YoungBudd and Young (1979), Reference Hamley, Smith and YoungHamley and others (1985), and Hamley (unpublished) in Figure 5a are compared with the balance velocities (Figs 5b and 6a) and the computed dynamics velocities in Figures 5c and 6b. Although there is considerable variation between the results, the general trend of increase of velocity towards the coast, particularly over the deep bedrock regions, is common to all three.

This is not a trivial result because it involves stresses which on some flow lines increase from 0.5 bar to over 1.5 bar then decrease to 0.5 bar while the velocity continues to increase (due to sliding) near the coast.

Since there is inherent error or uncertainty in each of the velocity distributions, statistical regression is useful to compare the observations with both the balance and dynamics velocities. The first conclusion is similar to that expressed by Reference Hamley, Smith and YoungHamley and others (1985) that the observations in this region do not differ substantially from balance. This involves also some uncertainty in the comparison for the ratio of the average column velocity to surface velocity which is expected (from the modelling) to be generally in the range of 80-90% until the high-sliding zones are reached where it may approach 100%. The second conclusion is that the present dynamics velocity formulation gives a reasonable representation of the observed velocities. The regions where large differences appear between balance and dynamics velocities may well be areas in which the ice sheet is not in balance. On the other hand, in those regions, there may be errors in the data used as input to the model which affect the dynamics velocities. Observed velocities in these regions, as well as improved observations for the input data, would be particularly valuable for clarifying both the state of balance and the dynamics.

In the interim, the comparison between observed and computed dynamics velocities provides some encouragement to use the computed dynamics velocities interactively for the time-varying past conditions through the ice-age changes.

Appendix I

Deformation Velocity and Flow Relation

For a coordinate system with axes horizontal in the line of flow, vertical, and across the flow, it is assumed for the large horizontal scales of the grid that the flow is governed primarily by the horizontal shear stress τxz, given by

(A1-1)

where ρ is the ice density, is the gravitational

Table II Ice-flow rates

at minimum for isotropic ice and stress τ0 = 1 bar (octahedral)

Note: The temperature used for flow in the model is the temperature difference from pressure-melting point at each depth.

acceleration, α is the surface slope, and is the depth below the surface.

The flow relation used is of the form which relates the shear strain-rate at depth z,

, to the shear stress by
(A1-2)

where = 3 and θ is dependent on the temperature and ice properties.

For this context, the dependence on temperature (θ) is given by the data of Budd and Jacka (in press) as shown in Table II, shown for octahedral shear stress τ0 = 1 bar, and strain-rate

, from laboratory measurements.

The relations between the octahedral and horizontal shear stress and strain-rate are

(A1-3)

Table II gives the strain-rate at minimum for isotropic ice tests. A further factor 1 is used with the flow law to allow for anisotropy, tertiary flow, other factors which may affect the flow of ice, or to match with observed surface velocities.

The velocity depth gradient d/d is then given by

(A1-4)

The horizontal velocity at depth is given by

(A1-5)

where s is the surface velocity minus the basal sliding and θ is determined at each depth using Table II from the temperature profile from the heat conduction equation as given in Appendix II.

The Vertical Velocity Profile

Following Reference Budd, Jenssen, Radok, Budd, Keage and BlundyBudd and others (1970), the usual vertical coordinate, , is replaced by the relative coordinate ξ = (Ez)/Z, where Z is the ice thickness and is the ice-surface elevation. The continuity equation for an ice element is then

(a)

where

is the horizontal velocity, and (= dξ/dt) is the relative vertical velocity.

The continuity equation for the column as a whole is

(b)

where

is the average horizontal velocity in the column, and is the surface-accumulation rate. It is assumed that
(c)

where ψ is a function of ξ (or ) which may be computed from the flow law used. Use of Equations (a) and (c), with substitution into Equation (b) yields, after some manipulation,

Replacing ξ by (1 - ζ) and integrating from bedrock (ζ = 0) to any general ζ gives, since ζ = 0 at the ice base,

(d)

where

(e)

Appendix II

Heat-Conduction Equation Including Internal Deformation Heating and Advection

The basic thermodynamic equation used in the numerical model is

(A2-1)

where θ is the temperature (°C), is the ice thickness (m),

is the horizontal velocity (m a−1), is the heat generation from deformation as given below, κ is the thermal diffusivity, κ′ is essentially the derivative of κ with respect to temperature (m2 a−1°C−1), ξ is a relative vertical coordinate with 0 at the surface and 1 at the base, and is the relative vertical velocity at ξ (w = dξ/dt).

The quantities Z and A are not functions of depth but other quantities are. In addition, κ, κ′, q, and are functions of temperature so that Equation (A2-l) is highly non-linear.

The horizontal velocity u is computed from the flow law as given in Appendix I, involving the temperature and stress at any point. The average velocity through the column is u .

The heating q is considered in two parts: q i which is internal deformational heating, and q b which is basal frictional heating from sliding.

(A2-2)

where is the thermal conductivity (kg m a−3 °C−1) and

are computed from the ice-flow law as given in Appendix I.

For the sliding component,

(A2-3)

where b is the basal sliding velocity. Note is a function of the base temperature.

The computation of basal sliding is based on the empirical sliding studies of Reference Budd and YoungBudd and others (1979) and the comparison of dynamics and balance velocities in West Antarctica by Reference Budd, Jenssen and SmithBudd and others (1984).

The current form used is given by

(A2-4)

where * is the ice thickness above buoyancy (i.e. that required for floating) and 2 3, and 4 are constant, 2 = 1.5 × 108 m3 year1 bar1, 3 = 50 m2, 4 = 0.0035 m−1, θ2 is the temperature relative to pressure-melting point at ν = 0.1 °C−1.

For cases in which the model is run with prescribed (surface or balance) velocities, the total velocity is normalized to the prescribed velocities but the fractions of sliding and internal deformation are determined by the above relations.

References

Reference016

Budd, W.F. Jacka, T.H.. 1987 Ice rheology relevant to ice sheet modelling. International Association of Hydrological Sciences Publication 170 (Symposium at Vancouver 1987 — The Physical Basis of Ice Sheet Modelling), 29. Google Scholar
Budd, W.F. Jacka, T.H.. In press. A review of the rheology of ice for ice sheet modelling. Cold Reg Sci. Technol. Google Scholar
Budd, W.F. Jenssen, D.. 1987 Numerical modelling of large–scale basal water flux under the West Antarctic ice sheet. In Veen van der, C.J.,Oeriemans, J. , eds. Dynamics of the West Antarctic Ice Sheet. Proceedings of a Workshop held in Utrecht, May 6–8, 1985 Dordrecht, etc., D. Reidel Publishing Company, 293320. Google Scholar
Budd, W.F. Smith, I.N.. 1982 Large–scale numerical modelling of the Antarctic ice sheet. Ann. Glacial., 3, 4249. CrossRefGoogle Scholar
Budd, W.F. Smith, I.N.. 1985 The state of balance of the Antarctic ice sheet — an updated assessment. In Glaciers, Ice Sheets and Sea Level: Effects of a CO2–induced Climatic Change. Report of a Workshop held in Seattle, Washington, September 13–15, 1984 Washington,DC .US Department of Energy, 172177. Google Scholar
Budd, W.F. Young, N.W.. 1979 Results from the I.A.G.P. flow–line study inland of Casey, Wilkes Land, Antarctica. J. Glacial 24(90), 89101. CrossRefGoogle Scholar
Budd, W.F. Young, N.W..1983a. Application of modelling techniques to measured profiles of temperatures and isotopes. In Robin G. climatic record in polar ice sheets. Cambridge, etc., Cambridge University Press, 150177. Google Scholar
Budd, W.F. Young, N.W. 1983b. Techniques for the analysis of temperature–depth profiles in ice sheets. In Robin, G. climatic record in polar ice sheets. Cambridge, etc., Cambridge University Press, 145150. Google Scholar
Budd, W.F. Jenssen, D. Radok, U.. 1970 Derived physical characteristics of the Antarctic ice sheet.ANARE Interim Rep. 120. Budd, W.F. Keage, P.L. Blundy, N.A.. 1979 Empirical studies of ice sliding. J. Glaciol 23(89), 157170. Google Scholar
Budd, W.F. Jenssen, D. Smith, I.N.. 1984 A three–dimensional time–dependent model of the West Antarctic ice sheet. Ann. Glaciol., 5, 2936. Google Scholar
Budd, W.F. Mclnnes, D Jenssen, B.J Smith, I.N.. 1987 Modelling the response of the West Antarctic ice sheet to a climatic warming. In Veen van der, C.J.,Oeriemans, J., eds. Dynamics of the West Antarctic Ice Sheet. Proceedings of a Workshop held in Utrecht, May 6–8, 1985 Dordrecht, etc., D. Reidel Publishing Company, 321358. Google Scholar
Drewry, D.J. ed. 1983.Antarctica: glaciological and geophysical folio. Cambridge, University of Cambridge Scott Polar Research Institute. Google Scholar
Hamley, T.C. Unpublished. The East Antarctic ice from ice sheet flow to iceberg dissolution. (M.S thesis , University of Melbourne, 1987) Google Scholar
Hamley, T.C. Smith, I.N. Young, N.W.. 1985 Mass–balance and ice flow–law parameters for East Antarctica. J. Glaciol., 31(109), 334339. Google Scholar
Jenssen, D. Budd, W.F. In press. A–examination of the Vostok temperature–depth profile. ANARE Res Notes Google Scholar
Jenssen, D. Budd, W.F. Smith, I.N. Radok, U.. 1985.On the surging potential of polar ice streams part II Ice streams and physical characteristics of the Ross Sea drainage basin west Antarctica. Melbourne, University of Melbourne Meteorology Department; Boulder, CO, University of Colorado Cooperative Institute for Research in Environmental Sciences. Google Scholar
Jouzel, J. and 6 others. 1987 Vostok ice core: a continuous isotope temperature record over the last climatic cycle (160,000 years). Nature, 329(6138), 403408. Google Scholar
Lorius, C. and 6 others. 1985 A 15,000–year climatic record from Antarctic ice. Nature. 315(6029), 591596. Google Scholar
Mclnnes, B.J. Budd, W.F.. 1984 A cross–sectional model for West Antarctica. Ann. Glaciol., 5, 9599. Google Scholar
Paterson, W.S.B. Budd, W.F.. 1982 Flow parameters for ice sheet modelling. Cold Reg. Sci. Technol., 6(2), 175177. Google Scholar
Radok, U. Barry, R.G Jenssen, D. Keen, R.A. Kiladis, G.N Mclnnes, B.. 1982.Climatic and physical characteristics of the Greenland ice sheet. Boulder, CO, University of Colorado Cooperative Institute for Research in Environmental Sciences. (Report CO 80309.) Google Scholar
Brown, D. Jenssen-Smith, I.N. Budd, W.F.. 1986.On the surging potential of polar ice streams Part I Antarctic ice accumulation basins and their main discharge regions. Washington,DC .US Department of Energy; Melbourne, University of Melbourne Meteorology Department; Boulder, CO, University of Colorado Cooperative Institute for Research in Environmental Sciences. Google Scholar
Radok, U. Jenssen, D. Mclnnes, B.J.. 1987.On the surging potential of polar ice streams. Antartic surges — a clear and present danger? Washington,DC .US Department of Energy; Melbourne, University of Melbourne Meteorology Department; Boulder, CO, University of Colorado Cooperative Institute for Research in Environmental Studies. Google Scholar
Yiou, F. Raisbeck, G.M. Bourles, D. Lorius, C. Barkov, N.I.. 1985 10Be in ice at Vostok Antarctica during the last climatic cycle. Nature, 316(6029), 616617. Google Scholar
Figure 0

Fig. 1. Balance velocities. The distribution of balance velocities is shown as required for steady state with the present accumulation, ice thickness, and elevation contours.

Figure 1

Table I Geothermal heat flux and temperature gradients in ice

Figure 2

Fig. 2. Basal temperature. The computed basal temperatures are shown derived for steady state based on the present regime for the following geothermal flux in mW m−2: (a) ϕG = 46.2, (b) ϕG = 51.7, (c) ϕG = 57.5. In (d) the corresponding results are for ϕG = 51.7 but for the accumulation rate and balance velocity distributions at one-half of their present-day values.

Figure 3

Fig. 3. Dynamics velocity. The velocity distribution computed from the temperature distribution obtained using ϕG = 50.4 mW m−2 and the flow and sliding relations given in Appendix I. There are broad similarities with the balance velocities but notable regions of difference such as up-stream of the Shackleton Range in Dronning Maud Land.

Figure 4

Fig. 4. Sliding velocity. The distribution of sliding velocity computed from the shear and effective normal stresses shows low values inland but high values near the coast particularly in the region of the deep ice streams and outlet glaciers.

Figure 5

Fig. 5. Comparison with observations in East Antarctica. For the region of the IAGP surveys, the observed velocities from Hamley and others are shown in (b). The computed balance velocities are shown in (c) and the dynamics velocities in (d). The model elevations are shown in (a).

Figure 6

Fig. 6. Balance and dynamics velocities versus observed. The balance velocities are plotted against the observed velocities from Figure 5 in (a). The corresponding dynamics velocities are plotted against the observed in (b). Although considerable scatter exists, both show reasonable trends which can be used to assess balance and the effective flow parameters.

Figure 7

Table II Ice-flow rates at minimum for isotropic ice and stress τ0 = 1 bar (octahedral)