Abstract
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.
1. Introduction
Volcanoes are underlain by vertically extensive and geometrically complex magmatic plumbing systems that are capable of dynamically reorganizing themselves over the timescales of individual eruptions (Marsh, 2004; Cashman et al., 2017; Sparks and Cashman, 2017; Magee et al., 2018; Maclennan, 2019; Sparks et al., 2019). It is therefore essential to quantify magma storage and transport timescales if we are to develop accurate models of volcanic behavior. Magma ascent rates are of particular importance in this respect because they play a central role in determining the style of volcanic eruptions, with eruptions fed by rapidly ascending magmas typically being more explosive than those fed by slowly ascending magmas (Gonnermann and Manga, 2007; Rutherford, 2008). Changes in magma ascent rate can also drive potentially hazardous transitions in eruption style (Woods and Koyaguchi, 1994). Importantly, global systematics in the eruptive behavior of volcanoes can now be resolved through the local complexity of individual events (Cassidy et al., 2018). However, most published estimates of magma ascent rate are from silicic systems despite the fact that basaltic systems are geographically widespread and dominate global magma budgets (Cassidy et al., 2018, and references therin). Although a number of studies have gone some way toward redressing this imbalance (e.g., Ruprecht and Plank, 2013; Lloyd et al., 2014; Peslier et al., 2015; Newcombe et al., 2020), estimates of magma ascent rate from basaltic systems remain few in number and somewhat uncertain in nature.
Magma ascent rates are commonly estimated using numerical modeling approaches of varying complexity (Wilson and Head, 1981; Mastin, 2002; La Spina et al., 2015, 2019). Analog modeling approaches have also been used to constrain magma ascent rates and inform the application of numerical models to natural systems (Namiki and Manga, 2006, 2008; Wright et al., 2012). However, ascent rates determined using modeling approaches are often subject to large uncertainties because conduit dimensions and eruptive mass fluxes are often very poorly constrained. Direct estimates of magma ascent rate have been obtained by modeling the diffusive re-equilibration of crystals entrained shortly before eruption (e.g., Ruprecht and Plank, 2013; Mollo et al., 2015; Mutch et al., 2019b; Newcombe et al., 2020), and the diffusive loss of volatiles from melt inclusions and melt embayments as a result of decompression-driven degassing (e.g., Humphreys et al., 2008; Lloyd et al., 2014; Hartley et al., 2018; Barth et al., 2019). Crystal textures have also been used to determine ascent rates in cases where crystallization kinetics are sufficiently well-understood (Armienti et al., 2013). However, ascent rates estimated from the re-equilibration of crystals, whether chemical or textural, convolve isobaric intervals of crystal entrainment with polybaric intervals of magma ascent, and are thus likely to underestimate true values. Although techniques that exploit decompression-driven volatile exsolution circumvent this particular problem, they can only be applied to magmas that degas significantly en route to the surface. Given that primitive basalts from ocean island and mid-ocean ridge settings often degas little during ascent because of their low volatile (i.e., H2O) contents (<0.2 wt.% H2O; Michael, 1988; Saal et al., 2002; Miller et al., 2019), the rates at which they ascend must be derived using different and volatile-independent approaches.
Basalts frequently carry crystals and nodules formed at depth to the surface, and exploiting the pressure-dependent stability of these crystalline cargoes provides a means of directly estimating the rates at which H2O-poor primitive basalts ascend (e.g., Brearley and Scarfe, 1986). Clinopyroxene is especially important in this regard because it commonly crystallizes from primitive basalts and its stability correlates strongly with pressure (Figure 1; Bender et al., 1978; Presnall et al., 1978; Fujii and Bougault, 1983; Chen and Zhang, 2009; Holland et al., 2018). Clinopyroxene crystals formed at depth may therefore dissolve during ascent, as evidenced by the resorbed nature of phenocrystic and xenocrystic clinopyroxenes erupted from some basaltic volcanoes (e.g., Tsuchiyama, 1986; Gurenko and Sobolev, 2006; Eason and Sinton, 2009). Indeed, the dissolution of clinopyroxene crystals during magma ascent may be responsible for the pyroxene paradox whereby many mid-ocean ridge basalts preserve geochemical signals of clinopyroxene crystallization but contain no clinopyroxene crystals (O'Donnell and Presnall, 1980; Thompson et al., 1980; Grove et al., 1992).
Figure 1
Here we exploit a quantitative but little used description of clinopyroxene stability in basaltic magmas to estimate the rates at which primitive Icelandic basalts ascend through the crust. We do this by modeling the dissolution of clinopyroxene crystals in a wehrlitic nodule from the Borgarhraun lava flow in North Iceland. Our findings suggest that H2O-poor basalts from ocean island and mid-ocean ridge systems can ascend as fast as many of their H2O-rich counterparts from arc systems. This has important implications for evaluating risks posed by basaltic volcanoes as primitive magmas appear capable of traversing the whole crust within a few days.
2. Nodules in Icelandic Basalts
Basaltic magmas often carry nodules of crystalline material to the surface. These nodules provide crucial insights into the chemical structure of Earth's interior and the dynamics of magmatic plumbing systems (e.g., Rudnick, 1992; Pearson et al., 2003; Holness et al., 2007; Cooper et al., 2016). Magmatic nodules occur in the products of numerous Icelandic eruptions and are variably wehrlitic, gabbroic, troctolitic and feldspathic in composition (Trønnes, 1990; Hansen and Grönvold, 2000; Maclennan et al., 2003a; Gurenko and Sobolev, 2006; Holness et al., 2007). Gabbroic, troctolitic and feldspathic nodules in lavas and tephras from the Eastern Rift Zone of Iceland are thought to constitute the disaggregated remnants of plagioclase-rich crystal mushes from the mid-crust (Hansen and Grönvold, 2000; Holness et al., 2007); macrocrysts and glomerocrysts of anorthitic plagioclase from the same samples are likewise considered to reflect more efficiently disaggregated portions of the same mushes (Halldórsson et al., 2008; Neave et al., 2014). In contrast, wehrlitic, gabbroic and troctolitic nodules from the Miðfell and Borgarhraun lava flows in the Western and Northern Volcanic Zones of Iceland, respectively, are thought to reflect the crystallization, mixing and reaction of primitive magmas in the lower crust and uppermost mantle (Trønnes, 1990; Maclennan et al., 2003a; Gurenko and Sobolev, 2006). Nodules from the Borgarhraun lava flow are of particular interest because of their extremely primitive compositions [olivine XFo up to 0.92, where XFo = molar Mg/(Mg+Fe); plagioclase XAn up to 0.92, where XAn = molar Ca/(Ca+Na); and clinopyroxene Mg#cpx up to 0.90, where Mg#cpx = Mg/(Mg+Fe)]. These crystals therefore preserve information about the earliest evolution of mantle-derived melts that would have otherwise been obscured by mixing and sampling bias (Maclennan, 2008; Winpenny and Maclennan, 2011; Neave et al., 2019b).
Nodules from the Borgarhraun lava flow are well-suited for estimating primitive basalt ascent rates because there are excellent geobarometric constraints on the pressure (i.e., depth) at which they were stored prior to eruption (Figure 1); the depth from which magma ascended is well-known. Clinopyroxene-liquid geobarometry using the model of Putirka et al. (1996) indicates that the Borgarhraun magma was stored at 810 ± 110(1σ) MPa in the uppermost mantle prior to eruption (Winpenny and Maclennan, 2011). The model uncertainty (standard error of estimate, SEE) associated with the geobarometer of Putirka et al. (1996) is 140 MPa. This pre-eruptive storage pressure is corroborated by experimental petrology and geobarometry exploiting the pressure dependence of olivine-plagioclase-augite-melt cotectic positions (Yang et al., 1996; Maclennan et al., 2003b, 2012). In contrast, clinopyroxene-liquid geobarometry using the model of Neave and Putirka (2017) suggests that the Borgarhraun magma was stored within the lower crust at a slightly lower pressure of 570 ± 120(1σ) MPa. This model is also associated with an SEE of 140 MPa. Despite this discrepancy in storage pressure estimates, which probably reflects the model of Neave and Putirka (2017) being optimized for lower pressures and more evolved compositions than those considered here, it is nevertheless clear that the Borgarhraun magma was stored near, and quite probably below, the Moho immediately before eruption.
A further advantage of investigating the Borgarhraun lava flow is that two recent studies into magma storage and ascent timescales provide an excellent framework within which we can interpret our new observations (Mutch et al., 2019a,b). Specifically, Mutch et al. (2019a) exploited the diffusive re-equilibration of spinel chadacrysts within clinopyroxene oikocrysts from wehrlitic nodules to estimate deep magma residence times on the order of 1,000 years. Exploiting the diffusive re-equilibration of olivine macrocrysts in the rims of wehrlitic nodules then allowed Mutch et al. (2019b) to estimate transcrustal magma transport times on the order of 10 days that correspond to minimum magma ascent rates of 0.02–0.1 m.s−1.
Here we investigate a wehrlitic nodule in sample JM27 collected in July 1999 from the eastern edge of the Borgarhraun lava flow (65°51.05′N, 16°59.93′W). Although nodules with variably dunitic, wehrlitic, gabbroic and troctolitic compositions have been described from the same location by Maclennan et al. (2003a), we focus on wehrlitic sample JM27 for three reasons: firstly, the abundance of clinopyroxene makes it possible to robustly constrain the minimum apparent dissolution lengths required to avoid dissolution time estimates being compromised by 3-D cutting effects; secondly, the presence of variably orientated crystals allows the effects of anisotropy on apparent dissolution lengths to be evaluated; and thirdly, the nodule's internal porosity likely prevented the dissipation of dissolution textures by convection or turbulent mingling. Nevertheless, we note that other nodules described by Maclennan et al. (2003a) show similar features to those we describe in sample JM27 (Supplementary Material), indicating that our findings are applicable to the Borgarhraun eruption as a whole.
3. Methods
QEMSCAN imaging of sample JM27 was performed on a FEI Quanta-650F instrument in the Department of Earth Sciences at the University of Cambridge, UK. Images were collected with a 4-μm pixel spacing, and 2000-count EDX spectra were used to produce phase and Al concentration maps following the principles discussed by Pirrie et al. (2004) and Neave et al. (2017a). Additional backscattered electron (BSE) imaging was performed on a FEI Quanta-650F instrument in the Department of Earth and Environmental Sciences at the University of Manchester, UK.
Mineral compositions were determined by electron probe microanalysis (EPMA) on a Cameca SX100 instrument in the Department of Earth Sciences at the University of Cambridge, UK. Silicon, Ti, Al, Cr, Fe, Mn, Mg, Ca, Na, K, and Ni were measured in minerals with an accelerating voltage of 15 kV. Olivine, clinopyroxene and spinel were measured with a current of 20 nA and a spot size of 1 μm. Plagioclase was measured with a current of 10 nA and a spot size of 5 μm. The following standards were used for calibration: diopside (Si and Ca), rutile (Ti), corundum (Al), Cr metal (Cr), fayalite (Fe), Mn metal (Mn), St. Johns olivine (Mg), jadeite (Na), orthoclase (K), and NiO (Ni). Peak and background counting times were typically 20 and 10 s, respectively for major elements and 40 and 20 s, respectively for minor elements, with the exception of Na that was counted on-peak for 10 s. Data quality was monitored by measuring the following secondary standards (representative analyses are provided in the Supplementary Material): San Carlos olivine (NMNH 111312-44), Lake County plagioclase (NMNH 115900), Kakanui augite (NMNH 122142), and Ney County Cr-augite (NMNH 164905) (Jarosewich et al., 1980, 1987). Accuracy and precision are both estimated as ≤ 1% relative for major elements (present at >1 wt.%) and ≤ 5% relative for minor elements (present at <1 wt.%).
4. Results
4.1. Petrography and QEMSCAN Imaging
The wehrlitic nodule in sample JM27 is composed of clinopyroxene with subordinate olivine and rare plagioclase (Figures 2A, 3A; Supplementary Material). Equant clinopyroxene crystals in the nodule are typically 500 μm to 1 mm in diameter, while equant olivine crystals are typically 200–500 μm in diameter. Rare tabular plagioclase crystals up to 250 μm in length also occur. The nodule hosts a network of microcrystalline pockets and embayments of groundmass that are contiguous with the host lava groundmass (Figure 3B). These pockets and channels are sometimes vesiculated, indicating that the nodule was porous to its host lava prior to lava flow emplacement. Both lava and nodule groundmasses are composed of high-aspect-ratio plagioclase crystals intergrown with small crystals of olivine, clinopyroxene and ilmenite (Figure 3B). The boundary between the wehrlitic nodule and the groundmass—defined by the limit of connected clinopyroxene crystals—is highly sinuous (Figure 2A). Crystals of olivine, plagioclase and clinopyroxene occur outside the nodule both as individual macrocrysts and as macrocrysts within monomineralic and polymineralic glomerocrysts (Figures 2A, 3D). Typical macrocrystic and glomerocrystic crystal diameters span the following ranges: 200 μm to 2 mm for olivine; 200–800 μm for plagioclase; and 200–800 μm for clinopyroxene. Glomerocrystic plagioclase and clinopyroxene crystals often occur in ophitic arrangements (Figure 3D).
Figure 2
Figure 3
QEMSCAN phase maps highlight abundant plagioclase inclusions within the cores of nodule clinopyroxene crystals that have not been described in previous studies (Figure 2B). These plagioclase inclusions are typically 20–100 μm in length and irregular to vermicular in form (Figure 3A). Many plagioclase inclusions are associated with small (10–50 μm) and irregular pockets of mesostasis that are texturally distinct from spheroidal melt inclusions (Figure 3C). Numerous small crystals of plagioclase also occur within the rims of nodule clinopyroxene crystals. Some of these crystals are separated from the groundmass by their host crystals and represent true inclusions. Others form continuous chains of crystals that extend from clinopyroxene rims deep into the interior of clinopyroxene cores. Within the limited compositional resolution of QEMSCAN phase identification, plagioclase inclusions have similar compositions to groundmass plagioclase crystals (bytownitic to labradoritic), while plagioclase macrocrysts and nodule plagioclase crystals are more primitive (anorthitic to bytownitic).
Semi-quantitative Al concentration maps indicate that plagioclase inclusions are located within domains of low-Al2O3 clinopyroxene that also appear to form inclusions within high-Al2O3 clinopyroxene cores (Figure 2C). Note that we refer to these low-Al2O3 clinopyroxene domains as inclusions without implying any genetic connotations. In other words, low-Al2O3 clinopyroxene inclusions are simply considered as domains of low-Al2O3 clinopyroxene that appear to be enclosed by high-Al2O3 clinopyroxene. Importantly, the presence and nature of these low-Al2O3 clinopyroxene inclusions was only revealed through Al concentration mapping; variations in clinopyroxene Al2O3 content were not apparent during optical microscopy or phase mapping. Nodule clinopyroxene rims have similarly low Al2O3 contents to clinopyroxene inclusions. Indeed, the chains of plagioclase crystals described in the preceding paragraph lie within low-Al2O3 channels that extend from clinopyroxene rims deep into their cores. Low-Al2O3 clinopyroxene inclusions are on the order of 50–250 μm in diameter, while the thinnest nodule clinopyroxene rims—those least affected by sectioning effects—are 50–100 μm wide.
High-Al2O3 clinopyroxene cores have relatively low BSE intensities while low-Al2O3 clinopyroxene rims and inclusions have relatively high BSE intensities (Figures 3A,B). Boundaries between nodule clinopyroxene cores and rims are often rounded (Figure 3B). As previously documented by Winpenny and Maclennan (2011), subtle variations in BSE intensity within nodule clinopyroxene cores delimit sector zones. Patchy variations in BSE intensity also occur at the edges of some clinopyroxene inclusions (Figure 3A), but appear to represent complex interfaces between inclusions and their host crystals rather than zoning intrinsic to low-Al2O3 clinopyroxene inclusions (Figures 3B,C). Importantly, we see no crystallographically controlled variations in BSE intensity that would reveal the presence of exsolution lamellae within nodule clinopyroxene cores (e.g., Holness et al., 2011).
Although BSE intensity is broadly constant within low-Al2O3 clinopyroxene inclusions, some inclusions contain films and pockets of very high BSE intensity (Figure 3C). These films and pockets typically occur next to plagioclase inclusions and along the inclusions' boundaries with high-Al2O3 clinopyroxene hosts, suggesting that they represent domains of late-stage crystallization. Further evidence for late-stage crystallization is provided by the small oxide crystals (5–20 μm in diameter, and subsequently identified as titanomagnetite) that occur along some interfaces between plagioclase inclusions and portions of clinopyroxene inclusions with very high BSE intensities (Figure 3C).
4.2. Electron Probe Microanalysis
Plagioclase and olivine compositions from sample JM27 are summarized in Figure 4. Plagioclase compositions range from XAn= 0.70 to XAn= 0.91 (Figure 4A), and different compositions are found in different textural associations. High XAn contents (>0.85) are only found in rare nodule plagioclase cores and the cores of some macrocrysts and glomerocrysts. Moderate XAn contents (0.78–0.85) are found in nodule plagioclase rims, some plagioclase inclusions and in the cores and rims of some macrocrysts and glomerocrysts. Low XAn contents (<0.78) are found in some plagioclase inclusions and the outermost rims of macrocrysts, glomerocrysts and nodule plagioclase crystals. These low XAn contents are typically located close to films and pockets of late-stage crystallization (Figure 3C). Olivine compositions range from XFo = 0.81 to XFo = 0.90 (Figure 4B). Relationships between composition and texture are less clear for olivine than plagioclase, though nodule, macrocryst, and glomerocryst olivine cores typically contain the highest XFo contents (>0.85). Lower XFo contents (0.82–0.85) are found in the cores and rims of some macrocrysts and glomerocrysts, and the lowest XFo content measured (0.81) is from a nodule olivine rim in contact with the groundmass.
Figure 4
Clinopyroxene compositions are summarized in Figures 5, 6, and range from Mg#cpx = 0.35 to Mg#cpx = 0.91, with the majority of compositions lying between Mg#cpx = 0.85 and Mg#cpx = 0.90. Different clinopyroxene compositions are closely linked with different textural associations. Nodule clinopyroxene cores are primitive (Mg#cpx = 0.88–0.91; Figures 5, 6) and rich in octahedral Al (Alvi = 0.08–0.11, where Alvi is calculated on a six-oxygen basis; Figure 6A). They are also somewhat enriched in jadeite component [XJd = 0.012–0.018, where XJd is calculated following Putirka (2008); Figure 6B]. Variability in TiO2 (0.11–0.20 wt.%; Figure 5B), Al2O3 (4.0–5.5 wt.%; Figure 5D) and wollastonite component [XWo = 0.41–0.44, where XWo = Ca/(Ca+Mg+Fe) on a molar basis; Figure 5F] over the range of Mg#cpx contents present in nodule clinopyroxene cores reflects sector zoning (Figure 3A). Some macrocryst and glomerocryst cores have similar compositions to nodule clinopyroxene cores (i.e., Mg#cpx = 0.88–0.90 and Alvi = 0.09–0.11), while others are slightly more evolved. Nodule clinopyroxene rims are consistently more evolved than nodule clinopyroxene cores (Mg#cpx = 0.72–0.88), a feature highlighted by the positive correlation between Mg#cpx and TiO2 defined by analyses from nodule clinopyroxene rims (Figure 5B). Rims are also considerably poorer in Alvi (<0.08; Figure 6A) and slightly poorer in XJd (<0.014; Figure 6B) than nodule clinopyroxene cores. Some macrocryst and glomerocryst compositions overlap with nodule clinopyroxene compositions, though a few analyses extend to much more evolved compositions (Mg#cpx reaches as low as 0.55 and TiO2 reaches as high as 1.0 wt.%), probably reflecting late-stage crystallization during lava flow emplacement.
Figure 5
Figure 6
Clinopyroxene inclusions are the most compositionally variable of all clinopyroxene textural types (Mg#cpx = 0.35–0.90). The majority of inclusion analyses fall within a compositional range shared with some macrocryst and glomerocryst analyses, as well as with analyses from nodule clinopyroxene rims (Mg#cpx = 0.85–0.87, Al2O3 = 2.0–3.0 wt.% and Alvi = 0.03–0.06; Figures 5B,D, 6A). These analyses are located within constant-BSE-intensity domains of clinopyroxene inclusions (Figure 3A). A primitive subpopulation of analyses from patchy zones at inclusion edges overlaps with analyses from nodule clinopyroxene core compositions (Mg#cpx = 0.88–0.90, Al2O3 = 4.5–5.0 wt.% and Alvi = 0.08–0.09), and probably constitutes analyses of complex, folded interfaces between inclusions and nodule clinopyroxene cores (Figure 3B). An evolved subpopulation of inclusion analyses overlaps partly with analyses from nodule clinopyroxene rims. These analyses are located within domains of very high BSE intensity that we have interpreted as pockets of late-stage crystallization (Figure 3C). The inflection of TiO2 content with decreasing Mg#cpx within this subpopulation reflects the onset of titanomagnetite crystallization.
5. Identifying Clinopyroxene Dissolution
High-Mg# clinopyroxene crystals in the Borgarhraun lava flow record high pre-eruptive magma storage pressures of 570–810 MPa (Winpenny and Maclennan, 2011; Neave and Putirka, 2017). These pressures are consistent with the high Al2O3, Alvi, and XJd contents of macrocryst, glomerocryst, and nodule clinopyroxene cores investigated here, confirming the ultimately deep origin of sample JM27 (Figures 3, 5B,C, 6; Aoki and Kushiro, 1968; Thompson, 1974; Blundy et al., 1995; Putirka et al., 1996). Given that the liquidus temperature of clinopyroxene is significantly higher at 570–810 MPa than 1 atm in Borgarhraun-like magmas [90–140°C higher according to calculations with the MELTS algorithm (Figure 1; Ghiorso and Sack, 1995; Smith and Asimow, 2005), and 70–120°C higher according to calculations on a compositionally analogous basalt from the Reykjanes Peninsula with THERMOCALC (RE46; Yang et al., 1996; Holland et al., 2018)], clinopyroxene crystals formed near the Moho are thus unlikely to have been in equilibrium with their carrier liquids upon eruption. Assuming that the ascending Borgarhraun magma cooled along an adiabatic gradient of <1°C.km−1 (Katz et al., 2003), it thus seems probable that deep-formed high-Al2O3 clinopyroxene crystals would have at least partly dissolved en route to the surface.
The rounded nature of high-Al2O3, high-Alvi, and high-XJd nodule clinopyroxene cores indicates that they underwent dissolution prior to their encapsulation within low-Al2O3, low-Alvi, and low-XJd rims at or near the surface (Figures 3B, 5B,C, 6; Aoki and Kushiro, 1968; Thompson, 1974; Blundy et al., 1995; Putirka et al., 1996). A low-pressure origin for clinopyroxene rims is indicated by comparisons with the products of 1-atm experiments performed by Yang et al. (1996), as well as by clinopyroxene-liquid geobarometry performed on nodule clinopyroxene rims using the model of Neave and Putirka (2017). Geobarometry was undertaken on clinopyroxene rims using the matrix glass composition of Sigurdsson et al. (2000), and returned a pressure of 210 ± 35(1σ) MPa with an associated model SEE of 140 MPa. All clinopyroxene-liquid pairs were checked for multicomponent equilibrium following the approach of Neave et al. (2019a). Although 210 MPa is appreciably greater than 1 atm, it is, model uncertainties notwithstanding, considerably lower than the 570–810 MPa range obtained from high-Al2O3 cores (Figure 1; Winpenny and Maclennan, 2011; Neave and Putirka, 2017). Moreover, that only a few analyses passed equilibrium filters (n = 3) suggests that clinopyroxene rims crystallized under disequilibrium conditions, feasibly resulting in an overestimation of crystallization pressures (Mollo et al., 2010).
The presence of abundant low-Al2O3 clinopyroxene inclusions within high-Al2O3 nodule clinopyroxene cores implies that dissolution was not only restricted to crystal faces, but also permeated crystal cores (Figure 3A). Similar textures have been described in clinopyroxene macrocrysts from Haleakala volcano in Hawaii by Welsch et al. (2016), but were attributed to rapid, kinetically controlled crystal growth rather than dissolution. These authors came to this interpretation because both low-Al2O3 and high-Al2O3 domains of clinopyroxene macrocrysts are in contact the groundmass in their samples, indicating that low- and high-Al2O3 clinopyroxene grew concurrently. However, only low-Al2O3 rims are in contact the groundmass in sample JM27, ruling out concurrent growth in our case. Low-Al2O3 clinopyroxene inclusions in the Borgarhraun lava flow also cut across boundaries between sector zones within high-Al2O3 nodule clinopyroxene cores, demonstrating that the inclusions formed by the resorption of pre-existing crystal cores. The occurrence of low-Al2O3 clinopyroxene embayments within high-Al2O3 cores also suggests that clinopyroxene inclusions may represent ramifying networks of low-Al2O3 clinopyroxene channels in 3-D rather than isolated inclusions of the kind visible in 2-D. Indeed, the compositional similarity between clinopyroxene inclusions and nodule clinopyroxene rims suggests that they share a common dissolution-reprecipitation origin that we discuss in detail below (Figure 5). Finally, domains of patchy zoning associated with complex interface geometries between inclusions and their host crystals are similar to cellular dissolution textures in other systems (Figure 3B; Streck, 2008), highlighting dissolution's central role in clinopyroxene inclusion genesis.
Our observations from the Borgarhraun lava flow are not unique; similar observations from lavas erupted elsewhere in Iceland validate our identification of substantial clinopyroxene dissolution. For example, Trønnes (1990) and Hansteen (1991) described resorbed “Al- and Cr-rich endiopside” crystals equivalent to our nodule clinopyroxene cores in primitive basalts from the Hengill volcanic system in the Western Rift Zone. Gurenko and Sobolev (2006) also describe spongy and embayed high-Al2O3 clinopyroxene crystals in glass-hosted nodules from the primitive Miðfell lava flow, also from the Hengill volcanic system, though these authors account for dissolution by melt-rock reaction rather than decompression during ascent. Clinopyroxene dissolution has also been proposed as an explanation for the chemical systematics of lavas from both the Eastern and Western Rift Zones (Halldórsson et al., 2008; Eason and Sinton, 2009), illustrating its potentially widespread role in basalt petrogenesis.
In contrast with the glassy nodules from the Miðfell lava flow described by Gurenko and Sobolev (2006), the wehrlitic nodule we describe here is fully crystalline; gabbroic and troctolitic nodules collected from the same location are similarly crystalline (Supplementary Material). Thus, while previously studied clinopyroxenes from the Miðfell lava flow are surrounded by glass, those we have studied from the Borgarhraun lava flow are surrounded by low-Al2O3 clinopyroxene rims. They are also permeated by low-Al2O3 clinopyroxene inclusions that host inclusions of relatively low-XAn plagioclase. This important textural difference—the presence or absence of low-Al2O3 clinopyroxene around and within high-Al2O3 clinopyroxene cores—reflects the contrasting emplacement histories of the Miðfell and Borgarhraun lava flows. Specifically, the Miðfell lava flow was erupted subglacially, making it possible for Gurenko and Sobolev (2006) to collect glassy material from rapidly quenched pillow basalt rims. In contrast, the Borgarhraun lava flow erupted subaerially, meaning that sample JM27 cooled slowly and crystallized in situ. The relatively slow cooling of the Borgarhraun lava flow did however allow important information about the extent of pre-eruptive high-Al2O3 dissolution to be captured. Namely, low-Al2O3 clinopyroxene rims and inclusions record the approximate dimensions of boundary layers formed by the dissolution of high-Al2O3 clinopyroxene during magma ascent which were then preserved during cooling by the high closure temperature of major element diffusion within clinopyroxene (Freer et al., 1982).
Fractional crystallization calculations performed on the mean composition of high-Al2O3 nodule clinopyroxene cores with the MELTS algorithm at 1 atm reproduce the phase assemblage and phase proportions observed in clinopyroxene rims and inclusions (Ghiorso and Sack, 1995; Smith and Asimow, 2005). Crystallization was assumed to be fractional because lava flow cooling rates probably outpaced mineral-mineral equilibration rates. Olivine crystallization was also suppressed in favor of epitactic clinopyroxene growth. Calculations were performed without imposing oxygen fugacity constraints; oxygen fugacity was imposed from the starting clinopyroxene Fe2+/Fe3+ estimated by charge balance following Lindsley and Andersen (1983). Once the calculated mass fraction of melt (F) reached 0.05, calculated mass fractions of clinopyroxene, plagioclase and spinel-structured oxide were 0.89, 0.05, and 0.01, respectively, broadly in line with BSE observations from clinopyroxene inclusions (Figures 3A,C). Melts with the mean composition of high-Al2O3 nodule clinopyroxene cores thus crystallize clinopyroxene and plagioclase in a ratio of ~18:1, a ratio much closer to that observed in clinopyroxene inclusions than the ~1:1 ratio calculated to result from matrix glass crystallization at 1 atm (Figures 3A,B). Moreover, pigeonite was predicted to crystallize at F <0.12, consistent with the low XWo content of low-Mg#cpx analyses from rims and inclusions associated with pockets of late-stage crystallization (Figures 3C, 5E).
Despite calculations with the MELTS algorithm reproducing the phase assemblage and phase proportions observed in clinopyroxene inclusions, they nevertheless return consistently higher Mg#cpx and XAn than observed in sample JM27 (Mg#cpx>0.9 until F <0.25 and XAn>0.8 until F <0.15). Although discrepancies between calculations with the MELTS algorithm and observations from both natural and experimental systems are well-documented (e.g., Neave et al., 2019b), they are probably exacerbated in this case by differences between the mean composition of high-Al2O3 nodule clinopyroxene cores and the compositions for which the MELTS algorithm has been calibrated (Ghiorso and Sack, 1995). However, a more fundamental explanation for these discrepancies stems from assuming that the boundary layers created by clinopyroxene dissolution were closed and compositionally homogeneous systems for their entire lifespans. Although these assumptions make performing illustrative calculations with the MELTS algorithm tractable, they are incorrect in detail. This is because our observations suggest Mg-rich and Na-poor high-Al2O3 clinopyroxene cannot be mass balanced against a Mg-poor and Na-rich assemblage containing low-Al2O3 clinopyroxene, plagioclase and titanomagnetite.
The dissolution rate of crystals is strongly controlled by the diffusivity of their constituent components in their surrounding melts, which results in the creation of concentration gradients around dissolving crystals (e.g., Zhang et al., 1989). Mean boundary layer compositions will therefore lie between those of dissolving crystals and their host melts. This may account for the slightly lower Mg# of most clinopyroxene inclusion and nodule clinopyroxene rim analyses (i.e., those not associated with pockets of late-stage crystallization) with respect to nodule clinopyroxene core analyses (Figure 5). Boundary layer compositions will also depend on the relative diffusivities of different elements, meaning that some exchange with the surrounding melt is highly likely to occur; even convectively stable boundary layers are unlikely to be fully isolated from the surrounding melt. Complex multicomponent effects notwithstanding (e.g., Liang et al., 1996), Na typically diffuses more quickly than other major species in silicate melts (Zhang et al., 2010b), providing a mechanism by which initially Na-poor boundary layers could have acquired sufficient Na to crystallize XAn ≤ 0.82 plagioclase. In contrast, Al diffuses more slowly than most other major species in silicate melts (Chen and Zhang, 2008; Yu et al., 2016), providing a mechanism by which Al-poor boundary layers formed by dissolution could have been preserved during magma ascent (cf. Neave et al., 2017b). However, evaluating the exact composition of boundary layers in sample JM27 would be challenging, and we suggest that such information would be best recovered from glassy nodules with feasibly intact boundary layers (e.g., Gurenko and Sobolev, 2006).
The narrowest nodule clinopyroxene rims we observed are typically 50–100 μm wide (Figures 3A,B). These minimum rim widths provide our best estimates of the distances that nodule clinopyroxene cores dissolved during ascent. This is because they are probably associated with crystal faces oriented normal to the plane of the thin section that are hence are minimally affected by 3-D cutting effects (e.g., Shea et al., 2015). Moreover, rim widths do not appear to vary substantially between different crystal faces, indicating that mineral anisotropy does not significantly affect clinopyroxene dissolution lengthscales in sample JM27. Clinopyroxene inclusions and embayments probably resulted from dissolution over greater distances than 100 μm, but may also reflect the exploitation of pre-existing fractures and cleavage planes (Chen and Zhang, 2009). It is also possible that the development of dissolution channels in 3-D may have affected the apparent dissolution lengths visible in 2-D. Importantly, the geometry of the wehrlitic nodule investigated here means that nodule clinopyroxene rims were largely protected from turbulence in the external magma during ascent. That is, we expect that nodule clinopyroxene rims record dissolution that was primarily diffusive rather than convective in nature (Kerr, 1995; Zhang and Xu, 2003). As a consequence, the thickness of nodule clinopyroxene rims encodes robust information about clinopyroxene dissolution times that can be used to estimate magma decompression and ascent rates.
6. Modeling Clinopyroxene Dissolution
Experimental investigations of clinopyroxene dissolution by Chen and Zhang (2009) provide a framework for modeling clinopyroxene dissolution recorded in the Borgarhraun lava flow. Specifically, Chen and Zhang (2009) characterized rates of diffusive diopside dissolution into a tholeiitic basalt by experimentally dissolving diopside crystals at a range of pressures and temperatures. They then used their experimental results to calibrate a model from which dissolution time could be predicted as a function of pressure, temperature, melt composition, clinopyroxene composition and dissolution length. We do however note the important limitation that Chen and Zhang (2009) did not incorporate any effects of mineral anisotropy into their simplified model. Nevertheless, coherence in their experimental results suggests that anisotropy exerts only second order effects over crystal dissolution rates. Importantly, dissolution experiments were carefully designed to avoid convection, making the resulting model well-suited for estimating dissolution times from the nodule clinopyroxene rims we describe here. Furthermore, the high-Mg#cpx compositions from the Borgarhraun lava flow are amongst the closest to pure diopside of any known compositions from Iceland (Maclennan et al., 2003a; Winpenny and Maclennan, 2011), making the Borgarhraun eruption an ideal candidate for estimating magma decompression and ascent rates using a model based on diopside dissolution. Moreover, clinopyroxene crystals in JM27 are free from exsolution lamellae that may complicate the interpretation of features created by diffusion and dissolution processes (e.g., Dohmen and Milke, 2010).
Chen and Zhang (2009) provided a series of five equations (their Equations 4, 5a–5d) that, in theory, can be solved iteratively to obtain the parameters required to quantitatively relate dissolution length and dissolution time in diopside-liquid systems. However, our attempts to iteratively solve these equations with the optim() function in the stats package of R often failed to converge (R Development Core Team, 2016), making it difficult to explore clinopyroxene stability across diverse pressure and temperature conditions. We therefore rearranged the equations of Chen and Zhang (2009) to reduce the number of unknown parameters from five to one, allowing the single remaining unknown parameter (the fitting parameter α) to be obtained through a more stable root-finding approach implemented with the uniroot() function in the stats package of R. These rearrangements are documented below.
As described by Chen and Zhang (2009) in their Equation (2a), chemical profiles in melts adjacent to dissolving crystals can be described by the 1-D diffusive dissolution equation (Crank, 1975):
where C is the concentration of a component in the melt, C0 is concentration in the melt at the crystal-melt interface and C∞ is the far-field concentration. x is the distance from the crystal-melt interface, D is the effective binary diffusivity of the element or oxide in question and t is duration of experiment. α is a fitting parameter related to the dissolution rate that, as per Equation (2b) of Chen and Zhang (2009), satisfies:
where Cc is the concentration of the element or oxide in question in the crystal.
As per Equation (2c) of Chen and Zhang (2009), the crystal diffusive dissolution distance can be expressed as:
where ρm and ρc are the densities of melt and crystal, respectively.
The stability of clinopyroxene crystals in basaltic melts is mediated, in part, by the diffusion of Mg and Ca cations within their host liquids (Chen and Zhang, 2009); the diffusion of Mg and Ca within clinopyroxene crystals is many orders of magnitude slower than within basaltic melts (Zhang et al., 2010a), and thus does not affect clinopyroxene stability over the timescales of interest here. As per Equations (3a) and (3b) of Chen and Zhang (2009), Mg and Ca diffusion within basaltic melts can be described by:
and:
where T is in K and D is in m2.s−1.
As per Equation (4) of Chen and Zhang (2009), the saturation condition at the clinopyroxene-melt interface can be expressed when:
where C0 values are expressed in oxide wt.%, P is pressure in GPa and a, b, e and g are fitting parameters with values of 22.85 ± 1.82(2σ), 2.10 ± 0.44(2σ), 26360 ± 2565(2σ), and 1.22 ± 0.43(2σ) × 106, respectively.
Rearranging Equation (6) and gathering the terms then allows us to define a function in which depends upon P and T alone:
and thus:
Given that L is equal for both MgO and CaO, Equation (3) leads to the expression:
If we define:
then rearranging Equation (2) provides:
Given that, as per Equation (8), , then:
and:
Rearranging these equations then gives:
and therefore:
A root finding approach can then be used to find values of αMg where the following expression is zero and substitutions are made using Equation (9):
and Equation (3) can then be used to calculate t.
The results of dissolution calculations performed assuming a fixed dissolution time of 1 s are presented in Figure 7 and demonstrate the strong effects of pressure and temperature on the stability of clinopyroxene with the mean composition of nodule clinopyroxene cores in the Borgarhraun matrix glass composition reported by Sigurdsson et al. (2000). Calculated clinopyroxene liquidus temperatures vary from ~1,190°C at 1 atm to ~1,320°C at 800 MPa and are encouragingly consistent with the results of calculations with the MELTS algorithm (Figure 1). These calculations also illustrate the strong dependence of dissolution rate on temperature and hence the vital importance of estimating magma eruption temperatures accurately when determining magma decompression and ascent rates from clinopyroxene dissolution lengths.
Figure 7
Indicative clinopyroxene dissolution times for the Borgarhraun lava flow were estimated by assuming that all dissolution took place at the surface following instantaneous magma ascent from near-Moho depths (Figure 8). A magma temperature of 1,230°C was taken from Mutch et al. (2019b) who performed calculations with the geothermometer of Sugawara (2000) that has an associated model SEE of 30°C. The temperature decrease associated with adiabatic cooling of the ascending magma over the pressure interval of relevance is well within the SEE of the geothermometer used (<8°C over 215 MPa; see below; Katz et al., 2003), and is thus not considered further. Uncertainties in dissolution time estimates were evaluated by resampling (n = 1,000) temperature and fitting parameters in Equation (6) for each dissolution length increment according to their respective uncertainties (Sugawara, 2000; Chen and Zhang, 2009). Resulting median dissolution times range from 30 min for a dissolution length of 50 μm to 150 min for a dissolution length of 100 μm—the range of minimum observed dissolution lengths (Figure 3). Lower 1σ bounds on dissolution times (i.e., 17th percentiles) are typically a factor of ten shorter than the median times and 1σ bounds (83rd percentiles) are typically a factor of 18 longer. Dissolution times are thus associated with 1σ uncertainties on the order of one order of magnitude which primarily reflect uncertainties in magma eruption temperatures that are challenging to overcome with currently available approaches (Putirka, 2008).
Figure 8
Mean magma decompression rates were estimated from clinopyroxene dissolution times by assuming that dissolution started as soon the ascending magma reached a pressure of 215 MPa, the pressure at which Equations (3) and (6) predict a change from clinopyroxene crystallization to clinopyroxene dissolution (Figure 7). Estimated mean decompression rates vary from 110 kPa.s−1 for a dissolution length of 50 μm to 24 kPa.s−1 for a dissolution length of 100 μm (Figure 9A). Estimated rates increase substantially to 210 MPa.s−1 for a dissolution length of 5 μm, and decrease modestly to 3.3 kPa.s−1 for a dissolution length of 300 μm. Mean magma ascent rates were estimated by converting the clinopyroxene saturation pressure of 215 MPa into a depth of 7.6 km by assuming a mean crustal density of 2.88 Mg.m−3 (Maclennan et al., 2001). Resulting mean ascent rate estimates vary from 3.9 m.s−1 for a dissolution length of 50 μm to 0.85 m.s−1 for a dissolution length of 100 μm (Figure 9B), and are subject to propagated uncertainties on order of one order of magnitude. However, assuming that all dissolution occurred at the surface following instantaneous ascent is likely to overestimate true magma decompression and ascent rates because it unrealistically implies that all dissolution occurred under the highest possible degrees of disequilibrium. In line with other studies (e.g., Lloyd et al., 2014; Mutch et al., 2019b; Newcombe et al., 2020), we therefore performed additional calculations with the more naturalistic assumption that magma ascended at a constant rate.
Figure 9
Assuming that clinopyroxene dissolution took place in response to decompression at a constant rate, dissolution lengths of 50 and 100 μm imply decompression rates of 15 and 3.0 kPa.s−1 respectively (Figure 9A). These values are approximately one order of magnitude lower than the mean decompression rates estimated by assuming that all dissolution took place following instantaneous ascent. Estimated magma ascent rates are also correspondingly lower: 0.53 and 0.11 m.s−1 for dissolution lengths of 50 and 100 μm respectively (Figure 9B). The difference between rates estimated with instantaneous and constant-rate assumptions primarily reflects the dependence of clinopyroxene dissolution rate on pressure (Figure 7). That is, clinopyroxene dissolution rates are much lower at high pressures where disequilibrium is negligible than at the surface where disequilibrium is greatest, meaning that very little clinopyroxene dissolves during the initial stages of clinopyroxene-undersaturated magma decompression.
Magma decompression rates estimated from clinopyroxene dissolution are compared with those estimated from the diffusive re-equilibration of primitive olivine crystals carried by the Borgarhraun magma in Figure 9C (Mutch et al., 2019b). Encouragingly, our best estimates (3.0–15 kPa.s−1) are broadly comparable with those obtained from a wholly independent diffusion chronometry approach for the same eruption (0.7–3.3 kPa.s−1). Our findings thus corroborate those of Mutch et al. (2019b) who estimated a mean magma ascent rate of 0.02–0.1 m.s−1, and confirm that transcrustal magma transport in Iceland may be considerably more rapid than reported in arc settings (0.5–1 mm.s−1; Ruprecht and Plank, 2013). On the basis of comparisons with transcrustal transport estimates from Mutch et al. (2019b), our best estimates also suggest that the Borgarhraun magma ascended more rapidly near the surface than it did at depth, consistent with either the upward acceleration of rising magma or the convolution of isobaric intervals of magma mixing with polybaric intervals of magma ascent by diffusion-based methods.
Magma decompression rates estimated for the Borgarhraun eruption are compared with published estimates from other basaltic systems in Figure 9C. Note that we compare decompression rather than ascent rates in order to avoid introducing additional uncertainties from converting pressures to depths across a range of different settings with different crustal density structures. Perhaps unsurprisingly, our best decompression rate estimates for the Borgarhraun eruption (3.0–15 kPa.s−1) are about two orders of magnitude lower than those reported for the sub-plinian 1974 eruption of Volcán de Fuego, Guatemala, based on the loss of volatiles from melt embayments and the loss of H2O from olivine crystals (240–600 and 380–450 kPa.s−1, respectively; Lloyd et al., 2014; Newcombe et al., 2020). Decompression rates estimated from the loss of volatiles from melt embayments for a series of Hawaiian to sub-plinian eruptions of Kı̄lauea, Hawaii, are also greater than those estimated for the Borgarhraun eruption (50–450 kPa.s−1), with the greatest Kı̄lauean rates being associated with the most explosive, sub-plinian eruptions (Ferguson et al., 2016). Indeed, decompression rates estimated for episode 1 of the Hawaiian 1959 Kı̄lauea Iki eruption from both the loss of volatiles from melt embayments and the loss of H2O from olivine crystals are only about one order of magnitude greater than our estimates for the Borgarhraun eruption (~50 kPa.s−1; Ferguson et al., 2016; Newcombe et al., 2020). Similar rates also been estimated for the Hawaiian 1977 eruption of Seguam, Alaska, using the same techniques (20–130 kPa.s−1 from embayments and 40–230 kPa.s−1 from olivine; Newcombe et al., 2020). Thus, magmas that have fed explosive basaltic eruptions ranging from Hawaiian to sub-plinian in style have typically decompressed more rapidly than the magma that fed the Borgarhraun eruption. However, this does not appear to have been the case for all explosively erupted basaltic magmas.
Decompression rates estimated for two moderately explosive eruptions of Cerro Negro, Nicaragua, in 1992 (VEI3, likely vulcanian) and 1995 (VEI2, likely Hawaiian) are broadly similar to those estimated for the Borgarhraun eruption, especially in the case of less explosive 1995 eruption (Barth et al., 2019; Newcombe et al., 2020). Specifically, rates estimated for the 1992 and 1995 eruptions from the loss of volatiles from melt embayments are 2.3–20 and 0.039–10 kPa.s−1 respectively (Barth et al., 2019), and those from the loss of H2O from olivine crystals are 7–12 and 3 kPa.s−1 respectively (Newcombe et al., 2020). Decompression rates comparable to those from the Borgarhraun eruption have also been reported for the vulcanian 2013 eruption of Chaparrastique, El Salvador, based on plagioclase dissolution textures (~3 kPa.s−1; Scarlato et al., 2017). Although deep magma decompression rates of ~16 and ~2.2 Pa.s−1 beneath Mount Etna and Campi Flegrei, Italy, estimated from clinopyroxene textures by Armienti et al. (2013) and Bonechi et al. (2020) respectively are about three orders of magnitude lower than our best estimates for the Borgarhraun eruption, shallow decompression rates estimated for Hawaiian eruptions of Mount Etna in 2011 on the basis of clinopyroxene-ilmenite re-equilibration are very similar to those we report here (1.3–13 kPa.s−1; Mollo et al., 2015). Finally, we note that our estimated decompression rates overlap with those obtained by modeling the diffusive loss of H2O from moderately evolved olivine-hosted melt inclusions in magmatic tephra collected during weakly Hawaiian phases of the 2014–2015 Holuhraun eruption, Iceland (3–8 kPa.s−1; Hartley et al., 2018). During its final ascent, the Borgarhraun magma therefore appears to have decompressed as rapidly as magmas that have fed various explosive basaltic eruptions ranging from Hawaiian to vulcanian in style.
The magma decompression rates that both we and Mutch et al. (2019b) estimate for the dominantly effusive and at most transiently explosive Borgarhraun eruption are broadly consistent with recently established relationships between explosivity and magma decompression rate (e.g., Barth et al., 2019). Namely, we find that the effusively erupted Borgarhraun magma decompressed more slowly than the magmas that fed basaltic sub-plinian eruptions at Fuego and Kı̄lauea. However, differences in magma decompression rate are modest or even absent when comparing the Borgarhraun eruption with Hawaiian and vulcanian eruptions at Cerro Negro, Chaparrastique and Mount Etna. Our calculations thus suggest that effusively erupted H2O-poor basalts from ocean island and mid-ocean ridge settings can decompress and ascend through the shallow crust just as rapidly as some explosively erupted H2O-rich basalts from arc settings. Moreover, it seems unlikely that the ascent rates of basaltic magmas are simply related to their H2O contents. Although H2O exsolution certainly plays a major role in driving the ascent of myriad magma types, other factors, such as melt buoyancy and CO2 exsolution must therefore drive the rapid ascent of primitive and H2O-poor basalts like the Borgarhraun magma.
Comparing rates of shallow, syn-eruptive magma decompression (including our estimate for the Borgarhraun eruption) with those for deep decompression associated with inter-eruptive magma assembly beneath Mount Etna and Campi Flegrei (Figure 9C; Armienti et al., 2013; Bonechi et al., 2020), suggests that magmas ascend much more rapidly during eruptions than during pre-eruptive magma recharge at depth. Magma transport leading to eruption may therefore operate in a fundamentally different way from magma transport out of the mantle and between magma reservoirs. It thus remains unclear how well seismic signals of inter-eruptive magma transport reflect likely rates of pre-eruptive magma ascent (White et al., 2011; Tarasewicz et al., 2012; Mutch et al., 2019b). Indeed, our results suggest that primitive and H2O-poor basalts may erupt with no more than a few days' precursory warning, making it challenging to evaluate the risks posed by some basaltic volcanoes.
7. Conclusions
Mineral textures and compositions in a wehrlitic nodule from the primitive Borgarhraun lava flow erupted in North Iceland record clinopyroxene dissolution during magma ascent from near-Moho depths. Mapping the distribution of Al distinguishes high-Al2O3 clinopyroxene cores from low-Al2O3 clinopyroxene rims and inclusions. The high Al2O3 content of clinopyroxene cores is consistent with crystallization at ~800 MPa, in line with published estimates of pre-eruptive magma storage pressures. In contrast, the low Al2O3 of clinopyroxene rims and inclusions is consistent with crystallization at or near the surface. Given that clinopyroxene rim and inclusion boundaries cut across sector zone boundaries within resorbed clinopyroxene cores, we infer that high-Al2O3 clinopyroxene cores dissolved during ascent as a consequence of the strong dependence of clinopyroxene stability on pressure. We thus interpret clinopyroxene rims and inclusions as recrystallized boundary layers protected from convective or turbulent dissipation by their host nodule. Corroborating evidence for clinopyroxene dissolution is provided by plagioclase and titanomagnetite crystals within clinopyroxene rims and inclusions that occur in proportions consistent with low-pressure crystallization of molten high-Al2O3 clinopyroxene cores during lava flow emplacement. Our findings are supported by textural and geochemical observations from across Iceland that indicate a potentially widespread role for clinopyroxene dissolution in basalt petrogenesis. It thus seems feasible that clinopyroxene dissolution could occur within any vertically extensive basaltic plumbing system and could hence offer a compelling explanation for the pyroxene paradox reported from mid-ocean ridge basalts.
Experimental calibrations of clinopyroxene stability suggest clinopyroxene began to dissolve once the Borgarhraun magma decompressed below the pressure of clinopyroxene saturation at 215 MPa. Our best decompression rate estimates of 3.0–15 kPa.s−1 were obtained from calculations assuming that decompression took place at a constant rate. These estimates could nonetheless be improved in future by quantifying the effects of mineral anisotropy on clinopyroxene dissolution kinetics with further experiments. Converting magma decompression rates into magma ascent rates results in values of 0.11–0.53 m.s−1 that are similar published estimates obtained from the diffusive re-equilibration of olivine (Mutch et al., 2019b). Nevertheless, our ascent rates are slightly faster than these published values, consistent with magma acceleration during ascent or the temporary stalling of magma after mush disaggregation. Although the Borgarhraun eruption broadly conforms to global correlations between eruption explosivities and magma decompression rates, the decompression rates we estimate are comparable to those obtained from some moderately explosive basaltic eruptions in arc settings, implying that there is no simple relationship between magma decompression rate and magma H2O content. Indeed, comparing magma ascent rates from Iceland with those from arc settings suggests that primitive and H2O-poor basalts can traverse the crust considerably faster than some H2O-rich basalts, though further work is required to evaluate how ascent rates are affected by the differing structures of mid-ocean ridge, ocean island and arc plumbing systems. It is, however, clear that magmas tapped from reservoirs near the Moho in ocean island and mid-ocean ridge settings may reach the surface within days. This represents a major challenge when evaluating the risks posed by infrequent but nonetheless potentially hazardous eruptions of primitive basalt.
Statements
Data availability statement
All datasets generated for this study are included in the article/Supplementary Material.
Author contributions
DN and JM co-designed the project. DN performed the BSE imaging and EPMA, implemented the clinopyroxene dissolution model, and wrote the manuscript. JM collected the sample, facilitated the QEMSCAN imaging, rearranged the equations of Chen and Zhang (2009), and contributed to writing the manuscript.
Funding
DN and this work were supported by a Presidential Fellowship from the University of Manchester.
Acknowledgments
We thank Lewis Hughes for assisting with BSE imaging at the University of Manchester and Iris Buisman for performing QEMSCAN imaging and assisting with EPMA at the University of Cambridge. We also thank CP and SM for their insightful and constructive reviews, and MP for his efficient editorial handling.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/feart.2020.00188/full#supplementary-material
Supplementary Data Sheet 1Raw BSE images.
Supplementary Table 1EMPA data.
Supplementary Figure 1JM21 4 μm QEMSCAN Al Map.
Supplementary Figure 2JM21 4 μm QEMSCAN Phase Map.
Supplementary Figure 3JM24 4 μm QEMSCAN Al Map.
Supplementary Figure 4JM24 4 μm QEMSCAN Phase Map.
Supplementary Figure 5JM27 4 μm QEMSCAN Al Map.
Supplementary Figure 6JM27 4 μm QEMSCAN Phase Map.
Supplementary Figure 7QEMSCAN legend.
References
1
AokiK.KushiroI. (1968). Some clinopyroxenes from ultramafic inclusions in Dreiser Weiher, Eifel. Contrib. Mineral. Petrol. 18, 326–337. 10.1007/BF00399694
2
ArmientiP.PerinelliC.PutirkaK. D. (2013). A new model to estimate deep-level magma ascent rates, with applications to Mt. Etna (Sicily, Italy). J. Petrol. 54, 795–813. 10.1093/petrology/egs085
3
BarthA.NewcombeM.PlankT.GonnermannH.HajimirzaS.SotoG. J.et al. (2019). Magma decompression rate correlates with explosivity at basaltic volcanoes–constraints from water diffusion in olivine. J. Volcanol. Geotherm. Res. 387:106664. 10.1016/j.jvolgeores.2019.106664
4
BenderJ. F.HodgesF. N.BenceA. E. (1978). Petrogenesis of basalts from the project FAMOUS area: experimental study from 0 to 15 kbars. Earth Planet. Sci. Lett. 41, 277–302. 10.1016/0012-821X(78)90184-X
5
BlundyJ. D.FalloonT. J.WoodB. J.DaltonJ. A. (1995). Sodium partitioning between clinopyroxene and silicate melts. J. Geophys. Res. Solid Earth100, 15501–15515. 10.1029/95JB00954
6
BonechiB.PerinelliC.GaetaM. (2020). Clinopyroxene growth rates at high pressure: constraints on magma recharge of the deep reservoir of the Campi Flegrei Volcanic District (south Italy). Bull. Volcanol. 82:5. 10.1007/s00445-019-1342-5
7
BrearleyM.ScarfeC. M. (1986). Dissolution rates of upper mantle minerals in an alkali basalt melt at high pressure: an experimental study and implications for ultramafic xenolith survival. J. Petrol. 27, 1157–1182. 10.1093/petrology/27.5.1157
8
CashmanK. V.SparksR. S. J.BlundyJ. D. (2017). Vertically extensive and unstable magmatic systems: a unified view of igneous processes. Science355:eaag3055. 10.1126/science.aag3055
9
CassidyM.MangaM.CashmanK. V.BachmannO. (2018). Controls on explosive-effusive volcanic eruption styles. Nat. Commun. 9:2839. 10.1038/s41467-018-05293-3
10
ChenY.ZhangY. (2008). Olivine dissolution in basaltic melt. Geochim. Cosmochim. Acta72, 4756–4777. 10.1016/j.gca.2008.07.014
11
ChenY.ZhangY. (2009). Clinopyroxene dissolution in basaltic melt. Geochim. Cosmochim. Acta73, 5730–5747. 10.1016/j.gca.2009.06.016
12
CooperG. F.DavidsonJ. P.BlundyJ. D. (2016). Plutonic xenoliths from Martinique, Lesser Antilles: evidence for open system processes and reactive melt flow in island arc crust. Contrib. Mineral. Petrol. 171:87. 10.1007/s00410-016-1299-8
13
CrankJ. (1975). The Mathematics of Diffusion. 2nd Edn. Oxford: Clarendon Press.
14
DohmenR.MilkeR. (2010). Diffusion in polycrystalline materials: grain boundaries, mathematical models, and experimental data. Rev. Mineral. Geochem. 72, 921–970. 10.2138/rmg.2010.72.21
15
EasonD. E.SintonJ. M. (2009). Lava shields and fissure eruptions of the Western Volcanic Zone, Iceland: evidence for magma chambers and crustal interaction. J. Volcanol. Geotherm. Res. 186, 331–348. 10.1016/j.jvolgeores.2009.06.009
16
FergusonD. J.GonnermannH. M.RuprechtP.PlankT. A.HauriE. H.HoughtonB. F.et al. (2016). Magma decompression rates during explosive eruptions of Klauea volcano, Hawaii, recorded by melt embayments. Bull. Volcanol. 78:71. 10.1007/s00445-016-1064-x
17
FreerR.CarpenterM. A.LongJ. V. P.ReedS. J. B. (1982). “Null result” diffusion experiments with diopside: implications for pyroxene equilibria. Earth Planet. Sci. Lett. 58, 285–292. 10.1016/0012-821X(82)90201-1
18
FujiiT.BougaultH. (1983). Melting relations of a magnesian abyssal tholeiite and the origin of MORBs. Earth Planet. Sci. Lett. 62, 283–295. 10.1016/0012-821X(83)90091-2
19
GhiorsoM. S.SackR. O. (1995). Chemical mass transfer in magmatic processes IV. A revised and internally consistent thermodynamic model for the interpolation and extrapolation of liquid-solid equilibria in magmatic systems at elevated temperatures and pressures. Contrib. Mineral. Petrol. 119, 197–212. 10.1007/s004100050036
20
GonnermannH. M.MangaM. (2007). The fluid mechanics inside a volcano. Annu. Rev. Fluid Mech. 39, 321–356. 10.1146/annurev.fluid.39.050905.110207
21
GroveT. L.KinzlerR. J.BryanW. B. (1992). Fractionation of mid-ocean ridge basalt (MORB), in Mantle Flow and Melt Generation at Mid-Ocean Ridges, Geophysical Monograph, Vol. 71, eds MorganJ. P.BlackmanD. K.SintonJ. M. (Washington DC: American Geophysical Union), 281–310. 10.1029/GM071p0281
22
GurenkoA. A.SobolevA. V. (2006). Crust-primitive magma interaction beneath neovolcanic rift zone of Iceland recorded in gabbro xenoliths from Midfell, SW Iceland. Contrib. Mineral. Petrol. 151, 495–520. 10.1007/s00410-006-0079-2
23
HalldórssonS. A.ÓskarssonN.GrönvoldK.SigurdssonG.SverrisdottirG.SteinthórssonS. (2008). Isotopic-heterogeneity of the Thjorsa lava-Implications for mantle sources and crustal processes within the Eastern Rift Zone, Iceland. Chem. Geol. 255, 305–316. 10.1016/j.chemgeo.2008.06.050
24
HansenH.GrönvoldK. (2000). Plagioclase ultraphyric basalts in Iceland: the mush of the rift. J. Volcanol. Geotherm. Res. 98, 1–32. 10.1016/S0377-0273(99)00189-4
25
HansteenT. H. (1991). Multi-stage evolution of the picritic Mælifell rocks, SW Iceland: constraints from mineralogy and inclusions of glass and fluid in olivine. Contrib. Mineral. Petrol. 109, 225–239. 10.1007/BF00306481
26
HartleyM. E.BaliE.MaclennanJ.NeaveD. A.HalldórssonS. A. (2018). Melt inclusion constraints on petrogenesis of the 2014–2015 Holuhraun eruption, Iceland. Contrib. Mineral. Petrol. 173:10. 10.1007/s00410-017-1435-0
27
HollandT. J. B.GreenE. C. R.PowellR. (2018). Melting of peridotites through to granites: a simple thermodynamic model in the system KNCFMASHTOCr. J. Petrol. 59, 881–900. 10.1093/petrology/egy048
28
HolnessM. B.AndersonA. T.MartinV. M.MaclennanJ.PassmoreE.SchwindingerK. R. (2007). Textures in partially solidified crystalline nodules: a window into the pore structure of slowly cooled mafic intrusions. J. Petrol. 48, 1243–1264. 10.1093/petrology/egm016
29
HolnessM. B.StrippG.HumphreysM. C. S.VekslerI. V.NielsenT. F. D.TegnerC. (2011). Silicate liquid immiscibility within the crystal mush: late-stage magmatic microstructures in the Skaergaard intrusion, East Greenland. J. Petrol. 52, 175–222. 10.1093/petrology/egq077
30
HumphreysM. C. S.MenandT.BlundyJ. D.KlimmK. (2008). Magma ascent rates in explosive eruptions: constraints from H2O diffusion in melt inclusions. Earth Planet. Sci. Lett. 270, 25–40. 10.1016/j.epsl.2008.02.041
31
JarosewichE.GooleyR.HuslerJ. (1987). Chromium augite–a new microprobe reference sample. Geostand. Geoanal. Res. 11, 197–198. 10.1111/j.1751-908X.1987.tb00027.x
32
JarosewichE.NelenJ. A.NorbergJ. A. (1980). Reference samples for electron microprobe analysis. Geostand. Newslett. 4, 43–47. 10.1111/j.1751-908X.1980.tb00273.x
33
KatzR. F.SpiegelmanM.LangmuirC. H. (2003). A new parameterization of hydrous mantle melting. Geochem. Geophys. Geosyst. 4, 1–19. 10.1029/2002GC000433
34
KerrR. C. (1995). Convective crystal dissolution. Contrib. Mineral. Petrol. 121, 237–246. 10.1007/BF02688239
35
La SpinaG.BurtonM.de' Michieli VitturiM. (2015). Temperature evolution during magma ascent in basaltic effusive eruptions: a numerical application to Stromboli volcano. Earth Planet. Sci. Lett. 426, 89–100. 10.1016/j.epsl.2015.06.015
36
La SpinaG.ClarkeA. B.de' Michieli VitturiM.BurtonM.AllisonC. M.RoggensackK.et al. (2019). Conduit dynamics of highly explosive basaltic eruptions: the 1085 CE sunset crater sub-Plinian events. J. Volcanol. Geotherm. Res. 387:106658. 10.1016/j.jvolgeores.2019.08.001
37
LiangY.RichterF. M.WatsonE. B. (1996). Diffusion in silicate melts: II. multicomponent diffusion in CaO-Al2O3-SiO2 at 1500°C and 1 GPa. Geochim. Cosmochim. Acta60, 5021–5035. 10.1016/S0016-7037(96)00352-3
38
LindsleyD. H.AndersenD. J. (1983). A two-pyroxene thermometer. J. Geophys. Res. 88, A887–A906. 10.1029/JB088iS02p0A887
39
LloydA. S.RuprechtP.HauriE. H.RoseW.GonnermannH. M.PlankT. A. (2014). NanoSIMS results from olivine-hosted melt embayments: magma ascent rate during explosive basaltic eruptions. J. Volcanol. Geotherm. Res. 283, 1–18. 10.1016/j.jvolgeores.2014.06.002
40
MaclennanJ. (2008). Concurrent mixing and cooling of melts under Iceland. J. Petrol. 49, 1931–1953. 10.1093/petrology/egn052
41
MaclennanJ. (2019). Mafic tiers and transient mushes: evidence from Iceland. Philos. Trans. R. Soc. A377:20180021. 10.1098/rsta.2018.0021
42
MaclennanJ.GaetaniG. A.HartleyM. E.NeaveD. A.WinpennyB. (2012). Petrological constraints on the crustal structure under rift zones, in AGU Fall Meeting Abstracts (San Francisco, CA), T41G.
43
MaclennanJ.McKenzieD.GrönvoldK.ShimizuN.EilerJ. M.KitchenN. (2003a). Melt mixing and crystallization under Theistareykir, northeast Iceland. Geochem. Geophys. Geosyst. 4, 1–40. 10.1029/2003GC000558
44
MaclennanJ.McKenzieD.GrönvoldK.SlaterL. (2001). Crustal accretion under Northern Iceland. Earth Planet. Sci. Lett. 191, 295–310. 10.1016/S0012-821X(01)00420-4
45
MaclennanJ.McKenzieD.HiltonF.GrönvoldK.ShimizuN. (2003b). Geochemical variability in a single flow from northern Iceland. J. Geophys. Res. 108, 1–21. 10.1029/2000JB000142
46
MageeC.StevensonC. T. E.EbmeierS. K.KeirD.HammondJ. O. S.GottsmannJ. H.et al. (2018). Magma plumbing systems: a geophysical perspective. J. Petrol. 59, 1217–1251. 10.1093/petrology/egy064
47
MarshB. D. (2004). A magmatic mush column Rosetta stone: the McMurdo dry valleys of Antarctica. EOS Trans. Am. Geophys. Union85, 497–508. 10.1029/2004EO470001
48
MastinL. G. (2002). Insights into volcanic conduit flow from an open-source numerical model. Geochem. Geophys. Geosyst. 3, 1–18. 10.1029/2001GC000192
49
MichaelP. J. (1988). The concentration, behavior and storage of H2O in the suboceanic upper mantle: implications for mantle metasomatism. Geochim. Cosmochim. Acta52, 555–566. 10.1016/0016-7037(88)90110-X
50
MillerW. G. R.MaclennanJ.ShorttleO.GaetaniG. A.Le RouxV.KleinF. (2019). Estimating the carbon content of the deep mantle with Icelandic melt inclusions. Earth Planet. Sci. Lett. 523:115699. 10.1016/j.epsl.2019.07.002
51
MolloS.Del GaudioP.VenturaG.IezziG.ScarlatoP. (2010). Dependence of clinopyroxene composition on cooling rate in basaltic magmas: implications for thermobarometry. Lithos118, 302–312. 10.1016/j.lithos.2010.05.006
52
MolloS.GiacomoniP. P.AndronicoD.ScarlatoP. (2015). Clinopyroxene and titanomagnetite cation redistributions at Mt. Etna volcano (Sicily, Italy): footprints of the final solidification history of lava fountains and lava flows. Chem. Geol. 406, 45–54. 10.1016/j.chemgeo.2015.04.017
53
MutchE. J. F.MaclennanJ.HollandT. J. B.BuismanI. (2019a). Millennial storage of near-Moho magma. Science264, 260–264. 10.1126/science.aax4092
54
MutchE. J. F.MaclennanJ.ShorttleO.EdmondsM.RudgeJ. F. (2019b). Rapid transcrustal magma movement under Iceland. Nat. Geosci. 12, 569–574. 10.1038/s41561-019-0376-9
55
NamikiA.MangaM. (2006). Influence of decompression rate on the expansion velocity and expansion style of bubbly fluids. J. Geophys. Res. Solid Earth111, 1–17. 10.1029/2005JB004132
56
NamikiA.MangaM. (2008). Transition between fragmentation and permeable outgassing of low viscosity magmas. J. Volcanol. Geotherm. Res. 169, 48–60. 10.1016/j.jvolgeores.2007.07.020
57
NeaveD. A.BaliE.GufinnssonG. H.HalldórssonS. A.KahlM.SchmidtA.-S.et al. (2019a). Clinopyroxene-liquid equilibria and geothermobarometry in natural and experimental tholeiites: the 2014–2015 Holuhraun eruption, Iceland. J. Petrol. 60, 1653–1680. 10.1093/petrology/egz042
58
NeaveD. A.BuismanI.MaclennanJ. (2017a). Continuous mush disaggregation during the long-lasting Laki fissure eruption, Iceland. Am. Mineral. 102, 2007–2021. 10.2138/am-2017-6015CCBY
59
NeaveD. A.HartleyM. E.MaclennanJ.EdmondsM.ThordarsonT. (2017b). Volatile and light lithophile elements in high-anorthite plagioclase-hosted melt inclusions from Iceland. Geochim. Cosmochim. Acta205, 100–118. 10.1016/j.gca.2017.02.009
60
NeaveD. A.MaclennanJ.HartleyM. E.EdmondsM.ThordarsonT. (2014). Crystal storage and transfer in basaltic systems: the Skuggafjöll eruption, Iceland. J. Petrol. 55, 2311–2346. 10.1093/petrology/egu058
61
NeaveD. A.NamurO.ShorttleO.HoltzF. (2019b). Magmatic evolution biases basaltic records of mantle chemistry towards melts from recycled sources. Earth Planet. Sci. Lett. 520, 199–211. 10.1016/j.epsl.2019.06.003
62
NeaveD. A.PutirkaK. D. (2017). A new clinopyroxene-liquid barometer, and implications for magma storage pressures under Icelandic rift zones. Am. Mineral. 102, 777–794. 10.2138/am-2017-5968
63
NewcombeM. E.PlankT.BarthA.AsimowP.HauriE. (2020). Water-in-olivine magma ascent chronometry: every crystal is a clock. J. Volcanol. Geotherm. Res. 398:106872. 10.1016/j.jvolgeores.2020.106872
64
O'DonnellT.PresnallD. C. (1980). Chemical variations of the glass and mineral phases in basalts dredged from 25–30 N along the Mid-Atlantic Ridge. Am. J. Sci. 280, 845–868.
65
PearsonD. G.CanilD.ShireyS. B. (2003). Mantle samples included in volcanic rocks: xenoliths and diamonds. Treat. Geochem. 2–9, 171–275. 10.1016/B0-08-043751-6/02005-3
66
PeslierA. H.BizimisM.MatneyM. (2015). Water disequilibrium in olivines from Hawaiian peridotites: recent metasomatism, H diffusion and magma ascent rates. Geochim. Cosmochim. Acta154, 98–117. 10.1016/j.gca.2015.01.030
67
PirrieD.ButcherA. R.PowerM. R.GottliebP.MillerG. L. (2004). Rapid quantitative mineral and phase analysis using automated scanning electron microscopy (QemSCAN); potential applications in forensic geoscience. Geol. Soc. Lond. Spec. Publ. 232, 123–136. 10.1144/GSL.SP.2004.232.01.12
68
PresnallD. C.DixonS. A.DixonJ. R.O'DonnellT. H.BrennerN. L.SchrockR. L.et al. (1978). Liquidus phase relations on the join diopside-forsterite-anorthite from 1 atm to 20 kbar: their bearing on the generation and crystallization of basaltic magma. Contrib. Mineral. Petrol. 66, 203–220. 10.1007/BF00372159
69
PutirkaK. D. (2008). Thermometers and barometers for volcanic systems. Rev. Mineral. Geochem. 69, 61–120. 10.2138/rmg.2008.69.3
70
PutirkaK. D.JohnsonM.KinzlerR. J.LonghiJ.WalkerD. (1996). Thermobarometry of mafic igneous rocks based on clinopyroxene-liquid equilibria, 0–30 kbar. Contrib. Mineral. Petrol. 123, 92–108. 10.1007/s004100050145
71
R Development Core Team (2016). R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing.
72
RudnickR. L. (1992). Xenoliths–samples of the lower continental crust. Contin. Lower Crust23, 269–316.
73
RuprechtP.PlankT. A. (2013). Feeding andesitic eruptions with a high-speed connection from the mantle. Nature500, 68–72. 10.1038/nature12342
74
RutherfordM. J. (2008). Magma ascent rates. Rev. Mineral. Geochem. 69, 241–271. 10.2138/rmg.2008.69.7
75
SaalA. E.HauriE. H.LangmuirC. H.PerfitM. R. (2002). Vapour undersaturation in primitive mid-ocean-ridge basalt and the volatile content of earth's upper mantle. Nature419, 451–455. 10.1038/nature01073
76
ScarlatoP.MolloS.Del BelloE.von QuadtA.BrownR. J.GutierrezE.et al. (2017). The 2013 eruption of Chaparrastique volcano (El Salvador): effects of magma storage, mixing, and decompression. Chem. Geol. 448, 110–122. 10.1016/j.chemgeo.2016.11.015
77
SheaT.LynnK. J.GarciaM. O. (2015). Cracking the olivine zoning code: distinguishing between crystal growth and diffusion. Geology43, 935–938. 10.1130/G37082.1
78
SigurdssonI. A.SteinthórssonS.GrönvoldK. (2000). Calcium-rich melt inclusions in Cr-spinels from Borgarhraun, northern Iceland. Earth Planet. Sci. Lett. 183, 15–26. 10.1016/S0012-821X(00)00269-7
79
SmithP. M.AsimowP. D. (2005). Adiabat-1ph: a new public front-end to the MELTS, pMELTS, and pHMELTS models. Geochem. Geophys. Geosyst. 6, 1–8. 10.1029/2004GC000816
80
SparksR. S. J.AnnenC.BlundyJ. D.CashmanK. V.RustA. C.JacksonM. D. (2019). Formation and dynamics of magma reservoirs. Philos. Trans. R. Soc. A377:20180019. 10.1098/rsta.2018.0019
81
SparksR. S. J.CashmanK. V. (2017). Dynamic magma systems: implications for forecasting volcanic activity. Elements13, 35–40. 10.2113/gselements.13.1.35
82
StaplesR. K.WhiteR. S.BrandsdóttirB.MenkeW.MaguireP. K. H.McBrideJ. H. (1997). Färoe-Iceland ridge experiment 1. Crustal structure of northeastern Iceland. J. Geophys. Res. Solid Earth102, 7849–7866. 10.1029/96JB03911
83
StreckM. J. (2008). Mineral textures and zoning as evidence for open system processes. Rev. Mineral. Geochem. 69, 595–622. 10.2138/rmg.2008.69.15
84
SugawaraT. (2000). Empirical relationships between temperature, pressure, and MgO content in olivine and pyroxene saturated liquid. J. Geophys. Res. 105, 8457–8472. 10.1029/2000JB900010
85
TarasewiczJ.BrandsdóttirB.WhiteR. S.HenschM.ThorbjarnardóttirB. (2012). Using microearthquakes to track repeated magma intrusions beneath the Eyjafjallajkull stratovolcano, Iceland. J. Geophys. Res. Solid Earth117, 1–13. 10.1029/2011JB008751
86
ThompsonG.BryanW. B.MelsonW. G. (1980). Geological and geophysical investigation of the mid-Cayman rise spreading center: geochemical variation and petrogenesis of basalt glasses. J. Geol. 88, 41–55. 10.1086/628472
87
ThompsonR. N. (1974). Some high-pressure pyroxenes. Mineral. Mag. 39, 768–787. 10.1180/minmag.1974.039.307.04
88
TrønnesR. G. (1990). Basaltic melt evolution of the Hengill volcanic system, SW Iceland, and evidence for clinopyroxene assimilation in primitive tholeiitic magmas. J. Geophys. Res. Solid Earth95, 15893–15910. 10.1029/JB095iB10p15893
89
TsuchiyamaA. (1986). Melting and dissolution kinetics: application to partial melting and dissolution of xenoliths. J. Geophys. Res. 91:9395. 10.1029/JB091iB09p09395
90
WelschB.HammerJ. E.BaronnetA.JacobS.HellebrandE.SintonJ. M. (2016). Clinopyroxene in postshield Haleakala ankaramite 2. Texture, compositional zoning, and supersaturation in the magma. Contrib. Mineral. Petrol. 171:6. 10.1007/s00410-015-1213-9
91
WhiteR. S.DrewJ.MartensH. R.KeyJ.SoosaluH.JakobsdóttirS. S. (2011). Dynamics of dyke intrusion in the mid-crust of Iceland. Earth Planet. Sci. Lett. 304, 300–312. 10.1016/j.epsl.2011.02.038
92
WilsonL.HeadJ. W. (1981). Ascent and eruption of basaltic magma on the earth and moon. J. Geophys. Res. Solid Earth86, 2971–3001. 10.1029/JB086iB04p02971
93
WinpennyB.MaclennanJ. (2011). A partial record of mixing of mantle melts preserved in Icelandic phenocrysts. J. Petrol. 52, 1791–1812. 10.1093/petrology/egr031
94
WoodsA. W.KoyaguchiT. (1994). Transitions between explosive and effusive eruptions of silicic magmas. Nature370, 641–644. 10.1038/370641a0
95
WrightH. M. N.CashmanK. V.MothesP. A.HallM. L.RuizA. G.Le PennecJ. L. (2012). Estimating rates of decompression from textures of erupted ash particles produced by 1999–2006 eruptions of Tungurahua volcano, Ecuador. Geology40, 619–622. 10.1130/G32948.1
96
YangH.-J.KinzlerR. J.GroveT. L. (1996). Experiments and models of anhydrous, basaltic olivine-plagioclase-augite saturated melts from 0.001 to 10 kbar. Contrib. Mineral. Petrol. 124, 1–18. 10.1007/s004100050169
97
YuY.ZhangY.ChenY.XuZ. (2016). Kinetics of anorthite dissolution in basaltic melt. Geochim. Cosmochim. Acta179, 257–274. 10.1016/j.gca.2016.02.002
98
ZhangX.GangulyJ.ItoM. (2010a). Ca-Mg diffusion in diopside: tracer and chemical inter-diffusion coefficients. Contrib. Mineral. Petrol. 159, 175–186. 10.1007/s00410-009-0422-5
99
ZhangY.NiH.ChenY. (2010b). Diffusion data in silicate melts. Rev. Mineral. Geochem. 72, 311–408. 10.2138/rmg.2010.72.8
100
ZhangY.WalkerD.LesherC. E. (1989). Diffusive crystal dissolution. Contrib. Mineral. Petrol. 102, 492–513. 10.1007/BF00371090
101
ZhangY.XuZ. (2003). Kinetics of convective crystal dissolution and melting, with applications to methane hydrate dissolution and dissociation in seawater. Earth Planet. Sci. Lett. 213, 133–148. 10.1016/S0012-821X(03)00297-8
Summary
Keywords
magma ascent rate, plumbing system, clinopyroxene dissolution, basalt, disequilibrium, Iceland
Citation
Neave DA and Maclennan J (2020) Clinopyroxene Dissolution Records Rapid Magma Ascent. Front. Earth Sci. 8:188. doi: 10.3389/feart.2020.00188
Received
09 March 2020
Accepted
11 May 2020
Published
10 June 2020
Volume
8 - 2020
Edited by
Mattia Pistone, University of Georgia, United States
Reviewed by
Silvio Mollo, Sapienza University of Rome, Italy; Chiara Maria Petrone, Natural History Museum, United Kingdom
Updates
Copyright
© 2020 Neave and Maclennan.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: David A. Neave david.neave@manchester.ac.uk
This article was submitted to Petrology, a section of the journal Frontiers in Earth Science
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.