Interlaboratory comparison of cosmogenic 21Ne in quartz

Pieter Vermeescha,b, Greg Balcoc, Pierre-Henri Blarde, Tibor J. Dunaif, Florian Kobera, Samuel Niedermanng, David L. Shusterd,c, Stefan Straskya, Finlay M. Stuarth, Rainer Wielera, Laurent Zimmermanne

aInstitute of Geochemistry and Petrology, ETH Zurich, Zürich, Switzerland bLondon Geochronology Centre, University College London, London, United Kingdom cBerkeley Geochronology Center, Berkeley, United States dDepartment of Earth and Planetary Science, University of California, Berkeley, United States eCentre de Recherches Pétrologiques et Géochimiques, Vandoeuvre-lès-Nancy, France fUniversity of Cologne, Köln, Germany gDeutsches GeoForschungsZentrum GFZ, Potsdam, Germany hScottish Universities Environmental Research Centre, Glasgow, United Kingdom

_________________________________________________________________

Abstract

We performed an interlaboratory comparison study with the aim to determine the accuracy of cosmogenic 21Ne measurements in quartz. CREU-1 is a natural quartz standard prepared from amalgamated vein clasts which were crushed, thoroughly mixed, and sieved into 125-250 μm and 250-500μm size fractions. 50 aliquots of CREU-1 were analyzed by five laboratories employing six different noble gas mass spectrometers. The released gas contained a mixture of 16-30% atmospheric and 70-84% non-atmospheric (predominantly cosmogenic) 21Ne, defining a linear array on the 22Ne/20Ne-21Ne/20Ne three isotope diagram with a slope of 1.108±0.014. The internal reproducibility of the measurements is in good agreement with the formal analytical precision for all participating labs. The external reproducibility of the 21Ne concentrations between labs, however, is significantly overdispersed with respect to the reported analytical precision. We report an average reference concentration for CREU-1 of 348±10×106at[21Ne]/g[SiO 2], and suggest that the 7.1% (2σ) overdispersion of our measurements may be representative of the current accuracy of cosmogenic 21Ne in quartz. CREU-1 was tied to CRONUS-A, which is a second reference material prepared from a sample of Antarctic sandstone. We propose a reference value of 320±11×106at/g for CRONUS-A. The CREU-1 and CRONUS-A intercalibration materials may be used to improve the consistency of cosmogenic 21Ne to the level of the analytical precision.

_________________________________________________________________

1. Introduction

Cosmogenic neon is a relatively little used tool for studying Earth surface processes. It is powerful for four reasons. First, it is produced and retained in quartz (Niedermann et al.19931994Shuster and Farley2005) as well as most other silicates, such as pyroxene (Schäfer et al.1999), olivine (Poreda and Cerling1992), sanidine (Kober et al.2005), hornblende and biotite (Amidon and Farley2012). Therefore, it is applicable to most rock types found on the Earth’s surface. Second, cosmogenic 21Ne is a stable nuclide. This gives it an age range limited essentially only by the erosion rate and allows exceptionally old landscapes to be dated (Schäfer et al.1999Dunai et al.2005). Third, neon has three isotopes (20Ne, 21Ne, and 22Ne), each of which have different abundances in the various reservoirs (atmospheric, nucleogenic, or magmatic) that may contribute to the natural 21Ne background (Niedermann2002). By simultaneously analyzing all three isotopes and verifying whether they plot on a mixing line between atmospheric and spallogenic components, the cosmogenic neon method provides an internal ‘reliability check’ which is absent from other commonly used nuclides. Fourth, neon can be measured using a standard sector field noble gas mass spectrometer. Sample requirements are modest (typically 100-200 mg) and sample preparation is relatively straightforward as it does not require extensive chemical purification or chromatography. This greatly increases sample throughput, which in turn opens up exciting opportunities for detrital work (Dunai et al.2005Codilean et al.2008).

Cosmogenic 21Ne is even more useful when combined with one or more cosmogenic radionuclides such as 10Be or 26Al. Such double- or triple-dating may be used for burial dating (Balco and Shuster2009aVermeesch et al.2010), for catchment-wide erosion studies with complex exposure histories (Kober et al.2009), or to measure the exposure age of old and very slowly eroding surfaces (Fujioka et al.2005). An implicit assumption of many of these studies is that the accuracy of the 21Ne method equals its analytical precision. Violation of this assumption may lead to erroneous results such as samples plotting in the ‘forbidden zone’ of the 21Ne/10Be two-nuclide diagram (Lal1991Kober et al.2011). An interlaboratory comparison study was set up in the framework of the CRONUS-EU initiative (Stuart and Dunai2009) with the aim to address this issue and provide the noble gas community with a well-characterized reference standard for the analysis of cosmogenic 21Ne in quartz. The CREU-1 standard is a mixture of natural quartz pebbles, rich in cosmogenic 21Ne, which were crushed and thoroughly homogenized to ensure optimal reproducibility (Section 2). Two size-fractions of CREU-1 were analyzed by five prominent cosmogenic noble gas laboratories, each of which used different experimental setups and data reduction protocols (Section 3). In total, 50 aliquots of CREU-1 were analyzed, with reported analytical precisions of 2-6%, but an external reproducibility of 7.1% (Section 4). These analyses were tied to a further 10 measurements of CRONUS-A, which is a second reference material prepared from an Antarctic quartzite analysed by two of the participating labs (Section 5).

2. Standard material

The CREU-1 standard material is pure quartz prepared from exposed vein-quartz clasts of a Miocene erosion surface (1933’53.4”S, 707’1.5”W, 930 m) in the Atacama desert near Pisagua, Chile (between sites B and C of Dunai et al.2005). The clasts were shed onto the surface from local sources after the main sedimentation episode at ~12-14 Ma (Dunai et al.2005). Approximately 400g of material was mixed from five clasts (sample name CH04/5, pebbles 5, 6, 7, 8 and 13, weighing 81g, 104g, 77g, 109g and 55g respectively) that had 21Ne excess concentrations within 5% of their mean value. After crushing in a W-carbide disk mill, five size fractions were prepared using stainless steel sieves:

  • 40-62μm: 10.45g, wet sieved and dried overnight at 50C
  • 63-125μm: 40.62g, wet sieved and dried overnight at 50C
  • 125-250μm: 74.4g, dry sieved
  • 250-500μm: 188g, dry sieved
  • >500μm: 15.8g, dry sieved

