Reuse & Permissions

It is not necessary to obtain permission to reuse this article or its components as it is available under the terms of the Creative Commons Attribution 4.0 International license. This license permits unrestricted use, distribution, and reproduction in any medium, provided attribution to the author(s) and the published article's title, journal citation, and DOI are maintained. Please note that some figures may have been included with permission from other third parties. It is your responsibility to obtain the proper permission from the rights holder directly for these figures.

Export citation

Export citation

Choose format for download:

Download Citation
  • Open Access

Consecutive shooting during adjoint optimization enables highly accurate estimation of unobserved state variables of high-dimensional chaotic systems

Kota Shiozawa1,*, Inga Kottlarz2,3,4,5,*, Isao T. Tokuda1, Ulrich Parlitz2,3,4,†, and Sebastian Herzog4,6

  • *These authors contributed equally to this work.
  • †Contact author: ulrich.parlitz@ds.mpg.de

Phys. Rev. Research 8, 013014 – Published 9 January, 2026

DOI: https://doi.org/10.1103/8c7w-dpdt

Abstract

Recent advancements in machine learning have led to the development of sophisticated algorithms for automatic differentiation. These software tools simplify the implementation of the adjoint method by automating the derivation of adjoint equations, enabling researchers to apply the adjoint method to a wide range of problems without requiring extensive expertise in the underlying mathematics. One particularly significant application is the state estimation of partially observable systems. We show that for this application the performance of the adjoint method can be significantly improved by a shooting protocol in which the end of the previous segment is used as the initial condition for estimating the next segment, independent of the success of the estimate for the previous segment. The superior performance of this adjoint optimization-based consecutive shooting method is demonstrated using Lorenz-96 models with dimensions from 9 to 300.

View figure in article

Physics Subject Headings (PhySH)

Article Text

