Please use this identifier to cite or link to this item:
Title: Evaluation of the effects of 3D diffusion, crystal geometry, and initial conditions on retrieved time-scales from Fe–Mg zoning in natural oriented orthopyroxene crystals
Authors: Krimer, Daniel
Costa, Fidel
Keywords: Volcano
Issue Date: 2017
Source: Krimer, D., & Costa, F. (2017). Evaluation of the effects of 3D diffusion, crystal geometry, and initial conditions on retrieved time-scales from Fe–Mg zoning in natural oriented orthopyroxene crystals. Geochimica et Cosmochimica Acta, 196, 271-288.
Series/Report no.: Geochimica et Cosmochimica Acta
Abstract: Volcano petrologists and geochemists increasingly use time-scale determinations of magmatic processes from modeling the chemical zoning patterns in crystals. Most determinations are done using one-dimensional traverses across a two-dimensional crystal section. However, crystals are three-dimensional objects with complex shapes, and diffusion and re-equilibration occurs in multiple dimensions. Given that we can mainly study the crystals in two-dimensional petrographic thin sections, the determined time-scales could be in error if multiple dimensional and geometrical effects are not identified and accounted for. Here we report the results of a numerical study where we investigate the role of multiple dimensions, geometry, and initial conditions of Fe–Mg diffusion in an orthopyroxene crystal with the view towards proper determinations of time scales from modeling natural crystals. We found that merging diffusion fronts (i.e. diffusion from multiple directions) causes ‘additional’ diffusion that has the greatest influence close to the crystal’s corners (i.e. where two crystal faces meet), and with longer times the affected area widens. We also found that the one dimensional traverses that can lead to the most accurate calculated time-scales from natural crystals are along the b- crystallographic axis on the ab-plane when model inputs (concentration and zoning geometry) are taken as measured (rather than inferred from other observations). More specifically, accurate time-scales are obtained if the compositional traverses are highly symmetrical and contain a concentration plateau measured through the crystal center. On the other hand, for two dimensional models the ab- and ac-planes are better suited if the initial (pre-diffusion) concentration and zoning geometry inputs are known or can be estimated, although these are a priory unknown, and thus, may be difficult to use in practical terms. We also found that under certain conditions, a combined one-dimensional and two-dimensional model performed on the ab-section can reveal the initial (pre-diffusional) concentration, and thus, offer a unique opportunity to recover lost petrologic information. The influence of three-dimensional diffusion on the one-dimensional two-dimensional model estimates is the combined result of the crystal shape, in particular the presence of facets, and somewhat of the aspect ratio. Our study focuses on Fe–Mg in orthopyroxene, but many of the effects we report are also applicable to other minerals and elements.
ISSN: 0016-7037
DOI: 10.1016/j.gca.2016.09.037
Schools: Asian School of the Environment 
Research Centres: Earth Observatory of Singapore 
Rights: © 2017 The Authors. Published by Elsevier Ltd. This is an open access article under the CC BY license (
Fulltext Permission: open
Fulltext Availability: With Fulltext
Appears in Collections:EOS Journal Articles

Citations 20

Updated on Jul 9, 2024

Web of ScienceTM
Citations 20

Updated on Oct 30, 2023

Page view(s) 50

Updated on Jul 13, 2024

Download(s) 50

Updated on Jul 13, 2024

Google ScholarTM




Items in DR-NTU are protected by copyright, with all rights reserved, unless otherwise indicated.