Thermohaline Circulation Strength Estimator

Scale the overturning transport of a thermohaline cell from its density contrast, pycnocline depth and latitude, and see how a freshening of the high latitudes weakens it.

🌐 Thermohaline Circulation Strength Estimator
Density difference (high minus low latitude)
kg/m³
0.13
Change in density difference
%
-50100
Temperature difference (low minus high latitude)
K
020
Salinity contrast (low minus high latitude)
psu
03
Change in salinity contrast (+ = fresher north)
%
-50100
Pycnocline depth (H)
m
2001500
Basin width (W)
km
10008000
Meridional distance (L)
km
100012000
Latitude (for f)
deg
570
Overturning transport
—
Transport in m³/s
—
After the percent change
—
Change in transport
—
Density difference used
—
Reduced gravity (g′)
—
Coriolis parameter (f)
—
Versus observed AMOC
—
Extra
—

🌐 What is the Thermohaline Circulation Strength Estimator?

Thermohaline circulation strength is the volume of water the density-driven overturning cell moves per second, usually quoted in Sverdrups (1 Sv = 1,000,000 m³/s). The name combines temperature and salinity, the two properties that set the density of seawater. Cold or salty water at high latitudes becomes dense, sinks and spreads at depth, and lighter water flows poleward near the surface to replace it. This estimator turns that picture into a number by scaling the transport with the density difference, the depth of the layer involved and the rotation of the Earth.

The calculation is a scaling estimate based on thermal wind balance, not an ocean model. A density difference between low and high latitudes creates a pressure difference, rotation turns that pressure difference into a vertical shear of the flow, and integrating over the pycnocline depth gives a transport that grows with the square of that depth and falls with the Coriolis parameter. Scaling arguments of this kind are used in textbooks and idealised model studies to understand why the overturning responds to buoyancy forcing, and to estimate how large the response might be when high latitudes become fresher.

It is useful for students working through physical oceanography, for climate communicators who want to translate a salinity or density change into a Sverdrup figure, and for modellers who want a quick sanity check on an idealised experiment. Temperature and salinity mode uses a linear equation of state with alpha = 2 x 10-4 per K and beta = 7.6 x 10-4 per psu, the same kind of coefficients used in Stommel's two-box model. That mode also reports the salinity contrast at which temperature and salt effects cancel exactly, so you can see how close a given freshening is to reversing the buoyancy forcing.

It differs from the climate-science AMOC Strength Estimator, which multiplies a basin cross section by a mean northward velocity and converts the result into heat transport. Here the driver is the density difference itself, so the page answers a different question: how much does the transport change if the density difference changes? The comparison with the roughly 17 Sv measured by the RAPID array at 26.5 degrees north is given only as context, because the observed value includes wind-driven and boundary current components that a pure density scaling does not contain.

Two cautions apply. The scaling constant C is of order one and cannot be fixed from first principles, and several studies have argued that pycnocline depth is not a universal control on overturning strength. Treat the result as a guide to sensitivity and orders of magnitude, not as a forecast of any real ocean circulation or of an AMOC collapse.

📐 Formula

Ψ  =  C × g′ × H² × W ÷ (f × L)    with    g′ = g × Δρ ÷ ρ0
Ψ = overturning transport in m³/s (divide by 106 for Sverdrups)
C = order-one scaling constant, taken as 0.5 here (the average of a linearly varying thermal wind profile)
g′ = reduced gravity in m/s², with g = 9.81 m/s² and ρ0 = 1027 kg/m³
Δρ = density difference between high and low latitudes (kg/m³)
H = pycnocline depth (m), W = basin width (m), L = meridional distance over which the density changes (m)
f = 2 Ω sin(φ), the Coriolis parameter, with Ω = 7.2921 × 10-5 rad/s and φ the latitude
Temperature and salinity mode: Δρ = ρ0 (α ΔT − β ΔS), with α = 2 × 10-4 per K and β = 7.6 × 10-4 per psu
Zero-flow salinity contrast: ΔScrit = α ΔT ÷ β. Beyond it Δρ changes sign.
Example: Δρ = 1 kg/m³, H = 600 m, W = 4000 km, L = 6000 km, 26.5° gives Ψ = 1.7614 × 107 m³/s = 17.61 Sv

The scaling is the standard thermal wind estimate (Psi proportional to g′ H² / f). The width W and distance L enter through the horizontal density gradient and the cross section. Stommel's 1961 box model has the same linear flow law, flow proportional to the density difference, but adds salt advection feedbacks that this page does not simulate.

📖 How to Use This Calculator

Steps

1
Choose a mode. Pick Density difference to enter the contrast directly, or Temperature and salinity to build it from alpha times the temperature difference minus beta times the salinity contrast.
2
Enter the geometry and latitude. Set the pycnocline depth in metres, the basin width and meridional distance in kilometres, and the latitude in degrees.
3
Set the percent change. Enter the percent change in the density difference, or in the salinity contrast, to see how the overturning responds to freshwater forcing.
4
Read the transport and chart. Review the transport in Sverdrups and m³/s, the changed value, the comparison with the roughly 17 Sv observed at 26.5 degrees north, and the chart.