References (49)

  1. A. Abur and A. G. Exposito, Power System State Estimation: Theory and Implementation (CRC Press, Boca Raton, FL, 2004).
  2. F. L. Chernousko, State Estimation for Dynamic Systems (CRC Press, Boca Raton, FL, 1993).
  3. R. E. Kalman, A new approach to linear filtering and prediction problems, J. Basic Eng. 82, 35 (1960).
  4. H. D. I. Abarbanel, The Statistical Physics of Data Assimilation and Machine Learning (Cambridge University Press, Cambridge, 2022).
  5. A. Isidori, Nonlinear Control Systems: An Introduction (Springer, Berlin, Heidelberg, 1985).
  6. H. Nijmeijer and A. Van der Schaft, Nonlinear Dynamical Control Systems (Springer, New York, NY, 1990), Vol. 175.
  7. E. Ott, B. R. Hunt, I. Szunyogh, A. V. Zimin, E. J. Kostelich, M. Corazza, E. Kalnay, D. Patil, and J. A. Yorke, A local ensemble Kalman filter for atmospheric data assimilation, Tellus A: Dyn. Meteorol. Oceanogr. 56, 415 (2004).
  8. G. Datseris and U. Parlitz, Nonlinear Dynamics: A Concise Introduction Interlaced with Code (Springer Nature, Cham, Switzerland, 2022).
  9. G. Evensen, Data Assimilation (Springer, Berlin, 2009).
  10. L. Junge and U. Parlitz, Synchronization and control of coupled Ginzburg-Landau equations using local coupling, Phys. Rev. E 61, 3736 (2000).
  11. S. Berg, S. Luther, and U. Parlitz, Synchronization based system identification of an extended excitable system, Chaos 21, 033104 (2011),.
  12. E. Baake, M. Baake, H. G. Bock, and K. M. Briggs, Fitting ordinary differential equations to chaotic data, Phys. Rev. A 45, 5524 (1992).
  13. H. D. Abarbanel, D. R. Creveling, R. Farsian, and M. Kostuk, Dynamical state and parameter estimation, SIAM J. Appl. Dyn. Syst. 8, 1341 (2009).
  14. J. Ye, N. Kadakia, P. J. Rozdeba, H. D. I. Abarbanel, and J. C. Quinn, Improved variational methods in statistical data assimilation, Nonlinear Processes Geophys. 22, 205 (2015).
  15. J. R. Stroud and T. Bengtsson, Sequential state and variance estimation within the ensemble Kalman filter, Mon. Weather Rev. 135, 3194 (2007).
  16. D. Rey, M. Eldridge, M. Kostuk, H. D. Abarbanel, J. Schumann-Bischoff, and U. Parlitz, Accurate state and parameter estimation in nonlinear systems with sparse observations, Phys. Lett. A 378, 869 (2014).
  17. D. Rey, M. Eldridge, U. Morone, H. D. I. Abarbanel, U. Parlitz, and J. Schumann-Bischoff, Using waveform information in nonlinear data assimilation, Phys. Rev. E 90, 062916 (2014).
  18. J. Schumann-Bischoff, U. Parlitz, H. D. I. Abarbanel, M. Kostuk, D. Rey, M. Eldridge, and S. Luther, Basin structure of optimization based state and parameter estimation, Chaos 25, 053108 (2015).
  19. K. J. H. Law, D. Sanz-Alonso, A. Shukla, and A. M. Stuart, Filter accuracy for the Lorenz 96 model: Fixed versus adaptive observation operators, Physica D 325, 1 (2016).
  20. S. G. Penny, T. A. Smith, T. Chen, J. A. Platt, H. Lin, M. Goodliff, and H. D. I. Abarbanel, Integrating recurrent neural networks with data assimilation for scalable data-driven state estimation, J. Adv. Model. Earth Syst. 14, e2021MS002843 (2022).
  21. J. Brajard, A. Carrassi, M. Bocquet, and L. Bertino, Combining data assimilation and machine learning to emulate a dynamical model from sparse and noisy observations: A case study with the Lorenz 96 model, J. Comput. Sci. 44, 101171 (2020).
  22. J. Guillot, P. Ailliot, E. Frénod, J. Ruiz, and P. Tandeo, State and stochastic parameters estimation with combined ensemble Kalman and particle filters, Mon. Weather Rev. 153, 1141 (2025).
  23. R. Fablet, B. Chapron, L. Drumetz, E. Mémin, O. Pannekoucke, and F. Rousseau, Learning variational data assimilation models and solvers, J. Adv. Model. Earth Syst. 13, e2021MS002572 (2021).
  24. A. Amemiya, M. Shlok, and T. Miyoshi, Application of recurrent neural networks to model bias correction: Idealized experiments with the Lorenz-96 Model, J. Adv. Model. Earth Syst. 15, e2022MS003164 (2023).
  25. T. J. Cocucci, M. Pulido, M. Lucini, and P. Tandeo, Model error covariance estimation in particle and ensemble Kalman filters using an online expectation–maximization algorithm, Q. J. R. Meteorol. Soc. 147, 526 (2021).
  26. L. Lettermann, A. Jurado, T. Betz, F. Wörgötter, and S. Herzog, Tutorial: A beginner’s guide to building a representative model of dynamical systems using the adjoint method, Commun. Phys. 7, 128 (2024).
  27. E. Lorenz, Predictability: A problem partly solved, in Seminar on Predictability, 4-8 September 1995 (ECMWF, Shinfield Park, Reading, 1995), Vol. 1, pp. 1–18.
  28. I. Kottlarz and U. Parlitz, Ordinal pattern-based complexity analysis of high-dimensional chaotic time series, Chaos 33, 053105 (2023).
  29. The Lyapunov spectra were calculated using the DynamicalSystems.jl library [8, 46], which implements the method “H2” in Ref. [30], originally stated in Ref. [31].
  30. K. Geist, U. Parlitz, and W. Lauterborn, Comparison of different methods for computing Lyapunov exponents, Prog. Theor. Phys. 83, 875 (1990).
  31. G. Benettin, L. Galgani, A. Giorgilli, and J.-M. Strelcyn, Lyapunov characteristic exponents for smooth dynamical systems and for Hamiltonian systems; A method for computing all of them. Part 1: Theory, Meccanica 15, 9 (1980).
  32. J. Bradbury, R. Frostig, P. Hawkins, M. J. Johnson, C. Leary, D. Maclaurin, G. Necula, A. Paszke, J. VanderPlas, S. Wanderman-Milne, and Q. Zhang, JAX: Composable transformations of Python+NumPy programs (2018), http://github.com/jax-ml/jax.
  33. J. S. Arora and E. J. Haug, Methods of design sensitivity analysis in structural optimization, AIAA J. 17, 970 (1979).
  34. A. D. Belegundu, Lagrangian approach to design sensitivity analysis, J. Eng. Mech. 111, 680 (1985).
  35. G. Bianconi, A. Arenas, J. Biamonte, L. D. Carr, B. Kahng, J. Kertesz, J. Kurths, L. Lü, C. Masoller, A. E. Motter, M. Perc, F. Radicchi, R. Ramaswamy, F. A. Rodrigues, M. Sales-Pardo, M. S. Miguel, S. Thurner, and T. Yasseri, Complex systems in the spotlight: Next steps after the 2021 Nobel prize in physics, J. Phys.: Complex. 4, 010201 (2023).
  36. R. M. Errico, What is an adjoint model? Bull. Am. Meteorol. Soc. 78, 2577 (1997).
  37. P. Ji, J. Ye, Y. Mu, W. Lin, Y. Tian, C. Hens, M. Perc, Y. Tang, J. Sun, and J. Kurths, Signal propagation in complex networks, Phys. Rep. 1017, 1 (2023).
  38. S. A. Jurovics and J. E. McIntyre, The adjoint method and its application to trajectory optimization, ARS J. 32, 1354 (1962).
  39. M. Lotkin and H. N. Browne, On the accuracy of the adjoint method of differential corrections, Am. Math. Monthly 63, 97 (1956),.
  40. R.-E. Plessix, A review of the adjoint-state method for computing the gradient of a functional with geophysical applications, Geophys. J. Int. 167, 495 (2006).
  41. R. Tomović, Sensitivity analysis of dynamic systems, Electronic Sciences Series (McGraw-Hill, New York, 1963).
  42. R. Tomović, The role of sensitivity analysis in engineering problems, in Sensitivity Methods in Control Theory, edited by L. Radanovic (Pergamon, Oxford, 1966), pp. 103–109.
  43. J. Schumann-Bischoff, S. Luther, and U. Parlitz, Nonlinear system identification employing automatic differentiation, Commun. Nonlinear Sci. Numer. Simul. 18, 2733 (2013).
  44. https://gitlab.gwdg.de/sherzog3/adoptODE.
  45. https://github.com/KotaShiozawa/adoptODE-lorenz96.
  46. G. Datseris, DynamicalSystems.jl: A Julia software library for chaos and nonlinear dynamics, J. Open Source Software 3, 598 (2018).
  47. Y. Kuramoto, Diffusion-induced chaos in reaction systems, Prog. Theor. Phys. Suppl. 64, 346 (1978).
  48. G. I. Sivashinsky, On flame propagation under conditions of stoichiometry, SIAM J. Appl. Math. 39, 67 (1980).
  49. A.-K. Kassam and L. N. Trefethen, Fourth-order time-stepping for stiff PDEs, SIAM J. Sci. Comput. 26, 1214 (2005).

Outline

Information

Sign In to Your Journals Account

Filter

Filter

Article Lookup

Enter a citation