Gradients in H2O, CO2, and exsolved gas in a large‐volume silicic magma system: Interpreting the record preserved in melt inclusions from the Bishop Tuff

Abstract

Infrared spectroscopic analyses of ∼140 melt inclusions in quartz phenocrysts from the zoned Bishop rhyolitic tuff demonstrate that systematic gradients in dissolved magmatic H 2 O and CO 2 concentrations were present during preemptive crystallization of the magma body. Melt inclusions from the earliest erupted samples contain lower H 2 O (5.3±0.4 wt %) and CO 2 (62±37 ppm) than inclusions from the middle of the eruption (5.7±0.2 wt % H 2 O; 120±60 ppm CO 2). Melt inclusions from late erupted samples have much lower H 2 O (4.1±0.3 wt %) and higher and variable CO 2 (150–1085 ppm). Trace element analyses of melt inclusions by ion microprobe show that inclusions within single pumice clasts from the early and middle Bishop Tuff have an inverse correlation between CO 2 and incompatible elements. This pattern indicates that the magma was gas‐saturated during crystallization, with CO 2 partitioning into a coexisting gas phase. Quantitative modeling using H 2 O‐CO 2 solubility relations reveals a preeruptive gradient in exsolved gas, with gas contents varying from ∼1 wt % in the deeper regions of the magma body to nearly 6 wt % near the top. Dissolved Cl, B, Li, and Be in melt inclusions correlate negatively with CO 2. Mass balance modeling of Cl loss to exsolving H 2 O‐rich gas during crystallization provides strong corroborating evidence for the mass fractions of exsolved gas estimated from H 2 O, CO 2, and trace element data. Pressures of quartz crystallization and melt inclusion entrapment calculated from inclusion H 2 O‐CO 2 data are consistent with progressive downward tapping of a zoned magma body during the eruption. Melt inclusion gas saturation pressures, magma volume estimates, and time‐stratigraphic‐compositional relations suggest that early erupted magma was stored at the top of a downward widening magma body. Melt inclusion data and the inferred gradients in dissolved H 2 O, CO 2 and exsolved gas in the Bishop magma body suggest that gas saturation plays an important role in the formation and subsequent preservation of compositional gradients in silicic magma reservoirs.