Repeat-pass Synthetic Aperture Radar Interferometry Research Articles

A High-Precision Baseline Calibration Method Based on Estimation of Azimuth Fringe Frequency with THz Interferometry SAR

In this study, repeat-pass synthetic aperture radar interferometry (repeat-pass THz InSAR) is first extended to the terahertz band, and it has tremendous potential in the application of high-resolution three-dimensional (3D) imaging due to its shorter wavelength, larger bandwidth, and greater sensitivity to elevation variation. The super-resolution and high sensitivity of THz InSAR pose greater demands on the baseline calibration for high-precision digital elevation model (DEM) generation. To meet the elevation accuracy requirement of THz InSAR, we propose a baseline calibration method relying on the estimation of the azimuth fringe frequency (EAFF) of the interferometric phase. Initially, a model for non-parallel sampling path errors within the squint SAR repeat-pass interferometry was established, and then, we conducted the theoretical analysis of the phase errors induced by the non-parallel errors. Following this, using a reference DEM, the relationship between the fringe frequency of the error phase and the bias in the repeat-path positioning was established. This allowed the estimation of the position errors to be transformed into the frequency spectrum estimation based on the FFT, which would mitigate the impact of unknown SAR sampling positions. Ultimately, we investigated the accuracy of the proposed EAFF calibration method, and the simulation showed that it can achieve the theoretical accuracy when the correlation coefficient exceeds 0.3. Furthermore, we configured the repeat-pass THz InSAR system with the 0.3 THz stepped-frequency radar. Compared to the conventional calibration based on ground control points (GCPs), the 3D reconstruction of both a knife and a terrain model, calibrated using the proposed EAFF algorithm, demonstrated that the elevation accuracy can achieve millimeter-level precision across the entire image swath. The above results also proved the great potential of THz InSAR in high-precision 3D imaging and remote sensing.

Assessing the Impact of Positioning Errors in Car-Borne Repeat-Pass SAR Interferometry With a Controlled Rail-Based Experiment

Agile synthetic aperture radar (SAR) platforms such as car-borne and UAV-borne SAR systems require combined inertial navigation systems (INS) and global navigation satellite systems (GNSS) to measure the radar sensor trajectories used for focusing and interferometric processing. Measurement inaccuracies from INS/GNSS systems lead to residual phase errors in the SAR products whose minimisation is crucial to derive accurate topographic and deformation information. In this work, we analyse the impact of residual positioning errors on car-borne repeat-pass SAR interferometry at L-band for different INS/GNSS measurement configurations and for the typical car-borne acquisition geometry. The positioning errors are evaluated both during single SAR acquisitions with long integration times and between different acquisitions as a function of the distance of the radar platform from the GNSS reference stations. We show the reduction of interferometric phase errors achievable by additionally using a GNSS receiver mounted in the vicinity of the SAR platform as compared to remote reference stations of the national network of permanent GNSS receivers. Test results obtained in a controlled setup with a rail-based SAR system equipped with a navigation-grade INS/GNSS system show maximum repeat-pass trajectory errors on the order of 1–2 cm using a local GNSS reference station and up to 10–15 cm using the remote reference stations, leading to azimuth and range phase trends in the interferometric products.

The potential of sentinel-1 InSAR coherence for grasslands monitoring in Eastern Cape, South Africa

This study evaluates the potential of the repeat-pass synthetic aperture radar interferometry (InSAR) coherence for grasslands monitoring in Eastern Cape Province, in South Africa. We used 12 pairs of InSAR data in the Interferometric Wide swath mode (IW), in VV-polarization (with different spatial baseline B⊥ and different temporal repeat intervals between the subsequent acquisitions), covering grasslands in the study area during two growing seasons; drought season (2015/2016) and non-drought season (2016/2017). Results show that non-drought season has higher variability in the coherence between the different regions of interest (ROIs) than the drought season and higher variability in the backscattering coefficient. Based on our results, we conclude that InSAR-based data have shown potential for grasslands monitoring, and the regularity of Sentinel-1 data is well suited for InSAR applications. However, the potential to separate grasslands from other vegetation categories appears to be limited with longer temporal baselines of 12 days intervals using only Sentinel-1A data. This potential will be improved with the 6 days intervals that will become available with Sentinel-1B in the study area.

Vegetated Target Decorrelation in SAR and Interferometry: Models, Simulation, and Performance Evaluation

