Vol48/04/2005def 647 ANNALS OF GEOPHYSICS, VOL. 48, N. 4/5, August/October 2005 Key words noble gases – solubility – degassing – silicate melts 1. Introduction When investigating processes involving magmas and coexisting vapor phases, a direct observation of the system is obviously impossi- ble. The information available comes from out- gassed vapors that are released at the surface and from fluid inclusions trapped in igneous products. This implies that the observation con- ditions are different from those at the source, and that the pristine gas content may be modi- fied by chemical reactions, unless unreactive species are considered. Furthermore, degassing of the dominant volatiles (H2O, CO2, S, Cl and F) strongly modifies the physico-chemical properties of the magmatic system and plays an active role in driving the degassing process, whereas components at trace concentrations passively follow the evolution of the system. Among these latter, poorly soluble species mark out the degassing process evolution very well, as degassing involves the exsolution of the vapor phase where they are preferentially partitioned. All the highlighted features can be Noble gas solubility in silicate melts: a review of experimentation and theory, and implications regarding magma degassing processes Antonio Paonita Istituto Nazionale di Geofisica e Vulcanologia, Sezione di Palermo, Italy Abstract Noble gas solubility in silicate melts and glasses has gained a crucial role in Earth Sciences investigations and in the studies of non-crystalline materials on a micro to a macro-scale. Due to their special geochemical features, noble gases are in fact ideal tracers of magma degassing. Their inert nature also allows them to be used to probe the structure of silicate melts. Owing to the development of modern high pressure and temperature technologies, a large number of experimental investigations have been performed on this subject in recent times. This paper reviews the related literature, and tries to define our present state of knowledge, the problems encountered in the experimental procedures and the theoretical questions which remain unresolved. Throughout the manuscript I will also try to show how the thermodynamic and structural interpretations of the growing experimental dataset are greatly improving our understanding of the dissolution mechanisms, although there are still several points under discussion. Our improved capability of predicting noble gas solubilities in conditions closer to those found in magma has allowed scientists to develop quantitative models of magma degassing, which provide constraints on a number of questions of geological impact. Despite these recent improvements, noble gas solubility in more complex systems involving the main volatiles in magmas, is poorly known and a lot of work must be done. Ex- pertise from other fields would be extremely valuable to upcoming research, thus focus should be placed on the structural aspects and the practical and commercial interests of the study of noble gas solubility. Mailing address: Dr. Antonio Paonita, Istituto Nazio- nale di Geofisica e Vulcanologia, Sezione di Palermo, Via Ugo La Malfa 153, 90146 Palermo (Italy); e-mail: paoni- ta@pa.ingv.it 648 Antonio Paonita recognized in the geochemical group of the no- ble gases: they are all chemically inert and volatile, and are present in trace concentrations in magmatic systems. Moreover, their physico- chemical properties gradually vary within the group with the increase in their atomic mass or size. Thus, although noble gases do not thermo- dynamically define the magmatic system, their abundance ratios are a very powerful tool in studying degassing processes, on a micro to a macro-scale, from magmatic intrusions to man- tle outgassing and the development of atmos- phere. For these studies, noble gas solubilities are a crucial parameter and a number of investiga- tions show a growing interest in the subject (Carroll and Draper, 1994; Carroll and Webster, 1994 and reference therein; Draper and Carroll, 1995; Chamorro-Perez et al., 1996, 1998; Shi- bata et al., 1996, 1998; Shackelford, 1999; Nuc- cio and Paonita, 2000; Paonita et al., 2000; Wal- ter et al., 2000; Schmidt and Keppler, 2002), al- so thanks to the development of modern experi- mental and analytical technologies. Experimen- tal devices make it possible to reproduce mag- matic temperatures from atmospheric to mantle pressures, and to measure the noble gas concen- trations down to parts per million (ppm). Never- theless, technical questions that mainly concern the distinction between the dissolved and ex- solved fractions of gas must be solved to obtain accurate vapor-melt partition coefficients. Moreover, available experimental data on noble gas solubility are far from exhaustive, particu- larly as regards helium. The existing studies on- ly partially cover the compositional ranges of magmatic systems and possible thermo-baric conditions, thus their use in modeling natural processes has still to be optimized. Due to their inert nature, noble gases devel- op van der Waals-type interactions with silicate melts and seem to display a preferentially «physical solubility» (i.e. the noble gas atom is similar to a fixed-volume sphere entering a rigid network). Nonetheless, their dissolution mechanism is not fully defined. Empirical or semi-empirical (Lux, 1987; Jambon, 1987; White et al., 1989; Chennaoui-Aoudjehane and Jambon, 1990; Carroll and Stolper, 1993; Shi- bata et al., 1998; Nuccio and Paonita, 2000) ap- proaches are still the most frequently used tools to describe the noble gas solubility and to ex- trapolate the experimental data towards com- plex magmatic systems, whereas statistical-me- chanics models (Doremus, 1966; Studt et al., 1970; Shacklford et al., 1972) work within a very narrow range of temperature, pressure (T, P) and composition. Current efforts also aim at investigating the structural environment of no- ble gas atoms inside the silicate network, prom- ising results coming from the use of X-ray tech- niques (Wulf et al., 1999). In this framework of growing high-quality scientific literature, this paper describes the sol- ubility of noble gas in silicate melts, which is mainly applied to geological problems involv- ing magma degassing in its wide range of pos- sible conditions. Throughout the text, this paper will try to draw attention to existing literature and the main problems encountered in experi- ments and modeling. While stressing the im- portance of noble gases as tools in studying magma degassing, the paper shows two things: a) what type of data is needed for geological ap- plications; b) which are the most promising guidelines of research for achieving these re- sults. 2. Experimental techniques The investigation of noble gas solubility in silicate materials at magmatic conditions needs to experimentally reproduce high temperature and atmospheric up to mantle pressures. Mod- ern technology has designed and optimized re- sistance furnaces and Externally Heated and In- ternally Heated Pressure Vessels (EHPV and IHPV) to perform experiments at atmospheric pressure up to 500-600 MPa at magmatic tem- peratures. In this equipment, the noble gas un- der investigation is normally used as the pres- sure medium in the high pressure line and en- ters the unsealed noble-metal capsules where the sample is contained. In EHPV and IHPV experiments, noble gas can also be contained in a sealed capsule, that is loaded prior to the run by using a custom-designed loading device (Boettcher et al., 1989) or a gas-bearing mate- rial (Paonita et al., 2000). Experiments at pres- 649 A review on noble gas solubility and degassing in magmas sure up to mantle values (several GPa) have been performed using piston cylinder apparatus or multi-anvil cells. In these instruments, a sealed capsule (single or double) containing the sample and the noble gas is inserted, and high pressures are reached using a solid pressure medium. Run duration has to be long enough to attain equilibrium. Based on the diffusion rates of no- ble gases, the appropriate duration ranges from hours to months, depending on the melt or glass composition and temperature (longer periods are necessary for heavy gases in silica glasses). The use of small grain sized sample powders yields a shorter diffusive walk of gas atoms in melt, because the gas penetrates through the grain interstices before melting, resulting in faster equilibration (Jambon et al., 1986; Lux, 1989; Roselieb et al., 1992). The noble gas content in quenched samples is analyzed using different methods. The elec- tron microprobe has been frequently used for in situ analyses of Ar or heavier noble gases (White et al., 1989; Carroll and Stolper, 1991, 1993; Montana et al., 1993; Draper and Carroll, 1995; Chamorro-Perez et al., 1996; Schmidt and Keppler, 2002), and has obtained concen- tration profiles in quenched glasses with a spa- tial resolution of around 10 µm. X-ray induced emission by synchrotron radiation and proton beams have been applied to the analysis of Kr in albite (Carroll et al., 1993b). Light gases (He, Ne) cannot be detected by these methods, and bulk extraction techniques have been ap- plied to extract the dissolved gases (Hayatsu and Waboso, 1985; Jambon et al., 1986; Lux, 1987; Broadhurst et al., 1992; Roselieb et al., 1992; Shibata et al., 1996, 1998; Paonita et al., 2000). Solid samples are melted in a vacuum line, in order to extract the dissolved volatiles. Various analytical methods are capable of measuring the amount of extracted gas, namely Gas-Chromatography (GC), Mass Spectrome- try (MS) and Knudsen cell Mass Spectrometry (KMS). The main problems consist in distin- guishing between gas fractions from surface ad- sorption or trapping in fractures, secondary flu- id inclusions and truly dissolved volatiles in the silicate melt (see Roselieb et al., 1992). In order to avoid the gas fraction released from incorpo- rated bubbles, powdering of quenched glasses has been used and grain sizes below 100 µm were found to be necessary to remove most of the physically trapped gas. However, an exces- sively small grain size could cause significant gas loss by diffusion before analysis, at least as regards helium and, in part, neon, if the samples are not analyzed within a matter of hours (Pao- nita et al., 2000). Stepwise degassing methods, coupled to KMS or quadrupole MS measure- ments, have been developed to overcome these difficulties, based on the fact that different gas contributions are extracted in different ranges of temperature (fig. 1) (Roselieb et al., 1992, 1995; Paonita et al., 2000). In situ extraction techniques (ultra-violet laser beam) have been combined with mass spectrometry (Kelly et al., 1994; Brooker et al., 1998; Chamorro et al., 2002), and the concen- tration of dissolved argon in glass volumes hav- ing around a 10 µm diameter size have been ob- tained. These methods are very promising as they are able to remove the contribution from physically trapped gas and they can be used for lighter noble gases. Fig. 1. Release spectrum of helium from basaltic glasses heated at a rate of 10°C/min. Data from Paonita et al. (2000). Peak 2 provides the amount of dissolved gas, whereas peak 1 has been attributed to the decrepitation of fluid inclusions. 650 Antonio Paonita 3. Noble gas solubility: experimental results and modeling The main parameters affecting the solubili- ty of noble gases in silicate melts are pressure, temperature and composition of the coexisting liquid and vapor phases. In general (see table I), the concentration of any dissolved noble gas displays a linear relation with the gas pressure up to several hundred MPa (figs. 2a,b and 3). At Table I. Noble gas solubility in melt and glasses, and parameters for pressure and temperature dependence. Argon (a) P° T° V° (b) ∆H° (b) Ref. cc(STP)g−1bar−1 (MPa) (°C) (cc/mol) (cal/mol) Mg2Si2O6 2.0E-05 0.1 1500 1 Basalt 2.5E-05 0.1 1200 2 Bas.-andes. 9.0E-05 0.1 1200 16700 3 Tholeiite 2.3E-05 0.1 1200 3 Alk.Ol-basalt 4.6E-05 0.1 1300 3 Tholeiite 5.9E-05 0.1 1250 5700 4 Basanite 1.0E-04 0.1 1200 4500 5 Tholeiite 6.7E-05 0.1 1200 6800 5 Alk.Ol-basalt 3.6E-05 0.1 1200 14200 5 Andesite 1.5E-04 0.1 1350 5 Ugandite 4.5E-05 0.1 1350 5 NaAlSi3O8 2.9E-04 1500 1600 22.8 1962 6 KAlSi3O8 1.9E-04 1500 1600 20.6 1722 6 CaAl2Si2O6 2.6E-05 1500 1600 22.5 3373 6 CaMg2Si2O6 6.0E-06 1500 1600 21.5 3421 6 K2Si4O9 3.2E-04 1500 1400 22.6 3062 6 Basalt 4.7E-05 1500 1600 23.1 2225 6 Granite 3.9E-04 1500 1600 23.7 3086 6 SiO2 1.2E-03 0.1 700 16.44 −4761 7,13 (d) NaAlSi3O8 3.0E-04 10 1000 20 (c) 8 Basalt (synt.) 4.6E-05 0.1 1300 9 Rhyolite 4.6E-04 0.1 700 15.6 −1457 10,13 (d) NaAlSi3O8 3.0E-04 0.1 700 22 −1807 10,13 (d) KAlSi3O8 3.0E-04 240 700 22 −3447 10 Ol-tholeiite 4.0E-05 0.1 1275 23 10 Basalt 9.9E-05 0.1 1200 23 10 CaAl2Si2O6 5.6E-04 5000-12000 (e) 14 SiO2 5.0E-04 5000-12000 (e) 14 NS1 (f) 1.4E-04 197 1300 5410 15 NS2 2.6E-04 197 1300 9690 15 NS3 3.4E-04 197 1300 10840 15 NCS1 1.4E-04 197 1300 11390 15 NCS2 1.1E-04 197 1300 11320 15 NCS3 6.0E-05 188 1400 15 CMS1 2.6E-05 199 1500 15 CMS2 2.7E-05 199 1500 15 CMS3 2.5E-05 199 1500 15 Haplogranite 4.8E-04 3800 1800 24.3 17 Tholeiite 1.0E-04 3800 1800 24.7 17 651 A review on noble gas solubility and degassing in magmas Table I (continued). Helium (a) P° T° V° (b) ∆H° (b) Ref. cc(STP)g−1bar−1 (MPa) (°C) (cc/mol) (cal/mol) Mg2Si2O6 1.2E-04 0.1 1500 1 Tholeiite 5.6E-04 0.1 1250 7100 4 Basanite 6.8E-04 0.1 1200 1500 5 Tholeiite 5.5E-04 0.1 1200 4900 5 Alk.Ol-basalt 5.7E-04 0.1 1200 540 5 Ugandite 4.8E-04 0.1 1350 5 SiO2 7.9E-03 0.1 25 5.5-10.8 −641 7 (g) Basalt 3.8E-03 (5.1) 215 1160 16 (h) 2.0E-03 (3.3) 125 1160 16 2.6E-03 (3.3) 112 1160 16 Rhyolite 8.6E-03 (6.0) 215 1140 16 5.4E-03 (5.0) 165 1140 16 4.4E-03 (3.9) 112 1140 16 Neon (a) P° T° V° (b) ∆H° (b) Ref. cc(STP)g−1bar−1 (MPa) (°C) (cc/mol) (cal/mol) Mg2Si2O6 7.0E-05 0.1 1500 1 Bas.-andes. 2.6E-04 0.1 1200 4200 3 Tholeiite 1.9E-04 0.1 1200 3 Alk.Ol-basalt 2.2E-04 0.1 1300 3 Tholeiite 2.5E-04 0.1 1250 5000 4 Basanite 3.6E-04 0.1 1200 5100 5 Tholeiite 2.9E-04 0.1 1200 1300 5 Alk.Ol-basalt 1.7E-04 0.1 1200 12300 5 Ugandite 2.1E-05 0.1 1350 5 SiO2 4.9E-03 0.1 150 10.3 −1833 7 (g) NaAlSi3O8 1.5E-03 100 1000 16.8 (d) 8 Basalt (synt.) 2.8E-05 0.1 1300 9 NS1 (f) 8.4E-05 197 1300 8230 15 NS2 1.0E-04 197 1300 9590 15 NS3 1.2E-04 197 1300 7250 15 NCS1 7.0E-04 197 1300 13160 15 NCS2 5.0E-04 197 1300 10410 15 NCS3 2.6E-04 188 1400 15 CMS1 1.4E-04 199 1500 15 CMS2 1.4E-04 199 1500 15 CMS3 1.9E-04 199 1500 15 Kripton (a) P° T° V° (b) ∆H° (b) Ref. cc(STP)g−1bar−1 (MPa) (°C) (cc/mol) (cal/mol) Basalt 1.3E-05 0.1 1200 2 Bas.-andes. 2.1E-05 0.1 1200 19300 3 Tholeiite 8.0E-06 0.1 1200 3 Alk.Ol-basalt 1.6E-05 0.1 1300 3 Tholeiite 3.0E-5 0.1 1250 19500 4 Andesite 1.1E-04 0.1 1350 5 Basanite 8.2E-05 0.1 1200 3700 5 652 Table I (continued). Tholeiite 4.8E-05 0.1 1200 9400 5 Alk.Ol-basalt 2.1E-05 0.1 1200 19000 5 Ugandite 3.0E-05 0.1 1350 5 NaAlSi3O8 1.8E-04 50 750 29 (d) 8 Basalt (synt.) 1.3E-05 0.1 1300 9 NaAlSi3O8 1.4E-04 0.1 800 25 11 SiO2 3.8E-4 0.1 800 11 NS1 (f) 5.7E-05 197 1300 2890 15 NS2 1.0E-04 197 1300 16200 15 NS3 1.5E-04 197 1300 19190 15 NCS1 6.3E-05 197 1300 11840 15 NCS2 4.5E-05 197 1300 15020 15 NCS3 2.4E-05 188 1400 15 CMS1 1.1E-05 199 1500 15 CMS2 1.1E-05 199 1500 15 CMS3 1.3E-05 199 1500 15 Xenon (a) P° T° V° (b) ∆H° (b) Ref. cc(STP)g−1bar−1 (MPa) (°C) (cc/mol) (cal/mol) Basalt 9.0E-06 0.1 1200 2 Tholeiite 9.0E-6 0.1 1200 3 Andesite 8.3E-05 0.1 1350 5 Basanite 3.1E-05 0.1 1200 4400 5 Tholeiite 2.5E-05 0.1 1200 2900 5 Alk.Ol-basalt 9.2E-06 0.1 1200 14400 5 Ugandite 1.0E-05 0.1 1350 5 Basalt (synt.) 1.1E-05 0.1 1300 9 NaAlSi3O8 9.7E-05 1500 1600 28.5 2488 12 KAlSi3O8 1.0E-04 1500 1600 28 2560 12 K2Si4O9 1.2E-04 1500 1400 29.4 3876 12 NS1 (f) 3.5E-05 197 1300 −6360 15 NS2 4.9E-05 197 1300 35650 15 NS3 7.5E-05 197 1300 30140 15 NCS1 2.6E-05 197 1300 11120 15 NCS2 2.2E-05 197 1300 23420 15 NCS3 1.2E-05 188 1400 15 CMS1 9.3E-06 199 1500 15 CMS2 6.2E-06 199 1500 15 CMS3 7.2E-06 199 1500 15 Tholeiite 2.8E-05 3500 1800 35.76 17 (a) Noble gas solubility at P° and T°; the dissolved concentrations are normalized by noble gas pressure. Uncer- tainties are normally within 10%. (b) V° and ∆H° are molar volume and enthalpy of solution of the reference state (see text). (c) Calculated from the isothermal data of 8 by fitting eq. (3.6). (d) Solubilities and related parameters are reliable for both pure Ar and Ar in He-Ar mixtures. (e) Pressure range of the experiments. (f) NS refers to the system Na2O-SiO2 (1, 2, 3 have about 35-65, 31-70, 24-77 wt% respectively). NCS is the system Na2O-CaO-SiO2 (1, 2, 3 have about 25-9-64, 19-19-60, 11-31-56 wt% respectively). CMS is the system CaO-MgO-SiO2 (1, 2, 3 have about 37-8-54, 26-16-56, 14-26-59 wt% respectively). (g) Calculated from data from Frank et al. (1961), Shackelford et al. (1972), Shelby (1972a,b, 1976). (h) Values in hydrous melts at Ptot ≈ PH2O and He concentration in vapour > 0.1 mol%. Numbers in parenthesis are Antonio Paonita 653 A review on noble gas solubility and degassing in magmas wt% of dissolved water. Solubilities at Ptot ≈ PH2O + PCO2 were obtained at PH2O very close to that given in table and do not display appreciable differences with respect to given values, therefore they were not reported. Ref.: 1 – Kirsten (1968); 2 – Fisher (1970); 3 – Hayatsu and Waboso (1985); 4 – Jambon et al. (1986); 5 – Lux (1987); 6 – White et al. (1989); 7 – Carroll and Stolper (1991); 8 – Roselieb et al. (1992); 9 – Broadhurst et al. (1992); 10 – Carroll and Stolper (1993); 11 – Carroll et al. (1993); 12 – Montana et al. (1993); 13 – Draper and Carroll (1995); 14 – Chamorro-Perez et al. (1996); 15 – Shibata et al. (1998); 16 – Paonita et al. (2000); 17 – Schmidt and Keppler (2002). Fig. 2a,b. Dissolved Ar as a function of Ar pressure in a) rhyolite (800-900°C) and basalt (1200°C); b) haplogranite and basalt (1500°C). Data are from Car- roll and Stolper (1993) and White et al. (1989). Curves were calculated by using eq. (3.6) with pa- rameters given in table I. The linear relation between concentration and pressure, as well as the higher sol- ubilities in the acidic compositions, should be noted. higher pressure, the solubility (hereafter de- fined as the ratio between gas concentration in melt and gas fugacity, unless differently speci- fied) decreases, and recent experiments indicate Fig. 3. Neon and krypton solubilities in albite glass at 1000 and 750°C, respectively. Data are taken from Roselieb et al. (1992). Fig. 4. Effect of very high pressure on Ar solubility in basaltic and haplogranitic melts between 1500 and 2000°C. Data from Schmidt and Keppler (2002); curves calculated by using eq. (3.6) with the param- eters given by the authors (see table I). At about 5 GPa, dissolved Ar in the melt reaches a threshold concentration which does not increase even when Ar pressure is increased by some GPa. a b 654 Antonio Paonita that a threshold concentration is reached, at least for Ar and Xe, at very high pressures (4-5 GPa and 5 GPa, respectively, fig. 4; Schmidt and Keppler, 2002). The dependence of solubil- ity on temperature is very modest and is often lower than the experimental uncertainties. 3.1. Dependence on the composition of melt and glass The parameters which most affect the solu- bility of noble gases in silicate melts and glasses are the composition of the solvent and the atom- ic radius of the noble gas (Doremus, 1966; Shackelford et al., 1972; Shelby, 1976; Hayatsu and Waboso, 1985; Jambon et al., 1986; Lux, 1987; Broadhurst et al., 1992; Roselieb et al., 1992; Carroll and Stolper, 1993; Shibata et al., 1996, 1998; Shackelford, 1999; Paonita et al., 2000; Schmidt and Keppler, 2002). In general, solubility increases with the increasing concen- tration of SiO2 in the melt while smaller noble gases display higher solubilities (fig. 5). These variations suggest that noble gas atoms dissolve in holes (free spaces) of the silicate melt struc- ture, following the so-called physical mecha- nism of dissolution (Doremus, 1966; Studt et al., 1970; Shackelford et al., 1972; Shelby, 1976). Five basic structural units form the main net- work of silicate melts and glasses, namely SiO2 (three-dimensional network), Si2O22− (sheet), Si2O64− (chain), Si2O76− (dimer), SiO44− (mono- mer) (Virgo et al., 1980). Silicon and other tetrahedrally coordinated cations (Al, Fe3+, Ti4+) act as network-formers, whereas Na+, K+, Mg2+, Ca2+, Fe2+ break the silicate polymers and thus are called network-modifiers. According to Shi- bata et al. (1998), the fully polymerized, three- dimensional SiO2 structure of silicate melts cre- ates a number of holes where gas atoms can be accommodated, whereas the less polymerized structures of the other units have fewer free spaces, therefore less volume is available to store gas atoms (i.e. the holes are too small). This is why silica-rich compositions display the highest noble gas solubilities. A quantitative expression of the degree of polymerization is given by the ratio NBO/T, namely the number of Non-Bridg- ing Oxygens per atom of a Tetrahedrally coordi- nated cation (Brawer and White, 1975; Mysen et al., 1985). The three-dimensional network has bridging oxygens only (NBO/T= 0, that is to say the highest polymerization), whereas the ratio grows when other units are present; if TO44− monomers are the only units in the melt, then NBO/T= 4. Shibata et al. (1998) found a direct relationship between noble gas solubility and de- gree of polymerization in the silicate melts, ac- cording to which gas atoms dissolve preferen- tially in more polymerized, low NBO/T melts (fig. 6). They proposed the following equation to calculate the vapor-melt equilibrium constant of noble gas i: (3.1) where ABO and ANBO are the fraction of the units with bridging oxygens and the units with bridg- ing and non-bridging oxygens, while KiBO and Ki NBO are the equilibrium constants that refer to each group of units. The model was calibrated for noble gases from Ne to Xe, but unfortunate- ly helium was not treated. From fig. 6 it can be seen that the relationship between solubility and the NBO/T ratio is also quite evident when plot- ting the available experimental data for helium, ln ln lnK A K A Ki i i BO BO NBO NBO= + Fig. 5. Solubility (expressed under form of Henry’s constant) of the noble gases as a function of their size at about 1200°C. Solubility exponentially increases from light to heavy noble gases. 655 A review on noble gas solubility and degassing in magmas thereby suggesting that the solubility mecha- nism is analogous. A similar approach has al- ready been proposed by Jambon (1987) and Chennaoui-Aoudjehane and Jambon (1990), al- though no clear interpretation in terms of melt structure was given by the authors. Further- more, even though the model by Shibata et al. (1998) does work for a large range of natural melt compositions, significant deviations from the relationship solubility-degree of polymer- ization occur in melts containing large amounts of Al (or network-formers other than Si). How- ever, the measured solubilities appear to be lower than those predicted from the NBO/T ra- tio (fig. 6). Although further data on noble gas solubilities in aluminosilicate melts are neces- sary to confirm this hypothesis, it is possible that cations (Na, K, Ca, Mg), added to balance the surplus of negative charges due to Si versus Al substitution, occupy structural vacancies otherwise available for noble gases. Fig. 6. Effect of melt polymerization (expressed by the NBO/T ratio, see text) on the solubility of some noble gases (modified from Shibata et al., 1998). Da- ta for Ar are a compilation selected by Shibata et al. (1998), He and Ne from Jambon et al. (1986) and Lux (1987). Ar data are divided into two classes: T/(T+Si) ratio higher (circles) than or lower (trian- gles) than 0.1, where T is the number of tetrahedral- ly coordinated cations other than silicon (i.e. Al, Ti). It is evident that the NBO/T ratio fails to predict the low Ar solubilities in aluminosilicate melts. Fig. 7a-c. Relationships between noble gas solubil- ity and a) ionic porosity of melt, b) melt density and c) melt molar volume (modified from Carroll and Stolper, 1993). In all the graphs, the squares are sol- ubilities in CaO-MgO-Al2O3-SiO2 (CAMS) melts from Broadhurst et al. (1992). Henry’s constant and molar volumes were computed on eight-oxygen ba- sis melts. Linear fits of data, with the related R 2 val- ues for Ar, are also displayed. Solubilities are much better correlated to ionic porosity than other parame- ters. The anomalous behavior of CAMS melts is un- clear (see text for discussion). a b c 656 Antonio Paonita Based on the mechanism of interstitial disso- lution in silicate networks, relationships have al- so been found between noble gas solubility and the physical properties of melt that give the gauge of the free space in the silicate network (fig. 7a- c). In particular melt density (Lux, 1987; White et al., 1989), melt molar volume (Broadhurst et al., 1992) and melt Ionic Porosity (IP, Carroll and Stolper, 1993). Carroll and Stolper (1993) clearly show that ionic porosity works better than the other properties (density, molar volume) in predicting solubility, and is thus a better expres- sion of the free space in the melt structure (see fig. 7a-c). In fact, ionic porosity is defined as (3.2) where vca is the total volume of cations plus an- ions in one gram of melt and vL is the specific volume of melt. The relation found by Carroll and Stolper (1993) between the solubility loga- rithm and ionic porosity is (3.3) where m and n are fitting parameters to be cali- brated by linear regression of the experimental- ly measured solubility of the noble gas in sev- eral melt compositions, and Henry’s constant Kh is defined as (3.4) where fi is the noble gas fugacity in the gas phase. Carroll and Stolper (1993) estimated these parameters for all stable noble gases, by using experimental data acquired at 1000- 1400°C and 0.1 MPa. In spite of its general ac- complishment, ionic porosity fails to predict solubility in CaO-MgO-Al2O3-SiO2 (CMAS) melts. Furthermore, although it estimates the overall free volume in melt, the same volume could be distributed either in a few large holes or in many small cavities. Following Shibata et al. (1998), the incorporation of network- modifier cations probably generates a kind of free space that is unable to accommodate noble gas atoms (i.e. too small interstitial sites). Nearly all of the available experimental data (Carroll and Webster, 1994 and reference there- /K xfh i i= ln K m nIPh- = + ( / )v vIP 100 1 Lca= - in; Draper and Carroll, 1994, 1995; Chamorro- Perez et al., 1996, 1998; Shibata et al., 1996, 1998; Shackelford, 1999; Walter et al., 2000; Schmidt and Keppler, 2002) have been obtained by equilibrating a silicate liquid with a gas phase composed of solely noble gases, which intro- duces a large uncertainty in the determination of the real solubility of noble gases in magmas. In fact, the gas phases coexisting with magmas are normally H2O-CO2 dominated vapors; more rarely, sulfur can reach several percents, where- as halogens and N2 remain minor components and the noble gases occur in trace concentrations (Gerlach and Nordlie, 1975; Giggenbach, 1996). Even at some GPa of total pressure, the partial pressures of noble gases are very low (a few bars; Carroll and Webster, 1994). In a recent ar- ticle, Paonita et al. (2000) experimentally inves- tigated the effect of H2O and CO2 on helium sol- ubility, in order to simulate data at more geolog- ically relevant conditions. They measured the Fig. 8. Effect of dissolved water concentration on the helium Henry’s constant, the latter expressed in MPa and calculated on 8-oxygen basis melt (modi- fied from Nuccio and Paonita, 2000). Experimental data (symbols) from Paonita et al. (2000) were ob- tained in the condition Ptot ≈ PH2O in the range 100- 200 MPa, and about 0.1 mol% He in the vapor. Curves were computed by using the model by Nuc- cio and Paonita (2000), which accounts for the effect of H2O and CO2 on both noble gas solubility in melt and fugacity in vapor (see text). 657 A review on noble gas solubility and degassing in magmas isothermal He solubility in basaltic and rhyolitic melts at pressures of 100 up to 200 MPa, in the presence of mixed H2O-CO2 vapors (0 to 50 mol% H2O) with ≈ 0.1 mol% He. The solubility increased by about a factor of three with 3 wt% H2O added to the melt, whereas a further addi- tion of water, up to 6 wt%, seemed to have very little effect on the solubility (fig. 8). The effect of CO2 addition (up to 0.05 wt%) was within the experimental uncertainties or slightly negative. Based on these data, Nuccio and Paonita (2000) found that the helium solubility is in very good agreement with the model of Carroll et al. (1993) (eq. (3.3) with the same parameters as those for dry melts, fig. 9), provided that the ionic porosi- ty melt is computed by taking into account the presence of dissolved H2O and CO2 (namely the H2O and CO2 partial molar volumes). It is obvi- ous that the H2O and CO2 concentrations in the coexisting gas and liquid phase must be con- strained by the solubility of these two volatiles. Nuccio and Paonita coupled their model to the thermodynamic model by Papale (1999) for sol- ubility of H2O-CO2 mixed vapors in silicate melts. They obtained the following equation for ionic porosity: (3.5) where r and i are ionic radius and molar fraction of each of the I ions forming the melt; x and v° are the molar fraction and the reference molar vol- ume of each of the n oxides in melt at P and T, and wi, CO2 are binary interaction coefficients be- tween dissolved CO2 and other oxides (wi, CO2 val- ues are taken from Papale, 1999; reference molar volumes of oxides are from Lange and Char- michael, 1987; and ionic radii from Shannon and Prewitt, 1969). By inserting the ionic porosity ++ j z x ( ) xx P x w r x x P w P x w i IP 100 1 1 3 4 1 1 , , ,P T i i i n I z z P T z z n j j z n CO CO CO CO ion ion ion CO CO CO CO 1 3 1 1 CO 1 2 2 2 2 2 2 2 2 2 = - + - + + - - o r o ! ! = = + z 1= j 1= c c ^ h < > F H Z [ \ ] ] ] ] ) _ ` a b b bb 4 / / / / value from eq. (3.5) into eq. (3.3), the helium Henry’s constant can be considered valid for H2O-CO2 bearing melts. The empirical parame- ters m and n were the same as those obtained for dry melts (fig. 9). In order to calculate the helium concentration in the melt by eq. (3.4), the authors computed the noble gas fugacity by using a mod- ified Redlich-Kwong equation of state for the H2O-CO2-noble gas (De Santis et al., 1974). With the Nuccio and Paonita (2000) model, you can theoretically calculate the solubility of other noble gases by assuming that m and n pa- rameters in eq. (3.3) are the same as those for dry melts. Although the last assumption has on- ly been experimentally confirmed for He, due to the lack of experimental data, the authors predicted that the solubilities would also in- crease with water addition, and become similar when 6-7 wt% H2O are dissolved in melt. The relation between He solubility and wa- ter content of the melt can be tentatively dis- Fig. 9. Relationship between helium solubility (Kh in MPa and on 8-oxygen basis melt) versus Ionic Porosity (IP) in H2O and H2O-CO2 bearing silicate melts (H2O and CO2 contents up to 6 and 0.05 wt%, respectively). IP was calculated by taking into ac- count the effect of the main volatiles (see text). Line was calculated by using eq. (3.3) with m and n pa- rameters in Carroll and Stolper (1993) for dry melts, and it displays a very good agreement with the ex- perimental data (data from Paonita et al., 2000). 658 Antonio Paonita cussed from a structural point of view. Water dissolution mainly occurs as OH− groups at low concentrations of dissolved H2O, whereas mo- lecular H2O becomes more and more important with increasing water content in basaltic to rhy- olitic melts (McMillan, 1994; Dixon et al., 1995; Nowak and Behrens, 1995, 2001; Ihinger et al., 1999; Withers et al., 1999; King and Hol- loway, 2002). Although it is generally accepted that water depolymerizes silica and alkali-sili- cate compositions, a similar effect is not clear in aluminosilicate melts (such as magmas). As first recognized by Sykes and Kubicki (1993), water could break Si-O-Al bonds and have a depolymerizing effect on the network. On the contrary, the Si-O-Al bond could be main- tained, but weakened by the formation of a hy- drogen bond, following which the silicate net- work would undergo merely modest deforma- tions (Kohn et al., 1989, 1998; Oglesby et al., 2001). Based on recent studies (Zeng et al., 1999, 2000; Schmidt et al., 2000, 2001), both mechanisms could occur in the melt. In fact, the significant effect of OH− formation on He solu- bility matches better the hypothesis of signifi- cant structural changes rather than small defor- mations, suggesting that noble gases can pro- vide valuable insights into the structure of alu- minosilicate melts. We are unable to evaluate whether or not Shibata et al.’s approach (1998) predicts the observed increase in the solubility of helium, because the uncertainties regarding the structural effect of water and its depolymer- izing role prevent us from obtaining a reliable estimation of the NBO/T ratio in water-bearing melts. On the other hand, the higher He solubil- ity is consistent with the calculated increase in the ionic porosity due to the addition of water. At higher concentrations of dissolved water, he- lium solubility seems to decrease with the in- crease of the molecular H2O/OH ratio, suggest- ing that the H2O molecules probably occupy free spaces in the silicate structure. 3.2. Pressure and temperature effects Equilibrium partitioning of noble gases be- tween vapor and silicate melt has been de- scribed using the following thermodynamic re- lation, which is based on the equality condition of the chemical potential of the gas species in vapor and liquid phases (3.6) where a and f refer to the activity of the noble gas in melts and its fugacity in vapor; P° and T° are the standard state conditions; ∆H 0P°, T is the enthalpy of dissolution of the noble gas at P° and T; V 0P, T is the molar volume of noble gas in its reference state at P and T, and R is the gas constant. Noble gas fugacity is normally com- puted, as a function of P and T, by using a mod- ified Redlich-Kwong equation of state (Hol- loway, 1977). The parameters a°/f °, V° and ∆H° must be calibrated using experimental data for each no- ble gas as a function of pressure, temperature and silicate melt composition. Use of eq. (3.6) and the consequent parameter estimation de- pends on the theoretical model adopted for the dissolution process. Based on the previous con- siderations concerning the physical dissolution of the noble gases, one can hypothesize that no- ble gas atoms dissolve into a fixed population of holes in the silicate melt network. One can also assume that the available sites for a given noble gas are identical, that their number is in- dependent of T and P, and that gas atoms do not influence each other reciprocally. This allows us to define an activity-composition relation, namely a=n/(N−n), where n is the number of dissolved atoms and N is the number of avail- able sites in the melt, referenced to a hypothet- ical reference state of pure noble gas. By fixing the temperature and assuming that the volume of any solubility site is unchanged after the dis- solution of a noble gas atom (namely, V°=0), eq. (3.6) and the previous definition of activity can be combined to obtain the well-known for- malism of the Langmuir adsorption isotherm (3.7) where KL(T) indicates that the equilibrium con- stant depends on temperature. The above for- ( ) / ( )K T n N n fL = - ( ) ln ln f a f a R H T T RT V P P 1 1 , , , , , , P T P T P T P T T P T = - - + - - ∆ P 0 0 c c c c c c c c c b l 659 A review on noble gas solubility and degassing in magmas malism has been widely applied to describe no- ble gas solubility in melts and glasses by con- sidering the number of sites as a parameter to be evaluated (Shelby, 1976; Roselieb et al., 1992; Carroll and Stolper, 1991; Walter et al., 2000). The very high-pressure experiments re- cently performed by Schmidt and Keppler (2002) reached a threshold Ar concentration at which the dissolved gas content did not in- crease with increasing Ar pressure in either sili- cic or tholeiitic melts (fig. 4). The authors sug- gested that the sites available to store the Ar atoms were fully occupied, therefore the N val- ue was accordingly constrained to 13.8×1020 sites/cm3 and the obtained Langmuir constant was able to fit the data very well. Approximat- ing the structure of the granitic liquid as a tridymite-like network, the number of intersti- tial holes can also be computed. When compar- ing such a number with the N value obtained by Schmidt and Keppler, the latter represents only 60% of the interstitial holes. This underlines the fact that only a fraction of the holes is available to accommodate Ar atoms. In accordance with its much more compact structure, the estimated N value in the tholeiitic melt was almost one or- der of magnitude lower (3.1×1020 sites/cm3). The number of holes available for the physical dissolution of Xe was also significantly lower than that for Ar, which is consistent with the larger size of the former. Although the Lang- muir isotherm was successfully applied, noble gas dissolution could violate one of its main as- sumptions, namely that gas atoms do not affect the structure of the solvent (V° ≠ 0). In fact, a re- cent X-ray absorption spectroscopy study in- vestigated the structural environment of dis- solved krypton in vitreous silica, revealing that Kr atoms are coordinated by oxygen atoms hav- ing a characteristic geometry (Wulf et al., 1999): thus, Kr atoms create or adjust to their own solubility sites. However, whether or not smaller gas atoms have a similar network-mod- ifier role must be verified. Models involving additional parameters which also contain information on the structure of the silicate network have been developed fol- lowing statistical-mechanical approaches (Dore- mus, 1966; Studt et al., 1970; Shackelford et al., 1972). In addition to the number of sites, these models include the vibrational frequency of the dissolved atom, which is treated as a harmonic oscillator in the structural cavity, and the binding energy to the interstitial hole. In silica glass, these two parameters depend upon the atomic size of the noble gas (lower frequency and more nega- tive binding energy for larger noble gases; Shel- by, 1976; Shackelford and Brown, 1980; Nakaya- ma and Shackelford, 1990). A sharp increase in the binding energy and decrease in the vibrational frequency were obtained along the join Na2O- SiO2 when moving from the silica-rich to the sodium-rich edge of the miscibility gap (Shelby, 1973). This is in agreement with both the stronger structural constraints and the necessity to spend energy in creating adequate holes in a less poly- merized solvent (see also below). A different theoretical approach in applying eq. (3.6) considers that the units of the silicate melt and noble gas atoms mix ideally. In this case, the activity of the dissolved gas is equal to its molar fraction when a Henrian reference state of pure noble gas is assumed (a hypothet- ical material composed by noble gas arranged like a silicate melt). The two parameters V° and ∆H° define the dependence of the solubility on P and T, respectively (see table I). Considering them as constants in the fitting procedure, the model is generally able to successfully repro- duce the experimental data in a large range of P and T conditions (Lux, 1987; White et al., 1989; Carroll and Stolper, 1991, 1993; Draper and Carroll, 1995; Shibata et al., 1998; Schmidt and Keppler, 2002). The calculated molar vol- umes for heavy noble gases (Ar, Xe) in melts are comparable to their respective co-volumes in MRK equations (the parameter named b), which reflects the atomic volume (White et al., 1989; Schmidt and Keppler, 2002). Slightly lower values than the co-volumes were deter- mined in glasses (Carroll and Stolper, 1991, 1993; Draper and Carroll, 1995), suggesting that the noble gases fit into pre-existing holes and only a small amount of new space must be created in the melt. Similar results were found for helium in glasses, which displays V° values lower than its co-volumes (Carroll and Stolper, 1991). The small V° values of the noble gases account for their Henrian behavior at low pres- sure (at least up to some hundreds of MPa), due 660 Antonio Paonita to the negligible weight of the V°(P−P°)/RT term in eq. (3.6). The Ar molar volumes com- puted for different melts (tholeiite to granitic melts; White et al. 1989, Schmidt and Keppler 2002) display a negligible dependence on the liquid composition, suggesting that a similar environment surrounds the dissolved atom in spite of the notable structural differences of the silicate structure. The heat of the solution (∆H°), in both glasses and melts, is generally small (see table I), and experimental data show that temperature dependence is close to that of the experimental uncertainties. Negative ∆H° values, namely de- creasing solubility at increasing temperature, have been estimated for silicate glasses (Shack- elford et al., 1972; Shelby, 1972a,b; Carroll and Stolper, 1991; Draper and Carroll, 1995) with compositions from albitic and rhyolitic up to fused silica (fig. 10). Increasing Na2O and K2O concentrations to above 20 mol% seem to invert this solubility-temperature relation (Shelby, 1973, 1974). On the other hand, a positive en- thalpy of solution and a positive direct tempera- ture dependence of solubility have normally been measured in basaltic melts (Lux, 1987; Jambon et al., 1986; Hayatsu and Waboso, 1985; see fig. 11), as well in Na2O-SiO2, Na2O-CaO- SiO2, Na2O-MgO-SiO2 melts (White et al., 1989; Shibata et al., 1988). All the described effects are normally more pronounced for heavier no- ble gases. The dissolution of a noble gas atom in the silicate melt can be usefully divided into two steps, each having its own enthalpy varia- tion: i) the creation of a hole to accommodate a gas atom in the melt structure, and ii) the trans- fer of the atom into the hole. The former step needs energy, namely its enthalpy variation is positive; instead, moving an atom from a gas into a more condensed phase involves the evo- lution of heat and a negative ∆H° sign. The whole enthalpy of solution is the sum of the two terms and its sign will depend on which contri- bution is dominant. Glasses and silica-rich melts have open structures, so we can expect that only a small contribution of energy is need- ed to create cavities inside. Moving atoms into the melt will be the energetically dominant process and, in view of the negative ∆H° sign of such a process, we can explain the negative dependence of solubility on temperature. Oppo- site considerations can be made to explain the direct solubility-temperature relationship of less polymerized, alkali-rich melts, where a few pre-existing large holes are available and ener- gy has to be spent to create or enlarge cavities. Fig. 10. Temperature effect on noble gas solubility in a basanite melt between 1200 and 1500°C at 0.1 MPa. Data are from Lux (1987) and uncertainties are in the order of the symbol size. Curves were comput- ed by using eq. (3.6) with parameters in table I. Sol- ubility slightly increases with the increase of temper- ature. Fig. 11. Dissolved Ar in a rhyolite glass at 120-130 MPa as a function of temperature (data from Carroll and Stolper, 1993). Curves were computed by using eq. (3.6) with parameters in table I. Solubility dis- plays a slightly negative relation with the increase of temperature. 661 A review on noble gas solubility and degassing in magmas 3.3. Size distribution of solubility sites Several researchers (Shackelford and Masa- ryk, 1978; Shackelford and Brown, 1980; Shack- elford, 1982; Nakayama and Shackelford, 1990; Carroll and Stolper, 1993) used the different sol- ubilities among the noble gases as a gas-probe to investigate the distribution of hole size in vitreous silica and natural silicate melts. The logarithm of noble gas solubility displays a linear relation with the atomic radius (fig. 5), meaning that hole sizes result as having a log-normal distribution. The mode shifts towards higher hole size with in- creasing ionic porosity, thus large holes become more common than small ones and the solubility of heavier noble gases increases more quickly with respect to that of light gases. In accordance with the above conclusions, a numerical simula- tion of the interstitial structure of silica glass con- firmed the log-normal distribution (Chan and El- liot, 1991). Nuccio and Paonita (2000) assumed that the m and n parameters in eq. (3.3) are composi- tion-independent for all noble gases, in accor- dance with what was experimentally verified for helium. In these conditions, their extended ionic porosity model for helium solubility in H2O-CO2 bearing melts is valid for other noble gases. They calculated that a low percent of wa- ter in a melt causes an exponential increase of the number of large holes, as a consequence of potential structural modifications of the melt. These conclusions need to be proved using sol- ubility data for noble gases having a large atomic radius. 4. Implications in modeling magma degassing Several elements in magmas are able to form chemical species which may be exsolved in va- por phases (H, C, O, S mainly, halogens, nitro- gen, noble gases etc.; Gerlach and Nordilie, 1975; Giggenbach, 1996). These elements can react with the liquid forming chemical bonds (i.e. OH groups), or they exist under the form of volatile atomic or molecular species dissolved in the silicate melt (i.e. molecular H2O, noble gas- es). Any dissolved volatile species (true or po- tential) is characterized by its own vapor pres- sure, which depends directly upon its concentra- tion in the melt. Magma reaches supersaturation when the sum of the vapor pressures of the indi- vidual volatiles becomes higher than the total pressure of the system, as a consequence of as- cent and depressurization of the melt (Carroll and Webster, 1994). In such conditions, the mag- ma is potentially able to exsolve a vapor phase. As H2O and CO2 are the most abundant volatiles in natural magmatic systems, the assessment of magma saturation has been modeled as a gas mixture of these two compounds (Papale, 1999), although sulfur could have significant effects when its concentrations reach several weight percent in the system (Moretti et al., 2003). When the vapor phase exsolves, minor volatiles are partitioned between the silicate liq- uid and the coexisting vapor, depending upon their solubilities. In fact, early vapor is strongly enriched in highly volatile species, which are consequently depleted in the melt (Giggenbach, 1996). The more soluble species are only able to exsolve into the gas phase during the late stages of degassing. As a consequence, less sol- Fig. 12. Fraction of noble gases remaining in a tholeiitic melt with respect to the amount of exsolved vapor during a closed system degassing. Vapor/melt ratio is the volume of gas phase exsolved per volume unit of melt. Adopted equations for calculations were those by Giggenbach (1996), as well as the H2O and CO2 solubilities. CO2 solubility causes CO2 to degas before light noble gases, whereas H2O is much more soluble and is exsolved during late degassing. 662 Antonio Paonita uble, heavy noble gases (Ar, Kr, Xe) are de- gassed earlier, whereas He and Ne degassing occurs later (fig. 12). This solubility-controlled fractionation of the noble gases is much more significant if in- finitesimal aliquots of volatiles are exsolved in equilibrium from the liquid phase and immedi- ately removed, thereby creating open system degassing conditions. However, this behavior is also evident during closed system degassing processes, during which the exsolved vapor stays in contact with the melt throughout the degassing process (Carroll and Webster, 1994). Noble gas fractionation in open system de- gassing can be modeled by applying the typical Rayleigh distillation process. By considering two noble gases i and j, the extent of fractiona- tion between two gases is related to the residual fraction of one of the two volatiles (4.1) where C is the volatile concentration in melt (m); subscript «0» refers to initial conditions (undegassed parental melt) and K is the Henry’s constant. Calculations by eq. (4.1) have been ap- plied to ratios between noble gases as a function of the total content of one noble gas in basalts (Taylor, 1986; Sarda and Graham, 1990; Carroll and Webster, 1994; Marty, 1995; Giggenbach, 1996; Burnard, 1999; Moreira and Sarda, 2000). By writing two Rayleigh equations for a ternary system, such as He-Ar-CO2, and coupling them, it can be shown that a plot of ln(Ar/CO2) versus ln(Ar/He) yields a straight line with a slope de- pendent upon the solubility ratios among the gases and an intercept dependent upon the Ar/CO2 and Ar/He ratios in the pristine unde- gassed melt (Burnard, 1999, 2001). In a closed system, vesicles of an exsolved vapor remain in the melt, so the total amount of gas does not vary. Models have been developed that take into account the mass balance among the total amount of noble gas, its dissolved con- centration and the volume of exsolved vapor per mass unit of melt (Jambon et al., 1986; Sar- da and Moriera, 2002). By coupling Henry’s Law to the mass balance, the fractionation be- C C C C C C ,j i j i i i K K 00 1 j i = - c c d m m n tween noble gases can be related to the vesicu- larity of the melt (the volume of vapor phase exsolved per volume unit of melt) (4.2) where T0 is equal to 273 K; Te is the tempera- ture of the gas-melt equilibrium; ρ is the densi- ty of the silicate melt (assumed constant), and V* is the melt vesicularity. The above models have been frequently ap- plied to study degassing on spatial scales from local to planetary. For example, investigations were performed on the gas phase trapped in vesi- cles and dissolved in glasses recovered from Middle-Oceanic Ridge Basalts (MORB) and Oceanic Island Basalts (OIB) (Jambon et al., 1986; Taylor, 1986; Sarda and Graham, 1990; Carroll and Webster, 1994; Giggenbach, 1996; Burnard, 1999, 2001; Honda and Patterson, 1999; Moreira and Sarda, 2000; Sarda and Mori- era, 2002). Furthermore, constraints were ob- tained on the outgassing history of the mantle and its primordial composition. The higher He/Ar ratios in MORB glasses with respect to their production/accumulation rate in the mantle (due to K and U decay), have normally been ex- plained by degassing processes (Jambon et al., 1985), because Ar is preferentially degassed in respect of He. On this basis, Moriera and Sarda (2000) focused on the different degassing modalities that occur from MORB and from mantle plumes, the former degassing as a closed system, and the latter as an open-system Rayleigh-type process. However, some ques- tions have been raised about this «degassing» hypothesis, namely that higher He/Ar ratios oc- cur at higher He contents (Honda and Patterson, 1999) and in basalts erupted at greater depths (Fisher, 1997). Solubility controlled degassing does not explain these constraints, whereas mechanisms of helium diffusion from mantle minerals or crystals in basalt melts have been hy- pothesized (Matsuda and Marty, 1995; Fisher, 1997; Honda and Patterson, 1999). In fact, Mor- eira and Sarda (2000) and Sarda and Moreira (2002) have suggested that the geochemical pat- tern of noble gas abundances in MORB is con- sistent with equilibrium vesiculation followed by / ( ) / ( ) C C C C K K V T T K V T T K , , * * j i j i i j e e i j 0 0 0 0 # #= + + t t 663 A review on noble gas solubility and degassing in magmas various episodes of bubble loss, in accordance with the early model by Sarda and Graham (1990). However, Sarda and Moreira (2002) have also suggested that helium solubility should be two times higher than that experimentally measured at 0.1 MPa, and that Ar and Ne solu- bilities should be lower by factors of 3.5 and 7, respectively. This opposite variation in solubili- ties makes their hypothesis about a pressure ef- fect on solubility rather unlikely. Burnard (1999) moves criticisms on the experimental values of solubilities too, providing evidence that the no- ble gas/CO2 solubility ratio has to be higher than experimentally measured. Nevertheless, a new view on noble gas data from MORB takes into account the size of the vesicles, as small bubbles should represent the previous extent of the de- gassing with respect to large vesicles (Burnard, 2001). Sequential crushing of basalts allows for the analysis of gas from progressively smaller bubbles, which can be used to determine the tra- jectory of the degassing. Preliminary data seem to reconcile with predicted solubility-controlled fractionation studies (Burnard, 2001). The foregoing discussion stresses how noble gas abundances from natural samples are used as a reference point and as a check of experi- mental results on noble gas solubilities. It also highlights the wide use of the noble gas solubil- ity data in Earth Sciences problems and the strong inferences they provide. Considering the fact that the experimental conditions do not match the geological ones, degassing models that use such data must be applied very careful- ly. Equations (4.1) and (4.2) describe how noble gases are partitioned in a vapor while they are progressively exsolved, regardless of: i) what type of volatiles dominate the vapor and if they behave non-ideally in such a gas mixture, ii) changes in concentration and, above all, in solu- bility due to the exsolution of the main volatiles. Application to MORB degassing could disre- gard the second conditions, at least for a first ap- proximation, because the amounts of H2O and CO2 in these melts are lower than 1 wt% at least for P < 1 GPa. However, MORB degassing should take into account the noble gas fugacity in CO2-dominated vapor. Bottinga and Javoy (1990) discussed the exsolution of a non-ideal H2O-CO2-He-Ar mixture from bas-altic melts, although the H2O and CO2 effects on noble gas solubility were not considered. Nuccio and Paonita (2001) started from the already dis- cussed model for noble gas solubility in H2O- CO2 bearing melts (Nuccio and Paonita, 2000; see Section 3.1). The authors obtained an equa- tion that allowed them to calculate the concen- tration of noble gas in an exsolved H2O-CO2 mixed vapor as a function of the dissolved con- centration of H2O, CO2 and noble gas in melt and their total amounts in the vapor+melt sys- tem. Thus, the authors numerically simulated the degassing of multicomponent H2O-CO2-no- Fig. 13. He/Ar fractionation during closed system degassing of basaltic melt at 1200°C. Initial He/Ar ratio was 3. Thick grey lines were computed by us- ing the model of Jambon et al. (1986) (eq. (4.2) in the text), assuming noble gas solubilities were not af- fected by the presence of the main volatiles (K = con- stant) and ideal gas behavior. Thin black lines were computed on the basis of the model by Nuccio and Paonita (2001), which accounts for these effects. In this model, the initial conditions were 85 mol% CO2 in the H2O-CO2 mixed vapor (noble gases are com- puted as minor species) at 700 MPa pressure. Black thin lines for melt and vesicles approach each other at low pressure, due to the effect of water that makes the solubilities of noble gases more similar. Dashed curves were computed for an open system (grey line by Jambon et al., 1986 model black line by Nuccio and Paonita, 2001). 664 Antonio Paonita ble gas mixtures from depressurizing magmas. As a result, remarkable differences were high- lighted between the comprehensive approach and simpler models (figs. 13-15). In spite of the constant ratio between He/Ar in melt and vesi- cles computed by the simpler model (namely, (He/Ar)gas / (He/Ar)melt = KHe / KAr), curves from the more complex model come closer at low pressure. In fact, taking into account the fact that high pressure causes very different solubility among noble gases and that dissolved water has the opposite effect, a few wt% of dissolved wa- ter at low pressure are able to significantly re- duce the solubility difference between He and Ar (fig. 13 and see Section 3.1). As a conse- quence, noble gases may fractionate less than predicted by simpler calculations, depending on the composition of the magma and on the main volatiles it contains. The application of the mod- el to very different contexts, from silicic to basaltic volcanoes (Nuccio and Paonita 2001; Caracausi et al., 2003), shows how the inert gas content of outgassed magmatic volatiles may produce information and constraints on deep magmatic conditions: magma pressure and depth over time, degassing and volatile content, melt composition and temperature. The current research aims at developing approaches to de- scribe saturation surfaces in more complex sys- tems, including sulfur and halogens. On this ba- sis, the solubility of noble gases in such systems progressively more similar to true geological conditions, will have to be theoretically and ex- perimentally investigated, allowing for the de- velopment of more refined degassing models. Finally, vapor-melt equilibrium for noble gases has been assumed in all the discussed models. Experimental observations and theoret- ical modeling have been performed to investi- gate the growth of steam bubbles in silicate liq- uids (Bottinga and Javoy, 1990; Proussevich and Sahagian, 1996, 1998; Gardner et al., 1999, 2000). Bubbles seem to maintain chemical equilibrium with the surrounding melt even in those silicic magmas whose ascent rates are Fig. 14. He/Ar fractionation during closed system degassing of rhyolitic melt at 1000°C. Curves as in fig. 13. In the calculation using the model by Nuccio and Paonita (2001), initial conditions were 40 mol% CO2 in H2O-CO2 mixed vapor (noble gases are com- puted as minor species) at 500 MPa pressure. The higher amount of dissolved water (6 wt%) in respect of the basaltic melt amplifies the discussed water ef- fect, causing very low He/Ar fractionation. Fig. 15. He/Ne fractionation during closed system degassing of basaltic melt at 1200°C. Curves and conditions as in fig. 13. The behavior is qualitatively similar to that discussed in the case of argon, al- though the fractionation is more modest. 665 A review on noble gas solubility and degassing in magmas three orders of magnitude higher than their common pre-eruptive values (∼ 0.01 m/s), whereas non-equilibrium exsolution could oc- cur during explosive rising. Nevertheless, there is a total lack both of experimentation and mod- els for the investigation of possible disequilib- rium fractionation among noble gases. In addi- tion, noble gas diffusion in melts containing H2O and/or CO2, which is crucial to describe such processes, is also poorly known. 5. Upcoming research Although our knowledge of noble gas solu- bility in magmas has significantly improved over the last twenty years, a lot of work has yet to be accomplished from both experimental and theo- retical points of view. Magma is a multicompo- nent system carrying several volatile compo- nents in its gas phase; nevertheless our studies on the solubility of noble gases in such complex volatile-bearing systems are still in their early stages. Few experimental data are available re- garding helium solubility in silicate melts coex- isting with H2O-CO2 mixed vapors. Although these studies are not comprehensive enough to fully describe and understand the effect of the main volatile components on noble gas solubili- ty, they do highlight that it cannot be neglected. The dissolution of noble gases is strictly related to the availability of solubility sites in the silicate network and, due to the significant structural changes in melts caused by the main volatiles, H2O and CO2 may well have a very significant control over noble gas solubility. On this basis, sulfur and chlorine (perhaps fluorine as well) may exhibit similar effects, however we lack critical information to evalu- ate them. As far as the gas phase is concerned, there are equations of state that can calculate the activity-composition relations for Ar in H- C-S-O mixed vapor in a wide range of pressure and temperature conditions (Belenosko et al., 1992). More recently, equations have become available for noble gases in H-C-S-O-halogens (Churakov and Gottschalk, 2003). Neverthe- less, the activities of noble gas in S-Cl bearing silicate melts are totally unknown, and there is a great need for experimentation in this field of research. In addition, the inadequate knowledge of the structural effect of sulfur and chlorine does not allow us to evaluate their influence on the distribution of the hole-size in melts and, as a consequence, on noble gas solubility. Regard- less of the analytical method, the experimental determination of noble gas solubility in com- plex systems requires the charging of known amounts of volatiles inside capsules, in order to obtain vapor phases having similar composi- tions to those of magmatic ones. In view of this, the method of loading noble gases by gas-bear- ing glasses (see Section 2) should be broadly applicable to the investigation of the solubility of all noble gases (even He and Ne) in melts up to very high pressures. Although the mass bal- ance between the total and dissolved noble gas allows us to calculate the concentration of no- ble gas in vapor and then the partition coeffi- cient, the ability to measure the content of no- ble gas in both the vapor and liquid phase be- comes ever more important. In this respect, coupling in situ extraction techniques and mass spectrometry seems the most powerful tool, es- pecially in investigating the solubility of light noble gases. In addition, X-ray absorption techniques have started to provide insights into how noble gas atoms are accommodated in the silicate melt structure and the short-range order. Although these data do not produce a direct measurement of solubility, they are still able to constrain the parameters used in the statistical-mechanics sol- ubility models. Such theoretical approaches give detailed information and explanations of the mechanism of noble gas dissolution in melt and thus they need to be further improved. Addition- al information could also come, in the future, from computer simulations of noble gas-bearing melts by using Molecular Dynamics (MD) meth- ods, although a lot of work must be done to de- fine pair potentials between the noble gas and the main particles of the melt by fitting experimen- tal data or by performing ab-initio calculations. For this work, the expertise of mineralogists and materials scientists, interested in the consider- able inferences of noble gas solubility on the structure of silicate glass and melt, could be very helpful. This co-operation should aim at devel- oping several aspects having technological and 666 Antonio Paonita commercial interests and involving the interac- tion of gas atoms with rigid glass materials. For instance, noble gases (He, Ar) are normally used in the growth of silicon-made materials for use in semiconductor devices, or with the aim of pre- venting gas inclusion formation during the man- ufacture of glass and ceramics. Ceramics are al- so porous enough to allow for noble gas perme- ation, allowing Magnetic Resonance Imaging (MRI) using laser-polarized noble-gas-testing of the structural integrity of the material. Large dif- ferences in the solubility and diffusion of the gases can be used to separate the components of gas mixtures. Finally, an interesting field to be developed in the future involves the possible effects of dis- equilibrium during noble gas degassing from magmas. In their recent paper, Caracausi et al. (2003) show that the helium isotopic composi- tion monitored in volcanic emissions at Mt. Etna displays variations over time, which cannot be referred to shallow secondary processes. Kinetic fractionation of helium isotopes during diffusion in growing bubbles has been suggested to ex- plain this type of signal (Nuccio and Valenza, 1998; Caracausi et al., 2003). In fact, light 3He atoms could diffuse at higher rates in growing bubbles with respect to heavier 4He atoms. In the same way, the ratios between noble gases both in volcanic gases and fluid inclusions could be af- fected by a similar disequilibrium process, espe- cially when taking into consideration that they have similar dissolution mechanisms and that there is a big difference between their atomic masses. Experiments on this topic are not simple and theoretical calculations, to be based on noble gas diffusion rates, require adequate knowledge of the effects of the main volatiles on noble gas diffusivity in silicate melts. A lot of effort still has to be made in this direction, which however seems very promising in its application to mag- matic systems. Acknowledgements The author wishes to thank P.M. Nuccio for useful comments and suggestions. The review by D. Baker and C. Romano greatly improved the manuscript. REFERENCES BELENOSHKO, A. and S.K. SAXENA (1992): A unified equa- tion of state for fluids of C-H-O-N-S-Ar composition and their mixtures up to very high temperatures and pressures, Geochim. Cosmochim. Acta, 56, 3611-3626. BOETTCHER, S.L., Q. GUO and A. MONTANA (1989): A sim- ple device for loading gases in high-pressure experi- ments, Am. Mineral., 74, 1383-1384. BOTTINGA, Y. and M. JAVOY (1990): MORB degassing: bub- ble growth and ascent, Chem. Geol., 81, 255-270. BRAWER, S.A. and W.B. WHITE (1975): Raman spec- troscopy investigation of the structure of silicate glass- es, I. The binary silicate glasses, J. Chem. Phys., 63, 2421-2432. BROADHURST, C.L., M.J. DRAKE, B.E. HAGEE and T.J. BERNATOWICZ (1992): Solubility and partitioning of Ne, Ar, Kr, and Xe in minerals and synthetic basalt melts, Geochim. Cosmochim. Acta, 56, 709-723. BROOKER, R.A., J.-A. WARTHO, M.R. CARROLL, S.P. KEL- LEY and D.S. DRAPER (1998): Prelimirary UVLAMP de- terminations of argon partition coefficients for olivine and clinopyroxene grown from silicate melts, Chem. Geol., 147, 185-200. BURNARD, P. (1999): The bubble-by-bubble volatile evolu- tion of two mid-ocean ridge basalts, Earth Planet. Sci. Lett., 174, 199-211. BURNARD, P. (2001): Correction for volatile fractionation in ascending magmas: noble gas abundances in primary mantle melts, Geochim. Cosmochim. Acta, 65, 2605- 2614. CARACAUSI, A., F. ITALIANO, P.M. NUCCIO, A. PAONITA and A. RIZZO (2003): Evidence of deep magma degassing and ascent by geochemistry of peripheral gas emis- sions at Mt. Etna (Italy): assessment of the magmatic reservoir pressure, J. Geophys. Res., 108, 2463-2484. CARROLL, M.R. and D.S. DRAPER (1994): Noble gases as trace elements in magmatic processes, Chem. Geol., 117, 37-56. CARROLL, M.R. and E.M. STOLPER (1991): Argon solubility and diffusion in silica glass: implications for the solu- tion behavior of molecular gases, Geochim. Cos- mochim. Acta, 55, 211-225. CARROLL, M.R. and E.M. STOLPER (1993): Noble gas solu- bilities in silicate melts and glasses: new experimental results for Ar and the relationship between solubility and ionic porosity, Geochim. Cosmochim. Acta, 57, 5039-5051. CARROLL, M.R. and J.D. WEBSTER (1994): Solubilities of sulfur, noble gases, nitrogen, chlorine, and fluorine in magmas, in Volatiles in Magmas, edited by M.R. CAR- ROLL and J.R. HOLLOWAY, Rev. Mineral., 30, 231-279. CARROLL, M.R., S.R. SUTTON, M.L. RIVERS and D.S. WOOLUM (1993): A experimental study of krypton dif- fusion and solubility in silicic glasses, Chem. Geol., 109, 9-28. CHAMORRO-PEREZ, E., P. GILLET and A. JAMBON (1996): Ar- gon solubility in silicate melts at very high pressures. Experimental set-up and preliminary results for silica and anorthite melts, Earth Planet. Sci. Lett., 145, 97- 107. CHAMORRO-PEREZ, E., P. GILLET, A. JAMBON, J. BARDO and P. MCMILLAN (1998): Low argon solubility in silicate 667 A review on noble gas solubility and degassing in magmas melts at high pressure, Nature, 393, 352-355. CHAMORRO, E.M., R.A. BROOKER, J.A. WARTHO, B.J. WOOD, S.P. KELLEY and J.D. BLUNDY (2002): Ar and K partitioning between clinopyroxene-silicate at 8 GPa, Geochim. Cosmochim. Acta, 66, 507-519. CHAN, S.L. and S.R. ELLIOT (1991): Theoretical study of the interstice statistics of the oxygen sublattice in vitre- ous SiO2, Phys. Rev. B, 43, 4423-4432. CHENNAOUI-AOUDJEANE, H. and A. JAMBON (1990): He sol- ubility in silicate glasses at 250°C: a model for calcu- lation, Eur. J. Mineral., 2, 539-545. CHURAKOV, S.V. and M. GOTTSCHALK (2003): Perturbation theory based equation of state for polar molecular flu- ids, II. Fluid mixtures, Geochim. Cosmochim. Acta, 67, 2415-2425. DE SANTIS, R., G.J.F. BREEDVELD and J.M. PRAUSNITZ (1974): Thermodynamic properties of aqueous gas mixtures at advanced pressures, Ind. Eng. Chem. Process. Des. Develop., 13, 374-377. DIXON, J.E., E.M. STOLPER and J.R. HOLLOWAY (1995): An experimental study of water and carbon dioxide solu- bilities in Mid-Ocean Ridge basaltic liquids, Part I. Calibration and Solubility models, J. Petrol., 36 (6), 1607-1631. DOREMUS, H. (1966): Physical solubility of gases in fused silica, J. Am. Ceram. Soc., 49, 461-462. DRAPER, D.S. and M.R. CARROLL (1995): Argon diffusion and solubility in silicic glasses exposed to an Ar-He mixtures, Earth. Planet. Sci. Lett., 132, 15-24. FRANK, R.C., D.E. SWETS and R.W. LEE (1961): Diffusion of neon isotopes in fused quartz, J. Chem. Phys., 35, 1451-1459. FISHER, D.E. (1970): Heavy rare gas in a Pacific seamount, Earth Planet. Sci. Lett., 9, 331-335. FISHER, D.E. (1997): Helium, argon and xenon in crushed and melted MORB, Geochim. Cosmochim. Acta, 61, 3003-3012. GARDNER, J.E., M. HILTON and M.R. CARROLL (1999): Exper- imental constrains on degassing of magma: isothermal bubble growth during continuous decompression from high pressure, Earth Planet. Sci. Lett., 168, 201-218. GARDNER, J.E., M. HILTON and M.R. CARROLL (2000): Bub- ble growth in highly viscous silicate melts during con- tinuous decompression from high pressure, Geochim. Cosmochim. Acta, 64, 1473-1483. GERLACH, T.M. and B.E. NORDILIE (1975): The C-O-H-S gaseous system, Part I. Composition limits and trends in basaltic glasses, Am. J. Sci., 275, 353-376. GIGGENBACH, W.F. (1996): Chemical composition of volcanic gases, in Monitoring and Mitigation of Volcanic Hazards, edited by R. SCARPA and R. TILLING (Springer), 221-256. HAYATSU, A. and C.E. WABOSO (1985): The solubility of rare gases in silicate melts and implications for K-Ar dating, in Terrestrial Noble Gases, edited by F.A. PO- DOSEK, Chem. Geol., 52, 97-102. HOLLOWAY, J.R. (1977): Fugacity and activity of molecular species in supercritical fluids, in Thermodynamics in Geology, edited by D.G. FRASER (Reidel Pub. Comp., Dordrecht-Holland, Boston). HONDA, M. and D.B. PATTERSON (1999): Systematic ele- mental fractionation of mantle-derived helium, neon and argon in mid-ocean ridge glasses, Geochim. Cos- mochim. Acta, 63, 2863-2874. IHINGER, P.D., Y. ZHANG and E.M. STOLPER (1999): The spe- ciation of dissolved water in rhyolitic melt, Geochim. Cosmochim. Acta, 63, 3567-3578. JAMBON, A. (1987): He solubility in silicate melts: a tenta- tive model of calculation, Chem. Geol., 62, 131-136. JAMBON, A., H.W. WEBER and F. BEGEMAN (1985): Helium and argon from an Atlantic MORB glass: concentra- tion, distribution and isotopic composition, Earth Planet. Sci. Lett., 73, 255-267. JAMBON, A., H.W. WEBER and O. BRAUN (1986): Solubility of He, Ne, Ar, Kr, Xe in a basalt melt in the range of 1250-1600°C: geochemical implications, Geochim. Cosmochim. Acta, 50, 401-408. KELLEY, S.P., N.O. ARNAUD and G. TURNER (1994): High spatial resolution 40Ar/39Ar investigation of quartz, bi- otite and plagioclase using a new ultra-violet laser probe extraction technique, Geochim. Cosmochim. Ac- ta, 58, 3519-3525. KING, P.L. and J.R. HOLLOWAY (2002): CO2 solubility and speciation in intermediate (andesitic) melts: the role of H2O and composition, Geochim. Cosmochim. Acta, 66, 1627-1640. KIRSTEN, T. (1968): Incorporation of rare gas in solidifying enstatite melts, J. Geophys. Res., 73, 2807-2810. KOHN, S.C., R. DUPREE and M.E. SMITH (1989): A multinu- clear magnetic resonance study of the structure of hy- drous albite glasses, Geochim. Cosmochim. Acta, 53, 2925-2935. KOHN, S.C., M.E. SMITH, P.J. DIRKEN, E.R.H. VAN ECK, A.P.M. KENTGENS and R. DUPREE (1998): Sodium en- vironment in hydrous and dry albite glasses: improved 23Na solid state NMR data and their implication for wa- ter dissolution mechanisms, Geochim. Cosmochim. Ac- ta, 62, 79-87. LANGE, R.M. and I.S.E. CHARMICHAEL (1987): Densities of Na2O-K2O-CaO-MgO-FeO-Fe2O3-Al2O3-TiO2-SiO2 liquids: new measurements and derived partial molar properties, Geochim. Cosmochim. Acta, 51, 2931-2946. LUX, G. (1987): The behavior of noble gases in silicate liq- uids: solution, diffusion, bubbles, and surface effects, with applications to natural samples, Geochim. Cos- mochim. Acta, 51, 1549-1560. MARTY, B. (1995): Nitrogen content of the mantle inferred from N2-Ar correlation in oceanic basalts, Nature, 377, 326-329. MATSUDA, J. and B. MARTY (1995): The 40Ar/36Ar ratio of the undepleted mantle: a reevaluation, Geophys. Res. Lett., 22, 1937-1940. MCMILLAN, P.F. (1994): Water solubility and speciation models, in Volatiles in Magmas, edited by M.R. CAR- ROLL and J.R. HOLLOWAY, Rev. Mineral., 30, 131-156. MONTANA, A., Q. GUO, S.L. BOETTCHER, B.S. WHITE and M. BREARLEY (1993): Xe and Ar in high-pressure sili- cate liquids, Am. Mineral., 78, 1135-1142. MOREIRA, M. and P. SARDA (2000): Noble gas constraints on degassing processes, Earth Planet. Sci. Lett., 176, 375-386. MORETTI, R., P. PAPALE and G. OTTONELLO (2003): A mod- el for the saturation of C-O-H-S fluids in silicate melts, in Volcanic Degassing, edited by C. OPPENHEIMER, D.M. PYLE and J. BARCLAY, Geol. Soc. Spec. Publ. 213, 81-101. MYSEN, B.O., D. VIRGO and F.A. SEIFERT (1985): Relation- 668 Antonio Paonita ships between properties and structure of aluminosili- cate melts, Am. Mineral., 70, 88-105. NAKAYAMA, G.S. and J.F. SHACKELFORD (1990): Solubility and diffusivity of argon in vitreous silica, J. Non-Cryst. Solids, 126, 249-260. NOWAK, M. and H. BEHRENS (1995): The speciation of wa- ter in haplogranitic glasses and melts determined by in situ near-infrared spectroscopy, Geochim. Cosmochim. Acta, 59, 3445-3450. NOWAK, M. and H. BEHRENS (2001): Water in rhyolitic magmas: getting a grip on a slippery problem, Earth Planet. Sci. Lett., 184, 515-522. NUCCIO, P.M. and A. PAONITA (2000): Investigation of the noble gas solubility in H2O-CO2 bearing silicate liquids at moderate pressure, II. The Extended Ionic Porosity (EIP) model, Earth Planet. Sci. Lett., 183, 499-512. NUCCIO, P.M. and A. PAONITA (2001): Magmatic degassing of multicomponent vapors and assessment of magma depth: application to Vulcano Island (Italy), Earth Planet. Sci. Lett., 193, 467-481. NUCCIO, P.M. and M. VALENZA (1998): Magma degassing and geochemical detection of its ascent, in Water-Rock Interaction, edited by G.B. AREHART and J.R. HULSTON (Balkema), 475-478. OGLESBY, J.V., S. KROEKER and J.F. STEBBINS (2001): Potas- sium hydrogen disilicate: a possible model compound for 17O NMR spectra of hydrous silicate glasses, Am. Mineral., 86, 341-347. PAONITA, A., G. GIGLI, D. GOZZI, P.M. NUCCIO and R. TRIG- ILA (2000): Investigation of He solubility in H2O-CO2 bearing silicate liquids at moderate pressure: an exper- imental method, Earth Planet. Sci. Lett., 181, 595-604. PAPALE, P. (1999): Modeling of the solubility of a two com- ponent H2O+CO2 fluid in silicate liquids, Am. Miner- al., 84 (4), 477-492. PROUSSEVITCH, A.A. and D.L. SAHAGIAN (1996): Dynamics of coupled diffusive and decompressive bubble growth in magmatic systems, J. Geophys. Res., 101, 17447-17456. PROUSSEVITCH, A.A. and D.L. SAHAGIAN (1998): Dynamics and energetics of bubble growth in magmas: analytical formulation and numerical modeling, J. Geophys. Res., 103, 18223-18251. ROSELIEB, K., W. RAMMENSEE, H. BUTTNER and M. ROSEN- HAUER (1992): Solubility and diffusion of noble gases in vitreous albite, Chem. Geol., 96, 241-266. ROSELIEB, K., W. RAMMENSEE, H. BUTTNER and M. ROSEN- HAUER (1995): Diffusion of noble gases in melts of the system SiO2-NaAlSi2O6, Chem. Geol., 120, 1-13. SARDA, P. and D. GRAHAM (1990): Mid-oceanic ridge pop- ping rocks: implications for degassing at ridge crests, Earth Planet. Sci. Lett., 97, 268-289. SARDA, P. and M. MORIEIRA (2002): Vesiculation and vesi- cle loss in mid-ocean ridge basalt glasses: He, Ne, Ar elemental fractionation and pressure influence, Geo- chim. Cosmochim. Acta, 66, 1449-1458. SCHMIDT, B.C. and H. KEPPLER (2002): Experimental evi- dence for high noble gas solubilities in silicate melts under pressures, Earth Planet. Sci. Lett., 195, 277-290. SCHMIDT, B.C., T. RIEMER, S.C. KOHN, H. BEHRENS and R. DUPREE (2000): Different water solubility mechanisms in hydrous glasses along the quartz-albite join. Evi- dence from NMR spectroscopy, Geochim. Cosmochim. Acta, 64, 513-526. SCHMIDT, B.C., T. RIEMER, S.C. KOHN, F. HOLTZ AND R. DUPREE (2001): Structural implications of water disso- lution in haplogranitic glasses from NMR spectroscopy: influence of total water content and mixed alkali effect, Geochim. Cosmochim. Acta, 65, 2949-2964. SHACKELFORD, J.F. (1982): A gas probe analysis of structure in bulk and surface layers of vitreous silica, J. Non- Cryst. Solids, 49, 299-307. SHACKELFORD, J.F. (1999): Gas solubility in glasses – Prin- ciples and structural implications, J. Non-Cryst. Solids, 253, 231-241. SHACKELFORD, J.F. and B.D. BROWN (1980): A gas probe analysis of structure in silicate glasses, J. Am. Ceram. Soc., 63, 562-565. SHACKELFORD, J.F. and J.S. MASARYK (1978): The intersti- tial structure of vitreous silica, J. Non-Cryst. Solids, 30, 127-139. SHACKELFORD, J.F., P.L. STUDT and R.M. FULRATH (1972): Solubility of gases in glass II. He, Ne, and H2 in fused silica, J. Appl. Phys., 43, 1619-1626. SHANNON, R.D. and C.T. PREWITT (1969): Effective ionic radii in oxides and fluorides, Acta Cryst., 25, 925-946. SHELBY, J.E. (1972a): Helium migration in natural and syn- thetic vitreous silica, J. Am. Ceram. Soc., 55, 61-64. SHELBY, J.E. (1972b): Neon migration in vitreous silica, J. Am. Ceram. Soc., 55, 167-170. SHELBY, J.E. (1973): Effects of phase separation on helium migration in sodium silicate glasses, J. Am. Ceram. Soc., 56, 263-266. SHELBY, J.E. (1974): Helium diffusion and solubility in K2O-SiO2 glasses, J. Am. Ceram. Soc., 57, 236-263. SHELBY, J.E. (1976): Pressure dependence of helium and neon in vitreous silica, J. Appl. Phys., 47, 135-139. SHEN, A. and H. KEPPLER (1995): Infrared spectroscopy of hydrous silicate melts to 1000°C and 10 kbar: direct observation of H2O speciation in a diamond-anvil cell, Am. Mineral., 80, 1335-1338. SHIBATA, T., E. TAKAHASHI and J. MATSUDA (1996): Noble gas solubility in binary CaO-SiO2 system, Geophys. Res. Lett., 23, 3139-3142. SHIBATA, T., E. TAKAHASHI and J. MATSUDA (1998): Solubil- ity of neon, argon, krypton and xenon in binary and ter- nary silicate system: a new view on noble gas solubili- ty, Geochim. Cosmochim. Acta, 62, 1241-1253. STUDT, P.L., J.F. SHACKELFORD and R.M. FULRATH (1970): Solubility of gases in glass – A monoatomic model, J. Appl. Phys., 44, 2777-2780. SYKES, D. and J. KUBICKI (1993): A model for H2O solubil- ity mechanisms in albite melts from infrared spec- troscopy and molecular orbital calculations, Geochim. Cosmochim. Acta, 57, 1039-1052. TAYLOR, B.E. (1986): Magmatic volatiles: isotopic varia- tions of C, H, and S, Rev. Mineral., 16, 185-225. VIRGO, D., B.O. MYSEN and I. KUSHIRO (1980): Anionic con- stitution of 1-atmosphere silicate melts: implications for the structure of igneous melts, Science, 208, 1371-1373. WALTER, H., K. ROSELIEB, H. BUTTNER and M. ROSENHAUER (2000): Pressure dependence of the solubility of Ar and Kr in melts of the system SiO2-NaAlSi2O6, Am. Miner- al., 85, 1117-1127. WHITE, B.S., M. BREARLEY and A. MONTANA (1989): Solu- bility of Argon in silicate liquids at high pressures, Am. Mineral., 74, 513-529. 669 A review on noble gas solubility and degassing in magmas WITHERS, A.C., Y. ZHANG and H. BEHRENS (1999): Recon- ciliation of experimental results on H2O speciation in rhyolitic glass using in situ and quenching techniques, Earth Planet. Sci. Lett., 173, 343-349. WULF, R., G. CALAS, A. RAMOS, H. BUTTNER, K. ROSELIEB and M. ROSENHAUER (1999): Structural environment of krypton dissolved in vitreous silica, Am. Mineral., 84, 1461-1463. ZENG, Q., H. NEKVASIL and C.P. GREY (1999): Proton envi- ronments in hydrous aluminosilicate glasses: a 1H MAS, 1H/27Al, and 1H/23Na TRAPDOR NMR study, J. Phys. Chem. B, 103, 7406-7415. ZENG, Q., H. NEKVASIL and C.P. GREY (2000): In support of a depolymeryzation model for water in sodium alumi- nosilicate glasses: information from NMR spec- troscopy, Geochim. Cosmochim. Acta, 64, 883-896.