Water solubility in trachytic and pantelleritic melts: an experimental study

Solubility experiments were performed on a trachyte and a pantellerite from Pantelleria. The trachyte has SiO2 = 65.2 wt%, Al2O3 = 15.2 wt% and a peralkaline index (P.I. = molar[(Na2O+ K2O)/Al2O3]) ∼ 1 while the pantellerite has SiO2 = 72.2 wt%, Al2O3 = 11 wt% and a P.I. = 1.3. Solubility experiments were performed in the pressure range of 50–300 MPa at T = 950 °C for the trachyte and 50–200 MPa at T = 850 °C for the pantellerite. The water content of experimental glasses was determined by Karl Fischer titration, elemental analyser and FT–IR spectroscopy. Water content appears similar in both compositions for analogous pressure conditions, varying from ∼2.5 wt% at 50 MPa to ∼5.8 at 200 MPa. Comparison of the experimental data with the thermodynamic models shows that the empirical model of Moore et al. [1998] better matches the experimental data for both compositions, while the thermodynamic models of Papale et al. [2006] and Ghiorso and Gualda [2015] tend to overestimate melt water content, probably due to the paucity of peralkaline melt compositions used for calibration. One inference of these new solubility data is on MI-derived depth of felsic reservoirs at Pantelleria which deepens from 2.4 to 3.5 km.