The paper addresses the temporal stability of distributed targets, particularly referring to vegetation, to evaluate the degradation affecting synthetic aperture radar (SAR) imaging and repeat-pass interferometry, and provide efficient SAR simulation schemes for generating big dataset from wide areas. The models that are mostly adopted in literature are critically reviewed, and aim to study decorrelation in a range of time (from hours to days), of interest for long-term SAR, such as ground-based or geosynchronous, or repeat-pass SAR interferometry. It is shown that none of them explicitly account for a decorrelation occurring in the short-term. An explanation is provided, and a novel temporal decorrelation model is proposed to account for that fast decorrelation. A formal method is developed to evaluate the performance of SAR focusing, and interferometry on a homogenous, stationary scene, in terms of Signal-to-Clutter Ratio (SCR), and interferometric coherence. Finally, an efficient implementation of an SAR simulator capable of handling the realistic case of heterogeneous decorrelation over a wide area is discussed. Examples are given by assuming two geostationary SAR missions in C and X band.

Baseline Estimation Using DEM Prior Knowledge and Capability Analysis for Gaofen-3 Repeat-Pass SAR Interferometry

In this paper, the study on the digital elevation model (DEM) reconstruction using Gaofen-3 satellite repeat-pass synthetic aperture radar (SAR) interferometry is presented. The coherence of the single-look complex (SLC) images and residuals of interferogram provide a powerful guarantee for interferometric processing for Gaofen-3. Baseline estimation is the key factor to generate accurate DEM with SAR interferometry. This paper uses the coarse orbit data and the precise orbit data of Gaofen-3 to carry on the baseline estimation and uses the prior coarse DEM information to refine the baseline estimation. Then, the topographic height is reconstructed with more precise spatial baseline. An experiment is conducted to verify the conclusion, and some mountain and flatten region are chosen as the study area. The results show that the DEM errors based on the prior DEM method are rapidly converged compared with the ones based on the orbit data method. Therefore, benefiting from the precise orbit data of the Gaofen-3 satellite, the accuracy of the baseline estimation can be improved further by using the method based on prior DEM information; thus, higher accuracy of DEM can be achieved so that Gaofen-3 repeat-pass data have some good characteristics, such as high imaging quality, stable spatial baseline, and high correlation, which has a good repeat-pass SAR interferometry application prospect.

Estimation of Residual Motion Errors in Airborne SAR Interferometry Based on Time-Domain Backprojection and Multisquint Techniques

For airborne repeat-pass synthetic aperture radar interferometry (InSAR), precise trajectory information is needed to compensate for deviations of the platform movement from a linear track. Using the trajectory information, motion compensation (MoCo) can be implemented within SAR data focusing. Due to the inaccuracy of current navigation systems, residual motion errors (RMEs) exist between the real and measured trajectory, causing phase undulations in the final interferograms. Up to now, MoCo and RME estimation have usually been combined in airborne InSAR to estimate ground deformation. Conventional MoCo methods generally involve azimuthal and range resampling and phase correction. Then frequency-domain focusing techniques can be used to generate the SAR images. After focusing SAR images with MoCo, both multisquint and autofocus approaches can be used to estimate RME. In addition to the MoCo-based frequency-domain focusing, the time-domain backprojection (BP) technique can also focus the SAR data obtained from highly nonlinear platform trajectories. In this paper, we present, for the first time, the combination of BP and multisquint techniques for RME estimation. A detailed derivation of the implementation of the multisquint approach using the BP-focusing images is presented. Repeat-pass data from the SlimSAR system over Slumgullion landslide are used to demonstrate the feasibility of RME estimation for both stationary and nonstationary scenes. We conclude that the proposed method can effectively remove the RME.

Application of Repeat-Pass TerraSAR-X Staring Spotlight Interferometric Coherence to Monitor Pasture Biophysical Parameters: Limitations and Sensitivity Analysis

This paper describes the potential and limitations of repeat-pass synthetic aperture radar interferometry (InSAR) to retrieve the biophysical parameters of intensively managed pastures. We used a time series of eight acquisitions from the TerraSAR-X Staring Spotlight (TSX-ST) mode. The ST mode is different from conventional Stripmap mode; therefore, we adjusted the Doppler phase correction for interferometric processing. We analyzed the three interferometric pairs with an 11-day temporal baseline, and among these three pairs found only one gives a high coherence. The results show that the high coherence in different paddocks is due to the cutting of the grass in the month of June, however the temporal decorrelation in other paddocks is mainly due to the grass growth and high sensitivity of the X-band SAR signals to the vegetation cover. The InSAR coherence (over coherent paddocks) shows a good correlation with SAR backscatter ( ${R}^{2}_{\rm dB}=0.65, p ) and grassland biophysical parameters ( ${R}^{2}_{\rm Height}=0.55, p , ${R}^{2}_{\rm Biomass}=0.75, p ). It is thus possible to detect different management practices (e.g., grazing, mowing/cutting) using SAR backscatter (dB) and coherence information from high spatial short baseline X-band imagery; however, the rate of decorrelation over vegetated areas is high. Initial findings from the June pair show the possibility of change detection due to the grass growth, grazing, and mowing events by using InSAR coherence information. However, it is not possible to automatically categorize different paddocks undergoing these changes based only on the SAR backscatter and coherence values, due to the ambiguity caused by tall grass flattened by the wind.

