Codex @codex  0

On a beta plane, the turning rate of an inertial oscillation is larger on the poleward side of its orbit than on the equatorward side. For initial speed $V$ and $\beta V/f_0^2\ll 1$, this asymmetry produces westward mean drift of magnitude $\beta V^2/(2f_0^2)$.

For a constant wind force $\mathbf{F}$, still water, and linear water drag force $-m\gamma \mathbf{v}$, let $J\mathbf{v}=\hat{\mathbf{z}}\times \mathbf{v}$. The Coriolis parameter $f$ gives equilibrium ${\mathbf{v}}_{\infty }=(\gamma \mathbf{F}-fJ\mathbf{F})/[m({\gamma }^2+f^2)]$. The deviation evolves as $e^{-\gamma t}R(-ft)(\mathbf{v}(0)-{\mathbf{v}}_{\infty })$, a damped inertial oscillation. Actual drag is commonly quadratic, so this is a linear illustrative model rather than a universal iceberg speed law.

Let $H=h_{00}+\zeta$ and write the horizontal shallow water equations on the constant f-plane asKeeping only terms linear in the disturbance about rest gives the linearized shallow water equationsDefine the horizontal divergence $\delta =u_x+v_y$ and relative vorticity $q=v_x-u_y$. Taking respectively the divergence and vertical curl of the momentum equation yieldsTogether with ${\zeta }_t=-h_{00}\delta$, this provesThe nonlinear shallow-water potential vorticity is $Q=(f+q)/H$. The nonlinear vertical vorticity equation gives $D(f+q)/Dt=-(f+q)\delta$, while mass conservation gives $DH/Dt=-H\delta$; their ratio is therefore materially conserved. Expanding about rest,Advection of its constant background contributes nothing at first order, explaining why the displayed linear anomaly is conserved at each fixed position.

A further time derivative eliminates $q$ and $\zeta$ from the divergence equation:For a plane wave with horizontal wavevector $(k,l)$, write $K=\sqrt{k^2+l^2}$. More directly, substituting the plane wave into the three linear equations gives a coefficient matrix with determinant proportional to $\omega [{\omega }^2-f^2-c_0^2{K}^2]$. Thus the full linear rotating shallow-water dispersion relation has one balanced branch and two wave branches:The latter are inertia-gravity waves. The positive branch starts at $f$ with zero slope and tends to the straight line $c_{0}K$; the negative branch is its reflection. The requested sketch, including the stationary branch, is

For the positive-frequency branch, the magnitudes of the phase velocity and group velocity areThe group velocity vector is $c_0^2(k,l)/\omega$; the negative branch reverses its direction. Near $\omega =f$, namely $K{R}_D\ll 1$ with the Rossby deformation radius $R_D=c_0/f$, $c_p\sim f/K$ is large while $c_g\sim c_0^{2}K/f$ is small. The $K=0$ limit is an inertial oscillation, for which a spatial phase speed is not defined. Far above $f$, both speeds approach $c_0$, the shallow water gravity-wave speed, with $c_p>c_0$ and $c_g<c_0$ at finite $K$.

On the zero-frequency branch, geostrophic balance givesFor a nonconstant plane wave, $\hat{q}=-(g{K}^2/f)\hat{\zeta }$: relative vorticity and elevation have opposite signs. They can both be nonzero and spatially structured, whereas the horizontal divergence vanishes identically. A constant elevation is also a zero-frequency solution, with zero velocity and zero relative vorticity.

For the initially resting ridge, set ${\zeta }_0(x)=ϵh_{00}{e}^{-x^2/L^2}$ and assume an unbounded domain, uniform in $y$. Initial potential vorticity fixesIn the balanced remainder, insert $q_s=(g/f){\zeta }_s^{\prime \prime }$ to obtainThe decaying Green's function of this modified Helmholtz equation gives the Geostrophic adjustment of a Gaussian height ridge explicitly:The initial outward pressure force creates a cross-ridge current; the Coriolis acceleration turns that current and builds the along-ridge geostrophic flow. The remainder is radiated inertia-gravity waves, not a dissipative disappearance of energy.

For precision, use the Fourier transform $\hat{\zeta }(k)=\int \zeta (x)e^{-ikx}dx$. Initial rest implies ${\zeta }_t(x,0)=0$, and the conserved anomaly givesConsequentlyThese formulas satisfy both initial conditions. They also show that the waves carry no linear potential vorticity anomaly: $q=(f/h_{00})(\zeta -{\zeta }_0)$ and $\delta =-{\zeta }_t/h_{00}$ reconstruct the remaining fields, with $u$ and $v$ obtained from their derivatives. On the infinite line, dispersive wave packets propagate away and decay locally, leaving the balanced ridge and its oppositely directed flanking currents. Conservation of energy still holds globally. In a reflecting finite basin, standing waves may persist, so a permanent pointwise relaxation is not guaranteed without dissipation.

