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.     

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.     

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.     

An absolute paleotemperature record from 10 to 6 Ka inferred from fluid inclusion D/H ratios of a stalagmite from Vancouver Island, British Columbia, Canada

By using continuous helium flow during the crushing of calcite speleothem samples, we are able to recover liberated inclusion waters without isotopic fractionation. A paleotemperature record for the Jacklah Jill Cave locality, Vancouver Island, BC, was obtained from a 30-cm tall stalagmite that grew 10.3-6.3 Ka ago, using δ¹⁸O values of the crushed calcite and of the inclusion water as inferred from its δD. It is found that the locality experienced mean annual temperature variations up to 11 °C over a 4-Ka period in the early Holocene. At the beginning of the period, local temperature quickly increased from a minimum of ∼1 °C to around 10 °C, but this early climate optimum, about 3 °C warmer than today, only lasted for ∼1200 years. About 8.6 Ka ago, temperature had declined to ∼7 °C, approximately the same as the modern cave temperature. Since then, the study area has experienced only minor temperature fluctuations, but there was a brief fall to ∼4 °C at around 7 Ka ago, which might be caused by a short lived expansion of local alpine glaciers. The long-term T-dependence of δD was 1.47‰/°C, identical to the value in modern precipitation.     

Low-temperature anisotropic diffusion of helium in zircon: Implications for zircon (U-Th)/He thermochronometry

In the last decade the zircon (U-Th)/He (ZHe) thermochronometer has been applied to a variety of geologic problems. Although bulk diffusion coefficients for He in zircon are available from laboratory step-heating experiments, little is known about the diffusion mechanism(s) and their dependence on the crystallographic structure of zircon. Here, we investigate the diffusion of He in perfectly crystalline zircon using atomistic simulation methods that provide insights into the structural pathways of He migration in zircon. Empirical force fields and quantum-mechanical calculations reveal that the energy barriers for He diffusion are strongly dependent on structure. The most favorable pathway for He diffusion is the [0 0 1] direction through the open channels parallel to the c-axis (, activation energy for tracer diffusion of a He atom along [0 0 1]). In contrast, energy barriers are higher in other directions where narrower channels for He diffusion are identified, such as [1 0 0], [1 0 1], and [1 1 0] (ΔE^∗ of 44.8, 101.7, and 421.3 kJ mol⁻¹, respectively). Molecular dynamics simulations are in agreement with these results and provide additional insight in the diffusion mechanisms along different crystallographic directions, as well as the temperature dependence. Below the closure temperature of He in zircon [T_c ∼ 180 °C, Reiners P. W., Spell T. L., Nicolescu S., and Zanetti K. A. (2004) Zircon (U-Th)/He thermochronometry: He diffusion and comparisons with Ar-40/Ar-39 dating. Geochim. Cosmochim. Acta68, 1857-1887], diffusion is anisotropic as He moves preferentially along the [0 0 1] direction, and calculated tracer diffusivities along the two most favorable directions differ by approximately five orders of magnitude (D_[001]/D_[100] ∼ 10⁵, at T = 25 °C). Above this temperature, He atoms start to hop between adjacent [0 0 1] channels, along [1 0 0] and [0 1 0] directions (perpendicular to the c-axis). The diffusion along [1 0 0] and [0 1 0] is thermally activated, such that at higher temperatures, He diffusion in zircon becomes nearly isotropic (D_[001]/D_[100] ∼ 10, at T = 580 °C). These results suggest that the anisotropic nature of He diffusion at temperatures near the closure temperature should be considered in future diffusivity experiments. Furthermore, care should be taken when making geologic interpretations (e.g., exhumation rates, timing of cooling, etc.) from this thermochronometer until the effects of anisotropic diffusion on bulk ages and closure temperature estimates are better quantified.     

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.     

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.     

Thermal history of the unsaturated zone at Yucca Mountain,Nevada, USA

