Solubility of CO2 in a Ca-rich leucitite: effects of pressure,temperature, and oxygen fugacity

The solubility of carbon dioxide in a Ca-rich leucitite has been investigated as a function of pressure (0.1–2.0 GPa), temperature (1200–1600°C), and oxygen fugacity. The experiments were done in a rapid-quench internally-heated pressure vessel (0.1 GPa) and a piston cylinder (0.5–2.0 GPa). The leucitite glass, previously equilibrated at NNO, and silver oxalate were loaded in Fe-doped Pt capsules (oxidized conditions) and graphitelined Pt capsules (reduced conditions). Secondary Ion Mass Spectrometry and bulk carbon analyses were used to determine the amount of dissolved carbon. Speciation of carbon was characterized by Fourier transform microinfrared spectroscopy. At oxidized conditions, only CO₃ ^2- is observed as a dissolved species. The solubility is high with CO₂ contents in the melt attaining 6.2 wt% at 2.0 GPa and 1350°C. The solubility increases with pressure and shows a significant negative temperature dependence. An excellent correlation is obtained when the data are fit to a model, based on the simplified solubility reaction CO₂ ^(vapor)+O^2-(melt)CO₃ ^2-(melt), which describes the solubility of CO₂ as a function of pressure, temperature and fCO₂. At reduced conditions, the amount of carbon dissolved is significantly lower, and CO₃ ^2- is still the only species present in the melt. If the solubility model established at oxidized conditions is applied, the carbon dissolved appears to be essentially a function of fCO₂ alone although divergence increases in a consistent manner with pressure and temperature. This could suggest a low but significant solubility of CO with a positive temperature dependence or a departure of the calculated fluid compositions determined by the equation of state from the actual ones. The strong preferential solubility of carbon in its oxidized C⁴⁺ form, even at reduced conditions, implies that ascending melts with high CO₂ solubility can experience significant oxidation through degassing. This could reconcile the oxidized nature of some Ca-rich alkaline magmas with more reduced mantle source regions.     

CO2 solubility in dacitic melts equilibrated with H2O-CO2 fluids: Implications for modeling the solubility of CO2 in silicic melts

The solubility of CO₂ in dacitic melts equilibrated with H₂O-CO₂ fluids was experimentally investigated at 1250°C and 100 to 500 MPa. CO₂ is dissolved in dacitic glasses as molecular CO₂ and carbonate. The quantification of total CO₂ in the glasses by mid-infrared (MIR) spectroscopy is difficult because the weak carbonate bands at 1430 and 1530 cm⁻¹ can not be reliably separated from background features in the spectra. Furthermore, the ratio of CO_2,mol/carbonate in the quenched glasses strongly decreases with increasing water content. Due to the difficulties in quantifying CO₂ species concentrations from the MIR spectra we have measured total CO₂ contents of dacitic glasses by secondary ion mass spectrometry (SIMS).At all pressures, the dependence of CO₂ solubility in dacitic melts on x^fluid_CO2,total shows a strong positive deviation from linearity with almost constant CO₂ solubility at x_CO2fluid > 0.8 (maximum CO₂ solubility of 795 ± 41, 1376 ± 73 and 2949 ± 166 ppm at 100, 200 and 500 MPa, respectively), indicating that dissolved water strongly enhances the solubility of CO₂. A similar nonlinear variation of CO₂ solubility with x_CO2fluid has been observed for rhyolitic melts in which carbon dioxide is incorporated exclusively as molecular CO₂ (Tamic et al., 2001). We infer that water species in the melt do not only stabilize carbonate groups as has been suggested earlier but also CO₂ molecules.A thermodynamic model describing the dependence of the CO₂ solubility in hydrous rhyolitic and dacitic melts on T, P, f_CO2 and the mol fraction of water in the melt (x_water) has been developed. An exponential variation of the equilibrium constant K₁ with x_water is proposed to account for the nonlinear dependence of x_CO2,totalmelt on x_CO2fluid. The model reproduces the CO₂ solubility data for dacitic melts within ±14% relative and the data for rhyolitic melts within 10% relative in the pressure range 100-500 MPa (except for six outliers at low x_CO2fluid). Data obtained for rhyolitic melts at 75 MPa and 850°C show a stronger deviation from the model, suggesting a change in the solubility behavior of CO₂ at low pressures (a Henrian behavior of the CO₂ solubility is observed at low pressure and low H₂O concentrations in the melt). We recommend to use our model only in the pressure range 100-500 MPa and in the x_CO2fluid range 0.1-0.95. The thermodynamic modeling indicates that the partial molar volume of total CO₂ is much lower in rhyolitic melts (31.7 cm³/mol) than in dacitic melts (46.6 cm³/mol). The dissolution enthalpy for CO₂ in hydrous rhyolitic melts was found to be negligible. This result suggests that temperature is of minor importance for CO₂ solubility in silicic melts.     

