Reaction Rates Control High-Temperature Chemistry of Volcanic Gases in Air

When volcanic gases enter the atmosphere, they encounter a drastically different chemical and physical environment, triggering a range of rapid processes including photochemistry, oxidation and aerosol formation. These processes are critical to understanding the reactivity and evolution of volcanic emissions in the atmosphere yet are typically challenging to observe directly at the lava-atmosphere interface due to the nature of volcanic activity. Inferences are instead drawn largely from observations of volcanic plumes as they drift across a crater’s edge and further downwind and the application of thermodynamic models that neglect reaction kinetics as gas and air mix and thermally equilibrate. Here, we foreground chemical kinetics in simulating this critical zone. Volcanic gases are injected into a chain-of-reactors model that simulates time-resolved high-temperature chemistry in the dispersing plume. Boundary conditions of decreasing temperature and increasing proportion of air interacting with volcanic gases are specified with time according to an offline plume dynamics model. In contrast to equilibrium calculations, our chemical kinetics model predicts that CO is only partially oxidised, consistent with observed CO in volcanic plumes downwind from source. Formation of sulfate precursor SO3 at SO3/SO2 = 10-3 mol/mol is consistent with the range of reported sulfate aerosol to SO2 ratios observed close to crater rims. High temperature chemistry also forms oxidants OH, HO2, and H2O2. The H2O2 will likely augment volcanic sulfate yields by reacting with SO2(aq) in the cooled-condensed plume. Calculations show that high-temperature OH will react with volcanic halogen halides (HBr, HCl) to yield reactive halogens (Br, Cl) in the young plume. Strikingly, high-temperature production of radical oxidants (including HOx) is enhanced by volcanic emissions of reduced gases (CO, H2, H2S) due to chemical feedback mechanisms, although the kinetics of some reactions are uncertain, especially regarding sulfur. Our findings argue strongly that the chemistry of the hot near-source plume cannot be captured by equilibrium model assumptions, and highlight the need for development of more sophisticated, kinetics-based, high-temperature CHONS-halogen reaction models.


INTRODUCTION
Volcanoes release gases and aerosols to the atmosphere through both quiescent (passive) degassing, and effusive and explosive eruptions. Volcanic emissions of SO 2 to the stratosphere are well known to impact climate through the atmospheric oxidation of SO 2 to form radiatively active sulfate aerosol particles. In fact, volcanoes release a wide range of highly reactive gases that can impact atmospheric composition and climate. Volcanic bromine and chlorine emissions undergo plume atmospheric chemistry that can destroy both tropospheric and stratospheric ozone. In descending order of typical abundance (which varies with volcano setting, magma redox conditions, and eruptive style) volcanic gas emissions consist of: H 2 O, CO 2 , SO 2 , HCl, HF, H 2 S, OCS, CO, and HBr, as well as other trace species e.g., metals. Observations of volcanic plumes identify several additional species (e.g., NO, NO 2 , HNO 3 , BrO, OClO, SO 4 2− , HO 2 NO 2 , and H 2 O 2 ), (Allen et al., 2000;Mather et al., 2004b;Bobrowski et al., 2007;Oppenheimer et al., 2010;Carn et al., 2011;Martin et al., 2012;Kern and Lyons, 2018). These species are formed by oxidizing chemical reactions as the magmatic gases mix with air, first at high temperatures near to the source and then at low temperatures as the cooled plume disperses further into the background atmosphere. Numerical models of low-temperature atmospheric chemistry of the plume have been developed (e.g., Bobrowski et al., 2007;Roberts et al., 2009;von Glasow, 2010;Jourdain et al., 2016) that are able to reproduce some -but not all -of these observed oxidized species. This points to an important role of high-temperature chemistry. Current approaches apply thermodynamic equilibrium models to represent the high-temperature chemistry of the near-source plume (e.g., Gerlach, 2004;Martin et al., 2006). However, thermodynamic equilibrium calculations may be misleading if the chemistry is limited by kinetics (rates of reactions) in the rapidly cooling volcanic gas and air mixture within seconds of emission. Here, a chemical kinetics approach is developed to simulate the time-varying chemical processing of C-H-O-S gases in the hot near-source plume on its release and mixing with background air.

