Hostname: page-component-586b7cd67f-rdxmf Total loading time: 0 Render date: 2024-11-22T13:08:11.125Z Has data issue: false hasContentIssue false

Geodetic Results of the Ross Ice Shelf Survey Expeditions, 1962–63 and 1965–66

Published online by Cambridge University Press:  30 January 2017

Egon Dorrer
Affiliation:
Lehrstuhl für Photogrammetrie und Kartographie, Technische Universität Braunschweig, 33 Braunschweig, West Germany
Walther Hofmann
Affiliation:
Lehrstuhl für Photogrammetrie und Kartographie, Technische Universität Braunschweig, 33 Braunschweig, West Germany
Wilfried Seufert
Affiliation:
Lehrstuhl für Photogrammetrie und Kartographie, Technische Universität Braunschweig, 33 Braunschweig, West Germany
Rights & Permissions [Opens in a new window]

Abstract

By means of modern geodetic observation techniques the ice movement along an east-west and a north-south profile across the Ross Ice Shelf, Antarctica, was measured during the two Antarctic summers, 1962–63 and 1965–66. 103 markers were placed on the 910 km long traverse. Distances were measured by tellurometer, and traverse angles by a precision theodolite between all consecutive markers, normally 8 to 9 km apart. For this type of observation method, six men distributed into three groups of two men each were necessary.

The main part of the paper deals with data processing and with the computation of the ice movement. As the ice moves, the geometrical configuration of the traverse changes during the epoch of observation. For this “reduction to epoch” problem two methods are described in detail: (1) time reduction of observations, and (2) time reduction of positions. Between the two field journeys, only linear ice movement can be assumed. It is possible, however, to determine acceleration and curvature of the ice flow at all traverse points where the traverse angles differ considerably from 180°.

The result of all computations is the field of velocity vectors along the traverse. Obvious characteristics are the rapid increase of velocity between the McMurdo Ice Shelf and Ross Ice Shelf, the uniform and nearly parallel movement in the middle of the ice shelf (maximum velocity 935 m year−1), the decrease of velocity along the north-south profile, and the systematic increase of divergence of the flow lines towards the ice margins. Careful study of the velocity vector field shows some deviations from an entirely uniform distribution.

Résumé

Résumé

Au moyen de techniques d’observations géodésiques modernes, le mouvement de la glace le long de profils est–ouest et nord–sud sur le Ross Ice Shelf, Antarctique, a été mesuré durant les deux étés antarctiques 1962–63 et 1965–66. Au total, 103 balises ont été placées sur une distance de 910 km. Les distances étaient mesurées à l’aide de telluromètres, et les angles à l’aide d’un théodolite de précision à chaque balise, normalement distantes de 8 à 9 km. Pour ce type de méthode d’observation, 6 personnes réparties en trois groupes de deux étaient nécessaires.

Le chapitre principal porte sur l’exploitation des données et la détermination du mouvement de la glace. Comme la glace s’écoule, la configuration géométrique de la polygonale change pendant la période d’observation. Pour ce problème de réduction dans le temps, deux méthodes sont décrites en détail: (1) réduction du temps des observations et (2) réduction du temps des positions. Entre les deux campagnes, un mouvement linéaire de la glace peut seul être admis. Il est cependant possible de déterminer l’accélération et la courbure de l’écoulement de la glace en chaque point de la polygonale où les angles diffèrent beaucoup de 180°.

Le résultat obtenu est le champ des vecteurs vitesse le long de la polygonale. Les caractéristiques marquantes sont l’augmentation de vitesse entre McMurdo et le Ross Ice Shelf, le mouvement uniforme et quasi parallèle au milieu de l’ice shelf (maximum 935 m an−1), la diminution de la vitesse le long du profil nord–sud, et la divergence systématique des lignes d’écoulement vers les bords de la glace. L’étude soignée du champ des vecteurs de vitesse montre quelques déviations d’une distribution entièrement uniforme.

Zusammenfassung

Zusammenfassung

Mit Hilfe moderner geodätischer Beobachtungsmethoden wurde während der zwei Südsommer 1962–63 und 1965–66 die Eisbewegung längs eines Ostwest- und eines Nordsüd-Profils auf dem Ross Ice Shelf, Antarktis, bestimmt. Insgesamt wurden auf der 910 km langen Strecke 103 Geschwindigkeitspegel errichtet. Die Abstände aufeinanderfolgender, normalerweise 8 bis 9 km voneinander entfernter Pegel wurden mittels Tellurometer, die Polygonwinkel mit einem Präzisionstheodolit gemessen. Für diese Beobachtungsmethode waren sechs Mann, verteilt auf drei Gruppen zu je zwei Mann, nötig.

Das Hauptkapitel behandelt die Verarbeitung der Messdaten und die Berechnung der Eisbewegung. Da sich das Eis ständig bewegt, ändert sich die geometrische Konfiguration des Polygonzuges während des Beobachtungszeitraumes. Zwei Methoden zur Reduktion auf einen bestimmten Zeitpunkt werden ausführlich beschrieben: (1) Zeitreduktion von Beobachtungen. (2) Zeitreduktion von Punktlagen. Zwischen den beiden Messkampagnen kann nur eine lineare Eisbewegung angenommen werden. Es ist jedoch möglich, Beschleunigung und Krümmung der Strömungslinien an all den Polygonpunkten zu bestimmen, deren Polygonwinkel sich beträchtlich von 180° unterscheiden.

