## Introduction

Renewed interest in the detailed movements of icebergs began with the advent of oil exploration on the continental shelf off eastern Canada, where their occurrence presents considerable hazards to drilling operations and shipping.

On the other hand, icebergs may possibly be a source of fresh water if they can be towed from the Antarctic to arid parts of the world. Whether an iceberg is to be used or avoided, it is important to learn about its drift under the influence of environmental forces such as winds and currents.

In this paper, an iceberg drift model is presented which considers the motion of an iceberg due to wind and water drag forces, the Coriolis force, the pressure-gradient force resulting from the sea-surface slope, and the water acceleration. The validity of this model can only be studied by comparing predicted and observed iceberg drift trajectories when data on wind and currents are collected simultaneously and used as input to the model.

Such an opportunity presented itself during a storm on 21–22 August 1972, near Saglek, Labrador, where the Memorial University of Newfoundland (MUN) was conducting an oceanographic study in collaboration with the Bedford Institute of Oceanography (BIO). During the study, iceberg positions were monitored with the help of a land-based radar at Saglek, while current and wind data were gathered by CSS *Dawson*, provided by BIO. The predicted iceberg drift trajectory was close to that observed for a specific iceberg in this area. A detailed account of activities during this study is given by Allen (1972).

## The Iceberg Drift Model

Iceberg motion is the net result of forces which vary with time and space. Some of these forces are due to gravity, pressure gradient, wind drag, water drag, Coriolis effects, waves, and swells. Since this study is mainly concerned with the horizontal movement of icebergs, only the horizontal components of these forces are considered. Wave and swell forces are neglected in this study as their magnitudes in the horizontal directions are unknown in comparison to other forces. The mathematical model takes into account the significant environmental forces due to water drag, wind drag, Coriolis acceleration, and sea-surface slope (pressure gradient).

The drag force due to the water is proportional to the square of the relative velocity of water with respect to the iceberg. The constant of proportionality depends upon the size and shape of the under-water portion of the iceberg. Ocean currents represent a complicated mixture of different types of aperiodic and periodic water movements, a few of which are the geostrophic current, the wind-driven current, the inertia current, and the tidal current. The distribution of the magnitude and direction of the various current components varies with time, depth, and location. For the purpose of modelling iceberg drift, the water column is divided into two layers, and the water drag force is then obtained as the vectorial sum of the drag forces which are expressed in terms of relative velocity of the current with respect to the iceberg in each layer.

Similarly, the magnitude and direction of the wind drag force depends on the size and shape of the above-water portion of the iceberg. The average ratio of the iceberg velocity to wind speed is about 0.03–0.04 (Murray 1969); thus in the expression for wind drag force, the relative velocity of wind with respect to the iceberg is taken to be the wind velocity itself.

The Coriolis force, due to the rotating frame of reference with the Earth, tends to move the iceberg and the water surrounding the iceberg to the right of their path (clockwise) in the northern hemisphere. In a geostrophic current, the pressure-gradient force on an iceberg due to a sloping sea surface balances the Coriolis force due to its movement. If the iceberg motion is not along a geostrophic current direction, there are two forces acting on the iceberg; the Coriolis force due to its movement and the pressure gradient due to the sea-surface slope. In the present study, the pressure-gradient force in each layer is expressed as the negative of the current’s Coriolis force (Neumann 1968, Mountain 1980). Thus the combined action of the Coriolis force and the pressure-gradient force on the iceberg is equivalent to expressing the Coriolis force of the iceberg in terms of the relative velocity of the iceberg with respect to the farfield current in a particular layer.

If the water around the iceberg is accelerating due to some potential field, the same potential field would also be acting on the iceberg to accelerate it. So the momentum equation must include a term which takes into account the forces accelerating the water mass and the iceberg at the same time, and this force on the iceberg will be equal to the product of mass of the iceberg and the acceleration of water surrounding the iceberg. This was first pointed out by Bayly (1971) in a discussion of a paper by Cochkanoff and others (1971) (also see Napoleoni, unpublished).

The concept of added water mass is introduced to account for inertial drag due to the acceleration of the water around an object, and the object behaves as if a mass were added to it. The added mass can be determined from the potential flow theory, and it has been assumed to be equal to half of the mass of the iceberg.

The equations of motion taking all the above-mentioned forces into account are written below in the component form:

where the symbols used are defined above.

If the parameters related to the iceberg are known, the above set of equations may be integrated to obtain the response of the model in the form of the iceberg velocity and position, provided the current and wind data are supplied as the forcing function (or input to the model). Since the set of differential equations are coupled and nonlinear, it is expeditious to integrate them with the help of a digital computer, and the fourth-order Runge Kutta method is employed for this purpose.

## The Iceberg Parameters

Since the mass, area, and drag coefficients of the icebergs under study are not known, these values are chosen from the range quoted in the literature. Thus, the cross-sectional area and mass data for the iceberg under study are given the values shown in Table 1, which was compiled by the U.S. Coast Guard from observations made over a number of years to establish seven classes of iceberg (Mountain 1980), and the values of the drag coefficients in water and air are assumed to be 1 and 1.5, respectively, regardless of the size and shape of the iceberg.

