One-atmosphere high-temperature CO–CO 2 –SO 2 gas-mixing furnace: design, operation, and applications

. In this paper, we present a new design for a 1 atm gas-mixing furnace using the gas mixture CO– CO 2 –SO 2 . This furnace can simulate disequilibrium processes such as magmatic and volcanic degassing. Here, we present the technical aspects of the design. The furnace can sustain temperatures of up to 1650 ◦ C and has a hot zone that spans 200 mm vertically, where the hotspot is determined to be ∼ 32 mm below the midpoint of the furnace enclosure. The four mass ﬂow controllers are individually calibrated and accurate to within 0.8 % of the speciﬁed value. The f O 2 is accurately reproduced in the furnace within ± 0 . 002 log units, as calibrated by the Fe–FeO reaction across the iron–wüstite (IW) buffer at 1300 ◦ C. The furnace can reliably simulate dynamic conditions, where the f O 2 can be modulated at a maximum rate of 2.0 log units min − 1 by varying the gas mixture. A delay of 40 s is observed to attain the f O 2 calculated from the gas mixture, at the hotspot. A series of safety measures to protect the user from exposure to the toxic gases are detailed. In our experiments, the furnace is used to determine sulfur isotope fractionation factors among melt, sulﬁde, and the gas phase, within a magmatic context, using either crystals of olivine or silica glass tubes. The furnace has the potential to investigate various other dynamic high-temperature reactions occurring on Earth.


