Estimating volcanic deformation source parameters with a finite element inversion: The 2001–2002 unrest at Cotopaxi volcano, Ecuador

Publisher： John Wiley & Sons Inc

E-ISSN： 2169-9356|120|3|1473-1486

ISSN： 2169-9313

Source： JOURNAL OF GEOPHYSICAL RESEARCH: SOLID EARTH, Vol.120, Iss.3, 2015-03, pp. : 1473-1486

Disclaimer: Any content in publications that violate the sovereignty, the constitution or regulations of the PRC is not accepted or approved by CNPIEC.

Previous Menu Next

Abstract

AbstractDeformation at Cotopaxi was observed between 2001 and 2002 along with recorded seismicity beneath the northeast (NE) flank, despite the fact that the last eruption occurred in 1942. We use electronic distance meter deformation data along with the patterns of recorded seismicity to constrain the cause of this unrest episode. To solve for the optimum deformation source parameters we employ inverse finite element (FE) models that account for material heterogeneities and surface topography. For a range of source shapes the models converge on a shallow reservoir beneath the southwest (SW) flank. The individual best fit model is a small oblate‐shaped source, approximately 4–5 km beneath the summit, with a volume increase of roughly 20 × 106 m3. This SW source location contrasts with the NE seismicity locations. Subsequently, further FE models that additionally account for temperature‐dependent viscoelasticity are used to reconcile the deformation and seismicity simultaneously. Comparisons of elastic and viscous timescales allude to aseismic pressurization of a small magma reservoir in the SW. Seismicity in the NE is then explained through a mechanism of fluid migration from the SW to the NE along fault systems. We extend our analyses to further show that if future unrest crises are accompanied by measurable seismicity around the deformation source, this could indicate a higher magma supply rate and increased likelihood of a forthcoming eruption.