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.     

α-Emitting mineral inclusions in apatite, their effect on (U-Th)/He ages, and how to reduce it

U-Th rich mineral inclusions in apatite are often held responsible for erroneously old (U-Th)/He ages, because they produce “parentless” He. Three aspects associated with this problem are discussed here. First, simple dimensional considerations indicate that for small mineral inclusions, the parentless helium problem might not be as serious as generally thought. For example, a mineral inclusion that is 10% the length, width and height of its host apatite needs to be a thousand times more concentrated in U and Th to produce an equal amount of He. Therefore, single isolated inclusions smaller than a few μm are unlikely to contribute significant helium. For larger or more abundant inclusions, the parentless helium problem can be solved by dissolution of the apatite and its inclusions in hot HF. Second, besides creating parentless helium, inclusions also complicate α-ejection corrections. Mathematical exploration of this latter problem for spherical geometries reveals that for randomly distributed inclusions, the probability distribution of single-grain ages is predicted to have a sharp mode at the mean age, with tails towards younger and older ages. Multiple-grain measurements will yield accurate and precise age estimates if 10 or more randomly distributed α-emitting mineral inclusions are present in a sample. Third, thermal modeling indicates that mineral inclusions have a non-trivial but minor (<5 °C) effect on the closure temperature. These predictions were tested on apatites from rapidly cooled migmatites of Naxos (Greece) which contain abundant U-rich zircon inclusions. Thirty-seven samples were subjected to two kinds of treatment. The “pooled” age (i.e., the synthetic multi-grain age computed from a number of single-grain analyses) of four inclusion-free samples (13 apatites), prepared in HNO₃ is 10.9 Ma, close to apatite and zircon fission-track ages from the same rock. (U-Th)/He ages of 14 inclusion-bearing samples dissolved in HNO₃ range between 9 and 45 Ma, with a pooled age of 22.6 Ma. The ages of 19 HF-treated samples range between 5 and 16 Ma, with 10 of 14 single-grain samples between 9 and 13 Ma and a pooled age of 10.9 Ma. These observations agree with the theoretical predictions and support the addition of HF-treated apatite (U-Th)/He dating to the thermochronological toolbox.     

Interpreting data dispersion and “inverted” dates in apatite (U-Th)/He and fission-track datasets: An example from the US midcontinent

New apatite (U-Th)/He (AHe) and apatite fission-track (AFT) data were acquired for cratonic basement samples from an 80 m span of drillcore in northeastern Kansas. The short depth interval over which the samples were collected indicates that they should have undergone thermal histories that would be indistinguishable using low temperature thermochronometry techniques. Individual AHe dates from four samples range from 99 to 464 Ma. Three samples yield dates <300 Ma that display a correlation with apatite eU (9-34 ppm) and a weaker correlation with grain size. eU concentration maps of apatites from these samples reveal low to moderate zonation in eU. Results for a fourth sample are characterized by dates >300 Ma, higher eU (39-113 ppm), and substantial data dispersion uncorrelated with eU and grain size. These apatites have strong and variable eU zonation. AFT dates for five samples range from 242 to 291 Ma. The sample with the highest eU apatites and oldest AHe dates yields the youngest AFT results. These results are “inverted”, with AHe dates distinctly older than the corresponding AFT date.We explore both the causes of data dispersion and the overall compatibility of this cratonic dataset. We find that geologically reasonable thermal histories can (1) explain the distribution of the moderate eU AHe data when accounting for the influence of radiation damage, grain size, and eU zonation on apatite He diffusivity, (2) reproduce the observed dispersion in the high eU AHe data when using a viable range of eU zonation and grain size, and (3) explain the AFT data for the same samples. The AHe and AFT data are mutually consistent, and viable thermal histories successfully predict the observed pattern of older AHe than AFT dates for the high eU apatites. Together these results suggest that appropriately accounting for the known controls on apatite He diffusivity can explain the observed dispersion and “inverted” AHe and AFT results in some thermochronometry datasets. A range of AHe dates should be especially common in cratonic data, because small differences in apatite He diffusivity are amplified by the thermal histories that typify cratonic settings. We use these results to develop some guidelines for interpreting dispersed AHe datasets. First, date-eU and date-grain size correlations should be evaluated, and if these patterns occur they can be used to better resolve the thermal history. Second, for samples that yield inexplicably large dispersion of AHe dates uncorrelated with eU and crystal size, the appropriate strategy is either to reject these samples from the suite used for thermal history interpretation or to acquire additional data to help decipher the significance of the age distribution.     

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.     

