We generated a 13‐year InSAR time series from 1995–2008 to investigate the spatiotemporal characteristics of two neighboring volcano's deformations for the Lazufre volcanic area, central Andes. The ...data reveal two scales of uplift initiating during the observation time: (1) a large‐scale uplift started in 1997 that shows an increase of the mean uplift rate of up to 3.2 cm/yr, now affecting several eruptive centers situated in an area larger than 1800 km2 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.
We integrate geologic, structural, leveling and Differential SAR Interferometry data to show that Vesuvius began to spread onto its sedimentary substratum about 3,600 years ago. Moreover, we model ...the detected deformation with a solution of the lubrication approximation of the Navier‐Stokes equations to show that spreading may continue for about 7,200 years more. Correlation of volcanic spreading with phases of the eruptive activity suggests that Plinian eruptions, which are thought to pose the major hazard, are less likely to occur in the near future.
The aim of this paper is to propose a methodology to perform inverse numerical modelling of slow landslides that combines the potentialities of both numerical approaches and well-known remote-sensing ...satellite techniques. In particular, through an optimization procedure based on a genetic algorithm, we minimize, with respect to a proper penalty function, the difference between the modelled displacement field and differential synthetic aperture radar interferometry (DInSAR) deformation time series. The proposed methodology allows us to automatically search for the physical parameters that characterize the landslide behaviour. To validate the presented approach, we focus our analysis on the slow Ivancich landslide (Assisi, central Italy). The kinematical evolution of the unstable slope is investigated via long-term DInSAR analysis, by exploiting about 20 years of ERS-1/2 and ENVISAT satellite acquisitions. The landslide is driven by the presence of a shear band, whose behaviour is simulated through a two-dimensional time-dependent finite element model, in two different physical scenarios, i.e. Newtonian viscous flow and a deviatoric creep model. Comparison between the model results and DInSAR measurements reveals that the deviatoric creep model is more suitable to describe the kinematical evolution of the landslide. This finding is also confirmed by comparing the model results with the available independent inclinometer measurements. Our analysis emphasizes that integration of different data, within inverse numerical models, allows deep investigation of the kinematical behaviour of slow active landslides and discrimination of the driving forces that govern their deformation processes.
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.
Precision synthetic aperture radar (SAR) data processing requires the compensation of the range-dependent range cell migration (RDRCM) phenomenon. This paper describes a new method that permits the ...compensation of the RDRCM effect using a nonstandard Fourier transform (FT). This operation is applied, in the two-dimensional Fourier domain, to the range signal spectrum and allows the compensation of the azimuth frequency-dependent scaling factor due to the RDRCM. The nonstandard FT is performed via a chirp z-transform that is carried out with one convolution operation and two phase multiplications.< >
Here we present the results of the inversion of a new geodetic data set covering the 2012 Emilia seismic sequence and the following 1 year of postseismic deformation. Modeling of the geodetic data ...together with the use of a catalog of 3‐D relocated aftershocks allows us to constrain the rupture geometries and the coseismic and postseismic slip distributions for the two main events (Mw 6.1 and 6.0) of the sequence and to explore how these thrust events have interacted with each other. Dislocation modeling reveals that the first event ruptured a slip patch located in the center of the Middle Ferrara thrust with up to 1 m of reverse slip. The modeling of the second event, located about 15 km to the southwest, indicates a main patch with up to 60 cm of slip initiated in the deeper and flatter portion of the Mirandola thrust and progressively propagated postseismically toward the top section of the rupture plane, where most of the aftershocks and afterslip occurred. Our results also indicate that between the two main events, a third thrust segment was activated releasing a pulse of aseismic slip equivalent to a Mw 5.8 event. Coulomb stress changes suggest that the aseismic event was likely triggered by the preceding main shock and that the aseismic slip event probably brought the second fault closer to failure. Our findings show significant correlations between static stress changes and seismicity and suggest that stress interaction between earthquakes plays a significant role among continental en echelon thrusts.
Key Points
We computed a new geodetic data set covering the 2012 Emilia (Northern Italy) seismic sequence
We defined the source models for the coseismic and postseismic phases of the two main events
We explored how these thrust events have interacted with each other during the sequence
Volcanoes deform dynamically due to changes in both their magmatic system and instability of their edifice. Mt. Etna features vigorous and almost continuous eruptive activity from its summit craters ...and periodic flank eruptions. Even though its shape is that of a large stratovolcano, its structure features two rift systems and a flank collapse structure similar to Hawaiian shield volcanoes. We analyze European remote sensing (ERS) satellite differential interferometric synthetic aperture radar (InSAR) data (1993–1996) for Mt. Etna spanning its quiescence from 1993 through the initiation of renewed eruptive activity in late 1995. We use synthetic aperture radar (SAR) data from both ascending and descending ERS satellite tracks. Comparison of independent interferograms covering the first 2 years of the inflationary period shows a pattern consistent with inflation of the volcano. Calculation of the tropospheric path delay based on meteorological data does not change this interpretation. Interferograms from late summer 1995–1996 show no significant deformation. Joint inversion of interferograms from ascending and descending satellite tracks require both inflation from a spheroidal magmatic source located beneath the summit at 5 km below sea level, and displacement of the east flank of Etna along a basal decollement. Both sources of deformation were contemporaneous within the resolution of our data and suggest that inflation of the central magma chamber acted to trigger slip of Etna's eastern flank. These results demonstrate that flank instability and recharge of a volcano's magma system must both be considered toward understanding how volcanoes work and in their hazard evaluation.
We present an efficient solution to mitigate phase unwrapping (PhU) errors that can affect a sequence of multi-temporal differential synthetic aperture radar (SAR) interferograms. To this aim, we ...propose a strategy that, starting from a properly chosen network of differential interferograms, complements PhU operations with an advanced multi-temporal region-growing (RG) procedure that exploits the space-time relationships among the computed interferograms. In particular, the proposed method implements an iterative procedure that, at each step, allows correcting a sequence of previously unwrapped interferograms at one selected pixel, namely candidate pixel, by exploiting the (unwrapped) phase values at its neighbouring 'seed' pixels (i.e. the ones already correctly unwrapped). Following their estimation, the unwrapped phases are then used to retrieve surface deformation products, such as mean deformation velocity maps and displacement time series, through (advanced) small baseline differential SAR interferometry (DInSAR) techniques. The effectiveness of the presented RG PhU algorithm is demonstrated by analysing a data set of SAR images acquired by the European Remote Sensing (ERS)-1/2 sensors over the megacity area of Istanbul, Turkey.
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. Die 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. Die 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.
We exploit the amplitude information of a sequence of synthetic aperture radar (SAR) images, acquired at different times, in order to generate displacement time-series in areas characterized by large ...and/or rapid deformation, the size of which is on the order of the image's pixel dimensions. We follow the same rationale of the Small BAseline Subset (SBAS) differential SAR interferometry (DInSAR) approach, by coupling the available SAR images into pairs characterized by a small separation between the acquisition orbits. We exploit the amplitudes of the selected image pairs in order to calculate the relative across-track (range) and along-track (azimuth) pixel-offsets (PO). Finally, we apply the SBAS inversion strategy to retrieve the range and azimuth displacement time-series. This approach, referred to as pixel-offset (PO-) SBAS technique, has been applied to a set of 25 ENVISAT SAR observations of the Sierra Negra caldera, Galápagos Islands, spanning the 2003-2007 time interval. The retrieved deformation time-series show the capability of the technique to detect and measure the large displacements affecting the inner part of the caldera that, in correspondence to the October 2005 eruption, reached several meters. Moreover, by comparing the PO-SBAS results to continuous GPS measurements, we estimate that the accuracy of the PO-SBAS time-series is on the order of 1/30th of a pixel for both range and azimuth directions.