Das Resultat der Berechnungen ist das Geschwindigkeitsvektorfeld längs des Polygonzugs. Auffällig ist der rapide Geschwindigkeitsanstieg zwischen McMurdo- und Ross Ice Shelf, die gleichförmige und nahezu parallele Bewegung in der Mitte des Eisschelfs (Maximalgeschwindigkeit 935 m pro Jahr), der Geschwindigkeitsabfall im Nord-Süd-Profil und die systematische Zunahme der Divergenz der Strömungslinien gegen die Eisränder. Sorgfáltiges Studium des Geschwindigkeitsvektorfeldes zeigt einige Abweichungen von einer völlig gleichförmigen Verteilung längs des Profiles.

Type
Research Article
Copyright
Copyright © International Glaciological Society 1969

Introduction

Since its discovery by Sir James Clark Ross in 1841, the Ross Ice Shelf, Antarctica, has attracted more and more sciengific attention. But it was not until the International Geophysical Year (I.G.Y.), 1957–58, that Reference CraryCrary and others (1962) began to undertake various glaciological and seismological traverses across this part of the white continent. A series of field operations was initiated by J. H. Zumberge, who after his deformation studies at “Camp Michigan” (Reference ZumbergeZumberge and others, 1960) realized with C. W. M. Swithinbank the importance of the so-called Dawson trailFootnote * for long-term glaciological investigations. His “Ross Ice Shelf Studies” consisted of the determination of ice deformation, movement and snow accumulation along this route. During the second traverse of the route (Ross Ice Shelf Traverse, RIST) in the Antarctic summer 1959–60, Swithinbank was able to measure the heights of 1 800 bamboo poles, and to set out simple deformation patterns each 20 miles (32 km) (Reference Zumberge and MellorZumberge, 1964, p. 71). He also made sun observations at points along the trail. A comparison of the geographical coordinates of these points with the corresponding values of Dawson suggested ice velocities (Reference Zumberge and MellorZumberge, 1964, p. 67) with errors in places twice as big as the ice velocity itself.

Therefore, when Zumberge was planning a repetition of RIST in 1962–63, W. Hofmann proposed improvements, originating from his geodetic field work with Expédition Glaciologique Internationale au Groenland (EGIG) on the Greenland ice sheet in 1959 (Reference HofmannHofmann, 1964). During the Antarctic summer 1962–63, a geodetic traverse was measured by means of modern tellurometers along the Dawson trail (Reference HofmannHofmann, 1963). Particulars of the instruments and their use were published by Reference Hofmann and MellorHofmann and others (1964). Due to the accumulation of observational errors in a geodetic traverse, the accuracy of traverse points decreases rapidly with increasing distance from the starting point. Therefore, traverse angles and distances have to be determined with high precision. Only then, and by complying with some special conditions for observing in polar regions, can individual traverse points be determined with a precision one to two orders of magnitude higher than by astronomical means. The traverse, therefore, has been called the “Ross Ice Shelf Survey” (RISS).

In the course of this RISS-I expedition, 114 marking poles were set out and surveyed along the Dawson trail, and also along a profile roughly parallel to the 168° W. meridian. They consisted of 12 ft (3.6 m) aluminium tubes which projected about 6 ft (1.8 m) above the snow surface. During the traverse, Reference Heap, Rundle and MellorHeap and Rundle (1964) re-measured the 1 800 bamboo poles that had been set up in 1958, thus giving for the first time a complete accumulation profile across the northern part of the Ross Ice Shelf.

The measurements were repeated during the Antarctic summer 1965–66 under the leadership of E. Dorrer. The purpose was to re-measure the RISS-I markers using the same geodetic methods. The differences in distance and angle between the measurements would then give the displacement, and hence the actual movement for each marker. The following sections deal with RISS-II and the data processing of both sets of measurements.

Geodetic Field Work

Personnel

Each expedition consisted of six men (for RISS-I see Reference Hofmann and MellorHofmann and others (1964)). The participants of RISS-II were E. Dorrer as geodesist and field leader, K. Nottarp as electronic specialist, O. Rcinwarth as glaciologist, N. O’Hara as geologist and navigator, W. Seufert as assistant geodesist and D. Stelling as assistant glaciologist.

Instruments

During both traverses, the same surveying methods were used. Experience in Greenland (Reference HofmannHofmann, 1964) showed that a geodetic traverse was most accurate if its distances were measured with tellurometers and its angles with a precision theodolite. Unlike the 1959 journey in Greenland, the aerial system was separated from the tellurometer and mounted on a light aluminium tube approximately 4 m high (Reference NottarpNottarp, 1963). It was therefore possible to measure distances of 8 to 11 km directly, almost independently of local topographie features on the ice shelf. Using experience from RISS-I, special extension legs for the theodolite tripod were constructed by Nottarp, thus enabling the observer to sight his angle targets during RISS-II from a standpoint up to 0.5 m higher. The targets were identical with the tellurometer reflectors (diameter 65 cm). Because of the high mounting of the theodolite, traverse angles could be measured even under abnormal atmospheric conditions. One tellurometer (MA 1-417) was borrowed from the Ohio State University and two others (MR II-3MV, MRA 11-4MV) from the U.S. Army Corps of Engineers at Fort Belvoir, Maryland. A Kern DKM3 precision theodolite had been made available to the party by the Institut fir Photogrammetric and Kartographie at the Technische Hochschule München, Germany, a Zeiss Th 40 theodolite by the Geodätisches Institut of the Technische Hochschule Braunschweig, Germany, and a Kern DKM2 theodolite with astronomical accessories by United States Antarctic Research Programs (USARP).

