Data assimilation into nonlinear stochastic models

Article English OPEN
Miller, Robert N. ; Carter, Everett F. ; Blue, Sally T. (2011)

With very few exceptions, data assimilation methods which have been used or proposed foruse with ocean models have been based on some assumption of linearity or near-linearity. Thegreat majority of these schemes have at their root some least-squares assumption. While onecan always perform least-squares analysis on any problem, direct application of least squaresmay not yield satisfactory results in cases in which the underlying distributions are significantlynon-Gaussian. In many cases in which the behavior of the system is governed by intrinsicallynonlinear dynamics, distributions of solutions which are initially Gaussian will not remain soas the system evolves. The presence of noise is an additional and inevitable complicating factor.Besides the imperfections in our models which result from physical or computational simplifyingassumptions, there is uncertainty in forcing fields such as wind stress and heat flux which willremain with us for the foreseeable future. The real world is a noisy place, and the effects ofnoise upon highly nonlinear systems can be complex. We therefore consider the problem ofdata assimilation into systems modeled as nonlinear stochastic differential equations. When themodels are described in this way, the general assimilation problem becomes that of estimatingthe probability density function of the system conditioned on the observations. The quantitywe choose as the solution to the problem can be a mean, a median, a mode, or some otherstatistic. In the fully general formulation, no assumptions about moments or near-linearity arerequired. We present a series of simulation experiments in which we demonstrate assimilationof data into simple nonlinear models in which least-squares methods such as the (Extended)Kalman filter or the weak-constraint variational methods will not perform well. We illustratethe basic method with three examples: a simple one-dimensional nonlinear stochastic differentialequation, the well known three-dimensional Lorenz model and a nonlinear quasigeostrophicchannel model. Comparisons to the extended Kalman filter and an extension to the extendedKalman filter are presented.DOI: 10.1034/j.1600-0870.1999.t01-2-00002.x
  • References (47)
    47 references, page 1 of 5

    Bennett, A. F. and Thorburn, M. A. 1992. The generalized inverse of a nonlinear quasigeostrophic ocean circulation model. J. Phys. Oceanogr. 22, 213-230.

    Bennett, A. F., Leslie, L. M., Hagelberg, C. R. and Powers, P. E. 1993. Tropical cyclone prediction using a barotropic model initialized by a generalized inverse method. Mon. Wea. Rev. 121, 1714-1729.

    Bergman, L. A., Wojtkiewicz, S. F., Johnson, E. A. and Spencer, B. F. 1996. A robust numerical solution of the Fokker-Planck equation for second order dynamical systems under parametric and external white noise excitation. Fields Inst. Comm. 9, 23-37.

    Bergman, L. A. and Spencer, B. F. Jr. 1992. Robust numerical solution of the transient Fokker-Planck equation for nonlinear dynamical systems. IUTAM Symposium, Turin. N. Bellomo and F. Casciati, eds., Springer-Verlag, Berlin, Heidelberg, 49-59.

    Bierman, G. J. 1977. Factorization methods for discrete sequential estimation. Academic Press, New York. 241 pp.

    Bouttier, F. 1994. A dynamical estimation of forecast error covariances in an assimilation system. Mon. Wea. Rev. 122, 2376-2390.

    Box, G. E. P and Muller, M. E. 1958. A note on the generation of random normal deviates. Annals Math. Stat. 29, 610-611.

    Budgell, W. P. 1986. Nonlinear data assimilation for shallow water equations in branched channels. J. Geophys. Res. 91, 10633-10644.

    Chang, J. S. and Cooper, G. 1970. A practical diVerence scheme for Fokker-Planck equations. J. Comput. Phys. 6, 1-16.

    Charney, J. G. and Devore, J. G. 1979. Multiple flow equilibria in the atmosphere and blocking. J. Atmos. Sci. 36, 1205-1216.

  • Metrics
    No metrics available
Share - Bookmark