Numerical Modeling of Surface Deformation due to Magma Chamber Inflation/Deflation in a Heterogeneous Viscoelastic Half-space

Abstract:

Interpreting surface deformation patterns in terms of deeper processes in regions of active magmatism is challenging and inherently non-unique. This study focuses on interpreting the unusual sombrero-shaped pattern of surface deformation in the Altiplano Puna region of South America, which has previously been modeled as the effect of an upwelling diapir of material in the lower crust. Our goal is to investigate other possible interpretations of the surface deformation feature using a suite of viscoelastic models with varying material heterogeneity. We use the finite-element code PyLith to study surface deformation due to a buried time-varying (periodic) overpressure source, a magma body, at depth within a viscoelastic half-space. In our models, the magma-body is a penny-shaped crack, with a cylindrical region above the crack that is weak relative to the surrounding material. We initially consider a magma body within a homogeneous viscoelastic half-space to determine the effect of the free surface upon deformation above and beneath the source region. We observe a complex depth-dependent phase relationship between stress and strain for elements that fall between the ground surface and the roof of the magma body. Next, we consider a volume of weak material (faster relaxation time relative to background) that is distributed with varying geometry around the magma body. We investigate how surface deformation is governed by the spatial distribution of the weak material and its rheologic parameters. We are able to reproduce a "sombrero" pattern of surface velocities for a range of models with material heterogeneity. The wavelength of the sombrero pattern is primarily controlled by the extent of the heterogeneous region, modulated by flexural effects. Our results also suggest an "optimum overpressure forcing frequency" where the lifetime of the sombrero pattern (a transient phenomenon due to the periodic nature of the overpressure forcing) reaches a maximum. Through further research we hope to better understand how the parameter space of our forward model controls the distribution of surface deformation and eventually develop a better understanding of the observed pattern of surface deformation in the Altiplano Puna.