The mass of volatiles degassed from volcanic eruptions is often estimated by comparing the volatile concentrations in undegassed glassy melt inclusions with the volatile concentrations in the ...degassed matrix glass. However, melt inclusions are prone to post-entrapment modification, including diffusive H+ loss through the host olivine crystal lattice which lowers the H2O content of the inclusion, and the degassing of CO2 into a bubble in response to cooling and crystallisation on the inclusion walls. Such bubbles are very common in olivine-hosted melt inclusions from the AD 1783–1784 Laki eruption, south–east Iceland. We have determined the CO2 content of these bubbles using micro-Raman spectroscopy, and the CO2 concentration in the glass by SIMS. Our results show that >90% of the total inclusion CO2 may be sequestered into the bubble, which demonstrates the importance of measuring the compositions of both vapour bubbles and the glass phase in melt inclusions. We reconstruct the deep degassing path of the Laki magma by using Nb as proxy for the undegassed CO2 content of the melt inclusions. The substantial CO2/Nb variation in the Laki melt inclusions (3.8–364) can be explained by concurrent crystallisation and CO2 degassing in the Laki magmatic system. We calculate the amount of CO2 lost from individual melt inclusions, assuming CO2/Nb ≈ 435 for enriched Icelandic mantle and CO2/Nb ≈ 171 for depleted mantle. Melt inclusions with the greatest saturation pressures have lost the least CO2 prior to inclusion trapping. At any given saturation pressure, the most enriched melt inclusions have lost the most CO2, while the most depleted inclusions have lost very little CO2. Enriched primary melts with high initial CO2 concentrations are therefore useful for investigating deep degassing behaviour in magmatic systems because a range of melt inclusion saturation pressures are recorded during crystallisation and degassing. Depleted melt inclusions with low initial CO2 concentrations remain vapour-undersaturated to shallow levels and cannot be used to constrain deep degassing behaviour. The cumulative CO2 mass release from the Laki magma is determined as a function of pressure and extent of crystallisation. Using an updated petrologic method that takes into account the diversity of primary melts and CO2 sequestration into vapour bubbles, we calculate the total mass of CO2 exsolved from the Laki magma to be 304 Mt.
•CO2 in shrinkage bubbles in melt inclusions can be measured by Raman spectroscopy.•>90% of the total melt inclusion CO2 may be contained in the bubble.•The deep degassing path of the Laki magma is reconstructed using CO2/Nb systematics.•Laki melt inclusions preserve a record of concurrent crystallisation and degassing.•The total mass of CO2 exsolved from the Laki magma is ∼304 Mt.
The Askja volcanic system, North Iceland, experienced a volcano‐tectonic episode between 1874 and 1876, the climax of which was a rhyolitic, phreatoplinian to Plinian eruption at Askja central ...volcano on 28–29 March 1875. Fissure eruptions also occurred in 1875, producing the Nýjahraun lava, 45–65 km north of Askja. The Nýjahraun basalt is indistinguishable, in terms of whole‐rock major elements, from the small‐volume basaltic eruptions that took place at Askja in the early 20th century. It has been suggested that all of these basalts originated from a shallow magma chamber beneath Askja, with the Nýjahraun eruptions being fed by northward‐propagating lateral dykes. It has also been conjectured that the Holuhraun lava, located at the southern tip of the Askja volcanic system 15–25 km south of Askja, was connected with the 1874–1876 Askja volcano‐tectonic episode. We re‐examine these interpretations in light of new whole‐rock, glass and melt inclusion analyses from samples collected along the length of the Askja volcanic system. Glasses from Nýjahraun and the Askja 20th century eruptions are geochemically distinct. We suggest that the Askja 20th century basalts mixed with evolved melts in the crust, while the Nýjahraun magma evolved without such interactions. The Holuhraun basalt is more similar to lavas erupted on the Bárðarbunga‐Veiðivötn volcanic system than to postglacial basalts from Askja, indicating that particular geochemical signatures are not necessarily confined to the tectonic or structural surface expression of single volcanic systems. This has important implications for the identification and delineation of individual volcanic systems beneath the northwest sector of Vatnajökull.
Key Points
We present new geochemical data for eruptions on the Askja volcanic system
Lateral flow is not required to explain fissure eruptions north of Askja in 1875
The Holuhraun lava is compositionally similar to Veidivotn basalts
Petrological constraints on the timescales of pre-eruptive crystal storage and magma degassing provide an important framework for the interpretation of seismic, geodetic and gas monitoring data in ...volcanically active regions. We have used Fe–Mg diffusion chronometry in 86 olivine macrocrysts from the AD 1783–1784 Laki eruption on Iceland's Eastern Volcanic Zone to characterise timescales of crystal storage and transport in the lead-up to this eruption. The majority of these olivines have core compositions of Fo<76, and rim compositions in the range Fo69–Fo74 that are close to equilibrium with the Laki melt. Diffusion modelling using the greyscale intensity of backscattered electron images as a proxy for olivine composition reveals that the most probable Fe–Mg diffusion timescale for Laki olivines is 7.8 days, which reflects the characteristic olivine residence time in the carrier melt prior to eruption. A small population of Fo>81 olivines record Fe–Mg diffusion timescales of ∼124 days; these crystals are likely to have formed in mid-crustal magma chambers, been transferred to storage at shallower levels and then entrained into the Laki melt prior to eruption. Typical Fe–Mg diffusion timescales of 6–10 days are shorter than the average time interval between discrete episodes of the Laki eruption, indicating variable or pulsed disaggregation of stored crystals into the carrier liquid prior to the onset of each episode. The diffusion timescales coincide with historical accounts of strong and frequent earthquakes in southeast Iceland, which we interpret as being associated with mush disaggregation related to melt withdrawal and the initiation of dyke propagation from a crustal magma reservoir at ∼6±3 km depth to the surface. We calculate pre-eruptive CO2 fluxes of 2–6 Mt d−1, assuming a pre-eruptive CO2 outgassing budget of 189.6 Mt for the Laki eruption and a constant rate of CO2 release in the 6–10 days preceding each eruptive episode. Our dataset indicates that petrological constraints on the timescales of magmatic processes occurring in the days leading up to historic eruptions may enhance our ability to forecast the onset of future large eruptions, both in Iceland and further afield.
•Fe–Mg diffusion in olivine is used to characterise timescales of crystal storage and transport.•Olivines were entrained into their carrier liquid 6–10 days prior to the Laki eruption.•Diffusion timescales are comparable to historical records of pre-eruptive seismicity.•CO2 fluxes in the 6–10 days preceding eruption were 2–6 Mt/d.•CO2 flux monitoring may enhance accuracy of eruption forecasting in Iceland.
The cycling of material from Earth's surface environment into its interior can couple mantle oxidation state to the evolution of the oceans and atmosphere. A major uncertainty in this exchange is ...whether altered oceanic crust entering subduction zones can carry the oxidised signal it inherits during alteration at the ridge into the deep mantle for long-term storage. Recycled oceanic crust may be entrained into mantle upwellings and melt under ocean islands, creating the potential for basalt chemistry to constrain solid Earth–hydrosphere redox coupling.
Numerous independent observations suggest that Iceland contains a significant recycled oceanic crustal component, making it an ideal locality to investigate links between redox proxies and geochemical indices of enrichment. We have interrogated the elemental, isotope and redox geochemistry of basalts from the Reykjanes Ridge, which forms a 700 km transect of the Iceland plume. Over this distance, geophysical and geochemical tracers of plume influence vary dramatically, with the basalts recording both long- and short-wavelength heterogeneity in the Iceland plume. We present new high-precision Fe-XANES measurements of Fe3+/∑Fe on a suite of 64 basalt glasses from the Reykjanes Ridge. These basalts exhibit positive correlations between Fe3+/∑Fe and trace element and isotopic signals of enrichment, and become progressively oxidised towards Iceland: fractionation-corrected Fe3+/∑Fe increases by ∼0.015 and ΔQFM by ∼0.2 log units. We rule out a role for sulfur degassing in creating this trend, and by considering various redox melting processes and metasomatic source enrichment mechanisms, conclude that an intrinsically oxidised component within the Icelandic mantle is required. Given the previous evidence for entrained oceanic crustal material within the Iceland plume, we consider this the most plausible carrier of the oxidised signal.
To determine the ferric iron content of the recycled component (Fe2O3source) we project observed liquid compositions to an estimate of Fe2O3 in the pure enriched endmember melt, and then apply simple fractional melting models, considering lherzolitic and pyroxenitic source mineralogies, to estimate Fe2O3(source) content. Propagating uncertainty through these steps, we obtain a range of Fe2O3(source) for the enriched melts (0.9–1.4 wt%) that is significantly greater than the ferric iron content of typical upper mantle lherzolites. This range of ferric iron contents is consistent with a hybridised lherzolite–basalt (pyroxenite) mantle component. The oxidised signal in enriched Icelandic basalts is therefore potential evidence for seafloor–hydrosphere interaction having oxidised ancient mid-ocean ridge crust, generating a return flux of oxygen into the deep mantle.
•Oxidation correlates with enrichment in Reykjanes Ridge basalts.•Oxidation and enrichment correlate on 1000 km and 10 km length scales.•The geochemistry of the enriched source is consistent with having been recycled.•Recycling of oxidised oceanic crust may contribute to the solid Earth oxygen cycle.
Clinopyroxene is a major rock forming mineral capable of incorporating diverse metal cations. As a consequence, clinopyroxene preserves valuable archives of magmatic processes. Understanding ...clinopyroxene is thus essential for understanding Earth’s wider chemical evolution. However, knowledge about the relative abundances of ferrous and ferric iron in magmatic clinopyroxene remains sparse because it is not currently possible to routinely measure the valence of iron in clinopyroxene crystals without either separating single crystals for bulk analysis or securing access to Mössbauer spectroscopy or a synchrotron radiation source to perform in-situ microanalysis. This is despite magmatic clinopyroxene crystals often containing appreciable quantities of ferric iron that affect its stability and behaviour in currently ill-constrained ways and limit our ability to exploit its chemistry to robustly reconstruct conditions of magma storage and evolution. Here we integrate optimised electron probe microanalysis and Mössbauer spectroscopy on endmember and single-crystal clinopyroxene samples to re-evaluate previously discredited approaches for estimating clinopyroxene ferric iron contents by stoichiometry. By ensuring that we measured all major and minor elements in clinopyroxene crystals with sufficient precision, we show that it is possible to readily obtain stoichiometric estimates of clinopyroxene ferric-to-total iron ratios with similar precisions to those derived from Mössbauer spectroscopy (1
σ
∼
3.5% absolute). Being able to robustly determine clinopyroxene ferric iron contents enables us to propose a new empirical scheme for assigning clinopyroxene components that explicitly accounts for ferric iron, which is primarily hosted within esseneite component (CaFe
3
+
AlSiO
6
) in clinopyroxenes dominated by quadrilateral components and aegirine component (NaFe
3
+
Si
2
O
6
) in alkali clinopyroxenes. Our new scheme provides a framework for documenting the full spectrum of clinopyroxene compositions in both natural and experimental systems when analyses are performed with sufficient precision.
The timescales of magmatic processes within a volcanic system may be variable over a volcano's geological history. Crystals reflect environmental perturbations under which they grew, and ...compositional gradients quenched inside crystals on eruption can be exploited to extract timescales of magmatic processes. Here, we use multi-element diffusion chronometry in olivine macrocrysts to recover their residence time in a melt that ultimately erupted at the surface. The macrocrysts were mobilized by the carrier melt from mushy layers in the magma reservoir, and diffusion timescales likely reflect the time interval between mush disaggregation, ascent, and eruption. To unravel the evolution of mush disaggregation timescales with time, we target early-Holocene, middle-Holocene, and historical magmatic units erupted in the Bardarbunga-Veidivotn volcanic system in Iceland's Eastern Volcanic Zone. Macrocryst contents vary between samples; early-Holocene samples are highly phyric (10-45 vol% macrocrysts) and contain gabbroic nodules, whereas middle-Holocene (5-15 vol%) and historical units (5-10 vol%) tend to be generally less phyric. Early-Holocene olivine macrocrysts have core compositions in the range Fo84-87, while middle-Holocene and historical samples record a wider range in core compositions from Fo80 to Fo86.5. Olivine rims are in chemical equilibrium with their carrier liquid and are slightly more evolved in early-Holocene units (Fo76-81) compared to middle-Holocene (Fo78-80) and historical (Fo81-83) units. Diffusion chronometry reveals that the timescale between mush dis-aggregation and eruption has changed over time, with timescales getting shorter approaching recent times. Early-Holocene olivine macrocrysts dominantly record Fe-Mg diffusion timescales between 200-400 days, while middle-Holocene and historical units typically record timescales of about 70 and 60 days, respectively. Barometric studies suggest that melts and crystals are likely stored and gradually transferred throughout an interconnected multi-tiered system that ultimately culminates in a mid-crustal reservoir(s) at about 6.8-7.5±2.5 km depth, where final disaggregation by the carrier liquid took place. We argue that, as a result of extensional processes enhanced by rifting events, well-mixed melts got drawn into mid-crustal reservoir(s), causing crystal mush loosening and mobilization. In addition, we propose that more energy in the form of heat and/or melt supply was required in the early-Holocene to break up the dense mush fabric and convert it into an eruptible magma. Conversely, as evidenced by the diverse macrocryst content of the historical units and by the lack of gabbroic nodules, the system has become characterized by a less compact mush fabric since at least the middle-Holocene, such that fresh injection of melt would easily loosen and mobilize the mush, resulting in an eruption within a couple of months. This study provides evidence that along axial rift settings, rifting-related processes can help to "pull the mush apart" with no requirement for primitive magma injection as an eruption trigger. Furthermore, we provide evidence that in the Bardarbunga-Veidivotn volcanic system specifically, the time between mush disaggregation and eruption has decreased considerably with time, indicating shorter warning times before imminent eruptions.
The 2014–2015 Holuhraun eruption, on the Bárðarbunga volcanic system in central Iceland, was one of the best-monitored basaltic fissure eruptions that has ever occurred, and presents a unique ...opportunity to link petrological and geochemical data with geophysical observations during a major rifting episode. We present major and trace element analyses of melt inclusions and matrix glasses from a suite of ten samples collected over the course of the Holuhraun eruption. The diversity of trace element ratios such as La/Yb in Holuhraun melt inclusions reveals that the magma evolved via concurrent mixing and crystallization of diverse primary melts in the mid-crust. Using olivine–plagioclase–augite–melt (OPAM) barometry, we calculate that the Holuhraun carrier melt equilibrated at 2.1 ± 0.7 kbar (7.5 ± 2.5 km), which is in agreement with the depths of earthquakes (6 ± 1 km) between Bárðarbunga central volcano and the eruption site in the days preceding eruption onset. Using the same approach, melt inclusions equilibrated at pressures between 0.5 and 8.0 kbar, with the most probable pressure being 3.2 kbar. Diffusion chronometry reveals minimum residence timescales of 1–12 days for melt inclusion-bearing macrocrysts in the Holuhraun carrier melt. By combining timescales of diffusive dehydration of melt inclusions with the calculated pressure of H
2
O saturation for the Holuhraun magma, we calculate indicative magma ascent rates of 0.12–0.29 m s
−1
. Our petrological and geochemical data are consistent with lateral magma transport from Bárðarbunga volcano to the eruption site in a shallow- to mid-crustal dyke, as has been suggested on the basis of seismic and geodetic datasets. This result is a significant step forward in reconciling petrological and geophysical interpretations of magma transport during volcano-tectonic episodes, and provides a critical framework for the interpretation of premonitory seismic and geodetic data in volcanically active regions.
The valence state of Fe plays a vital role in setting and recording the oxidation state of magmas, commonly expressed in terms of oxygen fugacity (
f
O
2
). However, our knowledge of how and why
f
O
...2
varies within and between magmatic systems remains patchy because of diverse challenges associated with estimating the valence state of Fe in glasses and minerals routinely. Here we investigate Fe valence systematics in magmatic clinopyroxene crystals from ocean island basalts (OIBs) erupted in Iceland and the Azores to explore whether they record information about magma Fe
3
+
contents and magmatic
f
O
2
conditions. Although many studies assume that all Fe in augitic clinopyroxene crystals from OIBs occurs as Fe
2
+
, we find that up to half of the total Fe present can occur as Fe
3
+
, with crystals from alkali systems typically containing more Fe
3
+
than those from tholeiitic systems. Thus, Fe
3
+
is a major if under-appreciated constituent of augitic clinopyroxene crystals erupted from ocean island volcanoes. Most Fe
3
+
in these crystals is hosted within esseneite component (CaFe
3
+
AlSiO
6
), though some may be hosted in aegirine component (NaFe
3
+
Si
2
O
6
) in crystals from alkali systems. Observations from samples containing quenched matrix glasses suggest that the incorporation of Fe
3
+
is related to the abundance of tetrahedrally coordinated Al (
IV
Al), implying some steric constraints over Fe
3
+
partitioning between clinopyroxene and liquid (i.e.,
D
Fe
2
O
3
cpx
-
liq
values), though this may not be an equilibrium relationship. For example,
IV
Al-rich
{
h
k
0
}
prism sectors of sector-zoned crystals contain more Fe
3
+
than
IV
Al-poor
{
1
¯
11
}
hourglass sectors. Moreover,
IV
Al-rich compositions formed during disequilibrium crystallisation are enriched in Fe
3
+
. Apparent clinopyroxene-liquid Fe
2
+
–Mg exchange equilibria (i.e.,
K
D
,
Fe
2
+
-
Mg
cpx
-
liq
values) are similarly affected by disequilibrium crystallisation in our samples. Nonetheless, it is possible to reconcile our observed clinopyroxene compositions with glass Fe valence systematics estimated from olivine-liquid equilibria if we assume that
K
D
,
Fe
2
+
-
Mg
cpx
-
liq
values lies closer to experimentally reported values of 0.24
-
0.26 than values of
∼
0.28 returned from a general model. In this case, olivine-liquid and clinopyroxene-liquid equilibria record equivalent narratives, with one of our glassy samples from Iceland recording evolution under
f
O
2
conditions about one log unit above fayalite-magnetite-quartz (FMQ) equilibrium (i.e.,
∼
FMQ+1) and our glassy Azorean sample recording evolution under significantly more oxidising conditions (
≥
FMQ+2.5) before experiencing syn-eruptive reduction, likely as a result of SO
2
degassing; our other glassy sample from Iceland was also affected by reductive SO
2
degassing. Overall, our findings demonstrate that the Fe valence systematics of clinopyroxene crystals can record information about the conditions under which OIBs evolve, but that further experimental work is required to properly disentangle the effects of magma composition, disequilibrium and
f
O
2
conditions on clinopyroxene-liquid equilibria involving Fe
2
+
and Fe
3
+
.
The majority of basaltic magmas stall in the Earth's crust as a result of the rheological evolution caused by crystallization during transport. However, the relationships between crystallinity, ...rheology and eruptibility remain uncertain because it is difficult to observe dynamic magma crystallization in real time. Here, we present in-situ 4D data for crystal growth kinetics and the textural evolution of pyroxene during crystallization of trachybasaltic magmas in high-temperature experiments under water-saturated conditions at crustal pressures. We observe dendritic growth of pyroxene on initially euhedral cores, and a surprisingly rapid increase in crystal fraction and aspect ratio at undercooling ≥30 °C. Rapid dendritic crystallization favours a rheological transition from Newtonian to non-Newtonian behaviour within minutes. We use a numerical model to quantify the impact of rapid dendritic crystallization on basaltic dike propagation, and demonstrate its dramatic effect on magma mobility and eruptibility. Our results provide insights into the processes that control whether intrusions lead to eruption or not.