Of these five fractions, the 125-250μm and 250-500μm fractions were taken to produce the standard material, while the remaining fractions were preserved, but not processed any further. The 125-250μm and 250-500μm fractions were soaked in concentrated sulfuric acid at 120C overnight, to remove all iron coatings and accessory minerals such as rutile, sphene and fluorite. After the acid treatment, the material was rinsed ten times in cold de-ionized water, followed by five times one hour ultrasonic rinsing in de-ionized water at 80C. Next, the quartz was dried overnight at 110C. Although the preparation steps outlined above probably already ensured a thoroughly mixed quartz sand, a FRITSCH© rotary cone sample divider laborette 27 was used to split the material into 16 equal fractions. Different aliquots of CREU-1 have been analyzed by five noble gas laboratories, at BGC (Berkeley), CRPG (Nancy), ETH (Zürich), GFZ (Potsdam) and SUERC (Glasgow).

3. Analytical methods

The five participating laboratories employed a variety of noble gas mass spectrometers and analytical procedures for cosmogenic 21Ne analysis. Rather than forcing all the participants to use the same heating schedules, gettering times and so forth, they were allowed to use their own measurement routines, so that the calibration exercise fully captured the diversity of approaches used for 21Ne analysis. The temperature steps and amount of material used are reported in Tables 1-3.

3.1. BGC

Neon extraction from quartz at BGC employed a 14-sample vacuum chamber with a 3 inch diameter sapphire viewport. Samples of up to 150 mg quartz were encapsulated in a Ta packet and heated through the viewport by a 150 W diode laser (λ = 810 nm) using a feedback control system in which the temperature of the packet was continuously monitored by an optical pyrometer coaxial with the laser delivery optic. Calibration of the pyrometer for the emissivity of the Ta packets was accomplished by placing a thermocouple in the same apparatus. Collateral heating of adjacent samples was prevented by completing one heating step for all samples before beginning the next heating step. This procedure was tested by interspersing blanks consisting of an empty Ta packet. After heating, sample gas was reacted with a SAES®; getter and adsorbed to a cryogenic trap at 20 K. Neon was then released into the mass spectrometer at 70 K. All sample heating, gas processing, and measurement operations were automatically controlled. Analyses were done with a MAP-215 mass spectrometer updated with modern ion-counting electronics. Under normal operating conditions, this machine had a relatively low Ar+/Ar++ ratio (2.5-5, depending on source tuning) and inadequate mass resolution to fully resolve 20Ne+ from 40Ar++, so a correction for background 40Ar++ was required. As described in Balco and Shuster (2009b), this was accomplished by introducing a 39Ar spike and monitoring the Ar charge ratio as well as the 40Ar+ signal throughout each analysis. The resulting correction on mass 20 varied between analyses, but was typically equivalent to 5.00±0.02×108 atoms 20Ne. Similarly, a correction for 12C16O 2++ on mass 22 was made by establishing a relationship between the Ar and CO2 charge ratios. Absolute calibration of Ne abundance was made by peak height comparison against an air standard processed in the same way as the samples and analyzed several times daily. Linearity of machine response was verified by varying the volume of the air standard. The pressure of the air standard reservoir was measured during loading with an MKS Baratron manometer, and corrected for atmospheric water vapor using three separate hygrometers at the time of air sample collection. Absolute volumes of the reservoir and pipette were determined by differential pressure measurements, again using the Baratron, against two separate reference glass ampules whose volumes were independently measured by Hg weighing. The amount of cosmogenic 21Ne was calculated by assuming two-component mixing of atmospheric and cosmogenic neon. Reported uncertainties include i) counting uncertainties on all masses, including those used to generate corrections for 40Ar++ and CO 2++; ii) uncertainty in blank subtraction (the 21Ne process blank was ~ 0.5 Hz or ~ 90,000 atoms, which was < 1% of typical signals on mass 21 for these measurements); and iii) the reproducibility of the air standards (~1% for 20Ne, ~3% for 21Ne).

3.2. CRPG

After 10 minutes cleaning in an acetone ultrasonic bath, quartz aliquots were wrapped in copper foils (Alfa Aesar®;, 0.025 mm thick, 99.8%). Samples were then loaded under high vacuum in a stainless steel carousel that had been baked during 10 h at 80C. Gas extraction from the quartz was realized by 25 minutes heating in a home-designed single vacuum resistance furnace with a boron nitride crucible (Zimmermann et al.in press). Sequential purification with charcoals in liquid nitrogen, titanium sponges (JohnsonMatthey®;, mesh m3N8 t2N8) and SAES®; getters (ST172/HI/20-10/650C) permitted gas cleaning by removal of H2O, Ar, Kr, Xe and hydrocarbons. Ne was not separated from He. The purified gas was finally analyzed using a VG5400 mass spectrometer. Corrections for isobaric interferences of 40Ar++ at m/e = 20 and 12C16O 2++ at m/e = 22 were negligible compared to the amount of analyzed neon. The mass spectrometer sensitivity was determined by peak height comparison against a 0.2 cm3 (~1.6×1010 atoms of 20Ne) pipette of a gas standard having an atmospheric composition. Typical furnace blanks at 1000-1300C (25 min) were 1.0±0.2×108, 3±1×105 and 1.63±0.06×107 atoms of 20Ne, 21Ne and 22Ne, respectively. Excess 21Ne (21Ne*) concentrations were calculated following:

21   *       20
  N e =  Rc ×   N em × (Rm  - Ra)∕(Rc -  Ra)
(1)

where 20Ne m is the measured 20Ne, R c is the cosmogenic 21Ne/20Ne-ratio (R c = 0.8; Niedermann2002), Rm is the measured 21Ne/20Ne-ratio, and R a is the atmospheric 21Ne/20Ne-ratio (R a = 0.00296).

3.3. ETH

