Crystal-melt partitioning of noble gases (helium, neon, argon, krypton, and xenon) for olivine and clinopyroxene

Mineral-melt partition coefficients of all noble gases (^min/meltD_i) have been obtained for olivine (ol) and clinopyroxene (cpx) by UV laser ablation (213 nm) of individual crystals grown from melts at 0.1 GPa mixed noble gas pressure. Experimental techniques were developed to grow crystals virtually free of melt and fluid inclusions since both have been found to cause profound problems in previous work. This is a particularly important issue for the analysis of noble gases in crystals that have very low partition coefficients relative to coexisting melt and fluid phases. The preferred partitioning values obtained for the ol-melt system for He, Ne, Ar, Kr, and Xe are 0.00017(13), 0.00007(7), 0.0011(6), 0.00026(16), and , respectively. The respective cpx-melt partition coefficients are 0.0002(2), 0.00041(35), 0.0011(7), 0.0002(2), and . The data confirm the incompatible behaviour of noble gases for both olivine and clinopyroxene but unlike other trace elements these values show little variation for a wide range of atomic radius. The lack of dependence of partitioning on atomic radius is, however, consistent with the partitioning behaviour of other trace elements which have been found to exhibit progressively lower dependence of ^min/meltD_i on radius as the charge decreases. As all noble gases appear to exhibit similar ^min/meltD_i values we deduce that noble gases are not significantly fractionated from each other by olivine and clinopyroxene during melting and fractional crystallisation. Although incompatible, the partitioning values for noble gases also suggest that significant amounts of primordial noble gases may well have been retained in the mantle despite intensive melting processes. The implication of our data is that high primordial/radiogenic noble gas ratios (³He/⁴He, ²²Ne/²¹Ne, and ³⁶Ar/⁴⁰Ar) characteristic of plume basalt sources can be achieved by recycling a previously melted (depleted) mantle source rather than reflecting an isolated, non-degassed primordial mantle region.     

Solubilities of nitrogen and noble gases in silicate melts under various oxygen fugacities: implications for the origin and degassing history of nitrogen and noble gases in the earth

Solubility experiments for nitrogen and noble gases (Ar and Ne) in silicate melts were conducted using two experimental configurations: one was conducted at 1 atmospheric pressure, T =1300°C and oxygen fugacity (fO₂) of IW + 0.9 (i.e., 0.9 log units higher than the iron-wüstite buffer) and the other at high pressures (P_total ∼ 2 × 10⁸ Pa), 1500°C and fO₂ ∼ IW + 6. For the former experiment, isotopically labeled-nitrogen (¹⁵N¹⁵N-enriched) was used to distinguish dissolved nitrogen from contaminating atmospheric or organic nitrogen and to examine dissolution mechanisms of nitrogen in silicate melts. The results obtained for the two series of experiments are consistent with each other, suggesting that Henry's law is satisfied for fN₂ of up to ∼250 atm (2.5 × 10⁷ Pa). The results are also consistent with our earlier results (Miyazaki et al., 1995) obtained at highly oxidizing conditions (fO₂ ∼ IW + 10). All these results support physical dissolution of nitrogen as N₂ molecules in silicate melts for fO₂ from ∼IW + 10 down to ∼IW. The observed solubility (Henry's constant) of nitrogen (3-5 × 10⁻⁹ mol/g/atm) is comparable to that of Ar (2-4 × 10⁻⁹ mol/g/atm), and much lower than that of Ne (11-14 × 10⁻⁹ mol/g/atm) at 1300°C. A preliminary experiment was also performed for partitioning of nitrogen and noble gases between clinopyroxene (cpx) and basaltic melt using a piston cylinder-type apparatus at 1.5 GPa and at 1270 to 1350°C. The obtained cpx/melt partition coefficient of nitrogen is 0.06, slightly lower than those of noble gases (∼0.1 for Ne to Xe), suggesting that nitrogen is as incompatible as or even slightly more incompatible than noble gases. The present results imply that a large nitrogen/Ar fractionation would not be produced by magmatic processes. Therefore, the two orders of magnitude difference between the N₂/³⁶Ar ratios in the Earth's atmosphere (∼10⁴) and that in the mantle (∼10⁶) must be explained by some other processes, such as incomplete segregation of metal blobs into the core and their later oxidation.     

