Comparison of thermally driven circulations from a depth-coordinate model and an isopycnal-layer model. Part II: The difference and structure of the circulations
Yg. Park et K. Bryan, Comparison of thermally driven circulations from a depth-coordinate model and an isopycnal-layer model. Part II: The difference and structure of the circulations, J PHYS OCEA, 31(9), 2001, pp. 2612-2624
Thermally driven ocean circulations in idealized basins are calculated with
two well-known model codes. one based on depth-level coordinates and the o
ther based on isopycnal coordinates. In addition, the two models have very
different representations of convection. In the level-coordinate model. con
vective adjustment is used, while in the isopycnal-coordinate model. convec
tion is simulated by a transformation of the surface layer to the layer bel
ow. Both models indicate a three-layer structure in the circulation. The lo
wer and middle layers have a flow structure that corresponds with the class
ical abyssal circulation models. The upper flow is strongly constrained by
the buoyancy flux field at the upper surface and the convective parameteriz
ation. The model with convective adjustment and level coordinates is domina
ted by an eastward flow, which sinks to subsurface level at the eastern bou
ndary. It lacks any indication of a surface cyclonic flow, even in the vici
nity of sinking at the northern wall. On the other hand, in the model based
on density coordinates the eastward surface flow turns to the north at the
eastern boundary and forms a pronounced cyclonic circulation at high latit
udes. Due to the cyclonic circulation, the coldest surface water is found n
ear the northwestern corner, while in the level model the coldest water is
near the northeastern corner. The isopycnal model appears to be a more real
istic representation of the real ocean since both wind and the thermohaline
circulation are thought to contribute to the North Atlantic subarctic cycl
onic gyre.
Although the zonally averaged buoyancy flux produced by the two model codes
is the same. the actual patterns of buoyancy flux at the surface are not s
imilar at high latitudes. This suggests that the two types of numerical mod
els would indicate very different air-sea interaction if coupled to atmosph
eric models and used to simulate climate. The application of the Gent-McWil
liams parameterization of mesoscale eddies to the model with z coordinates
and convective adjustment reduces the differences between the surface circu
lation of the two models by a small amount.