Noble gases were extracted by heating in a molybdenum crucible. Released gases were cleaned in a stainless steel extraction line equipped with Al/Zr-getters (SAES®;) and activated charcoal held at the temperature of liquid nitrogen before He and Ne were expanded to a cryogenic pump. Helium and neon were separated by adsorbing neon at 14 K on stainless steel frits and analyzing helium first. After pumping away the helium, neon was released from the cryotrap at 50 K. Noble gas analyses were performed in a custom-made, all-metal magnetic sector mass-spectrometer (90, 210 mm radius) equipped with a modified Baur-Signer ion source with essentially constant sensitivity over the pressure range relevant for this work (Baur1980). The ion source was equipped with a compressor device increasing the sensitivity by factors of 120 and 200 for 3He and 21Ne, respectively (Baur1999) compared to the sensitivities of the same spectrometer with the compressor turned off. The absolute sensitivity and mass discrimination of the mass spectrometer were determined by analysing known amounts of standard noble gas mixtures prepared from commercially available pure gases. The Ne isotopic composition of the standard gas was cross calibrated against two air standards (Heber et al.2009). Similarly, the Ne amounts delivered by the standard pipette were cross calibrated with air standards as well as with other independently filled standard gas bottles. The uncertainty of the Ne standard gas amounts is estimated to be 2% (Heber et al.2009). Full procedural blanks (45’ at 600oC + 20’ at 800 + 15’ at 1750C) were 1.211±0.006×108, 3.5±0.2×105, and 1.17±0.01×107 atoms of 20Ne, 21Ne and 22Ne, respectively. Corrections for isobaric interferences on mass 20 have been applied for 40Ar++ and H218O+ but were always less than 2%. No correction for CO 2++ on 22Ne was necessary. The low correction factors for doubly charged species were the results of a low electron acceleration voltage of 45V in the ion source. Excess 21Ne (21Ne*) concentrations were calculated with Equation 1.

3.4. GFZ

CREU-1 quartz samples were wrapped in aluminium foil and loaded in a sample carousel without further treatment, except for two aliquots of the 250-500μm fraction (GFZ-6-7) which were crushed to ~50μm grain size in an agate mortar before loading. Noble gases were extracted in a resistance-heated furnace equipped with a tantalum crucible and molybdenum liner and analyzed in either of two VG5400 noble gas mass spectrometers, with measurements GFZ1-7 being measured on one machine, and GFZ8-11 on the other (Tables 1 and 2). GFZ-8 was not heated, but instead crushed in vacuo between two hard metal jaws in order to test whether Ne trapped in fluid inclusions of CREU-1 has an atmospheric isotopic composition. Gas purification involved a dry ice trap, two titanium sponge and foil getters, and two SAES®; (Zr-Al) getters. The noble gases were trapped on stainless steel frits and/or activated charcoal in cryogenic adsorbers and sequentially released for He, Ne, and Ar-Kr-Xe analysis. Isobaric interferences of 40Ar++ at m/e=20 (up to 20% at 400C) and 12C16O 2++ at m/e=22 (up to 10% at 400C) were corrected according to the method described by Niedermann et al. (19931997). A correction for H218O+ at m/e=20 was not necessary due to the mass resolution of 600. Blanks had an atmospheric composition and contained 1-3×107 atoms of 20Ne, depending on temperature. Excess 21Ne was calculated without applying a blank correction, assuming an atmospheric origin of all the measured 20Ne:

21N e * =21 N em × (Rm  - Ra )∕Rm
(2)

with all abbreviations as in Equation 1. In some cases a high atmospheric Ne memory (i.e., rapid decay of non-atmospheric Ne isotope ratios) required the application of a special procedure to derive the Ne concentration and isotopic composition at the time of gas admission to the mass spectrometer (see Goethals et al.2009). Absolute noble gas concentrations were obtained by peak height comparison against a 0.1 cm3 pipette of calibration gas (an artificial mixture of the five noble gases in nitrogen provided by Linde company; Niedermann et al.1997), which was cross-calibrated in the 1990s against glass ampoule noble gas standards made available by O. Eugster (University of Bern) and whose noble gas concentrations are judged accurate to ~3% at 95% confidence level, and have been propagated into the overall uncertainty.

3.5. SUERC

The clean quartz was thoroughly rinsed in ultra-pure acetone and packed into aluminium foil cylinders. Cosmogenic Ne was extracted by heating each sample packet for 20 minutes. The active gases were removed by exposure to two hot SAES®; (Zr-Al) getters during heating, and for a further 20 minutes as the furnace cooled. The heavy noble gases and residual active gases were subsequently adsorbed on liquid nitrogen cooled activated charcoal for 10 minutes and exposed to a getter at room temperature to adsorb hydrogen. Neon was then adsorbed on activated charcoal in a cryostatic cold head at 30K. The helium was pumped for 1 minute, then the Ne was desorbed from the charcoal trap at 100K. Neon isotopes were analyzed statically in a MAP-215 magnetic sector mass spectrometer equipped with a modified Nier-type ion source, an axial electron multiplier (Burle Channeltron) operated in pulse-counting mode and a Faraday detector. A room temperature SAES®; G50 getter and a liquid nitrogen-cooled activated charcoal trap were used to minimize the contribution of interfering species during analysis. The data presented here were taken over a period of two years. Consequently source conditions changed to a small degree. Typically the source was tuned for Ne sensitivity prior to analytical periods; electron voltage of 88 V, trap current of 500 μA and an acceleration voltage of 3 kV. A slit in front of the electron multiplier was used to achieve a resolving power (m/Δm) of approximately 400. For all samples and calibrations the abundances of masses 18, 19, 20, 21, 22, 40 and 44 were determined by integrating counts recorded in 40-100 blocks of 5 seconds each. Peak heights of masses 2 and 16 were measured on the Faraday detector. Instrumental sensitivity was calculated from repeated analysis of aliquots of 2.2×1010 atoms 20Ne in air sampled from a 5 liter reservoir. Isotopic mass discrimination was approximately 0.50 ± 0.03 %/amu. The average high temperature 20Ne blank was 1×108 atoms. There was no observed increase when empty Al foil was heated. The Ne isotopic composition of blank measurements after correction for interfering species (see below) was indistinguishable from air ratios. Since it is likely that a significant amount of air-derived Ne is released from the quartz during heating, no blank correction has been made to the data. Excess 21Ne concentrations were calculated assuming an atmospheric origin of all the measured 20Ne according to Equation 2. Interference at m/e = 20 from H218O+ was calculated from measurement of H216O+ at mass 18. The contribution never exceeded 0.03%. No H19F+ signal was observed in blanks and mass spectrometer backgrounds. The dominant interference at m/e = 20 came from 40Ar++. The charge state ratio 40Ar+/40Ar++ is governed by the partial pressure of H in the mass spectrometer ionization region. A first-order relationship between 40Ar+/40Ar++ and H+ beam size was recorded. The partial pressure of H remained constant resulting in 40Ar+/40Ar++ = 2.30-2.32. The contribution of 40Ar++ to the measured 20Ne signal in CREU quartz samples was <1%. Correction for 12C16O 2++ at m/e = 22 was calculated from measured mass 44 (12C16O 2+) using a CO2+/CO 2++ = 50 to 58 (determined by repeated measurements interspersed with sample measurements). No pressure dependence on the CO2+/CO 2++ ratio was recorded for a 50-fold variation in the partial pressure of H and CO2. Correction for interfering 12C16O 2++ never exceeded 1%.

