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.     

He diffusion systematics in minerals: Evidence from synthetic monazite and zircon structure phosphates

Helium diffusivity was measured in synthetic rare-earth-element orthophosphates with systematically varying properties to evaluate potential controls on He transport in minerals. In the zircon structure phosphates (in this study, the phosphates of Tb, Dy, Ho, Er, Tm, Yb, and Lu as well as synthetic xenotime, YPO₄), He diffusion is strongly anisotropic. Transport apparently proceeds preferentially through channels aligned with the c-axis. The activation energy for diffusion is almost the same (122 ± 6 kJ/mol) in all members of this family, but there is a monotonic decrease in D_o with atomic number from TbPO₄ (∼10⁵ cm²/s) to LuPO₄ (∼10 cm²/s). The c-parallel channels become increasingly constricted in the same sequence, likely accounting for the systematically decreasing diffusivity. The He closure temperature (r = 1 cm, dT/dt = 10 °C/Myr) increases with atomic number from 44 °C for TbPO₄ to 88 °C for LuPO₄. Diffusion of radiogenic helium from natural zircon and xenotime is much slower than these synthetic analogs predict, suggesting that coupled substitution of REE and P for Zr and Si and/or radiation damage profoundly modify the energetics of interstitial He diffusion. In particular, α-recoil may play a key role by damaging the continuity and integrity of the channels.Monazite structure phosphates (here La, Ce, Pr, Nd, Sm, and Gd phosphate) are far more He retentive than those of the zircon structure. Activation energies increase smoothly with atomic number from LaPO₄ (183 kJ/mol) to NdPO₄ (224 kJ/mol) then decrease again to GdPO₄ (198 kJ/mol). D_o values mimic this pattern, spanning a range from ∼10⁻¹ cm²/s (GdPO₄) to 10⁴ cm²/s (NdPO₄). Nevertheless, He closure temperatures increase monotonically with atomic number, from 300 °C in LaPO₄ to 410 °C in GdPO₄. No evidence was obtained bearing on diffusion anisotropy, but the monazite structure lacks through-going channels so it is not expected. Diffusion parameters for radiogenic helium in natural monazite are similar to those obtained on the synthetic analogs.Ionic porosity is not the primary control on He diffusion in the orthophosphates. Within a given structure and with limited elemental substitution, ionic porosity and He closure temperature are negatively correlated, as predicted. However, differences between crystal structures are far more important than ion packing density: at comparable ionic porosity the monazite structure phosphates have He closure temperatures ∼300 °C higher than the xenotime structure phosphates. Modifications to the structures by radiation damage likely play a similarly significant role in controlling He diffusion.     

Zircon (U-Th)/He thermochronometry: He diffusion and comparisons with Ar/Ar dating

(U-Th)/He chronometry of zircon has a wide range of potential applications including thermochronometry, provided the temperature sensitivity (e.g., closure temperature) of the system be accurately constrained. We have examined the characteristics of He loss from zircon in a series of step-heating diffusion experiments, and compared zircon (U-Th)/He ages with other thermochronometric constraints from plutonic rocks. Diffusion experiments on zircons with varying ages and U-Th contents yield Arrhenius relationships which, after about 5% He release, indicate E_a = 163-173 kJ/mol (39-41 kcal/mol), and D₀ = 0.09-1.5 cm²/s, with an average E_a of 169 ± 3.8 kJ/mol (40.4 ± 0.9 kcal/mol) and average D₀ of 0.46^+0.87_−0.30 cm²/s. The experiments also suggest a correspondence between diffusion domain size and grain size. For effective grain radius of 60 μm and cooling rate of 10°C/myr, the diffusion data yield closure temperatures, T_c, of 171-196°C, with an average of 183°C. The early stages of step heating experiments show complications in the form of decreasing apparent diffusivity with successive heating steps, but these are essentially absent in later stages, after about 5-10% He release. These effects are independent of radiation dosage and are also unlikely to be due to intracrystalline He zonation. Regardless of the physical origin, this non-Arrhenius behavior is similar to predictions based on degassing of multiple diffusion domains, with only a small proportion (<2-4%) of gas residing in domains with a lower diffusivity than the bulk zircon crystal. Thus the features of zircon responsible for these non-Arrhenius trends in the early stages of diffusion experiments would have a negligible effect on the bulk thermal sensitivity and closure temperature of a zircon crystal.We have also measured single-grain zircon (U-Th)/He ages and obtained ⁴⁰Ar/³⁹Ar ages for several minerals, including K-feldspar, for a suite of slowly cooled samples with other thermochronologic constraints. Zircon He ages from most samples have 1 σ reproducibilities of about 1-5%, and agree well with K-feldspar ⁴⁰Ar/³⁹Ar multidomain cooling models for sample-specific closure temperatures (170-189°C). One sample has a relatively poor reproducibility of ∼24%, however, and a mean that falls to older ages than predicted by the K-feldspar model. Microimaging shows that trace element zonation of a variety of styles is most pronounced in this sample, which probably leads to poor reproducibility via inaccurate α-ejection corrections. We present preliminary results of a new method for characterizing U-Th zonation in dated grains by laser-ablation, which significantly improves zircon He age accuracy.In summary, the zircon (U-Th)/He thermochronometer has a closure temperature of 170-190°C for typical plutonic cooling rates and crystal sizes, it is not significantly affected by radiation damage except in relatively rare cases of high radiation dosage with long-term low-temperature histories, and most ages agree well with constraints provided by K-spar ⁴⁰Ar/³⁹Ar cooling models. In some cases, intracrystalline U-Th zonation can result in inaccurate ages, but depth-profiling characterization of zonation in dated grains can significantly improve accuracy and precision of single-grain ages.     

