InSAR Deformation Research Articles

Inversion of reservoir parameters in Shuguang Oil Production Plant of the Liaohe Oilfield based on InSAR deformation

Inversion of reservoir parameters in Shuguang Oil Production Plant of the Liaohe Oilfield based on InSAR deformation

Land subsidence monitoring based on InSAR and inversion of aquifer parameters

In order to accurately separate the elastic and inelastic deformation information caused by aquifer compression in the land subsidence signal, and to invert the hydrogeological parameters of high spatial and temporal resolution to better apply the groundwater-ground subsidence model, a CWT (Continuous Wavelet Transform) separation method for aquifer elastic and inelastic deformation signals based on CWT is adopted, and the deformation signal is extracted by InSAR technology. The large-scale synthetic aperture radar dataset obtained by Envisat satellite from 2007 to 2009 is collected to obtain the surface deformation characteristic of the area by SBAS-InSAR technology, and then the independence provided by the observation well is used. Using the independent water level data provided by the observation wells, combined with the vertical InSAR deformation component and the head data, the CWT method is used to separate the periodic deformation signal components and long-term trends. Finally, the isolated signal component is used to invert the elastic and inelastic storage coefficient based on the ground subsidence model. The settlement signal separation method used in this paper makes up for the shortcomings of the two kinds of information in the previous settlement signal that are difficult to separate, which allows for more accurate inversion of aquifer parameters and helps to understand the aquifer parameters and continuously manage groundwater resources.

Sequential Estimation of Dynamic Deformation Parameters for SBAS-InSAR

The synthetic aperture radar (SAR) interferometry (InSAR) has been developed for more than 20 years for historical surface deformation reconstruction. In particular, the onboard Sentinel-1/A/B satellite, newly planned NASA-ISRO SAR (NISAR), and Germany Tandem-L will continue to provide unprecedented SAR data with an increased number of acquisitions. However, processing of real-time SAR data has been experiencing challenges regarding the InSAR deformation parameter estimation over a long time with the small baseline subsets (SBAS) InSAR technology. We use sequential adjustment for the estimation of the deformation parameters, which uses Bayesian estimation theory under the least square criteria to inverse long time-series deformation dynamically. Finally, both simulated and real Sentinel-1A SAR data verify the performance of the sequential estimation. It can be regarded as an effective data processing tool in the coming era of SAR big data.

Investigation of a Small Landslide in the Qinghai-Tibet Plateau by InSAR and Absolute Deformation Model

Landslides are one of the major geohazards in the Qinghai-Tibet Plateau, and have recently increased in both frequency and size. SAR interferometry (InSAR) has been widely applied in landslide research, but studies on monitoring small-scale landslides are rare. In this study, we investigated the performance of Small Baseline Subsets method (SBAS) in monitoring small-scale landslide and further developed a new deformation model to obtain the absolute deformation time series. The results showed that SBAS could well capture the small-scale landslide characteristics including spatiotemporal abnormal displacement and progressive failure processes. The newly developed absolute deformation model further detected the process of landslide details, such as instances of noticeable creeps induced by rainfall and snowmelt. Finally, a conceptual model of the kinematics-based failure mechanism for small-scale landslide was proposed. This study extended the monitoring capability of InSAR and improved our knowledge on the deformation in the frozen ground regions.

Co-Seismic and Postseismic Fault Models of the 2018 Mw 6.4 Hualien Earthquake Occurred in the Junction of Collision and Subduction Boundaries Offshore Eastern Taiwan