Algorithm on the Estimation of Residual Motion Errors in Airborne SAR Images

Due to the lack of accuracy in the navigation system, deviations in the order of centimeters between the real trajectory and the measured one, called residual motion errors (RMEs), frequently appear in synthetic aperture radar (SAR) images. They usually cause azimuth defocusing and phase errors. For very high-resolution SAR imaging and repeat-pass SAR interferometry, the RME must be estimated and compensated. In the literature, the algorithms used to estimate such errors are mainly through interferometry. The error difference for the SAR image pair is calculated, and then the correction is applied either on one image, or alternatively fractions of it on both images. However, the accuracy of this approach is greatly deteriorated by the decorrelation and ground movements. In this paper, we associate the phase difference between two adjacent subaperture images with the second derivative of the RME, and propose a new algorithm to estimate the RME for individual SAR image. It exploits point-like targets distributed along the azimuth direction, and not only corrects the phase error, but also improves the azimuth focusing. Therefore, besides in airborne repeat-pass interferometry, the algorithm can also be applied in very high-resolution SAR imaging. Simulated and real SAR data are used to demonstrate its feasibility and accuracy. Its limitations and extensions are also discussed at the end of the paper.

Ground subsidence phenomena in the Delta municipality region (Northern Greece): Geotechnical modeling and validation with Persistent Scatterer Interferometry

Land subsidence is a common phenomenon occurring in several regions worldwide. Persistent subsidence causes strong consequences on the affected areas and related problems include environmental, economic and social aspects.A set of forty-two Synthetic Aperture Radar (SAR) images, acquired in 1995–2001 by the European Space Agency (ESA) satellites ERS1 and ERS2, were processed with Persistent Scatterer Interferometry (PSI) technique to investigate spatial and temporal patterns of deformation in the Delta municipality (Thessaloniki plain, Northern Greece), a deltaic area with a long history of land subsidence related to aquifer system compaction.Exploitation of output products of a PSI analysis, both average LOS (Line of Sight) deformation rates and displacement time series, revealed a large subsidence area due to intense groundwater withdrawal. Higher displacement velocities have been observed south-west of Kalochori and south of Sindos, from both sides of the Gallikos River. In those areas deformation rates of roughly 4.5 cm/yr have been recorded, during the period from 1995 to 2001.Increasing subsidence rates are measured moving toward the mouth of the Gallikos River, where the thickest sequence of compressible Quaternary sediments is observed. Displacement time series retrieved by PSI technique has been compared with the temporal evolution of the deformation as measured by pre-existing leveling surveys, showing a great agreement.A 2-D finite element model has been run along two representative cross sections in the Kalochori area, in order to simulate the observed temporal evolution of subsidence, coupling the geotechnical behavior of the formations and the piezometric surface level. Finally, results obtained by the subsidence model have been positively compared with the PSI-based information on displacement, providing accurate and perfectly verified results.Outcomes of this work demonstrated the potential of repeat-pass satellite SAR interferometry (InSAR) as suitable technique for increasing knowledge about the extent and the rate of the deformations in case of subsidence events. Moreover, InSAR turned out to be a valuable tool to validate subsidence models and represents a cost-efficiency method, alternative to ground-based measurements for investigating surface deformation phenomena.

Advanced interpretation of land subsidence by validating multi-interferometric SAR data: the case study of the Anthemountas basin (Northern Greece)

Abstract. The potential of repeat-pass space borne SAR (Synthetic Aperture Radar) interferometry has been exploited to investigate spatial patterns of land subsidence in the Anthemountas basin, in the northern part of Greece. The PSI (Persistent Scatterer Interferometry) approach, based on the processing of long series of SAR acquisitions, has been applied to forty-two images acquired in 1995–2001 by ERS1/2 satellites. Interferometric results have been analysed at a basin scale as support for land motion mapping and at a local scale for the characterisation of ground motion events affecting the village of Perea in the Thermaikos municipality and the "Macedonia" international airport. PSI results revealed a moderate subsidence phenomenon along the wider coastal zone of Anthemountas basin corresponding to intense groundwater extraction. Highest values, exceeding −20 mm yr−1, were measured in the airport area where the thickest sequence of compressible Quaternary sediments occurs. Intense subsidence has been detected also in the Perea village (maximum deformation of −10 to −15 mm yr−1), where a series of fractures, causing damages to both buildings and infrastructure, occurred in 2005–2006.

