Parallel Filter Algorithms for Data Assimilation in Oceanography
A consistent systematic comparison of filter algorithms based on the Kalman filter and intended for data assimilation with large-scale nonlinear models is presented. Considered are the EnsembleKalman Filter (EnKF), the Singular Evolutive Extended Kalman (SEEK) filter, and the Singular Evolutive Interpolated Kalman (SEIK) filter. Within the two parts of this thesis, the filter algorithms are compared with a focus on their mathematical properties as Error Subspace KalmanFilters (ESKF). Further, the filters are studied as parallel algorithms. This study includes the development of an efficient framework for parallel filtering. In the first part, the filters are motivated in the context of statistical estimation. The unified interpretation as ESKF algorithms provides the basis for the consistent comparison of the filters. Numerical data assimilation experiments with a model based on the shallow water equations show how choices of the filter schemeand particular state ensembles for the filter initialization lead to variations of the data assimilation performance.The application of the three filter algorithms on parallel computers is studied in the second part. The parallelization possibilities of the different phases of the algorithms are examined. Further, a framework for parallel filtering is developed which allows to combine filter algorithms with existing numerical models requiring only minimal changes to the source code of the model.The framework is used to combine the parallel filters with the 3D finite element ocean model FEOM. Numerical data assimilation experiments are utilized to assess the parallel efficiency of the filtering framework and the parallel filters. The experiments yield an excellent parallel efficiency for the filtering framework. Further, the framework and the filter algorithms are well suited for application to realistic large-scale data assimilation problems.
Helmholtz Research Programs > MARCOPOLI (2004-2008) > German community ocean model