The global pattern of trace-element distributions in ocean floor basalts

High Mg# of the continental crust explained by calc-alkaline differentiation

We used compiled geochemical data to investigate the mechanisms that control Mg# (molar ratio of Mg/(Mg + Fe_T)) in andesitic arc lavas. We find that andesites from mature continental arcs with crustal thickness of >45 km have systematically higher Mg# than those from oceanic arcs with crustal thickness of <30 km. The elevated Mg# in continental arc lavas results from strong Fe depletion during high-pressure differentiation favored in thick crusts. This proposal is reinforced by our compiled melting/crystallization experiment data. We show that the Mg# characteristics of continental arc lavas match that of the continental crust. These findings suggest that the formation of many high-Mg# andesites and the continental crust may not require slab-melt/peridotite interactions. Instead, the high Mg# of the continental crust can be explained by intracrustal calc-alkaline differentiation processes in magmatic orogens.

Basalt derived from highly refractory mantle sources during early Izu-Bonin-Mariana arc development

The magmatic character of early subduction zone and arc development is unlike mature systems. Low-Ti-K tholeiitic basalts and boninites dominate the early Izu-Bonin-Mariana (IBM) system. Basalts recovered from the Amami Sankaku Basin (ASB), underlying and located west of the IBM’s oldest remnant arc, erupted at ~49 Ma. This was 3 million years after subduction inception (51-52 Ma) represented by forearc basalt (FAB), at the tipping point between FAB-boninite and typical arc magmatism. We show ASB basalts are low-Ti-K, aluminous spinel-bearing tholeiites, distinct compared to mid-ocean ridge (MOR), backarc basin, island arc or ocean island basalts. Their upper mantle source was hot, reduced, refractory peridotite, indicating prior melt extraction. ASB basalts transferred rapidly from pressures (~0.7-2 GPa) at the plagioclase-spinel peridotite facies boundary to the surface. Vestiges of a polybaric-polythermal mineralogy are preserved in this basalt, and were not obliterated during persistent recharge-mix-tap-fractionate regimes typical of MOR or mature arcs.

Magmatism associated with early growth of subduction zones is unlike that of mature island arc systems. Here, the authors find basalts with distinct mineralogical and geochemical characteristics were erupted during this early stage, and derived from extremely refractory, hot mantle sources.

Iron isotopes trace primordial magma ocean cumulates melting in Earth's upper mantle

The differentiation of Earth ~4.5 billion years (Ga) ago is believed to have culminated in magma ocean crystallization, crystal-liquid separation, and the formation of mineralogically distinct mantle reservoirs. However, the magma ocean model remains difficult to validate because of the scarcity of geochemical tracers of lower mantle mineralogy. The Fe isotope compositions (δ⁵⁷Fe) of ancient mafic rocks can be used to reconstruct the mineralogy of their mantle source regions. We present Fe isotope data for 3.7-Ga metabasalts from the Isua Supracrustal Belt (Greenland). The δ⁵⁷Fe signatures of these samples extend to values elevated relative to modern equivalents and define strong correlations with fluid-immobile trace elements and tungsten isotope anomalies (μ¹⁸²W). Phase equilibria models demonstrate that these features can be explained by melting of a magma ocean cumulate component in the upper mantle. Similar processes may operate today, as evidenced by the δ⁵⁷Fe and μ¹⁸²W heterogeneity of modern oceanic basalts.

Competing effects of spreading rate, crystal fractionation and source variability on Fe isotope systematics in mid-ocean ridge lavas

