Search USGS⌕ Search

Geology topics

Terry Plank

Publications and source records attributed to Terry Plank.

5 recordsLinked to original sources
USGS2025

Volcanic gases reflect magma stalling and launching depths

Many open-vent arc volcanoes display two modes in their continuous gas emissions, one with a characteristic CO 2 / S T ratio typical of periods of quiescent degassing and another punctuated by high CO 2 / S T gas emitted in the weeks before eruption, a recently recognized eruption precursor. In this study we explore the origin of the two modes of degassing revealed by time-series gas data at Turrialba volcano (Costa Rica) in the context of new melt inclusion (MI) data. To reconstruct the c[CO 2 ] of undegassed magma, we developed a rapid-quench piston-cylinder assembly to rehomogenize the vapor bubble commonly contained in MIs. We focus on olivine-hosted MIs from a mafic scoria sample erupted from Turrialba in 1864–1866. The reconstructed CO 2 contents in MIs decrease from ∼4000 to <1000 ppmw as S contents decrease from 3500 to <1000 ppmw. The highest reconstructed S and CO 2 in the MIs resulted in an initial magmatic CO 2 / S T ratio (molar) of 0.83. Informed by the MI data, we modeled the decompression degassing of Turrialba magma and vapor composition using the Sulfur_X and EVo models. Instead of being controlled by initial magmatic CO 2 /S T ratio as suggested by previous studies, we find that the quiescent gas emitted from Turrialba during 2014–2018 (CO 2 / S T = 2.3 ± 0.8, molar) appears to reflectequilibrium with magmas stored at 4–8 km (Sulfur_X) or 2 km (EVo) depth, when H 2 O is degassing extensively from the magma. A magma storage region at 4–8 km is also supported by seismic tomography. The second gas mode is noted by spikes in CO 2 / S T ∼ 7.9 ± 2 in the weeks prior to eruption. This gas reflects equilibrium with a magma at 12–18 km (Sulfur_X) or 4–8 km (EVo), where the ascending magma is saturated with a CO 2 -rich vapor. Thus, there are two important trans crustal depths beneath the volcano: one where the rate of H 2 O loss from the magma and thus magma viscosity increases, and one at greater depths where high CO 2 /S T vapor forms and may facilitate dike propagation. We interpret the shallower, H 2 O-loss region as the main site of magma stalling and storage, where quiescent gas is generated continuously. We interpret the greater depth (12–18 km) as the source of the precursory gas that precedes eruption, and where the mafic melt lastly equilibrated with a mush zone before ascending and triggering eruption weeks later. This hypothesis is ripe for testing at other volcanoes that exhibit two modes in gas geochemistry.

Turrialba volcano↗
USGS2020

Vapor-bubble growth in olivine-hosted melt inclusions

Melt inclusions record the depth of magmatic processes, magma degassing paths, and volatile budgets of magmas. Extracting this information is a major challenge. It requires determining melt volatile contents at the time of entrapment when working with melt inclusions that have suffered post-entrapment modifications. Several processes decrease internal melt inclusion pressure, resulting in nucleation and growth of a vapor bubble and, time permitting, diffusion of volatiles (especially CO2) into the vapor bubble. Methods exist that attempt to reconstruct the entrapped CO2 contents, but they are difficult to apply and yield inconsistent results. Here, we explore bubble growth, evaluate CO2 reconstruction approaches, and develop improved experimental and computational approaches. Piston-cylinder experiments were conducted on olivine-hosted melt inclusions from Seguam (Alaska, USA) and Fuego (Guatemala) volcanoes at the following conditions: 500-800 MPa, 1140-1200 °C for Seguam and 1110-1140 °C for Fuego, 4-8 wt% H2O in the KBr brine, and run durations of 10-120 minutes. Heated melt inclusions form well-defined S-CO2 trends that can be described by degassing models. CO2 contents are enriched by a factor of ~2.5, on average, relative to those of the glasses within unheated melt inclusions, whereas S contents of heated and unheated melt inclusion glasses overlap, indicating insignificant amounts of S partition into the vapor bubble. Low closure temperatures enable CO2 diffusion into vapor bubbles during quench upon eruption, while a higher closure temperature for S limits its loss to vapor bubbles. We evaluate the timescales of post-entrapment processes and use the results to develop a new computational model to restore entrapped CO2 contents: MIMiC (Melt Inclusion Modification Corrections). Heated melt inclusion data are used as a benchmark to evaluate of the results from MIMiC and other published methods of CO2 reconstruction. The methods perform variably well. Key advantages to our experimental rehomogenization technique are that it enables accurate measurements of CO2 contents and allows for large quantities of melt inclusions to be rehomogenized efficiently. Our new computational model produces more accurate results than other computational methods, has similar accuracy to the Raman method of CO2 reconstruction in cases where Raman can be applied (i.e., no C-bearing phases in bubble), and can be applied to the vast body of published melt inclusion data. To obtain the most robust data on bubble-bearing melt inclusions, we recommend taking both experimental- and MIMiC-based approaches.

