The isotope effect for Mg-Fe interdiffusion in olivine and its dependence on crystal orientation, composition and temperature

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 Di/Dj = (mj/mi)β, 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. For the Fo83.4-Fo88.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 Fo88.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.