Infrared Hyperspectral and Ultraviolet Remote Measurements of Volcanic Gas Plume at MT Etna during IMAGETNA Campaign

: Quantiﬁcation of gaseous emission ﬂuxes from volcanoes can yield valuable insights on processes occurring in the Earth’s interior as part of hazard monitoring. It is also an important task in the framework of climate change, in order to reﬁne estimates of natural emissions. Passive open-path UltraViolet (UV) scattered observation by UV camera allows the imaging of volcanic plumes and evaluation of sulfur dioxide (SO 2 ) ﬂuxes at high temporal resolution during daytime. Another technique of imaging is now available in the InfraRed (IR) spectral domain. Infrared hyperspectral imagers have the potential to overcome the boundary of daytime sampling of the UV, providing measurements also during the night and giving access simultaneously to additional relevant gas species. In this context the IMAGETNA campaign of measurements took place at Mt Etna (Italy) in June 2015. Three di ﬀ erent IR imagers (commercial and under developments) were deployed, together with a Fourier Transform InfraRed spectrometer (FTIR) instrument, a UV camera, a Long Wavelength InfraRed (LWIR) camera and a radiometer. We present preliminary results obtained by the two IR cameras under development, and then the IR hyperspectral imager results, coming from full physics retrieval, are compared to those of the UV camera. The comparison points out an underestimation of the SO 2 Slant Column Densities (SCD) of the UV camera by a factor of 3.6. The detailed study of the retrieved SO 2 SCD highlights the promising application of IR imaging in volcanology for remotely volcanic plume gas measurements. It also provides a way to investigate uncertainties in the SO 2 SCD imaging in the UV and the IR. comparisons, highlights the need to quantify the scattering e ﬀ ect on UV camera remote sensing measurements for the characterization of volcanic plume emissions.


