Numerical error assessment and a temporal horizon for internal waves in a hydrostatic model

by Ben R. Hodges, Bernard Laval, and Bridget M. Wadzuk

Citation: Hodges, B.R., B. Laval, and B.M. Wadzuk, Numerical error assessment and a temporal horizon for internal waves in a hydrostatic model, Ocean Modelling, 13(1):44-64,doi:10.1016/j.ocemod.2005.09.005.

Abstract

Three-dimensional hydrostatic models are often used to represent stratified dynamics and large-scale circulations in lakes and coastal oceans. However, non-linear internal wave evolution is fundamentally non-hydrostatic, so any hydrostatic model will eventually fail if non-hydrostatic effects accumulate. In this paper, we demonstrate that the time scale of baroclinic wave steepening provides a temporal horizon for model skill when sufficiently fine model grid scales are applied. At coarser grid scales, model results may be dominated by either diffusion or dissipation prior to the time scale of steepening. Methods for estimating numerical viscosity and diffusivity are developed, along with approximate predictive scales for the Centre for Water Research Estuary and Lake Computer Model (CWR-ELCOM). Results indicate that numerical diffusion of the pycnocline (i.e. thermocline or halocline) may rapidly increase past the time scale of steepening for a well-resolved model, and that somewhat coarser horizontal grid scales may be preferable for integrations beyond the time scale of steepening. While the present application and analyses are for a specific model, the methodology is developed in a general formulation and could be applied to other models for insight into error accumulation..

< pdf of first page>

©2006 Ben R. Hodges • last updated April 19, 2006

UT College of EngineeringDept of Civil EngineeringCRWREFM Home