i
Estimating observation impact signals in NCEP GSI using the Lanczos method
-
2010
-
Details:
-
Personal Author:
-
Corporate Authors:
-
NOAA Program & Office:
-
Description:Despite the tremendous progress that has been made in data assimilation (DA) methodology, observing systems which reduce observation errors, and model improvements which reduce background errors, the analyses produced by the best available DA systems are still far from the truth. Analysis error and error covariance are important since they describe the accuracy of the analyses, and directed related to the future forecast errors, i.e. the forecast quality. In addition, analysis error covariance is critically important in building an efficient ensemble forecast system (EFS). Estimating analysis error covariance in an ensemble based Kalman filter DA is easier in theory, but it is always challenging in variational DA systems which have been in operation at most NWP centers. In this study, we use the Lanczos method in the NCEP Gridpoint Statistical Interpolation (GSI) DA system to look into other important aspects and properties of this method which were not exploited before. We apply this method to estimate the observation impact signals (OIS) which are directly related to the analysis error variances. It is found that the smallest eigenvalue of the transformed Hessian matrix converges to one as the number of minimization iterations increases. When more observations are assimilated, the converging speed becomes slower and more eigenvectors are needed to retrieve the observation impacts. It is also found that the OIS over data-rich regions can be picked up easily by the eigenvectors with dominant eigenvalues. Since only a limited number of eigenvectors can be computed due to computational expense, the OIS is severely underestimated, and the analysis error variance is consequently overestimated. It is found that the mean OIS values for temperature and wind components at typical model levels are increased by about 1.5 times when the number of eigenvectors is doubled. We have proposed four different calibration schemes to compensate for the missing trailing eigenvectors. Results show that the method with calibration for smaller number of eigenvectors cannot pick up the observation impacts over the regions with less observations as well as an 'ideal case' with a large number of eigenvectors. But proper calibrations do enhance and improve the impact signals over regions with more data. When compared with the observation locations, the method generally captures the OIS over regions with more observation data, including satellite data over the southern oceans. Over the tropics, some observation impacts may be missed due to the smaller background errors specified in the GSI, which is not related to the method. It is found that a lot more eigenvectors are needed to retrieve impact signals that resemble the banded structures from satellite observations, particularly over the tropics. Another benefit from the Lanczos method is that the dominant eigenvectors can be used in preconditioning the conjugate gradient algorithm in the GSI to speed up the convergence.
-
Content Notes:Mozheng Wei, Manuel S.F.V. De Pondeca, Zoltan Toth, David Parrish.
"October 25, 2010."
System requirements: Adobe Acrobat Reader.
Includes bibliographical references (pages 15-17).
-
Keywords:
-
Series:
-
Document Type:
-
Place as Subject:
-
Rights Information:Public Domain
-
Compliance:Library
-
Main Document Checksum:
-
Download URL:
-
File Type: