Advances in Meteorology

Volume 2013 (2013), Article ID 410812, 6 pages

http://dx.doi.org/10.1155/2013/410812

Research Article

A New Data Assimilation Scheme: The Space-Expanded Ensemble Localization Kalman Filter

Hongze Leng,¹ Junqiang Song,¹ Fengshun Lu,² and Xiaoqun Cao¹

¹College of Computer, National University of Defense Technology, Changsha 410073, China

²China Aerodynamics Research and Development Center, Mianyang, Sichuan 621000, China

Received 6 January 2013; Revised 23 February 2013; Accepted 9 March 2013

Academic Editor: Hann-Ming Henry Juang

Copyright © 2013 Hongze Leng et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

This study considers a new hybrid three-dimensional variational (3D-Var) and ensemble Kalman filter (EnKF) data assimilation (DA) method in a non-perfect-model framework, named space-expanded ensemble localization Kalman filter (SELKF). In this method, the localization operation is directly applied to the ensemble anomalies with a Schur Product, rather than to the full error covariance of the state in the EnKF. Meanwhile, the correction space of analysis increment is expanded to a space with larger dimension, and the rank of the forecast error covariance is significantly increased. This scheme can reduce the spurious correlations in the covariance and approximate the full-rank background error covariance well. Furthermore, a deterministic scheme is used to generate the analysis anomalies. The results show that the SELKF outperforms the perturbed EnKF given a relatively small ensemble size, especially when the length scale is relatively long or the observation error covariance is relatively small.

1. Introduction

Data assimilation (DA) combines the information from the short-term model forecast and the information from the observations to form an initial estimate of the state in a dynamic system. The three-dimensional variational (3D-Var) method [1, 2] is known as an effective scheme for DA system, and 4D-Var [3, 4] is an extension of 3D-Var with an extra time dimension. However, due to the use of a static background error covariance, the DA system with 3D-Var lacks flow-dependent error statistics, which will yield an inaccurate estimate of the state.

Another scheme, called Kalman filter (KF), as well as its variant extended Kalman filter (EKF) [5, 6], is then used, which can estimate both the state and its real-time uncertainty in a linear or weakly nonlinear dynamic system. The ensemble Kalman filter (EnKF) [7, 8] is a variant of the EKF, and it utilizes an ensemble-based covariance to approximate the forecast uncertainty. Hybrid schemes, regarded to absorb both the advantages of 3D-Var and EnKF, are attracting much attention these years, such as MLEF [9], En4DVar [10], incorporating ensemble covariance in GSI [11], PODEn4DVar [12], and IEnKF [13]. Localization is often used to reduce spurious correlations in ensemble-based covariance matrices [14–16]. However, research on directly localizing the ensemble anomalies and expanding the correction space of analysis increment is small, and the deterministic scheme of generating the analysis anomalies in a hybrid frame attracts little attention as well.

The aim of this paper is to present a new method, which can directly localize the ensemble anomalies, rather than the full forecast or the full analysis error covariance, and increase the rank of the ensemble-based covariance matrix. Therefore, a new hybrid 3D-Var and EnKF scheme, named space-expanded ensemble localization Kalman filter (SELKF), is proposed. In the new method, a localization operation is directly executed in the ensemble anomalies, rather than the full forecast error covariance or analysis error covariance. Meanwhile, the correction space of analysis increment is expanding from a space spanned by the ensemble perturbations to a new space with a larger dimension. Furthermore, a deterministic scheme is applied to form the analysis anomalies. We perform a number of experiments to compare the performance of the SELKF against the traditional (perturbed) EnKF, exploring its capabilities in different ensemble sizes, length scales of the correlation, inflations, observation time intervals, and observation error variances.

2. The New Scheme: SELKF

We firstly define some notations used in this paper as follows:

 I: number of ensemble members; N: dimension of state vector; M: dimension of observation vector; r: rank of the truncated square root of correlation matrix S.

In the EnKF, the ensemble-based covariances are defined as

where and denote the forecast error covariance and the analysis error covariance, respectively; and indicate the forecast ensemble anomalies and the analysis ensemble anomalies, where each column is the deviation of an ensemble member to the mean, respectively.

Unlike the perturbed EnKF, the ensemble transformation Kalman filter (ETKF) [14] used a deterministic scheme to generate the analysis anomalies, in which the transformation from forecast ensemble anomalies to analysis ensemble anomalies can be written as

where V and are the eigenvector matrix and the eigenvalue matrix of , respectively. ETKF is theoretically equivalent to a 3D-Var DA method in the linear framework, named MLEF [9], in which

