Monitoring deformation in rapidly subsiding funnels within mining areas remains a major challenge in the InSAR community. Most existing studies have primarily focused on the identification, filtering, and unwrapping of such subsidence funnels, while overlooking the issue of coherence underestimation in regions of rapid subsidence. This neglect often results in a sparse distribution of measurement points and difficulties in satisfying the continuity assumptions required for reliable phase unwrapping. To address this, this article presents a coherence-maximization algorithm to compensate for the underestimation bias caused by nonlinear phase gradients contributed from rapid subsidence in mining funnels. The algorithm begins with automatic detection of small objects by integrating the YOLOv11 model with the Slicing Aided Hyper Inference (SAHI) framework. Subsequently, each detected subsidence funnel is individually modeled from the wrapped interferometric phases using a two-dimensional mixed Gaussian model, optimized by a hybrid Simulated Annealing-Genetic Algorithm. The objective function is designed to maximize the interferometric coherence. Using both semi-simulated and real data from Sentinel-1 and LuTan-1 over Shanxi Province, China, we validated the effectiveness of the coherence-maximization algorithm under varying levels of decorrelation. The observed improvements in phase unwrapping further confirm the potential value of the proposed method for enhancing the reliability of InSAR measurements over mining-affected scenarios.

Phase linking is the crucial step in distributed scatterer InSAR processing, which determines the quality of time-series interferometric phases. The weight matrix is the key measure of phase linking method, which controls the participation of each interferometric pair for the single-master phase linking. Existing methods don’t consider the impacts of temporal-changed polarimetric scattering characteristics, leading to large phase closure errors. Based on the polarimetric stationarity, we propose a novel scattering consistency weight measure. Combined with the coherence weight, a joint weight is generated to improve three general phase linking methods. The methods are validated with time-series Radarsat-2 PolSAR data over Kilauea Volcano, Hawaii. The results show that the proposed methods obtain higher temporal posterior coherence and better equivalent single-master (ESM) interferometric phase than traditional methods.

DSM block adjustment is a crucial step to improve absolute height accuracy and consistency within overlap in large-scale topographic mapping using bistatic interferometric SAR. More objects, fixed error models, and more observations types with varying accuracies pose significant challenges for adjustment. To this end, this letter proposes an integration of iterative method, hypothesis testing and surveying adjustment into the solution. Preconditioning technology and -test are combined for adaptive correcting of the error model, while robust Helmert variance component estimation is employed for dynamic update of the stochastic model. Large-scale spaceborne bistatic SAR data containing diverse land cover types are utilized to validate the method. The 90% linear error (LE90) of DSM height against ICESat-2 decreases from 26.27 m to 4.90 m, corresponding to a decrease from 27.98 m to 3.70 m in LE90 of height difference at tie-points. The results demonstrate not only the ability to obtain the most suitable height error model for each DSM scene, but also the determination of appropriate weighting ratio between control point and tie-point while resisting gross error.

An advanced InSAR-based iterative tropospheric decomposition model for more effectively mitigating both vertically stratified and turbulent delays is proposed in this paper. The novelty of the model lies in two aspects: 1) incorporating turbulence variance, which characterizes the statistical properties of atmospheric turbulence, into a stochastic model when iteratively removing stratified delay, and 2) taking the distance and turbulence variance as the weight factors of samples, simultaneously, in spatially interpolating turbulence. In addition, our model is implemented interferogram by interferogram, so we do not need time-series InSAR datasets. The effectiveness of the proposed method is tested and validated using a set of simulated experiments and an in-site experiment over Hawaii Island, and the experimental results show that our method can mitigate tropospheric delays efficiently and robustly.

