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.

Fluid-assisted grain size reduction leads to strain localization in oceanic transform faults

Oceanic Transform Faults are major plate boundaries representing the most seismogenic part of the mid ocean ridge system. Nonetheless, their structure and deformation mechanisms at depth are largely unknown due to rare exposures of deep sections. Here we study the mineral fabric of deformed mantle peridotites - ultramafic mylonites - collected from the transpressive Atobá ridge, along the northern fault of the St. Paul transform system in the Equatorial Atlantic Ocean. We show that, at pressure and temperature conditions of the lower oceanic lithosphere, the dominant deformation mechanism is fluid-assisted dissolution-precipitation creep. Grain size reduction during deformation is enhanced by dissolution of coarser pyroxene grains in presence of fluid and contextual precipitation of small interstitial ones, leading to strain localization at lower stresses than dislocation creep. This mechanism potentially represents the dominant weakening factor in the oceanic lithosphere and a main driver for the onset and maintenance of oceanic transform faults.

Lower oceanic crust formed by in situ melt crystallisation revealed by seismic layering

Oceanic crust forms at mid-ocean spreading centres through a combination of magmatic and tectonic processes, with the magmatic processes creating two distinct layers: the upper and the lower crust. While the upper crust is known to form from lava flows and basaltic dikes based on geophysical and drilling results, the formation of the gabbroic lower crust is still debated. Here we perform a full waveform inversion of wide-angle seismic data from relatively young (7-12-million-year-old) crust formed at the slow spreading Mid-Atlantic Ridge. The seismic velocity model reveals alternating, 400-500 m thick, high and low velocity layers with ±200 m/s velocity variations, below ~2 km from the oceanic basement. The uppermost low-velocity layer is consistent with hydrothermal alteration, defining the base of extensive hydrothermal circulation near the ridge axis. The underlying layering supports that the lower crust is formed through the intrusion of melt as sills at different depths, that cool and crystallise in situ. The layering extends up to 5-15 km distance along the seismic profile, covering 300,000-800,000 years, suggesting that this form of lower crustal accretion is a stable process.

Magmatic water content controls the pre-eruptive depth of arc magmas

Vanguard efforts in forecasting volcanic eruptions are turning to physics-based models, which require quantitative estimates of magma conditions during pre-eruptive storage. Below active arc volcanoes, observed magma storage depths vary widely (~0 to 20 kilometers) and are commonly assumed to represent levels of neutral buoyancy. Here we show that geophysically observed magma depths (6 ± 3 kilometers) are greater than depths of neutral buoyancy, ruling out this commonly assumed control. Observed depths are instead consistent with predicted depths of water degassing. Intrinsically wetter magmas degas water and crystallize deeper than dry magmas, resulting in viscosity increases that lead to deeper stalling of ascending magma. The water-depth relationship provides a critical constraint for forecasting models by connecting depth of eruption initiation to its volatile fuel.

Three-dimensional magnetic stripes require slow cooling in fast-spread lower ocean crust

Earth's magnetic field is recorded as oceanic crust cools, generating lineated magnetic anomalies that provide the pattern of polarity reversals for the past 160 million years¹. In the lower (gabbroic) crust, polarity interval boundaries are proxies for isotherms that constrain cooling and hence crustal accretion. Seismic observations^2-4, geospeedometry^5-7 and thermal modelling^8-10 of fast-spread crust yield conflicting interpretations of where and how heat is lost near the ridge, a sensitive indicator of processes of melt transport and crystallization within the crust. Here we show that the magnetic structure of magmatically robust fast-spread crust requires that crustal temperatures near the dike-gabbro transition remain at approximately 500 degrees Celsius for 0.1 million years. Near-bottom magnetization solutions over two areas, separated by approximately 8 kilometres, highlight subhorizontal polarity boundaries within 200 metres of the dike-gabbro transition that extend 7-8 kilometres off-axis. Oriented samples with multiple polarity components provide direct confirmation of a corresponding horizontal polarity boundary across an area approximately one kilometre wide, and indicate slow cooling over three polarity intervals. Our results are incompatible with deep hydrothermal cooling within a few kilometres of the axis^2,7 and instead suggest a broad, hot axial zone that extends roughly 8 kilometres off-axis in magmatically robust fast-spread ocean crust.

