More Like this

1. Connecting Magma Chamber Mush Dynamics to Post Eruptive Reinflation at Cordón Caulle

   Patrick Phelps, Helge Gonnermann (January 2021, AGU Fall Meeting, New Orleans, LA & Online Everywhere, 13-17 December, 2021)

   Volcano inflation, for durations of months to years immediately following an eruption, has been observed at a number of volcanoes, including the 2011/12 eruption of Cordón Caulle, Chile. Such reinflation is often explained by replenishment of the magma reservoir from a deeper source. Whether and why that is the case remains uncertain in most instances, but the implications for renewed eruptive potential may be profound. Here, we posit redistribution of melt within a zoned magma reservoir consisting of a crystal rich mush overlain by an eruptible layer of crystal poor rhyolite as an alternate mechanism for reinflation. Such a zoned magma body is consistent with conceptual models for how crystal poor rhyolites form and with the presence of mafic enclaves within the Cordón Caulle rhyolite. The enclaves can be interpreted as pieces of mush entrained into the overlying rhyolite during its withdrawal from the reservoir. We test the hypothesis that melt from the inter-crystalline pores of the mush can redistribute by porous flow into the overlying crystal poor rhyolite, causing inflation after an eruption. We simulate the flow of melt within the zoned reservoir during and after eruption with a numerical model. As crystal poor rhyolite is erupted, magma pressure within the rhyolite layer above the mush decreases. Consequently, interstitial melt flows upward within the mush, toward the reduced pressure at the interface of mush and crystal poor rhyolite. The mush is treated as a poroelastic material, with interstitial melt flow governed by Darcy's law. Thus, the change in pressure caused by withdrawal from the overlying rhyolite diffuses downward into the mush as the interstitial melt flows upward. The change in pore pressure results in an elastic deformation of the mush matrix. Because pore pressure diffusivity is small, melt redistribution can persist for years after eruptive activity ends, leading to slow inflation compared to fast eruptive deflation. We predict a partial recovery of volume lost from the eruption. Reinflation occurs because the expansion of decompressing melt flowing from the mush into the crystal poor rhyolite exceeds compression of the poroelastic mush. For cases where the interstitial melt is moderately compressible due to exsolved volatiles, our model reproduces the deformation observed at Cordón Caulle. 

   more » « less
2. Effects of H2O–CO2 Fluids, Temperature, and Peridotite Fertility on Partial Melting in Mantle Wedges and Generation of Primary Arc Basalts

   https://doi.org/10.1093/petrology/egad047  

   Lara, Michael; Dasgupta, Rajdeep (July 2023, Journal of Petrology)

   Abstract Many lines of evidence from high P–T experiments, thermodynamic models, and natural observations suggest that slab-derived aqueous fluids, which flux mantle wedges contain variable amounts of dissolved carbon. However, constraints on the effects of H2O–CO2 fluids on mantle melting, particularly at mantle wedge P–T conditions, are limited. Here, we present new piston cylinder experiments on fertile and depleted peridotite compositions with 3.5 wt.% H2O and XCO2 [= molar CO2 / (CO2 + H2O)] of 0.04–0.17. Experiments were performed at 2–3 GPa and 1350°C to assess how temperature, peridotite fertility, and XCO2 of slab-derived fluid affects partial melting in mantle wedges. All experiments produce olivine + orthopyroxene +7 to 41 wt.% partial melt. Our new data, along with previous lower temperature data, show that as mantle wedge temperature increases, primary melts become richer in SiO2, FeO*, and MgO and poorer CaO, Al2O3, and alkalis when influenced by H2O–CO2 fluids. At constant P–T and bulk H2O content, the extent of melting in the mantle wedge is largely controlled by peridotite fertility and XCO2 of slab-fluid. High XCO2 depleted compositions generate ~7 wt.% melt, whereas, at identical P–T, low XCO2 fertile compositions generate ~30 to 40 wt.% melt. Additionally, peridotite fertility and XCO2 have significant effects on peridotite partial melt compositions. At a constant P–T–XCO2, fertile peridotites generate melts richer in CaO and Al2O3 and poorer in SiO2, MgO + FeO, and alkalis. Similar to previous experimental studies, at a constant P–T fertility condition, as XCO2 increases, SiO2 and CaO of melts systematically decrease and increase, respectively. Such distinctive effects of oxidized form of dissolved carbon on peridotite partial melt compositions are not observed if the carbon-bearing fluid is reduced, such as CH4-bearing. Considering the large effect of XCO2 on melt SiO2 and CaO concentrations and the relatively oxidized nature of arc magmas, we compare the SiO2/CaO of our experimental melts and melts from previous peridotite + H2O ± CO2 studies to the SiO2/CaO systematics of primitive arc basalts and ultra-calcic, silica-undersaturated arc melt inclusions. From this comparison, we demonstrate that across most P–T–fertility conditions predicted for mantle wedges, partial melts from bulk compositions with XCO2 ≥ 0.11 have lower SiO2/CaO than all primitive arc melts found globally, even when correcting for olivine fractionation, whereas partial melts from bulk compositions with XCO2 = 0.04 overlap the lower end of the SiO2/CaO field defined by natural data. These results suggest that the upper XCO2 limit of slab-fluids influencing primary arc magma formation is 0.04 < XCO2 < 0.11, and this upper limit is likely to apply globally. Lastly, we show that the anomalous SiO2/CaO and CaO/Al2O3 signatures observed in ultra-calcic arc melt inclusions can be reproduced by partial melting of either CO2-bearing hydrous fertile and depleted peridotites with 0 < XCO2 < 0.11 at 2–3 GPa, or from nominally CO2-free hydrous fertile peridotites at P > 3 GPa. 

   more » « less
