Graph Theory for Analyzing Pair-wise Data: Application to Interferometric Synthetic Aperture Radar Data ReinischEC_MS_thesis_20160728.pdf

Graph theory is useful for estimating time-dependent model parameters via weighted least-squares using interferometric synthetic aperture radar (InSAR) data. Plotting acquisition dates (epochs) as vertices and pair-wise interferometric combinations as edges defines an incidence graph. The edge-vertex incidence matrix and the normalized edge Laplacian matrix are factors in the covariance matrix for the pair-wise data. Using empirical measures of residual scatter in the pair-wise observations, we estimate the variance at each epoch by inverting the covariance of the pair-wise data. We evaluate the rank deficiency of the corresponding least-squares problem via the edge-vertex incidence matrix. We implement our method in a MATLAB software package called GraphTreeTA available on GitHub ( We apply temporal adjustment to the data set described in Lu et al. (2005) at Okmok volcano, Alaska, which erupted most recently in 1997 and 2008. The data set contains 44 differential volumetric changes and uncertainties estimated from interferograms between 1997 and 2004. Estimates show that approximately half of the magma volume lost during the 1997 eruption was recovered by the summer of 2003. Between June 2002 and September 2003, the estimated rate of volumetric increase is (6.2 +/- 0.6) x 10^6 m^3/yr. Our preferred model provides a reasonable fit that is compatible with viscoelastic relaxation in the five years following the 1997 eruption. Although we demonstrate the approach using volumetric rates of change, our formulation in terms of incidence graphs applies to any quantity derived from pair-wise differences, such as wrapped phase or wrapped residuals.

Date of final oral examination: 05/19/2016 This thesis is approved by the following members of the Final Oral Committee: Kurt L. Feigl, Professor, Geoscience Michael Cardiff, Assistant Professor, Geoscience Clifford H. Thurber, Vilas Distinguished Professor, Geoscience Thesis submitted by Elena C. Reinisch at the University of Wisconsin - Madison in partial fulfillment of the requirements for the degree of Master of Science (Geophysics). Methods developed and documented in this thesis are used by the PoroTomo project for analyzing InSAR data at Brady Hot Springs, Nevada.

Data and Resources

Metadata Source

Additional Info

Field Value
Citation Date 2016-07-28T00:00:00-06:00

Harvest Information

Original ID f0000000-58cc-4372-a567-000000001075
Index Date 2020-01-27T15:40:32-07:00
Original Format ISO-USGIN
Original Version 1.2


Name Kurt Feigl
Position primary contact
Organization University of Wisconsin

Geographic Extent

North Bound 39.9883
South Bound 39.5883
East Bound -118.8167
West Bound -119.2167