Introduction
Reference Meier, Meier, Rasmussen, Krimmel, Olsen and FrankMeier and others (1985) have recently published a remarkably comprehensive photogrammetric study of Columbia Glacier, Alaska, which includes (as plate 7) a detailed map of its lower part showing the trajectories of principal strain-rate on the surface, averaged over the measurement year. The pattern of trajectories displays two isotropic points and it is interesting to compare them with theoretical predictions made for glaciers in general (Reference NyeNye, 1983). The trajectories of the lesser (more compressive) principal strain-rate are shown schematically in Figure 1, with the isotropic points denoted P1 and P2.
The Upper Isotropic Point P1
This point is associated with a change from longitudinal extension to longitudinal compression in the lower part of the glacier. From the observations, Reference Meier, Meier, Rasmussen, Krimmel, Olsen and FrankMeier and others (1985) concluded that the pattern at this point is lemon (L) (Fig. 2), in contrast to the theoretical prediction that it would be monstar (M), The essential difference between the lemon and monstar patterns is that the latter has three straight-line trajectories running into the isotropic point, while the former has only one. Observationally, it can be very difficult to distinguish between the two patterns and I believe that, given the density of data points used in the survey of Columbia Glacier, either interpretation is possible.
To explain this in more detail let me take a coordinate system, as in my paper, with origin at P1, with x longitudinal and increasing down-glacier, and with y trans-verse, and put
where έ xx and έ xy are tensor strain-rate components, both derivatives being taken at the origin. The strain-rate tensor is then approximated in the neighbourhood of P1 by
where E is the common value of έ xx and έ yy at P1. Note that this approximation includes only the two most important linear variations of the components. The monstar prediction was based on the assumption that, in most glaciers, 0 < a < b. At P1, on Columbia Glacier there is approximate equality between the two dominant derivatives, 0 < a ≈ b, and consequently the point is borderline monstar–lemon.
If the three straight lines of the monstar pattern are opened out as far as the data allow, the two outer lines make angles of ±12° with the x-axis. From the formula tan-1 , this corresponds to a/b = 0.95. Thus, if the data are interpreted most favourably for the monstar pattern, the condition a < b is only just satisfied. As a/b passes through 1, the three straight lines close up and the pattern becomes lemon. As Reference Meier, Meier, Rasmussen, Krimmel, Olsen and FrankMeier and others (1985) indicated, the data can also be fitted with a lemon pattern, corresponding to a > b. My conclusion is that one cannot be certain about the pattern category. In a long, narrow glacier one would expect transverse variations, represented by b, to dominate longitudinal variations, represented by a, and hence a clear monstar pattern. The lower part of Columbia Glacier is evidently too wide for this.
The Lower Isotropic Point P2
Although much of Columbia Glacier has its bed below sea-level, it is grounded throughout its length, ending in an ice cliff standing in the sea and producing icebergs. The lower isotropic point P2 is associated with the change from longitudinal compression to extension in the terminal part of the glacier. Thus, taking the origin now at P2, we have a < 0 < b. The index discriminant given by Reference Thorndike, Thorndike, Cooley and NyeThorndike and others (1978) reduces to D I = ab, and is therefore negative, giving index -½. The line discriminant is D L = 12b(b–a)3 and is therefore positive, indicating three lines. The theoretical pattern is therefore star, in agreement with the observations. We see that the star pattern is inevitable if the arrangement of signs is a < 0 < b. Drag from the valley sides ensures that 0 < b in the centre of the glacier (although not necessarily near the edges), while the absence of restraint at the terminus must surely entail a < 0. Thus the star pattern should be universal for these conditions.
The Contour Classification
The third kind of classification for an isotropic point, in addition to index (±½) and line (1 or 3), relates to the shape of the contours showing the magnitudes of the two principal strain-rates near the point. These form two inter-secting families of curves. Generically, both families are ellipses or both families are hyperbolas. For the tensor in Equation (1) the contour discriminant reduces to D c = 0, showing that we have the borderline, non-generic, case where the contours are parabolas. This is structurally unstable, and the perturbations represented by linear terms omitted from Equation (1), which must be present, will change the contours into either very elongated ellipses or hyperbolas. The observations on Columbia Glacier cannot be expected to make such fine distinctions, either for P1 or P2.
However, as shown in Reference NyeNye (1983), the contour classifi-cation for stress, rather than strain-rate, is structurally stable. The contour discriminant has the same sign as –a 2 b 2, and is therefore unconditionally negative (unless a = b = 0). This indicates hyperbolic, rather than elliptic, contours for the magnitudes of the principal stresses around both P1 and P2.
We may note that P1 and P2 are special, as isotropic points, in that E ≈ 0; that is, έ xx = έ yy ≈ 0 at the points themselves. This means that they occur nearly, but not exactly, at the two places (on the line defined by έ xy = 0) where the longitudinal strain-rate έ xx changes sign. Thus, P1 is very close to the point where longitudinal extension changes to compression, and P2 is very close to the point where it changes back to extension.
Stationary Points for the Speed
Reference Meier, Meier, Rasmussen, Krimmel, Olsen and FrankMeier and others (1985) also showed, as plate 6, a contour map of the speed of the lower part of Columbia Glacier, averaged over the measurement year. It is notable that the two isotropic points P1 and P2 that they show in plate 7 are close to, but not coincident with, a maximum and a saddle point, respectively, in the contours of speed. If we take a local coordinate system with origin at the isotropic point and if we assume v = 0, then the condition for a stationary point in the speed |u |, namely
is identical with that for an isotropic point, which is
Since, if v were uniformly zero, the stationary points for speed would coincide precisely with the isotropic points, any small lack of coincidence must be attributed to the presence of a small v component. The stationary points for speed and the isotropic points are individually structurally stable, but the coincidences between them are only approximate because they are unstable to perturbations of the transverse velocity.
Plate 5 of Reference Meier, Meier, Rasmussen, Krimmel, Olsen and FrankMeier and others (1985) shows contours for the speed over the entire Columbia Glacier system, which contains several tributaries. There are about 34 points where the speed is stationary, some 14 of which are on the main ice stream, and they are all maxima or saddles, with no minima. This is explained by the fact that is almost always negative, while can be positive (saddle) or negative (maximum). If we neglect v, and The analysis of P1 and P2 then suggests that these 34 stationary points for speed are close to 34 isotropic points (on the measurement scale adopted), the maxima corresponding to lemons or monstars and the saddles to stars. On any one ice stream, maxima for speed must alternate with saddles (Fig. 3) and, likewise, lemons or monstars must alternate with stars. All the isotropic points will have hyperbolic, rather than elliptic, contours of magnitude of principal stress.