The ascending and descending InSAR deformations derived from ALOS-2 and Sentinel-1 satellite SAR images and GPS displacements are used to estimate the fault model of the 2018 Mw 6.4 Hualien earthquake. The sinistral strike-slip fault dipping to the west with a high dip angle of 89.4° and a rake angle of 201.7° is considered as the seismogenic fault of this event. This seismogenic fault also triggered the ruptures of the Milun fault, which dips to the east with a dip angle of ~72°, and an unknown west-dipping fault with a dip angle of 85.2°. Two predicted faulting models indicate that the InSAR deformation fields include more postseismic slip than those of the GPS data. The north segment of the Milun fault and west-dipping fault have been triggered by the rupture of the seismogenic fault, but the postseismic slip occurred only in the south segment of the Milun fault. The InSAR-derived co-seismic and postseismic faulting model suggests that the significant slip concentrates at depths of 2.4–15.0 km of the main fault, 0.0–14.0 km of the Milun fault. Only minor slip is detected on the west-dipping fault. The maximum fault slip of ca. 2.1 m is located at the depth of ca. 2.4 km under the Meilun Tableland. The Coulomb failure stress (CFS) change calculated by the co-seismic and postseismic faulting model shows that there is a significant CFS increase in the east of the south segment of the Milun fault, but few of the aftershocks occur in this area, which indicates a high risk of future seismic hazard.

MT-INSAR LANDSLIDE MONITORING WITH THE AID OF HOMOGENEOUS PIXELS FILTER

Abstract. SAR interferograms are often contaminated by random noises related to temporal decorrelation, geometrical decorrelation and thermal noises, which makes the fringes obscured and greatly decreases the density of the coherent target and the accuracy of InSAR deformation results, especially for the landslide monitoring in vegetated region and in rainy season. Two different SAR interferogram filtering methods, that is Goldstein filter and homogeneous pixels filter, for one specific landslide are compared. The results show that homogeneous pixels filter is better than Goldstein one for small-scale loess landslide monitoring, which can increase the density of monitoring points. Moreover, the precision of InSAR result can reach millimeter by comparing with GPS time series measurements.

INVESTIGATION ON INSAR TIME SERIES DEFORMATION MODEL CONSIDERING RHEOLOGICAL PARAMETERS FOR SOFT CLAY SUBGRADE MONITORING

Abstract. The stability control is one of the major technical difficulties in the field of highway subgrade construction engineering. Building deformation model is a crucial step for InSAR time series deformation monitoring. Most of the InSAR deformation models for deformation monitoring are pure empirical mathematical models, without considering the physical mechanism of the monitored object. In this study, we take rheology into consideration, inducing rheological parameters into traditional InSAR deformation models. To assess the feasibility and accuracy for our new model, both simulation and real deformation data over Lungui highway (a typical highway built on soft clay subgrade in Guangdong province, China) are investigated with TerraSAR-X satellite imagery. In order to solve the unknows of the non-linear rheological model, three algorithms: Gauss-Newton (GN), Levenberg-Marquarat (LM), and Genetic Algorithm (GA), are utilized and compared to estimate the unknown parameters. Considering both the calculation efficiency and accuracy, GA is chosen as the final choice for the new model in our case study. Preliminary real data experiment is conducted with use of 17 TerraSAR-X Stripmap images (with a 3-m resolution). With the new deformation model and GA aforementioned, the unknown rheological parameters over all the high coherence points are obtained and the LOS deformation (the low-pass component) sequences are generated.

GROUND SUBSIDENCE MONITORING WITH MT-InSAR AND MECHANISM INVERSION OVER XI’AN, CHINA

Abstract. The ancient Xi’an, China, has been suffering severe land subsidence and ground fissure hazards since the 1960s, which has affected the safety of Subways. Multi-sensor SAR data are conducted to monitor the latest complex ground deformation and its influence on subway line No.3 over Xi’an. Annual deformation rates have been retrieved to reveal the spatiotemporal evolution of ground subsidence in Xi’an city from 2013 to 2017. Meanwhile, the correlation between land subsidence and ground fissures are analyzed by retrieving the deformation differences in both sides of the fissures. Besides, the deformation along subway line No. 3 is analyzed, and the fast deformation section is quantitatively studied. Finally, a flat lying sill model with distributed contractions is implemented to model the InSAR deformation over YHZ subsidence center, which manifests that the ground deformation is mainly caused by groundwater withdrawal.

