Abstract - Pit craters and calderas are volcanic depressions produced by subsidence of a magma reservoir roof. To identify how geometric and mechanical factors may influence the structural evolution of this subsidence, we used 2D Distinct Element Method (DEM) numerical models. The reservoir host rock was represented as an assemblage of bonded circular particles that interact according to elastic-frictional laws. Varying particle and bond properties produced a range of bulk material properties characteristic of natural rock masses. Fracturing results when bonds break, once their shear or tensile strength is exceeded. The magma reservoir was represented as a region of non-bonded, low-friction particles. Withdrawal of magma was simulated by incrementally reducing the area of the reservoir particles. Resultant gravity-driven failure and subsidence of the reservoir roof were explicitly replicated. Interaction of the roof's strength, Young's modulus, thickness/diameter ratio (T/D), and of the reservoir's shape yields a variety of model structures and subsidence styles. In conceptual terms, four end-member subsidence styles developed: (1) 'central sagging', favored by low strength and low T/D; (2) 'central snapping', favored by high strength, low T/D and a sill-like reservoir shape; (3) 'single central block', favored by low to intermediate strength, high Young's modulus, and intermediate T/D; (4) 'multiple central blocks', favored by high strength, low Young's modulus and high T/D. Most model realizations incorporated some combination of each style, however. The models provide a framework for understanding natural pit crater or caldera structures, as at Nindiri (Nicaragua), Fernandina (Galapagos), Dolomieu (La Reunion), and Miyakejima (Japan).
Journal of Geophysical Research, 116, B07202, 2011.