An evaluation of the carbon-14 beta track technique: Implications for solubilities and partition coefficients determined by beta track mapping

The accuracy and precision of absolute C concentrations measured by beta track autoradiography (BTA) have been evaluated by synthesizing ¹⁴CO₂-bearing glasses (0.5–2.0 wt.% CO₂) of albite and diopside composition and comparing the amount of CO₂ introduced into the glass to the amount measured by BTA. Measured CO₂ contents in diopside glass are less than loaded CO₂ contents, but agree within experimental error; measured contents in albite glass are greater than loaded contents by 10–50%, which exceeds experimental error in many cases. Carbon isotope fractionation between calcite and CO₂ vapor is not significant above 1000°C at 0.1 or 3.0 GPa and cannot explain the discrepancies. The discrepancies represent an error in BTA which has not been previously recognized and is apparently related to beta absorption.Based on this study, CO₂ concentrations in albite measured by BTA should be corrected downward by 40%, which would bring such solubilities into better agreement with those determined by infrared spectroscopy. CO₂ solubilities in diopside melts measured by BTA are in agreement with phase equilibria and are accurate to within 20%. Partition coefficients measured by BTA may be subject to analytical errors inherent to autoradiography. Lack of awareness of these potential errors casts doubt on the accuracy of these measurements in many cases.     

Solubility of water and carbon dioxide in an icelandite at 1400 °C and 10 kilobars

The concentrations of water and carbon dissolved in an icelandite glass quenched from 1400 °C and 10 kbar were measured using Fourier transform infra-red spectroscopy and elemental analyses of carbon and hydrogen. Only carbon dioxide and water were observed in the fluid phase as analysed after quenching with a qudrupole mass analyser. The mole fraction of carbon dioxide in the fluid phase ranged from 0.36 to 0.95. Carbon is dissolved as carbonate except at the highest CO₂ fluid fugacity, where a small amount of molecular CO₂ is observed. Dissolved carbon in the glasses, calculated as CO₂, remained constant at approximately 1 wt %, in spite of the different CO₂ fluid fugacities. Water was dissolved as molecular water and as hydroxyl groups, the hydroxyl concentration in the quenched glasses remaining almost constant over the whole interval, whereas the molecular water dissolves in accordance with Henry's law. Molecular water peaks at 5200␣cm⁻¹ and 1630 cm⁻¹, the hydroxyl peak at 4500␣cm⁻¹, and the carbonate peaks at 1400 cm⁻¹–1550 cm⁻¹ have been calibrated using elemental analyses of C and H in the quenched glasses. As molecular water decreases in the melt the higher wavenumber carbonate peak is observed to move towards the molecular water peak at 1630 cm⁻¹ causing a split of the carbonate peaks, ranging from 45 cm⁻¹ to 100 cm⁻¹. Received: 15 November 1995 / Accepted: 21 September 1996     

CO2 solubility and solubility mechanisms in silicate melts at high pressures

CO₂ solubility has a slight negative temperature dependence in olivine melilitite at 30 kb with 9% CO₂ dissolved at 1,450 °C, 8.5% at 1,550 °C and 8.3% at 1,650° C. CO₂ is dissolved as the carbonate molecule (CO ₃ ^2– ) only. Feldspar melts (albite-anorthite) dissolve much less CO₂ at 30 kb (around 2%) with a slight increase with increasing anorthite content. A CO₂ absorption peak in infrared spectra of albite-rich glasses diappears in favour of the CO ₃ ^2– peak with increasing anorthite content. It is inferred that CO₂ was present as CO ₃ ^2– in albite-rich melts also, but reverts to CO₂ during quenching because of bonding differences related to Ca²⁺ and Na⁺ in the melts.     

The speciation of carbon dioxide in sodium aluminosilicate glasses

