Asymmetric magma plumbing system beneath Axial Seamount based on full waveform inversion of seismic data

The architecture of magma plumbing systems plays a fundamental role in volcano eruption and evolution. However, the precise configuration of crustal magma reservoirs and conduits responsible for supplying eruptions are difficult to explore across most active volcanic systems. Consequently, our understanding of their correlation with eruption dynamics is limited. Axial Seamount is an active submarine volcano located along the Juan de Fuca Ridge, with known eruptions in 1998, 2011, and 2015. Here we present high-resolution images of P-wave velocity, attenuation, and estimates of temperature and partial melt beneath the summit of Axial Seamount, derived from multi-parameter full waveform inversion of a 2D multi-channel seismic line. Multiple magma reservoirs, including a newly discovered western magma reservoir, are identified in the upper crust, with the maximum melt fraction of ~15-32% in the upper main magma reservoir (MMR) and lower fractions of 10% to 26% in other satellite reservoirs. In addition, a feeding conduit below the MMR with a melt fraction of ~4-11% and a low-velocity throat beneath the eastern caldera wall connecting the MMR roof with eruptive fissures are imaged. These findings delineate an asymmetric shallow plumbing system beneath Axial Seamount, providing insights into the magma pathways that fed recent eruptions.

Beyond spreading rate: Controls on the thermal regime of mid-ocean ridges

The thermal state of mid-ocean ridges exerts a crucial modulation on seafloor spreading processes that shape ~2/3 of our planet's surface. Standard thermal models treat the ridge axis as a steady-state boundary layer between the hydrosphere and asthenosphere, whose thermal structure primarily reflects the local spreading rate. This framework explains the deepening of axial melt lenses (AMLs)-a proxy for the basaltic solidus isotherm-from ~1 to ~3 km from fast- to intermediate-spreading ridges but fails to account for shallow crustal AMLs documented at slow-ultraslow spreading ridges. Here, we show that these can be explained by a numerical model that decouples the potentially transient ridge magma supply from spreading rate, captures the essential physics of hydrothermal convection, and considers multiple modes of melt emplacement. Our simulations show that melt flux is a better thermal predictor than spreading rate. While multiple combinations of melt/dike emplacement modes, permeability structure, and temporal fluctuations of melt supply can explain shallow crustal AMLs at slow-ultraslow ridges, they all require elevated melt fluxes compared to most ridge sections of comparable spreading rates. This highlights the importance of along-axis melt focusing at slow-ultraslow ridges and sheds light on the natural variability of their thermal regimes.

Porosity evolution of mafic crystal mush during reactive flow

The emergence of the "mush paradigm" has raised several questions for conventional models of magma storage and extraction: how are melts extracted to form eruptible liquid-rich domains? What mechanism controls melt transport in mush-rich systems? Recently, reactive flow has been proposed as a major contributing factor in the formation of high porosity, melt-rich regions. Yet, owing to the absence of accurate geochemical simulations, the influence of reactive flow on the porosity of natural mush systems remains under-constrained. Here, we use a thermodynamically constrained model of melt-mush reaction to simulate the chemical, mineralogical, and physical consequences of reactive flow in a multi-component mush system. Our results demonstrate that reactive flow within troctolitic to gabbroic mushes can drive large changes in mush porosity. For example, primitive magma recharge causes an increase in the system porosity and could trigger melt channelization or mush destabilization, aiding rapid melt transfer through low-porosity mush reservoirs.

Influence of Petrogenesis on the Engineering Properties of Ultramafic Aggregates and on Their Suitability in Concrete

