Search for: All records

As the volume of data collected at monitored volcanoes continues to expand, researchers will need quick, reliable, and automated methods of inverting those data into useful models of the underlying magma systems. Recently adapted from other fields for use in volcanology, the Ensemble Kalman Filter (EnKF) is one such inversion technique that has been used to produce several successful forecasts and hind‐casts of volcanic unrest, correlating geodetic deformation with mechanical stresses around the magma reservoir. However, given the similarity in which changes to a reservoir's size and pressure are expressed at the surface, the filter can have trouble fully resolving magmatic conditions. In this study, we therefore test several different published variations of the EnKF workflow to produce an optimal configuration for use in future forecasting efforts. By generating synthetic observations of ground deformation under known conditions and then assimilating them through different implementations of the EnKF, we find that many variants favored in other fields underperform for this specific application. We conclude that correlations between model parameters that develop within the EnKF's Monte Carlo ensemble distort the filter's ability to correctly update the model state, causing the filter to systematically favor changes in some parameters over others and ultimately converge to a partially inaccurate solution. This effect can be somewhat mitigated by interrupting these parameter correlations, and the filter remains sensitive to many aspects of the magma system regardless. However, further research and novel approaches will be needed to truly optimize the EnKF for use in volcanology.

 

Surface deformation and seismicity provide critical information to understand the dynamics of volcanic unrest. During 2006–2007, >80 mm/yr uplift was observed by interferometric synthetic aperture radar (InSAR) at the central Atka volcanic center, Alaska, coinciding with an increasing seismicity rate. On November 25, 2006, a phreatic eruption occurred at the Korovin volcanic vent, 5‐km north of the central Atka, following the drainage of its crater lake a month prior to the eruption. The InSAR data are assimilated into three‐dimensional finite element models using the Ensemble Kalman Filter to investigate: (1) the pressure source creating the surface deformation; (2) the triggering of the volcano‐tectonic (VT) earthquakes in the Atka volcanic center; and (3) the triggering of the phreatic eruption at Korovin. The models show that the pressure source required to create the surface deformation is a NE‐tilted, oblate ellipsoid, which rotated from steep to gentle dipping from June to November 2006 before the eruption. The modeled dilatancy in a pre‐existing weak zone, coinciding with the Amlia‐Amukta fault, driven by the pressure source has a spatial and temporal correlation with the evolution of the VT earthquakes during the unrest. The fault dilatancy may have increased the connected porosity and permeability of the fault zone allowing fluid injection which triggered the observed seismicity. In addition, the dilatated fault may have increased the fluid capacity of the fault zone by ∼10⁵ m³, causing the discharge of the crater lake at Korovin. Consequently, the phreatic eruption of the Korovin volcano may have been triggered.

 

Volcanic eruptions pose a significant and sometimes unpredictable hazard, especially at systems that display little to no precursory signals. For example, the 2008 eruption of Okmok volcano in Alaska notably lacked observable short‐term precursors despite years of low‐level unrest. This unpredictability highlights that direct monitoring alone is not always enough to reliably forecast eruptions. In this study, we use the Ensemble Kalman Filter (EnKF) to produce a successful hindcast of the Okmok magma system in the lead up to its 2008 eruption. By assimilating geodetic observations of ground deformation, finite element models track the evolving stress state of the magma system and evaluate its stability using mechanical failure criteria. The hindcast successfully indicates an increased eruption likelihood due to tensile failure weeks in advance of the 2008 eruption. The effectiveness of this hindcast illustrates that EnKF‐based forecasting methods may provide critical information on eruption probability in systems lacking obvious precursors.

 

Forecasting the onset of a volcanic eruption from a closed system requires understanding its stress state and failure potential, which can be investigated through numerical modeling. However, the lack of constraints on model parameters, especially rheology, may substantially impair the accuracy of failure forecasts. Therefore, it is essential to know whether large variations and uncertainties in rock properties will preclude the ability of models to predict reservoir failure. A series of two‐dimensional, axisymmetric models are used to investigate sensitivities of brittle failure initiation to assumed rock properties. The numerical experiments indicate that the deformation and overpressure at failure onset simulated by elastic models will be much lower than the viscoelastic models, when the timescale of pressurization exceeds the viscoelastic relaxation time of the host rock. Poisson's ratio and internal friction angle have much less effect on failure forecasts than Young's modulus. Variations in Young's modulus significantly affect the prediction of surface deformation before failure onset when Young's modulus is < 40 GPa. Longer precursory volcano‐tectonic events may occur in weak host rock (E< 40 GPa) due to well‐developed Coulomb failure prior to dike propagation. Thus, combining surface deformation with seismicity may enhance the accuracy of eruption forecast in these situations. Compared to large and oblate magma systems, small and prolate systems create far less surface uplift prior to failure initiation, suggesting that more frequent measurements are necessary.

 