Weathering geochronology by (U-Th)/He dating of goethite

Nine samples of supergene goethite (FeOOH) from Brazil and Australia were selected to test the suitability of this mineral for (U-Th)/He dating. Measured He ages ranged from 61 to 8 Ma and were reproducible to better than a few percent despite very large variations in [U] and [Th]. In all samples with internal stratigraphy or independent age constraints, the He ages corroborated the expected relationships. These data demonstrate that internally consistent He ages can be obtained on goethite, but do not prove quantitative ⁴He retention. To assess possible diffusive He loss, stepped-heating experiments were performed on two goethite samples that were subjected to proton irradiation to produce a homogeneous distribution of spallogenic ³He. The ³He release pattern indicates the presence of at least two diffusion domains, one with high helium retentivity and the other with very low retentivity at Earth surface conditions. The low retentivity domain, which accounts for ∼ 5% of ³He, contains no natural ⁴He and may represent poorly crystalline or intergranular material which has lost all radiogenic ⁴He by diffusion in nature. Diffusive loss of ³He from the high retentivity domain is independent of the macroscopic dimensions of the analyzed polycrystalline aggregate, so probably represents diffusion from individual micrometer-size goethite crystals. The ⁴He/³He evolution during the incremental heating experiments shows that the high retentivity domain has retained 90%-95% of its radiogenic helium. This degree of retentivity is in excellent agreement with that independently predicted from the helium diffusion coefficients extrapolated to Earth surface temperature and held for the appropriate duration. Considering both the high and low retentivity domains, these data indicate that one of the samples retained 90% of its radiogenic ⁴He over 47.5 Ma and the other retained 86% over 12.3 Ma. Thus while diffusive-loss corrections to supergene goethite He ages are required, these initial results indicate that the corrections are not extremely large and can be rigorously quantified using the proton-irradiation ⁴He/³He method.     

Laser depth profiling studies of helium diffusion in Durango fluorapatite

Ultraviolet lasers coupled with sensitive mass spectrometers provide a useful way to measure laboratory-induced noble gas diffusion profiles in minerals, thus enabling the calculation of diffusion parameters. We illustrate this laser ablation depth profiling (LADP) technique for a previously well-studied mineral-isotopic system: ⁴He in Durango fluorapatite. LADP studies were conducted on oriented, polished slabs from a single crystal that were heated under vacuum to a variety of temperatures between 300 and 450 °C for variable times. The resolved ⁴He profiles exhibited error-function loss as predicted by previous bulk ⁴He diffusion studies. All of the slabs, regardless of crystallographic orientation, yielded modeled diffusivities that are statistically co-linear on an Arrhenius diagram, suggesting no diffusional anisotropy of ⁴He in this material. The data indicate an activation energy of 142.2 ± 5.0 (2σ) kJ/mol and diffusivity at infinite temperature - reported as ln(D₀) - of −4.71 ± 0.94 (2σ) m²/s. These values imply a bulk closure temperature for ⁴He in Durango fluorapatite of 74 °C for a 50 μm radius grain, infinite cylinder geometry, and a cooling rate of 10 °C/Myr.     

Identification of He sources and estimation of He ages in groundwater of the North China Plain