Laser microprobe (U-Th)/He geochronology

A new analytical method had been developed to enable high-spatial-resolution (U-Th)/He dating of accessory minerals. It involves the use of a focused ArF excimer to ablate pits in a polished grain surface, with the evolved gases spiked for isotope-dilution measurement of radiogenic ⁴He. These data are converted to concentrations by precise measurement of each pit using an optical interferometric microscope. U, Th, and Sm concentration measurements are made using one of several alternative microanalytical techniques (e.g., wavelength-dispersive electron microprobe analysis or laser-ablation, inductively coupled plasma mass spectrometry). By way of illustration, we present both conventional and laser microprobe (U-Th)/He dating results for a Brazilian monazite sample. Laser microprobe data (28 measurements on two crystal fragments) yield a weighted mean (U-Th)/He date of 455.3 ± 3.7 Ma (2SE). This result is statistically indistinguishable from the mean conventional (U-Th)/He date for three separate grain fragments: 449.6 ± 9.8 Ma (2SE). The agreement of conventional and laser ablation dates should encourage a wide variety of applications of the technique, including: (1) detrital mineral dating for provenance and unroofing studies; (2) the dating of broken, included, highly zoned, or irregular grains which are not easily corrected for α-ejection; and (3) measuring ⁴He loss profiles that can be inverted to determine cooling histories.     

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.     

矿物中4He同位素提取和含量测定研究

(U-Th)/He定年是一种有效的低温热年代学定年技术,现已被广泛应用于地质研究的各个领域,而矿物中4He同位素的有效提取和含量准确测定是该技术的关键。磷灰石和锆石是(U-Th)/He定年最常用的矿物,其4He提取条件及铀钍含量测定方法都较为成熟;而其他矿物(如磁铁矿、橄榄石、针铁矿、石榴子石等)的研究则相对较少。文章介绍了当前国内外(U-Th)/He研究中采用的4He同位素提取方法———真空炉加热法和激光加热法,激光加热法因具有低4He背景值和耗时短的优点而成为主要的提取方法。以磷灰石样品测试为例,介绍了成都地质矿产研究所建立的采用激光加热法和四极杆质谱提取4He同位素及其含量测量过程、含量计算和校正方法。指出未来(U-Th)/He测试技术除继续改进现有分析方法外,应加强对更多不同矿物的测试研究。     

Combined U-Th/He and Ar/Ar geochronology of post-shield lavas from the Mauna Kea and Kohala volcanoes, Hawaii

Late Quaternary, post-shield lavas from the Mauna Kea and Kohala volcanoes on the Big Island of Hawaii have been dated using the ⁴⁰Ar/³⁹Ar and U-Th/He methods. The objective of the study is to compare the recently demonstrated U-Th/He age method, which uses basaltic olivine phenocrysts, with ⁴⁰Ar/³⁹Ar ages measured on groundmass from the same samples. As a corollary, the age data also increase the precision of the chronology of volcanism on the Big Island. For the U-Th/He ages, U, Th and He concentrations and isotopes were measured to account for U-series disequilibrium and initial He. Single analyses U-Th/He ages for Hamakua lavas from Mauna Kea are 87 ± 40 to 119 ± 23 ka (2σ uncertainties), which are in general equal to or younger than ⁴⁰Ar/³⁹Ar ages. Basalt from the Polulu sequence on Kohala gives a U-Th/He age of 354 ± 54 ka and a ⁴⁰Ar/³⁹Ar age of 450 ± 40 ka. All of the U-Th/He ages, and all but one spurious ⁴⁰Ar/³⁹Ar ages conform to the previously proposed stratigraphy and published ¹⁴C and K-Ar ages. The ages also compare favorably to U-Th whole rock-olivine ages calculated from ²³⁸U-²³⁰Th disequilibria. The U-Th/He and ⁴⁰Ar/³⁹Ar results agree best where there is a relatively large amount of radiogenic ⁴⁰Ar (>10%), and where the ⁴⁰Ar/³⁶Ar intercept calculated from the Ar isochron diagram is close to the atmospheric value. In two cases, it is not clear why U-Th/He and ⁴⁰Ar/³⁹Ar ages do not agree within uncertainty. U-Th/He and ⁴⁰Ar/³⁹Ar results diverge the most on a low-K transitional tholeiitic basalt with abundant olivine. For the most alkalic basalts with negligible olivine phenocrysts, U-Th/He ages were unattainable while ⁴⁰Ar/³⁹Ar results provide good precision even on ages as low as 19 ± 4 ka. Hence, the strengths and weaknesses of the U-Th/He and ⁴⁰Ar/³⁹Ar methods are complimentary for basalts with ages of order 100-500 ka.     

