Thermohaline Circulation — The Ocean Conveyor Belt That Regulates Earth's Climate

1. Seawater Density and the T-S Relationship

Thermohaline circulation is density-driven. Seawater density depends on both temperature (T) and salinity (S). Cold, salty water is denser and sinks; warm, fresh water is buoyant and rises.

The T-S diagram (temperature vs salinity plot) is the oceanographer's fingerprint for water masses. Each deep water mass — North Atlantic Deep Water, Antarctic Bottom Water, Mediterranean Outflow Water — occupies a distinct T-S region, allowing passive tracers to track their origin and mixing.

2. Deep Water Formation

New deep water is formed at only a handful of locations worldwide where surface water becomes dense enough to sink through the entire water column — a violent process called convection.

North Atlantic Deep Water (NADW)

In the Labrador Sea, Irminger Sea, and Nordic Seas (between Greenland and Norway), winter cooling and ice formation both act to densify surface water:

• Cooling: reduces temperature → increases density
• Brine rejection: when sea ice forms, dissolved salt is expelled into the remaining liquid, raising salinity (and density) of the surrounding water

The resulting plumes of dense water sink to 2,000–4,000 m depth and flow southward through the Atlantic as NADW — the main driver of the AMOC.

Antarctic Bottom Water (AABW)

The densest natural water in the global ocean forms on the Antarctic continental shelves, particularly in the Weddell and Ross Seas. Brine rejection under sea ice produces water cold enough (~−1.9°C) and saline enough (~34.7 psu) to sink all the way to the ocean floor (~5,000 m) and spread northward beneath NADW in all ocean basins.

3. The Global Ocean Conveyor (MOC)

Wallace Broecker popularized the image of the "ocean conveyor belt" in 1991. In his simplified picture, a single loop connects the North Atlantic sinking region with upwelling in the Southern Ocean and Pacific:

1. Warm surface water flows north in the Atlantic, losing heat to the atmosphere
2. Dense water sinks in the North Atlantic and Nordic Seas → forms NADW
3. NADW flows south along the ocean floor, beneath the equator
4. It joins Antarctic Circumpolar Current (ACC) — the only uninterrupted circumpolar flow
5. Upwelling returns deep water to the surface over decades in the Southern Ocean and Pacific
6. Shallow warm return flow carries heat back toward the Northern Hemisphere

4. The Atlantic Meridional Overturning Circulation (AMOC)

The AMOC is the Atlantic branch of the MOC. It carries roughly 17–19 Sverdrups (1 Sv = 10⁶ m³/s) of water northward in the upper ocean — equivalent to ~100 Amazon rivers — and returns the same volume at depth.

Heat Transport

The AMOC carries approximately 1.3 petawatts (PW) of heat northward across 26°N — about 25% of the total poleward energy transport at that latitude. This heat flux warms the North Atlantic region by 5–10°C above what orbital forcing alone would predict.

RAPID Array Monitoring

Since 2004, the RAPID–MOCHA array of moorings at 26.5°N has continuously measured AMOC strength. Key findings:

• Mean strength ~17 Sv; significant interannual variability ±3 Sv
• Notable slowdown in 2009–2010 (−30%): cold European winter of 2009/10 followed by ~6 months
• Long-term weakening trend: AMOC ~15% weaker in 2004–2023 than pre-industrial estimates from sediment proxies

5. Climate Impacts of the AMOC

Why Britain Is Warmer Than Labrador

London (51°N) and Calgary (51°N) share the same latitude. London's January mean is +5°C; Calgary's is −8°C. The AMOC transports heat equivalent to ~5 million power plants to the North Atlantic — replacing the radiative deficit of high-latitude Europe.

Sea Level

A weaker AMOC reduces the southward geostrophic pressure gradient that normally depresses sea level along the US East Coast. A 30% AMOC weakening would add ~20–30 cm of sea level rise to Boston, New York, and Miami — on top of thermal expansion from warming.

Monsoon and Tropical Rainfall

A cooler North Atlantic shifts the Intertropical Convergence Zone (ITCZ) southward, weakening the West African and Indian monsoons and intensifying drought in the Sahel. Paleoclimate evidence from Dansgaard-Oeschger events shows that past AMOC collapses triggered near-global climate reorganization within decades.

6. AMOC Slowdown and Collapse Risk

Why Freshwater Is the Threat

Accelerating Greenland ice sheet melt (currently ~280 Gt/yr) and increased Arctic river runoff inject fresh (low-density) water into the North Atlantic. This freshwater cap suppresses deep convection — the density excess that drives NADW formation is neutralized.

The Stommel Bifurcation — Two Stable States

Henry Stommel (1961) showed analytically that a simple 2-box ocean model has two stable equilibria: an "on" state (strong circulation) and an "off" state (weak or collapsed circulation). The system can flip between them catastrophically beyond a freshwater forcing threshold:

IPCC Assessment

The IPCC AR6 (2021) assessed an abrupt collapse this century as unlikely but not impossible — implying a 5–10% probability under high-emission scenarios. A 2023 study using early-warning statistical signals suggested the AMOC may be approaching a critical transition point, though this interpretation remains contested.

7. Simple Box Model in JavaScript

Stommel's 2-box model can be integrated as a pair of ODEs to visualize the two stable states and the bifurcation:

// Stommel two-box thermohaline model
const k = 1.5e-6; // hydraulic constant [m³/(s·kg)]
const alpha = 1.5e-4; // thermal expansion [K⁻¹]
const beta  = 8e-4;  // haline contraction [psu⁻¹]
const gamma = 1 / (60 * 86400); // restoring rate to atmosphere (1/month)

function stommelStep(T1, S1, T2, S2, H, dt) {
  // H = freshwater flux into polar box [psu/s]
  const q = k * (alpha * (T1 - T2) - beta * (S1 - S2));

  // Salinity equations (temperature relaxes quickly to atmospheric forcing)
  const dS1 = Math.abs(q) * (S2 - S1) + H; // salt import from polar box
  const dS2 = Math.abs(q) * (S1 - S2) - H; // freshwater + flow

  return {
    S1: S1 + dS1 * dt,
    S2: S2 + dS2 * dt,
    q_Sv: q * 1e-6 // convert to Sverdrups
  };
}

// Try H = 0 → stable ON state (~17 Sv)
// Try H = 3e-9 → AMOC collapses to near-zero / reversal