Equipment

The group was provided with four used Polaris “Sno Traveler” motor toboggans by USARP. One of these motor toboggans broke down at R18 and could not be repaired within reasonable time; it was center behind. Sciengific instruments, tents, food, fuel, aluminium tubes, bamboo poles and personal gear were distributed on eleven Nansen sledges, three or four of which had to be dragged behind each motor toboggan.

Method

To eliminate gross errors in the observations and to check the results, distances as well as angles must be measured at least twice and preferably repeatedly. This requires the arrangement of three groups of two men each, distributed over three neighbouring traverse stations (Fig. 1). The measuring process was directed by means of the voice communication channel in the tellurometer system.

Fig. 1. Arrangement of traverse measurement.

The task of the first group after leaving its preceding station was to find the next marker. Along the Dawson trail, between R1 and R30, most of the old bamboo poles could be seen and followed. Later on, a navigation technique similar to that of Reference BlackBlack (1962) was used. Sitting or standing and facing backwards on the motor toboggan, the driver sighted back along the track of his sledge train and to a fairly straight line of 12 ft (3.6 m) bamboo poles set out each kilometre by his companion. After having noted the bearing to the next station by using the RISS-I angle, the navigator never missed the new marker by more than 50 m. Immediately after arrival, the aerial mast was placed vertically on the marker tube and oriented towards the middle group; that is to say, backward by the first group and forward by the third group.

Meanwhile group two (Fig. 1) had set its theodolite vertically above the centre of the upper end of the aluminum tube. This centre was chosen as the reference point to which all measurements were to be related. Normally, the markers of the first traverse could be used directly. Sometimes, however, these old markers were higher than the tripod. In those cases, an auxiliary marker was used a few metres away, the eccentrically measured angles and distances later being reduced to the old marker. Each traverse angle was measured in five full sets as quickly as possible in order to avoid the effects of the tripod sinking into the snow and tripod torsion due to the sun’s radiation. In favourable observational conditions, the angle measurement of a single point could be completed in 15 to 20 min.

For the distance measurements, the theodolite and tripod had to be replaced by the tellurometer aerial mast. From each station, group two measured the distance to the forward and to the backward station with its master tellurometer (MA I-417). Ten fine readings were made between the initial and final coarse readings. It was also necessary to make meteorological observations and to measure the inclination of the mast. During the course of the work all traverse distances were measured twice, giving an important check against gross errors. No single distance measurement took more than 15 min.

On each traverse it took about 2 months to measure the 103 traverse stations on the Dawson trail and on the north–south profile (910 km). This gives a rough figure of the rate of measurement using light-weight equipment on a smooth level surface without crevasses such as the northern part of the Ross Ice Shelf.

Accuracy

From a total of 106 angle measurements during the second traverse, a mean square error of ±±1·9ec (dispersion range 0.3 to 5.2ec) was obtained (Fig. 2). Compared with the corresponding value of the first journey (±2.3ec), a remarkable increase in accuracy was noted, due to the experience gained on RISS-I.

Fig. 2. Frequency distribution of the standard errors of traverse angles.

Using 101 distances measured back and forth, a mean square error of ±6.7/2 = ±3.4 cm (the same as in RISS-I) resulted for every traverse distance measured twice. As indicated in Figure 3, a systematic error of 4.3/2 = 2.2 cm with undetermined sign reveals a possible constant deviation between the two remote tellurometers (Fig. 1), probably due to deformation of their antennae. During RISS-I, no similar effect had occurred.

Fig. 3. Frequency distribution of the distance differences between fiirward and backward measurements.

Determination of Surface Movement

Reduction of data

At first, all data observed during both field journeys had to be reduced with respect to deviations from the geometrical model adopted. The transmission time of the electromagnetic waves between two stations depends on the refractive index of air; changes in the pattern frequency of the tellurometer influence the scale of distance measurement. Therefore, meteorological reductions due to air pressure, humidity and temperature (Reference Hofmann and MellorHofmann and others, 1964), and frequency corrections are necessary. Due to the fact that the tellurometer antennae were not exactly vertical, all traverse distances and angles had to be reduced to the reference point defined above.

For each traverse distance measured twice, the arithmetic mean value was calculated. Since both observations were carried out at different instants, this is only correct to a degree depending on the amount of ice deformation. Theoretically, the mean of both distances and the mean of both corresponding observation instants are incoherent quantities.

All angles and distances observed eccentrically on auxiliary markers had to be reduced to the corresponding central marker. Since geographical coordinates refer to an earth ellipsoid coinciding approximately with sea-level, all traverse distances on the ice-shelf surface had to be reduced to sea-level. The corresponding altitudes were taken from Reference CraryCrary and others (1962) and the altitude of Observation Hill is based on the H.O. Chart No. 6712.

Methods

Normally, geodetic measurements are carried out on the rigid earth surface, whose points have fixed relative and absolute positions. However, for measurements over a glacier, the displacement of the surveying point due to the deformation of the surface during the measurement has to be considered. This means that any geometrical observation (angle, azimuth or distance) is valid only for a certain instant (Reference DorrerDorrer, 1967). Measurements combined in the same survey have to be reduced to the same reference time. Deformations and movements in the survey area can be detected and determined by repeated measurements at different times.

