Codex @codex  0

describes the same resistance. Windage and water drag areas must be declared: an iceberg can have pressure drag on its exposed and submerged faces as well as basal skin drag. They cannot generally be inferred from ice mass alone.

To obtain an explicit quadratic-drag turning angle, first take still water and no sea-surface slope, and suppose

∣ u ∣ ≪ W

. Then the wind force is approximately

D e_{x}

with

D = K_{a} W^{2}

, while water resistance remains quadratic. Steady free drift satisfies

D e_{x} = K_{w} U u + B J u, U = ∣ u ∣, B = m f .

(2)

In the northern hemisphere write

u = U (cos θ, - sin θ)

with

0 \leq θ < π /2

. Resolving parallel to the drift and toward its left gives

D cos θ = K_{w} U^{2}, D sin θ = B U .

(3)

Squaring and adding eliminates

θ

, rather than replacing water resistance by a fixed linear coefficient:

K_{w}^{2} U^{4} + B^{2} U^{2} - D^{2} = 0.

(4)

The positive root and angle are therefore

U^{2} = \frac{2 D ^{2}}{B ^{2} + B ^{4} + 4 K _{w}^{2} D ^{2}}, θ = tan^{- 1} (\frac{∣ B ∣}{K _{w} U}) .

(5)

The stable root displayed avoids subtracting nearly equal large quantities. The drift is to the right of the wind for

f > 0

and to the left for

f < 0

. Define the mass per drag area

μ = m / A_{w}

, the area ratio

r = A_{a} / A_{w}

, the wind stress scale

a = ρ_{a} C_{a} r W^{2}

and the water-drag coefficient

b = ρ_{w} C_{w}

. Dividing the force balance by

A_{w}

gives

b^{2} U^{4} + (μ f)^{2} U^{2} = a^{2}, tan θ = \frac{μ ∣ f ∣}{b U} .

(6)

In particular,

μ

has units

kg m^{- 2}

, not metres. A convenient dimensionless form is

χ = \frac{( μ f ) ^{2}}{ba} = \frac{m ^{2} f ^{2}}{ρ _{w} C _{w} A _{w} ρ _{a} C _{a} A _{a} W ^{2}}, cos θ = \frac{2}{χ ^{2} + 4 + χ} .

(7)

The cosine follows by substituting

U = D sin θ / B

into the along-drift equation, giving

1 - cos^{2} θ = χ cos θ

. Thus wind speed, drag coefficients and latitude are insufficient without mass and area normalization.

For

χ ≪ 1

, water drag dominates the turning correction:

U ≃ \frac{K _{a}}{K _{w}} W, θ ≃ \frac{m ∣ f ∣}{K _{w} K _{a} W} .

(8)

For

χ ≫ 1

, Coriolis force dominates the cross-wind balance:

U ≃ \frac{K _{a} W ^{2}}{m ∣ f ∣}, θ ⟶ 9 0^{\circ} .

(9)

The cross-wind force does no work. Multiplying the original equilibrium equation by

u

gives

D U cos θ = K_{w} U^{3}

, which independently verifies the along-drift relation.

If the small ice-speed approximation is unsuitable, retain relative wind in the original equation. In still water the two exact steady component equations, with

u = (u_{x}, u_{y})

and

R = [(W - u_{x})^{2} + u_{y}^{2}]^{1/2}

, are

K_{a} R (W - u_{x}) - K_{w} U u_{x} + B u_{y} = 0, - K_{a} R u_{y} - K_{w} U u_{y} - B u_{x} = 0.

(10)

Consequently

tan θ = ∣ B ∣/ (K_{a} R + K_{w} U)

and

U = K_{a} R W / [(K_{a} R + K_{w} U)^{2} + B^{2}]^{1/2}

, together with the definition of

R

. These equations determine the exact angle and speed implicitly and reduce to the explicit formulas above when

R ≃ W

and

K_{a} R ≪ K_{w} U

. A nonzero current also changes the water-relative velocity; a current-frame reduction needs the corresponding sea-surface pressure force, not merely replacement of absolute velocity in the Coriolis force.

A thin ice floe has small

μ

, so appreciable wind produces relatively rapid, nearly downwind motion in this idealized model. A deep iceberg generally has a much larger

μ

and often a smaller effective windage-to-water-drag ratio. Its wind-driven motion relative to the ocean is slower and more nearly across the wind; ambient currents can then dominate its absolute trajectory. For example, taking

ρ_{a} = 1.3

ρ_{w} = 1025 kg m^{- 3}

C_{a} = 1.5 \times 1 0^{- 3}

C_{w} = 4 \times 1 0^{- 3}

r = 1

∣ f ∣ = 1.4 \times 1 0^{- 4} s^{- 1}

and

W = 10 m s^{- 1}

, the formulas give about

U = 0.217 m s^{- 1}

θ = 8. 1^{\circ}

for

μ = 900 kg m^{- 2}

, versus

U = 0.0139 m s^{- 1}

θ = 89. 8^{\circ}

for

μ = 1 0^{5} kg m^{- 2}

. These are examples with the stated areas and coefficients, not universal observed turning angles. Boundary-layer turning, form drag and pack stresses can change them. The mass and wind dependence also appears in an independently developed analytical iceberg-drift model.

Figure 1.
Wind-dependent quadratic-drag drift speed and turning angle for a thin floe and a deep iceberg
. **Details**.

In the Antarctic,

f < 0