Secondary calcite, silica and minor amounts of fluorite deposited in fractures and cavities record the chemistry, temperatures, and timing of past fluid movement in the unsaturated zone at Yucca Mountain, Nevada, the proposed site of a high-level radioactive waste repository. The distribution and geochemistry of these deposits are consistent with low-temperature precipitation from meteoric waters that infiltrated at the surface and percolated down through the unsaturated zone. However, the discovery of fluid inclusions in calcite with homogenization temperatures (T_h) up to ∼80 °C was construed by some scientists as strong evidence for hydrothermal deposition. This paper reports the results of investigations to test the hypothesis of hydrothermal deposition and to determine the temperature and timing of secondary mineral deposition. Mineral precipitation temperatures in the unsaturated zone are estimated from calcite- and fluorite-hosted fluid inclusions and calcite δ¹⁸O values, and depositional timing is constrained by the ²⁰⁷Pb/²³⁵U ages of chalcedony or opal in the deposits. Fluid inclusion T_h from 50 samples of calcite and four samples of fluorite range from ∼35 to ∼90 °C. Calcite δ¹⁸O values range from ∼0 to ∼22‰ (SMOW) but most fall between 12 and 20‰. The highest T_h and the lowest δ¹⁸O values are found in the older calcite. Calcite T_h and δ¹⁸O values indicate that most calcite precipitated from water with δ¹⁸O values between −13 and −7‰, similar to modern meteoric waters.     

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.     

U and Th zonation in apatite observed by laser ablation ICPMS, and implications for the (U-Th)/He system

A laser-ablation inductively-coupled plasma mass spectrometry technique was developed to measure U, Th, and Ce zonation in polished sections of apatite for assessing the consequences of parent zonation for (U-Th)/He thermochronometry. The technique produces concentration maps with an averaging length-scale of ∼20 μm, comparable to the α-stopping distance, and a precision of ∼5% down to few ppm concentration levels. A model was developed to transform the measured concentration distribution into a simplified representation for use in spherical-geometry He production-diffusion models. To illustrate these methods, 30 sections of apatite from a single granite (GC863) were mapped. Every analyzed apatite from GC863 is zoned, with most grains having variable thickness rims and terminations that are enriched in U and Th by about a factor of three over the grain cores.Parent zonation has three independent effects on (U-Th)/He He ages: it influences the α ejection correction, the ⁴He concentration profile which governs diffusive loss, and, via radiation damage trap accumulation, spatial variability of diffusivity within the crystal. If the observed zonation is typical of the apatite population in GC863, use of the standard homogenous α ejection correction would cause He ages to be on average 3% too young, and with a large amount of grain-to-grain variability (9% too young in the most rim-enriched case to 6% too old in a core-enriched case). Independent of the ejection correction, the concentration profile modifies the effective closure temperature of the apatites by placing more (or less) ⁴He near the grain edge. The parent zonation in GC863 apatites causes closure temperatures to range from four degrees lower (rim-enriched case) to two degrees higher (core-enriched case) than applies in the homogenous case. Alpha ejection and concentration profile effects on He age are additive and of the same sense. In the case of typical grains in GC863 cooled between 1 and 10 °C/Ma, the two effects are roughly equal in magnitude. The effects of intracrystalline variations in radiation damage trap accumulation become apparent at slow cooling rates (1 °C/Ma). For example, in rim-enriched GC863 grains cooled at 1 °C/Ma, preferential accumulation of radiation damage traps near the grain rim almost compensates for the higher loss rate expected of ⁴He also located preferentially near the rim. Under some circumstances strong rim-enrichment may actually increase the effective closure temperature of an apatite. Zonation at the level observed in GC863 modifies the ⁴He/³He spectra substantially from that expected from a uniform distribution. Measured ⁴He/³He spectra are strikingly similar to predictions based on the mapped eU distributions of the very same crystals, supporting the overall validity of the analytical and interpretive approach presented here.The magnitude and style of U, Th zonation documented in GC863 is one possible source of frequently observed over-dispersion of apatite (U-Th)/He ages as well as anomalous ⁴He/³He spectra.     

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.     

Ar-Ar and I-Xe ages and thermal histories of three unusual metal-rich meteorites

