A dearth of intermediate melts at subduction zone volcanoes and the petrogenesis of arc andesites

The pristine precursor of Andean-type magmatism preserved in magma mingling zones

Intermediate magma compositions have been postulated to be parental to Andean-type magmatism in the recent years. Geochemical and experimental methods have allowed the modelling of a hypothetical parental composition that accounts for the major element trends displayed by Andean-type batholiths. However, natural plutonic examples matching the modelled composition remain lacking, likely due to the predominance of fractionated liquids and cumulates in the batholiths after protracted and large-scale differentiation. Contrary to this, magma mingling zones, a common feature in Andean-type batholiths, are characterised by quenching phenomena, minimising differentiation. In this paper, we present data from intermediate magmatism in the world-class Gerena magma mingling zone in the Seville Sierra Norte batholith (southern Iberia), compositionally equivalent to Andean-type magmatic series. Geochemical data from quenched dark globules of variable scale and the corresponding host granodiorites are contrasted with the bimodal trends displayed by the host batholith. Results suggest that the smaller-scale dark globules have not undergone any significant fractionation. Furthermore, after conducting geochemical modelling we conclude the dark globules represent a composition that could be parental to Andean-type magmas. We propose that magma mingling zones are an optimal place to probe for parental magmas of Andean-type magmatism, particularly those represented in pristine melanocratic, intermediate globules.

Alkaline magmas in shallow arc plutonic roots: a field and experimental investigation of hydrous cumulate melting in the southern Adamello batholith

Despite the first-order importance of crystallisation-differentiation for arc magma evolution, several other processes contribute to their compositional diversity. Among them is the remelting of partly crystallised magmas, also known as cumulate melting or 'petrological cannibalism'. The impact of this process on the plutonic record is poorly constrained. We investigate a nepheline-normative dyke suite close to the Blumone gabbros, a large amphibole-gabbro unit of the Tertiary Southern Alpine Adamello igneous complex. The compositions of the studied dykes are characterised by low SiO2 (43-46 wt. %), MgO (5.0-7.2 wt. %), Ni (18-40 μg/g), and high Al2O3 (20.2-22.0 wt. %) contents. Phenocrystic plagioclase in these dykes exhibits major, trace, and Sr isotope compositions similar to Blumone cumulate plagioclase, suggesting a genetic link between the nepheline-normative dykes and the amphibole-gabbro cumulates. We tested this hypothesis by performing saturation experiments on a nepheline-normative dyke composition in an externally heated pressure vessel at 200 MPa between 975 and 1100 °C at fO2 conditions close to the Ni-NiO buffer. Plagioclase and spinel are near-liquidus phases at and above 1050 °C, contrasting with the typical near-liquidus olivine ± spinel assemblage in hydrous calc-alkaline basalts. The alkaline nature of the dykes results from the abundance of amphibole in the protolith, consistent with melting of amphibole-gabbro cumulates. We modelled the heat budget from the repeated injection of basaltic andesite into a partly crystallised amphibole-gabbro cumulate. The results of this model show that no more than 7% of the cumulate pile reaches temperatures high enough to produce nepheline-normative melts. We propose that such nepheline-normative dykes are a hallmark of hydrous cumulate melting in subvolcanic plumbing systems. Therefore, ne-normative dykes in arc batholiths may indicate periods with high magma fluxes.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00410-023-02047-3.

Global volcanic rock classification of Holocene volcanoes

This data descriptor assigns the major and minor rock names from worldwide Holocene volcanoes of the Global Volcanism Program (GVP) using the Total Alkali-Silica diagram (TAS) for the chemical classification of volcanic rocks using the Geochemistry of Rocks of the Oceans and Continents (GEOROC) database. The precompiled files of the GEOROC database provide the chemical composition of volcanic rock samples, from which we computed major and minor rocks for global Holocene volcanoes reported in GVP. The combined dataset associates each volcano with the relative abundance of each volcanic sample type (whole rock, glass, melt inclusion) and provides the five major (more than 10% abundance) and minor rock names. In total, over 138,000 GEOROC volcanic rock samples were considered, for ~1000 Holocene volcanoes. The resulting major rock compositions are in general consistent with those given in GVP. The dataset provides a global panorama of rock composition for Holocene volcanoes.

