Mid-ocean-ridge seismicity reveals extreme types of ocean lithosphere

Highly variable magmatic accretion at the ultraslow-spreading Gakkel Ridge

Crustal accretion at mid-ocean ridges governs the creation and evolution of the oceanic lithosphere. Generally accepted models1-4 of passive mantle upwelling and melting predict notably decreased crustal thickness at a spreading rate of less than 20 mm year-1. We conducted the first, to our knowledge, high-resolution ocean-bottom seismometer (OBS) experiment at the Gakkel Ridge in the Arctic Ocean and imaged the crustal structure of the slowest-spreading ridge on the Earth. Unexpectedly, we find that crustal thickness ranges between 3.3 km and 8.9 km along the ridge axis and it increased from about 4.5 km to about 7.5 km over the past 5 Myr in an across-axis profile. The highly variable crustal thickness and relatively large average value does not align with the prediction of passive mantle upwelling models. Instead, it can be explained by a model of buoyant active mantle flow driven by thermal and compositional density changes owing to melt extraction. The influence of active versus passive upwelling is predicted to increase with decreasing spreading rate. The process of active mantle upwelling is anticipated to be primarily influenced by mantle temperature and composition. This implies that the observed variability in crustal accretion, which includes notably varied crustal thickness, is probably an inherent characteristic of ultraslow-spreading ridges.

Microseismicity and lithosphere thickness at a nearly-amagmatic oceanic detachment fault system

Oceanic detachment faults play a central role in accommodating the plate divergence at slow-ultraslow spreading mid-ocean ridges. Successive flip-flop detachment faults in a nearly-amagmatic region of the ultraslow spreading Southwest Indian Ridge (SWIR) at 64°30'E accommodate ~100% of plate divergence, with mostly ultramafic smooth seafloor. Here we present microseismicity data, recorded by ocean bottom seismometers, showing that the axial brittle lithosphere is on the order of 15 km thick under the nearly-amagmatic smooth seafloor, which is no thicker than under nearby volcanic seafloor or at more magmatic SWIR detachment systems. Our data reveal that microearthquakes with normal focal mechanisms are colocated with seismically-imaged damage zones of the active detachment fault and of antithetic hanging-wall faults. The level of the hanging-wall seismicity is significantly higher than that documented at more magmatic detachments of slow-ultraslow ridges, which may be a unique feature of nearly-amagmatic flip-flop detachment systems.

Giant sponge grounds of Central Arctic seamounts are associated with extinct seep life

The Central Arctic Ocean is one of the most oligotrophic oceans on Earth because of its sea-ice cover and short productive season. Nonetheless, across the peaks of extinct volcanic seamounts of the Langseth Ridge (87°N, 61°E), we observe a surprisingly dense benthic biomass. Bacteriosponges are the most abundant fauna within this community, with a mass of 460 g C m⁻² and an estimated carbon demand of around 110 g C m⁻² yr⁻¹, despite export fluxes from regional primary productivity only sufficient to provide <1% of this required carbon. Observed sponge distribution, bulk and compound-specific isotope data of fatty acids suggest that the sponge microbiome taps into refractory dissolved and particulate organic matter, including remnants of an extinct seep community. The metabolic profile of bacteriosponge fatty acids and expressed genes indicate that autotrophic symbionts contribute significantly to carbon assimilation. We suggest that this hotspot ecosystem is unique to the Central Arctic and associated with extinct seep biota, once fueled by degassing of the volcanic mounts.

This study reports the discovery of dense sponge gardens across the peaks of permanently ice-covered, extinct volcanic seamounts of the Langseth Ridge and on the remnants of a now extinct seep ecosystem. Using approaches to sample and infer food and energy sources to this ice-covered community, the authors suggest that the sponges use refractory organic matter trapped in the extinct seep community on which they sit.

Spatial and Temporal Variations in Crustal Structure in the Eurasian Basin

