Volatile composition of microinclusions in diamonds from the Panda kimberlite, Canada: Implications for chemical and isotopic heterogeneity in the mantle

In order to better investigate the compositions and the origins of fluids associated with diamond growth, we have carried-out combined noble gas (He and Ar), C and N isotope, K, Ca and halogen (Cl, Br, I) determinations on fragments of individual microinclusion-bearing diamonds from the Panda kimberlite, North West Territories, Canada. The fluid concentrations of halogens and noble gases in Panda diamonds are enriched by several orders of magnitude over typical upper mantle abundances. However, noble gas, C and N isotopic ratios (³He/⁴He = 4-6 Ra, ⁴⁰Ar/³⁶Ar = 20,000-30,000, δ¹³C = −4.5‰ to −6.9‰ and δ¹⁵N = −1.2‰ to −8.8‰) are within the worldwide range determined for fibrous diamonds and similar to the mid ocean ridge basalt (MORB) source value. The high ³⁶Ar content of the diamonds (>1 × 10⁻⁹ cm³/g) is at least an order of magnitude higher than any previously reported mantle sample and enables the ³⁶Ar content of the subcontinental lithospheric mantle to be estimated at ∼0.6 × 10⁻¹² cm³/g, again similar to estimates for the MORB source. Three fluid types distinguished on the basis of Ca-K-Cl compositions are consistent with carbonatitic, silicic and saline end-members identified in previous studies of diamonds from worldwide sources. These fluid end-members also have distinct halogen ratios (Br/Cl and I/Cl). The role of subducted seawater-derived halogens, originally invoked to explain some of the halogen ratio variations in diamonds, is not considered an essential component in the formation of the fluids. In contrast, it is considered that large halogen fractionation of a primitive mantle ratio occurs during fluid-melt partitioning in forming silicic fluids, and during separation of an immiscible saline fluid.     

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.     

Grain boundary partitioning of Ar and He

An experimental procedure has been developed that permits measurement of the partitioning of Ar and He between crystal interiors and the intergranular medium (ITM) that surrounds them in synthetic melt-free polycrystalline diopside aggregates. ³⁷Ar and ⁴He are introduced into the samples via neutron irradiation. As samples are crystallized under sub-solidus conditions from a pure diopside glass in a piston cylinder apparatus, noble gases diffusively equilibrate between the evolving crystal and intergranular reservoirs. After equilibration, ITM Ar and He is distinguished from that incorporated within the crystals by means of step heating analysis. An apparent equilibrium state (i.e., constant partitioning) is reached after about 20 h in the 1450 °C experiments. Data for longer durations show a systematic trend of decreasing ITM Ar (and He) with decreasing grain boundary (GB) interfacial area as would be predicted for partitioning controlled by the network of planar grain boundaries (as opposed to ITM gases distributed in discrete micro-bubbles or melt). These data yield values of GB-area-normalized partitioning, , with units of (Ar/m³ of solid)/(Ar/m² of GB) of 6.8 × 10³-2.4 × 10⁴ m⁻¹. Combined petrographic microscope, SEM, and limited TEM observation showed no evidence that a residual glass phase or grain boundary micro-bubbles dominated the ITM, though they may represent minor components. If a nominal GB thickness (δ) is assumed, and if the density of crystals and the grain boundaries are assumed equal, then a true grain boundary partition coefficient , may be determined. For reasonable values of δ, is at least an order of magnitude lower than the Ar partition coefficient between diopside and melt. Helium partitioning data provide a less robust constraint with between 4 × 10³ and 4 × 10⁴ cm⁻¹, similar to the Ar partitioning data. These data suggest that an ITM consisting of nominally melt free, bubble free, tight grain boundaries can constitute a significant but not infinite reservoir, and therefore bulk transport pathway, for noble gases in fine grained portions of the crust and mantle where aqueous or melt fluids are non-wetting and of very low abundance (i.e., <0.1% fluid). Heterogeneities in grain size within dry equilibrated systems will correspond to significant differences in bulk rock noble gas content.     

The noble gas systematics of late-orogenic H2O-CO2 fluids, Mt Isa, Australia