where is a control variable with dimension , which is the size of ensemble, and is the observation operator. MLEF constructs analysis ensemble anomalies with the same form as (2).

The modification to the forecast state in the MLEF (or the ETKF) is restricted to the ensemble space. If the ensemble size is relatively small, one can only obtain a suboptimal analysis due to the restricted subspace. In addition, in the original ETKF, there is no localization, which is known as a powerful method to reduce sample errors by lessening spurious correlation in an error covariance [15].

In the following, a space-expanded ensemble localization Kalman filter (SELKF) is presented. In SELKF, the correction space is expanded from to Ir, where is a spectral truncated dimension of a correlation matrix S, and the localization operation is mixed into the ensemble anomalies directly, as was done by Buehner [17], rather than being applied into the full error covariance.

Introduce a new control variable and define

where indicates a diagonal matrix with the elements of in the diagonal, is the deviation of from the mean , and denotes the truncated square root of correlation matrix . is an matrix, and w is a vector with dimension Ir

where is a vector with dimension , corresponding to .

The following identity can be derived [17]:

If is linear or weakly nonlinear, we define

where

Thus, we can obtain a new cost function and its gradient

The analysis solution can then be calculated

From (4), we can get a relationship between and as follows:

where is the analysis error covariance with respect to w, and it is equal to the inversion of the Hessian matrix of cost function (9).

We define

where .

Thus, the analysis ensemble anomalies can be calculated

in which

and is a vector with the corresponding elements of .

Finally, the analysis ensemble is deterministically generated as follows:

To reduce the computation cost, some transformations can be made as done in [9]

where and are the eigenvector matrix and the eigenvalue matrix of , respectively.

A similar scheme [18–20] (hereafter, we call the scheme LW) was proposed to localize in the ensemble perturbation space. However, there are two major differences between the LW and the SELKF. (1) The dimension of the control variable in the SELKF is Ir, which is much smaller than IN in the LW, due to the use of a truncated spectral expansion in the SELKF. Therefore, the computational cost is significantly reduced. (2) The way to generate the analysis anomalies is quite different. In Lorenc 2003, no deterministic scheme is utilized to update the analysis ensemble, and in Wang et al. 2008a, 2008b, an ETKF scheme is directly used to form the analysis anomalies. However in this paper, the analysis anomalies are generated from the localized square root, , of the analysis error covariance matrix, as is shown in functions (14) and (15). This scheme keeps the property of localization and explicitly reduces the spurious correlations in the analysis error covariance matrix.

3. Experiments and Results

The performance of the new method will be investigated in the Lorenz96 system, which is often used as an approximate to real atmospheric systems [12, 21]. It is a chaos system in which a small perturbation of the initial state can result in a much different forecast.

The system can be written as

where , and . The initial state of this simulation is defined as

The model used in this experiment is not perfect, with an error covariance of , where A is a matrix with 1 at the diagonal and 0.5 at the sub- and super-diagonals. The observation operator is defined as , in which is a perturbed value with variance of . We assume that all grid points are observed every time steps. The observations are generated by adding a random to the truth at observation times. We assume that each variable is observed independently; that is, the observation error covariance is diagonal.

Here, we use a compactly supported weighting function as a Schur Product function to reduce the spurious correlation as follows:

where , denotes the length scale, and indicates the distance between two points (observation points and grid points) [22]. For the approximation of the correlation matrix S in the SELKF, the truncation strategy is simply set to choose the eigenvalues equaling to or bigger than 1 and the eigenvectors corresponding to those eigenvectors, while the full S is used in the EnKF. The inflation strategy for both methods is simply defined as

where .

The comparison of the new method with the EnKF will be made after a 2000-step spin-up of the model. The performance with different length scales of the correlation matrix and inflation factors will be investigated firstly. Every RMSE (the root mean square error) is an average of 50 runs over 20000 simulation time steps. We use in the first experiment.

Figure 1 shows the RMSE versus the length scale and the covariance inflation factor . It can be seen that when the ensemble size is large enough, for example, , the SELKF performs as well as the EnKF; however, it performs better than the EnKF as the number of ensemble members decreases, and the smaller the ensemble size, the larger the performance difference. It is interesting that the optimal performance of the EnKF is achieved at a relatively small inflation factor, especially with a relatively large ensemble size, for example, for , whereas the best results of SELKF are obtained at a much larger inflation factor, for example, for . However, when the length scale is increasing, a little larger inflation factor could be used to obtain a better performance for the EnKF, as is shown in Figure 1(b). It should be noted that with the decreasing of ensemble members, the optimal inflation factor of the EnKF increases to a much larger one, for example, 2 through 3 when .

