USGS Science⌕ Search
Home/Archive/USGS

USGS · 70023294

Adiabatic temperature changes of magma-gas mixtures during ascent and eruption

L.G. Mastin·M.S. Ghiorso
View original source ↗Export metadata ↓

Abstract

Most quantitative studies of flow dynamics in eruptive conduits during volcanic eruptions use a simplified energy equation that ignores either temperature changes, or the thermal effects of gas exsolution. In this paper we assess the effects of those simplifications by analyzing the influence of equilibrium gas exsolution and expansion on final temperatures, velocities, and liquid viscosities of magma–gas mixtures during adiabatic decompression. For a given initial pressure ( p 1 ), temperature ( T 1 ) and melt composition, the final temperature ( T f ) and velocity ( u max ) will vary depending on the degree to which friction and other irreversible processes reduce mechanical energy within the conduit. The final conditions range between two thermodynamic end members: (1) constant enthalpy ( dh =0), in which T f is maximal and no energy goes into lifting or acceleration; and (2) constant entropy ( ds =0), in which T f is minimal and maximum energy goes into lifting and acceleration. For ds =0, T 1 =900 °C and p 1 =200 MPa, a water-saturated albitic melt cools by ~200 °C during decompression, but only about 250 °C of this temperature decrease can be attributed to the energy of gas exsolution per se: the remainder results from expansion of gas that has already exsolved. For the same T 1 and p 1 , and dh =0, T f is 10–15 °C hotter than T 1 but is about 10–25 °C cooler than T f in similar calculations that ignore the energy of gas exsolution. For ds =0, p 1 =200 MPa and T 1 =9,000 °C, assuming that all the enthalpy change of decompression goes into kinetic energy, a water-saturated albitic mixture can theoretically accelerate to ~800 m/s. Similar calculations that ignore gas exsolution (but take into account gas expansion) give velocities about 10–15% higher. For the same T 1 , p I =200 MPa, and ds =0, the cooling associated with gas expansion and exsolution increases final melt viscosity more than 2.5 orders of magnitude. For dh =0, isenthalpic heating decreases final melt viscosity by about 0.7 orders of magnitude. Thermal effects of gas exsolution are responsible for less than 10% of these viscosity changes. Isenthalpic heating could significantly reduce flow resistance in eruptive conduits if heat generation were concentrated along conduit walls, where shearing is greatest. Isentropic cooling could enhance clast fragmentation in near-surface vents in cases where extremely rapid pressure drops reduce gas temperatures and chill the margins of expanding pyroclasts.

Explore related subjects

Contributions to Mineralogy and PetrologyArticle

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

L.G. Mastin, M.S. Ghiorso. 2001. Adiabatic temperature changes of magma-gas mixtures during ascent and eruption. https://doi.org/10.1007/s004100000210

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

USGS2026

Storage timescales and the crystal growth history recorded in the layered reservoir of the Unit 3 Kulanaokuaiki Tephra Member of the Uēkahuna Ash from Kīlauea volcano

The Kulanaokuaiki Tephra Unit 3 (900 C.E.) was a subplinian basaltic scoria eruption, the largest known to have occurred at Kīlauea’s summit. The initiation mechanism and cause(s) of this highly hazardous eruption remain poorly understood, particularly as this style of activity has not been observed in the modern period. We use olivine diffusion timescales from Fe-Mg and lithium zoning to assess magma storage timescales and evaluate evidence for mixing prior to eruption. Of the 200 crystals assessed, 55% are normally zoned for forsterite and 45% are non-zoned. Zoned crystals (> 1 mm) ( n = 96) have average compositions of Fo 89 cores and Fo 81 rims, while small crystals 0.5–1 mm ( n = 14) have Fo 85 cores and Fo 82 rims. Modeled Fo diffusion timescales range from < 1 to 12.6 years (median of 3.0 years). Lithium maps in 32 of the forsterite-zoned crystals reveal a variety of lithium zoning types, most of which we attribute to crystal growth. Many lithium profiles are coupled to phosphorus within the outer ~ 100 μm of the crystal ( n = 17), reflecting growth processes and precluding their use for diffusion studies. Overall, our Unit 3 olivine dataset displays no evidence of magmatic priming shortly prior to eruption, unlike prior observations for the basaltic subplinian eruption of the Keanakāko‘i Tephra Unit E (1650 C.E.). Thus, explosive eruption magmas may be stored for considerably longer than previously thought, and magmatic priming is not a necessary precursor for Kīlauea’s explosive eruptions. We hypothesize that gas-driven filter pressing from crystallization of microlites and macrocrysts in the reservoir slurry contributed to the explosivity of the Unit 3 eruption.