（U-Th）/He定年——低温热年代学研究的一种新技术

近年来（U-Th）/He定年方法在技术和低温热年代学应用上的优势吸引人们以此开展构造、地形地貌演化等的研究；与裂变径迹、K Ar、Ar Ar等定年方法的比较，该方法是有效的；人们研究了温度及矿物粒径对He扩散的影响，确定了校正方法。（U-Th）/He的封闭温度低于其它体系（磷灰石的仅为75℃±）；有质谱仪只需数十毫克样品就能获得 U、Th、He含量；有较裂变径迹法快速、省样、省力的优点；该技术虽还有待完善，但仍不失为开展低温热年代学和低温热演化研究的潜在有效方法。     

Duluth Complex FC1 Apatite and Zircon: Reference Materials for (U-Th)/He Dating?

The accuracy and validation of geo- and thermochronological dating hinges on the availability of well-characterised age reference materials. The Mesoproterozoic gabbroic anorthosite FC1 from the Duluth Complex, Minnesota is a reference material for zircon U-Pb and a suggested reference material for apatite fission-track dating. We evaluate FC1 as (U-Th)/He reference material, and determine its apatite U-Pb, and zircon and apatite (U-Th)/He age. Our dating results constrain the thermal history of FC1, showing that fast cooling occurred between ~ 1099 and 1040 Ma from ≥ 600 °C to ~ 200 °C. The zircon (U-Th)/He data from air-abraded grains give a robust isochron age of 1037 ± 25 Ma (2s) without overdispersion. The within-grain homogeneity of U and Th, the availability of FC1 zircon, and the absence of radiation-damage effects on the (U-Th)/He age support its use as reference material. Unabraded zircon grains give lower and more dispersed ages, highlighting the usefulness of air abrasion to control for α-ejection in (U-Th)/He dating. Our apatite (U-Th-Sm)/He single-grain ages vary between 180 and 300 Ma. Their wide dispersion argues against the use of FC1 apatite as (U-Th-Sm)/He reference material and makes the interpretation of their low-temperature thermal history complicated.     

Petrogenesis of oligocene plutonic rocks in western Anatolia (NW Turkey): Insights from mineral and rock chemistry,Sr-Nd isotopes,and U-Pb,Ar-Ar and (U-Th)/He geochronology

