The Laguna del Maule (LdM) volcanic field comprises the greatest concentration of postglacial rhyolite in the Andes and includes the products of ~40 km3of explosive and effusive eruptions. Recent observations at LdM by interferometric synthetic aperture radar and global navigation satellite system geodesy have revealed inflation at rates exceeding 20 cm/year since 2007, capturing an ongoing period of growth of a potentially large upper crustal magma reservoir. Moreover, magnetotelluric and gravity studies indicate the presence of fluids and/or partial melt in the upper crust near the center of inflation. Petrologic observations imply repeated, rapid extraction of rhyolitic melt from crystal mush stored at depths of 4–6 km during at least the past 26 ka. We utilize multiple types of surface‐wave observations to constrain the location and geometry of low‐velocity domains beneath LdM. We present a three‐dimensional shear‐wave velocity model that delineates a ~450‐km3shallow magma reservoir ~2 to 8 km below surface with an average melt fraction of ~5%. Interpretation of the seismic tomography in light of existing gravity, magnetotelluric, and geodetic observations supports this model and reveals variations in melt content and a deeper magma system feeding the shallow reservoir in greater detail than any of the geophysical methods alone. Geophysical imaging of the LdM magma system today is consistent with the petrologic inferences of the reservoir structure and growth during the past 20–60 kyr. Taken together with the ongoing unrest, a future rhyolite eruption of at least the scale of those common during the Holocene is a reasonable possibility.
- Award ID(s):
- 1757495
- NSF-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
-
more » « less
Abstract -
more » « less
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
Abstract Imaging silicic systems using geophysics is challenging because many interrelated factors (e.g., temperature, melt fraction, melt composition, geometry) can contribute to the measured geophysical anomaly. Joint interpretation of models from multiple geophysical methods can better constrain interpretations of the subsurface structure. Previously published resistivity and shear wave velocity (Vs) models, derived separately from magnetotelluric (MT) and surface wave seismic data, respectively, have been used to model the restless Laguna del Maule Volcanic Field, central Chile. The Vs model contains a 450 km3low‐velocity zone (LVZ) interpreted as a region with an average melt fraction of 5–6%. The resistivity model contains a conductor (C3) interpreted as a region with a melt fraction >35%. The spatial extents of the LVZ and C3 overlap, but the geometries and interpretations of these features are different. To resolve these discrepancies, this study investigates the resolution of the MT data using hypothesis testing and constrained MT inversions. It is shown that the MT data are best fit with discrete conductors embedded within the larger LVZ. The differences between the MT and seismic models reflect resolution differences between the two data sets as well as varying sensitivities to physical properties. The MT data are sensitive to smaller volumes of extractable mush that contain well‐connected crystal‐poor melt (C3). The seismic data have lower spatial resolution but image the full extent of the poorly connected crystal‐rich magma storage system. The combined images suggest that the LdMVF magma plumbing system is thermally heterogeneous with coexisting zones of warm and cold storage.
-
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
-
more » « less
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.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.3389/feart.2021.725917