Portales Valley, Sombrerete, and Northwest Africa (NWA) 176 are three unrelated meteorites, which consist of silicate mixed with substantial amounts of metal and which likely formed at elevated temperatures as a consequence of early impacts on their parent bodies. Measured ³⁹Ar-⁴⁰Ar ages of these meteorites are 4477 ± 11 Ma and 4458 ± 16 Ma (two samples of Portales Valley), 4541 ± 12 Ma, and 4524 ± 13 Ma, respectively (Ma = million years; all one-sigma errors). The Ar-Ar data for Portales Valley show no evidence of later open system behavior suggested by some other chronometers. Measured ¹²⁹I-¹²⁹Xe ages of these three meteorites are 4559.9 ± 0.5 Ma, 4561.9 ± 1.0 Ma, and ∼4544 Ma, respectively (relative to Shallowater = 4562.3 ± 0.4 Ma). From stepwise temperature release data, we determined the diffusion characteristics for Ar and Xe in our samples and calculated approximate closure temperatures for the K-Ar and I-Xe chronometers. Adopting results and interpretations about these meteorites from some previous workers, we evaluated all these data against various thermal cooling models. We conclude that Portales Valley formed 4560 Ma ago, cooled quickly to below the I-Xe closure temperature, then cooled deep within the parent body at a rate of ∼4 °C/Ma through K-Ar closure. We conclude that Sombrerete formed 4562 Ma ago and cooled relatively quickly. NWA 176 likely formed and cooled quickly ∼4544 Ma ago, or later than formation times of most meteorite parent bodies. For all three meteorites, the Ar-Ar ages are in better agreement with I-Xe ages and preferred thermal models if we increase these Ar-Ar ages by ∼20 Ma. Such age corrections would be consistent with probable errors in ⁴⁰K decay parameters in current use, as suggested by others. The role of impact heating and possible disruption and partial reassembly of meteorite parent bodies to form some meteorites likely was an important process in the early solar system.     

The thermal history of equilibrated ordinary chondrites and the relationship between textural maturity and temperature

The compositions and textures of phases in eleven equilibrated ordinary chondrites from the H, L, and LL groups spanning petrographic types 4-6 were studied and used to constrain the thermal histories of their parent bodies. Based on Fe-Mg exchange between olivine and spinel, average equilibration temperatures for type 4-6 chondrites encompass a small range, 586-777 °C, relative to what is commonly assumed for peak temperatures (600-950 °C). The maximum temperatures recorded by individual chondrites, which are minima relative to peak metamorphic temperatures, increase subtly but systematically with metamorphic type and are tightly clustered for H4-6 (733-754 °C) and LL4-6 (670-777 °C). For the Ls, Ausson (L5) records a higher maximum olivine-spinel temperature (761 °C) than does the L4 chondrite Saratov (673 °C) or the L6 chondrite Glatton (712 °C). Our data combined with olivine-spinel equilibration temperatures calculated for other equilibrated ordinary chondrites using mineral compositions from the literature demonstrate that, in general, type 4 chondrites within each chemical group record temperatures lower than or equal to those of types 5-6 chondrites.For H chondrites, the olivine-spinel closure temperature is a function of spinel grain size, such that larger grains, abundant in types 5-6 chondrites, record temperatures of ∼740 °C or more while smaller grains, rare in types 5-6 but abundant in type 4 chondrites, record lower temperatures. Olivine-spinel temperatures in the type 6 chondrites Guareña and Glatton are consistent with rapid (50-100 °C/Myr) cooling from high temperatures in the ordinary chondrite parent bodies. With one exception (∼500 °C/Myr), olivine-spinel data for St.-Séverin (LL6) are consistent with similar cooling rates. Cooling rates of order 100 °C/Myr at ∼750 °C for type 6 chondrites are considerably higher than previously determined cooling rates for lower temperatures (?550 °C) based on metallography, fission tracks, and geochronology. For H chondrites, current thermal models of an “onion shell” parent body are inconsistent with a small range of peak temperatures based on olivine-spinel and two pyroxene thermometry combined with a wide dispersion of cooling rates at low temperatures. Equilibrated chondrites may have sampled regions near a major transition in physical properties such as near the base of a regolith pile.     

Geology, geochemistry and earthquake history of L?ihi Seamount, Hawaìi's youngest volcano