Low solubility of He and Ar in carbonatitic liquids: Implications for decoupling noble gas and lithophile isotope systems

In order to asses the importance of carbonatitic liquids in transporting noble gases in the mantle, the solubilities of He and Ar in carbonatitic liquids were estimated from analyses of calcium-potassium carbonate glasses that had been synthesized at 1 bar and temperatures between 850 and 950 °C under He or Ar enriched atmospheres. Despite poor reproducibility related to difficulties synthesizing carbonatite glass, we have been able to estimate He and Ar solubilities in carbonatite liquids to be 1 × 10⁻⁸ and 2 × 10⁻⁹ mol g⁻¹ at 1 bar respectively (with ?50% uncertainty). Despite the significant uncertainties on these estimates, it is clear that the noble gases are not massively soluble in carbonatite liquids (within error, these solubilities are identical to their equivalent solubilities in tholeiitic melts). Assuming the results of these low pressure experiments can be applied to mantle conditions, it seems unlikely that carbonatite metasomatism per se transports noble gases within the mantle. It is nevertheless possible that partitioning of lithophile trace elements (including the important radioelements, U, K and Th) and noble gases between a carbonatitic melt and a silicate melt could effectively decouple lithophile and noble gas isotope systematics because the carbonatitic melt expressedly does not transport noble gases, yet is known to efficiently transport incompatible trace elements.     

Solubilities of noble gases in magnetite: implications for planetary gases in meteorites

Solubilities of noble gases in magnetite were determined by growing magnetite in a noble-gas atmosphere between 450 and 700°K. Henry's law is obeyed at pressures up to 10^?2 atm for He, Ne, Ar and up to 10^?5 atm for Kr, Xe, with the following distribution coefficients at 500° (cc STP g^?1 atm^?5): He 0.042, Ne 0.016, Ar 3.6, Kr 1.3, Xe 0.88, some 10²–10⁵ times higher than previous determinations on silicate and fluoride melts. Apparent heats of solution in kcal/mole are: He ?2.42 ±0.12, Ne ?2.20 ±0.10, Ar ?15.25 ±0.25, Kr ?13.0 ±0.3, Xe ?12-5 ± 0.5. These values, too, stand in sharp contrast with earlier determinations on melts which were small and positive, but are comparable to the values for clathrates. Presumably the gases are held in anion vacancies.Extrapolation of the magnetite data to the formation temperature of C1 chondrites, 360°K, shows that the Ar_p³⁶ content of Orgueil magnetite could be acquired by equilibrium solubility at a total nebular pressure of 4 × 10^?6 atm. In the absence of data for silicates (the principal host phase of planetary gas), an attempt is made to estimate the solubilities required to account for planetary gases in meteorites. These values do not appear grossly unreasonable in the light of the magnetite data, when structural differences between the two minerals are taken into account. It seems that equilibrium solubility may be able to account for four features of planetary gas: elemental ratios, amounts, correlations with other volatiles and retentive siting. It cannot account for the isotopic fractionation of planetary gas, however.     

Dual speciation of nitrogen in silicate melts at high pressure and temperature: An experimental study

