Physical grounding of the energy-constraint theory of thermohaline circulation: threshold robustness beyond conceptual models, an energy-to-mixing parameterization, and reconciliation with buoyancy-driven dynamics

Conceptual models that replace Stommel's buoyancy constraint on the thermohaline circulation (THC) with an energy constraint — the overturning scale is set by external mechanical energy sustaining diapycnal mixing — predict that saline-mode multiple equilibria exist only between a freshwater-flux threshold and a mixing-energy threshold, with abrupt transitions to a thermal mode when either is crossed. These results rest entirely on two-box, three-box, and two-layer conceptual models with an imposed closure and fixed parameter choices. The open problem is to give this energy-constraint view a physical and model-level foundation: test whether the threshold and abrupt-transition structure survives in higher-complexity configurations, construct a physically grounded parameterization linking external mechanical energy (wind and tidal input) to mixing and overturning in stratified three-dimensional oceans, and determine under what realistic conditions the buoyancy-driven and energy-constrained viewpoints each dominate, formulating a closure that represents both jointly.

Background

The thermohaline circulation (THC) is the global-scale overturning of the ocean driven by differences in heat and freshwater at the sea surface; in the Atlantic it appears as the Atlantic Meridional Overturning Circulation (AMOC), with northward flow of warm surface water and southward return flow of cold deep water. Stommel (1961) introduced the canonical conceptual description: two well-mixed boxes exchanging water, heat, and salt, with the exchange velocity set by the density (buoyancy) difference between the boxes. This 'buoyancy constraint' embodies the assumption that surface thermohaline forcing drives the overturning. Depending on whether temperature or salinity dominates the density contrast, the circulation is in a thermal mode (temperature-dominated, heat transported poleward) or a saline mode (salinity-dominated), and Stommel's model admits two stable equilibria — the origin of the idea that the AMOC can undergo abrupt transitions, a recurring concern for climate dynamics (Weijer et al. 2019).

A distinct line of work emphasizes energetics. Because the ocean is stably stratified, mechanical energy is required to lift dense water and sustain diapycnal mixing (mixing across density surfaces) against gravity; the required power — on the order of 2 terawatts — is supplied not by surface heating but by winds and tides (Wunsch and Ferrari 2004). Sandström's-theorem-type arguments — Sandström's theorem is the classical result that heating and cooling applied at the same pressure level cannot by themselves drive a steady closed circulation, so an external source of mechanical work is required — and laboratory analogues (Wang and Huang 2005) support the view that the ocean is not a heat engine: without external mechanical stirring the overturning would spin down. Huang (1999) showed in a general circulation model that the poleward heat flux scales nearly linearly with the external energy flux driving vertical mixing. Guan and Huang (2008) recast Stommel's two-box model by replacing the buoyancy constraint with an energy constraint: the overturning velocity takes the form $w = e/[\rho_0\alpha(T_1-T_2)-\rho_0\beta(S_1-S_2)]$, where $e$ is the mixing-energy parameter (proportional to the external mechanical energy sustaining mixing), $\rho_0$ the mean density, $\alpha$ and $\beta$ the thermal and saline expansion coefficients, and $T_i, S_i$ the box temperatures and salinities — i.e. the overturning strengthens with available mixing energy and with weaker density contrast, opposite to the buoyancy-constraint scaling. In this framework the multiple equilibria of Stommel-type models can be eliminated by strong wind-driven gyres, and freshwater forcing can strengthen rather than slow the overturning (Nilsson and Walin 2001, in a two-layer model with fixed mixing energy).

