Prevalence of Seismic Rate Anomalies Preceding Volcanic Eruptions in Alaska

Seismic rate increases often precede eruptions at volcanoes worldwide. However, many eruptions occur without such precursors. Additionally, identifying seismic rate increases near volcanoes with high levels of background seismicity is non-trivial and many periods of elevated seismicity occur without ensuing eruptions, limiting their usefulness for forecasting in some cases. Although these issues are commonly known, efforts to quantify them are limited. In this study, we consistently apply a common statistical tool, the β-statistic, to seismically monitored eruptions in Alaska of various styles to determine the overall prevalence of seismic rate anomalies immediately preceding eruptions. We find that 6 out of 20 (30%) eruptions have statistically significant precursory seismic rate increases. Of these 6 eruptions, 3 of them occur at volcanoes with relatively felsic compositions, repose periods >15 years, and VEI ≥ 3. Overall, our results confirm that seismic rate increases are common prior to larger eruptions at long dormant, “closed-system” volcanoes, but uncommon preceding smaller eruptions at more frequently active, “open-system” volcanoes with more mafic magmas. We also explore the rate of other anomalies not precursory to eruptions and investigate their origins. Some of these non-eruptive anomalies can be explained by aftershocks of regional seismic events, magmatic activity that did not lead to eruption, or unrest at other nearby volcanoes. Some open-system volcanoes have high non-eruptive anomaly rates and low pre-eruptive anomaly rates and are thus not amenable to forecasting based on earthquake catalogs. In this study, we find that 31% of anomalies lead to eruption. With continued calibration at more volcanoes, the β-statistic that we apply may be used more broadly to analyze future periods of seismic unrest at other volcanoes, properly placing such episodes into the context of the long-term background rate. These results may be useful for informing future eruption forecasts around the world, and the statistical tool may aid volcano observatories in identifying future seismic rate anomalies under changing network conditions

Volcanic earthquake catalog enhancement using integrated detection, matched-filtering, and relocation tools

Volcanic earthquake catalogs are an essential data product used to interpret subsurface volcanic activity and forecast eruptions. Advances in detection techniques (e.g., matched-filtering, machine learning) and relative relocation tools have improved catalog completeness and refined event locations. However, most volcano observatories have yet to incorporate these techniques into their catalog-building workflows. This is due in part to complexities in operationalizing, automating, and calibrating these techniques in a satisfactory way for disparate volcano networks and their varied seismicity. In an effort to streamline the integration of catalog-enhancing tools at the Alaska Volcano Observatory (AVO), we have integrated four popular open-source tools: REDPy, EQcorrscan, HypoDD, and GrowClust. The combination of these tools offers the capability of adding seismic event detections and relocating events in a single workflow. The workflow relies on a combination of standard triggering and cross-correlation clustering (REDPy) to consolidate representative templates used in matched-filtering (EQcorrscan). The templates and their detections are then relocated using the differential time methods provided by HypoDD and/or GrowClust. Our workflow also provides codes to incorporate campaign data at appropriate junctures, and calculate magnitude and frequency index for valid events. We apply this workflow to three datasets: the 2012–2013 seismic swarm sequence at Mammoth Mountain (California), the 2009 eruption of Redoubt Volcano (Alaska), and the 2006 eruption of Augustine Volcano (Alaska); and compare our results with previous studies at each volcano. In general, our workflow provides a significant increase in the number of events and improved locations, and we relate the event clusters and temporal progressions to relevant volcanic activity. We also discuss workflow implementation best practices, particularly in applying these tools to sparse volcano seismic networks. We envision that our workflow and the datasets presented here will be useful for detailed volcano analyses in monitoring and research efforts

The 2013–2020 seismic activity at Sabancaya Volcano (Peru): Long lasting unrest and eruption