The noble gases (He, Ne, Ar, Kr and Xe) are powerful geochemical tracers because they have distinctive isotopic compositions in the atmosphere, crust and mantle. This study illustrates how noble gases can be used to trace fluid origins in high-temperature metamorphic and mineralising environments; and at the same time provides new information on the composition of noble gases in deeper parts of the crust than have been sampled previously.We report data for H₂O and CO₂ fluid inclusions trapped at greenschist to amphibolite facies metamorphic conditions associated with three different styles of mineralisation and alteration in the Proterozoic Mt Isa Inlier of Australia. Sulphide fluid inclusions are dominated by crustal ⁴He. However, co-variations in fluid inclusion ²⁰Ne/²²Ne, ²¹Ne/²²Ne, ⁴⁰Ar/³⁶Ar and ¹³⁶Xe/¹³⁰Xe indicate noble gases were derived from three or more reservoirs. In most cases, the fluid inclusions elemental noble gas ratios (e.g. Ne/Xe) are close to the ranges expected in sedimentary and crystalline rocks. However, the elemental ratios have been modified in some of the samples providing evidence for independent pulses of CO_2, and interaction of CO₂ with high-salinity aqueous fluids.Compositional variation is attributed to mixing of: (i) magmatic fluids (or deeply sourced metamorphic fluids) characterised by basement-derived noble gases with ²⁰Ne/²²Ne ∼ 8.4, ²¹Ne/²²Ne ∼ 0.4, ⁴⁰Ar/³⁶Ar ∼ 40,000 and ¹³⁶Xe/¹³⁰Xe ∼ 8; (ii) basinal-metamorphic fluids with a narrow range of compositions including near-atmospheric values and (iii) noble gases derived from the meta-sedimentary host-rocks with ²⁰Ne/²²Ne ∼ 8-9.8, ²¹Ne/²²Ne < 0.1, ⁴⁰Ar/³⁶Ar < 2500 and ¹³⁶Xe/¹³⁰Xe ∼ 2.2.These data provide the strongest geochemical evidence available for the involvement of fluids from two distinct geochemical reservoirs in Mt Isa’s largest ore deposits. In addition the data show how noble gases in fluid inclusions can provide information on fluid origins, the composition of the crust’s major lithologies, fluid-rock interactions and fluid-fluid mixing or immiscibility processes.     

Noble gas composition of the solar wind as collected by the Genesis mission

We present the elemental and isotopic composition of noble gases in the bulk solar wind collected by the NASA Genesis sample return mission. He, Ne, and Ar were analyzed in diamond-like carbon on a silicon substrate (DOS) and ^84,86Kr and ^129,132Xe in silicon targets by UV laser ablation noble gas mass spectrometry. Solar wind noble gases are quantitatively retained in DOS and with exception of He also in Si as shown by a stepwise heating experiment on a flown DOS target and analyses on other bulk solar wind collector materials. Solar wind data presented here are absolutely calibrated and the error of the standard gas composition is included in stated uncertainties. The isotopic composition of the light noble gases in the bulk solar wind is as follows: ³He/⁴He: (4.64 ± 0.09) × 10⁻⁴, ²⁰Ne/²²Ne: 13.78 ± 0.03, ²¹Ne/²²Ne: 0.0329 ± 0.0001, ³⁶Ar/³⁸Ar 5.47 ± 0.01. The elemental composition is: ⁴He/²⁰Ne: 656 ± 5, and ²⁰Ne/³⁶Ar 42.1 ± 0.3. Genesis provided the first Kr and Xe data on the contemporary bulk solar wind. The preliminary isotope and elemental composition is: ⁸⁶Kr/⁸⁴Kr: 0.302 ± 0.003, ¹²⁹Xe/¹³²Xe: 1.05 ± 0.02, ³⁶Ar/⁸⁴Kr 2390 ± 150, and ⁸⁴Kr/¹³²Xe 9.5 ± 1.0. The ³He/⁴He and the ⁴He/²⁰Ne ratios in the Genesis DOS target are the highest solar wind values measured in exposed natural and artificial targets. The isotopic composition of the other noble gases and the Kr/Xe ratio obtained in this work agree with data from lunar samples containing “young” (∼100 Ma) solar wind, indicating that solar wind composition has not changed within at least the last 100 Ma. Genesis could provide in many cases more precise data on solar wind composition than any previous experiment. Because of the controlled exposure conditions, Genesis data are also less prone to unrecognized systematic errors than, e.g., lunar sample analyses. The solar wind is the most authentic sample of the solar composition of noble gases, however, the derivation of solar noble gas abundances and isotopic composition using solar wind data requires a better understanding of fractionation processes acting upon solar wind formation.     

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.     

Apparent decoupling of the He and Ne isotope systematics of the Icelandic mantle: The role of He depletion, melt mixing, degassing fractionation and air interaction