To understand the tectonic–magmatic history, crustal structure and crustal accretion mode of the Eurasian Basin in the Arctic, we calculated the crustal thickness, residual bathymetry (RB) and non-isostatic topography of the Eurasian Basin by using the latest bathymetry, free-air gravity anomaly, crustal age and sediment thickness data. The tectonic–magmatic process of the Eurasian Basin can be divided into two stages. During magnetic isochrons 24-13 (C24-C13), the crustal structure of the basin was characterized by regional variations, while the crustal structure of the basin was characterized by local variations during magnetic isochrons 13-0 (C13-0). On the whole, the western part of the basin had a thicker crustal thickness and higher RB than the eastern part of the basin during C24-C13, which should result from the northward movement of Greenland. During magnetic isochrons 24-20 (C24-C20), the crustal structure of the eastern part of the basin had abnormally strong asymmetry. We speculate that there may be mantle upwelling beneath the Kara Sea Shelf in the south of the Eastern Eurasian Basin, which provides a large amount of melt for the crustal accretion of the southern part of the Eastern Eurasian Basin. The melt focusing supply could generate abnormally thick crust (>7 km) during magma enhancement period. The Western Eurasian Basin had stronger spatial variability and more frequency asymmetric polarity reversal than the Eastern Eurasian Basin during magnetic isochrons 6-0 (C6-0). We attribute this to the inflow of the North Atlantic mantle.

Thermochemical anomalies in the upper mantle control Gakkel Ridge accretion

Despite progress in understanding seafloor accretion at ultraslow spreading ridges, the ultimate driving force is still unknown. Here we use ⁴⁰Ar/³⁹Ar isotopic dating of mid-ocean ridge basalts recovered at variable distances from the axis of the Gakkel Ridge to provide new constraints on the spatial and temporal distribution of volcanic eruptions at various sections of an ultraslow spreading ridge. Our age data show that magmatic-dominated sections of the Gakkel Ridge spread at a steady rate of ~11.1 ± 0.9 mm/yr whereas amagmatic sections have a more widely distributed melt supply yielding ambiguous spreading rate information. These variations in spreading rate and crustal accretion correlate with locations of hotter thermochemical anomalies in the asthenosphere beneath the ridge. We conclude therefore that seafloor generation in ultra-slow spreading centres broadly reflects the distribution of thermochemical anomalies in the upper mantle.

13 million years of seafloor spreading throughout the Red Sea Basin

The crustal and tectonic structure of the Red Sea and especially the maximum northward extent of the (ultra)slow Red Sea spreading centre has been debated—mainly due to a lack of detailed data. Here, we use a compilation of earthquake and vertical gravity gradient data together with high-resolution bathymetry to show that ocean spreading is occurring throughout the entire basin and is similar in style to that at other (ultra)slow spreading mid-ocean ridges globally, with only one first-order offset along the axis. Off-axis traces of axial volcanic highs, typical features of (ultra)slow-spreading ridges, are clearly visible in gravity data although buried under thick salt and sediments. This allows us to define a minimum off-axis extent of oceanic crust of <55 km off the coast along the complete basin. Hence, the Red Sea is a mature ocean basin in which spreading began along its entire length 13 Ma ago.

Here, based on earthquake data, vertical gravity gradient data and high-resolution bathymetry, the authors show that the Red Sea is not in transition from rifting to spreading as previously proposed. They instead suggest it to be a mature ocean basin in which continuous seafloor spreading began quasi-instantaneously along its entire length around 13 Ma ago.

Semibrittle seismic deformation in high-temperature mantle mylonite shear zone along the Romanche transform fault

Oceanic transform faults, a key element of plate tectonics, represent the first-order discontinuities along mid-ocean ridges, host large earthquakes, and induce extreme thermal gradients in lithosphere. However, the thermal structure along transform faults and its effects on earthquake generation are poorly understood. Here we report the presence of a 10- to 15-kilometer-thick in-depth band of microseismicity in 10 to 34 kilometer depth range associated with a high-temperature (700° to 900°C) mantle below the brittle lithosphere along the Romanche mega transform fault in the equatorial Atlantic Ocean. The occurrence of the shallow 2016 moment magnitude 7.1 supershear rupture earthquake and these deep microearthquakes indicate that although large earthquakes occur in the upper brittle lithosphere, a substantial amount of deformation is accommodated in the semibrittle mylonitic mantle that resides at depths below the 600°C isotherm. We also observe a rapid westward deepening of this band of seismicity indicating a strong lateral heterogeneity.

Global variation of seismic energy release with oceanic lithosphere age

Variations in Mid Ocean Ridge seismicity with age provide a new tool to understand the thermal evolution of the oceanic lithosphere. The sum of seismic energy released by earthquakes during a time, and for an area, is proportional to its lithospheric age. Asthenospheric temperatures emerge on ridge centers with new crust resulting in high seismic activity; thus, the energy released sum is highest on the young lithosphere and decreases with age. We propose a general model that relates the systematic variation of seismic energy released with the lithospheric age. Our analysis evaluates the main physical factors involved in the changes of energy released sum with the oceanic lithosphere age in MOR systems of different spreading rates. These observations are substantiated based on three cross-sections of the East Pacific Rise, six sections in the Mid Atlantic Ridge, and three profiles in the Central Indian Ridge. Our global model provides an additional tool for understanding tectonic processes, including the effects of seismicity and mid-plate volcanism, and a better understanding of the thermal evolution for the young oceanic lithosphere.

Migration of Natural Hydrogen from Deep-Seated Sources in the São Francisco Basin, Brazil

Hydrogen gas is seeping from the sedimentary basin of São Franciso, Brazil. The seepages of H2 are accompanied by helium, whose isotopes reveal a strong crustal signature. Geophysical data indicates that this intra-cratonic basin is characterized by (i) a relatively high geothermal gradient, (ii) deep faults delineating a horst and graben structure and affecting the entire sedimentary sequence, (iii) archean to paleoproterozoïc basements enriched in radiogenic elements and displaying mafic and ultramafic units, and (iv) a possible karstic reservoir located 400 m below the surface. The high geothermal gradient could be due to a thin lithosphere enriched in radiogenic elements, which can also contribute to a massive radiolysis process of water at depth, releasing a significant amount of H2. Alternatively, ultramafic rocks that may have generated H2 during their serpentinization are also documented in the basement. The seismic profiles show that the faults seen at the surface are deeply rooted in the basement, and can drain deep fluids to shallow depths in a short time scale. The carbonate reservoirs within the Bambuí group which forms the main part of the sedimentary layers, are crossed by the fault system and represent good candidates for temporary H2 accumulation zones. The formation by chemical dissolution of sinkholes located at 400 m depth might explain the presence of sub-circular depressions seen at the surface. These sinkholes might control the migration of gas from temporary storage reservoirs in the upper layer of the Bambuí formation to the surface. The fluxes of H2 escaping out of these structures, which have been recently documented, are discussed in light of the newly developed H2 production model in the Precambrian continental crust.

Deep high-temperature hydrothermal circulation in a detachment faulting system on the ultra-slow spreading ridge

Coupled magmatic and tectonic activity plays an important role in high-temperature hydrothermal circulation at mid-ocean ridges. The circulation patterns for such systems have been elucidated by microearthquakes and geochemical data over a broad spectrum of spreading rates, but such data have not been generally available for ultra-slow spreading ridges. Here we report new geophysical and fluid geochemical data for high-temperature active hydrothermal venting at Dragon Horn area (49.7°E) on the Southwest Indian Ridge. Twin detachment faults penetrating to the depth of 13 ± 2 km below the seafloor were identified based on the microearthquakes. The geochemical composition of the hydrothermal fluids suggests a long reaction path involving both mafic and ultramafic lithologies. Combined with numerical simulations, our results demonstrate that these hydrothermal fluids could circulate ~ 6 km deeper than the Moho boundary and to much greater depths than those at Trans-Atlantic Geotraverse and Logachev-1 hydrothermal fields on the Mid-Atlantic Ridge.