High-Temperature Products in Volcanic Plumes: Observations and Modeling
Observational evidence for the high-temperature production of oxidants in volcanic plumes includes SO xy , HO xy , NO xy , and rapid formation of reactive halogens. Sulfate-rich aerosols have been observed at the crater-rim or in very neardownwind plumes of many volcanoes (exception: volcanic domes), Table 1. Sulfate is reported relative to co-measured SO 2 which serves as a convenient yet only approximate plume tracer, allowing observations to be normalized for dilution effects. The sulfate:SO 2 ratios typically range from 0.00002 up to around 0.01 mol/mol, Table 1 and references therein. Many of the ratios are too high to be explained by atmospheric oxidation of SO 2 at low-temperatures on the plume transport timescale of seconds to minutes (e.g., Galeazzo et al., 2018). Instead, the observations indicate a sulfate-rich aerosol formed very close to source. One hypothesized pathway to its formation is high-temperature oxidation of a small fraction of SO 2 in the near-source plume to produce sulfate precursor SO 3 (Allen et al., 2000). For volcanic domes, sulfate formation may be enhanced by circulation of air inside the dome (Zelenski et al., 2015). The range in observed SO 4 2− :SO 2 indicates variability in the near-source sulfate formation processes at different volcanoes, as well as possibly some measurement uncertainty. In all cases sulfate contributes a minor fraction to the total sulfur emission that is typically dominated by SO 2 and H 2 S (the proportion depending strongly on magma redox conditions; Oppenheimer et al., 2014). Nevertheless, at source (or so-called "primary") volcanic sulfate is critical to plume processes such as multi-phase halogen reactions (discussed further below) and acts as a condensational sink that can inhibit new particle formation (Sahyoun et al., 2019).
Frontiers in Earth Science | www.frontiersin.org 2 July 2019 | Volume 7 | Article 154 originated from the self-reaction of volcanic HO 2 . Oppenheimer et al. (2010) measured elevated HO 2 NO 2 in the (cooled and dilute) Mt Erebus plume where the cold Antarctic atmosphere likely contributed to its thermal stability. HO 2 NO 2 is formed by the reaction of HO 2 with NO 2 therefore its presence provides evidence for volcanogenic sources of both HO x and NO x .
Volcanogenic HO x has also been observed by a chance encounter of an instrumented aircraft with the Hekla 2000 eruption plume in the lower stratosphere (Rose et al., 2006). Observations showed OH and HO 2 were elevated at up to 0.5 and 15 pptv above background in 33-34 h aged volcanic cloud during night. These measurements have eluded any explanation given the lack of photochemistry, a cooled plume, and anticipated fast HO x destruction by reaction with volcanic SO 2 and halogens (both in excess abundances). In general, the reported elevated abundances of HO xy in volcanic plumes compared to the background atmosphere are surprising because low-temperature volcanic plume chemistry is expected to deplete HO xy (Galeazzo et al., 2018). Volcanic HO xy may originate from high-temperature chemistry of the near-source plume. Elevated concentrations of NO, NO 2 , and HNO 3 have been measured in several volcanic plumes in the troposphere in addition to the HO 2 NO 2 observed in the Mt Erebus plume: a review is given by Martin et al. (2012) and references therein, reporting a typical molar ratio for NO y /H 2 O of 10 −5 . Whilst filter pack sampling identified HNO 3 in crater-rim plumes from several volcanoes (e.g., Mather et al., 2004a), on other occasions there was no evidence for HNO 3 , even at the same volcano (e.g., Martin et al., 2010). This may point to variability in emissionsprocessing or reliability of measurement technique. Nevertheless, real-time aircraft-based measurements by Atmospheric Chemical Ionization Mass Spectrometer confirm volcanic HNO 3 in plumes downwind from Mt Etna at HNO 3 /SO 2 = 0.02 (Voigt et al., 2014), and in the Mt Erebus plume (Oppenheimer et al., 2010). Sources of volcanic NO xy are uncertain. Martin et al. (2012) present a kinetic model of CHON chemistry that demonstrates high-temperature formation of thermal NO x is too slow to occur in the near-source plume that rapidly disperses and cools over seconds-minutes. Volcanic lightning is a likely source of NO x in clouds from explosive eruptions. It has also been hypothesized that NO x may sometimes be formed by oxidation of ammonia, for example in the stratospheric cloud from the Hekla 2000 eruption (Rose et al., 2006).
The need to characterize high-temperature chemistry in the near-source plume is underlined by model studies of volcanic halogen transformations in the plume (Bobrowski et al., 2007;von Glasow, 2010;Roberts et al., 2014). Halogens are emitted as hydrogen halides (HBr and HCl) from volcanoes but are transformed into reactive halogens (BrO and OClO) in the tropospheric plume via a multi-phase autocatalytic "bromine-explosion" chemistry occurring at low-temperatures. However, numerical models of the plume atmospheric halogen chemistry can only reproduce the observed magnitude and rate of formation of BrO if the volcanic emission is modified to include radicals believed to form at high-temperatures in the near-source plume (the so-called "effective source region"). Radicals such as OH (that react with HBr, HCl to produce Br, and Cl) provide an initial source of reactive bromine to initiate the low-temperature bromine explosion cycles. Furthermore, high-temperature production of sulfateprecursor SO 3 generates near-source (so-called "primary") sulfate aerosols that promote heterogeneous reactions that drive the bromine explosion chemistry. Observations of BrO (and occasionally OClO) in tropospheric volcanic plumes globally (e.g., Bobrowski et al., 2007;Boichu et al., 2011;Hörmann et al., 2013;Donovan et al., 2014;Kern and Lyons, 2018) indicate pervasive emissions of volcanic halogens to the troposphere accompanied by both high-and low-temperature chemistry in volcanic plumes.
The high-temperature region of the plume also facilitates gasash interactions including adsorption of gases such as SO 2 and HCl onto ash particles that may then fall-out (Delmelle et al., 2018 and references therein). High-temperature SO 3 is a precursor to sulfate that can condense as the eruption column cools and could change ash properties (e.g., Hoshyaripour et al., 2015) and ashrelated atmospheric impacts e.g., regarding plume lightning or ice nuclei. Whilst such heterogeneous interactions are outside the scope of our study, quantifying the high-temperature oxidizing chemistry occurring in the gas-phase as volcanic gases mix with air is a fundamental first-step to characterising the nearsource plume composition and its possible influence on gasash interactions.
Thermodynamic models such as HSC have been applied to calculate the equilibrium compositions of mixtures of magmatic gases with air, and HSC output has been used in numerical model studies of the plume atmospheric chemistry and impacts (e.g., Gerlach, 2004;Martin et al., 2006Martin et al., , 2007Martin et al., , 2009Aiuppa et al., 2007;Bobrowski et al., 2007Bobrowski et al., , 2015Roberts et al., 2009Roberts et al., , 2014Roberts et al., , 2018von Glasow, 2010;Martin and Ilyinskaya, 2011;Hoshyaripour et al., 2012;Ilyinskaya et al., 2012;Surl et al., 2015;Jourdain et al., 2016). HSC is used to calculate the composition of the nearsource plume for temperatures above 600 • C according to Gibbs free energy minimization. The emission is released at magmatic temperatures and as the magmatic gases entrain air, the mixture cools toward background temperature. The calculation assumes thermodynamic equilibrium composition is reached at a specified temperature and degree of mixing of air with magmatic gases (often expressed as a volumetric or molar ratio, VA:VM where VA is air volume and VM is magmatic gas volume). Plume composition is then assumed to be "quenched" with respect to high-temperature chemistry upon further mixing with air. Choice of the VA:VM "quenching-point" for the HSC calculation is rather uncertain and may vary with eruption/emission type. For studies of passively degassing volcanoes it was found that a high VA:VM = 40:60 (Bobrowski et al., 2007) yielded extremely high SO 3 relative to SO 2 (much greater than observed sulfate/SO 2 ratios, see Table 1), therefore more recent studies chose lower VA:VM = 10:90 or 5:95 (e.g., Roberts et al., 2014) as the point where high-temperature chemistry is quenched. Conceivably, large eruptions might sustain high-temperature conditions over longer periods allowing the chemistry to reach equilibrium conditions at higher VA:VM. Hoshyaripour et al. (2012) simulated the high-temperature chemistry in an eruption column by co-varying VA:VM from 0 to 1 with temperature decreasing FIGURE 1 | Schematics of (A) high-temperature plume chemistry that generates radicals and oxidized products in the plume and (B) our chain-of-reactors approach to simulate the processes based on chemical kinetics under time-varying temperature (T) and proportion of mixing of air with magmatic gases (VA:VM).
from 1273 K to 873 K (1000 • C to 600 • C), predicting that the eruption column acts as a hot oxidizing reactor for S species.
In general, volcanic plume applications of thermodynamic equilibrium models such as HSC yield elevated abundances of SO x (dominated by SO 3 ), HO x (dominated by OH) and NO x (dominated by NO), and halogen radicals Br, Cl (e.g., Martin et al., 2006). The relative abundances of these oxidized products increase with increasing VA:VM. The equilibrium models also predict near-complete oxidation of reduced gases emitted from volcanoes (e.g., CO, H 2 S, and H 2 ), yet these gas species are often detected at the crater-rim or in downwind (cooled) plumes (e.g., Aiuppa et al., 2005;Schumann et al., 2011;Moussallam et al., 2012;Roberts et al., 2012;de Moor et al., 2016). This suggests that some gas species may not always be under equilibrium control . Plume observations of reduced gases relative to major gases (e.g., CO/CO 2 , H 2 /H 2 O, and H 2 S/SO 2 ) are used to infer emission properties such as temperature and oxidative state, but such inferences cannot hold true if there is substantial oxidation of the reduced gases following emission, as predicted by HSC. Here, a method is developed to evaluate how reaction rates control high-temperature chemistry of volcanic gases in air as a function of time-varying plume dynamics, Figure 1.
Our model uses a chemical kinetics rather than thermodynamics approach, and requires time-resolved boundary conditions in both temperature and degree of mixing of magmatic gases with air (VA:VM), as well as the composition of the volcanic gases emitted as model input.
The temperature of the volcanic gas emission is uncertain. It is typically assumed to be identical to petrological magmatic estimates (600-1200 • C), although analysis of gas compositions measured at Kilauea, Hawaii (Oppenheimer et al., 2018) and Mt Etna, Italy (Aiuppa et al., 2011) show that gases may cool prior to emission into the atmosphere. Introduction of surface waters may also cool the gas emission. Physical properties of the plume are poorly constrained for the few seconds following release from the eruptive vent (or other emission sources e.g., lava lake or dome), as there are few direct measurements or dedicated model studies. A study of gas-puffing dynamics at passively degassing volcano Masaya yielded estimates of mean dilutions of 0.09 and 0.24 between emission and detection by thermal radiometer about 5-10 m above the vent, for gas puffs with mean velocities 2.5 and 4.3 m/s (Branan et al., 2008).
One way to estimate temperature and VA:VM in the nearsource plume can be through dynamical model simulations of plume rise (e.g., Mastin, 2007;Aubry et al., 2017). This study relies on one such eruptive case to supply physical boundary conditions: the PlumeRise model applied to the 2010 Eyjafjallajökull eruption plume (Woodhouse et al., 2013), for which observations of volcanic gases including CO/SO 2 have also been reported both close to the source and far downwind (Allard et al., 2011;Schumann et al., 2011). This Eyjafjallajökull case forms the basis for an exploratory investigation into the chemical kinetics of gas-phase reactions occurring at high-temperatures in volcanic plumes.
Frontiers in Earth Science | www.frontiersin.org

