Pre-eruptive magmatic processes re-timed using a non-isothermal approach to magma chamber dynamics

Open Source PaperThis work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/. The attached file is the published version of the article

Melt Inclusion Vapour Bubbles: The Hidden Reservoir for Major and Volatile Elements

Olivine-hosted melt inclusions (MIs) provide samples of magmatic liquids and their dissolved volatiles from deep within the plumbing system. Inevitable post-entrapment modifications can lead to significant compositional changes in the glass and/or any contained bubbles. Re-heating is a common technique to reverse MI crystallisation; however, its effect on volatile contents has been assumed to be minor. We test this assumption using crystallised and glassy basaltic MIs, combined with Raman spectroscopy and 3D imaging, to investigate the changes in fluid and solid phases in the bubbles before and after re-heating. Before re-heating, the bubble contains CO2 gas and anhydrite (CaSO4) crystallites. The rapid diffusion of major and volatile elements from the melt during re-heating creates new phases within the bubble: SO2, gypsum, Fe-sulphides. Vapour bubbles hosted in naturally glassy MIs similarly contain a plethora of solid phases (carbonates, sulphates, and sulphides) that account for up to 84% of the total MI sulphur, 80% of CO2, and 14% of FeO. In both re-heated and naturally glassy MIs, bubbles sequester major and volatile elements that are components of the total magmatic budget and represent a “loss” from the glass. Analyses of the glass alone significantly underestimates the original magma composition and storage parameters

Crystal and melt inclusion timescales reveal the evolution of magma migration before eruption

Volatile element concentrations measured in melt inclusions are a key tool used to understand magma migration and degassing, although their original values may be affected by different re-equilibration processes. Additionally, the inclusion-bearing crystals can have a wide range of origins and ages, further complicating the interpretation of magmatic processes. To clarify some of these issues, here we combined olivine diffusion chronometry and melt inclusion data from the 2008 eruption of Llaima volcano (Chile). We found that magma intrusion occurred about 4 years before the eruption at a minimum depth of approximately 8 km. Magma migration and reaction became shallower with time, and about 6 months before the eruption magma reached 3–4 km depth. This can be linked to reported seismicity and ash emissions. Although some ambiguities of interpretation still remain, crystal zoning and melt inclusion studies allow a more complete understanding of magma ascent, degassing, and volcano monitoring data.NRF (Natl Research Foundation, S’pore)MOE (Min. of Education, S’pore)Published versio

The evolution and storage of primitive melts in the Eastern Volcanic Zone of Iceland: the 10 ka Grímsvötn tephra series (i.e. the Saksunarvatn ash)