4. Results


PIC

Figure 1: Neon three-isotope plots of (a) all the individual heating steps and (b) the total released gas for each analyzed CREU-1 aliquot. The data fit a spallation line with a slope of 1.108 ± 0.014 (2σ, MSWD = 3.4). Error symbols are 1σ.


All five labs reported data for the coarse fraction, while three labs measured the fine fraction as well. The results for both sets of analyses are reported in Tables 1 and 2. The 21Ne/20Ne and 22Ne/20Ne compositions of the individual heating steps and their sums are consistent with a predominantly spallogenic origin of the released 21Ne (Figure 1). The pooled analyses comprise 16-30% atmospheric and 70-84% excess 21Ne, with individual heating steps containing up to 98% excess 21Ne. Linear regression of the spallation line yields a slope of 1.108 ± 0.014 (2σ), which is in statistical agreement with previously published values (Table 8 of Niedermann2002).

The total excess 21Ne contents of all the aliquots are shown in Figure 2. The reported 2σ analytical uncertainties are between 2 and 6%. The MSWD (Mean Square of the Weighted Deviates, a.k.a. ‘reduced Chi-square’, McIntyre et al.1966) is reasonably close to unity for ETH, GFZ and BGC, indicating good agreement of the observed scatter with the measurement errors. The extremely low MSWD of 0.005 for CRPG may indicate overestimated analytical uncertainties, but could also be due to chance, as only two aliquots were analyzed. Finally, the coarse fraction of SUERC is characterized by an MSWD of 4.1, which may indicate underestimated analytical uncertainties. However, measurements of the fine fraction by the same lab have an MSWD of 1.3. There is no systematic difference between the fine and the coarse grain size fractions of ETH, SUERC and GFZ. Measurements GFZ-6-7 were performed on material from the coarse fraction that was crushed for 5 minutes in air with an agate mortar to ~50μm, resulting in some loss of excess 21Ne. GFZ-8 was crushed in vacuo, and the data shown are for that crushing extraction. The 21Ne excess of GFZ-8 is considerably less than the 21Ne deficit in GFZ-6-7, probably because the in-vacuo crusher was much less efficient than the mortar. Measurements GFZ-6-8 were not included in subsequent calculations and figures. Total 3He concentrations measured at GFZ were 204±10×106at/g for the coarse fraction (seven measurements) and 109±7×106at/g (a single analysis) for the fine fraction. The resulting 21Ne∕3He-ratios are significantly greater than the production-rate ratio. This is likely caused by a combination of helium loss due to hot acid etching during sample preparation, and the fact that helium is not quantitatively retained in quartz at surface temperatures (Shuster and Farley2005).

BGC analyzed material from two different vials of CREU-1, thus presenting an opportunity to verify the homogeneity of the standard. Measurements BGC-1-4 were performed on material from the same vial as ETH, whereas measurements BGC-5-8 were done on the same vial as CRPG. The observed difference between the two vials analyzed by BGC falls within the analytical uncertainty. The difference between the results of BGC and ETH/CRPG, however, falls well outside the statistically acceptable range. The error-weighted means of all the labs do not agree with each other within the analytical uncertainties, defined as the standard errors of those means. Therefore, in order to calculate a global average of all the data (using both the fine and the coarse grain fractions), we used a random effects model with two sources of uncertainty. We assume that the intra-laboratory averages xi (where i is an identifier for each participating lab) come from a normal distribution of the form:

xi ~ N (μ,σ2i + ζ2)
(3)

where μ is the global mean, σi2 the analytical uncertainty (variance) of the ith lab, and ζ2 is the amount of overdispersion, i.e. the excess scatter (variance) that cannot be explained by the analytical uncertainty alone. To understand this formula, consider the following two special cases. If σi = 0 (perfect reproducibility within each lab) then μ is the arithmetic mean of the laboratory averages. And if ζ = 0 (perfect reproducibility between all labs) then μ is the error-weighted mean of those same laboratory averages. In order to simultaneously take into account the finite analytical precision of each lab and the variance between the labs, Equation 3 was iteratively solved for both μ and ζ, yielding an average 21Ne concentration of 348±10×106at/g and an overdispersion (defined as 2ζ/μ) of 7.1%.


PIC

Figure 2: Overview of all the reported 21Ne concentrations and 2σ uncertainties, with indication of the error-weighted means for each participating laboratory. White bars are considered outliers and were not used to calculate the averages. Left and middle panels: coarse and fine fractions of CREU-1; right panel: CRONUS-A. Gray band marks the average and 2σ uncertainty of CREU-1.


5. Comparison with CRONUS-A

In addition to CREU-1, two of the participating labs also analyzed CRONUS-A as a second reference material. CRONUS-A was collected in Antarctica’s Arena Valley (77 52’ 58.9”S, 160 56’ 35.1”E, 1666m elevation), from a large (40kg) yet thin (~2cm) slab of Beacon sandstone. Quartz was purified at the University of Vermont by crushing, sieving and repeated etching in dilute HF, using procedures designed for cosmogenic 10Be-26Al analysis. CRPG reported one and BGC a further nine analyses of CRONUS-A, using the same protocols that were used for the CREU-1 measurements (Table 3). The average cosmogenic 21Ne content of the nine CRONUS-A samples measured by BGC was 338.9±3.8×106at/g, i.e. 7.6±3.7% lower than that of CREU-1. The single CRONUS-A analysis of CRPG is lower than its CREU-1 measurements by a similar amount (4±17%), although the reported analytical precision of the latter estimate is much poorer. Additionally, published CRONUS-A values have been reported by two laboratories which did not participate in the interlaboratory comparison, at Harvard University (330±3×106at/g, Middleton et al.2012) and the California Institute of Technology (338±10×106at/g, Amidon and Farley2012). Normalizing the average CRONUS-A value reported by BGC to the CREU-1 reference value results in a 21Ne concentration of 320±11×106at/g. We propose that when this value is used as a reference, CRONUS-A can serve as an alternative to CREU-1.