We propose an atmospheric turbulence statistical characteristic-enhanced stacking-interferometry synthetic aperture radar (InSAR) method to more accurately derive the average deformation velocity for slow or linear surface displacements. It can be used to estimate the variance and covariance matrix (VCM) for time-series InSAR observations pixel by pixel. Moreover, a minimum variance linear estimator is used to calculate optimal weights for interferogram stacking. Compared with traditional stacking, our method is an improvement, and its advantages include the following: 1) the relevance between time-series atmospheric delays of interferograms is considered; 2) the atmospheric turbulence is better suppressed; and 3) the performance is not limited by the quantity of interferograms. The effectiveness of the new method is verified through a series of simulated experiments and real Sentinel-1 experiments over Southern California. The results demonstrate that our proposed method is more effective and robust in removing turbulent atmospheres.

A novel distributed scatterer interferometric synthetic aperture radar (DS-InSAR) method is presented in which the sum of Kronecker product (SKP) decomposition method is applied to DS candidates. Unlike existing polarimetric optimization methods, the proposed method considers polarimetric and interferometric coherence information simultaneously, resulting in separation of the polarimetric scattering process for each target and the maximum diversity for the corresponding phase center locations. Physical reliability of the phase optimization solution is thereby ensured. The performance of the novel method is evaluated using 30 quad-polarized C-band Radarsat-2 synthetic aperture radar (SAR) images over Kilauea Volcano, Hawaii. The proposed method provides a higher density of measurement scatterer (MS) points and a higher quality of DS interferometric phase with a temporally stable phase center than single-polarization (HH) method and quad-polarization exhaustive search polarimetric optimization (ESPO) method. Thus, the proposed method shows good performance in phase quality improvement and point density increasement.

Homogeneous pixel selection (HPS) plays an important role in the application of multitemporal SAR interferometry. The statistical goodness-of-fit testing of the temporal samples has been widely used for HPS. However, the detection rates of the existing methods are unsatisfactory under small datasets. In this paper, a stacked auto-encoder (SAE) network based method is proposed for the selection of homogeneous pixels under the idea of deep learning image classification, as termed by SAEHPS. The SAE network is used to learn the spatial distribution behavior of the average intensity image. The deep network is trained and tested on different high-resolution SAR datasets of the Hong Kong Airport and the Fuzhou City, and three pixel-wise labels (i.e., high, medium, and low reflections) are regarded as outputs of model learning. The unsupervised training and supervised fine-tuning realize the class prediction. The results show that the SAE can achieve robust accuracies above 90% based on empirically labeled samples, especially in non-architectural areas where the distributed scatterers exist. The SAE results are devoted to the multitemporal PS/DS InSAR approach to identify homogeneous pixels. Both qualitative and quantitative experiments in HPS, phase optimization, and deformation monitoring have demonstrated the superiority of the novel method.

With the increasing amount of synthetic aperture radar (SAR) data with various imaging geometries (at least ascending/descending tracks), it is possible to obtain accurate multi-dimensional (MD) deformation time series with long time span. However, in most cases SAR data of different geometries are un-synchronously acquired over the same region, making it impossible to directly solve the underdetermined observation model (OSM) between the interferometric SAR (InSAR) measurements and the MD deformations. Kalman filter (KF), as one of the most famous dynamic estimators, can obtaina prioriinformation of the unknowns based on the preexisting time series, therefore it can be used to deal with this InSAR underdetermined problem. This article employs the KF to realize the dynamic estimation of MD deformations with short-baseline interferograms. The innovation lies in the establishment of the KF state transition model (STM) and OSM, which aims to make the InSAR monitoring problem better adapt to the KF. Particularly, by assuming a smooth deforming process, existing deformation time series are used to establish the STM and to predict the deformations at current moment. Besides, a strain model (SM) is employed to assist the establishment of the OSM. Simulation and real experiments in the Geysers geothermal field (GGF), U.S. demonstrate that, compared with the state-of-the-art methods, the proposed KF method allows more robust deformation estimation and achieves higher computational efficiency for dynamic estimation.