Synthesis of global satellite observations of magmatic and volcanic deformation: implications for volcano monitoring & the lateral extent of magmatic domains

Global Synthetic Aperture Radar (SAR) measurements made over the past decades provide insights into the lateral extent of magmatic domains, and capture volcanic process on scales useful for volcano monitoring. Satellite-based SAR imagery has great potential for monitoring topographic change, the distribution of eruptive products and surface displacements (InSAR) at subaerial volcanoes. However, there are challenges in applying it routinely, as would be required for the reliable operational assessment of hazard. The deformation detectable depends upon satellite repeat time and swath widths, relative to the spatial and temporal scales of volcanological processes. We describe the characteristics of InSAR-measured volcano deformation over the past two decades, highlighting both the technique’s capabilities and its limitations as a monitoring tool. To achieve this, we draw on two global datasets of volcano deformation: the Smithsonian Institution Volcanoes of the World database and the Centre for the Observation and Modelling of Earthquakes, Volcanoes and Tectonics volcano deformation catalogue, as well as compiling some measurement characteristics and interpretations from the primary literature.We find that a higher proportion of InSAR observations capture non-eruptive and non-magmatic processes than those from ground-based instrument networks, and that both transient (< month) and long-duration (> 5 years) deformation episodes are under-represented. However, satellite radar is already used to assess the development of extended periods of unrest and long-lasting eruptions, and improved spatial resolution and coverage have resulted in the detection of previously unrecognised deformation at both ends of the spatial scale (~ 10 to > 1000 km2). ‘Baseline’ records of past InSAR measurements, including ‘null’ results, are fundamental for any future interpretation of interferograms in terms of hazard‚ both by providing information about past deformation at an individual volcano, and for assessing the characteristics of deformation that are likely to be detectable (and undetectable) using InSAR.More than half of all InSAR deformation signals attributed to magmatic processes have sources in the shallow crust (< 5 km depth). While the depth distribution of InSAR-derived deformation sources is affected by measurement limitations, their lateral distribution provides information about the extent of active magmatic domains. Deformation is common (24% of all potentially magmatic events) at loci ≥5 km away from the nearest active volcanic vent. This demonstrates that laterally extensive active magmatic domains are not exceptional, but can comprise the shallowest part of trans-crustal magmatic systems in a range of volcanic settings.

Complex surface deformation monitoring and mechanism inversion over Qingxu-Jiaocheng, China with multi-sensor SAR images

Qingxu-Jiaocheng, China has been suffering severe land subsidence along with the development of ground fissure, which are controlled by local fault and triggered by groundwater withdrawal. With multi-sensor SAR images, we study the spatiotemporal evolution of ground deformation over Qingxu-Jiaocheng with an IPTA InSAR technique and assess the role of groundwater withdrawal to the observed deformation. Discrete GPS measurements are applied to verify the InSAR results. The RMSE of the differences between InSAR and GPS, i.e. ALOS and GPS and Envisat and GPS, are 5.7 mm and 6.3 mm in the LOS direction, respectively. The east-west and vertical components of the observed deformation from 2007 to 2010 are decomposed by using descending-track Envisat and ascending-track ALOS interferograms, indicating that the east-west component cannot be neglected when the deformation is large or the ground fissure is active. Four phases of land subsidence in the study region are successfully retrieved, and its spatiotemporal evolution is quantitatively analyzed. Lastly, a flat lying sill model with distributed contractions is implemented to model the InSAR deformation over Qingxu-Jiaocheng, which manifests that the ground deformation is mainly caused by groundwater withdrawal. This research provides new insights into the land subsidence monitoring and its mechanism inversion over Qingxu-Jiaocheng region.

Sequential combination of multi-source satellite observations for separation of surface deformation associated with serial seismic events