Extensive vertical deformation (>4.5 m) observed at Sierra Negra volcano Galápagos, Ecuador, between 1992 and the 2005 eruption led scientists to hypothesize that repeated faulting events relieved magma chamber overpressure and prevented eruption. To better understand the catalyst of the 2005 eruption, thermomechanical models are used to track the stress state and stability of the magma storage system during the 1992–2005 inflation events. Numerical experiments indicate that the host rock surrounding the Sierra Negra reservoir remained in compression with minimal changes in overpressure (~10 MPa) leading up to the 2005 eruption. The lack of tensile failure and minimal overpressure accumulation likely inhibited dike initiation and accommodated the significant inflation without the need for pressure relief through shallow trapdoor faulting events. The models indicate that static stress transfer due to the M_w5.4 earthquake 3 hr prior to the eruption most likely triggered tensile failure and catalyzed the 2005 eruption.

 

Numerical modeling framework tracks volcanic activity and forecasts the 2018 eruption of Sierra Negra Volcano, Galapagos. 

Monitoring the activity of subglacial volcanoes along the Aleutian Arc in Alaska is important to the safety of local populations, as well as air traffic flying through the region. However, observations of volcanic unrest are limited by accessibility and resources, particularly at glacier-covered systems, making investigations of their stability challenging. Westdahl Peak, a subglacial volcano on Unimak Island in the Aleutian Arc has experienced significant unrest and uplift since its most recent VEI three eruption in 1991-1992. Given the magnitude of observed uplift, previous investigations suggested the potential for eruption by 2010, but no such event has occurred. One hypothesis to explain this prolonged unrest is that the 1-km thick glacier may increase the stability of the magma system. However, the impact of ice caps and glaciers on the short-term stability of volcanoes is not well understood. In this study, thermomechanical finite element models are used to evaluate how the stability of a glaciated volcano is impacted by variations in ice cap thickness, magma chamber depth, geometry, magma flux rate, and seasonal changes in ice cover thickness. Our numerical experiments indicate that the presence of an ice cap (1–3 km thick) increases the average repose interval for a magma system. Among models with different magma chamber geometries, depths, and flux rates, the greatest increases in repose interval are observed in prolate systems where the increase is up to 57% for a chamber located at 5 km-depth. Spherical and oblate also experience smaller, yet significant, increases in repose interval. Additionally, the percentage increase in repose interval is not impacted by variations in magma flux rate for a given ice cap thickness and magma chamber geometry. However, flux rates do influence the timing of eruptions when the system is experiencing seasonal variations in ice thickness. Our results show that systems with low flux rates are more likely to fail when the ice thickness is at its lowest. The numerical estimates further suggest that the ice cap on Westdahl Peak, which is ∼1 km, may slightly increase the stability of the magma system. In general, given flux rates and magma chamber geometries estimated for the Westdahl system, the repose interval can increase by ∼7 years due to the Westdahl glacier. This increase is small on a geologic scale but is significant on human time scales and the impact of glaciers must be considered in future forecasting efforts. 

Ensemble based data assimilation approaches, such as the Ensemble Kalman Filter (EnKF), have been widely and successfully implemented to combine observations with dynamic forecast models. In this study the EnKF is adapted to assimilate ground deformation observations from interferometric synthetic-aperture radar (InSAR) and GPS into thermomechanical finite element models (FEM) to evaluate volcanic unrest. Two eruption hindcasts are investigated: the 2008 eruption of Okmok volcano, Alaska and the 2018 eruption of Sierra Negra volcano, Galápagos, Ecuador. At Okmok, EnKF forecasts tensile failure and the lateral movement of the magma from a central pressure source in the lead up to its 2008 eruption indicating potential for diking. Alternatively, at Sierra Negra, the EnKF forecasts significant shear failure coincident with a Mw 5.4 earthquake that preceded the 2018 eruption. These successful hindcasts highlight the flexibility and potential of the volcano EnKF approach for near real time monitoring and hazard assessment at active volcanoes worldwide.