Benchmarking and developing numerical Finite Element models of volcanic deformation
Journal of Volcanology and Geothermal Research
Open Access funded by Natural Environment Research Council, under a Creative Commons license- Attribution 3.0 Unported (CC BY 3.0)
Volcanic deformation during an unrest episode can provide valuable insights into potential magmatic plumbing system dynamics. Paramount to this is a model linking the recorded ground movement to the causative source at depth. Most models employ analytical techniques due to their simplicity, but these are limited in their approach due to a number of necessary assumptions, and restricted to crude subsurface representations. We address this problem by providing guidelines and example model files to benchmark against a simple, analytical model with a numerical Finite Element approach using COMSOL Multiphysics. The boundary conditions should also be applicable to other Finite Element modeling packages. Then, due to the flexibility of the Finite Element method, this allows a progression of adding increasing complexities to reproduce the likely intricacies of the subsurface. We thus provide further guidelines and accompanying model files to incorporate subsurface heterogeneity, benchmarked viscoelastic rheology and temperature-dependent mechanics. In doing so, we highlight that setting up more integrated geodetic models is not particularly difficult and can alter inferred source characteristics and dynamics. The models are applied to Uturuncu volcano in southern Bolivia to demonstrate the approach.
This work was supported by the European Union, Framework Program 7 (grant #282759, “VUELCO”, and grant #308665, “MEDSUV”), the Natural Environmental Research Council (NE/G01843X/1) and the Royal Society (UF090006). We thank Maurizio Battaglia, an anonymous reviewer and editor Joan Martí for their constructive reviews.
This is the final version of the article. Available from the publisher via the DOI in this record.
Vol. 280, pp. 126 - 130