Blue Glacier was subject to 5 to 10% increase in thickness during the period 1958 to 1978. A surface velocity increase of 15 to 50% has accompanied this change in ice thickness. A non-uniform distribution of thickening over the glacier produced a general decrease in surface slope. The flow law parameters in a power-law relation for ice can be obtained from the resulting increase in velocity.

In 1957–58 the surface velocities and surface elevation along several transverse profiles were determined by others. During the period 1977 to 1980 similar measurements were made of surface velocity, elevation, and ice thickness along profiles spatially equivalent to those of the 1957–58 survey. Using these results a detailed study of the flow response has been made. Preliminary results indicate that the controlling surface slope a is one that is averaged over approximately ten times the glacier thickness and also indicate a relatively high stress exponent n = 4 to 5.

Blue Glacier flows in a complex channel and any detailed interpretation of the flow response requires the inclusion of this geometry in a model. For this reason multi-dimensional finite-element methods using a non-linear constitutive relation were employed. Various degrees of model sophistication allowed for the specification of different parameters governing the flow of the glacier. Variations in the flow law parameters and the curvature of the channel were found to have significant effects on the velocity and stress fields, as well as on the characteristics of the flow response to the thickness perturbation. Fully three-dimensional models of the flow of Blue Glacier and its response to the thickness increase were developed. Important effects such as those arising from longitudinal thickness and surface slope variations were included in three-dimensional modeling. The three-dimensional effects were seen to play an important role in flow of a rheologically non-linear fluid such as ice.