Interpreting unrest at silicic volcanoes requires knowledge of the magma storage conditions and dynamics that precede eruptions. The Laguna del Maule volcanic field, Chile, has erupted ~40 km3of rhyolite over the last 20 ka. Astonishing rates of sustained surface inflation at >25 cm/year for >12 years reveal a large, restless system. Integration of geochronologic, petrologic, geomorphic, and geophysical observations provides an unusually rich context to interpret ongoing and prehistoric processes. We present new volatile (H2O, CO2, S, F, and Cl), trace element, and major element concentrations from 109 melt inclusions hosted in quartz, plagioclase, and olivine from seven eruptions. Silicic melts contain up to 8.0 wt. % H2O and 570 ppm CO2. In rhyolites melt inclusions track decompression‐driven fractional crystallization as magma ascended from ~14 to 4 km. This mirrors teleseismic tomography and magnetotelluric findings that reveal a domain containing partial melt spanning from 14 to 4 km. Ce and Cl contents of rhyolites support the generation of compositionally distinct domains of eruptible rhyolite within the larger reservoir. Heat, volatiles, and melt derived from episodic mafic recharge likely incubate and grow the shallow reservoir. Olivine‐hosted melt inclusions in mafic tephra contain up to 2.5 wt. % H2O and 1,140 ppm CO2and proxy for the volatile load delivered via recharge into the base of the silicic mush at ~14 to 8 km. We propose that mafic recharge flushes deeper reaches of the magma reservoir with CO2that propels H2O exsolution, upward accumulation of fluid, pressurization, and triggering of rhyolitic eruptions.
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.
more » « less- Award ID(s):
- 1411779
- NSF-PAR ID:
- 10371341
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Solid Earth
- Volume:
- 124
- Issue:
- 3
- ISSN:
- 2169-9313
- Page Range / eLocation ID:
- p. 2858-2872
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
more » « less
Abstract With substantial postglacial rhyolite eruptions and ongoing rapid uplift, the Laguna del Maule volcanic field in the southern Andes provides an exceptional opportunity to study the dynamics of an active silicic magmatic system. Using 4,093
P arrivals from 137 distant earthquakes recorded by 44 local stations over∼ 2.25 years, we conduct teleseismic tomography to image the crustal structure down to 40 km below the volcanic field. A prominent low‐velocity body at depths between∼ 0 and 12 km below sea level (b.s.l.), characterized by a volume of∼ 500 km3 and a peak anomaly of− 400 m/s (∼ 9%), overlaps the location of the upper‐crustal magma reservoir detected in recent gravity and surface wave tomography studies. Its estimated averageP wave velocity of∼ 4.6 km/s corresponds to an average melt fraction of about 14% and a melt volume of∼ 70 km3 . Petrologic observations are also consistent with generation and storage of rhyolitic melts at depths corresponding to the anomalous zone. Moreover, the tomographic results support a lower crust zone of MASH (melting, assimilation, storage, and homogenization) from a depth of∼ 25 km to the base of the model, which likely reflects a deep crustal source of magma that contributes to and incubates the shallow silicic reservoir. -
more » « less
Abstract The rhyolite‐producing Laguna del Maule volcanic field (LdMVF), Chile, has had numerous post‐glacial eruptions that produced large explosions and voluminous lava flows. During the Holocene ∼60 m of surface uplift is recorded by paleo‐shorelines of the fresh‐water Laguna del Maule, with an inflation source near the Barrancas volcanic complex. Rhyolites from the Barrancas complex erupted over ∼14 ka including some of the youngest (1.4 ± 0.6 ka) lava flows in the field. New gravity data collected on the Barrancas complex reveals a residual gravity low (−6 mGal, “Barrancas anomaly”) that is distinct from the pronounced gravity low (−19 mGal; “Lake anomaly”) associated with present‐day ground uplift to the northwest. Three‐dimensional inversion of the Barrancas anomaly indicates the presence of a magma body with a maximum density contrast with the host rock of −250 kg/m3centered at a depth of ∼3 km below surface. Nearby Miocene high‐silica granites represent frozen remnants of highly evolved rhyolitic magma. Comparison of the densities measured from samples of these plutons with the geophysical model densities, and integration of thermodynamic modeling of silicic melt evolution, provide constraints on our interpretation. We propose a magma body, containing <30% melt phase and low volatile content, exists beneath Barrancas. The Barrancas and Lake gravity lows represent magma in different physical states, associated with past and present‐day storage beneath LdMVF. The gravity model mirrors geochemical observations which independently indicate that at least two distinct rhyolites were generated and stored as discrete magma bodies within the broader LdMVF.
-
more » « less
SUMMARY The 2011–2012 eruption at Cordón Caulle in Chile produced crystal-poor rhyolitic magma with crystal-rich mafic enclaves whose interstitial glass is of identical composition to the host rhyolite. Eruptible rhyolites are thought to be genetically associated with crystal-rich magma mushes, and the enclaves within the Cordón Caulle rhyolite support the existence of a magma mush from which the erupted magma was derived. Moreover, towards the end of the 2011–2012 eruption, subsidence gave way to inflation that has on average been continuous through at least 2020. We hypothesize that magma segregation from a crystal mush could be the source of the observed inflation. Conceptually, magma withdrawal from a crystal-poor rhyolite reservoir caused its depressurization, which could have led to upward flow of interstitial melt within an underlying crystal mush, causing a new batch of magma to segregate and partially recharge the crystal-poor rhyolite body. Because the compressibility of the crystalline matrix of the mush is expected to be lower than that of the interstitial melt, which likely contains some fraction of volatile bubbles, this redistribution of melt would result in a net increase in volume of the system and in the observed inflation. We use numerical modelling of subsurface magma flow and storage to show under which conditions such a scenario is supported by geodetic and petrologic observations.
-
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.
