Source mechanism of Vulcanian degassing at Popocatepetl Volcano, Mexico, determined from waveform inversions of very long period signals

[1] The source mechanism of very long period (VLP) signals accompanying volcanic degassing bursts at Popocatepetl is analyzed in the 15 - 70 s band by minimizing the residual error between data and synthetics calculated for a point source embedded in a homogeneous medium. The waveforms of two eruptions ( 23 April and 23 May 2000) representative of mild Vulcanian activity are well reproduced by our inversion, which takes into account volcano topography. The source centroid is positioned 1500 m below the western perimeter of the summit crater, and the modeled source is composed of a shallow dipping crack ( sill with easterly dip of 10 degrees) intersecting a steeply dipping crack ( northeast striking dike dipping 83 degrees northwest), whose surface extension bisects the vent. Both cracks undergo a similar sequence of inflation, deflation, and reinflation, reflecting a cycle of pressurization, depressurization, and repressurization within a time interval of 3 - 5 min. The largest moment release occurs in the sill, showing a maximum volume change of 500 - 1000 m(3), pressure drop of 3 - 5 MPa, and amplitude of recovered pressure equal to 1.2 times the amplitude of the pressure drop. In contrast, the maximum volume change in the dike is less ( 200 - 300 m 3), with a corresponding pressure drop of 1 - 2 MPa and pressure recovery equal to the pressure drop. Accompanying these volumetric sources are single-force components with magnitudes of 10(8) N, consistent with melt advection in response to pressure transients. The source time histories of the volumetric components of the source indicate that significant mass movement starts within the sill and triggers a mass movement response in the dike within a few seconds. Such source behavior is consistent with the opening of a pathway for escape of pent-up gases from slow pressurization of the sill driven by magma crystallization. The opening of this pathway and associated rapid evacuation of volcanic gases induces the pressure drop. Pressure recovery in the magma filling the sill is driven by diffusion of gases from the resulting supersaturated melt into bubbles. Assuming a penny-shaped crack at ambient pressure of 40 MPa, the observed pressure and volume variations can be modeled with the following attributes: crack radius ( 100 m), crack aperture ( 5 m), bubble number density (10(10) - 10(12) m(-3)), initial bubble radius (10(-6) m), final bubble radius (similar to 10(-5) m), and net decrease of gas concentration in the melt (0.01 wt %).