Low magnitude earthquakes generating significant subsidence: the Lunigiana case study

We applied the Differential Interferometric Synthetic Aperture Radar (DInSAR) technique to investigate and measure surface displacements due to the ML 5.2, June 21, 2013, earthquake occurred in the Apuan Alps (NW Italy) at a depth of about 5 km. The Centroid Moment Tensor (CMT) solution from INGV indicates an almost pure normal fault mechanism. Two differential interferograms showing the coseismic displacement were generated using X- band and C-band data respectively. The X-Band interferogram was obtained from a Cosmo-SkyMed ascending pair (azimuth -7.9° and incidence angle 40°) with a time interval of one day (June 21 - June 22) and 139 m spatial baseline, covering an area of about 40x40 km around the epicenter. The topographic phase component was removed using the 90 m SRTM DEM. The C-Band interferferogram was computed from two RADARSAT-2 Standard-3 (S3) images, characterized by 24 days temporal and 69 m spatial baselines, acquired on June 18 and July 12, 2013 on ascending orbit (azimuth -10.8°) with an incidence angle of 34° and covering 100x100 km area around the epicenter. The topographic phase component was removed using 30 m ASTER DEM. Adaptive filtering, phase unwrapping with Minimum Cost Flow (MCF) algorithm and orbital refinement were also applied to both interferograms. We modeled the observed SAR deformation fields using the Okada analytical formulation within a nonlinear inversion scheme, and found them to be consistent with a fault plane dipping towards NW at an angle of about 45°. In spite of the small magnitude, this earthquake produces a surface subsidence of about 1.5 cm in the Line-Of-Sight (LOS) direction, corresponding to about 3 cm along the vertical axis, that can be observed in both interferograms and appears consistent with the normal fault mechanisms.