Table I Table I. Iceberg Characteristics (Mountain 1980 with acknowledgement to R.Q. Robe)

A more detailed discussion of this can be found in E1-Tahan (1980).

## Input to the Model

### Current data

As mentioned above, the icebergs were tracked by radar at Saglek while the currents and wind speeds were recorded simultaneously by current meters and C S S *Dawson*, respectively. Figure 1 shows the trajectories of five icebergs and the locations of the current meters A,B,C, and D. At each of these four locations, current meters were installed at depths of 13 m (A1,B1,CI,D1) and 75m (A2,B2,C2,D2) below the sea surface, and three additional meters were installed at locations A, B, and C at depths of 165 m (A3), 146 m (B3), and 176 m (C3), respectively. The data obtained from these current meters were presented and analysed by Holden (1974), who concluded that the current near the surface oscillated with the tidal period of 12.5 h before the storm of 22 August 1972 and with the inertial oscillation period after the storm.

Fig. 1 Drift pattern of icebergs offshore Saglek after the storm in August 1972, showing locations of current meters.

We shall pay attention to the iceberg drift trajectory 20C because it is close to the current meters at locations A and B. Since the iceberg drift model needs the average water velocity components in the two layers, the oceanic boundary layer (*U*
_{1}, *V*
_{1}) and the deeper water below (*U*
_{2}, *V*
_{2}), the following procedure is adopted to obtain the average velocities from the data of the current meters. The average current velocity of the deeper layer is taken to be equal to the current velocities measured by current meters at 75 m depth. All the current meters at 75 m depth recorded only the direction of the current and not the magnitude due to malfunction, and hence the magnitudes measured by current meters A3 and B3 at deeper depths are assumed to be the same as those at 75 m depth. This is justified because the geostrophic current magnitude is more or less constant below the oceanic boundary layer. Let us designate the components of the measured current at 75 m depth to be *U*_{g}
and *V*_{g}
, and the corresponding components at 13 m depth to be *U*
_{s} and *V*
_{s}. From the above discussion, we have

and

To calculate the average velocity component (*U*
_{1},*V*
_{1}) in the boundary layer, we assume that the wind-generated currents (*U*
_{w},*V*
_{w}) have the distribution with respect to depth in the form given by Ekman (Neumann 1968) . If the wind shearstress vector is acting in the positive *y* direction (i.e. to the north), the wind-generated velocity components are given below:

and

where

The latitude of Saglek is 58°9’ N. and the depth of the boundary layer *D* is calculated to be 39.7 m which is approximately equal to the depth of the mixing layer as is evident from the contour lines of measured salinity-temperature depth data (Allen 1972) . The net water mass transport, *S*_{x}
and *S*_{y}
, in the easterly and northerly directions are given by Neumann (1968) as :

and

This is a remarkable result for it shows that the total water-mass transport is to the right of wind shear-stress direction in the northern hemisphere and to the left in the southern hemisphere. The net effect of wind-driven currents on the motion of an iceberg is to integrate the drag forces at different levels in the boundary layer. In this paper, it is assumed that the boundary layer has a uniform velocity such that the net water-mass transport due to wind effects is equal to that given by Equations 9 and 10, and the average velocity components are *Ū*
_{w} = 0.225 *V*
_{0} and *V̄*
_{w} = 0. Using Equations 7 and 8, the velocity components at 13 m depth are *U*
_{w} = 0.358 *V*
_{0} cos (-13.9°) and *V*
_{w} = 0.358 *V*
_{0} sin(-l3.9°). Hence, we obtain a factor, equal to 0.628, which is the ratio of the magnitude of the average velocity in the boundary layer to that at 13 m depth, and the difference in their directions is 13.9° (clockwise). On the basis of this discussion, the average water velocity components (*U*
_{1},*V*
_{1}) in the boundary layer are obtained by the following expression which effectively subtracts the geostrophic current from the surface current, rotates the resultant by 13.9° in the anti-clockwise direction, multiplies it by a factor of 0.628, and adds the geostrophic velocity components.

and

Since the current velocity varies spatially as well as temporally, the current data near the iceberg are obtained by interpolating the current values at corresponding depths at locations A and B. A cubic interpolation function is used to give more weight to the current meter closer to the iceberg.

### Wind data

The surface wind velocity and direction were obtained from the log books of C S S *Dawson* which was in the area gathering other oceanographical and meteorological data.

## Results and Discussion

As mentioned previously, the mass and cross-section areas of the iceberg under study are not known, and the values of these parameters are assumed to be those given in Table 1. Figures 2 and 3 show the predicted trajectories of small, medium, and large non-tabular icebergs along with the observed trajectory of iceberg 20C during the storm. The predicted trajectory of a medium non-tabular iceberg (Fig. 2) is close to the observed trajectory, and thus iceberg 20C is assumed to be a medium non-tabular iceberg, an assumption also supported by field observations (personal communication from R.T. Dempster).