Major, trace and volatile elements were measured in a suite of primitive macrocrysts and melt inclusions from the thickest layer of the 10 ka Grímsvötn tephra series (i.e. Saksunarvatn ash) at Lake Hvítárvatn in central Iceland. In the absence of primitive tholeiitic eruptions (MgO > 7 wt.%) within the Eastern Volcanic Zone (EVZ) of Iceland, these crystal and inclusion compositions provide an important insight into magmatic processes in this volcanically productive region. Matrix glass compositions show strong similarities with glass compositions from the AD 1783–84 Laki eruption, confirming the affinity of the tephra series with the Grímsvötn volcanic system. Macrocrysts can be divided into a primitive assemblage of zoned macrocryst cores (An_78–An_92, Mg#_cpx = 82–87, Fo_79.5–Fo_87) and an evolved assemblage consisting of unzoned macrocrysts and the rims of zoned macrocrysts (An_60–An_68, Mg#_cpx = 71–78, Fo_70–Fo_76). Although the evolved assemblage is close to being in equilibrium with the matrix glass, trace element disequilibrium between primitive and evolved assemblages indicates that they were derived from different distributions of mantle melt compositions. Juxtaposition of disequilibrium assemblages probably occurred during disaggregation of incompatible trace element-depleted mushes (mean La/Yb_melt = 2.1) into aphyric and incompatible trace element-enriched liquids (La/Yb_melt = 3.6) shortly before the growth of the evolved macrocryst assemblage. Post-entrapment modification of plagioclase-hosted melt inclusions has been minimal and high-Mg# inclusions record differentiation and mixing of compositionally variable mantle melts that are amongst the most primitive liquids known from the EVZ. Coupled high field strength element (HFSE) depletion and incompatible trace element enrichment in a subset of primitive plagioclase-hosted melt inclusions can be accounted for by inclusion formation following plagioclase dissolution driven by interaction with plagioclase-undersaturated melts. Thermobarometric calculations indicate that final crystal-melt equilibration within the evolved assemblage occurred at ~1140°C and 0.0–1.5 kbar. Considering the large volume of the erupted tephra and textural evidence for rapid crystallisation of the evolved assemblage, 0.0–1.5 kbar is considered unlikely to represent a pressure of long-term magma accumulation and storage. Multiple thermometers indicate that the primitive assemblage crystallised at high temperatures of 1240–1300°C. Different barometers, however, return markedly different crystallisation depth estimates. Raw clinopyroxene-melt pressures of 5.5–7.5 kbar conflict with apparent melt inclusion entrapment pressures of 1.4 kbar. After applying a correction derived from published experimental data, clinopyroxene-melt equilibria return mid-crustal pressures of 4±1.5 kbar, which are consistent with pressures estimated from the major element content of primitive melt inclusions. Long-term storage of primitive magmas in the mid-crust implies that low CO_2 concentrations measured in primitive plagioclase-hosted inclusions (262–800 ppm) result from post-entrapment CO_2 loss during transport through the shallow crust. In order to reconstruct basaltic plumbing system geometries from petrological data with greater confidence, mineral-melt equilibrium models require refinement at pressures of magma storage in Iceland. Further basalt phase equilibria experiments are thus needed within the crucial 1–7 kbar range.D.A.N. was supported by a Natural Environment Research Council studentship (NE/1528277/1) at the start of this project. SIMS analyses were supported by Natural Environment Research Council Ion Microprobe Facility award (IMF508/1013).This is the final version of the article. It first appeared from Springer via http://dx.doi.org/10.1007/s00410-015-1170-

Fe-Mg interdiffusion rates in clinopyroxene: Experimental data and implications for Fe-Mg exchange geothermometers

Chemical interdiffusion of Fe-Mg along the c-axis [001] in natural diopside crystals (XDi = 0.93) was experimentally studied at ambient pressure, at temperatures ranging from 800 to 1,200 °C and oxygen fugacities from 10-11 to 10-17 bar. Diffusion couples were prepared by ablating an olivine (XFo = 0.3) target to deposit a thin film (20-100 nm) onto a polished surface of a natural, oriented diopside crystal using the pulsed laser deposition technique. After diffusion anneals, compositional depth profiles at the near surface region (~400 nm) were measured using Rutherford backscattering spectroscopy. In the experimental temperature and compositional range, no strong dependence of DFe-Mg on composition of clinopyroxene (Fe/Mg ratio between Di93-Di65) or oxygen fugacity could be detected within the resolution of the study. The lack of fO2-dependence may be related to the relatively high Al content of the crystals used in this study. Diffusion coefficients, DFe-Mg, can be described by a single Arrhenius relation with (Formula presented). DFe-Mg in clinopyroxene appears to be faster than diffusion involving Ca-species (e.g., DCa-Mg) while it is slower than DFe-Mg in other common mafic minerals (spinel, olivine, garnet, and orthopyroxene). As a consequence, diffusion in clinopyroxene may be the rate-limiting process for the freezing of many geothermometers, and compositional zoning in clinopyroxene may preserve records of a higher (compared to that preserved in other coexisting mafic minerals) temperature segment of the thermal history of a rock. In the absence of pervasive recrystallization, clinopyroxene grains will retain compositions from peak temperatures at their cores in most geological and planetary settings where peak temperatures did not exceed ~1,100 °C (e.g., resetting may be expected in slowly cooled mantle rocks, many plutonic mafic rocks, or ultra-high temperature metamorphic rocks)

Pre- and syn-eruptive degassing and crystallisation processes of the 2010 and 2006 eruptions of Merapi volcano, Indonesia

The 2010 eruption of Merapi (VEI 4) was the volcano’s largest since 1872. In contrast to the prolonged and effusive dome-forming eruptions typical of Merapi’s recent activity, the 2010 eruption began explosively, before a new dome was rapidly emplaced. This new dome was subsequently destroyed by explosions, generating pyroclastic density currents (PDCs), predominantly consisting of dark coloured, dense blocks of basaltic andesite dome lava. A shift towards open-vent conditions in the later stages of the eruption culminated in multiple explosions and the generation of PDCs with conspicuous grey scoria and white pumice clasts resulting from sub-plinian convective column collapse. This paper presents geochemical data for melt inclusions and their clinopyroxene hosts extracted from dense dome lava, grey scoria and white pumice generated during the peak of the 2010 eruption. These are compared with clinopyroxene-hosted melt inclusions from scoriaceous dome fragments from the prolonged dome-forming 2006 eruption, to elucidate any relationship between pre-eruptive degassing and crystallisation processes and eruptive style. Secondary ion mass spectrometry analysis of volatiles (H2O, CO2) and light lithophile elements (Li, B, Be) is augmented by electron microprobe analysis of major elements and volatiles (Cl, S, F) in melt inclusions and groundmass glass. Geobarometric analysis shows that the clinopyroxene phenocrysts crystallised at depths of up to 20 km, with the greatest calculated depths associated with phenocrysts from the white pumice. Based on their volatile contents, melt inclusions have re-equilibrated during shallower storage and/or ascent, at depths of ~0.6–9.7 km, where the Merapi magma system is interpreted to be highly interconnected and not formed of discrete magma reservoirs. Melt inclusions enriched in Li show uniform “buffered” Cl concentrations, indicating the presence of an exsolved brine phase. Boron-enriched inclusions also support the presence of a brine phase, which helped to stabilise B in the melt. Calculations based on S concentrations in melt inclusions and groundmass glass require a degassing melt volume of 0.36 km3 in order to produce the mass of SO2 emitted during the 2010 eruption. This volume is approximately an order of magnitude higher than the erupted magma (DRE) volume. The transition between the contrasting eruptive styles in 2010 and 2006 is linked to changes in magmatic flux and changes in degassing style, with the explosive activity in 2010 driven by an influx of deep magma, which overwhelmed the shallower magma system and ascended rapidly, accompanied by closed-system degassing