Comparison of overflow simulations on different vertical grids using the Finite Element Ocean circulation Model

Abstract The Finite Element Ocean circulation Model (FEOM) is applied to study the sensitivity of density driven overflows to the vertical discretization and bottom topography representation using the dynamics of overflow mixing and entrainment (DOME) setup. FEOM allows for hybrid grids combining σ, z + σ, full cell, partly shaved cell and fully shaved cell grids within the same numerical kernel thus isolating as far as possible effects of mesh geometry from those of model numerics. The sensitivity of diapycnal mixing, entrainment, plume thickness and plume meridional distribution to vertical discretization and partly to the subgrid process parameterization is explored. It is shown that simulations on pure σ grids or the combination of z + σ resolve the overflow processes best in terms of downslope plume propagation, plume thickness and dilution, and also have the least resolution dependence. Grids using z-levels generate excessive spurious mixing when resolution is insufficient. Applying partial cells improves the plume representation, but still requires higher horizontal and vertical resolution to converge to the σ grid results. It is demonstrated that increasing lateral viscosity causes the plume thickness to reduce whereas increasing lateral diffusivity has opposite effect. When keeping the Prandtl number constant, the increase in diffusivity and viscosity leads to an increase in mixing and plume thickness on z-level grids and also on σ-grids when lateral dissipation is oriented along geopotential surfaces. Using the along σ- diffusion helped to obtain correct plume thickness and entrainment on σ grids. Increasing the vertical mixing coefficients leads to an increase in diapycnal mixing and in downslope penetration as well.

[1]  Lin Jiang,et al.  Three-Dimensional Simulations of Overflows on Continental Slopes , 1996 .

[2]  Qiang Wang,et al.  Finite element ocean circulation model based on triangular prismatic elements, with application in studying the effect of topography representation , 2008 .

[3]  J. Girton,et al.  Descent and Modification of the Overflow Plume in the Denmark Strait , 2003 .

[4]  Qiang Wang,et al.  The finite element ocean model and its aspect of vertical discretization , 2007 .

[5]  Robert Hallberg,et al.  Simulation of Density-Driven Frictional Downslope Flow in Z-Coordinate Ocean Models , 1998 .

[6]  Jens Schröter,et al.  A finite-element ocean model: principles and evaluation , 2004 .

[7]  H. Goosse,et al.  Parameterization of density-driven downsloping flow for a coarse-resolution ocean model in z-coordinate , 1999 .

[8]  R. Hallberg,et al.  Numerical simulations of the Mediterranean sea outflow: impact of the entrainment parameterization in an isopycnic coordinate ocean model , 2003 .

[9]  P. Killworth,et al.  A Turbulent Bottom Boundary Layer Code for Use in Numerical Ocean Models , 1999 .

[10]  S. Legg,et al.  Regional simulations of the Faroe Bank Channel overflow in a level model , 2007 .

[11]  Robert Hallberg,et al.  Time Integration of Diapycnal Diffusion and Richardson Number–Dependent Mixing in Isopycnal Coordinate Ocean Models , 2000 .

[12]  J. Jungclaus,et al.  Cyclogenesis in the Denmark Strait Overflow Plume , 2001 .

[13]  R. Döscher,et al.  A Method for Improved Representation of Dense Water Spreading over Topography in Geopotential-Coordinate Models , 1997 .

[14]  George L. Mellor,et al.  A generalized coordinate ocean model and a comparison of the bottom boundary layer dynamics in terrain-following and in z-level grids , 2004 .

[15]  Anand Gnanadesikan,et al.  Transient Response in a Z-Level Ocean Model That Resolves Topography with Partial Cells , 1998 .

[16]  P. Killworth On the rate of descent of overflows , 2001 .

[17]  Tal Ezer,et al.  Entrainment, diapycnal mixing and transport in three-dimensional bottom gravity current simulations using the Mellor–Yamada turbulence scheme , 2005 .

[18]  Y. Chao,et al.  An Embedded Bottom Boundary Layer Formulation for Z-Coordinate Ocean Models , 2000 .

[19]  G. Mellor,et al.  A three-dimensional model study of the Mediterranean outflow , 2000 .

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

[21]  C. Adduce,et al.  Mixing induced in a dense current flowing down a sloping bottom in a rotating fluid , 2005 .

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

[23]  J. Verron,et al.  How Topographic Smoothing Contributes to Differences between the Eddy Flows Simulated by Sigma- and Geopotential-Coordinate Models , 2002 .

[24]  D. Dietrich,et al.  Dissipation effects in North Atlantic Ocean modeling , 2004 .

[25]  Molly O. Baringer,et al.  Outflows and deep water production by marginal seas , 1994 .

[26]  R. Pacanowski,et al.  Parameterization of Vertical Mixing in Numerical Models of Tropical Oceans , 1981 .

[27]  Rüdiger Gerdes,et al.  A primitive equation ocean circulation model using a general vertical coordinate transformation: 2. Application to an overflow problem , 1993 .

[28]  D. Dietrich,et al.  NOTES AND CORRESPONDENCE Entrainment and Transport in Idealized Three-Dimensional Gravity Current Simulation , 2005 .

[29]  D. Haidvogel,et al.  A numerical simulation of flow at Fieberling Guyot , 1997 .

[30]  G. Mellor,et al.  Development of a turbulence closure model for geophysical fluid problems , 1982 .

[31]  J. Peraire,et al.  Finite Element Flux-Corrected Transport (FEM-FCT) for the Euler and Navier-Stokes equations , 1987 .

[32]  R. Käse,et al.  Structure and variability of the Denmark Strait Overflow: Model and observations , 2003 .

[33]  Sonya Legg,et al.  Comparison of entrainment in overflows simulated by z-coordinate, isopycnal and non-hydrostatic models☆ , 2006 .

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