Northwestern Anatolia is characterized by voluminous Paleozoic to Cenozoic granitoid bodies with varying compositions. Most of them are composite plutons emplaced into western Anatolia orogenic crust during the Eocene, Oligocene and Miocene along the İzmir-Ankara-Erzincan suture zone. This paper reports systematic good quality mineral and bulk-rock chemistry, Sr-Nd isotope data, honblend Ar-Ar, zircon U-Pb and first apatite (U-Th)/He (AHe) ages to reveal possible source compositions of the Evciler and Eybek granitoids and petrogenetic/geodynamic processes involved during their genesis, and thermochronology of Oligocene magmatism in the NW Anatolia. The Evciler and Eybek granitoids are mainly granodiorite and composed of K-feldspar (usually orthoclase and rarely microcline), plagioclase (albite, oligoclase), hornblende, biotite, quartz and accessory minerals (e.g., titanite, zircon, apatite, opaque), and secondary minerals such as chlorite, sericite and clay minerals. Estimated temperature-pressure conditions are 690–770 ° C at 1.6–2.7 kbar for the Evciler granitoid and 690–760 ° C at 3.2–4.01 kbar for the Eybek granitoid. These two granitoids enriched in LILEs (e.g., U, Th, Rb, and K), LREEs and Pb, and depleted in HREEs (e.g., Nb, Ti) and Sr, Ba and P relative to LILEs, and display small negative Eu anomalies. They belong to calc-alkaline, high-K calc-alkaline and minor shoshonite series, and display metaluminous and I-type character. Their REE patterns show a large fractionation between LREE and HREE ((La/Yb)N = 4.6–21.4) and a small negative Eu anomaly (Eu* = 0.2–0.3). The Evciler granitoid has homogeneous ⁸⁷Sr/⁸⁶Sr = 0.7060−0.7063 and ¹⁴³Nd/¹⁴⁴Nd = 0.51259−0.51262, and the Eybek granitoid has ⁸⁷Sr/⁸⁶Sr = 0.7060−0.7080 and ¹⁴³Nd/¹⁴⁴Nd = 0.51243−0.51263. New precise ⁴⁰Ar/³⁹Ar age data of hornblende and ²⁰⁶Pb–²³⁸U ages of zircons and (U-Th)/He ages of apatites from the plutons allow a more accurate temporal reconstruction of the Cenozoic magmatism of the western Anatolia. ⁴⁰Ar/³⁹Ar dating of hornblendes from the Evciler and Eybek granitoids gave plateau ages of between ca. 28 Ma and 25 Ma. Laser ablation inductively coupled plasma mass spectrometer (LA-ICP-MS) ²⁰⁶Pb-²³⁸Pb ages of euhdral magmatic zircons from the samples of these granitoids yield between ca. 28 and 26 Ma. The new high-temperature age constraints indicate Oligocene emplacement ages for the two intrusive bodies. The closeness of the zircon U-Pb and the hornblende Ar-Ar ages show that they experienced quick post-crystallization cooling. However, the significant difference between the apatite (U-Th)/He ages of 19.8 Ma and 7.6 Ma obtained on the Evciler and Eybek granitoids warns that in the post-Oligocene times the two structural blocks had different exhumation histories.     

Polyphase movement on the Lavanttal Fault Zone (Eastern Alps): reconciling the evidence from different geochronological indicators

The Lavanttal Fault Zone (LFZ) is generally considered to be related to Miocene orogen-parallel escape tectonics in the Eastern Alps. By applying thermochronological methods with retention temperatures ranging from ~450 to ~40°C we have investigated the thermochronological evolution of the LFZ and the adjacent Koralm Complex (Eastern Alps). ⁴⁰Ar/³⁹Ar dating on white mica and zircon fission track (ZFT) thermochronology were carried out on host rocks (HRs) and fault-related rocks (cataclasites and fault gouges) directly adjacent to the unfaulted protolith. These data are interpreted together with recently published apatite fission track (AFT) and apatite (U-Th)/He ages. Sample material was taken from three drill cores transecting the LFZ. Ar release spectra in cataclastic shear zones partly show strongly rejuvenated incremental ages, indicating lattice distortion during cataclastic shearing or hydrothermal alteration. Integrated plateau ages from fault rocks (~76 Ma) are in parts slightly younger than plateau ages from HRs (>80 Ma). Incremental ages from fault rock samples are in part highly reduced (~43 Ma). ZFT ages within fault gouges (~65 Ma) are slightly reduced compared to the ages from HRs, and fission tracks show reduced lengths. Combining these results with AFT and apatite (U-Th)/He ages from fault rocks of the same fault zone allows the recognition of distinct faulting events along the LFZ from Miocene to Pliocene times. Contemporaneous with this faulting, the Koralm Complex experienced accelerated cooling in Late Miocene times. Late-Cretaceous to Palaeogene movement on the LFZ cannot be clearly proven. ⁴⁰Ar/³⁹Ar muscovite and ZFT ages were probably partly thermally affected along the LFZ during Miocene times.     

Electrodynamic Disaggregation: Does it Affect Apatite Fission-Track and (U-Th)/He Analyses?