6. Discussion

It is interesting to note that significant amounts of excess 21Ne remained trapped in the quartz after the second highest heating step, at temperatures of up to 820C. Total degassing was not achieved until the final temperature step at 1140C and more. This is significantly higher than the 800C release temperature for cosmogenic neon reported by Niedermann (2002). Nevertheless, for all samples of all labs, the data points of the higher temperature steps plot on the mixing line between atmospheric and cosmogenic neon (Figure 1), which strongly suggests that the non-atmospheric neon in all samples is essentially purely cosmogenic, although quartz occasionally also contains a nucleogenic neon component released at high temperature with a 21Ne/22Ne ratio of approximately unity (NeHT Niedermann et al.1994Niedermann2002). However, in view of the position of all data points in Figure 1 it seems very improbable that a sizeable fraction of the non-atmospheric 21Ne in our samples could be nucleogenic NeHT . Even in this unlikely case this would be largely irrelevant for the purpose of interlaboratory comparison, because for all samples we sum the non-atmospheric 21Ne from all temperature steps.

Despite the fact that CREU-1 is pure and highly enriched in spallogenic neon, the 21Ne concentrations reported by the participating labs are significantly overdispersed with respect to the formal analytical uncertainty. In theory, this overdispersion could be due to inhomogeneity of the standard material itself, as different labs analyzed aliquots from different vials of CREU-1. However, the analysis of two of these vials by BGC, and comparison with measurements of those same vials by ETH and CRPG, shows that this is not the case. Therefore, CREU-1 is homogenous. If the overdispersion cannot be attributed to the standard material itself, then it must be due to biases introduced by the different standard calibration bottles used (Heber et al.2009), or to differences in the neon sensitivity between samples and standards introduced by sample processing or tuning conditions.

7. Conclusion

Our calibration experiment has shown that, although the reported analytical precision of cosmogenic noble gas measurements may be as low as 2%, the accuracy is not quite as good. We suggest that the 7.1% dispersion observed in our study be used as a more realistic estimate of the accuracy of the 21Ne method at the present time. It should be borne in mind that this may even be an optimistic value, for a highly enriched and well behaved standard material. Using realistic and conservative analytical uncertainties is especially important for studies combining 21Ne with other (radio)nuclides, and to assess the resolving power of such studies. For single nuclide studies, CREU-1 or CRONUS-A measurements can be used to normalize 21Ne to the reference values reported in this paper, so that measurements from different labs can be compared on an equal footing and relative differences in 21Ne can be compared on the level of the analytical precision (Dunai and Stuart2009). Those interested in obtaining aliquots of these standards may contact T. Dunai (tdunai@uni-koeln.de) for CREU-1 or T. Jull (jull@email.arizona.edu) for CRONUS-A.

Acknowledgments

This research was funded by CRONUS-EU (Marie Curie RTN project 511927). We would like to thank Mark Kurz and an anonymous reviewer for positive and constructive feedback on the submitted manuscript.


Table 1: Summary table of the coarse fraction of CREU-1. ‘21/20’ and ‘22/20’ are the 21Ne∕20Ne and 22Ne∕20Ne ratios, 21Ne* is excess 21Ne. Temperatures of ETH-4-7 (marked by an asterisk) and SUERC-8-9 (omitted) were set when the crucible was really full causing samples to be degassed at positions where the temperature was lower than the nominal crucible temperature. GFZ-6-7 were crushed to small grain size (~50μm) before loading, while GFZ-8 was degassed by in vacuo crushing instead of heating. These samples were not included in Figures 1 and 2.













mass T 20Ne 2σ 21/20 2σ 22/20 2σ 21Ne* 2σ sum 2σ
[mg] [C][×109at/g] [×10-3] [×10-3] [×106at/g] [×106at/g]













ETH-1 66.37 600 13.370.27 24.230.15 126.610.83 285.5 6.1













