Imaging irregular magma reservoirs with InSAR and GPS observations: Application to Kilauea and Copahue volcanoes

The availability of synthetic aperture radar (SAR) interferometry (InSAR) data has increased our awareness of the complexity of volcano deformation sources. InSAR's spatial completeness helps identify or clarify source process mechanisms at volcanoes (i.e. Mt. Etna east flank motion; Lazufre crustal magma body; Kilauea dike complexity) and also improves potential model realism. In recent years, Bayesian inference methods have gained widespread use because of their ability to constrain not only source model parameters, but also their uncertainties. They are computationally intensive, however, which tends to limit them to a few geometrically rather simple source representations (for example, spheres). An alternative approach involves solving for irregular pressure and/or density sources from a three-dimensional (3-D) grid of source/density cells. This method has the ability to solve for arbitrarily shaped bodies of constant absolute pressure/density difference. We compare results for both Bayesian (a Markov chain Monte Carlo algorithm) and the irregular source methods for two volcanoes: Kilauea, Hawaii, and Copahue, Argentina-Chile border. Kilauea has extensive InSAR and GPS databases from which to explore the results for the irregular method with respect to the Bayesian approach, prior models, and an extensive set of ancillary data. One caveat, however, is the current restriction in the irregular model inversion to volume-pressure sources (and at a single excess pressure change), which limits its application in cases where sources such as faults or dikes are present. Preliminary results for Kilauea summit deflation during the March 2011 Kamoamoa eruption suggests a northeast-elongated magma body lying roughly 1-1.5 km below the surface. Copahue is a southern Andes volcano that has been inflating since early 2012, with intermittent summit eruptive activity since late 2012. We have an extensive InSAR time series from RADARSAT-2 and COSMO-SkyMed data, although both are from descending tracks. Preliminary modeling suggests a very irregular magma body that extends from the volcanic edifice to less than 5 km depth and located slightly north of the summit at shallow depths but to the ENE at greater depths. In our preliminary analysis, we find that there are potential limitations and trade-offs in the Bayesian results that suggest the simplicity of the assumed analytic source may generate systematic biases in source parameters. Instead, the irregular 3-D solution appears to provide greater realism, but is limited in the number and type of sources that can be modeled.