A spatially resolved surface kinetic model for forsterite dissolution

The development of complex alteration layers on silicate mineral surfaces undergoing dissolution is a widely observed phenomenon. Given the complexity of these layers, most kinetic models used to predict rates of mineral-fluid interactions do not explicitly consider their formation. As a result, the relationship between the development of the altered layers and the final dissolution rate is poorly understood. To improve our understanding of the relationship between the alteration layer and the dissolution rate, we developed a spatially resolved surface kinetic model for olivine dissolution and applied it to a series of closed-system experiments consisting of three-phases (water (+/- NaCl), olivine, and supercritical CO2) at conditions relevant to in situ mineral carbonation (i.e. 60 degrees C, 100 bar CO2). We also measured the corresponding delta Mg-26/24 of the dissolved Mg during early stages of dissolution. Analysis of the solid reaction products indicates the formation of Mg-depleted layers on the olivine surface as quickly as 2 days after the experiment was started and before the bulk solution reached saturation with respect to amorphous silica. The delta Mg-26/24 Mg of the dissolved Mg decreased by approximately 0.4 parts per thousand in the first stages of the experiment and then approached the value of the initial olivine (-0.35 parts per thousand) as the steady-state dissolution rate was approached. We attribute the preferential release of Mg-24 to a kinetic effect associated with the formation of a Mg-depleted layer that develops as protons exchange for Mg2+. We used experimental data to calibrate a surface kinetic model for olivine dissolution that includes crystalline olivine, a distinct active layer from which Mg can be preferentially removed, and secondary amorphous silica precipitation. By coupling the spatial arrangement of ions with the kinetics, this model is able to reproduce both the early and steady-state long-term dissolution rates, and the kinetic isotope fractionation. In the early stages of olivine dissolution the overall dissolution rate is controlled by exchange of protons for Mg, while the steady-state dissolution rate is controlled by the net removal of both Mg and Si from the active layer. Modeling results further indicate the importance of the spatial coupling of individual reactions that occur during olivine dissolution. The inclusion of Mg isotopes in this study demonstrates the utility of using isotopic variations to constrain interfacial mass transfer processes. Alternative kinetic frameworks, such as the one presented here, may provide new approaches for modeling fluid-rock interactions. (c) 2015 Elsevier Ltd. All rights reserved.