Sabancaya volcano is the youngest and second most active volcano in Peru. It is part of the Ampato-Sabancaya volcanic complex which sits to the south of the ancient Hualca Hualca volcano and several frequently active faults, thus resulting in complex volcano-tectonic interactions. After 15 years of repose, in 2013, a series of 4 earthquakes with magnitude >4.5 occurred within 24 h, marking the beginning of a new episode of unrest. Several additional swarms of earthquakes occurred in the following years until magmatic eruptive activity started on 6 November 2016. This activity is ongoing as of this writing, with an average of 50 explosions per day. In this study, we present results of multiparametric monitoring of Sabancaya's activity observed during 2013–2020. Seismic data are used to create a one-dimensional seismic velocity model, to catalog, locate, and characterize earthquakes, to detect repeating earthquake families, and to monitor seismic velocity variations by ambient noise cross-correlation. These analyses are complemented by visual and remote sensing observations and ground deformation measurements. All monitored parameters showed significant changes on 6 November 2016, the day of eruption onset, thus dividing the eruptive activity into pre-eruptive and eruptive stages. The unrest is characterized by high levels of seismic activity with hundreds of events detected per day. Volcano-tectonic (VT) earthquakes were dominant during the pre-eruptive period while long-period (LP) events and explosions have been most numerous since the eruption onset. Earthquake locations highlight long-lasting seismogenic zones along multiple previously active regional faults, as well as along newly identified faults. This VT seismicity is mainly distributed in a sector from the northwest to the east of the volcanic complex at distances of up to 30 km from the crater. We focus our analysis on two eruptive episodes: the eruption onset and subsequent crater migration from south to north, and the increase of lava dome extrusion rate in 2019. Both episodes are accompanied by seismic velocity decreases of up to 0.2% and are preceded by a few weeks by bursts of distal VT activity, including numerous repeating earthquakes. These repeated events were located on several remote tectonic faults (5–25 km from the vent). We suggest that these phenomena could be due to the injection of a batch of magma in the deep reservoir and/or conduit, which would generate 1) a pressure wave propagating in the hydrothermal system, triggering the bursts of seismic activity and 2) slow rising of magma by melting old material filling the conduit that eventually produced the eruptive and dome growth acceleration events.SCOPUS: ar.jinfo:eu-repo/semantics/publishe

The 2013–2020 seismic activity at Sabancaya Volcano (Peru): Long lasting unrest and eruption

International audienceSabancaya volcano is the youngest and second most active volcano in Peru. It is part of the Ampato-Sabancaya volcanic complex which sits to the south of the ancient Hualca Hualca volcano and several frequently active faults, thus resulting in complex volcano-tectonic interactions. After 15 years of repose, in 2013, a series of 4 earthquakes with magnitude >4.5 occurred within 24 h, marking the beginning of a new episode of unrest. Several additional swarms of earthquakes occurred in the following years until magmatic eruptive activity started on 6 November 2016. This activity is ongoing as of this writing, with an average of 50 explosions per day. In this study, we present results of multiparametric monitoring of Sabancaya's activity observed during 2013-2020. Seismic data are used to create a one-dimensional seismic velocity model, to catalog, locate, and characterize earthquakes, to detect repeating earthquake families, and to monitor seismic velocity variations by ambient noise cross-correlation. These analyses are complemented by visual and remote sensing observations and ground deformation measurements. All monitored parameters showed significant changes on 6 November 2016, the day of eruption onset, thus dividing the eruptive activity into pre-eruptive and eruptive stages. The unrest is characterized by high levels of seismic activity with hundreds of events detected per day. Volcano-tectonic (VT) earthquakes were dominant during the pre-eruptive period while long-period (LP) events and explosions have been most numerous since the eruption onset. Earthquake locations highlight long-lasting seismogenic zones along multiple previously active regional faults, as well as along newly identified faults. This VT seismicity is mainly distributed in a sector from the northwest to the east of the volcanic complex at distances of up to 30 km from the crater. We focus our analysis on two eruptive episodes: the eruption onset and subsequent crater migration from south to north, and the increase of lava dome extrusion rate in 2019. Both episodes are accompanied by seismic velocity decreases of up to 0.2% and are preceded by a few weeks by bursts of distal VT activity, including numerous repeating earthquakes. These repeated events were located on several remote tectonic faults (5-25 km from the vent). We suggest that these phenomena could be due to the injection of a batch of magma in the deep reservoir and/or conduit, which would generate 1) a pressure wave propagating in the hydrothermal system, triggering the bursts of seismic activity and 2) slow rising of magma by melting old material filling the conduit that eventually produced the eruptive and dome growth acceleration events