This study focuses on the use of petrology as a useful tool in construction applications (i.e., concrete). More specifically, this study investigates how the petrogenetic characteristics of ultramafic rocks derived from ophiolite complexes (Veria–Naousa, Gerania) can act as a key tool for the prediction of the final behaviour of ultramafic aggregates as concrete aggregates. For this reason, their petrographic, chemical and mineralogical characteristics were examined and correlated with their engineering properties for evaluating their suitability as concrete aggregates. This study had come up, for the first time, that the genesis environment of the ultramafic rocks is the determinant factor for their physico-mechanical characteristics. Their suitability is relevant to the impact of their mineralogical and structural characteristics both from the two different ophiolite complexes (Veria–Naousa and Gerania). Except serpentinization, the basic alteration process-index of ultramafic rocks, there are also other chemical indices which can be used for ultramafic rocks that may determine their properties. In this context the term ‘fertility rate’ (FR) was introduced which may characterize ultramafic rocks as fertile or not. Furthermore, the Ultramafic Rock Health Index (U.R.H.I.) as well as the Normalized U.R.H.I. (U.R.H.I.N) was also introduced and correlated with the engineering properties of the investigated aggregate rocks. The last index aims to assess and quantify the overall health conditions, encompassing the two major modifying factors that include removal of primary mineral phases, as well as the extent of serpentinization. The main conclusion of this paper is that the genesis environment of the ultramafic rocks is the critical factor that determines their mineralogical, petrographic and chemical characteristics which consequently determines the basic engineering properties of rocks that determine their suitability or not as concrete aggregates.

Volatile Content Implications of Increasing Explosivity of the Strombolian Eruptive Style along the Fracture Opening on the NE Villarrica Flank: Minor Eruptive Centers in the Los Nevados Group 2

Potential flank eruptions at the presently active Villarrica, Southern Andes Volcanic Zone (33.3–46 °S) require the drawing of a comprehensive scenario of eruptive style dynamics, which partially depends on the degassing process. The case we consider in this study is from the Los Nevados Subgroup 2 (LNG2) and constitutes post-glacial minor eruptive centers (MECs) of basaltic–andesitic and basaltic composition, associated with the northeastern Villarrica flank. Petrological studies of the melt inclusions volatile content in olivine determined the pre-eruptive conditions of the shallow magma feeding system (<249 Mpa saturation pressure, 927–1201 °C). The volatile saturation model on “pressure-dependent” volatile species, measured by Fourier Transform Infrared Microspectrometry (FTIR) (H2O of 0.4–3.0 wt.% and CO2 of 114–1586 ppm) and electron microprobe (EMP), revealed that fast cooling pyroclasts like vesicular scoria preserve a ~1.5 times larger amount of CO2, S, Cl, and volatile species contained in melt inclusions from primitive olivine (Fo76–86). Evidence from geological mapping and drone surveys demonstrated the eruption chronology and spatial changes in eruption style from all the local vents along a N45° corridor. The mechanism by which LNG2 is degassed plays a critical role in increasing the explosivity uphill on the Villarrica flank from volcanic vents in the NE sector (<9 km minimum saturation depth) to the SW sector (<8.1 km), where many crystalline ballistic bombs were expulsed, rather than vesicular and spatter scoria.

Evolution of the Crustal and Upper Mantle Seismic Structure From 0-27 Ma in the Equatorial Atlantic Ocean at 2° 43'S

We present seismic tomographic results from a unique seismic refraction and wide-angle survey along a 600 km long flow-line corridor of oceanic lithosphere ranging in age from 0 to 27 Ma in the equatorial Atlantic Ocean at 2° 43'S. The velocities in the crust near the ridge axis rapidly increase in the first 6 Myr and then change gradually with age. The upper crust (Layer 2) thickness varies between 2 and 2.4 km with an average thickness of 2.2 km and the crustal thickness varies from 5.6 to 6 km along the profile with an average crustal thickness of 5.8 km. At some locations, we observe negative velocity anomalies (∼-0.3 km/s) in the lower crust which could be either due to chemical heterogeneity in gabbroic rocks and/or the effects of fault related deformation zones leading to an increase in porosities up to 1.6% depending on the pore/crack geometry. The existence of a low velocity anomaly beneath the ridge axis suggests the presence of partial melt (∼1.3%) in the lower crust. Upper mantle velocities also remain low (∼7.8 km/s) from ridge axis up to 5 Ma, indicating a high temperature regime associated with mantle melting zone underneath. These results suggest that the evolution of the crust and uppermost mantle at this location occur in the first 10 Ma of its formation and then remains unchanged. Most of the structures in the older crust and upper mantle are fossilized structures and could provide information about past processes at ocean spreading centers.

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.

Cryptic evolved melts beneath monotonous basaltic shield volcanoes in the Galápagos Archipelago