Infrared spectroscopy has been used to study the speciation of CO₂ in glasses near the NaAlO₂-SiO₂ join quenched from melts held at high temperatures and pressures. Absorption bands resulting from the antisymmetric stretches of both molecular CO₂ (2,352 cm^–1) and CO ₃ ^2– (1,610 cm^–1 and 1,375 cm^–1) are observed in these glasses. The latter are attributed to distorted Na-carbonate ionic-complexes. Molar absorptivities of 945 liters/mole-cm for the molecular CO₂ band, 200 liters/mole-cm for the 1,610 cm^–1 band, and 235 liters/mole-cm for the 1,375 cm^–1 band have been determined. These molar absorptivities allow the quantitative determination of species concentrations in the glasses with a precision on the order of several percent of the amount present. The accuracy of the method is estimated to be ±15–20% at present.The ratio of molecular CO₂ to CO ₃ ^2– in sodium aluminosilicate glasses varies little for each silicate composition over the range of total dissolved CO₂ content (0–2%), pressure (15–33 kbar) and temperature (1,400–1,560° C) that we have studied. This ratio is, however, a strong function of silicate composition, increasing both with decreasing Na₂O content along the NaAlO₂-SiO₂ join and with decreasing Na₂O content in peraluminous compositions off the join.Infrared spectroscopic measurements of species concentrations in glasses provide insights into the molecular level processes accompanying CO₂ solution in melts and can be used to test and constrain thermodynamic models of CO₂-bearing melts. CO₂ speciation in silicate melts can be modelled by equilibria between molecular CO₂, CO ₃ ^2– , and oxygen species in the melts. Consideration of the thermodynamics of such equilibria can account for the observed linear relationship between molecular CO₂ and carbonate concentrations in glasses, the proposed linear relationship between total dissolved CO₂ content and the activity of CO₂ in melts, and observed variations in CO₂ solubility in melts.     

CO2, 13C/12C and H2O variability in natural basaltic glasses: a study comparing stepped heating and ftir spectroscopic techniques

A comparison of two independent techniques was used to assess the homogeneity of CO₂ and H₂O concentrations in a number of natural basaltic glasses. Variations in carbon concentration and isotopic ratio were determined by comparison of stepped heating data obtained in two different laboratories. Dissolved volatile concentrations were also obtained by stepped heating and Fourier Transform Infrared (FTIR) spectroscopy. Replicate stepped heating analyses of a mid-ocean ridge basaltic glass show that the concentration and ¹³C/¹²C of bulk magmatic and dissolved CO₂ vary by less than ±10% and ±0.5‰, respectively. A similar degree of correlation is observed for replicate stepped heating analyses of Mariana Trough glasses conducted in two different laboratories. Dissolved CO₂ concentrations determined by stepped heating also correlate well with concentrations measured by FTIR spectroscopy. The correspondence of results obtained in these experiments provide an upper limit to the degree of natural variation in concentrations and isotopic ratios of these volatiles in basaltic glasses and suggest that intrinsic, magmatic carbon has a relatively homogeneous distribution in these glasses. Water concentrations determined through extraction by heating and FTIR also show excellent agreement.     

Influence of composition and thermal history of volcanic glasses on water content as determined by micro-Raman spectrometry

Development of Raman spectrometry for quantification of water content in natural glasses requires the assessment of the dependence of the technique on glass composition and thermal history. In the low frequency domain, Raman spectra topology varies due to glass depolymerization and substitution in the framework of (Si⁴⁺)^IV by alkali-balanced (Al³⁺)^IV and (Fe³⁺)^IV in calcalkaline (rhyolite to basaltic andesite) and alkaline (trachyte, phonolite to alkali basalt) glasses. These processes result in strong dependence of previous analytical procedure (internal calibration) on glass composition. Here, we show that an analytical procedure based on calibration to an external standard is only faintly composition-dependent for Si-rich alkaline glasses (trachytes–phonolites). For a given glass composition, thermal history also plays a fundamental role in the choice of Raman procedure for water analysis. Repeated cycles of thermal annealing induce microcrystallization of hydrous trachyte glasses and modify cation distribution in the glass structure. Application of these concepts to analysis of banded obsidians suggests that small-scale heterogeneities in glasses are not simply related to magma degassing, but could depend on thermal history and consequent relaxation paths in the melt.     

Modeling of the solubility of a one-component H2O or CO2 fluid in silicate liquids

