Improving Finite Element Optimization of InSAR-Derived Deformation Source Using Integrated Multiscale Approach

Parametric optimization/inversion of Interferometric Synthetic Aperture Radar (InSAR) measurements enables the modeling of the volcanic deformation source by considering the approximation of the analytic formulations or by defining refined scenarios within a Finite Element (FE) framework. However, the geodetic data modeling can lead to ambiguous solutions when constraints are unavailable, turning out to be time-consuming. In this work, we use an integrated multiscale approach for retrieving the geometric parameters of volcanic deformation sources and then constraining a Monte Carlo optimization of FE parametric modeling. This approach allows for contemplating more physically complex scenarios and more robust statistical solutions, and significantly decreasing computing time. We propose the Campi Flegrei caldera (CFc) case study, considering the 2019–2022 uplift phenomenon observed using Sentinel-1 satellite images. The workflow firstly consists of applying the Multiridge and ScalFun methods, and Total Horizontal Derivative (THD) technique to determine the position and horizontal sizes of the deformation source. We then perform two independent cycles of parametric FE optimization by keeping (I) all the parameters unconstrained and (II) constraining the source geometric parameters. The results show that the innovative application of the integrated multiscale approach improves the performance of the FE parametric optimization in proposing a reliable interpretation of volcanic deformations, revealing that (II) yields statistically more reliable solutions than (I) in an extraordinary tenfold reduction in computing time. Finally, the retrieved solution at CFc is an oblate-like source at approximately 3 km b.s.l. embedded in a heterogeneous crust.