On the Hydrogranular Dynamics of Magmatic Gravity Currents

Magmatic processes are generally governed by multi-phase interactions of silicate liquid, crystals, and bubbles. However, the modes of dissipation and the manner that stress is transmitted are poorly understood. We use a model of a simple but widely applicable gravity current as a means to exemplify the hydrogranular dynamics in crystal-rich magmas. Viscous and lubrication forces are of special interest because they have a dual role in dispersal and mixing in a crystal-rich gravity current. For example, lubrication forces provide an initial apparent yield strength by inducing a negative pore pressure as crystals move apart. However, once the gravity current is underway, lubrication forces reduce the dissipation due to collision and frictional contact.The gravity current is initiated by a combination of toppling and sliding along a well-defined granular fault. This produces three distinct regimes: a quasi-static base, an overlying particle hump that translates in a quasi-plastic fashion by grain-passing and rolling until the angle of repose is reached, and a viscous particle current. The current initially forms a leading vortex at the head, but the loss of crystals by sedimentation-assisted granular capture by an upward growing particle front drains energy from the flow. The vortex is soon abandoned, but persists in the reservoir as a fossil feature of orphaned crystals in a smear of previous intercumulate fluid. The kinetic energy of the most active crystals decays in a dual fashion, initially linearly, then parabolically with a near symmetrical increase and loss of kinetic energy.There is very little entrainment and mixing between intercumulate and reservoir fluids from magmatic gravity currents. Only a thin seam of reservoir melt is captured by the base of the flow as it descends across the floor. Hence magmatic gravity currents, while producing modest amounts of crystal sorting, are not effective agents of mixing as lubrication and viscous forces inhibit interpenetration of reservoir fluid.