skip to main content


Title: Melt Migration in Crystal Mushes by Viscous Fingering: Insights From High‐Temperature, High‐Pressure Experiments
Abstract

We conducted experiments to study melt migration in crystal‐rich mushes, with application to magma ascent within transcrustal magma reservoirs. Mushes with crystal volume fractions of 0.59–0.83 were prepared by hot‐pressing crushed borosilicate glass mixed with different proportions of quartz sand particles. Each experimental sample comprises stacked disks of mush and soda‐lime glass, a proxy for crystal‐free magma. Samples were subjected to confining pressures of 100–300 MPa and a temperature of 900°C (above the glass transition temperatures of the borosilicate and soda‐lime glasses) for up to 6 h. The bottom and circumference of the mush and soda lime disks experience the confining pressure, but the top of the mush disks is at room pressure, resulting in a pore‐pressure gradient across the mush layer. Following cooling and decompression, we determined the area fraction and morphology of soda‐lime melt that migrated into the mush layer during experiments. Melt fraction is more strongly correlated to crystal fraction than pore‐pressure gradient, increasing with crystal fraction before sharply decreasing as crystal fractions exceed 0.8. This change at 0.8 coincides with the transition from crystals in the mush moving during soda‐lime migration to crystals forming a continuous rigid network. In our experiments, melt migration occurred by viscous fingering, but near the mobile‐to‐rigid transition, melt migration is enhanced by additional capillary action. Our results indicate that magma migration may peak when rigid mushes “unlock” to become mobile. This transition may mark an increase in magma migration, a potential precursor to volcanic unrest and eruption.

 
more » « less
Award ID(s):
2052599
NSF-PAR ID:
10372042
Author(s) / Creator(s):
 ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Solid Earth
