Numerical models of ionic diffusion in one and three dimensions: application to dehydration of mantle olivine
Résumé
The hydrogen content of nominally anhydrous minerals is of great interest, because it can influence many physical and mechanical properties of mantle rocks. Moreover, the hydrogen diffusion profiles can be used to constrain timescales related to magma eruptions. Here, we report models of ionic diffusion for trace elements in anisotropic crystals and apply them to hydrogen diffusing out of mantle-derived olivine. We first compare and discuss the characteristics of 1D and 3D models and show that only 3D anisotropic diffusion models can lead to diffusion profiles exhibiting non-equilibrium plateau at the center of the solid along the slowest axis, as measured in natural samples. In a second part, we discuss the differences between hydration and dehydration of olivine for diffusion that is linked to two different atomic sites involved in hydrogen mobility. Finally, we apply our 3D anisotropic model to previous results on mantle-derived olivine from Pali-aike to better characterize diffusion coefficients and their anisotropy that could be relevant for dehydration of olivine. Our results show that dehydration has to be strongly anisotropic, with a fast [100] axis and a significantly slower [001] axis.