Introduction
Volcanoes are a substantial source of gases and aerosols to the atmosphere. A major scientific challenge is the quantification of volcanic gas emissions for the understanding of eruptive processes and to assess the overall impact and climate change response [1][2][3][4]. It has been demonstrated that emissions of sulfur dioxide (SO 2 ) by volcanoes in quiescent non-eruptive stages (passive degassing) [5][6][7][8][9] may significantly impact the overall sulphur budget of the Earth's atmosphere compared to other natural sources or anthropogenic activities. Though passive SO 2 emission does not extend to the stratosphere domain respect to episodic SO 2 emission during eruption [10,11], this species is persistently released by volcanoes in the troposphere at different temporal and spatial variability representing a climate-forcing factor [12][13][14]. In addition, on open-system volcanoes, volcanic degassing exerts a first-order role in the style and timing of volcanic eruption. Changes in the flux rate or the chemical composition of gases emitted from volcanoes can accompany or precede changes in volcanic activity i.e., represent a signal of coming eruption [15,16]. This motivates the monitoring of volcanic volatile emissions at numerous volcanoes around the world [17][18][19]. Sulphur dioxide (SO 2 ) is typically one of the most abundant volatile components in volcanic plumes (after H 2 O and CO 2 ) and its strong UltraViolet (UV) absorption band (and low abundance in the background atmosphere) allows relatively straightforward detection [20,21]. In addition, SO 2 flux is a particularly useful parameter to monitor in the volcanic environments, as it can be considered a proxy for magma flux, assuming complete exsolution of gas and lossless transport to the point of measurement [22]. Volcanic SO 2 emission rate measurements have been carried out using ground-based UV remote sensing techniques for more than 30 years, becoming a standard practice in volcanic surveillance [17,18,23]. Emission fluxes have been documented by COSPEC [24,25] and networks of ultraviolet scanner spectrometers [19,26] based on differential optical absorption spectroscopy (DOAS). Another technique is UV camera technology [27][28][29] that uses band-pass filters to visualize the SO 2 plume as a 2D image. UV cameras can monitor SO 2 emission flux variability at much higher frequency (~2 Hz) than COSPEC/DOAS spectrometers, but similarly are subject to uncertainties due to time-varying wind field and possible ground effects for measurements near the point of gas emission.
Volcanic ash plumes and gas composition is also detected remotely in the infrared region by open-path Fourier Transform InfraRed (FTIR) [30][31][32][33]. Thermal infrared (8 to 14 µm) offers advantages, with respect to UV, such as the possibility to observe SO 2 flux during both day and night in various configurations (absorption and emission mode). Infrared cameras using band-pass filters, like the Cyclops camera developed by Prata and Bernardo [34], can offer the possibility to detect volcanic SO 2 . FTIR enables to retrieve simultaneously several volcanic species not observable by the UV, and has been applied in volcanology in absorption mode to detect volcanic gases (H 2 O, CO 2 , CO 2 , CH 4 , N 2 O, HCl, HF and CO), using hot rocks or lava as source [35,36]. However, the measurement is limited to 1D (single pixel). Hyperspectral InfraRed (IR) imagers offer the potential to image SO 2 and various other gas species in 2D in emission mode during day and night.
Whilst multi-spectral and hyperspectral sensing infrared imagers have been relatively common in space satellite application [37,38] but their ground-based counterparts are not widely used, though have prospective applications in volcanology [39,40]. The suitability of multi-pixel spectroscopy and hyperspectral imaging relies on the opportunity to simultaneously measure spectral radiance disturbed by gas signatures and the background of spectral radiance. The combined advantages of spectroscopy and imaging allow the acquisition of an image in which each pixel is associated with a determined portion of the spectrum, and where each image is a two-dimensional map of concentrations of gas compounds. Hyperspectral imagers adapted and/or optimized for gas detection in the atmosphere have mainly been developed for military applications [41] and secondarily for atmosphere characterization. For example, the Gimballed Limb Observer for Radiance Imaging of the Atmosphere (GLORIA) instrument (IMK-KIT, Germany) is used for airborne or balloon-borne measurement campaigns to characterize the chemical composition (H 2 O, O 3 , HNO 3 , CH 4 , N 2 O, CFC-11, CFC-12, ...) of the upper troposphere and lower stratosphere [42][43][44].
Remote Sens. 2019, 11, 1175 3 of 23 In this study, we used the commercial instrument Hyper-Cam Long Wavelength InfraRed (LWIR) from TELOPS Company and two multispectral IR cameras in the MidWave InfraRed (MWIR) under development at ONERA (Office National d'Etudes et de Recherches Aerospatiales) for remote observation of volcanic gas in the plume of Mt Etna. Measurements were performed side-by-side with UV imaging camera to measure SO 2 Slant Column Densities (SCDs) in the plume and to evaluate the discrepancy between measurements obtained in the IR and the UV.
The objectives of research are threefold: first is to achieve quantitative measurements of volcanic gas; second to compare the results achieved with an IR imaging system with those simultaneous measured by UV camera to explore uncertainty and to assess the suitability of the current performance of hyperspectral imaging for measuring volcanic plume gases; and third to test and quantify the performance of our hyperspectral data processing algorithm. A secondary goal of the campaign was to test new lightweight multi or hyperspectral instruments in operational conditions, in order to estimate the ability to these new instrumental concepts to deliver useful data for future campaigns.
In the following section, we present the campaign context and the characteristics of the different instruments deployed. We report preliminary results from the two IR multispectral cameras under development and make recommendations for future use. After checking the quality of the radiance spectrum obtained by Hyper-Cam LWIR hyperspectral imager by comparing with the widely used single-pixel FTIR instrument, we present the radiative transfer retrieval code and method applied to retrieve SO 2 slant column densities. We also present sensitivity tests on the IR retrieval parameters. Then, in the last part we discuss in detail the SO 2 slant column densities results obtained by the Hyper-Cam LWIR hyperspectral imager and compare them to those obtained simultaneously by the UV camera.

IMAGETNA Campaign
IMAGETNA campaign was carried out between 21 and 26 June 2015. Volcanic plume observations were remotely made from the Pizzi De Neri Volcano observatory located on the North side of Mt Etna at 2,850 m of altitude. Over the five days of field activity, Mt Etna was in a quiescent stage of volcanic activity. Emissions from the four summit craters of the volcano resulted in a near-source plume that was ash poor, but variable in its consistency and magnitude. Due to the wind-speed and direction on the days of field activity and the location of the sampling site, our measurements relate mainly to the volcanic plume issued from North East Craters ( Figure 1). The site of measurements allowed us to image the plume almost perpendicular to the axis of the wind and from a moderately far distance of 2 km from the main emission area such as to observe simultaneously the volcanic plume and the clear sky. This geometry configuration allows us to perform limb remote sensing measurements under a 15 • angle with respect to the horizontal.
Measurements have been performed during the five days of the campaign early in the morning (starting at 6:00 Local Time) to get the maximum of thermal contrast between the plume and the clear sky and also to prevent the presence of convective clouds, which develop in the afternoon.
The first day of the campaign was dedicated to the instrumentation deployment, tests and calibration. The collection of the data was properly carried out in four days. Nevertheless, among of recorded data, some sequences were rejected due to meteorological conditions such as the sporadic presence of clouds, which induce too strong opacity for infrared remote sensing measurements but also strong winds inducing too much vibration of the instruments. Good quality measurements sequences were recorded between 25 and 26 June. speed and direction on the days of field activity and the location of the sampling site, our measurements relate mainly to the volcanic plume issued from North East Craters (Figure 1). The site of measurements allowed us to image the plume almost perpendicular to the axis of the wind and from a moderately far distance of ~ 2 km from the main emission area such as to observe simultaneously the volcanic plume and the clear sky. This geometry configuration allows us to perform limb remote sensing measurements under a 15 ° angle with respect to the horizontal.  Data collection was performed deploying six instruments in order to observe the gas plume using different spectral regions, at different spectral resolution, and using the properties of gas emission under limb viewing (Table 1). Three IR imagers were deployed, two of them (SIBI, which stands for "Spectro-Imageur Biréfringent" and Filament) are under development at ONERA and were deployed out of the laboratory for the first time, the third one was provided by TELOPS Company (Hyper-Cam LWIR). Together with the three IR imaging systems, measurements were also carried out by a FTIR single pixel spectrometer in emission mode (OPAG 33: Open-Path Gas Analyzer, Brucker company), which covers the entire spectral range of the different IR imagers, and a UV camera for the SO 2 measurements. Each instrument differs by the physical technology described below. JAI CM-140GE-UV a 640 × 512 is the raw image size, a spatial scan is necessary to measure the hyperspectral cube; b more and newer filters were used for the IMAGETNA Campaign then in the proceeding from Sakat et al. [46]. Both, UV camera and FTIR, were placed side-by-side to capture the entire scene allowing to complement the IR devices. We thus have the possibility to compare SO 2 Slant Column Densities (SCD) images obtained simultaneously in the IR region from Hyper-Cam LWIR and in the UV region from the UV camera.
The FTIR spectrometer and the Filament imager were set up on the same platform assuring the instruments looking at the same scene. The SIBI, the Hyper-Cam LWIR, and the UV camera instruments were deployed close to each other looking in the same direction to cover at least a part of the same scene ( Figure 2 instruments were deployed close to each other looking in the same direction to cover at least a part of the same scene ( Figure 2).

Lightweight SIBI Imager
SIBI, is a novel compact hyperspectral camera in the mid-IR which uses Static Fourier Transform (SFT) to produce a hyperspectral data cube with a moderate spectral resolution and a full width at half maximum (FWHM) of 14 cm -1 . SIBI is a very compact concept in comparison with existing hyperspectral cameras [45]. This is possible by the use of a birefringent interferometer, and a very compact cryogenic setup. SIBI camera combines a static interferometer and an imaging lens to obtain an interference pattern on each raw image. Thanks to this interference pattern and to a spatial scan of the scene through a rotation of the whole instrument, we can recover the data cube [48,49]. The size reduction of SIBI, compared to other cameras, is firstly due to the use of a birefringent interferometer, based on polarization splitting. Indeed, such an interferometer amounts to a plate merely set in front of the camera (more precisely, this plate is a stack of two parallel plates between two polarizers). The second cause of size reduction comes from the use of a very compact cryogenic camera. Details can be found in PhD manuscript of Pola Fossi [50]. A first prototype of the spectral imager (12 × 12 × 20 cm 3 and 2.7 kg) has been involved in the IMAGETNA campaign for proof of concept. It was installed on a rotation stage ( Figure 2). This operation generated an interferometric cube which was post-processed through four steps to create the spectra of the points of the scene. The first step is to correct the image by subtracting the background, removing bad pixels and performing digital gain correction. Then, the optical image distortion is corrected using a pinhole camera model and homography for the camera panning distortion. A correlation algorithm is applied to register images and thus to obtain a full interferogram for each point of the scene. Spectrum is calculated by a Fourier transform from this interferogram and is radiometrically calibrated to deliver the spectral radiance at the entrance window level.

Filament imager
Filament is a multispectral mid-IR imager. Its concept is derived from the Tombo system [51]: the pupil is divided by an array of micro-lenses, to form multiple images (which we will call "small images") of the scene on the focal plane array. Besides, a spectral filter is set on each micro-lens, so that each small image of the scene is associated to one spectral content. We thus obtain a snapshot multispectral imager. It is a compact instrument fully cryogenic to reduce the instrumental

Lightweight SIBI Imager
SIBI, is a novel compact hyperspectral camera in the mid-IR which uses Static Fourier Transform (SFT) to produce a hyperspectral data cube with a moderate spectral resolution and a full width at half maximum (FWHM) of 14 cm −1 . SIBI is a very compact concept in comparison with existing hyperspectral cameras [45]. This is possible by the use of a birefringent interferometer, and a very compact cryogenic setup. SIBI camera combines a static interferometer and an imaging lens to obtain an interference pattern on each raw image. Thanks to this interference pattern and to a spatial scan of the scene through a rotation of the whole instrument, we can recover the data cube [48,49]. The size reduction of SIBI, compared to other cameras, is firstly due to the use of a birefringent interferometer, based on polarization splitting. Indeed, such an interferometer amounts to a plate merely set in front of the camera (more precisely, this plate is a stack of two parallel plates between two polarizers). The second cause of size reduction comes from the use of a very compact cryogenic camera. Details can be found in PhD manuscript of Pola Fossi [50]. A first prototype of the spectral imager (12 × 12 × 20 cm 3 and 2.7 kg) has been involved in the IMAGETNA campaign for proof of concept. It was installed on a rotation stage ( Figure 2). This operation generated an interferometric cube which was post-processed through four steps to create the spectra of the points of the scene. The first step is to correct the image by subtracting the background, removing bad pixels and performing digital gain correction. Then, the optical image distortion is corrected using a pinhole camera model and homography for the camera panning distortion. A correlation algorithm is applied to register images and thus to obtain a full interferogram for each point of the scene. Spectrum is calculated by a Fourier transform from this interferogram and is radiometrically calibrated to deliver the spectral radiance at the entrance window level.

Filament Imager
Filament is a multispectral mid-IR imager. Its concept is derived from the Tombo system [51]: the pupil is divided by an array of micro-lenses, to form multiple images (which we will call "small images") of the scene on the focal plane array. Besides, a spectral filter is set on each micro-lens, so that each Remote Sens. 2019, 11, 1175 6 of 23 small image of the scene is associated to one spectral content. We thus obtain a snapshot multispectral imager. It is a compact instrument fully cryogenic to reduce the instrumental background, currently cooled by liquid Nitrogen. The innovation in Filament lies in the spectral filters. Indeed, these filters are subwavelength gratings [52] made of narrow slits in a gold layer, deposited on a dielectric layer. Since the grating period is lower than the wavelength, the diffracted orders are evanescent, and the zero th order transmission is extremely low because of the reflective gold layer. However, at a specific wavelength, the first diffracted order in transmission is coupled into the waveguide formed by the dielectric layer and then re-diffracted in the zero order by the grating. Thus, the component exhibits a quite sharp transmission at this specific wavelength. Thanks to a 2D geometry for the grating, the transmitted light is unpolarized.

Instrument Description
The Hyper-Cam LWIR from TELOPS Company combines an IR camera as well as a FTIR which means that for every pixel in an image we have an infrared spectrum. It provides an excellent image quality with a 320 × 256 pixels Focal-Plane Array (FPA) detector. An external telescope for magnification of 0.25 was added to the camera during the campaign in order to widen the field of view (FOV) by 4 (25.6 • × 20.4 • ). The Hyper-Cam LWIR offers user-selectable spectral resolution from 0.25 cm −1 to 150 cm −1 . During the campaign, we have mainly used spectral resolution of 2 and 4 cm −1 . A radiometric calibration was performed before and at regular intervals between measurement sequence using two black bodies at different temperatures. The black bodies temperatures must be close to the scene temperature for higher accuracy (~1 • C). Over the campaign, data were recorded selecting the temperatures of 15 and 25 • C, since the use of lower temperatures would have generated condensation on the blackbodies.
The Hyper-Cam acquisition software "Reveal Pro 5", developed by TELOPS Company, was used to record the data. Each file in a given sequence of measurement is called datacube. Datacubes, which corresponds to an IR image with an IR spectrum for each pixel, as said above, are available in raw format (interferogram) or in calibrated radiance. The acquisition time for each datacube is fully dependent on both spectral and spatial resolutions. In order to capture the plume dynamics, we've adjusted the spectral and spatial resolutions of the Hyper-Cam LWIR to get an acquisition time between 1 to 2.5 seconds per datacube.

Retrieval Methodology
The LATMOS Atmospheric Retrieval Algorithm (LARA) used to retrieve SO 2 slant column densities from Hyper-Cam LWIR instrument measurements has been developed over the past decade [53,54] (LATMOS: Laboratoire Atmosphère, Milieux, Observations Spatiales). It has been widely used to analyze atmospheric spectra recorded from ground-based, balloon or satellite-borne experiments, both in absorption or emission mode, but also for the limb or the nadir geometry: LPMA-balloon [55,56], MIPAS-Envisat satellite [57], IASI-Metop satellite [58,59]. LARA allows the simultaneous retrieval of spectra in several spectral windows for the joint retrieval of slant column densities of various species. Surface temperature and emissivity, and if needed instrumental line shape or instrumental spectral shift, may be fitted together with the species. The retrieval algorithm chooses the optimal estimation method and includes an accurate line-by-line radiative transfer model and an efficient minimization algorithm of the Levenberg-Marquardt type [60]. The full error covariance matrix is calculated within the retrieval process (measurement noise and the smoothing error as described in Rodgers, 2000 [60]). The forward model (i.e., the radiative transfer model) uses molecular parameters which are mainly extracted from the HITRAN 2012 database [61]. Individual line shapes are calculated with a Voigt profile based on the Lorentzian parameters listed in the spectroscopic database and the line shifting coefficient can be used when non-zero in HITRAN 2012. The calculation takes into account the water vapor continuum [62] as well as water vapor self-broadening. The reflected downward flux and the reflected or diffused sunlight are modeled. For the present work, the algorithm was tailored to the specificities of Hyper-Cam LWIR spectra, Mt Etna plume characteristics and shape, and the limb geometry. As in Rix et al. [63], the plume shape is assumed with a predefined central plume height Z plume and a Gaussian distribution with a half-width at half-maximum equal to plume thickness Tk plume . Assumptions on the thickness of the plume and the altitude of the plume centre were made by combining visual observations and the characteristics of the relief of the Mt Etna.
The characteristics of the plume are defined according to the following 3 parameters: 1.
The additional temperature inside the plume compared to ambient air ∆ T plume = 1 K 2.
The thickness of the plume Tk plume = 400 m 3.
The altitude of the plume centre (where gas and particulate due to Etna emission is supposed to be maximum) Z plume = 3.2 km To take into account particles contribution to the IR spectrum, we parametrize that contribution considering a spectral dependency of the optical thickness (first order linear modeling) of the plume. The line of sight (limb view) considered is divided into 43 layers from the instrument altitude (2.8 km) up to the top of the atmosphere (80 km), with layers of 100 m in the plume. For each path in each layer, temperature, pressure, path length, air concentration and slant column densities (SCD) are estimated. The temperature and pressure have been derived from the combination of meteorological sounding data over Trapani performed the same day and ECMWF (European Centre of Medium-range Weather Forecast) analyses. Trapani is the closest radio-sounding station to Mt Etna. SO 2 is retrieved in the spectral window 1,100-1,200 cm −1 , considering H 2 O, CO 2 , O 3 , N 2 O, CO, CH 4 and SO 2 SCDs, as well as optical thickness due to particles in the plume contributing to the IR spectrum. Particles optical thickness and SO 2 SCD are retrieved in the plume whereas all the other species are distributed on the entire slant column. It is a full physic retrieval.
For each spectrum, the following parameters have been retrieved (state vector): • One scaling factor for the initial profiles of each molecular species (H 2 O, CO 2 , O 3 , N 2 O, CO, and CH 4 ) except SO 2 ; • One scaling factor of the SO 2 Volume Mixing Ratio (VMR) in the centre of the plume (SO 2 VMR outside the plume is fixed at the climatological value), with a linear variation of the VMR from the scaling factor from the centre to the edge; • Two parameters to consider the spectral dependence of the optical thickness due to particles; From these retrieved state vector variables, the slant column densities of SO 2 is then derived together with corresponding quantities like H 2 O SCD. The SCDs of SO 2 are then directly comparable to other measurements (see Section 3). The quality of the retrievals is evaluated by performing χ 2 tests. A reduced χ 2 value close to 1 indicates that the residual of the retrieval is at the level of the instrumental noise.

FTIR Instrument
The FTIR is set to a spectral resolution of 0.5 cm −1 . The non-linearity effects are processed by the Bruker software. The raw results are used in the spectra output of the capture software. Measurements of a black body reference temperature (5 • C, 10 • C, 15 • C, 20 • C) enable radiometric calibration (linear least squares) for each wavelength. As expected for ambient temperature emitter, the resulting spectra are much more intense between 8 and 12 µm than between 3 and 5 µm. The radiometric calibration is more accurate in the (8-12) µm spectral range which allows us to compare brightness temperature (BT) spectra with those of Hyper-Cam LWIR instrument.

UV Camera
Volcanic SO 2 observations in the UV were performed by a dual-SO 2 camera system. The configuration consists of two CCD cameras (JAY CM-140GE-UV) with Sony ICX407BLA sensors (1,380 × 1,040 pixels) with sensor quantum efficiency of~7%. The two cameras are mounted side-by-side and aligned such that they view identical scenes. Light is focused onto the sensors using 12 mm focal length quartz lenses.
Two interference band-pass optical filters, with Full Width at Half Maximum (FWHM) of 10 nm and central wavelengths of 310 and 330 nm, respectively, are installed to the lenses to enhance sensitivity to SO 2 and aerosol, and just aerosol, respectively (e.g. [64]).
Field of view is of 28.8 • width × 21.9 • height. The SO 2 camera is controlled by custom-built software, which manages the collection of images at a sampling rate between 0.1 and 0.03 Hz depending on scene luminosity. Shutter time is adjusted automatically to avoid over-exposure.
The system is coupled with a UV spectrometer (Ocean Optics usb2000+), which collects open-path in-plume ultraviolet spectra such as to calibrate the apparent SO 2 column amount observed by the camera system in absolute. Raw UV camera images are processed following the protocols outlined in Burton et al. [64] in which previously stored clear sky images from each camera are used to normalize the measured images from the 310 nm and 330 nm cameras, which are then divided by each other to produce an SO 2 transmittance image. After taking the negative natural logarithm of the transmittance image, it is converted in a differential optical thickness, which is then converted to SO 2 SCD by calibrating the amounts with SO 2 inverted by DOAS ultraviolet spectroscopy [26,65].  (Figure 3d). The noise level was estimated from the experimental hyperspectral cube to take into account both the temporal detector noise and the spatial impact of the image processing. We find a Noise Equivalent Spectral Radiance (NESR) of 11 nW/sr/cm 2 /cm −1 , which means a Noise Equivalent Temperature Difference (NETD) of 880 mK at 2100 cm −1 for a blackbody at 20 • C. Although the NESR is quite constant in the 3.7-4.8 µm spectral band, the NETD increases at higher wavenumbers (lower wavelengths), as it can be seen on Figure 3d: this is due to the high dynamic of the Planck's law for ambient temperature in the MWIR (see Figure 3c).

Preliminary Results of SIBI and Filament Imagers
The behavior of this very compact multi spectral IR camera was quite satisfactory [50]. However, the spectral characterization of the instrument (measurement of the map of the Optical Path Differences) was only done in the laboratory at a temperature about 20 • C, while in the field the temperature was below 10 • C. This temperature difference modified the birefringent interferometer, inducing errors on the spectral calculation (roughly, a dilatation of the spectral axis reaching few tens of cm −1 ). Such an issue could be corrected in future experiments either by a more exhaustive characterization of the instrument in the laboratory, or by implementing a thermal control of the birefringent plate, or by performing in situ spectral calibration (with a monochromatic source).
The Filament instrument has a simpler processing chain: it consists in calibrating (as in a classical camera) each small image, taking into account its spectral sensitivity, and then in interpolating the 24 small images on a common grid, with fixed registration parameter measured in the laboratory.
The spectral transmission of Filament filters is plotted on Figure 4a. The spectral resolution (FWHM) is about 200 cm −1 , much wider than the separation between filters. We can also see that the filters exhibit a secondary transmission peak at higher wavenumbers, and their extinction ratio is not zero. Nevertheless, this is not a problem for gas analysis, since the spectral transmission of the filters is fully characterized in the laboratory and could thus be considered when retrieving the data with Remote Sens. 2019, 11, 1175 9 of 23 the radiative transfer code. The main advantage of these filters is the ease of manufacturing an array of filters with different spectral transmission, especially because all filters have the same structure, with the same layer thicknesses, the only difference being on the grating period. In the case of Filament, there are 6 × 4 filters with a 640 × 512 pixels focal plane array. Nevertheless, due to undesirable cross-talk, the final spectral image is limited to 56 × 56 pixels, by 24 spectral bands (Figure 4b). Figure 3 presents different data from SIBI: a raw (interferometric) image of the scene (Figure 3a), a monochromatic image at 4.8 μm (2,060 cm -1 ) (Figure 3b), and few examples of calibrated spectra at entrance pupil plotted either in spectral radiance (Figure 3c) or in brightness temperature (Figure 3d). The noise level was estimated from the experimental hyperspectral cube to take into account both the temporal detector noise and the spatial impact of the image processing. We find a Noise Equivalent Spectral Radiance (NESR) of 11 nW/sr/cm 2 /cm -1 , which means a Noise Equivalent Temperature Difference (NETD) of 880 mK at 2,100 cm -1 for a blackbody at 20 °C. Although the NESR is quite constant in the 3.7-4.8 μm spectral band, the NETD increases at higher wavenumbers (lower wavelengths), as it can be seen on Figure 3d: this is due to the high dynamic of the Planck's law for ambient temperature in the MWIR (see Figure 3c).  The behavior of this very compact multi spectral IR camera was quite satisfactory [50]. However, the spectral characterization of the instrument (measurement of the map of the Optical Path Differences) was only done in the laboratory at a temperature about 20 °C, while in the field the temperature was below 10 °C. This temperature difference modified the birefringent interferometer, inducing errors on the spectral calculation (roughly, a dilatation of the spectral axis reaching few tens of cm -1 ). Such an issue could be corrected in future experiments either by a more exhaustive characterization of the instrument in the laboratory, or by implementing a thermal control of the birefringent plate, or by performing in situ spectral calibration (with a monochromatic source).

Preliminary results of SIBI and Filament imagers
The Filament instrument has a simpler processing chain: it consists in calibrating (as in a classical camera) each small image, taking into account its spectral sensitivity, and then in interpolating the 24 small images on a common grid, with fixed registration parameter measured in the laboratory.
The spectral transmission of Filament filters is plotted on Figure 4a. The spectral resolution (FWHM) is about 200 cm -1 , much wider than the separation between filters. We can also see that the filters exhibit a secondary transmission peak at higher wavenumbers, and their extinction ratio is not zero. Nevertheless, this is not a problem for gas analysis, since the spectral transmission of the filters is fully characterized in the laboratory and could thus be considered when retrieving the data with the radiative transfer code. The main advantage of these filters is the ease of manufacturing an array of filters with different spectral transmission, especially because all filters have the same structure, with the same layer thicknesses, the only difference being on the grating period. In the case of Filament, there are 6 × 4 filters with a 640 × 512 pixels focal plane array. Nevertheless, due to undesirable cross-talk, the final spectral image is limited to 56 × 56 pixels, by 24 spectral bands ( Figure  4b). The next steps for Filament would be to correct the cross-talk between the channels by inserting an array of separation walls between the micro-lens array and the image plane. It may also be useful for gas analysis to replace the filters by newly developed ones, with finer spectral resolution or better rejection ratio, as described by Tardieu   The next steps for Filament would be to correct the cross-talk between the channels by inserting an array of separation walls between the micro-lens array and the image plane. It may also be useful for gas analysis to replace the filters by newly developed ones, with finer spectral resolution or better rejection ratio, as described by Tardieu et al. (2015) [66].

Results of Hyper-Cam LWIR
Examples of several individual spectra obtained with Hyper-Cam LWIR at different locations in the FOV of one of the sequences acquired on 26 June 2015 are presented in Figure 5. In the clear sky spectrum, the ozone signature in the spectral band 1000-1100 cm −1 is clearly visible (blue) and the Mt Etna is visible as a 300 K black body emission (black). By looking at different pixels in the plume (purple; orange; red; green), we can see that the greater the distance from the volcano crater, the greater the atmospheric ozone signature is detectable. Close to the crater (purple spectrum) the plume is opaque, spectral gas emission is no more detectable.
Remote Sens. 2018, 10, x FOR PEER REVIEW 10 of 25 (purple; orange; red; green), we can see that the greater the distance from the volcano crater, the greater the atmospheric ozone signature is detectable. Close to the crater (purple spectrum) the plume is opaque, spectral gas emission is no more detectable. 3.2.1. IR brightness temperature spectra: Hyper-Cam LWIR versus FTIR As a first step, we compared FTIR OPAG 33 and Hyper-Cam LWIR brightness temperature spectra observed simultaneously in the plume. An example of a direct comparison is given in Figure  6b for diluted plume. The single pixel of FTIR instrument (blue circle in Figure 6a) is over plotted on the FOV of the Hyper-Cam LWIR sequence at 07:24 UTC on the 25 June 2015. This pixel corresponds roughly to 930 spectra collected by Hyper-Cam LWIR. A good agreement between the FTIR spectrum and the 930 mean Hyper-Cam LWIR BT spectra is observed between 880 cm -1 and 1,310 cm -1 , especially in the region of O3 1,000-1,100 cm -1 as well as the one of SO2 1,100-1,200 cm -1 . The high spatial resolution of Hyper-Cam LWIR, compared to the one of the FTIR single pixel, induces a standard deviation up to 15 K. This shows the large variability of the particle and gas content in the Mt Etna plume, but also the high sensitivity of the measurements of the Hyper-Cam LWIR instrument. Indeed, the spectral radiometric accuracy of Hyper-Cam LWIR is 0.5 K at a black body temperature of 30 °C [67].

IR Brightness Temperature Spectra: Hyper-Cam LWIR versus FTIR
As a first step, we compared FTIR OPAG 33 and Hyper-Cam LWIR brightness temperature spectra observed simultaneously in the plume. An example of a direct comparison is given in Figure 6b for diluted plume. The single pixel of FTIR instrument (blue circle in Figure 6a) is over plotted on the FOV of the Hyper-Cam LWIR sequence at 07:24 UTC on the 25 June 2015. This pixel corresponds roughly to 930 spectra collected by Hyper-Cam LWIR. A good agreement between the FTIR spectrum and the 930 mean Hyper-Cam LWIR BT spectra is observed between 880 cm −1 and 1310 cm −1 , especially in the region of O 3 1000-1100 cm −1 as well as the one of SO 2 1100-1200 cm −1 . The high spatial resolution of Hyper-Cam LWIR, compared to the one of the FTIR single pixel, induces a standard deviation up to 15 K. This shows the large variability of the particle and gas content in the Mt Etna plume, but also the high sensitivity of the measurements of the Hyper-Cam LWIR instrument. Indeed, the spectral radiometric accuracy of Hyper-Cam LWIR is 0.5 K at a black body temperature of 30 • C [67].

Radiance spectra retrieved with LARA
Examples of three radiance spectra measured and the retrieved results by LARA code for dense plume, diluted plume, and clear sky are displayed in Figure 7. Since the measurements have a spectral resolution of 2 cm -1 the water vapor bands are well defined and well simulated, which prevent us to attribute water vapor emission to SO2 emission, or the opposite.
A good quality of the retrievals is obtained with respectively normalized χ 2 values of 2, 1.6 and 4.6 for dense plume, diluted plume pixel, and clear sky. As a first approach, we reject retrievals with a normalized χ 2 greater than 10. This value corresponds to a mean retrieved error around 10 % on the SO2 SCD [60].

Radiance Spectra Retrieved with LARA
Examples of three radiance spectra measured and the retrieved results by LARA code for dense plume, diluted plume, and clear sky are displayed in Figure 7. Since the measurements have a spectral resolution of 2 cm −1 the water vapor bands are well defined and well simulated, which prevent us to attribute water vapor emission to SO 2 emission, or the opposite.

Radiance spectra retrieved with LARA
Examples of three radiance spectra measured and the retrieved results by LARA code for dense plume, diluted plume, and clear sky are displayed in Figure 7. Since the measurements have a spectral resolution of 2 cm -1 the water vapor bands are well defined and well simulated, which prevent us to attribute water vapor emission to SO2 emission, or the opposite.
A good quality of the retrievals is obtained with respectively normalized χ 2 values of 2, 1.6 and 4.6 for dense plume, diluted plume pixel, and clear sky. As a first approach, we reject retrievals with a normalized χ 2 greater than 10. This value corresponds to a mean retrieved error around 10 % on the SO2 SCD [60].  A good quality of the retrievals is obtained with respectively normalized χ 2 values of 2, 1.6 and 4.6 for dense plume, diluted plume pixel, and clear sky. As a first approach, we reject retrievals with a normalized χ 2 greater than 10. This value corresponds to a mean retrieved error around 10 % on the SO 2 SCD [60].

SO 2 SCD Retrieved from Hyper-Cam LWIR Measurements
The outputs of LARA give us the slant column of SO 2 in molecules/cm 2 in the plume which are then divided by the air density to obtain the slant column density of SO 2 in ppm m. The retrievals of the spectra have been performed for each pixel independently on a computer of the AERIS data centre (www.aeris-data.fr). The full physics retrieval of one image needs roughly a week of calculation considering the configuration of the selected sequence: 320 × 64 pixels, spectral resolution 2 cm −1 and SO 2 spectral window 1,100-1,200 cm −1 . This large computation time must be reduced to be able to use such instrument for monitoring purposes. A pixel classification methodology is under investigation at the moment to reduce the number of retrievals. Nevertheless, a resolution reduction of the images was applied for the sensitivity tests that will be presented in the next part. Figure 8 presents an example of SO 2 SCD (a) and normalized χ 2 (b) images obtained on the 26 June at 08:27 UTC. The selection of the pixels is made according to two threshold values: the error on the retrieved scaling factor on the SO 2 SCD calculated by the code that must be lower than 10% and a normalized χ 2 < 10. The error considered here is the accuracy determined from the full covariance error matrix that includes the measurement noise and the smoothing error as described in Rodgers, 2000 [60]. The instrumental noise is equal to 5E-8 W/cm 2 /sr/cm −1 . We clearly recognize the shape of

SO2 SCD retrieved from Hyper-Cam LWIR measurements
The outputs of LARA give us the slant column of SO2 in molecules / cm 2 in the plume which are then divided by the air density to obtain the slant column density of SO2 in ppm m. The retrievals of the spectra have been performed for each pixel independently on a computer of the AERIS data centre (www.aeris-data.fr). The full physics retrieval of one image needs roughly a week of calculation considering the configuration of the selected sequence: 320 × 64 pixels, spectral resolution 2 cm -1 and SO2 spectral window 1,100-1,200 cm -1 . This large computation time must be reduced to be able to use such instrument for monitoring purposes. A pixel classification methodology is under investigation at the moment to reduce the number of retrievals. Nevertheless, a resolution reduction of the images was applied for the sensitivity tests that will be presented in the next part. Figure 8 presents an example of SO2 SCD (a) and normalized χ 2 (b) images obtained on the 26 June at 08:27 UTC. The selection of the pixels is made according to two threshold values: the error on the retrieved scaling factor on the SO2 SCD calculated by the code that must be lower than 10% and a normalized χ 2 < 10. The error considered here is the accuracy determined from the full covariance error matrix that includes the measurement noise and the smoothing error as described in Rodgers, 2000 [60]. The instrumental noise is equal to 5E-8 W / cm ²/ sr / cm -1 . We clearly recognize the shape of the plume. Values up to 15,000 ppm m are encountered close to the crater (right part of the image) and values in the range of 1,000-4,000 ppm m in the diluted part of the plume (left part of the image).  Normalized χ 2 image of the plume (Figure 8b) shows that most of the spectra obtained are significantly relevant with normalized χ 2 close to 1. Very few pixels with normalized χ 2 greater than 10 exist in the plume. They probably denote emission regions perturbed by the dynamic of the plume and/or the presence of particles (water droplets, sulfate droplets, ash) which affect the IR spectra.
While the normalized χ 2 values obtained are good and the data processing is easy, the SO 2 SCD values we retrieved in dense and diluted plume regions are very high (4 to 10 times higher) compared to those previously published in the literature by measuring in the UV [16,23,68]. Such high values of SO 2 amounts were similarly encountered in the other measurement sequences of the campaign(some other examples are presented in Figure 9) To investigate the dependency of our results on the assumptions made (characteristic of plume parameters) for SO 2 SCDs retrievals, we performed sensitivity tests. Normalized χ 2 image of the plume (Figure 8b) shows that most of the spectra obtained are significantly relevant with normalized χ 2 close to 1. Very few pixels with normalized χ 2 greater than 10 exist in the plume. They probably denote emission regions perturbed by the dynamic of the plume and/or the presence of particles (water droplets, sulfate droplets, ash) which affect the IR spectra.
While the normalized χ 2 values obtained are good and the data processing is easy, the SO2 SCD values we retrieved in dense and diluted plume regions are very high (4 to 10 times higher) compared to those previously published in the literature by measuring in the UV [16,23,68]. Such high values of SO2 amounts were similarly encountered in the other measurement sequences of the campaign(some other examples are presented in Figure 9) To investigate the dependency of our results on the assumptions made (characteristic of plume parameters) for SO2 SCDs retrievals, we performed sensitivity tests.

Sensitivity tests on retrieval parameters
Smekens and Gouhier [69] used the same instrument (Hyper-Cam LWIR) deployed at Stromboli observing a great diversity of spectra within a collected scene. They show that simple geometry cannot be applied to the entire scene with considerable variation from the top to the bottom of the image due to viewing angle. While the geometry of each pixel is considered different in our retrievals, we took assumptions on the geometry of the plume (altitude of the plume centre Z plume and thickness of the plume Tk plume) and on the temperature difference between the plume and the ambient air ΔT plume. We tested the sensitivity of our results against these parameters. Table 2 resumes the different cases considered. For each of the nine sensitivity tests we performed, the values of two parameters were set with the reference configuration values defined in Section 2.4.2, and the third parameter was tested according to three different values presented in Table 2.

Sensitivity Tests on Retrieval Parameters
Smekens and Gouhier [69] used the same instrument (Hyper-Cam LWIR) deployed at Stromboli observing a great diversity of spectra within a collected scene. They show that simple geometry cannot be applied to the entire scene with considerable variation from the top to the bottom of the image due to viewing angle. While the geometry of each pixel is considered different in our retrievals, we took assumptions on the geometry of the plume (altitude of the plume centre Z plume and thickness of the plume Tk plume ) and on the temperature difference between the plume and the ambient air ∆T plume . We tested the sensitivity of our results against these parameters. Table 2 resumes the different cases considered. For each of the nine sensitivity tests we performed, the values of two parameters were set with the reference configuration values defined in Section 2.4.2, and the third parameter was tested according to three different values presented in Table 2. As mentioned in Section 3.2.2.2, the retrieval of one image takes a significant calculation time. To minimize it for the sensitivity tests of the different parameters, we decided to reduce the number of pixels to be retrieved. We have selected one pixel out of ten for the horizontal axis and one pixel out of five for the vertical axis (Figure 10), which represent 384 pixels and a calculation time of a few hours for each tested value.  As mentioned in 3.2.2.2, the retrieval of one image takes a significant calculation time. To minimize it for the sensitivity tests of the different parameters, we decided to reduce the number of pixels to be retrieved. We have selected one pixel out of ten for the horizontal axis and one pixel out of five for the vertical axis (Figure 10), which represent 384 pixels and a calculation time of a few hours for each tested value. For the interpretation of the results of the sensitivity tests, we chose to separate the image vertically in two parts, diluted plume part and dense plume part at horizontal pixel number 195 (white dashed line in Figure 10). The left part of the image contains normalized χ² ≤ 5 and SO2 SCD values up to 7,000 ppm m (diluted plume part) whereas, the right part of the image corresponds to the dense part of the plume where SO2 SCD values can attain 15,000 ppm m. The sensitivity tests correspond respectively to 225 pixels for diluted part and 65 pixels for dense part of the plume. For the interpretation of the results of the sensitivity tests, we chose to separate the image vertically in two parts, diluted plume part and dense plume part at horizontal pixel number 195 (white dashed line in Figure 10). The left part of the image contains normalized χ 2 ≤ 5 and SO 2 SCD values up to 7,000 ppm m (diluted plume part) whereas, the right part of the image corresponds to the dense part of the plume where SO 2 SCD values can attain 15,000 ppm m. The sensitivity tests correspond respectively to 225 pixels for diluted part and 65 pixels for dense part of the plume. Figure 11 shows, for each sensitivity test, the correlation points of SO 2 SCDs for tested parameters: SO 2 SCDs for the reference. All the correlation points obtained have been fitted by a linear regression and the coefficients of determination (R 2 ) were calculated.  respectively. An additional temperature of 5 K and 10 K in the plume would decrease the SO 2 SCD by 30% and 40%, respectively. However, for the diluted plume, it would be surprising that the temperature difference reaches these values because we are far (several hundred meters) from the craters. Regarding sensitivity tests on plume thickness parameter (Figure 11b)  The sensitivity of the plume parameters is significant in the dense plume retrieval of SO 2 SCD and highlights the complexity of this part of the image with the difficulty to perform a retrieval independent of the assumptions made on the plume characteristics.
The impact of plume parameters can be significant close to the crater, but it is small when we are far from the crater. The results of the sensitivity tests in the diluted plume part show the robustness of the retrieval, with a small dependency of the retrieval with respect to plume parameters and normalized χ 2 lower than 5. In the dense part of the plume we have a high dependency of the retrieval to the assumptions made on plume parameters and a higher normalized χ 2 . However, it is important to notice that plume parameters used for the retrieval do not explain the high values of SO 2 SCD highlighted in the previous part.

Discussion
In order to investigate the SO 2 SCD retrieved by Hyper-Cam LWIR and their consistence with data reported in the literature at Mt Etna, we've compared the results obtained in the IR to those by the UV camera collected simultaneously. Figure 12 presents the image of inverted SO 2 SCD obtained by the UV camera collocated and simultaneously to the IR imager. As the UV camera has a wider FOV (28.8 • × 27.9 • ) respect to Hyper-Cam LWIR and in order to compare the results, the UV image was cropped in relation to the IR hyperspectral imager FOV (25.6 • × 5.1 • ). Moreover, the spatial resolution of the UV camera was also resampled to the spatial resolution of the IR imager. Along vertical axis, one pixel in the IR image corresponds to 3.5 pixels in the UV image, and along horizontal axis, one pixel in the IR image corresponds to 3.8 pixels in the UV image. We also did not consider and masked out the pixels that did not satisfy our retrieval quality conditions of the IR SO 2 SCD results. Remote Sens. 2018, 10, x FOR PEER REVIEW 18 of 25 The order of magnitude of the maximum values obtained by the UV camera close to the crater is ~1,500 ppm m and is ~500 ppm m in the diluted plume region. It is well known that remote observation of passive open-path SO2 in volcanic plumes by UV spectroscopy for remote observation are affected by light dilution: not all the photons measured by the spectrometer travelled through the plume, and this effect may be very large due to in-plume scattering [70]. Unlike UV, IR is not affected by scattering. Figure 13 shows the correlation points between SO2 SCD in IR and SO2 SCD in UV where two compact structures are exhibited. Indeed, as previously mentioned in 3.2.3, we consider two parts of the image, corresponding to 12,419 correlation points in the diluted part (black circles) with low SO2 SCD values and 2,455 points in the dense part (blue circles) with high SO2 SCD. For the diluted plume, the SO2 SCD values are roughly in the range of 0-7,000 ppm m for the IR measurements while it corresponds to values in the range of 0-1,500 ppm m in the UV image. Since the correlation points present a compact shape it confirms that both instruments capture the same type of information in the plume. We derived a second order polynomial fit in the diluted plume (green line in Figure 13):

UV SO 2 SCD
1.36 10 1.55 501.27 (1) where y stands for the SO2 SCD in the IR image and x for the SO2 SCD in the UV image. We also derived a linear fit (red line in Figure 13) to roughly give us the coefficient existing between SO2 SCD in the IR and in the UV images. In the diluted plume the obtained slope is 3.57 with a good determination coefficient of 0.83. The SO 2 SCD inverted from UV images show a plume shape similar to that obtained in the IR image (previously presented in Figure 8a). It means that both instruments capture the same type of information.
The order of magnitude of the maximum values obtained by the UV camera close to the crater is 1,500 ppm m and is~500 ppm m in the diluted plume region. It is well known that remote observation of passive open-path SO 2 in volcanic plumes by UV spectroscopy for remote observation are affected by light dilution: not all the photons measured by the spectrometer travelled through the plume, and this effect may be very large due to in-plume scattering [70]. Unlike UV, IR is not affected by scattering. Figure 13 shows the correlation points between SO 2 SCD in IR and SO 2 SCD in UV where two compact structures are exhibited. Indeed, as previously mentioned in 3.2.3, we consider two parts of the image, corresponding to 12,419 correlation points in the diluted part (black circles) with low SO 2 SCD values and 2,455 points in the dense part (blue circles) with high SO 2 SCD. For the diluted plume, the SO 2 SCD values are roughly in the range of 0-7,000 ppm m for the IR measurements while it corresponds to values in the range of 0-1,500 ppm m in the UV image. Since the correlation points present a compact shape it confirms that both instruments capture the same type of information in the plume. We derived a second order polynomial fit in the diluted plume (green line in Figure 13): y = 1.36 × 10 −3 x 2 + 1.55x + 501.27 (1) where y stands for the SO 2 SCD in the IR image and x for the SO 2 SCD in the UV image. We also derived a linear fit (red line in Figure 13) to roughly give us the coefficient existing between SO 2 SCD in the IR and in the UV images. In the diluted plume the obtained slope is 3.57 with a good determination coefficient of 0.83. Results show that UV camera measurements roughly underestimate the SO2 SCD by approximately a factor of 3.6 (slope of the fit).
In the region of high SO2 SCD values, the correlation points highlight that when the IR imager measurements vary from 5,000 ppm m to 10,000 ppm m, the UV camera measurements present only low variation around 1,200 ppm m. In the dense plume, the UV camera cannot capture the variability of SO2, which can probably be explained by light dilution.
It is very difficult to know if all of the SO2 molecules present along the entire thickness of the plume contribute to the UV signal i.e., only the photons interacting with a partial layer of the plume contribute to the UV signal. Moreover, the UV camera measurements may be limited due to the opacity of the plume probably linked to particle concentrations, which may multiple light scattering within the plume increasing the effective path. As a result, the measured amount could either be reduced or enhanced, especially for low solar altitude angles [20,71,72]. The effective path of the UV scattered light through the plume and between the plume and the instrument can cause a significant attenuation of the measured absorbance, especially at greater distance from the plume [73][74][75]. Mori et al. [70] report that the larger effect of the UV scattering appears at the shorter wavelength bands and with increasing distance from a plume. The authors estimate SO2 attenuation was not significant with any wavelength band (< ± 10%) at 0.6 km but was 35-50 % with shorter wavelength band at 2.6 km distance.
Kern et al. [76] performed a detailed study on the accuracy of SO2 column densities and emission rates measured with a UV camera, by applying radiative transfer code associated with scattering effects due to aerosols. They show that taking into account scattering effects due to aerosols induces SO2 SCD values to possibly be 1.5 or 2 times higher, and emission rates to be 20 % and 90% higher than those obtained with a conventional DOAS retrieval, depending on the different SO2 and aerosol content present in the volcanic plume.
Smekens and Gouhier [69] show that clouds in the background could as well affect the measurements, suggesting the use of radiative transfer retrieval .ie. full physic retrievals. Our study follows this recommendation with a state vector in the SO2 spectral range considering H2O, CO2, O3, N2O, CO, CH4, SO2, and plume parameters including spectral dependency of the optical thickness of the plume and considering the contribution of all these parameters to signal along the line of sight from the instrument up to 80 km. Results show that UV camera measurements roughly underestimate the SO 2 SCD by approximately a factor of 3.6 (slope of the fit).
In the region of high SO 2 SCD values, the correlation points highlight that when the IR imager measurements vary from 5,000 ppm m to 10,000 ppm m, the UV camera measurements present only low variation around 1,200 ppm m. In the dense plume, the UV camera cannot capture the variability of SO 2 , which can probably be explained by light dilution.
It is very difficult to know if all of the SO 2 molecules present along the entire thickness of the plume contribute to the UV signal i.e., only the photons interacting with a partial layer of the plume contribute to the UV signal. Moreover, the UV camera measurements may be limited due to the opacity of the plume probably linked to particle concentrations, which may multiple light scattering within the plume increasing the effective path. As a result, the measured amount could either be reduced or enhanced, especially for low solar altitude angles [20,71,72]. The effective path of the UV scattered light through the plume and between the plume and the instrument can cause a significant attenuation of the measured absorbance, especially at greater distance from the plume [73][74][75]. Mori et al. [70] report that the larger effect of the UV scattering appears at the shorter wavelength bands and with increasing distance from a plume. The authors estimate SO 2 attenuation was not significant with any wavelength band (<±10%) at 0.6 km but was 35-50 % with shorter wavelength band at 2.6 km distance.
Kern et al. [76] performed a detailed study on the accuracy of SO 2 column densities and emission rates measured with a UV camera, by applying radiative transfer code associated with scattering effects due to aerosols. They show that taking into account scattering effects due to aerosols induces SO 2 SCD values to possibly be 1.5 or 2 times higher, and emission rates to be 20 % and 90% higher than those obtained with a conventional DOAS retrieval, depending on the different SO 2 and aerosol content present in the volcanic plume.
Smekens and Gouhier [69] show that clouds in the background could as well affect the measurements, suggesting the use of radiative transfer retrieval i.e., full physic retrievals. Our study follows this recommendation with a state vector in the SO 2 spectral range considering H 2 O, CO 2 , O 3 , N 2 O, CO, CH 4 , SO 2 , and plume parameters including spectral dependency of the optical thickness of the plume and considering the contribution of all these parameters to signal along the line of sight from the instrument up to 80 km.
Similarly, uncertainty exists for the IR imager measurements, whereas we parameterize particle emission to retrieved SO 2 SCD, the opacity of the dense plume could also affect the IR signal detected overestimating retrieved SO 2 SCD. The sensitivity tests we performed highlighted the complexity of the dense plume part and additional work is needed to give robust measurements of SO 2 SCD in this region. Joint explicit particle and gas radiative transfer model could be applied if the instrument could attain larger wavelength [77] to characterize the nature and the distribution of particles for dense plume.

Conclusions
During the IMAGETNA campaign at Mt Etna (Pizzi De Neri observatory) we have deployed three infrared imagers: SIBI, Filament and the IR hyperspectral camera Hyper-Cam LWIR from TELOPS Company. Combined with a FTIR single pixel and a UV camera we have investigated the ability of such instruments to measure SO 2 degassing from volcanic plume and to provide further constraints on the uncertainty in open-path passive UV scattering SO 2 observation.
SIBI IR imager operates in the spectral domain 3.7-4.8 µm with a spectral resolution about 14 cm −1 . It is easy to implement and the post processing to retrieve spectra have been performed. The limitation of SIBI during IMAGETNA campaign mainly comes from the stability of the interferometer with respect to the temperature. This could have been prevented by an in situ calibration, or a thermal control of the interferometer. Further efforts are consequently needed to evaluate its ability to use it for other gases measurements such as CO 2 .
Filament infrared multispectral imager has 24 bands between 3 µm and 5 µm and uses filters leading to a spectral resolution of about 200 cm −1 . The next developments would be to correct the cross-talk between the channels. An option is to replace the filters to have a finer spectral resolution or a better rejection ratio, and to improve the cryogenic system in order to facilitate use during campaign.
For the first time, an IR hyperspectral camera (Hyper-Cam LWIR) has been deployed with a data processing using a line-by-line radiative transfer model to retrieve as a first step SO 2 SCD. Combined with the widely used single pixel FTIR, we have checked the radiance spectra intensity obtained by the Hyper-Cam LWIR. A good agreement was observed. The LARA line by line radiative transfer model has been applied to retrieve SO 2 SCD images from the Hyper-Cam LWIR. The SO 2 SCD obtained are higher than those reported in the literature using UV camera at Mt Etna.
The sensitivity tests on the plume characteristics assumptions highlight a limited impact on SO 2 SCD in the diluted part of the plume.
A detailed comparison of simultaneous spatial and temporal measurements from the UV camera and the IR imager have been done. Not expected very high SO 2 SCD values are obtained, with strong differences whatever the pixel considered. While the IR retrievals were easily performed with good quality (most of the normalized χ 2 are lower than five, and only a few bad pixels were detected) SO 2 SCD IR measurements are greater than SO 2 SCD UV measurements. In diluted plume condition, we have obtained a compact UV:IR correlation. This result shows that both instruments capture the same plume shape and that a polynomial correlation fit can be deduced. The UV seems to underestimate the SO 2 SCD by a factor of roughly 3.6. In the dense plume; the opacity of the plume, probably due to numerous particles, leads to a very small SO 2 SCD variability in the UV in comparison to the one in the IR. While IR measurements are less bothered by particles, this dense plume region needs to be investigated in more detail.
The light dilution in passive scattered UV remote sensing of SO 2 in volcanic plumes is complex, and mainly relies on distance between plume and instruments and aerosol burden. The quantification of effective SO 2 amounts is still an opened question. The infrared limb measurements associated with a full physic retrieval could help to quantify this phenomenon.
As a conclusion, this study highlights the potential of IR application in volcanic gas measurements at active volcanoes to investigate plume emissions all day long and that additional species can also potentially be simultaneously accessible. This study, with detailed IR and UV SO 2 measurements comparisons, highlights the need to quantify the scattering effect on UV camera remote sensing measurements for the characterization of volcanic plume emissions.