Diffusion of helium in zircon and apatite

D. J. Cherniak, E. B. Watson, Jay Thomas

Research output: Contribution to journalArticle

67 Citations (Scopus)

Abstract

Diffusion of helium has been characterized in natural zircon and apatite. Polished slabs of zircon and apatite, oriented either normal or parallel to c were implanted with 100 keV 3He at a dose of 5 × 1015 3 He/cm2. Diffusion experiments on implanted zircon and apatite were run in Pt capsules in 1-atm furnaces. 3He distributions following experiments were measured with Nuclear Reaction Analysis using the reaction 3He(d,p)4He. For diffusion in zircon we obtain the following Arrhenius relations:D⊥ c = 2.3 × 10- 7 exp (- 146 ± 11 kJ mol- 1 / R T) m2 s- 1 .D∥ c = 1.7 × 10- 5 exp (- 148 ± 17 kJ mol- 1 / R T) m2 s- 1 . Although activation energies for diffusion normal and parallel to c are comparable, there is marked diffusional anisotropy, with diffusion parallel to c nearly 2 orders of magnitude faster than transport normal to c. These diffusivities bracket the range of values determined for He diffusion in zircon in bulk-release experiments, although the role of anisotropy could not be directly evaluated in those measurements. In apatite, the following Arrhenius relation was obtained over the temperature range of 148-449 °C for diffusion normal to c:D = 2.10 × 10- 6 exp (- 117 ± 6 kJ mol- 1 / R T) m2 s- 1 . In contrast to zircon, apatite shows little evidence of anisotropy. He diffusivities obtained in this study fall about an order of magnitude lower than diffusivities measured through bulk release of He through step-heating, and within an order of magnitude of determinations where ion implantation was used to introduce helium and He distributions measured with elastic recoil detection. Since the diffusion of He in zircon exhibits such pronounced anisotropy, helium diffusional loss and closure cannot be modeled with simple spherical geometries and the assumption of isotropic diffusion. A finite-element code (CYLMOD) has recently been created to simulate diffusion in cylindrical geometry with differing radial and axial diffusion coefficients. We present some applications of the code in evaluating helium lost from zircon grains as a function of grain size and length to diameter ratios, and consider the effects of "shape anisotropy", where diffusion is isotropic (as in the case of apatite) but shapes of crystal grains or fragments may depart significantly from spherical geometry.

Original languageEnglish (US)
Pages (from-to)155-166
Number of pages12
JournalChemical Geology
Volume268
Issue number1-2
DOIs
StatePublished - Oct 20 2009
Externally publishedYes

Fingerprint

Apatites
Helium
helium
apatite
zircon
anisotropy
Anisotropy
diffusivity
geometry
Geometry
experiment
activation energy
Nuclear reactions
Experiments
slab
grain size
Ion implantation

Keywords

  • Apatite
  • Diffusion
  • Helium
  • Nuclear reaction analysis
  • Thermochronology
  • Zircon

ASJC Scopus subject areas

  • Geochemistry and Petrology
  • Geology

Cite this

Diffusion of helium in zircon and apatite. / Cherniak, D. J.; Watson, E. B.; Thomas, Jay.

In: Chemical Geology, Vol. 268, No. 1-2, 20.10.2009, p. 155-166.

Research output: Contribution to journalArticle

Cherniak, D. J. ; Watson, E. B. ; Thomas, Jay. / Diffusion of helium in zircon and apatite. In: Chemical Geology. 2009 ; Vol. 268, No. 1-2. pp. 155-166.
@article{6169d569192c4e28b3c59115ef41f0cd,
title = "Diffusion of helium in zircon and apatite",
abstract = "Diffusion of helium has been characterized in natural zircon and apatite. Polished slabs of zircon and apatite, oriented either normal or parallel to c were implanted with 100 keV 3He at a dose of 5 × 1015 3 He/cm2. Diffusion experiments on implanted zircon and apatite were run in Pt capsules in 1-atm furnaces. 3He distributions following experiments were measured with Nuclear Reaction Analysis using the reaction 3He(d,p)4He. For diffusion in zircon we obtain the following Arrhenius relations:D⊥ c = 2.3 × 10- 7 exp (- 146 ± 11 kJ mol- 1 / R T) m2 s- 1 .D∥ c = 1.7 × 10- 5 exp (- 148 ± 17 kJ mol- 1 / R T) m2 s- 1 . Although activation energies for diffusion normal and parallel to c are comparable, there is marked diffusional anisotropy, with diffusion parallel to c nearly 2 orders of magnitude faster than transport normal to c. These diffusivities bracket the range of values determined for He diffusion in zircon in bulk-release experiments, although the role of anisotropy could not be directly evaluated in those measurements. In apatite, the following Arrhenius relation was obtained over the temperature range of 148-449 °C for diffusion normal to c:D = 2.10 × 10- 6 exp (- 117 ± 6 kJ mol- 1 / R T) m2 s- 1 . In contrast to zircon, apatite shows little evidence of anisotropy. He diffusivities obtained in this study fall about an order of magnitude lower than diffusivities measured through bulk release of He through step-heating, and within an order of magnitude of determinations where ion implantation was used to introduce helium and He distributions measured with elastic recoil detection. Since the diffusion of He in zircon exhibits such pronounced anisotropy, helium diffusional loss and closure cannot be modeled with simple spherical geometries and the assumption of isotropic diffusion. A finite-element code (CYLMOD) has recently been created to simulate diffusion in cylindrical geometry with differing radial and axial diffusion coefficients. We present some applications of the code in evaluating helium lost from zircon grains as a function of grain size and length to diameter ratios, and consider the effects of {"}shape anisotropy{"}, where diffusion is isotropic (as in the case of apatite) but shapes of crystal grains or fragments may depart significantly from spherical geometry.",
keywords = "Apatite, Diffusion, Helium, Nuclear reaction analysis, Thermochronology, Zircon",
author = "Cherniak, {D. J.} and Watson, {E. B.} and Jay Thomas",
year = "2009",
month = "10",
day = "20",
doi = "10.1016/j.chemgeo.2009.08.011",
language = "English (US)",
volume = "268",
pages = "155--166",
journal = "Chemical Geology",
issn = "0009-2541",
publisher = "Elsevier",
number = "1-2",

}

