Uncertainty assessment and applicability of an inversion method for volcanic ash forecasting

Abstract

Abstract. Significant improvements in the way we can observe and model volcanic ash clouds have been obtained since the 2010 Eyjafjallajökull eruption. One major development has been the application of data assimilation techniques, which combine models and satellite observations such that an optimal understanding of ash clouds can be gained. Still, questions remain regarding the degree to which the forecasting capabilities are improved by inclusion of such techniques and how these improvements depend on the data input. This study explores how different satellite data and different uncertainty assumptions of the satellite and a priori emissions affect the calculated volcanic ash emission estimate, which is computed by an inversion method that couples the satellite retrievals and a priori emissions with dispersion model data. Two major ash episodes over 4 days in April and May of the 2010 Eyjafjallajökull eruption are studied. Specifically, inversion calculations are done for four different satellite data sets with different size distribution assumptions in the retrieval. A reference satellite data set is chosen, and the range between the minimum and maximum 4-day average load of hourly retrieved ash is 121 % in April and 148 % in May, compared to the reference. The corresponding a posteriori maximum and minimum emission sum found for these four satellite retrievals is 26 and 47 % of the a posteriori reference estimate for the same two periods, respectively. Varying the assumptions made in the satellite retrieval is seen to affect the a posteriori emissions and modelled ash column loads, and modelled column loads therefore have uncertainties connected to them depending on the uncertainty in the satellite retrieval. By further exploring our uncertainty estimates connected to a priori emissions and the mass load uncertainties in the satellite data, the uncertainty in the a priori estimate is found in this case to have an order-of-magnitudegreater impact on the a posteriori solution than the mass load uncertainties in the satellite. Part of this is explained by a too-high a priori estimate used in this study that is reduced by around half in the a posteriori reference estimate. Setting large uncertainties connected to both a priori and satellite mass load shows that they compensate each other, but the a priori uncertainty is found to be most sensitive. Because of this, an inversion-based emission estimate in a forecasting setting needs well-tested and well-considered assumptions on uncertainties for the a priori emission and satellite data. The quality of using the inversion in a forecasting environment is tested by adding gradually, with time, more observations to improve the estimated height versus time evolution of Eyjafjallajökull ash emissions. We show that the initially too-high a priori emissions are reduced effectively when using just 12 h of satellite observations. More satellite observations (> 12 h), in the Eyjafjallajökull case, place the volcanic injection at higher altitudes. Adding additional satellite...

Investigating Volcanic Plumes from Mt. Etna Eruptions of December 2015 by Means of AVHRR and SEVIRI Data

“…Equation (3), which was derived by using a dataset of historical eruptions, such as other similar empirical formulations (e.g., [38]), provides rough estimates of the volumetric flow rate especially in the presence of ash plumes from small eruptions, which are more significantly affected by windy conditions (e.g., [39,40]). The volumetric flow rate, derived using Equation (3), was then converted into MER (mass eruption rate, kg/s) by multiplying its value by the assumed density of ash (2500 kg/m 3 ; e.g., [41]).…”

Modelling the volcanic ash plume from Eyjafjallajökull eruption (May 2010) over Europe: evaluation of the benefit of source term improvements and of the assimilation of aerosol measurements

“…"Resolved" source terms are the result of an explicit simulation of the thermodynamic and buoyancy processes in the plume (such as the steady 1D model Plumeria; Mastin, 2007) and even the microphysical aerosol processes, including aggregation (FPLUME; Folch et al, 2016). Source inversion of volcanic ash columns measured by satellites has also been developed in different institutes (Stohl et al, 2011;Steensen et al, 2017a;Beckett et al, 2020). The purpose of this work is generally to optimise a source term for which the model ash load columns match observed columns.…”

Volcanic ash forecast using ensemble-based data assimilation: an ensemble transform Kalman filter coupled with the FALL3D-7.2 model (ETKF–FALL3D version 1.0)

“…This allowed for the reconstruction of the full emission profile using data from different sensors. Stohl et al (2011), Kristiansen et al (2012), Denlinger et al (2012), Pelley et al (2015), and Steensen et al (2017) discussed further developments and evaluations of the proposed approach. In particular, Pelley et al (2015) describe the operational implementation of this algorithm at the London Volcanic Ash Advisory Centre (VAAC).…”