Introduction
The first person who conducted a petrology experiment using a furnace was Sir James Hall (1806). Since then, numerous high-temperature (HT) petrology studies (Roeder and Emslie, 1970;Lofgren, 1980;Rutherford, 1993;Hammer and Rutherford, 2002) have been conducted to determine the thermodynamic and kinetic parameters of various natural petrologic processes. Today, a vertical HT gas-mixing furnace is an essential instrument in an experimental petrology laboratory. It provides the ability to perform experiments at a high temperature (more than 2200 • C in some cases) within a controlled atmosphere. Several designs (i.e. temperature range, size, gas mixture) are available depending on the application. In general, these apparatuses are designed to maintain a steady temperature and internal atmosphere, at 1 atm, so that chemical and/or isotopic equilibrium can be attained between the condensed phases (liquid or solid) and the surrounding gas. In this contribution, we present a newly designed vertical 1 atm HT furnace, capable of dynamically changing the partial pressure of mixed gas species in the CO-CO 2 -SO 2 system during the experiment. Thermodynamic equilibrium between the gas and melt may not be obtained during degassing in the conduit prior to eruption. This is a process that can be simulated with the new furnace configuration (Pichavant et al., 2013;Le Gall and Pichavant, 2016;Moussallam et al., 2022).
The important characteristics of a gas-mixing furnace are the gas mixture, the size of the hotspot, and the quench mechanism. Here, the gas mixture of CO-CO 2 -SO 2 is well suited to studying sulfur isotope fractionation among melt, sulfide, 322 S. Prabha-Mohan et al.: One-atmosphere high-temperature CO-CO 2 -SO 2 gas-mixing furnace and gas. Given the strict safety rules that need to be followed, the design has eliminated any potential exposure to toxic gases and high temperatures. To our knowledge, no studies have reported gas-mixing furnaces that perform systematic changes in the gas mixture during an experimental run. Additionally, few designs use SO 2 in the gas mixture (Haughton et al., 1974;Katsura and Nagashima, 1974;Gaetani and Grove, 1997;Brenan and Caciagli, 2000;O'Neill and Mavrogenes, 2002). Experimental configurations using SO 2 are rare because of the increased toxicity and the potential for S precipitation during cooling to block the extraction pipes. The design by O'Neill and Mavrogenes (2002) is one of the few furnace designs to use a CO-CO 2 -SO 2 gas mixture. However, their design cannot be programmed to automatically change the composition of gases during an experiment. Also, their samples are quenched in water, whereas our design uses a cold hollowed-out stainless-steel bolt to quench the sample. While Brenan and Caciagli (2000) used a similar gas-mixing design, the sample was quenched by manually lowering it into an ice-water mixture, a process which is both slow and unsafe. However, the most crucial factors regarding the unique furnace design which is proposed here are the dynamically controlled atmosphere and a large volume. A much larger volume and hotspot allow for larger samples and a variety of sample holders. Increased safety measures also allow for longer experiments with minimal risk of accidental exposure.

Design and setup
The description of this 1 atm HT gas-mixing furnace is in three parts: (1) the furnace and sample holder, (2) the gasmixing unit, and (3) the safety systems.

Furnace and sample holders
The experimental setup is constructed around a 660 mm tall vertical tube furnace (Elite TSV 17/50/200), with a hot zone rated to be 200 mm vertically and with a maximum temperature of 1700 • C (sustained maximum temperature of 1650 • C). The six heating elements, Kanthal Super 1800, are set in a hexagon around the 60 mm outer diameter (50 mm inner diameter) and 840 mm long alumina tube. Heating is controlled with a Eurotherm 3216CP temperature controller. The top end of the alumina tube is cooled by a water circuit, assuring the electronic components are not exposed to heat and protecting the O-rings which seal the inside of the furnace from the atmosphere.
Our instrument can use one of two ceramic sample holders to set the sample in the hotspot of the furnace. One of them incorporates an oxygen fugacity sensor (ECONOX ContrOx 650) which can be used to calibrate the oxygen fugacity of the system. Both sample holders incorporate a thermocouple (Type B, Pt -30 % RH vs. Pt -6 % RH) and a pair of 1 mm thick platinum wires, used as quench electrodes (Fig. 1). The sample holder is introduced from the top of the furnace, and several O-rings prevent gas leaks and seal the inside of the furnace. The sample capsule (olivine crystal or silica tube) is suspended at the end of the sample holder, from the quench electrodes, using a thin platinum wire (diameter: 0.2 mm); this completes the quench current circuit. The sample is quenched by passing an electrical current of 20 A with a maximum DC voltage of 20 V, causing the thin platinum wire to melt, resulting in a drop of the sample to the cold end of the furnace (quench bowl). Typically, f O 2 is calibrated prior to the actual experiment using the sample holder with the oxygen sensor, whereas the simpler sample holder is used for the experiment. This is to safeguard and protect the f O 2 sensor from corrosion (potential reaction with S) during long experiments (6-24 h) as the ceramic sample holder is not perfectly sealed. The calibration also verifies the proper functioning of the mass flow controllers (MFCs) prior to each experiment.
Olivine crystals are used as sample capsules for experiments starting with basalt powder. We use high-Mg olivines (∼ Fo 90 , purchased from a seller in Arizona who sourced the olivine from peridotites in Pakistan), as they react minimally with basalt and impose a crystal-melt equilibrium which is appropriate for the natural systems investigated. The capsules are machined from whole crystals, rounded, and measure 5-10 mm. Silica glass tubes are used as sample capsules for experiments starting with other materials, such as rhyolite, and metals (Ni, Fe).

Gas-mixing unit/control interface LabVIEW
The gas-mixing unit was designed to ensure the safety of the operator; therefore, the gas flow is controlled by a computer placed in a separate room. Furthermore, the four gas tanks (Ar, CO, CO 2 , SO 2 ) are stored in a cabinet with an exhaust vent, next to the furnace. Four Alicat mass flow controllers for Ar (model MC, range up to 1000 cm 3 min −1 ), CO (model MC, range up to 500 cm 3 min −1 ), CO 2 (model MC, range up to 1500 cm 3 min −1 ), and SO 2 (model MCS, range up to 400 cm 3 min −1 ) are also placed inside the cabinet. Only a single gas line with the mixed gas leaves the cabinet and is connected to the alumina tube of the furnace. The mass flow controllers are chosen to achieve 3.4 mm s −1 flow velocity (400 cm 3 min −1 ) in the alumina tube, and Ar is used to flush the system at a rate of 8.5 mm s −1 (1000 cm 3 min −1 ). The mass flow controllers are controlled through a Lab-VIEW control interface, developed specifically for this furnace (Fig. 2), which enables dynamic control of the MFCs during an experiment and records temperature and gas flow rates. In brief, the control signal is sent to a mass flow controller from the independent unit (National Instruments cRIO-9066) which communicates with the control interface in the PC. This unit is also designed to control the security actions in case of component failures, such as overheating, inlet, which is pumped in with a pump built into the ECONOX controller. At the lower end: two 1 mm thick platinum wires with hooks (quench electrodes) and the thermocouple which extends out to roughly the height at which the sample is suspended. This sample holder has a total length of 68 cm. (b) A schematic drawing of the simple ceramic sample holder. At the top: the connector which transmits the current from the thermocouple to the computer. At the lower end: two 1 mm thick platinum wires with hooks (quench electrodes) from which the sample is suspended. This sample holder has a total length of 62 cm. malfunction of the air extraction system, and detection of toxic gases above safety thresholds. Gas set-point files are loaded into the control interface as .csv files and are used to control the gas flow rates during an experiment.

Safety systems
The safety of the user against exposure to toxic gases and high temperatures was a priority in this design; therefore multiple safety features were set in place.
The CO and SO 2 lines are connected to a security box that keeps the gas valves closed until the safety relay is activated. This system is connected to the gas sensors and the fume hood sensor in the room; if they exceed safety thresholds or if the fume hood stops functioning, the valves are automatically shut. The gas is cleaned downstream by passing it through two glass beakers with desiccants to get rid of moisture (that might come from the third beaker) and to precipitate sulfur. The dry gas then goes on to a third glass beaker and is bubbled in water with some carbonates to dissolve CO 2 , whereas the SO 2 reacts with carbonates to form sulfates. A minimal amount of CO is also dissolved. The remaining gas is then fed straight into the evacuation hood, all while keeping the entire system closed (Fig. 3). The furnace room is also evacuated by an additional fume hood to ensure maximum safety in case of a leak elsewhere. The fume hoods evacuate the room at a rate of 125 L s −1 . Using a Dräger Xam 2500 portable gas detector at the roof exit, it was determined that of the gas used, roughly 98 % of the CO and 1 % of the SO 2 are released into the atmosphere. Therefore, the CO concentration is within the 8 h exposure limit (< 50 ppm) at the roof exhaust. Additionally, the CO and SO 2 have extra gas lines for purging to ensure a safe release of the tank when they need to be replaced.
The safety of the operator is assured by a portable gas sensor which is typically placed near the furnace, in addition to the array of fixed gas detectors. Human presence in the furnace room is minimized. An Internet protocol (IP) camera is used to monitor the system remotely. In case of a leak, masks designed to filter out SO 2 and CO are available. In case the cooling water stops flowing for any reason, the safety mechanism automatically stops the furnace to prevent overheating.

Operation
The furnace is usually kept at 700 • C and is heated to the experimental temperature as the first step. The f O 2 sensor bearing sample holder is then introduced and lowered to the hotspot, and the furnace is flushed with Ar for approximately 5 min, which corresponds to a volume several times the size of the system. The desired gas mixture is passed through for ∼ 30 min until a stable f O 2 reading is recorded. Once again, the system is flushed with Ar. The f O 2 sensor bearing sample holder is then swapped out for the simple holder that suspends the prepared sample. At this point, the sample is placed near the top of the furnace, which is ∼ 600 • C cooler than the hotspot. The furnace is immediately flushed with argon to prevent the occurrence of any reaction with air. The sample is gradually lowered to the hotspot in about four to five steps over 30 min; this is done to avoid thermal shock to the sample capsule, as some olivine crystals have invisible micro-defects that might result in them breaking apart due to the thermal shock. Once the thermocouple reads a temperature that is 20-25 • C from the target, the experiment is  The stainless-steel quench bowl is heated to ∼ 160 • C with the external heating coil to prevent the build-up of sulfur. The gas leaves through two exits, each going into a beaker containing desiccants to precipitate sulfur. The gas is then collected from the bottom of the beaker and sent into the third beaker containing water with carbonates to react with CO 2 and SO 2 . Via another tube, the gas is then let into the fume hood which releases the gas into the atmosphere at the roof. started with the flow of the CO-CO 2 -SO 2 gas mixture. The flow of the gas mixture immediately raises the temperature by approximately 5 • C, likely due to the change in the thermal conductivity. The bubbling of the outgoing gas is monitored in the third beaker with water. The bubbling is typically intermittent and is considered a reliable indicator of normalcy. The intermittence is likely caused by the CO 2 dissolution into the water, and there are no such observations for the Ar gas flow. The gas mass flow controllers are also closely monitored for the first 30 min of the experiment to ensure stability. At the end of the experiment, the sample is quenched into the hollowed quench bolt (a piece in the centre, Fig. 4) in the bottom bowl of the furnace. The quench bowl is heated to a temperature of ∼ 160 • C during the experiment with a heating coil to prevent sulfur precipitation. This temperature was chosen as it is sufficiently hotter than the melting point of S (112.8 • C) and lower than the melting point of the O-rings. The quench bolt is cooled externally with a beaker of cold water prior to dry-quenching the sample, which falls into the cooled inner quench bolt. The inner bolt is clean, dry, and free of sulfur, as any sulfur liquid drips down pools outside the outer S-drain bolt. The gas flow and the heating coil are then stopped; the furnace is flushed with argon for 5 min before recovering the sample. The furnace is then cooled down to its idle temperature of 700 • C. The accuracy of the thermocouple (Type B) on the simple ceramic sample holder was determined by melting a pure gold metal strip. The strip was suspended across the platinum quench electrodes of the sample holder, and the height of the holder was adjusted so that the strip was at the hotspot. The temperature was then raised to 1100 • C on the temperature controller, and the system was left to equilibrate until the temperature measured by the thermocouple remained constant. The set temperature on the controller was then increased by 1 • C every 2 min until the Au strip melted. The resistance across the quench electrodes was measured using a voltmeter (∼ 1.5 ) to determine when the Au strip melted. A minute after the set temperature on the controller reached 1110 • C, the resistance jumped to 0.5 M , indicating that the circuit was open due to the gold melting. At that moment, the temperature indicated by the thermocouple, adjusted for its distance from the hotspot, was 1061.3 • C, which soon stabilized to 1062 • C. Accounting for hysteresis, the melting temperature of Au is 1064.43 • C (Haynes et al., 2016); therefore, the thermocouple records temperatures 2.43 • C colder than the actual temperature.

Hotspot position
Determination of the hotspot location was carried out at temperatures of 1200, 1300, and 1400 • C, which were set using the temperature controller, with a steady flow of Ar at 3.4 mm s −1 (400 cm 3 min −1 ), representative of the experimental gas flow rate and temperature. Results show that the hotspot is located at ∼ 32 ± 2 mm (Fig. 5), below the midpoint of the furnace enclosure. The location of the hotspot at each temperature is nearly the same, and the shift is insignificant. The hotspot location is also found to be independent of the gas flow rates for the range of 1.7-6.8 mm s −1 (200-800 cm 3 min −1 ). An offset of 46-52 • C is observed in the temperature range of 1200-1400 • C between the two thermocouples, due to their position. The thermocouple used by the temperature controller is outside the alumina tube, whereas the recorded temperature is from the thermocouple on the ceramic sample holder within the furnace and records a lower temperature. No change in this offset has been noted in over 200 experiments. The temperature variation at the hotspot is within 10 • C over 5 cm, giving the gas mixture approxi- Figure 5. Temperature measured vs. depth from the midpoint of the furnace. y-axis labels indicate the temperature set at the controller, which reads a higher temperature near the heating element. The curve was derived from a square fit and then maximized to find the location of the hotspot. mately 7 s to react completely at the desired temperature before reaching the hotspot and the sample.

Mass flow controller performance
The performance of the mass flow controllers (MFCs) was inspected from the log produced by the control interface. They are precise and send the requested flow of gas with an error of ±2 cm 3 min −1 . According to Alicat, the MFCs for CO, CO 2 , and Ar have a precision of ±0.6 % of the reading or ±0.1 % of the full scale (±0.8 % of the reading or ±0.2 % of the full scale for SO 2 ), whichever is greater. It is attempted to keep the values for individual gases above 20 cm 3 min −1 (0.17 mm s −1 ), at which errors are 5 % relative to the requested flow rate. This lower limit of the MFC sets the lower limit of attainable f S 2 in the furnace, indicated by the field in Fig. 8. The use of an MFC capable of regulating lower flow rates can allow further reduction of f S 2 in the gas mixture.

f O 2 -f S 2 calibrations
The atmosphere in the furnace is controlled by the thermodynamic equilibrium imposed by the mixture of the gases and the temperature. Here, the partial pressures of nine gas species are calculated (O 2 , CO, CO 2 , COS, CS 2 , S 2 , SO, SO 2 , SO 3 ) by Gibbs-free-energy minimization using iterative Lagrangian optimization (modified from an open-source code, Thermochemical Equilibrium Abundances (TEA), developed by Blecic et al., 2016, based on the methods of Eriksson et al., 1971, andWhite et al., 1958, using thermodynamic data from Chase, 1998). Based on the Gibbsfree-energy minimization mentioned above, Fig. 6 shows the Figure 6. The range of f O 2 -f S 2 conditions that can be attained within the furnace at 1300 • C, using the gas mixture CO 2 -CO-SO 2 . Calculated using Gibbs-free-energy minimization. Panel (a) shows the proportion of CO in the gas mixture, while (b) shows the proportion of SO 2 in the gas mixture. The redlines show the limits of the reaction boundaries as reference, sloped line: 3FeS + 2O 2 = Fe 3 O 4 + 3/2 S 2 , vertical line: range of f O 2 -f S 2 conditions that can be attained in the furnace at 1300 • C, using the gas mixture CO-CO 2 -SO 2 . The gas mixture begins to react as it enters the top of the furnace. When it reaches the 200 mm long hotspot (50 mm long within 10 • C), the mixture is expected to be fully reacted and consists of the equilibrium gas species as calculated above.
The f O 2 -f S 2 condition within the furnace was calibrated in steps. First, f O 2 control was verified using the Fe-FeO and Ni-NiO reactions at their respective buffers ( ironwüstite (IW) and NNO) followed by tests using the f O 2 sensor. Once calibrated for f O 2 , f S 2 was calibrated by the synthesis of FeS or magnetite at the FeS-magnetite buffer.
The f O 2 condition within the furnace can be monitored using a Y-doped zirconia sensor and mineral reactions. The mineral reactions used were the Fe-FeO and the Ni-NiO reactions. Both were tested using different methods. The f O 2 condition was verified at the NNO buffer by an oxidation reaction of Ni metal (Ni-NiO, Ni + 1/2 O 2 NiO), which was visually determined after quench. Two experiments (Table 2), starting with pure NiO powder, in a Si glass capsule were carried out at 1350 • C. In both cases, the result follows expectations; i.e. we have NiO at a positive NNO condition and Ni at a negative NNO condition. Based on the above calibration, the f O 2 condition within the furnace is accurate to within ±0.27 log f O 2 .
The f O 2 condition was verified at the IW buffer (Fe-FeO, Fe + 1/2 O 2 FeO) at 1300 • C by measuring the resistance across a 1 mm wide pure Fe strip suspended between the Pt quench electrodes. The log f O 2 was varied from −11.5 to −10 in 30 min. The log f O 2 according to the gas mixtures calculated at 1300 • C was varied from −11.5 to −10 in 30 min. The equilibrated T at the hotspot was 1288.3 • C. In these conditions, IW = 0 at log f O 2 = −10.8826. The initial resistance across the Fe strip was 2.53 . The resistance stayed constant for ∼ 5 min and then slowly started increasing at the rate of 0.02-0.04 min −1 . A minute after crossing the buffer (calculated from the gas flow values at the MFCs), the resistance starts increasing dramatically (> 0.2 min −1 ), and within 5 min, the strip breaks (likely due to fast conversion to FeO). Given a 40 s travel time for the gas from the MFC to the hotspot, the corrected log f O 2 at the time of crossing the buffer is −10.886. Based on this assessment criteria, the f O 2 is accurate to within ±0.002 log f O 2 .
The zirconia sensor works by measuring the potential difference generated by the difference in oxygen activity at both ends (gas mixture and reference air). The f O 2 at the sample location is calculated using the following formula (derived from Nernst's law): where mV is the potential difference detected by the sensor in millivolts, T is the temperature in kelvin, pO 2 ref is the reference pressure of oxygen in bars (0.21 is used), and pO 2 is the partial pressure of oxygen at the location of the sensor in bars. The constant value of 0.0496 is derived from the equation 2.303·R n·F , where R is the gas constant; F is the Faraday constant; 2.303 is the conversion factor from natural log to log 10 ; and n equals 4, since four electrons are exchanged in the reaction between the zirconia membrane, the reference air (atmosphere), and the atmosphere of the furnace. The location of the zirconia sensor is 3 cm above the sample; however the temperature of the thermocouple of the sample holder is applied to the equation to determine the f O 2 at the sample location. We tested attainment of target f O 2 conditions by running experiments under several gas mixtures for 30 min each at temperatures of 1200-1400 • C (Table 1). We observed that the atmosphere within the furnace appears to be more oxidizing (by about 0.4 log units) than expected from the gas mixture (i.e. the sensor reads higher f O 2 ). The shift observed is independent of the gas proportions introduced and flow rate. We conclude this is a shift in the zero point of the sensor and that the absolute readings of the sensor are not a reliable representation of the atmosphere within the furnace. However, the zirconia sensor accurately tracks the changes in the f O 2 conditions of the furnace. This is discussed further in the following section (Sect. 4.4). We therefore conclude that the f O 2 sensor in our possession is inaccurate but precise.
With the f O 2 calibrated, reversal experiments were run to calibrate f S 2 and ensure the target condition is accurately attained. Since an f S 2 sensor does not exist, calibrated redox conditions were used to verify the sulfide's formation and disappearance. These experiments follow the standard operating procedure outlined in Sect. 3. Here, the key reactions taking place in this system are as follows: 3FeS(s) + 2O 2 (g) Fe 3 O 4 (s) + 3/2S 2 (g).
Reaction (R3) describes the exchange of sulfide or oxide, and as f O 2 is calibrated and controlled by the gas mixture, f S 2 can be determined once the equilibrium boundary of the reaction is experimentally determined. Practically, the reaction boundary needs to be bracketed using Reaction (R3) by starting with FeS and Fe 3 O 4 . The precise details of these experiments are provided in Table 2. Figure 7a shows the sample (no. 179) that started out with pure Fe powder and is now an iron sulfide. For this starting condition, the iron metal is expected to quickly react to form metastable oxide (3Fe + 2O 2 Fe 3 O 4 ) and sulfide (2Fe+S 2 2FeS) before transforming into the equilibrium phase. Figure 7b is a BSE image of the sample (no. 163) that started out with pure FeS 2 powder and was oxidized to form an iron oxide, most likely Fe 3 O 4 . For this starting condition, the disequilibrium FeS 2 reacts to form iron oxide. Since this experiment was run for just 2 h and at a lower temperature, some cores (brighter part) of unreacted iron sulfide are visible in the image (Fig. 7b). Based on these observations, we conclude that the f S 2 imposed by the gas mixture reflects that of the calculation. Figure 8 indicates the conditions of all sulfide-forming experiments and one iron oxide formation (no. 163). The horizontal axis is f O 2 relative to the fayalite-magnetite-quartz buffer (Frost, 1991), and the vertical axis is f S 2 relative to Reaction (R3). All our experiments are consistent with the expected boundary of Reaction (R3), and the boundary is bracketed within ±0.6 log units of f S 2 . However, given the accuracy of the MFCs in reproducing the f O 2 values accurately, we see no reason why the f S 2 values should be significantly off. We can therefore safely assume that the f S 2 accuracy of the furnace is as good as the f O 2 accuracy.

Dynamic conditions
To perform experiments that simulate conditions of degassing, it is necessary to constrain which f O 2 conditions in the furnace can be dynamically controlled and the rates at which they can be changed. There exists an upper limit to how fast the atmosphere within the furnace can be changed, as fast changes in the gas mixture may be dampened during transport. The limit was determined by continually modifying the gas mixture in the furnace at fixed time intervals by creating steps in the gas set-point file. A constant flow rate of 400 cm 3 min −1 was used during this calibration. Figure 9 shows the result of the calibration, in which 328 S. Prabha-Mohan et al.: One-atmosphere high-temperature CO-CO 2 -SO 2 gas-mixing furnace  The tests with shorter periods at 2, 4, and 8 min reproduced the sinusoidal curves, but the amplitude of these curves was dampened, indicating diffusion and/or incomplete reaction within the gas mixture during transport (not shown here). The typical "degassing" experiment imposes a constant rate of change in f O 2 and f S 2 after the system reaches equilibrium. Such experiments simulate certain degassing conditions observed in nature, such as the ones modelled by Burgisser and Scaillet (2007). This constant change in the f O 2 and f S 2 of the system requires precise management of the gas flow rates as they are not linear or systematic. The degassing path is divided into multiple sub-steps, and the corresponding gas mixture is calculated for each point and interpolated between them. The points are closely spaced where the changes in the gas mixture are greatest. Based on the previous oscillation experiment, the duration of the change can be reduced to approximately 5 min, and it is possible to simulate a 1 log unit change in partial pressure of oxygen and sulfur. Based on a coupled chemical-physical model calculation, f S 2 in a magma can change by approximately 0.6 log units as it ascends from 1 kbar to the surface (Burgisser and Scaillet, 2007).

Potential applications
The creation of a gas-mixing furnace that can faithfully replicate equilibrium and disequilibrium conditions lends itself to plenty of potential applications in high-temperature, nearsurface geology, especially volcanology. In addition to the direct applications discussed in the Introduction, other potential applications are discussed below: 1. High-temperature and lower-pressure conditions are believed to cause significant sulfur isotope fractionation between gas species, melt, and sulfides (de Moor et al., 2013;Mandeville et al., 2009;Marini et al., 2011Marini et al., , 1998;   Sakai et al., 1982); such fractionation can be replicated and quantified.
2. Quantifying the behaviour of sulfur and the amount of oxygen in the sulfide, the Fe and Ni content in the associated melt in various oxygen and sulfur fugacity conditions is of great interest (Brenan and Caciagli, 2000;Mungall et al., 2005;Shima and Naldrett, 1975) and can tell us more about the accumulation of chalcophile elements in sulfide and their associated melts (Holzheid and Lodders, 2001;Parat et al., 2011;Rose and Brenan, 2001). This is of great importance as many chalcophile elements are economically important and understanding the behaviour of sulfur and its isotopes in volcanic and magmatic processes can improve our understanding of how and where deposits of these chalcophile elements might have formed or may form in the future (Kiseeva et al., 2017;Kiseeva and Wood, 2015). 3. Volcanic degassing processes are known to cause isotopic and elemental fractionation in a wide range of elements (Mandeville et al., 2009;Sossi et al., 2019;Vlastelic et al., 2021). This fractionation occurs in dynamic conditions where temperature and gas compositions change with time. Many questions also exist around the HT reactions that take place at volcanic craters, fumaroles, and vents (Kusakabe et al., 2000;Shinohara et al., 2015;Vlastelic et al., 2021). By determining the isotopic fractionation factors between gas, melt, and sulfide in equilibrium and disequilibrium conditions, we can estimate critical parameters related to degassing processes, such as the ascent rate of the magma or residence time of the magma in the reservoir. The known isotopic composition of the gas used in the experiment can be used to characterize the gasmelt fractionation. This is because the input gas floods the system and, given the minute size of the sample, the overall isotopic equilibrium of the system is controlled by the gas composition. To experiment with disequilibrium conditions, multiple experiments are carried out, and each is stopped at a different point to capture compositional and isotopic snapshots of the melt. The melt exchanges sulfur isotopes with a gas with a constant isotopic composition but varying molecular species, which imposes disequilibrium.
Data availability. No data sets were used in this article.
Author contributions. SPM prepared the manuscript and figures with contributions from KTK and DFN. The initial furnace design was created by KTK and constructed by AM and FP. Initial calibration and testing were carried out by DFN and KTK. Further modifications to the design were carried out by SPM, KTK, and FP.