Local Mass Conservation and Velocity Splitting in PV-Based Balanced Models. Part II: Numerical Results
Authors: Mohebalhojeh Ali R. , McIntyre Michael E. 
Published: June 2007

The effects of enforcing local mass conservation on the accuracy of non-Hamiltonian potential-vorticity- based balanced models (PBMs) are examined numerically for a set of chaotic shallow-water f-plane vortical flows in a doubly periodic square domain. The flows are spawned by an unstable jet and all have domain-maximum Froude and Rossby numbers Fr ∼0.5 and Ro ∼1, far from the usual asymptotic limits Ro → 0, Fr → 0, with Fr defined in the standard way as flow speed over gravity wave speed. The PBMs considered are the plain and hyperbalance PBMs defined in Part I. More precisely, they are the plain-δδ, plain-γγ, and plain-δγ PBMs and the corresponding hyperbalance PBMs, of various orders, where “order” is related to the number of time derivatives of the divergence equation used in defining balance and potential-vorticity inversion. For brevity the corresponding hyperbalance PBMs are called the hyper-δδ, hyper-γγ, and hyper-δγ PBMs, respectively. As proved in Part I, except for the leading-order plain-γγ each plain PBM violates local mass conservation. Each hyperbalance PBM results from enforcing local mass conservation on the corresponding plain PBM. The process of thus deriving a hyperbalance PBM from a plain PBM is referred to for brevity as plain-to-hyper conversion. The question is whether such conversion degrades the accuracy, as conjectured by McIntyre and Norton. Cumulative accuracy is tested by running each PBM alongside a suitably initialized primitive equation (PE) model for up to 30 days, corresponding to many vortex rotations. The accuracy is sensitively measured by the smallness of the ratio ϵ = ||QPBM − QPE||2/||QPE||2, where QPBM and QPE denote the potential vorticity fields of the PBM and the PEs, respectively, and || ||2 is the L2 norm. At 30 days the most accurate PBMs have ϵ ≈ 10−2 with PV fields hardly distinguishable visually from those of the PEs, even down to tiny details. Most accurate is defined by minimizing ϵ over all orders and truncation types δδ, γγ, and δγ. Contrary to McIntyre and Norton’s conjecture, the minimal ϵ values did not differ systematically or significantly between plain and hyperbalance PBMs. The smallness of ϵ suggests that the slow manifolds defined by the balance relations of the most accurate PBMs, both plain and hyperbalance, are astonishingly close to being invariant manifolds of the PEs, at least throughout those parts of phase space for which Ro ≲ 1 and Fr ≲ 0.5. As another way of quantifying the departures from such invariance, that is, of quantifying the fuzziness of the PEs’ slow quasimanifold, initialization experiments starting at days 1, 2, . . . 10 were carried out in which attention was focused on the amplitudes of inertia–gravity waves representing the imbalance arising in 1-day PE runs. With balance defined by the most accurate PBMs, and imbalance by departures therefrom, the results of the initialization experiments suggest a negative correlation between early imbalance and late cumulative error ϵ. In such near-optimal conditions the imbalance seems to be acting like weak background noise producing an effect analogous to so-called stochastic resonance, in that a slight increase in noise level brings PE behavior closer to the balanced behavior defined by the most accurate PBMs when measured cumulatively over 30 days.