Two-thirds of the Earth is covered by mid-ocean ridge basalts, which form along a network of divergent plate margins. Basalts along these margins display a chemical diversity, which is consequent to a complex interplay of partial mantle melting in the upper mantle and magmatic differentiation processes in lower crustal levels. Igneous differentiation (crystal fractionation, partial melting) and source heterogeneity, in general, are key drivers creating variable chemistry in mid-ocean ridge basalts. This variability is reflected in iron isotope systematics (expressed as δ⁵⁷Fe), showing a total range of 0.2 ‰ from δ⁵⁷Fe =  + 0.05 to + 0.25 ‰. Respective contributions of source heterogeneity and magma differentiation leading to this diversity, however, remain elusive. This study investigates the iron isotope systematics in basalts from the ultraslow spreading Gakkel Ridge in the Arctic Ocean and compares them to existing data from the fast spreading East Pacific Rise ridge. Results indicate that Gakkel lavas are driven to heavier iron isotope compositions through partial melting processes, whereas effects of igneous differentiation are minor. This is in stark contrast to fast spreading ridges showing reversed effects of near negligible partial melting effects followed by large isotope fractionation along the liquid line of descent. Gakkel lavas further reveal mantle heterogeneity that is superimposed on the igneous differentiation effects, showing that upper mantle Fe isotope heterogeneity can be transmitted into erupting basalts in the absence of homogenisation processes in sub-oceanic magma chambers.

An Early Cretaceous subduction-modified mantle underneath the ultraslow spreading Gakkel Ridge, Arctic Ocean

Earth's upper mantle, as sampled by mid-ocean ridge basalts (MORBs) at oceanic spreading centers, has developed chemical and isotopic heterogeneity over billions of years through focused melt extraction and re-enrichment by recycled crustal components. Chemical and isotopic heterogeneity of MORB is dwarfed by the large compositional spectrum of lavas at convergent margins, identifying subduction zones as the major site for crustal recycling into and modification of the mantle. The fate of subduction-modified mantle and if this heterogeneity transmits into MORB chemistry remains elusive. Here, we investigate the origin of upper mantle chemical heterogeneity underneath the Western Gakkel Ridge region in the Arctic Ocean through MORB geochemistry and tectonic plate reconstruction. We find that seafloor lavas from the Western Gakkel Ridge region mirror geochemical signatures of an Early Cretaceous, paleo-subduction zone, and conclude that the upper mantle can preserve a long-lived, stationary geochemical memory of past geodynamic processes.

Consequences of a crystal mush-dominated magma plumbing system: a mid-ocean ridge perspective

Crystal mush is rapidly emerging as a new paradigm for the evolution of igneous systems. Mid-ocean ridges provide a unique opportunity to study mush processes: geophysical data indicate that, even at the most magmatically robust fast-spreading ridges, the magma plumbing system typically comprises crystal mush. In this paper, we describe some of the consequences of crystal mush for the evolution of the mid-ocean ridge magmatic system. One of these is that melt migration by porous flow plays an important role, in addition to rapid, channelized flow. Facilitated by both buoyancy and (deformation-enhanced) compaction, porous flow leads to reactions between the mush and migrating melts. Reactions between melt and the surrounding crystal framework are also likely to occur upon emplacement of primitive melts into the mush. Furthermore, replenishment facilitates mixing between the replenishing melt and interstitial melts of the mush. Hence, crystal mushes facilitate reaction and mixing, which leads to significant homogenization, and which may account for the geochemical systematics of mid-ocean ridge basalt (MORB). A second consequence is cryptic fractionation. At mid-ocean ridges, a plagioclase framework may already have formed when clinopyroxene saturates. As a result, clinopyroxene phenocrysts are rare, despite the fact that the vast majority of MORB records clinopyroxene fractionation. Hence, melts extracted from crystal mush may show a cryptic fractionation signature. Another consequence of a mush-dominated plumbing system is that channelized flow of melts through the crystal mush leads to the occurrence of vertical magmatic fabrics in oceanic gabbros, as well as the entrainment of diverse populations of phenocrysts. Overall, we conclude that the occurrence of crystal mush has a number of fundamental implications for the behaviour and evolution of magmatic systems, and that mid-ocean ridges can serve as a useful template for trans-crustal mush columns elsewhere. This article is part of the Theo Murphy meeting issue 'Magma reservoir architecture and dynamics'.

Sulfide resorption during crustal ascent and degassing of oceanic plateau basalts

