VENUSS: a unified finite-element model of solidifying lava

Abstract

The development of a solid rind or carapace at the surface of lava flows and domes results in a transition in deformation mechanism from dominantly viscous to elastic or plastic. This transition has a significant impact on the rate and style of emplacement, including on the construction of channelized flows, over-steepened margins, and flow advance due to lava breakouts. These processes are particularly important in subaqueous, subglacial, and extraterrestrial environments in which cooling is accelerated, requiring models specifically calibrated for these environments. We present a new numerical model, Viscous-Elastic Numerically Unified Solver for Solidifying flows (VENUSS), for cooling and solidifying free surface flows. The model couples a viscous fluid interior with an elastic shell whose thickness grows in response to cooling. As a demonstration of the impact of including a solidified crust in the flow model, we show that a dome-like shape fed from below with an elastic shell coupled to the basal topography results in more lateral expansion and less vertical uplift than a comparable highly-viscous rind, demonstrating the need for lateral transfer of stress in solid layers to accurately interpret and predict dome deformation.