MATERIALS AND METHODS
Our approach to simulate high-temperature chemical reactions in the near-source plume includes three components: (i) the chemical kinetics model (chain-of-reactors), (ii) physical boundary conditions applied to the model, and (iii) the volcanic gas emissions injection as model input.

Chemical Kinetic Model (Chain-of-Reactors)
Chemkin is a fortran-based tool for incorporating complex chemical kinetics into model simulations of fluid dynamics (Kee et al., 2000 and related references from Reaction Design). Gas species and (reversible) reactions are first declared, as well as initial conditions (such as initial gas concentrations, temperature, pressure), and then the reactions are simulated in a plug-flow reactor (PFR) to evaluate temporal changes in gas composition.
Rates of progressions of reactions are quantified as the difference between the rate of the forward and reverse reactions, for example E1 for the reversible reaction A + B = C + D (where gases A, B, C, and D may also be involved in other reactions, each with its corresponding differential equation).
The rate constants for forward reactions, k f , are described in the form, E2, where R is the gas constant and T the gas temperature. The pre-exponential factor, A, the temperature exponent, β, and activation energy E, are specified in the mechanism file (see Glarborg and Marshall, 2013).
The reverse reaction rate constant, k r , is given by E3 where K is the equilibrium constant that is determined by thermodynamic properties (molar enthalpies and entropies as a function of temperature); these properties are specified in the thermodynamic data file (see Glarborg and Marshall, 2013).
The Chemkin model also accounts for third-body, unimolecular, chemically activated reactions, and pressuredependencies of reaction rates. The resulting system of ordinary differential equations and initial conditions is solved by Chemkin to quantitatively predict how the gas composition in the chemical reactor evolves as a function of time according to the kinetics-based reaction scheme and the thermodynamic conditions.
The reaction mechanism applied within ChemKin is userspecified. There are very few studies at high-temperature dealing with sulfur species oxidation and their interactions with carbon, nitrogen oxides, or halogens. Here, the mechanism proposed by Glarborg and Marshall (2013) was selected for its good performance in comparison to laboratory shock-tube experiments, as recently highlighted by Mathieu et al. (2017). It includes 38 species involved in 161 reversible reactions. This reaction mechanism, and corresponding thermodynamic properties, relies on earlier works on sulfur chemistry and was updated to model the oxidation of carbonyl sulfide (OCS). A full list of gas species and chemical reactions is given in Supplementary Tables S1, S2. This version does not include halogens or NO x chemistry. The interactions between sulfur, carbon, halogens, and nitrogen oxides are beyond the scope of this work, but will need to be addressed in future. Also, all the species considered were in the gas-phase; heterogeneous reactions or photochemical reactions were not included. This is a clear limitation of the present study which should be seen as a preliminary step needed to account for time evolution of chemicals in a near-source plume. Whilst the chemical scheme is one of the most complete regarding sulfur compounds, there exist some uncertainties in the high-temperature reaction kinetics, in particular related to H 2 S (Mathieu et al., 2017). A single PFR simulates the high-temperature chemistry for given starting conditions in gas composition and temperature. A timevarying temperature can be introduced by applying a temperature gradient in the PFR and imposing an initial velocity of the reactor gases. However, air cannot be added as a function of time into a single reactor, i.e., the proportion of air to magmatic gases in the simulation is held constant according to the initial composition.
Therefore, a chain-of-reactors approach was developed using the Chemkin Pro software to simulate the chemical evolution of a near-source plume that simultaneously cools (decreasing T) and entrains air (increasing VA:VM) as it dilutes, Figure 1. Thirty PFRs were placed in a chain. A continuous decreasing trend in plume temperature was imposed across the chain by defining the start and end temperatures for each PFR according to our chosen boundary conditions (details in the next section). The starting composition for each reactor was defined by taking the output from the preceding one and diluting it with air (78% N 2 , 21% O 2 , 1% Ar) according to time-varying VA:VM boundary conditions. The pressure decrease was also considered even if it was not found to be significant. This alternate simulation of the plume chemistry (within reactor boxes) and mixing-dilution steps (between reactor boxes) is similar to the approach used in kinetics-based modeling of the low-temperature atmospheric chemistry of volcanic plumes e.g., PlumeChem (Roberts et al., 2009).

