Sverdrup Balance

Sverdrup Balance is a foundational dynamical relation in physical oceanography and geophysical fluid dynamics derived by the Norwegian oceanographer and meteorologist Harald Ulrik Sverdrup in 1947. It states that across the vast, frictionless interior of an ocean basin—away from western boundary currents and topographic obstacles—the depth-integrated meridional (north–south) mass or volume transport V is governed directly by the vertical component of the curl of the surface wind stress τ divided by the meridional gradient of the planetary Coriolis effect parameter β = ∂f/∂y:

βV = (1 / ρ₀) curl_z(τ) = (1 / ρ₀) (∂τ_y / ∂x − ∂τ_x / ∂y)

Here, V = ∫_−H⁰ v dz is the vertically integrated northward volume transport per unit zonal width (with SI units of m²/s, and basin-integrated transport measured in Sverdrups, where 1 Sv = 10⁶ m³/s), ρ₀ ≈ 1,025 kg/m³ is the reference seawater density, f = 2Ω sin(φ) is the Coriolis parameter at latitude φ, β = (2Ω cos φ) / R_E ≈ 2 × 10⁻¹¹ m⁻¹ s⁻¹ is the Rossby parameter on a mid-latitude β-plane, and (τ_x, τ_y) are the zonal (eastward) and meridional (northward) components of wind stress acting on the sea surface.

By linking the curl of atmospheric winds directly to deep geostrophic meridional drift, Sverdrup balance provided the first quantitative explanation for the mid-latitude subtropical and subpolar ocean gyres, the existence of eastward equatorial countercurrents, and the necessity of narrow, intense western boundary currents such as the Gulf Stream and Kuroshio (completed by Henry Stommel in 1948 and Walter Munk in 1950 via the Munk layer).

Linearized Momentum Equations on the Beta-Plane

Consider a steady, Boussinesq, hydrostatic ocean on a mid-latitude β-plane (f = f₀ + βy) where horizontal friction is negligible and vertical turbulent stress ∂τ/∂z is confined to the upper wind-forced Ekman transport layer of thickness d_E ≈ 30–100 m. The linearized horizontal momentum equations and the incompressible continuity equation are:

−fv = −(1 / ρ₀)(∂p / ∂x) + (1 / ρ₀)(∂τ_x / ∂z)

fu = −(1 / ρ₀)(∂p / ∂y) + (1 / ρ₀)(∂τ_y / ∂z)

∂u / ∂x + ∂v / ∂y + ∂w / ∂z = 0

Cross-differentiating the two momentum equations—taking ∂/∂x of the meridional momentum equation and subtracting ∂/∂y of the zonal momentum equation—eliminates the horizontal pressure gradient terms ∂²p/(∂x∂y) completely. Substituting the horizontal divergence ∂u/∂x + ∂v/∂y = −∂w/∂z from the continuity equation and noting that ∂f/∂y = β while ∂f/∂x = 0 yields the linear planetary vorticity equation:

βv = f (∂w / ∂z) + (1 / ρ₀) (∂ / ∂z) [∂τ_y / ∂x − ∂τ_x / ∂y]

Depth Integration and Ekman Pumping Coupling

Integrating the vorticity equation vertically from a deep level of no motion z = −H (where both vertical velocity w(−H) = 0 and turbulent stress τ(−H) = 0 vanish, above the deep abyssal seafloor) up to the sea surface z = 0 (where rigid-lid vertical velocity w(0) = 0 and surface wind stress equals τ) collapses the two right-hand terms into the surface wind-stress curl:

β ∫_−H⁰ v dz = βV = (1 / ρ₀) curl_z(τ)

To understand the physical mechanism of Sverdrup balance, one splits the total meridional transport V = V_E + V_g into the surface Ekman transport (V_E = −τ_x / [ρ₀f]) and the underlying interior geostrophic current transport (V_g):