Hawaii↗
USGS2024

Shallow storage of the explosive Earthquake Flat Pyroclastics magma body, Okataina Volcanic Center, Taupo Volcanic Zone, New Zealand: Evidence from phase-equilibria experiments

Rhyolitic tuffs range widely in their crystal contents from nearly aphyric to crystal-rich, and their crystal cargoes inform concepts of upper crustal magma reservoirs. The Earthquake Flat pyroclastics (Okataina Volcanic Center, Taupo Volcanic Zone, New Zealand) are 10 km 3 of rhyolitic tuffs with abundant (~ 40 vol.%) plagioclase and quartz, minor biotite, hornblende, and orthopyroxene, and accessory Fe-Ti oxides, apatite, and zircon, set in high-silica rhyolitic glass. Major minerals form large, euhedral phenocrysts and abundant glomerocrysts with few disequilibrium textures excepting some faintly resorbed quartz. Plagioclase phenocrysts have thick rims of nearly constant composition near An 30 , and hornblende is weakly zoned or unzoned. The abundant and texturally complex mineral assemblage contrasts with the nearby (~ 25 km), nearly synchronous, but more voluminous and crystal-moderate rhyolite tuffs from Rotoiti caldera. New H 2 O-saturated phase-equilibria results on the erupted Earthquake Flat melt (glass) determine its co-saturation with the partial phenocryst assemblage of plagioclase, quartz, biotite, and Fe-Ti oxides at: 140 MPa, 755 ºC. These closely approximate the conditions of the pre-eruptive magma body assuming it was saturated with nearly pure H 2 O and at an f O 2 of ~ Ni–NiO. Absence of hornblende and orthopyroxene from the synthesized assemblages may result from those minerals being in a peritectic reaction relation with melt to produce biotite, so they would not grow from the liquid used as starting material. Experimental results on Rotoiti rhyolite (Nicholls et al. 1992) show that the two bodies resided at similar pressures, temperatures, and f O 2 s. Lower crystal abundance of the Rotoiti tuffs may result from slight compositional differences. We interpret that the Earthquake Flat pyroclastics were sourced from the crystal-rich periphery of a mushy reservoir system with the Rotoiti occupying a more melt-rich central location. Uncertain is whether this was a single intrusion zoned continuously in crystallinity, or discrete adjacent intrusions, but our results illustrate and quantify complexities of magma storage across relatively short distances.

North Island, Taupo Volcanic Zone↗
USGS2024

Constraining magma storage conditions of the Toba magmatic system: A plagioclase and amphibole perspective

Silicic magma reservoirs are responsible for producing the largest explosive eruptions in the geologic record. Petrologic and geochronological data provide evidence for these systems spending substantial periods of time (10 4 –10 5 yrs) within the upper crust prior to eruption; however, the long-term thermochemical evolution of these systems is not fully understood, as existing petrologic data make it challenging to quantify the time interval a magmatic system has spent at certain temperatures, or its “thermal history”. Here, we investigate the 74 ka Youngest Toba Tuff (YTT), one of the largest explosive eruptions in the geologic record, to better constrain the long-term thermal evolution of its magmatic system. We combine forward models of Sr diffusion in plagioclase and hornblende, mineral thermometry, and pre-existing trace-element evolution models to quantify the thermochemical evolution of the YTT magmatic system. We find that plagioclase crystals record decades to centuries of storage at temperatures 750 C, while hornblende records up to 6200 years at the same temperatures. Hornblende crystallizes at temperatures around 820 C and adjusting our diffusion modeling to this temperature results in no more than 900 years at initial crystallization conditions. Combined with previous trace-element modeling work, these results indicate that although there was chemical diversity for long durations in the YTT magma system sufficient to produce unique composition eruptive products, the entire system was experiencing a relatively similar thermal history that did not allow for large bodies of eruptible magma to be present for long periods ( 10 2 –10 3 years). Rather, we suggest that magmas within the YTT magmatic system were stored for long durations at thermal conditions where they were uneruptible and only remobilized within a few centuries prior to eruption.

Sumatra↗