Boundary Conditions to the Chain-of-Reactors Provided by PlumeRise Model
The chain-of-reactors model requires time-resolved boundary conditions for plume temperature and VA:VM as the emitted volcanic gases rapidly mix with air. Here, output from the 1D PlumeRise model is used. PlumeRise describes the ascent of a turbulent buoyant plume composed of solid pyroclasts, magmatic gases, and air and includes the effects of cross-wind on turbulent entrainment of air into the plume. The PlumeRise dynamics model approximates volcanic gases as entirely H 2 O (a more complex emission composition provides input to the chain-ofreactors model, see Section "Magmatic Gas Emissions" below).
Output from the PlumeRise model includes the plume temperature (T, in K), the total mass flux (Q in kg/s), the gas Frontiers in Earth Science | www.frontiersin.org mass fraction (n), the vertical velocity, U (m/s), the mass fraction of water (liquid+gas) in the plume (ϕ), specific humidity (C), and altitude (z in m), each provided at a resolution (dz) of a few meters. These variables can be used to derive temperature and VA:VM as a function of time. The time, in seconds is given by E4, and integrating dt.
E4 dt = dz/U The plume total mass (gases and solids), M in kg, in the discrete volume at each time-step is given by E5.
The mass of all gases (volcanic H 2 O and air) in the plume, M (g) is the product of M and the gas mass fraction, n, E6.
Under high temperature plume conditions, all H 2 O in the plume is in gaseous form, whose mass M (w) is the product of M and ϕ, E7. E7 The mass of dry air in the plume, M (da) is the difference between M (g) and M (w) , E8. The mass of background water in the plume, M (bw) is the product of M (da) with the specific humidity, C, E9. The mass of volcanic gases (here: entirely as H 2 O), M (vg) , is the difference between M (w) and M (bw) , E10.
The masses M (vg) , M (da) , and M (bw) in kg are converted into moles of volcanic gases, dry air, and background water, m (vg) , m (da) , m (bw) , by dividing by their respective molecular weights (18 × 10 −3 for H 2 O and 29 × 10 −3 for dry air containing N 2 , O 2 , and Ar, expressed as kg/mol), E11,12,13. The total moles of entrained wet air, m (air) , is calculated by the sum of m (da) and m (bw) , E14.
Finally, the ratio VA:VM is the moles of wet air divided by moles of volcanic gases, E15.
The calculation yields VA:VM as a function of dz or time when performed for each row of output from the 1D plume dynamics model. Assumptions that water exists purely in the gas phase and that water content is dominated by the volcanic emission are valid for the hot and concentrated plume region investigated in this study. The resulting trend in temperature with VA:VM is consistent with simple specific-heat calculations for the magmatic-air-ash mixture. Plume temperature and VA:VM depend strongly on the T, n, U starting conditions used for the PlumeRise model. Several PlumeRise 1D model results showing decreasing temperature and increasing VA:VM with time or distance from the source are shown in Supplementary Figure S1, for the Eyjafjallajökull eruption April 14-17, 2010. These are based on the model runs of Woodhouse et al. (2013) that used meteorological data from local balloon soundings 1 . The model input parameters are only partially constrained, but Woodhouse et al. (2013) found good predictive capability of the model to reproduce plume height for starting parameter ranges of temperature T = 600-1000 K, vertical velocity U = 50-100 m/s, gas mass fraction n = 0.030-0.086 (Woodhouse et al., 2013; Supplementary Table S3). Notably, the model runs are highly diverse in their starting temperatures and predicted VA:VM. The starting temperatures are all substantially lower than petrological estimates of the magmatic temperature, (Keiding and Sigmarsson, 2012) who reported 1443 K (1170 ± 25 • C) for the early flank eruption, and 1273-1333 K (1000-1060 • C) for benmoritic products related to the main summit eruption that are pertinent to this study. Reasons for this are unclear but could be due to decoupling of magmatic and gas temperatures or a cooling effect of incorporated water. The Introduction Section "High-Temperature Products in Volcanic Plumes: Observations and Modeling" highlighted considerable uncertainty in the physical conditions (temperature, VA:VM) of the near-source plume.
In this exploratory study we use output from the PlumeRise simulation with starting temperature at 1000 K as example boundary conditions in time-varying T and VA:VM as a basis for chemical kinetics modeling of the hot plume chemistry. This chosen PlumeRise dynamics scenario from Woodhouse et al. (2013) combined with magmatic gas emissions based on Allard et al. (2011; Table 2) yields an SO 2 flux of 5.5 kt/hr SO 2 , that is consistent with the range of SO 2 fluxes derived from satellite observations during the Eyjafjallajökull eruption (Boichu et al., 2013).

Magmatic Gas Emissions
A volcanic gas emission composition is input to the chain-ofreactors model, based on reported observations. It should be emphasized that available observations are an imperfect basis for the initialization of model of high-temperature chemistry in the near-source plume. First, no individual measurement technique can observe all gases in a volcanic plume. Available observations may not fully characterize all magmatic gases in the emission, nor any temporal evolution in their release as a function of volcanic activity. Second, observations of the cooled or partially cooled plume necessarily already incorporate any high-temperature modification of the emissions. Our use of such observations inherently assumes that preceding plume chemistry has only slightly modified the volcanic gas emission. This is indeed the prevalent assumption behind near-source monitoring of volcanic gas emissions (with the exception of volcanic BrO, e.g., Guttmann et al., 2018). Unfortunately, sources such as lava bodies are typically hazardous and inaccessible, therefore few direct measurements are available of the gas emission prior to any mixing with air.
In this study the composition of the Eyjafjallajökull magmatic gas emission is based on observations of H 2 O-CO 2 -SO 2 -CO-(HCl) in the young plume of the summit eruption reported by Allard et al. (2011) who performed FTIR spectroscopy of volcanic gases in early May 2010, using hot rocks as an IR source. Note that the volcanic emission of HCl was replaced by Ar in our modeling that does include halogens. The emission scenarios for the model studies are summarized in Table 2.
A first simulation emissions scenario (Eyja CO) is based entirely on observations of Allard et al. (2011). This emission contains CO but no other reduced gases such as H 2 or H 2 S, which typically cannot be resolved in this experimental set up. HF was also not reported, being below detection limits. The measurements by Allard et al. (2011) yield a CO/SO 2 molar ratio of 0.175 and CO/CO 2 = 0.0044. To evaluate how the presence of CO and SO 2 impact the high-temperature chemistry of the nearsource plume, sensitivity studies were undertaken that excluded CO and SO 2 from the emission (Eyja no-CO no-SO2).
A further simulation (Eyja CO-H 2 -H 2 S) was undertaken that included estimates for H 2 and H 2 S in the Eyjafjallajökull emission (based on thermodynamic equilibrium calculations) alongside CO and the other gases reported by Allard et al. (2011). It is assumed that the magmatic gases were in thermodynamic equilibrium and at magmatic temperature on emission. For the Eyjafjallajökull summit eruption a temperature of 1273 K is assumed (Keiding and Sigmarsson, 2012). The gases obey the equilibria CO + 0.5 O 2 = CO 2 , H 2 + 0.5 O 2 = H 2 O, H 2 S + 1.5 O 2 = SO 2 + H 2 O with corresponding thermodynamic equilibrium constants of 1.13 × 10 7 , 1.92 × 10 7 , and 1.54 × 10 17 at 1273 K (constants obtained from HSC). For these conditions the CO/CO 2 molar ratio 0.0044 observed by Allard et al. (2011) corresponds to an oxygen fugacity of 4.0 × 10 −10 bar (i.e., QFM+1.6 log units). The corresponding H 2 O/H 2 molar ratio is 387 and SO 2 /H 2 S molar ratio is 1379. Combining these theoretical molar ratios with the mole fraction composition for H 2 O and SO 2 in Table 2 yields estimates for the mole fractions of H 2 = 0.0024 and H 2 S = 1.45 × 10 −6 in the Eyjafjallajökull emission. Thus, it is predicted that the eruption emitted a similar quantity (by mole) of H 2 compared to SO 2 , but only a trace emission of H 2 S. Some volcanoes release a larger proportion of sulfur in reduced forms. To investigate the potential impact of H 2 S on the high-temperature chemistry of near-source volcanic plumes, we undertook a sensitivity study using a composition based on measurements made at Mt Etna (passive degassing) whose emission composition is still relatively oxidized, but includes H 2 S (Roberts et al., 2017 and references therein) and H 2 (Aiuppa et al., 2011;Supplementary Table S4). We emphasize that our combination of Eyjafjallajökull plume dynamics with Mt Etna emission composition for this sensitivity study is not representative of plume chemistry for Mt Etna passive degassing conditions. Rather it reflects a possible scenario for an ash-rich explosion from Mt Etna, assuming emission composition similar to that observed during passive degassing, and plume dynamics similar to Eyjafjallajökull. The Mt Etna composition used as input to our kinetics-model study is identical to the input used for a thermodynamic model study by  [where the H 2 emission was calculated theoretically and is somewhat higher than that observed by Aiuppa et al. (2011)]. Notably, the SO 2 content of the Mt Etna emission composition (0.029 mole fraction) is more than an order of magnitude higher than the SO 2 content of Eyjafjallajökull emission (0.002 mole fraction). The H 2 S content of the emission is also much higher for Mt Etna (also with a higher ratio to SO 2 ), whilst the H 2 , H 2 O, and CO 2 mole fractions are similar (although ratios to SO 2 , X/SO 2 , are lower). A sensitivity study was also performed without H 2 S in the Etna emission (Etna no-H 2 S).