We present new He-Ne data for geothermal fluids and He-Ne-Ar data for basalts from throughout the Icelandic neovolcanic zones and older parts of the Icelandic crust. Geothermal fluids, subglacial glasses, and mafic phenocrysts are characterized by a wide range in helium isotope ratios (³He/⁴He) encompassing typical MORB-like ratios through values as high as 36.8 R_A (where R_A = air ³He/⁴He). Although neon in geothermal fluids is dominated by an atmospheric component, samples from the northwest peninsula show a small excess of nucleogenic ²¹Ne, likely produced in-situ and released to circulating fluids. In contrast, geothermal fluids from the neovolcanic zones show evidence of a contribution of mantle-derived neon, as indicated by ²⁰Ne enrichments up to 3% compared to air. The neon isotope composition of subglacial glasses reveals that mantle neon is derived from both depleted MORB-mantle and a primordial, ‘solar’ mantle component. However, binary mixing between these two endmembers can account for the He-Ne isotope characteristics of the basalts only if the ³He/²²Ne ratio of the primordial mantle endmember is lower than in the MORB component. Indeed, the helium to neon elemental ratios (⁴He/²¹Ne∗ and ³He/²²Ne_s where ²¹Ne∗ = nucleogenic ²¹Ne and ²²Ne_s = ‘solar’-derived ²²Ne) of the majority of Icelandic subglacial glasses are lower than theoretical values for Earth’s mantle, as observed previously for other OIB samples. Helium may be depleted relative to neon in high-³He/⁴He ratio parental melts due to either more compatible behavior during low-degree partial melting or more extensive diffusive loss relative to the heavier noble gases. However, Icelandic glasses show higher ⁴He/⁴⁰Ar∗ (⁴⁰Ar∗ = radiogenic Ar) values for a given ⁴He/²¹Ne∗ value compared to the majority of other OIB samples: this observation is consistent with extensive open-system equilibrium degassing, likely promoted by lower confining pressures during subglacial eruptions of Icelandic lavas. Taken together, the He-Ne-Ar systematics of Icelandic subglacial glasses are imprinted with the overlapping effects of helium depletion in the high-³He/⁴He ratio parental melt, binary mixing of two distinct mantle components, degassing fractionation and interaction with atmospheric noble gases. However, it is still possible to discern differences in the noble gas characteristics of the Icelandic mantle source beneath the neovolcanic zones, with MORB-like He-Ne isotope features prevalent in the Northern Rift Zone and a sharp transition to more primitive ‘solar-like’ characteristics in central and southern Iceland.     

Refertilization of lithospheric mantle beneath the Yangtze craton in south-east China: Evidence from noble gases geochemistry

The Yangtze craton (YC), in eastern China, is one of the oldest cratons in the world and is characterized by a complex tectonic and geodynamic evolution. This evolution regards most of the eastern China craton, which since Mesozoic time has undergone significant thinning (> 200 km) of Archean lithosphere. This thinning favored the refertilization of the old refractory subcontinental lithospheric mantle (SCLM) by the upwelling of younger fertile asthenosphere. Whether this feature is localized only beneath certain areas of eastern China or is a more widespread characteristic of the mantle, including the YC, is a matter of debate.In order to constrain the history of the YC SCLM, we have measured the He- and Ar-isotopic compositions of fluid inclusions hosted in mantle xenoliths in the Lianshan area, which is part of the poorly investigated YC in south-east China. We also report new mineral chemistry and trace element compositions of clinopyroxenes from the same suite of samples, for comparison with noble gases. Two distinct types of xenoliths can be identified: Type 1, characterized by mantle-like He-isotopic (³He/⁴He) ratios (up to 9.1 Ra), represents fragments of a fertile lithospheric mantle; Type 2, showing ³He/⁴He values in the SCLM range (³He/⁴He < 7 Ra), represents shallow relicts of a refractory mantle. The patterns of rare-earth elements as well as the Y and Yb concentrations in the clinopyroxenes normalized to primitive mantle (Y_N and Yb_N, respectively) indicate that fractional partial melting might have affected the local mantle by < 3% in Type 1 and up to 20% in Type 2 xenoliths from Lianshan, respectively. The range of ⁴He/⁴⁰Ar* (⁴⁰Ar* is corrected for atmospheric contamination) ranges from 4.9 × 10^− 4 to 3.6 × 10^− 1, which is below the typical production ratio of the mantle (⁴He/⁴⁰Ar* = 1–5); this range is however compatible with this fractional partial melting. The variable ³He/⁴He and ⁴He/⁴⁰Ar* values in Lianshan xenoliths suggest that the local mantle source was also influenced by kinetic fractionation, possibly triggered by metasomatic melts. Metasomatism associated with carbonatitic melts, together with fluxing by CO₂-rich fluids, have permeated the mantle beneath Lianshan, generating the observed decoupling between noble gases and trace elements. The interpretative framework is also applicable for other mantle xenoliths from eastern China, indicating that the refertilization of the SCLM by ascending mantle-like melts is common also to YC, which can be identified using noble gases.     

Vesiculation and vesicle loss in mid-ocean ridge basalt glasses: He, Ne, Ar elemental fractionation and pressure influence

