Discrete element simulations of gravitational volcanic deformation

We have carried out two-dimensional particle dynamics simulations of granular piles subject to frictional Coulomb failure criteria to gain a first-order understanding of different modes of gravitational deformation within volcanoes. Under uniform basal and internal strength conditions, granular piles grow self-similarly, developing distinctive stratigraphies, morphologies, and structures. Piles constructed upon cohesive substrates exhibit particle avalanching, forming outward dipping strata and angle of repose slopes. Systematic decreases in basal strength lead to progressively deeper and steeper internal detachment faults and slip along a basal decollement; landslide forms grade from shallow slumps to deep-seated landslide and, finally, to axial subsidence and outward flank displacements, or volcanic spreading. Surface slopes decrease and develop concave up morphologies with decreasing decollement strength; depositional layers tilt progressively inward. Spatial variations in basal strength cause lateral transitions in pile structure, stratigraphy, and morphology. This approximation of volcanoes as Coulomb granular piles reproduces the richness of deformational structures and surface morphologies in many volcanic settings. The gentle slopes of Hawaiian volcanoes and Olympus Mons on Mars suggest weak basal decollements that enable volcanic spreading. High-angle normal faults, favored above weak decollements, are interpreted in both settings and may explain catastrophic sector collapse in Hawaii and broad aureole deposits surrounding Olympus Mons. In contrast, steeper slopes and shallow detachment faults predominate in the Canary Islands, thought to lack a weak decollement, favoring smaller, more frequent slope failures than predicted for Hawaii. The numerical results provide a useful predictive tool for interpreting dynamic behavior and associated geologic hazards of active volcanoes.