Equilibrium Modeling: HSC Thermodynamic Model
Our study compares output from the chemical kinetics model outlined above to an equilibrium model (HSC) that has been previously applied to investigate volcanic plume chemistry (e.g., Martin et al., 2006). HSC model version 7.1 was used to calculate equilibrium compositions of mixtures of air with magmatic gases under the same temperature and VA:VM boundary conditions as for the chain-of-reactors model. The input composition was determined by mixing air (composed of 78% N 2 , 21% O 2 and 1% Ar) with the Eyjafjallajökull magmatic gas composition ( Table 2) for each VA:VM, whilst temperature was specified according to the PlumeRise boundary conditions. Pure magmatic gases are usually assumed to have compositions at chemical equilibrium when they are released from the volcano (and there is good evidence for this to first order for major Cand S-bearing species, e.g., Oppenheimer et al., 2018); however, some repartitioning of Eyjafjallajökull gases may be expected for thermodynamic equilibrium calculations even at VA:VM = 0:100 because the reported gas emission composition is incomplete and the initial PlumeRise temperature is slightly lower than magmatic temperature (see section "Magmatic Gas Emissions").

Model Diagnostics: Cooling and Mixing of Volcanic Gases With Air
Physical parameters from the PlumeRise model (temperature T, VA:VM, distance from the source, and plume cross-sectional area) used as boundary conditions for the chain-of-reactors model vary as a function of time, Figures 2A,B. Temperature is varied smoothly across the chain-of-reactors and decreases from 1000 K to 890 K (over the 30 PFRs) whereas VA:VM increases stepwise from 0 to 5.3 due to air-entrainment between each reactor. This T and VA:VM scenario is not representative of all volcanoes. Smaller volcanic emissions would likely disperse faster, particularly if ash-poor. The presence of ash in Eyjafjallajökull plume maintains high temperatures at higher VA:VM than would occur in an ash-poor emission (that for example for VA:VM > 3 would yield T < 500 K according to simple gasphase mixing). This result is not surprising given 1D plumerise models of explosive eruptions typically assume an emission containing volcanic gases at a few weight-percent, with the remainder as pyroclasts. We highlight that this "temperature buffering" effect of ash may significantly prolong the period in which volcanic gases undergo high-temperature chemistry as they mix with air.
As VA:VM increases, the abundance of air-dominated gases such as O 2 increases step-wise toward an asymptote of 21% (pure air) whilst major magmatic gases such as H 2 O and CO 2 decrease step-wise due to dilution, Figure 2C. Plume chemistry can additionally cause gas abundances to rise or fall within each reactor. Figure 2C shows that there is negligible or only small impact of plume chemistry on O 2 , H 2 O, and CO 2 whose abundances are largely controlled by air entrainment/dilution. Also shown in Figure 2C is the product of the O 2 and H 2 O mole fractions that reaches a maximum at around 0.5 s, of relevance to the high-temperature reactions discussed in Section "Discussion."

High-Temperature Production of Oxidants and Radicals
The chemical evolution of the Eyjafjallajökull near-source plume according to the high-temperature chain-of-reactors model is shown in Figure 3A (Eyja CO run) over a time-period of about 2.4 s, during which VA:VM increases from 0 to 5.3, and temperature decreases from 1000 K to 890 K. The slight declines in major gases H 2 O, CO 2 , and SO 2 mixing ratios show they are largely under the control of dilution whilst air is entrained, i.e., are largely conserved. The production of oxidant radicals such as OH, HO 2 , and H 2 O 2 increases rapidly with plume time evolution during the first 0.5 s, after which their absolute mixing ratios decrease. Formation of SO 3 follows a similar pattern. The volcanic gas CO decreases only slightly faster than the major gases, indicating that it is not substantially oxidized by the plume chemistry. The general pattern is similar for the Eyja CO-H 2 -H 2 S simulation, Figure 3B: a higher maximum abundance of oxidant radicals is reached (slightly earlier) and slightly more CO is oxidized, but still not completely depleted. In the Eyja CO-H 2 -H 2 S simulation both volcanic gases H 2 and H 2 S become oxidized at around 0.5 s. To analyze the chemical kinetics model results in more detail and compare to field-observations and thermodynamic calculations, ratios are presented of the minor volcanic gases relative to the major (largely conserved) gases SO 2 , CO 2 , and H 2 O.