Dissolved helium concentrations and ³He/⁴He ratios were measured for 18 groundwater samples collected from the Quaternary confined aquifers in the North China Plain (NCP). The dissolved helium concentrations ranged from 1 × 10⁻⁷ to 1 × 10⁻⁶ cm³STP·g⁻¹ in the 14 samples from the central plain, but was approximately two orders of magnitude higher, between 6 × 10⁻⁶ and 9 × 10⁻⁵ cm³STP·g⁻¹, in 4 samples from the coastal plain. Based on these concentrations and the corresponding ³He/⁴He ratios varying from 0.09 to 0.55 R_a (where R_a is the ³He/⁴He ratio of air), the dissolved helium in groundwater in the central plain was identified to be primarily a mixture of atmospheric helium with radiogenic helium and a representative radiogenic helium ratio was estimated to be 0.035 R_a. Despite the high fraction of terrigenic ⁴He in the samples from the coastal plain, their ³He/⁴He ratios were found to be significantly above this radiogenic value, ranging between 0.20 and 0.37 R_a, indicating the presence of a mantle-derived He component in this area. About 2–4% mantle helium was estimated to be present in the groundwater of the coastal plain, which probably is associated with the regional Cangdong fault and tectonic activities. Based on the radiogenic He component, ⁴He ages of the groundwater in the central plain were calculated by assuming either pure in situ production or an external helium flux J₀ of 4.7 × 10⁻⁸ cm³STPcm⁻²a⁻¹. The estimated ⁴He ages fall between 9.5 and 51.4 ka and are comparable to the ¹⁴C ages, suggesting that the results of ⁴He dating are reasonable and can be an effective tool to estimate groundwater residence times under suitable conditions.     

4He diffusion in specular hematite

In order to test the chronometer qualities of speculante for the (U + Th)/He dating method, ⁴He release experiments by stepwise heating of two specularites from the Rimbach mineralization locality in the southern Vosgues (France) have been carried out. The diffusion coefficients define linear Arrhenius plots within a temperature interval of 250 to 830 °C, which is suggestive of volume diffusion. Extrapolation of the diffusion behavior to 20° C yields diffusion coefficients (D₂₀ values) smaller than 10^?26 [cm² s^?1] for both hematites with activation energies at 116 [kJ/mole]. The results of our study suggest that specularite is a very helium retentive hematite variety which is capable of quantitatively retaining radiogenic helium over geologic periods of time.     

Rare gas isotopes and parent trace elements in ultrabasic-alkaline-carbonatite complexes, Kola Peninsula: identification of lower mantle plume component

During the Devonian magmatism (370 Ma ago) ∼20 ultrabasic-alkaline-carbonatite complexes (UACC) were formed in the Kola Peninsula (north-east of the Baltic Shield). In order to understand mantle and crust sources and processes having set these complexes, rare gases were studied in ∼300 rocks and mineral separates from 9 UACC, and concentrations of parent Li, K, U, and Th were measured in ∼70 samples. ⁴He/³He ratios in He released by fusion vary from pure radiogenic values ∼10⁸ down to 6 × 10⁴. The cosmogenic and extraterrestrial sources as well as the radiogenic production are unable to account for the extremely high abundances of ³He, up to 4 × 10⁻⁹ cc/g, indicating a mantle-derived fluid in the Kola rocks. In some samples helium extracted by crushing shows quite low ⁴He/³He = 3 × 10⁴, well below the mean ratio in mid ocean ridge basalts (MORB), (8.9 ± 1.0) × 10⁴, indicating the contribution of ³He-rich plume component. Magnetites are principal carriers of this component. Trapped ³He is extracted from these minerals at high temperatures 1100°C to 1600°C which may correspond to decrepitation or annealing primary fluid inclusions, whereas radiogenic ⁴He is manly released at a temperature range of 500°C to 1200°C, probably corresponding to activation of ⁴He sites degraded by U, Th decay.Similar ⁴He/³He ratios were observed in Oligocene flood basalts from the Ethiopian plume. According to a paleo-plate-tectonic reconstruction, 450 Ma ago the Baltica (including the Kola Peninsula) continent drifted not far from the present-day site of that plume. It appears that both magmatic provinces could relate to one and the same deep-seated mantle source.The neon isotopic compositions confirm the occurrence of a plume component since, within a conventional ²⁰Ne/²²Ne versus ²¹Ne/²²Ne diagram, the regression line for Kola samples is indistinguishable from those typical of plumes, such as Loihi (Hawaii). ²⁰Ne/²²Ne ratios (up to 12.1) correlate well with ⁴⁰Ar/³⁶Ar ones, allowing to infer a source ⁴⁰Ar/³⁶Ar ratio of about 4000 for the mantle end-member, which is 10 times lower than that of the MORB source end-member. In (³He/²²Ne)_PRIM versus (⁴He/²¹Ne)_RAD plot the Kola samples are within array established for plume and MORB samples; almost constant production ratio of (⁴He/²¹Ne)_RAD ≅ 2 × 10⁷ is translated via this array into (³He/²²Ne)_PRIM ∼ 10. The latter value approaches the solar ratio implying the non-fractionated solar-like rare gas pattern in a plume source.The Kola UACC show systematic variations in the respective contributions of in situ-produced radiogenic isotopes and mantle-derived isotopes. Since these complexes were essentially plutonic, we propose that the depth of emplacement exerted a primary control on the retention of both trapped and radiogenic species, which is consistent with geological observations. The available data allow to infer the following sequence of processes for the emplacement and evolution of Kola Devonian UACC: 1) Ascent of the plume from the lower mantle to the subcontinental lithosphere; the plume triggered mantle metasomatism not later than ∼700 to 400 Ma ago. 2) Metasomatism of the lithosphere (beneath the central part of the Kola Peninsula), including enrichment in volatile (e.g., He, Ne) and in incompatible (e.g., U, Th) elements. 3) Multistage intrusions of parental melts, their degassing, and crystallisation differentiation ∼370 Ma ago. 4) Postcrystallisation migration of fluids, including loss of radiogenic and of trapped helium. Based on model compositions of the principle terrestrial reservoirs we estimate the contributions (by mass) of the plume material, the upper mantle material, and the atmosphere (air-saturated groundwater), into the source of parent melt at ∼2%, 97.95%, and ∼0.05%, respectively.     

