Abstract We investigate the origin of diversity of eruption styles in silicic volcanoes on the basis of a 1-dimensional steady conduit flow model that considers vertical relative motion between gas and liquid (i.e., vertical gas escape). The relationship between the assemblage of steady-state solutions in the conduit flow model and magma properties or geological conditions is expressed by a regime map in the parameter space of the ratio of liquid-wall friction force to liquid–gas interaction force (non-dimensional number ε), and a normalized conduit length Λ. The regime map developed in the companion paper shows that when ε is smaller than a critical value ε cr, a solution of explosive eruption exists for a wide range of Λ, whereas an effusive solution exists only when Λ ~ 1. On the other hand, when ε > ε cr, an effusive solution exists for a wide range of Λ. Diversity of eruption styles observed in nature is explained by the change in ε accompanied by the change in magma viscosity during magma ascent. As magma ascends, the magma viscosity increases because of gas exsolution and crystallization, leading to the increase in ε. For the viscosity of hydrous silicic magma at magma chamber, ε is estimated to be smaller than ε cr, indicating that an explosive solution exists for wide ranges of geological parameters. When magma flow rate is small, the viscosity of silicic magma drastically increases because of extensive crystallization at a shallow level in the conduit. In this case, ε can be greater than ε cr; as a result, a stable effusive solution co-exists with an explosive solution.