Clinopyroxene–Liquid Equilibria and Geothermobarometry in Natural and Experimental Tholeiites: the 2014–2015 Holuhraun Eruption, Iceland

Abstract

Abstract Clinopyroxene–liquid geothermobarometry is a widely used tool for estimating the conditions under which mafic magmas are stored before they erupt. However, redox variability, sector zoning and disequilibrium crystallization present major challenges to the robust estimation of magma storage conditions. Moreover, most recent studies seeking to address these challenges have focused on clinopyroxenes from alkalic systems and are thus of limited use for understanding clinopyroxenes from the tholeiitic systems that dominate global magma budgets. Here we combine observations on natural clinopyroxenes from the 2014–2015 Holuhraun lava in Iceland with observations on experimental clinopyroxenes synthesized during high-pressure, high-temperature experiments on the same lava in order to investigate clinopyroxene–liquid equilibria in tholeiitic systems and optimize of geothermobarometric strategies. Natural clinopyroxenes from the 2014–2015 Holuhraun lava are sector zoned, with {1-11} hourglass sectors being enriched in the enstatite–ferrosillite component at the expense of all other components with respect to {hk0} prism sectors. In contrast with observations on clinopyroxenes from alkalic systems, sector zoning in clinopyroxenes from the 2014–2015 Holuhraun lava is characterized by differences in Ca and Na contents as well as in Ti and Al contents. The products of crystallization experiments performed at 100–600 MPa and 1140–1220 °C on a powdered starting glass at two sets of melt H2O content–oxygen fugacity conditions (∼0·1 wt % H2O and close to the graphite-oxygen redox buffer, and 0·5–1·0 wt % H2O and approximately one and half log units above the quartz–fayalite–magnetite redox buffer) demonstrate that clinopyroxene crystals from nominally equilibrium experiments can preserve strongly disequilibrium compositions. The compositional systematics of experimental clinopyroxenes are consistent with the presence of sector zoning. Furthermore, the magnitude of compositional variability increases with decreasing melt H2O content and increasing deviations of experimental temperatures below clinopyroxene liquidus temperatures (i.e. degrees of undercooling sensu lato), indicating that kinetic processes play a key role in controlling clinopyroxene compositions, even under notionally equilibrium conditions. Few published analyses of experimental clinopyroxene crystals may thus represent truly equilibrium compositions. Stoichiometric calculations on natural and experimental clinopyroxenes show that Fe3+ is a major constituent of clinopyroxenes from tholeiitic magmas under naturally relevant oxygen fugacity conditions. They also show that Fe3+ is most likely incorporated as Ca- and Al- bearing Ca–Fe-Tschermak’s component rather than Na-bearing aegirine component at oxygen fugacities up to one and a half log units above the quartz–fayalite–magnetite buffer. Elevated oxygen fugacities are thus less likely to compromise clinopyroxene–liquid geothermobarometry than previously thought. Guided by our experimental results, we combined published descriptions of clinopyroxene–liquid equilibria with geothermobarometric equations to develop an internally consistent and widely applicable method for performing geothermobarometry on tholeiitic magmas that does not require equilibrium zones to be selected a priori. Applying this method to natural clinopyroxene crystals from the 2014–2015 Holuhraun lava that formed under low but variable degrees of undercooling (perhaps 25 °C or less) returns values in excellent agreement with those from independent methods (232 ± 86 MPa, 1161 ± 11 °C). Robust estimates of magma storage conditions can thus be obtained by performing clinopyroxene–liquid geothermobarometry on tholeiitic magmas when disequilibrium is suitably accounted for.