Irradiation records in regolith materials, II: Solar wind and solar energetic particle components in helium, neon, and argon extracted from single lunar mineral grains and from the Kapoeta howardite by stepwise pulse heating

High-resolution stepped heating has been used to extract light noble gases implanted in a suite of 13 individual lunar ilmenite and iron grains and in the Kapoeta howardite by solar wind (SW) and solar energetic particle (SEP) irradiation. Isotopic analyses of gases evolved at low temperatures from the lunar grains confirm the neon and argon compositions obtained by Pepin et al. (Pepin R. O., Becker R. H., and Schlutter D. J., “Irradiation records in regolith materials, I: Isotopic compositions of solar-wind neon and argon in single lunar regolith grains”, Geochim. Cosmochim. Acta63, 2145-2162, 1999) in an initial study of 11 regolith grains, primarily ilmenites. Combination of the data sets from both investigations yields ²⁰Ne/²²Ne = 13.85 ± 0.04, ²¹Ne/²²Ne = 0.0334 ± 0.0003, and ³⁶Ar/³⁸Ar = 5.80 ± 0.06 for the lunar samples; the corresponding ³⁶Ar/³⁸Ar ratio in Kapoeta is 5.74 ± 0.06. The neon ratios agree well with those measured by Benkert et al. (Benkert J.-P., Baur H., Signer P., and Wieler R., “He, Ne, and Ar from the solar wind and solar energetic particles in lunar ilmenites and pyroxenes”, J. Geophys. Res. (Planets)98, 13147-13162, 1993) in gases extracted from bulk lunar ilmenite samples by stepped acid etching and attributed by them to the SW. The ³⁶Ar/³⁸Ar ratios, however, are significantly above both Benkert et al.’s (1993) proposed SW value of 5.48 ± 0.05 and a later estimate of 5.58 ± 0.03 from an acid-etch analysis of Kapoeta (Becker R. H., Schlutter D. J., Rider P. E., and Pepin R. O., “An acid-etch study of the Kapoeta achondrite: Implications for the argon-36/argon-38 ratio in the solar wind”, Meteorit. Planet. Sci.33, 109-113, 1998). We believe, for reasons discussed here and in our earlier report, that 5.80 ± 0.06 ratio most nearly represents the wind composition. The ³He/⁴He ratio in low-temperature gas releases, not measured in the first particle suite, is found in several grains to be indistinguishable from Benkert et al.’s (1993) SW estimate. Elemental ratios of He, Ne, and Ar initially released from grain-surface SW implantation zones are solar-like, as found earlier by Pepin et al. (1999). Gases evolved from these reservoirs at higher temperatures show evidence for perturbations from solar elemental compositions by prior He loss, thermal mobilization of excess Ne from fractionated SW components, or both.Attention in this second investigation was focused on estimating the isotopic compositions of both the SW and the more deeply sited SEP components in regolith grains. Several high-temperature “isotopic plateaus”—approximately constant isotopic ratios in gas fractions released over a number of consecutive heating steps—were observed in the close vicinities of the SEP ratios for He, Ne, and Ar reported by Benkert et al. (1993). Arguments presented in the text suggest that these plateaus are relatively free of interferences from multicomponent mixing artifacts that can mimic pure component signatures. Average SEP compositions derived from the stepped-heating plateau measurements are in remarkable agreement with the Zürich acid-etch values for all three gases.     

A New Single Vacuum Furnace Design for Cosmogenic 3He Dating

