Using a resolution function to regulate parameterizations of oceanic mesoscale eddy effects

Abstract Mesoscale eddies play a substantial role in the dynamics of the ocean, but the dominant length-scale of these eddies varies greatly with latitude, stratification and ocean depth. Global numerical ocean models with spatial resolutions ranging from 1° down to just a few kilometers include both regions where the dominant eddy scales are well resolved and regions where the model’s resolution is too coarse for the eddies to form, and hence eddy effects need to be parameterized. However, common parameterizations of eddy effects via a Laplacian diffusion of the height of isopycnal surfaces (a Gent–McWilliams diffusivity) are much more effective at suppressing resolved eddies than in replicating their effects. A variant of the Phillips model of baroclinic instability illustrates how eddy effects might be represented in ocean models. The ratio of the first baroclinic deformation radius to the horizontal grid spacing indicates where an ocean model could explicitly simulate eddy effects; a function of this ratio can be used to specify where eddy effects are parameterized and where they are explicitly modeled. One viable approach is to abruptly disable all the eddy parameterizations once the deformation radius is adequately resolved; at the discontinuity where the parameterization is disabled, isopycnal heights are locally flattened on the one side while eddies grow rapidly off of the enhanced slopes on the other side, such that the total parameterized and eddy fluxes vary continuously at the discontinuity in the diffusivity. This approach should work well with various specifications for the magnitude of the eddy diffusivities.

[1]  Alistair Adcroft,et al.  Simulations of underwater plumes of dissolved oil in the Gulf of Mexico , 2010 .

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

[3]  Mehmet Ilicak,et al.  Spurious dianeutral mixing and the role of momentum closure , 2012 .

[4]  Peter R. Gent,et al.  Response to Increasing Southern Hemisphere Winds in CCSM4 , 2011 .

[5]  A. E. Gill.  Atmosphere-Ocean Dynamics , 1982 .

[6]  Robert Hallberg,et al.  The Role of Eddies in Determining the Structure and Response of the Wind-Driven Southern Hemisphere Overturning: Results from the Modeling Eddies in the Southern Ocean (MESO) Project , 2006 .

[7]  M. Spall,et al.  Specification of eddy transfer coefficients in coarse resolution ocean circulation models , 1997 .

[8]  R. Hallberg,et al.  Reconciling estimates of the free surface height in Lagrangian vertical coordinate ocean models with mode-split time stepping , 2009 .

[9]  G. Flierl,et al.  Baroclinically Unstable Geostrophic Turbulence in the Limits of Strong and Weak Bottom Ekman Friction: Application to Midocean Eddies , 2004 .

[10]  Isaac M. Held,et al.  A Scaling Theory for Horizontally Homogeneous, Baroclinically Unstable Flow on a Beta Plane , 1996 .

[11]  David P. Marshall,et al.  Do We Require Adiabatic Dissipation Schemes in Eddy-Resolving Ocean Models? , 1998 .

[12]  D. Stammer.  Global Characteristics of Ocean Variability Estimated from Regional TOPEX/POSEIDON Altimeter Measurements , 1997 .

[13]  S. Jayne,et al.  Eddy–Mean Flow Interaction in the Kuroshio Extension Region , 2011 .

[14]  A. Adcroft,et al.  Parameterization of ocean eddies: Potential vorticity mixing, energetics and Arnold’s first stability theorem , 2010 .

[15]  R. Hallberg,et al.  A Parameterization of Shear-Driven Turbulence for Ocean Climate Models , 2008 .

[16]  Stephen M. Griffies,et al.  Biharmonic Friction with a Smagorinsky-Like Viscosity for Use in Large-Scale Eddy-Permitting Ocean Models , 2000 .

[17]  T. Schneider,et al.  The Surface Branch of the Zonally Averaged Mass Transport Circulation in the Troposphere , 1999 .

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

[19]  D. Chelton,et al.  Geographical Variability of the First Baroclinic Rossby Radius of Deformation , 1998 .

[20]  Stephen M. Griffies,et al.  Spurious Diapycnal Mixing Associated with Advection in a z-Coordinate Ocean Model , 2000 .

[21]  J. Pedlosky,et al.  A model of multiple zonal jets in the oceans : dynamical and kinematical analysis , 2009 .

[22]  Ronald,et al.  GFDL’s ESM2 Global Coupled Climate–Carbon Earth System Models. Part I: Physical Formulation and Baseline Simulation Characteristics , 2012 .

[23]  G. Danabasoglu,et al.  Effects of vertical variations of thickness diffusivity in an ocean general circulation model , 2007 .

[24]  Stephen M. Griffies,et al.  A boundary-value problem for the parameterized mesoscale eddy transport , 2010 .

[25]  K. S. Smith.  The geography of linear baroclinic instability in Earth's oceans , 2007 .

[26]  G. Vallis.  Atmospheric and Oceanic Fluid Dynamics: Fundamentals and Large-Scale Circulation , 2017 .

[27]  J. Pedlosky.  Geophysical Fluid Dynamics , 1979 .

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

[29]  P. Cessi.  An Energy-Constrained Parameterization of Eddy Buoyancy Flux , 2008 .

[30]  Carsten Eden,et al.  Effects of different closures for thickness diffusivity , 2009 .

[31]  C. Eden,et al.  Towards a mesoscale eddy closure , 2008 .

[32]  Norman A. Phillips,et al.  Energy Transformations and Meridional Circulations associated with simple Baroclinic Waves in a two-level, Quasi-geostrophic Model , 1954 .

[33]  Eric P. Chassignet,et al.  Gulf Stream Separation in Numerical Ocean Models , 2013 .

[34]  B. Fox‐Kemper,et al.  Parameterization of Mixed Layer Eddies. Part I. Theory and Diagnosis , 2008 .