Implications of Multiple Disequilibrium Textures in Quartz-Hosted Embayments

The faces of volcanic phenocrysts may be marked by imperfections occurring as holes that penetrate the crystal interior. When filled with glass these features, called embayments or reentrants, have been used to petrologically constrain magmatic ascent rate. Embayment ascent speedometry relies on the record of disequilibrium preserved as diffusion-limited volatile concentration gradients in the embayment glass. Clear, glassy embayments are carefully selected for speedometry studies. The use and subsequent descriptions of pristine embayments overrepresent their actual abundance. Here, we provide a textural analysis of the number, morphology, and filling characteristics of quartz-hosted embayments. We target a collection of large (i.e., >20 km3 erupted volume) silicic eruptions, including the Bishop Tuff, Tuff of Bluff Point, Bandelier Tuff, Mesa Falls Tuff, and Huckleberry Ridge Tuff in the United States, Oruanui Tuff in New Zealand, Younger Toba Tuff in Indonesia, the Kos Plateau Tuff in Greece, and the Giant Pumice from La Primavera caldera in Mexico. For each unit, hundreds of quartz crystals were picked and the total number of embayment-hosting crystals were counted and categorized into classifications based on the vesicularity and morphology. We observed significant variability in embayment abundance, form, and vesicularity across different eruptions. Simple, cylindrical forms are the most common, as are dense glassy embayments. Increasingly complex shapes and a range of bubble textures are also common. Embayments may crosscut or deflect prominent internal cathodoluminescence banding in the host quartz, indicating that embayments form by both dissolution and growth. We propose potential additional timescales recorded by embayment disequilibrium textures, namely, faceting, bubbles, and the lack thereof. Embayment formation likely occurs tens to hundreds of years before eruption because embayment surfaces are rounded instead of faceted. Bubble textures in embayments are far from those predicted by equilibrium solubility. Homogenous nucleation conditions likely allow preservation of pressures much greater than magmastatic inside embayments. Our textural observations lend insight into embayment occurrence and formation and guide further embayment studies.

The Bishop Tuff, California: New Insights into Magmatic Timescales and Processes from Micro-Analytical Approaches

