Numerical simulation of plagioclase rim growth during magma ascent at Bezymianny Volcano, Kamchatka

Slow CaAl-NaSi interdiffusion in plagioclase crystals preserves chemical zoning of plagioclase in detail, which, along with strong dependence of anorthite content in plagioclase on melt composition, pressure, and temperature, make this mineral an important source of information on magma processes. A numerical model of zoned crystal growth is developed in the paper. The model is based on equations of multicomponent diffusion with diagonal cross-component diffusion terms and accounts for mass conservation on the melt-crystal interface and growth rate controlled by undercooling. The model is applied to the data of plagioclase rim zoning from several recent Bezymianny Volcano (Kamchatka) eruptions. We show that an equilibrium growth model cannot explain crystallization of naturally observed plagioclase during magma ascent. The developed non-equilibrium model reproduced natural plagioclase zoning and allowed magma ascent rates to be constrained. Matching of natural and simulated zoning suggests ascent from 100 to 50 MPa during 15-20 days. Magma ascent rate from 50 MPa to the surface varies from eruption to eruption: plagioclase zoning from the December 2006 eruption suggests ascent to the surface in less than 1 day, whereas plagioclase zoning from March 2000 and May 2007 eruptions are better explained by magma ascent over periods of more than 30 days). Based on comparison of diffusion coefficients for individual elements a mechanism of atomic diffusion during plagioclase crystallization is proposed.