The “reduction to epoch” problem has been treated by Reference SwithinbankSwithinbank (1958). Reference HofmannHofmann (1964) has also discussed the reduction of all observations in a glacial survey to a reference time. As small parts of a glacier do not generally change their form very much within short periods, this problem was neglected by Dittrich and Schwarz (196Reference Dittrich and Schwarz6) for the loo km long trigonometric chain near the Soviet station Mirny in Antarctica. The geodetic traverse over the Ross Ice Shelf, however, extends over a distance of more than 900 km and is tied to a fixed point on rock only at its beginning (it is a so-called “end-free traverse”). Owing to the cumulative propagation of errors of traverse angles, time reduction in this case is essential.

If only two measurements are available, the reduction must be based on the assumption that the movement of the markers took place along straight lines and with constant speed. Neither the curvature of the flow lines nor an acceleration of the trajectory can be derived from such a solution. At the ice margin, however, flow lines curve as the ice conforms to the outline of the land. The start of the RISS traverse (cf. Fig. 11) lies on the McMurdo Ice Shelf in a zone subject to rapid local deformation. Neglecting this would affect at least the directions of the traverse. It is therefore worth considering whether other glaciological surveys in the area might throw light on this problem.

Assuming the traverse points to be moving along straight lines and at constant speed, the time reduction can be made in two distinct ways:

  1. The field observations of the two journeys can be reduced to two reference dates in order to determine the displacements in the period between them. The displacement of all points during the period is computed by comparison between the two series of time-reduced observations. The traverses I and II, reduced to their reference dates t z and t z are represented in Figure 4 by thin lines.

  2. The original observations can he used if applied to points whose positions correspond to the actual dates of measurement. These positions are obtained by time reduction along the observed displacement vector. In this case the displacement of all points results from the comparison of two terraced traverses (Fig. 4).

Both methods, if correctly applied, lead to the same final displacement vectors for a certain period, normally one year. As they have both been used on the RISS traverse, they will be described in detail.

Fig. 4. Terraced traverses caused by time dependency of measurements.

Time reduction of observations (W.H. and W.S.)

The fundamental equations can be derived from Figure 5 in which o is a fixed point whereas p moves along the straight line p′ p″ with constant velocity. The traverse angles ω′, ω″ and distances s′, s″ were measured at times t′, t″, respectively. The displacement vector for the time interval (t″−t′) is r. The displacement vectors r t , or r t , for any time intervals (tt′) or (t″−t) are determined by linear interpolation:

(1a)

(1b)

Fig. 5. Time-reduction of observations.

With the denominations in Figure 5 the corresponding values ω t and s t can be determined with the following formulae:

(2)

(3a)

(3b)

or by exchanging s′ with s″ and c′ with c

(4a)

(4b)

(5)

Formulae (3)(5) show that the time reduction of the observations ω and s depends only on the observations themselves and their changes with time; it can be done without knowledge of the displacement vector r.

In order to reduce errors caused by deviations from the assumed linear movement of p, formulae (3a) and (3b) will be used if (tt′) is small, i.e. for the first traverse, whereas for the same reason formulae (4a) and (4b) are suitable for the second traverse where (t″−t) is small.

For the application of formulae (3)(5) along the whole traverse, Figure 5 has to be adjusted to the general case where not only p but also o is movable.

Figure 6 shows the straight flow lines of three adjacent points p i−1, p i , p i+1 of the traverse. The positions p′ refer to the first traverse, p″ to the second. The positions p′ i−1 and p″ i−1 represent the case that all measurements were executed at the same moments t i−1 and t i−1 on the two traverses.

Fig. 6. Time-reduction of observations.

By parallel shift of the distance

into the position

we originate the triangle

.

is the vector of the relative displacement between p i−1 and p i in the time interval (t i−1t i−1). As all movements are taken to be straight and proportional to time, the length of this vector is proportional to the absolute displacement

. Therefore, triangle

corresponds to the reduction triangle op′ p″ in Figure 5 and the time reductions of the traverse elements ω and s can be derived in this triangle. The angle

is the sum of the differences ω′−ω″ of the traverse angles from ω 0 to ω i−1 between the two traverses; for the present, it may be taken as known.

If the measurements in p′ i and p″ i took place later at t i and t i the markers will have moved to the positions p′ i and p″ i . The corresponding reduction triangles are

and

. The triangle

permits the time reduction of the measurements in

and

to the times t i−1 and t i−1 in the following way.

The angle

differs from

by the difference of the measured traverse angles ω i and

:

(6)

As Figure 6 shows, angle

is connected with angle

by the equation:

(7)

The generally small angles δ i−1 and δ i−1 transform the reduction triangle

into the triangle

and can be computed in the triangle

by application of Equations (3a) and (4a), using the known time intervals (t i t i−1), (t i t i−1) and (t i−1t i−1) between the measurements in p i−1, and p i . Combining Equation (7) with (6) gives

(8)

With

, all necessary elements of the reduction triangle

are known.

The time reduction of triangle

to triangle

according to formulae (3) and (4) gives the reduced distances

and

and the generally small angles d i and d i .

Finally, as Figure 6 shows, the reduced traverse angles

and

can be computed with the formulae:

(8a)

(8b)

In this way the traverse elements—angles ω and distances s—of the leg p i p i+1, measured at time t i , can be reduced to any time t 0 if the corresponding elements of the preceding leg p i−1p i are known for the time t i−1. As the elements of the first triangle p0 p′1 p″1 between the starting point p0 (Observation Hill) and the first marker p1 are known as measurements in p0 at times t0 and t0, all following measurements can be reduced to the dates t0 and t0. For that purpose, the angles d i and d i and the reduced distances

