If model parameterizations of unresolved physics, such as the variety of upper ocean mixing processes, are to hold over the large range of time and space scales of importance to climate, they must be strongly physically based. Observations, theories, and models of oceanic vertical mixing are surveyed. Two distinct regimes are identified: ocean mixing in the boundary layer near the surface under a variety of surface forcing conditions (stabilizing, destabilizing, and wind driven), and mixing in the ocean interior due to internal waves, shear instability, and double diffusion (arising from the different molecular diffusion rates of heat and salt). Mixing schemes commonly applied to the upper ocean are shown not to contain some potentially important boundary layer physics. Therefore a new parameterization of oceanic boundary layer mixing is developed to accommodate some of this physics. It includes a scheme for determining the boundary layer depth h, where the turbulent contribution to the vertical shear of a bulk Richardson number is parameterized. Expressions for diffusivity and nonlocal transport throughout the boundary layer are given. The diffusivity is formulated to agree with similarity theory of turbulence in the surface layer and is subject to the conditions that both it and its vertical gradient match the interior values at h. This nonlocal “K profile parameterization” (KPP) is then verified and compared to alternatives, including its atmospheric counterparts. Its most important feature is shown to be the capability of the boundary layer to penetrate well into a stable thermocline in both convective and wind-driven situations. The diffusivities of the aforementioned three interior mixing processes are modeled as constants, functions of a gradient Richardson number (a measure of the relative importance of stratification to destabilizing shear), and functions of the double-diffusion density ratio, Rρ. Oceanic simulations of convective penetration, wind deepening, and diurnal cycling are used to determine appropriate values for various model parameters as weak functions of vertical resolution. Annual cycle simulations at ocean weather station Papa for 1961 and 1969–1974 are used to test the complete suite of parameterizations. Model and observed temperatures at all depths are shown to agree very well into September, after which systematic advective cooling in the ocean produces expected differences. It is argued that this cooling and a steady salt advection into the model are needed to balance the net annual surface heating and freshwater input. With these advections, good multiyear simulations of temperature and salinity can be achieved. These results and KPP simulations of the diurnal cycle at the Long-Term Upper Ocean Study (LOTUS) site are compared with the results of other models. It is demonstrated that the KPP model exchanges properties between the mixed layer and thermocline in a manner consistent with observations, and at least as well or better than alternatives.

Oceanic vertical mixing: a review and a model with a nonlocal boundary layer parameterization

We examine observations of turbulence in the geophysical environment, primarily from oceans but also from lakes, in light of theory and experimental studies undertaken in the laboratory and with numerical simulation. Our focus is on turbulence in density-stratified environments and on the irreversible fluxes of tracers that actively contribute to the density field. Our understanding to date has come from focusing on physical problems characterized by high Reynolds number flows with no spatial or temporal variability, and we examine the applicability of these results to the natural or geophysical-scale problems. We conclude that our sampling and interpretation of the results remain a first-order issue, and despite decades of ship-based observations we do not begin to approach a reliable sampling of the overall turbulent structure of the ocean interior.

Density Stratification, Turbulence, but How Much Mixing?

Second-order turbulence closure models for geophysical boundary layers. A review of recent work

Laboratory experiments on stably stratified grid turbulence have suggested that turbulent diffusivity , suggesting that such Reynolds–Richardson number or Reynolds–Froude number aggregates are more descriptive of stratified turbulent flow conditions than the conventional reliance on Richardson number alone.

Parameterization of turbulent fluxes and scales using homogeneous sheared stably stratified turbulence simulations

Direct numerical simulations of stably and strongly stratified turbulent flows with Reynolds number Re >> 1 and horizontal Froude number F-h > 1, viscous forces are unimportant and l(v) scales as l ...

/pdf/scaling-analysis-and-simulation-of-strongly-stratified-2rq6ootada.pdf

Scaling analysis and simulation of strongly stratified turbulent flows

Homogeneous sheared stratified turbulence was simulated using a DNS code. The initial turbulent Reynolds numbers (Re) were 22, 44, and 89, and the initial dimensionless shear rate (S*) varied from 2 to 16. We found (similarly to Rogers (1986) for unstratified flows) the final value of S* at high Re to be ∼ 11, independent of initial S*. The final S* varies at low Re, in agreement with Jacobitz et al. (1997). At low Re, the stationary Richardson number (Ris) depends on both Re and S*, but at higher Re, it varies only with Re. A scaling based on the turbulent kinetic energy equation which suggests this result employs instantaneous rather than initial values of flow parameters.At high Re the dissipation increases with applied shear, allowing a constant final S*. The increased dissipation occurs primarily at high wavenumbers due to the stretching of eddies by stronger shear. For the high-Re stationary flows, the turbulent Froude number (Frt) is a constant independent of S*. An Frt-based scaling predicts the final value of S* well over a range of Re. Therefore Frt is a more appropriate parameter for describing the state of developed stratified turbulence than the gradient Richardson number.

/pdf/scaling-and-parameterization-of-stratified-homogeneous-3zakzegqw5.pdf

Lucinda H. Shih

Papers

Parameterization of turbulent fluxes and scales using homogeneous sheared stably stratified turbulence simulations

Scaling and parameterization of stratified homogeneous turbulent shear flow