Welcome to the new version of CaltechAUTHORS. Login is currently restricted to library staff. If you notice any issues, please email coda@library.caltech.edu
Published December 1994 | metadata_only
Journal Article

Fast grain boundary: A FORTRAN-77 program for calculating the effects of retrograde interdiffusion of stable isotopes


Exchange of stable isotopes between coexisting minerals is recognized widely as an important factor in the interpretation of stable isotope geochemistry of plutonic and high-grade metamorphic rocks. Where retrogression has occurred without major recrystallization events, the rate limiting step for stable isotope exchange will be diffusion. The mathematics of diffusion are well known for many problems, but no analytical solution, including that for closure temperature, adequately describes the complex and highly variable controls of rate and mass balance that will dominate many diffusion processes in rocks. We have implemented a model describing diffusional exchange for rocks in which grain boundary diffusion is sufficiently rapid that a representative volume of rock (typically millimeter to centimeter) is able to have mutual equilibration of all grain boundaries for the time scale of cooling. This Fast Grain Boundary model explicitly links intracrystalline diffusion rates and abundances of all minerals in a rock, and allows study of the impact of rock type on stable isotope thermometry, retrogression, and zonation. The FORTRAN-77 program for the Fast Grain Boundary model presented here can be used with a personal computer to solve typical problems in minutes. Input includes the grain size(s), model abundance(s), diffusion coefficient, and fractionation factor for each constituent mineral, and a cooling rate for the rock. Output includes the diffusion profile and integrated (bulk) composition of every mineral in a rock, as well as the apparent temperatures that would be observed by applying bulk-mineral stable isotope thermometry to such a rock.

Additional Information

© 1994 Elsevier Science Ltd. Received 20 July 1993; accepted 20 May 1994; Available online 31 March 2003. This code was developed and run on a Macintosh 11-ci using an MPW compiler and language systems FORTRA 3.0.1. This work was supported by an NSF Graduate Fellowship to JME and NSF grants EAR- 9105709 (JWV) and EAR-9106313 (LPB). Thanks are due to Jean Bahr, Herb Wang, and Martha Gerdes for helpful discussions on this work. This manuscript and the program described were improved by two thorough and helpful anonymous reviews.

Additional details

August 22, 2023
August 22, 2023