skip to main content


Title: The Impact of Ice Caps on the Mechanical Stability of Magmatic Systems: Implications for Forecasting on Human Timescales
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.  more » « less
Award ID(s):
1752477
NSF-PAR ID:
10342367
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Frontiers in Earth Science
Volume:
10
ISSN:
2296-6463
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. The recent identification of unrest at multiple volcanoes that have not erupted in over 10 kyr presents an intriguing scientific problem. How can we distinguish between unrest signaling impending eruption after kyr of repose and non-magmatic unrest at a waning volcanic system? After ca. 250 kyr without a known eruption, in recent decades Uturuncu volcano in Bolivia has exhibited multiple signs of unrest, making the classification of this system as “active”, “dormant”, or “extinct” a complex question. Previous work identified anomalous low resistivity zones at <10 km depth with ambiguous interpretations. We investigate subsurface structure at Uturuncu with new gravity data and analysis, and compare these data with existing geophysical data sets. We collected new gravity data on the edifice in November 2018 with 1.5 km spacing, ±15 μ Gal precision, and ±5 cm positioning precision, improving the resolution of existing gravity data at Uturuncu. This high quality data set permitted both gradient analysis and full 3-D geophysical inversion, revealing a 5 km diameter, positive density anomaly beneath the summit of Uturuncu (1.5–3.5 km depth) and a 20 km diameter arc-shaped negative density anomaly around the volcano (0.5–7.5 depth). These structures often align with resistivity anomalies previously detected beneath Uturuncu, although the relationship is complex, with the two models highlighting different components of a common structure. Based on a joint analysis of the density and resistivity models, we interpret the positive density anomaly as a zone of sulfide deposition with connected brines, and the negative density arc as a surrounding zone of hydrothermal alteration. Based on this analysis we suggest that the unrest at Uturuncu is unlikely to be pre-eruptive. This study shows the value of joint analysis of multiple types of geophysical data in evaluating volcanic subsurface structure at a waning volcanic center. 
    more » « less
  2. Abstract

    The along‐strike variations of the velocity, thickness, and dip of subducting slabs and the volcano distribution have been observed globally. It is, however, unclear what controls the distribution of volcanoes and the associated magma generation. With the presence of nonuniform volcanism, the Aleutian‐Alaska subduction zone (AASZ) is an ideal place to investigate subduction segmentation and its relationship with volcanism. Using full‐wave ambient noise tomography, we present a high‐resolution 3‐D shear wave velocity model of the AASZ for the depths of 15–110 km. The velocity model reveals the distinct high‐velocity Pacific slab, the thicker, flatter, and more heterogeneous Yakutat slab, and the northeasterly dipping Wrangell slab. We observe low velocities within the uppermost mantle (at depth <60 km) below the Aleutian arc volcanoes, representing partial melt accumulation. The large crustal low‐velocity anomaly beneath the Wrangell volcanic field suggests a large magma reservoir, likely responsible for the clustering of volcanoes. The Denali volcanic gap is above an average‐velocity crust but an extremely fast mantle wedge, suggesting the lack of subsurface melt. This is in contrast with the lower‐velocity back‐arc mantle beneath the adjacent Buzzard Creek‐Jumbo Dome volcanoes to the east. The back‐arc low velocities associated with the Pacific, the eastern Yakutat, and the Wrangell slabs may reflect subduction‐driven mantle upwelling. The structural variation of the downgoing slabs and the overriding plate explains the change of volcanic activity along the AASZ. Our findings demonstrate the combined role of the subducting slab and the overriding plate in controlling the characteristics of arc magmatism.

     
    more » « less
  3. Abstract

    Since the 1919 foundation of the International Association of Volcanology and Chemistry of the Earth’s Interior (IAVCEI), the fields of volcano seismology and acoustics have seen dramatic advances in instrumentation and techniques, and have undergone paradigm shifts in the understanding of volcanic seismo-acoustic source processes and internal volcanic structure. Some early twentieth-century volcanological studies gave equal emphasis to barograph (infrasound and acoustic-gravity wave) and seismograph observations, but volcano seismology rapidly outpaced volcano acoustics and became the standard geophysical volcano-monitoring tool. Permanent seismic networks were established on volcanoes (for example) in Japan, the Philippines, Russia, and Hawai‘i by the 1950s, and in Alaska by the 1970s. Large eruptions with societal consequences generally catalyzed the implementation of new seismic instrumentation and led to operationalization of research methodologies. Seismic data now form the backbone of most local ground-based volcano monitoring networks worldwide and play a critical role in understanding how volcanoes work. The computer revolution enabled increasingly sophisticated data processing and source modeling, and facilitated the transition to continuous digital waveform recording by about the 1990s. In the 1970s and 1980s, quantitative models emerged for long-period (LP) event and tremor sources in fluid-driven cracks and conduits. Beginning in the 1970s, early models for volcano-tectonic (VT) earthquake swarms invoking crack tip stresses expanded to involve stress transfer into the wall rocks of pressurized dikes. The first deployments of broadband seismic instrumentation and infrasound sensors on volcanoes in the 1990s led to discoveries of new signals and phenomena. Rapid advances in infrasound technology; signal processing, analysis, and inversion; and atmospheric propagation modeling have now established the role of regional (15–250 km) and remote (> 250 km) ground-based acoustic systems in volcano monitoring. Long-term records of volcano-seismic unrest through full eruptive cycles are providing insight into magma transport and eruption processes and increasingly sophisticated forecasts. Laboratory and numerical experiments are elucidating seismo-acoustic source processes in volcanic fluid systems, and are observationally constrained by increasingly dense geophysical field deployments taking advantage of low-power, compact broadband, and nodal technologies. In recent years, the fields of volcano geodesy, seismology, and acoustics (both atmospheric infrasound and ocean hydroacoustics) are increasingly merging. Despite vast progress over the past century, major questions remain regarding source processes, patterns of volcano-seismic unrest, internal volcanic structure, and the relationship between seismic unrest and volcanic processes.

     
    more » « less
  4. Abstract

    Improvement of our understanding of the role of ground deformation due to viscoelastic relaxation following eruptions is important, as the generated signals can resemble renewed magma inflow. We study post‐eruptive unrest at the subglacial Bárðarbunga volcano, Iceland, after a caldera collapse and major magma drainage in 2014–2015. Elevated seismicity began about 6 months after the eruption ended, including nine Mlw> 4.5 earthquakes. Global Navigation Satellite System and Sentinel‐1 Interferometric Synthetic Aperture Radar geodesy are applied to evaluate post‐eruptive ground deformation from 2015 to 2018. Horizontal velocities locally exceed 10 cm/year and rapidly decay with distance away from the caldera. We explore two end‐member models and their combination to explain the post‐eruptive deformation field: 1) viscoelastic relaxation caused by the co‐eruptive caldera collapse and magma withdrawal, and 2) renewed magma inflow. We find parameter combinations for each model that explain the observed ground deformation. The purely viscoelastic relaxation model, consisting of a half‐space composed of a 7‐km thick elastic layer on top of a viscoelastic layer with a viscosity of 3.0 × 1018 Pa s reproduces broadly the observations. A simple magma inflow model consisting of a single point source with an inflow rate of 1 × 107 m3/year at 0.7 km depth broadly fits the observations, but may be unrealistic. A more elaborate model of magma inflow into a 10‐km deep sill combined with slip on the caldera ring fault explains the observations well. Our results suggest that the co‐eruptive deformation field is likely influenced by viscoelastic relaxation, renewed magma inflow, or a combination of both processes.

     
    more » « less
  5. null (Ed.)
    Almost two-thirds of Earth’s submarine volcanism is expressed in the ocean basins along the ~65,000 km long mid-ocean-ridge (MOR) system, but the remaining third takes place along intraoceanic arcs and seamounts. The 22,000 km long intraoceanic arc system appears to surpass the MOR system in terms of both the frequency of hydrothermal activity, with several sites per 100 km of arc, and the range in chemical composition of the fluids being discharged. Hydrothermal systems hosted by submarine arc volcanoes differ substantially from those in spreading environments in that they commonly contain a large component of magmatic fluid. Our primary scientific goal is to discover the fundamental underlying processes that develop as a consequence of this difference. A magmatic-hydrothermal signature, coupled with the shallow depths of arc volcanoes and their high-volatile contents, strongly influences the chemistry of fluids and the resulting mineralization and likely has important consequences for the biota associated with these systems. Because of the high metal contents and very acidic fluids, these hydrothermal systems are also thought to be important analogs of many porphyry copper and epithermal gold deposits mined today on land. Drilling at Brothers volcano on the Kermadec arc will provide the missing link (i.e., the third dimension) in our understanding of mineral deposit formation along arcs, the subseafloor architecture of these volcanoes and their related permeability, as well as the relationship between the discharge of magmatic fluids and the deep biosphere. Expedition 376 will drill and log two caldera sites, one on the rim of the caldera and a second inside the caldera, and a third site at the summit of a volcanic cone, located inside the caldera at Brothers. These sites represent discharge zones of geochemically distinct fluids that are variably affected by magmatic volatile input, allowing us to directly address the consequences of magma degassing for metal transport to the seafloor and its effect on the functioning of microbial communities. Drilling will provide access to critical zones dominated by magma degassing and high-temperature hydrothermal circulation over depth intervals regarded as crucial, not only in the development of multiphase mineralizing systems but also in identifying subsurface microbially habitable environments. The specific objectives of Expedition 376 are • To characterize the subvolcano, magma chamber–derived volatile phase to test model-based predictions that this is either a single-phase gas or two-phase brine-vapor; • To determine the subseafloor distribution of base and precious metals and metalloids and the reactions that have taken place along pathways to the seafloor; • To quantify the mechanisms and extent of fluid-rock interaction and consequences for mass transfer of metals and metalloids into the ocean and the role of magmatically derived carbon and sulfur species in mediating these fluxes; and • To assess the diversity, extent, and metabolic pathways of microbial life in an extreme, metal-toxic, and acidic volcanic environment. 
    more » « less