“…In simple form, these models consist of the laminar 1D steady state ascent of magma through a rigid cylinder [e.g., Jaupart, 1996;Sparks, 1999, 2002;Mastin and Ghiorso, 2000]. Even these simple models may display strongly nonlinear behavior and have multiple steady state solutions due to nonlinearity in constitutive laws and dependencies between different parameters, and additional complexity is introduced by consideration of lateral variation of magma properties in the conduit [Massol et al, 2001;Collier and Neuberg, 2006;Mastin, 2005], non-Newtonian rheologies [Melnik and Sparks, 2005], noncylindrical or depth-varying conduit geometries [Costa et al, 2007a[Costa et al, , 2007bde' Michieli Vitturia et al, 2008], and elastic or viscoelastic behavior of the conduit walls and/or inclusion of the magma chamber feeding the eruption [Maeda, 2000;Barmin et al, 2002;Costa et al, 2007a].…”