B-preconditioned minimization algorithms for variational data assimilation with the dual formulation

Quarterly Journal of the Royal Meteorological Society

Volumn 140, Issue 679, 2014, Pages 539-556

  Gurol S ^a   Weaver, A T ^a   Moore, A M ^b   Piacentini, A ^a   Arango, H G ^c   Gratton, S ^a,d  

^a CERFACS   (France)

^b UNIVERSITY OF CALIFORNIA   (United States)

^c RUTGERS UNIVERSITY   (United States)

^d UNIVERSITÉ DE TOULOUSE   (France)

Author keywords

3D Var; 4D Var; Conjugate gradient method; Dual approach; Lanczos method; Ocean data assimilation; PSAS

Indexed keywords

ALGORITHMS; CONJUGATE GRADIENT METHOD; COVARIANCE MATRIX; VALUE ENGINEERING; VECTOR SPACES;

3D-VAR; 4D-VAR; DUAL APPROACH; LANCZOS METHODS; OCEAN DATA ASSIMILATION; PSAS;

ITERATIVE METHODS;

ALGORITHM; COVARIANCE ANALYSIS; DATA ASSIMILATION; GRADIENT ANALYSIS; LEAST SQUARES METHOD; MATRIX; OCEANOGRAPHY; VECTOR;

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

References (74)

