3D multi-source model of elastic volcanic ground deformation

Developments in Interferometric Synthetic Aperture Radar (InSAR) and GNSS (Global Navigation Satellite System) during the past decades have promoted significant advances in geosciences, providing high-resolution ground deformation data with dense spatio-temporal coverage. This large dataset can be exploited to produce accurate assessments of the primary processes occurring in geologically active areas. We present a new, original methodology to carry out a multi-source inversion of ground deformation data to better understand the subsurface causative processes. A nonlinear approach permits the determination of location, size and three-dimensional configuration, without any a priori assumption as to the number, nature or shape of the potential sources. The proposed method identifies a combination of pressure bodies and different types of dislocation sources (dip-slip, strike-slip and tensile) that represent magmatic sources and other processes such as earthquakes, landslides or groundwater-induced subsidence through the aggregation of elemental cells. This approach has the following features: (1) simultaneous inversion of the deformation components and/or line-of-sight (LOS) data; (2) simultaneous determination of diverse structures such as pressure bodies or dislocation sources, representing local and regional effects; (3) a fully 3D context; and (4) no initial hypothesis about the number, geometry or types of the causative sources is necessary. This methodology is applied to Mt. Etna (Southern Italy). We analyze the ground deformation field derived from a large InSAR dataset acquired during the January 2009 - June 2013 time period. The application of the inversion approach models several interesting buried structures as well as processes related to the volcano magmatic plumbing system, local subsidence within the Valle del Bove and seaward motion of eastern flank of the volcano. (C) 2020 Elsevier B.V. All rights reserved.