Electrical resistivity model of Newberry Volcano in the central Oregon Cascades. Model created by inverting magnetotelluric (MT) data with the ModEM inversion software. All support files used in the inversion are present. Model is centered at 43.7261, -121.3156.

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 (https://github.com/feigl/gipht). 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

In 2002 and 2003 a collaborative effort was undertaken between Lawrence Berkeley National Laboratory, Sandia National Laboratories, the USGS Menlo Park, the USGS Hawaiian Volcano Observatory, and Electromagnetic Instruments Inc. to study the Kilauea volcano in Hawaii using the magnetotelluric (MT) technique. The work was motivated by a desire to improve understanding of the magma reservoirs and conduits within Kilauea and the East and Southwest Rift zones, which has implications for understanding Kilaueas plumbing system. An improved understanding of the rift zones has implications in understanding large-scale landslides that are generated in the Hilina Slump, which produce significant impacts on coastal communities. Up to eight stations operated simultaneously, with multiple remote reference sites, and data were processed using multi-station robust processing techniques. In total, data were acquired at 70 sites over the Southwest and East rift zones. Good to excellent quality data were obtained even in the harshest conditions, such as those encountered on the fresh lava flows of the East Rift Zone (ERZ), where electrical contact resistances are on the order of 100 kOhm. This data supports the continuing efforts to increase geothermal power on the island of Hawaii. Each of the 70 EDI files are the MT impedance tensors for 1 site. There is also a description of the processing of the data and a site map showing the locations of each site.

This submission includes all magnetotelluric (MT) transfer functions acquired during the 2014 EGS stimulation at Newberry Volcano in central Oregon as well as previously acquired MT data for the overall volcano. Plots of all data are provided (including forward response from a model used in a publication now in review in G-cubed). Also included is a kmz file giving all station locations.