and

have to be computed with the time intervals (t i t0) and t i t0).

The displacement vectors follow from a comparison of the marker coordinates in both traverses.

Time reduction of positions (E.D.)

The use of original observations requires that they are applied at the points where they were actually measured. The position of points on a surface can be described by differential geometry. The mathematical relations are simple if a reference plane is used. If, however, we consider that all movements take place on the earth ellipsoid (ellipsoid of revolution), the reduction becomes rather complicated. Any points (markers) p on the surface of this ellipsoid can be indicated by their geographical coordinates p = (λ,ϕ), where λ = longitude (positive eastward, negative westward), and ϕ = latitude (positive in the Northern Hemisphere, negative in the Southern Hemisphere).Footnote *

If at instants t′, t″ a marker p is located at the corresponding positions p′ = p(t′), p″ = p(t″), the position of the ice flow line is determined but not its shape. On this line the position of the marker itself is fixed at an arbitrary instant t if the velocity vector v is known.

The location of a marker p as a function of time can be expanded into a Taylor series in powers of (tt 0):

(9)

where

and

are the velocity and acceleration vectors respectively at a reference time t 0. Applying t′ and t″ to Equation (9), subtracting the resulting equations from each other, and putting t 0 = (t′+t″)/2 gives

(10)

Neglecting all residual terms of higher than second order leads to an average velocity vector

(11)

Since no acceleration vector can be deduced from Equation (10), and therefore no force system assumed to act upon the marker, the flow line has to coincide with the geodesic between p′ and p″. Denoting the metric tensor of the surface by G, and the magnitude of the velocity by v, then

where

(12)

represents the geodetic distance between the two points p′, p″. The variable surface vector p has to meet the differential equations of the geodesic. This means that Equation (12) is a differential equation which may be solved by any numerical method (see for example Reference DorrerDorrer (1966)).

The reduction procedure will be described by means of Figure 7, which shows the beginning of a terraced traverse repeated after a certain time interval. Starting from a fixed point p0 and a fixed azimuth from p0 to a second fixed point p−1, the positions of the first marker p1 can be computed at instants t0, t0 at which the angles ω0, ω0 and distances s1, s1, respectively, have been measured.

Fig. 7. Geometrical relations at the beginning of a repeated traverse on moving ground (simplified).

On the earth ellipsoid this computation leads to the solution of two systems of differential equations of the geodesics between p0 and

,

, respectively. This is known in geodetic literature as the “problem of transferring geographical coordinates over long distances” (in the following denoted by “problem I”), for the numerical treatment of which any appropriate method may be used (see for example Reference DorrerDorrer (1964)).Footnote The displacement vector r 1 is now fixed between

and

, Footnote and can be determined from Equation (12) which represents the reversal of problem I, viz. the problem of determining distance and azimuth of a geodesic between two known points (problem II Footnote §).

Since measurements of the angles ω1, ω1 in p1 and distances s2, s2 have been carried out at different instants t1, t1, respectively, the corresponding positions

,

have to be determined on the geodesic through

and

. The small displacements from

,

to

,

, denoted by

,

, respectively, are defined by

(13)

,

can therefore be found by using problem I. To transfer the azimuths between p0 and

,

to the azimuths between

,

and

,

, respectively, at first problem II has to be applied to the geodesics between p0 and

,

respectively.

The two resulting azimuths

,

simply lead to the tranferred ones:

(14a)

(14b)

By means of these azimuths and the distances s2, s2 the positions

and

of marker p2 can be determined by problem I, and the procedure described above has to be repeated. From marker p2 on, however, the transfer of azimuths is possible only if the corresponding positions

,

of marker p1 were computed as analogues to

,

, respectively. The corresponding displacements

,

Footnote * can be determined by using Equation (14):

(15)

and problem I.

If required, the positions of all markers at their different observation times (two terraced traverses in Figure 4) can now be reduced to a certain reference date which may be chosen arbitrarily unless there are important reasons for a particular time. For the RISS traverse, this applies in the McMurdo Ice Shelf region where great changes of ice velocity occur. Therefore the reference date should be chosen within the period of observations in this region. If this is done, deviations of the supposed straight flow lines from the actual ones will have least effect on the whole traverse. For the sake of simplicity, the dates t z = t0 and t z = t0 of observations on Observation Hill (p0) have been taken as reference dates for both traverses. Reducing the markers to these dates analogue to Equations (13) or (15) leads to the positions

,

, i = 1, 2, …, n (Fig. 7) of the two reduced traverses in Figure 4. The geographical coordinates of all traverse markers, listed in the Appendix, are identical with these reduced positions.

The process of computation from a marker p i−1 to the next marker p i requires the solution of problem I six times, and of problem II three times, and, if reduced traverses are wanted, an additional two solutions of problem I.

As the velocity vectors still refer to sea-level, they have to be transferred in reverse to the ice-shelf surface.

Improvement of the method. The method described above can be refined with respect to time reduction and with respect to the acceleration and curvature of the ice movement.

Since at a marker p i−1 the traverse angle ω i−1 and the traverse distance s i towards the next marker p i , i = 1, 2, …, n were normally measured at slightly different times, the two markers changed their positions relative to each other during the time interval. The computation of p i by azimuth and distance (problem I) is therefore theoretically incorrect unless s i is reduced to the date of angle measurement. Figure 8 illustrates the geometrical conditions of the traverse between neighbouring stations.

