Empirical constraints on closure temperatures from a single diffusion coefficient

Abstract

The elucidation of thermal histories by geochronological and isotopic means is based fundamentally on solid-state diffusion and the concept of closure temperatures. Because diffusion is thermally activated, an analytical solution of the closure temperature (Tc*) can only be obtained if the diffusion coefficient D of the diffusion process is measured at two or more different temperatures. If the diffusion coefficient is known at only one temperature, however, the true closure temperature (Tc*) cannot be calculated analytically because there exist an infinite number of possible (apparent) closure temperatures (T¯c) which can be generated by this single datum. By introducing further empirical constraints to limit the range of possible closure temperatures, however, mathematical analysis of a modified form of the closure temperature equation shows that it is possible to make both qualitative and quantitative estimates of Tc* given knowledge of only one diffusion coefficient DM measured at one temperature TM. Qualitative constraints of the true closure temperature Tc* are obtained from the shapes of curves on a graph of the apparent Tc (T¯c) vs. activation energy E, in which each curve is based on a single diffusion coefficient measurement DM at temperature TM. Using a realistic range of E, the concavity of the curve shows whether TM is less than, approximately equal to, or greater than Tc*. Quantitative estimates are obtained by considering two dimensionless parameters [lnERT^c vs. Tc*/TM] derived from these curves. When these parameters are plotted for known argon diffusion data and for a given diffusion size and cooling rate, it is found that the resultant curves are almost identical for all of the commonly dated K–Ar minerals – biotite, phlogopite, muscovite, hornblende and orthoclase – in spite of differences in their diffusion parameters. A common curve for Ar diffusion can be derived by least-squares fitting of all the Ar diffusion data and provides a way of predicting a “model” closure temperature Tcm from a single diffusion coefficient DM at temperature TM. Preliminary diffusion data for a labradorite lead to a Tcm of 507 ± 17 °C and a corresponding activation energy of about 65 kcal/mol, given a grain size of 200 μm and a cooling rate of 5 °C/Ma. Curves for He diffusion in silicates (augite, quartz and sanidine) also overlap to a significant degree, both among themselves and with the Ar model curve, suggesting that a single model curve may be a good representation of noble gas closure temperatures in silicates. An analogous model curve for a selection of 18O data can also be constructed, but this curve differs from the Ar model curve. A single model curve for cationic species does not appear to exist, however, suggesting that chemical bonding relationships between the ionic size/charge and crystal structure may influence the closure temperatures of diffusing cations. An indication of the degree of overlap among the various curves for Ar, He, 18O and cations is also obtained by considering the dimensionless parameter E/RTc*; for the noble gases and 18O, E/RTc* values for the respective minerals are very similar, whereas for cations, there is significant dispersion. Given these constraints, this may be a potential method of estimating closure temperatures for certain diffusing species when there are limited diffusion data.