Estimation of the internal pressure gradient in σ-coordinate ocean models: comparison of second-, fourth-, and sixth-order schemes

Sigma-coordinate ocean models are attractive because of their abilities to resolve bottom and surface boundary layers. However, these models can have large internal pressure gradient (IPG) errors. In this paper, two classes of methods for the estimation of the IPGs are assessed. The first is based on the integral approach used in the Princeton Ocean Model (POM). The second is suggested by Shchepetkin and McWilliams (2003) based on Green’s theorem; thus, area integrals of the pressure forces are transformed into line integrals. Numerical tests on the seamount problem, as well as on a northwestern Atlantic grid using both classes of methods, are presented. For each class, second-, fourth-, and sixth-order approximations are tested. Results produced with a fourth-order compact method and with cubic spline methods are also given. The results show that the methods based on the POM approach in general give smaller errors than the corresponding methods given in Shchepetkin and McWilliams (2003). The POM approach also is more robust when noise is added to the topography. In particular, the IPG errors may be substantially reduced by using the computationally simple fourth-order method from McCalpin (1994).

[1]  L. Slørdal,et al.  THE PRESSURE GRADIENT FORCE IN SIGMA-CO-ORDINATE OCEAN MODELS , 1997 .

[2]  Dale B. Haidvogel,et al.  Numerical Simulation of Flow around a Tall Isolated Seamount. Part I: Problem Formulation and Model Accuracy , 1993 .

[3]  P. Chu,et al.  Hydrostatic correction for sigma coordinate ocean models , 2003 .

[4]  Robert L. Haney,et al.  On the Pressure Gradient Force over Steep Topography in Sigma Coordinate Ocean Models , 1991 .

[5]  Julie D. Pietrzak,et al.  On the pressure gradient error in sigma coordinate ocean models: A comparison with a laboratory experiment , 1999 .

[6]  G. Stelling,et al.  On the approximation of horizontal gradients in sigma co‐ordinates for bathymetry with steep bottom slopes , 1994 .

[7]  Y. Song,et al.  A General Pressure Gradient Formulation for Ocean Models. Part II: Energy, Momentum, and Bottom Torque Consistency , 1998 .

[8]  A. Schuster.  On electric currents induced by rotating magnets, and their application to some phenomena of terrestrial magnetism , 1896 .

[9]  J. Berntsen,et al.  Estimating the internal pressure gradient errors in a σ-coordinate ocean model for the Nordic Seas , 2007 .

[10]  Peter C. Chu,et al.  Sixth-order difference scheme for sigma coordinate ocean models , 1997 .

[11]  Hyun-chul Lee,et al.  Effects of winds and Caribbean eddies on the frequency of Loop Current eddy shedding: A numerical model study , 2003 .

[12]  George L. Mellor,et al.  The Pressure Gradient Conundrum of Sigma Coordinate Ocean Models , 1994 .

[13]  G. Mellor.  USERS GUIDE for A THREE-DIMENSIONAL, PRIMITIVE EQUATION, NUMERICAL OCEAN MODEL , 1998 .

[14]  A. M. Davies,et al.  Numerical studies of internal waves at a sill: Sensitivity to horizontal grid size and subgrid scale closure , 2008 .

[15]  N. Heaps,et al.  Three-dimensional coastal ocean models , 1987 .

[16]  J. Berntsen.  Internal Pressure Errors in Sigma-Coordinate Ocean Models , 2002 .

[17]  M. Levine,et al.  A Correction to the Baroclinic Pressure Gradient Term in the Princeton Ocean Model , 2001 .

[18]  A. Ciappa.  An operational comparative test of z-levels PGF schemes to reduce pressure gradient errors of the ocean model POM , 2006 .

[19]  Hernan G. Arango,et al.  Developments in terrain-following ocean models: intercomparisons of numerical aspects , 2002 .

[20]  André B. Fortunato,et al.  Evaluation of horizontal gradients in sigma-coordinate shallow water models , 1996 .

[21]  L. Kantha,et al.  Numerical models of oceans and oceanic processes , 2000 .

[22]  John M. Gary,et al.  Estimate of Truncation Error in Transformed Coordinate, Primitive Equation Atmospheric Models , 1973 .

[23]  John D. McCalpin.  A comparison of second‐order and fourth‐order pressure gradient algorithms in a σ‐co‐ordinate ocean model , 1994 .

[24]  Jarle Berntsen,et al.  Internal pressure errors in sigma-coordinate ocean models—sensitivity of the growth of the flow to the time stepping method and possible non-hydrostatic effects , 2005 .

[25]  L. Oey,et al.  Sigma Coordinate Pressure Gradient Errors and the Seamount Problem , 1998 .

[26]  M Israeli,et al.  Numerical Simulation of Viscous Incompressible Flows , 1974 .

[27]  S. G. L. Smith,et al.  Numerical and Analytical Estimates of M2 Tidal Conversion at Steep Oceanic Ridges , 2006 .

[28]  James C. McWilliams,et al.  A method for computing horizontal pressure‐gradient force in an oceanic model with a nonaligned vertical coordinate , 2003 .

[29]  M. D. Sikirić,et al.  A new approach to bathymetry smoothing in sigma-coordinate ocean models , 2009 .

[30]  J. Molines,et al.  A sigma-coordinate primitive equation model for studying the circulation in the South Atlantic. Part I: Model configuration with error estimates , 1998 .

[31]  Y. Song,et al.  A General Pressure Gradient Formulation for Ocean Models. Part I: Scheme Design and Diagnostic Analysis , 1998 .