Understanding the Equation

dζ/dt + J(ψ,ζ) + βv = curl(τ)/ρH − rζ + A∇⁴ψ

This single equation governs the large-scale circulation of the ocean. It explains why the Gulf Stream exists, why it's on the western side of the Atlantic, and why ocean currents form the patterns they do.

Let's break it apart, piece by piece.

The Streamfunction ψ

The pressure map of the ocean

ψ — streamfunction

ψ is a map of ocean pressure. Water flows along lines of constant ψ, just like contour lines on a topographic map. Higher ψ means higher pressure.

Click to place pressure highs (left click) and lows (right click). Watch how flow arrows appear perpendicular to the gradient — water doesn't flow downhill, it flows along the contours.

Think of ψ like a topographic map of the ocean surface. Gravity wants water to flow downhill, but the Coriolis effect deflects it to the right (in the Northern Hemisphere), so it ends up flowing ALONG the contours, not across them. This is called geostrophic balance.

Vorticity ζ = ∇²ψ

Spin in the fluid

ζ = ∇²ψ — relative vorticity

Vorticity is spin. Positive ζ means counterclockwise rotation, negative means clockwise. Mathematically, it's the Laplacian of ψ — the curvature of the pressure field.

Click and drag to create spinning fluid. Blue = counterclockwise (positive), red = clockwise (negative).

Drop a leaf in a river. If it rotates, there's vorticity. Hurricanes have enormous vorticity. So do bathtub drains. Any time fluid is spinning, vorticity is at work.

The β Effect: βv

Why western boundary currents exist

βv — planetary vorticity gradient

β measures how the Coriolis force changes with latitude. As water moves north, it gains planetary vorticity because the Earth's rotation projects more strongly onto the local vertical at higher latitudes.

Watch the column of water: as it moves north (drag up), it spins up. Moving south, it spins down. This asymmetry between east and west is everything.

This single term — β — is why the Gulf Stream is on the WEST side of the Atlantic, not the east. Henry Stommel showed this in 1948. Remove β from the equation and you get a symmetric gyre with no western intensification. It's the most important result in physical oceanography.

Wind Forcing: curl(τ)/ρH

The engine that drives ocean gyres

curl(τ)/ρH — wind stress curl

Wind pushes the ocean surface. But it's not the wind itself that drives the gyres — it's the curl of the wind stress. The spatial variation in wind is what matters.

Drag to change the wind pattern. The curl (shown as color) appears where wind speed or direction varies in space. Trade winds blow west near the equator, westerlies blow east at mid-latitudes — the curl between them drives a clockwise subtropical gyre.

Imagine rubbing the top of a book to the right and the bottom to the left. The book spins. That rotation comes not from either push alone, but from the difference between them — the curl.

Friction: −rζ

The brake on ocean spin

−rζ — bottom friction

Bottom friction removes vorticity from the system. Without it, wind would inject vorticity forever and the ocean would spin up without limit.

Watch the spinning disk slow down as friction acts. The key: friction is proportional to vorticity itself, so it's strongest where the flow is most intense — in the western boundary current.

In Stommel's model, friction is what closes the vorticity budget. Wind injects vorticity everywhere across the gyre. That vorticity must be removed somewhere. Friction does it — but only in the narrow western boundary layer where vorticity is intense enough. The Gulf Stream IS the western boundary friction layer.

Viscosity: A∇⁴ψ

Smoothing sharp gradients

A∇⁴ψ — lateral viscosity (biharmonic)

Lateral viscosity acts like diffusion — it smooths out sharp gradients. Without it, the western boundary current would be infinitely thin, a mathematical discontinuity.

Watch as a sharp gradient gets smoothed over time. Munk (1950) showed that this viscous term creates a boundary layer of realistic width (~100 km).

Stir cream into coffee. At first there's a sharp boundary between cream and coffee. Molecular diffusion (and turbulent mixing) gradually smooth that boundary until everything is uniform. Viscosity does the same thing to velocity gradients in the ocean.

The Nonlinear Term: J(ψ,ζ)

The chaos term

J(ψ,ζ) — advection of vorticity

The Jacobian J(ψ,ζ) means the flow carries its own vorticity. The ocean current transports the very spin that defines the current. This is what makes the equation nonlinear.

Watch particles get advected by a flow that they themselves influence. Small perturbations grow into eddies and meanders.

This is why ocean currents meander and shed eddies. Without this term, the equation is linear and the solution is smooth and steady. With it, you get instabilities, turbulence, and the chaotic beauty of real ocean circulation. It's also why we need computers — no analytical solution exists.

Putting It Together

The complete vorticity balance

dζ/dt + J(ψ,ζ) + βv = curl(τ)/ρH − rζ + A∇⁴ψ

The equilibrium of the ocean is a vorticity budget:

Wind injects vorticity everywhere across the basin (the curl of the wind stress).

β redistributes that vorticity, pushing it westward.

Friction removes it — but only where vorticity is concentrated enough, which is in the narrow western boundary layer.

Viscosity gives that boundary layer a finite width.

Nonlinearity creates the eddies, meanders, and instabilities that make the real ocean complex and beautiful.

The Gulf Stream, the Kuroshio, the Agulhas — every western boundary current on Earth exists because of this balance. The wind pushes, β concentrates, friction dissipates.

Now go play with the simulation →