Carbon dioxide in silicate melts: A molecular dynamics simulation study

We have performed a series of molecular dynamics simulations aimed at the evaluation of the solubility of CO₂ in silicate melts of natural composition (from felsic to ultramafic). In making in contact within the simulation cell a supercritical CO₂ phase with a silicate melt of a given composition, we have been able to evaluate (i) the solubility of CO₂ in the P-T range 1473-2273 K and 20-150 kbar, (ii) the density change experienced by the CO₂-bearing melt, (iii) the respective concentrations of CO₂ and species in the melt, (iv) the lifetime and the diffusivity of these species and (v) the structure of the melt around the carbonate groups. The main results are the following:(1) The solubility of CO₂ increases markedly with the pressure in the three investigated melts (a rhyolite, a mid-ocean ridge basalt and a kimberlite) from about ∼2 wt% CO₂ at 20 kbar to ∼25 wt% at 100 kbar and 2273 K. The solubility is found to be weakly dependent on the melt composition (as far as the present compositions are concerned) and it is only at very high pressure (above ∼100 kbar) that a clear hierarchy between solubilities occurs (rhyolite < MORB < kimberlite). Furthermore at a given pressure the calculated solubility is negatively correlated with the temperature.(2) In CO₂-saturated melts, the proportion of carbonate ions is positively correlated with the pressure at isothermal condition and is negatively correlated with the temperature at isobaric condition (and vice versa for molecular CO₂). Furthermore, at fixed (P, T) conditions the proportion of carbonate ions is higher in CO₂-undersaturated melts than in the CO₂-saturated melt. Although the proportion of molecular CO₂ decreases when the degree of depolymerization of the melt increases, it is still significant in CO₂-saturated basic and ultrabasic compositions at high temperatures. This finding is at variance with experimental data on CO₂-bearing glasses which show no evidence of molecular CO₂ as soon as the degree of depolymerization of the melt is high (e.g. basalt). These conflicting results can be reconciled with each other by noticing that a simple low temperature extrapolation of the simulation data predicts that the proportion of molecular CO₂ in basaltic melts might be negligible in the glass at room temperature.(3) The carbonate ions are found to be transient species in the liquid phase, with a lifetime increasing exponentially with the inverse of the temperature. Contrarily to a usual assumption, the diffusivity of carbonate ions into the liquid silicate is not vanishingly small with respect to that of CO₂ molecules: in MORB they differ from each other by a factor of ∼6 at 1473 K and only a factor of ∼2 at 2273 K. Although the bulk diffusivity of CO₂ is governed primarily by the diffusivity of CO₂ molecules, the carbonate ions contribute significantly to the diffusivity of CO₂ in depolymerized melts.(4) Concerning the structure of the CO₂-bearing silicate melt, the carbonate ions are found to be preferentially associated with NBO’s of the melt, with an affinity for NBOs which exceeds that for BOs by almost one order of magnitude. This result explains why the concentration in carbonate ions is positively correlated with the degree of depolymerization of the melt and diminishes drastically in fully polymerized melts where the number of NBO’s is close to zero. Furthermore, the network modifier cations are not randomly distributed in the close vicinity of carbonate groups but exhibit a preferential ordering which depends at once on the nature of the cation and on the melt composition. However at the high temperatures investigated here, there is no evidence of long lived complexes between carbonate groups and metal cations.     

The effect of bulk composition on the solidus of carbonated eclogite from partial melting experiments at 3 GPa

