A Near-Real-Time Method for Estimating Volcanic Ash Emissions Using Satellite Retrievals

We present a Bayesian inversion method for estimating volcanic ash emissions using satellite retrievals of ash column load and an atmospheric dispersion model. An a priori description of the emissions is used based on observations of the rise height of the volcanic plume and a stochastic model of th...

Full description

Bibliographic Details
Published in:Atmosphere
Main Authors: Rachel E. Pelley, David J. Thomson, Helen N. Webster, Michael C. Cooke, Alistair J. Manning, Claire S. Witham, Matthew C. Hort
Format: Article in Journal/Newspaper
Language:English
Published: MDPI AG 2021
Subjects:
Online Access:https://doi.org/10.3390/atmos12121573
https://doaj.org/article/3dee3478807341638f576b43e34ddab2
Description
Summary:We present a Bayesian inversion method for estimating volcanic ash emissions using satellite retrievals of ash column load and an atmospheric dispersion model. An a priori description of the emissions is used based on observations of the rise height of the volcanic plume and a stochastic model of the possible emissions. Satellite data are processed to give column loads where ash is detected and to give information on where we have high confidence that there is negligible ash. An atmospheric dispersion model is used to relate emissions and column loads. Gaussian distributions are assumed for the a priori emissions and for the errors in the satellite retrievals. The optimal emissions estimate is obtained by finding the peak of the a posteriori probability density under the constraint that the emissions are non-negative. We apply this inversion method within a framework designed for use during an eruption with the emission estimates (for any given emission time) being revised over time as more information becomes available. We demonstrate the approach for the 2010 Eyjafjallajökull and 2011 Grímsvötn eruptions. We apply the approach in two ways, using only the ash retrievals and using both the ash and clear sky retrievals. For Eyjafjallajökull we have compared with an independent dataset not used in the inversion and have found that the inversion-derived emissions lead to improved predictions.