Petrological geodynamics of mantle melting II. AlphaMELTS + multiphase flow: Dynamic fractional melting

Abstract

In this second installment of a series that aims to investigate the dynamic interaction between the composition and abundance of the solid mantle and its melt products, the classic interpretation of fractional melting is extended to account for the dynamic nature of the process. A multiphase numerical flow model is coupled with the program AlphaMELTS, which provides at the moment possibly the most accurate petrological description of melting based on thermodynamic principles. The conceptual idea of this study is based on a description of the melting process taking place along a 1-D vertical ideal column where chemical equilibrium is assumed to apply in two local sub-systems separately on some spatial and temporal scale. The solid mantle belongs to a local sub-system (ss1) that does not interact chemically with the melt reservoir which forms a second sub-system(ss2). The local melt products are transferred in the melt sub-system ss2 where the melt phase eventually can also crystallize into a different solid assemblage and will evolve dynamically. The main difference with the usual interpretation of fractional melting is that melt is not arbitrarily and instantaneously extracted from the mantle, but instead remains a dynamic component of the model, hence the process is named dynamic fractional melting (DFM). Some of the conditions that may affect the DFM model are investigated in this study, in particular the effect of temperature, mantle velocity at the boundary of the mantle column. A comparison is made with the dynamic equilibrium melting (DEM) model discussed in the first installment. The implications of assuming passive flowor active floware also considered to some extent. Complete data files of most of the DFM simulations, four animations and two new DEM simulations (passive/active flow) are available following the instructions in the Supplementary Material.