Helium‐3 is a stable cosmogenic isotope that can be used to determine the time interval during which a rock sample has been at or close to the Earth’s surface. As a result of the high production rate of ‘cosmogenic’³He (≈ 130 at g^?1 year^?1) and the low detection limit of modern mass spectrometers, it is possible to date exceptionally young surfaces (≈ 1000 years). The precision and accuracy of cosmogenic ³He measurements depend critically on the passive helium blank (produced by the metalwork of the extraction furnace) which can be significant relative to the sample signals. We have developed and constructed, at the CRPG (Nancy, France), a new high temperature furnace (< 1500 °C) to extract helium in minerals such as apatite, pyroxene and olivine at 1050, 1350 and 1450 °C, respectively. The furnace demonstrated an excellent helium extraction yield (> 99% for olivine and pyroxene for heating times of 20–30 min and temperatures in the range 1050–1450 °C) and low residual helium contributions (the blank, obtained under the same analytical conditions as for sample extraction: 1 × 10^?15 mole ⁴He and < 4 × 10^?21 mole ³He). This is approximately an order of magnitude lower than those reported by other laboratories using conventional furnaces.     

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.     

The influence of artificial radiation damage and thermal annealing on helium diffusion kinetics in apatite

Recent work [Shuster D. L., Flowers R. M. and Farley K. A. (2006) The influence of natural radiation damage on helium diffusion kinetics in apatite. Earth Planet. Sci. Lett.249(3-4), 148-161] revealing a correlation between radiogenic ⁴He concentration and He diffusivity in natural apatites suggests that helium migration is retarded by radiation-induced damage to the crystal structure. If so, the He diffusion kinetics of an apatite is an evolving function of time and the effective uranium concentration in a cooling sample, a fact which must be considered when interpreting apatite (U-Th)/He ages. Here we report the results of experiments designed to investigate and quantify this phenomenon by determining He diffusivities in apatites after systematically adding or removing radiation damage.Radiation damage was added to a suite of synthetic and natural apatites by exposure to between 1 and 100 h of neutron irradiation in a nuclear reactor. The samples were then irradiated with a 220 MeV proton beam and the resulting spallogenic ³He used as a diffusant in step-heating diffusion experiments. In every sample, irradiation increased the activation energy (E_a) and the frequency factor (D_o/a²) of diffusion and yielded a higher He closure temperature (T_c) than the starting material. For example, 100 h in the reactor caused the He closure temperature to increase by as much as 36 °C. For a given neutron fluence the magnitude of increase in closure temperature scales negatively with the initial closure temperature. This is consistent with a logarithmic response in which the neutron damage is additive to the initial damage present. In detail, the irradiations introduce correlated increases in E_a and ln(D_o/a²) that lie on the same array as found in natural apatites. This strongly suggests that neutron-induced damage mimics the damage produced by U and Th decay in natural apatites.To investigate the potential consequences of annealing of radiation damage, samples of Durango apatite were heated in vacuum to temperatures up to 550 °C for between 1 and 350 h. After this treatment the samples were step-heated using the remaining natural ⁴He as the diffusant. At temperatures above 290 °C a systematic change in T_c was observed, with values becoming lower with increasing temperature and time. For example, reduction of T_c from the starting value of 71 to ∼52 °C occurred in 1 h at 375 °C or 10 h at 330 °C. The observed variations in T_c are strongly correlated with the fission track length reduction predicted from the initial holding time and temperature. Furthermore, like the neutron irradiated apatites, these samples plot on the same E_a − ln(D_o/a²) array as natural samples, suggesting that damage annealing is simply undoing the consequences of damage accumulation in terms of He diffusivity.Taken together these data provide unequivocal evidence that at these levels, radiation damage acts to retard He diffusion in apatite, and that thermal annealing reverses the process. The data provide support for the previously described radiation damage trapping kinetic model of Shuster et al. (2006) and can be used to define a model which fully accommodates damage production and annealing.     

Diffusion of lead in plagioclase and K-feldspar: an investigation using Rutherford Backscattering and Resonant Nuclear Reaction Analysis

