Crystal-bearing magmas are often assumed to be intrinsically non-Newtonian, but the mechanism underlying suspension shear thinning remains unconstrained. Here, we test the hypothesis that shear thinning in these suspensions arises solely from unrelaxed shear thinning in the melt phase or from viscous heating or both. We compile existing data for the rheology of high-temperature crystal-bearing silicate magmas across a wide range of conditions. In order to test this, we define a ‘lever’ function L(ϕ) which scales the strain rate in the melt phase between crystals relative to the bulk strain rate of the suspension as a function of crystal volume fraction ϕ. We show that, in general, the parameterisation of the amplification of strain rates via the L(ϕ) factor coupled with a shear-thinning law for the melt phase fully accounts for the non-Newtonian behaviour of the suspension. Next, we use the Brinkman number Br to demonstrate that some existing data derive from experiments that were subject to viscous heating, resulting in a further manifestation of apparent shear thinning. Taken together, our results provide a theoretical framework that predicts microphysical regimes for shear thinning. This implies that crystal-bearing magmas are Newtonian except where the scaled conditions for viscous heating are met, or when the strain rates in the melt phase approach the inverse of the structural relaxation timescale. This shear thinning due to non-Newtonian behaviour of the melt between crystals occurs very close to the brittle threshold, yielding an extension of the brittle failure criterion to crystal-bearing magmas, with direct implications for volcanic eruption dynamics.
Read full abstract