Events on Tuesday, December 9th, 2008
- Chaos & Complex Systems Seminar
- A technique for modeling radar interferograms without phase unwrapping: Application to the M 5 Fawnskin, California earthquake of 4 December 1992
- Time: 12:05 pm
- Place: 4274 Chamberlin (Refreshments will be served)
- Speaker: Kurt Feigl, UW Department of Geology and Geophysics
- Abstract: Interferometric analysis of synthetic aperture radar images (InSAR) measures the phase shifts between two images acquired at two distinct times. These ambiguous "wrapped" phase values range from 1/2 to +1/2 cycles. The standard approach interprets the phase values in terms of the change in distance between the ground and the radar instrument by resolving the integer ambiguities in a process known as "unwrapping". To avoid unwrapping, we have developed, validated, and applied a new method for modeling the wrapped phase data directly. The method defines a cost function in terms of wrapped phase to measure the misfit between the observed and modeled values of phase. By minimizing the cost function with a simulated annealing algorithm, the method estimates parameters in a nonlinear model. Since the wrapped phase residuals are compatible with a von Mises distribution, several parametric statistical tests can be used to evaluate the fit of the model to the data. The method, named General Inversion for Phase Technique (GIPhT), can handle noisy, wrapped phase data. Applying GIPhT to two interferograms in the area of Fawnskin, California, we estimate a set of model parameters describing a magnitude 5 aftershock of the 1992 Landers earthquake. The resulting simulation fits the data well. The phase final residuals have a circular mean deviation less than 0.15 cycles/datum. Sampling the final residuals, we find the circular standard deviation of a phase measurement to be approximately 0.2 cycle, corresponding to 6 mm in range.