Ekman Pumping

Ekman Pumping is the vertical fluid motion—comprising downward Ekman pumping (downwelling) and upward Ekman suction (upwelling)—generated at the boundary between a turbulent, friction-dominated planetary boundary layer and the nearly inviscid geostrophic interior of the ocean or atmosphere. First formulated quantitatively by the Swedish oceanographer Vagn Walfrid Ekman in 1905 and extended to basin-scale wind-driven circulation by Harald Sverdrup (1947) and Jule G. Charney and Arnt Eliassen (1949), the mechanism translates horizontal frictional stresses into vertical mass flux through the three-dimensional continuity equation. Because the Coriolis effect deflects depth-integrated Ekman transport by 90° to the right of the surface wind stress in the Northern Hemisphere (and 90° to the left in the Southern Hemisphere), spatial variations in the wind field—specifically a non-zero vertical component of wind-stress curl—drive horizontal convergence or divergence of water within the upper 20–100 meters of the ocean. Conservation of mass requires this converging or diverging boundary-layer fluid to exit or enter the boundary layer vertically at a characteristic velocity w_E = k · ∇ × (τ / ρ₀f), typically on the order of 10⁻⁶ to 10⁻⁵ m s⁻¹ (30 to 300 meters per year) across open-ocean gyres. Despite its modest magnitude, Ekman pumping is the primary mechanical driver of the large-scale ocean thermocline circulation, subtropical and subpolar gyres, marine nutrient upwelling, and atmospheric boundary-layer spin-down of synoptic cyclones and anticyclones.

Horizontal Momentum Balance in the Surface Ekman Layer

Consider a horizontally broad, rotating fluid on a midlatitude β-plane or f-plane with reference density ρ₀ and Coriolis parameter f = 2Ω sin φ, where Ω = 7.292 × 10⁻⁵ rad s⁻¹ is Earth's angular rotation rate and φ is latitude. Within the oceanic surface boundary layer (extending from the sea surface z = 0 downward to a depth z = −d_E below which turbulent shear stress vanishes), the steady, linear, ageostrophic horizontal momentum equations balance the Coriolis acceleration against the vertical divergence of the turbulent Reynolds shear stress τ(x, y, z) = (τ_x, τ_y):

−f v_E = (1 / ρ₀) ∂τ_x / ∂z , +f u_E = (1 / ρ₀) ∂τ_y / ∂z

Here (u_E, v_E) = (u − u_g, v − v_g) denotes the ageostrophic Ekman velocity departure from the interior geostrophic flow (u_g, v_g). Integrating these equations vertically across the entire Ekman layer from z = −d_E (where τ(−d_E) ≈ 0) to the free surface z = 0 (where the atmosphere applies the surface wind stress τ₀ = (τ_0x, τ_0y)) yields the depth-integrated Ekman volume transport U_E = (U_E, V_E) = ∫_−dE⁰ (u_E, v_E) dz:

U_E = + τ_0y / (ρ₀ f) , V_E = − τ_0x / (ρ₀ f) ⇒ U_E = (τ₀ × k) / (ρ₀ f)

where k is the upward unit vector. Remarkably, this integrated mass transport depends only on the surface wind stress τ₀ and the local Coriolis parameter f, and is independent of the specific vertical profile of eddy viscosity A_z(z) within the mixed layer.

Mass Continuity and the Ekman Pumping Velocity Formula

To determine the vertical velocity w_E = w(x, y, −d_E) at the base of the surface Ekman layer, one integrates the incompressible three-dimensional continuity equation ∂u/∂x + ∂v/∂y + ∂w/∂z = 0 vertically from z = −d_E to the rigid-lid sea surface z = 0 (where w(x, y, 0) = 0):

w(0) − w(−d_E) = − ∫_−dE⁰ ( ∂u_E / ∂x + ∂v_E / ∂y ) dz = − ( ∂U_E / ∂x + ∂V_E / ∂y ) = − ∇_h · U_E

Solving for w_E ≡ w(−d_E) and substituting the Ekman transport components (U_E, V_E) gives the general Ekman pumping equation on a spherical or β-plane Earth where f(y) varies with latitude (β = df/dy):

w_E = + ∇_h · U_E = ∂ / ∂x ( τ_0y / ρ₀f ) − ∂ / ∂y ( τ_0x / ρ₀f ) = k · ∇ × ( τ₀ / ρ₀f )

