Numerical simulations of the Mediterranean sea outflow: impact of the entrainment parameterization in an isopycnic coordinate ocean model

Gravity current entrainment is essential in determining the properties of the interior ocean water masses that result from marginal sea overflows. Although the individual entraining billows will be unresolvable in large-scale ocean models for the foreseeable future, some large-scale simulations are now being carried out that do resolve the intermediate scale environment which may control the rate of entrainment. Hallberg Mon. Wea. Rev. 128 (2000) 1402] has recently developed an implicit diapycnal mixing scheme for isopycnic coordinate ocean models that includes the Richardson number dependent entrainment parameterization of Turner J. Fluid Mech. 173 (1986) 431], and which may be capable of representing the gravity current evolution in large-scale ocean models. The present work uses realistic regional simulations with the Miami Isopycnic Coordinate Ocean Model (MICOM) to evaluate ability of this scheme to simulate the entrainment that is observed to occur in the bottom boundary currents downstream of the Mediterranean outflow. These simulations are strikingly similar to the observations, indicating that this scheme does produce realistic mixing between the Mediterranean outflow and the North Atlantic Central Water. Sensitivity studies identify the critical Richardson number below which vigorous entrainment occurs as a particularly important parameter. Some of these experiments also show meddies detaching from the Mediterranean undercurrent at locations that appear to be highly influenced by topographic features.