the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Quantitative imaging of volcanic SO2 plumes using Fabry–Pérot interferometer correlation spectroscopy
Christopher Fuchs
Nicole Bobrowski
Ulrich Platt
We present first measurements with a novel imaging technique for atmospheric trace gases in the UV spectral range. Imaging Fabry–Pérot interferometer correlation spectroscopy (IFPICS) employs a Fabry–Pérot interferometer (FPI) as the wavelength-selective element. Matching the FPI's distinct, periodic transmission features to the characteristic differential absorption structures of the investigated trace gas allows us to measure differential atmospheric column density (CD) distributions of numerous trace gases with high spatial and temporal resolution. Here we demonstrate measurements of sulfur dioxide (SO2), while earlier model calculations show that bromine monoxide (BrO) and nitrogen dioxide (NO2) are also possible. The high specificity in the spectral detection of IFPICS minimises cross-interferences to other trace gases and aerosol extinction, allowing precise determination of gas fluxes. Furthermore, the instrument response can be modelled using absorption cross sections and a solar atlas spectrum from the literature, thereby avoiding additional calibration procedures, e.g. using gas cells. In a field campaign, we recorded the temporal CD evolution of SO2 in the volcanic plume of Mt. Etna, with an exposure time of 1 s and 400×400 pixel spatial resolution. The temporal resolution of the time series was limited by the available non-ideal prototype hardware to about 5.5 s. Nevertheless, a detection limit of 2.1×1017 molec cm−2 could be reached, which is comparable to traditional and much less selective volcanic SO2 imaging techniques.
- Article
(10506 KB) - Full-text XML
-
Supplement
(17388 KB) - BibTeX
- EndNote
Ground-based imaging of atmospheric trace gas distributions has a great potential to give new insights into mixing processes and chemical conversion of atmospheric trace gases by allowing their observation at high spatio-temporal resolution. Whereas present space-borne trace gas imaging provides daily global coverage with a spatial resolution of a few kilometres (e.g. Veefkind et al., 2012), ground-based observation can potentially reach a spatial resolution in the order of metres and a temporal resolution in the single digit Hz range. Such techniques in particular allow the investigation of trace gas distributions with strong gradients and short timescale chemical conversions.
There are several approaches for imaging trace gas distributions using scattered sunlight in the UV–Vis wavelength range (see, for example, Platt and Stutz, 2008; Platt et al., 2015): an image can be scanned pixel by pixel with a telescope, and recorded spectra are evaluated to determine the trace gas column density (whiskbroom approach). Alternatively, with a more complex optics and a two-dimensional detector, one detector dimension of the spectrograph can be used for spatially resolving an image column. Column by column (or pushbroom) scanning then resolves an image. The high spectral resolution of the spectrograph-based techniques allows the accurate and simultaneous identification of several trace gases; however, the light throughput and the scanning process severely limit the temporal resolution. A third approach applies tunable filters to resolve the trace gas spectral features, e.g. acousto-optical tunable filters (Dekemper et al., 2016), as wavelength-selective elements for an entire image frame. The application of tunable filters can have high spectral resolution and hence high trace gas selectivity; however, due to limited light throughput the temporal resolution lies in the order of minutes. A fourth imaging technique uses a small number (typically two) of wavelength channels selected by static filters, e.g. interference filters (Mori and Burton, 2006). This approach can be quite fast with a temporal resolution in the order of seconds; the trace gas selectivity, however, strongly depends on the correlation of trace gas absorption with the wavelength-selective elements and usually is rather marginal.
Fabry–Pérot interferometers (FPIs) exhibit a periodic spectral transmission pattern which can be matched to periodic spectral features (typically due to rotational or vibrational structures of electronic transitions) of the trace gas absorption, thereby yielding very high correlation for some trace gases. Imaging Fabry–Pérot interferometer correlation spectroscopy (IFPICS) thus essentially combines the advantage of fast image acquisition with selective spectral identification of the target trace gas. IFPICS was proposed by Kuhn et al. (2014) and discussed in Platt et al. (2015) for volcanic SO2. Kuhn et al. (2019) demonstrated the feasibility with a 1 pixel prototype for volcanic SO2 and evaluated its applicability to other trace gases.
Here we present first imaging measurements (at a resolution of 400×400 pixels, 1 s exposure time) performed with IFPICS and confirm its high selectivity and sensitivity. A prototype instrument for SO2 was tested at Mt. Etna volcano, Italy, showing a noise-equivalent signal between 2.1×1017–5.5×1017 . Furthermore, we show that the instrument response can be modelled and thereby intrinsically calibrated, using a solar atlas spectrum and literature trace gas absorption cross sections.
Existing interference-filter-based SO2 cameras used, for example, for the quantification of volcanic trace gas emission fluxes into the atmosphere (Mori and Burton, 2006; Bluth et al., 2007; Kern et al., 2015a), exhibit strong cross-interferences to aerosol scattering extinction and other trace gases (Lübcke et al., 2013; Kuhn et al., 2014). Furthermore, these techniques require in-field calibration. Besides the thereby induced systematic errors that propagate into the emission flux quantification, the detection limit is mostly determined by these cross-interferences. Thus, the applicability of the technique is limited to strong emitters with respective plume and weather conditions. The much higher selectivity of IFPICS largely extends the range of applicable conditions (e.g. to ship emissions and weaker emitting volcanoes) and significantly reduces the systematic errors. Furthermore, the extension of the technique to other trace gases, e.g. bromine monoxide (BrO), formaldehyde (HCHO), or nitrogen dioxide (NO2), can give new important insights into short-scale chemical conversion processes in the atmosphere.
Similarly to the SO2 camera principle (e.g. Mori and Burton, 2006; Bluth et al., 2007), IFPICS uses an apparent absorbance (AA) , i.e. the difference between two measured optical densities τA and τB, to quantify the column density (CD) , i.e. the integrated concentration c of the measured gas along a light path L for each pixel of the image. The AA is calculated from two (or more) spectral settings that yield a maximum correlation difference to the gas absorption spectrum. Ideally the periodicity of the FPI fringes is matched to periodic spectral absorption features as shown in Fig. 1 for SO2. For IFPICS we use two spectral settings A and B. Setting A exhibits on-band absorption, where the FPI transmission maxima coincide with the SO2 absorption maxima and hence correlate with the differential absorption structures of SO2. Setting B uses an off-band position, where the FPI transmission maxima anti-correlate with the differential SO2 absorption structures (see Fig. 1). The spectral separation between setting A and B is thereby reduced by a factor of ≈30 (in the case of SO2) to only ≈0.5 nm in contrast to ≈10–15 nm for traditional SO2 cameras (see Lübcke et al., 2013; Kern et al., 2015a), which minimises broadband interferences due to, for example, scattering and extinction by aerosols or other absorbing gases. This application of an FPI is similar to approaches reported by Wilson et al. (2007) and Vargas-Rodríguez and Rutt (2009), for the detection of carbon monoxide, carbon dioxide, and methane in the infrared spectral range.
By measuring the optical density and in both spectral settings A and B respectively, the relation between the AA with the CD S is given by
where IA and IB denote the radiances with the presence of the target trace gas in the absorption light path and I0,A and I0,B the radiance without. The absorber-free reference radiances I0,A and I0,B can be determined from, for example, a reference region within the image. The differential weighted effective trace gas absorption cross section becomes independent of S for small AAs (). At higher AAs, saturation effects occur due to the non-linearity of Lambert–Beer's law; however knowledge of the absorption cross sections, the background radiation spectrum, and the instrument transmission allows to be calculated for arbitrary CDs S using a numerical model.
2.1 Instrument model
The AA is modelled for given target trace gas CDs S by simulating the incoming radiances IA and IB and I0,A and I0,B. As incident radiation, a high-resolution, top of atmosphere (TOA) solar atlas spectrum I0,TOA(λ) is used according to Chance and Kurucz (2010). The TOA spectrum is scaled by the wavelength λ−4 approximating a Rayleigh scattering atmosphere. Since our measurement wavelength range, of 304–313 nm for SO2, overlaps with absorption of ozone (O3), the TOA spectrum is corrected for the stratospheric O3 absorption by multiplying all intensities by Lambert–Beer's term , where denotes the total atmospheric ozone slant column density, e.g. according to the TEMIS database (Veefkind et al., 2006), and the O3 absorption cross section according to Serdyuchenko et al. (2014). This yields the scattered skylight radiance I0(λ) as follows:
indicated in Fig. 2. Based on I0(λ), the radiances measured by the instrument for the two respective spectral settings are calculated with the absorption of trace gases and the spectral instrument transfer function Tinstr(λ). The investigated target trace gas j (in this work SO2) and potentially interfering trace gas species k (in this work O3) are added according to Lambert–Beer's law. In the following we use the index i, denoting the FPI settings A and B, respectively. The quantity I0,i thereby denotes the reference radiance excluding the target trace gas j from the light path (see Fig. 2).
The spectral instrument transfer functions Tinstr,i(λ) for the two spectral settings,
consists of the measured band pass filter (BPF) transmission spectrum TBPF(λ), the spectral (i.e. wavelength-dependent) quantum efficiency Q(λ) of the detector, and a spectral loss factor η(λ) of the employed optical components (e.g. by reflections). Considering only a single, parallel beam of light traversing the instrument, the FPI transmission spectrum TFPI,i(λ) is defined by the Airy function (Perot and Fabry, 1899)
with the light beam incidence angle αi for the two spectral settings onto the FPI, the FPI mirror separation d, the refractive index n of the medium inside the FPI, and the FPI reflectivity R (see Table 1, Figs. 1 and 3d). The FPI used in this work is static and air-spaced, meaning d, n, and R are fixed. Hence, the incidence angle αi is the exclusive free parameter available to tune the FPIs transmission spectrum TFPI,i between settings i=A and i=B respectively. The change in αi is achieved by tilting the FPI optical axis with respect to the imaging optical axis (see Sect. 2.2).
However, in reality a spectral setting will always impinge with a range of incidence angles onto the FPI. In this work we assume cone-shaped light beams, with half cone opening angles ωc, where the entire cone can be tilted by αi relative to the normal of the FPI mirrors (see Fig. 3d). From this assumption, it follows that the incidence angles αi are distributed over a cone with the incidence angle distribution , where ϑ and φ are the polar and azimuth angles, respectively. Hence, the single-beam FPI transmission spectrum TFPI,i(λ) of Eq. (6) is extended by a weighted average over , giving the effective FPI transmission spectrum :
Thereby, N(γ(αi,ωc)) denotes the weighting function with given by the integral in Eq. (7), excluding the integrand TFPI,i itself, ϑ the polar angle, and φ the azimuth angle of the spherical integration within boundaries defined by the tilted cone-shaped light beams. For example, for a non-tilted FPI (αi=0) the integration boundaries are and ; for a tilted FPI, however, the transformation of is more complex and requires several case analyses.
The incidence angle distribution γ(αi,ωc) will affect the shape of the FPI transmission spectrum by decreasing the effective finesse F of the FPI, leading to a blurring of the FPI fringes (see Fig. 1).
2.2 The IFPICS prototype
The IFPICS prototype is a newly developed instrument, designed to function under harsh environmental conditions in remote locations like, for example, in proximity to volcanoes. Hence, the prototype is designed to be small, with dimensions of , lightweight at 4.8 kg (see Fig. 3a), and a power consumption <10 W; thus it can be battery-operated for several hours. A 2D UV-sensitive CMOS sensor with 2048×2048 pixel resolution (SCM2020-UV provided by EHD imaging) is used to acquire images. The sensor is operated in 4×4 binning mode, yielding a final image resolution of 512×512 pixels. However, we found that the software of the SCM2020-UV image sensor does not allow sufficiently precise triggering. Therefore ≈0.6 s is lost in each image acquisition, which severely limits the operation of the IFPICS camera. Replacement of the sensor by a scientific-grade UV detector array will solve this problem in future studies.
The internal camera optics is highly modular and easily adjustable. The IFPICS prototype employs an image side telecentric optical setup as proposed in Kuhn et al. (2014, 2019). A photograph and a schematic drawing are shown in Fig. 3. An aperture and a lens (lens 1) parallelise incoming light from the imaging field of view (FOV) before it traverses the FPI and the BPF. A second lens (lens 2) focusses the light onto the 2D UV-sensitive sensor. Thereby, in good approximation, all the pixels of the image experience the same spectral instrument transfer function Tinstr,i(λ) for the two wavelength settings.
The FPI is the central optical element of the IFPICS prototype and is implemented as static air-spaced etalon with fixed d, n, and R (provided by SLS Optics Ltd.). The mirrors are separated using ultra-low expansion glass spacers to maintain a constant mirror separation d and parallelism over the large clear aperture of 20 mm, even in variable environmental conditions. In order to tune the spectral transmission between setting A and B, a variation of the incidence angle α is applied. The FPI can be tilted within the parallelised light path using a stepper motor. The stepper motor has a resolution of 0.9∘ per step, is equipped with a planetary gearbox (reduction rate 1∕9), and is operated in micro-stepping mode (1∕16), resulting in a resolution of 0.00625∘ per motor step. The time required for tilting between our settings A and B is ≈0.15 s. We favour the approach of tilting the FPI over changing internal physical properties like, for example, the mirror separation d by piezoelectric actuators, as it keeps simplicity, robustness, and accuracy high for measurements under non-laboratory conditions. However it needs to be considered that the tilting of the FPI will generate a linear shift between the respective images acquired in setting A and B, requiring an alignment in the evaluation process.
The half cone opening angle ωc is determined by the entrance aperture a and the focal length f of lens 1 and can be calculated by . The physical properties of the optical components and the instrument are listed in Table 1 and were mostly chosen according to the dimensioning assumed in the calculations of Kuhn et al. (2019).
The FPI design with fixed d, n, and R (see Fig. 3c) in particular is chosen to inherently generate a transmission spectrum matching the differential absorption structures of SO2. This includes the basic idea that the untilted FPI () already matches the on-band position A. In our case, however, the manufacturing accuracy of d lies within one free spectral range (≈2 nm for SO2) so that corresponds to an off-band (B) position (), and the on-band (A) position () is reached by a small tilt of . The basic advantages of using small incident angles αi are that they keep the spread of the incidence angle distribution γ(αi,ωc) (see Sect. 2.1) low and thereby retain the FPI's effective finesse F high (since the reflectivity R of the FPI mirror coating is somewhat dependent on the angle of incidence so is the finesse F). This leads to a much weaker blurring of the FPI fringes in the FPI transmission spectrum resulting in a higher sensitivity of the instrument (see Fig. 1). With the prototype setup, however, we encountered disturbing reflections for low FPI incidence angles. For that reason we used the subsequent correlating order of the FPI transmission with for an on-band and for an off-band setting (see Table 1), thereby making a compromise between sensitivity and accurate evaluable images.
3.1 Measurements at Mt. Etna, Italy
First measurements with the prototype described above were performed at the Osservatorio Vulcanologico Pizzi Deneri (lat. 37.766, long. 15.017; 2800 ) at Mt. Etna, on 21 and 22 July 2019. The physical properties of the IFPICS prototype and the FPI tilt angles αi for tuning between on-band i=A and off-band setting i=B were selected according to Table 1. The tilt of the FPI generates a linear shift between the recorded on-band and off-band images on the detector and accounts for 6 pixels using tilt angles αi. This shift needs to be corrected before cross-evaluating images recorded in setting A and B. The exposure time was set to 1 s for all measurements, and 4×4 binning (total spatial resolution of 512×512 pixels) was applied for all acquired images.
3.2 Validation of the instrument model
To quantify the accuracy of our model, two SO2 gas cells were measured with the IFPICS prototype and by differential optical absorption spectroscopy (DOAS; see Platt and Stutz, 2008), on 21 July 2019, 11:10–11:20 CET. The sky was used as light source with a constant viewing angle (10∘ elevation, 270∘ N azimuth) in a plume-free part of the sky. To enhance the image quality, a flat-field correction is used, compensating pixel-to-pixel variations in sensitivity. The flat-field correction requires the acquisition of dark and flat-field images. The dark images are determined by the arithmetic mean over five images with no light entering the IFPICS instrument, and the flat-field images are obtained by the arithmetic mean over five images acquired in a plume-free sky region. The flat-field images thereby directly include the reference measurement I0,i, making a later correction for the atmospheric background unnecessary. In the same viewing direction, Ii is measured for each gas cell and FPI setting i in order to calculate the AA according to Eq. (1). Figure 4 shows the gas cell measurements (red) including uncertainties (error bars, 1σ). The uncertainties directly arise from the errors of the DOAS measurement and due to variations in optomechanical settings of the IFPICS prototype.
The instrument model (Eqs. 2–7) was used to calculate the IFPICS AA from a given SO2 CD . The model parameters are mostly fixed by the IFPICS prototype optics, as given in Table 1. The remaining parameter, the atmospheric O3 slant column density (see Eq. 2), is calculated in a geometric approximation using the solar zenith angle (SZA) and vertical O3 column density (VCD), which both are location-, date-, and time-dependent. They were (according to the solar geometry calculator by NOAA, 2020) and VCD (according to the TEMIS database; Veefkind et al., 2006). The VCD can be treated to be approximately constant over the period of a day.
The output of the instrument model (drawn, black) for an SZA of 53∘ is shown in Fig. 4. The model uncertainty (shaded grey) is determined by a root mean square over the errors in the output by individually varying the input parameters within their stated uncertainties. The thus calculated calibration function using the instrument model matches the SO2 gas cells' validation measurement within the range of confidence. The model nicely describes the flattening of the AA–CD relation for high CDs (up to ), which originates from the CD dependence of (see Eq. 1).
To show the impacts of the SZA on the instrument model, the model output is also calculated for three other SZAs while keeping the other parameters constant. The model output is shown in Fig. 4 for an SZA of 80∘ (dashed–dotted black) for early morning and late afternoon conditions, an SZA of 70∘ (dashed black) for morning and afternoon conditions, and an SZA of 25∘ (dotted black) for noon conditions. High SZAs lead to an increase of stratospheric O3 absorption, which alters the spectral shape of the scattered skylight radiance I0(λ) (see Eq. 2) which is used in the forward model. In other words, for high O3 absorption, lower wavelength radiance, where the differential SO2 absorption features are stronger, will contribute less to the integrated radiances Ii, I0,i (Eqs. 3, 4). The thereby induced SZA dependence of the sensitivity can easily be accounted for in the model. Note that this influence of strong O3 absorption only occurs at our chosen wavelength range for the SO2 measurement. When applying IFPICS to other trace gases, e.g. BrO or NO2, at higher wavelength, this effect will be negligible.
3.3 Results of the field measurements
Volcanic plume measurements were performed on 22 July 2019 at 08:50–09:10 CET. The instrument was pointing towards the plume of Mt. Etna's South-East Crater with a constant viewing direction (azimuth 204∘ N, elevation 5∘; see Fig. 5). The wind direction was ≈5∘ N, with a velocity of (wind data from UWYO, 2020). Hence, the plume was partly covered by the crater flank. The frame rate during the measurement was 0.2 Hz for a pair (IA and IB) of images.
The flat-field correction was performed as described in Sect. 3.2, using the arithmetic mean over 10 dark images and five flat-field images, obtained in a plume-free sky region. An exemplary set of volcanic plume SO2 images, obtained with the IFPICS instrument in on-band setting IA and off-band setting IB, are shown in Fig. 6. Further images of IA and IB are shown in Appendix A. The circular shape of the retrieved image arises from the FPI's circular clear aperture limiting the imaging FOV.
The IFPICS SO2 AA is calculated pixel-wise according to Eq. (1) from IA and IB. For the conversion into SO2 CD , the forward instrument model (Eqs. 2–7) is inverted by least-squares fitting of a fourth-order polynomial to the calculated CD relation . The model input parameters of the instrument are shown in Table 1. The SZA during the time of the measurement is (NOAA, 2020), with a VCD of 335±5 DU (according to the TEMIS database; Veefkind et al., 2006). The retrieved calibration function is
with , , , and in units of molec cm−2 respectively, with x0 fixed to zero. This approximation yields an average relative deviation of 0.007 % for from the modelled value, with a maximum relative deviation of 0.08 % for small SO2 CDs.
An evaluated image of the volcanic plume SO2 CD distribution corresponding to the intensities shown in Fig. 6 is shown in Fig. 7. Further evaluated CD distribution images of the same time series are presented in Appendix A. A time series of the plume evolution is visualised in a flip-book in the Supplement.
The volcanic plume of Mt. Etna's South-East Crater is clearly visibly and reaches SO2 CDs higher than . The atmospheric background is and was determined by the arithmetic mean over a plume-free area within the evaluated image (white square, 100×100 pixels, in Fig. 7a). Since the is determined from an evaluated CD distribution image, it accounts for the residual signal in between the direction of the volcanic plume and the direction of the flat-field images used in the evaluation. The was subtracted from the displayed image in the final step of the evaluation. The similar plume-free area (white square, 100×100 pixels, in Fig. 7a) is further used to give an estimation for the SO2 detection limit of the IFPICS prototype by calculating the 1σ pixel–pixel standard deviation. The obtained detection limit for an exposure time of 1 s is 5.5×1017 molec cm−2, given by the noise-equivalent signal. The measurements were performed in the morning with an SZA of 78∘ and therefore reduced sensitivity and under relatively low light conditions. For decreasing SZA, the sensitivity will increase according to Fig. 4, and the increasing sky radiance will reduce the photon shot noise. In other words, the gas cell measurements (taken at SZA of 53∘, with approximately twice the sky radiance compared to SZA of 78∘) show a detection limit of for an exposure time of 1 s. For ideal measurement conditions (lowest SZA, highest sky radiance), the detection limit will be further improved.
After the proof of concept, showing the capability of IFPICS to determine SO2 CD images, it is possible to determine fluxes from a CD image time series. In particular, if the series allows us to trace back individual features in consecutively recorded images, it can be used to directly determine the plume velocity using the approach of cross-correlation (e.g. McGonigle et al., 2005; Mori and Burton, 2006; Dekemper et al., 2016) or optical flow algorithms (e.g. Kern et al., 2015b) and to determine the plume propagation direction (e.g. Klein et al., 2017). However, the viewing geometry on the day of our measurement was unfavourable as it was not possible to reach another measurement location due to a lack of infrastructure. The plume propagation direction and central line of sight show an inclination of 19∘ only, resulting in high pixel contortions, especially for pixels close to the edges of the FOV. Further, significant parts of the plume are covered by the crater flank due to its propagation direction. For the sake of completeness, we would like to give a rough estimate of the SO2 flux obtained from our data.
The SO2 flux is determined by integrating the SO2 CD along a transect through the volcanic plume and subsequent multiplication by the wind velocity perpendicular to the FOV direction; however due to the viewing geometry issues, we will use external wind data (direction: 5∘; velocity data from UWYO, 2020) for the calculation. As the camera pixel size is finite, the integral is replaced by a discrete summation over the pixel n:
including the perpendicular wind velocity v⟂, the SO2 CD , and the pixel extent hn. The perpendicular wind velocity can directly be calculated from geometric considerations (see Fig. 5a), accounting to . To determine the pixel extent, the distance between the volcanic plume and the location of measurement is required. In the centre of the FOV, this distance is ≈3500 m, yielding hn≈2.7 m. To keep the impact of pixel contortions low, the plume transect is located centrally in the FOV at column 250 and ranges from rows n=230 to 330. Using these quantities, we retrieve a mean SO2 mass flux for the measurement of for the investigated plume of the South-East Crater, which is comparable to previous flux measurements of the South-East Crater (Aiuppa et al., 2008; D′Aleo et al., 2016). Nevertheless, the flux should be regarded as a lower limit, since the plume was covered by the crater flank to an unknown extent.
By imaging and quantifying the SO2 distribution in the volcanic plume of Mt. Etna, we successfully demonstrate the feasibility of the IFPICS technique proposed by Kuhn et al. (2014). We were able to unequivocally resolve the dynamical evolution of SO2 in a volcanic plume with a high spatial and temporal resolution (400×400 pixels, 1 s integration time, 4×4 binning). The retrieved detection limit for the SO2 measurement is . The detection limit however varies with the SZA and can reach values below under ideal conditions, comparable to traditional SO2 imaging techniques (see Kern et al., 2015a). Also, the imaging technique lends itself to the determination of gas fluxes, and we obtained an SO2 mass flux of for Mt. Etna's South-East Crater plume. However, due to unfavourable conditions in the viewing geometry, the retrieved flux should be treated as a lower limit. In general, it is possible to apply optical flow algorithms on image series acquired under more ideal viewing geometry conditions (e.g. Kern et al., 2015b). These allow the plume velocity and angle to be determined between the observation direction and plume propagation direction in order to retrieve accurate SO2 fluxes (e.g. Klein et al., 2017).
The specific spectral detection scheme of IFPICS allows a numerical instrument model to be used to directly convert the measured AA into CD S distributions. This inherent calibration method makes in-field calibrations methods, e.g. by gas cells, unnecessary. The accuracy of the instrument model could be demonstrated using SO2 cells with a known CD, determined by simultaneous DOAS measurements.
Our IFPICS instrument is still an early stage prototype. The optics employed is highly modular, allowing easy adjustments even outside a laboratory. The physical dimensions of <10 L and <5 kg and the low power consumption of <10 W, combined with the fact that no maintenance and in-field calibration are needed, make it already a close to ideal field instrument. Furthermore, the temporal resolution of the instrument can further be increased by replacing the sensor employed as it does not allow for time-optimised control of image acquisition.
Compared to traditional SO2 cameras, the minimised cross-interferences to broadband plume extinction increase the selectivity and thus should allow to apply the IFPICS technique to much weaker SO2 sources. Furthermore, the expected smaller interference to broadband effects in comparison to traditional SO2 imaging techniques should allow the range of meteorological conditions acceptable for field measurement to be extended (see Kuhn et al., 2014).
The demonstrated IFPICS technique is not limited to the detection of SO2. In general the technique is applicable to numerous further trace gases which show a distinct pattern (ideally periodic) in their absorption spectrum (see Kuhn et al., 2019). In the case of volcanic emissions, detectable trace species are, for example, bromine monoxide (BrO) or chlorine dioxide (OClO). Beyond volcanic applications IFPICS could be used to investigate, for example, air pollution by measuring nitrogen dioxide (NO2) or formaldehyde (HCHO).
Further evaluated images of the time series acquired on 22 July 2019 at 08:50–09:10 CET at Mt. Etna, Italy, are shown in Fig. A1. The evaluation procedure is analogous to the routine explained in Sect. 3.3. The time difference between a set (1–4) of images accounts for ≈120 s and allows us to trace back plume dynamics.
The data can be obtained from the authors upon request.
The supplement related to this article is available online at: https://doi.org/10.5194/amt-14-295-2021-supplement.
JK, NB, and UP developed the question of research. JK, NB, and CF conducted the field campaign. JK and CF developed the instrument model. CF designed, constructed, and characterised the instrument, evaluated the data, and wrote the manuscript, with all authors contributing by revising it within several iterations.
The authors declare that they have no conflict of interest.
We would like to thank SLS Optics Ltd. for sharing their expertise in designing and manufacturing the etalons. Support from the Deutsche Forschungsgemeinschaft (project DFG PL 193/23-1) is gratefully acknowledged. We also thank Emmanuel Dekemper and Toshiya Mori for their valuable and constructive reviews.
The article processing charges for this open-access publication were covered by the Max Planck Society.
This paper was edited by Andreas Hofzumahaus and reviewed by Emmanuel Dekemper and Toshiya MORI.
Aiuppa, A., Giudice, G., Gurrieri, S., Liuzzo, M., Burton, M., Caltabiano, T., McGonigle, A. J. S., Salerno, G., Shinohara, H., and Valenza, M.: Total volatile flux from Mount Etna, Geophys. Res. Lett., 35, L24302, https://doi.org/10.1029/2008gl035871, 2008. a
Bluth, G., Shannon, J., Watson, I., Prata, A., and Realmuto, V.: Development of an ultra-violet digital camera for volcanic SO2 imaging, Journal of Volcanol. Geoth. Res., 161, 47–56, 2007. a, b
Bogumil, K., Orphal, J., Homann, T., Voigt, S., Spietz, P., Fleischmann, O., Vogel, A., Hartmann, M., Kromminga, H., Bovensmann, H., Frerick, J., and Burrows, J.: Measurements of molecular absorption spectra with the {SCIAMACHY} pre-flight model: instrument characterization and reference data for atmospheric remote-sensing in the 230–2380 nm region, J. Photoch. Photobio. A, 157, 167–184, https://doi.org/10.1016/S1010-6030(03)00062-5, 2003. a
Chance, K. and Kurucz, R.: An improved high-resolution solar reference spectrum for earth's atmosphere measurements in the ultraviolet, visible, and near infrared, J. Quant. Spectrosc. Ra., 111, 1289–1295, https://doi.org/10.1016/j.jqsrt.2010.01.036, 2010. a
D′Aleo, R., Bitetto, M., Donne, D. D., Tamburello, G., Battaglia, A., Coltelli, M., Patanè, D., Prestifilippo, M., Sciotto, M., and Aiuppa, A.: Spatially resolved SO2 flux emissions from Mt Etna, Geophys. Res. Lett., 43, 7511–7519, https://doi.org/10.1002/2016gl069938, 2016. a
Dekemper, E., Vanhamel, J., Van Opstal, B., and Fussen, D.: The AOTF-based NO2 camera, Atmos. Meas. Tech., 9, 6025–6034, https://doi.org/10.5194/amt-9-6025-2016, 2016. a, b
Kern, C., Lübcke, P., Bobrowski, N., Campion, R., Mori, T., Smekens, J.-F., Stebel, K., Tamburello, G., Burton, M., Platt, U., and Prata, F.: Intercomparison of SO2 camera systems for imaging volcanic gas plumes, J. Volcanol. Geoth. Res., 300, 22–36, https://doi.org/10.1016/j.jvolgeores.2014.08.026, 2015a. a, b, c
Kern, C., Sutton, J., Elias, T., Lee, L., Kamibayashi, K., Antolik, L., and Werner, C.: An automated SO2 camera system for continuous, real-time monitoring of gas emissions from Kīlauea Volcano′s summit Overlook Crater, J. Volcanol. Geoth. Res., 300, 81–94, https://doi.org/10.1016/j.jvolgeores.2014.12.004, 2015b. a, b
Klein, A., Lübcke, P., Bobrowski, N., Kuhn, J., and Platt, U.: Plume propagation direction determination with SO2 cameras, Atmos. Meas. Tech., 10, 979–987, https://doi.org/10.5194/amt-10-979-2017, 2017. a, b
Kuhn, J., Bobrowski, N., Lübcke, P., Vogel, L., and Platt, U.: A Fabry–Perot interferometer-based camera for two-dimensional mapping of SO2 distributions, Atmos. Meas. Tech., 7, 3705–3715, https://doi.org/10.5194/amt-7-3705-2014, 2014. a, b, c, d, e
Kuhn, J., Platt, U., Bobrowski, N., and Wagner, T.: Towards imaging of atmospheric trace gases using Fabry–Pérot interferometer correlation spectroscopy in the UV and visible spectral range, Atmos. Meas. Tech., 12, 735–747, https://doi.org/10.5194/amt-12-735-2019, 2019. a, b, c, d
Lübcke, P., Bobrowski, N., Illing, S., Kern, C., Alvarez Nieves, J. M., Vogel, L., Zielcke, J., Delgado Granados, H., and Platt, U.: On the absolute calibration of SO2 cameras, Atmos. Meas. Tech., 6, 677–696, https://doi.org/10.5194/amt-6-677-2013, 2013. a, b
McGonigle, A. J. S., Hilton, D. R., Fischer, T. P., and Oppenheimer, C.: Plume velocity determination for volcanic SO2 flux measurements, Geophys. Res. Lett., 32, L11302, https://doi.org/10.1029/2005GL022470, 2005. a
Mori, T. and Burton, M.: The SO2 camera: A simple, fast and cheap method for ground-based imaging of SO2 in volcanic plumes, Geophys. Res. Lett., 33, L24804, https://doi.org/10.1029/2006GL027916, 2006. a, b, c, d
NOAA: NOAA, Global Monitoring Division, Solar Geometry Calculator, available at: https://www.esrl.noaa.gov/gmd/grad/antuv/SolarCalc.jsp, last access: April 2020. a, b
Perot, A. and Fabry, C.: On the Application of Interference Phenomena to the Solution of Various Problems of Spectroscopy and Metrology, Astrophys. J., 9, 87, https://doi.org/10.1086/140557, 1899. a
Platt, U. and Stutz, J.: Differential optical absorption spectroscopy, Springer Verlag, Berlin, Heidelberg, 2008. a, b
Platt, U., Lübcke, P., Kuhn, J., Bobrowski, N., Prata, F., Burton, M., and Kern, C.: Quantitative imaging of volcanic plumes – Results, needs, and future trends, J. Volcanol. Geoth. Res., 300, 7–21, https://doi.org/10.1016/j.jvolgeores.2014.10.006, 2015. a, b
Serdyuchenko, A., Gorshelev, V., Weber, M., Chehade, W., and Burrows, J. P.: High spectral resolution ozone absorption cross-sections – Part 2: Temperature dependence, Atmos. Meas. Tech., 7, 625–636, https://doi.org/10.5194/amt-7-625-2014, 2014. a
UWYO: University of Wyoming, Department of Atmopsheric Science, atmospheric soundings; Station data: Trapani (LICT), Italy 2019-07-22 12:00, available at: http://weather.uwyo.edu/upperair/sounding.html, last access: April 2020. a, b, c
Vargas-Rodríguez, E. and Rutt, H.: Design of CO, CO2 and CH4 gas sensors based on correlation spectroscopy using a Fabry–Perot interferometer, Sensor. Actuat. B-Chem., 137, 410–419, https://doi.org/10.1016/j.snb.2009.01.013, 2009. a
Veefkind, J., de Haan, J., Brinksma, E., Kroon, M., and Levelt, P.: Total ozone from the ozone monitoring instrument (OMI) using the DOAS technique, IEEE T. Geosci. Remote, 44, 1239–1244, https://doi.org/10.1109/tgrs.2006.871204, 2006. a, b, c
Veefkind, J., Aben, I., McMullan, K., Förster, H., de Vries, J., Otter, G., Claas, J., Eskes, H., de Haan, J., Kleipool, Q., van Weele, M., Hasekamp, O., Hoogeveen, R., Landgraf, J., Snel, R., Tol, P., Ingmann, P., Voors, R., Kruizinga, B., Vink, R., Visser, H., and Levelt, P.: TROPOMI on the ESA Sentinel-5 Precursor: A GMES mission for global observations of the atmospheric composition for climate, air quality and ozone layer applications, Remote Sens. Environ., 120, 70–83, https://doi.org/10.1016/j.rse.2011.09.027, 2012. a
Wilson, E. L., Georgieva, E. M., and Heaps, W. S.: Development of a Fabry–Perot interferometer for ultra-precise measurements of column CO2, Meas. Sci. Technol., 18, 1495–1502, https://doi.org/10.1088/0957-0233/18/5/040, 2007. a