Expanding the quotient derivative using ∂f/∂x = 0 and ∂f/∂y = β separates the vertical velocity into two distinct physical contributions:

w_E = ( 1 / ρ₀f ) ( ∂τ_0y / ∂x − ∂τ_0x / ∂y ) + ( β τ_0x ) / ( ρ₀ f² ) = ( k · ∇ × τ₀ ) / ( ρ₀ f ) + ( β τ_0x ) / ( ρ₀ f² )

1. Wind-Stress Curl Term [(k** · ∇ × τ₀) / (ρ₀f)]:** Arises from horizontal shear or curvature in the surface wind stress. In the Northern Hemisphere (f > 0), anticyclonic wind stress (negative curl, ∇ × τ₀ · k < 0) causes horizontal convergence of Ekman transport (∇_h · U_E < 0) and downward Ekman pumping (w_E < 0). Cyclonic wind stress (positive curl, ∇ × τ₀ · k > 0) causes horizontal divergence of Ekman transport (∇_h · U_E > 0) and upward Ekman suction (w_E > 0).
2. Planetary β-Divergence Term [(β τ_0x) / (ρ₀ f²)]: Even under a spatially uniform zonal wind stress τ_0x with zero curl, the meridional variation of the Coriolis parameter f(y) alters the magnitude of the meridional Ekman transport V_E = −τ_0x/(ρ₀f) across latitude circles, generating a non-zero divergence ∂V_E/∂y = +β τ_0x / (ρ₀f²). At midlatitudes, the wind-stress curl term is typically 5 to 10 times larger than the β-term, whereas in the tropics (small f), both terms become dynamically comparable.

Subtropical Downwelling vs. Subpolar and Equatorial Upwelling

Across the world's major ocean basins, the meridional profile of atmospheric surface winds—easterly trade winds in the tropics (0°–30° latitude), prevailing westerlies at midlatitudes (30°–60° latitude), and polar easterlies at high latitudes (>60°)—establishes permanent belts of Ekman pumping and suction:

• Subtropical Gyres (15°–45° N/S — Ekman Downwelling): Between the tropical easterlies (τ_0x < 0) and midlatitude westerlies (τ_0x > 0), the meridional shear −∂τ_0x/∂y imparts a negative (anticyclonic) wind-stress curl in the Northern Hemisphere (and positive anticyclonic curl in the Southern Hemisphere). In both hemispheres, Ekman transport converges toward the center of the subtropical gyre (the Sargasso Sea in the North Atlantic and the North Pacific Subtropical Gyre), driving downward Ekman pumping of w_E ≈ −1 × 10⁻⁶ to −2 × 10⁻⁶ m s⁻¹ (−30 to −60 m yr⁻¹). This downwelling depresses the main thermocline to depths of 600–1,000 m, creating nutrient-depleted (oligotrophic) surface waters and storing warm water in the upper ocean.
• Subpolar Gyres (45°–65° N and the Southern Ocean — Ekman Suction): Between the midlatitude westerlies and polar easterlies (or across the northern flank of the Antarctic Circumpolar Current along the Antarctic Divergence at ~60°–65° S), cyclonic wind-stress curl drives horizontal Ekman divergence and upward Ekman suction of w_E ≈ +1 × 10⁻⁶ to +4 × 10⁻⁶ m s⁻¹ (+30 to +125 m yr⁻¹). This upward flux domes the pycnocline toward the surface and brings cold, carbon- and nutrient-rich deep water into the euphotic zone.
• Equatorial Upwelling Band (5° S–5° N): Because the Coriolis parameter f = 2Ω sin φ changes sign across the equator, easterly trade winds (τ_0x < 0) drive northward Ekman transport (V_E > 0) just north of the equator and southward Ekman transport (V_E < 0) just south of the equator. This sharp meridional divergence drives vigorous equatorial upwelling (w_E ≈ 10⁻⁵ m s⁻¹, or ~1 m day⁻¹) from the Equatorial Undercurrent, producing the cold tongue and high primary productivity of the eastern equatorial Pacific and Atlantic.

Coupling to the Geostrophic Interior and Sverdrup Vorticity Balance

