A coupled ocean-atmosphere climate model: temperature versus salinity effects on the thermohaline circulation

A simple nonlinear three-box ocean model of the North Atlantic Ocean including the rudiments of eddy mixing, vertical stratification and thermohaline circulation is first presented. It is subject to uniform latitudinal differential heating, q, and net evaporation m sub(e) , and includes a linear equation of state. Two quite different limiting steady-state solutions exist. The first has a warm saline surface water and a cold, low-salinity deep ocean; deep water is primarily formed in higher latitudes by the prevalence of differential heating. A second limiting solution consists of a warm saline deep ocean underlying a cool, low-salinity surface ocean; deep water is formed primarily in lower latitudes as a consequence of large differential evaporation. A coupled ocean-atmosphere model, in which the oceanic surface heat fluxes are determined internally but with differential evaporation at the ocean surface m sub(e) remaining an external parameter, is next presented. The atmosphere component is a simple energy balance model that emphasizes the vertical fluxes of radiative, sensible and latent heat fluxes but does not include temperature-albedo feedback. Model response depends on the external parameters m sub(e) and mu , controlling the magnitude of the thermohaline-driven circulation, and on the magnitudes of the eddy mixing coefficients and the solar constant. For small m sub(e) , a steady-state solution corresponding to a cold fresh deep ocean is found, qualitatively similar to the modern ocean. For large m sub(e) , a steady-state solution with a warm saline deep ocean occurs; this solution resembles conceptual models that have been proposed for the warm saline Cretaceous ocean. There exists an intermediate region of values of m sub(e) for which the solutions are more complex. On the lower end of this region, both the cold fresh deep-ocean and warm saline deep-ocean circulations coexist as stable equilibria. On the upper end, the cold-deep ocean becomes unstable, manifesting oscillations with growing amplitude, and ultimately reaches the warm saline deep-ocean solution. In the neighborhood of a `cusp' on the mu , m sub(e) plane, that is, for relatively small mu , more complex behaviour occurs, which has not yet been fully analyzed. The model response in the region of complexity is not sensitive to changes in the solar constant but is sensitive to the eddy mixing coefficients.