Sarda and Graham (1990) proposed that in mid-ocean ridge basalts (MORBs), degassing occurs through equilibrium vesiculation followed by various extents of vesicle loss. This model predicts that in a bulk sample of MORB glass with vesicles, the rare gases represent a binary mixture between a vesicle component and a component dissolved in the melt. As vesiculation is expected to produce very different rare gas concentrations and elemental ratios in gas and melt, binary mixing systematics should be recorded in the MORB rare gas abundance data. Indeed, a large range of ⁴He/⁴⁰Ar∗ ratios was known to exist, but these binary mixing systematics remained elusive because helium was used as a proxy for rare gas abundance because helium is not affected by air addition. Here we show that using Ar instead of He, the ⁴He/⁴⁰Ar∗ ratio is higher where the Ar concentration is lower, as expected from simple binary mixing systematics.Taking advantage of the growing Ne database, we further show that the predicted binary mixing is recorded by the He-Ar and He-Ne couples, provided He concentration is not used to trace vesicle abundance. This is because a significant part of helium remains in the melt due to its higher solubility. In contrast, Ar or Ne concentrations, which can both be corrected for air addition, clearly trace vesicles and yield binary mixing patterns that hold for ridges worldwide. The model of vesiculation and vesicle loss thereby finds geochemical support in the rare gas abundance data.The He-Ne-Ar concentration data is best explained by assuming the ratio of helium to neon or argon solubility is about 5 to 15 times higher than values measured in 1 bar laboratory experiments, due to higher He and lower Ne and Ar solubilities. We propose that this is a pressure effect, and vesiculation mainly occurs during magma ascent in the mantle after melting.     

Geochemical constraints on the genesis of the Bayan Obo Fe-Nb-REE deposit in Inner Mongolia, China

Trace element and isotopic compositions of carbonate from ore bodies, country rock which hosts the ore bodies (H8 dolomite), a carbonatite dyke exposed in Dulahala near Bayan Obo, and rare earth element (REE)-rich dolomite in Bayan Obo have been determined to understand the genesis of the Bayan Obo Fe-Nb-REE ore deposit, the world’s largest resource of REE. The REE and trace element distribution patterns of samples from the REE-rich carbonatite dykes are identical to those of mineralized carbonate rocks, indicating a genetic linkage between the REE-rich carbonatite and mineralization in this region. By contrast, carbon and oxygen isotopes in the mineralized carbonate varied significantly, δ¹³C = −7.98‰ to −1.12‰, δ¹⁸O = 8.60-25.69‰, which are distinctively different from those in mantle-derived carbonatite. Abnormal isotopic fractionations between dolomite and calcite suggest that these two minerals are in disequilibrium in the carbonatite dyke, ore bodies, and H8 marble from Bayan Obo. This isotopic characteristic is also found in mineralized sedimentary marine micrite from Heinaobao, ∼25 km southeast of the Bayan Obo Fe-Nb-REE ore deposit. These facts imply that the carbonate minerals in the Bayan Obo deposit have resulted from sedimentary carbonate rocks being metasomatised by mantle-derived fluids, likely derived from a REE-enriched carbonatitic magma. The initial Nd isotope values of ore bodies and carbonatite dykes are identical, indicating that ore bodies, carbonatite dykes and veins may have a similar REE source.     

Trace element partitioning between carbonatitic melts and mantle transition zone minerals: Implications for the source of carbonatites

The occurrence of CO₂-rich lavas (carbonatites, kimberlites) and carbonate-rich xenoliths provide evidence for the existence of carbonatitic melts in the mantle. To model the chemical composition of such melts in the deep mantle, we experimentally determined partition coefficients for 23 trace elements (including REE, U-Th, HFSE, LILE) between deep mantle minerals and carbonatite liquids at 20 and 25 GPa and 1600 °C. Under these conditions, majoritic garnet and CaSiO₃ perovskite are the main reservoirs for trace elements. This study used both femtosecond LA-ICP-MS and SIMS techniques to measure reliable trace element concentrations. Comparison of the two techniques shows a general agreement, except for Sc and Ba. Our experimentally determined partition coefficients are consistent with the lattice strain model. The data suggest an effect of melt structure on partition coefficients in this pressure range. For instance, strain-free partition coefficient (D₀) for majorite-carbonatite melts do not follow the order of cation valence, , observed for majorite-CO₂-free silicate melts. The newly determined partition coefficients were combined with trace element composition of majoritic garnets found as inclusions in diamond to model trace element patterns of deep-seated carbonatites. The result compares favorably with natural carbonatites. This suggests that carbonatites can originate from the mantle transition zone.     

Noble gas study of HIMU and EM ocean island basalts in the Polynesian region

