Monitoring an earthfill dam using differential SAR interferometry: La Pedrera dam, Alicante, Spain

Surface displacement at the dykes of La Pedrera reservoir (SE Spain) has been measured by satellite differential Synthetic Aperture Radar (SAR) interferometry. At the main dyke, a displacement of about 13 cm along the satellite line of sight has been estimated between August 1995 and May 2010, from a dataset composed by ERS-1, ERS-2 and Envisat-ASAR images. Two independent short-term processing tasks were also carried out with ERS-2/Envisat-ASAR (from June 2008 to May 2010) and TerraSAR-X (from August 2008 to June 2010) images which have shown similar spatial and temporal displacement patterns. The joint analysis of historical instrument surveys and DInSAR-derived data has allowed the identification of a long-term deformation process which is reflected at the dam's surface and is also clearly recognizable in the inspection gallery. The plausible causes of the displacements measured by DInSAR are also discussed in the paper. Finally, DInSAR data have been used to compute the long-term settlement of La Pedrera dam, showing a good agreement with external studies. Consequently, this work demonstrates the integration of DInSAR with in-situ techniques which helps provide a complete spatial vision of the displacements in the dam thereby helping to differentiate the causal mechanisms.

Monitorización de la subsidencia del terreno en la Vega Media del río Segura mediante Interferometría SAR diferencial (DInSAR)

Ground subsidence caused by aquifer exploitation is a geotechnical hazard that affects wide areas, causing important economic damages. This phenomenon is due to soil consolidation produced by the increase of effective stress caused by piezometric depletion. In this work a Remote Sensing Technique called Coherent Pixel (CPT) is applied to monitor subsidence in the Vega Media of the Segura River during a period of eleven years.

Radar interferometry techniques for the study of ground subsidence phenomena: a review of practical issues through cases in Spain

Subsidence related to multiple natural and human-induced processes affects an increasing number of areas worldwide. Although this phenomenon may involve surface deformation with 3D displacement components, negative vertical movement, either progressive or episodic, tends to dominate. Over the last decades, differential SAR interferometry (DInSAR) has become a very useful remote sensing tool for accurately measuring the spatial and temporal evolution of surface displacements over broad areas. This work discusses the main advantages and limitations of addressing active subsidence phenomena by means of DInSAR techniques from an end-user point of view. Special attention is paid to the spatial and temporal resolution, the precision of the measurements, and the usefulness of the data. The presented analysis is focused on DInSAR results exploitation of various ground subsidence phenomena (groundwater withdrawal, soil compaction, mining subsidence, evaporite dissolution subsidence, and volcanic deformation) with different displacement patterns in a selection of subsidence areas in Spain. Finally, a cost comparative study is performed for the different techniques applied.

InSAR deformation measurements of the earthquake cycle in transcurrent tectonic domains, analytical and analog modeling

I applied the SBAS-DInSAR method to the Mattinata Fault (MF) (Southern Italy) and to the Doruneh Fault System (DFS) (Central Iran). In the first case, I observed limited internal deformation and determined the right lateral kinematic pattern with a compressional pattern in the northern sector of the fault. Using the Okada model I inverted the observed velocities defining a right lateral strike slip solution for the MF. Even if it fits the data within the uncertainties, the modeled slip rate of 13-15 mm yr-1 seems too high with respect to the geological record. Concerning the Western termination of DFS, SAR data confirms the main left lateral transcurrent kinematics of this fault segment, but reveal a compressional component. My analytical model fits successfully the observed data and quantifies the slip in ~4 mm yr-1 and ~2.5 mm yr-1 of pure horizontal and vertical displacement respectively. The horizontal velocity is compatible with geological record. I applied classic SAR interferometry to the October–December 2008 Balochistan (Central Pakistan) seismic swarm; I discerned the different contributions of the three Mw > 5.7 earthquakes determining fault positions, lengths, widths, depths and slip distributions, constraining the other source parameters using different Global CMT solutions. A well constrained solution has been obtained for the 09/12/2008 aftershock, whereas I tested two possible fault solutions for the 28-29/10/08 mainshocks. It is not possible to favor one of the solutions without independent constraints derived from geological data. Finally I approached the study of the earthquake-cycle in transcurrent tectonic domains using analog modeling, with alimentary gelatins like crust analog material. I successfully joined the study of finite deformation with the earthquake cycle study and sudden dislocation. A lot of seismic cycles were reproduced in which a characteristic earthquake is recognizable in terms of displacement, coseismic velocity and recurrence time.

Digital elevation model (DEM) of the ice sheet in western Dronning Maud Land, Antarctica

In this paper, a new digital elevation model (DEM) is derived for the ice sheet in western Dronning Maud Land, Antarctica. It is based on differential interferometric synthetic aperture radar (SAR) from the European Remote Sensing 1/2 (ERS-1/2) satellites, in combination with ICESat's Geoscience Laser Altimeter System (GLAS). A DEM mosaic is compiled out of 116 scenes from the ERS-1 ice phase in 1994 and the ERS-1/2 tandem mission between 1996 and 1997 with the GLAS data acquired in 2003 that served as ground control. Using three different SAR processors, uncertainties in phase stability and baseline model, resulting in height errors of up to 20 m, are exemplified. Atmospheric influences at the same order of magnitude are demonstrated, and corresponding scenes are excluded. For validation of the DEM mosaic, covering an area of about 130,000 km**2 on a 50-m grid, independent ICESat heights (2004-2007), ground-based kinematic GPS (2005), and airborne laser scanner data (ALS, 2007) are used. Excluding small areas with low phase coherence, the DEM differs in mean and standard deviation by 0.5 +/- 10.1, 1.1 +/- 6.4, and 3.1 +/- 4.0 m from ICESat, GPS, and ALS, respectively. The excluded data points may deviate by more than 50 m. In order to suppress the spatially variable noise below a 5-m threshold, 18% of the DEM area is selectively averaged to a final product at varying horizontal spatial resolution. Apart from mountainous areas, the new DEM outperforms other currently available DEMs and may serve as a benchmark for future elevation models such as from the TanDEM-X mission to spatially monitor ice sheet elevation.