Apatite fission-track and (U-Th)/He analyses require the liberation of intact idiomorphic apatite grains from rock samples. While routinely being carried out by mechanical methods, electrodynamic disaggregation (ED) offers an alternative approach. The high-voltage discharges produced during the ED process create localised temperature peaks (10000 K) along a narrow plasma channel. In apatite, such high temperatures could potentially reduce the length of fission tracks, which start to anneal at temperatures > 60 °C, and could also enhance He diffusion, which becomes significant at 30–40 °C over geological time scales. A comparison of fission-track analyses and (U-Th)/He ages of apatites prepared both by mechanical (jaw crusher, disk mill) and ED processing provides a way of determining whether heating during the latter method has any significant effect. Apatites from three samples of different geological settings (an orthogneiss from Madagascar, the Fish Canyon Tuff, and a muscovite-gneiss from Greece) yielded statistically identical track length distributions compared to samples prepared mechanically. Additionally, (U-Th)/He ages of apatites from a leucogranite from Morocco prepared by both methods were indistinguishable. These first results indicated that during electrodynamic disaggregation apatite crystals were not heated enough to partially anneal the fission tracks or induce significant diffusive loss of He.     

Oxygen and carbon isotope composition of structural carbonate in weathering apatites from laterites, southern Brazil and western Senegal

The oxygen (δ¹⁸O_c) and carbon (δ¹³C_c) isotope compositions of the structural carbonate group (CO₃) in apatites from lateritic profiles were investigated. The weathering profiles, located in southern Brazil and in western Senegal, are developed on three different types of apatite-rich parent rock: carbonatite, metamorphosed marine phosphorite and sedimentary marine phosphorite. The parent rock apatites are of magmatic, hydrothermal, metamorphic and sedimentary origins. The in situ formation of apatite of weathering origin in the profiles is well documented petrographically and geochemically.The overall range of measured δ¹⁸O_c and δ¹³C_c values of apatites of weathering origin (22 to 27 SMOW for δ¹⁸O_c and −15 to −10 PDB for δ¹³C_c) is much smaller than the range of measured and/or published isotope compositions of parent rock apatites (4–35 for δ¹⁸O_c and −11 to +1 for δ¹³C_c). In any profile, the apatites of weathering origin can exhibit lower, similar or higher δ¹⁸O_c values than parent rock apatites. In contrast, their δ¹³C_c values are systematically and significantly lower than those of the parent rock apatites. Apatites formed as a result of weathering in laterites can therefore be readily distinguished from apatites of other origin on the basis of their isotope composition.Assuming that apatite CO₃ fractionates O in a way similar to calcite CO₃, the structural carbonate group of the apatites of weathering origin appears to form in approximate isotopic equilibrium with the weathering solutions. The very low δ¹³C_c values exhibited by these apatites indicate that the dominant sources of dissolved CO₂ in the soil water are organic. The isotope composition of structural carbonate in apatite of weathering origin in lateritic profiles may provide useful information for paleoenvironmental studies.     

Apatite Reference Materials for SIMS Microanalysis of Isotopes and Trace Elements

Twelve apatite samples have been tested as secondary ion mass spectrometry (SIMS) reference materials. Laser ablation-inductively coupled plasma-mass spectrometry (LA-ICP-MS) analysis shows that the SLAP, NUAN and GR40 apatite gems are internally homogeneous, with most trace element mass fractions having 2 standard deviations (2s) ≤ 2.0%. BR2, BR5, OL2, AFG2 and AFB1, which have U > 63 μg g^-1, ²⁰⁶Pb/²⁰⁴Pb > 283, and homogeneous SIMS U-Pb data, have respective isotope dilution thermal ionisation mass spectrometry (ID-TIMS) ages of 2053.83 ± 0.21 Ma, 2040.34 ± 0.09 Ma, 868.87 ± 0.25 Ma, 478.71 ± 0.22 Ma and 473.25 ± 0.09 Ma. Minor U-Pb heterogeneity exists and accurate SIMS results require correction with the 3D Concordia-constrained common Pb composition. Among the studied samples, AFG2 and BR5 are the most homogeneous U-Pb reference materials. The SIMS sulfur isotopic compositions of eight of the apatites shows they are homogeneous, with 2s for both 10³δ³⁴S and 10³δ³³S < 0.55‰. One apatite, BR96, has Δ³³S = -0.36 ± 0.2‰. The apatite samples have ID-TIMS ⁸⁷Sr/⁸⁶Sr between 0.704214 ± 0.000030 and 0.723134 ± 0.000035.     