<p>The Bishop Tuff is the product of one of the largest eruptions on Earth in the last 1 Myr. This thesis studies the Bishop Tuff in order to better understand the nature of the pre-eruptive magma body, with an emphasis on the processes that occurred within it and the timescales over which they operated. In situ geochemical analyses of crystals and glass from samples collected throughout the Bishop Tuff stratigraphic succession yields insights into the nature of zoning and mixing within this supervolcanic system. Timescales for zircon growth (inferred to represent longevity of the magma chamber) are investigated using U-Pb dating of zircons. Zircon textural and trace element data obtained by SIMS (SHRIMP-RG) are presented from 15 stratigraphically controlled Bishop Tuff samples and two older Glass Mountain (GM) lava samples. The resulting eruption age estimate derived from the weighted mean of 166 rim ages of 766.6±3.1 ka (95% confidence) is identical within uncertainty to published values from ID-TIMS and 40Ar/39Ar techniques. An eruption age is also derived for GM dome YA (the youngest GM dome) of 862±23 ka (95% confidence), significantly older than the widely used 790±20 ka K-Ar age. The oldest zircon cores from late-erupted Bishop material (including those with GM-type textures) have a weighted mean of 838.5±8.8 ka (95% confidence), implying that the Bishop Tuff system was only active for ~80 kyr, and had effectively no temporal overlap with the GM system. Bishop zircon textures are divided into four suites whose proportions change systematically through the eruptive sequence. Trace element variations in Bishop zircons are influenced strongly by sector zoning for many elements, and thus restrict the value of trace element variations in discerning compositional stratification within the magma chamber. In later-erupted units, bright-rim overgrowths are common, and are inferred to have crystallized from the same „bright-rim‟ magma as generated the contrasting rims seen in CL or BSE imaging on quartz, feldspar and orthopyroxene. From zircon zonation patterns, this less-evolved, slightly hotter magma invaded deeper parts of the chamber represented in the late-erupted northern units possibly up to ~10 kyr prior to eruption. In order to better quantify the timescales of interaction with the „bright-rim‟ magma, two-feldspar thermometry data are presented on multiple Bishop Tuff samples to constrain temperature variations within the pre-eruptive magma body and yield values for diffusion modelling. Two-feldspar thermometry agrees well with published Fe–Ti-oxide thermometry and reveals a ~80 °C uniform thermal gradient between the upper and lower regions of the magma chamber. Using this thermometry, diffusion of Ti in quartz, Ba in sanidine, Sr in sanidine and Fe-Mg interdiffusion in orthopyroxene are modelled to estimate timescales for the formation of overgrowth rims on crystals. Ti in quartz and Fe-Mg in orthopyroxene diffusion both yield timescales of <150 years for the formation of overgrowth rims, although differing by about an order of magnitude in their timing. However, Ba and Sr diffusion modelling in sanidine yields disparate timescales 1-2 orders of magnitude longer than for Ti in quartz. The main cause for this discrepancy is inferred to be an incorrect assumption for the initial profile shape for Ba and Sr diffusion modelling (i.e. the profile is influenced by growth zoning). Using the comparison with Sr, constraints are placed on the initial width of the core-rim interface and the initial conditions can be refined, bringing Ba and Sr diffusion timescales into mutual alignment and closer to the values from Ti in quartz. This modelling shows that piecemeal rejuvenation of lower Bishop Tuff magma chamber occurred over a period of ~500 years leading up to eruption. In situ major and trace element analyses of sanidine, plagioclase, biotite, orthopyroxene, clinopyroxene, zircon and matrix glass from the Bishop Tuff and two GM lavas are presented to investigate the pre-eruptive stratification of the Bishop magma chamber and its chemical relationship to the GM system. Analyses of samples from the entire Bishop stratigraphy confirm that the magma chamber was thermally and compositionally zoned prior to growth of crystals and the intrusion of the „bright-rim‟ forming magma. Study of rare mixed swirly and dacitic pumice samples shows enrichments in Ba, Sr and Ti (the elements responsible for bright-rim overgrowths in phenocryst phases) and identifies these pumices as possible representatives of the „bright-rim‟ magma. This integrated study of phenocrysts and glass from the Bishop Tuff leads to development of a revised magma chamber model, in which there is a unitary chamber with a stepped or sloping roof. The chamber has an upper, volumetrically dominant (~2/3) part showing no evidence for convection and with unzoned crystals, and a lower part which had experienced mixing of crystals and interaction with the „bright-rim‟ magma. Intrusion of the „bright-rim‟ magma introduced orthopyroxene and dominantly bright zircon crystals, and caused overgrowth of bright rims enriched in Ti, Sr and Ba on sanidine and quartz phenocrysts. Chemical compositions of GM and Bishop Tuff materials show a shared consanguinity, implying common modes of magma generation, yet the generation of GM and Bishop eruptible magma bodies were physically and temporally separate events.</p>

The Bishop Tuff, California: New Insights into Magmatic Timescales and Processes from Micro-Analytical Approaches

