No abstract

No abstract

No abstract

Explosive eruptions of large-volume rhyolitic magma systems are common in the geologic record and pose a major potential threat to society. Unlike other natural hazards, such as earthquakes and tsunamis, a large rhyolitic volcano may provide warning signs long before a caldera-forming eruption occurs. Yet, these signs-and what they imply about magma-crust dynamics-are not well known. This is because we have learned how these systems form, grow, and erupt mainly from the study of ash flow tuffs deposited tens to hundreds of thousands of years ago or more, or from the geophysical imaging of the unerupted portions of the reservoirs beneath the associated calderas. The Laguna del Maule Volcanic Field, Chile, includes an unusually large and recent concentration of silicic eruptions. Since 2007, the crust there has been inflating at an astonishing rate of at least 25 cm/yr. This unique opportunity to investigate the dynamics of a large rhyolitic system while magma migration, reservoir growth, and crustal deformation are actively under way is stimulating a new international collaboration. Findings thus far lead to the hypothesis that the silicic vents have tapped an extensive layer of crystal-poor, rhyolitic melt that began to form atop a magmatic mush zone that was established by ca. 20 ka with a renewed phase of rhyolite eruptions during the Holocene. Modeling of surface deformation, magnetotelluric data, and gravity changes suggest that magma is currently intruding at a depth of ~5 km. The next phase of this investigation seeks to enlarge the sets of geophysical and geochemical data and to use these observations in numerical models of system dynamics.

The Laguna del Maule (LdM) volcanic field comprises the greatest concentration of postglacial rhyolite in the Andes and includes the products of ~40 km3 of 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‐km3 shallow 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.

A warped paleoshoreline records 10,000 years of magma-driven surface deformation above an active rhyolite-producing reservoir.

Accurate and precise ages of large silicic eruptions are critical to calibrating the geologic timescale and gauging the tempo of changes in climate, biologic evolution, and magmatic processes throughout Earth history. The conventional approach to dating these eruptive products using the Ar/Ar method is to fuse dozens of individual feldspar crystals. However, dispersion of fusion dates is common and interpretation is complicated by increasingly precise data obtained via multicollector mass spectrometry. Incremental heating of 49 individual Bishop Tuff (BT) sanidine crystals produces Ar/Ar dates with reduced dispersion, yet we find a 16-ky range of plateau dates that is not attributable to excess Ar. We interpret this dispersion to reflect cooling of the magma reservoir margins below ∼475 °C, accumulation of radiogenic Ar, and rapid preeruption remobilization. Accordingly, these data elucidate the recycling of subsolidus material into voluminous rhyolite magma reservoirs and the effect of preeruptive magmatic processes on the Ar/Ar system. The youngest sanidine dates, likely the most representative of the BT eruption age, yield a weighted mean of 764.8 ± 0.3/0.6 ka (2σ analytical/full uncertainty) indicating eruption only ∼7 ky following the Matuyama-Brunhes magnetic polarity reversal. Single-crystal incremental heating provides leverage with which to interpret complex populations of Ar/Ar sanidine and U-Pb zircon dates and a substantially improved capability to resolve the timing and causal relationship of events in the geologic record.

Despite the hazard posed by explosive silicic eruptions, the magma storage conditions and dynamics that precede these events remain controversial. The Laguna del Maule volcanic field, central Chile, is an exceptional example of postglacial (younger than ca. 20,000 years) rhyolite volcanism and sustained unrest driven by a large, shallow, active silicic magma system. New zircon petrochronologic data reveal that compositionally distinct domains developed concurrently within the Laguna del Maule magma reservoir, which produced two episodes of concentrated rhyolitic eruptions at 23–19 and 8–2 ka. Zircon crystallization ages record 160 kyr of magma emplacement resulting in a several hundreds of cubic kilometers reservoir that has been imaged geophysically. The average magma emplacement rate inferred from the zircon geochronology and tomographically defined magma volume is consistent with those required by thermal models to maintain a shallow silicic system. Ti‐in‐zircon temperatures of crystal cores and rims and hiatuses in crystal growth indicates most of this volume persisted in a near‐solidus state. However, consistent patterns of trace element zoning in crystal interiors and crystallization rates derived from a model of diffusion‐limited zircon growth suggest the erupted rhyolite magma batches originated from long‐lived hot zones of extractable mush embedded within the larger, cool reservoir of rigid mush. These contrasting, coeval magma storage conditions obviate a simple hot versus cold storage dichotomy for large silicic magma systems.