800 9.850.20 6.960.11 105.590.67 39.6 1.0
1750 8.700.19 6.220.10 106.571.00 28.4 0.8 353.5 7.6
ETH-2 66 600 11.450.24 27.780.26 129.170.74 285.1 6.5
800 8.420.17 7.520.13 106.921.30 38.6 1.0
1750 6.620.18 6.700.14 105.900.98 24.9 0.8 348.5 7.7
ETH-3 49.23 800 22.130.45 17.150.13 117.110.29 315.3 6.9
1750 11.160.22 5.990.12 104.811.15 34.0 0.9 349.2 7.6
ETH-4 82.13 600* 1.340.03 74.660.76 180.3 1.9 96.3 2.4
800* 1.390.15 26.910.56 125.9 3.1 33.4 3.6
1750* 25.180.50 11.160.11 109.0 1.2 207.3 4.7 335.8 7.3
ETH-5 48.56 600* 1.890.04 47.5 1.0 155.0 2.7 84.7 2.6
800* 1.580.28 22.440.42 119.2 2.3 30.9 5.5
1750* 29.420.59 10.740.17 108.4 1.8 229.8 5.8 344.2 8.3
ETH-6 48.73 600* 1.430.03 66.6 1.0 171.9 4.2 91.5 2.6
800* 2.210.15 24.180.55 124.2 1.7 47.1 3.3
1750* 28.540.57 10.310.10 108.6 1.0 210.5 4.7 347.8 7.6
ETH-7 81.3 600* 1.390.03 75.2 1.5 187.1 3.5 100.8 3.0
800* 1.750.13 27.370.42 126.4 2.1 42.8 3.3
1750* 28.560.57 10.240.09 108.5 0.7 208.8 4.5 351.2 7.7
BGC-1 109.9 370 4.1 1.5 24.4 9.1 117 49 86.9 6.0
740 14.0 0.9 17.3 1.2 119 11 200.0 7.7
1140 20.3 1.3 7.0 0.5 105 8 81.5 7.4 368 12
BGC-2 129.4 370 5.3 0.8 22.4 3.6 125 25 101.6 5.2
740 15.0 0.5 15.9 0.6 114 7 194.9 7.4
1140 19.0 0.6 6.8 0.3 105 5 71.1 4.2 368 10
BGC-3 115.7 390 4.7 0.6 28.2 3.6 132 32 121.3 6.9
780 19.8 1.6 13.5 0.7 114 7 210 12
1140 12.7 0.5 6.7 0.4 109 6 46.5 4.1 378 15
BGC-4 104.8 390 5.8 0.7 23.4 3.0 124 26 119.6 8.1
780 16.5 0.9 15.4 0.9 113 8 203.4 8.6
1140 14.9 1.0 6.6 0.5 105 10 50.1 5.7 373 13
BGC-5 103.5 370 3.1 0.6 32.3 6.0 137 67 90.7 6.9
740 14.2 0.7 17.1 0.8 121 12 203.4 9.6
1140 22.7 1.2 6.5 0.3 108 5 79.6 7.5 374 14
BGC-6 83.2 370 3.5 1.7 28 14 141 78 87.2 6.9
740 12.9 0.8 18.3 1.2 122 11 201.9 7.8
1140 18.2 1.0 7.3 0.5 109 7 80.3 6.5 369 12
BGC-7 75 390 6.3 2.0 17.8 5.7 115 43 92.4 8.6
780 15.3 1.0 16.2 1.1 118 10 199 13
1140 17.6 0.9 6.6 0.4 110 8 61.8 6.1 353 17
BGC-8 144.4 390 5.9 1.0 20.7 3.3 121 24 103.4 5.6
780 15.3 0.5 16.0 0.6 115 5 196.3 4.4
1140 17.2 0.6 6.5 0.2 106 4 59.7 3.5 359 8
BGC-9 110.2 370 3.9 0.6 28.9 4.8 139 42 101.7 5.8
740 14.0 0.8 17.2 1.0 117 10 202 11
1140 18.9 0.6 6.6 0.3 108 5 65.5 5.3 369 13
SUERC-1 165.1 1350 30.2 1.8 14.500.38 111.5 3.5 343 23 343 23
SUERC-2 239.6 480 1.820.11 24.120.80 127.8 5.8 37.9 2.7
550 1.020.06 40.3 1.6 143.3 6.4 37.8 3.0
650 3.670.22 20.210.61 118.2 5.2 62.4 4.4
800 8.550.51 15.880.45 116.0 4.8 108.9 7.4
1400 18.0 1.1 6.020.14 105.2 3.9 54.2 3.5 301 10
SUERC-3 258.4 480 4.910.29 25.100.67 125.4 5.5 107.3 7.2
550 5.310.32 25.200.71 127.6 5.1 116.5 7.9
650 2.790.17 16.400.63 113.5 5.4 37.0 2.9
800 11.330.68 8.340.27 106.9 4.4 60.1 4.3
1200 14.730.88 6.200.17 104.7 4.4 47.0 3.2
1350 5.240.31 5.240.19 103.7 4.8 11.8 0.9 380 12
SUERC-4 203 1350 36.5 2.2 12.570.32 112.8 2.3 350 23 350 23
SUERC-5 204.6 1350 33.8 2.0 13.320.30 114.4 2.3 349 23 349 23
SUERC-6 261.3 1350 32.9 2.0 13.060.33 115.1 2.3 332 22 332 22
SUERC-7 160.2 1350 40.6 2.4 11.590.23 110.8 1.6 350 23 350 23
SUERC-8 237.3 - 17.7 1.1 19.330.33 119.5 1.5 289 18
- 4.950.30 6.270.26 106.3 3.8 16.3 1.2
- 3.170.19 5.120.16 104.1 1.3 6.8 0.5
- 2.530.15 4.860.17 106.6 1.8 4.8 0.4 317 18
SUERC-9 196.1 - 21.1 1.3 16.720.28 117.2 0.9 289 18
- 6.530.39 5.750.10 105.0 0.8 18.1 1.1
- 6.160.37 5.140.14 105.4 1.3 13.4 0.9
- 3.450.21 4.400.20 105.8 1.5 4.9 0.4 325 18
SUERC-10 175.4 1350 35.4 2.1 12.790.27 113.1 0.9 346 23
1350 0.810.05 2.600.85 105.8 3.8 - - 346 23
SUERC-11 209.2 1350 11.880.71 25.000.44 126.3 1.1 260 16
1350 23.8 1.4 6.700.13 106.5 0.8 88.7 5.7 349 17
SUERC-12 187.8 1350 32.5 1.9 13.330.21 113.0 0.8 335 21
1350 0.880.05 4.450.43 105.5 3.4 1.3 0.2 336 21
SUERC-13 202.3 1350 33.5 2.0 13.590.22 114.5 0.8 355 22 355 22
SUERC-14 92.8 1350 39.7 2.4 12.220.12 111.5 0.5 357 22 357 22
GFZ-1 50.56 400 0.70.13 149 25 277 34 99.4 9.5
800 18.6 1.3 15.040.47 115.6 2.1 224 17
1200 5.70.47 6.700.41 107.4 4.5 21.2 2.6 345 20
GFZ-2 102.5 400 1.10.14 95.2 9.6 205 13 105.1 8.3
800 16.2 1.0 15.880.34 115.8 1.2 209 13
1200 11.30.73 6.320.12 105.9 1.4 38.1 2.6 352 16
GFZ-3 99.58 400 0.40.14 181 54 300 61 72.8 6.7
800 15.2 1.1 18.720.95 117.1 1.2 239 22
1200 11.80.84 6.330.19 103.4 1.2 39.6 3.5 351 23
GFZ-4 100.4 400 0.70.08 113 10 220 14 76.2 4.9
800 21.1 1.1 13.820.14 113.1 3.3 229 12
1200 10.80.58 6.220.20 104.5 0.8 35.2 2.9 340 13
GFZ-5104.52 400 0.50.09 174 31 282 33 79.8 7.2
800 16.7 1.0 16.770.30 115.3 0.8 231 14
1200 9.40.56 6.420.18 107.5 1.3 32.5 2.5 343 16
GFZ-6101.26 400 1.10.10 216 14 335 13 239 19
800 2.60.20 25.210.93 126.2 4.8 57.0 4.0
1200 0.10.06 4.204.40 90 35 0.1 0.2 296 19
GFZ-7110.74 400 0.90.13 229 29 349 33 194 14
800 2.10.21 28.9 2.0 130.5 3.4 54.6 3.8
1200 0.00.10 - - - - 0.1 0.2 249 15
GFZ-8 502.7 20 10.20.53 3.960.07 102.9 0.9 10.160.89 10.160.89
GFZ-9201.19 400 1.10.11 97.7 5.9 207.3 7.2 99.1 7.5
600 3.70.28 32.690.61 133.7 2.7 111.2 8.1
800 12.70.92 10.000.16 110.3 0.7 89.2 6.7
1200 10.1 1.1 6.210.16 105.7 1.6 32.8 3.8 332 13
GFZ-10 201.1 400 2.40.19 57.3 3.1 162.8 3.4 129.2 7.3
600 5.50.32 24.600.69 126.5 1.4 118.2 6.7
800 10.30.61 7.790.15 106.6 0.6 50.0 3.2
1200 6.60.41 6.610.20 108.5 0.9 24.0 1.7 321 11
CRPG-1 149.1 820 14.150.54 21.7 1.2 122.4 6.5 267 19
1260 16.950.64 6.220.45 104.5 5.4 55.4 7.9 322 21
CRPG-2 83.3 1180 28.0 1.1 14.460.79 114.3 6.0 323 25
1260 0.030.16 12 79 - - 0.3 2.6 323 26















