Glacier motion estimation using SAR offset-tracking procedures

Two image-to-image patch offset techniques for estimating feature motion between satellite synthetic aperture radar (SAR) images are discussed. Intensity tracking, based on patch intensity cross-correlation optimization, and coherence tracking, based on patch coherence optimization, are used to estimate the movement of glacier surfaces between two SAR images in both slant-range and azimuth direction. The accuracy and application range of the two methods are examined in the case of the surge of Monacobreen in Northern Svalbard between 1992 and 1996. Offset-tracking procedures of SAR images are an alternative to differential SAR interferometry for the estimation of glacier motion when differential SAR interferometry is limited by loss of coherence, i.e., in the case of rapid and incoherent flow and of large acquisition time intervals between the two SAR images. In addition, an offset-tracking procedure in the azimuth direction may be combined with differential SAR interferometry in the slant-range direction in order to retrieve a two-dimensional displacement map when SAR data of only one orbit configuration are available

Application of Differential and Polarimetric Synthetic Aperture Radar (SAR) Interferometry for Studying Natural Hazards

In the following work, I address the problem of coherence loss in standard Differential Interferometric SAR (DInSAR) processing, which can result in incomplete or poor quality deformation measurements in some areas. I incorporate polarimetric information with DInSAR in a technique called Polarimetric SAR Interferometry (PolInSAR) in order to acquire more accurate and detailed maps of surface deformation. In Chapter 2, I present a standard DInSAR study of the Ahar double earthquakes (Mw=6.4 and 6.2) which occurred in northwest Iran, August 11, 2012. The DInSAR coseismic deformation map was affected by decorrelation noise. Despite this, I employed an advanced inversion technique, in combination with a Coulomb stress analysis, to find the geometry and the slip distribution on the ruptured fault plane. The analysis shows that the two earthquakes most likely occurred on a single fault, not on conjugate fault planes. This further implies that the minor strike-slip faults play more significant role in accommodating convergence stress accumulation in the northwest part of Iran. Chapter 3 presents results from the application of PolInSAR coherence optimization on quad-pol RADARSAT-2 images. The optimized solution results in the identification of a larger number of reliable measurement points, which otherwise are not recognized by the standard DInSAR technique. I further assess the quality of the optimized interferometric phase, which demonstrates an increased phase quality with respect to those phases recovered by applying standard DInSAR alone. Chapter 4 discusses results from the application of PolInSAR coherence optimization from different geometries to the study of creep on the Hayward fault and landslide motions near Berkeley, CA. The results show that the deformation rates resolved by PolInSAR are in agreement with those of standard DInSAR. I also infer that there is potential motion on a secondary fault, northeast and parallel to the Hayward fault, which may be creeping with a lower velocity

Parameters affecting interferometric coherence and implications for long-term operational monitoring of mining-induced surface deformation

Includes abstract.Includes bibliographical references.Surface deformation due to underground mining poses risks to health and safety as well as infrastructure and the environment. Consequently, the need for long-term operational monitoring systems exists. Traditional field-based measurements are point-based meaning that the full extent of deforming areas is poorly understood. Field-based techniques are also labour intensive if large areas are to be monitored on a regular basis. To overcome these limitations, this investigation considered traditional and advanced differential radar interferometry techniques for their ability to monitor large areas over time, remotely. An area known to be experiencing mining induced surface deformation was used as test case. The agricultural nature of the area implied that signal decorrelation effects were expected. Consequently, four sources of data, captured at three wavelengths by earth-orbiting satellites were obtained. This provided the opportunity to investigate different phase decorrelation effects on data from standard imaging platforms using real-world deformation phenomenon as test-case. The data were processed using standard dInSAR and polInSAR techniques. The deformation measurement results together with an analysis of parameters most detrimental to long-term monitoring were presented. The results revealed that, contrary to the hypothesis, polInSAR techniques did not provide an enhanced ability to monitor surface deformation compared to dInSAR techniques. Although significant improvements in coherence values were obtained, the spatial heterogeneity of phase measurements could not be improved. Consequently, polInSAR could not overcome ecorrelation associated with vegetation cover and evolving land surfaces. However, polarimetric information could be used to assess the scattering behaviour of the surface, thereby guiding the definition of optimal sensor configuration for long-term monitoring. Despite temporal and geometric decorrelation, the results presented demonstrated that mining-induced deformation could be measured and monitored using dInSAR techniques. Large areas could be monitored remotely and the areal extent of deforming areas could be assessed, effectively overcoming the limitations of field-based techniques. Consequently, guidelines for the optimal sensor configuration and image acquisition strategy for long-term operational monitoring of mining-induced surface deformation were provided

Application of DInSAR-GPS optimization for derivation of fine-scale surface motion maps of Southern California

