Qijie Wang

InSAR analysis of surface deformation over permafrost to estimate active layer thickness based on one-dimensional heat transfer model of soils.

This paper presents a novel method to estimate active layer thickness (ALT) over permafrost based on InSAR (Interferometric Synthetic Aperture Radar) observation and the heat transfer model of soils. The time lags between the periodic feature of InSAR-observed surface deformation over permafrost and the meteorologically recorded temperatures are assumed to be the time intervals that the temperature maximum to diffuse from the ground surface downward to the bottom of the active layer. By exploiting the time lags and the one-dimensional heat transfer model of soils, we estimate the ALTs. Using the frozen soil region in southern Qinghai-Tibet Plateau (QTP) as examples, we provided a conceptual demonstration of the estimation of the InSAR pixel-wise ALTs. In the case study, the ALTs are ranging from 1.02 to 3.14u2009m and with an average of 1.95u2009m. The results are compatible with those sparse ALT observations/estimations by traditional methods, while with extraordinary high spatial resolution at pixel level (~40u2009meter). The presented method is simple, and can potentially be used for deriving high-resolution ALTs in other remote areas similar to QTP, where only sparse observations are available now.

A Refined Strategy for Removing Composite Errors of SAR Interferogram

In standard differential synthetic aperture radar interferometry, there could still be a residual tilt (orbital error) in the interferometric phase due to inaccurate baseline estimation. We demonstrated theoretically that the orbital errors were partially elevation dependent. On the basis of this, we introduced an elevation-dependent item to the conventional polynomial model to simulate, and therefore, compensate the orbital errors, as well as the small scale topographic and/or topography-related phase errors. Robust regression approach was suggested to determine the parameters of the proposed model. The model was validated with both synthetic and real ALOS PALSAR data of the Zhouqu, China mudslide. The synthetic test indicated that upon applying the refined model, the accuracies of phase measurements were improved by nearly two times, compared to those using conventional linear and quadratic models. The real data experiment indicated that after utilizing the refined model, the correlation between the interferogram and the digital elevation model of Zhouqu reduced to about 1/5 of those using linear and quadratic models. This demonstrates that the elevation-dependent phase components have been largely removed by the new model. More importantly, the interferogram corrected by the new model visibly disclosed the deformation area affected by the Zhouqu mudslide.

Investigating the Ground Deformation and Source Model of the Yangbajing Geothermal Field in Tibet, China with the WLS InSAR Technique

Ground deformation contains important information that can be exploited to look into the dynamics of a geothermal system. In recent years, InSAR has manifested its strong power in the monitoring of ground deformation. In this paper, a multi-temporal InSAR algorithm, WLS InSAR, is employed to monitor and characterize the Yangbajing geothermal field in Tibet, China, using 51 ENVISAT/ASAR images acquired from two overlapping descending tracks. The results reveal that the WLS InSAR algorithm can suppress the adverse effects of seasonal oscillations, associated with the freezing-thawing cycle of the permafrost in the Qinghai-Tibet Plateau. Deformations of up to 2 cm/yr resulting from the exploitation of the geothermal resource have been detected in the southern part of the Yangbajing field between 2006 and 2010. A source model inversion of the subsurface geothermal fluids was carried out based on the elastic half-space theory using the accumulated deformations. It was found that most geothermal fluid loss has occurred in the southern part of the shallow reservoir as the pore space beneath the northern part of field was recharged by the ascending flow from the deep layers of the reservoir through well-developed faults in the region.

Coastal Subsidence Monitoring Associated with Land Reclamation Using the Point Target Based SBAS-InSAR Method: A Case Study of Shenzhen, China