1. Axelsson O. 1996. Iterative Solution Methods. Cambridge University Press: Cambridge, UK.
2. Balmaseda MA, Mogensen K, Weaver AT. 2013. Evaluation of the ECMWF ocean reanalysis ORAS4. Q. J. R. Meteorol. Soc. DOI: 10.1002/qj.2063.
3. Bennett AF. 2002. Inverse Modelling of the Ocean and Atmosphere. Cambridge University Press: Cambridge, UK.
4. Broquet G, Edwards CA, Moore AM, Powell BS, Veneziani M, Doyle JD. 2009a. Application of 4D-variational data assimilation to the California Current System. Dyn. Atmos. Oceans 48: 69-91.
5. Broquet G, Moore AM, Arango HG, Edwards CA, Powell BS. 2009b. Ocean state and surface forcing correction using the ROMS-IS 4DVAR data assimilation system. Mercator Ocean Quarterly Newsletter 34: 5-13.
6. Buehner M. 2005. Ensemble-derived stationary and flow-dependent background-error covariances: evaluation in a quasi-operational NWP setting. Q. J. R. Meteorol. Soc. 131: 1013-1043.
7. Cardinali C, Pezzulli S, Andersson E. 2004. Influence-matrix diagnostic of a data assimilation system. Q. J. R. Meteorol. Soc. 130: 2767-2786.
8. Chan TF, Chow E, Saad Y, Yeung MC. 1999. Preserving symmetry in preconditioned Krylov subspace methods. SIAM J. Sci. Comput. 20: 568-581.
9. Chapman DC. 1985. Numerical treatment of cross-shelf open boundaries in a barotropic coastal ocean model. J. Phys. Oceanogr. 15: 1060-1075.
10. Checkley DM, Barth JA. 2009. Patterns and process in the California Current system. Prog. Oceanogr. 83: 49-64.
11. Chien CS, Chang SL. 2003. Application of the Lanczos algorithm for solving the linear systems that occur in continuation problems. Numer. Linear Algebra Appl. 10: 335-355. DOI: 10.1002/nla.306.
12. Chua BS, Bennett AF. 2001. An inverse ocean modeling system. Ocean Model. 3: 137-165.
13. Cohn S, Da Silva A, Guo J, Sienkiewicz M, Lamich D. 1998. Assessing the effects of data selection with the DAO physical-space statistical analysis system. Mon. Weather. Rev. 126: 2913-2926.
14. Courtier P. 1997. Dual formulation of four-dimensional variational assimilation. Q. J. R. Meteorol. Soc. 123: 2449-2461.
15. Courtier P, Thépaut J-N, Hollingsworth A. 1994. A strategy for operational implementation of 4D-Var using an incremental approach. Q. J. R. Meteorol. Soc. 120: 1367-1388.
16. Da Silva A, Pfaendtner J, Guo J, Sienkiewicz M, Cohn S. 1995. Assessing the effects of data selection with the DAO physical-space statistical analysis system. In: Proceedings of the 2nd WMO Symposium on assimilation of observations in meteorology and oceanography. World Meteorological Organization: Geneva, Switzerland. 273-278.
17. Daley R, Barker E. 2001. NAVDAS: Formulation and diagnostics. Mon. Weather Rev. 120: 869-883.
18. Derber J, Rosati A. 1989. A global oceanic data assimilation system. J. Phys. Oceanogr. 19: 1333-1347.
19. Desroziers G, Berre L. 2012. Accelerating and parallelizing minimizations in ensemble and deterministic variational assimilations. Q. J. R. Meteorol. Soc. DOI: 10.1002/qj.1886 In press.
20. Doyle JD, Jiang Q, Chao Y, Farrara J. 2009. High-resolution atmospheric modeling over the Monterey Bay during AOSN II. Deep Sea Res. II 56: 87-99.
21. Egbert GD, Bennett AF, Foreman MCG. 1994. TOPEX/POSEIDON tides estimated using a global inverse method. J. Geophys. Res. 99: 24821-24852.
22. 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.
23. 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.
24. El Akkraoui A, Trémolet Y, Todling R. 2012. Preconditioning of variational data assimilation and the use of a bi-conjugate gradient method. Q. J. R. Meteorol. Soc. DOI:10.1002/qj.1997. In press.
25. Evensen G. 2009. Data Assimilation: The Ensemble Kalman Filter. Springer Verlag: Berlin.
26. Fairall CW, Bradley EF, Godfrey JS, Wick GA, Ebson JB, Young GS. 1996a. Cool-skin and warm layer effects on the sea surface temperature. J. Geophys. Res. 101: 1295-1308.
27. Fairall CW, Bradley EF, Rogers DP, Ebson JB, Young GS. 1996b. Bulk parameterization of air-sea fluxes for tropical ocean global atmosphere coupled-ocean atmosphere response experiment. J. Geophys. Res. 101: 3747-3764.
28. Fisher M. 2003. Background-error covariance modelling. In: Proceedings of Seminar on Recent Developments in Data Assimilation for Atmosphere and Ocean. ECMWF: Reading, UK. 87-101.
29. Fisher M, Nocedal J, Trémolet Y, Wright SJ. 2009. Data assimilation in weather forecasting: A case study in PDE-constrained optimization. Optim. Eng. 10: 409-426.
30. Flather RA. 1976. A tidal model of the northwest European continental shelf. Mem. Soc. R. Sci. Liege 10: 141-164.
31. Gauthier P, Thépaut JN. 2001. Impact of the digital filter as a weak constraint in the pre-operational 4D-VAR assimilation system of Météo-France. Mon. Weather Rev. 129: 2089-2102.
32. Gelaro R, Zhu Y. 2009. Examination of observation impacts derived from observing system experiments (OSEs) and adjoint models. Tellus 61A: 179-193.
33. Golub GH, Van Loan CF. 1996. Matrix Computations. (third ed.) Johns Hopkins University Press: Baltimore, MA.
34. Gratton S, Tshimanga J. 2009. An observation-space formulation of variational assimilation using a Restricted Preconditioned Conjugate-Gradient algorithm. Q. J. R. Meteorol. Soc. 135: 1573-1585.
35. Gratton S, Lawless A, Nichols NK. 2007. Approximate Gauss-Newton methods for nonlinear least-squares problems. SIAM J. Optim. 18: 106-132.
36. Gratton S, Toint PL, Tshimanga J. 2009. 'Inexact range-space Krylov solvers for linear systems arising from inverse problems'. Technical Report 09/20, Department of Mathematics, FUNDP, University of Namur: Belgium.
37. Gratton S, Laloyaux P, Sartenaer A, Tshimanga J. 2011a. A reduced- and limited-memory preconditioned approach for the 4D-Var data-assimilation problem. Q. J. R. Meteorol. Soc. 137: 452-466.
38. Gratton S, Sartenaer A, Tshimanga J. 2011b. On a class of limited-memory preconditioners for large-scale linear systems with multiple right-hand sides. SIAM J. Optim. 21: 912-935.
39. Gratton S, Gürol S, Toint PL. 2013. Preconditioning and globalizing conjugate gradients in dual space for quadratically penalized nonlinear least-squares problems. Comput. Optim. Appl. 54: 1-25.
40. Haidvogel DB, Arango HG, Budgell WP, Cornuelle BD, Curchitser E, Di Lorenzo E, Fennel K, Geyer WR, Hermann AJ, Lanerolle L, Levin J, McWilliams JC, Miller AJ, Moore AM, Powell TM, Shchepetkin AF, Sherwood CR, Signell RP, Warner JC, Wilkin J. 2008. Ocean forecasting in terrain-following coordinates: Formulation and skill assessment of the Regional Ocean Modeling System. J. Comput. Phys. 227: 3595-3624.
41. Hickey BM. 1998. Coastal oceanography of western North America from the tip of Baja California to Vancouver Island. The Sea 11: 345-393.
42. Horn RA, Johnson CR. 1999. Matrix Analysis. Cambridge University Press: Cambridge, UK.
43. Lawless AS, Gratton S, Nichols NK. 2005. Approximate iterative methods for variational data assimilation. Int. J. Numer. Meth. Fluids. 47: 1129-1135.
44. Liu WT, Katsaros KB, Businger JA. 1979. Bulk parameterization of the air-sea exchange of heat and water vapor including the molecular constraints at the interface. J. Atmos. Sci. 36: 1722-1735.
45. Lorenc AC. 1988. Optimal nonlinear objective analysis. Q. J. R. Meteorol. Soc. 114: 205-240.
46. Lorenc AC. 1997. Development of an operational variational assimilation scheme. J. Meteorol. Soc. Japan. 75: 339-346.
47. Madec G. 2008. 'NEMO ocean engine'. Technical Report 27, Note du Pôle de modélisation, Institut Pierre-Simon Laplace (IPSL), France. http://www.nemo-ocean.eu.
48. Marchesiello P, McWilliams JC, Shchepetkin AF. 2001. Open boundary conditions for long-term integration of regional oceanic models. Ocean Model. 3: 1-20.
49. Mogensen K, Balmaseda MA, Weaver AT, Martin M, Vidard A. 2009. 'NEMOVAR: A variational data assimilation system for the NEMO model'. ECMWF Newsletter 120: 17-22.
50. Mogensen K, Balmaseda MA, Weaver AT. 2012. 'The NEMOVAR ocean data assimilation system as implemented in the ECMWF ocean analysis for System 4'. Tech. Report 668, ECMWF: Reading, UK.
51. Moore AM, Arango HG, Broquet G, Powell BS, Zavala-Garay J, Weaver AT. 2011a. The Regional Ocean Modeling System (ROMS) 4-dimensional variational data assimilation systems. Part I: System overview and formulation. Prog. Oceanogr. 91: 34-49.
52. Moore AM, Arango HG, Broquet G, Edwards C, Veneziani M, Powell B, Foley D, Costa D, Robinson P. 2011b. The Regional Ocean Modeling system (ROMS) 4-dimensional variational data assimilation systems. Part II: Performance and application to the California Current System. Prog. Oceanogr. 91: 50-73.
53. Moore AM, Arango HG, Broquet G, Edwards C, Veneziani M, Powell B, Foley D, Doyle JD, Costa D, Robinson P. 2011c. The Regional Ocean Modeling system (ROMS) 4-dimensional variational data assimilation systems. Part III: Observation impact and observation sensitivity in the California Current System. Prog. Oceanogr. 91: 74-94.
54. Nocedal J, Wright SJ. 2006. Numerical Optimization. Series in Operations Research, Springer Verlag: Berlin.
55. Nour-Omid B, Parlett BN, Raefsky A. 1988. Comparison of Lanczos with conjugate gradient using element preconditioning. In: Proceedings of First International Symposium on Domain Decomposition Methods for Partial Differential Equations. SIAM: Philadelphia, PA. 250-260.
56. Paige CC, Saunders MA. 1975. Solution of sparse indefinite systems of linear equations. SIAM J. Numer. Anal. 12: 617-629.
57. Papadrakakis M, Smerou S. 1990. A new implementation of the Lanczos method in linear problems. Int. J. Numer. Meth. Eng. 29: 141-159. DOI: 10.1002/nme.1620290110.
58. Parlett BN. 1980. The symmetric eigenvalue problem. Prentice-Hall: Englewood Cliffs, NJ, USA. Reprinted (1998) as Classics in Applied Mathematics 20: SIAM: Philadelphia, PA.
59. Purser RJ, Wu WW, Parrish DF, Roberts NM. 2003. Numerical aspects of the application of recursive filters to variational statistical analysis. Part I: Spatially homogeneous and isotropic Gaussian covariances. Mon. Weather Rev. 131: 1524-1535.
60. Saad Y. 1996. Iterative Methods for Sparse Linear Systems. PWS Publishing Company: Boston, MA.
61. Sasaki YK. 1970. Some basic formalisms in numerical variational analysis. Mon. Weather. Rev. 98: 875-883.
62. Shchepetkin AF, McWilliams JC. 2003. A method for computing horizontal pressure-gradient force in an oceanic model with a nonaligned vertical grid. J. Geophys. Res. 108C3 DOI: 10.1029/2001/JC001047.
63. Shchepetkin AF, McWilliams JC. 2005. The Regional Oceanic Modeling System (ROMS): a split explicit, free-surface, topography-following-coordinate oceanic model. Ocean Model. 9: 347-404.
64. Tarantola A. 2005. Inverse Problem Theory. Methods for Data Fitting and Model Parameter Estimation. SIAM: Philadelphia, PA.
65. Trémolet Y. 2006. Accounting for an imperfect model in 4D-Var. Q. J. R. Meteorol. Soc. 132: 2483-2504.
66. Tshimanga J, Gratton S, Weaver AT, Sartenaer A. 2008. Limited-memory preconditioners with application to incremental four-dimensional variational data assimilation. Q. J. R. Meteorol. Soc. 134: 751-769.
67. Veneziani M, Edwards CA, Doyle JD, Foley D. 2009a. A central California coastal ocean modeling study: 1. Forward model and the influence of realistic versus climatological forcing. J. Geophys. Res. 114C04015 DOI: 10.1029/2008JC004774.
68. Veneziani M, Edwards CA, Moore AM. 2009b. A central California coastal ocean modeling study: 2. Adjoint sensitivities to local and remote forcing mechanisms. J. Geophys. Res. 114C04020 DOI: 10.1029/2008JC004775.
69. Wang X, Synder C, Hamill TM. 2007. On the theoretical equivalence of differently proposed ensemble-3DVAR hybrid analysis schemes. Mon. Weather Rev. 135: 222-227.
70. Wang X, Barker DM, Synder C, Hamill TM. 2008. A hybrid ETKF-3DVAR data assimilation scheme for the WRF model. Part I: Observing system simulation experiment. Mon. Weather Rev. 136: 5116-5131.
71. Weaver AT, Mirouze I. 2013. On the diffusion equation and its application to isotropic and anisotropic correlation modelling in variational assimilation. Q. J. R. Meteorol. Soc. 139: 242-260.
72. Weaver AT, Deltel C, Machu E, Ricci S, Daget N. 2005. A multivariate balance operator for variational ocean data assimilation. Q. J. R. Meteorol4. Soc. 131: 3605-3625.
73. Weaver AT, Vialard J, Anderson DLT. 2003. Three- and four-dimensional variational assimilation with a general circulation model of the tropical Pacific ocean. Part 1: Formulation, internal diagnostics and consistency checks. Mon. Weather Rev. 131: 1360-1378.
74. Wunsch C, Heimbach P. 2007. Practical global ocean state estimation. Physica D 230: 197-208.