Paper detail

A Multiresolution Ensemble Kalman Filter using Wavelet Decomposition

We present a method of using classical wavelet based multiresolution analysis to separate scales in model and observations during data assimilation with the ensemble Kalman filter. In many applications, the underlying physics of a phenomena involve the interaction of features at multiple scales. Blending of observational and model error across scales can result in large forecast inaccuracies since large errors at one scale are interpreted as inexact data at all scales. Our method uses a transformation of the observation operator in order to separate the information from different scales of the observations. This naturally induces a transformation of the observation covariance and we put forward several algorithms to efficiently compute the transformed covariance. Another advantage of our multiresolution ensemble Kalman filter is that scales can be weighted independently to adjust each scale's effect on the forecast. To demonstrate feasibility we present applications to a one dimensional Kuramoto-Sivashinsky (K-S) model with scale dependent observation noise and an application involving the forecasting of solar photospheric flux. The latter example demonstrates the multiresolution ensemble Kalman filter's ability to account for scale dependent model error. Modeling of photospheric magnetic flux transport is accomplished by the Air Force Data Assimilative Photospheric Transport (ADAPT) model.