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

Jarle Berntsen, Lie Yauw Oey

Research output: Contribution to journalArticlepeer-review

32 Scopus citations

Abstract

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).

Original languageEnglish
Pages (from-to)317-330
Number of pages14
JournalOcean Dynamics
Volume60
Issue number2
DOIs
StatePublished - Apr 2010

Keywords

  • Internal pressure
  • Ocean modeling
  • Sigma-coordinates

Fingerprint

Dive into the research topics of 'Estimation of the internal pressure gradient in σ-coordinate ocean models: Comparison of second-, fourth-, and sixth-order schemes'. Together they form a unique fingerprint.

Cite this