Preconditioning of variational data assimilation and the use of a bi-conjugate gradient method

Quarterly Journal of the Royal Meteorological Society

Volumn 139, Issue 672, 2013, Pages 731-741

  El Akkraoui, Amal ^a,b   Trémolet, Yannick ^c   Todling, Ricardo ^a  

^a NASA HEADQUARTERS   (United States)

^b SCIENCE SYSTEMS AND APPLICATIONS INC   (United States)

^c EUROPEAN CENTRE FOR MEDIUM RANGE WEATHER FORECASTS   (United Kingdom)

Author keywords

4D Var; GSI; Lanczos method; Minimization algorithms; Preconditioned conjugate gradient methods

Indexed keywords

COMPUTATIONAL EFFICIENCY; COVARIANCE MATRIX; EIGENVALUES AND EIGENFUNCTIONS; NUMBER THEORY; VALUE ENGINEERING;

4D-VAR; BI-CONJUGATE GRADIENT METHODS; CONJUGATE GRADIENT ALGORITHMS; FOUR DIMENSIONAL VARIATIONAL; LANCZOS METHODS; MINIMIZATION ALGORITHMS; PRECONDITIONED CONJUGATE GRADIENT METHOD; VARIATIONAL DATA ASSIMILATION;

CONJUGATE GRADIENT METHOD;

ALGORITHM; COVARIANCE ANALYSIS; DATA ASSIMILATION; EIGENVALUE; ERROR ANALYSIS; FOUR-DIMENSIONAL MODELING; GRADIENT ANALYSIS; METEOROLOGY;

EID: 84897957200     PISSN: 00359009     EISSN: 1477870X     Source Type: Journal    
DOI: 10.1002/qj.1997     Document Type: Article

References (50)