1. Subtropical Gyre Convergence (Ekman Downwelling w_E < 0): Between the easterly Trade Winds (~15°N) and the mid-latitude Westerlies (~45°N), the surface wind stress exerts a clockwise (anticyclonic) torque on the ocean surface, so curl_z(τ) < 0 in the Northern Hemisphere. Surface Ekman transports converge toward the center of the subtropical gyre, driving downward Ekman pumping (w_E = curl_z[τ] / [ρ₀f] < 0) at the base of the Ekman layer (z = −d_E).
2. Vorticity Compression and Equatorward Drift: In the frictionless geostrophic interior below z = −d_E, the vorticity balance reduces to βv_g = f ∂w/∂z. Because downward Ekman pumping (w_E < 0 at the top, w = 0 at z = −H) vertically squashes interior water columns (∂w/∂z < 0), conservation of angular momentum (potential vorticity f/h) requires the squashed water columns to move toward lower planetary vorticity f—that is, equatorward (v_g < 0, V < 0).
3. Subpolar Gyre Divergence (Ekman Upwelling w_E > 0): At subpolar latitudes (50°N–65°N), cyclonic wind stress curl (curl_z[τ] > 0) drives surface Ekman divergence and upward Ekman suction (w_E > 0), stretching interior water columns (∂w/∂z > 0) and forcing poleward Sverdrup transport (V > 0).

Defining a horizontal volume transport streamfunction Ψ(x, y) such that U = −∂Ψ/∂y and V = ∂Ψ/∂x, Sverdrup balance becomes a first-order partial differential equation in x:

∂Ψ / ∂x = V(x, y) = [1 / (ρ₀β)] curl_z(τ)

Because this equation is only first-order in x, it can satisfy a no-normal-flow boundary condition (Ψ = 0) at only one of the two meridional continental boundaries (either the eastern boundary x = x_E or the western boundary x = x_W, but not both). Because planetary Rossby waves carry long-wavelength energy westward at group velocity c_gx = −βR_d², the Sverdrup interior solution integrates westward from the eastern boundary x = x_E:

Ψ_Sverdrup(x, y) = − ∫_x^x_E [1 / (ρ₀β)] curl_z(τ(x′, y)) dx′

Across the 5,000-km width of the subtropical North Atlantic or North Pacific, this westward integration accumulates ~30 Sv of southward interior transport. Mass conservation requires an equal and opposite +30 Sv northward return flow inside a frictional western boundary layer (Stommel's bottom-friction layer or the lateral-viscosity Munk layer), where frictional vorticity torque balances the planetary vorticity gain βV.

Why does water in the subtropical gyre interior flow southward even where the winds blow eastward or westward?

Sverdrup balance is a statement of angular momentum (vorticity) conservation rather than direct downwind drag. Clockwise winds over the subtropical gyre push surface Ekman water inward toward the gyre center, where it sinks (Ekman downwelling) and vertically compresses the underlying water column. On a rotating spherical Earth where the Coriolis parameter f increases northward (β > 0), a compressed water column can only conserve its potential vorticity (f/H) by drifting southward toward the equator where f is smaller.

Why must the Sverdrup streamfunction be integrated from the eastern boundary rather than the western boundary?

Baroclinic and barotropic long Rossby waves—which adjust the ocean's density and pressure field toward steady Sverdrup balance—always have a westward phase and group velocity. Consequently, boundary information propagates from the eastern continental coast westward across the interior, leaving the western boundary to close the mass budget through a fast, narrow frictional boundary current.

Where in the global ocean does Sverdrup balance break down?

Sverdrup balance holds remarkably well across the upper 1,000–2,000 m of subtropical and tropical ocean interiors (confirmed by Argo float arrays and satellite scatterometry), but breaks down in western boundary currents (where lateral eddy viscosity and nonlinear advection dominate), in the Antarctic Circumpolar Current (where deep flow interacts strongly with seafloor topography via bottom pressure torque), and in the deep abyssal thermohaline circulation.

1. 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
2. Stommel, H. (1948). "The Westward Intensification of Wind-Driven Ocean Currents." Transactions, American Geophysical Union, 29(2), 202–206. https://doi.org/10.1029/TR029i002p00202
3. Wunsch, C. (2011). "The Decadal Mean Ocean Circulation and Sverdrup Balance." Journal of Marine Research, 69(2–3), 417–434. https://doi.org/10.1357/002224011798765303
4. Gray, A. R., & Riser, S. C. (2014). "A Global Analysis of Sverdrup Balance Using Absolute Geostrophic Velocities from Argo." Journal of Physical Oceanography, 44(4), 1213–1229. https://doi.org/10.1175/JPO-D-12-0206.1