Projected shadowing-based data assimilation

Abstract

In this article we develop algorithms for data assimilation based upon a computational time dependent stable/unstable splitting. Our particular method is based upon shadowing refinement and synchronization techniques and is motivated by work on assimilation in the unstable subspace [Carrassi et al., Chaos, 18 (2008), 023112; Trevisan, D'Isidoro, and Talagrand, Q. J. R. Meteorol. Soc., 136 (2010), pp. 487-496; Palatella, Carrassi, and Trevisan, J. Phys. A, 46 (2013), 254020] and pseudo-orbit data assimilation [Judd and Smith, Phys. D, 151 (2001), pp. 125-141; Judd et al., J. Atmos. Sci., 65 (2008), pp. 1749-1772; Du and Smith, J. Atmos. Sci., 71 (2014), pp. 469-482]. The algorithm utilizes time dependent projections onto the nonstable subspace determined by employing computational techniques for Lyapunov exponents/vectors. The method is extended to parameter estimation without changing the problem dynamics and we address techniques for adapting the method when (as is commonly the case) observations are not available in the full model state space. We use a combination of analysis and numerical experiments (with the Lorenz 63 and Lorenz 96 models) to illustrate the efficacy of the techniques and show that the results compare favorably with other variational techniques.