Single satellite geodetic technique has weakness for mapping sequence of ground deformation associated with serial seismic events, like InSAR with long revisiting period readily leading to mixed complex deformation signals from multiple events. It challenges the observation capability of single satellite geodetic technique for accurate recognition of individual surface deformation and earthquake model. The rapidly increasing availability of various satellite observations provides good solution for overcoming the issue. In this study, we explore a sequential combination of multiple overlapping datasets from ALOS/PALSAR, ENVISAT/ASAR and GPS observations to separate surface deformation associated with the 2011 Mw 9.0 Tohoku-Oki major quake and two strong aftershocks including the Mw 6.6 Iwaki and Mw 5.8 Ibaraki events. We first estimate the fault slip model of major shock with ASAR interferometry and GPS displacements as constraints. Due to the used PALSAR interferogram spanning the period of all the events, we then remove the surface deformation of major shock through forward calculated prediction thus obtaining PALSAR InSAR deformation associated with the two strong aftershocks. The inversion for source parameters of Iwaki aftershock is conducted using the refined PALSAR deformation considering that the higher magnitude Iwaki quake has dominant deformation contribution than the Ibaraki event. After removal of deformation component of Iwaki event, we determine the fault slip distribution of Ibaraki shock using the remained PALSAR InSAR deformation. Finally, the complete source models for the serial seismic events are clearly identified from the sequential combination of multi-source satellite observations, which suggest that the major quake is a predominant mega-thrust rupture, whereas the two aftershocks are normal faulting motion. The estimated seismic moment magnitude for the Tohoku-Oki, Iwaki and Ibaraki evens are Mw 9.0, Mw 6.85 and Mw 6.11, respectively.

Groundwater and Subsidence Modeling Combining Geological and Multi-Satellite SAR Data over the Alto Guadalentín Aquifer (SE Spain)

In the current context of climate change, improving groundwater monitoring and management is an important issue for human communities in arid environments. The exploitation of groundwater resources can trigger land subsidence producing damage in urban structures and infrastructures. Alto Guadalentín aquifer system in SE Spain has been exploited since 1960 producing an average piezometric level drop of 150 m. This work presents a groundwater model that reproduces groundwater evolution during 52 years with an average error below 10%. The geometry of the model was improved introducing a layer of less permeable and deformable soft soils derived from InSAR deformation and borehole data. The resulting aquifer system history of the piezometric level has been compared with ENVISAT deformation data to calculate a first-order relationship between groundwater changes, soft soil thickness, and surface deformation. This relationship has been validated with the displacement data from ERS and Cosmo-SkyMed satellites. The resulting regression function is then used as an empirical subsidence model to estimate a first approximation of the deformation of the aquifer system since the beginning of the groundwater extraction, reaching 1 to 5.5 m in 52 years. These rough estimations highlight the limitations of the proposed empirical model, requiring the implementation of a coupled hydrogeomechanical model.

Fusing Meter-Resolution 4-D InSAR Point Clouds and Optical Images for Semantic Urban Infrastructure Monitoring

Using synthetic aperture radar (SAR) interferometry to monitor long-term millimeter-level deformation of urban infrastructures, such as individual buildings and bridges, is an emerging and important field in remote sensing. In the state-of-the-art methods, deformation parameters are retrieved and monitored on a pixel basis solely in the SAR image domain. However, the inevitable side-looking imaging geometry of SAR results in undesired occlusion and layover in urban area, rendering the current method less competent for a semantic-level monitoring of different urban infrastructures. This paper presents a framework of a semantic-level deformation monitoring by linking the precise deformation estimates of SAR interferometry and the semantic classification labels of optical images via a 3-D geometric fusion and semantic texturing. The proposed approach provides the first “SARptical” point cloud of an urban area, which is the SAR tomography point cloud textured with attributes from optical images. This opens a new perspective of InSAR deformation monitoring. Interesting examples on bridge and railway monitoring are demonstrated.

Interpolation of GPS and Geological Data Using InSAR Deformation Maps: Method and Application to Land Subsidence in the Alto Guadalentín Aquifer (SE Spain)

