Constraints on time-dependent volcanic source models at Long Valley Caldera from 1996 to 2009 using InSAR and geodetic measurements

SUMMARY Continuous monitoring of Long Valley Caldera since the late 1970s, including data from seismic and geodetic networks has shown renewed episodic unrest with accelerated uplift separated by reduced uplift, no activity, or slow deflation. We examine the time-dependent behaviour at Long Valley Caldera from 1996 to 2009 by integrating InSAR and continuous GPS (CGPS) measurements. The ERS-1/2 radar data between 1996 and 2008 and reprocessed three-component CGPS data from the Long Valley GPS network from 1996 to 2009 were combined to invert for source geometry and volume change for the following deformation episodes: 1997–1998 uplift, 2002–2003 uplift, 2004–2007 slow subsidence, and 2007–2009 slow uplift. We employed non-linear Monte-Carlo random search approaches (random cost and simulated annealing) in our inversion and examined models including spherical and finite sources (dipping prolate spheroid). Our results show that the sources of all post-2000 events (i.e. 2002–2003, 2007–2009 uplift and 2004–2007 subsidence) locate at shallow depths of ∼6–8 km and have nearly identical surface locations, suggesting that these events may be caused by the same source in the mid-crust, possibly a mixture of hydrothermal and partial-melt magma. All three events, 2002–2003, 2004–2007, 2007–2009, are characterized by the low total volume change: ∼0.01, ∼–0.003, ∼0.006 km3, respectively, with corresponding volume change rates of ∼0.007, ∼–0.001, ∼0.002 km3 yr–1. The 1997–1998 inflationary event has a steeper source geometry and much greater volume change rate (∼0.03 km3 yr–1) than the other events, in agreement with previous studies. This suggests it is possibly driven by magma intrusion beneath the resurgent dome from a deeper source. If we regard post-2000 events as proxies for future eruption hazard, the inferred source dynamics (e.g. mid-crustal location and subdued volume change) and accompanied decrease in earthquake activity from these post-2000 events suggest that the probability for near-term eruption is low. Our study demonstrates that CGPS, along with InSAR, provide a valuable tool for monitoring time-dependent source processes at active volcanic regions.

[1]  Thomas R. Walter,et al.  Time‐dependent volcano source monitoring using interferometric synthetic aperture radar time series: A combined genetic algorithm and Kalman filter approach , 2010 .

[2]  S. Owen,et al.  Estimation of interplate coupling in the Nankai trough, Japan using GPS data from 1996 to 2006 , 2010 .

[3]  M. Bonafede,et al.  A numerical model of dyke propagation in layered elastic media , 2010 .

[4]  Thomas R. Walter,et al.  Randomly iterated search and statistical competency as powerful inversion tools for deformation source modeling: Application to volcano interferometric synthetic aperture radar data , 2009 .

[5]  David P. Hill,et al.  Analytical modeling of gravity changes and crustal deformation at volcanoes: The Long Valley caldera, California, case study , 2009 .

[6]  A. Newman,et al.  Constraints on continued episodic inflation at Long Valley Caldera, based on seismic and geodetic observations , 2009 .

[7]  Henriette Sudhaus,et al.  Improved source modelling through combined use of InSAR and GPS under consideration of correlated data errors: application to the June 2000 Kleifarvatn earthquake, Iceland , 2009 .

[8]  Jeffrey T. Freymueller,et al.  Tracking magma volume recovery at Okmok volcano using GPS and an unscented Kalman filter , 2009 .

[9]  Paolo Berardino,et al.  Uplift and magma intrusion at Long Valley caldera from InSAR and gravity measurements , 2009 .

[10]  Jean Virieux,et al.  Seismic reflections reveal a massive melt layer feeding Campi Flegrei caldera , 2008 .

[11]  Antonio Pepe,et al.  The 2004–2006 uplift episode at Campi Flegrei caldera (Italy): Constraints from SBAS‐DInSAR ENVISAT data and Bayesian source inference , 2008 .

[12]  T. Menand The mechanics and dynamics of sills in layered elastic rocks and their implications for the growth of laccoliths and other igneous complexes , 2008 .

[13]  Paolo Berardino,et al.  Surface deformation of Long Valley caldera and Mono Basin, California, investigated with the SBAS-InSAR approach , 2007 .

[14]  Daniel Dzurisin,et al.  Volcano deformation : geodetic monitoring techniques , 2007 .