<p>The Bishop Tuff is the product of one of the largest eruptions on Earth in the last 1 Myr. This thesis studies the Bishop Tuff in order to better understand the nature of the pre-eruptive magma body, with an emphasis on the processes that occurred within it and the timescales over which they operated. In situ geochemical analyses of crystals and glass from samples collected throughout the Bishop Tuff stratigraphic succession yields insights into the nature of zoning and mixing within this supervolcanic system. Timescales for zircon growth (inferred to represent longevity of the magma chamber) are investigated using U-Pb dating of zircons. Zircon textural and trace element data obtained by SIMS (SHRIMP-RG) are presented from 15 stratigraphically controlled Bishop Tuff samples and two older Glass Mountain (GM) lava samples. The resulting eruption age estimate derived from the weighted mean of 166 rim ages of 766.6±3.1 ka (95% confidence) is identical within uncertainty to published values from ID-TIMS and 40Ar/39Ar techniques. An eruption age is also derived for GM dome YA (the youngest GM dome) of 862±23 ka (95% confidence), significantly older than the widely used 790±20 ka K-Ar age. The oldest zircon cores from late-erupted Bishop material (including those with GM-type textures) have a weighted mean of 838.5±8.8 ka (95% confidence), implying that the Bishop Tuff system was only active for ~80 kyr, and had effectively no temporal overlap with the GM system. Bishop zircon textures are divided into four suites whose proportions change systematically through the eruptive sequence. Trace element variations in Bishop zircons are influenced strongly by sector zoning for many elements, and thus restrict the value of trace element variations in discerning compositional stratification within the magma chamber. In later-erupted units, bright-rim overgrowths are common, and are inferred to have crystallized from the same „bright-rim‟ magma as generated the contrasting rims seen in CL or BSE imaging on quartz, feldspar and orthopyroxene. From zircon zonation patterns, this less-evolved, slightly hotter magma invaded deeper parts of the chamber represented in the late-erupted northern units possibly up to ~10 kyr prior to eruption. In order to better quantify the timescales of interaction with the „bright-rim‟ magma, two-feldspar thermometry data are presented on multiple Bishop Tuff samples to constrain temperature variations within the pre-eruptive magma body and yield values for diffusion modelling. Two-feldspar thermometry agrees well with published Fe–Ti-oxide thermometry and reveals a ~80 °C uniform thermal gradient between the upper and lower regions of the magma chamber. Using this thermometry, diffusion of Ti in quartz, Ba in sanidine, Sr in sanidine and Fe-Mg interdiffusion in orthopyroxene are modelled to estimate timescales for the formation of overgrowth rims on crystals. Ti in quartz and Fe-Mg in orthopyroxene diffusion both yield timescales of <150 years for the formation of overgrowth rims, although differing by about an order of magnitude in their timing. However, Ba and Sr diffusion modelling in sanidine yields disparate timescales 1-2 orders of magnitude longer than for Ti in quartz. The main cause for this discrepancy is inferred to be an incorrect assumption for the initial profile shape for Ba and Sr diffusion modelling (i.e. the profile is influenced by growth zoning). Using the comparison with Sr, constraints are placed on the initial width of the core-rim interface and the initial conditions can be refined, bringing Ba and Sr diffusion timescales into mutual alignment and closer to the values from Ti in quartz. This modelling shows that piecemeal rejuvenation of lower Bishop Tuff magma chamber occurred over a period of ~500 years leading up to eruption. In situ major and trace element analyses of sanidine, plagioclase, biotite, orthopyroxene, clinopyroxene, zircon and matrix glass from the Bishop Tuff and two GM lavas are presented to investigate the pre-eruptive stratification of the Bishop magma chamber and its chemical relationship to the GM system. Analyses of samples from the entire Bishop stratigraphy confirm that the magma chamber was thermally and compositionally zoned prior to growth of crystals and the intrusion of the „bright-rim‟ forming magma. Study of rare mixed swirly and dacitic pumice samples shows enrichments in Ba, Sr and Ti (the elements responsible for bright-rim overgrowths in phenocryst phases) and identifies these pumices as possible representatives of the „bright-rim‟ magma. This integrated study of phenocrysts and glass from the Bishop Tuff leads to development of a revised magma chamber model, in which there is a unitary chamber with a stepped or sloping roof. The chamber has an upper, volumetrically dominant (~2/3) part showing no evidence for convection and with unzoned crystals, and a lower part which had experienced mixing of crystals and interaction with the „bright-rim‟ magma. Intrusion of the „bright-rim‟ magma introduced orthopyroxene and dominantly bright zircon crystals, and caused overgrowth of bright rims enriched in Ti, Sr and Ba on sanidine and quartz phenocrysts. Chemical compositions of GM and Bishop Tuff materials show a shared consanguinity, implying common modes of magma generation, yet the generation of GM and Bishop eruptible magma bodies were physically and temporally separate events.</p>

Fumarolic Pathways Were Structurally Controlled by a Strike-Slip Fault System Beneath the Bishop Tuff, Bishop, California

