Century-scale variability in a randomly forced,two-dimensional thermohaline ocean circulation model

The response of a two-dimensional thermohaline ocean circulation model to a random freshwater flux superimposed on the usual mixed boundary conditions for temperature and salinity is considered. It is shown that for a wide range of vertical and horizontal diffusivities and a box geometry that approximates the Atlantic Ocean, 200–300 yr period oscillations exist in the basic-state, interhemispheric meridional overturning circulation with deep convection in the north. These fluctuations can also be described in terms of propagating salinity anomalies which travel in the direction of the thermohaline flow. For large horizontal (K _h = 15 × 103 m²/s) and small vertical (K _v = 0.5 × 10^–4 m²/s) diffusivities, the random forcing also excites deca-millennial oscillations in the basic structure of the thermohaline circulation. In this case, the meridional circulation pattern slowly oscillates between three different stages: a large positive cell, with deep convection in the North Atlantic and upwelling in the south; a symmetric two-cell circulation, with deep convection in both polar regions and upwelling near the equator; and a large negative cell, with deep convection in the South Atlantic and upwelling in the north. Each state can persist for 0 (10 kyr).