3. Feasibility of melt segregation from a crystal mush in response to the 2011–2012 eruption at Cordón Caulle, Chile

   https://doi.org/10.1093/gji/ggad259  

   Phelps, Patrick_R; Gonnermann, Helge_M; Winslow, Heather; Ruprecht, Philipp; Pritchard, Matthew_E; Delgado, Francisco; Liao, Yang (June 2023, Geophysical Journal International)

   SUMMARY The 2011–2012 eruption at Cordón Caulle in Chile produced crystal-poor rhyolitic magma with crystal-rich mafic enclaves whose interstitial glass is of identical composition to the host rhyolite. Eruptible rhyolites are thought to be genetically associated with crystal-rich magma mushes, and the enclaves within the Cordón Caulle rhyolite support the existence of a magma mush from which the erupted magma was derived. Moreover, towards the end of the 2011–2012 eruption, subsidence gave way to inflation that has on average been continuous through at least 2020. We hypothesize that magma segregation from a crystal mush could be the source of the observed inflation. Conceptually, magma withdrawal from a crystal-poor rhyolite reservoir caused its depressurization, which could have led to upward flow of interstitial melt within an underlying crystal mush, causing a new batch of magma to segregate and partially recharge the crystal-poor rhyolite body. Because the compressibility of the crystalline matrix of the mush is expected to be lower than that of the interstitial melt, which likely contains some fraction of volatile bubbles, this redistribution of melt would result in a net increase in volume of the system and in the observed inflation. We use numerical modelling of subsurface magma flow and storage to show under which conditions such a scenario is supported by geodetic and petrologic observations. 

   more » « less
