Being aligned with the sloping seabed, the transverse flow transports less dense water down in the southern flank of the channel. Therefore the salinity/density contours bend downwards, Bcl-2 cleavage displaying a tendency to become vertical and eventually produce inverted, hydrostatically unstable stratification. However, when the density contours approach the vertical, the density stratification weakens and the stratified shear gravity current becomes hydrodynamically unstable, producing turbulent mixing together with vertical homogenization of BBL, thereby establishing a pure horizontal density gradient. This was demonstrated in the POM simulation (Figure 4), where the instability of the stratified shear current is plausibly

parameterized by the 21/2 moment turbulence closure (Mellor & Yamada 1982). The

parameterization explicitly describes the effect of stratification on vertical mixing, since the vertical turbulent viscosity KM   and heat/salt diffusivity KH   are expressed as equation(5) KM=lqSM(Rit),KH=lqSH(Rit),where q   is the root mean square velocity fluctuation (so that q  2 is the specific kinetic energy of turbulence), Z-VAD-FMK mouse l   is the external length scale of turbulence, and SM   and SH   are functions of the Richardson number Rit   equation(6) Rit=l2q2gρ0∂ρpot∂z,where ρpot   is the potential density and ρ  0 is the reference density. Note that Rit   < 0 when stratification is hydrostatically stable (in this case −(g/ρ0)(∂ρpot/∂z)≡N2−(g/ρ0)(∂ρpot/∂z)≡N2 is the squared buoyancy frequency), Rit = 0 for neutral stratification, and Rit > 0 for hydrostatically unstable stratification. For neutral stratification (Rit = 0) SM = 0.8 SH = 0.39 and for stable stratification SM and SH are infinitesimally Liothyronine Sodium small with |Rit| (i.e. SM ≈ SH → 0 at Rit → –∞, and, for example, SM ≈ SH = 0.014 at Rit = –1). And finally, for unstable stratification, SM and SH increase rapidly with the growth of an unstable

(inverted) potential density gradient, achieving in the POM code a practical limit of SM = 0.75 SH = 12.7 at Rit = 0.028 and further retaining the same limiting value at Rit > 0.028. Therefore, even when an inverted density gradient was formed as a result of differential transverse advection, the above described drastic increase of vertical eddy diffusivity/viscosity at unstable density stratification would mix up the inversion and establish vertical quasi-homogeneity, so that the residual inverted gradients would be strongly depressed. Unlike POM, the MIKE 3 simulation is based on the Smagorinsky subgrid scale model turbulent closure, which does not explicitly allow for stratification. The Smagorinsky subgrid diffusivity is simply taken to be proportional to the product of the squared vertical grid size and velocity gradients, implying that the model is able to resolve the instability of shear stratified flow and the related intensification of vertical mixing.