(a)   SELKF,

(b) EnKF,

(c) SELKF,

(d) EnKF,

(e) SELKF,

(f) EnKF,

Figure 1: Parameter dependence of the analysis RMSE of two methods. The horizontal and the vertical axes indicate the inflation factor and the length scale, respectively.

It is also noted that when the length scale is relatively large, for example, through 10, the EnKF diverges severely at a relatively small ensemble size. In contrast, the SELKF can get a good performance with an appropriate inflation factor, even though the length scale is very large. Smaller length scales will reduce more spurious correlation in the covariance; however, it destroys physical balance at the same time [23]. It is therefore shown that the SELKF is more applicable than the EnKF especially when the ensemble size is limited and the length scale is relatively large.

Figures 1(a) and 1(b) demonstrate the comparable performance between the SELKF and the EnKF when the ensemble size is large enough. In the following, we will investigate the sensitivity of the two methods to varying observation time intervals and observation errors. Here, we choose the parameters that result in best performance for the two methods, respectively: for the SELKF and for the EnKF, and for both.

It is shown in Table 1 that the shorter the observation time interval, the better the performance for both methods, as expected. When the observation error variance is relatively large, the EnKF performs a litter better than the SELKF in a relatively small observation interval. If becomes larger, the SELKF can get a better performance than the EnKF even when is relatively large. However, when is small enough, the improvement of performance for the SELKF is significantly great, whereas the EnKF performs much worse than the SELKF, even worse than itself in the case of relatively large . The reason is partly because the EnKF takes use of the observation perturbations to generate analysis ensemble. Once the observation error variance is not large enough, the space spanned by the analysis ensemble cannot accommodate the truth. In contrast, the deterministic way of constructing the analysis anomalies in the SELKF can avoid that problem and benefit from the better accuracy of .

Table 1: Sensitivity to the observation interval and the observation error. The first value of each pair indicates the RMSE of SELKF, while the second one is the RMSE of EnKF.

0.384, 1.019
0.474, 0.746
0.607, 0.581
0.864, 0.737

0.480, 1.087
0.584, 0.794
0.732, 0.723
0.962, 0.925

0.562, 1.224
0.674, 0.867
0.834, 0.839
1.072, 1.086

0.641, 1.423
0.764, 0.961
0.942, 0.946
1.202, 1.238

4. Summary and Discussion

In order to reduce the spurious correlation in the ensemble forecast covariance and expand the limited correction space of the analysis increment in the EnKF, a new method called SELKF is presented.

This new method is compared with the traditional EnKF in the Lorenz96 system with model errors, and the results demonstrate that when the ensemble size is relatively large, the EnKF performs as well as the SELKF if the length scale is relatively small and the observation error covariance stays in a proper range. With decreasing ensemble size, the SELKF will outperform the EnKF, especially when the length scale is relatively large, where the EnKF diverges rapidly. This may be partly because the SELKF uses an expanded correction space, which can accommodate the truth of the analysis increment within the space spanned by the expanded ensemble anomalies, even given a relatively small ensemble size. It is noted that the SELKF performs much better than the EnKF when the observation error covariance is very small. This is partly due to the use of a deterministic scheme of forming the analysis anomalies. The SELKF uses exact factorization of the analysis error covariance and thus avoids the suboptimality resulting from the stochastic nature of the perturbed EnKF.

However, there are some shortcomings in the SELKF. If the length scale is very small, the truncation of the correlation matrix will yield a relatively large truncated rank , which results in a relatively expensive computation. This indicates that the SELKF would better be applied to systems with relatively long length scales. In addition, if the ensemble size is relatively large, the best performance of the SELKF is achieved at a much larger inflation factor than in the EnKF, the reason of which is still an open question.

Acknowledgments

This work was partly supported by the National Science Foundation of China (Grant no. 41105063). The authors are grateful to the anonymous reviewers for helpful comments.

References

N. Gustafsson, L. Berre, S. Hörnquist et al., “Three-dimensional variational data assimilation for a limited area model, part I: general formulation and the background error constraint,” Tellus A, vol. 53, no. 4, pp. 425–446, 2001.

S. A. Haben, A. S. Lawless, and N. K. Nichols, “Conditioning and preconditioning of the variational data assimilation problem,” Computers and Fluids, vol. 46, no. 1, pp. 252–256, 2011.