NAVIGATION AND REMOTE SENSING PAYLOADS AND METHODS OF THE SARVANT UNMANNED AERIAL SYSTEM

Abstract. In a large number of scenarios and missions, the technical, operational and economical advantages of UAS-based photogrammetry and remote sensing over traditional airborne and satellite platforms are apparent. Airborne Synthetic Aperture Radar (SAR) or combined optical/SAR operation in remote areas might be a case of a typical "dull, dirty, dangerous" mission suitable for unmanned operation – in harsh environments such as for example rain forest areas in Brazil, topographic mapping of small to medium sparsely inhabited remote areas with UAS-based photogrammetry and remote sensing seems to be a reasonable paradigm. An example of such a system is the SARVANT platform, a fixed-wing aerial vehicle with a six-meter wingspan and a maximumtake- of-weight of 140 kilograms, able to carry a fifty-kilogram payload. SARVANT includes a multi-band (X and P) interferometric SAR payload, as the P-band enables the topographic mapping of densely tree-covered areas, providing terrain profile information. Moreover, the combination of X- and P-band measurements can be used to extract biomass estimations. Finally, long-term plan entails to incorporate surveying capabilities also at optical bands and deliver real-time imagery to a control station. This paper focuses on the remote-sensing concept in SARVANT, composed by the aforementioned SAR sensor and envisioning a double optical camera configuration to cover the visible and the near-infrared spectrum. The flexibility on the optical payload election, ranging from professional, medium-format cameras to mass-market, small-format cameras, is discussed as a driver in the SARVANT development. The paper also focuses on the navigation and orientation payloads, including the sensors (IMU and GNSS), the measurement acquisition system and the proposed navigation and orientation methods. The latter includes the Fast AT procedure, which performs close to traditional Integrated Sensor Orientation (ISO) and better than Direct Sensor Orientation (DiSO), and features the advantage of not requiring the massive image processing load for the generation of tie points, although it does require some Ground Control Points (GCPs). This technique is further supported by the availability of a high quality INS/GNSS trajectory, motivated by single-pass and repeat-pass SAR interferometry requirements.

Residual motion estimation with point targets and its application to airborne repeat-pass SAR interferometry

Due to the lack of accuracy in the navigation system, airborne repeat-pass synthetic aperture radar (SAR) interferometry (InSAR) suffers residual motion errors (RMEs). Previously, we proposed the multisquint technique with point targets (MTPT) to estimate the errors, but its implementation needs improvement and its accuracy has not been assessed by real airborne repeat-pass InSAR data. In this article, the modified MTPT is introduced first. The modification is mainly embodied in two aspects: automatic target selection and noise removing. Because the multisquint (MS) technique is also capable of estimating the RMEs in SAR interferogram and its accuracy has been verified to be high, here it is used as a comparison. In addition, from the viewpoint of MTPT, MS can be understood from another perspective. Using real X-band airborne repeat-pass SAR data, the performance of MTPT is evaluated. The experiment shows that MTPT is able to achieve high accuracy when a large number of point targets are distributed in the observed scene. The limitations of MTPT are also discussed at the end of this article.

Tree height influence on ERS interferometric phase in boreal forest

The European Remote Sensing 1/2 (ERS-1/2) "tandem" coherence has been shown to provide estimates of stem volume in boreal forest in agreement with in situ data. Tree height estimation from ERS interferometric phase represents a further step in the investigations concerning the retrieval of biophysical parameters using repeat-pass synthetic aperture radar (SAR) interferometry. At two test sites located in Sweden and Finland, sets of respectively nine and eight ERS "tandem" interferograms were available. Images acquired under stable winter weather conditions and during nighttime were found to be less affected by atmospheric artifacts. Reduction of atmospheric artifacts in interferograms was performed with a phase screen estimated over a dense grid of open areas. Nonetheless, at each test site, only a limited set of pairs was useful for tree height investigations. Under stable winter conditions, the interferometric tree height obtained from an inversion of the differential interferometric phase at stand level was found to be much lower than the true tree height. Spread and uncertainty in the interferometric tree height measurements were caused by phase noise and residual atmospheric artifacts. Using the semiempirical interferometric water cloud model (IWCM), the modeled interferometric tree height was generally in reasonable agreement with the measurements, showing the need of a phase term in interferometric modeling of forests. The inversion of the IWCM for tree height retrieval showed the strong effect of phase noise and atmospheric artifacts on the estimates. Hence, tree height retrieval from ERS repeat-pass SAR interferometry seems to have limited forestry applications. The results also indicate under what conditions the forest influence is small on digital elevation models derived from repeat-pass interferometry.