Accurately predicting ground surface deformation is a crucial task in mining-related geohazards control. Interferometric synthetic aperture radar (InSAR) technique can widely detect historical line-of-sight (LOS) displacements with a high spatial resolution. By incorporating with spatio-temporal deformation models, InSAR can predict kinematic 3-D displacements due to underground mining. However, this method depends on the geometric parameters (at least seven generally) of underground mined-out areas, hindering its practical applications especially over a large area. To circumvent this, we proposed a new method for predicting kinematic 3-D mining displacements by incorporating InSAR with a temporal model, rather than spatio-temporal models used before, in this article. In doing so, much less prior parameters (only three and can be empirically given) are required, with respect to the previous InSAR-based methods. To achieve this, we first revealed that InSAR LOS displacements caused by underground longwall mining at a single point temporally follow an S-shaped growth pattern. Meanwhile, we also showed that a Weibull model can describe the temporal evolution well. Based on these findings, the proposed method first predicts, in a point-wise manner, the kinematic LOS mining displacements from historical InSAR measurements using the Weibull model and a Kalman filter. The kinematic 3-D displacements are then resolved from the predicted LOS displacements with the help of a common prior information relating to mining deformation. The proposed method was tested in the Datong coal mining area of north China. The results show an averaged accuracy of about 0.007 m of the resolved kinematic 3-D displacements.

It is a common method to resolve three-dimensional (3-D) deformation components associated with underground mining by incorporating Single-track interferometric synthetic aperture radar (InSAR) with a prior deformation model termed linear proportion model (LPM) (hereinafter referred to as Sin-LPM). Nevertheless, the Sin-LPM method relies on three model parameters that are needed to bein situcollected, and it neglects their dynamic changes during the period of underground extraction, narrowing the practical applications of the Sin-LPM method, and degrading the accuracy of the estimated 3-D displacements. In this article we propose a new method to resolve 3-D mining displacements from multi-track InSAR observations by incorporating with the LPM. In which, the model parameters are first considered as dynamic and further adaptively estimated from the multi-track InSAR observations using a robust solver. Following that, 3-D mining displacements are resolved from the multi-track InSAR using the conjugate gradient method (CGM). The proposed method was tested in Datong coalfield, China. The results suggest that the proposed method can well estimate 3-D mining displacements with a mean error of about 1.8 cm. Compared with the previous Sin-LPM, the proposed method can effectively improve the accuracy of the estimated 3-D displacements (e.g., 69% in this study), and can work well even over a large area where the model parameters are unknown. The proposed method offers a new insight to improve the InSAR-based retrieval of 3-D displacements induced by other anthropologic or geophysical activities.

In this letter, we propose a forest height inversion model based on three-order Fourier-Legendre (FL) polynomials from the multibaseline polarimetric synthetic aperture radar interferometry (PolInSAR) data. The proposed model expresses the vertical structure of the volume layer as three-order FL polynomials. Meanwhile, the forest height is treated as an unknown parameter, rather than a priori information, as adopted in polarization coherence tomography technology. On the other hand, to be more realistic, the proposed model uses multipolarization PolInSAR data and considers that the synthetic aperture radar (SAR) signals in different polarizations describe the forest vertical structure in different ways so that we can obtain a more comprehensive forest vertical structure. Airborne P-band PolInSAR data acquired over the boreal and tropical forest areas were selected for testing the forest height inversion method. The results show that, compared to random volume over ground (RVoG) model-based inversion, the accuracy of the proposed model is improved by 28.20% and 17.30%, respectively, for the boreal and tropical forest scenes.

Deformation monitoring of a converter station of an electrical power transmission system will help to prevent potential damages to power facilities and properties. Corner reflector InSAR (CRInSAR) enables local deformation measurements in low coherence areas like construction-engineering projects. However, the accuracy of CRInSAR will be degraded by the phase unwrapping errors, especially when the deformation is large or dominated by nonlinear component. This letter reports a study that employs CRInSAR technique to monitor the deformations of ten CRs installed in Xiangtan converter station, China, with seven TerraSAR-X Spotlight images. A two-step phase unwrapping tactics is proposed based on the least squares ambiguity decorrelation adjustment (LAMBDA) algorithm to focus on estimating nonlinear deformation without being affected by unwrapping errors. The results reveal that the two-step LAMBDA algorithm can achieve an accuracy of less than 2 mm for CRInSAR deformation monitoring regardless of small- or large-scale deformations, as validated by the trigonometric leveling measurements.