Oxidation of CO Is Kinetics Limited in the Near-Source Plume
CO is only slightly oxidized by high-temperature chemistry, as shown by CO/SO 2 ratio, Figure 4A (as well as CO/CO 2 , Figure 4B). At the end of the simulations (VA:VM = 5.3 after 2.4 s), CO/SO 2 has only declined by about 15% of the initial ratios for Eyja CO, and by 37% of the initial ratios for Eyja CO-H 2 -H 2 S. In contrast, HSC thermodynamic calculations predict near-complete oxidation of CO as soon as air is mixed with the Frontiers in Earth Science | www.frontiersin.org 8 July 2019 | Volume 7 | Article 154 and indeed in other volcanic gas emissions from more reduced magmas, e.g., at Mt Erebus (Oppenheimer and Kyle, 2008). Additional observations of volcanic CO in the Eyjafjallajökull plume are reported by Schumann et al. (2011) who made instrumented aircraft measurements of the aged (10-100 h) plume over Europe during mid-April to mid-May. Both CO and SO 2 were observed above background levels. The molar ratio of excess CO and excess SO 2 (as averages) in transects of the aged plume is around ∼1, i.e., somewhat higher than the ratio reported by Allard et al. (2011) in the young plume. Reasons for this discrepancy are not clear. It is unlikely that CO was formed in the plume in an oxidizing atmosphere. Plume structure in the aircraft measurements as reported by Schumann et al. (2011) is smoother for SO 2 compared to CO. This may in part be due to the vastly different response times of the aircraft-based instruments for SO 2 (T95 response time to reach 95% signal is 80 s) and CO (response time < 2 s), which might cause a high bias in the derived CO/SO 2 . There is some scatter in the CO and SO 2 measurements as well as variability in background CO. In addition, SO 2 could undergo atmospheric oxidation processes, resulting in an increase in the CO/SO 2 ratio in the downwind plume. Unusually, at Erebus volcano a rapid loss of SO 2 has been identified in the very young plume possibly aided by cloud processing (Oppenheimer et al., 2010). However, oxidation of SO 2 in volcanic plumes is typically slow (e.g., Galeazzo et al., 2018). The measurements of CO/SO 2 in the downwind Eyjafjallajökull plume do not show any clear dependency on plume age. In any case, the aircraft measurements of Schumann et al. (2011) detected volcanic plume CO alongside SO 2 , confirming that CO was emitted during the 2010 Eyjafjallajökull eruption and that CO was not substantially oxidized by high-temperature nearsource chemistry as is commonly predicted by thermodynamic modeling. Oxidation of CO was evidently kinetics-limited in the transient near-source plume. A summary of the model and observed CO/SO 2 is given in Table 3.
The Eyja CO-H 2 -H 2 S simulation that includes theoretical estimates for volcanic emission of H 2 and H 2 S predicts nearcomplete oxidation of both gases, Figure 4B, as well as partial oxidation of CO. This is due to an oxidant enhancing feedback mechanism (see section "Discussion"). No specific observations of H 2 or H 2 S are available for comparison in the Eyjafjallajökull near-source or far downwind plume. However, observations of H 2 S in plumes from both passively degassing volcanoes (e.g., Roberts et al., 2017) and explosive eruption plumes (e.g., de Moor et al., 2016) indicate that H 2 S is (largely) conserved, as for CO. Measurements of H 2 have been reported in plumes from some passively degassing volcanoes (e.g., Aiuppa et al., 2011;Moussallam et al., 2012). In general, observations of reduced gases in volcanic plumes (specifically CO/CO 2 , H 2 O/H 2 , and SO 2 /H 2 S gas ratios) may be used to infer magmatic properties at emission such as oxygen fugacity (see section "Magmatic Gas Emissions"), on the condition that modification of the emitted composition by near-source plume chemistry is either negligible or well-constrained by plume chemistry models. Here the small decrease in CO/CO 2 predicted by the kinetics model (Eyja CO) is equivalent to a difference in oxygen fugacity of 1.5 × 10 −10 bars or about 0.14 log unit QFM. We highlight two important areas of model uncertainty regarding the near-source oxidation of reduced gases. First, the reactions are temperature-dependent yet the temperature of the volcanic gas emission as well as subsequent plume physical evolution are both rather poorly known. Second, the kinetics of certain reactions, particularly for H 2 S are poorly constrained by laboratory experiments. Thus, the presence of H 2 and H 2 S observed in certain volcanic plumes may suggest either a lower emission temperature or faster cooling rate than the Eyjafjallajökull case presented here, or may reflect uncertainties in the model chemical schemes.  High-Temperature Production of SO 3 as a Sulfate Aerosol Precursor SO 3 abundance in magmatic gas is negligible but is predicted to form through high-temperature plume chemistry. SO 3 is a precursor to sulfate aerosol (as is well-known e.g., from studies of power station emissions). Upon further plume cooling, volcanic SO 3 will react readily with H 2 O (g) to form H 2 SO 4(g) that is highly hygroscopic and will form sulfate particles. The widespread observation of near-source volcanic sulfate (see section "Introduction, " Table 1), motivates our model investigation of high-temperature SO 3 as a possible sulfate source. The SO 3 abundance increases rapidly at around 0.5 s, Figure 5A. The SO 3 /SO 2 molar ratio increases throughout the model simulations, indicating a continual production of SO 3 that is similar for Eyja CO and Eyja CO-H 2 -H 2 S. By the end of the simulations (2.4 s, VA:VM = 5.3, T = 890 K), the SO 3 /SO 2 molar ratio reaches 10 −3 . The fact that SO 3 /SO 2 ratio is still increasing at the end of the model simulation indicates that additional  high-temperature production of SO 3 is expected as the plume cools and dilutes further. Extrapolation of Figure 5A suggests SO 3 /SO 2 may potentially reach up to ∼2 × 10 −3 . Whilst there are no reported sulfate/SO 2 measurements for the Eyjafjallajökull summit eruption for specific comparison, this modeled ratio of SO 3 /SO 2 of around 10 −3 lies within the range of observed near-source sulfate/SO 2 ratios, including the Eyjafjallajökull flank eruption ( Figure 5B and Table 1). In contrast, calculations using the HSC thermodynamic equilibrium model at VA:VM = 5.3, T = 890 K predict SO 3 /SO 2 of around 3.1, an order of magnitude above the highest ratio of sulfate/SO 2 in the field observations, and several orders of magnitude higher than the Eyjafjallajökull observations. We conclude that the chemistry producing SO 3 in the near-source volcanic plume is under kinetic rather than thermodynamic control.

High-Temperature Production of Oxidants HO 2 , OH, and H 2 O 2
High-temperature oxidants HO xy (HO 2 , OH, and H 2 O 2 ) are predicted to form in the Eyjafjallajökull plume through high-temperature plume chemistry, Figures 6A,B (absolute concentrations shown). Also shown for scale is the SO 2 concentration that decreases largely due to dilution. Hightemperature HO xy reaches a maximum around 0.5 s. HO xy production is greater for Eyja CO-H 2 -H 2 S than for the Eyja CO simulation. In both cases HO x exists predominantly as HO 2 . This contrasts with thermodynamic calculations that predict HO x to exit predominantly as OH (e.g., Gerlach, 2004). Another important feature of the chemical kinetics modeling is that H 2 O 2 is predicted to form (at greatest abundance for Eyja CO-H 2 -H 2 S), whilst H 2 O 2 is essentially not formed in HSC thermodynamic calculations (e.g., Martin et al., 2006). High-temperature SO 3 was discussed above as a precursor to sulfate aerosols. The H 2 O 2 formed by high-temperature chemistry may promote further sulfate formation upon plume cooling, due to its reaction with dissolved SO 2(aq) . If a condensed phase is formed already in the crater zone this could yield an additional source of "primary" or near-source sulfate. Alternatively, if plumes become condensed clouds further  Frontiers in Earth Science | www.frontiersin.org FIGURE 7 | (A) schematic of key high-to-low temperature processes in the volcanic plume, including high-temperature production of SO 3 as sulfate aerosol precursor, HOx oxidants and halogen radicals, and low-temperature multi-phase chemistry producing BrO. (B) Calculated lifetimes of HCl and HBr in the plume with respect to their reaction with high-temperature OH as a function of plume evolution time. Solid lines: Eyja CO simulation, Dashed lines: Eyja CO-H 2 -H 2 S simulation. Lifetime minima are <1 s, see Table 4. downwind, the volcanic H 2 O 2 may form secondary sulfates. In this respect, Carn et al. (2011) measured H 2 O 2 in the plume emitted by Huila volcano (elevated above background by about 0.5 ppbv) that was rapidly titrated when the plume entered a cloudy phase. We are unaware of other observations of H 2 O 2 in volcanic plumes.
We highlight two further potential impacts of volcanic H 2 O 2 . Studies have shown that iron on the surface of volcanic ash has the potential to generate substantial quantities of hydroxyl radicals that may affect the toxicity of inhaled ash particles (Horwell et al., 2003;Hillman et al., 2012). This occurs via the Fenton reaction involving H 2 O 2 . The elevated abundance of H 2 O 2 in the volcanic plume predicted by our model study would likely enhance the production of OH radicals and their potential health hazard. Indeed, an observation that has puzzled the volcanology-atmospheric community to date is the presence of elevated OH (and HO 2 ) measured in the 2010 Mt Hekla eruption cloud in the stratosphere (Rose et al., 2006). This elevated OH is all the more surprising given the measurement was made during night-time and in the presence of major OH reactive sinks (e.g., HCl and SO 2 ). We propose the reaction of volcanic H 2 O 2 on ash particles as a possible mechanism that could generate OH radicals in volcanic plumes even during night-time.