InSAR observations of the 1993–95 Bering Glacier (Alaska, U.S.A.) surge and a surge hypothesis

AbstractTime-varying accelerations were observed on Bagley Icefield during the 1993–95 surge of Bering Glacier, Alaska, U.S.A., using repeat-pass synthetic aperture radar interferometry. Observations were from datasets acquired during winter 1991/92 (pre-surge), winter 1993/94 (during the surge) and winter 1995/96 (post-surge). The surge is shown to have extended 110 km up the icefield from Bering Glacier to within 15 km or less of the flow divide. Acceleration and step-like velocity profiles are strongly associated with an along-glacier series of central phase bull’s-eyes with diameters of 0.5–4 km. These bull’s-eyes are interpreted to represent glacier surface rise/fall events of ∼3–30 cm during 1–3 day observation intervals and indicate possible migrating pockets of subglacial water. We present a surge hypothesis that relates late-summer climate to englacial water storage and thence to the subglacial water dynamics — pressurization, hydraulic jacking, depressurization and migration — suggested by our observations.

Post-emplacement lava subsidence and the accuracy of ERS InSAR digital elevation models of volcanoes

Repeat-pass synthetic aperture radar interferometry (InSAR) using data acquired by the ERS platforms is an attractive method for acquiring topographic data of volcanoes. Caution is advised, however, when using this technique in regions covered by young, thick lava flows. In this study, the magnitude of post-emplacement subsidence associated with the 1991-93 lava flow at Mount Etna, Sicily, was measured using differential radar interferometric techniques, and it was found that the rates of subsidence are large enough to contribute a significant component to the measured phase shift, even in ERS data acquired on consecutive orbits. It demonstrates the detrimental effect that such phase shifts have on the accuracy of digital elevation models derived by repeat-pass radar interferometry.

Mapping wet snowcovers with SAR interferometry

The capability of mapping wet snowcovers by means of Synthetic Aperture Radar (SAR) is well known. However, using only the backscattering intensity at C-band, there still remain some incertitudes regarding the effect of either rough snow surfaces or very small liquid water contents for measurements performed at small incidence angles. Since the European remote sensing (ERS) satellites usually operate at 23 incidence angle, it is important to consider these effects. This paper shows that these difficulties can be overcome using repeatpass SAR interferometry. ERS data over a part of the Netherlands and a part of Switzerland are used as illustrative examples. In the situations examined the occurrence of wet snow could not be directly observed by means of the backscattering intensity, but the use of the degree of coherence allowed an easy discrimination. The results are validated with in situ measurements.

ENVISAT ASAR data reduction: impact on SAR interferometry

An onboard synthetic aperture radar (SAR) signal-data reduction algorithm called Flexible Block Adaptive Quantization (FBAQ) was developed in 1994 for the Advanced SAR (ASAR) on the ENVISAT satellite. This paper presents work done in a follow-on study that examined the impact of the data-reduction algorithm on the accuracy of two digital elevation models (DEM) produced by using the technique of repeat-pass SAR interferometry. All three allowable data compression ratios were investigated to determine the maximum compression ratio appropriate for SAR interferometry. Based on the scenes studied it was concluded that a reduction from 8 to 4-bits/sample was the maximum data reduction ratio appropriate for precision SAR interferometry, while 8 to 3-bits/sample and 8 to 2-bits/sample encoding were only appropriate for less-demanding wide-swath applications.

Repeat-pass SAR interferometry over forested terrain

Repeat-pass synthetic aperture radar (SAR) interferometry provides the possibility of producing topographic maps and geocoded as well as radiometrically calibrated radar images. However, the usefulness of such maps and images depends on our understanding of how different types of terrain affect the radar measurements. It is essential that the scene coherence between passes is sdcient. In this paper, we derive a general system model including both radar system and scene scattering properties. The model is used to interpret measurements over a forested area where the scene coherence varies between 0.2 and 0.5. The coherence is found to be sensitive to temperature changes around 0°C but surprisingly insensitive to wind speed. The interferometric height discontinuity at the forest to openfield boundary shows good agreement with in situ tree height measurements for a dense boreal forest, but is observed to decrease for a less dense forest. This suggests the possibility of estimating bole volume from the interferometric tree height and a ground DEM. The decrease of scene coherence over a dense forest with increasing baseline is also used to estimate the effective scattering layer thickness.