skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Modeling Magma System Evolution During 2006–2007 Volcanic Unrest of Atka Volcanic Center, Alaska
Abstract 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 ∼105 m3, causing the discharge of the crater lake at Korovin. Consequently, the phreatic eruption of the Korovin volcano may have been triggered.  more » « less
Award ID(s):
1752477
PAR ID:
10374527
Author(s) / Creator(s):
 ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Solid Earth
Volume:
126
Issue:
2
ISSN:
2169-9313
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    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
  2. Abstract The Kolumbo submarine volcano in the southern Aegean (Greece) is associated with repeated seismic unrest since at least two decades and the causes of this unrest are poorly understood. We present a ten‐month long microseismicity data set for the period 2006–2007. The majority of earthquakes cluster in a cone‐shaped portion of the crust below Kolumbo. The tip of this cone coincides with a low Vp‐anomaly at 2–4 km depth, which is interpreted as a crustal melt reservoir. Our data set includes several earthquake swarms, of which we analyze the four with the highest events numbers in detail. Together the swarms form a zone of fracturing elongated in the SW‐NE direction, parallel to major regional faults. All four swarms show a general upward migration of hypocenters and the cracking front propagates unusually fast, compared to swarms in other volcanic areas. We conclude that the swarm seismicity is most likely triggered by a combination of pore‐pressure perturbations and the re‐distribution of elastic stresses. Fluid pressure perturbations are induced likely by obstructions in the melt conduits in a rheologically strong layer between 6 and 9 km depth. We conclude that the zone of fractures below Kolumbo is exploited by melts ascending from the mantle and filling the crustal melt reservoir. Together with the recurring seismic unrest, our study suggests that a future eruption is probable and monitoring of the Kolumbo volcanic system is highly advisable. 
    more » « less
  3. Abstract Volcanic summit craters are typically noted to form by roof collapse into a depressurized magma chamber or by explosive excavation. Recent examples of effusive activity (e.g., Kilauea Volcano, Hawai'i) allowed specifically for quantification of the collapse process. However, small spatiotemporal morphologic change related to background mass wasting and low‐level explosive activity has not been well quantified in volcanic craters. Telica volcano, Nicaragua, is a persistently restless basaltic‐andesite stratovolcano. Telica's persistent restlessness is caused by a long‐lived magmatic‐hydrothermal system with high‐temperature crater fumaroles and low‐frequency seismicity, punctuated by subdecadal, low‐explosivity (VEI 1–2) phreatic eruptions. We use photographic observations (1994 to 2017) and structure‐from‐motion point cloud construction and differencing (2011 to 2017) to analyze changes at Telica in the context of summit crater formation and eruptive precursors. Crater wall retreat (up to 40 m) spatially correlates with long‐lived high‐temperature fumaroles in the crater walls, whereas eruptions eject material (>5 m) from the crater floor through vent formation and/or clearing. These processes sustain a morphology similar to that of pit craters but without a shallow depressurized magma chamber. Our observations indicate system‐wide sealing prior to eruption by viscous magma in the conduit and eruption of a dome in 2017 and hydrothermal mineralization, not from vent covering talus; though, vent covering talus can redirect the shallow conduit. This study shows promise for photogrammetric techniques in correlating morphologic change with summit crater formation and volcanic activity and the power of long‐term visual observations in understanding active volcanic processes. 
    more » « less
  4. Rock strength has long been linked to lithospheric deformation and seismicity. However, independent constraints on the related elastic heterogeneity are missing, yet could provide key information for solid Earth dynamics. Using coseismic Global Navigation Satellite Systems (GNSS) data for the 2011 M9 Tohoku-oki earthquake in Japan, we apply an inverse method to infer elastic structure and fault slip simultaneously. We find compliant material beneath the volcanic arc and in the mantle wedge within the partial melt generation zone inferred to lie above ~100 km slab depth. We also identify low-rigidity material closer to the trench matching seismicity patterns, likely associated with accretionary wedge structure. Along with traditional seismic and electromagnetic methods, our approach opens up avenues for multiphysics inversions. Those have the potential to advance earthquake and volcano science, and in particular once expanded to InSAR type constraints, may lead to a better understanding of transient lithospheric deformation across scales. 
    more » « less
  5. Volcano monitoring and eruption forecasting require accurate characterization of transcrustal magmatic structures to place volcanic unrest in context within the system where it occurs. Structural imaging using local seismicity is limited to seismogenic depths. Here, we exploit arrivals in teleseismic receiver functions that change polarity with backazimuth to image two surfaces beneath Akutan volcano in the Aleutian arc. The two surfaces delineate an upper to midcrustal inverted conical volume that deepens and thickens away from the volcanic center, with thicknesses of 3–13 km. The top of the volume is at depths of 2–3 km below sea level at distances of ∼5–15 km from the caldera center. The bottom is at depths of 7–15 km at the same distances, and the cone’s thickness increases outward from ∼5 to ∼10 km. The signal is best fit by a volume with anisotropy with fast symmetry planes that dip outward from the center and downward increases in shear velocity at both interfaces. The upper boundary coincides with the top of Akutan’s volcanotectonic (VT) seismogenic zone, with the VT seismicity exhibiting outward dipping planar features that match the anisotropic fast plane orientation within the volume. The bottom of the anisotropic volume is below the termination depth of the majority of the VT seismicity and is therefore likely associated with the brittle–ductile transition. Long-period (LP) events associated previously with magma movement are concentrated below the anisotropic VT volume. Because of the strong spatial association with VT seismicity, we interpret the volume as consisting of concentric outward dipping faults and dikes that align the seismogenic response to stress changes from magmatic processes. Our observations map this volume independent of the present-day seismicity distribution and thus provide a spatially more complete image of the magmatic system. 
    more » « less