Eddy Mixing of Potential Vorticity versus Thickness in an Isopycnic Ocean Model

Sybren S. Drijfhout Royal Netherlands Meteorological Institute, De Bilt, Netherlands

Search for other papers by Sybren S. Drijfhout in
Current site
Google Scholar
PubMed
Close
and
Wilco Hazeleger Royal Netherlands Meteorological Institute, De Bilt, Netherlands

Search for other papers by Wilco Hazeleger in
Current site
Google Scholar
PubMed
Close
Restricted access

Abstract

Parameterizations of the eddy-induced velocity that advects tracers in addition to the Eulerian mean flow are traditionally expressed as a downgradient Fickian diffusion of either isopycnal layer thickness or large-scale potential vorticity (PV). There is an ongoing debate on which of the two closures is better and how the spatial dependence of the eddy diffusivity should look like. To increase the physical reasoning on which these closures are based, the authors present a systematic assessment of eddy fluxes of thickness and PV and their relation to mean-flow gradients in an isopycnic eddy-resolving model of an idealized double-gyre circulation in a flat bottom, closed basin. The simulated flow features strong nonlinearities, such as tight inertial recirculations, a meandering midlatitude jet, pools of homogenized PV, and regions of weak flow where β/h dominates the PV gradient. It is found that the zonally averaged eddy flux of thickness scales better with the zonally averaged meridional thickness gradient than the eddy flux of PV with the PV gradient. The reason for this is that the two-scale approximation, which is often invoked to derive a balance between the downgradient eddy flux of PV and enstrophy dissipation, does not hold. It is obscured by advection of perturbation enstrophy, which is multisigned and weakly related to mean-flow gradients. On the other hand, forcing by vertical motions, which enters the balance between the downgradient eddy flux of thickness and dissipation in most cases, acts to dissipate thickness variance. It is dominated by the conversion from potential to kinetic energy and the subsequent downgradient transport of thickness. Also, advection of perturbation thickness variance tends to be more single-signed than advection of perturbation enstrophy, forcing the eddy flux of thickness to be more often down the mean gradient. As a result, in the present configuration a downgradient diffusive closure for thickness seems more appropriate to simulate the divergent eddy fluxes than a downgradient diffusive closure for PV, especially in dynamically active regions where the eddy fluxes are large and in regions of nearly uniform PV.

Corresponding author address: Sybren S. Drijfhout, Royal Netherlands Meteorological Institute (KNMI), P.O. Box 201, 3730 AE, De Bilt, Netherlands.

Email: drijfhou@knmi.nl

Abstract

Parameterizations of the eddy-induced velocity that advects tracers in addition to the Eulerian mean flow are traditionally expressed as a downgradient Fickian diffusion of either isopycnal layer thickness or large-scale potential vorticity (PV). There is an ongoing debate on which of the two closures is better and how the spatial dependence of the eddy diffusivity should look like. To increase the physical reasoning on which these closures are based, the authors present a systematic assessment of eddy fluxes of thickness and PV and their relation to mean-flow gradients in an isopycnic eddy-resolving model of an idealized double-gyre circulation in a flat bottom, closed basin. The simulated flow features strong nonlinearities, such as tight inertial recirculations, a meandering midlatitude jet, pools of homogenized PV, and regions of weak flow where β/h dominates the PV gradient. It is found that the zonally averaged eddy flux of thickness scales better with the zonally averaged meridional thickness gradient than the eddy flux of PV with the PV gradient. The reason for this is that the two-scale approximation, which is often invoked to derive a balance between the downgradient eddy flux of PV and enstrophy dissipation, does not hold. It is obscured by advection of perturbation enstrophy, which is multisigned and weakly related to mean-flow gradients. On the other hand, forcing by vertical motions, which enters the balance between the downgradient eddy flux of thickness and dissipation in most cases, acts to dissipate thickness variance. It is dominated by the conversion from potential to kinetic energy and the subsequent downgradient transport of thickness. Also, advection of perturbation thickness variance tends to be more single-signed than advection of perturbation enstrophy, forcing the eddy flux of thickness to be more often down the mean gradient. As a result, in the present configuration a downgradient diffusive closure for thickness seems more appropriate to simulate the divergent eddy fluxes than a downgradient diffusive closure for PV, especially in dynamically active regions where the eddy fluxes are large and in regions of nearly uniform PV.