Fig. 8. Measurement conditions between neighbouring stations of a traverse on moving ground.

At a certain time t i−1, the angle ω i−1 was measured in

. The azimuth

,i towards

can then be fixed provided the position of the preceding marker

is known at the same time (connection azimuth

, i−2). Hence, p i must lie on the geodetic line that goes through

with an azimuth

,i. At another instant t i−1;s , distance

was measured between

and

. Since the movement vectors of p i−1 and p i are not parallel, and since t i−1;s t i−1, the distances s i and

are not equal. The new position of the marker

could be determined exactly by transferring the geographical coordinates

into

over the distance

(problem I), if

were known. This distance, however, can only be obtained with mathematical rigour from the solution of an algebraic equation of fourth degree, which numerically is inconvenient. Although

and

could be determined in the case of RISS with sufficient accuracy by assuming a linear dependence of the distance with respect to time, a more general methodFootnote * may be considered.

Figure 9 illustrates how

and

can be obtained by means of successive approximations, provided that the time interval t i−1;s t i−1 is much smaller than the time difference between both traverses. Proceeding on the assumption that p i−1 is known as a function of time, the positions

,

and

,

of the marker are fixed (Equation (15)). By means of the known azimuths

,

,

, i and distances

,

approximate positions

,

for

,

, respectively, can be transferred from

,

(problem I). The greater the difference

, the more different are the exact and approximate positions. This leads (problem II) to an approximate displacement vector

, and to a velocity vector

according to Equation (II). Approximate positions

,

can be determined from Equation (9):

Fig. 9. Iterative solution for the elimination of time differences between consecutive angle and distance measurements.

Now, the geodetic distances

,

, obtained by problem II, have to be compared with the corresponding measured distances s i , s i . From the differences

,

thus found, the corresponding differences

,

can be found by linear interpolation from

Addition of these differences to the measured distances provides better approximations for

,

. The procedure has to be repeated until the differences δ s i and δ s i are equal to zero. The convergence of the recursive process is evident from Figure 9. At the end, the final marker positions

,

, as well as a definite velocity vector v i , will be obtained.

The second improvement considers that, if the curvature of the flow lines and the change of velocity are known at least approximately, the accuracy of the results can be greatly increased. Especially at the beginning of the traverse in the McMurdo Ice Shelf area, the flow lines depart considerably from straight lines (see Fig. 11). As such data could be obtained from a plot of velocity vectors on the McMurdo Ice Shclf,Footnote it was worth considering them in the computations. However, reliable values for curvature and acceleration could be determined only for the first three RISS markers (Table I).

Table I Additional Values for Some Flow Lines

Knowledge of the curvature c allows us to replace an assumed straight line by an are of a circle with radius c −1. The straight flow lines shown in the Figures 7, 8 and 9 can therefore be replaced by curves. It is evident that the mathematical relations become still more complicated, even if we neglect earth curvature and earth rotation (Coriolis force).

From Figure 10 there follows a relationship between the length of the are r c and the geodetic length of the chord r g:

(16)

Fig. 10. Flow line with constant curvature.

The whole computation procedure has to be altered with these refinements, i.e. if the positions

,

of a marker p i are determined (Figs. 7, 8 and 9), the geodetic length

II not only gives the arc length rc by Equation (16) but also the angle

(17)

between chord and arc at the endpoints

,

. The velocity v at these two points results from Equation (9), slightly changed however from its vectorial form into a representation along the flow line. Denoting the arc length by r and the tangent acceleration by a, the differentiation of

with respect to t then leads to (v 0 from Equation (11))

hence to

(18a)

(18b)

Thus the velocity vectors v i , v i are determined by Equations (17) and (18), and the positions

,

,

,

(Fig. 7), or additionally

,

,

,

(Fig. 8) can be computed by a procedure similar to that described above.

The effect of using these external data becomes evident if we consider the values for c and a of the RISS marker R2 in Table I. By Equations (17) and (18), an angle γ = 2° 02′ and a velocity difference v″−v′ = −2.2 m year−1 can be deduced.

Discussion of results

Application of the method described above by means of the reduced field data and Reference HeineHeine’s (1967) data (Table I) gives the results listed in the AppendixFootnote * and represented graphically in Figure 11.Footnote This figure shows the field of velocity vectors along the RISS traverse. Obvious characteristics are the rapid increase of velocity between McMurdo station and the Ross Ice Shelf, the uniform, nearly parallel movement in the middle of the ice shelf between the markers R33 and R37, and the systematically increasing divergence of the flow lines towards the ice margins.

Fig. 11. Map of the northern part of the Ross Ice Shelf, Antarctica,.showing the positions of the RISS aluminium mackers. the movement vectors (in year−1) and the strain-rates (10−2 year−1).

On the southern leg, a distinct decrease of velocity can be noted as well as a clear change of movement direction due to the turning of the ice flow lines around the western side of Roosevelt Island. The uniform decrease of velocity and the increasing divergence of movement direction along the nearly straight part of the traverse between r57 and r69 makes it possible to extrapolate these values beyond r69 to r76. We find that the region around r77 moves at a velocity of about 300 m year−1 almost perpendicular to the ice front west of the Bay of Whales.

Moreover, Figure 11 gives a graphical representation of the strain-rates between each pair of adjacent points on the traverse. These values, which are also listed in the appendix (column 6), give a measure of the relative change of distance per year. Two characteristics are evident. First, a very homogeneous strain field along the Dawson trail with significant variations only at the traverse vertices. Secondly, a high dependence of the strain-rates on the azimuth of the traverse, especially on the southern part of the north–south leg. This is evidently due to the compression caused by a major ice stream coming from Marie Byrd Land and another flowing from the polar plateau through the eastern part of the Queen Maud Range (Reference GiovinettoGiovinetto and others, 1966).