Magmatic and tectonic activity at mid-oceanic ridges can give detailed insights into high-temperature hydrothermal circulation of fluids. The authors here present geochemical and geophysical datasets that suggest a hydrothermal system penetrating the upper lithospheric mantle at an ultra-slow spreading mid-oceanic ridge.

Simulating Serpentinization as It Could Apply to the Emergence of Life Using the JPL Hydrothermal Reactor

The molecules feeding life's emergence are thought to have been provided through the hydrothermal interactions of convecting carbonic ocean waters with minerals comprising the early Hadean oceanic crust. Few laboratory experiments have simulated ancient hydrothermal conditions to test this conjecture. We used the JPL hydrothermal flow reactor to investigate CO₂ reduction in simulated ancient alkaline convective systems over 3 days (T = 120°C, P = 100 bar, pH = 11). H₂-rich hydrothermal simulant and CO₂-rich ocean simulant solutions were periodically driven in 4-h cycles through synthetic mafic and ultramafic substrates and Fe>Ni sulfides. The resulting reductants included micromoles of HS^- and formate accompanied possibly by micromoles of acetate and intermittent minor bursts of methane as ascertained by isotopic labeling. The formate concentrations directly correlated with the CO₂ input as well as with millimoles of Mg²⁺ ions, whereas the acetate did not. Also, tens of micromoles of methane were drawn continuously from the reactor materials during what appeared to be the onset of serpentinization. These results support the hypothesis that formate may have been delivered directly to a branch of an emerging acetyl coenzyme-A pathway, thus obviating the need for the very first hydrogenation of CO₂ to be made in a hydrothermal mound. Another feed to early metabolism could have been methane, likely mostly leached from primary CH₄ present in the original Hadean crust or emanating from the mantle. That a small volume of methane was produced sporadically from the ¹³CO₂-feed, perhaps from transient occlusions, echoes the mixed results and interpretations from other laboratories. As serpentinization and hydrothermal leaching can occur wherever an ocean convects within anhydrous olivine- and sulfide-rich crust, these results may be generalized to other wet rocky planets and moons in our solar system and beyond.

Abiotic methane synthesis and serpentinization in olivine-hosted fluid inclusions

The conditions of methane (CH4) formation in olivine-hosted secondary fluid inclusions and their prevalence in peridotite and gabbroic rocks from a wide range of geological settings were assessed using confocal Raman spectroscopy, optical and scanning electron microscopy, electron microprobe analysis, and thermodynamic modeling. Detailed examination of 160 samples from ultraslow- to fast-spreading midocean ridges, subduction zones, and ophiolites revealed that hydrogen (H2) and CH4 formation linked to serpentinization within olivine-hosted secondary fluid inclusions is a widespread process. Fluid inclusion contents are dominated by serpentine, brucite, and magnetite, as well as CH4(g) and H2(g) in varying proportions, consistent with serpentinization under strongly reducing, closed-system conditions. Thermodynamic constraints indicate that aqueous fluids entering the upper mantle or lower oceanic crust are trapped in olivine as secondary fluid inclusions at temperatures higher than ∼400 °C. When temperatures decrease below ∼340 °C, serpentinization of olivine lining the walls of the fluid inclusions leads to a near-quantitative consumption of trapped liquid H2O. The generation of molecular H2 through precipitation of Fe(III)-rich daughter minerals results in conditions that are conducive to the reduction of inorganic carbon and the formation of CH4 Once formed, CH4(g) and H2(g) can be stored over geological timescales until extracted by dissolution or fracturing of the olivine host. Fluid inclusions represent a widespread and significant source of abiotic CH4 and H2 in submarine and subaerial vent systems on Earth, and possibly elsewhere in the solar system.

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.

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.

Deep electrical imaging of the ultraslow-spreading Mohns Ridge

