skip to main content

Title: Experimental Constraints on Dacite Magma Storage beneath Volcán Quizapu, Chile
Abstract Volcán Quizapu, Chile, is an under-monitored volcano that was the site of two historical eruptions: an effusive eruption in 1846–1847 and a Plinian eruption in 1932, both of which discharged ∼5 km3 (dense rock equivalent) of lava and/or tephra. The majority of material erupted in both cases is trachydacite, nearly identical for each event. We present H2O-saturated, phase equilibrium experiments on this end-member dacite magma, using a pumice sample from the 1932 eruption as the main starting material. At an oxygen fugacity (fO2) of ∼NNO + 0·2 (where NNO is the nickel–nickel oxide buffer), the phase assemblage of An25–30 plagioclase + amphibole + orthopyroxene, without biotite, is stable at 865 ± 10 °C and 110 ± 20 MPa H2O pressure (PH2O), corresponding to ∼4 km depth. At these conditions, experiments also reproduce the quenched glass composition of the starting pumice. At slightly higher PH2O and below 860 °C, biotite joins the equilibrium assemblage. Because biotite is not part of the observed Quizapu phase assemblage, its presence places an upper limit on PH2O. At the determined storage PH2O of ∼110 MPa, H2O undersaturation of the magma with XH2Ofluid = 0·87 would align Ptotal to mineral-based geobarometry estimates of ∼130 MPa. However, XH2Ofluid < 1 is not required to reproduce the Quizapu dacite phase more » assemblage and compositions. A second suite of experiments at lower fO2 shows that the stability fields of the hydrous silicates (amphibole and biotite) are significantly restricted at NNO – 2 relative to NNO + 0·2. Additional observations of Quizapu lava and pumice samples support the existing hypothesis that rapid pre-eruptive heating drove the effusive 1846–1847 eruption, with important refinements. We demonstrate that microlites in the end-member dacite lavas are consistent with in situ crystallization (during ascent), rather than transfer from an andesite. In one end-member dacite lava, newly identified reverse zoning in orthopyroxene and incipient destabilization of amphibole are consistent with small degrees of heating. Our work articulates a clear direction for future Quizapu studies, which are warranted given the active nature of the Cerro Azul–Descabezado Grande volcanic axis. « less
Authors:
; ; ;
Award ID(s):
1347880
Publication Date:
NSF-PAR ID:
10273230
Journal Name:
Journal of Petrology
Volume:
62
Issue:
5
ISSN:
0022-3530
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract The Bishop Tuff (BT), erupted from the Long Valley caldera in California, displays two types of geochemical gradients with temperature: one is related to magma mixing, whereas the other is found in the high-SiO2 rhyolite portion of the Bishop Tuff and is characterized by twofold or lower concentration variations in minor and trace elements that are strongly correlated with temperature. It is proposed that the latter zonation, which preceded phenocryst growth, developed as a result of mineral–melt partitioning between interstitial melt and surrounding crystals in a parental mush, from which variable melt fractions were segregated. To test this hypothesis, trends of increasing vs decreasing element concentrations with temperature (as a proxy for melt fraction), obtained from published data on single-clast pumice samples from the high-SiO2 rhyolite portion of the Bishop Tuff, were used to infer their relative degrees of incompatibility vs compatibility between crystals and melt in the parental mush. Relative compatibility values (RCVi) for all elements i, defined as the concentration slope with temperature divided by average concentration, are shown to be linearly correlated with their respective bulk partition coefficients (bulk Di). Mineral–melt partition coefficients from the literature were used to constrain the average stoichiometry of the crystallization/meltingmore »reaction in the parental mush: 32 % quartz + 34 % plagioclase + 31 % K-feldspar + 1·60 % biotite + 0·42 % titanomagnetite + 0·34 % ilmenite + 0·093 % allanite + 0·024 % zircon + 0·025 % apatite = 100 % liquid. The proportions of tectosilicates in the reaction (i.e. location of eutectic) are consistent with depths of melt segregation of ~400–550 MPa and an activity of H2O of ~0·4–0·6. Temperatures of <770–780 °C are constrained by allanite in the reaction. Evidence that a fluid phase was present in the parental mush is seen in the decreasing versus increasing H2O and CO2 contents with temperature in the segregated interstitial melt that formed the high-SiO2 rhyolite portion of the Bishop Tuff. The presence of an excess fluid phase, which strongly partitions CO2 relative to the melt, is required to explain the compatible behavior of CO2, whereas the fluid abundance must have been low to explain the incompatible behavior of H2O. Calculated degassing paths for interstitial melts, which segregated from the parental mush and ascended to shallower depths to grow phenocrysts, match published volatile analyses in quartz-hosted melt inclusions and constrain fluid abundances in the mush to be ≤1 wt%. The source of volatiles in the parental mush, irrespective of whether it formed by crystallization or partial melting, must have been primarily from associated basalts, as granitoid crust is too volatile poor. Approximately twice as much basalt as rhyolite is needed to provide the requisite volatiles. The determination of bulk Di for several elements gives the bulk composition of the parental leucogranitic mush and shows that it is distinct from Mesozoic Sierran arc granitoids, as expected. Collectively, the results from this study provide new constraints for models of the complex, multi-stage processes throughout the Plio-Quaternary, involving both mantle-derived basalt and pre-existing crust, that led to the origin of the parental body to the Bishop Tuff.« less
  2. Volcanic eruptions of rhyolitic magma often show shifts from powerful (Vulcanian to Plinian) explosive episodes to a more gentle effusion of viscous lava forming obsidian flows. Another prevailing characteris-tic of these eruptions is the presence of pyroclastic obsidians intermingled with the explosive tephra. This dense, juvenile product is similar to the tephra and obsidian flow in composition, but is generally less degassed than its flow counterpart. The formation mechanism(s) of pyroclastic obsidians and the information they can provide concerning the extent to which magma degassing modulates the eruptive style of rhyolitic eruptions are currently subject to active research. Porous tephra and pyroclastic and flow obsidians from the 1060CE Glass Mountain rhyolitic eruption at Medicine Lake Volcano (California) were analyzed for their porosity, φ, water content, H2O, and hydrogen isotopic composition, δD. H2O in porous pyroclasts is correlated negatively with δD and positively with φ, indicating that the samples were affected by post-eruptive rehydration. Numerical modeling suggests that this rehydration occurred at an average rate of 10−23.5±0.5m2s−1during the ∼960 years since the eruption, causing some pyroclasts to gain up to 1 wt%of meteoric water. Pyroclastic and flow obsidians were not affected by rehydration due to their very low porosity. Comparison betweenmore »modeled δD-H2O relationships in degassing magma and values measured in the Glass Mountain samples supports the idea that rhyolitic magma degasses in closed-system until its porosity reaches a value of about 65±5%, beyond which degassing occurs in open-system until quench. During the explosive phase, rapidly ascending magma fragments soon after it becomes permeable, creating porous lapilli and ash that continue to degas in open-system within an expanding gas phase. As suggested by recent studies, some ash may aggregate and sinter on the conduit sides at different depths above the fragmentation level, partly equilibrating with the continuously fluxing heavier magmatic vapor, explaining the wide range of H2O contents and high variability in δD measured in the pyroclastic obsidians. Using only H2O and δD, it is impossible to rule out the possibility that pyroclastic obsidians may also form by permeable foam collapse, either syn-explosively near the conduit sides below the fragmentation level or during more effusive periods interspersed in the explosive phase. During the final effusive phase of the eruption, slowly ascending magma degasses in open-system until it reaches the surface, creating flows with low H2O and δD. This study shows that H2O measured in highly porous pyroclasts of a few hundred years or more cannot be used to infer syn-eruptive magma degassing pathways, unless careful assessment of post-eruptive rehydration is first carried out. If their mechanism of formation can be better understood, detailed analysis of the variations in texture and volatile content of pyroclastic obsidians throughout the explosive phase may help decipher the reasons why rhyolitic eruptions commonly shift from explosive to effusive phases.« less
  3. Hydrogen is a rapidly diffusing monovalent cation in nominally anhydrous minerals (NAMs, such as olivine, orthopyroxene, and clinopyroxene), which is potentially re-equilibrated during silicate melt-rock and aqueous fluid-rock interactions in massif and abyssal peridotites. We apply a 3D numerical diffusion modeling technique to provide first-order timescales of complete hydrogen re-equilibration in olivine, clinopyroxene, and orthopyroxene over the temperature range 600-1200°C. Model crystals are 1-3 mm along the c-axis and utilize H+ diffusion coefficients appropriate for Fe-bearing systems. Two sets of models were run with different boundary compositions: 1) “low-H models” are constrained by mineral-melt equilibrium partitioning with a basaltic melt that has 0.75 wt% H2O and 2) “high-H models,” which utilize the upper end of the estimated range of mantle water solubility for each phase. Both sets of models yield re-equilibration timescales that are identical and are fast for all phases at a given temperature. These timescales have strong log-linear trends as a function of temperature (R2 from 0.97 to 0.99) that can be used to calculate expected re-equilibration time at a given temperature and grain size. At the high end of the model temperatures (1000-1200°C), H+ completely re-equilibrates in olivine, orthopyroxene, and clinopyroxene within minutes to hours, consistent withmore »previous studies. These short timescales indicate that xenolith NAM mantle water contents are likely to be overprinted prior to eruption. The models also resolve the decoupled water-trace element relationship in Southwest Indian Ridge peridotites, in which peridotite REE abundances are reproduced by partial melting models whereas the relatively high NAM H2O contents require later re-equilibration with melt. At temperatures of 600-800°C, which correspond to conditions of hydrothermal alteration of pyroxene to amphibole and talc, H+ re-equilibration typically occurs over a range of timescales spanning days to years. These durations are well within existing estimates for the duration of fluid flow in oceanic hydrothermal systems, suggesting that peridotite NAM water contents are susceptible to diffusive overprinting during higher temperature hydrothermal alteration. Thus, diffusion during aqueous fluid-rock interactions may also explain NAM H2O contents that are too high to reflect residues of melting. These relatively short timescales at low temperatures suggest that the origin of water contents measured in peridotite NAMs requires additional constraints on sample petrogenesis, including petrographic and trace element analyses. Our 3D model results also hint that H+ may diffuse appreciably during peridotite serpentinization, but diffusion coefficients at low temperature are unconstrained and additional experimental investigations are needed.« less
  4. Abstract We present microbeam major- and trace-element data from 14 monzodiorites collected from the Malaspina Pluton (Fiordland, New Zealand) with the goal of evaluating processes involved in the production of andesites in lower arc crust. We focus on relict igneous assemblages consisting of plagioclase and amphibole with lesser amounts of clinopyroxene, orthopyroxene, biotite and quartz. These relict igneous assemblages are heterogeneously preserved in the lower crust within sheeted intrusions that display hypersolidus fabrics defined by alignment of unstrained plagioclase and amphibole. Trace-element data from relict igneous amphiboles in these rocks reveal two distinct groups: one relatively enriched in high field strength element concentrations and one relatively depleted. The enriched amphibole group has Zr values in the range of ∼25–110 ppm, Nb values of ∼5–32 ppm, and Th values up to 2·4 ppm. The depleted group, in contrast, shows Zr values <35 ppm and Nb values <0·25 ppm, and Th is generally below the level of detection. Amphibole crystallization temperatures calculated from major elements range from ∼960 to 830 °C for all samples in the pluton; however, we do not observe significant differences in the range of crystallization temperatures between enriched (∼960–840 °C) and depleted groups (∼940–830 °C). Bulk-rock Sr and Nd isotopes are also remarkably homogeneous and showmore »no apparent difference between enriched (εNdi = 0·1 to –0·1; 87Sr/86Sri = 0·70420–0·70413) and depleted groups (εNdi = 0·3 to –0·4; 87Sr/86Sri = 0·70424–0·70411). Calculated amphibole-equilibrium melt compositions using chemometric equations indicate that melts were highly fractionated (molar Mg# <50), andesitic to dacitic in composition, and were much more evolved than bulk lower continental crust or primitive basalts and andesites predicted to have formed from hydrous melting of mantle-wedge peridotite beneath an arc. We suggest that melts originated from a common, isotopically homogeneous source beneath the Malaspina Pluton, and differences between enriched and depleted trace-element groups reflect varying contributions from subducted sediment-derived melt and sediment-derived fluid, respectively. Our data demonstrate that andesites and dacites were the dominant melts that intruded the lower crust, and their compositions mirror middle and upper bulk-continental crust estimates. Continental crust-like geochemical signatures were acquired in the source region from interaction between hydrous mantle-wedge melts and recycled subducted sediment rather than assimilation and/or remelting of pre-existing lower continental crust.« less
  5. Abstract Drilling related to development of the platinum-group element deposit of the J-M Reef of the Stillwater Complex returned samples of a rare chromitite seam between anorthosite and norite in a discordant anorthositic body. Plagioclase core An concentrations are marginally higher and modestly reversely zoned on the norite side (average Ancore = 83·8; average Ancore – Anrim = –1·1) as compared with the anorthosite side (average Ancore 82·5; average Ancore – Anrim = +1·0). The anorthosites are also characterized by a slightly smaller average plagioclase grain size than plagioclase in the norite (1·41 mm and 1·54 mm, respectively). The chromite can contain single and polyphase inclusions of orthopyroxene, plagioclase, amphibole, biotite and Cl-rich apatite. These and other compositional and textural features, as well as inference from other discordant anorthositic bodies in the Banded series, are all consistent with a chromatographic model of chromite precipitation at a reaction front as a norite protolith reacts with a Cl-rich aqueous fluid saturated in plagioclase alone. Chromitite seam formation is modeled using an infiltration metasomatic model, in which a fluid becomes progressively undersaturated in pyroxene as it rises into the hotter part of the crystal pile. As this pyroxene-undersaturated fluid moves through a noritic protolith, itmore »dissolves the Cr-bearing orthopyroxene to produce an anorthosite. Chromite precipitates at the reaction front between the anorthosite and the norite owing to liberation of Mg and Cr from pyroxene. Continuous redissolution and reprecipitation of chromite occurs as the pyroxene dissolution front moves in the direction of fluid flow, collecting the Cr lost from the anorthosite. Owing to Cr dissolved mainly as a neutral divalent cation complex, CrCl(OH)0, in the solution, but incorporated as a trivalent cation in chromite, the required redox reaction can involve concurrent precipitation of sulfide with chromite. This mechanism differs from some recent models in that the anorthosites are themselves replacement bodies and are not original precipitates from a magma nor formed by loss of mafic material by partial melting. The results show the need for experimental mineral solubility data at T and P conditions appropriate to upper crustal mafic–ultramafic intrusions.« less