Geomorphic expression of rapid Holocene silicic magma reservoir growth beneath Laguna del Maule, Chile.

Large rhyolitic volcanoes pose a hazard, yet the processes and signals foretelling an eruption are obscure. Satellite geodesy has revealed surface inflation signaling unrest within magma reservoirs underlying a few rhyolitic volcanoes. Although seismic, electrical, and potential field methods may illuminate the current configuration and state of these reservoirs, they cannot fully address the processes by which they grow and evolve on geologic time scales. We combine measurement of a deformed paleoshore surface, isotopic dating of volcanism and surface exposure, and modeling to determine the rate of growth of a rhyolite-producing magma reservoir. The numerical approach builds on a magma intrusion model developed to explain the current, decade-long, surface inflation at >20 cm/year. Assuming that the observed 62-m uplift reflects several non-eruptive intrusions of magma, each similar to the unrest over the past decade, we find that ~13 km3 of magma recharged the reservoir at a depth of ~7 km during the Holocene, accompanied by the eruption of ~9 km3 of rhyolite. The long-term rate of magma input is consistent with reservoir freezing and pluton formation. Yet, the unique set of observations considered here implies that large reservoirs can be incubated and grow at shallow depth via episodic high-flux magma injections. These replenishment episodes likely drive rapid inflation, destabilize cooling systems, propel rhyolitic eruptions, and thus should be carefully monitored.

Incremental heating of Bishop Tuff sanidine reveals preeruptive radiogenic Ar and rapid remobilization from cold storage.

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 40Ar/39Ar 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 40Ar/39Ar 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 40Ar/39Ar 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 40Ar/39Ar sanidine and U-Pb zircon dates and a substantially improved capability to resolve the timing and causal relationship of events in the geologic record.