Redetermination of the Zalantun Group in the ARong Qi Area of Da Hinggan Mountains (Northeastern China): Evidence from Petrology, Geochronology and Geochemistry

The Da Hinggan Mountains are an important area in the tectonic evolution of the Central Asian Orogenic Belt (CAOB), and there are disputes over the closure time of the Paleo-Asian Ocean and the amalgamation spatiotemporal relationship between the Ergun-Hinggan Massif and the Songnun Massif. The geochronology and geochemistry of the Late Cambrian-Late Silurian volcanic rock assemblages in the ARong Qi area at the eastern margin of the Da Hinggan Mountains are studied in this paper. The results suggest that the U-Pb zircon ages of the Late Cambrian, Late Ordovician and Late Silurian volcanic rock assemblages are 507.5 ± 1.0 Ma, 456.2 ± 1.0 Ma, 446.1 ± 0.95 Ma and 423.3 ± 1.4 Ma, respectively. The Late Cambrian-Late Silurian volcanic rocks are quasi-aluminous-peraluminous, belonging to calc-alkaline-shoshonite series, which is rich in HREE but has insignificant europium anomalies. There are abundant large ion lithophile elements (LILE) in the rock, and remarkable Nb, Ta and Ti negative anomalies. The previous data and the current study indicate that a continental margin arc tectonic setting existed in the ARong Qi-Zalantun region during the Early Paleozoic, which is inferred to be the product of the subduction-accretion-amalgamation of the plates along the eastern margin of the Ergun Massif during the Early Paleozoic.

Zinc systematics quantify crustal thickness control on fractionating assemblages of arc magmas

Understanding the processes leading to the broad chemical variability of arc magmas is an essential, yet not fully elucidated, issue in Earth Sciences. Here, I show that Zn–MgO–SiO₂ systematics of magmatic arc rocks correlate significantly with arc thickness. Because Zn–MgO–SiO₂ systematics are mostly controlled by fractionation of different mineral phases, this suggests a systematic change in the proportions of fractionating mineral assemblages depending on arc thickness. Using a mass balance model with a Monte Carlo approach, I show that Zn–MgO–SiO₂ systematics can be quantitatively explained by a continuous transition from plagioclase-dominated fractionating assemblages in thin arcs to amphibole-garnet-magnetite-dominated assemblages in increasingly thicker arcs. Most likely, such a systematic change results from the increase of average depth of magma differentiation that is ultimately controlled by arc thickness. Results presented have implications on the causes of different geochemical trends in arcs, the role of arcs as H₂O filters, and their association with porphyry deposits.

Quantitative chemical mapping of plagioclase as a tool for the interpretation of volcanic stratigraphy: an example from Saint Kitts, Lesser Antilles

Establishing a quantitative link between magmatic processes occurring at depth and volcanic eruption dynamics is essential to forecast the future behaviour of volcanoes, and to correctly interpret monitoring signals at active centres. Chemical zoning in minerals, which captures successive events or states within a magmatic system, can be exploited for such a purpose. However, to develop a quantitative understanding of magmatic systems requires an unbiased, reproducible method for characterising zoned crystals. We use image segmentation on thin section scale chemical maps to segment textural zones in plagioclase phenocrysts. These zones are then correlated throughout a stratigraphic sequence from Saint Kitts (Lesser Antilles), composed of a basal pyroclastic flow deposit and a series of fall deposits. Both segmented phenocrysts and unsegmented matrix plagioclase are chemically decoupled from whole rock geochemical trends, with the latter showing a systematic temporal progression towards less chemically evolved magma (more anorthitic plagioclase). By working on a stratigraphic sequence, it is possible to track the chemical and textural complexity of segmented plagioclase in time, in this case on the order of millennia. In doing so, we find a relationship between the number of crystal populations, deposit thickness and time. Thicker deposits contain a larger number of crystal populations, alongside an overall reduction in this number towards the top of the deposit. Our approach provides quantitative textural parameters for volcanic and plutonic rocks, including the ability to measure the amount of crystal fracturing. In combination with mineral chemistry, these parameters can strengthen the link between petrology and volcanology, paving the way towards a deeper understanding of the magmatic processes controlling eruptive dynamics.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s00445-021-01476-x.

