Ensemble-type Kalman filter algorithm conserving mass, total energy and enstrophy

Yuefei Zeng, Tijana Janjic, Yvonne Ruckstuhl, Martin Verlaan

Research output: Contribution to journalArticleScientificpeer-review

6 Citations (Scopus)

Abstract

For numerical discretization schemes, the violation of enstrophy conservation causes a systematic and unrealistic energy cascade towards high wave numbers. The same occurs in data assimilation schemes, where the total energy, enstrophy and divergence could be strongly affected. In this article, we construct an ensemble data assimilation algorithm that conserves mass, total energy and enstrophy. The algorithm uses B-spline functions for localization and sequential quadratic programming to solve nonlinear constrained minimization problem. Idealized experiments are performed using a 2D shallow-water model, with selected contraints derived from the nature run. It is found that all experiments exhibit comparable root-mean-square errors, with a slight advantage for those that include the conservation constraint on the globally integrated enstrophy. However, the kinetic energy and enstrophy spectra in experiments with the enstrophy constraint are considerably closer to the true spectra, in particular at the smallest resolvable scales. Therefore, imposing conservation of enstrophy within the data assimilation algorithm effectively avoids the spurious energy cascade of the rotational part and thereby successfully suppresses the noise generated by the data assimilation algorithm. The 14 day deterministic free forecast, starting from the initial condition enforced by both total energy and enstrophy constraints, produces the best prediction. The same holds for the ensemble free forecasts.

Original languageEnglish
Pages (from-to)2902-2914
Number of pages13
JournalRoyal Meteorological Society. Quarterly Journal (online)
Volume143
Issue number708 Part A
DOIs
Publication statusPublished - 2017

Keywords

  • B-spline
  • conservation laws
  • enstrophy
  • kinetic energy spectrum
  • localization
  • noise
  • SQP
  • total energy

Fingerprint

Dive into the research topics of 'Ensemble-type Kalman filter algorithm conserving mass, total energy and enstrophy'. Together they form a unique fingerprint.

Cite this