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

A subgrid-scale form for mesoscale eddy mixing on isopycnal surfaces is proposed for use in non-eddy-resolving ocean circulation models. The mixing is applied in isopycnal coordinates to isopycnal layer thickness, or inverse density gradient, as well as to passive scalars, temperature and salinity. The transformation of these mixing forms to physical coordinates is also presented.

Isopycnal mixing in ocean circulation models

Communities of microscopic plant life, or phytoplankton, dominate the Earth's aquatic ecosystems. This important new book by Colin Reynolds covers the adaptations, physiology and population dynamics of phytoplankton communities in lakes and rivers and oceans. It provides basic information on composition, morphology and physiology of the main phyletic groups represented in marine and freshwater systems and in addition reviews recent advances in community ecology, developing an appreciation of assembly processes, co-existence and competition, disturbance and diversity. Although focussed on one group of organisms, the book develops many concepts relevant to ecology in the broadest sense, and as such will appeal to graduate students and researchers in ecology, limnology and oceanography.

/pdf/the-ecology-of-phytoplankton-4ev70n3jsq.pdf

The Ecology of Phytoplankton

We present the global general circulation model IPSL-CM5 developed to study the long-term response of the climate system to natural and anthropogenic forcings as part of the 5th Phase of the Coupled Model Intercomparison Project (CMIP5). This model includes an interactive carbon cycle, a representation of tropospheric and stratospheric chemistry, and a comprehensive representation of aerosols. As it represents the principal dynamical, physical, and bio-geochemical processes relevant to the climate system, it may be referred to as an Earth System Model. However, the IPSL-CM5 model may be used in a multitude of configurations associated with different boundary conditions and with a range of complexities in terms of processes and interactions. This paper presents an overview of the different model components and explains how they were coupled and used to simulate historical climate changes over the past 150 years and different scenarios of future climate change. A single version of the IPSL-CM5 model (IPSL-CM5A-LR) was used to provide climate projections associated with different socio-economic scenarios, including the different Representative Concentration Pathways considered by CMIP5 and several scenarios from the Special Report on Emission Scenarios considered by CMIP3. Results suggest that the magnitude of global warming projections primarily depends on the socio-economic scenario considered, that there is potential for an aggressive mitigation policy to limit global warming to about two degrees, and that the behavior of some components of the climate system such as the Arctic sea ice and the Atlantic Meridional Overturning Circulation may change drastically by the end of the twenty-first century in the case of a no climate policy scenario. Although the magnitude of regional temperature and precipitation changes depends fairly linearly on the magnitude of the projected global warming (and thus on the scenario considered), the geographical pattern of these changes is strikingly similar for the different scenarios. The representation of atmospheric physical processes in the model is shown to strongly influence the simulated climate variability and both the magnitude and pattern of the projected climate changes.

/pdf/climate-change-projections-using-the-ipsl-cm5-earth-system-2gm93qrdy6.pdf

Climate change projections using the IPSL-CM5 Earth System Model: From CMIP3 to CMIP5

▪ 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

The dependence of phytoplankton photosynthesis on light intensity may be altered by the range and frequency of variations in light intensity recentlv experienced by the organisms. A major source of the fluctuations in light intensity experienced by phytoplankton in the upper ocean is vertical motion. We estimate time and space scales for \\Tertical displacements of phytoplankton caused by turbulent mixing, internal waves, Langmuir circulations, and double diffusive processes. In the surface layer, depending on windspeed, current shear and stratification, we find that time scales for cycling of phytoplankton by turbulent eddies and mixing vary from about 0.5 h to hundreds of hours for vertical displacements of the order of 10 m. In the seasonal thermocline, turbulent diffusive time scales for displacements as small as several meters are weeks to months, whereas similar displacements by internal waves occur over periods of several minutes to several hours, according to the strength of the density stratification, and are then dominant. Langmuir cells seem to scale as the large turbulent eddies and need not be treated separately, and double diffusive processes seem to be of minor importance. The formulation used here of a vertical turbulent diffusion coefficient K, as a function of observable quantities-e the rate of dissipation of turbulent kinetic energy, and N the local buoyancy frequency-should also be us&d for estimating vertical fluxes of nutrients. In addition, this formulation is reversible in time and can be used to estimate the recent depth and light history of phytoplankton taken from the upper ocean. Traditionally the dependence of the photosynthetic production rate of phytoplankton on light intensity has been described by “P vs. I” curves relating the rate of photosynthetic production, P, to increasing light intensity, I, determined usually from experiments where portions of the same sample are incubated at several different constant light intensities for the same period. The uptake of labeled carbon over that period is measured and considered to be an index of the photosynthetic production (Peterson 1980). Division by the time interval gives the photosynthetic rate. A smooth curve fitted through the data resembles one of the family of curves shown in Fig. 1 (which, for illustration, has been generated from the model equation of Platt et al. 1980): P(I) = aem”‘( 1 eeVr). (1) Several problems exist: different curves are often found at different times of day (provided the incubation period is sufficiently brief, say 2-4 h), at different depths, during different seasons, and with different organisms (Harris 1973; MacCaull and Platt 1977; Platt et al. 1980). Furthermore, time-course measurements of varying photosynthetic production P(t)

