Clinopyroxene Dissolution Records Rapid Magma Ascent

Magma ascent rates control volcanic eruption styles. However, the rates at which basaltic magmas ascend through the crust remain highly uncertain. Although recent studies have successfully exploited records of decompression driven degassing to estimate the rates at which H2O-rich basalts ascend, such approaches cannot readily be applied to primitive and H2O-poor basalts that erupt in ocean island and mid-ocean ridge settings. Here we present magma ascent rates obtained by modeling the dissolution of clinopyroxene crystals in a wehrlitic nodule from the primitive Borgarhraun lava flow in North Iceland. High-Al2O3 clinopyroxene core compositions are consistent with crystallization near the Moho (~800 MPa), whereas low-Al2O3 clinopyroxene rims and inclusion compositions are consistent with crystallization at or near the surface. We interpret low-Al2O3 rims and inclusions as the crystallized remnants of boundary layers formed by the dissolution of high-Al2O3 clinopyroxene during magma ascent. By combining characteristic rim dissolution lengths of 50–100 μm with published experimental calibrations of clinopyroxene dissolution behavior, we estimate that the Borgarhraun magma most likely decompressed and ascended at rates of 3.0–15 kPa.s−1 and 0.11–0.53 m.s−1, respectively. These rates are slightly faster than published estimates obtained by modeling the diffusive re-equilibration of olivine crystals, suggesting that the Borgarhraun magma either accelerated upwards or that it stalled briefly at depth prior to final ascent. Comparisons with other basaltic eruptions indicate that the H2O-poor magma that fed the dominantly effusive Borgarhraun eruption ascended at a similar rate to some H2O-rich magmas that have fed explosive eruptions in arc settings. Thus, magma ascent rates do not appear to correlate simply with magma H2O contents. Overall, our findings confirm that primitive and H2O-poor basalts can traverse the crust within days, and may erupt with little precursory warning of magma ascent

Mantle-derived trace element variability in olivines and their melt inclusions

