Abstract. The accurate determination of the location, height, and loading of sulfur dioxide (SO2) plumes emitted by volcanic eruptions is essential for aviation safety. The SO2 layer height is also one of the most critical parameters with respect to determining the impact on the climate. Retrievals of SO2 plume height have been carried out using satellite UV backscatter measurements, but, until now, such algorithms are very time-consuming. We have developed an extremely fast yet accurate SO2 layer height retrieval using the Full-Physics Inverse Learning Machine (FP_ILM) algorithm. This is the first time the algorithm has been applied to measurements from the TROPOMI instrument onboard the Sentinel-5 Precursor platform. In this paper, we demonstrate the ability of the FP_ILM algorithm to retrieve SO2 plume layer heights in near-real-time applications with an accuracy of better than 2 km for SO2 total columns larger than 20 DU. We present SO2 layer height results for the volcanic eruptions of Sinabung in February 2018, Sierra Negra in June 2018, and Raikoke in June 2019, observed by TROPOMI.

No abstract

In this paper we analyze the accuracy and efficiency of several radiative transfer models for inferring cloud parameters from radiances measured by the Earth Polychromatic Imaging Camera (EPIC) on board the Deep Space Climate Observatory (DSCOVR). The radiative transfer models are the exact discrete ordinate and matrix operator methods with matrix exponential, and the approximate asymptotic and equivalent Lambertian cloud models. To deal with the computationally expensive radiative transfer calculations, several acceleration techniques such as, for example, the telescoping technique, the method of false discrete ordinate, the correlated k-distribution method and the principal component analysis (PCA) are used. We found that, for the EPIC oxygen A-band absorption channel at 764 nm, the exact models using the correlated k-distribution in conjunction with PCA yield an accuracy better than 1.5 % and a computation time of 18 s for radiance calculations at 5 viewing zenith angles.

Precise knowledge of the location and height of the volcanic sulphur dioxide (SO 2) plume is essential for accurate determination of SO 2 emitted by volcanic eruptions. Current SO 2 plume height retrieval algorithms based on ultraviolet (UV) satellite measurements are very time-consuming and therefore not suitable for near-real-time applications. In this work we present a novel method called the full-physics inverse learning machine (FP-ILM) algorithm for extremely fast and accurate retrieval of the SO 2 plume height. FP-ILM creates a mapping between the spectral radiance and the geophysical parameters of interest using supervised learning methods. The FP-ILM combines smart sampling methods, dimensionality reduction techniques, and various linear and non-linear regression analysis schemes based on principal component analysis and neural networks. The computationally expensive operations in FP-ILM are the radiative transfer model computations of a training dataset and the determination of the inversion operatorthese operations are performed off-line. The application of the resulting inversion operator to real measurements is extremely fast since it is based on calculations of simple regression functions. Retrieval of the SO 2 plume height is demonstrated for the volcanic eruptions of Mt. Kasatochi (in 2008) and Eyjafjallajökull (in 2010), measured by the GOME-2 (Global Ozone Monitoring Instrument-2) UV instrument on-board MetOp-A.

No abstract

No abstract

No abstract

The new generation of atmospheric composition sensors such as TROPOMI is capable of providing spectra of high spatial and spectral resolution. To process this vast amount of spectral information, fast radiative transfer models (RTMs) are required. In this regard, we analyzed the efficiency of two acceleration techniques based on the principal component analysis (PCA) for simulating the Hartley-Huggins band spectra. In the first one, the PCA is used to map the data set of optical properties of the atmosphere to a lower-dimensional subspace, in which the correction function for an approximate but fast RTM is derived. The second technique is based on the dimensionality reduction of the data set of spectral radiances. Once the empirical orthogonal functions are found, the whole spectrum can be reconstructed by performing radiative transfer computations only for a specific subset of spectral points. We considered a clear-sky atmosphere where the optical properties are defined by Rayleigh scattering and trace gas absorption. Clouds can be integrated into the model as Lambertian reflectors. High computational performance is achieved by combining both techniques without losing accuracy. We found that for the Hartley-Huggins band, the combined use of these techniques yields an accuracy better than 0.05% while the speedup factor is about 20. This innovative combination of both PCA-based techniques can be applied in future works as an efficient approach for simulating the spectral radiances in other spectral regions.