Ekman Transport

Ekman transport is the net vertically integrated mass or volume transport of fluid within a rotating, frictionally coupled boundary layer (the Ekman layer, spanning the upper ~20–150 m of the ocean) directed at a 90° angle to the surface wind stress vector—to the right of the wind in the Northern Hemisphere and to the left of the wind in the Southern Hemisphere. Formulated mathematically by the Swedish oceanographer Vagn Walfrid Ekman in 1902 and 1905, the phenomenon arises from the steady mechanical balance between downward momentum flux via vertical turbulent eddy viscosity and the planetary Coriolis effect. While the horizontal velocity vector at any discrete depth z rotates and decays downward in a three-dimensional Ekman spiral, integrating across the full boundary layer eliminates the turbulent viscosity parameter and locks the net transport perpendicular to the applied wind stress (ME = τ × k̂ / f). Across the global ocean, spatial gradients in Ekman transport govern coastal and equatorial upwelling, primary biological productivity by phytoplankton, and the wind-driven subtropical gyre circulation in physical oceanography.

During the 1893–1896 Norwegian Fram Polar Expedition, oceanographer Fridtjof Nansen froze his wooden schooner Fram into the Arctic pack ice to drift across the polar basin. Comparing daily astronomical positions against wind records, Nansen observed that Arctic sea ice never drifted downwind, but moved 20° to 40° to the right of the surface wind.

Nansen deduced that Earth's rotation deflected the drifting ice and that the ice dragged the underlying water, which was deflected further to the right with depth. Nansen posed the problem to Vilhelm Bjerknes at the University of Stockholm, whose student Vagn Walfrid Ekman solved it in 1902 and 1905 (On the Influence of the Earth's Rotation on Ocean-Currents). Replacing molecular kinematic viscosity (ν ≈ 10−6 m2/s) with turbulent vertical eddy viscosity (Az ≈ 10−3 to 10−1 m2/s), Ekman derived both the vertical velocity spiral and the 90° depth-integrated transport law.

---

Governing Momentum Balance in a Rotating Boundary Layer

Consider a steady (∂/∂t = 0), horizontally homogeneous (∂/∂x = ∂/∂y = 0), unstratified ocean of constant density ρ on an f-plane where horizontal pressure gradients vanish (∇hp = 0). Wind blowing across the sea surface exerts a horizontal wind stress vector τ = (τx, τy) (units of N/m2 or Pa, parameterized by τ = ρair CD |U10| U10, where ρair ≈ 1.22 kg/m3, CD ≈ 1.3 × 10−3, and U10 is 10 m wind speed). Within the upper oceanic boundary layer, the horizontal momentum equations reduce to a balance between the Coriolis force and the vertical divergence of turbulent shear stress:

−f v = (1 / ρ) (∂τxz / ∂z) = Az (∂2u / ∂z2)

+f u = (1 / ρ) (∂τyz / ∂z) = Az (∂2v / ∂z2)

where u(z) and v(z) are the eastward and northward current velocities at depth z (z ≤ 0), f = 2 Ω sin φ is the Coriolis parameter (Ω = 7.2921 × 10−5 rad/s and φ is latitude; f > 0 in the Northern Hemisphere and f < 0 in the Southern Hemisphere, with f ≈ 1.0 × 10−4 s−1 at 43° N), and Az is the kinematic vertical eddy viscosity (m2/s).

Analytical Solution: The Ekman Spiral and Ekman Depth

Aligning the coordinate system so that the surface wind stress blows northward along the +y axis (τx = 0, τy = τ) and applying the surface shear boundary condition ρ Az (∂u/∂z, ∂v/∂z)|z=0 = (0, τ) along with (u, v) → (0, 0) as z → −∞ yields the Ekman spiral solution (for f > 0, using positive depth d = −z ≥ 0 below the surface):

u(d) = V0 e−π d / DE cos(π/4 − π d / DE)

v(d) = V0 e−π d / DE sin(π/4 − π d / DE)

The total horizontal speed decays exponentially with depth d according to V(d) = √(u2 + v2) = V0 e−π d / DE, governed by two boundary-layer parameters:

1. Surface Ekman Speed (V0): At the ocean surface (d = 0), the current speed is V0 = τ / [ρ √(|f| Az)] (~2–3% of U10), and its direction is deflected 45° to the right of the wind vector in the Northern Hemisphere (45° to the left in the Southern Hemisphere).
2. Ekman Layer Depth (DE): Defined as the depth at which the current vector has rotated 180° (π radians) opposite to the surface current V0:

DE = π √(2 Az / |f|) At depth d = DE, the current speed has decayed to e−π ≈ 0.043 (4.3% of V0). For |f| = 10−4 s−1 and Az = 0.01–0.05 m2/s, the Ekman depth DE spans ~45 m to 100 m.

---

Integrating the horizontal momentum equations vertically from the base of the Ekman layer (z = −DE, where shear stress vanishes) to the sea surface (z = 0, where shear stress equals the surface wind stress τ) eliminates Az:

MEx = ∫−∞0 ρ u(z) dz = +τy / f and MEy = ∫−∞0 ρ v(z) dz = −τx / f

Or in vector notation (with unit vertical vector k̂):

ME = (τ × k̂) / f (kg · m−1 · s−1) and UE = ME / ρ = (τ × k̂) / (ρ f) (m2/s)

Because this vertical integral depends only on surface stress τ and planetary rotation f, the 90° net Ekman transport relation holds even when turbulent viscosity Az(z) varies arbitrarily with depth. For a 10 m/s wind (τ ≈ 0.16 N/m2) at 43° N (f = 10−4 s−1), net Ekman mass transport is ME = 1,600 kg/(m·s), or volume flux UE ≈ 1.56 m2/s.

---

Ekman Pumping, Sverdrup Gyre Dynamics, and Planetary Boundary Layers

By incompressible continuity (∂u/∂x + ∂v/∂y + ∂w/∂z = 0), horizontal divergence of Ekman volume transport sets the vertical Ekman pumping/suction velocity (wE) at the base of the Ekman layer:

wE = −∇h · UE = −(1 / ρ) k̂ · [∇ × (τ / f)]

In subtropical gyres (15°–45° latitude), negative wind stress curl (∇ × τ) drives Ekman pumping (wE < 0, ~30–50 m/yr downward), building a ~1.0 m sea-surface dome and compressing the underlying water column to drive equatorward interior Sverdrup transport (β Vg = (f / ρ) wE), returned poleward in western boundary currents (Gulf Stream, Kuroshio). Analogously, in the Atmospheric Planetary Boundary Layer (PBL, ~500–1,500 m thick), surface drag weakens the Coriolis force relative to the pressure gradient (−(1/ρ)∇p), turning surface winds 15°–35° across isobars toward low pressure.

---

Why does net Ekman transport occur at 90 degrees to the wind if the surface current moves at 45 degrees?

At the sea surface (z = 0), the water parcel is pulled forward by wind stress from above, deflected to the right (in the Northern Hemisphere) by the Coriolis force, and dragged backward and left by the slower water layer beneath it, producing a 45° surface deflection. Each deeper sub-layer is driven by the viscous shear of the layer above it, rotating further clockwise while decaying in speed as e−πd/DE. When velocity vectors from z = 0 to z = −DE are integrated vertically, downwind components cancel against deep reverse-angled components, leaving net mass transport at 90° to the surface wind stress.

What happens to Ekman transport at the Equator where the Coriolis parameter f goes to zero?

Because the steady Ekman transport equations MEx = τy / f and MEy = −τx / f contain f = 2 Ω sin φ in the denominator, the classical steady Ekman balance breaks down within ~1°–2° of the Equator (φ = 0°, where f = 0). Directly on the Equator, easterly trade winds push surface water downwind (westward) without Coriolis deflection, piling water up in the western Pacific warm pool and creating an eastward pressure gradient (∂p/∂x > 0) that drives the subsurface Equatorial Undercurrent (EUC). Simultaneously, just north (+1° to +3° N, f > 0) and just south (−1° to −3° S, f < 0) of the Equator, the same westward trade winds drive northward and southward Ekman transport, creating equatorial divergence and upwelling along 0° latitude.

Why is the classical 45-degree Ekman spiral rarely observed in its textbook form in the upper ocean?

Ekman's 1905 analytical spiral assumes a constant vertical eddy viscosity Az, a steady wind field, and an unstratified ocean. In the upper ocean, Az(z) varies with depth due to diurnal stratification, wave breaking, and Langmuir circulation (helical vortices aligned with the wind), which mix momentum across the upper layer and compress the surface deflection angle from 45° down to 10°–30°. Nevertheless, time-averaged ADCP and drifter arrays confirm that once diurnal fluctuations are averaged out, the depth-integrated net transport remains 90° to the wind stress (ME = τ × k̂ / f), because the vertical integral does not depend on Az(z).

---

1. Ekman, V. W. (1905). "On the influence of the Earth's rotation on ocean-currents." Arkiv för Matematik, Astronomi och Fysik, 2(11), 1–52. https://doi.org/10.1038/173590a0
2. Price, J. F., Weller, R. A., & Schudlich, R. R. (1987). "Wind-driven ocean currents and Ekman transport." Science, 238(4833), 1534–1538. https://doi.org/10.1126/science.238.4833.1534
3. Vallis, G. K. (2017). Atmospheric and Oceanic Fluid Dynamics: Fundamentals and Large-Scale Circulation (2nd ed.). Cambridge University Press. https://doi.org/10.1017/9781107588417
4. Chereskin, T. K. (1995). "Direct evidence for an Ekman balance in the California Current." Journal of Geophysical Research: Oceans, 100(C9), 18261–18269. https://doi.org/10.1029/95JC02182

---