Noble gas isotopes of HIMU and EM ocean island basalts from the Cook-Austral and Society Islands were investigated to constrain their origins. Separated olivine and clinopyroxene (cpx) phenocrysts were used for noble gas analyses. Since samples are relatively old, obtained from the oceanic area and showing chemical zoning in cpx phenocrysts, several tests on sample preparation and gas extraction methods were performed. First, by comparing heating and crushing methods, it has been confirmed that the crushing method is suitable to obtain inherent magmatic noble gases without radiogenic and cosmogenic components which were yielded after eruption, especially for He and Ne analyses. Second, noble gas compositions in the core and the rim of cpx phenocrysts were measured to evaluate the zoning effect on noble gases. The result has been that noble gas concentrations and He and Ne isotope ratios are different between them. The enrichment of noble gases in the rim compared to the core is probably due to fractional crystallization. Difference of He and Ne isotope ratios is explained by cosmogenic effect, and isotope ratios of the trapped component seem to be similar between the rim and the core. Third, leaching test reveals no systematic differences in noble gas compositions between leached and unleached samples.³He/⁴He ratios of HIMU samples in the Cook-Austral Islands are uniform irrespective of phenocryst type (olivine and cpx) and age of samples (10–18 Ma), and lower (average 6.8 R_A) than those of the Pacific MORB. On the other hand, ³He/⁴He of EM samples in the Cook-Austral Islands are similar to MORB values. EM samples in the Society Islands show rather higher ³He/⁴He than MORB. Ne, Kr and Xe isotope ratios are almost atmospheric within analytical uncertainties. ⁴⁰Ar/³⁶Ar are not so high as those of MORB. Anomalous noble gas abundance pattern such as He and Ne depletion and Kr and Xe enrichment relative to atmospheric abundances was observed. Furthermore, Ne/Ar and Kr/Ar show correlation with some trace elemental ratios like La/Yb.Lower ³He/⁴He of HIMU than MORB values requires relatively high time-integrated (U + Th)/³He for the HIMU source, which suggests that the HIMU source was produced from recycled materials which had been once located near the Earth’s surface. Moreover, extreme noble gas abundance pattern and strong correlation of Ne/Ar and Kr/Ar with La/Yb indicate that the HIMU endmember is highly depleted in light noble gases and enriched in heavy noble gases. Such feature is not common to mantle materials and is rather similar to the noble gas abundance patterns of the old oceanic crust and sediment, which supports the model that the HIMU source originates from subducted oceanic crust and/or sediment.If the HIMU source corresponds to the oceanic crust which subducted at 1–2 Ga as suggested by Pb isotope studies, however, the characteristic ³He/⁴He of HIMU (6.8 R_A) would be too high because radiogenic ⁴He produced by U and Th decay should dramatically decrease ³He/⁴He. To overcome this problem, the He open system model is introduced which includes the effects of ⁴He production and diffusion between the HIMU source material and the surrounding mantle. This model favors that the HIMU source resides in the upper mantle, rather than in the lower mantle. Furthermore, this model predicts the thickness of the HIMU source to be in the order of 1 km.In contrast to low and uniform ³He/⁴He character of HIMU, ³He/⁴He of EM are rather variable. Entrainment of upper mantle material and/or a less-degassed component are required to explain the observed ³He/⁴He of EM in the Polynesian area. Participation of the less-degassed component would be related to the “superplume” below the Polynesian region.     

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.     

U-Th-Pb and Lu-Hf isotopic constraints on the evolution of sub-continental lithospheric mantle, French Massif Central

