Ocean Turbulence. Part II: Vertical Diffusivities of Momentum, Heat, Salt, Mass, and Passive Scalars
V. M. Canuto,a m howard,Y. Cheng,M. S. Dubovikov
DOI: https://doi.org/10.1175/1520-0485(2002)032<0240:OTPIVD>2.0.CO;2
2010-01-01
Journal of Physical Oceanography
Abstract:A Reynolds stress-based model is used to derive algebraic expressions for the vertical diffusivities K-alpha (alpha = m, h, s) for momentum, heat, and salt. The diffusivities are expressed as K-alpha(R-rho, N, Ri(T), epsilon) in terms of the density ratio R-rho = alpha(s)partial derivativeS/partial derivativez(alpha(T)partial derivativeT/partial derivativez)(-1), the Brunt-Vaisala frequency N-2 = -grho(o)(-1)partial derivativerho/partial derivativez, the Richardson number Ri(T) = N-2/Sigma(2) (Sigma is the shear), and the dissipation rate of kinetic energy epsilon. The model is valid both in the mixed layer (ML) and below it. Here R-rho and N are computed everywhere using the large-scale fields from an ocean general circulation model while Ri(T) is contributed by resolved and unresolved shear. In the ML, the wind-generated large-scale shear dominates and can be computed within an OGCM. Below the ML, the wind is no longer felt and small-scale shear dominates. In this region, the model provides a new relation Ri(T) = cf (R-rho) with c approximate to 1 in lieu of Munk's suggestion Ri(T) approximate to c. Thus, below the ML, the K-alpha become functions of R-rho, N, and epsilon. The dissipation epsilon representing the physical processes responsible for the mixing, which are different in different parts of the ocean, must also be expressed in terms of the large-scale fields. In the ML, the main source of stirring is the wind but below the ML there is more than one possible source of stirring. For regions away from topography, one can compute epsilon using a model for internal waves. On the other hand, near topography, one must employ different expressions for epsilon. In agreement with the data, the resulting diffusivities are location dependent rather than universal values. Using North Atlantic Tracer Release Experiment (NATRE) data, the authors test the new diffusivities with and without an OGCM. The measured diffusivities are well reproduced. Also, a set of global T and S profiles is computed using this model and the KPP model. The profiles are compared with Levitus data. In the North Atlantic, at 24degreesN, the meridional overturning is close to the measured values of 17 +/- 4 Sv and 16 +/- 5 Sv (Sv = 10(6) m(3) s(-1)). The polar heat transport for the North Atlantic Ocean, the Indo-Pacific Ocean, and the global ocean are generally lowered by double diffusion. The freshwater budget is computed and compared with available data.