Formulation and implementation of a ''residual-mean'' ocean circulation model

A parameterization of mesoscale eddies in coarse-resolution ocean general circulation models (GCM) is formulated and implemented using a residual-mean formalism. In that framework, mean buoyancy is advected by the residual velocity (the sum of the Eulerian and eddy-induced velocities) and modified by a residual flux which accounts for the diabatic effects of mesoscale eddies. The residual velocity is obtained by stepping forward a residual-mean momentum equation in which eddy stresses appear as forcing terms. Study of the spatial distribution of eddy stresses, derived by using them as control parameters to ‘‘fit’’ the residual-mean model to observations, supports the idea that eddy stresses can be likened to a vertical down-gradient flux of momentum with a coefficient which is constant in the vertical. The residual eddy flux is set to zero in the ocean interior, where mesoscale eddies are assumed to be quasi-adiabatic, but is parameterized by a horizontal down-gradient diffusivity near the surface where eddies develop a diabatic component as they stir properties horizontally across steep isopycnals. The residual-mean model is implemented and tested in the MIT general circulation model. It is shown that the resulting model (1) has a climatology that is superior to that obtained using the Gent and McWilliams parameterization scheme with a spatially uniform diffusivity and (2) allows one to significantly reduce the (spurious) horizontal viscosity used in coarseresolution GCMs. � 2005 Elsevier Ltd. All rights reserved.

[1]  C. Reddi.  MIDDLE ATMOSPHERE DYNAMICS , 1998 .

[2]  P. Gent,et al.  Parameterizing eddy-induced tracer transports in ocean circulation models , 1995 .

[3]  L. Perelman,et al.  Hydrostatic, quasi‐hydrostatic, and nonhydrostatic ocean modeling , 1997 .

[4]  M. Redi.  Oceanic Isopycnal Mixing by Coordinate Rotation , 1982 .

[5]  L. Perelman,et al.  A finite-volume, incompressible Navier Stokes model for studies of the ocean on parallel computers , 1997 .

[6]  P. Gent,et al.  Isopycnal mixing in ocean circulation models , 1990 .

[7]  K. Trenberth,et al.  The mean annual cycle in global ocean wind stress , 1990 .

[8]  H. Hasumi,et al.  Developments in ocean climate modelling , 2000 .

[9]  G. Danabasoglu,et al.  Eulerian and Eddy-Induced Meridional Overturning Circulations in the Tropics , 2002 .

[10]  J. Marshall,et al.  Transformed Eulerian-Mean Theory. Part II: Potential Vorticity Homogenization and the Equilibrium of a Wind- and Buoyancy-Driven Zonal Flow , 2005 .

[11]  T. Keffer.  The Ventilation of the World's Oceans: Maps of the Potential vorticity Field , 1985 .

[12]  J. Green,et al.  Transfer properties of the large‐scale eddies and the general circulation of the atmosphere , 1970 .

[13]  D. Roemmich,et al.  Is the North Pacific in Sverdrup balance along 24°N? , 1994 .

[14]  William R. Young,et al.  Homogenization of potential vorticity in planetary gyres , 1982, Journal of Fluid Mechanics.

[15]  A. E. Gill,et al.  Energy partition in the large-scale ocean circulation and the production of mid-ocean eddies , 1974 .

[16]  James C. McWilliams,et al.  Sensitivity to Surface Forcing and Boundary Layer Mixing in a Global Ocean Model: Annual-Mean Climatology , 1997 .

[17]  Carl Wunsch,et al.  Global ocean circulation during 1992-1997, estimated from ocean observations and a general circulation model , 2002 .

[18]  John Marshall,et al.  Estimates and Implications of Surface Eddy Diffusivity in the Southern Ocean Derived from Tracer Transport , 2006 .

[19]  R. Greatbatch,et al.  On Parameterizing Vertical Mixing of Momentum in Non-eddy Resolving Ocean Models , 1990 .

[20]  J. Marshall,et al.  Equilibration of a warm pumped lens on a β plane , 2003 .

[21]  J. Marshall,et al.  Residual-Mean Solutions for the Antarctic Circumpolar Current and Its Associated Overturning Circulation , 2003 .

[22]  J. Marshall,et al.  Representation of Eddies in Primitive Equation Models by a PV Flux , 2000 .

[23]  Isaac M. Held,et al.  Parameterization of Quasigeostrophic Eddies in Primitive Equation Ocean Models. , 1997 .

[24]  R. A. Plumb.  Eddy Fluxes of Conserved Quantities by Small-Amplitude Waves , 1979 .

[25]  A. Adcroft,et al.  Representation of Topography by Shaved Cells in a Height Coordinate Ocean Model , 1997 .

[26]  Rüdiger Gerdes,et al.  The influence of numerical advection schemes on the results of ocean general circulation models , 1991 .

[27]  Carl Wunsch,et al.  Improved estimates of global ocean circulation, heat transport and mixing from hydrographic data , 2000, Nature.

[28]  John Marshall,et al.  On the Parameterization of Geostrophic Eddies in the Ocean , 1981 .

[29]  D. G. Andrews,et al.  Planetary Waves in Horizontal and Vertical Shear: The Generalized Eliassen-Palm Relation and the Mean Zonal Acceleration , 1976 .

[30]  Patrick Heimbach,et al.  Estimating Eddy Stresses by Fitting Dynamics to Observations Using a Residual-Mean Ocean Circulation Model and Its Adjoint , 2005 .

[31]  Gokhan Danabasoglu,et al.  Sensitivity of the global ocean circulation to parameterizations of mesoscale tracer transports , 1995 .

[32]  Dean Roemmich,et al.  Eddy Transport of Heat and Thermocline Waters in the North Pacific: A Key to Interannual/Decadal Climate Variability? , 2001 .

[33]  Peter H. Stone,et al.  An assessment of the Geophysical Fluid Dynamics Laboratory ocean model with coarse resolution: Annual‐mean climatology , 1999 .

[34]  Richard J. Greatbatch,et al.  Exploring the Relationship between Eddy-Induced Transport Velocity, Vertical Momentum Transfer, and the Isopycnal Flux of Potential Vorticity , 1998 .

[35]  K. Bryan,et al.  A Numerical Investigation of a Nonlinear Model of a Wind-Driven Ocean , 1963 .