Introduction
Water solubility in silicate melts has been widely investigated because of the important effects it has on magma properties, such as viscosity and density [e.g. Dingwell et al., 1996, Richet et al., 2 Pierangelo Romano et al. Moore et al., 1998, Behrens and Jantos, 2001, Tamic et al., 2001, Liu et al., 2005, Lesne et al., 2010 for which the general effects of composition, particularly that of alkalis and silica contents [e.g. Johannes and Holtz, 1996, Dingwell et al., 1997, Holtz et al., 2000, Schmidt and Behrens, 2008, Stabile et al., 2018, and pressure-temperature are well established [e.g. Holtz et al., 1992. More recently, several studies have also focused their attention on intermediate compositions such as andesite-dacite [Mandeville et al., 2002, Botcharnikov et al., 2006, phonotephrite [Behrens et al., 2009] and phonolite [Carroll and Blank, 1997, Larsen and Gardner, 2004, Iacono-Marziano et al., 2007, Schmidt and Behrens, 2008 melts. In contrast, less attention has been paid to compositions such as trachytic melts, which usually represents the parental composition to peralkaline rhyolite/phonolite melts: only Di Matteo et al. [2004] and Fanara et al. [2015] have obtained solubility data on such compositions. Trachyte and pantellerite magmas are involved in explosive eruptions, some of which can be of large magnitude (e.g. Tambora 1815, Campanian Ignimbrite 39 kyr, Green Tuff 45 kyr), and it is therefore vital to accurately constrain water solubility in such melts. Enlarging the spectrum of available solubility data is also crucial to further improve already existing solubility models [e.g. Papale et al., 2006, Ghiorso andGualda, 2015].
Here, we present new experimental H 2 O solubility data obtained on trachytic and pantelleritic melts over the pressure range of 50-300 MPa. The trachyte and pantellerite used in this study are representative of the most abundant magmas on Pantelleria island. The water contents dissolved in experimental glasses have been determined by FT-IR spectroscopy, Karl Fischer titration (KFT) and elemental analyser (EA). The necessity of performing additional experiments arises from the observed dependence of water solubility on the composition of silicate melts, as several studies quoted above have illustrated. In particular, a small variation in melt chemical composition influences not only water solubility but also the extrapolation of molar absorption coefficients. FT-IR spectroscopy is a common tool used to determine the water content of natural volcanic products (melt inclusions and glassy groundmasses) but its accuracy is strongly related to the use of proper absorp-tion coefficients, calibrated on compositions similar to that of natural glasses investigated. Otherwise, the extrapolation of absorption coefficients outside the calibration range may lead to a significant error in water determination. For this reason, in this work, we also used the results of FT-IR spectroscopy, KFT and EA analyses to determine the IR molar absorption coefficients for 4500 and 5200 cm −1 bands by using the peak intensity and a straight line correction procedure , and reference therein].

Starting material and charges preparation
We used, as starting materials, natural samples representative of the trachytic and peralkaline rhyolite (pantellerite) of Pantelleria Island (Table 1). Pantelleria, located in the Sicily Channel, is the type locality for pantellerite, an iron-rich rhyolite characterized by alkali/alumina molar ratio [(Na 2 O + K 2 O)/Al 2 O 3 ] > 1 [Macdonald, 1974]. Pantellerites and trachytes are the most abundant rocks outcropping in the island and also in other localities, usually representing the felsic end member in continental rift zones (Pantelleria, Ethopian Rift valley, etc.) or in oceanic islands (Mayor Island, Azores, etc.). The trachyte (sample GTT) and pantellerite (sample Pan113) specimen used in this study have been the subject of phase equilibrium experiments by Di Carlo et al. [2010] and Romano et al. [2018Romano et al. [ , 2020. To prepare the starting materials, the selected rocks were finely crushed and the resulting powders were melted twice at 1300°C for 3 h, then rapidly quenched. The glasses obtained were grinded in an agate mortar obtaining a powder of 10-40 µm mesh size. This technique, used in several previous experimental works [Di Carlo et al., 2010, Romano et al., 2020, allows to produce a crystal-free and homogeneous glass with a composition similar to the starting bulk rock. Experimental charges were prepared by loading 50-100 mg of dry glass powder and 5-10 mg distilled water in gold capsules, ensuring volatile saturation at experimental conditions. At 50 MPa, the amount of water loaded was slightly lower, to avoid the risk of capsule rupture. Once prepared, capsules were weighed and welded, then weighed again to verify any water loss during welding. Gold capsules were preferred to platinum or gold-palladium containers in order to minimize iron loss to the capsule walls, since the duration of the experiments exceeded the 24 h. The applied temperature was higher than 950°C for trachytes and 850°C for pantellerites, conditions that in both cases ensured to be above the liquidus in our experimental runs [Di Carlo et al., 2010, Romano et al., 2018.
Prior to experiments, capsules were left at least for 24 h in an oven at 120°C to allow the homogenisation of water throughout the capsule.

Experimental strategy
Experiments were performed in internally heated pressure vessels (at Institut des Sciences de la Terre d'Orléans) equipped with a molybdenum furnace and a fast quench device. As pressurizing medium, a mixture of argon and hydrogen loaded in the vessel at room temperature was used: the Ar/H 2 ratio was fixed to reach redox condition around the FMQ buffer [e.g. Scaillet et al., 1992, Romano et al., 2020. The pressure was recorded by a transducer calibrated against an Heise-Bourdon tube gauge (uncertainty ± 20 bar) while the temperature was continuously controlled through S-type thermocouples (accuracy ±5°C). Experiments were performed in a T range 950-985°C for trachytes and at 850°C for pantellerites. The pressure range explored varied from 50 to 300 MPa (at 50 MPa intervals) and the duration of each experiment was between 96 and 125 h. Experiments were terminated using the dropquench technique, which ensures a fast quench of the capsules from the top to the bottom of the vessel (quench rate >100°C·s −1 ). After the experiments, capsules were weighed to check for leaks and if the post run weight difference was >0.4 mg, capsules were discarded. The excess of H 2 O within the capsule was confirmed by the occurrence of small water drops during the opening of the capsules, establishing that saturation conditions were reached during the experiment. For each capsule, around onethird of the run product was embedded in epoxy resin for qualitative analysis by a scanning electron microscope (SEM-EDS) and chemical analysis by an electron microprobe (EMP) while the remaining part was used to determine water contents.

Scanning electron microscope and electron microprobe analyses
Glass fragments mounted in epoxy resin were first observed with a SEM and then analysed by an EMP to check glass for composition and homogeneity. EMP analyses were performed with a Cameca SX-Five located at CNRS-ISTO laboratory (France). Analytical conditions were an accelerating voltage of 15 kV, a beam current of 6 nA and a defocused beam of 10 × 10 µm. Na and K were analysed first to minimize volatile loss during the analyses. Microprobe analyses of alkali-rich hydrous glasses in some cases resulted problematic due to the alkali (mostly sodium) migration out of the excited volume. As suggested by Morgan and London [2005], we have used a current density of 0.006 nA/µm 2 that limits the Na migration under the beam and also provides adequate excitation of all elements.

Karl Fischer titration and elemental analyser
The total water contents of experimental glasses were determined by KFT and EA (EA Thermo Scientific Flash 2000) both housed at CNRS-ISTO laboratory.
In KFT analysis, 10-20 mg of glass fragments are placed within a Pt crucible in an induction furnace and heated up to 1300°C. Water liberated during the heating of glass fragments reacts quantitatively in a titration cell with iodine. The EA allows determining the concentration of H-C-N in experimental samples. In our case, only H 2 was measured to infer water content. The glass in the EA is first heated to 1800°C upon which the substance oxidizes into simple compounds, which are in turn quantified by thermal conductivity detection. Moussallam et al. [2015] used this technique to determine water and CO 2 contents in kimberlite melts. The EA was calibrated before each analytical session using pure pyrophyllite, which is known to contain 5 wt% of stoichiometric H 2 O in its crystalline structure. The experimental samples for the EA were prepared by grinding 1-10 mg of experimental glass into a fine powder and then loading it in a tin capsule folder. The instrument was tested repeatedly using an international standard with known amounts of carbon, hydrogen and analyses on experimental glasses synthetized in this study, haplogranite (HPG8) composition [Dingwell et al., 1997] and Eburru obsidian [Scaillet and Macdonald, 2006].
nitrogen and also using other hydrated glass samples (as internal standards) prepared for this purpose and for which water concentration was already known (haplogranite) or already analysed with KFT. In Figure 1 (Table 2), the water content determined by EA is compared with literature data for haplogranite composition [Dingwell et al., 1997], Eburru obsidian [Scaillet and Macdonald, 2006] and KFT analyses of this study; a good correlation is observed between H 2 O contents retrieved by these different techniques.

Density determination
Densities of anhydrous trachytic and pantelleritic starting glasses were measured by weighing single pieces of glass in air and ethanol. The determination of hydrous glass densities was more difficult because of the limited amount of run products; a minimum weight of 15 mg is indeed required to perform an accurate Archimedean density measure. Only for three hydrous trachytic glasses were reliable density data obtained ( Figure 2a, Table 3), which was used to derive an empirical relationship of the form ρ = −4.4017(H 2 O melt ) wt% + 2329. Standard deviation for  Dingwell et al. [1997].
Pan 113 and Eburru natural sample were analysed by Karl Fischer titration.
density values was estimated from repeated measures of glass fragments. For comparison, the densities of hydrous experimental glasses were also calculated using the procedure of Ochs and Lange [1999], calculated values being 4 to 11% higher than those measured. A less clear relationship has been found between density and melt water content for pantelleritic experimental glasses ( Figure 2b, Table 3) and calculated densities are 1 to 6% higher relative to measured ones. To maintain internal consistency, we preferred to use the density measured or calculated from the above empirical relationship, knowing that density measurements are affected by relative errors in the range 0.1-2.5%.

FT-IR spectroscopy
Fourier Transform Infrared analyses were performed using a Nicolet 760 Magna spectrometer connected to an IR microscope (fluxed with H 2 O-free compressed air) housed at ISTO laboratory. The spectrometer was equipped with an MCT detector, a Globar light source and a KBr beam splitter. Absorption spectra were acquired for each sample in the range 1000-6000 cm −1 with 128 scans and a resolution of 4 cm −1 . Experimental glass chip samples were placed on an NaCl crystal plate, which is transparent to the infrared beam; for each glass chip 3-5 points (50 × 50 µm spot size) were analysed. Experimental glasses were double-polished to obtain 100-300 µm thin wafers, and then carefully washed under acetone and ethanol in order to remove any impurities. The areas of the sample analysed were always checked optically to avoid bubble or impurities. The thickness of each glass chip was measured microscopically during the analysis, with an uncertainty of ±1-3 µm. Absorption spectra were converted in H 2 O concentration according to the Beer-Lambert equation: where C (wt%) is the water content dissolved in the melt, MW the molecular weight of water, A the height of the absorbance peak, t the thickness (cm) of the glass wafer, ρ the density (g/l) of the glass and ε the molar absorption extinction coefficient (L·mol −1 ·cm −1 ). To calculate the water content, we considered the absorption bands at 5200 cm −1 (molecular H 2 O) and 4500 cm −1 (OH − concentration) of the IR spectra. Only for a few samples the fundamental OH − stretching vibration at 3530 cm −1 was also measured, which required glass wafers of extremely thin thickness (<30 µm); thickness measures in the range 40-60 ± 2 µm are affected by a relative error of 2-5% while for a higher range of 6 Pierangelo Romano et al. thickness measures (∼200 µm) an incertitude in ±10 correspond to a relative error of 5%.

Microscopic observation and major element analyses
The run products consisted of microlite-free glasses.
Most of the quenched glasses present also bubbles at the melt/capsule interface indicating fluid saturation during the experiment; however, bubble content never exceeded the 0.1 vol.%. Major element composition of experimental glasses (Table 1)

Water solubility in trachyte and pantelleritic melts
The total H 2 O content dissolved in experimental glasses is reported in

Analyses of IR spectra and molar absorption coefficients calibration
The bands of IR spectra at 5200 cm −1 and 4500 cm −1 for water species are the most used to determine water content in geological samples (Figure 4). The former is assigned to the combination of stretching and bending mode of H 2 O molecules (H 2 O m ) while the second to the combination of stretching and bending of hydroxyls (OH − ) groups [Stolper, 1982].
To avoid an incorrect estimation of absorption band peak height, a baseline correction is required. Following the previous studies of Withers and Behrens [1999], Ohlhorst et al. [2001] and Mandeville et al. [2002], we used a simple straight baseline [TT baseline, Ohlhorst et al., 2001] connecting 4500 and 5200 backgrounds. We preferred to use this simple procedure knowing that a possible small underestimation of OH − band cannot be excluded, owing to the presence of the 4000 cm −1 band. As suggested by Stabile et al. [2020], depending on sample compo- sition and spectral variation, a critical evaluation of the approach used for the baseline correction has to be taken into account. Absorbance (peak intensity), density and thickness of the glass allow us to determine molar absorption coefficients (ε), in our case for the bands of molecular H 2 O (5200 cm −1 ) and OH − groups (4500 cm −1 ). The total water content dissolved in the experimental glasses is the sum of water species (H 2 O m and OH − ) and the Beer-Lambert equation can be rewritten as: C (wt%) = C H2O + C OH− = 18.01 * A 4500 /ρ * t * ε 4500 + 18.01 * A 5200 /ρ * t * ε 5200 .
Knowing C (wt%) from KFT or EA, Equation (2) can be solved to determine the molar absorption coefficients by plotting the absorbance of 5200 and 4500 bands normalized to density, thickness and water content. Using a linear regression the intercepts on X and Y axes are the ε 4500 and ε 5200 values (Figure 5a,b). Molar absorption coefficients were calculated assuming a linear relationship between the water species concentration and intensities of IR spectra bands. We obtained, for trachytic glasses, absorption coefficients of 0.88 L/mol·cm −1 (ε 4500 ) and 1.6 L/mol·cm −1 (ε 5200 ) while for pantellerite glasses ε 4500 is 1.21 L/mol·cm −1 and (ε 5200 ) 1.89 L/mol·cm −1 . Molar absorption coefficients are calibrated for a melt water content ranging between 2.6 and 7.4 in trachytic glasses and between 2.6 and 5.8 for pantelleritic ones. Regarding the uncertainties on the thickness and density measures, we can consider that the epsilon values are affected by relative errors in the range 2-7%.

Comparison with previous works
(i) Trachyte compositions. Di Matteo et al. [2004] and Fanara et al. [2015] presented solubility experi-ments on Campi Flegrei trachytes, obtained between 25 and 50 MPa (Figure 6a,b). Campi Flegrei trachytes are broadly similar in composition to Pantelleria trachytes, having a lower content in calcium, aluminium, potassium together with similar or slightly higher contents in silica, iron and sodium. Our solubility curve differs significantly from that of Di  and Moore et al. [1998]. The solubility curves show notable differences, the water content of our pantelleritic melt at 100 MPa being ∼1 wt% lower and up to ∼2 wt% lower at higher pressure relative to the data of Stabile et al. [2018]; our solubility data are very simi- Figure 6. Water contents solubilized in the melt as a function of pressure obtained in this study for trachyte and pantellerite compared with (a) Di Matteo et al. [2004], Fanara et al. [2015] experimental data and solubility models of Moore et al. [1998], Papale et al. [2006] and Ghiorso and Gualda [2015]. Experimental results on pantellerite compared with (b) Moore et al. [1998], Stabile et al. [2018] experimental data and solubility models of Moore et al. [1998], Papale et al. [2006] and Ghiorso and Gualda [2015]. lar to the experimental results of Moore et al. [1998] at 150 and 200 MPa (Figure 6b) remembering that Stabile et al. data are based on synthetic pantellerite. The water solubility model of Papale et al. [2006] predicts water content higher than 1 wt% to 1.5 wt% than those of our experiments while considering the Ghiorso and Gualda [2015]'s model, the difference between the calculated solubility and the experimental data ranges between 0.2 wt% and 1.2 wt%. Over- Table 4. Chemical composition of trachyte [Di Matteo et al., 2004] and pantellerite [Moore et al., 1998, Stabile et al., 2018 used in other experimental studies Di Matteo et al. [2004] Fanara et al. [2015 This study Moore et al. [1998] Stabile et al. [2018]  all, the model of Moore et al. [1998] better reproduces our experimental solubility data, these latter resulting lower by 0.1 to 0.5 wt% of the calculated ones.
In contrast, the model of Papale et al. [2006] tends to overestimate significantly the melt water contents of our compositions.
(iii) The effect of small variations in bulk composition. To a first approximation, these differences in water solubility can be ascribed to small differences in chemical composition between the starting materials considered (Table 4). As stated previously, the dependence of water solubility in trachyte and pantelleritic melts with respect to small variation in chemical composition has not yet been fully investigated. For instance, the trachyte of this study has aluminium and potassium content almost 3 wt% lower with respect to the trachyte and trachyphonolite of Di Matteo et al. [2004] and Fanara et al. [2015] which has a higher K/Na ratio (0.71 respect to 0.46 of GTT trachyte); on the contrary, the GTT trachyte has higher silica and iron contents. Considering pantelleritic composition [Moore et al., 1998, Stabile et al., 2018, the differences in water solubility are also related to slight differences in melt composition (Table 4). Our pantellerite and that of Moore et al. [1998] indeed present similar oxides concentration with respect to that of Stabile et al. [2018], which has lower iron, aluminium, sodium and higher silica content. Water solubility has been also shown to correlate positively with metal ionization potential and negatively with Al 2 O 3 content Jantos, 2001, Mysen, 2002]. Di Matteo et al. [2004] and Stabile et al. [2018] explored the effect of alkalis on water solubility evidencing that for both trachyte and pantellerite melts, higher Na content favour H 2 O incorporation in the melt, in accord with other studies on silicic compositions [Holtz et al., 1995, Dingwell et al., 1997, Carroll and Blank, 1997]. All these effects may explain the range in water contents obtained at the same P-T between the different sets of experiments. We note, however, that our comparison with previous experimental studies does not take into account the difference in experimental temperature, yet a clear relationship between temperature and water solubility has been evidenced for rhyolite [Yamashita, 1999] and phonolite melts [Schmidt and Behrens, 2008]. Such an effect still needs to be determined for peralkaline felsic compositions.

Effect of molar absorption coefficients
The molar absorption coefficients determined for our experimental glasses differ from those determined previously for trachytic by Di Matteo et al. [2004] and Fanara et al. [2015] and pantelleritic glasses by Stabile et al. [2020]. As stated above, molar absorption coefficients are determined assuming that they are constant over the range of water contents considered [Silver et al., 1990]. The general agreement between water contents obtained with EA ( Figure 7a) and KFT (Figure 7b) versus the total water content obtained by FT-IR spectroscopy, using the calibrated absorption coefficients, supports such an assumption (Figure 7a,b). There is a general consensus on the fact that absorption coefficients are composition dependent. Dixon et al. [1995] evidenced that absorption coefficients (ε 4500 , ε 5200 ) correlate positively with the sum of cations Si and Al, suggesting that such a linear relationship can be used to calculate molar absorptivities of silicate melts. Other studies have shown a relation of ε with Na/K ratio, silica content and excess of aluminium [Stolper, 1982, Silver et al., 1990 while Ohlhorst et al. [2001] reported a systematic increase of absorption coefficients from basalt to rhyolite and proposed a parabolic equation to calculate ε as a function of silica content. The calculated absorption coefficients proposed for glasses having similar SiO 2 content to our trachyte and pantellerite [Ohlhorst et al., 2001] are broadly similar to that determined in our study; yet, several authors [e.g. Dixon et al., 1995, Mandeville et al., 2002 suggest that a direct calibration for a specific composition yield to most accurate results. With respect to the published epsilon values on trachyte, we note that our epsilon values are 14-45% higher for ε (H 2 O) and 26-44% lower for ε (OH); for pantelleric glasses, ε (H 2 O) is 2% lower and ε (OH) is 20% higher with respect to the value of Stabile et al. [2020]. Anyway, we used the absorption coefficients obtained in the studies on experimental glasses of Di Matteo et al. [2004], Fanara et al. [2015] for trachyte and Stabile et al. [2020] for pantellerite (Table 5). Using the absorption coefficients obtained by Di Matteo et al. [2004] to determine the water content from our IR spectra decreases the retrieved values; for instance, at 100 MPa the melt water content decrease from 5 to 4.6, i.e. values are 8% lower. On the other hand, if our absorption coefficients are used on the data of Di Matteo et al. [2004], their water contents increase by 0.1-0.6 wt%. Doing the same exercise on pantelleritic glasses, and using the absorption coefficients of Stabile et al. [2020], the water content increase by 0.04-0.4 wt%. Using our absorption coefficients for the data of Stabile et al. [2020], the water content decrease by 0.2-0.5 wt%. Differences in absorption coefficients influence above all the estimation of water species that could reflect either an incorrect estimation of OH − relative to molecular species or the opposite. Taking as an example our results on pantellerite glasses, the OH/H 2 O m ratio increases from 0.07 to 0.2 when the absorption coefficients of Stabile et al. [2020] are used (i.e. ε 4500 is 0.98 L/(mol·cm −1 ) instead of 1.21 L/(mol·cm −1 ) as determined in our study).

Water speciation
Water species concentration (OH − and H 2 O m ) in the experimental glasses as a function of total water content are shown in Figure 8a,b, for trachyte and pantelleritic glasses. As already shown for other quenched hydrous glasses, the OH − group concentration is higher at low water content whereas molecular H 2 O dominates for water contents above ∼3 wt% [e.g. Stolper, 1982, Xue andKanzaki, 2004, Dixon 92 et al., 1995, Carroll and Blank, 1997, Withers and Behrens, 1999. Our experimental glasses have water contents higher than 2.5 wt% so that we cannot appreciate the water speciation at low water concentrations. Anyway, it is evident in Figure 8a that the experimental glass with ∼2.5 wt% represents the point at which OH − and H 2 O m are present at equivalent concentrations. This appears less clear in the experimental glasses of pantelleritic composition, in which molecular water is slightly higher than OH − at total water contents around 2.5 wt%. However, the speciation model of Dixon et al. [1995] and the experimental results on water speciation seem to be in good agreement within errors. Molecular water and hydroxyl species proportion in the experimental glasses at room temperature may reflect the structural re-equilibration during quench; this leads to an increase of molecular H 2 O over hydroxyls with respect to the proportions at high temperature and pressure conditions [Silver et al., 1990, Zhang et al., 1995, Nowak and Behrens, 1995, Carroll and Blank, 1997. For instance, Withers and Behrens [1999] investigated the effect of quenching on rhyolitic glasses containing 3 and 5 wt% H 2 O. Glasses with initial OH/H 2 O ratios of 1 and 0.6, after being reheated at 800°C and at 300 MPa and rapidly quenched, showed a slight increase of their OH/H 2 O ratios (1.17 and 0.72, respectively); they inferred that effect of quenching rate on water speciation was relatively minor at that melt water content. This gives information on the apparent equilibrium temperature (T ae ) [Zhang, 1994[Zhang, , 1999 considered the equilibrium temperature corresponding to the final speciation. However, we do not have sufficient data to appreciate if such a conclusion equally applies to other melt water contents.

Volcanological implications
Trachyte and pantellerite of Pantelleria show a similar capability to solubilize water content in the melt. We note only a subtle difference between the two compositions, water content in pantellerites being slightly higher (but always <1 wt%) with respect to trachytes. As evidenced by several works performed on melt inclusions [Lanzo et al., 2013, and references therein], water contents characteristic of natural pantellerites of Pantelleria are in the range 2-4.9 wt%. To date, the difference in H 2 O content has not yet been related with changes in eruptive dynamic of pantelleritic products (lava flow, strombolian and sub-plinian eruptions). However, a part of water, other factors can play an important role in controlling the eruptive styles [Cassidy et al., 2018] as already identified for mafic volcanism in a typical basaltic volcano as Mt. Etna [Moretti et al., 2018, and reference therein]. The study of water solubility in silicate melts has several implications for the study of hydrous natural glasses (melt inclusions and glass matrices). Ideally, melt inclusion studies allow one to assess the pre-eruptive volatile content of a magma, which can be then used to infer the (minimum) entrapment pressure (i.e. storage conditions if entrapment occurs in the reservoir) considering adequate solubility models and absence of diffusion towards the host crystal, post entrapment crystallization, etc. In the case of Pantelleria, several recent studies have determined the volatile content of melt inclusions trapped in phenocrysts of trachytic and pantelleritic rocks [Neave et al., 2012, Lanzo et al., 2013, Romano et al., 2019 of some key explosive eruptions. The maximum water content measured in melt inclusion with pantelleritic composition was 4.5 wt% [Lanzo et al., 2013] while the water content estimated for trachytic melt through thermodynamical modelling and mass balance by White et al. [2009] ranges between 3.4 and 4 wt%. Romano et al. [2019] report water contents in melt inclusions of trachytic composition (0.15 to 1.05 wt%), suggesting that melt inclusions experienced volatile loss through the crystals. Considering 4 wt% of water dissolved in the trachytic melt, the protracted crystallization to generate a pantellerite liquid derivate [∼80% of crystallization; White et al., 2009, Romano et al., 2018 will lead to an excess of volatiles in the magmatic system and a consequent outgassing of volatiles, which could be the deep source of the diffusive degassing observed at Pantelleria nowadays.
Considering now the results of phase equilibria obtained on representative trachyte and pantellerite magmas, the magmatic felsic system of Pantelleria is inferred to evolve close to water saturation conditions, water and chlorine (during the late stage of evolution) being the dominant volatile species [Lanzo et al., 2013]. Evolution at volatile saturation conditions has been recently proposed for similar peralkaline magmatic systems elsewhere [Iddon and Edmonds, 2020] with reservoirs located at pressures <200 MPa. Using the highest water content measured in trachytic and pantelleritic melt inclusions at Pantelleria (water is the main volatile component, with no detectable CO 2 ), the model of Papale et al. [2006] gives storage pressures between 55 and 75 MPa for trachyte and pantellerite, respectively. Our experimental data constrain a pressure range between 70 and 110 MPa, or values which are 30% higher. Evolution of peralkaline felsic magmas in shallow reservoirs around 100 MPa is suggested by phase equilibrium and geophysical constraints [e.g. Di Carlo et al., 2010, Mattia et al., 2007. Although such a difference might appear small, it corresponds to a difference of the top-depth of magma storage region (calculated assuming an average crustal density of 2600 kg·m −3 ) of ∼1 km. Underestimating by ∼1 km the depth of magma storage region in a shallow volcanic system such as Pantelleria might have important consequences with regard to volcanic hazard issues; for instance, it will affect the estimated travel time needed for the magma to breach the surface. The difference in solubilities may also affect the magma porosity evolution during magma ascent as simulated in numerical models, with an obvious impact on the explosive-effusive transition: the critical porosity of 70% beyond which fragmentation of an ascending magma might occur, will happen at deeper levels in a model using our solubility data instead of that of Di Matteo et al. [2004]. This underscores the importance of using accurate volatile solubility laws, hence the need to carry out specific hydrothermal experiments on representative compositions, if those laws are not available or if models are used too far from the compositional calibrated range.

Conclusions
This study was aimed at constraining the water solubility in trachyte and pantellerite melts from Pantelleria island. It is largely accepted that the solubility of water in silicate melt is composition-specific not only at the broad scale (silica-rich, silica-poor magmas) but also for variations in major element abundance within the same magma type (for instance, Pantelleria trachyte and Campi Flegrei trachyte). We thus performed this experimental study in order to add to the existing experimental data set on water solubility in felsic magmas, the missing data for trachytes straddling the chemical divide that separates metaluminous from peralkaline magmas. Our results place tighter constraints on water solubility-related arguments specific to eruptive phenomena (eruption style, pre-eruptive conditions, magmatic outgassing) but also on fundamental petrogenetic processes in which dissolved water plays a role (i.e. liquid line of descent oriented towards pantellerite melts).
Experiments were carried out in a pressure range 50-300 MPa in an IHPV, with only H 2 O as the sole volatile component. The FT-IR study of experimental glasses resulted in the definition of the following molar absorptivity coefficients ε 4500 = 0.88 L/mol·cm −1 and ε 5200 = 1.6 L/mol·cm −1 for trachyte and ε 4500 = 1.21 L/mol·cm −1 and ε 5200 = 1.89 L/mol·cm −1 for pantellerite glasses. Absorption coefficients are different from those estimated in other trachyte [Di Matteo et al., 2004, Fanara et al., 2015 or phonolite [Carroll and Blank, 1997] melt compositions, and hence contribute to define an important composition-specific variable needed for FT-IR studies of such melt compositions.
Our experimental data were tested against the models of Moore et al. [1998], Papale et al. [2006] and Ghiorso and Gualda [2015]. The Papale et al. [2006] model overestimates the H 2 O melt content for Na-rich magma, such as Pantelleria felsic magmas. This effect has significant implications when applied to the study of melt inclusions, whose volatile abundances are used to derive the depth of magma storage region. In the specific case of Pantelleria, trachyte and pantellerite storage pressure derived from the model of Papale et al. [2006] are about 50-70 MPa, for a dissolved H 2 O content in the range 3.4-4.5 wt%, which implies a magma storage depth of ∼2.4 km. Using our experimental solubility model the calculated depth increases to ∼3.5 km, a 30% difference, not trivial for such shallow magma bodies.