Abstract
For dealing with dynamical instability in predictions, numerical models should be provided with accurate initial values on the attractor of the dynamical system they generate. A discrete control scheme is presented to this end for trailing variables of an evolutive system of ordinary differential equations. The Influence Sampling (IS) scheme adapts sample values of the trailing variables to input values of the determining variables in the attractor. The optimal IS scheme has affordable cost for large systems. In discrete data assimilation runs conducted with the Lorenz 1963 equations and a nonautonomous perturbation of the Lorenz equations whose dynamics shows on-off intermittency the optimal IS was compared to the straightforward insertion method and the Ensemble Kalman Filter (EnKF). With these unstable systems the optimal IS increases by one order of magnitude the maximum spacing between insertion times that the insertion method can handle and performs comparably to the EnKF when the EnKF converges. While the EnKF converges for sample sizes greater than or equal to 10, the optimal IS scheme does so fromsample size 1. This occurs because the optimal IS scheme stabilizes the individual paths of the Lorenz 1963 equations within data assimilation processes.
References
[1] A.J. Majda, Challenges in climate science and contemporary applied mathematics, Comm. Pure App. Math., LXV (2012) 920-948.10.1002/cpa.21401Search in Google Scholar
[2] O.A. Ladyzhenskaya, The Mathematical Theory of Viscous Incompressible Flow, Gordon and Breach Sci. Pub., 1963.Search in Google Scholar
[3] V.K. Kalantarov, E.S. Titi, Global attractors and determining modes for the 3D Navier-Stokes-Voight equations, Chin. Ann. Math., 30B(6) (2009) 697-714.10.1007/s11401-009-0205-3Search in Google Scholar
[4] Z. Zhang, T.N. Krishnamurti, A perturbation method for hurricane ensemble predictions, Mon. Wea. Rev., 127 (1999) 447-469.Search in Google Scholar
[5] E. Kalnay, Z. Toth, Ensemble forecasting at NMC: the generation of perturbations, Bull. Amer. Met. Soc., 74:12 (1993) 2317-2330.Search in Google Scholar
[6] R. Buizza, The singular vector approach to the analysis of perturbation growth in the atmosphere, Ph.D. Thesis, University of London, 1997.Search in Google Scholar
[7] J.S. Rosenthal, Optimal proposal distributions and adaptive MCMC, in: S. Brooks, A. Gelman, G.L. Jones, X.L. Meng (Eds.), Handbook of Markov Chain Monte Carlo, Chapman&Hall/CRC Handbooks of Modern Statistical Methods, 2011, Chap. 4. 10.1201/b10905-5Search in Google Scholar
[8] P.J. Huber, Robust Statistics, Wiley, 1981.10.1002/0471725250Search in Google Scholar
[9] F.R. Hampel, E.M. Ronchetti, P.J. Rousseeuw, W.A. Stahel, Robust Statistics: the Approach Based on Influence Functions, Wiley, 1986.Search in Google Scholar
[10] A.J. Chorin, O.H. Hald, R. Kupferman, Optimal prediction with memory, Physica D, 166 (2002) 239-257.Search in Google Scholar
[11] J. Charney, M. Halem, R. Jastrow, Use of incomplete historical data to infer the present state of the atmosphere, J. Atm. Sci., 26 (1969) 1160-1163.Search in Google Scholar
[12] K. Hayden, E. Olson, E.S. Titi, Discrete data assimilation in the Lorenz and 2D Navier-Stokes equations, Physica D, 240 (2011) 1416-1425.Search in Google Scholar
[13] K.J.H. Law, A.M. Stuart, K.C. Zygalakis, Data Assimilation: A Mathematical Introduction, Springer, 2015.10.1007/978-3-319-20325-6Search in Google Scholar
[14] F. John, Partial Differential Equations, Springer, 1982.10.1007/978-1-4684-9333-7Search in Google Scholar
[15] D.J. Evans, G.P. Morriss, Statistical Mechanics of Nonequilibrium Liquids, Academic Press, 1990.Search in Google Scholar
[16] E.N. Lorenz, Deterministic nonperiodic flow, J. Atm. Sci., 20 (1963) 130-141.Search in Google Scholar
[17] E. Ott, J.C. Sommerer, Blowout bifurcations: the occurrence of riddled basins and on-off intermittency, Phys. Lett. A, 188 (1994) 39-47.Search in Google Scholar
[18] S.C. Venkataramani, T.M. Antonsen Jr, E. Ott, J.C. Sommerer, On-off intermittency: power spectrum and fractal properties of time series, Physica D, 96 (1996) 66-99.Search in Google Scholar
[19] L.M. Pecora, T.L. Carroll, Synchronization in chaotic systems, Phys. Rev. Lett., 64:8 (1990) 821-824.Search in Google Scholar
[20] C.H. Bishop, B.J. Etherton, S.J. Majumdar, Adaptive sampling with the Ensemble Transform Kalman Filter. Part I: theoretical aspects, Mon. Wea. Rev., 129 (2001) 420-436.10.1175/1520-0493(2001)129<0420:ASWTET>2.0.CO;2Search in Google Scholar
[21] A. Farhat, E. Lunasin, E.S. Titi, On the Charney conjecture of data assimilation employing temperature measurements alone: the paradigm of 3D planetary geostrophic model, Math. Clim. Wea. Forecasting, 2 (2016) 61-74.Search in Google Scholar
[22] J.E. Hoke, R.A. Anthes, The initialization of numerical models by a dynamic-initialization technique, Mon. Wea. Rev., 104 (1976) 1551-1556.Search in Google Scholar
[23] G.H. Golub, J.M. Ortega, Scientific Computing and Differential Equations, Academic Press, 1992.10.1016/B978-0-08-051669-1.50005-8Search in Google Scholar
© 2017
This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 License.