Trace element variability in oceanic basalts is commonly used to constrain the physics of mantle melting and the chemistry of Earth's deep interior. However, the geochemical properties of mantle melts are often overprinted by mixing and crystallisation processes during ascent and storage. Studying primitive melt inclusions offers one solution to this problem, but the fidelity of the melt-inclusion archive to bulk magma chemistry has been repeatedly questioned. To provide a novel check of the melt inclusion record, we present new major and trace element analyses from olivine macrocrysts in the products of two geographically proximal, yet compositionally distinct, primitive eruptions from the Reykjanes Peninsula of Iceland. By combining these macrocryst analyses with new and published melt inclusion analyses we demonstrate that olivines have similar patterns of incompatible trace element (ITE) variability to the inclusions they host, capturing chemical systematics on intra- and inter-eruption scales. ITE variability (element concentrations, ratios, variances and variance ratios) in olivines from the ITE-enriched Stapafell eruption is best accounted for by olivine-dominated fractional crystallisation. In contrast, ITE variability in olivines and inclusions from the ITE-depleted Háleyjabunga eruption cannot be explained by crystallisation alone, and must have originated in the mantle. Compatible trace element (CTE) variability is best described by crystallisation processes in both eruptions. Modest correlations between host and inclusion ITE contents in samples from Háleyjabunga suggest that melt inclusions can be faithful archives of melting and magmatic processes. It also indicates that degrees of ITE enrichment can be estimated from olivines directly when melt inclusion and matrix glass records of geochemical variability are poor or absent. Inter-eruption differences in olivine ITE systematics between Stapafell and Háleyjabunga mirror differences in melt inclusion suites, and confirm that the Stapafell eruption was fed by lower degree melts from greater depths within the melting region than the Háleyjabunga eruption. Although olivine macrocrysts from Stapafell are slightly richer in Ni than those from Háleyjabunga, their overall CTE systematics (e.g., Ni/(Mg/Fe), Fe/Mn and Zn/Fe) are inconsistent with being derived from olivine-free pyroxenites. However, the major element systematics of Icelandic basalts require lithological heterogeneity in their mantle source in the form of Fe-rich and hence fusible domains. We thus conclude that enriched heterogeneities in the Icelandic mantle are composed of modally enriched, yet nonetheless olivine-bearing, lithologies and that olivine CTE contents provide an incomplete record of lithological heterogeneity in the mantle. Modally enriched peridotites may therefore play a more important role in oceanic magma genesis than previously inferred.</p

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–1784 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 (An78–An92, Mg#cpx = 82–87, Fo79.5–Fo87) and an evolved assemblage consisting of unzoned macrocrysts and the rims of zoned macrocrysts (An60–An68, Mg#cpx = 71–78, Fo70–Fo76). 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/Ybmelt = 2.1) into aphyric and incompatible trace element-enriched liquids (La/Ybmelt = 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 CO2 concentrations measured in primitive plagioclase-hosted inclusions (262–800 ppm) result from post-entrapment CO2 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

Geochemical constraints on basalt petrogenesis in the Strait of Sicily Rift Zone (Italy): Insights into the importance of short lengthscale mantle heterogeneity

Igneous activity from the late Miocene to historic time (most recently 1891 CE) in the Strait of Sicily has created two volcanic islands (Pantelleria and Linosa) and several seamounts. These volcanoes are dominated by transitional (ol + hy-normative) to alkaline (ne-normative) basaltic lavas and scoriae; volcanic felsic rocks (peralkaline trachyte-rhyolite) crop out only on Pantelleria. Although most likely erupted through continental crust, basalts demonstrate no evidence of crustal contamination and are geochemically similar to oceanic island basalts (OIB). Despite their isotopic similarities, there are considerable compositional differences with respect to major and trace element geochemistry both between and within the two islands that are due to short-length scale mantle heterogeneity beneath the region as well as variability in partial melting and magma storage conditions. Published geophysical surveys suggest that lithospheric thickness beneath both islands is ~60 km; this is consistent with the results of our geochemical modelling (59\u201360 km), which also suggest mantle potential temperatures between 1415 and 1435 \ub0C, similar to those documented in other continental passive rifts. Trace element and isotopic data reveal that the asthenosphere beneath the Strait of Sicily is heterogenous at both interisland (100s of km) and intra-island (10s of km) scales. Although there is some compositional overlap between the two major synthems at Linosa, in general the older magmas (Arena Bianca, 700 ka) formed as a result of ~5% partial melting of a depleted MORB mantle (DMM) source enriched with a relatively small amount of recycled MORB material, whereas the younger magmas (Monte Bandiera, 530 ka) formed as a result of ~2% partial melting of a similar mantle source. Pantelleria magmas formed from a higher degree (~6%) of partial melting of a DMM source with a relatively greater amount of recycled MORB material and possibly other components. Geochemical modelling also suggests the older magmas on Linosa differentiated at a much shallower level (~8 km) than the younger magmas (~25 km, at or below the base of the crust) prior to eruption. Magmas stored in higher-level reservoirs were effectively homogenized and preserve a narrower compositional range than magmas sourced from depth. Data for the seamounts are scarce and compromised by significant seawater alteration; thus, these volcanic centers cannot be modelled but based on comparative geochemistry with the islands are likely the result of even smaller (&lt; 2%) degrees of partial melting beneath thicker (&gt; 60 km) lithosphere. Despite the geophysical similarities between the two islands in terms of lithospheric thickness and crustal thinning, melt productivity has been greater at Pantelleria, producing a much larger island and sustaining felsic magmatism, which we hypothesize may ultimately be entirely due to the local occurrence of much more fusible mantle

Geochemical Constraints on Mantle Sources and Basalt Petrogenesis in the Strait of Sicily Rift Zone (Italy): Insights into the Importance of Short Lengthscale Mantle Heterogeneity.

Igneous activity from the late Miocene to historic time (most recently 1891 ce) in the Strait of Sicily has created two islands (Pantelleria and Linosa) and several seamounts. These volcanoes are dominated by transitional (ol+hy-normative) to alkaline (ne-normative) basaltic lavas and scoriae; peralkaline felsic rocks (trachyte-rhyolite) crop out only on Pantelleria. Although most likely erupted through continental crust, basalts demonstrate no evidence of crustal contamination and are geochemically similar to oceanic island basalts (OIB). Despite their isotopic similarities, there are considerable compositional differences with respect to major and trace element geochemistry both between and within the two islands that are due to short-length scale mantle heterogeneity beneath the region as well as variability in partial melting and magma storage conditions. Published geophysical surveys suggest that lithospheric thickness beneath both islands is ~60 km; this is consistent with the results of our geochemical modelling (59-60 km), which also suggest mantle potential temperatures between 1415-1435°C, similar to other documented continental passive rifts. Although there is some compositional overlap between the three synthems at Linosa, in general the older magmas (Arena Bianca, 700 ka) formed as a result of ~5% partial melting of a depleted MORB mantle (DMM) source enriched with a relatively small amount of recycled MORB material, which differentiated in a shallow-level (~8 km) magma chamber prior to eruption whereas the younger magmas (Monte Bandiera, 530 ka) formed as a result of ~2% partial melting of a similar mantle source, which differentiated in a magma chamber at or below the base of the crust (~25 km). Pantelleria magmas formed from a higher degree (~6%) of partial melting of a DMM source enriched with a relatively greater amount of recycled MORB material with possibly other components. Data for the seamounts are scarce and compromised by significant seawater alteration; thus, these volcanic centers cannot be modelled but based on comparative geochemistry with the islands are likely the result of even smaller (60 km) lithosphere. Magmas stored in the higher-level chamber were more effectively homogenized and preserve a narrower compositional range. Despite the geophysical similarities between the two islands in terms of lithospheric thickness and crustal thinning, melt productivity has been greater at Pantelleria, producing a much larger island and sustaining felsic magmatism, which may ultimately be entirely due to the local occurrence of much more fusible mantle

Clinopyroxene dissolution records rapid magma ascent

Magma ascent rates control volcanic eruption styles. However, the rates at which basaltic magmas ascend through the crust remain highly uncertain. Although recent studies have successfully exploited records of decompression driven degassing to estimate the rates at which H2O-rich basalts ascend, such approaches cannot readily be applied to primitive and H2O-poor basalts that erupt in ocean island and mid-ocean ridge settings. Here we present magma ascent rates obtained by modelling the dissolution of clinopyroxene crystals in a wehrlitic nodule from the primitive Borgarhraun lava flow in North Iceland. High-Al2O3 clinopyroxene core compositions are consistent with crystallisation near the Moho (_800 MPa), whereas low-Al2O3 clinopyroxene rims and inclusion compositions are consistent with crystallisation at or near the surface. We interpret low-Al2O3 rims and inclusions as the crystallised remnants of boundary layers formed by the dissolution of high-Al2O3 clinopyroxene during magma ascent. By combining characteristic rim dissolution lengths of 50–100 m with published experimental calibrations of clinopyroxene dissolution behaviour, we estimate that the Borgarhraun magma most likely decompressed and ascended at rates of 3.0–15 kPa.s-1 and 0.11–0.53 m.s-1 respectively. These rates are slightly faster than published estimates obtained by modelling the diffusive re-equilibration of olivine crystals, suggesting that the Borgarhraun magma either accelerated upwards or that it stalled briefly at depth prior to final ascent. Comparisons with other basaltic eruptions indicate that the H2O-poor magma that fed the dominantly effusive Borgarhraun eruption ascended at a similar rate to some H2O-rich magmas that have fed explosive eruptions in arc settings. Thus, magma ascent rates do not appear to correlate simply with magma H2O contents. Overall, our findings confirm that primitive and H2O-poor basalts can traverse the crust within days, and may erupt with little precursory warning of magma ascent

Melt inclusion constraints on petrogenesis of the 2014–2015 Holuhraun eruption, Iceland

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 H2O 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

Diffusive over-hydration of olivine-hosted melt inclusions

The pre-eruptive water content of magma is often estimated using crystal-hosted melt inclusions. However, olivine-hosted melt inclusions are prone to post-entrapment modification by H+ diffusion as they re-equilibrate with their external environment. This effect is well established for the case of H+ loss from olivine-hosted inclusions that have cooled slowly in degassed magma. Here we present evidence for the opposite effect: the addition of H+ into inclusions that are held in melts that are enriched in H2O with respect to the trapped melts. The compositional variability in a suite of 211 olivine-hosted inclusions from the Laki and Skuggafjöll eruptions in Iceland's Eastern Volcanic Zone indicates that diffusive H+ gain governs the H2O content of incompatible trace element depleted inclusions. Individual eruptive units contain olivine-hosted inclusions with widely varying incompatible element concentrations but near-constant H2O. Furthermore, over 40% of the inclusions have H2O/Ce>380H2O/Ce>380, significantly higher than the H2O/Ce expected in primary Icelandic melts or mid-ocean ridge basalts (150–280). The fact that the highest H2O/Ce ratios are found in the most incompatible element depleted inclusions indicates that hydration is a consequence of the concurrent mixing and crystallisation of compositionally diverse primary melts. Hydration occurs when olivines containing depleted inclusions with low H2O contents are juxtaposed against more hydrous melts during mixing. Melt inclusions from a single eruption may preserve evidence of both diffusive H+ loss and H+ gain. Trace element data are therefore vital for determining H2O contents of melt inclusions at the time of inclusion trapping and, ultimately, the H2O content of the mantle source regions

Mixing between chemically variable primitive basalts creates and modifies crystal cargoes

From Springer Nature via Jisc Publications RouterHistory: received 2021-01-13, accepted 2021-07-23, registration 2021-09-03, pub-electronic 2021-09-17, online 2021-09-17, collection 2021-12Publication status: PublishedFunder: Deutsche Forschungsgemeinschaft (German Research Foundation); doi: https://doi.org/10.13039/501100001659; Grant(s): NE2097/1-1Funder: University of Manchester; doi: https://doi.org/10.13039/501100000770Funder: RCUK | Natural Environment Research Council (NERC); doi: https://doi.org/10.13039/501100000270; Grant(s): NE/T011106/1Abstract: Basaltic crystal cargoes often preserve records of mantle-derived chemical variability that have been erased from their carrier liquids by magma mixing. However, the consequences of mixing between similarly primitive but otherwise chemically variable magmas remain poorly understood despite ubiquitous evidence of chemical variability in primary melt compositions and mixing-induced disequilibrium within erupted crystal cargoes. Here we report observations from magma–magma reaction experiments performed on analogues of primitive Icelandic lavas derived from distinct mantle sources to determine how their crystal cargoes respond to mixing-induced chemical disequilibrium. Chemical variability in our experimental products is controlled dominantly by major element diffusion in the melt that alters phase equilibria and triggers plagioclase resorption within regions that were initially plagioclase saturated. Isothermal mixing between chemically variable basaltic magmas may therefore play important but previously underappreciated roles in creating and modifying crystal cargoes by unlocking plagioclase-rich mushes and driving resorption, (re-)crystallisation and solid-state diffusion