Synthetic aperture radar (SAR) interferometry is a technique that provides high‐resolution measurements of the ground displacement associated with many geophysical processes. Advanced techniques ...involving the simultaneous processing of multiple SAR acquisitions in time increase the number of locations where a deformation signal can be extracted and reduce associated error. Currently there are two broad categories of algorithms for processing multiple acquisitions, persistent scatterer and small baseline methods, which are optimized for different models of scattering. However, the scattering characteristics of real terrains usually lay between these two end‐member models. I present here a new method that combines both approaches, to extract the deformation signal at more points and with higher overall signal‐to‐noise ratio than can either approach alone. I apply the combined method to data acquired over Eyjafjallajökull volcano in Iceland, and detect time‐varying ground displacements associated with two intrusion events.
New satellite missions (e.g., the European Space Agency's Sentinel‐1 constellation), advances in data downlinking, and rapid product generation now provide us with the ability to access ...space‐geodetic data within hours of their acquisition. To truly take advantage of this opportunity, we need to be able to interpret geodetic data in a prompt and robust manner. Here we present a Bayesian approach for the inversion of multiple geodetic data sets that allows a rapid characterization of posterior probability density functions (PDFs) of source model parameters. The inversion algorithm efficiently samples posterior PDFs through a Markov chain Monte Carlo method, incorporating the Metropolis‐Hastings algorithm, with automatic step size selection. We apply our approach to synthetic geodetic data simulating deformation of magmatic origin and demonstrate its ability to retrieve known source parameters. We also apply the inversion algorithm to interferometric synthetic aperture radar data measuring co‐seismic displacements for a thrust‐faulting earthquake (2015 Mw 6.4 Pishan earthquake, China) and retrieve optimal source parameters and associated uncertainties. Given its robustness and rapidity in estimating deformation source parameters and uncertainties, our Bayesian framework is capable of taking advantage of real‐time geodetic measurements. Thus, our approach can be applied to geodetic data to study magmatic, tectonic, and other geophysical processes, especially in rapid‐response operational settings (e.g., volcano observatories). Our algorithm is fully implemented in a MATLAB®‐based software package (Geodetic Bayesian Inversion Software) that we make freely available to the scientific community.
Key Points
We present a Bayesian approach for the inversion of geodetic data and demonstrate successful applications to synthetic and real data
Our approach allows rapid estimates of source parameters and uncertainties and is well suited for rapid‐response and operational settings
We have implemented our approach in a MATLAB®‐based software package (GBIS) that is made freely available to the scientific community
Synthetic aperture radar (SAR) interferometry is a technique that permits remote detection of deformation at the Earth's surface, and has been used extensively to measure displacements associated ...with earthquakes, volcanic activity and many other crustal deformation phenomena. Analysis of a time series of SAR images extends the area where interferometry can be successfully applied, and also allows detection of smaller displacements, through the reduction of error sources. Here, we review recent advances in time series SAR interferometry methods that further improve accuracy. This is particularly important when constraining displacements due to processes with low strain rates, such as interseismic deformation. We include examples of improved algorithms applied to image deformation associated with the 2010 eruption of Eyjafjallajökull volcano in Iceland, slow slip on the Guerrero subduction zone in Mexico, and tectonic deformation in western Anatolia, Turkey.
Interferometric Synthetic Aperture Radar (InSAR) is widely used to measure deformation of the Earth's surface over large areas and long time periods. A common strategy to overcome coherence loss in ...long-term interferograms is to use multiple multilooked shorter interferograms, which can cover the same time period but maintain coherence. However, it has recently been shown that using this strategy can introduce a bias (also referred to as a “fading signal”) in the interferometric phase. We isolate the signature of the phase bias by constructing “daisy chain” sums of short-term interferograms of different length covering identical 1-year time intervals. This shows that the shorter interferograms are more affected by this phenomenon and the degree of the effect depends on ground cover types; cropland and forested pixels have significantly larger bias than urban pixels and the bias for cropland mimics subsidence throughout the year, whereas forests mimics subsidence in the spring and heave in the autumn. We, propose a method for correcting the phase bias, based on the assumption, borne out by our observations, that the bias in an interferogram is linearly related to the sum of the bias in shorter interferograms spanning the same time. We tested the algorithm over a study area in western Turkey by comparing average velocities against results from a phase linking approach, which estimates the single primary phases from all the interferometric pairs, and has been shown to be almost insensitive to the phase bias. Our corrected velocities agree well with those from a phase linking approach. Our approach can be applied to global compilations of short-term interferograms and provides accurate long-term velocity estimation without a requirement for coherence in long-term interferograms.
•The characteristics of the phase bias in InSAR data are investigated.•A mitigation strategy is proposed to estimate the correction in each interferogram.•The proposed strategy is simple and effective in addressing the phase bias.•The proposed method provides a close performance to a phase linking approach.
Resolving changes in topography through time using accurate high‐resolution digital elevation models (DEMs) is key to understanding active volcanic processes. For the first time in a volcanic ...environment, we utilize very high‐resolution tri‐stereo optical imagery acquired by the Pleiades‐1 satellite constellation and generate a 1 m resolution DEM of Fogo Volcano, Cape Verde—the most active volcano in the Eastern Atlantic region. Point cloud density is increased by a factor of 6.5 compared to conventional stereo imagery, and the number of 1 m2 pixels with no height measurements is reduced by 43%. We use the DEM to quantify topographic changes associated with the 2014–2015 eruption at Fogo. Height differences between the posteruptive Pleiades‐1 DEM and the preeruptive topography from TanDEM‐X give a lava flow volume of 45.83 ± 0.02 × 106 m3, emplaced over an area of 4.8 km2 at a mean rate of 6.8 m3 s−1.
Key Points
Tri‐stereo photogrammetry at Fogo increases point cloud density by a factor of 6.5 and reduces area with no height measurements by 43%
Estimated accuracy of heights in generated tri‐stereo Pleiades‐1 DEM is <0.51 m
Lava flow volume and mean output rate of 2014–2015 Fogo eruption are 45.83 ± 0.02 × 106 m3 and 6.8 m3 s−1
In contrast to the glacier mass losses observed at other locations around the world, some glaciers in the High Mountains of Asia appear to have gained mass in recent decades. However, changes in ...digital elevation models indicate that glaciers in Karakoram and Pamir have gained mass, while recent laser altimetry data indicate mass gain centred on West Kunlun. Here, we obtain results that are essentially consistent with those from altimetry, but with two-dimensional observations and higher resolution. We produced elevation models using radar interferometry applied to bistatic data gathered between 2011 and 2014 and compared them to a model produced from bistatic data collected in 2000. The glaciers in West Kunlun, Eastern Pamir and the northern part of Karakoram experienced a clear mass gain of 0.043 ± 0.078~0.363 ± 0.065 m w.e. yr
. The Karakoram showed a near-stable mass balance in its western part (-0.020 ± 0.064 m w.e. yr
), while the Eastern Karakoram showed mass loss (-0.101 ± 0.058 m w.e. yr
). Significant positive glacier mass balances are noted along the edge of the Upper Tarim Basin and indicate a decreasing gradient from northeast to southwest.
Land subsidence in urban environments is an increasingly prominent aspect in the monitoring and maintenance of urban infrastructures. In this study we update the subsidence information over Rome and ...its surroundings (already the subject of past research with other sensors) for the first time using Copernicus Sentinel-1 data and open source tools. With this aim, we have developed a fully automatic processing chain for land deformation monitoring using the European Space Agency (ESA) SentiNel Application Platform (SNAP) and Stanford Method for Persistent Scatterers (StaMPS). We have applied this automatic processing chain to more than 160 Sentinel-1A images over ascending and descending orbits to depict primarily the Line-Of-Sight ground deformation rates. Results of both geometries were then combined to compute the actual vertical motion component, which resulted in more than 2 million point targets, over their common area. Deformation measurements are in agreement with past studies over the city of Rome, identifying main subsidence areas in: (i) Fiumicino; (ii) along the Tiber River; (iii) Ostia and coastal area; (iv) Ostiense quarter; and (v) Tivoli area. Finally, post-processing of Persistent Scatterer Inteferometry (PSI) results, in a Geographical Information System (GIS) environment, for the extraction of ground displacements on urban infrastructures (including road networks, buildings and bridges) is considered.
High-precision monitoring of infrastructure using artificial reflectors is possible with freely available Sentinel-1 data, but large reflectors are needed. We find that a triangular trihedral corner ...reflector should typically have at least 1-m inner leg length. As such large reflectors are often not feasible for use in urban areas for essential infrastructure monitoring, we designed a multiple corner-reflector array to replace a single corner reflector with an inner leg length of 1 m. In this case, we use four reflectors where each of them is a truncated triangular trihedral with an inner leg length of 0.33 m. We measured interferometric synthetic aperture radar (InSAR) amplitude, phase, and coherence of this reflector array with various configurations of alignments of the array. We find that as long as great care is taken in the relative positioning of the four corner reflectors, so that they constructively interfere, each horizontal or vertical configuration provides the expected amplitude, coherence, and phase stability. Applications of multiple small corner reflectors in urban areas range from essential infrastructure monitoring (e.g., bridges, overpasses, and tunnel constructions), through assessment of structural health of buildings, to monitoring highway and railway embankments. We show that the multiple corner array works when placed in a single InSAR resolution cell, but depending on the application, the number and projection of corner reflectors can be varied, as long as sufficient signal-to-clutter ratio is achieved in the area of interest.
We present here a new InSAR persistent scatterer (PS) method for analyzing episodic crustal deformation in non‐urban environments, with application to volcanic settings. Our method for identifying PS ...pixels in a series of interferograms is based primarily on phase characteristics and finds low‐amplitude pixels with phase stability that are not identified by the existing amplitude‐based algorithm. Our method also uses the spatial correlation of the phases rather than a well‐defined phase history so that we can observe temporally‐variable processes, e.g., volcanic deformation. The algorithm involves removing the residual topographic component of flattened interferogram phase for each PS, then unwrapping the PS phases both spatially and temporally. Our method finds scatterers with stable phase characteristics independent of amplitudes associated with man‐made objects, and is applicable to areas where conventional InSAR fails due to complete decorrelation of the majority of scatterers, yet a few stable scatterers are present.
After 20 years of quiescence, Fogo volcano erupted in November 2014. The eruption produced fast‐moving lava flows that traveled for several kilometers and destroyed two villages. This event ...represents the first episode of significant surface deformation imaged by the new European Space Agency's Sentinel‐1 satellite in its standard acquisition mode, Terrain Observation by Progressive Scans (TOPS), which differs from that of previous synthetic aperture radar (SAR) missions. We perform a Bayesian inversion of Sentinel‐1 TOPS SAR interferograms spanning the eruption and accurately account for variations in the TOPS line‐of‐sight vector when modeling displacements. Our results show that magma ascended beneath the Pico do Fogo cone and then moved laterally toward its southwestern flank, where the eruptive fissure opened. This study provides important insights into the inner workings of Fogo volcano and shows the potential of Sentinel‐1 TOPS interferometry for geophysical (e.g., volcano monitoring) applications.
Key Points
The 2014–2015 Fogo eruption is first event imaged with Sentinel‐1 TOPS
Eruption was fed by a single but complex dike
Surface deformation modeling of Sentinel‐1 TOPS‐mode InSAR data