Summary
The adjoint method application in variational data assimilation provides a way of obtaining the exact gradient of the cost functionj with respect to the control variables. Additional information may be obtained by using second order information. This paper presents a second order adjoint model (SOA) for a shallow-water equation model on a limited-area domain. One integration of such a model yields a value of the Hessian (the matrix of second partial derivatives, ∇2 J) multiplied by a vector or a column of the Hessian of the cost function with respect to the initial conditions. The SOA model was then used to conduct a sensitivity analysis of the cost function with respect to distributed observations and to study the evolution of the condition number (the ratio of the largest to smallest eigenvalues) of the Hessian during the course of the minimization. The condition number is strongly related to the convergence rate of the minimization. It is proved that the Hessian is positive definite during the process of the minimization, which in turn proves the uniqueness of the optimal solution for the test problem. Numerical results show that the sensitivity of the response increases with time and that the sensitivity to the geopotential field is larger by an order of magnitude than that to theu andv components of the velocity field. Experiments using data from an ECMWF analysis of the First Global Geophysical Experiment (FGGE) show that the cost functionJ is more sensitive to observations at points where meteorologically intensive events occur. Using the second order adjoint shows that most changes in the value of the condition number of the Hessian occur during the first few iterations of the minimization and are strongly correlated to major large-scale changes in the reconstructed initial conditions fields.
Similar content being viewed by others
References
Cacuci, Dan G., 1981: Sensitivity theory for nonlinear systems. II. Extensions to additional classes of responses.J. Math. Phy.,22(12), 2803–2812.
Courtier, P., Talagrand, O., 1990: Variational assimilation of meteorological observations with the direct and adjoint shallow-water equations.Tellus,42A, 531–549.
Chao, W. C., Chang, L., 1992: Development of a 4-dimensional analysis system using the adjoint method at GLA. Part 1: Dynamics.Mon. Wea. Rev. 120, 1661–1673.
Courtier, P., 1985: Experiments in data assimilation using the adjoint model technique. Workshop on High-Resolution Analysis ECMWF (UK) June 1985.
Courtier, P., Talagrand, O., 1987: Variational assimilation of meteorological observations with the adjoint equations Part 2. Numerical results.Quart. J. Roy. Meteor. Soc.,113, 1329–1347.
Derber, J. C., 1985: The variational four dimensional assimilation of analysis using filtered models as constraints. Ph.D. Thesis, Univ. of Wisconsin-Madison, 141 pp.
Derber, J. C., 1987: Variational four dimensional analysis using the quasigeostrophic constraint.Mon. Wea. Rev.,115, 998–1008.
Fletcher, R., 1987:Practical Methods of Optimization, Vol. 1, 2nd ed. Chichester, New York: John Wiley, 436 pp.
Gill, P. E., Murray, W., Wright, M. H., 1981:Practical Optimization. London, Academic Press, 401 pp.
Grammeltvedt, A., 1969: A survey of finite-difference schemes for the primitive equations for a barotropic fluid.Mon. Wea. Rev.,97, 387–404.
Hall, M. C. G., Cacuci, D. G., Schlesinger, M. E., 1982: Sensitivity analysis of a radiative-convective model by the adjoint method.J. Atmos. Sci.,39, 2038–2050.
Hall, M. C. G., Cacuci, D. G., 1983: Physical interpretation of the adjoint functions for sensitivity analysis of atmospheric models.J. Atmos. Sci.,40, 2537–2546.
Hoffmann, R. N., 1986: A four dimensional analysis exactly satisfying equations of motion.Mon. Wea. Rev.,114, 388–397.
Kontarev, G., 1980: The adjoint equation technique applied to meteorological problem.EC MWF Tech. Rep., Reading, U. K., 21, 1–21 pp.
Le Dimet, F. X., 1982: A general formalism of variational analysis. CIMMS Report, Norman, OK 73091, 22, 1–34 pp.
Le Dimet, F. X., Talagrand, O., 1986: Variational algorithms for analysis and assimilation of meteorological observations: Theoretical aspects.Tellus,38A, 97–110.
Le Dimet, F. X., Nouailler, A., 1986: In: Sasaki, Y. K. (ed.)Assimilation of Dynamics in Geosciences. Amsterdam: Elsevier, 1986, 181–198.
Lewis, J. M., Derber, J. C., 1985: The use of adjoint equations to solve a variational adjustment problem with advective constraints.Tellus,37A, 309–322.
Liu, D. C., Jorge Nocedal, 1989: On the limited memory BFGS method for large scale minimization.Mathematical Programming,45, 503–528.
Lions, J. L., 1971:Optimal Control of Systems Governed by Partial Differential Equations (Translated by Mitter, S.K.). Berlin, Heidelberg: Springer, 396 pp
Lorenc, A. C., 1988a: Optimal nonlinear objective analysis.Quart. J. Roy. Meteor. Soc.,114, 205–240.
Lorenc, A. C., 1988b: A practical approximation to optimal four dimensional objective analysis.Mon. Wea. Rev.,116, 730–745.
Luenberger, D. G., 1984:Linear and Nonlinear Programming, 2nd ed. Reading, MA: Addison-Wesley, 491 pp.
Marchuk, G. I., 1974:Numerical Solution of the Problem of Dynamics of Atmosphere and Ocean (in Russian). Leningrad: Gidrometeoizdat, 303 pp.
Nash, S. G., 1985: Preconditioning of truncated-Newton methods.SIAM J. Sci. Stat. Comput.,6(3), 599–616.
Navon, I. M., Zou, X. L., Johnson, K., Derber, J., Sela, J., 1990: The adjoint of NMC spectral model. Tech. Report, FSU, Tallahassee, FSU-SCRI-REPORT-90-22, 71 pp.
Navon, I. M., Zou, X. L., Derber, J., Sela, J., 1992: Variational data assimilation with an adiabatic version of the NMC spectral model.Mon. Wea. Rev.,120, 1433–1466.
Penenko, V., Obraztsov, N. N., 1976: A variational initialization method for the fields of the meteorological elements.Meteorol. Gidrol. (English translation),11, 1–11.
Strang, G. 1986:Introduction to Applied Mathematics. Wellesley: Wellesley-Cambridge Press, 758 pp.
Sykes, J. F., Wilson, J. L., R. W. Andrews, 1985: Sensitivity analysis for steady state groundwater flow using adjoint operators.Water Resources Research,21 (3), 359–371.
Talagrand, O., Courtier, P., 1987: Variational assimilation of meteorological observations with the adjoint vorticity equation-Part 1. Theory.Quart. J. Roy. Meteor. Soc.,113, 1311–1328.
Thacker, W. C., Long, R. B., 1988: Fitting dynamics to data.J. Geophys. Res.,93, 1227–1240.
Thacker, W. C., 1989: The Role of Hessian Matrix in Fitting models to Measurements.J. Geophys. Res.,94, 6177–6196.
Thépaut, J. N., Courtier, 1991: Four-dimensional variational data assimilation using the adjoint of multilevel primitive equation model.Quart. J. Roy. Meteor. Soc.,117, 1225–1254.
Zou, X. L., Navon, I. M., Berger, M., Paul K. H. Phua, Schlick, T., Le Dimet, F. X., 1992: Numerical experience with limited-Memory Quasi-Newton methods and Truncated Newton methods.SIAM J. Optimization (in press).
Zou, X. L., Navon, I. M., Le Dimet, F. X., 1992: Incomplete observations and control of gravity waves in variational data assimilation.Tellus,44A, 273–296.
Zou, X. L., Navon, I. M., Sela, J., 1992: Control of gravity oscillations in variational data assimilation.Mon. Wea. Rev. (in press).
Author information
Authors and Affiliations
Additional information
With 17 Figures
Rights and permissions
About this article
Cite this article
Wang, Z., Navon, I.M., Le Dimet, F.X. et al. The second order adjoint analysis: Theory and applications. Meteorl. Atmos. Phys. 50, 3–20 (1992). https://doi.org/10.1007/BF01025501
Received:
Revised:
Issue Date:
DOI: https://doi.org/10.1007/BF01025501