摘要

We present an iterative algorithm for estimating the observation-noise covariance matrix R-t for ensemble-based data assimilation. The proposed algorithm enables us to estimate parameters in Rt, the number of which may become exceedingly large for data assimilation. The optimal Rt is estimated based on the maximum likelihood and the algorithm can be described as an expectation-maximization (EM) algorithm. The key to constructing the algorithm is to obtain the optimality condition of Rt in a closed form. The proposed method estimates R-t, which is adaptive to the current observation y(t). We present an application of the proposed method to a coupled atmosphere-ocean model and find that only a few iterations are necessary for convergence. We also find that R-t needs to have some kind of structure to obtain a reasonable covariance matrix using the proposed algorithm.

  • 出版日期2014-1