A new scheme for sulphur dioxide retrieval from IASI measurements: application to the Eyjafjallajökull eruption of April and May 2010

A new optimal estimation algorithm for the retrieval of sulphur dioxide (SO<sub>2</sub>) has been developed for the Infrared Atmospheric Sounding Interferometer (IASI) using the channels between 1000–1200 and 1300–1410 cm<sup>−1</sup>. These regions include the two SO<sub>2</sub> absorption bands centred at about 8.7 and 7.3 μm (the ν<sub>1</sub> and ν<sub>3</sub> bands respectively). The retrieval assumes a Gaussian SO<sub>2</sub> profile and returns the SO<sub>2</sub> column amount in Dobson units and the altitude of the plume in millibars (mb). Forward modelled spectra (against which the measurements are compared) are based on the Radiative Transfer for TOVS (RTTOV) code. In our implementation RTTOV uses atmospheric profiles from European Centre for Medium-Range Weather Forecasts (ECMWF) meteorological data. The retrieval includes a comprehensive error budget for every pixel derived from an error covariance matrix that is based on the SO<sub>2</sub>-free climatology of the differences between the IASI and forward modelled spectra. The IASI forward model includes the ability to simulate a cloud or ash layer in the atmosphere. This feature is used to illustrate that: (1) the SO<sub>2</sub> retrieval is not affected by underlying cloud but is affected if the SO<sub>2</sub> is within or below a cloud layer; (2) it is possible to discern if ash (or other atmospheric constituents not considered in the error covariance matrix) affects the retrieval using quality control based on the fit of the measured spectrum by the forward modelled spectrum. In this work, the algorithm is applied to follow the behaviour of SO<sub>2</sub> plumes from the Eyjafjallajökull eruption during April and May 2010. From 14 April to 4 May (during Phase I and II of the eruption) the total amount of SO<sub>2</sub> present in the atmosphere, estimated by IASI measurements, is generally below 0.02 Tg. During the last part of the eruption (Phase III) the values are an order of magnitude higher, with a maximum of 0.18 Tg measured on the afternoon of 7 May