Shen, Guan, Liang and Chen (2011) extended the energy-constraint two-box model to three boxes including a wind-driven gyre, and Shen and Guan (2015) formulated a two-layer (vertically stratified) conceptual model under the energy constraint, with a freshwater flux $p$ between the layers and the mixing-energy parameter $e$. By computing equilibria and classifying their linear stability via the eigenvalues of the Jacobian matrix (the matrix of first derivatives of the layer evolution equations) at fixed parameter choices, they found: (i) in thermal mode the overturning rate and meridional heat transport increase with $p$ and with $e$; (ii) in saline mode, bistability (the stable saline mode coexists with the stable thermal mode, separated by an unstable equilibrium) exists only for $p$ above a threshold and $e$ below a threshold; (iii) crossing either threshold produces an abrupt transition from the stable saline to the stable thermal mode; and (iv) the two thresholds interact — larger $p$ raises the critical $e$, larger $e$ raises the critical $p$. All of this, however, is confined to a conceptual model whose energy closure was derived for a horizontally partitioned box geometry and then imposed on the layered configuration, with thresholds computed numerically for specific parameter sets. Whether this structure is a robust property of the real ocean — or an artifact of the idealization — depends on three linked questions that the paper's own framework leaves open.

Problem Statement

Establish whether the energy-constraint theory of the thermohaline circulation is a robust property of the ocean rather than an artifact of conceptual models, by addressing three linked questions arising from the demonstrated limitations of Shen and Guan (2015) — its conceptual two-layer scope, its energy closure imported from a horizontally partitioned box geometry and never validated for a stratified three-dimensional ocean, and its thresholds computed for specific parameter choices: (1) Threshold robustness — determine whether the multiple-equilibria structure predicted by the energy-constraint closure (saline-mode bistability existing only between a freshwater-flux threshold $p_c$ and a mixing-energy threshold $e_c$, with abrupt transitions to a stable thermal mode when either threshold is crossed, and each threshold rising with the other control parameter) survives in a configuration of higher complexity than a small number of well-mixed boxes — one that resolves vertical stratification or basin geometry. (2) Energy-to-mixing parameterization — construct a physically grounded, spatially resolved closure that determines the effective mixing-energy parameter from the external mechanical-energy sources (wind and tidal input, of order 2 TW total) and their conversion to diapycnal mixing in a stratified three-dimensional ocean, and test it in higher-complexity models and against observations, replacing the current purely conceptual, approximately estimated parameter. (3) Reconciliation of constraints — establish under what realistic conditions the Stommel-type buoyancy constraint and the energy constraint each control the overturning, and formulate a closure that represents both influences jointly rather than treating them as mutually exclusive model choices. These deliverables jointly serve one scientific objective: a physically validated account of what drives and sustains the overturning, and of whether the abrupt saline–thermal mode transitions it predicts under changing freshwater forcing and mixing energy are credible features of the real ocean circulation. The three linked questions are logically independent: none presupposes the answer to another.

References

