Dynamic Model of Mesoscale Eddies

Abstract

Oceanic mesoscale eddies which are analogs of well known synoptic eddies (cyclones and anticyclones), are studied on the basis of the turbulence model originated by Dubovikov (Dubovikov, M.S., "Dynamical model of turbulent eddies", Int. J. Mod. Phys. B7, 4631-4645 (1993).) and further developed by Canuto and Dubovikov (Canuto, V.M. and Dubovikov, M.S., "A dynamical model for turbulence: I. General formalism", Phys. Fluids 8, 571-586 (1996a) (CD96a); Canuto, V.M. and Dubovikov, M.S., "A dynamical model for turbulence: II. Sheardriven flows", Phys. Fluids 8, 587-598 (1996b) (CD96b); Canuto, V.M., Dubovikov, M.S., Cheng, Y. and Dienstfrey, A., "A dynamical model for turbulence: III. Numerical results", Phys. Fluids 8, 599-613 (1996c)(CD96c); Canuto, V.M., Dubovikov, M.S. and Dienstfrey, A., "A dynamical model for turbulence: IV. Buoyancy-driven flows", Phys. Fluids 9, 2118-2131 (1997a) (CD97a); Canuto, V.M. and Dubovikov, M.S., "A dynamical model for turbulence: V. The effect of rotation", Phys. Fluids 9, 2132-2140 (1997b) (CD97b); Canuto, V.M., Dubovikov, M.S. and Wielaard, D.J., "A dynamical model for turbulence: VI. Two dimensional turbulence", Phys. Fluids 9, 2141-2147 (1997c) (CD97c); Canuto, V.M. and Dubovikov, M.S., "Physical regimes and dimensional structure of rotating turbulence", Phys. Rev. Lett. 78, 666-669 (1997d) (CD97d); Canuto, V.M., Dubovikov, M.S. and Dienstfrey, A., "Turbulent convection in a spectral model", Phys. Rev. Lett. 78, 662-665 (1997e) (CD97e); Canuto, V.M. and Dubovikov, M.S., "A new approach to turbulence", Int. J. Mod. Phys. 12, 3121-3152 (1997f) (CD97f); Canuto, V.M. and Dubovikov, M.S., "Two scaling regimes for rotating Raleigh-Benard convection", Phys. Rev. Letters 78, 281-284, (1998) (CD98); Canuto, V.M. and Dubovikov, M.S., "A dynamical model for turbulence: VII. The five invariants for shear driven flows", Phys. Fluids 11, 659-664 (1999a) (CD99a); Canuto, V.M., Dubovikov, M.S. and Yu, G., "A dynamical model for turbulence: VIII. IR and UV Reynolds stress spectra for shear driven flows", Phys. Fluids 11, 656-677 (1999b) (CD99b); Canuto, V.M., Dubovikov, M.S. and Yu, G., "A dynamical model for turbulence: IX. The Reynolds stress for shear driven flows", Phys. Fluids 11, 678-694 (1999c) (CD99c).). The CD model derives from general principles and does not resort to any free parameters. Yet, it successfully describes a wide variety of quite different turbulent flows. In the present work we apply CD model to the compressible ocean. The model yields mesoscale eddies generated by the baroclinic instability. The latter, in turn, arises from the nonhorizontal orientation of the surfaces of the constant potential density (isopycnals). The obtained dynamic equations for eddy fields reduce to a vertical eigen value problem, an eigen value real part yielding an eddy radius, while an imaginary part - an eddy drift velocity. The size of the eddy is about 3rd (where rd is the Rossby deformation radius). The eddy dynamics has the following distinctive features: (1) the large scale potential energy feeds the eddy potential energy (EPE) at scales ∼ rd , (2) from rd EPE cascades to the smaller scales down to ∼ l 1 determined from the condition that the spectral Rossby number Ro(q) ≡ qU'(q)f−1 becomes ∼ 1 (q is two-dimensional wave number within an isopycnal surface), (3) at scales ∼ l 1 EPE transforms into eddy kinetic energy (EKE) which cascades backwards to the larger scales up to ∼ rd , where it transforms back into EPE, thereby closing the energy flux circulation in a wavenumber space, (4) dissipation of the eddy energy (EE) occurs at scales ∼ l 1 since at those scales the fluctuating component of the vertical shear is maximal and equals to the Brunt-Vaisala frequency. The latter equality is the well known condition for generating the vertical turbulence which dissipates EE. The model enables to determine all turbulence characteristics, including the horizontal (isopycnal) diffusivity κ h in terms of the large scale mean fields. From the typical values of the latter follow estimates for the parameters of an eddy which agree well with the observational and simulational data: kh ∼ 103m2s−1, EKE K ∼ 103m2s−1, rd ∼ 3 × 104m, lI ∼ 10. In what concerns the bolus velocity, it contains additional terms (as compared to the model of Gent and McWilliams (Gent, P.R. and McWilliams, J.C., "Isopycnal mixing in ocean circulation models", J. Phys. Oceanogr. 20, 150-155 (1990)) which result from the eddy fields advection by a mean velocity ū. Since the latter varies with depth, it is inevitable to differ from the eddy drift velocity that produces a shearing force eroding the eddy coherent structures and, therefore, contributing negatively to EE production. This is in contrast with the positive contribution from the GM term (which is due to the baroclinic instability). In those regions where the disruptive action is stronger, there is no eddy generation.