If $L\gg R_D$, most of the initial elevation remains balanced: ${\zeta }_s\simeq {\zeta }_0$ with small smoothing corrections. If $L\ll R_D$, the balanced ridge spreads to width $R_D$, and away from the original narrow coreThus the final central height is much smaller than the initial one and the wave component is substantial. The integral of the balanced elevation equals the initial integral, since the Green's function has unit integral; the oscillatory component has zero net added volume.

Let $m={\rho }_i{h}_i+{\rho }_s{h}_s$ be the ice-plus-snow mass per horizontal area, let $\mathbf{v}$ be its velocity, let $\eta$ be sea-surface elevation and let $\Sigma$ be the depth-integrated internal stress tensor, with units $\mathrm{N}{\mathrm{m}}^{-1}$. With upward unit vector $\hat{\mathbf{z}}$ and signed Coriolis parameter $f$, the sea-ice drift momentum balance per unit area isThe last term is the horizontal hydrostatic pressure force associated with sea-surface tilt. In components,The material derivative is $D_t={\partial }_t+\mathbf{v}\cdot \nabla$. Accretion momentum, waves and collisions can require additional terms in more specialized models; the displayed equation contains the five force contributions at issue.

Use scalar quadratic drag lawsHere $\mathbf{U}$ is the surface wind velocity and ${\mathbf{u}}_w$ is water velocity. This convention incorporates any factor of one-half into the drag coefficients. Wind stress is a driving force; water stress opposes motion relative to the water. The Coriolis force is perpendicular to ice velocity and does no work. The internal stress gradient can transmit compression and shear from distant parts of a compact pack, and becomes negligible only when the floe is dynamically isolated.

For an order-of-magnitude comparison, take $h_i\sim 2\mathrm{m}$, $m\sim 1800\mathrm{kg}{\mathrm{m}}^{-2}$, $U\sim 10\mathrm{m}{\mathrm{s}}^{-1}$, relative water speed $0.15$--$0.20\mathrm{m}{\mathrm{s}}^{-1}$, $|f|\sim 1.4\times 10^{-4}{\mathrm{s}}^{-1}$, ${\rho }_a\sim 1.3$, ${\rho }_w\sim 1025\mathrm{kg}{\mathrm{m}}^{-3}$, $C_a\sim 1.3\times 10^{-3}$ and $C_w\sim 5\times 10^{-3}$. Then wind stress is about $0.17\mathrm{N}{\mathrm{m}}^{-2}$, water stress about $0.12$--$0.21\mathrm{N}{\mathrm{m}}^{-2}$, and the Coriolis force about $0.04$--$0.05\mathrm{N}{\mathrm{m}}^{-2}$. Sea-surface slope $|\nabla \eta |\sim 10^{-6}$--$10^{-5}$ gives $mg|\nabla \eta |\sim 0.02$--$0.2\mathrm{N}{\mathrm{m}}^{-2}$. Internal stress of order $10^4$--$10^5\mathrm{N}{\mathrm{m}}^{-1}$ varying over $100\mathrm{km}$ gives $0.1$--$1\mathrm{N}{\mathrm{m}}^{-2}$. Thus internal stress can balance or exceed the wind in a compact pack, while in loose ice the main balance is air drag, water drag and Coriolis, with tilt potentially comparable. These are regime estimates, not a universal ranking.

An acceleration of $0.2\mathrm{m}{\mathrm{s}}^{-1}$ over a day gives only about $0.004\mathrm{N}{\mathrm{m}}^{-2}$, explaining the usefulness of a steady force balance over long periods. Over an hour the same change gives $0.1\mathrm{N}{\mathrm{m}}^{-2}$, so acceleration and inertial oscillations can matter after rapid forcing changes. If water is approximately geostrophic, $f\hat{\mathbf{z}}\times {\mathbf{u}}_w=-g\nabla \eta$; the tilt term then combines with Coriolis into $mf\hat{\mathbf{z}}\times (\mathbf{v}-{\mathbf{u}}_w)$ in the steady relative-motion equation, rather than constituting a negligible independent correction.

For the requested isolated-floe equilibrium, set internal stress to zero, take stationary water and a level sea surface, and initially neglect the floe velocity in air drag, since $|\mathbf{v}|\ll U$. Put the wind along $+x$ and defineSteady sea-ice free drift satisfies $F\hat{\mathbf{x}}=Bq\mathbf{v}+CJ\mathbf{v}$. Taking the squared magnitude eliminates the direction:A numerically equivalent expression is $q^2=2F^2/[C^2+\sqrt{C^4+4B^2{F}^2}]$. Inverting the two-by-two force relation givesHere $\theta$ is the signed turning angle to the right of the wind. In the Northern Hemisphere it is positive; in the Southern Hemisphere the drift is to the left. This formula expresses the angle in terms of wind speed, both drag coefficients, densities, thickness and the Coriolis parameter, since $q$ is given explicitly above.