4. Experimentally derived F, Cl, and Br fluid/melt partitioning of intermediate to silicic melts in shallow magmatic systems

   https://doi.org/10.2138/am-2022-8109  

   Cassidy, Mike; Iveson, Alexander A; Humphreys, Madeleine CS; Mather, Tamsin A; Helo, Christoph; Castro, Jonathan M; Ruprecht, Philipp; Pyle, David M (October 2022, American Mineralogist)

   Abstract The conditions under which halogens partition in favor of an exsolved fluid relative to the coexisting melt are key for understanding many magmatic processes, including volcanic degassing, evolution of crustal melt bodies, and ore formation. We report new F, Cl, and Br fluid/melt partition coefficients for intermediate to silicic melts, for which F and Br data are particularly lacking; and for varying CO2-H2O contents to assess the effects of changing fluid composition (XH2O) on Br fluid/melt partitioning for the first time. The experiments were conducted at pressures 50–120 MPa, temperatures 800–1100 °C, and volatile compositions [molar XH2O = H2O/(H2O +CO2)] of 0.55 to 1, with redox conditions around the Nickel-Nickel Oxygen buffer (fO2 ≈ NNO). Experiments were not doped with Cl, Br, or F and were conducted on natural crystal-bearing volcanic products at conditions close to their respective pre-eruptive state. The experiments therefore provide realistic constraints on halogen partitioning at naturally occurring, brine-undersaturated conditions. Measurements of Br, Cl, and F were made by Secondary Ion Mass Spectrometry (SIMS) on 13 experimental glass products spanning andesite to rhyolitic compositions, together with their natural starting materials from Kelud volcano, Indonesia, and Quizapu volcano, Chile. Fluid compositions were constrained by mass balance. Average bulk halogen fluid/melt partition coefficients and standard deviations are: DClfluid/melt = 3.4 (±3.7 1 s.d.), DFfluid/melt = 1.7 (±1.7), and DBrfluid/melt = 7.1 (±6.4) for the Kelud starting material (bulk basaltic andesite), and DClfluid/melt = 11.1 (±3.5), DFfluid/melt = 0.8 (±0.8), and DBrfluid/melt = 31.3 (±20.9) for Quizapu starting material (bulk dacite). The large range in average partition coefficients is a product of changing XH2O, pressure and temperature. In agreement with studies on synthetic melts, our data show an exponential increase of halogen Dfluid/melt with increasing ionic radius, with partitioning behavior controlled by melt composition according to the nature of the complexes forming in the melt (e.g., SiF4, NaCl, KBr). The fundamental chemistry of the different halogens (differing ionic size and electronegativities) controls the way in which partitioning responds to changes in melt composition and other variables. Experimental results confirm that more Cl partitions into the fluid at higher bulk Cl contents, higher melt Na, higher fluid XH2O ratios, and lower temperatures. Bromine shows similar behavior, though it seems to be more sensitive to temperature and less sensitive to Na content and XH2O. In contrast, F partitioning into the fluid increases as the melt silica content decreases (from 72 to 56 wt% SiO2), which we attribute to the lower abundance of Si available to form F complexes in the melt. These new data provide more insights into the conditions and processes that control halogen degassing from magmas and may help to inform the collection and interpretation of melt inclusions and volcano gas data. 

   more » « less
5. Constraints on deep, CO2-rich degassing at arc volcanoes from solubility experiments on hydrous basaltic andesite of Pavlof Volcano, Alaska Peninsula, at 300 to 1200 MPa

   https://doi.org/10.2138/am-2021-7531  

   Mangan, Margaret T.; Sisson, Thomas W.; Hankins, W. Ben; Shimizu, Nobumichi; Vennemann, Torsten (May 2021, American Mineralogist)

   Abstract The solubility of CO2 in hydrous basaltic andesite was examined in fO2-controlled experiments at a temperature of 1125 °C and pressures between 310–1200 MPa. Concentrations of dissolved H2O and CO2 in experimental glasses were determined by ion microprobe calibrated on a subset of run glasses analyzed by high-temperature vacuum manometry. Assuming that the solubility of H2O in mafic melt is relatively well known, estimates of XH2Ofluid and PH2Ofluid in the saturating fluid were modeled, and by difference, values for XCO2fluid and PCO2fluid were obtained (XCO2 ~0.5–0.9); fCO2 could be then calculated from the fluid composition, temperature, and pressure. Dissolved H2O over a range of 2.3–5.5 wt% had no unequivocal influence on the dissolution of CO2 at the pressures and fluid compositions examined. For these H2O concentrations, dissolved CO2 increases with fCO2 following an empirical power-law relation: dissolved CO2 (ppmw) = 14.9−3.5+4.5[fCO2 (MPa)]0.7±0.03. The highest-pressure results plot farthest from this equation but are within its 1 standard-error uncertainty envelope. We compare our experimental data with three recent CO2-H2O solubility models: Papale et al. (2006); Iacono-Marziano et al. (2012); and Ghiorso and Gualda (2015). The Papale et al. (2006) and Iacono-Marizano et al. (2012) models give similar results, both over-predicting the solubility of CO2 in a melt of the Pavlof basaltic andesite composition across the fCO2 range, whereas the Ghiorso and Gualda (2015) model under-predicts CO2 solubility. All three solubility models would indicate a strong enhancement of CO2 solubility with increasing dissolved H2O not apparent in our results. We also examine our results in the context of previous high-pressure CO2 solubility experiments on basaltic melts. Dissolved CO2 correlates positively with mole fraction (Na+K+Ca)/Al across a compositional spectrum of trachybasalt-alkali basalt-tholeiite-icelandite-basaltic andesite. Shortcomings of current solubility models for a widespread arc magma type indicate that our understanding of degassing in the deep crust and uppermost mantle remains semi-quantitative. Experimental studies systematically varying concentrations of melt components (Mg, Ca, Na, K, Al, Si) may be necessary to identify solubility reactions, quantify their equilibrium constants, and thereby build an accurate and generally applicable solubility model. 

   more » « less