In a recent Journal of Glaciology article, “Ice and ocean tilt measurements in the Beaufort Sea”, Reference Weber and ErdelyiWeber and Erdelyi {1976) postulate that wind-induced tilt could cause ice sheets to break at right angles to the direction of drift. They consider the floes to be rigid and to be acted on by a horizontal-surface wind force and by an opposing bottom-surface drag force; the resulting moment tilts the floe until it is equilibrated by buoyancy forces of the water. Then for a given angle of tilt, the maximum bending moment depends on the cube of the floe length; hence, the maximum floe length can be found if the tensile strength of the ice is known. For 3 m ice and a tilt of 30 µrad, Weber and Erdelyi find maximum floe lengths to be 328 m for a strength of 10⁵ N m^-2 and 890 m for a strength of 2 × 10⁶ N m^-2 Footnote *.

The assumption that the ice is rigid, however, is unreasonable. A first approximation to its actual behavior is that of a floating elastic plate, and there is overwhelming evidence to support such an assumed behavior. As will now be shown, the incorporation of elastic flexural behavior into the analysis of the ice considerably changes the conclusions to be drawn.

Referring to the infinitesimal element shown in Figure 1, the differential equations of equilibrium are

where p is the density of the ice, M is the bending moment, q is the applied vertical load, v is the displacement of the ice sheet normal to its surface, V is the shearing force, T_a is the wind stress, T_W is the water stress, and a is the thickness of the ice. In the above, motion in the x-direction is assumed to be uniform and T_a and T_w are assumed to be constant; they must therefore be equal unless some horizontal force exists due to action of an adjacent floe.

Using the plate-bending constitutive law

where

is the flexural stiffness (Eis Young’s modulus, v is Poisson’s ratio), the equation of equilibrium in terms of displacement becomes

Fig. 1.

Equilibrium element.

For q = o, the solution to Equation (5) is

Where

(1/λ is the characteristic plate length).

Consider first a long ice sheet, floating free at x = o. Regularity conditions require that

while boundary conditions at x = o require the moment and shear there to be zero, i.e.

The corresponding solution for the displacement is

The maximum bending moment occurs at x = λπ/4 and is given by

The maximum bending stress σ_max is obtained directly from the maximum bending moment and is given by

Of importance here are the wind stress values required to cause fracture. Using a tensile strength value of 10⁵ Nm^-2, E/(l −,v ²) = 10⁹ N m^-2, and p = l Mg m^-3 and setting τ_a = τ_w one obtains, from Equations (ll) and (12), a wind stress to cause fracture of 5040 N m^-2 for 3 m ice and 2150 N m^-2 for 10 cm ice. These values are so much larger than those occurring (on the order of 0.5 N m^-2 that fracture of ice sheets due to wind loading must be ruled out. Although the above calculations were made for long (semi-infinite) plates, the results for plates of finite length do not differ significantly. The response of the ice sheet to the surface load given by Equation (10) amounts to an edge disturbance which rapidly dies out (as do flexural stresses); subsequently, the wind stress is equilibrated by the shearing force in the plate. The addition of an axial force at one end of the ice sheet resulting from the action of a contiguous floe has virtually no effect on the above solution; wind reaches of several hundred kilometers are required before interaction of axial and flexural effects becomes significant.

Finally, it is of interest to compare typical tilt values with those measured by Reference Weber and ErdelyiWeber and Erdelyi (1976). According to the elastic theory used here, tilt will be significant only near the end of the ice sheet. The maximum slope occurs at the free end and is given by

For a wind stress of 0.5 N m^-2 and other data as above, the till is 3.16.µrad for 3 m ice and 1 7.3 µrad for 10 cm ice. If tilt is interpreted as the difference in level between two points (as measured by Weber and Erdelyi) and the points are chosen to be x = o and x = π/2lambda;, then the tilt is 2 µrad for 3 m ice and 11 µrad for 10 cm ice. These values are of the same order of magnitude as those measured by Weber and Erdelyi using points 120 m apart. It would be interesting if tilt measurements could somehow be used to check Equation (10).