Four-dimensional ensemble variational data assimilation and the unstable subspace
The performance of (ensemble) Kalman filters used for data assimilation in the geosciences critically depends on the dynamical properties of the evolution model. A key aspect is that the error covariance matrix is asymptotically supported by the unstable–neutral subspace only, i.e. it is spanned by...
| Main Authors: | , |
|---|---|
| Format: | Article |
| Language: | English |
| Published: |
Stockholm University Press
2017-01-01
|
| Series: | Tellus: Series A, Dynamic Meteorology and Oceanography |
| Subjects: | |
| Online Access: | http://dx.doi.org/10.1080/16000870.2017.1304504 |
| _version_ | 1818216445652238336 |
|---|---|
| author | Marc Bocquet Alberto Carrassi |
| author_facet | Marc Bocquet Alberto Carrassi |
| author_sort | Marc Bocquet |
| collection | DOAJ |
| description | The performance of (ensemble) Kalman filters used for data assimilation in the geosciences critically depends on the dynamical properties of the evolution model. A key aspect is that the error covariance matrix is asymptotically supported by the unstable–neutral subspace only, i.e. it is spanned by the backward Lyapunov vectors with non-negative exponents. The analytic proof of such a property for the Kalman filter error covariance has been recently given, and in particular that of its confinement to the unstable–neutral subspace. In this paper, we first generalize those results to the case of the Kalman smoother in a linear, Gaussian and perfect model scenario. We also provide square-root formulae for the filter and smoother that make the connection with ensemble formulations of the Kalman filter and smoother, where the span of the error covariance is described in terms of the ensemble deviations from the mean. We then discuss how this neat picture is modified when the dynamics are nonlinear and chaotic, and for which analytic results are precluded or difficult to obtain. A numerical investigation is carried out to study the approximate confinement of the anomalies for both a deterministic ensemble Kalman filter (EnKF) and a four-dimensional ensemble variational method, the iterative ensemble Kalman smoother (IEnKS), in a perfect model scenario. The confinement is characterized using geometrical angles that determine the relative position of the anomalies with respect to the unstable–neutral subspace. The alignment of the anomalies and of the unstable–neutral subspace is more pronounced when observation precision or frequency, as well as the data assimilation window length for the IEnKS, are increased. These results also suggest that the IEnKS and the deterministic EnKF realize in practice (albeit implicitly) the paradigm behind the approach of Anna Trevisan and co-authors known as the assimilation in the unstable subspace. |
| first_indexed | 2024-12-12T06:52:05Z |
| format | Article |
| id | doaj.art-446df398baa044b9aaf4107d5028dc08 |
| institution | Directory Open Access Journal |
| issn | 1600-0870 |
| language | English |
| last_indexed | 2024-12-12T06:52:05Z |
| publishDate | 2017-01-01 |
| publisher | Stockholm University Press |
| record_format | Article |
| series | Tellus: Series A, Dynamic Meteorology and Oceanography |
| spelling | doaj.art-446df398baa044b9aaf4107d5028dc082022-12-22T00:34:03ZengStockholm University PressTellus: Series A, Dynamic Meteorology and Oceanography1600-08702017-01-0169110.1080/16000870.2017.13045041304504Four-dimensional ensemble variational data assimilation and the unstable subspaceMarc Bocquet0Alberto Carrassi1CEREA, joint laboratory École des Ponts ParisTech and EDF R&D, Université Paris-EstNansen Environmental and Remote Sensing CenterThe performance of (ensemble) Kalman filters used for data assimilation in the geosciences critically depends on the dynamical properties of the evolution model. A key aspect is that the error covariance matrix is asymptotically supported by the unstable–neutral subspace only, i.e. it is spanned by the backward Lyapunov vectors with non-negative exponents. The analytic proof of such a property for the Kalman filter error covariance has been recently given, and in particular that of its confinement to the unstable–neutral subspace. In this paper, we first generalize those results to the case of the Kalman smoother in a linear, Gaussian and perfect model scenario. We also provide square-root formulae for the filter and smoother that make the connection with ensemble formulations of the Kalman filter and smoother, where the span of the error covariance is described in terms of the ensemble deviations from the mean. We then discuss how this neat picture is modified when the dynamics are nonlinear and chaotic, and for which analytic results are precluded or difficult to obtain. A numerical investigation is carried out to study the approximate confinement of the anomalies for both a deterministic ensemble Kalman filter (EnKF) and a four-dimensional ensemble variational method, the iterative ensemble Kalman smoother (IEnKS), in a perfect model scenario. The confinement is characterized using geometrical angles that determine the relative position of the anomalies with respect to the unstable–neutral subspace. The alignment of the anomalies and of the unstable–neutral subspace is more pronounced when observation precision or frequency, as well as the data assimilation window length for the IEnKS, are increased. These results also suggest that the IEnKS and the deterministic EnKF realize in practice (albeit implicitly) the paradigm behind the approach of Anna Trevisan and co-authors known as the assimilation in the unstable subspace.http://dx.doi.org/10.1080/16000870.2017.1304504assimilation in the unstable subspaceensemble Kalman filteriterative ensemble Kalman smoother four-dimensional ensemble variational methodLyapunov vectors |
| spellingShingle | Marc Bocquet Alberto Carrassi Four-dimensional ensemble variational data assimilation and the unstable subspace Tellus: Series A, Dynamic Meteorology and Oceanography assimilation in the unstable subspace ensemble Kalman filter iterative ensemble Kalman smoother four-dimensional ensemble variational method Lyapunov vectors |
| title | Four-dimensional ensemble variational data assimilation and the unstable subspace |
| title_full | Four-dimensional ensemble variational data assimilation and the unstable subspace |
| title_fullStr | Four-dimensional ensemble variational data assimilation and the unstable subspace |
| title_full_unstemmed | Four-dimensional ensemble variational data assimilation and the unstable subspace |
| title_short | Four-dimensional ensemble variational data assimilation and the unstable subspace |
| title_sort | four dimensional ensemble variational data assimilation and the unstable subspace |
| topic | assimilation in the unstable subspace ensemble Kalman filter iterative ensemble Kalman smoother four-dimensional ensemble variational method Lyapunov vectors |
| url | http://dx.doi.org/10.1080/16000870.2017.1304504 |
| work_keys_str_mv | AT marcbocquet fourdimensionalensemblevariationaldataassimilationandtheunstablesubspace AT albertocarrassi fourdimensionalensemblevariationaldataassimilationandtheunstablesubspace |