As discovered by Harald Sverdrup (1947) and expanded in the ventilated thermocline theory of James R. Luyten, Joseph Pedlosky, and Henry Stommel (1983), vertical Ekman velocity w_E acts as a moving upper piston on the underlying geostrophic ocean interior. In the inviscid interior of thickness H beneath the Ekman layer, the linearized steady planetary-geostrophic vorticity equation balances advection of planetary vorticity β v_g against vortex-tube stretching f (∂w/∂z):

β v_g = f ( ∂w / ∂z )

Integrating vertically from a level of no motion at depth z = −H (where w(−H) = 0) up to the base of the Ekman layer z = −d_E (where w(−d_E) = w_E) yields the interior geostrophic meridional transport V_g = ∫_−H^−d_E v_g dz:

β V_g = f w_E = ( k · ∇ × τ₀ ) / ρ₀ + ( β τ_0x ) / ( ρ₀ f )

In a subtropical gyre where w_E < 0, vertical compression of interior water columns (∂w/∂z < 0) decreases their height; to conserve barotropic potential vorticity q ≈ f / H, the entire interior column must migrate equatorward (V_g < 0 in the Northern Hemisphere) toward smaller planetary vorticity f. Adding the meridional Ekman transport V_E = −τ_0x/(ρ₀f) to the interior geostrophic transport V_g cancels the β-divergence term and recovers the classic Sverdrup transport equation:

V_total = V_g + V_E = ( k · ∇ × τ₀ ) / ( ρ₀ β )

Quantitative Comparison of Ekman Pumping Regimes in the Ocean and Atmosphere

Charney–Eliassen Formula at a Solid Lower Boundary

At the bottom of the ocean (above the seafloor) or at the bottom of the atmosphere (above land or ocean surfaces), the boundary condition is a no-slip solid boundary (u = 0 at z = 0) rather than a prescribed external wind stress. Across this bottom Ekman layer of characteristic thickness d_E = (2 A_z / f)^1/2 (where A_z is the kinematic eddy viscosity), the horizontal velocity rotates with height along an Ekman spiral from u = 0 at the ground to the interior geostrophic wind or geostrophic bottom current u_g = (u_g, v_g) at the top of the boundary layer (z ≥ d_E).

Integrating the linear Ekman spiral profile vertically from z = 0 to z = ∞ and applying the continuity equation yields the classic Charney–Eliassen (1949) bottom Ekman pumping formula:

w_E = ( d_E / 2 ) ζ_g = ( A_z / 2f )^1/2 ( ∂v_g / ∂x − ∂u_g / ∂y ) = ( A_z / 2f )^1/2 ∇²ψ_g

where ζ_g = ∇²ψ_g is the vertical component of the interior geostrophic relative vorticity immediately above the boundary layer. Physically:

• Beneath a Cyclonic Vortex (ζ_g > 0 in the Northern Hemisphere): Surface friction slows the rotating parcel inside the boundary layer, weakening its outward Coriolis deflection relative to the inward radial pressure gradient force. Fluid spirals inward toward the low-pressure center (∇_h · U_E < 0) and is pumped upward out of the boundary layer (w_E > 0). In the atmosphere, this upward Ekman pumping lifts warm, moist air to its lifting condensation level, generating clouds and precipitation in synoptic low-pressure systems and tropical cyclones.
• Beneath an Anticyclonic Vortex (ζ_g < 0 in the Northern Hemisphere): Frictional cross-isobar flow spirals outward away from the high-pressure center (∇_h · U_E > 0), drawing dry air downward from the free troposphere into the boundary layer (w_E < 0, subsidence) and suppressing cloud formation.

Ekman Spin-Down Time Scale

Because upward Ekman pumping (w_E > 0) beneath a cyclonic vortex of depth H forces vertical compression of the overlying interior fluid column (∂w/∂z ≈ −w_E / H < 0), the interior vortex tubes shrink vertically and spin down via conservation of potential vorticity:

∂ζ_g / ∂t = f ( ∂w / ∂z ) = − ( f w_E ) / H = − [ ( f A_z / 2 )^1/2 / H ] ζ_g = − ζ_g / τ_E

Thus, the interior geostrophic vorticity decays as ζ_g(t) = ζ_g(0) e^−t/τ_E with the Ekman spin-down time scale:

τ_E = H / ( f A_z / 2 )^1/2 = ( 2 H ) / ( f d_E ) = E_V^−1/2 f⁻¹

where E_V = 2A_z/(fH²) = (d_E/H)² is the vertical Ekman number. Notice that τ_E is vastly shorter than the molecular or turbulent viscous diffusion time across the full fluid depth τ_diff = H² / A_z by a factor of d_E / H ≪ 1: secondary circulation via Ekman pumping flushes the entire interior fluid volume through the boundary layer in a fraction of the diffusive time. In the midlatitude troposphere (H ≈ 10 km, d_E ≈ 1 km, f ≈ 10⁻⁴ s⁻¹), the Ekman spin-down time is τ_E ≈ 2 × 10⁵ s ≈ 2.3 to 4 days, governing the decay lifetime of barotropic cyclones after baroclinic instability ceases.

What is the difference between Ekman pumping and Ekman suction?

Both terms refer to the vertical velocity w_E at the edge of an Ekman boundary layer, distinguished by the direction of vertical transport relative to the boundary layer. Ekman pumping refers to downward vertical velocity (w_E < 0 at the base of the oceanic surface layer) driven by horizontal Ekman convergence under anticyclonic wind-stress curl, which pumps surface water down into the ocean interior (or upward out of a bottom Ekman layer into the interior). Ekman suction refers to upward vertical velocity (w_E > 0 at the base of the oceanic surface layer) driven by horizontal Ekman divergence under cyclonic wind-stress curl, which draws deep water upward into the surface mixed layer.

Why does an anticyclonic wind field cause downward motion in the ocean in both hemispheres?

Although the Coriolis effect reverses sign between the Northern Hemisphere (f > 0) and Southern Hemisphere (f < 0), the sense of rotation of an anticyclonic atmospheric high-pressure system also reverses (clockwise in the Northern Hemisphere, with ∇ × τ₀ · k < 0; counterclockwise in the Southern Hemisphere, with ∇ × τ₀ · k > 0). Because Ekman transport is directed 90° to the right of the wind in the Northern Hemisphere and 90° to the left of the wind in the Southern Hemisphere, the ratio (k · ∇ × τ₀) / (ρ₀f) is negative in both hemispheres for an anticyclonic gyre, driving inward mass convergence and downward Ekman pumping (w_E < 0) in all five subtropical ocean gyres.

How do oceanic mesoscale eddies modify classical Ekman pumping?

Classical Ekman theory neglects both the surface ocean current velocity u_o in the bulk wind-stress formula τ₀ = ρ_a C_D |u_a − u_o|(u_a − u_o) and the interior relative vorticity ζ_o compared to f. As shown by Stern (1965), Niiler (1969), and Gaube et al. (2015), across strong oceanic eddies (100–200 km diameter), the curl of the relative wind stress τ(u_a − u_o) and the interaction of wind stress with the horizontal gradient of absolute vorticity (f + ζ_o) generate localized eddy-induced Ekman pumping anomalies of ±0.1 to ±1.0 m day⁻¹—up to an order of magnitude stronger than the background gyre-scale Ekman pumping—that systematically attenuate eddy amplitudes and modify marine chlorophyll distributions.

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.1007/s10236-005-0006-9
2. Sverdrup, H. U. (1947). "Wind-driven currents in a baroclinic ocean; with application to the equatorial currents of the eastern Pacific." Proceedings of the National Academy of Sciences, 33(11), 318–326. https://doi.org/10.1073/pnas.33.11.318
3. Charney, J. G., & Eliassen, A. (1949). "A numerical method for predicting the perturbations of the middle latitude westerlies." Tellus, 1(2), 38–54. https://doi.org/10.3402/tellusa.v1i2.8500
4. Luyten, J. R., Pedlosky, J., & Stommel, H. (1983). "The ventilated thermocline." Journal of Physical Oceanography, 13(2), 292–309. https://doi.org/10.1175/1520-0485(1983)013<0292:TVT>2.0.CO;2013<0292:TVT>2.0.CO;2)
5. Gaube, P., Chelton, D. B., Samelson, R. M., Schlax, M. G., & O'Neill, L. W. (2015). "Satellite observations of mesoscale eddy-induced Ekman pumping." Journal of Physical Oceanography, 45(1), 104–132. https://doi.org/10.1175/JPO-D-14-0032.1