Owing to the zigzag course of the north–south profile, due to primitive navigation techniques during the first traverse (Reference Hofmann and MellorHofmann and others, 1964), every station p i , the traverse angle ω i of which differs considerably from 180° (Fig. 12), allows us to determine the deformation tensor. For this purpose, the two adjacent strain-rates

,

and the rate of angle change

are sufficient, although the accuracy of the eigenvalues and eigenvectors decrease rapidly if ω i approaches 180°. The maximum and minimum strain-rates, together with their corresponding azimuths, have been computed for many points on the north–south profile but also for a few points on the Dawson trail. It must be noted, however, that these values appear to be significant only for the stations r5, r17, r53, r57 and r69. Though many traverse angles on the north–south profile differ more than 10° from 180°, the errors suggest a greater accuracy than can be expected. The conspicuous major changes of strain-rate (Appendix, column 6) and angle would require much more detailed information about the deformation field than is available from our observations. But this was beyond the scope of RISS.

Fig. 12. Adjacent strain-rates and angle change-rate at a traverse point Pi.

Although a careful study of the velocity vectors shows some deviations from an entirely uniform distribution along the four straight sections of the traverse between r5 and r69, velocity changes reveal a more detailed picture. The changes of velocity vectors along all traverse distances (except between r1 and r5), reduced to unit length, are shown in Figure 13. These changes are equivalent to the relative movement vectors between consecutive markers and they depend mainly on the corresponding direction of progress in the velocity vector field. This is evident at the traverse vertices r17, r53, r57, r69 and in the whole north–south profile. For each velocity-change vector the component parallel to its traverse distance is identical to the strain-rate, shown in Figure 11. Notable points lie between r8 and r10, where the border effect of Minna Bluff and Ross Island probably comes to an end, and between r120 and r133.

Fig. 13. Velocity vector changes along the traverse and reduced to unit length.

The velocity-change vectors around the traverse vertices, where the traverse angles differ considerably from 180°, enable us to determine the acceleration and curvature of the ice flow line at these points. The results are shown in Table II. According to this, r5 is slowed down considerably, but r17 still has a negative acceleration.

Table II Curvature and Acceleration Values for the Traverse Vertices

Conclusions

In spite of the geodetic character of the field work and data processing in this study, the results obtained are of major glaciological significance. The small expenditure in terms of personnel, equipment and money shows the efficiency of an operation of this kind extending over two Antarctic summer seasons. Use of the same method over inland ice sheets rather than ice shelves would of course need a different type of organization and different vehicles. Although further simplification might be considered (for instance, a reduction of personnel from three groups to two groups of two men each), such changes would be possible only at the expense of accuracy and of reliability of the observations and their results. To achieve a greater accuracy also for points farthest from the fixed starting points, there are two possibilities: (1) Measurement of auxiliary azimuths systematically distributed over the whole traverse, and (2) tying both ends of the traverse to fixed points. The rather sophisticated surveying technique used requires at least one specialist in electronics and one for surveying. In order to encourage further movement studies on the Ross Ice Shelf, all RISS markers were extended to a height above the snow level which ought to keep them visible at least until the 1969–70 summer.

Acknowledgements

The field work and data analysis were supported by grants from the National Science Foundation to the University of Michigan and to Grand Valley State College. President J. H. Zumberge was Principal Investigator. Logistic support was provided by Air Development Squadron Six of the U.S. Navy. Help and advice by Professor Dr E. Gotthardt, München, is gratefully acknowledged. Individuals whose help in the field work or in data reduction contributed to this publication are: K. Nottarp, J. A. Heap, A. S. Rundle, W. C. Campbell, O. Reinwarth, N. O’Hara, D. Stelling, A. J. Heine, D. Weber, M. Stephani. The authors wish to express their thanks to C. W. M. Swithinbank for carefully reading and improving this manuscript, and for his suggestions.

Appendix

Appendix

Footnotes

* In December 1958, Major M. Dawson had led a U.S. Army tractor train from “Little America V” almost to McMurdo, and had marked the route by 3 m high bamboo stakes about 350 m apart, and by empty fuel-drum cairns every 20 miles (32 km).

* p is a surface vector and stands for (λ,ϕ).

In the German language: “Erste geodätische Hauptaufgabe”.

The following explanations deal with subscripted variables. Their lower index always denotes the marker’s number in the traverse, the upper one distinguishes between several instants of observation.

§ In the German language: “Zweite geodätische Hauptaufgabe”.

*

may be denoted as the geodetic norm of a vector x.

* If, for instance, curvature and acceleration of the flow lines are considered.

This plot is the result of an intensive survey of a 10 km marker grid carried out during several Antarctic summers by New Zealand glaciologists (Heine, 1967).

* All computations were carried out by a TR 4 computer of the Kommission für Elektronisches Rechnen at the Bayerische Akademie der Wissenschaften, München.

Traverse points, velocity vectors and strain-rate “rectangles” have been plotted automatically by a ZUSE Z64 graphomat.

References