Potential for High-Temperature Formation of Halogen Radicals
High-temperature oxidants such as OH and radicals such as Br, Cl are believed to have an important role in accelerating the onset of low-temperature atmospheric chemistry "bromine explosion" cycles that form BrO and destroy ozone in the downwind plume (e.g., Bobrowski et al., 2007;Roberts et al., 2014;Surl et al., 2015; Figure 7A). The bromine explosion is autocatalytic so can be efficiently "kick-started" by a small quantity of high-temperature halogen radicals. Observations show that halogens were emitted during the eruption (young plume HCl/SO 2 ≈ 0.65, Allard et al., 2011) and that reactive halogen chemistry was active in the plume (downwind plume BrO/SO 2 ≈ 1.3 × 10 −4 , Heue et al., 2011). The Eyjafjallajökull volcanic halogen chemistry caused a depletion of tropospheric ozone (Vance et al., 2010;Schumann et al., 2011). Here offline calculations show the potential for formation of high-temperature halogen radicals to initiate this process.
The OH formed in the near-source volcanic plume is a powerful oxidant that can react with volcanic halogens (emitted as HCl and HBr) to produce halogen radicals (R1 and R2). The high-temperature chemical kinetics model scheme used in this study does not include halogen chemistry. Nevertheless, offline calculations of the lifetime of volcanic HCl, HBr emissions can be undertaken to estimate the potential for halogen radical formation, shown here for the Eyjafjallajökull simulations.  Góger et al. (2018) as k = (9.86 ± 2.38) × 10 −16 × T (1.23 ± 0.03) × exp[(5.93 ± 0.33) kJ mol −1 /RT] cm 3 molecule −1 s −1 (valid for 600-3200 K).

R1 HCl
Lifetimes of HBr and HCl in the near-source plume are calculated by the inverse of the first-order rate constant, i.e., the product of the 2nd order rate constants for R1 and R2 with the OH concentration. The temperature dependence of the rate constants is only slight (variation of up to about 15%) over our studied temperature range 890-1000 K. Rather, Figure 7B shows that the calculated HCl and HBr lifetimes during the Eyja 4 | Simulated maximum OH/SO 2 ratio and OH concentration in the high-temperature near-source plume according to the chain-of-reactors model, and minimum lifetime of HX (X = Br or Cl) with respect to reaction with OH by offline calculation (see section "Potential for High-Temperature Formation of Halogen Radicals" and Figure 7B). For comparison, the observed Eyjafjallajökull plume HCl/SO 2 and BrO/SO 2 molar ratios are 0.6 and ≈1.3 × 10 −4 , respectively, and the total time-duration of the high-temperature simulation is 2.4 s.
CO and Eyja CO-H 2 -H 2 S simulations are a strong function of OH concentration. The halogen halide lifetimes rapidly decrease as OH is formed in the near-source plume, and then increase as OH is consumed (and diluted). Lifetimes of HBr and HCl reach minima of <1 s in the plume, Table 4. Lifetimes are shortest for HBr that has the greater rate constant of reaction and for the Eyja CO-H 2 -H 2 S simulation that produces highest OH concentration. In all cases, the calculated HBr and HCl lifetime minima are shorter than the timescale of the plume simulation (2.4 s) indicating that significant halogen radical production will occur. Reactions R1 and R2 will likely cause a negative feedback on OH concentrations. Once formed the halogen radicals will react with CHOS and some may reform HBr or HCl. A more comprehensive understanding of the halogen interactions will require inclusion of halogen reactions explicitly in the hightemperature model chemical mechanism.

Chemical Reaction Mechanisms and Feedbacks
The chemical kinetics model identifies key reactions causing oxidation and oxidant production in the volcanic plume. In this study, the model simulations predict abundances of oxidants, and oxidized species increase rapidly around 0.5 s, at In high-temperature chemistry (generally temperatures above 800-850 K), formation of H 2 O 2 is usually considered as a sink of oxidant radicals. However, in the case of a near-source plume the removal of OH by CO (R5, followed by R6 and R7 leading to H 2 O 2 ) drives the formation of more OH and HO 2 from the reaction of H 2 O with O 2 (more specifically, the balance of this reaction, R3, with its reverse reaction). This feedback mechanism promotes the production of HO xy in volcanic plumes containing CO, where only a fraction of the emitted CO needs to be oxidized to substantially enhance HO xy abundances and promote the chain reaction. Similarly, the H 2 emission in the Eyja CO-H 2 -H 2 S simulation can also react with OH, R9 (followed by R6 and R7) to cause a similar positive feedback on HO x production.

R9 H 2 + OH H 2 O + H
These feedbacks involving reduced gases may explain the greater production of oxidants in Eyja CO-H 2 -H 2 S than Eyja CO simulation. This hypothesis is supported by a sensitivity study without CO and SO 2 in the emission (Eyja no-CO no-SO 2 , Supplementary Figure S2) that produces a lower OH abundance but where more OH accumulates more in the first hundreds of milliseconds, reflecting the balance of OH production and consumption.
We now consider the role of H 2 S. The very low (theoretically estimated) H 2 S emission from Eyjafjallajökull had a rather negligible impact on the near-source plume chemistry. The role of H 2 S is investigated in a simulation of the high-temperature chemistry of Mt Etna emissions, that include H 2 and H 2 S as well as SO 2 , H 2 O, and CO 2 , see Supplementary Table S4 (as well as Frontiers in Earth Science | www.frontiersin.org a control run that excludes H 2 S). The Etna simulation assumes an eruption with the same plume physics as for Eyjafjallajökull. A very rapid production of oxidants is predicted within the first 0.1 s of plume evolution, Supplementary Figure S3, coincident with a rapid decrease in H 2 S (and a slower decrease in H 2 ). Abundances of oxidant radicals initially substantially exceed those predicted for Eyja CO and Eyja CO-H 2 -H 2 S, and later decline to lower levels. According to our model, the oxidation of emitted H 2 S causes an extra production of O atoms (which is not observed in the case "without H 2 S") with in the first milliseconds of the plume's mixing with air, following the sequence: H 2 S + H (or OH) = SH + H 2 (or H 2 O), 2SH = H 2 S + S, S + O 2 = SO + O, and SO + O 2 = SO 2 + O. The production of sulfur atoms from H 2 S and their subsequent reactions with O 2 are mostly responsible for this production of O, which, in turn, converts SO 2 into SO 3 , resulting in the SO 3 peak observed at 0.1 s in the plume. Beyond this point, the conversion of SO 2 into SO 3 is mainly due to the abundance of HO 2 through SO 2 + HO 2 = SO 3 + OH. A significant part of SO 3 is then used to recycle SO 2 through the formation of HOSO 2 (SO 3 +HO 2 = HOSO 2 +O 2 and HOSO 2 +M = SO 2 +OH+M). Indeed, the formation of SO 3 in the Mt Etna scenario proceeds through two different mechanisms (with O, and then with HO 2 ) at different locations in the plume whereas in the Eyja CO and Eyja CO-H 2 -H 2 S or Etna no-H 2 S scenarios, the formation of SO 3 is only linked to HO 2 . It must be emphasized, however, that model predictions are limited by uncertainty in high-temperature H 2 S chemistry.
This sulfur cycling regulates SO 3 /SO 2 to reach around 10 −3 after 2.4 s in all of the simulations (Eyja CO, Eyja CO-H 2 -H 2 S, Etna and Etna no-H 2 S, Figure 5A and Supplementary Figure S3), despite their differences in CO-H 2 -H 2 S emissions that lead to vast differences in the production and abundances of HO xy . Temperature and VA:VM evolution with time may be important controls on SO 3 /SO 2 , therefore production of near-source sulfate could vary for different physical conditions other than the Eyjafjallajökull case considered here.

