Accelerating volcanic ash data assimilation using a mask-state algorithm based on an ensemble Kalman filter: a case study with the LOTOS-EUROS model (version 1.10)
Guangliang Fu1,Hai Xiang Lin1,Arnold Heemink1,Sha Lu1,Arjo Segers2,Nils van Velzen1,3,Tongchao Lu4,and Shiming Xu5Guangliang Fu et al. Guangliang Fu1,Hai Xiang Lin1,Arnold Heemink1,Sha Lu1,Arjo Segers2,Nils van Velzen1,3,Tongchao Lu4,and Shiming Xu5
1Delft University of Technology, Delft Institute of Applied Mathematics,
Mekelweg 4, 2628 CD Delft, the Netherlands
2TNO, Department of Climate, Air and Sustainability, P.O. Box 80015,
3508 TA Utrecht, the Netherlands
3VORtech, P.O. Box 260, 2600 AG Delft, the Netherlands.
4School of Mathematics, Shandong University, Jinan, Shandong, China
5Department of Earth System Science, Tsinghua University, Beijing,
China
1Delft University of Technology, Delft Institute of Applied Mathematics,
Mekelweg 4, 2628 CD Delft, the Netherlands
2TNO, Department of Climate, Air and Sustainability, P.O. Box 80015,
3508 TA Utrecht, the Netherlands
3VORtech, P.O. Box 260, 2600 AG Delft, the Netherlands.
4School of Mathematics, Shandong University, Jinan, Shandong, China
5Department of Earth System Science, Tsinghua University, Beijing,
China
Received: 01 Aug 2016 – Discussion started: 24 Aug 2016 – Revised: 07 Feb 2017 – Accepted: 03 Apr 2017 – Published: 24 Apr 2017
Abstract. In this study, we investigate a strategy to accelerate the data assimilation (DA) algorithm. Based on evaluations of the computational time, the analysis step of the assimilation turns out to be the most expensive part. After a study of the characteristics of the ensemble ash state, we propose a mask-state algorithm which records the sparsity information of the full ensemble state matrix and transforms the full matrix into a relatively small one. This will reduce the computational cost in the analysis step. Experimental results show the mask-state algorithm significantly speeds up the analysis step. Subsequently, the total amount of computing time for volcanic ash DA is reduced to an acceptable level. The mask-state algorithm is generic and thus can be embedded in any ensemble-based DA framework. Moreover, ensemble-based DA with the mask-state algorithm is promising and flexible, because it implements exactly the standard DA without any approximation and it realizes the satisfying performance without any change in the full model.
We propose a mask-state algorithm (MS) which records the sparsity information of the full ensemble state matrix and transforms the full matrix into a relatively small one. It will reduce the computational cost in the analysis step for plume assimilation applications. Ensemble-based DA with the mask-state algorithm is generic and flexible, because it implements exactly the standard DA without any approximation and it realizes the satisfying performance without any change of the full model.
We propose a mask-state algorithm (MS) which records the sparsity information of the full...
