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

▪ AbstractThe coexistence in the deep ocean of a finite, stable stratification, a strong meridional overturning circulation, and mesoscale eddies raises complex questions concerning the circulation energetics. In particular, small-scale mixing processes are necessary to resupply the potential energy removed in the interior by the overturning and eddy-generating process. A number of lines of evidence, none complete, suggest that the oceanic general circulation, far from being a heat engine, is almost wholly governed by the forcing of the wind field and secondarily by deep water tides. In detail however, the budget of mechanical energy input into the ocean is poorly constrained. The now inescapable conclusion that over most of the ocean significant “vertical” mixing is confined to topographically complex boundary areas implies a potentially radically different interior circulation than is possible with uniform mixing. Whether ocean circulation models, either simple box or full numerical ones, neither explic...

/pdf/vertical-mixing-energy-and-the-general-circulation-of-the-530xn4nnsb.pdf

Vertical mixing, energy, and the general circulation of the oceans

A simple eddy kinetic energy parameterization of the oceanic vertical mixing is presented. The parameterization scheme is based on recent works on atmospheric turbulence modeling. It is designed to simulate vertical mixing at all depths, from the upper boundary layer down to the abyss. This scheme includes a single prognostic equation for the turbulent kinetic energy. The computation of the turbulent length scales is diagnostic, rather than prognostic. In weakly turbulent regions the simulated vertical diffusivity is inversely proportional to the Brunt-Vaisala frequency. In the first validation experiments presented here, the vertical mixing scheme is embedded into a simple one-dimensional model and used for upper ocean simulations at two very different test sites: the station Papa in the Gulf of Alaska and the Long-Term Upper Ocean Study (LOTUS) mooring in the Sargasso Sea. At station Papa the model successfully simulates the seasonal evolution of the upper ocean temperature field. At LOTUS the focus is on a short 2-week period. A detailed analysis of the oceanic heat budget during that period reveals a large bias in the bulk-derived surface heat fluxes. After correction of the fluxes the model does well in simulating the evolution of the temperature and wind-driven current. In particular, the large observed diurnal cycles of the sea surface temperature are well reproduced. During the second (windy) week of the selected period the model accounts for about two thirds of the kinetic energy of the observed upper ocean currents at periods larger than 6 hours. The local wind forcing thus appears to be the dominant generation mechanism for the near-inertial motions, which are the most energetic. The velocity simulation is especially good at the low frequencies. During the second simulated week the model accounts for as much as 78% of the kinetic energy at subinertial frequencies. The simulated mean velocity profile is reminiscent of an Ekman spiral, in agreement with the observations.

A simple eddy kinetic energy model for simulations of the oceanic vertical mixing: Tests at Station Papa and long-term upper ocean study site

An improved mixed layer model, based on second-moment closure of turbulence and suitable for application to oceanic and atmospheric mixed layers, is described. The model is tested against observational data from different locations in the global oceans, including high latitudes and tropics. The model belongs to the Mellor-Yamada hierarchy but incorporates recent findings from research on large eddy simulations and second-moment closure. The modified expansion of Galperin, Kantha, Hassid and Rosati (1988) that leads to a much simpler and more robust quasi-equilibrium turbulence model is employed instead of the original Mellor and Yamada (1974) model. Findings from ongoing research at the National Center for Atmospheric Research on large eddy simulations of the atmospheric boundary layer are utilized to improve parameterizations of pressure covariance terms in the second-moment closure. Shortwave solar radiation penetration is given careful treatment in the model so that the model is applicable to investigations of biological and photochemical processes in the upper ocean. But by far the major improvement is in the inclusion of the shear instability-induced mixing in the strongly stratified region below the oceanic mixed layer that leads to a more realistic and reliable mixed layer model that is suitable for application to a variety of geophysical mixed layers and circulation problems. The model appears to predict the mixing in the upper ocean well on a variety of time scales, from event scale storm-induced deepening and diurnal scale variability to seasonal time scales. With proper attention to the heat and salt balances in the upper ocean, it should be possible to use it for simulations of interannual variability as well. While the model validation has been primarily against oceanic mixed layer data sets, it is believed that the improvements can be readily incorporated into a model of the atmospheric boundary layer as well.

An improved mixed layer model for geophysical applications

This book is written for college juniors and seniors and new graduate students in meteorology, ocean engineering, and oceanography. It begins with a brief overview of what is known about the ocean. This is followed by a description of the ocean basins, for the shape of the seas influences the physical processes in the water. Next, students will study the external forces, wind and heat, acting on the ocean, and the ocean's response. It also includes the equations describing dynamic response of the ocean. For example, the equations of motion, the influence of earth's rotation, and viscosity. Finally, students consider some particular examples: the deep circulation, the equatorial ocean and El NiEœno, and the circulation of particular areas of the ocean. Contents: 1) A Voyage of Discovery. 2) The Historical Setting. 3) The Physical Setting. 4) Atmospheric Influences. 5) The Oceanic Heat Budget. 6) Temperature, Salinity and Density. 7) The Equations of Motion. 8) Equations of Motion with Viscosity. 9) Response of the Upper Ocean to Winds. 10) Geostrophic Currents. 11) Wind Driven Ocean Circulation. 12) Vorticity in the Ocean. 13) Deep Circulation in the Ocean. 14) Equatorial Processes. 15) Numerical Models. 16) Ocean Waves. 17) Coastal Processes and Tides.

Introduction to Physical Oceanography

Diapycnal fluxes of momentum and heat produced by three-dimensional turbulence play important, but poorly understood, roles in the dynamics of the main thermocline and of the equatorial undercurrent. Diverse approaches—involving inferences, measurements, and process studies—are being pursued to define these fluxes and their significance. These are reviewed to improve the coordination of future mixing studies.

https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/JC092iC05p05249%4010.1002/%28ISSN%292169-9291.MSTRATF1

M. C. Gregg

Papers

Diapycnal mixing in the thermocline: A review