1. Andersson E, Fisher M, Munro R, McNally AP. 2000. Diagnosis of background errors for radiances and other observable quantities in a variational data assimilation scheme, and the explanation of a case of poor convergence. Q. J. R. Meteorol. Soc. 126: 1455-1472.
2. Courtier P, Talagrand O. 1987. Variational assimilation ofmeteorological observations with the adjoint vorticity equation. II: Numerical results. Q. J. R. Meteorol. Soc. 113: 1329-1347.
3. Courtier P, Thépaut J-N, Hollingsworth A. 1994. A strategy operational implementation of 4D-Var using an incremental approach. Q. J. R. Meteorol. Soc. 120: 1367-1387.
4. Dennis JE, Schnabel RB. 1983. Numerical methods for unconstrained optimization and nonlinear equations. Prentice-Hall: Englewood Cliffs, NJ. Reprinted 1996 as Classics in Applied Mathematics 16: SIAM: Philadelphia, U.S.A.
5. Derber JC, Rosati A. 1989. A global oceanic data assimilation system. J. Phys. Oceanogr. 19: 1333-1347.
6. Derber JC, Wu W-S. 1998. The use of TOVS cloud-cleared radiances in the NCEP SSI analysis system. Mon. Weather Rev. 126: 2287-2299.
7. Desroziers G, Hello G, Thépaut J-N. 2003. A 4D-Var reanalysis of the FASTEX. Q. J. R. Meteorol. Soc. 129: 1301-1325.
8. El Akkraoui A, Gauthier P. 2010. Convergence properties of the primal and dual forms of variational data assimilation. Q. J. R. Meteorol. Soc. 136: 107-115.
9. El Akkraoui A, Gauthier P, Pellerin S, Buis S. 2008. Intercomparison of the primal and dual formulations of variational data assimilation. Q. J. R. Meteorol. Soc. 134: 1015-1025.
10. Fisher M. 2003. 'Background-error covariancemodelling'. In Proceedings of Workshop on Recent developments in data assimilation for atmosphere and ocean, ECMWF: Reading, U.K.
11. Fisher M, Courtier P. 1995. 'Estimating the covariance matrices of analysis and forecast error in variational data assimilation'. Technical Memorandum 220, ECMWF: Reading, U.K.
12. Fletcher R. 1976. Conjugate gradient methods for indefinite linear systems. In Proc. Dundee Biennial Conf. Numer. Anal. Watson GA (ed.) Lecture Notes in Mathematics 506: 73-89. Springer-Verlag: Berlin.
13. Fletcher R, Reeves C. 1964. Function minimization by conjugate gradients. Computer J. 7: 81-84.
14. Freund RW, Golub GH, Nachtigal NM. 1992. Iterative solution of linear systems. Acta Numerica 1: 57-100.
15. Gauthier P, Thépaut J-N. 2001. Impact of the digital filter as a weak constraint in pre-operational 4D-Var assimilation system of Météo- France. Mon. Weather Rev. 129: 2089-2102.
16. Gauthier P, Tanguay M, Laroche S, Pellerin S, Morneau J. 2007. Extension of 3D-Var to 4D-Var: implementation of 4D-Var at the Meteorological Service of Canada. Mon. Weather Rev. 135: 2339-2354.
17. Gejadze I, LeDimet F-X, Shutyaev V. 2008. On analysis error covariances in variational data assimilation. SIAM J. Sci. Comput. 30: 1847-1874.
18. Gilbert JC, Lemaréchal C. 1989. Some numerical experiments with variable-storage quasi-Newton algorithms. Math. Program. 45: 407-435.
19. Gill PE, MurrayW,Wright MH. 1981. Practical optimization. Academic Press: San Diego, California.
20. Golub GH, Van Loan CF. 1996. Matrix computations. Johns Hopkins University Press: Baltimore and London.
21. Greenbaum A. 1997. Iterative methods for solving linear systems. SIAM: Philadelphia.
22. Hestenes MR, Stiefel E. 1952. Methods of conjugate gradients for solving linear systems. J. Res. National Bureau Standards 49: 409-436.
23. Honda Y, Nishijima M, Koizumi K, Ohta Y, Tamiya K, Kawabata T, Tsuyuki T. 2005. A pre-operational variational data assimilation system for a non-hydrostatic model at the Japan Meteorological Agency: Formulation and preliminary results. Q. J. R. Meteorol. Soc. 131: 3465-3475.
24. Janisková M, Thépaut J-N, Geleyn J-F. 1999. Simplified and regular physical parametrizations for incremental four-dimensional variational assimilation. Mon. Weather Rev. 127: 26-45.
25. Kleespies TJ, van Delst P, McMillin LM, Derber JC. 2004. Atmospheric transmittance of an absorbing gas. 6. OPTRAN status report and introduction to the NESDIS/NCEP Community Radiative Transfer Model. Appl. Opt. 43: 3103-3109.
26. Klinker E,Rabier F,Kelly G, Mahfouf J-F. 2000. The ECMWFoperational implementation of four-dimensional variational assimilation. III: Experimental results and diagnostics with operational configuration. Q. J. R. Meteorol. Soc. 126: 1191-1215.
27. Lanczos C. 1950. An iteration method for the solution of the eigenvalue problem of linear differential and integral operators. J. Res. National Bureau Standards 45: 255-282.
28. LeDimet FX, Talagrand O. 1986. Variational algorithms for analysis and assimilation ofmeteorological observations - Theoretical aspects. Tellus 38A: 97-110.
29. Lewis J,Derber J. 1985. The use of adjoint equations to solve a variational adjustment problem with advection constraints. Tellus 37A: 309-322.
30. Lewis J, Lakshmivarahan S, Dhall SK. 2006. Dynamic data assimilation: a least squares approach. In Encyclopedia of Mathematics. Cambridge University Press: Cambridge, U.K.
31. Lorenc AC. 1988. Optimal nonlinear objective analysis. Q. J. R. Meteorol. Soc. 114: 205-240.
32. Lorenc AC. 1997. Development of an operational variational assimilation scheme. J. Meteorol. Soc. Japan 75: 339-346.
33. Mahfouf J-F, Rabier F. 2000. The ECMWF operational implementation of four-dimensional variational assimilation. II: Experimental results with improved physics. Q. J. R. Meteorol. Soc. 126: 1171-1190.
34. Navon IM, LeglerDM. 1987. Conjugate-gradient methods for large-scale minimization in meteorology. Mon. Weather Rev. 115: 1479-1502.
35. Paige CC, Saunders MA. 1975. Solution of sparse indefinite systems of linear equations. SIAM J. Numer. Anal. 12: 617-629.
36. Parish DF, Derber JC. 1992. The National Meteorological Center's spectral statistical analysis system. Mon. Weather Rev. 120: 1747-1763.
37. Parlett BN, Taylor DR, Liu ZA. 1985. A look-ahead Lanczos algorithm for unsymmetric matrices. Math. Comput. 44: 105-124.
38. Rabier F, Järvinen H, Klinker E, Mahfouf J-F, Simmons A. 2000. The ECMWFoperational implementation of four-dimensional variational assimilation. I: Experimental results with simplified physics. Q. J. R. Meteorol. Soc. 126: 1143-1170.
39. Rawlins F, Ballard SP, Bovis KJ, Clayton AM, Li D, Inverarity GW, Lorenc AC, Payne TJ. 2007. The Met Office global four-dimensional variational data assimilation scheme. Q. J. R. Meteorol. Soc. 133: 347-362.
40. Rienecker MM, Suarez MJ, Gelaro R, Todling R, Bacmeister J, Liu E, Bosilovich MG, Schubert SD, Takacs L, Kim G-K, Bloom S, Chen J, Collins D, Conaty A, da Silva A, Gu W, Joiner J, Koster RD, Lucchesi R, Molod A, Owens T, Pawson S, Pegion P, Redder CR, Reichle R, Robertson FR, Ruddick AG, Sienkiewicz M, Woollen J. 2011. MERRA - NASA's Modern-Era Retrospective Analysis for Research and Applications. J. Climate 24: 3624-3648.
41. Saad Y. 1996. Iterative Methods for Sparse Linear Systems. PWSPublishing Co: Boston, Mass.
42. Saad Y, Schultz MH. 1986. GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems. SIAM J. Sci. Statist. Comput. 7: 856-869.
43. Sasaki Y. 1970. Some basic formalisms in numerical variational analysis. Mon. Weather Rev. 98: 875-883.
44. Sonneveld P. 1989. CGS: A fast Lanczos-type solver for nonsymmetric linear systems. SIAM J. Sci. Statist. Comput. 10: 36-52.
45. Thacker WC. 1989. The role of the Hessian matrix in fitting models to measurements. J. Geophys. Res. 94: C5, 6177-6196.
46. Trémolet Y. 2006. Accounting for an imperfect model in 4D-Var. Q. J. R. Meteorol. Soc. 132: 2483-2504.
47. Trémolet Y. 2007. Incremental 4D-Var convergence study. Tellus 59A: 706-718.
48. Van der Vorst HA. 1992. BiCGSTAB: A fast and smoothly converging variant of BiCG for the solution of nonsymmetric linear systems. SIAM J. Sci. Statist. Comput. 13: 631-644.
49. Wu W, Purser RJ, Parrish DF. 2002. Three-dimensional variational analysis with spacially inhomogeneous covariances. Mon. Weather Rev. 130: 2905-2916.
50. Zou X, Liu H, Derber J, Sela JG, Treadon R, Navon IM, Wang B. 2001. Four-dimensional variational data assimilation with a diabatic version of the NCEP global spectral model: System development and preliminary results. Q. J. R. Meteorol. Soc. 127: 1095-1122.