Chemical differentiation, cold storage and remobilization of magma in the Earth's crust.

The formation, storage and chemical differentiation of magma in the Earth's crust is of fundamental importance in igneous geology and volcanology. Recent data are challenging the high-melt-fraction 'magma chamber' paradigm that has underpinned models of crustal magmatism for over a century, suggesting instead that magma is normally stored in low-melt-fraction 'mush reservoirs'1-9. A mush reservoir comprises a porous and permeable framework of closely packed crystals with melt present in the pore space1,10. However, many common features of crustal magmatism have not yet been explained by either the 'chamber' or 'mush reservoir' concepts1,11. Here we show that reactive melt flow is a critical, but hitherto neglected, process in crustal mush reservoirs, caused by buoyant melt percolating upwards through, and reacting with, the crystals10. Reactive melt flow in mush reservoirs produces the low-crystallinity, chemically differentiated (silicic) magmas that ascend to form shallower intrusions or erupt to the surface11-13. These magmas can host much older crystals, stored at low and even sub-solidus temperatures, consistent with crystal chemistry data6-9. Changes in local bulk composition caused by reactive melt flow, rather than large increases in temperature, produce the rapid increase in melt fraction that remobilizes these cool- or cold-stored crystals. Reactive flow can also produce bimodality in magma compositions sourced from mid- to lower-crustal reservoirs14,15. Trace-element profiles generated by reactive flow are similar to those observed in a well studied reservoir now exposed at the surface16. We propose that magma storage and differentiation primarily occurs by reactive melt flow in long-lived mush reservoirs, rather than by the commonly invoked process of fractional crystallization in magma chambers14.

Trace element partitioning in basaltic systems as a function of oxygen fugacity.

Along with temperature, pressure and melt chemistry, magmatic oxygen fugacity (fO2) has an important influence on liquid and solid differentiation trends and melt structure. To explore the effect of redox conditions on mineral stability and mineral-melt partitioning in basaltic systems we performed equilibrium, one-atmosphere experiments on a picrite at 1200-1110 °C with fO2 ranging from NNO-4 log units to air. Clinopyroxene crystallizes from 1180 °C to near-solidus, along with plagioclase, olivine and spinel. Olivine Mg# increases with increasing fO2, eventually reacting to pigeonite. Spinel is absent under strongly reducing conditions. Mineral-melt partition coefficients (D) of redox-sensitive elements (Cr, Eu, V, Fe) vary systematically with fO2 and, in some cases, temperature (e.g. DCr in clinopyroxene). Clinopyroxene sector zoning is common; sectors along a- and b-axes have higher AlIV, AlVI, Cr and Ti and lower Mg than c-axis sectors. In terms of coupled substitutions, clinopyroxene CaTs (MgSi = AlVIAlIV) prevails under oxidized conditions (≥ NNO), where Fe3+ balances the charge, but is limited under reduced conditions. Overall, AlIV is maximised under high temperature, oxidizing conditions and in slowly grown (a-b) sectors. High AlIV facilitates incorporation of REE (REEAlIV = CaSi), but DREE (except DEu) show no systematic dependence on fO2 across the experimental suite. In sector zoned clinopyroxenes enrichment in REE3+ in Al-rich sectors is quantitatively consistent with the greater availability of suitably-charged M2 lattice sites and the electrostatic energy penalty required to insert REE3+ onto unsuitably-charged M2 sites. By combining our experimental results with published data, we explore the potential for trace element oxybarometry. We show that olivine-melt DV, clinopyroxene-melt DV/DSc and plagioclase-melt DEu/DSr all have potential as oxybarometers and we present expressions for these as a function of fO2 relative to NNO. The crystal chemical sensitivity of heterovalent cation incorporation into clinopyroxene and the melt compositional sensitivity of the Eu2+-Eu3+ redox potential limit the use of clinopyroxene-melt and plagioclase-melt, however, olivine-melt DV affords considerable precision and accuracy as an oxybarometer that is independent of temperature, and crystal and melt composition. Variation of DV and DV/DSc with fO2 for olivine and clinopyroxene contains information on redox speciation of V in coexisting melt. By comparing the redox speciation constraints from partitioning to data from Fe-free synthetic systems and XANES spectroscopy of quenched glasses, we show that homogenous equilibria involving Fe and V species modify V speciation on quench, leading to a net overall reduction in the average vanadium valence. Mineral-melt partitioning of polyvalent species can be a useful probe of redox speciation in Fe-bearing systems that is unaffected by quench effects. Supplementary Information: The online version contains supplementary material available at 10.1007/s00410-023-02069-x.

Formation and dynamics of magma reservoirs.

The emerging concept of a magma reservoir is one in which regions containing melt extend from the source of magma generation to the surface. The reservoir may contain regions of very low fraction intergranular melt, partially molten rock (mush) and melt lenses (or magma chambers) containing high melt fraction eruptible magma, as well as pockets of exsolved magmatic fluids. The various parts of the system may be separated by a sub-solidus rock or be connected and continuous. Magma reservoirs and their wall rocks span a vast array of rheological properties, covering as much as 25 orders of magnitude from high viscosity, sub-solidus crustal rocks to magmatic fluids. Time scales of processes within magma reservoirs range from very slow melt and fluid segregation within mush and magma chambers and deformation of surrounding host rocks to very rapid development of magma and fluid instability, transport and eruption. Developing a comprehensive model of these systems is a grand challenge that will require close collaboration between modellers, geophysicists, geochemists, geologists, volcanologists and petrologists. This article is part of the Theo Murphy meeting issue 'Magma reservoir architecture and dynamics'.

