Porosity evolution of mafic crystal mush during reactive flow

The emergence of the “mush paradigm” has raised several questions for conventional models of magma storage and extraction: how are melts extracted to form eruptible liquid-rich domains? What mechanism controls melt transport in mush-rich systems? Recently, reactive flow has been proposed as a major contributing factor in the formation of high porosity, melt-rich regions. Yet, owing to the absence of accurate geochemical simulations, the influence of reactive flow on the porosity of natural mush systems remains under-constrained. Here, we use a thermodynamically constrained model of melt-mush reaction to simulate the chemical, mineralogical, and physical consequences of reactive flow in a multi-component mush system. Our results demonstrate that reactive flow within troctolitic to gabbroic mushes can drive large changes in mush porosity. For example, primitive magma recharge causes an increase in the system porosity and could trigger melt channelization or mush destabilization, aiding rapid melt transfer through low-porosity mush reservoirs.