We have carried out a Pb double-spike and Lu-Hf isotope study of clinopyroxenes from spinel-facies mantle xenoliths entrained in Cenozoic intraplate continental volcanism of the French Massif Central (FMC). U-Th-Pb and Lu-Hf isotope systematics verify the existence of different lithospheric domains beneath the northern and southern FMC. Northern FMC clinopyroxenes have extreme Lu/Hf ratios and ultra-radiogenic Hf (ε_Hf = +39.6 to +2586) that reflect ∼15-25% partial melting in Variscan times (depleted mantle model ages ∼360 Ma). Zr, Hf and Th abundances in these clinopyroxenes are low and unaffected by hydrous/carbonatitic metasomatism that overprinted LILE and light REE abundances and caused decoupling of Lu/Hf-Sm/Nd ratios and Nd-Hf isotopes (ε_Nd = +2.1 to +91.2). Pb isotopes of northern FMC clinopyroxenes are radiogenic (²⁰⁶Pb/²⁰⁴Pb > 19), and typically more so than the host intraplate volcanic rocks. ²³⁸U/²⁰⁴Pb ratios range from 17 to 68, and most samples have distinctively low ²³²Th/²³⁸U (<1) and ²³²Th/²⁰⁴Pb (3-22). Clinopyroxenes from southern FMC lherzolites are generally marked by overall incompatible trace element enrichment including Zr, Hf and Th abundances, and have Pb isotopes that are similar to or less radiogenic than the host volcanic rocks. Hf isotope ratios are less radiogenic (ε_Hf = +5.4 to +41.5) than northern FMC mantle and have been overprinted by silicate-melt-dominated metasomatism that affected this part of FMC mantle. Major element and Lu concentrations of clinopyroxenes from southern FMC harzburgites are broadly similar to northern FMC clinopyroxenes and suggest they experienced similar degrees of melt extraction as northern FMC mantle. ²³⁸U/²⁰⁴Pb (53-111) and ²³²Th/²⁰⁴Pb ratios (157-355) of enriched clinopyroxenes from the southern FMC are extreme and significantly higher than the intraplate volcanic rocks. In summary, mantle peridotites from different parts of the FMC record depletion at ∼360 Ma during Variscan subduction, followed by differing styles of enrichment. Northern FMC mantle was overprinted by a fluid/carbonatitic metasomatic agent that carried elements like U, Pb, Sr and light REE. In contrast, much of the southern FMC mantle was metasomatised by a small-degree partial silicate melt resulting in enrichment of all incompatible trace elements. The extreme mantle ²³⁸U/²⁰⁴Pb (northern and southern FMC), ²³²Th/²³⁸U (northern FMC) and ²³²Th/²⁰⁴Pb ratios (southern FMC), coupled with unremarkable present-day Pb isotope ratios, constrain the timing of enrichment. Mantle metasomatism is a young feature related to melting of the upwelling mantle responsible for Cenozoic FMC volcanism, rather than subduction-related metasomatism intimately associated with mantle depletion during the Variscan orogeny. The varying metasomatic styles relate to pre-existing variations in the thickness of the continental lithospheric lid, which controlled the extent to which upwelling mantle could ascend and melt. In the northern FMC, a thicker and more refractory lithospheric lid (?80 km) only allowed incipient degrees of melting resulting in fluid/carbonatitic metasomatism of the overlying sub-continental lithospheric mantle. The thinner lithospheric lid of the southern FMC (?70 km) allowed larger degrees of melting and resulted in silicate-melt-dominated metasomatism, and also focused the location of the volcanic fields of the FMC above this region.     

Oxygen fugacity dependence of Os solubility in haplobasaltic melt

Os equilibrium solubilities were determined at 1350 °C over a wide range of oxygen fugacities (−12 < log fO₂ < −7) applying the mechanically assisted equilibration technique (MAE) at 10⁵ Pa (= 1 bar). Os concentrations in the glass samples were analysed using ID-NTIMS. Additional LA-ICP-MS and SEM analyses were performed to detect, visualize and analyse the nature and chemistry of “nanonuggets.” Os solubilities determined range at a constant temperature of 1350 °C from 0.63 ± 0.04 to 37.4 ± 1.16 ppb depending on oxygen fugacity. At the highest oxygen fugacities, Os³⁺ can be confirmed as the main oxidation state of Os. At low oxygen fugacities (below log fO₂ = −8), samples are contaminated by nanonuggets which, despite the MAE technique, were still not removed entirely from the melt. However, the present results indicate that applying MAE technology does reduce the amount of nanonuggets present significantly, resulting in the lowest Os solubility results reported to date under these experimental conditions, and extending the experimentally accessible range of fO₂ for these studies to lower values. Calculated metal/silicate melt partition coefficients are therefore higher compared to previous studies, making Os more siderophile. Neglecting the as yet unknown temperature dependence of the Os metal/silicate melt partition coefficient, extrapolation of the obtained Os solubilities to conditions for core-mantle equilibrium, results in a , while metallic alloy/silicate melt partition coefficients range from 1.4 × 10⁶ to 8.6 × 10⁷, in agreement with earlier findings. Therefore remains too high by 2-4 orders of magnitude to explain the Os abundance in the Earth’s mantle as result of core-mantle equilibrium during core formation.     

Isotopic fractionation of noble gases by diffusion in liquid water: Molecular dynamics simulations and hydrologic applications

Despite their great importance in low-temperature geochemistry, diffusion coefficients of noble gas isotopes in liquid water (D) have been measured only for the major isotopes of helium, neon, krypton and xenon. Data on the diffusion coefficients of minor noble gas isotopes are essentially non-existent and so typically have been estimated by a kinetic-theory model in which D varies as the inverse square root of the isotopic mass (m): D ∝ m^−0.5. To examine the validity of the kinetic-theory model, we performed molecular dynamics (MD) simulations of the diffusion of noble gases in ambient liquid water. Our simulation results agree with available experimental data on the solvation structure and diffusion coefficients of the major noble gas isotopes and reveal for the first time that the isotopic mass-dependence of all noble gas self-diffusion coefficients has the power-law form D ∝ m^−β with 0 < β < 0.2. Thus our results call into serious question the widespread assumption that the ‘square-root’ model can be applied to estimate the kinetic fractionation of noble gas isotopes caused by diffusion in ambient liquid water. To illustrate the importance of this finding, we used the diffusion coefficients determined in our MD simulations to reconsider the geochemical modeling of ²⁰Ne/²²Ne and ³⁶Ar/⁴⁰Ar isotopic ratios in three representative hydrologic studies. Our new modeling results indicate that kinetic isotopic fractionation by diffusion may play a significant role in noble gas transport processes in groundwater.     