The modeling of the solubility of water and carbon dioxide in silicate liquids (flash problem) is performed by assuming mechanical, thermal, and chemical equilibrium between the liquid magma and the gas phase. The liquid phase is treated as a mixture of ten silicate components + H₂O or CO₂, and the gas phase as a pure H₂O or CO₂. A general model for the solubility of a volatile component in a liquid is adopted. This requires the definition of a mixing equation for the excess Gibbs free energy of the liquid phase and an appropriate reference state for the dissolved volatile. To constrain the model parameters and identify the most appropriate form of the solubility equations for each dissolved volatile, a large number of experimental solubility determinations (640 for H₂O and 263 for CO₂) have been used. These determinations cover a large region of the P-T-composition space of interest. The resultant water and carbon dioxide solubility models differ in that the water model is regular and isometric, and the carbon dioxide model is regular and non-isometric. This difference is consistent with the different speciation modalities of the two volatiles in the silicate liquids, producing a composition-independent partial molar volume of dissolved water and a composition-dependent partial molar volume of dissolved carbon dioxide. The H₂O solubility model may be applied to natural magmas of virtually any composition in the P-T range 0.1 MPa–1 GPa and > 1000 K, whereas the CO₂ solubility model may be applied to several GPa pressures. The general consistency of the water solubility data and their relatively large number as compared to the calibrated model parameters (11) contrast with the large inconsistencies of the carbon dioxide solubility determinations and their low number with respect to the CO₂ model parameters (22). As a result, most of the solubility data in the database are reproduced within 10% of approximation in the case of water, and 30% in the case of carbon dioxide. When compared with the experimental data, the H₂O and CO₂ solubility models correctly predict many features of the saturation surface in the P-T-composition space, including the change from retrograde to prograde H₂O solubility in albitic liquids with increasing pressure, the so-called alkali effect, the increasing CO₂ solubility with increasing degree of silica undersaturation, the Henrian behavior of CO₂ in most silicate liquids up to about 30–50 MPa, and the proportionality between the fugacity in the gas phase, or the saturation activity in the liquid phase, and the square of the mole fraction of the dissolved volatile found in some unrelated silicate liquid compositions. Received: 21 August 1995 / Accepted: 8 July 1996     

Influence of glass polymerisation and oxidation on micro-Raman water analysis in alumino-silicate glasses

The development of an accurate analytical procedure for determination of dissolved water in complex alumino-silicate glasses via micro-Raman analysis requires the assessment of the spectra topology dependence on glass composition. We report here a detailed study of the respective influence of bulk composition, iron oxidation state and total water content on the absolute and relative intensities of the main Raman bands related to glass network vibrations (LF: ∼490 cm⁻¹; HF: ∼960 cm⁻¹) and total water stretching (H₂O_T: ∼3550 cm⁻¹) in natural glasses. The evolution of spectra topology was examined in (i) 33 anhydrous glasses produced by the re-melting of natural rock samples, which span a very large range of polymerisation degree (NBO/T from 0.00 to 1.16), (ii) 2 sets of synthetic anhydrous basaltic glasses with variable iron oxidation state (Fe³⁺/Fe_T from 0.05 to 0.87), and (iii) 6 sets of natural hydrous glasses (C_H2OT from 0.4 to 7.0 wt%) with NBO/T varying from 0.01 to 0.76.In the explored domain of water concentration, external calibration procedure based on the H₂O_T band height is matrix-independent but its accuracy relies on precise control of the focusing depth and beam energy on the sample. Matrix-dependence strongly affects the internal calibrations based on H₂O_T height scaled to that of LF or HF bands but its effect decreases from acid (low NBO/T, SM) to basic (high NBO/T, SM) glasses. Structural parameters such as NBO/T (non-bridging oxygen per tetrahedron) and SM (sum of structural modifiers) describe the matrix-dependence better than simple compositional parameters (e.g. SiO₂, Na₂O + K₂O). Iron oxidation state has only a minor influence on band topology in basalts and is thus not expected to significantly affect the Raman determinations of water in mafic (e.g. low SiO₂, iron-rich) glasses. Modelling the evolution of the relative band height with polymerisation degree allows us to propose a general equation to predict the dissolved water content in natural glasses:

     

Carbon isotopic fractionation between CO2 vapour,silicate and carbonate melts: an experimental study to 30 kbar