Deep roots for mid-ocean-ridge volcanoes revealed by plagioclase-hosted melt inclusions

The global mid-ocean ridge system is the most extensive magmatic system on our planet and is the site of 75 per cent of Earth's volcanism¹. The vertical extent of mid-ocean-ridge magmatic systems has been considered to be restricted: even at the ultraslow-spreading Gakkel mid-ocean ridge under the Arctic Ocean, where the lithosphere is thickest, crystallization depths of magmas that feed eruptions are thought to be less than nine kilometres². These depths were determined using the volatile-element contents of melt inclusions, which are small volumes of magma that become trapped within crystallizing minerals. In studies of basaltic magmatic systems, olivine is the mineral of choice for this approach^2-6. However, pressures derived from olivine-hosted melt inclusions are at odds with pressures derived from basalt major-element barometers⁷ and geophysical measurements of lithospheric thickness⁸. Here we present a comparative study of olivine- and plagioclase-hosted melt inclusions from the Gakkel mid-ocean ridge. We show that the volatile contents of plagioclase-hosted melt inclusions correspond to much higher crystallization pressures (with a mean value of 270 megapascals) than olivine-hosted melt inclusions (with a mean value of 145 megapascals). The highest recorded pressure that we find equates to a depth 16.4 kilometres below the seafloor. Such higher depths are consistent with both the thickness of the Gakkel mid-ocean ridge lithosphere and with pressures reconstructed from glass compositions. In contrast to previous studies using olivine-hosted melt inclusions, our results demonstrate that mid-ocean-ridge volcanoes may have magmatic roots deep in the lithospheric mantle, at least at ultraslow-spreading ridges.

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

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.

Hybrid shallow on-axis and deep off-axis hydrothermal circulation at fast-spreading ridges

Hydrothermal flow at oceanic spreading centres accounts for about ten per cent of all heat flux in the oceans and controls the thermal structure of young oceanic plates. It also influences ocean and crustal chemistry, provides a basis for chemosynthetic ecosystems, and has formed massive sulphide ore deposits throughout Earth's history. Despite this, how and under what conditions heat is extracted, in particular from the lower crust, remains largely unclear. Here we present high-resolution, whole-crust, two- and three-dimensional simulations of hydrothermal flow beneath fast-spreading ridges that predict the existence of two interacting flow components, controlled by different physical mechanisms, that merge above the melt lens to feed ridge-centred vent sites. Shallow on-axis flow structures develop owing to the thermodynamic properties of water, whereas deeper off-axis flow is strongly shaped by crustal permeability, particularly the brittle-ductile transition. About 60 per cent of the discharging fluid mass is replenished on-axis by warm (up to 300 degrees Celsius) recharge flow surrounding the hot thermal plumes, and the remaining 40 per cent or so occurs as colder and broader recharge up to several kilometres away from the axis that feeds hot (500-700 degrees Celsius) deep-rooted off-axis flow towards the ridge. Despite its lower contribution to the total mass flux, this deep off-axis flow carries about 70 per cent of the thermal energy released at the ridge axis. This combination of two flow components explains the seismically determined thermal structure of the crust and reconciles previously incompatible models favouring either shallower on-axis or deeper off-axis hydrothermal circulation.

Seismic reflection images of a near-axis melt sill within the lower crust at the Juan de Fuca ridge

The oceanic crust extends over two-thirds of the Earth's solid surface, and is generated along mid-ocean ridges from melts derived from the upwelling mantle. The upper and middle crust are constructed by dyking and sea-floor eruptions originating from magma accumulated in mid-crustal lenses at the spreading axis, but the style of accretion of the lower oceanic crust is actively debated. Models based on geological and petrological data from ophiolites propose that the lower oceanic crust is accreted from melt sills intruded at multiple levels between the Moho transition zone (MTZ) and the mid-crustal lens, consistent with geophysical studies that suggest the presence of melt within the lower crust. However, seismic images of molten sills within the lower crust have been elusive. Until now, only seismic reflections from mid-crustal melt lenses and sills within the MTZ have been described, suggesting that melt is efficiently transported through the lower crust. Here we report deep crustal seismic reflections off the southern Juan de Fuca ridge that we interpret as originating from a molten sill at present accreting the lower oceanic crust. The sill sits 5-6 km beneath the sea floor and 850-900 m above the MTZ, and is located 1.4-3.2 km off the spreading axis. Our results provide evidence for the existence of low-permeability barriers to melt migration within the lower section of modern oceanic crust forming at intermediate-to-fast spreading rates, as inferred from ophiolite studies.