More than a third of mid-ocean ridges have a spreading rate of less than 20 millimetres a year¹. The lack of deep imaging data means that factors controlling melting and mantle upwelling^2,3, the depth to the lithosphere-asthenosphere boundary (LAB)^4,5, crustal thickness^6-9 and hydrothermal venting are not well understood for ultraslow-spreading ridges^10,11. Modern electromagnetic data have greatly improved our understanding of fast-spreading ridges^12,13, but have not been available for the ultraslow-spreading ridges. Here we present a detailed 120-kilometre-deep electromagnetic joint inversion model for the ultraslow-spreading Mohns Ridge, combining controlled source electromagnetic and magnetotelluric data. Inversion images show mantle upwelling focused along a narrow, oblique and strongly asymmetric zone coinciding with asymmetric surface uplift. Although the upwelling pattern shows several of the characteristics of a dynamic system^3,12-14, it probably reflects passive upwelling controlled by slow and asymmetric plate movements instead. Upwelling asthenosphere and melt can be traced to the inferred depth of the Mohorovičić discontinuity and are enveloped by the resistivity (100 ohm metres) contour denoted the electrical LAB (eLAB). The eLAB may represent a rheological boundary defined by a minimum melt content. We also find that neither the melt-suppression model⁷ nor the inhibited-migration model¹⁵, which explain the correlation between spreading rate and crustal thickness^6,16-19, can explain the thin crust below the ridge. A model in which crustal thickness is directly controlled by the melt-producing rock volumes created by the separating plates is more likely. Active melt emplacement into oceanic crust about three kilometres thick culminates in an inferred crustal magma chamber draped by fluid convection cells emanating at the Loki's Castle hydrothermal black smoker field. Fluid convection extends for long lateral distances, exploiting high porosity at mid-crustal levels. The magnitude and long-lived nature of such plumbing systems could promote venting 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'.

Tracking the Serpentinite Feet of the Mediterranean Salt Giant

Interpretation of seismic profiles and results of scientific drillings in the Mediterranean subseafloor provided indication of gigantic salt deposits which rarely crop out on land, such as in Sicily. The salt giants were ascribed to the desiccation, driven by the solar energy, of the entire basin. Nevertheless, the evaporite model hardly explains deep-sea salt deposits. This paper considers a different hypothesis suggesting that seawater reached NaCl saturation during serpentinization of ultramafic rocks. Solid salts and brine pockets were buried within the serpentinite bodies being later (e.g., in the Messinian) released, due to serpentinite breakdown, and discharged at seafloor as hydrothermal heavy brines. Therefore, sea-bottom layers of brine at gypsum and halite saturation were formed. The model is applicable to the Mediterranean area since geophysical data revealed relicts of an aged (hence serpentinized) oceanic lithosphere, of Tethyan affinity, both in its western “Atlantic” extension (Gulf of Cádiz) and in eastern basins, and xenoliths from Hyblean diatremes (Sicily) provided evidence of buried serpentinites in the central area. In addition, the buoyant behavior of muddled serpentinite and salts (and hydrocarbons) gave rise to many composite diapirs throughout the Mediterranean area. Thus, the Mediterranean “salt giant” consists of several independent geobodies of serpentinite and salts.

Magma plumbing system and seismicity of an active mid-ocean ridge volcano

At mid-ocean ridges volcanism generally decreases with spreading rate but surprisingly massive volcanic centres occur at the slowest spreading ridges. These volcanoes can host unexpectedly strong earthquakes and vigorous, explosive submarine eruptions. Our understanding of the geodynamic processes forming these volcanic centres is still incomplete due to a lack of geophysical data and the difficulty to capture their rare phases of magmatic activity. We present a local earthquake tomographic image of the magma plumbing system beneath the Segment 8 volcano at the ultraslow-spreading Southwest Indian Ridge. The tomography shows a confined domain of partial melt under the volcano. We infer that from there melt is horizontally transported to a neighbouring ridge segment at 35 km distance where microearthquake swarms and intrusion tremor occur that suggest ongoing magmatic activity. Teleseismic earthquakes around the Segment 8 volcano, prior to our study, indicate that the current magmatic spreading episode may already have lasted over a decade and hence its temporal extent greatly exceeds the frequent short-lived spreading episodes at faster opening mid-ocean ridges.