TY - JOUR

T1 - Diffusion of helium in zircon and apatite

AU - Cherniak, D. J.

AU - Watson, E. B.

AU - Thomas, Jay

PY - 2009/10/20

Y1 - 2009/10/20

N2 - Diffusion of helium has been characterized in natural zircon and apatite. Polished slabs of zircon and apatite, oriented either normal or parallel to c were implanted with 100 keV 3He at a dose of 5 × 1015 3 He/cm2. Diffusion experiments on implanted zircon and apatite were run in Pt capsules in 1-atm furnaces. 3He distributions following experiments were measured with Nuclear Reaction Analysis using the reaction 3He(d,p)4He. For diffusion in zircon we obtain the following Arrhenius relations:D⊥ c = 2.3 × 10- 7 exp (- 146 ± 11 kJ mol- 1 / R T) m2 s- 1 .D∥ c = 1.7 × 10- 5 exp (- 148 ± 17 kJ mol- 1 / R T) m2 s- 1 . Although activation energies for diffusion normal and parallel to c are comparable, there is marked diffusional anisotropy, with diffusion parallel to c nearly 2 orders of magnitude faster than transport normal to c. These diffusivities bracket the range of values determined for He diffusion in zircon in bulk-release experiments, although the role of anisotropy could not be directly evaluated in those measurements. In apatite, the following Arrhenius relation was obtained over the temperature range of 148-449 °C for diffusion normal to c:D = 2.10 × 10- 6 exp (- 117 ± 6 kJ mol- 1 / R T) m2 s- 1 . In contrast to zircon, apatite shows little evidence of anisotropy. He diffusivities obtained in this study fall about an order of magnitude lower than diffusivities measured through bulk release of He through step-heating, and within an order of magnitude of determinations where ion implantation was used to introduce helium and He distributions measured with elastic recoil detection. Since the diffusion of He in zircon exhibits such pronounced anisotropy, helium diffusional loss and closure cannot be modeled with simple spherical geometries and the assumption of isotropic diffusion. A finite-element code (CYLMOD) has recently been created to simulate diffusion in cylindrical geometry with differing radial and axial diffusion coefficients. We present some applications of the code in evaluating helium lost from zircon grains as a function of grain size and length to diameter ratios, and consider the effects of "shape anisotropy", where diffusion is isotropic (as in the case of apatite) but shapes of crystal grains or fragments may depart significantly from spherical geometry.

AB - Diffusion of helium has been characterized in natural zircon and apatite. Polished slabs of zircon and apatite, oriented either normal or parallel to c were implanted with 100 keV 3He at a dose of 5 × 1015 3 He/cm2. Diffusion experiments on implanted zircon and apatite were run in Pt capsules in 1-atm furnaces. 3He distributions following experiments were measured with Nuclear Reaction Analysis using the reaction 3He(d,p)4He. For diffusion in zircon we obtain the following Arrhenius relations:D⊥ c = 2.3 × 10- 7 exp (- 146 ± 11 kJ mol- 1 / R T) m2 s- 1 .D∥ c = 1.7 × 10- 5 exp (- 148 ± 17 kJ mol- 1 / R T) m2 s- 1 . Although activation energies for diffusion normal and parallel to c are comparable, there is marked diffusional anisotropy, with diffusion parallel to c nearly 2 orders of magnitude faster than transport normal to c. These diffusivities bracket the range of values determined for He diffusion in zircon in bulk-release experiments, although the role of anisotropy could not be directly evaluated in those measurements. In apatite, the following Arrhenius relation was obtained over the temperature range of 148-449 °C for diffusion normal to c:D = 2.10 × 10- 6 exp (- 117 ± 6 kJ mol- 1 / R T) m2 s- 1 . In contrast to zircon, apatite shows little evidence of anisotropy. He diffusivities obtained in this study fall about an order of magnitude lower than diffusivities measured through bulk release of He through step-heating, and within an order of magnitude of determinations where ion implantation was used to introduce helium and He distributions measured with elastic recoil detection. Since the diffusion of He in zircon exhibits such pronounced anisotropy, helium diffusional loss and closure cannot be modeled with simple spherical geometries and the assumption of isotropic diffusion. A finite-element code (CYLMOD) has recently been created to simulate diffusion in cylindrical geometry with differing radial and axial diffusion coefficients. We present some applications of the code in evaluating helium lost from zircon grains as a function of grain size and length to diameter ratios, and consider the effects of "shape anisotropy", where diffusion is isotropic (as in the case of apatite) but shapes of crystal grains or fragments may depart significantly from spherical geometry.

KW - Apatite

KW - Diffusion

KW - Helium

KW - Nuclear reaction analysis

KW - Thermochronology

KW - Zircon

UR - http://www.scopus.com/inward/record.url?scp=70349383203&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=70349383203&partnerID=8YFLogxK

U2 - 10.1016/j.chemgeo.2009.08.011

DO - 10.1016/j.chemgeo.2009.08.011

M3 - Article

VL - 268

SP - 155

EP - 166

JO - Chemical Geology

JF - Chemical Geology

SN - 0009-2541

IS - 1-2

ER -