In this paper, we demonstrate that the spaceborne dual-platform TerraSAR-X (TSX) and TanDEM-X (TDX) pursuit monostatic mode full polarimetric (full-pol) synthetic aperture radar (SAR) data with a time lag can be used to monitor maritime traffic. For single polarization (single-pol) SAR data, the performance of vessel velocity estimation is mainly determined by 2-D cross correlation of SAR intensity data. As the sea clutter is changing dynamically during the TSX/TDX data acquisition, the correlation between two dual-platform images decreases significantly. We may get unstable or incorrect estimations of vessel velocity, especially under a higher wind condition. For solving this problem, we propose an object-oriented polarimetric likelihood ratio test (PolLRT) method based on the complex Wishart distribution. The proposed method makes PolLRT statistics of the detected target pixels for eliminating the effect of varied sea clutter. Two pairs of full-pol SAR data sets covering the Strait of Gibraltar acquired by dual-platform TSX/TDX in pursuit monostatic mode with a time lag of approximately 10 s are selected for the experiments. The experimental results demonstrate that the proposed PolLRT method has a better performance than that of the classical normalized cross correlation (NCC) method with VV polarization SAR data and the mutual information (MI) method with full-pol SAR data. Specifically, under the lower wind condition, the correct estimation rate of the NCC, the MI, and the proposed PolLRT methods are 85.7%, 57.1%, and 100%, respectively; under the relatively higher wind condition, the correct estimation rate of the above three methods are 48.8%, 23.2%, and 90.1%, respectively.

Efficient and accurate simulation of ground penetrating radar (GPR) in the open region helps immensely in both grasping the features of echoes and facilitating the interpretation of real GPR data. Due to the limitation of the computer model, however, the strong artificial boundary reflections, especially the low-frequency propagating waves encountered at the late stage of simulation greatly affect the simulation accuracy of GPR. This paper presents an innovative optimized unsplit-field convolutional perfectly matched layer (UCPML) based on rotated staggered grid (RSG) scheme to truncate the finite-difference time-domain (FDTD) lattices. Rather than obey the sharp variation based on an m th-order polynomial, the optimized approach employs a novel optimized term and an adjustment factor to seek a gentle variation on optimal constitutive coefficients. This guarantees that the determination of optimal constitutive coefficients can be less influenced by the order of polynomial and especially, to improve the absorptive performance on low-frequency propagating waves. The calculating efficiency and accuracy of the RSG-FDTD scheme, as well as the absorbing performance of the optimized UCPML, are verified by two numerical examples. In particular, the analysis of the amplitude-frequency features of low-frequency clutters at steady state of the electromagnetic (EM) field and the corresponding global reflection error in the time-frequency domain is also presented.

Interferometric synthetic aperture radar (InSAR) technique is a proven technique for measuring 3-D surface deformations by combining InSAR measurements from different techniques (i.e., differential InSAR, multiaperture InSAR, and pixel offset-tracking) and different tracks (i.e., ascending and descending) on a pixel-by-pixel basis. However, it is difficult to obtain the exact a priori variances or weights for such different kinds of InSAR measurements, resulting in inaccurate estimations of 3-D deformations. This paper proposes a method to retrieve 3-D deformations with InSAR by integrating the strain model and variance component estimation algorithm, which can exploit the spatial correlation of the adjacent points' deformations and produce accurate weights for multiple InSAR measurements. The proposed method is assessed with both simulated and real data sets. The results have shown that the proposed method can accurately measure 3-D surface deformations associated with geohazards, and even those occurring in a transient or short-term period (e.g., earthquake and volcanic eruption). In the case study of the 2007 eruption of Kilauea Volcano (Hawai'i), improvements of 51.2%, 22.4%, and 18.5% have been achieved for the derived east, north, and up displacements, respectively, with respect to those derived from the classical weighted least squares method.

