Which Is Better, an Ensemble of Positive-Negative Pairs or a Centered Spherical Simplex Ensemble?

摘要

New methods to center the initial ensemble perturbations on the analysis are introduced and compared with the commonly used centering method of positive-negative paired perturbations. In the new method, one linearly dependent perturbation is added to a set of linearly independent initial perturbations to ensure that the sum of the new initial perturbations equals zero; the covariance calculated from the new initial perturbations is equal to the analysis error covariance estimated by the independent initial perturbations, and all of the new initial perturbations are equally likely. The new method is illustrated by applying it to the ensemble transform Kalman filter (ETKF) ensemble forecast scheme, and the resulting ensemble is called the spherical simplex ETKF ensemble. It is shown from a multidimensional Taylor expansion that the symmetric positive-negative paired centering would yield a more accurate forecast ensemble mean and covariance than the spherical simplex centering if the ensemble were large enough to span all initial uncertain directions and thus the analysis error covariance was modeled precisely. However, when the number of uncertain directions is larger than the ensemble size, the spherical simplex centering has the advantage of allowing almost twice as many uncertain directions to be spanned as the symmetric positive-negative paired centering. The performances of the spherical simplex ETKF and symmetric positive-negative paired ETKF ensembles are compared by using the Community Climate Model Version 3 (CCM3). Each ensemble contains 1 control forecast and 16 perturbed forecasts. The NCEP-NCAR reanalysis data for the boreal summer in 2000 are used for the initialization of the control forecast and the verifications of the ensemble forecasts. The accuracy of the ensemble means, the accuracy of predictions of forecast error variance, and the ability of the ETKF ensembles to resolve inhomogeneities in the observation distribution were all tested. In all of these test categories, the spherical simplex ETKF ensemble was found to be superior to the symmetric positive-negative paired ETKF ensemble. The computational expense for generating spherical simplex ETKF initial perturbations is about as small as that for the symmetric positive-negative paired ETKF. Also shown is that the seemingly straightforward centering method, in which centered perturbations are obtained by subtracting the average of the perturbations from each individual perturbation, is unsatisfactory because the covariance estimated by the uncentered perturbations is not necessarily conserved after centering.