Interpretation of (U–Th)/He single grain ages from slowly cooled crustal terranes: A case study from the Transantarctic Mountains of southern Victoria Land

Low temperature thermochronologic techniques (e.g. apatite fission track (AFT) thermochronology and (U–Th)/He dating) constrain near-surface T–t paths and are often applied to uplift/denudation and landscape evolution studies. Samples collected in vertical profiles from granitic walls on either side of the Ferrar Glacier, southern Victoria Land, Antarctica were analyzed using AFT thermochronology and apatite (U–Th)/He dating to further constrain the lowest temperature thermal history of this portion of the Transantarctic Mountains. AFT central ages vary systematically with elevation and together with track length information define a multi-stage cooling/denudation history in the Cretaceous and early Tertiary. Apatite (U–Th)/He single grain age variation with elevation is not as systematic with considerable intra-sample age variation. Although many complicating factors (e.g., U- and Th-rich (micro)inclusions, fluid inclusions, variation in crystal size, α-particle ejection correction, zonation and α-particle ejection correction, implantation of He into a crystal or impediment of He diffusion out of a crystal, and ¹⁴⁷Sm-derived α-particles) may contribute to age dispersion, we found that variation in single grain ages correlated with cooling rate. Samples that cooled relatively quickly have less variation in single grain ages, whereas samples that cooled relatively slowly (< 3 °C/m.y.) or resided within an (U–Th)/He partial retention zone (HePRZ) prior to more rapid cooling have a comparatively greater variation in ages.Decay of U and Th via α-particle emission creates a ⁴He concentration profile dependent upon the initial parent [U,Th] within a crystal. Variation of single grain ages for samples with non-homogeneous [U,Th] distributions will be enhanced with long residence time in the partial retention zone (i.e., slow cooling) because of the relative importance of loss via volume diffusion and loss via α-particle ejection with respect to the [U,Th] zonation and the grain boundary. Correction of ages for α-particle ejection (F_T correction factor) typically assumes uniform U and Th distribution within the crystal and when applied to a population of crystals with different U and Th distributions will enhance the variation in ages. Most complicating factors (listed above) for apatite (U–Th)/He ages result in ages that are “too old”. We propose that if considerable variation in (U–Th)/He single grain ages exists, that a weighted mean age is determined once outlier single crystal ages are excluded using the criterion of Chauvenet or a similar approach. We suggest that the “true age” or most representative age for that age population lies between the minimum (U–Th)/He age and the weighted mean age. We apply this approach, coupled with composite age profiles to better constrain the T–t history of the profiles along the Ferrar Glacier. Significant intra-sample variation in single crystal apatite (U–Th)/He ages and other minerals dated by the (U–Th)/He method should be expected, especially when the cooling rate is slow. The variation of (U–Th)/He single crystal ages is therefore another parameter that can be used to constrain low-temperature thermal histories.     

Apatite fission track and (U-Th)/He thermochronology from the Archean Tanzania Craton: Contributions to cooling histories of Tanzanian basement rocks

Borehole and surface samples from the Archean Tanzania Craton were analysed for apatite fission track(AFT) and(U-Th)/He data with the aim of deciphering cooling histories of the basement rocks. Fission track dates from borehole and outcrop samples are Carboniferous-Permian(345± 33.3 Ma to271±31.7 Ma) whereas(U-Th)/He dates are Carboniferous-Triassic(336±45.8 Ma to 213±29 Ma) for outcrop grains and are consistently younger than corresponding AFT dates. Single grain(U-Th)/He dates from the borehole are likely to be flawed by excessive helium implantation due to their very low effective uranium contents, radiation damage and grain sizes. All AFT and(U-Th)/He dates are significantly younger than the stratigraphic ages of their host rocks, implying that the samples have experienced Phanerozoic elevated paleo-temperatures. Considerations of the data indicate removal of up to 9 km overburden since the Palaeozoic.Thermal modelling reveals a protracted rapid cooling event commencing during the early Carboniferous(ca. 350 Ma) at rates of 46 m/Ma ending in the Triassic(ca. 220 Ma). The model also suggests minor cooling during the Cretaceous of the samples to surface temperatures. The suggested later cooling event remains to be tested. The major cooling phase during the Carboniferous is interpreted to be associated with compressional tectonics during the Variscan Orogeny sensu far field induced stresses. Coeval sedimentation in the Karoo basins in the region suggests that most of the cooling of cratonic rocks during the Carboniferous was associated with denudation.     