A method based on random field theory and Gibbs-Markov random fields equivalency within Bayesian statistical framework is used to derive 3-D surface motion maps from sparse global positioning system (GPS) measurements and differential interferometric synthetic aperture radar (DInSAR) interferogram in the southern California region. The minimization of the Gibbs energy function is performed analytically, which is possible in the case when neighboring pixels are considered independent. The problem is well posed and the solution is unique and stable and not biased by the continuity condition. The technique produces a 3-D field containing estimates of surface motion on the spatial scale of the DInSAR image, over a given time period, complete with error estimates. Significant improvement in the accuracy of the vertical component and moderate improvement in the accuracy of the horizontal components of velocity are achieved in comparison with the GPS data alone. The method can be expanded to account for other available data sets, such as additional interferograms, lidar, or leveling data, in order to achieve even higher accuracy

Earthquake source parameters of the 2009 Mw 7.8 Fiordland (New Zealand) earthquake from L-band InSAR observations

The 2009 MW7.8 Fiordland (New Zealand) earthquake is the largest to have occurred in New Zealand since the 1931 Mw 7.8 Hawke’s Bay earthquake, 1 000 km to the northwest. In this paper two tracks of ALOS PALSAR interferograms (one ascending and one descending) are used to determine fault geometry and slip distribution of this large earthquake. Modeling the event as dislocation in an elastic half-space suggests that the earthquake resulted from slip on a SSW-NNE orientated thrust fault that is associated with the subduction between the Pacific and Australian Plates, with oblique displacement of up to 6.3 m. This finding is consistent with the preliminary studies undertaken by the USGS using seismic data

The 2010 MW 6.8 Yushu (Qinghai, China) earthquake: constraints provided by InSAR and body wave seismology

By combining observations from satellite radar, body wave seismology and optical imagery, we have determined the fault segmentation and sequence of ruptures for the 2010 Mw 6.8 Yushu (China) earthquake. We have mapped the fault trace using displacements from SAR image matching, interferometric phase and coherence, and 2.5 m SPOT-5 satellite images. Modeling the event as an elastic dislocation with three segments fitted to the fault trace suggests that the southeast and northwest segments are near vertical, with the central segment dipping 70° to the southwest; slip occurs mainly in the upper 10 km, with a maximum slip of 1.5 m at a depth of 4 km on the southeastern segment. The maximum slip in the top 1 km (i.e., near surface) is up to 1.2 m, and inferred locations of significant surface rupture are consistent with displacements from SAR image matching and field observations. The radar interferograms show rupture over a distance of almost 80 km, much larger than initial seismological and field estimates of the length of the fault. Part of this difference can be attributed to slip on the northwestern segment of the fault being due to an Mw 6.1 aftershock two hours after the main event. The remaining difference can be explained by a non-uniform slip distribution with much of the moment release occurring at depths of less than 10 km. The rupture on the central and southeastern segments of the fault in the main shock propagated at a speed of 2.5 km/s southeastward toward the town of Yushu located at the end of this segment, accounting for the considerable building damage. Strain accumulation since the last earthquake on the fault segment beyond Yushu is equivalent to an Mw 6.5 earthquake

SAR interferometry at Venus for topography and change detection

AbstractSince the Magellan radar mapping of Venus in the early 1990’s, techniques of synthetic aperture radar interferometry (InSAR) have become the standard approach to mapping topography and topographic change on Earth. Here we investigate a hypothetical radar mission to Venus that exploits these new methods. We focus on a single spacecraft repeat-pass InSAR mission and investigate the radar and mission parameters that would provide both high spatial resolution topography as well as the ability to detect subtle variations in the surface. Our preferred scenario is a longer-wavelength radar (S or L-band) placed in a near-circular orbit at 600km altitude. Using longer wavelengths minimizes the required radar bandwidth and thus the amount of data that will be transmitted back to earth; it relaxes orbital control and knowledge requirements. During the first mapping cycle a global topography map would be assembled from interferograms taken from adjacent orbits. This approach is viable due to the slow rotation rate of Venus, causing the interferometric baseline between adjacent orbits to vary from only 11km at the equator to zero at the inclination latitude. To overcome baseline decorrelation at lower latitudes, the center frequency of a repeated pass will be adjusted relative to the center frequency of its reference pass. During subsequent mapping cycles, small baseline SAR acquisitions will be used to search for surface decorrelation due to lava flows. While InSAR methods are used routinely on Earth, their application to Venus could be complicated by phase distortions caused by the thick Venus atmosphere

Improvement of the Accuracy of InSAR Image Co-Registration Based On Tie Points – A Review

Interferometric Synthetic Aperture Radar (InSAR) is a new measurement technology, making use of the phase information contained in the Synthetic Aperture Radar (SAR) images. InSAR has been recognized as a potential tool for the generation of digital elevation models (DEMs) and the measurement of ground surface deformations. However, many critical factors affect the quality of InSAR data and limit its applications. One of the factors is InSAR data processing, which consists of image co-registration, interferogram generation, phase unwrapping and geocoding. The co-registration of InSAR images is the first step and dramatically influences the accuracy of InSAR products. In this paper, the principle and processing procedures of InSAR techniques are reviewed. One of important factors, tie points, to be considered in the improvement of the accuracy of InSAR image co-registration are emphatically reviewed, such as interval of tie points, extraction of feature points, window size for tie point matching and the measurement for the quality of an interferogram