Solubility and speciation of nitrogen in silicate melts have been investigated between 1400 and 1700 °C and at pressures ranging from 10 to 30 kbar for six different binary alkali and alkaline-earth silicate liquids and a Ca-Mg-alumino silicate. Experiments were performed in a piston-cylinder apparatus. The nitrogen source is silver azide, which breaks down to Ag and molecular N₂ below 300 °C. At high pressure and temperature, the nitrogen content may be as high as 0.7 wt% depending on the melt composition, pressure, and temperature. It increases with T, P and the polymerization state of the liquid. Characterization by Raman spectroscopy and ¹⁵N solid state MAS NMR indicates that nitrogen is not only physically dissolved as N₂ within the melt structure like noble gases, but a fraction of nitrogen interacts strongly with the silicate network. The most likely nitrogen-bearing species that can account for Raman and NMR results is nitrosyl group. Solubility data follow an apparent Henry’s law behavior and are in good agreement with previous studies when the nitrosyl content is low. On the other hand, a significant departure from a Henry’s law behavior is observed for highly depolymerized melts, which contain more nitrosyl than polymerized melts. Possible solubility mechanisms are also discussed. Finally, a multi-variant empirical relation is given to predict the relative content of nitrosyl and molecular nitrogen as a function of P, T, and melt composition and structure. This complex speciation of nitrogen in melts under high pressure may have significant implication concerning crystal-melt partitioning of nitrogen as well as for potential elemental and isotopic fractionation of nitrogen in the deep Earth.     

Sorption of noble gases by solids,with reference to meteorites. III. Sulfides,spinels, and other substances; on the origin of planetary gases

To simulate trapping of noble gases by meteorites, we reacted 15 FeCr or FeCrNi alloy samples with CO, H₂O or H₂S at 350–720 K, in the presence of noble gases. The reaction products, including (Fe,Cr)₂O₃, FeCr₂S₄, FeS, C, and Fe₃C, were analyzed by mass spectrometry, usually after chemical separation by selective solvents. Three carbon samples were prepared by catalytic decomposition of CO or by dehydration of carbohydrates with H₂SO₄.The spinel and carbon samples were similar to those of earlier studies (Yang et al., 1982 and Yang and Anders, 1982), with only minor effects attributable to the presence of Ni. All samples sorted substantial amounts of noble gases, with distribution coefficients of 10^?1–10^?2 cm³ STP/g atm for Xe. On the basis of release temperature three gas components were distinguished: a generally dominant physisorbed component (20–80% of total), and two more strongly bound, chemisorbed and trapped components. Judging from the elemental pattern, the adsorbed components were acquired at the highest noble gas partial pressure encountered by the sample—atmosphere or synthesis vessel.Sulfides, particularly daubréelite, showed three distinctive trends relative to chromite or magnetite: the high-T component was larger, 30–70% of the total; $\text{Ne}\text{Xe}$ ratios were higher, by up to 10², possibly due to preferential diffusion of Ne during synthesis. In one synthesis, at relatively high P, the gases were sorbed with only minimal elemental fractionation, presumably by occlusion.Most of the features of primordial noble gases can be explained in terms of the data and concepts presented in the three papers of this series. The elemental fractionation pattern of Ar, Kr, Xe in meteorites, terrestrial rocks, and planets resembles the adsorption pattern on the solids studied: carbon, spinels, Sulfides, etc. The variation in $\text{Ne}\text{Ar}$ ratio may be explained by preferential diffusion of Ne. The high release temperature of meteoritic noble gases may be explained by transformation of physisorbed to chemisorbed gas, as observed in some experiments. The ready loss of meteoritic heavy gases on surficial oxidation (“Phase Q”) is consistent with adsorption, as is the high abundance. Extrapolation of the limited laboratory data suggests that the observed amounts of noble gases could have been adsorbed from a solar gas at 160–170 K and 10^?6–10^?5 atm, i.e. in the early contraction stages of the solar nebula. The principal unsolved problem is the origin of isotopically anomalous, apparently mass-fractionated noble gases in the Earth's atmosphere and in meteoritic carbon and chromite.     

The solubility of noble gases in water and in NaCl brine

A direct-sampling, mass-spectrometric technique has been used to measure simultaneously the solubilities of He, Ne, Ar, Kr, and Xe in fresh water and NaCl brine (0 to 5.2 molar) from 0° to 65 °C, and at 1 atm total pressure of moist air. The argon solubility in the most concentrated brines is 4 to 7 times less than in fresh water at 65 °C and 0°C, respectively. The salt effect is parameterized using the Setschenow equation.

