TY - JOUR

T1 - Solubility of H2O in rhyolitic melts at low pressures and a new empirical model for mixed H2O-CO2 solubility in rhyolitic melts

AU - Liu, Yang

AU - Zhang, Youxue

AU - Behrens, Harald

N1 - Funding Information: This study was supported by NSF grants EAR-9972937, EAR-0125506, EAR-0228752. We thank Otto Diedrich for preparation of polished samples after experiments, Sha Zhu and Daming Wang with the help of TEM, and P. Papale for providing calculated solubility using his solubility model. Discussion with Eric Essene was of great help. We would also like to thank M. Rutherford, D. Sahagian, and an anonymous reviewer for their constructive reviews that shaped the final manuscript.

PY - 2005/5/1

Y1 - 2005/5/1

N2 - Experiments of H2O solubility in synthetic haplogranitic and natural rhyolitic melts were conducted at 700-1200 °C and 0.098-25 MPa. Attainment of equilibrium was proven by reversals using starting materials with different initial H2O contents. Dissolved H2O contents were determined using Fourier transform infrared spectroscopy. When temperature decreases from 1000 to 700 °C, the solubility of H2O increases from 0.100 to 0.124 wt.% at 0.098 MPa, from 0.99 to 1.36 wt.% at 11 MPa, and from 1.46 to 2.17 wt.% at 25 MPa. At 6 MPa, the solubility of H2O increases from 0.63 to 0.77 wt.% from 1200 to 850 °C. Combining our data with those from the literature on pure H2O and pure CO2 solubility, mixed H2O-CO2 solubility, and H2O solubility for H2O-H2 gas mixture (240 measurements for H2O and 59 points for CO2), we develop a single empirical model for solubility of pure H2 O vapor, that of pure CO2 vapor, as well as that of mixed H2O-CO2 vapor in metaluminous high-silica rhyolitic melt. This model applies to the conditions of 700-1200 °C and 0-500 MPa, and consists of two equations: H2Ot = (354.94 Pw0.5 + 9.623 Pw - 1.5223 Pw1.5)/T + 0.0012439Pw1.5 + PCO2 (- 1.084 × 10-4 Pw0.5 - 1.362 × 10-5 Pw) CO2 = PCO2 (5668 - 55.99 Pw)/T + PCO2 (0.4133 Pw0.5 + 2.041 × 10-3 Pw1.5) where H2Ot is total dissolved H2O content in wt.%, CO2 content is in ppm by mass, and T is temperature in Kelvin, Pw=XwfP and PCO2 = XCO2fP (in MPa), where Xwf and XCO2f are the mole fraction of water and CO2 in the fluid. The 2σ relative uncertainty is 15% for the H2O equation, and 20% for the CO2 equation. The above H2O equation also applies to H2O solubility in mixed H2O-H2 fluid for XH2f<0.54 at 100 MPa, for XH2 f<0.49 at 200 MPa, and XH2f≥0.44 at 300 MPa (that is, the presence of H2 in the fluid insignificantly affects H2O solubility). This empirical model is recommended for the modeling of explosive volcanic eruptions and magma chamber dynamics.

AB - Experiments of H2O solubility in synthetic haplogranitic and natural rhyolitic melts were conducted at 700-1200 °C and 0.098-25 MPa. Attainment of equilibrium was proven by reversals using starting materials with different initial H2O contents. Dissolved H2O contents were determined using Fourier transform infrared spectroscopy. When temperature decreases from 1000 to 700 °C, the solubility of H2O increases from 0.100 to 0.124 wt.% at 0.098 MPa, from 0.99 to 1.36 wt.% at 11 MPa, and from 1.46 to 2.17 wt.% at 25 MPa. At 6 MPa, the solubility of H2O increases from 0.63 to 0.77 wt.% from 1200 to 850 °C. Combining our data with those from the literature on pure H2O and pure CO2 solubility, mixed H2O-CO2 solubility, and H2O solubility for H2O-H2 gas mixture (240 measurements for H2O and 59 points for CO2), we develop a single empirical model for solubility of pure H2 O vapor, that of pure CO2 vapor, as well as that of mixed H2O-CO2 vapor in metaluminous high-silica rhyolitic melt. This model applies to the conditions of 700-1200 °C and 0-500 MPa, and consists of two equations: H2Ot = (354.94 Pw0.5 + 9.623 Pw - 1.5223 Pw1.5)/T + 0.0012439Pw1.5 + PCO2 (- 1.084 × 10-4 Pw0.5 - 1.362 × 10-5 Pw) CO2 = PCO2 (5668 - 55.99 Pw)/T + PCO2 (0.4133 Pw0.5 + 2.041 × 10-3 Pw1.5) where H2Ot is total dissolved H2O content in wt.%, CO2 content is in ppm by mass, and T is temperature in Kelvin, Pw=XwfP and PCO2 = XCO2fP (in MPa), where Xwf and XCO2f are the mole fraction of water and CO2 in the fluid. The 2σ relative uncertainty is 15% for the H2O equation, and 20% for the CO2 equation. The above H2O equation also applies to H2O solubility in mixed H2O-H2 fluid for XH2f<0.54 at 100 MPa, for XH2 f<0.49 at 200 MPa, and XH2f≥0.44 at 300 MPa (that is, the presence of H2 in the fluid insignificantly affects H2O solubility). This empirical model is recommended for the modeling of explosive volcanic eruptions and magma chamber dynamics.

KW - Exsolution enthalpy

KW - HO solubility

KW - HO-CO solubility

KW - Metaluminous rhyolite

KW - Solubility model

KW - Volcanic eruption

UR - http://www.scopus.com/inward/record.url?scp=19944364251&partnerID=8YFLogxK

U2 - 10.1016/j.jvolgeores.2004.09.019

DO - 10.1016/j.jvolgeores.2004.09.019

M3 - Article

AN - SCOPUS:19944364251

VL - 143

SP - 219

EP - 235

JO - Journal of Volcanology and Geothermal Research

JF - Journal of Volcanology and Geothermal Research

SN - 0377-0273

IS - 1-3

ER -