Numerical simulations of Rossby–Haurwitz waves

Simulations of Rossby–Haurwitz waves have been carried out using four different high‐resolution numerical shallow water models: a spectral model, two semi‐Langrangian models predicting wind components and potential vorticity respectively, and a finite‐volume model on a hexagonal–icosahedral grid. The simulations show that (i) unlike the nondivergent case, the shallow water Rossby–Haurwitz wave locally generates small‐scale features and so has a potential enstrophy cascade, and (ii) contrary to common belief, the zonal wavenumber 4 Rossby–Haurwitz wave is dynamically unstable and will eventually break down if initially perturbed. Implications of these results for the use of the Rossby–Haurwitz wave as a numerical model test case are discussed. The four models tested give very similar results, giving confidence in the accuracy and robustness of the results. The most noticeable difference between the models is that truncation errors in the hexagonal–icosahedral grid model excite the Rossby–Haurwitz wave instability, causing the wave to break down quickly, whereas for the other models in the configurations tested the instability is excited only by roundoff error at worst, and the Rossby–Haurwitz wave breaks down much more slowly or not at all.