Investigating subsidence at volcanoes in northern California using InSAR

Both Medicine Lake Volcano (MLV) and Lassen Volcanic Center (LVC), northern CA, show signs of subsidence at rates of ~1 cm/yr. Leveling and campaign GPS measurements show that MLV has subsided at a constant rate for over 50 years, making the geodetic history of this volcano unique in both its duration and continuity. Here, we summarise and build upon the existing geodetic records at MLV and LVC, using interferometric synthetic aperture radar (InSAR) to extend the time-series of deformation measurements to 2011. We also use the improved spatial resolution of InSAR measurements to investigate causes of long-term subsidence, providing new insight into magmatic storage conditions at MLV and the timescales of deformation due to cooling and crystallization. A large InSAR dataset has been acquired for the volcanoes of northern CA, but application of the data has been limited by extensive noise and incoherence. We analyse multiple datasets from MLV and LVC and, with the use of multi-temporal InSAR analysis methods (noise-based stacking, π-RATE and StaMPS), demonstrate how InSAR may be used more successfully as a monitoring tool in this region. By comparing InSAR results for MLV to past geodetic studies, we demonstrate that subsidence is on going at ~1 cm/yr with no detectable change in rate. We find that the best fitting source geometry to InSAR data is a sill approximated by a horizontal penny-shaped crack, with radius 2 km and depth 11 km, undergoing volume loss at a rate of -0.0022 km³/yr. We discuss possible source mechanisms of long-term subsidence, investigating volume loss due to cooling and crystallization of an intrusion. We calculate the temperature, melt fraction and volume loss of an intrusion over time using petrological information and a numerical thermal model of heat loss by conduction. The geometry of the intrusion is based upon the depth and radius of the penny-shaped crack model. We run simulations for a range of thicknesses between that of a single intrusion (~50 m) and that of the larger column of intrusive material thought to exist beneath the edifice (~7000 m). Using constraints from the geodetic record, we identify a range of sills with volumes < 10 km³ that can account for the deformation recorded at MLV. We use these models to discuss the timing of intrusion and forecast the total duration of cooling. These processes are also significant at LVC and other Cascade volcanoes, where hydrothermal activity is likely to be driven by heat from magmatic intrusions and the exsolution of volatiles that occurs during cooling and crystallization.