Volume:
127
Issue:
8
ISSN:
2169-9313
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. BSE mosaics of mushes and experimental products. Abstract: "We conducted experiments to study melt migration in crystal-rich mushes, with application to magma ascent within transcrustal magma reservoirs. Mushes with crystal volume fractions of 0.59 to 0.83 were prepared by hot-pressing crushed borosilicate glass mixed with different amounts quartz sand particles. Each experimental sample comprises stacked disks of mush and soda-lime glass, a proxy for crystal-free magma. Samples were subjected to confining pressures of 100 to 300 MPa and a temperature of 900°C (above the glass transition temperatures of the borosilicate and soda-lime glasses) for up to 6 h. The bottom and circumference of the mush and soda lime disks experience the confining pressure, but the top of the mush disks are at room pressure, resulting in a pore-pressure gradient across the mush layer. Following cooling and decompression, we determined the area fraction and morphology of soda-lime melt that migrated into the mush layer during experiments. Melt fraction is more strongly correlated to crystal fraction than pore-pressure gradient, increasing with crystal fraction before sharply decreasing as crystal fractions exceed 0.8. This change at 0.8 coincides with the transition from crystals in the mush moving during soda-lime migration to crystals forming a continuous rigid network. In our experiments, melt migration occurred by viscous fingering, but near the mobile-to-rigid transition, melt migration is enhanced by additional capillary action. Our results indicate that magma migration may peak when rigid mushes “unlock” to become mobile. This transition may mark an increase in magma migration, a potential precursor to volcanic unrest and eruption." Imaging: "Transverse sections cut from the top and/or bottom of the vacuum hot-pressed mushes were polished, carbon-coated, and imaged in BSE mode using the JEOL JXA-8530FPlus Electron Probe Microanalyzer (EPMA) at UMN (15 kV, 10 nA). About ten 50x magnification images were taken per sample and then compiled into BSE mosaics using Affinity Designer. The different compositions of the borosilicate glass and crystalline materials are distinguishable by greyscale in BSE images. [...] Following each experiment, sample assemblies were cut longitudinally along the cylindrical axis to produce sections for microstructural analysis. Scored samples (pHi-19s, Int-20s, Lo-21s) were cut again to produce sections tangential to the sample cylinder. Cut sections were vacuum impregnated with EpoFix resin and hand-polished on diamond lapping film from 30 to 0.5 μm grit. Polished and carbon-coated samples were imaged in BSE mode in the EPMA at UMN (15 kV, 10 nA). The different compositions of the soda-lime glass, borosilicate glass, and crystalline materials are distinguishable by greyscale in BSE images. Twenty to forty 50x magnification images were taken per sample and then compiled into sample-scale BSE mosaics using Affinity Designer." 
    more » « less
  2. Abstract We investigate the shallow plumbing system of the Deccan Traps Large Igneous Province using rock and mineral data from Giant Plagioclase Basalt (GPB) lava flows from around the entire province, but with a focus on the Saurashtra Peninsula, the Malwa Plateau, and the base and top of the Western Ghats (WG) lava pile. GPB lavas in the WG typically occur at the transition between chemically distinct basalt formations. Most GPB samples are evolved basalts, with high Fe and Ti contents, and show major and trace elements and Sr-Nd-Pb isotopic compositions generally similar to those of previously studied Deccan basalts. Major element modeling suggests that high-Fe, evolved melts typical of GPB basalts may derive from less evolved Deccan basalts by low-pressure fractional crystallization in a generally dry magmatic plumbing system. The basalts are strongly porphyritic, with 6–25% of mm- to cm-sized plagioclase megacrysts, frequently occurring as crystal clots, plus relatively rare olivine and clinopyroxene. The plagioclase crystals are mostly labradoritic, but some show bytownitic cores (general range of anorthite mol%: 78–55). A common feature is a strong Fe enrichment at the plagioclase rims, indicating interaction with an Fe-rich melt similar to that represented by the matrix compositions (FeOt up to 16–17 wt%). Plagioclase minor and trace elements and Sr isotopic compositions analyzed by laser ablation inductively coupled plasma mass spectrometry show evidence of a hybrid and magma mixing origin. In particular, several plagioclase crystals show variable 87Sr/86Sri, which only partially overlaps with the 87Sr/86Sri of the surrounding matrix. Diffusion modeling suggests residence times of decades to centuries for most plagioclase megacrysts. Notably, some plagioclase crystal clots show textural evidence of deformation as recorded by electron back-scatter diffraction analyses and chemical maps, which suggest that the plagioclase megacrysts were deformed in a crystal-rich environment in the presence of melt. We interpret the plagioclase megacrysts as remnants of a crystal mush originally formed in the shallow plumbing system of the Deccan basalts. In this environment, plagioclase acquired a zoned composition due to the arrival of chemically distinct basaltic magmas. Prior to eruption, a rapidly rising but dense Fe-rich magma was capable of disrupting the shallow level crystal mush, remobilizing part of it and carrying a cargo of buoyant plagioclase megacrysts. Our findings suggest that basaltic magmas from the Deccan Traps, and possibly from LIPs in general, are produced within complex transcrustal magmatic plumbing systems with widespread crystal mushes developed in the shallow crust. 
    more » « less
  3. 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. Seismic imaging methods have provided detailed three-dimensional constraints on the physical properties of magmatic systems leading to invaluable insight into the storage, differentiation and dynamics of magma. These constraints have been crucial to the development of our modern understanding of magmatic systems. However, there are still outstanding knowledge gaps resulting from the challenges inherent in seismic imaging of volcanoes. These challenges stem from the complex physics of wave propagation across highly heterogeneous low-velocity anomalies associated with magma reservoirs. Ray-based seismic imaging methods such as travel-time and surface-wave tomography lead to under-recovery of such velocity anomalies and to under-estimation of melt fractions. This review aims to help the volcanologist to fully utilize the insights gained from seismic imaging and account for the resolution limits. We summarize the advantages and limitations of the most common imaging methods and propose best practices for their implementation and the quantitative interpretation of low-velocity anomalies. We constructed and analysed a database of 277 seismic imaging studies at 78 arc, hotspot and continental rift volcanoes. Each study is accompanied by information about the seismic source, part of the wavefield used, imaging method, any detected low-velocity zones, and estimated melt fraction. Thirty nine studies attempted to estimate melt fractions at 22 different volcanoes. Only five studies have found evidence of melt storage at melt fractions above the critical porosity that separates crystal mush from mobile magma. The median reported melt fraction is 13% suggesting that magma storage is dominated by low-melt fraction crystal mush. However, due to the limits of seismic resolution, the seismological evidence does not rule out the presence of small (<10 km 3 ) and medium-sized (<100 km 3 ) high-melt fraction magma chambers at many of the studied volcanoes. The combination of multiple tomographic imaging methods and the wider adoption of methods that use more of the seismic wavefield than the first arriving travel-times, promise to overcome some of the limitations of seismic tomography and provide more reliable constraints on melt fractions. Wider adoption of these new methods and advances in data collection are needed to enable a revolution in imaging magma reservoirs. 
    more » « less
  5. 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