On the pressure gradient force over steep topography in sigma coordinate ocean models

Abstract

The error in computing the pressure gradient force near steep topography using terrain following (o) coordinates is investigated in an ocean model using the family of vertical differencing schemes proposed by Arakawa and Suarez. The truncation error is estimated by substituting known buoyancy profiles in the finite difference hydrostatic and pressure gradient terms. The error due to "hydrostatic inconsistency," which is not simply a space truncation error, is also documented. The results show that the pressure gradient error is spread throughout the water column, and it is sensitive to the vertical resolution and to the placement of the grid points relative to the vertical structure of the buoyancy field being modeled. Removing a reference state, as suggested for the atmosphere by Gary, reduces the truncation error associated with the two lowest vertical modes by a factor of 2 to 3. As an example, the error in computing the pressure gradient using a standard 10-level primitive equation model applied to buoyancy profiles and topographic slopes typical of the California Current region corresponds to a false geostrophic current of the order of 10-12 cm s(-1). The analogous error is a hydrostatically consistent 30-level model with the reference state removed is about an order of magnitude smaller.