The carbon isotopic fractionation between CO₂ vapour and sodamelilite (NaCaAlSi₂O₇) melt over a range of pressures and temperatures has been investigated using solid-media piston-cylinder high pressure apparatus. Ag₂C₂O₄ was the source of CO₂ and experimental oxygen fugacity was buffered at hematite-magnetite by the double capsule technique. The abundance and isotopic composition of carbon dissolved in sodamelilite (SM) glass were determined by stepped heating and the ¹³C of coexisting vapour was determined directly by capsule piercing. CO₂ solubility in SM displays a complex behavior with temperature. At pressures up to 10 kbars CO₂ dissolves in SM to form carbonate ion complexes and the solubility data suggest slight negative temperature dependence. Above 20 kbars CO₂ reacts with SM to form immiscible Na-rich silicate and Ca-rich carbonate melts and CO₂ solubility in Na-enriched silicate melt rises with increasing temperature above the liquidus. Measured values for carbon isotopic fractionation between CO₂ vapour and carbonate ions dissoived in sodamelilite melt at 1200°–1400° C and 5–30 kbars average 2.4±0.2, favouring¹³C enrichment in CO₂ vapour. The results are maxima and are independent of pressure and temperature. Similar values of 2 are obtained for the carbon isotopic fractionation between CO₂ vapour and carbonate melts at 1300°–1400° C and 20–30 kbars.     

On leakage and seepage of CO2 from geologic storage sites into surface water

Geologic carbon sequestration is the capture of anthropogenic carbon dioxide (CO₂) and its storage in deep geologic formations. The processes of CO₂ seepage into surface water after migration through water-saturated sediments are reviewed. Natural CO₂ and CH₄ fluxes are pervasive in surface-water environments and are good analogues to potential leakage and seepage of CO₂. Buoyancy-driven bubble rise in surface water reaches a maximum velocity of approximately 30 cm s⁻¹. CO₂ rise in saturated porous media tends to occur as channel flow rather than bubble flow. A comparison of ebullition versus dispersive gas transport for CO₂ and CH₄ shows that bubble flow will dominate over dispersion in surface water. Gaseous CO₂ solubility in variable-salinity waters decreases as pressure decreases leading to greater likelihood of ebullition and bubble flow in surface water as CO₂ migrates upward.     

CO2 solubility and speciation in intermediate (andesitic) melts: the role of H2O and composition

We determined total CO₂ solubilities in andesite melts with a range of compositions. Melts were equilibrated with excess C-O(-H) fluid at 1 GPa and 1300°C then quenched to glasses. Samples were analyzed using an electron microprobe for major elements, ion microprobe for C-O-H volatiles, and Fourier transform infrared spectroscopy for molecular H₂O, OH⁻, molecular CO₂, and CO₃²⁻. CO₂ solubility was determined in hydrous andesite glasses and we found that H₂O content has a strong influence on C-O speciation and total CO₂ solubility. In anhydrous andesite melts with ∼60 wt.% SiO₂, total CO₂ solubility is ∼0.3 wt.% at 1300°C and 1 GPa and total CO₂ solubility increases by about 0.06 wt.% per wt.% of total H₂O. As total H₂O increases from ∼0 to ∼3.4 wt.%, molecular CO₂ decreases (from 0.07 ± 0.01 wt.% to ∼0.01 wt.%) and CO₃²⁻ increases (from 0.24 ± 0.04 wt.% to 0.57 ± 0.09 wt.%). Molecular CO₂ increases as the calculated mole fraction of CO₂ in the fluid increases, showing Henrian behavior. In contrast, CO₃²⁻ decreases as the calculated mole fraction of CO₂ in the fluid increases, indicating that CO₃²⁻ solubility is strongly dependent on the availability of reactive oxygens in the melt. These findings have implications for CO₂ degassing. If substantial H₂O is present, total CO₂ solubility is higher and CO₂ will degas at relatively shallow levels compared to a drier melt. Total CO₂ solubility was also examined in andesitic glasses with additional Ca, K, or Mg and low H₂O contents (<1 wt.%). We found that total CO₂ solubility is negatively correlated with (Si + Al) cation mole fraction and positively correlated with cations with large Gibbs free energy of decarbonation or high charge-to-radius ratios (e.g., Ca). Combining our andesite data with data from the literature, we find that molecular CO₂ is more abundant in highly polymerized melts with high ionic porosities (>∼48.3%), and low nonbridging oxygen/tetrahedral oxygen (<∼0.3). Carbonate dominates most silicate melts and is most abundant in depolymerized melts with low ionic porosities, high nonbridging oxygen/tetrahedral oxygen (>∼0.3), and abundant cations with large Gibbs free energy of decarbonation or high charge-to-radius ratio. In natural silicate melt, the oxygens in the carbonate are likely associated with tetrahedral and network-modifying cations (including Ca, H, or H-bonds) or a combinations of those cations.     