$\text{ln}\text{[}\text{β}{}^{\mathrm{i}}{}_{\mathrm{o}}\text{(T)}\text{β}{}^{\mathrm{i}}\text{(T)}\text{= MK}{}^{\mathrm{i}}{}_{\mathrm{M}}\text{(T)}$

where M is NaCl moiarity, βⁱ_o(T) and βⁱ(T) the Bunsen solubility coefficients for gas i in fresh water and brine, and Kⁱ_M(T) the empirical salting coefficient. Values of Kⁱ_M(T) are calculated using volumetric concentration units for noble gas and NaCl content and are independent of NaCl molarity. Below about 40°C, temperature coefficients of all Kⁱ_M are negative. The value of K^He_M is a minimum at 40°C. K^Ar_M decreases from about 0.40 at 0°C to 0.28 at 65 °C. The absolute magnitudes of the differences in salting coefficients (relative to K^Ar_M) decrease from 0° to 65°C. Over the range of conditions studied, all noble gases are salted out, and K^He_M ? K^Ne_M < K^Ar_M < K^Kr_M < K^Xe_M.From the solubility data, we calculated $\text{Δ}\text{G}{}^0{}_{\mathrm{tr}}$, $\text{Δ}\text{S}{}^0{}_{\mathrm{tr}}$, $\text{Δ}\text{H}{}^0{}_{\mathrm{tr}}$ and $\text{Δ}\text{C}{}^{\mathrm{O}}{}_{\mathrm{p,tr}}$ for the transfer of noble gases from fresh water to 1 molar NaCl solutions. At low temperatures $\text{Δ}\text{S}{}^0{}_{\mathrm{tr}}$, is positive, but decreases and becomes negative at temperatures ranging from about 25°C for He to 45°C for Xe. At low temperatures, the dissolved electrolyte apparently interferes with the formation of a cage of solvent molecules about the noble gas atom. At higher temperatures, the local environment of the gas atom in the brine appears to be slightly more ordered than in pure water, possibly reflecting the longer effective range of the ionic fields at higher temperature.The measured solubilities can be used to model noble gas partitioning in two-phase geothermal systems at low temperatures. The data can also be used to estimate the temperature and concentration dependence of the salt effect for other alkali halides. Extrapolation of the measured data is not possible due to the incompletely-characterized minima in the temperature dependence of the salting coefficients. The regularities in the data observed at low temperatures suggest relatively few high-temperature data will be required to model the behavior of noble gases in high-temperature geothermal brines.     

Gold solubility in oxidized and reduced, water-saturated mafic melt

We have performed experiments to evaluate Au solubility in natural, water-saturated basaltic melts as a function of oxygen fugacity. Experiments were carried out at 1000 °C and 200 MPa, and oxygen fugacity was controlled at the fayalite-magnetite-quartz (FMQ) oxygen fugacity buffer and FMQ + 4. All experiments were saturated with a metal-chloride aqueous solution loaded initially as a 10 wt% NaCl eq. fluid. The stable phase assemblage at FMQ consists of basalt melt, olivine, clinopyroxene, a single-phase aqueous fluid, and metallic Au. The stable phase assemblage at FMQ + 4 consists of basalt melt, clinopyroxene, magnetite-spinel solid solution, a single-phase aqueous fluid, and metallic Au. Silicate glasses (i.e., quenched melt) and their contained crystalline material were analyzed by using both electron probe microanalysis (EPMA) and laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS). Measured Au concentrations in the quenched melt range from 4.8 μg g⁻¹ to 0.64 μg g⁻¹ at FMQ + 4, and 0.54 μg g⁻¹ to 0.1 μg g⁻¹ at FMQ. The measured solubility of Au in olivine and clinopyroxene was consistently below the LA-ICP-MS limit of detection (i.e., 0.1 μg g⁻¹). These melt solubility data place important limitations on the dissolved Au content of water-saturated, Cl- and S-bearing basaltic liquids at geologically relevant fO₂ values. The new data are compared to published, experimentally-determined values for Au solubility in dry and hydrous silicate liquids spanning the compositional range from basalt to rhyolite, and the effects of melt composition, oxygen fugacity, pressure and temperature are discussed.     

