Ensemble‐type Kalman filter algorithm conserving mass, total energy and enstrophy. (21st November 2017)

• Record Type: Journal Article
• Authors: Zeng, Yuefei; Janjić, Tijana; Ruckstuhl, Yvonne; Verlaan, Martin
• Is Part Of: Quarterly journal of the Royal Meteorological Society. Volume 143:Number 708(2017)
• Page Start: 2902
• Page End: 2914
• Publication Date: 2017-11-21
• Subjects: conservation laws -- total energy -- enstrophy -- noise -- localization -- B‐spline -- SQP -- kinetic energy spectrum; Meteorology -- Periodicals; 551.5
• DOI: 10.1002/qj.3142
• Languages: English
• ISSNs: 0035-9009

Abstract

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. Abstract : Unphysical noise can be generated by assimilating data. We construct an ensemble‐type Kalman filter algorithm that enforces the conservation of mass, total energy and enstrophy as equality constraints. Experiments with the enstrophy constraints E_BSP_Es and E_BSP_EnEs (the latter additionally imposing total energy conservation) avoid the spurious energy cascade of the rotational part, thereby suppressing noise and improving the quality of a 14 day free forecast, compared with experiments without the enstrophy constraints E_BSP_NO and E_BSP_En.