Fluid infiltration and regional metamorphism of the Waits River Formation, north-east Vermont, USA

Abstract The Siluro-Devonian Waits River Formation of north-east Vermont was deformed, intruded by plutons and regionally metamorphosed during the Devonian Acadian Orogeny. Five metamorphic zones were mapped based on the mineralogy of carbonate rocks. From low to high grade, these are: (1) ankerite-albite, (2) ankerite-oligoclase, (3) biotite, (4) amphibole and (5) diopside zones. Pressure was near 4.5kbar and temperature varied from c. 450° C in the ankerite-albite zone to c. 525° C in the diopside zone. Fluid composition for all metamorphic zones was estimated from mineral equilibria. Average calculated χ_co2[= CO₂/(CO₂+ H₂O)] of fluid in equilibrium with the marls increases with increasing grade from 0.05 in the ankerite-oligoclase zone, to 0.25 in the biotite zone and to 0.44 in the amphibole zone. In the diopside zone, χ_CO2 decreases to 0.06. Model prograde metamorphic reactions were derived from measured modes, mineral chemistry, and whole-rock chemistry. Prograde reactions involved decarbonation with an evolved volatile mixture of χ_CO2 > 0.50. The χ_CO2 of fluid in equilibrium with rocks from all zones, however, was generally <0.40. This difference attests to the infiltration of a reactive H₂O-rich fluid during metamorphism. Metamorphosed carbonate rocks from the formation suggests that both heat flow and pervasive infiltration of a reactive H₂O-rich fluid drove mineral reactions during metamorphism. Average time-integrated volume fluxes (cm³ fluid/cm² rock), calculated from the standard equation for coupled fluid flow and reaction in porous media, are (1) ankerite-oligoclase zone: c. 1 × 10⁴; (2) biotite zone: c. 3 × 10⁴; (3) amphibole zone: c. 10 × 10⁴; and diopside zone: c. 60 × 10⁴. The increase in calculated flux with increasing grade is at least in part the result of internal production of volatiles from prograde reactions in pelitic schists and metacarbonate rocks within the Waits River Formation. The mapped pattern of time-integrated fluxes indicates that the Strafford-Willoughby Arch and the numerous igneous intrusions in the field area focused fluid flow during metamorphism. Many rock specimens in the diopside zone experienced extreme alkali depletion and also record low χ_CO2. Metamorphic fluids in equilibrium with diopside zone rocks may therefore represent a mixture of acid, H₂O-rich fluids given off by the crystallizing magmas, and CO₂-H₂O fluids produced by devolatilization reactions in the host marls. Higher fluxes and different fluid compositions recorded near the plutons suggest that pluton-driven hydrothermal cells were local highs in the larger regional metamorphic hydrothermal system.     

CO_2 Density-Raman Shift Relation Derived from Synthetic Inclusions in Fused Silica Capillaries and Its Application

Abstract: The densities of CO2 inclusions in minerals are commonly used to determine the crystallizing conditions of the host minerals. However, conventional microthermometry is difficult to apply for inclusions of small size (< 5–10 μm) or low density. Raman analysis is an alternative method for determining CO2 density, provided that the CO2 density–Raman shift relation is known. This study aims to establish this CO2 density–Raman shift relation by using CO2 inclusions synthesized in fused silica capillaries. By using this newly-developed synthetic technique, we formed pure CO2 inclusions, and their densities were determined by microthermometry. The Raman analysis showed that the relation between CO2 density (D in g/cm3) and the separations (Δ in cm?1) between the two main bands (i.e. Fermi diad bands) in CO2 Raman spectra can be represented by a cubic equation: D (g/cm3)=0.74203(?0.019Δ3+5.90332Δ2?610.79472Δ+21050.30165)?3.54278 (r2=0.99920). Our calculated D value for a given Δ is between those obtained from two previously-reported equations, which were derived from different experimental methods. An example was given in this study to demonstrate that the densities of natural CO2 inclusions that could not be derived from microthermometry could be determined by using our method.     

CO2 in haplo-phonolite melt: solubility, speciation and carbonate complexation

