R. Lanari

Satellite radar interferometry time series analysis of surface deformation for Los Angeles, California

[1] The Los Angeles, California, metropolitan area is a tectonically active region with surface deformation that is a combination of fault related tectonics plus a variety of natural and anthropogenic signals. We apply the small baseline subset (SBAS) algorithm to produce an interferometric synthetic aperture radar (InSAR) time series analysis for the Los Angeles area using data acquired by the ERS satellites from late 1995 into 2002. The result is a space-time deformation product that can be exploited to view not only the smoothly varying long-term surface motion, but also its time varying patterns. Large seasonal oscillations of the Santa Ana aquifer observed in Southern California Integrated GPS Network (SCIGN) data are accurately matched in the InSAR time series, moreover, correlations of the InSAR time series with an annual sinusoid allows us to investigate the dynamics of the hydrologic system. INDEX TERMS: 1243 Geodesy and Gravity: Space geodetic surveys; 1294 Geodesy and Gravity: Instruments and techniques; 1803 Hydrology: Anthropogenic effects. Citation: Lanari, R., P. Lundgren, M. Manzo, and F. Casu (2004), Satellite radar interferometry time series analysis of surface deformation for Los Angeles, California, Geophys. Res. Lett., 31, L23613, doi:10.1029/2004GL021294.

Modeling surface deformation observed with synthetic aperture radar interferometry at Campi Flegrei caldera

Satellite radar interferometry of Campi Flegrei caldera, Italy, reveals a pattern of subsidence during the period 1993–1998. Interferograms spanning the first half of the observation period (1993–1995) have a lower amplitude and average rate of subsidence than those spanning either the second half (1995–1998) or the entire period (1993–1998), consistent with observations of a slowing down or reversal of subsidence during the first half of the observation period. We calculate a time series of deformation images relative to a reference image on the basis of a least squares inversion. During the observation period the maximum subsidence progresses at a rate of roughly 38 ± 2 mm/yr, with periods of no apparent subsidence in late 1996 to early 1997. To understand the characteristics of the source, we jointly invert pairs of ascending and descending differential interferograms spanning similar time intervals (first half, second half, or entire interval) of the period 1993–1998. In each case the joint inversion fits the two unwrapped interferograms with a similar subhorizontal rectangular contracting tensile dislocation striking roughly N98°E with dimensions ? 4 × 2 km and located beneath the city of Pozzuoli at a depth of 2.5–3 km. Inversion for a spheroidal or Mogi point source also produced reasonable fits but with progressively poorer overall fits to the data, respectively. Our inversion assuming a simple source in an elastic half?space does not include the possible effects of local structure on the surface deformation, a factor that may also reduce the need for an asymmetric source. The solution we find is consistent with other studies that suggest subsidence due to hydrothermal diffusion as the primary deformation mechanism during this phase of caldera deflation.

Dynamic deformation of Etna Volcano observed by satellite radar interferometry

Satellite radar interferometry of Mt. Etna volcano, Sicily, Italy, reveals a sequence of deformation characterized by deflation during the end of the 1993 eruption, inflation from 1993–1995 with an increase in the inflation rate immediately before its resumed eruptive activity in late 1995. This was followed by very low deformation levels during the following year. The source of the deformation changed from a depth of 9 km during deflation to more than 11–14 km during the subsequent inflation, consistent with a model in which deflation at shallower levels is followed by inflation at greater depth as the volcano system recharges from below before its next eruption. This study demonstrates that radar interferometry provides an important contribution towards understanding the dynamic deformation of volcanoes. By revealing large scale changes in their pre-eruption deformation rates, radar interferometry could play an important role in volcano eruption monitoring.

Urban subsidence inside the city of Napoli (Italy) Observed by satellite radar interferometry

Satellite radar interferometry reveals subsidence occurring in an area of about 2 km by 2 km inside the city of Napoli, Italy. The observations show a maximum (vertical) displacement of about 5 cm between 1992/93 and 1996, while the deformation signal decreases between 1996 and 1998. The investigated area is densely populated and has been influenced, between 1992 and 1995, by intensive construction of a new underground railway. These interferometric results demonstrate the high correlation between the observed deformation and the excavation of the tunnels for the development of the underground.

Ground deformation and source geometry of the 24 August 2016 Amatrice earthquake (Central Italy) investigated through analytical and numerical modeling of DInSAR measurements and structural-geological data

We investigate the ground deformation and source geometry of the 2016 Amatrice earthquake (Central Italy) by exploiting ALOS2 and Sentinel-1 coseismic differential interferometric synthetic aperture radar (DInSAR) measurements. They reveal two NNW-SSE striking surface deformation lobes, which could be the effect of two distinct faults or the rupture propagation of a single fault. We examine both cases through a single and a double dislocation planar source. Subsequently, we extend our analysis by applying a 3-D finite elements approach jointly exploiting DInSAR measurements and an independent, structurally constrained, 3-D fault model. This model is based on a double fault system including the two northern Gorzano and Redentore-Vettoretto faults (NGF and RVF) which merge into a single WSW dipping fault surface at the hypocentral depth (8 km). The retrieved best fit coseismic surface deformation pattern well supports the exploited structural model. The maximum displacements occur at 5–7 km depth, reaching 90 cm on the RVF footwall and 80 cm on the NGF hanging wall. The von Mises stress field confirms the retrieved seismogenic scenario.

Stress transfer in the Lazufre volcanic area, central Andes

