Engineering topics

S. Matthews

Publications and source records attributed to S. Matthews.

3 records

NASA NTRS2022

VESIcal: A Critical Approach to Volatile Solubility Modelling Using the Open-Source Engine Vesical

Accurate models of H(2)O and CO(2) solubility in silicate melts are vital for understanding volcanic plumbing systems. These models are used to estimate the depths of magma storage regions from melt inclusion volatile contents, investigate the role of volatile exsolution as a driver of volcanic eruptions, and track the degassing path followed by a magma ascending to the surface. However, despite the large increase in the number of experimental constraints over the last two decades, many recent studies still utilize an earlier generation of models which were calibrated on experimental datasets with restricted compositional ranges. This may be because many of the available tools for more recent models require large numbers of input parameters to be hand-typed (e.g., temperature, concentrations of H(2)O, CO(2), and 8–14 oxides), making them difficult to implement on large datasets. Here, we use a new open-source Python3 tool, VESIcal, to critically evaluate the behaviors and sensitivities of different solubility models for a range of melt compositions. Using literature datasets of andesitic-dacitic experimental products and melt inclusions as case studies, we illustrate the importance of evaluating the calibration dataset of each model. Finally, we highlight the limitations of particular data presentation methods, such as isobar diagrams, and provide suggestions for alternatives, and best practices regarding the presentation and archiving of data. This review will aid the selection of the most applicable solubility model for different melt compositions, and identifies areas where additional experimental constraints on volatile solubility are required.

magma↗

NASA NTRS2021

VESIcal Part I: An open-source thermodynamic model engine for mixed volatile (H2O-CO2) solubility in silicate melts

Thermodynamics has been fundamental to the interpretation of geologic data and modeling of geologic systems for decades. However, more recent advancements in computational capabilities and a marked increase in researchers’ accessibility to computing tools has outpaced the functionality and extensibility of currently available modeling tools. Here we present VESIcal (Volatile Equilibria and Saturation Identification calculator): the first comprehensive modeling tool for H 2 O, CO 2 , and mixed (H 2 O-CO 2 ) solubility in silicate melts that: a) allows users access to seven commonly used models, plus easy inter-comparison between models; b) provides universal functionality for all models (e.g., functions for calculating saturation pressures, degassing paths, etc.); c) can process large datasets (1,000’s of samples) automatically; d) can output computed data into an excel spreadsheet for simple post-modeling analysis; e) integrates advanced plotting capabilities directly within the tool; and f) provides all of these within the framework of a python library, making the tool extensible by the user and allowing any of the model functions to be incorporated into any other code capable of calling python. The tool is presented within this manuscript, which is a Jupyter notebook containing worked examples accessible to python users with a range of skill levels. The basic functions of VESIcal can also be access via a web app (https://vesical.anvil.app). The VESIcal python library is open-source and available for download at https://github.com/kaylai/VESIcal.

K. Iacovino↗

NASA NTRS2020

Reconstructing Magma Storage Depths for the 2018 Kilauean Eruption from melt inclusion CO2 Contents: The importance of Vapor Bubbles

The 2018 Lower East Rift Zone (LERZ) eruption of Kīlauea Volcano and the accompanying collapse of the summit caldera marked the most destructive phase of activity on Hawai’i in the last 200 years. The integration of petrological data extracted from lava samples collected throughout the eruption with geodetic data examining the caldera collapse event, and estimates of the co-erupted flux of SO2 from the main eruptive fissure (Fissure 8), provides an exceptional opportunity to determine the reservoir geometry and magma transport paths supplying Kīlauea’s LERZ. The forsterite contents of erupted olivines and the degree of disequilibrium with their carrier melts indicate that two distinct olivine populations were erupted from Fissure 8. Melt inclusion entrapment pressures reveal that more evolved olivines (Fo<81.5) crystallized at ~2 km depth within the shallower Halema’uma’u reservoir, while more primitive olivines (Fo>81.5)crystallized within the deeper South Caldera reservoir at ~3–5 km depth. Crucially, primitive olivines experienced extensive post-entrapment crystallization, driving the growth of a vapor bubble. Raman spectroscopy reveals that this bubble contains up to 99% of the total inclusionCO2 budget (median=93%). Measurements of CO2 in only the glass phase would have underestimated entrapment depths by up to 60× (median=11×), and the importance of the SC reservoir as a source of magma to Fissure 8 would have been overlooked. Overall, we demonstrate that Raman measurements of bubbles, along with careful choice of suitably-calibrated H2O-CO2 solubility model, is vital to place accurate constraints on the depths of magma storage regions supplying volcanic eruptions.

SIMS↗