Accumulation, velocities, fluxes and mass balance of glacial tributaries to the Lambert Glacier-Amery Ice Shelf system, Antarctica

By incorporating recently available remote sensing data, we investigated the mass balance for all individual tributary glacial basins of the Lambert Glacier-Amery Ice Shelf system, East Antarctica. On the basis of the ice flow information derived from SAR interferometry and ICESat laser altimetry, we have determined the spatial configuration of eight tributary drainage basins of the Lambert-Amery glacial system. By combining the coherence information from SAR interferometry and the texture information from SAR and MODIS images, we have interpreted and refined the grounding line position. We calculated ice volume flux of each tributary glacial basin based on the ice velocity field derived from Radarsat three-pass interferometry together with ice thickness data interpolated from Australian and Russian airborne radio echo sounding (RES) surveys and inferred from ICESat laser altimetry data. Our analysis reveals that three tributary basins have a significant net positive imbalance, while five other subbasins are slightly positive or close to zero balance. Overall, in contrast to previous studies, we find that the grounded ice in Lambert Glacier-Amery Ice Shelf system has a positive mass imbalance of 22.9 ± 4.4 Gt/a. The net basal melting for the entire Amery Ice Shelf is estimated to be 27.0 ± 7.0 Gt/a. The melting rate decreases rapidly from the grounding zone to the ice shelf front. Significant basal refreezing is detected in the downstream section of the ice shelf. The mass balance estimates for both the grounded ice sheet and the ice shelf mass differ substantially from other recent estimates.

Advanced interpretation of subsidence in Murcia (SE Spain) using A-DInSAR data – modelling and validation

Subsidence is a natural hazard that affects wide areas in the world causing important economic costs annually. This phenomenon has occurred in the metropolitan area of Murcia City (SE Spain) as a result of groundwater overexploitation. In this work aquifer system subsidence is investigated using an advanced differential SAR interferometry remote sensing technique (A-DInSAR) called Stable Point Network (SPN). The SPN derived displacement results, mainly the velocity displacement maps and the time series of the displacement, reveal that in the period 2004–2008 the rate of subsidence in Murcia metropolitan area doubled with respect to the previous period from 1995 to 2005. The acceleration of the deformation phenomenon is explained by the drought period started in 2006. The comparison of the temporal evolution of the displacements measured with the extensometers and the SPN technique shows an average absolute error of 3.9±3.8 mm. Finally, results from a finite element model developed to simulate the recorded time history subsidence from known water table height changes compares well with the SPN displacement time series estimations. This result demonstrates the potential of A-DInSAR techniques to validate subsidence prediction models as an alternative to using instrumental ground based techniques for validation.

Relationship between piezometric level and ground deformations measured by means of DInSAR in the Vega Media of the Segura River (Spain)

Differential SAR Interferometry (DInSAR) is a remote sensing method with the well demonstrated ability to monitor geological hazards like earthquakes, landslides and subsidence. Among all these hazards, subsidence involves the settlement of the ground surface affecting wide areas. Frequently, subsidence is induced by overexploitation of aquifers and constitutes a common problem that affects developed societies. The excessive pumping of underground water decreases the piezometric level in the subsoil and, as a consequence, increases the effective stresses with depth causing a consolidation of the soil column. This consolidation originates a settlement of ground surface that must be withstood by civil structures built on these areas. In this paper we make use of an advanced DInSAR approach - the Coherent Pixels Technique (CPT) [1] - to monitor subsidence induced by aquifer overexploitation in the Vega Media of the Segura River (SE Spain) from 1993 to the present. 28 ERS-1/2 scenes covering a time interval of about 10 years were used to study this phenomenon. The deformation map retrieved with CPT technique shows settlements of up to 80 mm at some points of the studied zone. These values agree with data obtained by means of borehole extensometers, but not with the distribution of damaged buildings, well points and basements, because the occurrence of damages also depends on the structural quality of the buildings and their foundations. The most interesting relationship observed is the one existing between piezometric changes, settlement evolution and local geology. Three main patterns of ground surface and piezometric level behaviour have been distinguished for the study zone during this period: 1) areas where deformation occurs while ground conditions remain altered (recent deformable sediments), 2) areas with no deformation (old and non-deformable materials), and 3) areas where ground deformation mimics piezometric level changes (expansive soils). The temporal relationship between deformation patterns and soil characteristics has been analysed in this work, showing a delay between them. Moreover, this technique has allowed the measurement of ground subsidence for a period (1993-1995) where no instrument information was available.

King George Island TanDEM-X DEM, link to GeoTIFF

We have generated a new digital elevation model for entire King George Island, Antarctica, using summer TanDEM-X bistatic SAR satellite data. The data was processed using differential SAR interferometry with an older DEM as reference. 4 TanDEM-X scenes from January 2012 were used as input. The new DEM was referenced to and validated against DGPS measurements. Height values are given in reference to ellipsoid (WGS84).