Land subsidence resulting from groundwater extractions is a global phenomenon adversely affecting many regions worldwide. Understanding the governing processes and mitigating associated hazards require knowing the spatial distribution of the implicated factors (piezometric levels, lithology, ground deformation), usually only known at discrete locations. Here, we propose a methodology based on the Kriging with External Drift (KED) approach to interpolate sparse point measurements of variables influencing land subsidence using high density InSAR measurements. In our study, located in the Alto Guadalentín basin, SE Spain, these variables are GPS vertical velocities and the thickness of compressible soils. First, we estimate InSAR and GPS rates of subsidence covering the periods 2003–2010 and 2004–2013, respectively. Then, we apply the KED method to the discrete variables. The resulting continuous GPS velocity map shows maximum subsidence rates of 13 cm/year in the center of the basin, in agreement with previous studies. The compressible deposits thickness map is significantly improved. We also test the coherence of Sentinel-1 data in the study region and evaluate the applicability of this methodology with the new satellite, which will improve the monitoring of aquifer-related subsidence and the mapping of variables governing this phenomenon.

Influence of ocean tidal loading on InSAR offshore areas deformation monitoring

The ocean tide can cause the redistribution of the seawater mass, resulting in earth's surface deformation, namely ocean tidal loading (OTL). OTL vertical displacement may reach several centimeters, especially in coastal areas, so its effect in the field of high precision geodesy must be considered. This study concentrates on the influences of OTL on InSAR deformation measurements. We improve the osu.chinasea.2010 regional model and then compare the improved regional model with other regional models. It turns out that the improved regional model can achieve higher precision. Then we use it to replace the offshore part of the global model to generate the present model. We find that the displacement observed by the present model is 2–3 mm larger than that of other models on some sites. Finally, the present model is used to correct the deformation observed by InSAR of Shanghai and Los Angeles. A comparison between the displacements of IGS station with the corrected data shows that the OTL correction can improve the accuracy of InSAR deformation results by about 20%.

Use of a GPS-Derived Troposphere Model to Improve InSAR Deformation Estimates in the San Gabriel Valley, California

We evaluate the potential of troposphere models derived from ground meteorological data (pressure, temperature, and relative humidity) and Global Positioning System (GPS) data to improve InSAR measurements and models derived from them. We test this approach on an ERS-2/Envisat data set collected during a transient surface deformation episode that occurred from January to July 2005 in the San Gabriel Valley, southern California, USA. We find that the interferometric phase change observed over the corresponding period cannot be solely attributed to hydrological uplift associated with rising groundwater levels but also includes a significant contribution from differential tropospheric delay due to differing quantities of water vapor in the troposphere on the two SAR observation dates. We show that, if the tropospheric phase contribution is mistakenly interpreted as the range change associated with changes in groundwater storage, both the surface displacement and the groundwater storage coefficient may be overestimated by up to 30%. This method could be applied in real time where meteorological measurements are available near one or more GPS permanent site(s).

Confined aquifer head measurements and storage properties in the San Luis Valley, Colorado, from spaceborne InSAR observations

AbstractInterferometric Synthetic Aperture Radar (InSAR), a remote sensing technique for measuring centimeter‐level surface deformation, is used to estimate hydraulic head in the confined aquifer of the San Luis Valley (SLV), Colorado. Reconstructing head measurements from InSAR in agricultural regions can be difficult, as InSAR phase data are often decorrelated due to vegetation growth. Analysis of 17 L‐band ALOS PALSAR scenes, acquired between January 2007 and March 2011, demonstrates that comprehensive InSAR deformation measurements can be recovered over the vegetated groundwater basin with an improved processing strategy. Local skeletal storage coefficients and time delays between the head change and deformation are estimated through a joint InSAR‐well data analysis. InSAR subsidence estimates are transformed to head changes with finer temporal and spatial resolution than is possible using existing well records alone. Both InSAR and well data suggest that little long‐term water‐storage loss occurred in the SLV over the study period and that inelastic compaction was negligible. The seasonal head variations derived from InSAR are consistent with the existing well data at most locations where confined aquifer pumping activity dominates. Our results demonstrate the advantages of InSAR measurements for basin‐wide characterization of aquifer storage properties and groundwater levels over agricultural regions.