Table 2: Same as Table 1, but for the fine fraction of CREU-1.













mass T 20Ne 2σ 21/20 2σ 22/20 2σ 21Ne* 2σ sum2σ
[mg] [C][×109at/g] [×10-3] [×10-3] [×106at/g] [×106at/g]













ETH-8 73.38 600 8.590.20 40.610.34 147.1 1.5 323.4 7.9













800 6.780.21 6.210.18 107.6 1.3 22.1 0.9
1750 8.390.22 5.200.11 108.3 1.2 18.8 0.6 36412
ETH-9 70.42 600 9.910.21 35.200.31 142.6 1.1 319.3 7.4
800 6.240.14 6.000.08 106.8 1.2 19.0 0.5
1750 8.110.16 5.130.15 105.7 1.5 17.7 0.6 35611
SUERC-15 378.6 450 3.170.19 58.9 1.5 164.5 5.3 174.611.4
550 3.840.23 28.390.75 130.9 4.2 96.4 6.4
650 3.880.23 10.030.32 107.4 3.8 27.0 1.9
800 10.410.62 6.050.17 107.2 3.4 31.8 2.2
1350 9.360.56 4.620.16 104.2 3.3 15.3 1.1 34513
SUERC-16 237.1 550 0.860.05 52.0 1.9 158.5 8.8 41.8 3.2
650 2.340.14 54.8 1.6 161.7 7.1 119.9 8.3
800 2.110.13 26.900.62 121.5 4.5 49.9 3.2
1350 24.5 1.5 8.850.20 108.0 4.0 142.2 9.1
480 0.730.04 10.140.81 106.4 5.6 5.2 0.6 35913
SUERC-17 331350 31.6 1.9 15.040.29 115.9 1.3 374.924.5 37525
SUERC-18 156.21350 26.7 1.6 16.560.27 116.2 1.2 357.722.5 35822
SUERC-19 1661350 25.6 1.5 16.990.17 116.8 0.9 352.721.6 35322
GFZ-11100.24 400 0.960.13 152.016.0 263.018.0 142.913.7
600 3.300.27 38.3 1.2 139.6 3.1 116.6 8.6
800 12.570.94 7.870.16 108.2 1.0 61.8 4.8
1200 3.330.30 5.660.34 104.3 4.1 9.0 1.2 33017















Table 3: Same as Tables 1 and 2, but for CRONUS-A.













mass T 20Ne 2σ 21/20 2σ 22/202σ 21Ne* 2σ sum2σ
[mg] [C][×109at/g] [×10-3] [×10-3] [×106at/g] [×106at/g]













BGC-10137.7 390 2.210.72 69.2822.43 17259 144.4 8.2













780 11.481.10 19.13 1.63 11911 183.011.3
1140 1.650.35 9.53 2.16 11036 10.6 1.8 33814
BGC-11105.9 390 2.190.52 72.4816.99 17949 150.2 8.1
780 9.720.77 21.06 1.67 12612 173.2 8.4
1140 1.330.65 11.39 5.70 10169 11.1 2.6 33412
BGC-12122.5 370 2.681.24 50.5923.39 14072 124.7 6.5
740 9.860.33 22.09 0.76 124 8 192.0 4.5
1140 2.110.63 12.35 3.72 12243 20.0 2.5 337 8
BGC-13107.5 370 4.261.59 33.4412.55 12450 126.3 6.9
740 11.271.11 20.37 2.07 11815 195.1 7.7
1140 3.701.79 6.44 3.13 9754 12.8 5.6 33412
BGC-14 66.7 370 4.180.89 32.74 7.03 13846 123.7 5.5
740 10.430.46 21.66 1.10 12316 195.2 8.5
1140 2.130.89 11.16 4.81 12877 17.2 3.6 33611
BGC-15138.3 370 3.200.73 42.97 9.58 13941 127.110.0
740 10.240.42 22.06 0.68 123 8 193.1 8.7
1140 2.330.57 9.84 2.47 12639 15.7 2.3 33613
BGC-16167.8 390 4.530.63 33.39 4.36 13521 138.3 8.1
780 12.320.68 19.17 0.79 118 7 199.310.2
1140 1.930.52 10.52 2.91 12246 14.0 2.1 35213
BGC-17138.1 370 2.570.61 42.6810.09 15550 102.0 5.4
740 11.970.64 21.61 0.99 123 9 226.312.9
1140 2.960.56 10.70 2.10 11032 22.3 2.8 35114
BGC-18144.8 370 2.470.39 42.38 6.57 15356 97.0 8.0
740 11.290.45 22.40 0.93 12611 219.1 7.2
1140 3.930.50 9.24 1.26 11426 24.8 2.7 34111
CRPG-3 37.21200 13.210.64 26.43 1.75 138 9 311.227.6
1280 0.320.35 2.33 7.29 < DL 0 < DL 2.3 31128