Tracing nitrogen in volcanic and geothermal volatiles from the Nicaraguan volcanic front

We report new chemical and isotopic data from 26 volcanic and geothermal gases, vapor condensates, and thermal water samples, collected along the Nicaraguan volcanic front. The samples were analyzed for chemical abundances and stable isotope compositions, with a focus on nitrogen abundances and isotope ratios. These data are used to evaluate samples for volatile contributions from magma, air, air-saturated water, and the crust. Samples devoid of crustal contamination (based upon He isotope composition) but slightly contaminated by air or air-saturated water are corrected using N₂/Ar ratios in order to obtain primary magmatic values, composed of contributions from upper mantle and subducted hemipelagic sediment on the down-going plate. Using a mantle endmember with δ¹⁵N = −5‰ and N₂/He = 100 and a subducted sediment component with δ¹⁵N = +7‰ and N₂/He = 10,500, the average sediment contribution to Nicaraguan volcanic and geothermal gases was determined to be 71%. Most of the gases were dominated by sediment-derived nitrogen, but gas from Volcán Mombacho, the southernmost sampling location, had a mantle signature (46% from subducted sediment, or 54% from the mantle) and an affinity with mantle-dominated gases discharging from Costa Rica localities to the south. High CO₂/N_{2 exc.} ratios (N_{2 exc.} is the N₂ abundance corrected for contributions from air) in the south are similar to those in Costa Rica, and reflect the predominant mantle wedge input, whereas low ratios in the north indicate contribution by altered oceanic crust and/or preferential release of nitrogen over carbon from the subducting slab. Sediment-derived nitrogen fluxes at the Nicaraguan volcanic front, estimated by three methods, are 7.8 × 10⁸ mol N/a from ³He flux, 6.9 × 10⁸ mol/a from SO₂ flux, and 2.1 × 10⁸ and 1.3 × 10⁹ mol/a from CO₂ fluxes calculated from ³He and SO₂, respectively. These flux results are higher than previous estimates for Central America, reflecting the high sediment-derived volatile contribution and the high nitrogen content of geothermal and volcanic gases in Nicaragua. The fluxes are also similar to but higher than estimated hemipelagic nitrogen inputs at the trench, suggesting addition of N from altered oceanic basement is needed to satisfy these flux estimates. The similarity of the calculated input of N via the trench to our calculated outputs suggests that little or none of the subducted nitrogen is being recycled into the deeper mantle, and that it is, instead, returned to the surface via arc volcanism.     

Noble gas isotope sites and mobility in mafic rocks and olivine

Important He and Ar isotope studies on rocks and minerals, relevant to the geochemical and degassing history of the Earth, are often hampered by insufficient knowledge of the retentivity of different types of sites in minerals (inclusions, matrix) for these species, and of the relative importance of radiogenic and trapped components and possible differences in their behavior.To identify sites of noble gas isotopes, shed some light on their origin and estimate their residence times in olivine, which is a mineral considered as a good natural sampler, we investigated 2.5 Ga old ultramafic rocks from the Monche Pluton (Kola Peninsula, north-east part of the Baltic shield) using several extraction methods: crushing, fusion, slow step-wise and rapid incremental heating. Previous studies indicated that these rocks contain mainly trapped noble gases; however, to constrain the possible contribution of in-situ generated radiogenic helium, U and Th concentrations were also measured in the samples.The helium release pattern obtained by relatively fast (∼1.5 h long) incremental heating of olivine includes three distinct release peaks for helium: a low-temperature (600 °C) l-peak, a middle (800-1100 °C) m-peak and a high-temperature (∼1400 °C) h-peak. However, helium extraction from a powdered aliquot of the same olivine yields mainly the middle m-peak indicating that gases released in the l- and h-peaks occupy gas-liquid inclusions opened in the course of crushing and grinding. Moreover, slow step-wise heating (14 h) also results in a broad He release peak but in two well-separated l- and h-peaks of non-atmospheric ⁴⁰Ar∗. This feature implies helium migration from l- and h-vesicles into the matrix m during long step-wise heating experiments, whereas less movable Ar remains in inclusions at even relatively high almost-magmatic temperatures.Using a simple phenomenological model envisaging the three different residence sites for noble gases, both fast- and slow-heating release patterns for ⁴⁰Ar∗ and He, including those for the crushed sample, could be reproduced. The diffusion parameters inferred from the modeling of olivine (D₀ = 2.4 × 10⁻² cm² s⁻¹ and Ea = 133 kJ mol⁻¹) are similar to those published by Shuster et al. (2003) and Blard et al. (2008). The high matrix/fluid solubility coefficient for helium, H_He ∼ 0.01, exceeds estimates reported by Trull and Kurz (1993); however, the product D_He(T) × H_He, the “permeability” (that governs He migration in vesicles + matrix composed materials), is very similar to their value. Extrapolation to the ambient temperature (0 °C) gives long and similar helium residence times in l- and h-vesicles, exceeding 10¹⁰ yrs, and even longer time scales ∼10¹⁶ yrs are obtained for the helium residence in the matrix. Therefore, at low temperatures our samples may be considered as excellent samplers of trapped volatile species, including helium.     

