Search for: All records

Abstract Crystal‐hosted melt embayments and melt inclusions partially record magmatic processes at depth, but it is not always obvious how to interpret this record. One impediment is our incomplete understanding of how embayments and melt inclusions form. In this study, we investigate the formation mechanism of embayments and melt inclusions during quartz growth to quantify the relationship between the compositions of the entrapped and average melt. We study the growth of embayments and inclusions through direct numerical simulations that couple the growth of a crystal surface with the evolution of the concentrations of incompatible components in the surrounding melt. We find that H₂O is more enriched in the interior of defects on crystal surface compared to the exterior. The resultant lower disequilibrium in the defect interior causes lower growth rate than in the exterior, elongating the defect into an embayment. If crystal growth stops, the composition in the embayment equilibrates with the average melt within days to months. If crystal growth continues until the embayment neck closes, a melt inclusion forms. The melt entrapped by both embayments and melt inclusions is enriched in incompatible components, such as H₂O and CO₂. In addition to inclusion size, the enrichment of incompatible components in melt inclusions also depends on component diffusivity and the crystal growth regime. High‐diffusivity components like H₂O have similar enrichment levels in all scenarios, while lower‐diffusivity components like CO₂are more enriched in melt inclusions with smaller sizes or formed in continuous crystal growth. 

Abstract Crystals retain an imprint of the dynamic changes within a magma reservoir and hence contain invaluable information about the evolving conditions inside volcanic plumbing systems. However, instead of telling a single, simple story, they comprise overprinted evidence of numerous processes relating to temperature, pressure and composition that drive crystal precipitation and dissolution in magmatic systems. To decipher these different elements in the story that crystals tell, we attempt to identify the observational signatures of a simple, yet ubiquitous process: crystal precipitation and dissolution during magma cooling. To isolate this process in a complex magmatic system with intricate dynamic feedbacks, we assume that synthetic crystals precipitate and dissolve rapidly in response to deviations from thermodynamic equilibrium. In our crystalline‐scale simulations, synthetic crystals drag along the cooler‐than‐ambient melt in which they precipitated and can drive a temperature‐dependent, crystal‐driven convection. We analyze the non‐dimensional conditions for this coupled convection and record the heterogeneous thermal histories that synthetic crystals in this flow regime experience. We show that many synthetic crystals dissolve, loosing their thermal record of the convection. Based on our findings, we suggest that heterogeneity in the thermal history of crystals is more indicative of local, crystal‐scale processes than the overall, system‐wide cooling trend. 

Abstract Persistent volcanic activity is thought to be linked to degassing, but volatile transport at depth cannot be observed directly. Instead, we rely on indirect constraints, such as CO₂‐H₂O concentrations in melt inclusions trapped at different depth, but this data is rarely straight‐forward to interpret. In this study, we integrate a multiscale conduit‐flow model for non‐eruptive conditions and a volatile‐concentration model to compute synthetic profiles of volatile concentrations for different flow conditions and CO₂fluxing. We find that actively segregating bubbles in the flow enhance the mixing of volatile‐poor and volatile‐rich magma in vertical conduit segments, even if the radius of these bubbles is several orders of magnitude smaller than the width of the conduit. This finding suggests that magma mixing is common in volcanic systems when magma viscosities are low enough to allow for bubble segregation as born out by our comparison with melt‐inclusion data: Our simulations show that even a small degree of mixing leads to volatile concentration profiles that are much more comparable to observations than either open‐ or closed‐system degassing trends for both Stromboli and Mount Erebus. Our results also show that two of the main processes affecting observed volatile concentrations, magma mixing and CO₂fluxing, leave distinct observational signatures, suggesting that tracking them jointly could help better constrain changes in conduit flow. We argue that disaggregating melt‐inclusion data based on the eruptive behavior at the time could advance our understanding of how conduit flow changes with eruptive regimes. 

Abstract Many volcanoes around the world are persistently active with continuous degassing for years or even centuries, sometimes exceeding historic records. Such long‐term stability contrasts with short‐term instability, reflected in eruptive episodes that punctuate passive degassing. These two aspects of persistent activity, long‐term stability as opposed to short‐term instability, are often conceptualized through two distinct model frameworks: Exchange‐flow in volcanic conduits is commonly invoked to explain the long‐term thermal balance and sustained passive degassing, while the ascent of large gas slugs is called upon to understand explosive eruptions. While typically considered separately, we propose here that both flow processes could occur jointly in the conduits of persistently active volcanoes and in transient connections between subvolcanic melt lenses. To understand the dynamic interplay between exchange flow and slug ascent, we link analogue laboratory experiments with direct numerical simulations. We find that the two flows superimpose without creating major disruptions when only considering the ascent of a single gas slug. However, the sequential ascent of multiple gas slugs is disruptive to the ambient exchange flow, because it may entail continual buildup of buoyant magma at depth. While our study focuses on the laboratory scale, we propose that the dependence of exchange‐flow stability on sequential slug ascent is relevant for understanding why explosive sequences are sometimes followed by effusive eruptions. Taken together, our work suggests that integrating exchange flow and slug ascent could provide a more complete understanding of persistently active volcanoes than either model framework offers in isolation. 

Many volcanoes emit a significant portion of the gas they transport to the atmosphere during continual passive degassing rather than during eruptions. To maintain a high gas and thermal flux without erupting magma, the flow field in the volcanic conduit must be approximately balanced with gas-rich, buoyant magma ascending and degassed, heavy magma descending. In vertical conduits, this exchange flow takes the form of core–annular flow, where the gas-rich magma forms a core enclosed by an annulus of degassed magma. The flow dynamics of core–annular flow have been studied extensively in fluid dynamics, but mostly for constant material properties. Our study aims to advance our understanding of how core–annular flow responds to volatile exsolution – a simple, yet ubiquitous disruption in volcanic conduits, which alters both the density and the viscosity of the core fluid. By deriving an evolution equation for the core–annular interface based on a generalized exchange-flow condition using a lubrication approximation, we find that the response of the system to volatile exsolution depends on the conduit flow regime. The same nucleation event can lead to a flow adjustment only in the upper, only in the lower or in both portions of the volcanic conduit. Our results emphasize that the thermodynamic evolution of magma properties and volcanic conduit flow are intricately linked, which may help understand the observed variability of eruptive behaviour at persistently degassing volcanoes.