Constraining crustal silica on ancient Earth

Accurately quantifying the composition of continental crust on Hadean and Archean Earth is critical to our understanding of the physiography, tectonics, and climate of our planet at the dawn of life. One longstanding paradigm involves the growth of a relatively mafic planetary crust over the first 1 to 2 billion years of Earth history, implying a lack of modern plate tectonics and a paucity of subaerial crust, and consequently lacking an efficient mechanism to regulate climate. Others have proposed a more uniformitarian view in which Archean and Hadean continents were only slightly more mafic than at present. Apart from complications in assessing early crustal composition introduced by crustal preservation and sampling biases, effects such as the secular cooling of Earth's mantle and the biologically driven oxidation of Earth's atmosphere have not been fully investigated. We find that the former complicates efforts to infer crustal silica from compatible or incompatible element abundances, while the latter undermines estimates of crustal silica content inferred from terrigenous sediments. Accounting for these complications, we find that the data are most parsimoniously explained by a model with nearly constant crustal silica since at least the early Archean.

Evidence from plutonic xenoliths for magma differentiation, mixing and storage in a volatile-rich crystal mush beneath St. Eustatius, Lesser Antilles

Quantifying the storage conditions and evolution of different magmatic components within sub-volcanic plumbing systems is key to our understanding of igneous processes and products. Whereas erupted magmas represent a portion of the eruptible volcanic system, plutonic xenoliths provide a complementary record of the mushy roots of the plumbing system that cannot be mobilised easily to form lavas and consequently offer a unique record of magma diversity within the sub-volcanic plumbing system. Here, we present a detailed petrological and geochemical study of erupted plutonic xenoliths from the island of Sint Eustatius (Statia), in the northern Lesser Antilles volcanic arc. The plutonic xenoliths are predominantly gabbroic, but vary in texture, mineral assemblage and crystallisation sequence. We report major, trace and volatile (H2O and CO2) concentrations of xenolith-hosted melt inclusions (MIs) and interstitial glass. The MIs have a very large range in major element (49-78 wt% SiO2 and 0.1-6.1 wt% MgO) and trace element concentration (72-377 ppm Sr, 32-686 ppm Ba, 39-211 ppm Zr). Their chemistry varies systematically with host phase and sample type. Significantly, it shows that (1) plutonic xenoliths record a complete differentiation sequence from basalt to rhyolite (2) apatite, but not zircon, saturation was reached during crystallisation, (3) amphibole breakdown reactions play a role in the genesis of shallow gabbronorite assemblages, and (4) mixing between crystal cargos and multiple discrete bodies occurred. Residual melt volatile contents are high (≤ 9.1 wt% H2O and ≤ 1350 ppm CO2), returning volatile saturation pressures of 0-426 MPa. Multiple reaction geobarometry and experimental comparisons indicate that equilibration took place in the upper-middle crust (0-15 km). We infer that the Statia plutonic xenoliths represent portions of a large heterogeneous crystal mush within which a great diversity of melts was stored and mixed prior to eruption. Our data show that compositional variations in magmatic plumbing systems exceed those observed in volcanic products, a likely consequence of the blending that occurs prior to and during eruption.

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 space^1,10. However, many common features of crustal magmatism have not yet been explained by either the 'chamber' or 'mush reservoir' concepts^1,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 crystals¹⁰. Reactive melt flow in mush reservoirs produces the low-crystallinity, chemically differentiated (silicic) magmas that ascend to form shallower intrusions or erupt to the surface^11-13. These magmas can host much older crystals, stored at low and even sub-solidus temperatures, consistent with crystal chemistry data^6-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 reservoirs^14,15. Trace-element profiles generated by reactive flow are similar to those observed in a well studied reservoir now exposed at the surface¹⁶. 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 chambers¹⁴.

A Geochemical Overview of Mid-Archaean Metavolcanic Rocks from Southwest Greenland

The present contribution reviews bulk-rock geochemical data for mid-Archaean (ca. 3075–2840 Ma) metavolcanic rocks from the North Atlantic Craton of southwest Greenland. The data set includes the most recent high quality major and trace element geochemical analyses for ten different supracrustal/greenstone belts in the region. When distilling the data set to only include the least altered metavolcanic rocks, by filtering out obviously altered samples, mafic/ultramafic cumulate rocks, late-stage intrusive sheets (dolerites) and migmatites, the remaining data (N = 427) reveal two fundamentally distinct geochemical suites. The contrasting trends that emerge from the filtered geochemical data set, which best represents the melt compositions for these mid-Archaean metavolcanic rocks are: (1) tholeiitic (mainly basaltic) versus (2) calc-alkaline (mainly andesitic). These two rock suites are effectively separated by their La/Sm ratios (below or above three, respectively). It is demonstrated by geochemical modelling that the two contrasting suites cannot be related by either fractional crystallization or crustal assimilation processes, despite occurring within the same metavolcanic sequences. The tholeiitic basaltic rocks were directly mantle-derived, whereas the petrogenesis of the calc-alkaline andesitic rocks involve a significant (>50%) felsic component. The felsic contribution in the calc-alkaline suite could either represent slab-melt metasomatism of their mantle source, mafic-felsic magma mixing, or very large degrees of partial melting of mafic lower crust. At face value, the occurrence of andesites, and the negative Nb-Ta-Ti-anomalies of both suites, is consistent with a subduction zone setting for the origin of these metavolcanic rocks. However, the latter geochemical feature is inherent to processes involving crustal partial melts, and therefore independent lines of evidence are needed to substantiate the hypothesis that plate tectonic processes were already operating by the mid-Archaean.

Volcanic crystals as time capsules of eruption history

Crystals formed prior to a volcanic event can provide evidence of processes leading to and timing of eruptions. Clinopyroxene is common in basaltic to intermediate volcanoes, however, its ability as a recorder of pre-eruptive histories has remained comparatively underexplored. Here we show that novel high-resolution trace element images of clinopyroxene track eruption triggers and timescales at Mount Etna (Sicily, Italy). Chromium (Cr) distribution in clinopyroxene from 1974 to 2014 eruptions reveals punctuated episodes of intrusion of primitive magma at depth. Magma mixing efficiently triggered volcanism (success rate up to 90%), within only 2 weeks of arrival of mafic intrusions. Clinopyroxene zonations distinguish between injections of mafic magma and regular recharges with more evolved magma, which often fail to tip the system to erupt. High Cr zonations can therefore be used to reconstruct past eruptions and inform responses to geophysical signals of volcano unrest, potentially offering an additional approach to volcano hazard monitoring.

Crystals can record the timing of volcanic eruptive triggers at depth by examining their zoning histories. Here, the authors analyse clinopyroxene crystal zoning from eruptions at Mount Etna, Italy from 1974-2014 and show that the intrusion of magma can trigger an eruption within 2 weeks of arrival.

Petrological and experimental evidence for differentiation of water-rich magmas beneath St. Kitts, Lesser Antilles

