FEM inversion of DInSAR data of the 2009 L’Aquila Earthquake (Italy)

Abstract

Slip distribution on fault planes is usually retrieved by geodetic data assuming the local crust as an elastic, homogeneous and isotropic half-space. Realistic complexities such as non-homogeneous elastic structure and topographic relief can be handled only by numerical methods. However, such elaborated models are computationally expensive and are usually implemented for forward modelling rather than for inversion purposes. On the other hand, spatially dense geodetic data (e.g. DInSAR displacement maps) often reveal complex patterns of coseismic deformation, pointing out the oversimplification of the analytical models. We develop a procedure to perform inversion of geodetic data based on Finite Element (FE) method, accounting for a more realistic description of the Earth crust, e.g. medium heterogeneity, anisotropy, topographic relief. FE computed Green Functions are implemented in an inversion framework to constrain the fault slip distribution in complex media. The method is applied to the 2009 L’Aquila earthquake. Three DInSAR maps of coseismic displacement from ENVISAT and COSMO-Skymed satellites are inverted simultaneously. The FE model includes the heterogeneities of the local crust. Results from inversions highlight the non-negligible influence of the medium structure: homogeneous and heterogeneous models show discrepancies in the fault slip distribution values up to 20%.