1. Stommel, H. (1961). Thermohaline convection with two stable regimes of flow. Tellus 13:224–230. DOI 10.1111/j.2153-3490.1961.tb00079.x, https://doi.org/10.1111/j.2153-3490.1961.tb00079.x
2. Nilsson, J., and Walin, G. (2001). Freshwater forcing as a booster of thermohaline circulation. Tellus A 53:629–641. DOI 10.3402/tellusa.v53i5.12232, https://doi.org/10.3402/tellusa.v53i5.12232
3. Huang, R. X. (1999). Mixing and energetics of the oceanic thermohaline circulation. Journal of Physical Oceanography 29:727–746. DOI 10.1175/1520-0485(1999)029<0727:MAEOTO>2.0.CO;2, https://doi.org/10.1175/1520-0485(1999)029%3C0727%3AMAEOTO%3E2.0.CO%3B2
4. Wunsch, C., and Ferrari, R. (2004). Vertical mixing, energy, and the general circulation of the oceans. Annual Review of Fluid Mechanics 36:281–314. DOI 10.1146/annurev.fluid.36.050802.122121, https://doi.org/10.1146/annurev.fluid.36.050802.122121
5. Wang, W., and Huang, R. X. (2005). An experimental study on thermal circulation driven by horizontal differential heating. Journal of Fluid Mechanics 540:49–73. DOI 10.1017/S002211200500577X, https://doi.org/10.1017/S002211200500577X
6. Guan, Y. P., and Huang, R. X. (2008). Stommel's box model of thermohaline circulation revisited — The role of mechanical energy supporting mixing and the wind-driven gyration. Journal of Physical Oceanography 38:909–917. DOI 10.1175/2007JPO3535.1, https://doi.org/10.1175/2007JPO3535.1
7. Fürst, J. J., and Levermann, A. (2011). A minimal model for wind- and mixing-driven overturning: threshold behavior for both driving mechanisms. Climate Dynamics 38:239–260. DOI 10.1007/s00382-011-1003-7, https://doi.org/10.1007/s00382-011-1003-7
8. Shen, Y., Guan, Y. P., Liang, C. J., and Chen, D. K. (2011). A three-box model of thermohaline circulation under the energy constraint. Chinese Physics Letters 28:059201. DOI 10.1088/0256-307X/28/5/059201, https://doi.org/10.1088/0256-307X/28/5/059201
9. Shen, Y., and Guan, Y. P. (2015). Feature of thermohaline circulation in two-layer conceptual model based on energy constraint. Science China Earth Sciences 58:1397–1403. DOI 10.1007/s11430-015-5092-8, https://doi.org/10.1007/s11430-015-5092-8
10. Weijer, W., Cheng, W., Drijfhout, S. S., Fedorov, A. V., Hu, A., Jackson, L. C., Liu, W., McDonagh, E. L., Mecking, J. V., and Zhang, J. (2019). Stability of the Atlantic Meridional Overturning Circulation: A Review and Synthesis. Journal of Geophysical Research: Oceans 124:5336–5375. DOI 10.1029/2019JC015083, https://doi.org/10.1029/2019JC015083
11. MacKinnon, J. A., Zhao, Z., Whalen, C. B., Waterhouse, A. F., Trossman, D. S., Sun, O. M., St. Laurent, L. C., Simmons, H. L., Polzin, K., Pinkel, R., Pickering, A., Norton, N. J., Nash, J. D., Musgrave, R., Merchant, L. M., Melet, A. V., Mater, B., Legg, S., Large, W. G., Kunze, E., Klymak, J. M., Jochum, M., Jayne, S. R., Hallberg, R. W., Griffies, S. M., Diggs, S., Danabasoglu, G., Chassignet, E. P., Buijsman, M. C., Bryan, F. O., Briegleb, B. P., Barna, A., Arbic, B. K., Ansong, J. K., and Alford, M. H. (2017). Climate Process Team on Internal Wave–Driven Ocean Mixing. Bulletin of the American Meteorological Society 98:2429–2454. DOI 10.1175/BAMS-D-16-0030.1, https://doi.org/10.1175/BAMS-D-16-0030.1
12. Yu, Y., Liu, H., Lin, P., and Lan, J. (2020). The impact of oceanic processes on the transient climate response: a tidal forcing experiment. Acta Oceanologica Sinica 39:52–62. DOI 10.1007/s13131-019-1466-0, https://doi.org/10.1007/s13131-019-1466-0
13. Liu, Z., Zhang, S., Shen, Y., Guan, Y. P., and Deng, X. (2021). A study of capturing Atlantic meridional overturning circulation (AMOC) regime transition through observation-constrained model parameters. Nonlinear Processes in Geophysics 28:481–500. DOI 10.5194/npg-28-481-2021, https://doi.org/10.5194/npg-28-481-2021
14. Zhou, X., Yang, K., Yang, H., and Zhang, Q. (2026). Self-sustained multicentennial oscillation of the Atlantic meridional overturning circulation in two-hemisphere box models. Climate Dynamics 64:30. DOI 10.1007/s00382-025-07947-7, https://doi.org/10.1007/s00382-025-07947-7