Shenzhen, the first special economic zone of China, has witnessed earth-shaking changes since the late 1980s. In the past 35 years, about 80 km2 of land has been reclaimed from the sea in Shenzhen. In order to investigate coastal vertical land motions associated with land reclamation, we proposed an elaborated Point Target (PT) based Small Baseline Subset InSAR (SBAS-InSAR) strategy to process an ENVISAT ASAR ascending and descending orbits dataset both acquired from 2007 to 2010. This new strategy can not only select high density PTs but also generate a reliable InSAR measurement with full spatial resolution. The inter-comparison between InSAR-derived deformation velocities from different orbits shows a good self-consistency of the results extracted by the elaborated PT-based SBAS-InSAR strategy. The InSAR measurements show that the reclaimed land is undergoing remarkable coastal subsidence (up to 25 mm/year), especially at the Shenzhen Airport, Bao’an Center, Qianhai Bay and Shenzhen Bay. By analyzing the results together with land reclamation evolution, we conclude that the ground deformation is expected to continue in the near future, which will amplify the regional sea level rise.

Continent-Wide 2-D Co-Seismic Deformation of the 2015 Mw 8.3 Illapel, Chile Earthquake Derived from Sentinel-1A Data: Correction of Azimuth Co-Registration Error

In this study, we mapped the co-seismic deformation of the 2015 Mw 8.3 Illapel, Chile earthquake with multiple Sentinel-1A TOPS data frames from both ascending and descending geometries. To meet the requirement of very high co-registration precision, an improved spectral diversity method was proposed to correct the co-registration slope error in the azimuth direction induced by multiple Sentinel-1A TOPS data frames. All phase jumps that appear in the conventional processing method have been corrected after applying the proposed method. The 2D deformation fields in the east-west and vertical directions are also resolved by combing D-InSAR and Offset Tracking measurements. The results reveal that the east-west component dominated the 2D displacement, where up to 2 m displacement towards the west was measured in the coastal area. Vertical deformations ranging between −0.25 and 0.25 m were found. The 2D displacements imply the collision of the Nazca plate squeezed the coast, which shows good accordance with the geological background of the region.

Landslide Displacement Monitoring by a Fully Polarimetric SAR Offset Tracking Method

Landslide monitoring is important for geological disaster prevention, where Synthetic Aperture Radar (SAR) images have been widely used. Compared with the Interferometric SAR (InSAR) technique, intensity-based offset tracking methods (e.g., Normalized Cross-Correlation method) can overcome the limitation of InSAR’s maximum detectable displacement. The normalized cross-correlation (NCC) method, based on single-channel SAR images, estimates azimuth and range displacement by using statistical correlation between the matching windows of two SAR images. However, the matching windows—especially for the boundary area of landslide—always contain pixels with different moving characteristics, affecting the precision of displacement estimation. Based on the advantages of polarimetric scattering properties, this paper proposes a fully polarimetric SAR (PolSAR) offset tracking method for improvement of the precision of landslide displacement estimation. The proposed method uses the normalized inner product (NIP) of the two temporal PolSAR Pauli scattering vectors to evaluate their similarity, then retrieve the surface displacement of the Slumgullion landslide located in southwestern Colorado, USA. A pair of L-band fully polarimetric SAR images acquired by the Jet Propulsion Laboratory’s Uninhabited Aerial Vehicle Synthetic Aperture Radar (UAVSAR) system are selected for experiment. The results show that the Slumgullion landslide’s moving velocity during the monitoring time ranges between 1.6–10.9 mm/d, with an average velocity of 6.3 mm/d. Compared with the classical NCC method, results of the proposed method present better performance in the sub-pixel estimation. Furthermore, it performs better when estimating displacement in the area around the landslide boundaries.

Retrieving three-dimensional coseismic displacements of the 2008 Gaize, Tibet earthquake from multi-path interferometric phase analysis