To explore the effect of bulk composition on the solidus of carbonated eclogite, we determined near-solidus phase relations at 3 GPa for four different nominally anhydrous, carbonated eclogites. Starting materials (SLEC1, SLEC2, SLEC3, and SLEC4) were prepared by adding variable proportions and compositions of carbonate to a natural eclogite xenolith (66039B) from Salt Lake crater, Hawaii. Near-solidus partial melts for all bulk compositions are Fe–Na calcio-dolomitic and coexist with garnet + clinopyroxene + ilmenite ± calcio-dolomitic solid solution. The solidus for SLEC1 (Ca#=100 × molar Ca/(Ca + Mg + Fe_T)=32, 1.63 wt% Na₂O, and 5 wt% CO₂) is bracketed between 1,050°C and 1,075°C (Dasgupta et al. in Earth Planet Sci Lett 227:73–85, 2004), whereas initial melting for SLEC3 (Ca# 41, 1.4 wt% Na₂O, and 4.4 wt% CO₂) is between 1,175°C and 1,200°C. The solidus for SLEC2 (Ca# 33, 1.75 wt% Na₂O, and 15 wt% CO₂) is estimated to be near 1,100°C and the solidus for SLEC3 (Ca# 37, 1.47 wt% Na₂O, and 2.2 wt% CO₂) is between 1,100°C and 1,125°C. Solidus temperatures increase with increasing Ca# of the bulk, owing to the strong influence of the calcite–magnesite binary solidus-minimum on the solidus of carbonate bearing eclogite. Bulk compositions that produce near-solidus crystalline carbonate closer in composition to the minimum along the CaCO₃-MgCO₃ join have lower solidus temperatures. Variations in total CO₂ have significant effect on the solidus if CO₂ is added as CaCO₃, but not if CO₂ is added as a complex mixture that maintains the cationic ratios of the bulk-rock. Thus, as partial melting experiments necessarily have more CO₂ than that likely to be found in natural carbonated eclogites, care must be taken to assure that the compositional shifts associated with excess CO₂ do not unduly influence melting behavior. Near-solidus dolomite and calcite solid solutions have higher Ca/(Ca + Mg) than bulk eclogite compositions, owing to Ca–Mg exchange equilibrium between carbonates and silicates. Carbonates in natural mantle eclogite, which have low bulk CO₂ concentration, will have Ca/Mg buffered by reactions with silicates. Consequently, experiments with high bulk CO₂ may not mimic natural carbonated eclogite phase equilibria unless care is taken to ensure that CO₂ enrichment does not result in inappropriate equilibrium carbonate compositions. Compositions of eclogite-derived carbonate melt span the range of natural carbonatites from oceanic and continental settings. Ca#s of carbonatitic partial melts of eclogite vary significantly and overlap those of partial melts of carbonated lherzolite, however, for a constant Ca-content, Mg# of carbonatites derived from eclogitic sources are likely to be lower than the Mg# of those generated from peridotite.     

Carbon dioxide and argon diffusion in silicate melts: Insights into the CO2 speciation in magmas

To investigate the influence of temperature and composition on the diffusivities of dissolved carbon dioxide and argon in silicate melts, diffusion experiments were performed at magmatic pressure and temperature conditions in (a) albite melts with excess Na₂O (0-8.6 wt%) and a constant Si/Al ratio of 3, and (b) albite₇₀quartz₃₀ to jadeite melts with decreasing SiO₂ content and a constant Na/Al ratio of 1. We obtained diffusion coefficients at 500 MPa and 1323-1673 K. In the fully polymerized system Ab₇₀Qz₃₀ - Jd, the change in composition only has a weak effect on bulk CO₂ diffusivity, but Ar diffusivity increases clearly with decreasing SiO₂ content. In the system Ab + Na₂O, bulk CO₂ and Ar diffusivity increase significantly with gradual depolymerisation. The relatively small change in composition on molar basis in the depolymerized system leads to a significantly larger change in diffusivities compared to the fully polymerized Ab₇₀Qz₃₀-Jd join. Within error, activation energies for bulk CO₂ and Ar diffusion in both systems are identical with decreasing silica content (Ab + Na₂O: 159 ± 25 kJ mol⁻¹ for bulk CO₂ and 130 ± 8 kJ mol⁻¹ for Ar; Ab₇₀Qz₃₀-Jd: 163 ± 16 kJ mol⁻¹ for bulk CO₂ and 148 ± 15 kJ mol⁻¹ for Ar) even though this results in depolymerisation in one system and not the other.Although there is a variation in CO₂ speciation with changing composition as observed in quenched glasses, it has previously established that this is not a true representation of the species present in the melt, with the ratio of molecular CO₂ to carbonate decreasing during quenching. Thus, diffusion coefficients for the individual CO₂ species cannot be directly derived by measuring molecular CO₂ and CO₃^2- concentration-distance profiles in the glasses. To obtain diffusivities of individual CO₂ species, we have made two assumptions that (1) inert Ar can be used as a proxy for molecular CO₂ diffusion characteristics as shown by our previous work and (2) the diffusivity of CO₃²⁻ can be calculated assuming it is identical to network forming components (Si⁴⁺ and Al³⁺). This is derived from viscosity data (Eyring eqn.) and suggests that CO₃²⁻ diffusion would be several orders of magnitude slower than molecular CO₂ diffusion.The systematics of measured bulk CO₂ diffusivity rates and comparison with the Ar proxy all suggest that the faster molecular CO₂ species is much more dominant in melts than measurements on resulting quenched glasses would suggest. This study has confirmed an observation of surprisingly consistent bulk CO₂ diffusivity across a range of natural compositions were Ar diffusivity significantly increases. This is consistent with an actual increase in molecular CO₂ mobility (similar to Ar) that is combined with an increase in the proportion of the slower carbonate in the melt.These results demonstrate that the CO₂ diffusion and speciation model provides an insight into the transport processes in the melt and is promising and an alternative tool to in situ speciation measurements at magmatic conditions, which at the moment are technically extremely difficult. We present the first high pressure high temperature in situ MIR spectra of a CO₂ bearing albitic glass/melt suggesting that molecular CO₂ is a stable species at high temperature, which is qualitatively consistent with the modelled CO₂ speciation data.     

Defining an absolute reference frame for ‘clumped’ isotope studies of CO2

We present a revised approach for standardizing and reporting analyses of multiply substituted isotopologues of CO₂ (i.e., ‘clumped’ isotopic species, especially the mass-47 isotopologues). Our approach standardizes such data to an absolute reference frame based on theoretical predictions of the abundances of multiply-substituted isotopologues in gaseous CO₂ at thermodynamic equilibrium. This reference frame is preferred over an inter-laboratory calibration of carbonates because it enables all laboratories measuring mass 47 CO₂ to use a common scale that is tied directly to theoretical predictions of clumping in CO₂, regardless of the laboratory’s primary research field (carbonate thermometry or CO₂ biogeochemistry); it explicitly accounts for mass spectrometric artifacts rather than convolving (and potentially confusing) them with chemical fractionations associated with sample preparation; and it is based on a thermodynamic equilibrium that can be experimentally established in any suitably equipped laboratory using commonly available materials.By analyzing CO₂ gases that have been subjected to established laboratory procedures known to promote isotopic equilibrium (i.e., heated gases and water-equilibrated CO₂), and by reference to thermodynamic predictions of equilibrium isotopic distributions, it is possible to construct an empirical transfer function that is applicable to data with unknown clumped isotope signatures. This transfer function empirically accounts for the fragmentation and recombination reactions that occur in electron impact ionization sources and other mass spectrometric artifacts. We describe the protocol necessary to construct such a reference frame, the method for converting gases with unknown clumped isotope compositions to this reference frame, and suggest a protocol for ensuring that all reported isotopic compositions (e.g., Δ₄₇ values; [Eiler and Schauble, 2004] and [Eiler, 2007]) can be compared among different laboratories and instruments, independent of laboratory-specific analytical or methodological differences. We then discuss the use of intra-laboratory secondary reference frames (e.g., based on carbonate standards) that can be more easily used to track the evolution of each laboratory’s empirical transfer function. Finally, we show inter-laboratory reproducibility on the order of ±0.010 (1σ) for four carbonate standards, and present revised paleotemperature scales that should be used to convert carbonate clumped isotope signatures to temperature when using the absolute reference frame described here. Even when using the reference frame, small discrepancies remain between two previously published synthetic carbonate calibrations. We discuss possible reasons for these discrepancies, and highlight the need for additional low temperature (<15 °C) synthetic carbonate experiments.     

Rare earth partitioning between immiscible carbonate and silicate liquids and CO2 vapor: Results and implications for the formation of light rare earth-enriched rocks

Melting relations at 5 and 20 kbar on the composition join sanidine-potassium carbonate are dominated by a two-liquid region that covers over 60% of the join at 1,300 ° C. At this temperature, the silicate melt contains approximately 19 wt% carbonate component at 5 kbar and 32 wt% carbonate component at 20 kbar. The conjugate carbonate melt contains less than 5 wt% silicate component, and it varies less as a function of temperature than does the silicate melt.Partition coefficients for Ce, Sm, and Tm between the immiscible carbonate and silicate melts at 1,200 ° and 1,300 ° C at 5 and 20 kbar are in favor of the carbonate melt by a factor of 2–3 for light REE and 5–8 for heavy REE. The effect of pressure on partitioning cannot be evaluated independently because of complementary changes in melt compositions.Minimum REE partition coefficients for CO₂ vapor/carbonate melt and CO₂ vapor/silicate melt can be calculated from the carbonate melt/silicate melt partition coefficients, the known proportions of melt, and maximum estimates of the proportion of CO₂ vapor. The vapor phase is enriched in light REE relative to both melts at 20 kbar and enriched in all REE, especially the light elements, at 5 kbar. The enrichment of REE in CO₂ vapor relative to both melts is 3–4 orders of magnitude in excess of that in water vapor (Mysen, 1979) at 5 kbar and is approximately the same as that in water vapor at 20 kbar.Mantle metasomatism by a CO₂-rich vapor enriched in light REE, occurring as a precursor to magma genesis, may explain the enhanced REE contents and light REE enrichment of carbonatites, alkali-rich silicate melts, and kimberlites. Light REE enrichment in fenites and the granular suite of nodules from kimberlites attests to the mobility of REE in CO₂-rich fluids under both mantle and crustal conditions.     

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     

Carbonate xenoliths in La Palma: Carbonatite or alteration product?

Carbonate xenoliths containing olivine and rimmed by kaersutitic amphibole were collected in basaltic rocks of the Basal Complex of La Palma. The mineralogical composition and microscopic appearance may suggest a relationship with carbonatites in general, thus a major element, trace element and stable isotope study was conducted to investigate the origin of the carbonate formation. Based on electron microprobe analyses, the carbonate is calcite with up to 6.3 wt% MgO and 7.2 wt% SiO₂. The elevated SiO₂ content may suggest a melt origin for the carbonate. However, the C and O isotope compositions of the carbonate xenoliths (δ¹³C and δ¹⁸O around −1‰ and 13‰, respectively) are similar to those of calcite veins and amygdales in basaltic rocks of the Basal Complexes of La Palma and Fuerteventura and are interpreted as produced by fluid degassing and metasomatism by CO₂-H₂O fluid derived from mobilization of sedimentary material. Trace element contents determined by laser-ablation ICP-MS analyses support the assumed origin, thus, the relationship with carbonatitic melts can be excluded. Based on trace element compositions, the amphibole surrounding the xenoliths is not related genetically to the carbonate. The elevated SiO₂ content of the calcite can be attributed to submicron relics of pyroxene, thus, the use of this feature as an evidence for melt origin is questionable.     

Martian regolith in Elephant Moraine 79001 shock melts? Evidence from major element composition and sulfur speciation

Shock veins and melt pockets in Lithology A of Martian meteorite Elephant Moraine (EETA) 79001 have been investigated using electron microprobe (EM) analysis, petrography and X-ray Absorption Near Edge Structure (XANES) spectroscopy to determine elemental abundances and sulfur speciation (S²⁻ versus S⁶⁺). The results constrain the materials that melted to form the shock glasses and identify the source of their high sulfur abundances. The XANES spectra for EETA79001 glasses show a sharp peak at 2.471 keV characteristic of crystalline sulfides and a broad peak centered at 2.477 keV similar to that obtained for sulfide-saturated glass standards analyzed in this study. Sulfate peaks at 2.482 keV were not observed. Bulk compositions of EETA79001 shock melts were estimated by averaging defocused EM analyses. Vein and melt pocket glasses are enriched in Al, Ca, Na and S, and depleted in Fe, Mg and Cr compared to the whole rock. Petrographic observations show preferential melting and mobilization of plagioclase and pyrrhotite associated with melt pocket and vein margins, contributing to the enrichments. Estimates of shock melt bulk compositions obtained from glass analyses are biased towards Fe- and Mg- depletions because, in general, basaltic melts produced from groundmass minerals (plagioclase and clinopyroxene) will quench to a glass, whereas ultramafic melts produced from olivine and low-Ca pyroxene megacrysts crystallize during the quench. We also note that the bulk composition of the shock melt pocket cannot be determined from the average composition of the glass but must also include the crystals that grew from the melt - pyroxene (En_72-75Fs_20-21Wo_5-7) and olivine (Fo_75-80). Reconstruction of glass + crystal analyses gives a bulk composition for the melt pocket that approaches that of lithology A of the meteorite, reflecting bulk melting of everything except xenolith chromite.Our results show that EETA79001 shock veins and melt pockets represent local mineral melts formed by shock impedance contrasts, which can account for the observed compositional anomalies compared to the whole rock sample. The observation that melts produced during shock commonly deviate from the bulk composition of the host rock has been well documented from chondrites, rocks from terrestrial impact structures and other Martian meteorites. The bulk composition of shock melts reflects the proportions of minerals melted; large melt pockets encompass more minerals and approach the whole rock whereas small melt pockets and thin veins reflect local mineralogy. In the latter, the modal abundance of sulfide globules may reach up to 15 vol%. We conclude the shock melt pockets in EETA79001 lithology A contain no significant proportion of Martian regolith.     

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.     

Origin of CO2 and carbonate veins in mantle-derived xenoliths in the Pannonian Basin

The origin and evolution of CO₂ inclusions and calcite veins in peridotite xenoliths of the Pannonian Basin, Hungary, were investigated by means of petrographic investigation and stable isotope analyses. The fluid inclusions recovered in paragenetic olivine and clinopyroxene belong to distinct populations: type A (texturally early) inclusions with regular shapes (often with negative crystal forms) forming intragranular trails; type B (texturally late) inclusions defining randomly oriented trails that reach grain boundaries. Type B inclusions are often associated with silicate melt (type C) inclusions. Stable carbon isotope compositions in inclusion-hosted CO₂ were obtained by vacuum crushing followed by conventional dual inlet as well as continuous flow mass spectrometry in order to eliminate possible lab artifacts. Olivines, clino- and orthopyroxenes of the host peridotite have oxygen isotope compositions from 5.3 to 6.0‰ (relative to V-SMOW), without any relationship with xenolith texture. Some of the xenoliths contained calcite in various forms: veins and infillings in silicate globules in veins, secondary carbonate veins filling cracks and metasomatic veins with diffuse margins. The former two carbonate types have δ¹³C values around –13‰ (relative to V-PDB) and low Sr contents (< 0.5 wt.%), whereas the third type,veins with high-temperature metasomatic features have a δ¹³C value of –5.0‰ and high Sr contents up to 3.4 wt.%. In spite of the mantle-like δ¹³C value and the unusually high Sr content typical for mantle-derived carbonate, trace element compositions have proven a crustal origin. This observation supports the conclusions of earlier studies that the carbonate melt droplets found on peridotite xenoliths in the alkaline basalts represent mobilized sedimentary carbonate. The large δ¹³C range and the ¹²C-enrichment in the carbonates can be attributed to devolatilization of the migrating carbonate or infiltration of surficial fluids containing ¹²C-rich dissolved carbon.Carbon isotope compositions of inclusion-hosted CO₂ range from –17.8 to –4.8‰ (relative to V-PDB) with no relation to the amount of CO₂ released by vacuum crushing. Low-δ¹³C values measured by stepwise heating under vacuum suggest that the carbon component is pristine and not related to surficial contamination, and that primary mantle fluids with δ¹³C values around –5‰ were at least partly preserved in the xenoliths. Tectonic reworking and heating by the basaltic magma resulted in partial CO₂ release and local ¹³C-depletion.     

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.     

Carbonatite and highly peralkaline nephelinite melts from Oldoinyo Lengai Volcano,Tanzania: The role of natrite-normative fluid degassing

Oldoinyo Lengai, located in the Gregory Rift in Tanzania, is a world-famous volcano owing to its uniqueness in producing natrocarbonatite melts and because of its extremely high CO₂ flux. The volcano is constructed of highly peralkaline [PI = molar (Na₂O + K₂O)/Al₂O₃ > 2–3] nephelinite and phonolites, both of which likely coexisted with carbonate melt and a CO₂-rich fluid before eruption. Results of a detailed melt inclusion study of the Oldoinyo Lengai nephelinite provide insights into the important role of degassing of CO₂-rich vapor in the formation of natrocarbonatite and highly peralkaline nephelinites. Nepheline phenocrysts trapped primary melt inclusions at 750–800 °C, representing an evolved state of the magmas beneath Oldoinyo Lengai. Raman spectroscopy, heating-quenching experiments, low current EDS and EPMA analyses of quenched melt inclusions suggest that at this temperature, a dominantly natrite_ss-normative, F-rich (7–14 wt%) carbonate melt and an extremely peralkaline (PI = 3.2–7.9), iron-rich nephelinite melt coexisted following degassing of a CO₂ + H₂O-vapor. We furthermore hypothesize that the degassing led to re-equilibration between the melt and liquid phases that remained and involved 1/ mixing between the residual (after degassing) alkali carbonate liquid and an F-rich carbonate melt and 2/ enrichment of the coexisting nephelinite melt in alkalis. We suggest that in the geological past similar processes were responsible for generating highly peralkaline silicate melts in continental rift tectonic settings worldwide.     

Experimental determination of oxygen isotope fractionations between CO2 vapor and soda-melilite melt

We report results of experiments constraining oxygen isotope fractionations between CO₂ vapor and Na-rich melilitic melt at 1 bar and 1250 and 1400°C. The fractionation factor constrained by bracketed experiments, 1000^.lnα_{CO2-Na melilitic melt}, is 2.65±0.25 ‰ (±2σ; n=92) at 1250°C and 2.16±0.16 ‰ (2σ; n=16) at 1400°C. These values are independent of Na content over the range investigated (7.5 to 13.0 wt. % Na₂O). We combine these data with the known reduced partition function ratio of CO₂ to obtain an equation describing the reduced partition function ratio of Na-rich melilite melt as a function of temperature. We also fit previously measured CO₂-melt or -glass fractionations to obtain temperature-dependent reduced partition function ratios for all experimentally studied melts and glasses (including silica, rhyolite, albite, anorthite, Na-rich melilite, and basalt). The systematics of these data suggest that reduced partition function ratios of silicate melts can be approximated either by using the Garlick index (a measure of the polymerization of the melt) or by describing melts as mixtures of normative minerals or equivalent melt compositions. These systematics suggest oxygen isotope fractionation between basalt and olivine at 1300°C of approximately 0.4 to 0.5‰, consistent with most (but not all) basalt glass-olivine fractionations measured in terrestrial and lunar basalts.     

Searching for parental kimberlite melt

Constraining the composition of primitive kimberlite magma is not trivial. This study reconstructs a kimberlite melt composition using vesicular, quenched kimberlite found at the contact of a thin hypabyssal dyke. We examined the 4 mm selvage of the dyke where the most elongate shapes of the smallest calcite laths suggest the strongest undercooling. The analyzed bulk compositions of several 0.09-1.1 mm² areas of the kimberlite free from macrocrysts were considered to be representative of the melt. The bulk analyses conducted with a new “chemical point-counting” technique were supplemented by modal estimates, studies of mineral compositions, and FTIR analysis of olivine phenocrysts. The melt was estimated to contain 26-29.5 wt% SiO₂, ∼7 wt% of FeO_T, 25.7-28.7 wt% MgO, 11.3-15 wt% CaO, 8.3-11.3 wt% CO₂, and 7.6-9.4 wt% H₂O. Like many other estimates of primitive kimberlite magma, the melt is too magnesian (Mg# = 0.87) to be in equilibrium with the mantle and thus cannot be primary. The observed dyke contact and the chemistry of the melt implies it is highly fluid (η = 10¹-10³ Pa s at 1100-1000 °C) and depolymerized (NBO/T = 2.3-3.2), but entrains with 40-50% of olivine crystals increasing its viscosity. The olivine phenocrysts contain 190-350 ppm of water suggesting crystallization from a low SiO₂ magma (a_SiO2 below the olivine-orthopyroxene equilibrium) at 30-50 kb. Crystallization continued until the final emplacement at depths of few hundred meters which led to progressively more Ca- and CO₂-rich residual liquids. The melt crystallised phlogopite (6-10%), monticellite (replaced by serpentine, ∼10%), calcite rich in Sr, Mg and Fe (19-27%), serpentine (29-31%) and minor amounts of apatite, ulvöspinel-magnetite, picroilmenite and perovskite. The observed content of H₂O can be fully dissolved in the primitive melt at pressures greater than 0.8-1.2 kbar, whereas the amount of primary CO₂ in the kimberlite exceeds CO₂ soluble in the primitive kimberlite melt. A mechanism for retaining CO₂ in the melt may require a separate fluid phase accompanying kimberlite ascent and later dissolution in residual carbonatitic melt. Deep fragmentation of the melt as a result of volatile supersaturation is not inevitable if kimberlite magma has an opportunity to evolve.     

Solubility of CO2 in olivine melilitite at high pressures and role of CO2 in the earth's upper mantle

The positions of the liquidi and the near-liquidus phases of olivine-melilitite+CO₂ have been determined under MH-buffered and furnace-buffered conditions up to 40 kb. It is found that CO₂ alone lowers the liquidus compared to dry conditions, yet its influence is minor compared to H₂O. The major role of CO₂ is to favour the growth of orthopyroxene and garnet over that of olivine at least at high pressures. CO₂-contents of glasses from experiments just above the liquidus (MH-buffered) were determined as 5.1 % at 10kb; 7.5 % at 20kb, 9.3 % at 30kb and 10–11 % (estimated) at 40 kb. Experiments on (pyrolite –40 % olivine)+H₂O+CO₂ show that CO₂ occurs under mantle conditions as carbonate under subsolidus conditions and dissolved in a melt above the solidus. At 30kb, the solidus lies between 1,000 ° C and 1,050 ° C for vapour-saturated conditions, at and at .     

Carbonate clumped isotope thermometry of deep-sea corals and implications for vital effects

Here we calibrate the carbonate clumped isotope thermometer in modern deep-sea corals. We examined 11 specimens of three species of deep-sea corals and one species of a surface coral spanning a total range in growth temperature of 2-25 °C. External standard errors for individual measurements ranged from 0.005‰ to 0.011‰ (average: 0.0074‰) which corresponds to ∼1-2 °C. External standard errors for replicate measurements of Δ₄₇ in corals ranged from 0.002‰ to 0.014‰ (average: 0.0072‰) which corresponds to 0.4-2.8 °C. We find that skeletal carbonate from deep-sea corals shows the same relationship of Δ₄₇ (the measure of ¹³C-¹⁸O ordering) to temperature as does inorganic calcite. In contrast, the δ¹³ C and δ¹⁸O values of these carbonates (measured simultaneously with Δ₄₇ for every sample) differ markedly from equilibrium with seawater; i.e., these samples exhibit pronounced ‘vital effects’ in their bulk isotopic compositions. We explore several reasons why the clumped isotope compositions of deep-sea coral skeletons exhibit no evidence of a vital effect despite having large conventional isotopic vital effects.     

Argon and CO2 on the race track in silicate melts: A tool for the development of a CO2 speciation and diffusion model

We have analysed the kinetics of Argon and CO₂ diffusion in simplified iron free rhyolitic to hawaiitic melts using the diffusion couple technique. The concentration distance profiles of Ar and CO₂ were measured with electron microprobe analysis and Fourier Transform Infrared Spectroscopy, respectively. Error functions were fitted to the symmetrical concentration distance profiles to extract the diffusion coefficients.In the temperature range 1373 to 1773 K the activation energies for Ar diffusion range from 169 ± 20 to 257 ± 62 kJ mol⁻¹. Ar diffusivity increases exponentially with the degree of depolymerisation. In contrast, the mobility of total CO₂, that is identical to Ar mobility in rhyolitic melt, keeps constant with changing bulk composition from rhyolite to hawaiite. CO₂ speciation at 1623 K and 500 MPa was modeled for the range of compositions studied using the diffusion data of Ar and total CO₂ in combination with network former diffusion calculated from viscosity data. Within error this model is in excellent agreement with CO₂ speciation data extrapolated from temperatures near the glass transition temperature for dacitic melt composition. This model shows that even in highly depolymerised hawaiitic and tholeiitic melts molecular CO₂ is a stable species and contributes 70 to 80% to the total CO₂ diffusion, respectively.     

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.     

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.     

Generation and differentiation of group II kimberlites: constraints from a high-pressure experimental study to 10 GPa

Experiments have been performed in the multicomponent (natural) bulk system to constrain the conditions of generation and differentiation of a K-rich group II kimberlite (now also referred to as orangeite). The group II composition examined was saturated in olivine, orthopyroxene, and garnet at near liquidus conditions in the pressure range 4 to 10 GPa. In the range 2 to 3 GPa, the liquidus phase was olivine only. The potassic nature of the melts in the bulk compositions studied was ensured by the absence of any K-bearing phase in the residual assemblage at P > 4 GPa. Phlogopite is destabilized toward higher pressures by a carbonation reaction of the type phlogopite + CO₂ = enstatite + garnet + K₂CO₃ (liquid) + H₂O leading to alkalic, carbonatitic liquids coexisting with a garnet-peridotite (harzburgite or lherzolite) residue over a wide pressure-temperature space at pressures in excess of 4 GPa. Evidently, CO₂-bearing systems do not favor the stability of phlogopite and/or K-richterite amphibole at pressures in excess of 4 to 5 GPa, and it is suggested that the carbonate-bearing and potassic character of any mantle melt originating from this depth is most likely the product of a two-stage process: either a carbonate-bearing protolith is invaded by a potassic melt or fluid (probably supercritical), or a potassic protolith (after metasomatism) has been invaded by a carbonatite melt.