Model Uncertainties and Future Directions
Our chain-of-reactors model provides a framework for quantitative studies into the high-temperature reactions occurring in the near-source plume to form oxidized products such as sulfate-precursors, OH, H 2 O 2 that can impact subsequent physico-chemical processing in the cooled plume. We find that the high-temperature chemistry is under kinetic control; plume composition does not follow thermodynamic equilibrium (Figures 4A, 5B, with example model output in Supplementary  Tables S5, S6). Our exploratory study identifies several areas for further investigation.

Improved Model Chemistry at High-Temperatures and Link to Low-Temperature Plume Processes
Uncertainties in the high-temperature gas chemistry mechanism need to be addressed such as the chemical kinetics of sulfur species reactions, notably the oxidation of H 2 S, that are poorly quantified from laboratory experiments. Future model developments should also incorporate the chemistry of halogens and nitrogen species coupled to CHOS. Simulations of longer duration should cover the whole regime where high-temperature chemical reactions may be important (e.g., down to 800/850 K). The gas-phase mechanism should be extended to consider photolysis reactions for plumes that are not too optically thick, and, especially, heterogeneous chemical reactions, for example reactions on ash surfaces (Delmelle et al., 2018), noting the oxidized products from high-temperature gas chemistry could act to modify ash properties. Models should ultimately be developed that encompass both high-temperature chemistry and also the transition to low-temperature plume chemistry as simulated by atmospheric models (e.g., Jourdain et al., 2016).

Influence of Emitted Magmatic Gas Composition on High-Temperature Chemistry
The composition of magmatic gases released is specific to each volcano and can vary with volcanic activity/unrest. The chemical feedback mechanisms identified by our model suggest that emissions of reduced gases can enhance the high-temperature production of oxidants. This might lead to a greater production of oxidant radicals in the near-source plume of volcanoes with a relatively reduced emissions composition, such as Mt Erebus, Antarctica, whose emissions are rich in CO and H 2 . Few field-measurements exist to evaluate model predictions: the high-temperature region is difficult to access directly and few studies have quantified oxidized products in the cooled downwind plume, e.g., by instrumented aircraft. Another critical uncertainty is that the emission temperature of volcanic gases is poorly known. In some situations, according to the dynamics of the degassing, magmatic gas emissions at the surface may be neither in chemical nor thermodynamic equilibrium with the melt (Oppenheimer et al., 2018).

Influence of Plume Dynamics on High-Temperature Chemistry
The chemical kinetics model requires time-resolved boundary conditions in plume dynamics, notably temperature and VA:VM (degree of mixing of air with magmatic gases). This study used output from a 1D PlumeRise model for the Eyjafjallajökull 2010 eruption (Woodhouse et al., 2013) for this purpose. Variability in plume dynamics was not investigated but is expected to exert a strong influence on the near-source chemistry through the rate of entrainment of air and cooling. For example, smaller emission sources will typically disperse and cool faster in the troposphere. Ash-poor emissions will tend to cool more rapidly than ash-rich emissions where the solid pyroclasts buffer the temperature of the volcanic gas and air mixture. Large eruptions that inject gases and ash into the stratosphere should also be investigated. Spatial heterogeneity in the turbulent mixing of air into the plume might nonlinearly impact the plume chemistry, potentially requiring the development of 3D model approaches to investigate the plume's high-temperature chemistry as a function of a more detailed turbulent mixing.

CONCLUSIONS
We present chemical kinetics model simulations of the hightemperature gas chemistry that occurs as magmatic gases mix with air and cool in the near-source volcanic plume. Volcanic gases released into the atmosphere undergo rapid processing including photochemistry, oxidation, and aerosol formation. Quantifying these processes in volcanic plumes is essential to assess atmospheric and environmental impacts of volcanic emissions. During the first few seconds after release, the hot magmatic gases mix with background air containing oxidants and undergo a high-temperature processing. The gas chemistry in this high-temperature region of the plume is poorly characterized, yet can critically modify the composition of the emission to produce oxidized products and radicals, such as near-source sulfates, HOxy, and halogen radicals. These high-temperature products influence the subsequent physico-chemical processing of the plume at low-temperatures as it disperses in the atmosphere. Studies to date have largely represented the chemical composition of the high-temperature region using thermal equilibrium calculations (e.g., Martin et al., 2006;Hoshyaripour et al., 2012). A source of uncertainty in thermodynamic modeling of the hightemperature chemistry of volcanic plumes is the assumption of equilibrium conditions, that may not always be valid for transient plumes .
The chain-of-reactors approach of this study is timeresolved and based on chemical kinetics thus differs markedly from thermodynamic equilibrium calculations used to date. We performed an exploratory study of the high-temperature chemistry in a near-source volcanic plume over the first 2.4 s, using as boundary conditions for the time-varying temperature and mixing of air output from 1D PlumeRise model simulations of Eyjafjallajökull plume (Woodhouse et al., 2013). The magmatic gas emission includes H 2 O, CO 2 , SO 2 , and CO, based on observations of the young plume (Allard et al., 2011), with theoretical estimates for H 2 S and H 2 gas emissions added in a sensitivity study and offline calculations performed for halogens.
A key finding is that the plume chemistry is under kinetic control and that thermodynamic equilibrium cannot be assumed. The simulation predicts that CO emission from Eyjafjallajökull is not fully oxidized (largely conserved) by high-temperature chemistry. This is consistent with field-observations of CO in the volcanic plume (Allard et al., 2011;Schumann et al., 2011) and contrasts to thermal equilibrium calculations of near-complete CO oxidation, Table 3. Instead, our study shows CO oxidation is under chemical kinetics control. The model predicts that sulfateprecursor SO 3 is formed at levels (SO 3 /SO 2 = 10 −3 ) in general agreement to observations of near-source volcanic sulfate/SO 2 in volcanic plumes, Table 1. Conversely, thermodynamic calculations predict a much greater formation of SO 3 relative to SO 2 for the same plume dynamics conditions. The kinetics model predicts high-temperature formation of HO 2 , OH, and H 2 O 2. Formation of HO xy is predicted to be enhanced when additional reduced gases (H 2 S and H 2 ) besides CO are included in the emission.
Offline calculations demonstrate that volcanic halogens (HBr and HCl) will react rapidly with high-temperature OH, highlighting the potential for halogen radicals to form in the high-temperature plume. A key chemical mechanism is that HO x oxidants are formed by the (reversible) reaction of H 2 O with O 2 at high-temperatures to produce OH and HO 2 . Chemical feedbacks occur when reduced gases in the volcanic emission (CO and H 2 ) react with OH to produce HO 2 that can selfreact to form H 2 O 2 . Thus reduced gases in the volcanic emission act to promote formation of high-temperature HO xy . Another chemical feedback occurs via the oxidation pathway for volcanic H 2 S that rapidly generates additional oxidants in a chainreaction. However, there is large uncertainty in some hightemperature reactions, particularly for sulfur species. Quantifying the chemical processes occurring at high-temperatures in the near-source plume is an essential step in understanding the reactivity and evolution of volcanic emissions and their impacts in the troposphere and stratosphere. Improved characterization of this high-temperature region of the plume will require the combination of model, field-observation and laboratory studies.

AUTHOR CONTRIBUTIONS
All authors contributed to the study design, analysis and interpretation of data, and manuscript writing.