In this paper, synthetic aperture radar (SAR) data from ENVISAT ASAR ascending, descending and ALOS PALSAR ascending orbits are collected to investigate the coseismic displacements of the Mw 6.4 earthquake occurred in Gaize, Tibet on January 9, 2008 and the Mw 5.9 aftershock on January 16, 2008. Two interferometric phase analysis techniques, i.e., D-InSAR and multi-aperture InSAR, are employed to process the SAR data, with which the displacement measurements along three different line-of-sight (LOS) and three different azimuth directions are retrieved, respectively. Complete three-dimensional (3-D) coseismic displacement fields caused by the earthquake are then resolved by integrating the obtained LOS and azimuth displacement measurements with a weighted least squares adjustment, whose distributions are conformed to the two north-northeast trending northwest-dipping normal faults detected in previous studies. Ground subsidence and uplift are observed in the hanging wall and footwall of the main fault, respectively, and the subsidence reaches its maximum in the hanging wall of the second fault as a superimposed result of the Gaize earthquake and its aftershock. Anti-symmetric horizontal movements are also detected during the seismic events, which move inward in the focal region, but outward at the marginal. The left-lateral motions near the main fault indicate a small striking slip component caused by the Gaize earthquake. Finally, we discuss the potential of applying the derived spatially continuous 3-D displacement fields to determine the high-resolution 3-D strain fields of the Gaize earthquake, which provide important knowledge for assessing the source mechanism.

An Adaptive Offset Tracking Method with SAR Images for Landslide Displacement Monitoring

With the development of high-resolution Synthetic Aperture Radar (SAR) systems, researchers are increasingly paying attention to the application of SAR offset tracking methods in ground deformation estimation. The traditional normalized cross correlation (NCC) tracking method is based on regular matching windows. For areas with different moving characteristics, especially the landslide boundary areas, the NCC method will produce incorrect results. This is because in landslide boundary areas, the pixels of the regular matching window include two or more types of moving characteristics: some pixels with large displacement, and others with small or no displacement. These two kinds of pixels are uncorrelated, which result in inaccurate estimations. This paper proposes a new offset tracking method with SAR images based on the adaptive matching window to improve the accuracy of landslide displacement estimation. The proposed method generates an adaptive matching window that only contains pixels with similar moving characteristics. Three SAR images acquired by the Jet Propulsion Laboratory’s Uninhabited Aerial Vehicle Synthetic Aperture Radar (UAVSAR) system are selected to estimate the surface deformation of the Slumgullion landslide located in the southwestern Colorado, USA. The results show that the proposed method has higher accuracy than the traditional NCC method, especially in landslide boundary areas. Furthermore, it can obtain more detailed displacement information in landslide boundary areas.

Generalized functional model of maximum and minimum detectable deformation gradient for PALSAR interferometry

Abstract Empirical functional models for the maximum and minimum detectable deformation gradient of PALSAR interferometry were established based on coherence and discrete look numbers. Then, a least square regression method was used to fit the model coefficients and thus obtain the generalized functional models for both coherence and look numbers. The experimental results with ALOS PALSAR data of Wenchuan earthquake of China show that the new model works well for judging whether the deformation gradient can be detected by the D-InSAR technology or not. The results can help researchers to choose PALSAR data and to configure processing parameters, and also benefit the interpretation of the measured surface deformation.

Deriving time-series three-dimensional displacements of mining areas from a single-geometry InSAR dataset

This paper presents a method for deriving time-series three-dimensional (3-D) displacements of mining areas from a single-geometry interferometric synthetic aperture radar (InSAR) dataset (hereafter referred to as the SGI-based method). This is mainly aimed at overcoming the limitation of the traditional multi-temporal InSAR techniques that require SAR data from at least three significantly different imaging geometries to fully retrieve time-series 3-D displacements of mining areas. The SGI-based method first generates the multi-temporal observations of the mining-induced vertical subsidence from the single-geometry InSAR data, using a previously developed method for retrieving 3-D mining-related displacements from a single InSAR pair (SIP). The weighted least-squares solutions of the time series of vertical subsidence are estimated from these generated multi-temporal observations of vertical subsidence. Finally, the time series of horizontal motions in the east and north directions are estimated using the proportional relationship between the horizontal motion and the subsidence gradient of the mining area, on the basis of the SGI-derived time series of vertical subsidence. Seven ascending ALOS PALSAR images from the Datong mining area of China were used to test the proposed SGI-based method. The results suggest that the SGI-based method is effective. The SGI-based method not only extends the SIP-based method to time-series 3-D displacement retrieval from a single-geometry InSAR dataset, but also limits the uncertainty propagation from InSAR-derived deformation to the estimated 3-D displacements.