DOE OSTI2018
Isotopic fractionation associated with diffusion in crystals is the most reliable means of understanding the origin of mineral zoning in igneous and metamorphic rocks. We have experimentally determined the relative diffusivities of iron isotopes in olivine as a function of crystallographic orientation, composition, and temperature. For two isotopes i and j of an element, the isotope effect for diffusion is parameterized as D i /D j = (m j /m i ) β , where β is a dimensionless parameter, and D and m stand for diffusivity and mass, respectively. A series of single crystal diffusion couple experiments were conducted at an oxygen fugacity of QFM – 1.5 at temperatures of 1200, 1300, and 1400 °C. Here, for the Fo 83.4 -Fo 88.8 composition pair, β Fe is isotropic and a value of 0.16 ± 0.09 can be used to describe diffusion along all major crystallographic axes in olivine. Based on our experiments and previously reported coupled Mg-Fe isotopic data, we also estimate β Mg = 0.09 ± 0.05 for this range of olivine composition. For the Fo 88.8 -Fo100 composition pair, β Fe becomes anisotropic with β Fe [100] = 0.11 ± 0.03, β Fe [010 ] = 0.14 ± 0.03 (both within error of the value measured for the Fo83.4-Fo88.8 pair), and β Fe [001] = 0.03 ± 0.03. For Fo# between 83.4 and 100, β Fe [100] and β Fe [010] are thus independent of composition. The reason why βFe [001] transitions from ~0.16 to ~0.03 close to the Mg-endmember is unclear. Over the temperature range studied, a dependence of βFe on temperature was not resolved. General analytical expressions are introduced to calculate isotopic fractionation as a function of distance, time, β, and the concentration contrast between the diffusing media for spherical, cylindrical, and planar geometries.