Chemical diffusion of Pb has been measured in K-feldspar (Or₉₃) and plagioclase of 4 compositions ranging from An₂₃ to An₉₃ under anhydrous, 0.101 MPa conditions. The source of diffusant for the experiments was a mixture of PbS powder and ground feldspar of the same composition as the sample. Rutherford Backscattering (RBS) was used to measure Pb diffusion profiles. Over the temperature range 700–1050°C, the following Arrhenius relations were obtained (diffusivities in m²s^-1):Oligoclase (An₂₃): Diffusion normal to (001): log D = ( – 6.84 ± 0.59) – [(261 ± 13 kJ mol ^–1)/2.303RT]Diffusion normal to (010): log D = ( – 3.40 ± 0.50) – [(335 ± 11 kJ mol ^–1)/2.303RT]Andesine (An₄₃): Diffusion normal to (001): log D = ( – 6.73 ± 0.54) – [(266 ± 12 kJ mol ^–1)/2.303RT] Diffusion normal to (010) appears to be only slightly slower than diffusion normal to (001) in andesine.Labradorite (An₆₇): Diffusion normal to (001): log D = ( – 7.16 ± 0.61) – [(267 ± 13 kJ mol ^–1)/2.303RT] Diffusion normal to (010) is slower by 0.7 log units on average.Anorthite Diffusion normal to (010): log D = ( – 5.43 ± 0.48) – [(327 ± 11 kJ mol ^–1)/2.303RT]K-feldspar (Or₉₃): Diffusion normal to (001): log D = ( – 4.74 ± 0.52) – [(309 ± 16 kJ mol ^–1)/2.303RT] Diffusion normal to (010): log D = ( – 5.99 ± 0.51) – [(302 ± 11 kJ mol ^–1)/2.303RT]In calcic plagioclase, Pb uptake is correlated with a reduction of Ca, indicating the involvement of PbCa exchange in chemical diffusion. Decreases of Na and K concentrations in sodic plagioclase and K-feldspar, respectively, are correlated with Pb uptake and increase in Al concentration (measured by resonant nuclear reaction analysis), suggesting a coupled process for Pb exchange in these feldspars. These results have important implications for common Pb corrections and Pb isotope systematics. Pb diffusion in apatite is faster than in the investigated feldspar compositions, and Pb diffusion rates in titanite are comparable to both K-feldspar and labradorite. Given these diffusion data and typical effective diffusion radii for feldspars and accessory minerals, we may conclude that feldspars used in common Pb corrections are in general less inclined to experience diffusion-controlled Pb isotope exchange than minerals used in U-Pb dating that require a common Pb correction.     

Activation energy for diffusion of helium in water-saturated, compacted Na-montmorillonite

It is important to clarify the migration behavior of hydrogen gas dissolved in water-saturated, compacted bentonite, which is a promising material for geologic disposal of high-level waste and TRU waste disposal. The diffusion coefficients of helium, which can be detected under extremely low background conditions, in water-saturated, compacted Na-montmorillonite were determined as a function of temperature by a transient diffusion method. The activation energies for diffusion of helium were then obtained. The activation energies were from 6.9 ± 4.8 to 19 ± 2.8 kJ mol^− 1 and were regarded to be independent of dry density. The activation energies of helium in water-saturated Na-montmorillonite were roughly equal to those in bulk water, 14.9 kJ mol^− 1, and in ice, from 11 to 13 kJ mol^− 1. It is possible that helium diffuses not only in pore water but also in interlayer water.     

Geothermochronology based on noble gases: II. Stability of the (U-Th)/He isotope system in zircon

The kinetics of He migration from zircon of variable degree of metamictization was investigated. The migration parameters of He were experimentally determined, the influence of radiation damage and the degree of metamictization on the stability of the (U-Th)/He isotope system was evaluated, the mechanisms of noble gas escape from zircon were investigated, new data on the kinetics of He migration were obtained and compared with previous results for the kinetics of Xe migration from zircon of the same geologic objects. It was shown that He occurs in two energy positions in the zircon lattice: the main position (more than 80% He) with an activation energy of ∼39 kcal/mol and k ₀ = 10¹¹ yr⁻¹ and the second position with an activation energy for migration of 5–10 kcal/mol and k ₀ ∼ 10⁶ yr⁻¹. It was concluded that He migration from the main energy position is better described by a single-jump mechanism. The migration of He from the second energy position is consistent with the diffusion mechanism. It was shown that deviations from the linear dependence in the lnln(He₀/He_t)-1/T coordinates are probably related to the destruction of volume defects containing He atoms at high temperatures (more than 1000°C on the experimental time scale) resulting in a change from the single-jump to diffusion mechanism and the presence of atoms migrating via the diffusion mechanism. It was shown that the peak width in the spectrum of radiogenic He release and the appearance of a second peak also depend on the fraction of atoms migrating in accordance with the diffusion mechanism. It was found that the low activation energy for He release from the second energy position indicates the existence of continuous He loss from the zircon lattice.     

Mantle heat drives hydrothermal fluids responsible for carbonate-hosted base metal deposits: evidence from 3He/4He of ore fluids in the Irish Pb-Zn ore district

