Title: Geodetic Data Assimilation for Evaluating Volcanic Unrest
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.  more » « less
Award ID(s):
1634995 1752477
PAR ID:
10272179
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
IEEE IGARSS
Page Range / eLocation ID:
3618 to 3621
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Numerical modeling framework tracks volcanic activity and forecasts the 2018 eruption of Sierra Negra Volcano, Galapagos. 
    more » « less
  2. Abstract The 2018 eruption of Sierra Negra volcano, Galápagos, Ecuador has provided new insights into the mechanisms of caldera resurgence, subsidence, and fissuring at basaltic shield volcanoes. Here, we integrate local (∼0.4 km) seismo‐acoustic records and regional (∼85 km) infrasound array data to present new observations of the 2018 Sierra Negra eruption with improved time and spatial resolutions. These observations include: air‐to‐ground coupling ∼2 hr before the time of the eruption onset, migration of the infrasound tremor from 22:54 June 26 to 12:31 June 27 UT (all times in UT), and persistent infrasound detections during the weeks between 5 July and 18 August from an area that does not coincide with the previously documented eruptive fissures. We interpret air‐to‐ground coupling as infrasound tremor generated in the nearby fissures before the main eruptive phase started, although ambiguity remains in interpreting a single seismic‐infrasonic sensor pair. The progressive location change of the infrasound tremor agrees with the migration of the eruption down the north flank of Sierra Negra at a rate of ∼0.15 ± 0.04 m/s. The weeks‐long persistent detections coincide with a region that has thermal anomalies, co‐eruptive deformation, lava fields, and geological features that could be interpreted as multiple lava tube skylights. Our observations and interpretations provide constraints on the mechanisms underlying fissure formation and magma emplacement at Sierra Negra. 
    more » « less
  3. 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 Mw 5.4 earthquake 3 hr prior to the eruption most likely triggered tensile failure and catalyzed the 2005 eruption. 
    more » « less
  4. Abstract In their article entitled “Trapdoor Fault Activation:A Step Toward Caldera Collapse at Sierra Negra,Galapagos,Ecuador” Shreve and Delgado (2023,https://doi.org/10.1029/2023jb026437) examine co‐eruptive deformation during the 2018 eruption of Sierra Negra Volcano. One of their major conclusions is that the 2018 eruption, and specifically co‐eruptive faulting, represents the initial stages of caldera collapse. They reach this conclusion because they focus their analysis solely on co‐eruptive deformation, and do not investigate the total (net) deformation for the 2005 to 2018 eruption cycle. Bell, La Famina, et al. (2021,https://doi.org/10.1038/s41467‐021‐21596‐4) investigated both the pre‐ and co‐eruptive phases of the 2018 eruption and showed that net deformation was one of caldera resurgence, not subsidence. In this comment, we demonstrate that the conclusion of collapse, or even initiation of collapse, is attributable to not accounting for pre‐eruptive deformation on the intra‐caldera Trapdoor Fault system and incorrectly assuming that the volcano‐tectonic dynamics of Sierra Negra mimic those of other basaltic calderas. 
    more » « less
  5. Cinder cones are a common feature at many volcanic eruptions. Their shapes and volumes can reveal information about eruption conditions, and their geomorphological evolution shapes them and their surrounding environment. It is thus important to quantify the rate and patterns of erosion of young cinder cones. In this study, we examine the Ahmanilix cone, which formed during the 2008 eruption of Okmok volcano in the Aleutian islands region of Alaska. Ahmanilix, located on the eastern side of Okmok’s large caldera, is >250 meters tall and characterized by dramatic gullies formed by the harsh wind, snow and rain conditions typical of the Aleutians. We usd photogrammetry to create 3D models of Ahmanilix using aerial photographic surveys taken from a helicopter in 2021, 2022, 2023 and 2024. We utilize Agisoft Metashape to build point clouds, Cloud Compare to align the point clouds and build raster Digital Elevation Models (DEMs), and QGIS and Python to visualize and analyze these products. By subtracting DEM rasters we quantify year-to-year erosion. We compare our results with erosion rates estimated from satellite observations (Dai et al., 2020), identify regions dominated by erosion or deposition and correlate them with slopes and cinder lithology. Our observations can be extended to other cinder cones and help predict their geomorphological evolution. 
    more » « less