[15]  Yehuda Bock,et al.  Spatiotemporal filtering using principal component analysis and Karhunen-Loeve expansion approaches for regional GPS network analysis , 2006 .

[16]  T. Dixon,et al.  A four-dimensional viscoelastic deformation model for Long Valley Caldera, California, between 1995 and 2000 , 2006 .

[17]  James Jackson,et al.  The 1994 Sefidabeh (eastern Iran) earthquakes revisited: new evidence from satellite radar interferometry and carbonate dating about the growth of an active fold above a blind thrust fault , 2006 .

[18]  D. Hill Unrest in Long Valley Caldera, California, 1978–2004 , 2006, Geological Society, London, Special Publications.

[19]  Rowena B. Lohman,et al.  Some thoughts on the use of InSAR data to constrain models of surface deformation: Noise structure and data downsampling , 2005 .

[20]  Zhi-Xun Shen,et al.  Contemporary Crustal Deformation Around Southeast Borderland of Tibetan Plateau , 2004 .

[21]  Yehuda Bock,et al.  Error analysis of continuous GPS position time series , 2004 .

[22]  J. Langbein Deformation of the Long Valley Caldera, California: inferences from measurements from 1988 to 2001 , 2003 .

[23]  C. W. Roberts,et al.  The mechanics of unrest at Long Valley caldera, California. 2. Constraining the nature of the source using geodetic and micro-gravity data , 2003 .

[24]  D. Hill,et al.  Three-dimensional crustal structure of Long Valley caldera, California, and evidence for the migration of CO2 under Mammoth Mountain , 2003 .

[25]  Gianfranco Fornaro,et al.  A new algorithm for surface deformation monitoring based on small baseline differential SAR interferograms , 2002, IEEE Trans. Geosci. Remote. Sens..

[26]  H. Zebker,et al.  Fault Slip Distribution of the 1999 Mw 7.1 Hector Mine, California, Earthquake, Estimated from Satellite Radar and GPS Measurements , 2002 .

[27]  G. Fornaro,et al.  Modeling surface deformation observed with synthetic aperture radar interferometry at Campi Flegrei caldera , 2001 .

[28]  Mark Simons,et al.  Finite source modelling of magmatic unrest in Socorro, New Mexico, and Long Valley, California , 2001 .

[29]  Paul Segall,et al.  Estimating source parameters from deformation data, with an application to the March 1997 earthquake swarm off the Izu Peninsula, Japan , 2001 .

[30]  R. Hanssen Radar Interferometry: Data Interpretation and Error Analysis , 2001 .

[31]  C. W. Chen,et al.  Two-dimensional phase unwrapping with use of statistical models for cost functions in nonlinear optimization. , 2001, Journal of the Optical Society of America. A, Optics, image science, and vision.

[32]  Kristy F. Tiampo,et al.  Spherical and ellipsoidal volcanic sources at Long Valley caldera, California, using a genetic algorithm inversion technique , 2000 .

[33]  G. Wadge,et al.  An accurate and efficient method for including the effects of topography in three‐dimensional elastic models of ground deformation with applications to radar interferometry , 2000 .

[34]  Fuk K. Li,et al.  Synthetic aperture radar interferometry , 2000, Proceedings of the IEEE.

[35]  Roberts,et al.  Magma intrusion beneath long valley caldera confirmed by temporal changes in gravity , 1999, Science.

[36]  Charles A. Williams,et al.  The effects of topography on magma chamber deformation models: Application to Mt. Etna and radar interferometry , 1998 .

[37]  Robert W. King,et al.  Estimating regional deformation from a combination of space and terrestrial geodetic data , 1998 .

[38]  P. Dawson,et al.  Nonlinear teleseismic tomography at Long Valley Caldera, using three‐dimensional minimum travel time ray tracing , 1995 .

[39]  Bernd A. Berg,et al.  Locating global minima in optimization problems by a random-cost approach , 1993, Nature.

[40]  J. C. Savage Principal Component Analysis of Geodetically Measured Deformation in Long Valley Caldera, Eastern California, 1983–1987 , 1988 .

[41]  James H. Dieterich,et al.  Deformation from Inflation of a Dipping Finite Prolate Spheroid in an Elastic Half‐Space as a Model for Volcanic Stressing , 1988 .

[42]  Harold E. Lockwood,et al.  Remote sensing and digital image processing , 1982 .

[43]  K. Mogi Relations between the Eruptions of Various Volcanoes and the Deformations of the Ground Surfaces around them , 1958 .