Many volcanoes erupt compositionally homogeneous magmas over timescales ranging from decades to millennia. This monotonous activity is thought to reflect a high degree of chemical homogeneity in their magmatic systems, leading to predictable eruptive behaviour. We combine petrological analyses of erupted crystals with new thermodynamic models to characterise the diversity of melts in magmatic systems beneath monotonous shield volcanoes in the Galápagos Archipelago (Wolf and Fernandina). In contrast with the uniform basaltic magmas erupted at the surface over long timescales, we find that the sub-volcanic systems contain extreme heterogeneity, with melts extending to rhyolitic compositions. Evolved melts are in low abundance and large volumes of basalt flushing through the crust from depth overprint their chemical signatures. This process will only maintain monotonous activity while the volume of melt entering the crust is high, raising the possibility of transitions to more silicic activity given a decrease in the crustal melt flux.

In this study the authors show that monotonous basaltic volcanoes can host a range of melts in their sub-volcanic systems, extending to rhyolitic compositions. The study implies that volcanoes which have produced monotonous basaltic lavas on long timescales could transition to more explosive, silica-rich eruptions in the future.

The significance of plagioclase textures in mid-ocean ridge basalt (Gakkel Ridge, Arctic Ocean)

Textures and compositions of minerals can be used to infer the physiochemical conditions present within magmatic systems. Given that plagioclase is an abundant phase in many magmatic systems, understanding the link between texture and process is vital. Here, we present a database of textural and compositional data for > 1800 plagioclase crystals in mid-ocean ridge basalt from the Gakkel Ridge (Arctic Ocean) to investigate the physiochemical conditions and processes that govern the formation of plagioclase textures and compositions. The Gakkel basalts have high modal crystal contents (up to 50%). The crystal cargo is complex, with both individual plagioclase and glomerocrysts showing large variations in crystal habit, zoning and resorption. The most common types of zoning are reverse and patchy; we attribute patchy zoning to infilling following either skeletal growth or resorption. Resorption is abundant, with multiple resorption events commonly present in a single crystal, and results from both magmatic recharge and decompression. Periods of strong undercooling, distinct to quench crystallisation, are indicated by matured skeletal crystals and thin normally zoned melt inclusion-rich bands following resorption. Individual samples often contain diverse textural and compositional plagioclase groups. Furthermore, most plagioclase is not in equilibrium with its host melt. Finally, the porous open structures of some glomerocrysts suggest that they represent pieces of entrained disaggregated mush. We interpret this to indicate that the crystal cargo is not generally phenocrystic in origin. Instead, plagioclase crystals that formed in different parts of a mush-dominated plumbing system were entrained into ascending melts. The textures of individual crystals are a function of their respective histories of (under)cooling, magma mixing and decompression. The morphologies of melt inclusion trapped in the plagioclase crystals are associated with specific host crystal textures, suggesting a link between plagioclase crystallisation processes and melt inclusion entrapment. The database of plagioclase presented herein may serve as a template for the interpretation of plagioclase textures in magmatic systems elsewhere.

Architecture and dynamics of magma reservoirs

This introductory article provides a synopsis of our current understanding of the form and dynamics of magma reservoirs in the crust. This knowledge is based on a range of experimental, observational and theoretical approaches, some of which are multidisclipinary and pioneering. We introduce and provide a contextual background for the papers in this issue, which cover a wide range of topics, encompassing magma storage, transport, behaviour and rheology, as well as the timescales on which magma reservoirs operate. We summarize the key findings that emerged from the meeting and the challenges that remain. The study of magma reservoirs has wide implications not only for understanding geothermal and magmatic systems, but also for natural oil and gas reservoirs and for ore deposit formation. This article is part of the Theo Murphy meeting issue 'Magma reservoir architecture and dynamics'.

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'.

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'.

Mafic tiers and transient mushes: evidence from Iceland

