A Multiresolution Ensemble Kalman Filter using Wavelet Decomposition

TL;DR

This paper introduces a multiresolution ensemble Kalman filter that uses wavelet analysis to separate scales in data assimilation, improving forecast accuracy by handling scale-dependent errors and weights.

Contribution

It proposes a novel wavelet-based multiresolution approach for ensemble Kalman filtering, allowing independent scale weighting and efficient covariance computation.

Findings

Successfully applied to a 1D Kuramoto-Sivashinsky model with scale-dependent noise.

Demonstrated effectiveness in forecasting solar photospheric flux.

Enhanced handling of scale-dependent model errors.

Abstract

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.