Black, H. P. 1962. A navigation device for steering vehicles in featureless terrain. Polar Record, Vol. 11, No. 72 p. 29698.Google Scholar
Crary, A. P., and others. 1962. Glaciological studies of the Ross Ice Shelf, Antarctica, 1957–1960, by A. P. Crary, E. S. Robinson, H. F. Bennett and W. W. Boyd, Jr. IGϒ Glaciological Report Series (New York), No. 6.Google Scholar
Dittrich, G. Schwarz, G. 1966. Die geodätischen Arbeiten der deutschen Gruppe während der 7. sowjetischen Antarktis-Expedition 1962. Nationalkomitee für Geodäsie und Geophysik der Deutschen Demokratischen Republik bei der Deutschen Akademie der Wissenschaften zu Berlin. Fachgruppe Geodäsie, Reihe 3, Ht. 5.Google Scholar
Dorrer, E. 1964. Algorithm 240 coordinates on an ellipsoid (Z). Communications of the Association for Computing Machinery, Vol. 7, No. 9, p. 546.Google Scholar
Dorrer, E. 1966. Direkte numerische Lösung der geodätischen Hauptaufgaben auf Rotationsflächen. Deutsche Geodätische Kommission bei der Bayerischen Akademie der Wissenschaften, Reihe C, Nr. 90.Google Scholar
Dorrer, E. 1967. Die Bestimmung der Oberflächengeschwindigkeit ausgedehnter Gletschergebiete. Ein Beitrag zur Erfassung horizontaler Bodenbewegungen mit geodätischen Methoden. Zeitschrift für Vermessungswesen, 1967, p. 18389.Google Scholar
Giovinetto, M. B., and others. 1966. The regime of the western part of the Ross Ice Shelf drainage system, by M. [B.J Giovinetto, E. S. Robinson and C. W. M. Swithinbank. Journal of Glaciology, Vol. 6, No. 43, p. 5568.Google Scholar
Heap, J. A. Rundle, A. S. 1964. Snow accumulation on the Ross Ice Shelf, Antarctica. (In Mellor, M., ed. Antarctic snow and ice studies. Washington, D.C., American Geophysical Union, p. 11925. (Antarctic Research Series, Vol. 2.))Google Scholar
Heine, A. J. 1967. The McMurdo Ice Shelf, Antarctica: a preliminary report. New Zealand Journal of Geology and Geophysics, Vol. 10, No. 2, p. 47478.Google Scholar
Hofmann, W. 1963. Geodätisch-glaziologische Arbeiten in der Antarktis. Zeitschrift für Vermessungswesen, 1963, p. 25565.Google Scholar
Hofmann, W. 1964. Die geodätische Lagemessung über das grönländische Inlandeis der Internationalen Glaziologischen Grönland-Expedition (EGIG) 1959. Meddelelser om Gronland, Bd. 173, Nr. 6.Google Scholar
Hofmann, W., and others. 1964. The Ross Ice Shelf Survey (RISS), 1962–1963, by W. Hofmann, E. Dorrer and K. Nottarp. (In Mellor, M., ed. Antarctic snow and ice studies. Washington, D.C., American Geophysical Union, p. 83118. (Antarctic Research Series, Vol. 2.))Google Scholar
Nottarp, K. 1963. Hochantennen für Tellurometermessungen in der Antarktis. Allgemeine Vermessungsnachrichten. 1963, Ht. 8, p. 30407.Google Scholar
Swithinbank, C. W. M., 1958. Glaciology. I. The movement of the ice shelf at Maudheim. Norwegian–British–Swedish Antarctic Expedition, 1949–52. Scientific Results, Vol. 3, C.Google Scholar
Zumberge, J. H. 1964. Horizontal strain and absolute movement of the Ross Ice Shelf between Ross Island and Roosevelt Island, Antarctica, 1958–1963. (In Mellor, M., ed. Antarctic snow and ice studies. Washington, D.C., American Geophysical Union, p. 6581. (Antarctic Research Series, Vol. 2.))Google Scholar
Zumberge, J. H., and others. 1960. Deformation of the Ross Ice Shelf near the Bay of Whales, Antarctica, by J. H. Zumberge, M. [B.] Giovinetto, R. Kehle and J. [R.] Reid. IGϒ Glaciological Report Series (New York), No. 3.Google Scholar
Figure 0

Fig. 1. Arrangement of traverse measurement.

Figure 1

Fig. 2. Frequency distribution of the standard errors of traverse angles.

Figure 2

Fig. 3. Frequency distribution of the distance differences between fiirward and backward measurements.

Figure 3

Fig. 4. Terraced traverses caused by time dependency of measurements.

Figure 4

Fig. 5. Time-reduction of observations.

Figure 5

Fig. 6. Time-reduction of observations.

Figure 6

Fig. 7. Geometrical relations at the beginning of a repeated traverse on moving ground (simplified).

Figure 7

Fig. 8. Measurement conditions between neighbouring stations of a traverse on moving ground.

Figure 8

Fig. 9. Iterative solution for the elimination of time differences between consecutive angle and distance measurements.

Figure 9

Table I Additional Values for Some Flow Lines

Figure 10

Fig. 10. Flow line with constant curvature.

Figure 11

Fig. 11. Map of the northern part of the Ross Ice Shelf, Antarctica,.showing the positions of the RISS aluminium mackers. the movement vectors (in year−1) and the strain-rates (10−2 year−1).

Figure 12

Fig. 12. Adjacent strain-rates and angle change-rate at a traverse point Pi.

Figure 13

Fig. 13. Velocity vector changes along the traverse and reduced to unit length.

Figure 14

Table II Curvature and Acceleration Values for the Traverse Vertices