CO₂ solubility was measured in a synthetic iron-free phonolite (haplo-phonolite) by equilibrating melt with excess CO₂ fluid in a piston cylinder apparatus for a range of pressures (1.0- 2.5 GPa) and temperatures (1300 to 1550°C). The quenched glasses were then analysed using a bulk carbon analytical method (LECO). The measured solubilities are between 0.65 and 2.77 wt.% for the range of conditions studied and show a negative correlation with temperature as reported for most other silicate melt compositions.A range of carbonate species are present within the glass, as well as minor amounts of molecular CO₂. FTIR and NMR analyses suggest that carbonate is present as both ‘network’ and ‘depolymerised’ units as shown for relatively highly polymerised compositions in the model of Brooker et al. (2001b). The bulk CO₂ analyses were used to calibrate the IR extinction coefficient for the carbonate groups. However, the results show that the values obtained for the glasses vary with the melt equilibration conditions, presumably because the ratio of the different carbonate species changes as a complex function of run pressure, temperature and quench rate. Thus the use of IR may not be a reliable method for the quantification of dissolved CO₂ concentrations in natural glasses of ‘intermediate’ composition.     

C�CO�CH�CS fluids and granitic magma: how S partitions and modifies CO2 concentrations of fluid-saturated felsic melt at 200?MPa

Hydrothermal volatile-solubility and partitioning experiments were conducted with fluid-saturated haplogranitic melt, H₂O, CO₂, and S in an internally heated pressure vessel at 900°C and 200?MPa; three additional experiments were conducted with iron-bearing melt. The run-product glasses were analyzed by electron microprobe, FTIR, and SIMS; and they contain ??0.12 wt% S, ??0.097 wt% CO₂, and ??6.4 wt% H₂O. Apparent values of log f _O2 for the experiments at run conditions were computed from the [(S⁶⁺)/(S⁶⁺+S^2?)] ratio of the glasses, and they range from NNO ?0.4 to NNO?+?1.4. The C?CO?CH?CS fluid compositions at run conditions were computed by mass balance, and they contained 22?C99?mol% H₂O, 0?C78?mol% CO₂, 0?C12?mol% S, and <3 wt% alkalis. Eight S-free experiments were conducted to determine the H₂O and CO₂ concentrations of melt and fluid compositions and to compare them with prior experimental results for C?CO?CH fluid-saturated rhyolite melt, and the agreement is excellent. Sulfur partitions very strongly in favor of fluid in all experiments, and the presence of S modifies the fluid compositions, and hence, the CO₂ solubilities in coexisting felsic melt. The square of the mole fraction of H₂O in melt increases in a linear fashion, from 0.05 to 0.25, with the H₂O concentration of the fluid. The mole fraction of CO₂ in melt increases linearly, from 0.0003 to 0.0045, with the CO₂ concentration of C?CO?CH?CS fluids. Interestingly, the CO₂ concentration in melts, involving relatively reduced runs (log f _O2????NNO?+?0.3) that contain 2.5?C7?mol% S in the fluid, decreases significantly with increasing S in the system. This response to the changing fluid composition causes the H₂O and CO₂ solubility curve for C?CO?CH?CS fluid-saturated haplogranitic melts at 200?MPa to shift to values near that modeled for C?CO?CH fluid-saturated, S-free rhyolite melt at 150?MPa. The concentration of S in haplogranitic melt increases in a linear fashion with increasing S in C?CO?CH?CS fluids, but these data show significant dispersion that likely reflects the strong influence of f _O2 on S speciation in melt and fluid. Importantly, the partitioning of S between fluid and melt does not vary with the (H₂O/H₂O?+?CO₂) ratio of the fluid. The fluid-melt partition coefficients for H₂O, CO₂, and S and the atomic (C/S) ratios of the run-product fluids are virtually identical to thermodynamic constraints on volatile partitioning and the H, S, and C contents of pre-eruptive magmatic fluids and volcanic gases for subduction-related magmatic systems thus confirming our experiments are relevant to natural eruptive systems.     

Experimental determination of argon solubility in silicate melts: An assessment of the effects of liquid composition and temperature