Corresponding author address: Sybren S. Drijfhout, Royal Netherlands Meteorological Institute (KNMI), P.O. Box 201, 3730 AE, De Bilt, Netherlands.

Email: drijfhou@knmi.nl

Save
  • Andrews, D. G., and M. E. McIntyre, 1976: Planetary waves in horizontal and vertical shear: The generalized Eliassen–Palm relation and the zonal mean acceleration. J. Atmos. Sci.,33, 2031–2048.

  • ———, T. R. Holton, and C. B. Leovy, 1987: Middle Atmosphere Dynamics. Academic Press, 489 pp.

  • Bleck, R., and D. B. Boudra, 1981: Initial testing of a numerical ocean circulation model using a hybrid (quasi-isopycnic) vertical coordinate. J. Phys. Oceanogr.,11, 755–770.

  • ——, and ——, 1986: Wind driven spin-up in eddy-resolving ocean models formulated in isopycnic and isobaric coordinates. J. Geophys. Res.,91, 7611–7621.

  • Bryan, K., J. K. Dukowicz, and R. D. Smith, 1999: On the mixing coefficient in the parameterization of bolus velocity. J. Phys. Oceanogr.,29, 2442–2456.

  • Danabasoglu, G., and J. C. McWilliams, 1995: Sensitivity of the global ocean circulation to parameterizations of mesoscale tracer transports. J. Climate,8, 2967–2987.

  • Drijfhout, S. S., 1994: On the heat transport by mesoscale eddies in an ocean circulation model. J. Phys. Oceanogr.,24, 429–442.

  • Duffy, P. B., P. Eltgroth, A. J. Bourgeois, and K. Caldeira, 1995: Effect of improved subgrid scale transport of tracers on uptake of bomb radiocarbon in the GFDL ocean general circulation model. Geophys. Res. Lett.,22, 1065–1068.

  • England, M. H., 1995: Using chlorofluorocarbons to assess ocean climate models. Geophys. Res. Lett.,22, 3051–3054.

  • Gent, P. R., and J. C. McWilliams, 1990: Isopycnic mixing in ocean circulation models. J. Phys. Oceanogr.,20, 150–155.

  • ——, J. Willebrand, T. J. McDougall, and J. C. McWilliams, 1995: Parameterizing eddy-induced tracer transports in ocean circulation models. J. Phys. Oceanogr.,25, 463–474.

  • Gille, S. T., and R. E. Davis, 1999: The influence of mesoscale eddies on coarsely resolved density: An examination of subgrid-scale parameterization. J. Phys. Oceanogr.,29, 1109–1123.

  • Greatbatch, R. J., 1998: Exploring the relationship between eddy-induced transport velocity, vertical momentum transfer, and the isopycnal flux of potential vorticity. J. Phys. Oceanogr.,28, 422–432.

  • Hirst, A. C., and T. J. McDougall, 1996: Deep-water properties and surface buoyancy flux as simulated by a z coordinate model including eddy-induced advection. J. Phys. Oceanogr.,26, 1320–1343.

  • Holland, W. R., and P. B. Rhines, 1980: An example of eddy-induced ocean circulation. J. Phys. Oceanogr.,10, 1010–1031.

  • ——, T. Keffer, and P. B. Rhines, 1984: Dynamics of the oceanic general circulation: The potential vorticity field. Nature,308, 698–705.

  • Jiang, S., F.-F. Jin, and M. Ghil, 1995: Multiple equilibria, periodic and aperiodic solutions in a wind-driven, double-gyre, shallow-water model. J. Phys. Oceanogr.,25, 764–786.

  • Killworth, P. D., 1997: On the parameterization of eddy transfer. Part I: Theory. J. Mar. Res.,55, 1171–1197.

  • ——, 1998: On the parameterization of eddy transfer. Part II: Tests with a channel model. J. Mar. Res.,56, 349–374.

  • Lee, M.-M., D. P. Marshall, and R. G. Williams, 1997: On the eddy transfer of tracers: Advective or diffusive? J. Mar. Res.,55, 483–505.

  • Marshall, D. P., R. G. Williams, and M.-M. Lee, 1999: The relation between eddy-induced transport and isopycnic gradients of potential vorticity. J. Phys. Oceanogr.,29, 1571–1578.

  • Marshall, J. C., and G. Shutts, 1981: A note on rotational and divergent eddy fluxes. J. Phys. Oceanogr.,11, 1677–1680.

  • McDougall, T. J., 1998: Three-dimensional residual-mean theory. Ocean Modeling and Parameterization, E. P. Chassignet and J. Verron, Eds., NATO ASI Series C, Kluwer Academic, 95–122.

  • ——, and P. C. McIntosh, 1996: The temporal-residual-mean velocity. Part I: Derivation and the scalar conservation equations. J. Phys. Oceanogr.,26, 2653–2665.

  • McDowell, S., P. B. Rhines, and T. Keffer, 1982: North Atlantic potential vorticity and its relation to the general circulation. J. Phys. Oceanogr.,12, 1417–1436.

  • McWilliams, J. C., W. R. Holland, and J. H. S. Chow, 1978: A description of numerical Antarctic Circumpolar Currents. Dyn. Atmos. Oceans,2, 213–291.

  • O’Dwyer, J., and R. G. Williams, 1997: The climatological distribution of potential vorticity over the abyssal ocean. J. Phys. Oceanogr.,27, 2488–2506.

  • Plumb, R. A., and J. D. Mahlman, 1987: The zonally-averaged transport characteristics of the GFDL General Circulation/Transport Model. J. Atmos. Sci.,44, 298–327.

  • Press, W. H., B. P. Flannery, S. A. Teukolsky, and W. T. Vetterling, 1986: Numerical Recipes: The Art of Scientific Computing. Cambridge University Press, 818 pp.

  • Rhines, P. B., 1977: The dynamics of unsteady currents. The Sea. Vol 6: Marine Modelling, E. Goldberg, Ed., Wiley, 189–318.

  • ——, and W. R. Holland, 1979: A theoretical discussion of eddy-driven flows. Dyn. Atmos. Oceans,3, 289–325.

  • Rix, N. H., and J. Willebrand, 1996: Parameterization of mesoscale eddies as inferred from a high-resolution circulation model. J. Phys. Oceanogr.,26, 2281–2285.

  • Schmitz, W. J., 1977: On the deep general circulation in the western North Atlantic. J. Mar. Res.,35, 21–28.

  • Smagorinsky, J. S., 1963: General circulation experiments with the primitive equations. I: The basic experiment. Mon. Wea. Rev.,91, 99–164.

  • Stammer, D., 1998: On eddy characteristics, eddy transports, and mean flow properties. J. Phys. Oceanogr.,28, 727–739.

  • Treguier, A. M., 1999: Evaluating eddy mixing coefficients from eddy-resolving ocean models: A case study. J. Mar. Res.,57, 89–108.

  • ——, I. Held, and V. Larichev, 1997: On the parameterization of quasi-geostrophic eddies in primitive equation ocean models. J. Phys. Oceanogr.,27, 567–580.

  • Visbeck, M., J. Marshall, and T. Haine, 1997: Specification of eddy transfer coefficients in coarse-resolution ocean circulation models. J. Phys. Oceanogr.,27, 381–402.

  • Wunsch, C., 1999: Where do ocean eddy heat fluxes matter? J. Geophys. Res.,104, 13 235–13 249.

All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 231 75 1
PDF Downloads 65 14 1