A half-century of investigations are summarized here on the youngest Hawaiian volcano, L?ihi Seamount. It was discovered in 1952 following an earthquake swarm. Surveying in 1954 determined it has an elongate shape, which is the meaning of its Hawaiian name. L?ihi was mostly forgotten until two earthquake swarms in the 1970s led to a dredging expedition in 1978, which recovered young lavas. The recovery of young lavas motivated numerous expeditions to investigate the geology, geophysics, and geochemistry of this active volcano. Geophysical monitoring, including a real-time submarine observatory that continuously monitored L?ihi's seismic activity for 3 months, captured some of the volcano's earthquake swarms. The 1996 swarm, the largest recorded in Hawaìi, was preceded earlier in the year by at least one eruption and accompanied by the formation of a ∼300-m deep pit crater, Pele's Pit. Seismic and petrologic data indicate that magma was stored in a ∼8-9 km deep reservoir prior to the 1996 eruption.Studies on L?ihi have altered conceptual models for the growth of Hawaiian and other oceanic island volcanoes, and refined our understanding of mantle plumes. Petrologic and geochemical studies of L?ihi lavas showed that the volcano taps a relatively primitive part of the Hawaiian plume, producing a wide range of magma compositions. These compositions have become progressively more silica-saturated with time, reflecting higher degrees of partial melting as the volcano drifts toward the center of the hotspot. Helium and neon isotopes in L?ihi glasses are among the least radiogenic found at ocean islands, and may indicate a relatively deep and undegassed mantle source for the volcano. The north-south orientation of L?ihi rift zones indicates that they may have formed beyond the gravitational influence of the adjacent older volcanoes. A new growth model indicates that L?ihi is older, taller and more voluminous than previously thought. Seismic and bathymetric data have clarified the importance of landsliding in the early formation of ocean island volcanoes. However, a fuller understanding of L?ihi's internal structure and eruptive behavior awaits installation of monitoring equipment on the volcano.The presence of hydrothermal activity at L?ihi was initially proposed based on nontronite deposits on dredged samples that indicated elevated temperatures (31 °C), water temperature, methane and ³He anomalies, and clumps of benthic micro-organisms in the water column above the volcano in 1982. Submersible observations in 1987 confirmed a low temperature geothermal system (15-30 °C) prior to the 1996 formation of Pele's Pit. The sulfide mineral assemblage (wurtzite, pyrrhotite, and chalcopyrite) deposited after the pit crater collapsed are consistent with hydrothermal fluids with temperatures >250 °C, although the highest measured temperature was ∼200 °C. Vent temperatures decreased to ∼60 °C during the 2004 dive season indicating a waning of the current phase of hydrothermal activity.     

Assessment of grain-scale homogeneity and equilibration of carbon and oxygen isotope compositions of minerals in carbonate-bearing metamorphic rocks by ion microprobe

Nineteen samples of metamorphosed carbonate-bearing rocks were analyzed for carbon and oxygen isotope ratios by ion microprobe with a ∼5-15 μm spot, three from a regional terrain and 16 from five different contact aureoles. Contact metamorphic rocks further represent four groups: calc-silicate marble and hornfels (6), brucite marble (2), samples that contain a reaction front (4), and samples with a pervasive distribution of reactants and products of a decarbonation reaction (4). The average spot-to-spot reproducibility of standard calcite analyses is ±0.37‰ (2 standard deviations, SD) for δ¹⁸O and ±0.71‰ for δ¹³C. Ten or more measurements of a mineral in a sample that has uniform isotope composition within error of measurement can routinely return a weighted mean with a 95% confidence interval of 0.09-0.16‰ for δ¹⁸O and 0.10-0.29‰ for δ¹³C. Using a difference of >6SD as the criterion, only four of 19 analyzed samples exhibit significant intracrystalline and/or intercrystalline inhomogeneity in δ¹³C at the 100-500 μm scale, with differences within individual grains up to 3.7‰. Measurements are consistent with carbon isotope exchange equilibrium between calcite and dolomite in five of six analyzed samples at the same scale. Because of relatively slow carbon isotope diffusion in calcite and dolomite, differences in δ¹³C can survive intracrystalline homogenization by diffusion during cooling after peak metamorphism and likely represent the effects of prograde decarbonation and infiltration. All but 2 of 11 analyzed samples exhibit intracrystalline differences in δ¹⁸O (up to 9.4‰), intercrystalline inhomogeneity in δ¹⁸O (up to 12.5‰), and/or disequilibrium oxygen isotope fractionations among calcite-dolomite, calcite-quartz, and calcite-forsterite pairs at the 100-500 μm scale. Inhomogeneities in δ¹⁸O and δ¹³C are poorly correlated with only a single mineral (dolomite) in a single sample exhibiting both. Because of relatively rapid oxygen isotope diffusion in calcite, intracrystalline inhomogeneities in δ¹⁸O likely represent partial equilibration between calcite and fluid during retrograde metamorphism. Calcite is in oxygen isotope exchange equilibrium with forsterite in one of four analyzed samples, in equilibrium with dolomite in none of six analyzed samples, and in equilibrium with quartz in neither of two analyzed samples. There are no samples of contact metamorphic rock with analyzed reactants and products of an arrested metamorphic reaction that are in oxygen isotope equilibrium with each other. The degree of departure from equilibrium in analyzed samples is variable and is often related, at least in part, to alteration of δ¹⁸O of calcite during retrograde fluid-rock reaction. In situ sub-grain-scale carbon and oxygen isotope analyses of minerals are advisable in the common applications of stable isotope geochemistry to metamorphic petrology. Correlation of sub-mm scale stable isotope data with imaging will lead to improved understanding of reaction kinetics, reactive fluid flow, and thermal histories during metamorphism.     

