Particle-filter-based volcanic ash emission inversion applied to a hypothetical sub-Plinian Eyjafjallajökull eruption using the Ensemble for Stochastic Integration of Atmospheric Simulations (ESIAS-chem) version 1.0

• Published in: Geoscientific Model Development
• Main Authors: Franke, Philipp, Lange, Anne Caroline, Elbern, Hendrik
• Format: Article in Journal/Newspaper
• Language: English
• Published: Copernicus Publications 2022
• Subjects: article; Verlagsveröffentlichung; Eyjafjallajökull; Iceland
• Online Access: https://doi.org/10.5194/gmd-15-1037-2022; https://noa.gwlb.de/receive/cop_mods_00060121; https://noa.gwlb.de/servlets/MCRFileNodeServlet/cop_derivate_00059770/gmd-15-1037-2022.pdf; https://gmd.copernicus.org/articles/15/1037/2022/gmd-15-1037-2022.pdf

A particle-filter-based inversion system is presented, which enables us to derive time- and altitude-resolved volcanic ash emission fluxes along with its uncertainty. The system assimilates observations of volcanic ash column mass loading as retrieved from geostationary satellites. It aims to estimate the temporally varying emission profile endowed with error margins. In addition, we analyze the dependency of our estimate on wind field characteristics, notably vertical shear, within variable observation intervals. Thus, the proposed system addresses the special challenge of analyzing the vertical profile of volcanic ash clouds given only 2D high temporal-resolution column mass loading data as retrieved by geostationary satellites. The underlying method rests on a linear combination of height–time emission finite elements of arbitrary resolution, each of which is assigned to a model run subject to ensemble-based space–time source inversion. Employing a modular concept, this setup builds the Ensemble for Stochastic Integration of Atmospheric Simulations (ESIAS-chem). It comprises a particle smoother in combination with a discrete-grid ensemble extension of the Nelder–Mead minimization method. The ensemble version of the EURopean Air pollution Dispersion – Inverse Model (EURAD-IM) is integrated into ESIAS-chem but can be replaced by other models. As initial validation of ESIAS-chem, the system is applied to simulated artificial observations of both ash-contaminated and ash-free atmospheric columns using identical-twin experiments. Thus, in this idealized initial performance test the underlying meteorological uncertainty is neglected. The inversion system is applied to two notional sub-Plinian eruptions of the Eyjafjallajökull volcano, Iceland, with strong ash emission changes with time and injection heights. It demonstrates the ability of ESIAS-chem to retrieve the volcanic ash emission fluxes from the assimilation of column mass loading data only. However, the analyzed emission profiles strongly differ in their ...