CLIMATE AND THE OCEAN CIRCULATION 111 . THE OCEAN MODEL

The ocean model used in a calculation of the earth's climate is described in detail. Compared with earlier numerical models used in ocean circulation studies, the present model includes several new features. Temperature and salinity are treated separately. Density is calculated with an accurate equation of state for sea water. The model also includes a method for calculating the growth and movement of sea ice. Due to the very slow adjustment of the deep water in the ocean model, a numerical integration extending over the equivalent of a century fails t o reach a climatic equilibrium. At the termination of the run, the surface layers of the ocean show little change with respect t o time, but the average heating rate for the ocean as a whole is 2' per century. The salinity patterns at the termination of the run are highly realistic compared t o observations. A halocline forms in the Arctic Zone and a surface salinity maximum is present in the subtropics. A weak salinity minimum a t a depth of 1 km indicates an extensive water mass very similar to the Antarctic intermediate water of the Southern Hemisphere. Poleward heat transport is found to be closely related to the intensity of the thermohaline circulation. A vertical mixing coefficient, K, of 1.5 cm2 sec-* leads to very reasonable heat exchange with the atmospherc based on estimates of the heat balance of the North Atlantic. The calculation indicates that the thermal " relaxation " time of the ocean is too long for a numerical integration of the time-dependent equations to be a practical method of finding an equilibrium solution, and new methods should 806 807 808 808 809 8 10 810 813 814 819 822 822 823 826 826