Fig. 2 Predicted and observed trajectories of iceberg 20C during the storm which passed over the Labrador Sea on 21–22 August 1972.

Fig. 3 Effect of iceberg size on the predicted trajectory.

Figure 4 shows two predicted trajectories when the iceberg is driven by the wind alone and by the currents alone. Figure 5 shows two predicted trajectories when the iceberg is driven by geostrophic currents only and by geostrophic currents and wind together. In these two figures, the effect of excluding a particular environmental force can be seen as the observed trajectory is also plotted there. It is evident from these results that currents and wind have a significant, effect on iceberg drift.

Fig. 4 Effect of wind forces on iceberg drift.

Fig. 5 Iceberg drift in geostrophic current with inclusion and exclusion of wind.

Figure 6 shows the effect of including and excluding the Coriolis force on the predicted iceberg drift trajectory. There is a considerable difference in the trajectories, and hence inclusion of Coriolis force is essential for a good prediction of trajectories. Figure 7 depicts the effect of including and excluding the water-acceleration term in the drift model. Though the net effect of excluding this term is not large, a better correlation is obtained in this ease by including it in the model.

Fig. 6 Effect of Coriolis force on iceberg drift.

Fig. 7 Predicted iceberg drift with inclusion and exclusion of water acceleration.

Finally, Figure 8 shows the difference between the trajectories when current data at 13 m depth are used for the average current velocity in the boundary layer as opposed to using Equations 11 and 12. As is evident from the results, the current data at 13 m depth as input to the model do not give a good prediction of the iceberg drift, and these data along with the data at greater depths have to be used to estimate the wind-generated currents in the boundary layer.

Fig. 8 Effect of considering the distribution of the wind-generated current velocity.

## Conclusions

A numerical model for the iceberg drift is presented, and the model is used to predict the drift trajectory of an iceberg using current and wind data from the field. The trajectory of only one iceberg is predicted because this particular iceberg (20C) drifted close to the installed current meters. The predicted trajectory compares well with the observed trajectory when the iceberg parameters are assumed to be those of a medium non-tabular iceberg (Table I). The trajectories of other icebergs are disturbed by the storm in a similar manner as shown in Figure 1, and the deviations from the regular drift trajectories start and stop at about the same time. The direct action of the wind is evident from the fact that the icebergs are moved back and forth as the wind direction changes due to the passage of the storm centre through the area. The indirect action of wind on the drift of icebergs is important because of the wind-generated currents and inertial oscillations in the boundary layer.

In this paper, we have shown that the physics of the iceberg drift model are known to the extent that a good prediction of an iceberg trajectory during a storm can be made provided detailed current and wind data are available as input to the model.

## Acknowledgements

The authors are grateful to Dean R.T. Dempster, Faculty of Engineering and Applied Science, Memorial University of Newfoundland, St. John’s, Newfoundland, for making available the iceberg data and to Dr Miles McPhee of the U.S. Army Cold Regions Research and Engineering Laboratory, Hanover, N.H., for his discussion. The help of many others is also gratefully acknowledged in gathering all the data which led to successful completion of this work.

## References

Allen, J H and others 1973 Iceberg study, Saglek, Labrador, including cruise report C.S.S. *Dawson*, August 7 to August 26, 1972. Memorial University of Newfoundland. Faculty of Engineering and Applied Science. Report

Bayly, I M
1971
Contribution on the inclusion of certain terms in the equations used in the simulation for the prediction of ice movement by Dr. O. Cochkanoff et al. In Proceedings of the Canadian seminar on icebergs, Halifax, Nova Scotia, 1971. Ottawa, Dept of National Defence: 153–155

Cochkanoff, O, Graham, J W, Warner, J L
1971
Simulation techniques in the prediction of iceberg motion. In Proceedings of the Canadian seminar on icebergs, Halifax, Nova Scotia, 1971. Ottawa, Dept of National Defence: 135–152

El-Tahan, M
Unpublished Modelling of iceberg drift. M.Eng. thesis, Memorial University of Newfoundland, Faculty of Engineering and Applied Science, 1980

Holden, B J
Unpublished Some observations on the Labrador Current at Saglek, Labrador. M.Eng. thesis, Memorial University of Newfoundland, Faculty of Engineering and Applied Science, 1974

Mountain, D
1980
On predicting iceberg drift. Cold Regions Science and Technology
1(3–4) : 273–282

Murray, J E
1969
The drift, deterioration and distribution of icebergs in the North Atlantic Ocean. In Ice seminar … sponsored by the Petroleum Society of CIM, Calgary, Alberta, 1968. Montreal, Canadian Institute of Mining and Metallurgy: 3–18 (Special Volume 10)

Napoleoni, J G P
Unpublished The dynamics of iceberg drift. M.Sc. thesis, University of British Columbia, Department of Geophysics and Astronomy, 1979

Neumann, G
1968
Ocean currents. Amsterdam, Elsevier