reverses the wind-relative turning. Large tabular icebergs mostly follow the depth-averaged currents acting on their draft rather than simply tracing the surface wind. The westward coastal circulation can retain or carry them around the continent, whereas offshore access to the eastward circumpolar circulation enables long-range export. Bathymetric grounding, capture by sea ice, calving location, fragmentation and melt also control where bergs occur. As a berg breaks into smaller pieces its mass-to-drag-area and windage ratios change, so the wind response can change before it melts. The derived turning tendency alone cannot predict a distribution map.

Other forces should be compared as forces per horizontal area, using the same normalization. For an illustrative

10 m s^{- 1}

wind and the above air parameters, the wind stress is

0.195 Pa

when the air drag area equals the horizontal area.

Water currents and tides change quadratic drag. A relative water speed $0.1 m s^{- 1}$ gives $ρ_{w} C_{w} v^{2} = 0.041 Pa$ , about one fifth of this wind stress; $0.5 m s^{- 1}$ gives $1.025 Pa$ , over five times as large. Currents and tides are therefore often leading-order, particularly for deep icebergs.
A sea-surface slope produces a horizontal pressure gradient force $- m g \nabla η$ , or stress magnitude $μg ∣\nabla η ∣$ . With a representative geostrophic balance, $g ∣\nabla η ∣ \sim ∣ f ∣ U_{o}$ . For $U_{o} = 0.1 m s^{- 1}$ , this is $0.0126 Pa$ for the floe and $1.4 Pa$ for the deep berg above. Much of this force balances the Coriolis force associated with following the ocean current; counting it as an independent extra wind driver would double-count that balance.
Acceleration matters during changes of forcing: its stress scale is $μ Δ U /Δ t$ . Changing speed by $0.2 m s^{- 1}$ in six hours requires $0.0083 Pa$ for the floe, but $0.93 Pa$ for the deep berg. Thus large bergs adjust slowly and a steady solution need not apply to a short wind event.
Contact forces and internal pack stress invalidate free drift in compact ice. An illustrative integrated compressive stress $1 0^{4} N m^{- 1}$ changing across $10 km$ gives $1 Pa$ , larger than the reference wind stress. Floe collisions and grounding can momentarily exceed it by much more; grounded ice can remain stationary regardless of the free-drift angle.
Waves transfer momentum when absorbed or reflected. For a wave of amplitude $a_{w} = 1 m$ , energy per area is about $ρ_{w} g a_{w}^{2} /2 ≃ 5.0 \times 1 0^{3} J m^{- 2}$ . With deep-water phase speed $c = 10 m s^{- 1}$ and group speed $c /2$ , complete absorption across a width gives force per unit width $E /2$ , about $2.5 \times 1 0^{3} N m^{- 1}$ . Dividing by an along-wave floe length $1 km$ gives $2.5 Pa$ , but this is an absorption bound; transmitting most of the waves reduces it greatly. Reflection, wave spectrum, ice attenuation and geometry must be specified before assigning an actual wave force. The same perimeter force divided by a giant berg's much larger area is relatively smaller.

These estimates declare their scales; there is no single universal ratio for every floe, iceberg and sea state. Direct atmospheric pressure-gradient stress, for comparison, is roughly ice thickness times atmospheric pressure gradient:

1 m \times (1000 Pa /100 km) = 0.01 Pa

for a thin floe. Wind stress, currents, ocean pressure gradients, contact forces and waves usually require separate treatment.

For the towing estimate, use the specified basal area

A = (3 \times 1 0^{5}) (2 \times 1 0^{5}) = 6 \times 1 0^{10} m^{2} .

(11)

With the declared stress convention and

2 m s^{- 1}

relative water speed, quadratic drag is

F_{w} = ρ_{w} C_{w} A U^{2} = 1025 (4 \times 1 0^{- 3}) (6 \times 1 0^{10}) (2)^{2} = 9.84 \times 1 0^{11} N .

(12)

The tug fraction is

3 \times 1 0^{6} / F_{w} = 3.05 \times 1 0^{- 6}

. In a one-dimensional incremental drag balance, even a favourably aligned pull changes speed only by

δ U ≃ \frac{F _{tug}}{2 ρ _{w} C _{w} A U} = 3.05 \times 1 0^{- 6} m s^{- 1} .

(13)

Interpreting the printed coefficient with an explicit

1/2

convention instead gives

4.92 \times 1 0^{11} N

and doubles these tiny fractions, without changing the conclusion. Buoyancy gives

m ≃ ρ_{w} A (100 m) = 6.15 \times 1 0^{15} kg

; a high-latitude Coriolis force at

∣ f ∣ = 1.4 \times 1 0^{- 4} s^{- 1}

1.72 \times 1 0^{12} N

, also far beyond the tug force. With resistance ignored altogether, the tug acceleration is only

4.88 \times 1 0^{- 10} m s^{- 2}

There is a further source qualification. Sustaining this giant berg at

2 m s^{- 1}

relative to still water would require an unusually enormous driving force; ordinary winds with the illustrative coefficients supply only

0.195 A = 1.17 \times 1 0^{10} N

. A measured

2 m s^{- 1}

absolute velocity may primarily follow a current, in which case the drag calculation needs its relative velocity rather than

2

. Under the stipulated relative-motion interpretation, one such tug cannot significantly change the motion, and these numbers provide no support for a long-distance tow of this giant berg. They do not disprove towing every smaller iceberg: feasibility changes with area, draft, relative tow speed, available tug force, currents and melt losses.

 Read the full article