A diapycnal diffusion algorithm for isopycnal ocean circulation models with special application to mixed layers

Abstract

An algorithm is presented for solving the one-dimensional diffusion equation for density, written in terms of density (or a like surrogate) as the independent variable. The algorithm maintains nonnegative layer thicknesses, the premise of the transformation to density as the independent coordinate, under certain restrictions. Near-zero thickness layers can be maintained at the boundaries to accommodate future inflation in response to heating from the boundary. Layers can shrink to near-zero thickness in response to cooling from the boundary. A slight modification of the algorithm permits layers to have diffusion coefficients which differ by orders of magnitude. This provides a natural framework for a surface mixed layer in an isopycnal model, in which the mixed layer is distinguished as a zone of very high turbulent diffusivity overlying an ocean interior of much smaller turbulent diffusivity. The “mixed layer” may be an aggregation of several isopycnal layers rather than just one. A substantial jump in density at the mixed layer base can be represented by several near-zero thickness isopycnal layers. The specification of the thickness of the mixing zone, i.e., the mixed layer depth, is external to the algorithm. An illustration is given using a Kraus–Turner-type specification.