Noble gas distribution in oceanic basalt glasses

The distribution of noble gases has been investigated in six MORB glass samples using stepwise heating, vacuum crushing and the analysis of grain size fractions. These experiments indicate a strong noble gas partitioning into CO₂-filled vesicles. An inhomogeneous distribution of argon isotopes within individual glasses is observed in several cases and is believed to result from the combined effects of the deep-seated component partitioning into vesicles and the contamination of the residual component dissolved in the glass by atmospheric noble gases. Using a mixing diagram, we are able to discriminate against atmospheric contamination and estimate the Ar partition coefficients. These coefficients are in qualitative agreement with the Henry's law constant for silicate melts.     

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.     

The effect of compression on noble gas solubility in silicate melts and consequences for degassing at mid-ocean ridges

The effect of compression on noble gas solubility in silicate melts is still badly understood due to a lack of theoretical guidance. In the experimental literature, noble gases dissolving in liquid silicates are found to concentrate almost linearly with increasing pressure up to several tens of kbar, suggesting that Henry’s law could be valid up to very high pressures, although this law stipulates that the gaseous phase in contact with the liquid must be ideal. Recently, new experiments dealing with the dissolution of argon in synthetic and natural silicate melts have pointed out that the evolution of concentration with pressure exhibits a departure from linearity in the 50-100 kbar range, leading either to a levelling off or to a sudden collapse of the argon concentration above 50 kbar. Here, we investigate by means of liquid state physics how volatile species dissolve into silicate melts under pressure. We use a hard sphere model (the reference fluid in liquid state physics) to describe silicate melts and gas at high pressures. One of our main results is that, when pressure increases, the concurrent compaction of gas and melt explains the almost-linear behaviour of the noble gas concentration up to several tens of kbars, before melt compaction dominates and concentration either levels off or decreases gradually in the 50-100 kbar range. In spite of the existence of a quasi-linear regime over a large pressure range, our work disqualifies the use of the Henry law when dealing with high pressures. The implication of these findings to provide an understanding of degassing at mid-ocean ridges is next investigated. Applying our model to the scenario where CO₂ vesicle generation occurs in the magma at mantle depths during its ascent from melting regions, we evaluate magma vesicularity as well as noble gas concentrations in the basalt melt and in vesicles as a function of pressure at depth. It is stressed that the variable and usually strong noble gas elemental fractionation observed in mid-ocean ridge basalts can be explained by assuming a sequence of several vesiculation stages interrupted by vesicle loss during magma ascent.     

Heavy metal and sulfur transport during subcritical and supercritical hydrothermal alteration of basalt: Influence of fluid pressure and basalt composition and crystallinity

