The characteristics of magma reservoir failure beneath a volcanic edifice

Eruptions fed from subsurface reservoirs commonly construct volcanic edifices at the surface, and the growth of an edifice will in turn modify the subsurface stress state that dictates the conditions under which subsequent rupture of the inflating reservoir can occur. We re-examine this problem using axisymmetric finite element models of ellipsoidal reservoirs beneath conical edifices, explicitly incorporating factors (e.g., full gravitational loading conditions, an elastic edifice instead of a surface load, reservoir pressures sufficient to induce tensile rupture) that compromise previous solutions to illustrate why variations in rupture behavior can occur. Relative to half-space model results, the presence of an edifice generally rotates rupture toward the crest of a spherical reservoir, with increasing flank slope (for an edifice of constant volume) and larger edifices (or greater reservoir scaled depths) normally serving to enhance this trend. When non-spherical reservoirs are considered, the presence of an edifice amplifies previously identified half-space failure characteristics, shifting rupture to the crest more rapidly for prolate reservoirs while forcing rupture closer to the midpoint of oblate reservoirs. Rupture is always observed to occur in the σt orientation, and depending on where initial failure occurs rupture favors the initial emplacement of either lateral sills, circumferential intrusions or vertically ascending dikes. Ultimately, integration of our numerical model results with other information, for instance the sequence of intrusion/eruption events observed at a given volcano, can provide useful new insight into how a volcanoʹs subsurface magma plumbing system evolved. We demonstrate this process through application of our model to Summer Coon, a well-studied stratocone on Earth, and Ilithyia Mons, a large conical shield volcano on Venus.