The noble gas geochemistry of natural CO2 gas reservoirs from the Colorado Plateau and Rocky Mountain provinces, USA

Identification of the source of CO₂ in natural reservoirs and development of physical models to account for the migration and interaction of this CO₂ with the groundwater is essential for developing a quantitative understanding of the long term storage potential of CO₂ in the subsurface. We present the results of 57 noble gas determinations in CO₂ rich fields (>82%) from three natural reservoirs to the east of the Colorado Plateau uplift province, USA (Bravo Dome, NM., Sheep Mountain, CO. and McCallum Dome, CO.), and from two reservoirs from within the uplift area (St. John’s Dome, AZ., and McElmo Dome, CO.). We demonstrate that all fields have CO₂/³He ratios consistent with a dominantly magmatic source. The most recent volcanics in the province date from 8 to 10 ka and are associated with the Bravo Dome field. The oldest magmatic activity dates from 42 to 70 Ma and is associated with the McElmo Dome field, located in the tectonically stable centre of the Colorado Plateau: CO₂ can be stored within the subsurface on a millennia timescale.The manner and extent of contact of the CO₂ phase with the groundwater system is a critical parameter in using these systems as natural analogues for geological storage of anthropogenic CO₂. We show that coherent fractionation of groundwater ²⁰Ne/³⁶Ar with crustal radiogenic noble gases (⁴He, ²¹Ne, ⁴⁰Ar) is explained by a two stage re-dissolution model: Stage 1: Magmatic CO₂ injection into the groundwater system strips dissolved air-derived noble gases (ASW) and accumulated crustal/radiogenic noble gas by CO₂/water phase partitioning. The CO₂ containing the groundwater stripped gases provides the first reservoir fluid charge. Subsequent charges of CO₂ provide no more ASW or crustal noble gases, and serve only to dilute the original ASW and crustal noble gas rich CO₂. Reservoir scale preservation of concentration gradients in ASW-derived noble gases thus provide CO₂ filling direction. This is seen in the Bravo Dome and St. John’s Dome fields. Stage 2: The noble gases re-dissolve into any available gas stripped groundwater. This is modeled as a Rayleigh distillation process and enables us to quantify for each sample: (1) the volume of groundwater originally ‘stripped’ on reservoir filling; and (2) the volume of groundwater involved in subsequent interaction. The original water volume that is gas stripped varies from as low as 0.0005 cm³ groundwater/cm³ gas (STP) in one Bravo Dome sample, to 2.56 cm³ groundwater/cm³ gas (STP) in a St. John’s Dome sample. Subsequent gas/groundwater equilibration varies within all fields, each showing a similar range, from zero to ∼100 cm³ water/cm³ gas (at reservoir pressure and temperature).     

Helium isotopic fingerprints of the heavy land subsidence left in groundwater of the Saga Plain,near Beppu–Shimabara graben,Kyushu, Japan

The Saga Plain is near Beppu–Shimabara graben, a region of potential active volcanism. In the graben, mantle He, which has a high ³He/⁴He ratio of 1.1 × 10⁻⁵, escapes easily from the underlying subduction zone. In groundwater of the Saga Plain, except in the Shiroishi district, this ratio gradually increased as the dissolved He content increased, to a maximum of 5 × 10⁻⁶. In central Shiroishi, however, the ratio reached a minimum of 8.7 × 10⁻⁷ with increasing dissolved He content, suggesting that groundwater in central Shiroishi has selectively accumulated radiogenic He, which has a very low ratio of 1 × 10⁻⁸, rather than reflecting the regional He, which is rich in mantle He. This can be explained if groundwater in Shiroishi has become mixed with fossil pore water drawn from impermeable marine clay aquitard layers. The withdrawal of pore water has also caused severe land subsidence in central Shiroishi.