| Popis: |
State estimation in multi-layer turbulent flow fields with only a single layer of partial observation remains a challenging yet practically important task. Applications include inferring the state of the deep ocean by exploiting surface observations. Directly implementing an ensemble Kalman filter based on the full forecast model is usually expensive. One widely used method in practice projects the information of the observed layer to other layers via linear regression. However, when nonlinearity in the highly turbulent flow field becomes dominant, the regression solution will suffer from large uncertainty errors. In this paper, we develop a multi-step nonlinear data assimilation method. A sequence of nonlinear assimilation steps is applied from layer to layer recurrently. Fundamentally different from the traditional linear regression approaches, a conditional Gaussian nonlinear system is adopted as the approximate forecast model to characterize the nonlinear dependence between adjacent layers. The estimated posterior is a Gaussian mixture, which can be highly non-Gaussian. Therefore, the multi-step nonlinear data assimilation method can capture strongly turbulent features, especially intermittency and extreme events, and better quantify the inherent uncertainty. Another notable advantage of the multi-step data assimilation method is that the posterior distribution can be solved using closed-form formulae under the conditional Gaussian framework. Applications to the two-layer quasi-geostrophic system with Lagrangian data assimilation show that the multi-step method outperforms the one-step method with linear stochastic flow models, especially as the tracer number and ensemble size increase. |