Intra-shell oxygen isotope ratios in the benthic foraminifera genus Amphistegina and the influence of seawater carbonate chemistry and temperature on this ratio

Using secondary ion mass spectrometry (SIMS) we looked at the natural variability in the oxygen isotope ratio of the shallow water, symbionts-bearing foraminiferan Amphistegina lobifera. Live foraminifera were collected in February 2005 in the Gulf of Eilat, Israel. Vertical section exposing the knob area of this species represents the growth history of this species from August 2004 to February 2005. SIMS profile at a resolution of ∼15 μm (representing about 2 weeks considering the size of the knob area and the life span of ≈6 months of this foraminifera species) yielded δ¹⁸O changes of ∼1.5‰ that are compatible with the known temperature changes for the Gulf of Eilat for this period (21-27 °C). Natural variability between primary and secondary calcite at the knob area were obtained on horizontal section of the upper knob area. This section is semi-tangential to the growth lines and exposes relatively wide belts of the primary calcite which could be analysed using the SIMS (beam size of 10 × 20 μm). The primary calcite δ¹⁸O value is on average more than 3‰ lower than the secondary calcite that represents the bulk of the skeleton (more than 95% by weight). A vertical profile at the knob was obtained by rastering an area of 50 × 50 μm at vertical steps of roughly 1 μm. The profile revealed a narrow zone of lower δ¹⁸O compared to the higher values above and below it. The difference between the lowest δ¹⁸O and the highest one was also close to 2‰. The δ¹⁸O in the margin - keel area of A. lobifera is also lower compared to the bulk secondary calcite. Specimens that were cultured in the laboratory at a constant temperature and inorganic carbon but at different pH have increased their CaCO₃ weight by roughly a factor of 8. Single specimen from each pH (ranging between 7.90 and 8.45) were investigated with the SIMS at the knob area. While there is some variability within each specimen (perhaps related to the primary calcite), the general trend was a decrease in δ¹⁸O with increasing pH (or CO₃²⁻ concentration), in agreement with previous studies on planktonic foraminifera. Some other specimens grown at different temperatures (between 21 and 33 °C) were also measured with the SIMS at the knob area. For each temperature, we observed also some variability, nevertheless the trend of −0.2‰/°C in δ¹⁸O is observed.     

Apatite (U-Th)/He thermochronometry using a radiation damage accumulation and annealing model

Helium diffusion from apatite is a sensitive function of the volume fraction of radiation damage to the crystal, a quantity that varies over the lifetime of the apatite. Using recently published laboratory data we develop and investigate a new kinetic model, the radiation damage accumulation and annealing model (RDAAM), that adopts the effective fission-track density as a proxy for accumulated radiation damage. This proxy incorporates creation of crystal damage proportional to α-production from U and Th decay, and the elimination of that damage governed by the kinetics of fission-track annealing. The RDAAM is a version of the helium trapping model (HeTM; 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, 148-161), calibrated by helium diffusion data in natural and partially annealed apatites. The chief limitation of the HeTM, now addressed by RDAAM, is its use of He concentration as the radiation damage proxy for circumstances in which radiation damage and He are not accumulated and lost proportionately from the crystal.By incorporating the RDAAM into the HeFTy computer program, we explore its implications for apatite (U-Th)/He thermochronometry. We show how (U-Th)/He dates predicted from the model are sensitive to both effective U concentration (eU) and details of the temperature history. The RDAAM predicts an effective He closure temperature of 62 °C for a 28 ppm eU apatite of 60 μm radius that experienced a 10 °C/Ma monotonic cooling rate; this is 8 °C lower than the 70 °C effective closure temperature predicted using commonly assumed Durango diffusion kinetics. Use of the RDAAM is most important for accurate interpretation of (U-Th)/He data for apatite suites that experienced moderate to slow monotonic cooling (1-0.1 °C/Ma), prolonged residence in the helium partial retention zone, or a duration at temperatures appropriate for radiation damage accumulation followed by reheating and partial helium loss. Under common circumstances the RDAAM predicts (U-Th)/He dates that are older, sometimes much older, than corresponding fission-track dates. Nonlinear positive correlations between apatite (U-Th)/He date and eU in apatites subjected to the same temperature history are a diagnostic signature of the RDAAM for many but not all thermal histories.Observed date-eU correlations in four different localities can be explained with the RDAAM using geologically reasonable thermal histories consistent with independent fission-track datasets. The existence of date-eU correlations not only supports a radiation damage based kinetic model, but can significantly limit the range of acceptable time-temperature paths that account for the data. In contrast, these datasets are inexplicable using the Durango diffusion model. The RDAAM helps reconcile enigmatic data in which apatite (U-Th)/He dates are older than expected using the Durango model when compared with thermal histories based on apatite fission-track data or other geological constraints. It also has the potential to explain at least some cases in which (U-Th)/He dates are actually older than the corresponding fission-track dates.     