American Mineralogist↗
USGS2016

Magma decompression rates during explosive eruptions of Kīlauea volcano, Hawaii, recorded by melt embayments

The decompression rate of magma as it ascends during volcanic eruptions is an important but poorly constrained parameter that controls many of the processes that influence eruptive behavior. In this study, we quantify decompression rates for basaltic magmas using volatile diffusion in olivine-hosted melt tubes (embayments) for three contrasting eruptions of Kīlauea volcano, Hawaii. Incomplete exsolution of H 2 O, CO 2 , and S from the embayment melts during eruptive ascent creates diffusion profiles that can be measured using microanalytical techniques, and then modeled to infer the average decompression rate. We obtain average rates of ~0.05–0.45 MPa s −1 for eruptions ranging from Hawaiian style fountains to basaltic subplinian, with the more intense eruptions having higher rates. The ascent timescales for these magmas vary from around ~5 to ~36 min from depths of ~2 to ~4 km, respectively. Decompression-exsolution models based on the embayment data also allow for an estimate of the mass fraction of pre-existing exsolved volatiles within the magma body. In the eruptions studied, this varies from 0.1 to 3.2 wt% but does not appear to be the key control on eruptive intensity. Our results do not support a direct link between the concentration of pre-eruptive volatiles and eruptive intensity; rather, they suggest that for these eruptions, decompression rates are proportional to independent estimates of mass discharge rate. Although the intensity of eruptions is defined by the discharge rate, based on the currently available dataset of embayment analyses, it does not appear to scale linearly with average decompression rate. This study demonstrates the utility of the embayment method for providing quantitative constraints on magma ascent during explosive basaltic eruptions.

Hawaii↗
USGS2010

The role of water in generating the calc-alkaline trend: New volatile data for aleutian magmas and a new tholeiitic index

The origin of tholeiitic (TH) versus calc-alkaline (CA) magmatic trends has long been debated. Part of the problem stems from the lack of a quantitative measure for the way in which a magma evolves. Recognizing that the salient feature in many TH–CA discrimination diagrams is enrichment in Fe during magma evolution, we have developed a quantitative index of Fe enrichment, the Tholeiitic Index (THI): THI = Fe 4 ·0 /Fe 8 ·0 , where Fe 4 ·0 is the average FeO* concentration of samples with 4 ± 1 wt % MgO, and Fe 8 ·0 is the average FeO* at 8 ± 1 wt % MgO. Magmas with THI > 1 have enriched in FeO* during differentiation from basalts to andesites and are tholeiitic; magmas with THI < 1 are calc-alkaline. Most subduction zone volcanism is CA, but to varying extents; the THI expresses the continuum of Fe enrichment observed in magmatic suites in all tectonic settings. To test various controls on the development of CA trends, we present new magmatic water measurements in melt inclusions from eight volcanoes from the Aleutian volcanic arc (Augustine, Emmons, Shishaldin, Akutan, Unalaska, Okmok, Seguam, and Korovin). Least degassed H 2 O contents vary from ∼2 wt % (Shishaldin) to >7 wt % (Augustine), spanning the global range in arc mafic magmas. Within the Aleutian data, H 2 O correlates negatively with THI, from strongly calc-alkaline (Augustine, THI = 0·65) to moderately tholeiitic (Shishaldin, THI = 1·16). The relationship between THI and magmatic water is maintained when data are included from additional arc volcanoes, back-arc basins, ocean islands, and mid-ocean ridge basalts (MORBs), supporting a dominant role of magmatic water in generating CA trends. An effective break between TH and CA trends occurs at ∼2 wt % H 2 O. Both pMELTs calculations and laboratory experiments demonstrate that the observed co-variation of H 2 O and THI in arcs can be generated by the effect of H 2 O on the suppression of plagioclase and the relative enhancement of Fe-oxides on the liquid line of descent. The full THI–H 2 O array requires an increase in fO 2 with H 2 O, from ≤FMQ (where FMQ is the fayalite–magnetite–quartz buffer) in MORB to ∼ΔFMQ +0·5 to +2 in arcs, consistent with inferences from measured Fe and S species in glasses and melt inclusions. A curve fit to the data, H 2 O (wt % ± 1·2) = exp[(1·26 – THI)/0·32], may provide a useful tool for estimating the H 2 O content of magmas that are inaccessible to melt inclusion study.

Journal of Petrology↗