This paper presents a novel method for modeling and robustly estimating the time-series 3-D mining-induced large displacements from a single imaging geometry (SIG) synthetic aperture radar (SAR) amplitude data set using the offset-tracking (OT) technique (hereafter referred to as the OT-SIG). It first generates multitemporal observations of 3-D mining-induced displacements from the single-geometry SAR amplitude data set with the assistance of a prior model. Then, a functional relationship between mining-induced time-series 3-D displacements and the multitemporal 3-D deformation observations generated is constructed. Finally, the time-series 3-D displacements are robustly estimated based on the constructed function model using the M-estimator. The proposed OT-SIG provides a robust and cost-effective tool for retrieving time-series 3-D mining-induced large displacements, relaxing the basic requirement of the traditional method that at least two different viewing geometries' SAR data are needed. Finally, we tested the proposed OT-SIG with descending TerraSAR-X SAR amplitude data set over the Daliuta coal mining area in China. The results show that the root-mean-square errors (RMSEs) of OT-SIG-estimated time-series displacements are about 0.22 and 0.11 m in the vertical and horizontal directions, respectively. These RMSEs are around 5.7% and 10.9% of the maximum in situ deformation measurements in the corresponding directions, which can meet the accuracy requirements of practical applications.

Measuring 3-D mining-induced displacements is essential to understand mining deformation mechanisms and assess mining-related geohazards. In our previous work, we proposed a method for estimating 3-D mining-induced large displacements with the surface deformation along the radar line-of-sight (LOS) direction derived from a single amplitude pair (SAP) of synthetic aperture radar (SAR) using the offset tracking (OT) procedure (hereafter referred to as OT-SAP). The OT-SAP method effectively reduces the strict requirements on SAR data of the previous OT-based methods for 3-D mining-induced displacement retrieval. However, OT-SAP is not robust to errors in the LOS deformation, due to the lack of redundant observations. In this paper, we present an alternative approach (hereafter called AOT-SAP) to OT-SAP. The AOT-SAP method involves estimating the 3-D mining-induced large displacements with OT-derived 2-D deformation observations along the LOS and azimuth directions from an SAP of SAR, instead of just the LOS deformation in the OT-SAP method. Consequently, more redundant observations are incorporated in the AOT-SAP method compared with the previous OT-SAP method. The theoretical analysis and experiments based on both simulated and real data sets suggest that AOT-SAP can effectively improve the accuracies of the estimated 3-D displacements compared with the OT-SAP-estimated ones.

A new approach is presented for mapping 3-D surface displacement caused by subsurface fluid volumetric change based on 1-D interferometric synthetic aperture radar (InSAR) line-of-sight measurements and surface deformation modeling. The relationship between surface deformation and source fluid volumetric change is modeled according to elastic half-space theory. A distinctive advantage of the proposed approach is that it effectively extends the capability of the sun-synchronous orbit side-looking synthetic aperture radar that has been essentially only able to measure 1-D displacements accurately or at most 2-D displacements when InSAR measurements from more than one orbit or platform are combined. Experimental studies are carried out with both simulated and real data sets to test the performance of the method. The results have demonstrated that the approach works very well.

Underground extraction can be roughly classified into three types, i.e., subcritical, critical, and supercritical extraction, in accordance with the geological conditions in the overburden and the geometry of mined-out areas. In 2016, we proposed an approach based on the interferometric synthetic aperture radar (InSAR) technique and the probability integral method (PIM) for the cost-effective prediction of 3-D mining-induced displacements (abbreviated as InSAR-PIM). Due to the inherent assumption of critical extraction in the PIM, the InSAR-PIM method performs well in predicting the 3-D displacements caused by critical and/or supercritical extraction, but poorly for subcritical extraction. In this paper, we first propose a generalized PIM (GPIM) by modifying the traditional PIM with a simplified Boltzmann function. We then replace the PIM of the InSAR-PIM with the proposed GPIM to develop an extension of InSAR-PIM (referred as to InSAR-GPIM). The InSAR-GPIM was tested in the Qianyingzi coal mining area, China. The results show that the InSAR-GPIM-predicted horizontal and vertical displacements caused by subcritical, critical, and supercritical extraction agree well with the in situ observations, with average root-mean-square errors of about 0.032 and 0.050 m, respectively. These accuracies represent improvements of 60.9% and 59% when compared with the accuracies predicted by the InSAR-PIM in the horizontal and vertical directions. The results indicate that the InSAR-GPIM is capable of accurately predicting 3-D mining-induced displacements under different extraction conditions (i.e., subcritical, critical, and supercritical extraction), and it performs much better than the InSAR-PIM in the case of subcritical extraction. It is therefore believed that InSAR-GPIM will have a wider scope of applications than the previous InSAR-PIM.