Mantle plume-related magmas typically have higher chalcophile and siderophile element (CSE) contents than mid-ocean ridge basalts (MORB). These differences are often attributed to sulfide-under-saturation of plume-related melts. However, because of eruption-related degassing of sulfur (S) and the compositional, pressure, temperature and redox effects on S-solubility, understanding the magmatic behavior of S is challenging. Using CSE data for oceanic plateau basalts (OPB), which rarely degas S, we show that many OPB are sulfide-saturated. Differences in the timing of sulfide-saturation between individual OPB suites can be explained by pressure effects on sulfur solubility associated with ascent through over-thickened crust. Importantly, where S-degassing does occur, OPB have higher CSE contents than S-undegassed melts at similar stages of differentiation. This can be explained by resorption of earlier-formed sulfides, which might play an important role in enriching degassed melts in sulfide-compatible CSE and potentially contributes to anomalous enrichments of CSE in the crust.

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.

Growing magma chambers control the distribution of small-scale flood basalts

Small-scale continental flood basalts are a global phenomenon characterized by regular spatio-temporal distributions. However, no genetic mechanism has been proposed to explain the visible but overlooked distribution patterns of these continental basaltic volcanism. Here we present a case study from eastern China, combining major and trace element analyses with Ar-Ar and K-Ar dating to show that the spatio-temporal distribution of small-scale flood basalts is controlled by the growth of long-lived magma chambers. Evolved basalts (SiO2 > 47.5 wt.%) from Xinchang-Shengzhou, a small-scale Cenozoic flood basalt field in Zhejiang province, eastern China, show a northward younging trend over the period 9.4-3.0 Ma. With northward migration, the magmas evolved only slightly ((Na2O + K2O)/MgO = 0.40-0.66; TiO2/MgO = 0.23-0.35) during about 6 Myr (9.4-3.3 Ma). When the flood basalts reached the northern end of the province, the magmas evolved rapidly (3.3-3.0 Ma) through a broad range of compositions ((Na2O + K2O)/MgO = 0.60-1.28; TiO2/MgO = 0.30-0.57). The distribution and two-stage compositional evolution of the migrating flood basalts record continuous magma replenishment that buffered against magmatic evolution and induced magma chamber growth. Our results demonstrate that the magma replenishment-magma chamber growth model explains the spatio-temporal distribution of small-scale flood basalts.

Primitive layered gabbros from fast-spreading lower oceanic crust

Three-quarters of the oceanic crust formed at fast-spreading ridges is composed of plutonic rocks whose mineral assemblages, textures and compositions record the history of melt transport and crystallization between the mantle and the sea floor. Despite the importance of these rocks, sampling them in situ is extremely challenging owing to the overlying dykes and lavas. This means that models for understanding the formation of the lower crust are based largely on geophysical studies and ancient analogues (ophiolites) that did not form at typical mid-ocean ridges. Here we describe cored intervals of primitive, modally layered gabbroic rocks from the lower plutonic crust formed at a fast-spreading ridge, sampled by the Integrated Ocean Drilling Program at the Hess Deep rift. Centimetre-scale, modally layered rocks, some of which have a strong layering-parallel foliation, confirm a long-held belief that such rocks are a key constituent of the lower oceanic crust formed at fast-spreading ridges. Geochemical analysis of these primitive lower plutonic rocks--in combination with previous geochemical data for shallow-level plutonic rocks, sheeted dykes and lavas--provides the most completely constrained estimate of the bulk composition of fast-spreading oceanic crust so far. Simple crystallization models using this bulk crustal composition as the parental melt accurately predict the bulk composition of both the lavas and the plutonic rocks. However, the recovered plutonic rocks show early crystallization of orthopyroxene, which is not predicted by current models of melt extraction from the mantle and mid-ocean-ridge basalt differentiation. The simplest explanation of this observation is that compositionally diverse melts are extracted from the mantle and partly crystallize before mixing to produce the more homogeneous magmas that erupt.