Crystalline basalt, diabase and basalt glass have been reacted with a Na-Ca-K-Cl fluid of seawater ionic strength at 350–425°C, 375–400 bars pressure and fluid/rock mass ratios of 0.5–1.0, to assess the role of temperature, basalt/diabase chemistry and texture on heavy metal and sulfur mobility during hydrothermal alteration.Alteration of basalt/diabase is characterized by cation fixation and hydrolysis reactions which show increased reaction progress with increasing temperature at constant pressure. Correspondingly, pH in a series of 400 bar experiments ranges from 4.8 to 2.7 at 350 and 425°C, respectively and is typically lower for alteration of a SiO₂-rich crystalline basalt than for other rock types, due, in part, to relatively high SiO₂ concentrations in solution. High SiO₂ concentrations stabilize hydrous Na- and Ca-rich alteration phases, causing pH to decrease according to reactions such as: 3.0 CaAl₂Si₂O₈ + 1.0 Ca⁺⁺ + 2.0 H₂O = 2.0 Ca₂Al₃Si₃O₁₂(OH) + 2.0 H⁺Phases experimentally produced include: mixed layer chlorite/smectite, Ca-rich amphibole and clinozoisite. Clinozoisite was identified as a replacement product of plagioclase from diabase-solution interaction experiments.In direct response to H⁺ production, dissolved Fe, Mn and H₂S concentrations increase dramatically. For early-stage reaction, H₂S typically exceeds Fe and Mn. However, at 425°C and after long-term reaction at 400°C, H₂S is lost from solution, apparently in response to pyrite replacement of oxide and silicate phases.Pyrrhotite formed at temperatures ≤ 375°C, whereas magnetite was identified in all run products, except from basalt glass alteration.Cu and Zn concentrations in solution are not simple functions of pH. These metals achieve greatest solubility in fluids from experiments at 375–400°C, except when basalt glass is used as a reactant. The relatively low concentrations of these species in solution during basalt glass reaction may be due to adsorption by fine grained alteration phases.     

Methane solubilities in multisalt solutions

Solubilities of methane in multisalt solutions at 550 psia and 25°C can be predicted from single-salt salting coefficients. The ionic strength contribution of the ith salt, I_i, is multiplied by its molal salting coefficient, k_mi, in the following summation over all salts:

$\text{log}\text{M}{}_{\mathrm{o}}\text{M}{}_{\mathrm{s}}\text{= ∑}{}_{\mathrm{i}}\text{k}{}_{\mathrm{m}}{}_{\mathrm{i}}\text{I}{}_{\mathrm{i}}$

where m_o and m_s are molal methane solubilities in distilled water and the salt solution, respectively, at the T, P and methane fugacity of interest.This equation predicts methane solubility in multisalt brines containing Na⁺, K⁺, Mg⁺², Ca⁺², Cl^?, SO₄^?2 and CO₃^?2 ions. k_mi values reported by Stoessell and Byrne (1982b) can be used in solubility predictions in brines at earth surface conditions. Prediction in reservoir brines would require determination of k_mi, for the different salts at reservoir temperatures and pressures.     

Noble gases in ultramafic xenoliths from San Carlos,Arizona

This work reports the results of noble gas (Ne, Ar, Kr, Xe) analyses of accidental mantle xenoliths from San Carlos, Arizona. Except for the addition of radiogenic ⁴⁰Ar and mass fractionation effects, the isotopic structures of these gases are indistinguishable from atmospheric composition. The absence of ¹²⁹Xe excesses in these rocks may reflect indirect mixing of atmospheric gases with the source region of the xenoliths. The dominant influence on the noble gas abundances in the San Carlos xenoliths appears to have been diffusive gas loss, which may have occurred in a mantle metamorphic event or during contact with the host basanite magma. Evidence is presented for the partitioning of significant amounts of the heavy noble gases into fluid inclusions in the xenolith minerals; the proportion of each gas in the inclusions increases with increasing atomic weight of the gas, possibly reflecting solubility effects. The noble gases are present in greater concentration in pyroxenes than in olivine, similar to the behavior of other incompatible elements.     

Chronology and shock history of the Bencubbin meteorite: A nitrogen, noble gas, and Ar-Ar investigation of silicates, metal and fluid inclusions