https://aslopubs.onlinelibrary.wiley.com/doi/pdf/10.4319/lo.1983.28.5.0801

Time and space scales of vertical mixing and advection of phytoplankton in the upper ocean

Vertical eddy diffusivity in the ocean interior

The validity of the assumption of local isotropy is investigated using measurements of three orthogonal components of the turbulent velocity fields associated with initially high-Reynolds-number geophysical turbulence The turbulent fields, generated by various large-scale internal motions caused by tidal flows over an estuarine sill, decay under the influence of stable mean density gradients With measurements from sensors mounted on a submersible, we examine the evolution of spectral shapes and of ratios of cross-stream to streamwise components, as well as the degree of high-wavenumber universality, for the observational range of the parameter I≡ ks/kb = lb/ls This ratio is a measure of separation between the Kolmogoroff wavenumber ks≡ (e/ν3)¼ ≡ 2π/ls typical of scales by which turbulent kinetic energy has been dissipated (at rate e), and the buoyancy wavenumber kb ≡ (N3/e)½ ≡ 2π/lb typical of scales at which the ambient stratification parameter N ≡ (−gρz/ρ0)½ becomes important For values of I larger than ∼ 3000, inertial subranges are observed in all spectra, and the spectral ratio ϕ22/ϕ11 of cross-stream to streamwise spectral densities reaches the isotropic value of 4/3 for about a decade in wavenumber As ks/kb decreases, inertial subranges vanish, but spectra of the cross-stream and streamwise components continue to satisfy isotropic relationships at dissipation wavenumbers We provide a criterion for when e may safely be estimated from a single measured component of the dissipation tensor, and also explore questions of appropriate low-wavenumber normalization for buoyancy-modified turbulence

https://www.cambridge.org/core/services/aop-cambridge-core/content/view/S0022112084001592

Local isotropy and the decay of turbulence in a stratified fluid

The Romanche Fracture Zone is a major gap in the Mid-Atlantic Ridge at the equator, which is deep enough to allow significant eastward flows of Antarctic Bottom Water from the Brazil Basin to the Sierra Leone and Guinea Abyssal Plains. While flowing through the Romanche Fracture Zone, bottom-water properties are strongly modified due to intense vertical mixing. The diapycnal mixing coefficient in the bottom water of the Romanche Fracture Zone is estimated by using the finestructure of CTD profiles, the microstructure of high-resolution profiler data, and by constructing a heat budget from current meter data. The finestructure of density profiles is described using the Thorpe scalesLT. It is shown from microstructure data taken in the bottom water that the Ozmidov scale LO is related to LT by the linear relationship LO 5 0.95LT, similar to other studies, which allows an estimate of the diapycnal mixing coefficient using the Osborn relation. The Thorpe scale and the diapycnal mixing coefficient estimates show enhanced mixing downstream (eastward) of the main sill of the Romanche Fracture Zone. In this region, a mean diapycnal mixing coefficient of about 1000 3 1024 m2 s21 is found for the bottom water. Estimates of cross-isothermal mixing coefficient derived from the heat budgets constructed downstream of the current meter arrays deployed in the Romanche Fracture Zone and the nearby Chain Fracture Zone are in agreement with the finestructure estimates of the diapycnal mixing coefficient within the Romanche Fracture Zone. Although the two fracture zones occupy only 0.4% of the area covered by the Sierra Leone and Guinea Abyssal Plains, the diffusive heat fluxes across the 1.4 8C isotherm in the Romanche and Chain Fracture Zones are half that found over the abyssal plains across the 1.88C isotherm, emphasizing the role of these passages for bottom-water property modifications.

/pdf/mixing-in-the-romanche-fracture-zone-3ig7fj6oz9.pdf

Mixing in the Romanche Fracture Zone

While recent evidence suggests that North Pacific salmon stocks are influenced by decadal variability in atmospheric forcing of the ocean, the actual combination of physical and biological processes that determines this linkage has not been identified. This paper describes a possible scenario in which water column stability is the primary factor by which the physical environment influences phytoplankton production, the basis for production at higher trophic levels. Variation in the strength of the wintertime Aleutian Low pressure area affects water column stabilities, hence primary production, along the entire eastern boundary of the North Pacific. The `optimal stability window' explains the qualitative relationship between fish stocks and the strength of the winter Aleutian Low, as well as the observed out-of-phase variation between northern and southern salmon stocks.

Ann E. Gargett

Papers

Time and space scales of vertical mixing and advection of phytoplankton in the upper ocean

Vertical eddy diffusivity in the ocean interior

Local isotropy and the decay of turbulence in a stratified fluid

Mixing in the Romanche Fracture Zone

The optimal stability `window': a mechanism underlying decadal fluctuations in North Pacific salmon stocks?