Melting during late-stage rifting in Afar is hot and deep.

Investigations of a variety of continental rifts and margins worldwide have revealed that a considerable volume of melt can intrude into the crust during continental breakup, modifying its composition and thermal structure. However, it is unclear whether the cause of voluminous melt production at volcanic rifts is primarily increased mantle temperature or plate thinning. Also disputed is the extent to which plate stretching or thinning is uniform or varies with depth with the entire continental lithospheric mantle potentially being removed before plate rupture. Here we show that the extensive magmatism during rifting along the southern Red Sea rift in Afar, a unique region of sub-aerial transition from continental to oceanic rifting, is driven by deep melting of hotter-than-normal asthenosphere. Petrogenetic modelling shows that melts are predominantly generated at depths greater than 80 kilometres, implying the existence of a thick upper thermo-mechanical boundary layer in a rift system approaching the point of plate rupture. Numerical modelling of rift development shows that when breakup occurs at the slow extension rates observed in Afar, the survival of a thick plate is an inevitable consequence of conductive cooling of the lithosphere, even when the underlying asthenosphere is hot. Sustained magmatic activity during rifting in Afar thus requires persistently high mantle temperatures, which would allow melting at high pressure beneath the thick plate. If extensive plate thinning does occur during breakup it must do so abruptly at a late stage, immediately before the formation of the new ocean basin.

Primary carbonatite melt from deeply subducted oceanic crust.

Partial melting in the Earth's mantle plays an important part in generating the geochemical and isotopic diversity observed in volcanic rocks at the surface. Identifying the composition of these primary melts in the mantle is crucial for establishing links between mantle geochemical 'reservoirs' and fundamental geodynamic processes. Mineral inclusions in natural diamonds have provided a unique window into such deep mantle processes. Here we provide experimental and geochemical evidence that silicate mineral inclusions in diamonds from Juina, Brazil, crystallized from primary and evolved carbonatite melts in the mantle transition zone and deep upper mantle. The incompatible trace element abundances calculated for a melt coexisting with a calcium-titanium-silicate perovskite inclusion indicate deep melting of carbonated oceanic crust, probably at transition-zone depths. Further to perovskite, calcic-majorite garnet inclusions record crystallization in the deep upper mantle from an evolved melt that closely resembles estimates of primitive carbonatite on the basis of volcanic rocks. Small-degree melts of subducted crust can be viewed as agents of chemical mass-transfer in the upper mantle and transition zone, leaving a chemical imprint of ocean crust that can possibly endure for billions of years.

Chemical feedbacks during magma degassing control chlorine partitioning and metal extraction in volcanic arcs.

Hydrous fluids released from subducting oceanic lithosphere fuel arc magmatism and associated hydrothermal mineralization, including formation of porphyry copper deposits. Critical magma degassing parameters are the depth, chemistry and style of fluid release during magma ascent, notably the behaviour of chlorine, a key metal-transporting ligand. Currently, understanding is limited by restricted data on fluid-melt partitioning of chlorine as a function of pressure and magma chemistry, and the complex interplay between the two that occurs in polybaric magmatic systems. Here we present experimental determinations of chlorine partitioning as a function of fluid and melt composition at pressures from 50 to 800 MPa. We provide, for the first time, a quantitative understanding of chlorine and copper evolution that is valid for shallow, deep or transcrustal differentiation and degassing. Monte Carlo simulations using our new data reproduce the chemical evolution of melt inclusions from arc volcanoes and fluid inclusions from upper crustal intrusions and porphyry copper deposits. Our results not only provide a novel chemical framework for understanding magma degassing, but quantify the primacy of magmatic chlorine concentration at the point of fluid saturation in promoting efficient copper extraction from magmas.

Dopant incorporation into garnet solid solutions--a breakdown of Goldschmidt's first rule.

Atomistic simulations suggest trace elements are more soluble in a 50:50 pyrope (Mg3Al2Si3O12)-grossular (Ca3Al2Si3O12) garnet mixture than in either end-member; consistent with partitioning experiments, and, contrary to Goldschmidt's first rule, large trace element cations may substitute for Mg2+, small trace elements for Ca2+.

The 'zero charge' partitioning behaviour of noble gases during mantle melting.

Noble-gas geochemistry is an important tool for understanding planetary processes from accretion to mantle dynamics and atmospheric formation. Central to much of the modelling of such processes is the crystal-melt partitioning of noble gases during mantle melting, magma ascent and near-surface degassing. Geochemists have traditionally considered the 'inert' noble gases to be extremely incompatible elements, with almost 100 per cent extraction efficiency from the solid phase during melting processes. Previously published experimental data on partitioning between crystalline silicates and melts has, however, suggested that noble gases approach compatible behaviour, and a significant proportion should therefore remain in the mantle during melt extraction. Here we present experimental data to show that noble gases are more incompatible than previously demonstrated, but not necessarily to the extent assumed or required by geochemical models. Independent atomistic computer simulations indicate that noble gases can be considered as species of 'zero charge' incorporated at crystal lattice sites. Together with the lattice strain model, this provides a theoretical framework with which to model noble-gas geochemistry as a function of residual mantle mineralogy.