This paper presents a novel method for estimating the model parameters of the probability integral method (PIM) based on the line-of-sight deformation derived from the interferometric synthetic aperture radar. Then, it applies the settled PIM to forward predict the horizontal and vertical displacements induced by the extraction of a new working panel. The proposed method first constructed the functional relationship between the InSAR-derived LOS deformation and the model parameters of PIM. Subsequently, an improved genetic algorithm (GA), in which gross error elimination was imposed, was proposed, and used to estimate the model parameters of PIM with a large number of LOS deformation measurements. The estimated model parameters and PIM were then employed to forward predict the horizontal and vertical displacements induced by the extraction of a working panel. Simulated experiments show that the rmses of the predicted displacements along the up-down, west-east, and north-south directions are 1.5, 0.9, and 2.5 mm, respectively. Real data experiments over the Qianyingzi coal mining area of China indicate that the predicted displacements are highly consistent with those by field surveys, with rmses of 4.1 and 3 cm for the vertical and horizontal directions, respectively. These imply that the proposed approach can be a very promising tool for predicting the mining-induced displacements and will potentially contribute to the assessing and forecasting of possible geological hazards in the mining area.

Image segmentation plays a crucial role in object‐based remote sensing information extraction. This study improves the existing mean shift (MS) algorithm for segmenting high resolution remote sensing imagery by adopting two strategies. First, a pixel‐based, fixed bandwidth and weighted MS algorithm is applied to cluster the image. In this process, the space bandwidth is selected according to the resolution of remote sensing images, and the range bandwidths of each band are calculated based on grey feature and the plug‐in rule. Gaussian kernels are used for clustering. Second, a region‐based MS algorithm is applied to globally merge modes which are obtained in the first step. The spatial and range bandwidths are adaptively adjusted based on the clustering result of the first step. Experimental results with two Quickbird images show that the improved algorithm is superior to the typical MS algorithm, producing high precision and requiring less operation time.

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.

The coherence of radar echoes is a fundamental observable in interferometric synthetic aperture radar (InSAR) measurements. It provides a quantitative measure of the scattering properties of imaged surfaces and therefore is widely applied to study the physical processes of the Earth. However, unfortunately, the estimated coherence values are often biased due to various reasons such as radar signal nonstationarity and the bias in the estimators used. In this paper, we focus on multitemporal InSAR coherence estimation and present a hybrid approach that mitigates effectively the errors in the estimation. The proposed approach is almost completely self-adaptive and workable for both Gaussian and non-Gaussian SAR scenes. Moreover, the bias of the sample coherence can be mitigated with even only several samples included for a given pixel. Therefore, it is a more pragmatic method for accurate coherence estimation and can be applied actually. Different data sets are used to test the proposed method and demonstrate its advantages.

A novel coherence estimation method for small data sets is presented for interferometric synthetic aperture radar (SAR) (InSAR) data processing and geoscience applications. The method selects homogeneous pixels in both the spatial and temporal spaces by means of local and nonlocal adaptive techniques. Reliable coherence estimation is carried out by using such pixels and by correcting the bias in the estimated coherence caused by the non-Gaussianity in high-resolution SAR scenes. As an example, the proposed method together with coherence decomposition is applied to extract the temporal decorrelation component over an area in Macao. The results show that the proposed algorithms work well over various types of land cover. Moreover, the coherence change with time can be more accurately detected compared to other conventional methods.