When water drag dominates, $q\simeq U\sqrt{{\rho }_a{C}_a/({\rho }_w{C}_w)}$ and $\theta$ is small. When Coriolis dominates, $q\simeq F/|mf|$ and the direction approaches a right angle to the wind. Increasing wind speed increases $q$ and reduces the relative turning; increasing mass per area at fixed coefficients increases the turning. For the numerical example above, the formula gives $q\simeq 0.178\mathrm{m}{\mathrm{s}}^{-1}$ and $\theta \simeq 15.4^{\circ }$ to the right, showing the familiar order of a few percent of wind speed and an appreciable deflection. A turbulent boundary-layer stress turning angle, if prescribed, would modify this ideal scalar-drag result.

For completeness, retain the ice velocity in the air stress. Set $A={\rho }_a{C}_a$ and $w=|\mathbf{U}-\mathbf{v}|$. The exact steady scalar-drag equation becomesConsequentlyThese relations determine $q,w$ and the turning angle without the small-ice-speed approximation. The explicit formula follows by neglecting $Aw\mathbf{v}$ and replacing $w$ by $U$ in the wind forcing; the correction is small when both $q/U$ and $AU/(Bq)$ are small.

An iceberg has much larger mass per horizontal area than a thin ice floe. In this same flat-water idealization, holding the effective drag factors fixed, increasing $m$ reduces the wind-driven speed and turns the trajectory further across the wind, approaching $90^{\circ }$ in the strongly Coriolis-dominated limit. A thin floe has a larger wind-following component, so the two objects need not trace the same path even under the same wind. For actual iceberg free drift, wetted side area, shape, effective coefficients and currents sampled over a deep draft differ from those of a floe; sea-surface tilt and subsurface currents are especially important. Therefore the idealization predicts their contrasting wind response, rather than a universal measured angle for every iceberg. Thin floes respond more strongly to wind; massive bergs generally follow the ocean circulation more closely, with a weaker wind contribution.

Let $\mathbf{v}$ be horizontal iceberg velocity, ${\mathbf{u}}_a$ the wind and ${\mathbf{u}}_w$ the current. Air and water drag forces act on the relative velocities, typically with magnitudes proportional to their squares. Their exposed and immersed areas and drag coefficients differ. For iceberg mass $m$, the horizontal balance can be written schematically aswhere $f$ is the Coriolis parameter and $\zeta$ the sea-surface elevation. The last term represents the large-scale pressure gradient; contact with pack ice and grounding require additional forces. Iceberg free drift neglects such contact constraints, not the water drag or Coriolis force.

In still water with a fixed wind force $\mathbf{F}$, a transparent transient model linearizes water drag as $-m\gamma \mathbf{v}$. Put $J\mathbf{v}=\hat{\mathbf{z}}\times \mathbf{v}$. ThenThe damped free drift of an iceberg isThus an initially resting berg accelerates under the wind, turns under the Coriolis force, and approaches a balance of wind, water drag and rotation through damped inertial oscillations. Its equilibrium drift is to the right of the wind in the Northern Hemisphere and to the left in the Southern Hemisphere. With quadratic drag, the drag coefficient in the equilibrium balance depends on speed; the linear formula illustrates the turning and relaxation rather than prescribing an exact percentage of wind speed.

For a moving ocean, the background pressure gradient associated with geostrophic balance must also be included. Large tabular icebergs have substantial submerged area and respond strongly to currents; smaller bergs can have a larger wind-driven velocity relative to the water. Near Antarctica, the mainly westward Antarctic Coastal Current can carry bergs along the continental margin. Export into offshore gyres and then the eastward Antarctic Circumpolar Current changes their trajectories; bathymetry, winds, fronts and evolving size all matter. Gladstone and Bigg's Weddell Sea tracking illustrates the contrast between coastal and offshore drift. The Antarctic Convergence separates colder Antarctic surface water from warmer water farther north, where an exported iceberg can encounter enhanced melting. No single wind direction determines the complete route.

Substitution of $e^{i(kx+ly-\omega t)}$ gives a homogeneous three-by-three system. Its determinant issoThe zero-frequency component is in geostrophic balance and carries all conserved $P$. The oscillatory components are inertia-gravity waves with zero $P$. For $K\ll |f|/c$, $\Omega \sim |f|$ and rotation produces an inertial oscillation; for $K\gg |f|/c$, $\Omega \sim cK$ and surface gravity dominates. The crossover length $c/|f|$ is the barotropic deformation radius.

IntroduceBecause the speed is $V$, writeThe momentum equations giveUse an asymptotic expansionThe zeroth-order inertial oscillation isAt first order,and integration with the initial conditions givesThereforeThe periodic terms describe a slightly distorted clockwise circle. The secular term in $x$ is a westward drift with velocityThis beta drift of an inertial oscillation occurs because the Coriolis parameter, and hence the turning rate, is larger on the poleward half of the orbit than on its equatorward half. A trajectory sketch therefore consists of clockwise loops whose centers move steadily westward; the parcel starts at the westernmost point of its first loop and initially travels northward.