On the use of the Boussinesq equations, the reduced system, and the primitive equations for the computation of geophysical flows

Abstract

Motivated in part by the mathematical problems associated with the application of open boundary conditions to the hydrostatic primitive equations (PE), Browning et al. (1990, Dyn. Atmos. Oceans , 14: 303–332) proposed the use of the reduced system (RS) of equations to replace PE for oceanographic problems. The RS are essentially the Boussinesq equations (BO) with the non-hydrostatic terms in the vertical momentum equation multiplied by a constant δ 2 ⪢ 1. This artificially alters the physics (e.g. changing the intemal-inertial wave properties) to facilitate numerical integration, but the changes are assumed to have negligible effects on the dynamics of interest. We assess the accuracy and utility of the RS (following the guidelines for the choice of δ) by comparing numerical finite difference solutions of RS, PE and BO for initial-value problems involving three-dimensional instability of an ocean front and atmospheric frontal development in a two-dimensional Eady wave. Both explicit (BO, PE) and semi-implicit (BOSI, PESI) time-difference schemes are used for the Boussinesq and primitive equations. For RS, the same explicit scheme as for BO is used where δ ⪢ 1 allows larger time steps than with the other explicit models. It is found that relative to BO solutions, the errors for RS are small but increase rapidly and monotonically with increasing δ (over a range consistent with the guidelines) and are greater than the errors for the other models. The use of BOSI allows time steps at least as large as those for RS and results in smaller errors than RS. For these problems, BOSI is the preferable model to replace PE.