Geology Reports⌕ Search

Geology topics

Terry Plank

Publications and source records attributed to Terry Plank.

8 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↗
USGS2023

Tracking carbon from subduction to outgassing along the Aleutian-Alaska Volcanic Arc

Subduction transports volatiles between Earth’s mantle, crust, and atmosphere, ultimately creating a habitable Earth. We use isotopes to track carbon from subduction to outgassing along the Aleutian-Alaska Arc. We find substantial along-strike variations in the isotopic composition of volcanic gases, explained by different recycling efficiencies of subducting carbon to the atmosphere via arc volcanism and modulated by subduction character. Fast and cool subduction facilitates recycling of ~43 to 61% sediment-derived organic carbon to the atmosphere through degassing of central Aleutian volcanoes, while slow and warm subduction favors forearc sediment removal, leading to recycling of ~6 to 9% altered oceanic crust carbon to the atmosphere through degassing of western Aleutian volcanoes. These results indicate that less carbon is returned to the deep mantle than previously thought and that subducting organic carbon is not a reliable atmospheric carbon sink over subduction time scales.

Alaska↗
USGS2021

Volcanic seismicity beneath Chuginadak Island, Alaska (Cleveland and Tana volcanoes): Implications for magma dynamics and eruption forecasting

Cleveland and Tana are remote volcanoes located in the central Aleutian volcanic arc on the eastern end of the Islands of Four Mountains (IFM). The persistently active Mount Cleveland volcano, on the western side of Chuginadak Island, is surrounded by several closely spaced Quaternary volcanic centers including Carlisle, Herbert, Kagamil, Tana, and Uliaga, and numerous small satellite vents on Chiginadak between Cleveland and Tana. The Alaska Volcano Observatory (AVO) installed two permanent broadband seismometers on Chuginadak Island in 2014, and we operated a temporary broadband network focused on the western side of the island in 2015–2016. Collectively, these stations provided the first seismic observations of this frequently active volcano and the surrounding Holocene-aged volcanic vents. During the study period (July 2014–January 2019), eruptive activity at Cleveland was characterized by small explosions separated by periods of lava effusion that formed small domes in the volcano's summit crater. We characterize seismicity beneath Chuginadak Island through automated analysis of event waveform frequency content, development of a one-dimensional P-wave velocity model, calculation of earthquake hypocenters, magnitudes, focal mechanisms, and identification of earthquake families. This analysis reveals the full range of seismic event types expected in a highly active volcanic environment and includes Volcano-Tectonic (VT) earthquakes, Long-Period (LP) events, and explosion signals. LP events appear to cluster at shallow depth beneath the active crater of Mount Cleveland and almost all of the explosions occur without identifiable short-term (hours to days) seismic precursors. VT earthquakes beneath Mount Cleveland occur at depths of 2 to 8 km below sea level (BSL) and range in magnitude from −0.2 to 1.8. VT focal mechanisms have horizontal P-axes that align with the regional axis of maximum stress. These observations, and a relatively slow one-dimensional seismic velocity model, are consistent with a shallow body of magma that is fed through a deeper conduit system. The time-history of VT earthquakes and shallow LP events suggest their occurrence may track the transfer of magma and fluids from the mid-crust to the shallow portions of the conduit system and may provide a means to anticipate future explosions and periods of dome growth. VT hypocenters also extend ~7 km northeast of Cleveland's summit at depths of 5 to 10 km BSL, under a group of Holocene-aged vents between Mount Cleveland and Tana. These earthquakes have vertically-oriented P-axes and a greater percentage occur in families. These observations, combined with observations of vent orientation and morphology and gas flux, suggest the area between Cleveland and Tana represents a zone of complicated volcano-tectonic interaction, similar to calderas elsewhere in the Aleutian arc. The presence of a larger volcanic system in the eastern IFM could influence magmatism and account for the multiple closely spaced volcanic centers in this region.

Alaska↗
USGS2020

Linking subsurface to surface using gas emission and melt inclusion data at Mount Cleveland volcano, Alaska

Mount Cleveland is one of Alaska's most active volcanoes, yet little is known about the magmatic system driving persistent and dynamic volcanic activity. Volcanic gas and melt inclusion (MI) data from 2016 were combined to investigate shallow magmatic processes. SO 2 emission rates were between 166 and 324 t/day and the H 2 O/SO 2 was 600 ± 53, whereas CO 2 and H 2 S were below detection. Olivine‐, clinopyroxene‐, and plagioclase‐hosted MIs have up to 3.8 wt.% H 2 O, 514 ppm CO 2 , and 2,320 ppm S. Equilibration depths, based on MI H 2 O contents, suggest that a magmatic column extended from 0.5 to 3.0 km (~10–60 MPa). We used MI data to empirically model open‐system H‐C‐S degassing from 0 to 12 km and found that a column of magma between 0.5 and 3 km could produce the measured gas H 2 O/SO 2 ratio. However, additional magma deeper than 3 km is required to sustain emissions over periods greater than days to weeks, if the observed vent dimension is a valid proxy for the conduit. Assuming an initial S content of 2,320 ppm, the total magma supply needed to sustain the annual SO 2 flux was 5 to 9.8 Mm 3 /yr, suggesting a maximum intrusive‐to‐extrusive ratio of 13:1. The model predicts degassing of <50 t/day CO 2 for July 2016, which corresponds to a maximum predicted CO 2 /SO 2 of 0.2. Ultimately, frequent recharge from deeper, less degassed magma is required to drive the continuous activity observed over multiple years. During periods of recharge we would expect lower H 2 O/SO 2 and measurable volcanic CO 2 .

Alaska↗
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↗