References

   Amidon, W. H., Farley, K. A., 2012. Cosmogenic 3He and 21Ne dating of biotite and hornblende. Earth and Planetary Science Letters 313-314 (0), 86 – 94.

   Balco, G., Shuster, D. L., 2009a. 26Al-10Be-21Ne burial dating. Earth and Planetary Science Letters 286, 570–575.

   Balco, G., Shuster, D. L., 2009b. Production rate of cosmogenic 21Ne in quartz estimated from 10Be, 26Al, and 21Ne concentrations in slowly eroding Antarctic bedrock surfaces. Earth and Planetary Science Letters 281, 48–58.

   Baur, H., 1980. Numerische Simulation und praktische Erprobung einer rotationssymmetrischen Ionenquelle für Gasmassenspektrometer. Ph.D. Thesis, ETH-Zürich No. 6596.

   Baur, H., 1999. A noble-gas mass spectrometer compressor source with two orders of magnitude improvement in sensitivity. Eos, Transactions of the American Geophysical Union 80, F1118.

   Codilean, A. T., Bishop, P., Stuart, F. M., Hoey, T. B., Fabel, D., Freeman, S. P. H. T., 2008. Single-grain cosmogenic 21Ne concentrations in fluvial sediments reveal spatially variable erosion rates. Geology 36, 159–162.

   Dunai, T. J., González López, G. A., Juez-Larré, J., apr 2005. Oligocene Miocene age of aridity in the Atacama Desert revealed by exposure dating of erosion-sensitive landforms. Geology 33, 321–324.

   Dunai, T. J., Stuart, F. M., 2009. Reporting of cosmogenic nuclide data for exposure age and erosion rate determinations. Quaternary Geochronology 4 (6), 437–440.

   Fujioka, T., Chappell, J., Honda, M., Yatsevich, I., Fifield, K., Fabel, D., 2005. Global cooling initiated stony deserts in central Australia 2-4 Ma, dated by cosmogenic 21Ne-10Be. Geology 33, 993.

   Goethals, M. M., Niedermann, S., Hetzel, R., Fenton, C. R., 2009. Determining the impact of faulting on the rate of erosion in a low-relief landscape: A case study using in situ produced 21Ne on active normal faults in the Bishop Tuff, California. Geomorphology 103, 401–413.

   Heber, V. S., Wieler, R., Baur, H., Olinger, C., Friedmann, T. A., Burnett, D. S., 2009. Noble gas composition of the solar wind as collected by the Genesis mission. Geochimica et Cosmochimica Acta 73, 7414–7432.

   Kober, F., Alfimov, V., Ivy-Ochs, S., Kubik, P. W., Wieler, R., 2011. The cosmogenic 21Ne production rate in quartz evaluated on a large set of existing 21Ne-10Be data. Earth and Planetary Science Letters 302, 163–171.

   Kober, F., Ivy-Ochs, S., Leya, I., Baur, H., Magna, T., Wieler, R., Kubik, P. W., 2005. In situ cosmogenic 10Be and 21Ne in sanidine and in situ cosmogenic 3He in Fe Ti-oxide minerals. Earth and Planetary Science Letters 236, 404–418.

   Kober, F., Ivy-Ochs, S., Zeilinger, G., Schlunegger, F., Kubik, P. W., Baur, H., Wieler, R., 2009. Complex multiple cosmogenic nuclide concentration and histories in the arid Rio Lluta catchment, northern Chile. Earth Surface Processes and Landforms 34 (3), 398–412.

   Lal, D., 1991. Cosmic ray labeling of erosion surfaces: in situ nuclide production rates and erosion models. Earth and Planetary Science Letters 104, 424–439.

   McIntyre, G. A., Brooks, C., Compston, W., Turek, A., 1966. The Statistical Assessment of Rb-Sr Isochrons. Journal of Geophysical Research 71, 5459–5468.

   Middleton, J., Ackert Jr., R., Mukhopadhyay, S., 2012. Pothole and channel system formation in the mcmurdo dry valleys of antarctica: New insights from cosmogenic nuclides. Earth and Planetary Science Letters 355-356, 341 – 350.

   Niedermann, S., 2002. Cosmic-Ray-Produced Noble Gases in Terrestrial Rocks: Dating Tools for Surface Processes. In: Porcelli, D., Ballentine, C. J., Wieler, R. (Eds.), Noble Gases in Geochemistry and Cosmochemistry. Vol. 47 of Reviews in Mineralogy and Geochemistry. Mineralogical Society of America, pp. 731–784.

   Niedermann, S., Bach, W., Erzinger, J., 1997. Noble gas evidence for a lower mantle component in MORBs from the southern East Pacific Rise: Decoupling of helium and neon isotope systematics. Geochimica et Cosmochimica Acta 61, 2697–2715.

   Niedermann, S., Graf, T., Kim, J., Kohl, C., Marti, K., Nishiizumi, K., 1994. Cosmic-ray-produced 21Ne in terrestrial quartz: the neon inventory of Sierra Nevada quartz separates. Earth and Planetary Science Letters 125 (1-4), 341–355.

   Niedermann, S., Graf, T., Marti, K., 1993. Mass spectrometric identification of cosmic-ray-produced neon in terrestrial rocks with multiple neon components. Earth and Planetary Science Letters 118, 65–73.

   Poreda, R. J., Cerling, T. E., 1992. Cosmogenic neon in recent lavas from the western United States. Geophysical Research Letters 19, 1863–1866.

   Schäfer, J. M., Ivy-Ochs, S., Wieler, R., Leya, I., Baur, H., Denton, G. H., Schlüchter, C., 1999. Cosmogenic noble gas studies in the oldest landscape on earth: surface exposure ages of the Dry Valleys, Antarctica. Earth and Planetary Science Letters 167, 215–226.

   Shuster, D., Farley, K., 2005. Diffusion kinetics of proton-induced 21Ne, 3He, and 4He in quartz. Geochimica et Cosmochimica Acta 69, 2349–2359.

   Stuart, F. M., Dunai, T. J., 2009. Editorial. Quaternary Geochronology 4 (6), 435 – 436, advances in Cosmogenic Isotope Research from CRONUS-EU.

   Vermeesch, P., Fenton, C. R., Kober, F., Wiggs, G. F. S., Bristow, C. S., Xu, S., 2010. Sand residence times of one million years in the Namib Sand Sea from cosmogenic nuclides. Nature Geoscience 3, 862–865.

   Zimmermann, L., Blard, P.-H., Burnard, P., Medynski, S., Pik, R., Puchol., N., in press. A new single vacuum furnace design for cosmogenic 3He dating. Geostandards and Geoanalytical Research.