A Kalman-filter-based approach is presented for resolving 3-D surface displacements using multisensor, multitrack, and multitemporal interferometric synthetic aperture radar (SAR) measurements. Measurements from each interferogram are projected into the three reference directions and combined in the Kalman filter model with displacements determined from previous interferograms to produce updated displacement measurements. Both simulated and real data sets are used to test the proposed approach. It is found that the method works well when the measurement noise is low. The displacements in the north direction, however, are much lower in accuracy than those in the other two directions and even become unstable when the measurement noise is high due to the polar-orbiting imaging geometries of the current satellite SAR sensors.

We map the coseismic deformation of the 2011 Tohoku-Oki earthquake with data from three descending Envisat/ASAR tracks and six ascending ALOS/PALSAR tracks that cover most of northeastern Japan. Due to the inaccurate estimation of the satellite status, orbital ramps commonly exist in the coseismic interferograms, which resulted in inconsistency among the deformation maps released by several research groups. In this letter, calibration has been performed to accurately remove these ramps by a 2-D quadratic-phase model derived based on GPS measurements from the ARIA team at the Jet Propulsion Laboratory and Caltech. The average RMS of the interferometric synthetic aperture radar (InSAR) measurements, as compared with GPS measurements at the validation stations, has decreased from 17.8 to 7.7 cm after the orbital ramp correction is made, indicating that much more accurate InSAR measurements are achieved. The corrected coseismic deformation from the InSAR measurements is consistent with not only the GPS observations at the individual GPS stations but also with the coseismic deformation interferogram from interpolated GPS observation in the SAR viewing directions. The corrected coseismic deformation measurement results show a maximum line-of-sight displacement of up to 3.7 m from the ascending PALSAR tracks and 2.4 m from the descending ASAR tracks, respectively.

Previous approaches that integrate interferometric synthetic aperture radar (InSAR) and GPS measurements for 3-D surface displacement mapping require statistically estimating the variances of the measurements to yield optimal results. We present a variance component estimation approach to weigh the InSAR and GPS measurements in deriving 3-D surface displacements. The approach exploits the observations themselves for determining the weighting scheme, and therefore the a priori information on the stochastic model of the observations is not required. This is of great importance as accurate knowledge on the stochastic model is often unavailable. The performance of the proposed method is validated with both simulated and real datasets.

The Vondrak filter is a unique technique for smoothing data. The filter aims to achieve a balance between the fidelity and the smoothness of the filtered results. It can therefore preserve the original attributes of the observational data while, at the same time, smooth out the noise. We reformulate the 1-D Vondrak filter that has been widely used in data processing in fields such as astronomy and geophysics and then extend it into two dimensions. The method of conjugate gradients is used to solve the least squares optimization problem. The proposed 2-D filter is a powerful tool for enhancing the quality of various geoscience and remote sensing data such as satellite images. Various tests with simulated and real synthetic aperture radar interferograms show that the new filter is very effective in removing the noise.

Differential interferometric synthetic aperture radar (SAR) (DInSAR) has proven to be very useful in mapping and monitoring land subsidence in many regions of the world. Shanghai, China's largest city, is one of such areas suffering from land subsidence as a result of severe withdrawal of groundwater for different usages. DInSAR application in Shanghai with the C-band European Remote Sensing 1 & 2 (ERS-1/2) SAR data has been difficult mainly due to the problem of decorrelation of InSAR pairs with temporal baselines larger than 10 months. To overcome the coherence loss of C-band InSAR data, we used eight L-band Japanese Earth Resource Satellite (JERS-1) SAR data acquired during 2 October 1992 to 15 July 1998 to study land subsidence phenomenon in Shanghai. Three of the images were used to produce two separate digital elevation models (DEMs) of the study area to remove topographic fringes from the interferograms used for subsidence mapping. Six interferograms were used to generate 2 different time series of deformation maps over Shanghai. The cumulative subsidence map generated from each of the time series is in agreement with the land subsidence measurements of Shanghai city from 1990 - 1998, produced from other survey methods.

2005 IEEE International Geoscience and Remote Sensing Symposium, IGARSS 2005, Seoul, 25-29 July 2005