On 12 January 2020, Taal volcano, Philippines, erupted after 43 years of repose, affecting more than 500,000 people. Using interferometric synthetic aperture radar (InSAR) data, we present the complete pre- to post-eruption analyses of the deformation of Taal. We find that: 1) prior to eruption, the volcano experienced long-term deflation followed by short-term inflation, reflecting the depressurization-pressurization of its ~5 km depth magma reservoir; 2) during the eruption, the magma reservoir lost a volume of 0.531 +/- 0.004 km^3 while a 0.643 +/- 0.001 km^3 lateral dike was emplaced; and 3) post-eruption analyses reveal that the magma reservoir started recovery approximately 3 weeks after the main eruptive phase. We propose a conceptual analysis explaining the eruption and address why, despite the large volume of magma emplaced, the dike remained at depth. We also report the unique and significant contribution of InSAR data during the peak of the crisis.

We present evidence of volcano-tectonic interactions at Sabancaya volcano that we relate to episodic magma injection and high regional fluid pore pressures. We present a surface deformation time series at Sabancaya including observations from ERS-1/2, Envisat, Sentinel-1, COSMO-SkyMed, and TerraSAR-X that spans June 1992 - February 2019. These data show deep seated inflation northwest of Sabancaya from 1992-1997 and 2013-2019, as well as creep and rupture on multiple faults. Afterslip on the Mojopampa fault following a 2013 Mw 5.9 earthquake is anomalously long-lived, continuing for at least six years. The best fit fault plane for the afterslip is right-lateral motion on an EW striking fault at 1 km depth. We also model surface deformation from two 2017 earthquakes (Mw 4.4 and Mw 5.2) on unnamed faults, for which the best fit models are NW striking normal faults at 1-2 km depth. Our best fit model for a magmatic inflation source (13 km depth, volume change of 0.04 to 0.05 km^3 yr^-1), induces positive Coulomb static stress changes on these modeled fault planes. Comparing these deformation results with evidence from satellite thermal and degassing data, field observations, and seismic records, we interpret strong pre-eruptive seismicity at Sabancaya as a consequence of magmatic intrusions destabilizing tectonic faults critically stressed by regionally high fluid pressures. High fluid pressure likely also promotes fault creep driven by static stress transfer from the inflation source. We speculate that strong seismicity near volcanoes will be most likely with high pore fluid pressures and significant, offset magmatic inflation.

A combination of magmatic and tectonic processes occur on the western branch of the East African Rift System (EARS) driven by active volcanoes adjacent to active rift faults. Mt. Nyiragongo and Mt. Nyamuragira have the most recent eruptive histories of the 8 volcanoes in the Virunga Volcanic Zone (VVZ) located in a region between Rwanda, Uganda and Democratic Republic of Congo (DRC). On May 22nd 2021, Mt.Nyiragongo erupted the first major eruption following its 2002 eruption. This eruption didn’t have the common precursory seismic activity expected before an eruption as was observed in the 2002 eruption seismic record. Rather, there were numerous post-eruption earthquake events with the largest of those events being a magnitude ML 5.1. Around the region of the earthquake swarm, there was observable ground deformation in the city of Goma and Rubavu where surface fissures destroyed houses and split roads apart. This deformation appears to be related to a N-S striking dike intrusion from the volcano trending south towards and under Lake Kivu, according to observed seismicity. In collaboration with the Government of Rwanda, following the May Mt. Nyiragongo eruption, we established a network of 6 seismometers (2 Meridian Compact PH and 4 Trillium Compact PH) operating at 100 sps and two complimentary raspberry Shake and Booms (SBS) around Lake Kivu. This study will focus on characterizing deformation associated with the eruption and the subsequent seismic swarm. Here we present model results based on deformation during the May 2021 eruption as recorded through ALOS InSAR scenes to understand slip concentration during the dike intrusion. Using GTDef, a set of algorithms developed in MATLAB that can incorporate a wide range of geodetic data types to model deformation observed on the Earth’s surface, we model the slip distribution in this region based on the current hypothesis that the observed seismicity was a result of a dike intrusion defined by the southward propagation of the seismic swarm from Mt. Nyiragongo. Given an approximate source, we determine a preferred GNSS/GPS network design based on resolution-cost of additional stations at given locations and discuss first order characterization of the observed deformation.