Drilling to Gabbro in Intact Ocean Crust

Sampling an intact sequence of oceanic crust through lavas, dikes, and gabbros is necessary to advance the understanding of the formation and evolution of crust formed at mid-ocean ridges, but it has been an elusive goal of scientific ocean drilling for decades. Recent drilling in the eastern Pacific Ocean in Hole 1256D reached gabbro within seismic layer 2, 1157 meters into crust formed at a superfast spreading rate. The gabbros are the crystallized melt lenses that formed beneath a mid-ocean ridge. The depth at which gabbro was reached confirms predictions extrapolated from seismic experiments at modern mid-ocean ridges: Melt lenses occur at shallower depths at faster spreading rates. The gabbros intrude metamorphosed sheeted dikes and have compositions similar to the overlying lavas, precluding formation of the cumulate lower oceanic crust from melt lenses so far penetrated by Hole 1256D.

Frozen magma lenses below the oceanic crust

The Earth's oceanic crust crystallizes from magmatic systems generated at mid-ocean ridges. Whereas a single magma body residing within the mid-crust is thought to be responsible for the generation of the upper oceanic crust, it remains unclear if the lower crust is formed from the same magma body, or if it mainly crystallizes from magma lenses located at the base of the crust. Thermal modelling, tomography, compliance and wide-angle seismic studies, supported by geological evidence, suggest the presence of gabbroic-melt accumulations within the Moho transition zone in the vicinity of fast- to intermediate-spreading centres. Until now, however, no reflection images have been obtained of such a structure within the Moho transition zone. Here we show images of groups of Moho transition zone reflection events that resulted from the analysis of approximately 1,500 km of multichannel seismic data collected across the intermediate-spreading-rate Juan de Fuca ridge. From our observations we suggest that gabbro lenses and melt accumulations embedded within dunite or residual mantle peridotite are the most probable cause for the observed reflectivity, thus providing support for the hypothesis that the crust is generated from multiple magma bodies.

Modes of faulting at mid-ocean ridges

Abyssal-hill-bounding faults that pervade the oceanic crust are the most common tectonic feature on the surface of the Earth. The recognition that these faults form at plate spreading centres came with the plate tectonic revolution. Recent observations reveal a large range of fault sizes and orientations; numerical models of plate separation, dyke intrusion and faulting require at least two distinct mechanisms of fault formation at ridges to explain these observations. Plate unbending with distance from the top of an axial high reproduces the observed dip directions and offsets of faults formed at fast-spreading centres. Conversely, plate stretching, with differing amounts of constant-rate magmatic dyke intrusion, can explain the great variety of fault offset seen at slow-spreading ridges. Very-large-offset normal faults only form when about half the plate separation at a ridge is accommodated by dyke intrusion.

Bending-related faulting and mantle serpentinization at the Middle America trench

The dehydration of subducting oceanic crust and upper mantle has been inferred both to promote the partial melting leading to arc magmatism and to induce intraslab intermediate-depth earthquakes, at depths of 50-300 km. Yet there is still no consensus about how slab hydration occurs or where and how much chemically bound water is stored within the crust and mantle of the incoming plate. Here we document that bending-related faulting of the incoming plate at the Middle America trench creates a pervasive tectonic fabric that cuts across the crust, penetrating deep into the mantle. Faulting is active across the entire ocean trench slope, promoting hydration of the cold crust and upper mantle surrounding these deep active faults. The along-strike length and depth of penetration of these faults are also similar to the dimensions of the rupture area of intermediate-depth earthquakes.