Interferometric synthetic aperture radar (InSAR) technology has been widely applied to measure Earth surface motions related to natural and anthropogenic crustal deformation phenomena. With the widespread uptake of data captured by the European Space Agency’s Sentinel-1 mission and other recently launched or planned space-borne SAR missions, the usage of the InSAR technique to detect and monitor Earth surface displacements will increase even more in the coming years. However, InSAR can only measure a one-dimensional motion along the radar line of sight (LOS), which makes interpretation and communication of InSAR measurements challenging, and can add ambiguity to the modelling process. Within this paper, we investigate the implications of the InSAR LOS geometry using simulated and observed deformation phenomena and describe a methodology for multi-geometry data fusion of LOS InSAR measurements from many viewing geometries. We find that projecting LOS measurements to the vertical direction using the incidence angle of the satellite sensor (and implicitly assuming no horizontal motions are present) may result in large errors depending on the magnitude of horizontal motion and on the steepness of the incidence angle. We quantify these errors as the maximum expected error from simulated LOS observations based on a Mogi deformation model. However, we recommend to use LOS observations from several image geometries wherever data are available, in order to solve for vertical and E–W oriented horizontal motion. For an anthropogenic deformation phenomenon observed in seven independent InSAR analyses of Envisat SAR data from the Sydney region, Australia, we find that the strong horizontal motion present could lead to misinterpretation of the actual motion direction when projecting LOS measurements to vertical (uplift instead of subsidence). In this example, the difference between multi-geometry data fusion and vertical projection of LOS measurements (at an incidence angle of 33.8°) reach up to 67% of the maximum vertical displacement rate. Furthermore, the position of maximum vertical motion is displaced horizontally by several hundred metres when the LOS measurements are projected.

Trihedral corner reflectors are being increasingly used as point targets in deformation monitoring studies using interferometric synthetic aperture radar (InSAR) techniques. The frequency and size dependence of the corner reflector Radar Cross Section (RCS) means that no single design can perform equally in all the possible imaging modes and radar frequencies available on the currently orbiting Synthetic Aperture Radar (SAR) satellites. Therefore, either a corner reflector design tailored to a specific data type or a compromise design for multiple data types is required. In this paper, I outline the practical and theoretical considerations that need to be made when designing appropriate radar targets, with a focus on supporting multi-frequency SAR data. These considerations are tested by performing field experiments on targets of different size using SAR images from TerraSAR-X, COSMO-SkyMed and RADARSAT-2. Phase noise behaviour in SAR images can be estimated by measuring the Signal-to-Clutter ratio (SCR) in individual SAR images. The measured SCR of a point target is dependent on its RCS performance and the influence of clutter near to the deployed target. The SCR is used as a metric to estimate the expected InSAR displacement error incurred by the design of each target and to validate these observations against theoretical expectations. I find that triangular trihedral corner reflectors as small as 1 m in dimension can achieve a displacement error magnitude of a tenth of a millimetre or less in medium-resolution X-band data. Much larger corner reflectors (2.5 m or greater) are required to achieve the same displacement error magnitude in medium-resolution C-band data. Compromise designs should aim to satisfy the requirements of the lowest SAR frequency to be used, providing that these targets will not saturate the sensor of the highest frequency to be used. Finally, accurate boresight alignment of the corner reflector can be critical to the overall target performance. Alignment accuracies better than 4 • in azimuth and elevation will incur a minimal impact on the displacement error in X and C-band data.

This article shows how the array of corner reflectors (CRs) in Queensland, Australia, together with highly accurate geodetic synthetic aperture radar (SAR) techniques-also called imaging geodesy-can be used to measure the absolute and relative geometric fidelity of SAR missions. We describe, in detail, the end-to-end methodology and apply it to TerraSAR-X Stripmap (SM) and ScanSAR (SC) data and to Sentinel-1 interferometric wide swath (IW) data. Geometric distortions within images that are caused by commonly used SAR processor approximations are explained, and we show how to correct them during postprocessing. Our results, supported by the analysis of 140 images across the different SAR modes and using the 40 reflectors of the array, confirm our methodology and achieve the limits predicted by theory for both Sentinel-1 and TerraSAR-X. After our corrections, the Sentinel-1 residual errors are 6 cm in range and 26 cm in azimuth, including all error sources. The findings are confirmed by the mutual independent processing carried out at University of Zurich (UZH) and German Aerospace Center (DLR). This represents an improvement of the geolocation accuracy by approximately a factor of four in range and a factor of two in azimuth compared with the standard Sentinel-1 products. The TerraSAR-X results are even better. The achieved geolocation accuracy now approaches that of the global navigation satellite system (GNSS)-based survey of the CRs positions, which highlights the potential of the end-to-end SAR methodology for imaging geodesy.

[1] Knowledge of the number, location, and interseismic slip rate of crustal faults is key to our understanding of the bulk mechanics of continental deformation. However, at present, the sparsity of geodetic velocity observations is a major limiting factor to improving our knowledge of the fault network. In this study, we combine long-swath interferometric synthetic aperture radar (InSAR) data to form an interseismic displacement rate map of over 1000 km length spanning the central Tibetan Plateau. Our dense geodetic observations indicate that strain varies smoothly in space without major concentration at previously mapped faults. Discrepancy between line-of-sight InSAR and GPS observations between 33 ı N-35 ı N could indicate 4 mm/yr of surface uplift in this region. Derived strain rates are less than 5 10 -8 yr -1 across the plateau. At the E-W left-lateral Kunlun fault, postseismic deformation following the 2001 Kokoxili earthquake results in an observed strain rate of 18˙4 10 -8 yr -1 across a 100 km wide zone. Following removal of a viscoelastic model, we estimate an interseismic slip rate of 12˙2 mm/yr for the Kunlun fault. In the central Tibetan conjugate fault zone, the NE-SW striking Amdo-Sewa and NW-SE striking Beng Co faults have right-lateral slip rates of 6˙1 mm/yr and 1-4˙1 mm/yr, respectively. Meanwhile, the NE-SW striking Dongqiao fault has a left-lateral slip rate of 1-2˙1 mm/yr. The low slip rates we estimate for these faults indicate that tectonic strain is more broadly distributed across the plateau interior, perhaps due to the lithosphere being weak and without the abrupt strength contrasts that exist at the plateau margins.Citation: Garthwaite, M. C., H. Wang, and T. J. Wright (2013), Broadscale interseismic deformation and fault slip rates in the central Tibetan Plateau observed using InSAR,