It is well established that magmatism is trans-crustal, with melt storage and processing occurring over a range of depths. Development of this conceptual model was based on observations of the products of magmatism at spreading ridges, including Iceland. Petrological barometry and tracking of the solidification process has been used to show that the Icelandic crust is built by crystallization over a range of depths. The available petrological evidence indicates that most of the active rift zones are not underlain by extensive and pervasive crystal mush. Instead, the microanalytical observations from Iceland are consistent with a model where magmatic processing in the lower crust occurs in sills of decimetric vertical thickness. This stacked sills mode of crustal accretion corresponds to that proposed for the oceanic crust on the basis of ophiolite studies. A key feature of these models is that the country rock for the sills is hot but subsolidus. This condition can be met if the porosity in thin crystal mushes at the margins of the sills is occluded by primitive phases, a contention that is consistent with observations from cumulate nodules in Icelandic basalts. The conditions required for the stabilization of trans-crustal mushes may not be present in magmatic systems at spreading ridges. This article is part of the Theo Murphy meeting issue 'Magma reservoir architecture and dynamics'.

The prevalence of kilometer-scale heterogeneity in the source region of MORB upper mantle

The source regions of mid-ocean ridge basalts (MORB) are heterogeneous, consisting of chemically and lithologically distinct domains of variable size. Partial melting of such heterogeneous mantle sources gives rise to diverse isotopic compositions of MORB and abyssal peridotites. Variations in radiogenic isotope ratios in MORB are attributed to mixing of melts derived from enriched and depleted mantle components. However, melt mixing alone cannot fully account for the difference between the average ¹⁴³Nd/¹⁴⁴Nd in abyssal peridotites and their spatially associated MORB. We show that the more depleted Nd isotope composition in abyssal peridotites is a natural consequence of melt migration-induced mixing or smearing in the melting column. Sub-kilometer scale enriched mantle components or heterogeneities are significantly damped or homogenized in both the residue and erupted melt during their transit through the melting region. Heterogeneities with larger size and higher incompatible trace element abundance are more resistive to the mixing processes. The size-sensitive mixing depends on a parameter called the enrichment strength, which is the product of the heterogeneity size and the ratio between incompatible trace element abundance in the enriched and depleted mantle sources. Observed Nd-Hf isotope variations in MORB and abyssal peridotites can be reproduced if the enrichment strength is 20 to 60 km. These heterogeneities could be on the kilometer scale and have similar isotope ratios to but less incompatible trace element abundances than recycled oceanic crust.

Sensitivity of seafloor bathymetry to climate-driven fluctuations in mid-ocean ridge magma supply

Recent studies have proposed that the bathymetric fabric of the seafloor formed at mid-ocean ridges records rapid (23,000 to 100,000 years) fluctuations in ridge magma supply caused by sealevel changes that modulate melt production in the underlying mantle. Using quantitative models of faulting and magma emplacement, we demonstrate that, in fact, seafloor-shaping processes act as a low-pass filter on variations in magma supply, strongly damping fluctuations shorter than about 100,000 years. We show that the systematic decrease in dominant seafloor wavelengths with increasing spreading rate is best explained by a model of fault growth and abandonment under a steady magma input. This provides a robust framework for deciphering the footprint of mantle melting in the fabric of abyssal hills, the most common topographic feature on Earth.

Rapid hydrothermal cooling above the axial melt lens at fast-spreading mid-ocean ridge

Axial melt lenses sandwiched between the lower oceanic crust and the sheeted dike sequences at fast-spreading mid-ocean ridges are assumed to be the major magma source of oceanic crust accretion. According to the widely discussed “gabbro glacier” model, the formation of the lower oceanic crust requires efficient cooling of the axial melt lens, leading to partial crystallization and crystal-melt mush subsiding down to lower crust. These processes are believed to be controlled by periodical magma replenishment and hydrothermal circulation above the melt lens. Here we quantify the cooling rate above melt lens using chemical zoning of plagioclase from hornfelsic recrystallized sheeted dikes drilled from the East Pacific at the Integrated Ocean Drilling Program Hole 1256D. We estimate the cooling rate using a forward modelling approach based on CaAl-NaSi interdiffusion in plagioclase. The results show that cooling from the peak thermal overprint at 1000–1050°C to 600°C are yielded within about 10–30 years as a result of hydrothermal circulation above melt lens during magma starvation. The estimated rapid hydrothermal cooling explains how the effective heat extraction from melt lens is achieved at fast-spreading mid-ocean ridges.