Calcium isotope fractionation in calcite and aragonite

Calcium isotope fractionation was measured on skeletal aragonite and calcite from different marine biota and on inorganic calcite. Precipitation temperatures ranged from 0 to 28°C. Calcium isotope fractionation shows a temperature dependence in accordance with previous observations: 1000 · ln(α_cc) = −1.4 + 0.021 · T (°C) for calcite and 1000 · ln(α_ar) = −1.9 + 0.017 · T (°C) for aragonite. Within uncertainty the temperature slopes are identical for the two polymorphs. However, at all temperatures calcium isotopes are more fractionated in aragonite than in calcite. The offset in δ^44/40Ca is about 0.6‰. The underlying mechanism for this offset may be related to the different coordination numbers and bond strengths of the calcium ions in calcite and aragonite crystals, or to different Ca reaction behavior at the solid-liquid interface. Recently, the observed temperature dependence of the Ca isotope fractionation was explained quantitatively by the temperature control on precipitation rates of calcium carbonates in an experimental setting (Lemarchand et al., 2004). We show that this mechanism can in principle also be applied to CaCO₃ precipitation in natural environments in normal marine settings. Following this model, Ca isotope fractionation in marine Ca carbonates is primarily controlled by precipitation rates. On the other hand the larger Ca isotope fractionation of aragonite compared to calcite can not be explained by different precipitation rates. The rate control model of Ca isotope fractionation predicts a strong dependence of the Ca isotopic composition of carbonates on ambient CO₃²⁻ concentration. While this model is in general accordance with our observations in marine carbonates, cultured specimens of the planktic foraminifer Orbulina universa show no dependence of Ca-isotope fractionation on the ambient CO₃²⁻ concentration. The latter observation implies that the carbonate chemistry in the calcifying vesicles of the foraminifer is independent from the ambient carbonate ion concentration of the surrounding water.     

Comparison between grain size and multi-mineral Ar/Ar thermochronology

In the case of volume diffusion, the closure temperature of a mineral is function of, among other factors, the characteristic diffusion dimension, which can be approximated by the grain size of the mineral analysed for grains smaller than or similar in size to the diffusion domains. The theoretical possibility of single mineral grain size thermochronology had been demonstrated empirically in earlier studies, mostly using biotite. In order to examine the potential of this method, it was tested alongside the widely used multi-mineral ⁴⁰Ar/³⁹Ar thermochronology. The sample comes from the granitic McLean pluton, in the south section of the Grenville orogeny. Seven grain size separates of biotite (ranging between 90 and 1000 μm), eight size fractions of amphibole (between 63 and 1000 μm), and three size fractions of K-feldspar (250-600 μm) were extracted and dated by the laser step-heating ⁴⁰Ar/³⁹Ar method. The total gas ages obtained behave as theoretically predicted, with increasing ages for increasing grain sizes, including for K-feldspar, but with the exception of the smallest and the largest grains for biotite and amphibole. The calculated cooling rates are ca. 0.7 °C/Ma for K-feldspar, ca. 2.5 °C/Ma for biotite, and ca. 11 °C/Ma for amphibole, corresponding very well to a monotonic cooling of the McLean pluton. A quick initial thermal re-equilibration with the cooler host-rocks is followed by a much slower cooling on a thermal path parallel to that of the Frontenac Terrain situated immediately to the southeast. The validity of the single mineral grain size thermochronology is demonstrated by comparison with the thermal evolution of the adjacent units and with the cooling history derived from a multi-mineral thermochronology, suggesting that it can be routinely used. The application of this method can be hampered by insufficiently low analytical uncertainties.     

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.