D. N. Daescu and I. M. Navon, “A dual-weighted approach to order reduction in 4DVAR data assimilation,” Monthly Weather Review, vol. 136, no. 3, pp. 1026–1041, 2008.

X. Y. Huang, Q. Xiao, D. M. Barker et al., “Four-dimensional variational data assimilation for WRF: formulation and preliminary results,” Monthly Weather Review, vol. 137, no. 1, pp. 299–314, 2009.

R. E. Kalman, “A new approach to linear filtering and prediction problems,” Journal of Basic Engineering D, vol. 82, pp. 35–45, 1960.

M. L. Psiaki, “Backward-smoothing extended Kalman filter,” Journal of Guidance, Control, and Dynamics, vol. 28, no. 5, pp. 885–894, 2005.

G. Evensen, “Sequential data assimilation with a nonlinear quasi-geostrophic model using Monte Carlo methods to forecast error statistics,” Journal of Geophysical Research, vol. 99, no. C5, pp. 10143–10162, 1994.

K. Bergemann and S. Reich, “Amollified ensemble Kalman filter,” Quarterly Journal of the Royal Meteorological Society, vol. 136, no. 651, pp. 1636–1643, 2010.

M. Zupanski, “Maximum likelihood ensemble filter: theoretical aspects,” Monthly Weather Review, vol. 133, no. 6, pp. 1710–1726, 2005.

C. Liu, Q. Xiao, and B. Wang, “An ensemble-based four-dimensional variational data assimilation scheme, part I: technical formulation and preliminary test,” Monthly Weather Review, vol. 136, no. 9, pp. 3363–3373, 2008.

X. Wang, “Incorporating ensemble covariance in the gridpoint statistical interpolation variational minimization: a mathematical framework,” Monthly Weather Review, vol. 138, no. 7, pp. 2990–2995, 2010.

X. Tian, Z. Xie, and Q. Sun, “A POD-based ensemble four-dimensional variational assimilation method,” Tellus A, vol. 63, no. 4, pp. 805–816, 2011.

P. Sakov, D. S. Oliver, and L. Bertino, “An iterative EnKF for strongly nonlinear systems,” Monthly Weather Review, vol. 140, pp. 1988–2004, 2012.

C. H. Bishop, B. J. Etherton, and S. J. Majumdar, “Adaptive sampling with the ensemble transform Kalman filter, part I: theoretical aspects,” Monthly Weather Review, vol. 129, no. 3, pp. 420–436, 2001.

P. L. Houtekamer and H. L. Mitchell, “A sequential ensemble Kalman filter for atmospheric data assimilation,” Monthly Weather Review, vol. 129, no. 1, pp. 123–137, 2001.

J. L. Anderson, “A local least squares framework for ensemble filtering,” Monthly Weather Review, vol. 131, pp. 634–642, 2003.

M. Buehner, “Ensemble-derived stationary and flow-dependent background-error covariances: evaluation in a quasi-operational NWP setting,” Quarterly Journal of the Royal Meteorological Society, vol. 131, no. 607, pp. 1013–1043, 2005.

A. C. Lorenc, “The potential of the ensemble Kalman filter for NWP—a comparison with 4D-Var,” Quarterly Journal of the Royal Meteorological Society, vol. 129, no. 595, pp. 3183–3203, 2003.

X. Wang, D. M. Barker, C. Snyder, and T. M. Hamill, “A hybrid ETKF-3DVAR data assimilation scheme for the WRF model, part I: observation system simulation experiment,” Monthly Weather Review, vol. 136, no. 12, pp. 5116–5131, 2008.

X. Wang, D. M. Barker, C. Snyder, and T. M. Hamill, “A hybrid ETKF-3DVAR data assimilation scheme for the WRF model, part II: real observing experiments,” Monthly Weather Review, vol. 136, no. 12, pp. 5132–5147, 2008.

P. J. van Leeuwen, “Nonlinear data assimilation in geosciences: an extremely efficient particle filter,” Quarterly Journal of the Royal Meteorological Society, vol. 136, no. 653, pp. 1991–1999, 2010.

T. M. Hamill, J. S. Whitaker, and C. Snyder, “Distance-dependent filtering of background error covariance estimates in an ensemble Kalman filter,” Monthly Weather Review, vol. 129, no. 11, pp. 2776–2790, 2001.

G. Evensen, Ed., Data Assimilation: The Ensemble Kalman Filter, Springer, Berlin, Germany, 2nd edition, 2009.