The Volcanic Tableland, a plateau at the northern end of Owens Valley, CA, is capped by the rhyolitic Bishop Tuff. It hosts many tectonic and volcanic landforms, including hundreds of fault scarps, large joint sets, and inactive fumarolic mounds and ridges. The 1986 Chalfant Valley earthquake sequence shed light on a blind strike-slip fault system beneath the Bishop Tuff. The spatial relationships of the volcanic and tectonic structures have previously been well documented, however, the mechanisms of formation of structures and their enhancement as fumarolic pathways remain largely unexplored. We collected fault kinematic indicators, joint orientations, and documented fumarolic alterations of microcrystalline quartz in the Bishop Tuff and combined those field observations with fault response modeling to assess whether strike-slip activity played a key role in the development of fumarolic pathways. We found field evidence of dip-slip and strike-slip faulting that are consistent with the overall transtensional regional tectonics. Our modeling indicates that a blind strike-slip fault system would dilate joints in the overlying Bishop Tuff with preferred orientations that match observed orientations of joints along which fumarolic activity occurred. Our results imply that the localization of fumaroles was tectonically controlled and that fault activity in the valley floor likely initiated prior to tuff emplacement.

New Ti-in-quartz diffusivities reconcile natural Ti zoning with time scales and temperatures of upper crustal magma reservoirs

Titanium-in-quartz thermometry and diffusion chronometry are routinely applied to felsic magmatic systems. These techniques can be used to determine for how long, and at what temperatures, shallow crustal magmatic systems remain partially molten, both of which are fundamental for assessing volcanic hazards. We have conducted new Ti-in-quartz diffusion experiments at 1 bar, in air, between 900 and 1490 °C, and analyzed the products by secondary ion mass spectrometry (SIMS) depth profiling. The results show that Ti diffusivity is two to three orders of magnitude lower than previously determined {log10D = –8.3 ± 0.4 m2 s–1 – [311 ± 12 kJ mol–1/(2.303RT)]}, where R is the universal gas constant (kJK–1 mol–1) and T is the temperature in Kelvin. Application of these new diffusivities brings time scales determined by Ti-in-quartz diffusion chronometry, using quartz primarily from ignimbrites, into agreement with those determined from zircon U-Pb ages from the Bishop Tuff system (California, USA). This indicates that quartz crystallized early and recorded all, or much of, the thermal history of this magmatic system. These new data also show that sharp Ti zoning profiles can be maintained in quartz within slowly cooled rocks without necessitating that the quartz crystallization temperature is significantly lower than the experimentally determined H2O-saturated granite solidus, or that such samples underwent ultrafast cooling, as has recently been proposed for the granitoids from the Tuolumne Intrusive Suite (California, USA). Finally, our data also indicate that, at least regarding the Bishop Tuff, temperatures must have remained at near-solidus conditions for the entire pre-eruptive evolution of the system, thus relaxing interpretations of “cold storage” for this magmatic system.

A low-aluminum clinopyroxene-liquid geothermometer for high-silica magmatic systems

Several geothermobarometric tools have focused on clinopyroxene due to its prevalence in igneous rocks, however clinopyroxene produced in high-silica igneous systems is high in iron and low in aluminum, causing existing geothermometers that depend on aluminum exchange to fail or yield overestimated temperatures. Here we present a new clinopyroxene-liquid geothermometer recommended for use in natural igneous systems with bulk SiO2 ≥ 70 wt%, which contain clinopyroxene with Mg# ≤ 65 and Al2O3 ≤ 7 wt%. (1) T ( ∘ C ) = 300 [ − 1.89 − 0.601 ( X CaTs Cpx ) − 0.186 ( X DiHd 2003 Cpx ) + 4.71 ( X SiO 2 liq ) + 77.6 ( X TiO 2 liq ) + 10.9 ( X FeO liq ) + 33.6 ( X MgO liq ) + 15.5 ( X CaO liq ) + 15.6 ( X KO 0.5 liq ) ] The new geothermometer lowers calculated temperatures by ~85 °C on average relative to Putirka (2008, Eq. 33) and reduces the uncertainty by a factor of two (standard error of estimate ±20 °C). When applied to natural systems, we find this new clinopyroxene-liquid geothermometer reconciles many inconsistencies between experimental phase equilibria and preexisting geothermometry results for silicic volcanism, including those from the Bishop Tuff and Yellowstone caldera-forming and post-caldera rhyolites. We also demonstrate that clinopyroxene is not restricted to near-liquidus temperatures in rhyolitic systems; clinopyroxene can be stable over a broad temperature range, often down to the solidus. An Excel spreadsheet and Python notebook for calculating temperature with this new geothermometer may be downloaded from GitHub at http://bit.ly/cpxrhyotherm.