eruptive centers situated in an area larger than 1800 km 2 and (2) a small-scale uplift located at Lastarria volcano, which is the only volcano to show strong fumarolic activity in decades, with most of the clear deformation apparently not observed before 2000. Both the large and small uplift signals can be explained by magmatic or hydrothermal sources located at about 13 km and 1 km deep, respectively. To test a possible relationship, we use numerical modeling and estimate that the depth inflating source increased the tensile stress close to the shallow source. We discuss how the deep inflating source may have disturbed the shallow one and triggered the observed deformation at Lastarria. Citation: Ruch, J., A. Manconi, G. Zeni, G. Solaro, A. Pepe, M. Shirzaei, T. R. Walter, and R. Lanari (2009), Stress transfer in the Lazufre volcanic area, central Andes, Geophys. Res. Lett., 36, L22303, doi:10.1029/2009GL041276.

Long-term deformation analysis of historical buildings through the advanced SBAS-DInSAR technique: the case study of the city of Rome, Italy

Monitoring of deformation phenomena affecting urban areas and man-made structures is of key relevance for the preservation of the artistic, archaeological and architectural heritage. The differential SAR interferometry (DInSAR) technique has already been demonstrated to be an effective tool for non-invasive deformation analyses over large areas by producing spatially dense deformation maps with centimetre to millimetre accuracy. Moreover, by exploiting long sequences of SAR data acquired by different sensors, the advanced DInSAR technique referred to as the small baseline subset (SBAS) approach allows providing long-term deformation time series, which are strategic for guaranteeing the monitoring of urban area displacements. In this work, we investigate the effectiveness of the two-scale multi-sensor SBAS-DInSAR approach to detect and monitor displacements affecting historical and artistic monuments. The presented results, achieved by applying the full resolution SBAS technique to a huge set of ERS-1/2 and ENVISAT data, spanning the 1992–2010 time interval and relevant to the city of Rome (Italy), show the capability of this approach to detect and analyse the temporal evolution of possible deformation phenomena affecting historical buildings and archaeological sites. Accordingly, our analysis demonstrates the effectiveness of the full resolution multi-sensor SBAS approach to operate as a surface deformation tool for supporting the study and conservation strategies of the historical, cultural and artistic heritage.

Cloud Computing for Earth Surface Deformation Analysis via Spaceborne Radar Imaging: A Case Study

We present a case study on the migration to a Cloud Computing environment of the advanced differential synthetic aperture radar interferometry (DInSAR) technique, referred to as Small BAseline Subset (SBAS), which is widely used for the investigation of Earth surface deformation phenomena. In particular, we focus on the SBAS parallel algorithmic solution, namely P-SBAS, that allows the production of mean deformation velocity maps and the corresponding displacement time-series from a temporal sequence of radar images by exploiting distributed computing architectures. The Cloud migration is carried out by encapsulating the overall P-SBAS application in virtual machines running on the Cloud; moreover, the Cloud resources provisioning and configuration phases are implemented in an automatic way. Such an approach allows us to preserve the P-SBAS parallelization strategy and to straightforwardly evaluate its performance within a Cloud environment by comparing it with those achieved on a HPC in-house cluster. The results we present were achieved by using the Amazon Elastic Compute Cloud (EC2) of the Amazon Web Services (AWS) to process SAR datasets collected by the ENVISAT satellite and show that, thanks to the Cloud resources availability and flexibility, large DInSAR data volumes can be processed through the P-SBAS algorithm in short time frames and at reduced costs. As a case study, the mean deformation velocity map of the southern California area has been generated by processing 172 ENVISAT images. By exploiting 32 EC2 instances this processing took less than 17 hours to complete, with a cost of USD 850. Considering the available PB-scale archives of SAR data and the upcoming huge SAR data flow relevant to the recently launched (April 2014) Sentinel-1A and the forthcoming Sentinel-1B satellites, the exploitation of Cloud Computing solutions is particularly relevant because of the possibility to provide Cloud-based multi-user services allowing worldwide scientists to quickly process SAR data and to manage and access the achieved DInSAR results.

A new algorithm for monitoring localized deformation phenomena based on small baseline differential SAR interferograms

This paper presents a new solution for detecting and following the temporal evolution of small scale deformation phenomena; in particular our approach extends the capability of the SBAS technique, presented in P. Berardino et al. (2001), which is mainly focused on investigating large scale deformations with spatial resolutions of about 100 m/spl times/100 m. The proposed technique relies on small baseline differential SAR (DIFSAR) interferograms only, but it is implemented by using two different sets of data generated at low (multi-look data) and high spatial resolution (single-look data), respectively. The former are used to identify and estimate, via the SBAS technique or O. Mora et al. (2001, 2002), possible atmospheric phase artifacts and large scale deformation patterns; the latter to detect, on the high resolution residual phase components, structures highly coherent in time (buildings, rocks, lava structures, etc.), identified jointly to their heights and displacements. In particular the estimation of the temporal evolution of these local deformations is easily implemented by applying the SVD technique. The presented algorithm has been tested with data acquired by the European Remote Sensing (ERS) satellites relative to the Campania area (Italy).

A new approach for analyzing the temporal evolution of Earth surface deformations based on the combination of DIFSAR interferograms

We present a new approach for the evaluation of the Earth surface deformation evolution based on the combination of several differential interferograms spanning a time interval of interest. In particular we present a method that extends the least squares combination technique presented by S. Usai et al. (2000) by applying the singular value decomposition (SVD) method. Experiments carried out on ERS data validate the proposed approach.