There is little consensus on whether carbonate-hosted base metal deposits, such as the world-class Irish Zn + Pb ore field, formed in collisional or extensional tectonic settings. Helium isotopes have been analysed in ore fluids trapped in sulphides samples from the major base metal deposits of the Irish Zn-Pb ore field in order to quantify the involvement of mantle-derived volatiles that require melting to be realised, as well as test prevailing models for the genesis of the ore fields. ³He/⁴He ratios range up to 0.2 R _a, indicating that a small but clear mantle helium contribution is present in the mineralising fluids trapped in galena and marcasite. Sulphides from ore deposits with the highest fluid inclusion temperatures (~200 °C) also have the highest ³He/⁴He (>0.15 R _a). Similar ³He/⁴He are recorded in fluids from modern continental regions that are undergoing active extension. By analogy, we consider that the hydrothermal fluids responsible for the carbonate-hosted Irish base metal mineralization circulated in thinned continental crust undergoing extension and demonstrate that enhanced mantle heat flow is ultimately responsible for driving fluid convection.     

Argon diffusion in plagioclase and implications for thermochronometry: A case study from the Bushveld Complex, South Africa

Plagioclase is not only the most abundant mineral in the Earth’s crust, but is present in almost all terrestrial tectonic settings and is widespread in most extraterrestrial material. Applying the K-Ar system to this common mineral would provide a powerful tool for quantifying thermal histories in a wide variety of settings. Nonetheless, plagioclase has rarely been used for thermochronometry, largely due to difficulties in simultaneously acquiring precise geochronologic data and quantifying argon diffusion kinetics from a mineral with low-K concentration. Here we describe an analytical technique that generates high-precision ⁴⁰Ar/³⁹Ar data and quantifies Ar diffusion kinetics of low-K minerals. We present results of five diffusion experiments conducted on single crystals of plagioclase from the Bushveld Complex, South Africa. The observed diffusion kinetics yield internally consistent thermochronological constraints, indicating that plagioclase is a reliable thermochronometer. Individual grains have activation energies of 155-178 kJ/mol and ln(D₀/a²) varies between 3.5 and 6.5. These diffusion parameters correspond to closure temperatures of 225-300 °C, for a 10 °C/Ma cooling rate. Age spectra generally conform to single-domain diffusive loss profiles, suggesting that grain-scale diffusion dominates argon transport in this fairly simple plagioclase. Conjointly examining several single-grain analyses enables us to distinguish episodic reheating from slow cooling and indicates that the Bushveld Complex cooled rapidly and monotonically from magmatic temperature to <300 °C over 3 Ma, followed by protracted cooling to ambient crustal temperatures of 150-200 °C over ∼600 Ma.     

New insights into the carrier phase(s) of extraterrestrial He in geologically old sediments

To better understand the composition, characteristics of helium diffusion, and size distribution of interplanetary dust particles (IDPs) responsible for the long-term retention of extraterrestrial ³He, we carried out leaching, stepped heating, and sieving experiments on pelagic clays that varied in age from 0.5 Ma to ∼90 Myr. The leaching experiments suggest that the host phase(s) of ³He in geologically old sediments are neither organic matter nor refractory phases, such as diamond, graphite, Al₂O₃, and SiC, but are consistent with extraterrestrial silicates, Fe-Ni sulfides, and possibly magnetite. Stepped heating experiments demonstrate that the ³He release profiles from the magnetic and non-magnetic components of the pelagic clays are remarkably similar. Because helium diffusion is likely to be controlled by mineral chemistry and structure, the stepped heating results suggest a single carrier that may be magnetite, or more probably a phase associated with magnetite. Furthermore, the stepped outgassing experiments indicate that about 20% of the ³He will be lost through diffusion at seafloor temperatures after 50 Myrs, while sedimentary rocks exposed on the Earth’s surface for the same amount of time would lose up to 60%. The absolute magnitude of the ³He loss is, however, likely to depend upon the ³He concentration profile within the IDPs, which is not well known. Contrary to previous suggestions that micrometeorites in the size range of 50-100 μm in diameter are responsible for the extraterrestrial ³He in geologically old sediments [Stuart, F.M., Harrop, P.J., Knott, S., Turner, G., 1999. Laser extraction of helium isotopes from Antarctic micrometeorites: source of He and implications for the flux of extraterrestrial ³He flux to earth. Geochim. Cosmochim. Acta63, 2653-2665], our sieving experiment demonstrates that at most 20% of the ³He is carried by particles greater than 50 μm in diameter. The size-distribution of the ³He-bearing particles implies that extraterrestrial ³He in sediments record the IDP flux rather than the micrometeorite flux.     

Tracer diffusion of Mg, Ca, Sr, and Ba in Na-aluminosilicate melts

