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
Dissecting a Zombie: Joint Analysis of Density and Resistivity Models Reveals Shallow Structure and Possible Sulfide Deposition at Uturuncu Volcano, Bolivia
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
- Award ID(s):
- 1757495
- PAR ID:
- 10323277
- Date Published:
- Journal Name:
- Frontiers in Earth Science
- Volume:
- 9
- ISSN:
- 2296-6463
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
SUMMARY The ability to accurately and reliably obtain images of shallow subsurface anomalies within the Earth is important for hazard monitoring and a fundamental understanding of many geologic structures, such as volcanic edifices. In recent years, machine learning (ML) has gained increasing attention as a novel approach for addressing complex problems in the geosciences. Here we present an ML-based inversion method to integrate cosmic-ray muon and gravity data sets for shallow subsurface density imaging at a volcano. Starting with an ensemble of random density anomalies, we use physics-based forward calculations to find the corresponding set of expected gravity and muon attenuation observations. Given a large enough ensemble of synthetic density patterns and observations, the ML algorithm is trained to recognize the expected spatial relations within the synthetic input–output pairs, learning the inherent physical relationships between them. Once trained, the ML algorithm can then interpolate the best-fitting anomalous pattern given data that were not used in training, such as those obtained from field measurements. We test the validity of our ML algorithm using field data from the Showa-Shinzan lava dome (Mt Usu, Japan) and show that our model produces results consistent with those obtained using a more traditional Bayesian joint inversion. Our results are similar to the previously published inversion, and suggest that the Showa-Shinzan lava dome consists of a relatively high-density (2200–2400 km m–3) cylindrical anomaly, about 300 m in diameter. Adding noise to synthetic training and testing data sets shows that, as expected, the ML algorithm is most robust in areas of high sensitivity, as determined by the forward kernels. Overall, we discover that ML offers a viable alternate method to a Bayesian joint inversion when used with gravity and muon data sets for subsurface density imaging.more » « less
-
Abstract Enigmatic large‐scale (>150 km wide) ground deformation in southern Bolivia has been ongoing for more than 50 year. Concurrent changes in gravity recorded between 2010 and 2018 imply minor changes in subsurface density in the absence of significant mass changes. Numerical modeling of the gravity changes and concurrent InSAR LOS displacements gives annual bulk density changes of 0.002 kg m−3in the Altiplano‐Puna Magma Body (APMB) and −0.03 kg m−3in a vertical bulge‐column ensemble beneath Uturuncu volcano. We propose that the transcrustal migration of fluids from the APMB to shallower crustal levels by compressible flow is the source of ground deformation. Localized ground subsidence south of Uturuncu can be best explained by a density decrease of 20 ± 5 kg m−3between 2011 and 2013 in a hydrothermal reservoir. Our findings contribute to the growing recognition of transcrustal fluid migration as a source of volcanic unrest.more » « less
-
Santorini volcano in the South Aegean Volcanic Arc has a detailed history of ongoing volcanic and seismic activity, making it a prime location for studying magma storage and transport at arc volcanoes. The shallow magmatic system (<5 km depth) is well constrained by geophysical studies, but the deeper crustal structure is not. Located 15 km NE of Santorini, the Kolumbo seamount is also an active edifice, with consistently more seismicity and hydrothermal venting than Santorini. Geochemical studies indicate that Santorini and Kolumbo are fed by separate mantle and crustal magma sources, but prior seismic studies suggest otherwise (Dimitriadis et al, 2010; McVey et al, 2020). This study addresses the nature of lower-crustal magma structure beneath arc volcanoes and whether evolved volcanoes and nearby vents are connected through their plumbing. Tomographic inversion of P-wave Moho reflection (PmP) and turning P-wave (Pg) traveltimes is used to create 3-D models of Moho depth and P-wave velocity (Vp) down to depths of ~25 km. The PROTEUS experiment provides an exceptionally dense and large aperture traveltime dataset from an amphibious array of ~150 seismometers and ~14,000 active marine sources. The data are ~33,000 manually picked PmP arrivals and ~256,000 Pg arrivals from existing studies. Results show a low Vp anomaly extending from the Moho to the surface. This anomaly starts at the base of the crust under the NW Santorini caldera and extends up to the east. It is most pronounced at 10-15 km depth, where it is offset from both Santorini and Kolumbo. Limited resolution prevents imaging of a connection between this mid-crustal anomaly and the known shallow magma storage region under the Santorini caldera. A high-velocity core beneath Santorini is not found, a feature interpreted at other volcanoes as a cooled intrusive complex. Because no additional low Vp anomalies are found in the lower crust, we infer that a common mantle source and mid-crustal plumbing system is actively feeding both Santorini and Kolumbo. The spatial offset and elongated nature of magma storage implies a complex relationship between evolving magmatic structures and tectonics.more » « less
-
null (Ed.)Santorini is located in the central part of the Hellenic Volcanic Arc (South Aegean Sea) and is well known for the Late-Bronze-Age “Minoan” eruption that may have been responsible for the decline of the great Minoan civilization on the island of Crete. To use gravity to probe the internal structure of the volcano and to determine whether there are temporal variations in gravity due to near surface changes, we construct two gravity maps. Dionysos Satellite Observatory (DSO) of the National Technical University of Athens (NTUA) carried out terrestrial gravity measurements in December 2012 and in September 2014 at selected locations on Thera, Nea Kameni, Palea Kameni, Therasia, Aspronisi and Christiana islands. Absolute gravity values were calculated using raw gravity data at every station for all datasets. The results were compared with gravity measurements performed in July 1976 by DSO/NTUA and absolute gravity values derived from the Hellenic Military Geographical Service (HMGS) and other sources. Marine gravity data that were collected during the PROTEUS project in November and December 2015 fill between the land gravity datasets. An appropriate Digital Elevation Model (DEM) with topographic and bathymetric data was also produced. Finally, based on the two combined datasets (one for 2012–2014 and one for the 1970s), Free air and complete Bouguer gravity anomaly maps were produced following the appropriate data corrections and reductions. The pattern of complete Bouguer gravity anomaly maps was consistent with seismological results within the caldera. Finally from the comparison of the measurements made at the same place, we found that, within the caldera, the inner process of the volcano is ongoing both before, and after, the unrest period of 2011–2012.more » « less
An official website of the United States government

