Geology topics

Erik H. Hauri

Publications and source records attributed to Erik H. Hauri.

4 recordsLinked to original sources

USGS2019

A high carbon content of the Hawaiian mantle from olivine-hosted melt inclusions

The deep mantle carbon content and flux are fundamental quantities in understanding global volatile cycles and distributions. Here, we present CO2 concentrations measured in 407 olivine-hosted melt inclusions from Hualalai, Kilauea, Koolau, Loihi, and Mauna Loa to constrain the Hawaiian mantle CO2 content and flux. Quantification of melt inclusion CO2 is complicated by the ubiquitous presence of vapor or “shrinkage” bubbles. The contribution from exsolved shrinkage bubble CO2 was determined from the measured bubble size and a CO2 equation of state, and added to the dissolved CO2 to reconstruct total melt inclusion CO2 concentrations. Bubbles typically contain ~90% of melt inclusion C, much of which may be sequestered in precipitated phases on bubble walls, and thus not amenable to measurement by Raman spectroscopy. Based on our dataset of total (dissolved + bubble) CO2 concentrations, we estimate that parental melts from the five Hawaiian volcanoes have CO2 concentrations ranging from 3900 to 10,000 ppm CO2. Among the active volcanoes, CO2 concentrations decrease to the northwest, likely reflecting mantle source heterogeneity, although differences in the depths of the magma chambers may also play a role. Mantle sources of the Hawaiian volcanoes range from 380 to 480 ppm CO2 suggesting that the Hawaiian plume is at least a factor of ~4 more C-rich than the MORB mantle. This enrichment is likely due to the presence of recycled surficial C and/or C-rich primitive material in the Hawaiian mantle.

Hawaii↗

USGS2018

When does eruption run-up begin? Multidisciplinary insight from the 1999 eruption of Shishaldin volcano

During the run-up to eruption, volcanoes often show geophysically detectable signs of unrest. However, there are long-standing challenges in interpreting the signals and evaluating the likelihood of eruption, especially during the early stages of volcanic unrest. Considerable insight can be gained from combined geochemical and geophysical studies. Here we take such an approach to better understand the beginning of eruption run-up, viewed through the lens of the 1999 sub-Plinian basaltic eruption of Shishaldin volcano, Alaska. The eruption is of interest due to its lack of observed deformation and its apparent long run-up time (9 months), following a deep long-period earthquake swarm. We evaluate the nature and timing of recharge by examining the composition of 138 olivine macrocrysts and 53 olivine-hosted melt inclusions and through shear-wave splitting analysis of regional earthquakes. Magma mixing is recorded in three crystal populations: a dominant population of evolved olivines (Fo 60-69 ) that are mostly reversely zoned, an intermediate population (Fo 69-76 ) with mixed zonation, and a small population of normally zoned more primitive olivines (Fo 76-80 ). Mixing-to-eruption timescales are obtained through modeling of Fe-Mg interdiffusion in 78 olivines. The large number of resultant timescales provides a thorough record of mixing, demonstrating at least three mixing events: a minor event ~9 months prior to eruption, coincident with the onset of deep long-period seismicity; a major event ~50 days before eruption, coincident with a large (M5.2) shallow earthquake; and a final event about a week prior to eruption. Shear-wave splitting analysis shows a change in the orientation of the local stress field about a month after the deep long-period swarm and around the time of the M5.2 event. Earthquake depths and vapor saturation pressures of Raman-reconstructed melt inclusions indicate that the recharge magma originated from depths of at least 20 km, and that mixing with a shallow magma or olivine cumulates occurred in or just below the edifice (<3 km depth). Deformation was likely outside the spatial and temporal resolution of the satellite measurements. Prior to eruption magma was stored over a large range of depths (~0-2.5 km below the summit), suggesting a shallow, vertical reservoir that could provide another explanation for the lack of detectable deformation. The earliest sign of unrest (deep long-period seismicity) coincides temporally with magmatic activity (magma mixing and a change in the local stress state), possibly indicating the beginning of eruption run-up. The more immediate run-up began with the major recharge event ~50 days prior to eruption, after which the signs of unrest became continuous. This timescale is long compared to the seismic run-up to other basaltic eruptions (typically hours to days). Other volcanoes classified as open-system, based on their lack of precursory deformation, also tend to have relatively long run-up durations, which may be related to the time required to fill the shallow reservoir with magmas sourced from greater depth.

Alaska↗

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↗