The argon solubility of 38 liquids in the system Na₂O-CaO-MgO-Al₂O₃-SiO₂ (NCMAS) has been determined at 1873 K and 1 bar, the argon concentration of presaturated glasses being measured using a static mass spectrometer. For compositions in the subsystem diopside (CaMgSi₂O₆), nepheline (NaAlSiO₄), albite (NaAlSi₃O₈), anorthite (CaAl₂Si₂O₈), argon solubility is generally a linear function of the relative proportion of each end member, solubility being lowest in diopside melt (1.53 10⁻⁵ cm³ STP · g⁻¹ · bar⁻¹) and highest in albite melt (2.88 10⁻⁴ cm³ STP · g⁻¹ · bar⁻¹). For the tectosilicate joins studied (SiO₂-Na₂Al₂O₄, SiO₂-CaAl₂O₄, SiO₂-MgAl₂O₄) solubility decreases with decreasing silica content in all cases, being highest for Na-bearing liquids and lowest for Mg-bearing liquids at constant molar silica content. Where comparison is possible our results are in good agreement with data from the literature. When our data are considered in isolation we find that argon solubility shows an excellent correlation with calculated ionic porosity. The covariation of argon solubility and liquid density is also reasonable, that with molar volume less convincing and that with polymerization state (as defined by the ratio of the number of nonbridging oxygens and tetrahedral network forming cations; NBO/T) nonexistent. However, when our data are combined with those from the literature no well constrained correlation between argon solubility and ionic porosity is apparent. Based upon this observation and consideration of the temperature dependence of noble gas solubility it is concluded that ionic porosity is not a universally applicable parameter which may be used to predict noble gas solubility as a function of composition, temperature and pressure. Two new models for calculating argon solubility are proposed, both employing the notion of partial molar argon solubilities. The first uses oxide components, for which partial molar argon solubility is directly proportional to partial molar ionic porosity calculated at 1873 K, irrespective of the temperature of experimental equilibration. The second model, which offers the best fit to the available data, employs tetrahedral units rather than oxides as the proposed melt components. This latter model successfully accounts for reported argon solubilities in simple Al-free systems, in simple Al-bearing systems and in natural liquids. This is interpreted to infer that argon is incorporated in large sites in the liquid structure (such as the space within rings of n-tetrahedra) although further work is required to understand the quantitative links between melt structure and noble gas solubility.     

Post-emplacement serpentinization and related hydrothermal metamorphism in a kimberlite from Venetia, South Africa

The common serpentine–diopside matrix assemblage in volcaniclastic kimberlite (VK) at the Venetia Mine, South Africa is ascribed to a secondary origin, because of post‐emplacement serpentinization and associated hydrothermal metamorphism. Volcaniclastic deposits with 20–30% porosity infill kimberlite pipes in the waning stages of kimberlite eruptions. Olivine macrocrysts are typically rimmed by talc and are pseudomorphed by lizardite, with minor magnetite. The fine matrix consists of mixtures of lizardite, chlorite, smectite, brucite, calcite, titanite and andradite, an assemblage which either pseudomorphed microcrysts or in‐filled voids. Locally we recognize microcryst pseudomorphs rich in sub‐microscopic mixtures of lizardite with smectite, and other microcryst pseudomorphs and void‐filling matrix rich in chlorite and lizardite. Interstitial lizardite and associated phyllosilicates (brucite, smectite and chlorite) crystallized progressively from meteoric or hydrothermally derived pore waters, and Si⁴⁺ and Mg²⁺ released into the fluid phase during serpentinization of olivine macrocrysts. Radial‐fibrous fringes of diopside microlites around crystals display void‐filling textures because of unrestricted growth into pore spaces. Secondary diopside is attributed to Si⁴⁺, Mg²⁺ and Ca²⁺ cations released into the fluid phase by interaction with olivine, calcite and plagioclase in siliceous xenoliths. The paucity of primary, fine‐grained groundmass phases resistant to alteration, for example, perovskite and spinel, precludes an origin for the intergrain matrix as altered interstitial ash, glass or a late‐stage kimberlite melt. Isovolumetric replacement of olivine results in a volume increase of 60% so that pore spaces in the original deposit can be easily filled up with serpentine. The source of Al³⁺ to form chlorite and smectite is attributed to alteration of plagioclase in xenoliths which comprise 20–30 vol.% of the deposit. Titanite, hydro‐andradite and second‐generation diopside precipitate as hydrothermal minerals from calcium‐bearing serpentinizing fluids in replacement reactions and as void‐filling minerals. Consideration of mineral equilibria in the CaO‐MgO‐SiO₂‐H₂O‐CO₂ system constrains the common matrix assemblage of lizardite and diopside in X_CO2)–T space. At 300 bar, the assemblage is stable only at temperatures below 370 °C and X_CO2 < 0.01. This upper limit on temperature is well below the plausible solidus of ultrabasic magmas. Furthermore, the requirement of trace CO₂ in the fluid phase implies a post‐emplacement external source rather than ‘autometamorphism’ from kimberlite‐derived fluids, because of high P_CO2 commonly inferred for kimberlite magmas.