We have investigated the distribution and isotopic composition of nitrogen and noble gases, and the Ar-Ar chronology of the Bencubbin meteorite. Gases were extracted from different lithologies by both stepwise heating and vacuum crushing. Significant amounts of gases were found to be trapped within vesicles present in silicate clasts. Results indicate a global redistribution of volatile elements during a shock event caused by an impactor that collided with a planetary regolith. A transient atmosphere was created that interacted with partially or totally melted silicates and metal clasts. This atmosphere contained ¹⁵N-rich nitrogen with a pressure ?3 × 10⁵ hPa, noble gases, and probably, although not analyzed here, other volatile species. Nitrogen and noble gases were re-distributed among bubbles, metal, and partly or totally melted silicates, according to their partition coefficients among these different phases. The occurrence of N₂ trapped in vesicles and dissolved in silicates indicates that the oxygen fugacity (f_O2) was greater than the iron-wüstite buffer during the shock event. Ar-Ar dating of Bencubbin glass gives an age of 4.20 ± 0.05 Ga, which probably dates this impact event. The cosmic-ray exposure age is estimated at ∼40 Ma with two different methods. Noble gases present isotopic signatures similar to those of “phase Q” (the major host of noble gases trapped in chondrites) but elemental patterns enriched in light noble gases (He, Ne and Ar) relative to Kr and Xe, normalized to the phase Q composition. Nitrogen isotopic data together with ⁴⁰Ar/³⁶Ar ratios indicate mixing between a ¹⁵N-rich component (δ¹⁵N = +1000‰), terrestrial N, and an isotopically normal, chondritic N.Bencubbin and related ¹⁵N-rich meteorites of the CR clan do not show stable isotope (H and C) anomalies, precluding contribution of a nucleosynthetic component as the source of ¹⁵N enrichments. This leaves two possibilities, trapping of an ancient, highly fractionated atmosphere, or degassing of a primitive, isotopically unequilibrated, nitrogen component. Although the first possibility cannot be excluded, we favor the contribution of primitive material in the light of the recent finding of extremely ¹⁵N-rich anhydrous clasts in the CB/CH Isheyevo meteorite. This unequilibrated material, probably carried by the impactor, could have been insoluble organic matter extremely rich in ¹⁵N and hosting isotopically Q-like noble gases, possibly from the outer solar system.     

Partitioning of trace elements between carbonate-silicate melts and mantle minerals: Experiment and petrological consequences

The partitioning of a number of trace elements (Ba, Nb, Zr, Y, REE, etc.) between orthopyroxene, garnet, and carbonate-silicate melt was experimentally studied using a belt apparatus at pressures of 3.5–4.2 GPa and temperatures of 1300–1500°C. The experimental products were investigated by electron microprobe analysis and laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS). The experimental melts varied from carbonatitic (～5 wt % SiO₂) at low temperatures (1300–1350°C) to kimberlitic compositions (30 wt % SiO₂) at high temperatures (1500°C). The partition coefficients of most elements between orthopyroxene and melt (D _i ^Opx/L ) and garnet and melt (D _i ^Grt/L ) were almost independent of melt composition (temperature). The D _i ^Opx/L values ranged from <0.01 for the most incompatible Ba and light REE to 0.02–0.08 for moderately incompatible Zr, Y, and heavy REE. The D _i ^Grt/L values were approximately an order of magnitude higher, ～0.07 for light REE, 0.7 for Y, and 1.5 for Yb. The character of D _i ^Grt/L variations in the systems studied is in general similar to that established for silicate melts without volatile components. However, the differences in the behavior of moderately incompatible and compatible elements (e.g., light and heavy REE) in the experimental systems are less pronounced compared with CO₂-free systems. Considering carbonate-silicate and silicate melts as possible agents of mantle metasomatism, it can be concluded that the former can efficiently transport heavy REE, and the latter have a greater affinity for Nb, Ba, and light REE. A characteristic feature of mantle rocks enriched by carbonate-silicate melts is high Ba/La ratio coupled with relatively weakly fractionated REE distribution patterns. It was shown that the high degrees of enrichment observed in natural kimberlites can be explained by a two-stage scenario, including a preliminary invasion of carbonate-silicate melt into depleted harzburgites in the lower parts of the lithosphere and subsequent very low degree melting.     

Solubility of sulfur in some magmas at 1 atmosphere

