A Method for Adjoint Variational Data Assimilation with Physical “On–Off” Processes

Mu Mu LASG, Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing, China

Search for other papers by Mu Mu in
Current site
Google Scholar
PubMed
Close
and
Jiafeng Wang LASG, Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing, China

Search for other papers by Jiafeng Wang in
Current site
Google Scholar
PubMed
Close
Restricted access

Abstract

Since the accuracy of the tangent linear approximation of moist physics in a mesoscale model is case dependent, the problem related to the variational data assimilation with physical “on–off” processes is studied further in both time-continuous and discrete circumstances. Two kinds of typical on–off switches represented in idealized simple models are investigated: the zero-order discontinuous switch (Type I) and the first-order discontinuous switch (Type II). The main results are as follows. For Type I: 1) For the case in which the model is time continuous, the gradient of the cost function with respect to the initial condition exists except for the threshold. 2) In the time-discrete case, there are zigzag discontinuities in the cost function, and the method that keeps the switches in the tangent linear model the same as in the forward model (called Zou's method) is able to compute the correct gradient where it exists. An optimization with this gradient might yield a local minimum rather than the global minimum if the cost function has multiple minima, however. 3) A method based on the nonlinear perturbation equation is proposed that can give the accurate gradient in the time-continuous case. 4) In the discrete case, the method of this paper is useful to obtain the global descent direction of the cost function in optimization and is helpful to find the global minimum. In addition, it still employs the adjoint model constructed by Zou's method. For Type II, Zou's method can be used for both the time-continuous and discrete cases. The importance of reducing the model error in the context of variational data assimilation with discontinuous physics is also indicated.

Corresponding author address: Prof. Mu Mu, LASG, P.O. Box 9804, Beijing, 100029, China. Email: mumu@lasg.iap.ac.cn

Abstract

Since the accuracy of the tangent linear approximation of moist physics in a mesoscale model is case dependent, the problem related to the variational data assimilation with physical “on–off” processes is studied further in both time-continuous and discrete circumstances. Two kinds of typical on–off switches represented in idealized simple models are investigated: the zero-order discontinuous switch (Type I) and the first-order discontinuous switch (Type II). The main results are as follows. For Type I: 1) For the case in which the model is time continuous, the gradient of the cost function with respect to the initial condition exists except for the threshold. 2) In the time-discrete case, there are zigzag discontinuities in the cost function, and the method that keeps the switches in the tangent linear model the same as in the forward model (called Zou's method) is able to compute the correct gradient where it exists. An optimization with this gradient might yield a local minimum rather than the global minimum if the cost function has multiple minima, however. 3) A method based on the nonlinear perturbation equation is proposed that can give the accurate gradient in the time-continuous case. 4) In the discrete case, the method of this paper is useful to obtain the global descent direction of the cost function in optimization and is helpful to find the global minimum. In addition, it still employs the adjoint model constructed by Zou's method. For Type II, Zou's method can be used for both the time-continuous and discrete cases. The importance of reducing the model error in the context of variational data assimilation with discontinuous physics is also indicated.

Corresponding author address: Prof. Mu Mu, LASG, P.O. Box 9804, Beijing, 100029, China. Email: mumu@lasg.iap.ac.cn

Save
  • Errico, R. M., and K. D. Raeder, 1999: An examination of the accuracy of the linearization of a mesoscale model with moist physics. Quart. J. Roy. Meteor. Soc., 125 , 169195.

    • Search Google Scholar
    • Export Citation
  • Gilbert, J. C., and C. Lemarechal, 1989: Some numerical experiments with variable-storage quasi-Newton algorithms. Math. Prog., 45 , 407435.

    • Search Google Scholar
    • Export Citation
  • Janisková, M., J-F. Thepaut, and J-F. Geleyn, 1999: Simplified and regular physical parameterizations for incremental four-dimensional variational assimilation. Mon. Wea. Rev., 127 , 2645.

    • Search Google Scholar
    • Export Citation
  • Le Dimet, F-X., and O. Talagrand, 1986: Variational algorithm for analysis and assimilation of meteorological observations: Theoretical aspects. Tellus, 38A , 97110.

    • Search Google Scholar
    • Export Citation
  • Lu, C., and G. L. Browning, 1998: The impact of observational and model errors on four-dimensional variational data assimilation. J. Atmos. Sci., 55 , 9951011.

    • Search Google Scholar
    • Export Citation
  • Lu, C., and G. L. Browning, 2000: Discontinuous forcing generating rough initial conditions in 4DVAR data assimilation. J. Atmos. Sci., 57 , 16461656.

    • Search Google Scholar
    • Export Citation
  • Mahfouf, J-F., and F. Rabier, 2000: The ECMWF operational implementation of four-dimensional variational assimilation. Part II: Experimental results with improved physics. Quart. J. Roy. Meteor. Soc., 126 , 11711190.

    • Search Google Scholar
    • Export Citation
  • Navon, I. M., X. Zou, J. Derber, and J. Sela, 1992: Variational data assimilation with an adiabatic version of the NMC spectral model. Mon. Wea. Rev., 120 , 14331446.

    • Search Google Scholar
    • Export Citation
  • Rostaing, N., S. Dalmas, and A. Galligo, 1993: Automatic differentiation in Odyssee. Tellus, 45A , 558568.

  • Verlinde, J., and W. R. Cotton, 1993: Fitting microphysical observations of nonsteady convective clouds to a numerical model: An application of the adjoint technique of data assimilation to a kinematic model. Mon. Wea. Rev., 121 , 27762793.

    • Search Google Scholar
    • Export Citation
  • Xu, Q., 1996: Generalized adjoint for physical processes with parameterized discontinuities. Part I: Basic issues and heuristic examples. J. Atmos. Sci., 53 , 11231142.

    • Search Google Scholar
    • Export Citation
  • Xu, Q., 1997: Generalized adjoint for physical processes with parameterized discontinuities. Part IV: Problems and treatment in time discretization. J. Atmos. Sci., 54 , 27222728.

    • Search Google Scholar
    • Export Citation
  • Xu, Q., and J. Gao, 1999: Generalized adjoint for physical processes with parameterized discontinuities. Part VI: Minimization problems in multidimensional space. J. Atmos. Sci., 56 , 9941002.

    • Search Google Scholar
    • Export Citation
  • Xu, Q., J. Gao, and W. Gu, 1998: Generalized adjoint for physical processes with parameterized discontinuities. Part V: Coarse-grain adjoint and problems in gradient check. J. Atmos. Sci., 55 , 21302135.

    • Search Google Scholar
    • Export Citation
  • Zou, X., 1997: Tangent linear and adjoint of “on–off” processes and their feasibility for use in 4-dimensional variational data assimilation. Tellus, 49A , 331.

    • Search Google Scholar
    • Export Citation
  • Zou, X., I. M. Navon, and J. G. Sela, 1993: Variational data assimilation with moist threshold processes using the NMC spectral model. Tellus, 45A , 370387.

    • Search Google Scholar
    • Export Citation
All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 277 74 20
PDF Downloads 97 37 0