Covariance localisation and balance in an Ensemble Kalman Filter

Quarterly Journal of the Royal Meteorological Society - Tập 135 Số 642 - Trang 1157-1176 - 2009
Jeffrey D. Kepert1
1Centre for Australian Weather and Climate Research, Bureau of Meteorology, Melbourne, Australia

Tóm tắt

AbstractA major limitation of the Ensemble Kalman Filter (EnKF) is that the finite ensemble size introduces sampling error into the background covariances, with severe consequences for atmospheric and oceanographic applications. The negative effects of sampling error are customarily limited by covariance localisation, which earlier studies have suggested may introduce imbalance into the system. The deleterious effects of localisation upon balance are confirmed and detailed here, with localisation producing analyses with weaker geostrophic balance and stronger divergence than are obtained using the unlocalised covariances. These imbalances reduce as the localisation radius is increased, but are argued to be large for typical settings. An improved method for calculating local covariances from an ensemble is presented, in which the localisation is performed in streamfunction–velocity potential (ψ–χ), rather than wind component, space. Analyses using this method better preserve the balances contained within the unlocalised covariance model. This transformation further allows the option of intervariable localisation, in which the cross‐covariances involving χ, which are weak and therefore particularly subject to sampling error, are set to 0 instead of being calculated from the ensemble. The various localisations are compared in a series of identical‐twin experiments, with the new localisations producing analyses that are better balanced and significantly more accurate than the usual approach. The localisation with the χ cross‐covariances set to 0 is shown to be superior for the smaller ensemble sizes but not for the larger, implying that the larger ensembles are capable of resolving some of the true χ cross‐covariance in the test system. Copyright © 2009 Royal Meteorological Society

Từ khóa


Tài liệu tham khảo

10.1016/j.physd.2006.02.011

10.1175/JAS3510.1

10.1175/1520-0493(2001)129<0420:ASWTET>2.0.CO;2

10.1175/1520-0493(1972)100<0683:AEOPSM>2.3.CO;2

10.1175/1520-0493(1998)126<1719:ASITEK>2.0.CO;2

10.2151/jmsj1965.75.1B_257

10.1175/1520-0493(1998)126<2913:ATEODS>2.0.CO;2

10.1175/1520-0493(1985)113<1066:TAOSSD>2.0.CO;2

Daley R, 1986, Variational Methods in Geosciences

Daley R, 1991, Atmospheric Data Analysis

10.1034/j.1600-0870.1999.t01-2-00003.x

Evensen G, 1994, Sequential data assimilation with a nonlinear quasi‐geostrophic model using Monte Carlo methods to forecast error statistics, J. Geophys. Res., 99, 10143, 10.1029/94JC00572

10.1002/qj.49712555417

10.1175/1520-0493(2001)129<2776:DDFOBE>2.0.CO;2

10.1111/j.1600-0870.1986.tb00460.x

10.1175/1520-0493(1998)126<0796:DAUAEK>2.0.CO;2

10.1175/1520-0493(2001)129<0123:ASEKFF>2.0.CO;2

10.1256/qj.05.135

Kalman RE, 1960, A new approach to linear filtering and prediction problems, Trans. ASME—J. Basic Eng., 82, 35, 10.1115/1.3662552

10.1111/j.1600-0870.1986.tb00461.x

10.1175/1520-0493(1981)109<0701:AGTDMS>2.0.CO;2

10.1256/qj.02.132

10.1175/1520-0493(2002)130<2791:ESBAME>2.0.CO;2

Neef LJ, 2008, A low‐order model investigation of the analysis of gravity waves in the Ensemble Kalman Filter, J. Atmos. Sci.

SteinlePJ.2005.Generalized Statistical Interpolation. In proceedings of 4th Symposium on Assimilation of Observations for Meteorology and Oceanography Prague Czech Republic 18–22 April. (CD‐ROM) World Meteorological Organisation: Geneva.

10.1175/1520-0493(1999)127<1374:CODAUA>2.0.CO;2

10.1175/1520-0493(2002)130<1913:EDAWPO>2.0.CO;2

10.1175/1520-0493(1981)109<2357:TBBMAW>2.0.CO;2