Discussion

J. F. Nye: When computing non-steady distributions, what do you do about specifying initial conditions?

W. F. Budd: For the non-steady-state situation one can start from steady state and then apply changes to the input. Alternatively, the ice masses can be grown from zero thickness.

Nye: HOW do you deal with the necessary boundary condition for stress, etc., at the up-stream end of the field?

Budd: We have not worried about the tension at a divide which in non-steady state is unknown. However, our grid points are widely spaced to avoid the effects of longitudinal stresses and so the first point away from the divide is no longer affected by the unknown boundary stress.

L. Lliboutry: I wonder whether you have not forgotten some equation, since the rate of change of the surface altitude ought not to have to be put into the calculation. It should be found as an output if all the appropriate boundary conditions and equations are used.

Budd: If all the data described for the input were available, then the model would be overprescribed as you suggest. However, we do not know the data well and have to guess or try various possible outcomes. For example, if the flow law were prescribed the velocity could be calculated. Alternatively, if a velocity profile is prescribed a flow-law parameter could be calculated. Similarly, if the accumulation is prescribed the balance (or rate of change) could be calculated or vice versa.

Lliboutry: Duval’s experiments have shown that the normal ice fabric in moving ice (4 or 5 maxima forming a girdle around the compressive axis), although it allows only a quarter of the strain-rate compared with the simple-shear situation, is a very stable one. Nevertheless, in the time scale valid for ice sheets there may be a slow change in the fabric from the first to the second one. In this case the ice flow law should change with time. (A more perfect icc lattice would also lead to this result.)

Budd: From laboratory studies it seems that the girdle with 4 or 5 maxima can be caused by random sample clustering (Reference BuddBudd, 1972). We now have a reasonable picture of the crystal fabric distribution on the Law Dome flow line which seems to be consistent with the stress situation. If steady-state applied then the spatial distribution of fabrics through the icc sheet would be appropriate for calculating the flow. For changing icc sheets I agree that the development and change of the crystal fabrics should be considered.

G. DE Q. Robin: Your paper “Derived physical characteristics” may be described as determining physical characteristics (temperature, deformation, etc.) of internal layers from given boundary conditions. Your excellent curve-fitting techniques may be described, however, as determining boundary conditions from internal physical characteristics. Would you agree that the top 800 m of “Byrd” core and temperature profile are little affected, if at all, by basal and near-basal temperature gradients, and hence provide the best test to date of the agreement of isotopically derived temperatures with a measured temperature depth profile?

Budd: Yes, for the “Byrd” profile the indications are that the temperature variations assumed from the isotope profile give rise to the measured temperature profile provided the other parameters assumed, such as the accumulation rates and horizontal velocity, have remained at their current values. However, it needs to be noted that the temperature change could be either due to elevation, or climate change, or some as yet unknown combination of both. The effect of an ice-sheet surge in the region should also not be overlooked.

M. W. Mahaffy: In the graph showing the relation between surface and bed undulations along a cross-section, what was the relative size of the noise amplitude to the undulations?

Budd: The errors in the surface elevation were about ± 1 m relative and in the ice thickness about ±1 %, i.e. about 20 m. In both cases these errors were more than an order of magnitude below the amplitude of the undulations.