💡 Example Calculations

Example 1 — A Subtropical Cell Weakened by a Fifth

Density difference 1 kg/m³, H = 600 m, W = 4000 km, L = 6000 km, latitude 26.5°, change -20%

1
Reduced gravity g′ = 9.81 × 1 ÷ 1027 = 0.009552 m/s², and f = 6.5075 × 10-5 s-1
2
Ψ = 0.5 × 0.009552 × 600² × 4.0 × 106 ÷ (6.5075 × 10-5 × 6.0 × 106) = 1.7614 × 107 m³/s = 17.61 Sv
3
With Δρ = 0.8 kg/m³ the transport is 14.09 Sv, a change of -3.52 Sv (-20.0%)
Transport = 17.61 Sv (104% of the 17 Sv observed), falling to 14.09 Sv
Try this example →

Example 2 — A Deeper, Stronger Mid-Latitude Cell

Density difference 1.4 kg/m³, H = 800 m, W = 3000 km, L = 5000 km, latitude 40°, change -10%

1
g′ = 9.81 × 1.4 ÷ 1027 = 0.013373 m/s², and f = 9.3746 × 10-5 s-1
2
Ψ = 0.5 × 0.013373 × 800² × 3.0 × 106 ÷ (9.3746 × 10-5 × 5.0 × 106) = 2.7389 × 107 m³/s = 27.39 Sv
3
With Δρ = 1.26 kg/m³ the transport is 24.65 Sv, a change of -2.74 Sv (-10.0%)
Transport = 27.39 Sv, about 161% of the observed 17 Sv, falling to 24.65 Sv
Try this example →

Example 3 — Freshening the High Latitudes by 20 Percent

ΔT = 8 K, ΔS = 0.9 psu, H = 600 m, W = 4000 km, L = 6000 km, latitude 26.5°, salinity contrast +20%

1
Δρ = 1027 × (2 × 10-4 × 8 − 7.6 × 10-4 × 0.9) = 0.9407 kg/m³, so g′ = 0.008986 m/s²
2
Ψ = 1.6570 × 107 m³/s = 16.57 Sv, about 97% of the observed 17 Sv
3
The salinity contrast rises to 1.08 psu, Δρ falls to 0.8002 kg/m³ and the transport becomes 14.10 Sv, a change of -2.47 Sv (-14.9%)
A 20% larger salinity contrast cuts the transport by 14.9%, from 16.57 Sv to 14.10 Sv. The zero-flow contrast is 2.105 psu.
Try this example →

Example 4 — Freshening Past the Reversal Point

ΔT = 5 K, ΔS = 0.8 psu, H = 500 m, W = 3500 km, L = 5500 km, latitude 30°, salinity contrast +100%

1
Δρ = 1027 × (2 × 10-4 × 5 − 7.6 × 10-4 × 0.8) = 0.4026 kg/m³, so g′ = 0.003846 m/s² and f = 7.2921 × 10-5 s-1
2
Ψ = 4.1949 × 106 m³/s = 4.19 Sv
3
Doubling the contrast to 1.6 psu exceeds the zero-flow value of 1.316 psu, so Δρ = -0.2218 kg/m³ and Ψ = -2.31 Sv, a change of -6.51 Sv
The transport falls from 4.19 Sv to -2.31 Sv, meaning the buoyancy forcing has reversed sign
Try this example →

❓ Frequently Asked Questions