We employed the thin source technique to investigate tracer diffusion of Mg, Ca, Sr, and Ba in glasses and supercooled melts of albite (NaAlSi₃O₈) and jadeite (NaAlSi₂O₆) compositions. The experiments were conducted at 1 bar and at temperatures between 645 and 1025°C. Typical run durations ranged between 30 min and 35 days. The analysis of the diffusion profiles was performed with the electron microprobe. Diffusivities of Ca, Sr, and Ba were found to be independent of either duration t of the experiment or tracer concentration M, initially introduced into the sample. Mg exhibits a diffusivity depending on run time and concentration and tracer diffusivity is derived by extrapolation to M/√t = 0. Temperature dependence of the diffusivity D can be represented by an Arrhenius equation D = D_o exp(−E_a/RT), yielding the following least-squares fit parameters (with D in m²/s and E_a in kJ/mol): D_Mg = 1.8 · 10⁻⁵ exp(−234 ± 20/RT), D_Ca = 3.5 · 10⁻⁶ exp(−159 ± 6/RT), D_Sr = 3.6 · 10⁻⁶ exp(−160 ± 6/RT), and D_Ba = 6.0 · 10⁻⁶ exp(−188 ± 12/RT) for albite; and D_Mg = 8.3 · 10⁻⁶ exp(−207 ± 18/RT), D_Ca = 3.8 · 10⁻⁶ exp(−153 ± 4/RT), D_Sr = 2.3 · 10⁻⁶ exp(−150 ± 4/RT), and D_Ba = 3.7 · 10⁻⁵ exp(−198 ± 4/RT) for jadeite composition. Ca and Sr diffusivities agree within error in both compositions and exhibit the fastest diffusivities, whereas Mg reveals the lowest diffusivity. The relationship between activation energy and radius shows a minimum at Ca and Sr for albite and jadeite compositions extending the relationship already observed elsewhere for alkalies. With increasing substitution of Si by (Na + Al), diffusivities increase, whereas activation energies decrease. Furthermore, a simple model modified from that of Anderson and Stuart (Anderson O. L. and Stuart D. A., “Calculation of activation energy of ionic conductivity in silica glasses by classical methods,” J. Am. Ceram. Soc.37, 573-580, 1954) is discussed for calculating the activation energies.     

He,Ne, Ar,and C isotope systematics of geothermal emanations in the Lesser Antilles Islands Arc

We present He, Ne, Ar, and C isotope analyses of hydrothermal brines and gases from fumaroles, hot springs, mofettes and hydrothermal exploration drillings on the major islands of the Lesser Antilles Arc. The origin of hydrothermal brines, which have been analyzed also for O and H isotopes, is essentially meteoric-hydrothermal. Air-corrected isotope compositions of helium (2.2 Rc/Ra < ³He/⁴He < 8.6 Rc/Ra) and carbon (−20 < δ¹³C_PDB < +0.5) are variable and require a variety of crustal and magmatic sources. The diversity of δ¹³C_PDB and ³He/CO₂ ratios within individual volcanic centres suggests that crustal sources (e.g., limestone) contaminate magmatic CO₂ en route from high-level magma reservoirs (depth < 15 km) to the surface. A similar contamination may be found for magmatic helium on distal springs. The ³He/⁴He signature of summit fumaroles, thought to reflect the ³He/⁴He signature of high-level magmas, shows a remarkable systematic variation along the arc. In addition, there is a correlation throughout the arc between published Sr, Pb, and Nd isotope signatures of lavas and the ³He/⁴He signatures of summit fumaroles. On the northern islands (Nevis, Montserrat, Guadeloupe, and Dominica) summit fumaroles have the N-MORB signature (³He/⁴He = 8 ± 1 R/Ra), and the isotope signature of lavas is not dissimilar from comparable intra-oceanic arc tholeiites elsewhere. Variable enrichments in radiogenic Sr and Pb have been reported for lavas of individual volcanic centres of the Southern Islands (Martinique, St.Lucia, and Grenada), and summit fumaroles on these centres match these variations by variable radiogenic He-enrichments, i.e., lower ³He/⁴He ratios. This correlation suggests that radiogenic Sr and Pb enrichments of lavas and low ³He/⁴He signatures on summit fumaroles have a common origin, i.e., a terrigenous contaminant derived from the Orinoco depositionary fan. Crustal assimilation is thought to decouple the He isotope system from any other radiogenic isotope system and, therefore, we argue that the observed correlation of He, Sr, Pb, and Nd isotope systems is related to a terrigenous contaminant derived from subducted sediments. Support for this scenario also comes from the matching of low ³He/⁴He ratios and tectonic features of the forearc thought to favor the subduction of forearc sediments.The present study offers a first clue that, under suitable conditions, crustal helium from oceanic sediments might be subducted to the depth of arc magma sources and, possibly, even recycled into the deeper mantle.