In order to understand the distribution of sulfur in igneous rooks, we determined the solubility of sulfur in volcanic rock melts (tholeiite basalt, hawaiite and rhyodacite from Hawaii) at various gas compositions and at 1250° and 1300°C and 1 atm total pressure. The solubility of sulfur in the melt passes through a minimum with change in oxygen partial pressure, if other factors are held constant. For the basaltic liquid at 1200°C, most sulfur in the melt is as dissolved sulfide (S^?2) at oxygen partial pressures below 10^?8 atm and as dissolved sulfate at oxygen partial pressures above 10^?8 atm. Based on the present solubility data, 5 per cent is inferred for volcanic gas at 1 atm total pressure in equilibrium with subaerially extruded Hawaiian tholeiite basalt (Pele's hair with 180 ppm S) at 1200°C and 10^?8 atm P_O2.     

Noble-gas evidence for geothermal activity in a karstic terrain: Rocky Mountains,Canada.

Geothermally heated fluids are identified as a component in warm (up to 54°C) springs in the southern Canadian Rocky Mountains. High concentrations of radiogenic He (10^?3 ccSTP/cc gas) and atmospheric Ne, Ar, Kr and Xe in the gases that vigorously bubble at the Fairmont Hot Spring, assign the latter to be a “drowned” fumarole, fed by geothermal steam. Up to 75% depletions in the atmospheric noble gases in several warm springs indicate contributions of residual geothermal water. On the other hand, in a few cases noble gases were found in excess (W) over the expected concentrations in airequilibrated recharge water. The observed “reversed” pattern of w_nc >W_Ar >W_kr >W_xe is interpreted as excess air, incorporated during recharge through karstic conduits. The mixing ratios of geothermal and karstic waters can be deduced.     

Atomistic simulation of the mechanisms of noble gas incorporation in minerals

Atomistic simulations have been carried out to investigate the mechanisms of noble gas incorporation in minerals using both the traditional two-region approach and the “supercell” method. The traditional two-region approach has been used to calculate defect energies for Ne, Ar, Kr and Xe incorporation in MgO, CaO, diopside and forsterite in the static limit and at one atmosphere pressure. The possibilities of noble gas incorporation via both substitution and interstitial mechanisms are studied. The favored mechanism varies from mineral to mineral and from noble gas to noble gas. In all minerals studied, the variation of the solution energies of noble gas substitution with atomic radius appears approximately parabolic, analogous to those for 1⁺, 2⁺, 3⁺ and 4⁺ trace element incorporation on crystal lattice sites. Noble gas solution energies thus also fall on a curve, similar to those previously observed for cations with different charges, but with much lower curvature.The “supercell” method has been used to investigate the pressure dependence of noble gas incorporation in the same systems. Results indicate a large variation of the solubility of the larger noble gases, Kr and Xe with pressure. In addition, explicit simulation of incorporation at the (0 0 1) surface of MgO shows that the solubility of the heavier noble gases may be considerably enhanced by the presence of interfaces.     

A model for partitioning gases among brines and hydrocarbons in oil reservoirs: Examples from the IEA-GHG Weyburn CO2 Monitoring and Storage Project, Saskatchewan, Canada

We derive a thermodynamic model useful for calculating the partitioning of gases in brine–oil closed systems. We assumed real binary mixtures of all components, and solved a set of equilibrium equations relating the solubility of gases in brines to the solubility of gases in petroleum liquids. We applied the model to understand the distribution of gases in the IEA-GHG Weyburn CO₂ Monitoring and Storage Project, in Saskatchewan, Canada, throughout the first and second year after commencement of CO₂ injection. Results indicate that a free-gas phase formed in the south-east of the field 363 days after injection, which increased in size and volume with time, but did not correlate with zones showing excess CO₂ production. We also observed overlaps between areas of ¹³C depleted HCO₃⁻ and areas with free-vapor phase. We conclude that gas injection leads to short-term CO₂ dissolution in liquids, followed by the creation of a CO₂-rich vapor phase.