Monazite, iron oxide and barite exsolutions in apatite aggregates from CCSD drillhole eclogites and their geological implications

We have identified abundant exsolutions in apatite aggregates from eclogitic drillhole samples of the Chinese Continental Scientific Drilling (CCSD) project. Electron microscope and laser Raman spectroscopy analyses show that the apatite is fluorapatite, whereas exsolutions that can be classified into four types: (A) platy to rhombic monazite exsolutions; (B) needle-like hematite exsolutions; (C) irregular magnetite and hematite intergrowths; and (D) needle-like strontian barite exsolutions. The widths and lengths of type A monazite exsolutions range from about 6-10 μm (mostly 6 μm) and about 50-75 μm, respectively. Type B exsolutions are parallel with the C axis of apatite, with widths ranging from 0.5 to 2 μm, with most around 1.5 μm, and lengths that vary dramatically from 6 to 50 μm. Type C exsolutions are also parallel with the C axis of apatite, with lengths of ∼30-150 μm and widths of ∼10 to 50 μm. Type D strontian barite exsolutions coexist mostly with type B hematite exsolutions, with widths of about 9 μm and lengths of about 60-70 μm. Exsolutions of types B, C and D have never been reported in apatites before. Most of the exsolutions are parallel with the C axis of apatite, implying that they were probably exsolved at roughly the same time. Dating by the chemical Th-U-total Pb isochron method (CHIME) yields an U-Pb isochron age of 202 ± 28 Ma for monazite exsolutions, suggesting that these exsolutions were formed during recrystallization and retrograde metamorphism of the exhumed ultrahigh pressure (UHP) rocks. Quartz veins hosting apatite aggregates were probably formed slightly earlier than 202 Ma. Abundant hematite exsolutions, as well as coexistence of magnetite/hematite and barite/hematite in the apatite, suggest that the oxygen fugacity of apatite aggregates is well above the sulfide-sulfur oxide buffer (SSO). Given that quartz veins host these apatite aggregates, they were probably deposited from SiO₂-rich hydrous fluids formed during retrogression of the subducted slab. Such SiO₂-rich hydrous fluids may act as an oxidizing agent, a feasible explanation for the high oxygen fugacity in convergent margin systems.     

Mapping a geothermal anomaly using apatite (U‐Th)/He thermochronology in the Têt fault damage zone,eastern Pyrenees,France

(U‐Th)/He ages on apatite obtained in the vicinity of the Têt fault hydrothermal system show a large variability. In the inner damage zone adjacent to the fault core, where fluid flows are concentrated, AHe ages display a large scatter (3–41 Ma) and apatite ageing. Samples from the outer damage zone show young ages with less dispersion (0.9–21.1 Ma) and apatite rejuvenation. Outside the damage zone, ages are consistent with the regional exhumation history between 20 and 12 Ma. The important age dispersion found in the damage zone is interpreted as the result of ⁴He mobility during fluid infiltration. Our results show that thermochronological data close to a fault should be interpreted with caution, but may offer a new tool for geothermal exploration.     

Timing and mechanisms of Central Himalayan exhumation: discriminating between tectonic and erosion processes

The ability to deduce exhumation mechanisms from thermochronological data is hampered by the fact that assumptions on the thermal state of the lithosphere have to be made. Additional argumentation is generally required to discriminate between erosion-controlled and tectonically induced exhumation. This problem can be overcome by studying the spatial distribution of zircon and apatite (U-Th)/He and fission track data. In this work the variation of four different low temperature isotopic systems generating age trends along a sampling line is used to infer mechanisms of Quaternary exhumation in the Central High Himalayan Metamorphic Belt. Observed zircon age trends with southwards increasing cooling ages (from 0.5 to 1.7 Ma) are attributed to tectonically induced exhumation. The uniform apatite cooling ages clustered c. 0.5 Ma are attributed to erosion.