Multi-sensor InSAR Deformation Monitoring over Urban Area of Bratislava (Slovakia)

The integrated use of multiple Synthetic Aperture Radar (SAR) platforms for the deformation monitoring via satellite radar interferometry offers several perspectives for investigation of the behaviour of new and ageing structures, such as buildings and infrastructures, under varying or hazardous environment. Spanning almost 24 years of space-borne radar observations, this study aims to perform classical PSInSAR (Persistent Scatterer Interferometric Synthetic Aperture Radar) analysis incorporating measurements of ERS, Envisat, TerraSAR-X, Sentinel-1A and Radarsat-2 satellites. The results from the processing of different sensing geometries over Bratislava (Slovakia) urban area are presented, focusing on the description of characteristics associated with the specifics of every satellite platform in use. The discussion over technical feasibility of infrastructure monitoring is accompanied by the outline of possible future needs for the utilisation of the wealth source of information provided by the satellite radar imagery.

A persistent scatterer interpolation for retrieving accurate ground deformation over InSAR‐decorrelated agricultural fields

AbstractInterferometric synthetic aperture radar (InSAR) is a radar remote sensing technique for measuring surface deformation to millimeter‐level accuracy at meter‐scale resolution. Obtaining accurate deformation measurements in agricultural regions is difficult because the signal is often decorrelated due to vegetation growth. We present here a new algorithm for retrieving InSAR deformation measurements over areas with severe vegetation decorrelation using adaptive phase interpolation between persistent scatterer (PS) pixels, those points at which surface scattering properties do not change much over time and thus decorrelation artifacts are minimal. We apply this algorithm to L‐band ALOS interferograms acquired over the San Luis Valley, Colorado, and the Tulare Basin, California. In both areas, the pumping of groundwater for irrigation results in deformation of the land that can be detected using InSAR. We show that the PS‐based algorithm can significantly reduce the artifacts due to vegetation decorrelation while preserving the deformation signature.

Coseismic Slip Distribution of the 24 March 2011 Tarlay (Myanmar) Mw 6.8 Earthquake from ALOS PALSAR Interferometry

Abstract We investigate coseismic deformation of the 24 March 2011 M w 6.8 Tarlay, Myanmar, earthquake using ALOS PALSAR data from both descending and ascending passes. Using high‐quality synthetic aperture radar interferograms and amplitude‐offset images, the nearly linear surface rupture is well traced along the western end of the Nam Ma fault and strikes ∼69°. From both descending and ascending pass Interferometric Synthetic Aperture Radar data and a rigorous maximum a posteriori probabilistic inversion method, we infer that the event involved mostly a pure left‐lateral strike‐slip rupture with near‐vertical geometry. Our one‐segment model shows that the maximum slip of ∼4.1 m occurred at ∼4 km depth, much larger than the slip at the surface. Both interferograms also reveal a small segment to the east of the main rupture, in a densely populated farming area. Our inversion of a two‐segment model shows a similar slip distribution on the main fault, in addition to ∼0.1–0.3 m left‐lateral slip with normal component on a 58° north‐dipping segment. The total seismic moment from the two‐segment model is 1.95×10 19 N·m, equivalent to an M w 6.79 earthquake, which is comparable to the U.S. Geological Survey seismic inversion estimate of 2.10×10 19 N·m ( M w 6.84). The earthquake occurred within a group of east‐northeast‐striking left‐lateral strike‐slip faults near the Myanmar–Laos border, which are seismically active and reflect a system of actively clockwise rotating blocks. Online Material: Figures of InSAR observations, amplitude‐offset maps, InSAR deformation decomposition, one segment model inversion residuals, single data set inversion results, and tectonic setting.