Ground-Based Remote Sensing and Imaging of Volcanic Gases and Quantitative Determination of Multi-Species Emission Fluxes

The physical and chemical structure and the spatial evolution of volcanic plumes are of great interest since they influence the Earth’s atmospheric composition and the climate. Equally important is the monitoring of the abundance and emission patterns of volcanic gases, which gives insight into processes in the Earth’s interior that are difficult to access otherwise. Here, we review spectroscopic approaches (from ultra-violet to thermal infra-red) to determine multi-species emissions and to quantify gas fluxes. Particular attention is given to the emerging field of plume imaging and quantitative image interpretation. Here UV SO2 cameras paved the way but several other promising techniques are under study and development. We also give a brief summary of a series of initial applications of fast imaging techniques for volcanological research.


Introduction
The physical and chemical structure and the spatial evolution of volcanic plumes is of great interest for a number of reasons beyond scientific curiosity: (A) Volcanic gas emissions influence the atmosphere and therefore also the climate and other Earth system parameters in a number of ways and on different temporal and spatial scales (e.g., [1][2][3]).Investigations of plume chemistry and plume dispersal will help constrain these influences (see e.g., [4,5]).(B) The composition and emission rate of volcanic gases are linked to processes occurring in the Earth's interior, therefore measuring volcanic gases provides insights into these otherwise largely inaccessible processes.For instance, already Noguchi and Kamiya [6] showed that eruptions of Mt.Asama (Japan) could be forecast with some degree of accuracy by measuring the variable partitioning of sulphur species, chlorine species, and carbon dioxide in the emissions from the active crater.Malinconico [7] showed first that also the amount of gas, in particular the amount of emitted SO 2 varies when the volcanic activity changes.Several authors (e.g., [8,9]) used the measured SO 2 emission also to calculate the amount of magma involved in the simultaneously observed volcanic activity.
These remote sensing techniques have a number of decisive advantages over in-situ observations: The largest being that under most conditions the total amount of gas in the plume can be determined, rather than the gas concentration at the plume edge (usually probed by ground-based in-situ instruments).Since measurements can be made from distances of typically a few kilometres remote sensing is also much safer than in-situ sampling.Moreover, the technology allows easy automation and thus continuous measurements in real time are readily possible, even during periods of explosive activity.In fact continuous measurements of SO 2 and BrO fluxes are realised by the global Network of Observation of Volcanic and Atmospheric Change (NOVAC) [10][11][12].
Recent advances in technology, in particular for spectroscopic techniques, allow remote analysis of many species in volcanic plumes.Specifically, the two-dimensional or even three-dimensional distribution of gases within volcanic plumes and their temporal evolution can now be determined in real time.In other words, recent technology allows "imaging" of trace gas distributions and their motion in plumes.Compared to earlier techniques, which were only capable of measuring the total column density of a particular gas along a single line of sight within the plume (e.g., [13][14][15]), new imaging techniques (see e.g., [5,16]) give much better insight into transport and mixing processes, as well as into chemical transformations within plumes.Moreover, much more accurate, quantitative determinations of trace gas fluxes are now possible in most cases.If the imaging is performed at adequate time resolution, the plume speed can be derived directly from image series, either by correlation techniques (e.g., [13,[17][18][19]) or by more advanced image processing techniques (e.g., [20][21][22]).
The most abundant volcanic gases, i.e., H 2 O, CO 2 , SO 2 , HCl, HF, H 2 S and many other trace species like BrO, ClO, OClO, NO 2 CO, COS, SiF 4 , can be measured either by UV/vis spectroscopy (i.e., H 2 O, SO 2 , BrO, ClO, OClO, NO 2 are readily measured by DOAS) or IR spectroscopy (e.g., H 2 O, CO 2 , SO 2 , HCl, HF, H 2 S, CO, COS, SiF 4 ).Note, however, that the two species, which are usually most abundant in volcanic plumes (i.e., H 2 O and CO 2 ) also occur in high concentrations in the atmosphere surrounding the plume.Therefore it is frequently difficult to determine these gases by passive spectroscopy (see Section 2.3.3),since the contrast between the background atmospheric signal and background atmosphere plus plume is very small (sometimes in the 10 −3 range or below).
In the following we give an overview of the principles of optical remote sensing and of plume imaging techniques.These techniques can be applied from the ground, from (manned or unmanned) aircraft (e.g., [23]), or from satellite platforms.Although satellite imaging constitutes an interesting and important branch of these techniques (e.g., [24][25][26]) we restrict ourselves to the first type of application.
At present ground-based remote sensing techniques are in various stages of development, ranging from methods that are already routinely being applied to study volcanic degassing, to those which are so new that only first field test have been done.In addition we discuss applications to flux measurements.

Remote Sensing of Volcanic Gases
Remote sensing of volcanic plume composition relies on measuring the absorption, scattering or emission of radiation by the gases (and particulates) in the plume.At this point it is convenient to make a distinction between active (using an artificial light source) and passive (using natural light sources, e.g., the sun or thermal emission) remote sensing.In general, the radiative transfer equation (see e.g., [27]) underlying all these techniques would need to be solved for a remote sensing observer collecting radiation from a certain direction.The remote sensing observer could be located on a satellite in space, on an aircraft or another airborne platform, or on ground.The complexity of the remote sensing problem strongly depends on the chosen vantage point and viewing geometry, the spectral ranges covered, and the abundances of particulates and the target gases in the plume and in the background atmosphere.Figure 1 shows some basic ground based plume sensing schemes.In the following we give an overview of the principles of optical remote sensing and of plume imaging techniques.We limit ourselves to showcasing ground-based techniques that are regularly used for volcanic plume remote sensing.We recapitulate the basic aspects of absorption (Section 2.1) and thermal emission (Section 2.2) spectroscopic techniques and highlight some applications (Section 2.3).Imaging techniques are discussed in Section 3, below.

Absorption Spectroscopy
The classic absorption experiment, illustrated in Figure 2, puts a light source at a distance L from the observer such that the absorbing medium is in-between.Neglecting light-scattering and thermal emission by the medium itself, Beer-Lambert-Bouguer's law describes the change of the radiance spectrum I(λ, x) due to absorption along the light path x, with wavelength λ and absorption coefficient k a (λ, x).For a gaseous medium, the latter is given by where n i (x) is the number density of absorbing gas species i (of which there might be multiple), and σ i (λ, p(x), T(x), . . . is the respective absorption cross section that generally depends on the path via ambient pressure p, temperature T, or other variables like water vapour that causes spectroscopic effects such as foreign line broadening.
Assuming that the light source emits a background spectrum I(λ, 0) = I 0 (λ) at location x = 0, integration of Equation (1) yields the radiance spectrum at the location x = L of the remote sensing observer, where we define absorption optical density τ(λ, If the absorption cross sections can be assumed independent of ambient conditions along the path (p(x), T(x), ... = const.),the absorption optical density is approximately given by where S i denotes the slant column density of trace gas species i defined as the concentration integrated along the light path.For atmospheric remote sensing, this assumption typically holds for absorption cross sections in the UV and visible spectral range and thus, the Differential Optical Absorption Spectroscopy (DOAS) technique described in Section 2.3.1 adopts Equation ( 5).In the infra-red, the assumption is not generally applicable and absorption spectroscopic techniques need to start out from Equation (4) (Section 2.3.2),unless the absorbing gases can be safely assumed to only exist in a single plume layer with homogeneous pressure, temperature etc.Assuming that Equation ( 5) holds and that there is only a single absorbing gas in the volcanic plume, the respective gas slant column density S can, in principle, be determined by ratioing a background spectrum I 0 (λ) by the plume spectrum I(λ) according to, from which the plume average absorber concentration n is straightforward to calculate through if the species only exists inside the plume (and not in the background atmosphere) and if the geometric path length L p through the plume can be estimated e.g., using the geometric plume extent.In practice, various complications render this approach too simplistic and real-world solutions need to be found (Section 2.3).

Thermal Emission Spectroscopy
Thermal emission spectroscopy makes use of the infra-red emission of the atmospheric gases themselves i.e., it collects radiation emitted along the lines-of-sight through the atmosphere in the infra-red spectral range without targeting a specific light source as illustrated by Figure 3.Typically, Schwarzschild's approximation to the radiative transfer equation is used to describe the emitted radiance spectra, where k a (λ, x) is the absorption coefficient defined in Equation ( 2) and B(λ, T(x)) is Planck's thermal emission spectrum (c = speed of light, h = Planck constant, k = Boltzmann constant, T = temperature).Solving Schwarzschild's equation for an observer at location x = L pointing toward a cold background at location = 0 (I(λ, 0) = 0), with the volcanic plume in-between, yields the observed radiance with τ(λ, x, L) the absorption optical density defined in Equation (5).Equation (10) essentially tells that thermal radiation accumulates along the path toward the observer taking into account that radiation emitted at location x gets reabsorbed along the rest of the way to the observer Dividing the light path into N homogeneous layers k at constant temperature T k with layer boundaries x k and x k+1 , Equation (10) becomes which after calculating the integral reduces to Introducing the transmittance t k = exp(−λ(λ, x k , x k+1 ) of layer k for short-hand notation, Equation (11) reads A common view on Equation ( 13) for volcanic applications is the one shown in Figure 3.The emission along the line-of-sight is assumed to consist of three contributions, the atmospheric background radiance I behind (λ) entering the plume from behind, thermal emission by the plume layer (2), and thermal emission by the layer (3) in front of the plume which hosts the observer.Then, Equation (13) simplifies to Equation ( 14) carries the target gas concentrations in the plume transmittance t 2 , and thus can be used to setup plume remote sensing experiments (Section 2.3.3).
For further conceptual insight, we assume that a spectral range can be found where the target gas is the only absorber and that the target gas only exists in the volcanic plume.Then, all layer transmittances except t 2 become unity and the background radiance I behind (λ) vanishes, reducing Equation (14) to Following the discussion of Equation ( 5), for a homogeneous plume layer (p(x), T(x), ... = const), the optical density can be approximated by τ(λ) = σ(λ) × S.Then, Equation (15) readily yields the slant column density S of the target gas via from which the plume average concentration n calculates through Equation ( 7).Thus, under the above approximations, measuring a thermal emission spectrum I(λ, L) and measuring (or knowing from external sources) the plume layer temperature T 2 would allow for estimating the gas concentrations.
In practice, infra-red spectral windows are typically packed with overlapping absorption of various background and plume gases and thus, Equation ( 16) is too simplistic.
Figure 3. Schematic illustration of thermal emission sounding for a simplified 3-layer atmosphere consisting of a layer behind the volcanic plume, the plume-layer itself, and the layer between the plume and the observer.

Applications to Remote Sensing of Volcanic Plumes
There is a great variety of ways implement the above theoretical considerations for the purpose of volcanic research and surveillance.The currently available techniques can be broadly grouped into the following categories (see also Figure 4):

•
Active or passive spectroscopy (i.e., will there be an artificial light source or a natural one)  Optical absorption spectroscopy (see e.g., [27]) is a technique following the principles outlined in Section 2.1.Absorption cross sections can be assumed to be approximately independent of pressure and temperature in the UV/visible spectral range, and thus, Equation (3) can be combined with Equation (5).However, there is the challenge how to disentangle the minute target gas absorption from overlapping interfering absorption, from spectrally broad extinction of radiation due to light-scattering, and from solar Fraunhofer lines, if the sun is the light source (which is a particularly popular case).Facing these challenges the COSPEC (COrrelation SPECtrometer, e.g., [28]) techniques pioneered UV/vis spectroscopy for volcanic applications from the 1970s.As another approach to address the above challenge, the DOAS (Differential Optical Absorption Spectroscopy) techniques (e.g., [27]) measure absorption spectra in a contiguous range of wavelengths for which the molecular absorption cross sections show several differential structures that are unique fingerprints for the absorbing gases.This implies that absorption cross sections must show strong variations with wavelength to make a gas measurable by DOAS.Today, modern instruments are mostly based on the DOAS technique (e.g., [27]), which our discussions focus on in the following.DOAS can be performed with miniaturized spectrometers and, it is well suited to detect even smallest gas amounts.
Building on the overall idea of using differential absorption structures, DOAS techniques are applicable to direct light-source (sun, moon, or artificial sources, see below) arrangements as well as to applications using scattered (solar) light.In the former case, e.g., when directly pointing into sufficiently bright light sources such as the sun, the path x is the well-defined geometric path between the source and the observer.In the latter case, e.g., when pointing into the sky away from the sun, the path x is the average geometric path that photons travel between the light source and the observer along various individual trajectories such as suggested by Figure 5.This interpretation makes the path x, its length L and in consequence the slant column density S i in Equation (3) depend on wavelength since the average scattering light path depends on atmospheric light-scattering properties such as the particle scattering cross sections and phase functions as well as on the molecular absorption itself, which are all wavelength dependent.However, the narrow spectral ranges and the frequently small absorption optical densities (τ(λ) = 10 −4 ... 10 −1 ) encountered in DOAS applications, typically (but not always, see below) allow for dropping wavelength dependencies in x, L, and S i .
Although solar radiation has considerable intensity in the near IR (up to 2 to 3 µm) typically scattered sunlight absorption spectroscopy is only used in the UV/vis spectral region since it usually relies on Rayleigh scattering, where the scattering efficiency falls off strongly with wavelength (approximately proportional to λ −4 ).There is usually also a contribution of Mie (aerosol) scattering, however the intensity of this component also falls off with increasing wavelength (typically proportional to λ −1.3 ).A great advantage of using scattered solar light is its universal applicability with observation from only one point, which greatly simplifies the logistics.Dispersive or non-dispersive detectors can be used, up to now the former are employed for one dimensional work and the latter for imaging applications.Thus, conceptually, DOAS is fit to use Equations ( 6) and ( 7) to infer gas concentrations in volcanic plumes.As illustrated by Figure 5, the plume spectra I(λ) in Equation ( 6) are the measurements for which the lines-of-sight cross the volcanic plume.The background spectra I 0 (λ) are obtained by adjusting the observing geometry such that lines-of-sight avoid the volcanic plume or, if such an adjustment is not possible, by resorting to external sources such as solar spectral atlases [29].In practice, this simple approach faces various complications which need to be overcome.
First and foremost, the absorption optical densities of the target species are typically small (τ(λ) = 10 −4 10 −1 ) in the UV/visible spectral range and the target absorption usually interferes with strong Fraunhofer lines in the solar spectrum and with absorption by other gases.As explained above, DOAS relies on fingerprinting the target absorption by measuring contiguous spectral ranges covering a sufficient number of differential absorption structures of the target gas.For accessing the minute target absorption with high accuracy, however, it is further crucial to exclude systematic errors originating from uncertainties in the Fraunhofer spectrum or in the absorption cross sections causing spurious effects on the subjacent gas signature.Therefore, DOAS prefers using background spectra (and, if possible, absorption cross sections) that are measured by the actual remote sensing instrument.Then, the aforementioned error sources efficiently cancel, making minute absorption signals accessible.
Strictly, however, Equation ( 6), assumes the spectra I 0 (λ) and I(λ) to be available at an 'infinite' spectral resolution i.e., at a spectral resolution much better than the width of the atmospheric spectral features.This is generally not the case when using remote sensing measurements conducted by a field spectrometer.The instrument spectral response function (ISRF) smoothes the atmospheric spectra according to where the symbol ⊗ indicates the convolution operation and F(λ) is the ISRF that characterizes the observing instrument.The ISRF can be measured in the laboratory or calculated based on the instrument's optical and detector properties.Laser based approaches typically have spectrally narrow ISRF and thus, Equation ( 17) can indeed be dropped.For DOAS applications, replacing I, I 0 , and σ in Equation ( 6) by the measured (or ISRF convolved) quantities I * , I * 0 and σ * holds approximately in optically thin conditions.The approximation is applicable for most UV/visible absorbing gases and work-around's can be found for moderately thick absorption optical densities such that DOAS is widely applicable to trace gas remote sensing in this spectral range.
In practice, DOAS is applied to the retrieval of UV/visible absorbing gases such as H 2 O, SO 2 , BrO, OClO, NO 2 and O 3 using scattered sunlight as a light source or by directly pointing toward the sun or artificial light sources.In the latter case a broad-band (relative to the absorption band or line) artificial light source, typically in a search light arrangement, supplies a light beam with intensity I 0 (λ), which traverses the plume.At the other side of the plume the light is collected and the column density is calculated form the absorption spectrum.Frequently a variant of this set-up is used, where light source and detector are at one end of the light path, while at the other end there is only a reflector (typically a retro-reflector).This arrangement has two advantages: First, the light beam crosses the plume twice, thus doubling the sensitivity and second (and usually more important), all parts of the instrument, which require power are at one end of the light path, while only a passive reflector is at far end.Prerequisite for this approach is the availability of suitable sites on either side of the plume (e.g., two points at the crater rim) for mounting light source and detector or reflector (see e.g., [30,31]).In addition, the technique largely avoids the problems of the large background column density of e.g., H 2 O.In principle dispersive or non-dispersive detectors could be used, however, up to now in all applications only the former type is employed.
If the gas under investigation only exists inside the plume and not in the background atmosphere and if scattering effects are not too severe, estimating n from S can be as simple as Equation ( 7) i.e., dividing the slant column by the length of the line-of-sight transect through the plume.In practice, for measurements using scattered sunlight estimating the intra-plume concentrations from the slant columns S might be complicated due to several effects: (1) the intra-plume average light path has contributions from multiple scattering inside the plume due to large abundances of plume condensate or aerosols, (2) the overall average light path has contributions from photons that do not transect the plume at all and thus, they do not see any target gas absorption (commonly termed light dilution), (3) the target gas has non-negligible abundances in the background atmosphere and thus, background removal is required, or (4) the absorption optical density is not thin and the average scattering light path x cannot be assumed independent of wavelength and the convolution operation does not commute with the exponential in Beer-Lambert-Bouguer's law.
Following the early introduction of COSPEC (COrrelation SPECtrometer, e.g., [28]), today's instruments are based on the DOAS technique (e.g., [27]), which can be performed with miniaturized spectrometers.There are also variants with simplified spectral evaluation e.g., Flyspec (e.g., [32]).In principle reflection of sunlight on topographic objects or by clouds could be used, but this has not been attempted yet.

Absorption Spectroscopy in the Infra-Red
In the infra-red (from about 1 to 20 µm wavelength), the absorption cross sections usually depend strongly on pressure and temperature and the absorption optical densities are frequently large, this is particularly true for the thermal IR.Thus, the DOAS approach is not readily applicable.So, generally, Equation (3) in combination with Equations ( 4) and ( 17) must be used to forward model the measured absorption spectra based on an initial assumption of gas concentration profiles, plume geometry, and ambient meteorological conditions.An inverse estimation technique iteratively adjusts the target gas concentrations to statistically best match the measurements until some convergence criterion is met.
Absorption spectroscopy in the infra-red mostly relies on direct-light approaches i.e., techniques that have a detector which only accepts radiation from a narrow solid angle centred directly on the light source-be it the sun, a hot telluric source (e.g., lava) or an artificial light source.Thus, Figure 2 is directly applicable.The incoming radiance I 0 (λ) in Equation ( 3) would be the solar Fraunhofer spectrum for direct sun observations, it would be Planck's thermal emission spectrum B(λ, T), Equation ( 9), for a hot opaque source with temperature T such as lava, or it would be a silicon carbide lamp (Globar), incandescent lamp, LED (light emitting diode), or laser emission spectrum for artificial light source applications.
As described in Section 2.3.1 direct-light techniques have the great advantage that the light path is well-defined.As long as the light source is sufficiently bright to exclude contributions from light scattering and thermal emission along the atmospheric path, the light path computes easily from the positions of the source and the observer and, if necessary, from the atmosphere's refractive properties.The major disadvantage of direct-light techniques (as in the case of active UV/vis spectroscopy) is that one has little flexibility in choosing the observing geometry.The line-of-sight is strictly defined by the position of the light source and the receiving system, which have to be arranged such that the volcanic plume is in-between.
IR spectroscopy is implemented in four popular variants: (a) With broad-band, thermal light sources, e.g., globars in combination with dispersive detection systems, typically Fourier-Transform interferometers of the Michelson type ( [33][34][35]).At volcanoes hot lava can be used as source of radiation (e.g., [36], see below), which-according to our definition-would be classified as passive absorption spectroscopy.(b) Tunable Diode laser spectroscopy (TDLS), which (e.g., [37][38][39]) is a variety of dispersive spectroscopy where a narrow-band emitting light source (i.e., a semiconductor laser) is rapidly wavelength modulated in order to sweep across an absorption line of the gas to be measured.The original, rather unreliable lead salt diodes are now replaced by much more stable (though still expensive) quantum cascade laser diodes [40,41].Since the light source is wavelength modulated there is no need for dispersive detection.
All variants of the technique share the prerequisites (suitable sites) and advantages (low background) with the active UV/vis spectroscopy.Also the options for the arrangement of light paths (single path, retro reflector) are similar, in addition 'topographic targets' i.e., back-reflection of the transmitted radiation at terrain surfaces (see below) are in use.
Future developments for active IR spectroscopy could use Dual frequency comb spectroscopy (e.g., [42]), which-although it works differently-has some similarities to a combination of TDLS plus FT spectrometer.
(c) Passive IR spectroscopy using the sun (or the moon) as a direct light source is commonly referred to as the solar (or lunar) occultation technique which has been used for remote sensing of volcanic gases such as SO 2 , HF, HCl, and SiF 4 , e.g., [43][44][45].The technique simplifies substantially if the target gas only exists (at relevant amounts) in the volcanic plume and not in the background atmosphere which is, however, not the case for the major volcanic plume constituents CO 2 and H 2 O. Due to the rapid downwind dispersion of the volcanic plume, the volcanic enhancements (on top of the large background) become small and thus, increasingly difficult to measure the farther downwind the plume is sampled.Only recently, Butz et al. [46], demonstrated safe-distance remote sensing of volcanic CO 2 in Mt.Etna's plume during passive degassing conditions.They operated a sun-viewing, portable Fourier Transform Spectrometer (FTIR) on a truck in stop-and-go patterns underneath Mt.Etna's plume such that the lines-of-sight to the sun sampled the plume in 5-10 km distance from the crater.Co-measuring O 2 columns helped calibrating spurious variations in the targeted CO 2 columns which were merely due changes in observer position.Sequentially measuring intra-plume and extra-plume spectra and using co-measured HF, HCl and SO 2 as intra-plume tracers helped removing the atmospheric CO 2 background.These current generation instruments were able to discriminate the volcanic CO 2 signal out of a 300-1000 times larger atmospheric background path.(d) Hot volcanic material such as lava or volcanic rocks have been used heavily in open-path spectroscopic techniques, e.g., [33,35,36,43,[47][48][49][50][51]  Generally, the technique requires that the hot material or the lamp locates behind the plume and that it can be sighted by the observer.For many volcanoes, this requirement implies deploying instrumentation in the proximity of the crater or at the crater rim which spoils the general remote sensing advantage of avoiding hazards and hostile environments to operators and instrumentation.Recently, laser-based techniques for CO 2 have been developed using topographic reflection targets, e.g., [52,53], promising greater deployment flexibility but still requiring plume sampling close to the source where CO 2 enhancements are large compared to the atmospheric background.

Thermal Emission Spectroscopic Techniques
Typically, Equation ( 14) forms the basis for volcanic plume remote sensing by thermal emission spectroscopy, e.g., [54,55].Using a forward modelling approach similar to the one used for direct-light absorption spectroscopic techniques, one starts out with an initial guess for the ambient and plume temperatures as well as for the transmittances.The transmittances contain the layer-wise optical densities and thus, the information on the layer-wise gas composition.Iteratively, the targeted plume gas abundances are adjusted to yield a best match of measured and forward modelled spectra.
Thermal emission spectroscopy allows for the detection of other species than UV/vis absorption spectroscopy (see Section 2.3.1) and for observation independent of sunlight, i.e., also at night time.
Disadvantages include complicated set-up (FT-IR), usually cryogenic cooling of the detector is required.However recently also non-dispersive detection with detectors at room temperature was demonstrated [56,57].The classic approach is the differential one sketched in Figure 6, using the difference between two measurements, one pointing into the plume and one pointing next to the plume into the background sky, e.g., [54,58].For these two measurements, the recorded measurements approximately differ by temperature T 2 and transmittance t 2 in Equation ( 14) either referring to intra-plume or background conditions.Thus, the difference spectra isolate the information on the targeted plume composition which is straightforward to retrieve as outlined in Goff et al. [54], Stremme et al. [55], and Krueger et al. ( 2013) [59], used a thermal emission FTIR to sequentially scan the SO 2 and SiF 4 plume (and the adjacent background sky) of Popocatepetl volcano from 12 km distance providing a series of two-dimensional images of the volcanic gas columns.Then, they used the series of images to simultaneously estimate the plume average wind speed and the constituent outflux from the crater.Current developments toward imaging FTIRs (e.g., [57,60]) are promising for use in volcano plume monitoring since they allow for rapidly imaging the thermal emission spectrum emerging from two-dimensional scenes without the need for scanning a telescope.

Range Resolved (LIDAR) Techniques
A further, well known approach to active remote sensing trace gases is the LIght Detection and Ranging (LIDAR) technique, e.g., [61,62].The usual set-up is the mono-static LIDAR where transmitter and receiving system are at the same location.The (laser) light source emits short (<few 100 ns) pulses into the atmosphere.By analyzing the time profile of the intensity back-scattered from the atmosphere the spatial distribution of scattering and extinction along the direction of the emitted (and received) radiation can be deduced.With single-wavelength lasers aerosol backscatter and extinction are measured.
In order to determine trace gas distributions a more complicated set-up is required.These devices have become known as Differential Absorption Lidar (DIAL) using at least two wavelengths, one "on absorption" of the trace gas and one "off absorption" [13,14,63,64].The approach is similar to the SO 2 -camera (see Section 3.2.1).While DIAL systems have been bulky, relatively complicated and power consuming recent developments including micro-pulse diode lasers (e.g., [65]) promise smaller, lower power, and more portable systems.
Alternatively to the customary mono-static approach a bi-static LIDAR scheme has been developed (e.g., [4,66,67]).Here, radiation source and radiation detector are separated by a small distance (on the order of 10 m for a few km range) and the range resolution is accomplished by changing the elevation angle of either the radiation source or of the radiation detector.Advantages are the possibility to use continuous light sources (like LEDs) and simpler detection.
In principle two-dimensional or even three-dimensional push-broom imaging of volcanic plumes is possible with LIDAR instruments and to some extent has been demonstrated at volcanic plumes (e.g., [13,14,63]).

Imaging of Plumes
Many of the developments for volcanic plume remote sensing have started with collecting a series of individual plume (column density) measurements in order to sample one-dimensional plume cross sections.Modern technology allows to 'image' plumes, i.e., to provide two-dimensional arrays of optical densities (or 'images') of a scene containing volcanic gases.This approach offers many advantages over the simpler one-dimensional plume scanning:

•
Redundant measurements can be made by e.g., making trace gas flux determinations at several planes along the plume propagation direction allowing e.g., internal consistency checks (see e.g., [22]) Redundant flux measurements can be used to determine the exact plume propagation direction (see [22]) Last not least: the human visual system has powerful analysis capacities which can be used once images are available ('seeing is believing') For these (and more) reasons a series of imaging systems for volcanic plumes have been developed and applied (see e.g., [5,16]).In this context the requirements for plume imaging instruments include:

•
Capability to specifically detect the desired gas (or parameter)

•
To provide sufficient sensitivity (e.g., for SO 2 measurements a detection limit for SO 2 -column densities of the order of 10 17 molecules/cm 2 or ca.40 ppmm is required for volcanic emissions observations)

•
To provide sufficient spatial resolution to allow discrimination of the relevant features within the plume(s) • Sufficient time resolution (typically of the order of seconds) of the measurement is further required to be able to resolve the motion of the plume and variations in the volcanic source strength Some of these requirements are in conflict with each other.For example enhanced specificity may be achieved by recording of the intensity at several wavelengths, since this typically can only be done sequentially the time resolution will be reduced.Likewise, lower detection limits (i.e., higher sensitivity) are achieved by collecting a larger number of photons per pixel, which comes at the expense of time resolution (or, alternatively, spatial resolution).In particular, the number of gases which can be detected by imaging instruments is limited in comparison to one dimensional instruments.
It should be noted that there is a connection between spatial and temporal resolution in that a given spatial resolution requires sufficient temporal resolution to be useful.The motion of the plume ∆x = v•∆t within the acquisition time ∆t of an image (v = wind speed) should be at least of the order of the size of one picture element in the plane of the plume.In other words, high spatial resolution makes sense only if the technique also offers sufficient temporal resolution (for a given wind speed).Note that the opposite is not necessarily true: A temporal resolution exceeding the above limit would still result in useful (though partly oversampled) data.
As described in Section 2 several remote sensing techniques for the imaging of plumes have been developed during recent years.As described in Sections 2.3.2 and 2.3.3 the various techniques can be grouped into active and passive as well as in dispersive and non-dispersive approaches.For dispersive approaches several scanning schemes are possible as detailed in Section 3.1, below.
The characteristics of dispersive techniques is the determination of a spectrum (i.e., a series of intensities as function of wavelength) for each pixel of the image.The trace gas column density for each pixel is then derived from the measured spectrum using the algorithms described in Section 2, above.Non-dispersive approaches, on the other hand, derive the column density from the ratio of just a few (typically two) intensity measurements integrated over suitable (narrow) wavelength intervals.Frequently, a reference intensity ratio is recorded with the instrument pointed away from the plume.The trace gas column densities are then derived from the ratio of the sample and reference intensity ratios.
Dispersive approaches are typically more complex and since usually scanning is required slower than non-dispersive techniques.However, they have the advantage of being less vulnerable to spectroscopic interferences and to be able to also measure minor constituents like BrO due to higher spectral resolution enabling higher specificity for gas detection.Furthermore, analysis using traceable, well quantified absorption cross-sections increases the rigour of the technique.Many gas absorption features are naturally narrow (of the order of 1-2 nm), particularly for small molecules, and this limits the effectiveness (i.e., the sensitivity as well as the capability to specifically detect a particular species) of non-dispersive techniques.
Note that any combination of dispersive/non-dispersive techniques is in principle possible and several combinations have been realized, such as the combined use of a dispersive UV spectrometer alongside a non-dispersive UV imaging system (see e.g., [71]).In fact, an important aim of present research is the development of techniques combining the speed of non-dispersive techniques with the rigor of dispersive techniques.Overall, despite impressive accomplishments, techniques for quantitative imaging of volcanic plumes are still in an early stage of development and much progress can be expected in the near future.

Categories of Plume Imaging
Two dimensional images of trace gas column amounts are derived from three-dimensional data: two spatial dimensions and one wavelength (or interferogram) dimension are recorded.At present three different imaging techniques are known (see Figure 7 and Table 1): (1) Pixel at a time scanning ("whiskbroom" imaging, see Figure 7A): In this approach all pixels of an image scene are scanned sequentially according to a particular scheme (e.g., line by line as in early TV cameras), for each pixel a spectrum is determined.In this approach all (of e.g., 10 5 pixels) have to be scanned individually, such that it is potentially a rather slow approach.
Michelson interferometers have been combined with whisk-broom scanners to obtain 2-D images (e.g., [34,55,72]).( 2) Column (or row) at a time scanning ("push-broom" imaging, see Figure 7B): Here all pixels of an image column are scanned simultaneously while the image columns are scanned sequentially.
Because each column of the image is recorded at once, only of the order of several hundred columns have to be measured sequentially.Since only one scan-dimension is required the scanning mechanisms can become simpler though it is not necessarily faster.In fact, as will be explained below the amount of radiation collected by the entrance optics has to be split between all pixels of a column, thus, when the time to acquire an image is determined by the available number of photons the technique will not generally be faster than whiskbroom imaging.Push broom scanners have been realized with DOAS instruments as e.g., described by Lohberger et al. [73], Bobrowski et al. [74], Louban et al. [75] and Lee et al. [76], see Section 3.3, below.Michelson interferometers have also been combined with push-broom scanners (i.e., moving platform) to obtain 2-D images (e.g., [77]).A recent development based on a special type of interferometer (Sagnac interferometer) is the Thermal Hyperspectral Imager, which produces a spatial interferogram across the field of view, which is scanned across the image.After Fourier transformation a high resolution spectrum for each image pixel is obtained [57] (Gabrieli et al. 2016), which can be analysed for spectral signatures of volcanic gases.Gabrieli et al. [57,78] used such a device to produce images of the SO 2 distribution derived from spectra around 8.6 µm at a spectral resolution of about 0.25 µm.Explorative measurements were made at Kilauea Halema'uma'u crater (Hawaii) with a scan duration of 1 s.(3) Frame at a time scanning ("full frame" imaging, see Figure 7C): Here the entire frame is recorded at once (or in a sequence of steps in time).While in principle there could be large arrays of spectrometers determining the spectrum of each pixel (and in the future arrays of integrated micro-spectrometers could become viable), in practice the spectral information is usually determined by collecting sequential images with different wavelength selective elements (e.g., suitable filters, see Figure 8) in front of the camera sensors.
Current research also aims to employ an array of detectors for two-dimensional Michelson interferometers, where each pixel can effectively be thought of having its own interferometer.An example is GLORIA (e.g., [60]), which uses a 256 × 256 element Mercury Cadmium Telluride focal plane array (FPA) cooled to 60 K.The spectral coverage of the instrument is from 7.1 µm to 12.8 µm.

Gas correlation
Presently not used A A For volcanic plume imaging.

Non-Dispersive Plume Imaging
Non-dispersive imaging usually uses one or several two-dimensional image sensor(s), being sensitive in the desired spectral range in combination with a device offering (limited) spectral selectivity.This compromise favours spatial resolution and imaging speed over (usually) sensitivity and trace gas selectivity.

Quantifying Column Densities Using Classic (UV) SO 2 Cameras
The theory developed in Section 2.1 is best illustrated by briefly inspecting a "classic" SO 2 Camera.Such an instrument consists of an UV-sensitive camera viewing e.g., the area around an actively degassing volcanic vent through a UV band-pass filter that selectively transmits radiation at a wavelength λ A which is strongly absorbed by SO 2 (about λ A ≈ 310 nm).Thus the sections of the image receiving radiation which has passed through the volcanic plume (assumed to contain high amounts of SO 2 ) receive a lower radiation intensity I(λ A ) than the background sky (i.e., sections of the image where the radiation did not pass the plume), which ideally has the intensity I 0 (λ A ).The attenuation depends on the SO 2 column density and is quantitatively described by Equations ( 3) to ( 5) in Section 2.1.While this simple approach works in principle, in reality a higher sensitivity is reached by using two filters (usually referred to as "Channel A" and "Channel B") transmitting radiation at different wavelengths as described e.g., by Kern et al. [82,83], Platt et al. [5], Smekens et al. [84], McGonigle et al. [16].
In order to improve the absolute accuracy of SO 2 -cameras it has become customary to add a 1-D (DOAS) spectrometer pointing at a point within the field of view of the camera (typically the centre).This approach has been pioneered by Lübcke et al. [71].

Non-dispersive IR Imaging of Plumes
Plume imaging by thermal emission is a passive technique relying on the thermal emission from the plume (at ambient temperature) itself, see Section 2.3.3.Unlike the case for the SO 2 camera described above, the equivalent infra-red system has received far less attention, although several approaches are being studied: (1) IR-cameras with two or more filters similar to the SO 2 camera principle (e.g., Prata and Bernardo 2009) have been developed for SO 2 retrieval and ash detection ( [56,85]).A four filter IR camera was actually used by Lopez et al. [86] to simultaneously determine plume the SO 2 SCD, temperature, and ash content of the plume of Stromboli (Italy), Karymsky (Russia), and Láscar (Chile) volcanoes.(2) Gas-correlation spectroscopy has also been applied in the IR for real time imaging of ammonia and ethylene (e.g., [87,88]).
In the following we give a brief description of an IR-camera system already in use for volcanic imaging.The absorption spectrum in the thermal infra-red region (7-13 µm) exhibits many features attributable to gases present in volcanic emissions, including SO 2 , CO 2 , H 2 O, HCl among others.Of these gases SO 2 is the easiest to detect because of its generally low background concentration and broad absorption peaks at 7.3 µm and 8.6 µm.Atmospheric volcanic ash also absorbs and scatters infra-red radiation in this region of the electromagnetic spectrum and its characteristic features have been exploited by using infra-red satellite measurements for many years (e.g., [89]).There are some important differences when sensing gas and ash clouds from the ground compared to measurements from space.The most important of these is the thermal contrast between target (gas or ash cloud) and background, for example the warm surface below the plume in the case of satellite sensing, and the cold clear sky behind the plume in the case of ground-based sensing.A typical infra-red camera system consists of a cooled or uncooled detector array, focusing optics, infra-red interference filters and data recording electronics.For instance, Prata and Bernardo [56] describe the development of an uncooled multispectral imaging camera system based on a 320 by 240 bolometric detector and F/1 Germanium (Ge) optics.The system incorporates four narrowband interference filters with selectable central wavelengths and ~1 µm bandwidths.The usual configuration is to use 8.6 µm for SO 2 ; 10, 11 and or 12 µm for ash detection and quantification; and a broadband channel (8-12 µm) for plume and background scene temperature measurements.

Dispersive Imaging
Dispersive imaging instruments have been realized for the UV/visible as well as for the IR spectral ranges.
The Imaging DOAS (I-DOAS) technique brings advantages from different techniques together: a good spatial resolution, like the SO 2 camera, and high spectral information similar to the one used in conventional DOAS [73].I-DOAS measurements result in a three dimensional data set, including two-dimensions of spatial information (i.e., each picture element of an image (pixel) corresponds to a defined solid angle of space), and a third dimension containing highly resolved spectral information within each pixel.For this reason imaging spectroscopic instruments require both imaging and dispersive optical components.Today I-DOAS has been applied using various platforms: satellites (e.g., OMI [90]), airplane (e.g., [23,91]) and ground based (e.g., [73][74][75][76]).An imaging satellite instrument operating in the thermal IR was proposed by Wright et al. [77].
The major advantage of I-DOAS is the capability of the technique to measure several trace gases simultaneously, enabling plume dispersal and chemical transformations within the plume to be studied.Also, I-DOAS is less dependent on meteorological conditions compared to SO 2 cameras.Furthermore, no calibration is needed and radiative transfer corrections can be calculated [83,92].Drawbacks in comparison to the classic SO 2 camera might be the higher complexity of hardware, as well as the fact that it usually takes much longer (minutes rather than seconds) to acquire a full image of the plume, also the computation time for spectral evaluation is longer.
Both the push-broom and whiskbroom applications can be employed for the I-DOAS technique.The Whiskbroom approach is usually implemented by adding 2-D scanning entrance optics to a DOAS instrument consisting of either two motors turning two mirrors (or prisms) or one motor and moving the entire instrument in one direction.The push-broom approach (see Figures 7B and 9) uses a 2 dimensional CCD while the second spatial direction can be implemented by a scanning mirror or by moving the instrument (e.g., [73][74][75][76]).

Combining both Approaches
Obviously, it would be very desirable to combine the good specificity of dispersive techniques (see Section 3.3) with the speed (and ideally simplicity) of the non-dispersive approaches (see Section 3.2).In fact, there exist several technical solutions which are essentially non-dispersive and thus fast and simple while making use of the details of the spectral features of the gas to be measured.Two relatively popular (although not in the area of volcanic gas measurements) solutions are Fabry-Pérot Interferometers and gas correlation sensors.Both show much superior discrimination power than simple filter cameras while being rather simple and fast devices.In particular both approaches allow high-speed (seconds) two-dimensional imaging of trace gas distributions.
Fabry-Pérot Interferometers can be seen as filters with a periodic transmission pattern, i.e., transmission maxima are regularly spaced at given wavelength intervals ∆λ, while the width of each transmission maximum is δλ (with δλ < ∆λ).By manufacturing a Fabry-Pérot interferometer the quantities ∆λ and δλ can be tailored within rather wide limits, in particular the can be chosen such as to match the periodic structure of the absorption cross section of many trace gases.Thus, such a device (see Figure 8, case b) is most sensitive for a gas with a given separation between the maxima of the absorption cross section.This is described in detail by Kuhn et al. [79].A demonstration of an actual device built on this principle is reported by Kuhn et al. [80], the same authors describe an extension of the principle for the detection of BrO (and other gases) in volcanic plumes.
A similar general idea is behind the gas correlation sensor (see e.g., [93]), which consists of a cuvette containing the gas to be measured (e.g., SO 2 ) mounted within a camera (see Figure 8, case d).The cell can be moved in and out the optics and the intensity ratio of two images recorded with the cell in the camera and the cell removed contains the trace gas signal.The principle relies on the fact that, while the intensity is always reduced when the cell is moved in, the intensity reduction is smaller when the gas is present in the observed scene.This is due to the fact that the intensity I(λ) reaching the camera is already attenuated at certain wavelengths (i.e., where the gas has its absorption maxima) and thus cannot be attenuated much more by the gas in the cuvette.Although the technique is essentially limited to gases that are sufficiently stable to be sealed into the absorption cuvette, which is part of the instrument, it is technically simple, and lends itself to imaging applications.Imaging by the gas correlation method applied in the IR are reported e.g., by Sandsten et al. [88].

Volcanic Gas Flux Determination
Usually the gas flux from a volcanic source is more important than the gas concentration (or column density through the plume), since the former is an indicator of the activity while the latter quantities are influenced by (varying) dilution and other processes which are not related to volcanic activity.

The Principle of Volcanic Gas Flux Determination
As described above, conventional one-dimensional (1-D) spectroscopic measurements typically measure the trace gas column density of the species of interest along a single line of sight.Integration of the gas concentration along this line occurs intrinsically (see Equation ( 6) in Section 2.1), and the measured quantity is the number of molecules along the line of sight.In order to derive an emission rate, the number of molecules Q in a cross-section of the plume perpendicular to its propagation direction is needed.Conventional 1D instruments measure this quantity by traversing under the plume or scanning the viewing direction across the plume and integrating the derived column densities (e.g., [15]).
The procedure is sketched in Figure 10 it consists of measuring a series of column densities S i through the plume (ideally an infinite number of column densities) and integrating over the derived columns within the width w of the plume (Equation ( 18)).
From the quantity Q the emission rate J can be readily determined by multiplication with the plume speed v P .
While determining the position y across the plume is no problem when traversing under the plume (e.g., with a spectrometer mounted on a car) for instruments scanning the plume from a fixed point (e.g., all NOVAC stations work like this, see [10]) only the observation angle is known.In this case an estimate of the distance L P to the plume is required and the perpendicular position y can be derived from geometrical considerations (in the simplest case as ∆y = L P •∆γ, where ∆γ is the change in scan angle).
The plume speed (i.e., wind speed v P at the position of the plume) and the wind direction can be determined in a number of ways: (1) From local measurements (2) From large scale wind fields, which are available from regional or global data bases (e.g., ECMWF or MERRA).(3) From measured correlation data within the plume itself (see e.g., [18,64,94,95]) Approaches ( 1) and ( 2) suffer from various shortcomings: Local wind speed measurements at the crater rim are usually not available and wind speed measurements further away are frequently disturbed by orographic influences.Wind speeds from large scale wind fields usually are on a very coarse grid (e.g., 1 • by 1 • ), while they are typically less affected by local orography these data may just not be representative for the site of the volcano.In contrast to that approach (3) is based on the fact that the gas (e.g., SO 2 ) emission strength from a volcano usually varies considerably (typically several 10%) with time on time scales of seconds to minutes as sketched in Figure 11.Thus, the column densities S 1 (t) and S 2 (t) are measured at two points within the plume which are located at different distances d 1 and d 2 downwind from the volcanic source (i.e., usually the crater).From the two time series the time lag ∆t P for maximum correlation between the two time series (which typically may encompass a time period of several minutes) is calculated.From ∆t P and the difference d 2 -d 1 in downwind distance of the two measurement points the plume speed is given by: For the above reasons approach (3) is the preferred method of determining the wind speed in the plume within the NOVAC network.Imaging techniques, however, have the advantage that they already record two spatial dimensions.Therefore, scanning is only required if the imaging technique itself is based on scanning technology (e.g., Imaging DOAS, see Section 3).Trace gas cameras (e.g., SO 2 cameras described above and in Section 3) based on band-pass filters do not require scanning.Instead, the number of SO 2 molecules Q in the plume cross-section can be derived by spatially integrating the recorded column images across the plume width w, perpendicular to the plume direction (see Figure 10).In case the plume cross-section is not taken perpendicular to its propagation direction a trigonometric correction has to be applied as described by several authors (e.g., [10,82,88]) and is discussed in detail in Section 4.2 below.
In this case the plume speed can be determined directly from the camera data through analysis of sequentially recorded images (e.g., [68,69,84,96]), providing a significant advantage over non-imaging techniques, which typically employ independent measurements of the wind speed at the appropriate plume height.Special software packages were developed e.g., by Valade et al. [97] and a direct validation of thus derived SO 2 emissions were reported by Smekens et al. [84].

A More Detailed View-Determination of the Wind Direction
The simple approach sketched in Section 4.1-while applicable directly in certain cases-must be extended to compensate for particular conditions of the measurement: (1) Light dilution may affect the accuracy of the column density retrieval (2) Multiple scattering inside the plume can also affect the column density retrieval (3) The effect of the plume propagation direction being non-perpendicular to the viewing direction must be corrected The first two effects are discussed in several publications, e.g., [30,82,83,92,99,100] with the bottom line being that corrections are possible on the basis of a detailed analysis of the recorded trace gas (e.g., SO 2 ) spectra.Correction of the third effect is addressed in the following.The problem lies in the correct determination of the plume propagation direction.While this is usually not a problem when the trace gas column density is determined by a zenith looking instrument during traverses, since in these cases the position of the source (e.g., the crater) and the position of the measurement are usually well known.Consequently the angle α between the plume propagation direction (i.e., the wind direction) and the direction of the traverse is also known with high precision.
However, if measurements are made from a fixed point, either by a scanning system or a two dimensional trace gas camera (e.g., a SO 2 Camera), then the angle α between the plume propagation direction and the scanning direction or the image plane of the camera are not a priori known.
As sketched in Figure 12 there are several effects if the plume propagation direction is not perpendicular to the scanning direction or not in the image plane, respectively.Following the discussion in Klein et al. [22], we consider data derived by a plume imaging system.For simplicity we initially also assume a cylindrical plume (i.e., a plume where the cross section does not change with distance from the origin and with time).Later we will generalize for an arbitrary plume cross section and its evolution.As can be seen in Figure 12 top view (lower left panel) in the centre of the image plane the determined flux will be independent of the "tilt" angle α between image plane and plume propagation direction.This is due to the fact that two effects cancel (as noted e.g., by Mori and Burton [96]): (1) The length of the light path through the plume increases as 1/cos(α).
(2) The determined apparent wind speed is reduced by the factor cos(α) since the determined d 2 -d 1 appears shorter by this factor (see bottom right panel of Figure 12) while the determined time lag ∆t P stays the same.
However, this only holds for one point (e.g., the centre) within the field of view.Across the field of view there are two effects: (1) Still the light path through the plume is larger towards the edges of the image, while the determined velocity v P stays the same.Thus the flux appears to increase somewhat towards the edges of the image compared to the centre.(2) A further effect of at "tilt" is due to geometry in that the closer part of the plume appears larger than the part which is further away from the camera (see top right panel of Figure 12) thus the integral (Equation ( 18)) will extend over a larger extent and thus be larger.
When image series ("movies") of the plume are available these effects can be fortunately not only corrected, but rather the data can be used to determine the "tilt" angle α of the plume.This is done by (1) determining the flux in each image column of the image (where the plume is covered) and ( 2) making use of the usually well justified assumption that the amount of trace gas (e.g., SO 2 ) within the plume is conserved over a time period of a few minutes.Therefore the flux must be the same across the entire field of view (once the effect of the centre to edge enhancement (see above) is corrected for).If the thus corrected flux is not constant under the initial assumption of α = 0 then α is varied until the flux is constant across the field of view (see Figure 13).

Sample Applications
Technical advances in remote sensing and volcanic gas imaging have led to high temporal and spatial resolution data, which now allow the investigation of volcanic source behaviour taking place below the Earth surface by combining the surface observation of gas release with geophysical data and fluid dynamic models for melt and atmosphere including numerical models and laboratory data.Apart from the volcanological application also the volcanic plume chemistry and therefore the impact on the atmosphere can be investigated in more detail for instance by studying the heterogeneity of volcanic plumes and also by visualizing their mixing with the atmosphere.However, to our knowledge this possible application has not been realised so far.In the following section we will illustrate the new possibilities for volcanological investigations using some recent example studies.

UV Spectroscopy
Measurements of gases emitted from volcanoes, in particular the amount of SO 2 , have traditionally taken place using UV-light applications, starting with COSPEC [28] in the 70s and after 2002 measurements were more widely undertaken by applying DOAS technique.However, both approaches were rarely used for two-dimensional imaging and SO 2 fluxes were determined with a time resolution of >5 min.An exception are the wide field of view UV spectrometer applications [101,102].Highly improved time resolution and reduced error for plume velocities are possible since the onset of the SO 2 camera application [82,96].Using an SO 2 camera leads to the determination of SO 2 emissions with time resolution of the order of a second and therefore it opens the possibility to correlate gas emission measurements directly with geophysical signals (tremor, very long period events etc.) and thus adding information to test current models of physical processes below the Earth surface.
A recent example is the work of Moussallam et al. [103], describing different degassing patterns of lava lakes in relation with their underlying conduit flow processes (See Figure 14).The authors point out that there is no similarity between gas emissions of three out of four volcanoes (Erebus, Ambrym, Kilauea, Villarica) containing lava lakes, which were all studied with a comparable time resolution.The reason proposed are differences in magma viscosity and conduit geometry, which cause the different character of emission flux behaviour.No regular periodicity of gas fluxes as well as composition could be determined in data of Villarica volcano, Chile.Those surficial gas measurements, which were combined with thermal infra-red data, point to flow instabilities inside the conduit and a turbulent mixing between the ascending gas rich and the descending gas poor magma.Part of the observed high frequency variability of the emission strength of Villarica volcano is suggested to originate from atmospheric transport processes.Instead, periodic oscillations of flux and composition were found at Erebus and Ambrym (e.g., [101,104]) indicating, combined with geophysical data, a bidirectional magma flow, whereas at Kilauea [105], oscillation of degassing amounts has been interpreted as shallow degassing patterns due to gas percolation.However, not only has lava lake degassing been investigated so far, high time resolution SO 2 camera measurements were also done to investigate puffing and strombolian explosions at various volcanoes, showing similarly large differences in degassing characteristics of different volcanoes [19,106,107], came to the conclusion that explosions are only second order degassing processes compared with quiescent degassing.Those studies were undertaken at Stromboli and Assama volcano, respectively, showing that less than 16% of the gas amount is emitted by explosions.
In contrast studies at Fuego volcano by Waite et al. [108] suggest that 95% of the total degassing amount is contributed by explosions.
Table 2 is an expanded overview of the one of Burton et al. [109].Summarizing comparison studies between SO 2 camera-derived SO 2 fluxes and some geophysical parameters.+: Correlation is reported; −: no correlation has been found; 0: the parameter were independently reported and used together for a volcanological interpretation no investigation of the correlation between the parameter has been presented.

IR Spectroscopy
Although IR techniques have complex analytical requirements and still need relatively bulky and cost intensive instruments, IR techniques are highly desirable to investigate volcanic gas compositions, because nearly all major volcanic gas compounds show characteristic absorption bands in the infra-red.
Unfortunately, the two main volcanic gases, water vapour and CO 2 , although they have strong absorption features in the near IR are quite challenging to detect, mainly due to their high atmospheric background and its spatial variation.Nevertheless, the value of such information has been already demonstrated by using light path arrangements which avoid large atmospheric distances, i.e., light paths close to the source and using either artificial lamps or hot ground surfaces (e.g., lava lakes).For instance, at Erebus (e.g., [117]), Etna (e.g., [36,118,119]), Ambrym ( [104]), and Stromboli (e.g., [120,121]), high frequency variations of major plume composition could be successfully determined over time frames of hours to months.Further, progress has been done during the last few years (e.g., [46,53]) and even from space volcanic CO 2 emissions of at Yasur (on Tanna Island, Vanuatu), Aoba and Ambrym volcanoes (Vanuatu) have been detected by OCO-2 [122].Although these last examples are far from routinely usable measurements they show new technical developments for future research.
An interesting example for the advanced understanding of physical processes inside the feeding system of a volcano is the work by La Spina et al. [121].Only with spatial and high time resolved investigations by FT-IR spectroscopy the authors could show that the previously proposed explanation of CO 2 flushing, a commonly proposed process, is unlikely in the case of Stromboli.The authors could come to this conclusion by observing the differences in the composition of the three investigated craters at Stromboli (CC, NEC, SWC) and the high H 2 O/CO 2 ratios.They therefore suggested an alternative model (see Figure 15) to explain the enhanced CO 2 /SO 2 ratios weeks before major explosions and the decreasing CO 2 /SO 2 ratios during the hours just before this activity (without the need of CO 2 flushing) and also give an alternative explanation how a crystal pure magma may be transformed into a crystal rich one.This later fact could be explained by H 2 O exsolution and crystallization driven by depressurization when magma rises to the surface.La Spina et al. [121] assume that the increase of the CO 2 /SO 2 ratio is caused by more magma input into the shallow system.This should lead to a drop of permeability in the deep system which then causes the decrease of CO 2 /SO 2 ratios as the mixture of deep degassing and shallow degassing is shifted to a larger percentage of shallow degassing.According to the calculation of the authors this seems far more realistic than the incredible large amounts of CO 2 needed to produce significant isobaric dehydration to reproduce the measured H 2 O/CO 2 values.Further, La Spina et al. [121] utilize the observed spatial and temporal heterogeneity to produce a new model describing the geometrical structure of the shallow plumbing system of Stromboli (see Figure 15).carries the largest contribution of gas produced in closed system degassing, and is therefore richest in CO 2 , while the NEC crater has the weakest connection to the closed-system gas source, and is therefore relatively poor in CO 2 .The branching depth of 30-40 MPa is the minimum depth for which such an CO 2 /SO 2 fractionation could occur between the SWC and other craters.The final gas observed from each crater is therefore the superposition of different relative amounts of the three gas sources: CO 2 -rich gas from depth, H 2 O-SO 2 rich gas released from the ascending magma in the conduit and H 2 O-rich gas released from the shallow accumulation zone.Adapted from Figure 6 of [121].
In contrast to the observation of [121] at Stromboli, in another study of Oppenheimer et al. [123] did not find a significant variation of the CO 2 /SO 2 ratio at Yasur volcano between explosive and passive degassing activity.However, large variations (more than a magnitude) are observed in the HCl/SO 2 ratio with a ratio of about two during quiescent degassing and about 30 during explosions.These results are explained with gas accumulation below the HCl exsolution.
In addition, less abundant gases (e.g., SiF 4 ) where recently successfully studied using FTIR.After the first measurements of SiF 4 in a volcanic plume by Francis et al. [124], two dimensional SiF 4 distributions have been presented by [45] using thermal emission spectroscopy.Taquet et al. [125] also used thermal emission based spectroscopy to study the SiF 4 /SO 2 ratio over a period of 6 month during dome growth at Popocateptl.The authors found high variability of this ratio in coincidence with seismic data.SiF 4 is a secondary gas formed due to the interaction of HF with rocks or ash.Nevertheless, there is a high potential to study certain types of activity (e.g., dome growth) and eventually using this approach for monitoring in the future.In the article by Taquet et al. [125] the variation of the SiF 4 , SO 2 and the ratio SiF 4 /SO 2 helped to decipher the different processes related to the dome permeability.

Summary and Outlook
There is a clear development from occasional volcanic gas concentration measurements (typically by "wet" chemical techniques or gas chromatography) to continuous observations as pioneered by the COSPEC technique (e.g., [28]) and brought to a large scale by the NOVAC network (e.g., [10]) and other local networks (e.g., [109,126]).In addition gas concentration (or column density) measurements have been extended to flux measurements (see e.g., NOVAC).Another dimension is marked by the extension of SO 2 measurements to continuous multi-species observations.In this context also extremely affordable spectrometers based on mass produced smart phone imagers [127] can play a role.
From here a next quantitative jump (or paradigm shift, see [16]) is marked by gas imaging techniques like the SO 2 camera (see e.g., [109]), which allow-compared to scanning approaches-a much more detailed analysis of volcanic emissions and plume behaviour.Besides more precise flux measurements for instance short term periodicity in emission flux [16] or the angle between plume direction and image plane can be derived.Moreover, the spatial inhomogeneity of volcanic sources can be investigated.
High Frequency imaging of volcanic plumes was first realized with UV cameras and instruments based on this technique can be produced rather cost effectively and-as was recently shown by Wilkes et al. [128] have the potential to become very cheap instruments which may be ubiquitous in volcanic research.However, the approach has two drawbacks.
(1) The relatively crude spectroscopic approach, basically just ratioing the intensities in two wavelength intervals leads to interferences by aerosol, (stratospheric) ozone, and possibly other species.Moreover, radiation transport issues may also influence the accuracy of the technique (see e.g., [71,82,129]).(2) UV cameras rely on sunlight and can thus only operate during daylight hours.
There are techniques under development (or being adapted for the purpose of volcanic observations) which promise to overcome these weaknesses.These include gas correlation spectroscopy (see Section 3.4), titled filter imaging, and Fabry Pérot imaging (see [79,80]), which both make better use of the spectral finger print of the target gas than the classic UV camera, thus reducing interferences and potentially improving the accuracy of the measurement (see also [5]).
Infra-red thermal emission spectroscopy is independent of daylight and may be the technology of the future in either scanning or imaging applications.At present cost and logistic requirements (cooling of the IR detector, use of interferometers with moving parts) are still high.Here technological development aimed at lowering complexity and cost of the instruments would be very welcome.
At present, imaging Michelson interferometers exist and are used for certain purposes, but the instruments are still extremely expensive.Alternatively, the UV-camera approach can be extended to the thermal IR region (see e.g., [56]), while this approach suffers from similar problems as the UV camera (see above) the rapidly falling prices of thermal IR cameras may make this approach feasible for volcanological applications.
With all spectroscopic techniques radiation transport issues are important and may be limiting the achievable accuracy (see Section 2).Here, better radiation transport models and new approaches like using polarisation of the radiation and ratioing techniques based on the absorption due to "known" absorbers like O 2 or O 4 (oxygen dimers, see e.g., [130,131]) will be helpful.
Overall, it becomes clear that spectroscopic plume analysis is just at the beginning of its development and it constitutes a new and rapidly evolving technique to analyze magmatic degassing, which is (coupled with viscosity and crystallisation) the main driver of volcanic processes.

Figure 1 .
Figure 1.Overview of the different plume sensing schemes: (A) Thermal emission is recorded by passive spectroscopy in the thermal infra-red; (B) Direct light absorption spectroscopy in the UV, visible, or near IR; (C) Scattered sunlight spectroscopy in the UV or visible spectral ranges; (D) Active absorption spectroscopy in the UV, visible, or near IR.

Figure 2 .
Figure 2. Schematic illustration of the classic absorption experiment with a light source at a distance L from the observer and the absorbing volcano plume in-between.

Figure 4 .
Figure 4. Family tree of spectroscopic remote sensing schemes.

Figure 5 .
Figure 5. Schematic illustration of scattered skylight DOAS measurements with the sun as a light source.

Figure 6 .
Figure 6.Schematic illustration of thermal emission techniques exploiting the difference between lines-of-sight pointing through the volcanic plume and those pointing into the background sky.

Figure 8 .
Figure 8. Trace gas camera schemes: The basic set-up consists of a (UV) lens, aperture, wavelength selective element (WSE) and suitable 2-D detector.The WSE can either be mounted in front of the lens (dotted box) or between lens and detector (dashed box), the relative merits of the two approaches are discussed by Kern et al. [82].The different approaches are distinguished by their WSE, which can be one of the following: (a) two narrow band filters alternatingly being brought into the beam; (b) a Fabry-Pérot interferometer with adjustable transmission wavelength; (c) a narrow-band interference filter, which can be rotated; (d) a cell (cuvette) containing the gas to be measured, periodically introduced into the light beam.

Figure 9 .
Figure 9. Schematics of an Imaging-DOAS (I-DOAS) whisk-broom instrument and sketch of the evaluation procedure: The I-DOAS instrument simultaneously creates spectra for each pixel in a column of the image (left), these are recorded by a 2-dimensional CCD (centre).Upon DOAS evaluation trace-gas column density values for each pixel in the image column are derived.After completion of a scan all image columns are combined into a two-dimensional image (right), which is typically presented as "false colour image".

Figure 10 .
Figure 10.Schematics of flux determination from a plume scan.

Figure 11 .
Figure 11.Sketch of the determination the wind speed v P inside a plume.Left panel: Two time series (continuous red and black lines) are recorded at two different positions (difference d 2 − d 1 ≈ 900 m) downwind the plume of Mt.Etna.The dotted red lines indicate the shifting time series 2 versus series 1.In the graph, both series are low-pass filtered, to make the relevant structures clearly visible.Right panel: Correlation coefficient between the two time series as a function of time lag.The highest correlation is found for a time shift of ∆t = 59 ± 6.5 s resulting in a wind speed of 15.2 ± 1.6 m/s.Adapted from [98].

Figure 12 .
Figure 12.Schematics of flux determination from a plume scan.See Section 4.2 for details.

Figure 13 .
Figure 13.Sample application of the plume "tilt" correction on a set taken at Mount Etna on the 9th of July 2014.Top panel: Deviations of the apparent SO2 fluxes in three different vertical cross sections through the plume.These apparent deviations are caused by wrongly assuming an inclination ("tilt") of the plume with respect to the image plane of α = 0. Using the conservative flux condition an inclination angle of α = 38 • was determined.Figure from [22].

Figure 14 .
Figure 14.Three suggested degassing mechanisms (Kilauea, Villarica, and Erebus) with different conduit dynamics caused by the difference in magma input and viscosity.After [103], Figure 7.

Figure 15 .
Figure15.Illustrating a schematic model of the steady-state degassing system on Stromboli (described in the text).Bottom of figure: Point of transition to open system degassing, the gas exsolution is colour coded as explained in the box.The SWC crater (top left of figure) carries the largest contribution of gas produced in closed system degassing, and is therefore richest in CO 2 , while the NEC crater has the weakest connection to the closed-system gas source, and is therefore relatively poor in CO 2 .The branching depth of 30-40 MPa is the minimum depth for which such an CO 2 /SO 2 fractionation could occur between the SWC and other craters.The final gas observed from each crater is therefore the superposition of different relative amounts of the three gas sources: CO 2 -rich gas from depth, H 2 O-SO 2 rich gas released from the ascending magma in the conduit and H 2 O-rich gas released from the shallow accumulation zone.Adapted from Figure6of[121].

Table 1 .
Characteristics of the three basic imaging techniques.

Table 2 .
Comparison studies between SO 2 camera-derived SO 2 fluxes and geophysical parameters.