St. Kitts lies in the northern Lesser Antilles, a subduction-related intraoceanic volcanic arc known for its magmatic diversity and unusually abundant cognate xenoliths. We combine the geochemistry of xenoliths, melt inclusions and lavas with high pressure-temperature experiments to explore magma differentiation processes beneath St. Kitts. Lavas range from basalt to rhyolite, with predominant andesites and basaltic andesites. Xenoliths, dominated by calcic plagioclase and amphibole, typically in reaction relationship with pyroxenes and olivine, can be divided into plutonic and cumulate varieties based on mineral textures and compositions. Cumulate varieties, formed primarily by the accumulation of liquidus phases, comprise ensembles that represent instantaneous solid compositions from one or more magma batches; plutonic varieties have mineralogy and textures consistent with protracted solidification of magmatic mush. Mineral chemistry in lavas and xenoliths is subtly different. For example, plagioclase with unusually high anorthite content (An≤100) occurs in some plutonic xenoliths, whereas the most calcic plagioclase in cumulate xenoliths and lavas are An97 and An95, respectively. Fluid-saturated, equilibrium crystallisation experiments were performed on a St. Kitts basaltic andesite, with three different fluid compositions (XH2O = 1.0, 0.66 and 0.33) at 2.4 kbar, 950-1025 °C, and fO2 = NNO - 0.6 to NNO + 1.2 log units. Experiments reproduce lava liquid lines of descent and many xenolith assemblages, but fail to match xenolith and lava phenocryst mineral compositions, notably the very An-rich plagioclase. The strong positive correlation between experimentally determined plagioclase-melt KdCa-Na and dissolved H2O in the melt, together with the occurrence of Al-rich mafic lavas, suggests that parental magmas were water-rich (> 9 wt% H2O) basaltic andesites that crystallised over a wide pressure range (1.5-6 kbar). Comparison of experimental and natural (lava, xenolith) mafic mineral composition reveals that whereas olivine in lavas is predominantly primocrysts precipitated at low-pressure, pyroxenes and spinel are predominantly xenocrysts formed by disaggregation of plutonic mushes. Overall, St. Kitts xenoliths and lavas testify to mid-crustal differentiation of low-MgO basalt and basaltic andesite magmas within a trans-crustal, magmatic mush system. Lower crustal ultramafic cumulates that relate parental low-MgO basalts to primary, mantle -derived melts are absent on St. Kitts.

Vertically extensive and unstable magmatic systems: A unified view of igneous processes

Volcanoes are an expression of their underlying magmatic systems. Over the past three decades, the classical focus on upper crustal magma chambers has expanded to consider magmatic processes throughout the crust. A transcrustal perspective must balance slow (plate tectonic) rates of melt generation and segregation in the lower crust with new evidence for rapid melt accumulation in the upper crust before many volcanic eruptions. Reconciling these observations is engendering active debate about the physical state, spatial distribution, and longevity of melt in the crust. Here we review evidence for transcrustal magmatic systems and highlight physical processes that might affect the growth and stability of melt-rich layers, focusing particularly on conditions that cause them to destabilize, ascend, and accumulate in voluminous but ephemeral shallow magma chambers.

On the conditions of magma mixing and its bearing on andesite production in the crust

Mixing between magmas is thought to affect a variety of processes, from the growth of continental crust to the triggering of volcanic eruptions, but its thermophysical viability remains unclear. Here, by using high-pressure mixing experiments and thermal calculations, we show that hybridization during single-intrusive events requires injection of high proportions of the replenishing magma during short periods, producing magmas with 55-58 wt% SiO2 when the mafic end-member is basaltic. High strain rates and gas-rich conditions may produce more felsic hybrids. The incremental growth of crustal reservoirs limits the production of hybrids to the waning stage of pluton assembly and to small portions of it. Large-scale mixing appears to be more efficient at lower crustal conditions, but requires higher proportions of mafic melt, producing more mafic hybrids than in shallow reservoirs. Altogether, our results show that hybrid arc magmas correspond to periods of enhanced magma production at depth.