How do you estimate thermohaline circulation strength?+
A common scaling comes from thermal wind balance. The meridional density difference sets the vertical shear of the flow, and integrating over the pycnocline depth gives a transport Psi = C g′ H² W / (f L), where g′ is g times the density difference divided by a reference density. This calculator uses C = 0.5 and returns the result in Sverdrups.
What is a Sverdrup?+
One Sverdrup (Sv) equals one million cubic metres of water per second. It is the unit used for ocean volume transport because the Gulf Stream, the overturning circulation and the Antarctic Circumpolar Current move volumes far too large to quote conveniently in cubic metres per second.
How strong is the observed Atlantic overturning?+
The RAPID array at 26.5 degrees north reported a mean of about 17 Sv for April 2004 to October 2012, with individual months ranging far above and below that. The calculator uses 17 Sv as a dashed reference line on the chart and in the plain-language note.
What thermal and haline coefficients does the calculator use?+
Temperature and salinity mode uses a linear equation of state with a thermal expansion coefficient alpha = 2e-4 per kelvin, a haline contraction coefficient beta = 7.6e-4 per psu and a reference density of 1027 kg per cubic metre. Real values vary with temperature, pressure and salinity, so this is a simplification.
How does freshening the North Atlantic weaken the overturning?+
Fresher surface water is lighter, so it reduces the density difference between low and high latitudes. Because the scaling is proportional to that density difference, a smaller difference gives a weaker overturning. The calculator shows this by increasing the low-to-high latitude salinity contrast by a percent you choose.
What is the zero-flow salinity contrast?+
It is the salinity contrast at which the thermal and haline effects cancel, alpha times the temperature difference divided by beta. Beyond it the net density difference changes sign, which in Stommel's two-box model corresponds to the reversed, salt-driven circulation branch.
What is Stommel's two-box model?+
Stommel (1961) represented the ocean as a warm low-latitude box and a cold high-latitude box connected by pipes, with flow proportional to the density difference. Because the flow also moves salt, the system can have two stable states, one temperature driven and one salt driven. This page uses only its linear flow law and does not simulate the switching.
Why does the estimate depend on latitude?+
Geostrophic and thermal wind balance involve the Coriolis parameter f = 2 Omega sin(latitude). A larger f means a larger pressure difference is needed to carry the same flow, so the transport scales as 1 / f. Near the equator the balance fails, so latitudes below 5 degrees are rejected.
Why does pycnocline depth enter as a square?+
The thermal wind velocity grows in proportion to the depth of the layer, and the transport is that velocity times the same depth, so it scales as H squared. Doubling the pycnocline depth at a fixed density difference therefore quadruples the estimate, which is why the depth input matters so much.
Is this a climate model or a prediction?+
No. It is a one-line scaling estimate that shows how transport responds to density, depth and latitude. It cannot predict an AMOC collapse, which depends on feedbacks, mixing, wind forcing and ice melt that full models try to represent.
How is this different from the AMOC Strength Estimator?+
The AMOC Strength Estimator multiplies a basin cross section by a mean velocity and converts the result to heat transport. This estimator works from the physical driver instead, the density difference and the pycnocline depth, and shows how transport responds to freshening.
What does a negative transport mean?+
A negative value appears when the changed salinity contrast is larger than the zero-flow value, so the net density difference reverses sign. The scaling then gives flow in the opposite direction, the salt-driven branch of Stommel's picture. It is a warning that the linear estimate has left its range of validity, not a forecast.

How do you estimate thermohaline circulation strength?

A common scaling comes from thermal wind balance. The meridional density difference sets the vertical shear of the flow, and integrating over the pycnocline depth gives a transport Psi = C g' H^2 W / (f L), where g' = g times the density difference divided by a reference density. This calculator uses C = 0.5 and returns the result in Sverdrups.

What is a Sverdrup?

One Sverdrup (Sv) equals one million cubic metres of water per second. It is the unit used for ocean volume transport because the Gulf Stream, the overturning circulation and the Antarctic Circumpolar Current move volumes far too large to quote conveniently in cubic metres per second.

How strong is the observed Atlantic overturning?

The RAPID array at 26.5 degrees north reported a mean of about 17 Sv for April 2004 to October 2012, with individual months ranging far above and below that. The calculator uses 17 Sv as a dashed reference line on the chart and in the plain-language note.

What are the thermal and haline coefficients used?

Temperature and salinity mode uses a linear equation of state with a thermal expansion coefficient alpha = 2e-4 per kelvin, a haline contraction coefficient beta = 7.6e-4 per psu and a reference density of 1027 kg per cubic metre. Real values vary with temperature, pressure and salinity, so this is a simplification.

How does freshening the North Atlantic weaken the overturning?

Fresher surface water is lighter, so it reduces the density difference between low and high latitudes. Because the scaling is proportional to that density difference, a smaller difference gives a weaker overturning. The calculator shows this by increasing the low-to-high latitude salinity contrast by a percent you choose.

What is the zero-flow salinity contrast?

It is the salinity contrast at which the thermal and haline effects cancel, alpha times the temperature difference divided by beta. Beyond it the net density difference changes sign, which in Stommel's two-box model corresponds to the reversed, salt-driven circulation branch.

What is Stommel's two-box model?

Stommel (1961) represented the ocean as a warm low-latitude box and a cold high-latitude box connected by pipes, with flow proportional to the density difference. Because the flow also moves salt, the system can have two stable states, one temperature driven and one salt driven. This page uses only its linear flow law and does not simulate the switching.

Why does the estimate depend on latitude?

Geostrophic and thermal wind balance involve the Coriolis parameter f = 2 Omega sin(latitude). A larger f means a larger pressure difference is needed to carry the same flow, so the transport scales as 1 / f. Near the equator the balance fails, so latitudes below 5 degrees are rejected.

Is this a climate model or a prediction?

No. It is a one-line scaling estimate that shows how transport responds to density, depth and latitude. It cannot predict an AMOC collapse, which depends on feedbacks, mixing, wind forcing and ice melt that full models try to represent.

How is this different from the AMOC Strength Estimator?

The AMOC Strength Estimator multiplies a basin cross section by a mean velocity and converts the result to heat transport. This estimator works from the physical driver instead, the density difference and the pycnocline depth, and shows how transport responds to freshening.