The importance of water to oceanic mantle melting regimes

The formation of basaltic crust at mid-ocean ridges and ocean islands provides a window into the compositional and thermal state of the Earth's upper mantle. But the interpretation of geochemical and crustal-thickness data in terms of magma source parameters depends on our understanding of the melting, melt-extraction and differentiation processes that intervene between the magma source and the crust. Much of the quantitative theory developed to model these processes has neglected the role of water in the mantle and in magma, despite the observed presence of water in ocean-floor basalts. Here we extend two quantitative models of ridge melting, mixing and fractionation to show that the addition of water can cause an increase in total melt production and crustal thickness while causing a decrease in mean extent of melting. This may help to resolve several enigmatic observations in the major- and trace-element chemistry of both normal and hotspot-affected ridge basalts.     

Non-equilibrium degassing and a primordial source for helium in ocean-island volcanism

Radioactive decay of uranium and thorium produces 4He, whereas 3He in the Earth's mantle is not produced by radioactive decay and was only incorporated during accretion-that is, it is primordial. 3He/4He ratios in many ocean-island basalts (OIBs) that erupt at hotspot volcanoes, such as Hawaii and Iceland, can be up to sixfold higher than in mid-ocean ridge basalts (MORBs). This is inferred to be the result of outgassing by melt production at mid-ocean ridges in conjunction with radiogenic ingrowth of 4He, which has led to a volatile-depleted upper mantle (MORB source) with low 3He concentrations and low 3He/4He ratios. Consequently, high 3He/4He ratios in OIBs are conventionally viewed as evidence for an undegassed, primitive mantle source, which is sampled by hot, buoyantly upwelling deep-mantle plumes. However, this conventional model provides no viable explanation of why helium concentrations and elemental ratios of He/Ne and He/Ar in OIBs are an order of magnitude lower than in MORBs. This has been described as the 'helium concentration paradox' and has contributed to a long-standing controversy about the structure and dynamics of the Earth's mantle. Here we show that the helium concentration paradox, as well as the full range of noble-gas concentrations observed in MORB and OIB glasses, can self-consistently be explained by disequilibrium open-system degassing of the erupting magma. We show that a higher CO2 content in OIBs than in MORBs leads to more extensive degassing of helium in OIB magmas and that noble gases in OIB lavas can be derived from a largely undegassed primitive mantle source.     

Explosive volcanism on the ultraslow-spreading Gakkel ridge, Arctic Ocean

Roughly 60% of the Earth's outer surface is composed of oceanic crust formed by volcanic processes at mid-ocean ridges. Although only a small fraction of this vast volcanic terrain has been visually surveyed or sampled, the available evidence suggests that explosive eruptions are rare on mid-ocean ridges, particularly at depths below the critical point for seawater (3,000 m). A pyroclastic deposit has never been observed on the sea floor below 3,000 m, presumably because the volatile content of mid-ocean-ridge basalts is generally too low to produce the gas fractions required for fragmenting a magma at such high hydrostatic pressure. We employed new deep submergence technologies during an International Polar Year expedition to the Gakkel ridge in the Arctic Basin at 85 degrees E, to acquire photographic and video images of 'zero-age' volcanic terrain on this remote, ice-covered ridge. Here we present images revealing that the axial valley at 4,000 m water depth is blanketed with unconsolidated pyroclastic deposits, including bubble wall fragments (limu o Pele), covering a large (>10 km(2)) area. At least 13.5 wt% CO(2) is necessary to fragment magma at these depths, which is about tenfold the highest values previously measured in a mid-ocean-ridge basalt. These observations raise important questions about the accumulation and discharge of magmatic volatiles at ultraslow spreading rates on the Gakkel ridge and demonstrate that large-scale pyroclastic activity is possible along even the deepest portions of the global mid-ocean ridge volcanic system.     

Minimum speed limit for ocean ridge magmatism from 210Pb-226Ra-230Th disequilibria

Although 70 per cent of global crustal magmatism occurs at mid-ocean ridges-where the heat budget controls crustal structure, hydrothermal activity and a vibrant biosphere-the tempo of magmatic inputs in these regions remains poorly understood. Such timescales can be assessed, however, with natural radioactive-decay-chain nuclides, because chemical disruption to secular equilibrium systems initiates parent-daughter disequilibria, which re-equilibrate by the shorter half-life in a pair. Here we use 210Pb-226Ra-230Th radioactive disequilibria and other geochemical attributes in oceanic basalts less than 20 years old to infer that melts of the Earth's mantle can be transported, accumulated and erupted in a few decades. This implies that magmatic conditions can fluctuate rapidly at ridge volcanoes. 210Pb deficits of up to 15 per cent relative to 226Ra occur in normal mid-ocean ridge basalts, with the largest deficits in the most magnesium-rich lavas. The 22-year half-life of 210Pb requires very recent fractionation of these two uranium-series nuclides. Relationships between 210Pb-deficits, (226Ra/230Th) activity ratios and compatible trace-element ratios preclude crustal-magma differentiation or daughter-isotope degassing as the main causes for the signal. A mantle-melting model can simulate observed disequilibria but preservation requires a subsequent mechanism to transport melt rapidly. The likelihood of magmatic disequilibria occurring before melt enters shallow crustal magma bodies also limits differentiation and heat replenishment timescales to decades at the localities studied.     

The influence of ridge migration on the magmatic segmentation of mid-ocean ridges

The Earth's mid-ocean ridges display systematic changes in depth and shape, which subdivide the ridges into discrete spreading segments bounded by transform faults and smaller non-transform offsets of the axis. These morphological changes have been attributed to spatial variations in the supply of magma from the mantle, although the origin of the variations is poorly understood. Here we show that magmatic segmentation of ridges with fast and intermediate spreading rates is directly related to the migration velocity of the spreading axis over the mantle. For over 9,500 km of mid-ocean ridge examined, leading ridge segments in the 'hotspot' reference frame coincide with the shallow magmatically robust segments across 86 per cent of all transform faults and 73 per cent of all second-order discontinuities. We attribute this relationship to asymmetric mantle upwelling and melt production due to ridge migration, with focusing of melt towards ridge segments across discontinuities. The model is consistent with variations in crustal structure across discontinuities of the East Pacific Rise, and may explain variations in depth of melting and the distribution of enriched lavas.     

Evidence from gabbro of the Troodos ophiolite for lateral magma transport along a slow-spreading mid-ocean ridge

The lateral flow of magma and ductile deformation of the lower crust along oceanic spreading axes has been thought to play a significant role in suppressing both mid-ocean ridge segmentation and variations in crustal thickness. Direct investigation of such flow patterns is hampered by the kilometres of water that cover the oceanic crust, but such studies can be made on ophiolites (fragments of oceanic crust accreted to a continent). In the Oman ophiolite, small-scale radial patterns of flow have been mapped along what is thought to be the relict of a fast-spreading mid-ocean ridge. Here we present evidence for broad-scale along-axis flow that has been frozen into the gabbro of the Troodos ophiolite in Cyprus (thought to be representative of a slow-spreading ridge axis). The gabbro suite of Troodos spans nearly 20 km of a segment of a fossil spreading axis, near a ridge-transform intersection. We mapped the pattern of magma flow by analysing the rocks' magnetic fabric at 20 sites widely distributed in the gabbro suite, and by examining the petrographic fabric at 9 sites. We infer an along-axis magma flow for much of the gabbro suite, which indicates that redistribution of melt occurred towards the segment edge in a large depth range of the oceanic crust. Our results support the magma plumbing structure that has been inferred indirectly from a seismic tomography experiment on the slow-spreading Mid-Atlantic Ridge.     

The dynamics of melt and shear localization in partially molten aggregates

The volcanoes that lie along the Earth's tectonic boundaries are fed by melt generated in the mantle. How this melt is extracted and focused to the volcanoes, however, remains an unresolved question. Here we present new theoretical results with implications for melt focusing beneath mid-ocean ridges. By modelling laboratory experiments, we test a formulation for magma dynamics and provide an explanation for localized bands of high-porosity and concentrated shear deformation observed in experiments. These bands emerge and persist at 15 degrees-25 degrees to the plane of shear. Past theoretical work on this system predicted the emergence of melt bands but at an angle inconsistent with experiments. Our results suggest that the observed band angle results from a balance of porosity-weakening and strain-rate-weakening deformation mechanisms. Lower band angles are predicted for greater strain-rate weakening. From these lower band angles, we estimate the orientation of melt bands beneath mid-ocean ridges and show that they may enhance magma focusing toward the ridge axis.     

Geophysical evidence for reduced melt production on the Arctic ultraslow Gakkel mid-ocean ridge

Most models of melt generation beneath mid-ocean ridges predict significant reduction of melt production at ultraslow spreading rates (full spreading rates &<20 mm x yr(-1)) and consequently they predict thinned oceanic crust. The 1,800-km-long Arctic Gakkel mid-ocean ridge is an ideal location to test such models, as it is by far the slowest portion of the global mid-ocean-ridge spreading system, with a full spreading rate ranging from 6 to 13 mm x yr(-1) (refs 4, 5). Furthermore, in contrast to some other ridge systems, the spreading direction on the Gakkel ridge is not oblique and the rift valley is not offset by major transform faults. Here we present seismic evidence for the presence of exceptionally thin crust along the Gakkel ridge rift valley with crustal thicknesses varying between 1.9 and 3.3 km (compared to the more usual value of 7 km found on medium- to fast-spreading mid-ocean ridges). Almost 8,300 km of closely spaced aeromagnetic profiles across the rift valley show the presence of discrete volcanic centres along the ridge, which we interpret as evidence for strongly focused, three-dimensional magma supply. The traces of these eruptive centres can be followed to crustal ages of approximately 25 Myr off-axis, implying that these magma production and transport systems have been stable over this timescale.     

Aluminium control of argon solubility in silicate melts under pressure

Understanding of the crystal chemistry of the Earth's deep mantle has evolved rapidly recently with the gradual acceptance of the importance of the effect of minor elements such as aluminium on the properties of major phases such as perovskite. In the early Earth, during its formation and segregation into rocky mantle and iron-rich core, it is likely that silicate liquids played a large part in the transport of volatiles to or from the deep interior. The importance of aluminium on solubility mechanisms at high pressure has so far received little attention, even though aluminium has long been recognized as exerting strong control on liquid structures at ambient conditions. Here we present constraints on the solubility of argon in aluminosilicate melt compositions up to 25 GPa and 3,000 K, using a laser-heated diamond-anvil cell. The argon contents reach a maximum that persists to pressures as high as 17 GPa (up to 500 km deep in an early magma ocean), well above that expected on the basis of Al-free melt experiments. A distinct drop in argon solubility observed over a narrow pressure range correlates well with the expected void loss in the melt structure predicted by recent molecular dynamics simulations. These results provide a process for noble gas sequestration in the mantle at various depths in a cooling magma ocean. The concept of shallow partial melting as a unique process for extracting noble gases from the early Earth, thereby defining the initial atmospheric abundance, may therefore be oversimplified.     

Skew of mantle upwelling beneath the East Pacific Rise governs segmentation

Mantle upwelling is essential to the generation of new oceanic crust at mid-ocean ridges, and it is generally assumed that such upwelling is symmetric beneath active ridges. Here, however, we use seismic imaging to show that the isotropic and anisotropic structure of the mantle is rotated beneath the East Pacific Rise. The isotropic structure defines the pattern of magma delivery from the mantle to the crust. We find that the segmentation of the rise crest between transform faults correlates well with the distribution of mantle melt. The azimuth of seismic anisotropy constrains the direction of mantle flow, which is rotated nearly 10 degrees anticlockwise from the plate-spreading direction. The mismatch between the locus of mantle melt delivery and the morphologic ridge axis results in systematic differences between areas of on-axis and off-axis melt supply. We conclude that the skew of asthenospheric upwelling and transport governs segmentation of the East Pacific Rise and variations in the intensity of ridge crest processes.     

Isotopic evidence for Late Cretaceous plume-ridge interaction at the Hawaiian hotspot

When a mantle plume interacts with a mid-ocean ridge, both are noticeably affected. The mid-ocean ridge can display anomalously shallow bathymetry, excess volcanism, thickened crust, asymmetric sea-floor spreading and a plume component in the composition of the ridge basalts. The hotspot-related volcanism can be drawn closer to the ridge, and its geochemical composition can also be affected. Here we present Sr-Nd-Pb isotopic analyses of samples from the next-to-oldest seamount in the Hawaiian hotspot track, the Detroit seamount at 51 degrees N, which show that, 81 Myr ago, the Hawaiian hotspot produced volcanism with an isotopic signature indistinguishable from mid-ocean ridge basalt. This composition is unprecedented in the known volcanism from the Hawaiian hotspot, but is consistent with the interpretation from plate reconstructions that the hotspot was located close to a mid-ocean ridge about 80 Myr ago. As the rising mantle plume encountered the hot, low-viscosity asthenosphere and hot, thin lithosphere near the spreading centre, it appears to have entrained enough of the isotopically depleted upper mantle to overwhelm the chemical characteristics of the plume itself. The Hawaiian hotspot thus joins the growing list of hotspots that have interacted with a rift early in their history.     

Evidence from three-dimensional seismic reflectivity images for enhanced melt supply beneath mid-ocean-ridge discontinuities

Quantifying the melt distribution and crustal structure across ridge-axis discontinuities is essential for understanding the relationship between magmatic, tectonic and petrologic segmentation of mid-ocean-ridge spreading centres. The geometry and continuity of magma bodies beneath features such as overlapping spreading centres can strongly influence the composition of erupted lavas and may give insight into the underlying pattern of mantle flow. Here we present three-dimensional images of seismic reflectivity beneath a mid-ocean ridge to investigate the nature of melt distribution across a ridge-axis discontinuity. Reflectivity slices through the 9 degrees 03' N overlapping spreading centre on East Pacific Rise suggest that it has a robust magma supply, with melt bodies underlying both limbs and ponding of melt beneath large areas of the overlap basin. The geometry of melt distribution beneath this offset is inconsistent with large-scale, crustal redistribution of melt away from centres of upwelling. The complex distribution of melt seems instead to be caused by a combination of vertical melt transport from the underlying mantle and subsequent focusing of melt beneath a magma freezing boundary in the mid-crust.     

Cryptic striations in the upper mantle revealed by hafnium isotopes in southeast Indian ridge basalts

The Earth's mantle is isotopically heterogeneous on length scales ranging from centimetres to more than 10(4) kilometres. This heterogeneity originates from partial melt extraction and plate tectonic recycling, whereas stirring during mantle convection tends to reduce it. Here we show that mid-ocean ridge basalts from 2,000 km along the southeast Indian ridge (SEIR) display a bimodal hafnium isotopic distribution. This bimodality reveals the presence of ancient compositional striations (streaks) in the Indian Ocean upper mantle. The number density of the streaks is described by a Poisson distribution, with an average thickness of approximately 40 km. Such a distribution is anticipated for a well-stirred upper mantle, in which heterogeneity is continually introduced by plate tectonic recycling, and redistributed by viscous stretching and convective refolding.     

Mineralogy of the mid-ocean-ridge basalt source from neodymium isotopic composition of abyssal peridotites

Inferring the melting process at mid-ocean ridges, and the physical conditions under which melting takes place, usually relies on the assumption of compositional similarity between all mid-ocean-ridge basalt sources. Models of mantle melting therefore tend to be restricted to those that consider the presence of only one lithology in the mantle, peridotite. Evidence from xenoliths and peridotite massifs show that after peridotite, pyroxenite and eclogite are the most abundant rock types in the mantle. But at mid-ocean ridges, where most of the melting takes place, and in ophiolites, pyroxenite is rarely found. Here we present neodymium isotopic compositions of abyssal peridotites to investigate whether peridotite can indeed be the sole source for mid-ocean-ridge basalts. By comparing the isotopic compositions of basalts and peridotites at two segments of the southwest Indian ridge, we show that a component other than peridotite is required to explain the low end of the (143)Nd/(144)Nd variations of the basalts. This component is likely to have a lower melting temperature than peridotite, such as pyroxenite or eclogite, which could explain why it is not observed at mid-ocean ridges.     

Mantle segmentation along the Oman ophiolite fossil mid-ocean ridge

It has been difficult to relate the segmentation of mid-ocean ridges to processes occurring in the Earth's underlying mantle, as the mantle is rarely sampled directly and chemical variations observed in lavas at the surface are heavily influenced by details of their production as melt extracted from the mantle. Our understanding of such mantle processes has therefore relied on the analysis of pieces of fossil oceanic lithosphere now exposed at the Earth's surface, known as ophiolites. Here we present the phase chemistry and whole-rock major- and trace-element contents of 174 samples of the mantle collected along over 400 km of the Oman Sultanate ophiolite. We show that, when analysed along the fossil ridge, variations of elemental ratios sensitive to the melting process define a three-dimensional geometry of mantle upwellings, which can be related to the segmentation observed in modern mid-ocean ridge environments.     

Spreading-rate dependence of melt extraction at mid-ocean ridges from mantle seismic refraction data

A variety of observations indicate that mid-ocean ridges produce less crust at spreading rates below 20 mm yr(-1) (refs 1-3), reflecting changes in fundamental ridge processes with decreasing spreading rate. The nature of these changes, however, remains uncertain, with end-member explanations being decreasing shallow melting or incomplete melt extraction, each due to the influence of a thicker thermal lid. Here we present results of a seismic refraction experiment designed to study mid-ocean ridge processes by imaging residual mantle structure. Our results reveal an abrupt lateral change in bulk mantle seismic properties associated with a change from slow to ultraslow palaeo-spreading rate. Changes in mantle velocity gradient, basement topography and crustal thickness all correlate with this spreading-rate change. These observations can be explained by variations in melt extraction at the ridge, with a gabbroic phase preferentially retained in the mantle at slower spreading rates. The estimated volume of retained melt balances the approximately 1.5-km difference in crustal thickness, suggesting that changes in spreading rate affect melt-extraction processes rather than total melting.     

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.     

Seismic reflection images of the Moho underlying melt sills at the East Pacific Rise

The determination of melt distribution in the crust and the nature of the crust-mantle boundary (the 'Moho') is fundamental to the understanding of crustal accretion processes at oceanic spreading centres. Upper-crustal magma chambers have been imaged beneath fast- and intermediate-spreading centres but it has been difficult to image structures beneath these magma sills. Using three-dimensional seismic reflection images, here we report the presence of Moho reflections beneath a crustal magma chamber at the 9 degrees 03' N overlapping spreading centre, East Pacific Rise. Our observations highlight the formation of the Moho at zero-aged crust. Over a distance of less than 7 km along the ridge crest, a rapid increase in two-way travel time of seismic waves between the magma chamber and Moho reflections is observed, which we suggest is due to a melt anomaly in the lower crust. The amplitude versus offset variation of reflections from the magma chamber shows a coincident region of higher melt fraction overlying this anomalous region, supporting the conclusion of additional melt at depth.     

Re-visiting barium isotope compositions of mid-ocean ridge basalts and the implications

Barium (Ba) isotopes can be used as potential tracers for crustal material recycling in the mantle. Determination of the Ba isotope composition of the depleted mantle is essential for such applications. However, Ba isotope data for mantle-derived basalts are still rare. In this study, we reported high-precision Ba isotope data of 30 oceanic basalts including 25 mid-ocean ridge basalts (MORBs) from geochemically and geologically diverse mid-ocean ridge segments and five back-arc basin basalts. The δ138/134Ba values of these samples varied from ?0.06‰ to +0.11‰, with no systematic cross-region variation. Together with published data, we constrained the average δ138/134Ba of global MORBs to +0.05‰±0.09‰ (2 standard deviation, n = 51). Based on depleted MORBs that have (La/Sm)_N < 0.8, low 87Sr/86Sr (< 0.70263), and low Ba/Th < 71.3, we estimated the average δ138/134Ba of the depleted MORB mantle (DMM) as + 0.05‰ ± 0.05‰ (2SD, n = 16) that is significantly lower than the DMM (≈ 0.14‰) reported previously. If a new estimation of the DMM is applied, it is unreasonable to infer that the Ba isotope signatures of the “enriched-type” MORBs (E-MORBs) could be attributed to pervasive sediment recycling in the upper mantle. We, therefore, conclude that the Ba isotope compositions of the E-MORBs could be sourced from the incorporation of subducted altered oceanic crust and/or sediments depending on the Ba isotope composition and other geochemical information of the local mantle.     

Ancient, highly heterogeneous mantle beneath Gakkel ridge, Arctic Ocean

The Earth's mantle beneath ocean ridges is widely thought to be depleted by previous melt extraction, but well homogenized by convective stirring. This inference of homogeneity has been complicated by the occurrence of portions enriched in incompatible elements. Here we show that some refractory abyssal peridotites from the ultraslow-spreading Gakkel ridge (Arctic Ocean) have very depleted 187Os/188Os ratios with model ages up to 2 billion years, implying the long-term preservation of refractory domains in the asthenospheric mantle rather than their erasure by mantle convection. The refractory domains would not be sampled by mid-ocean-ridge basalts because they contribute little to the genesis of magmas. We thus suggest that the upwelling mantle beneath mid-ocean ridges is highly heterogeneous, which makes it difficult to constrain its composition by mid-ocean-ridge basalts alone. Furthermore, the existence of ancient domains in oceanic mantle suggests that using osmium model ages to constrain the evolution of continental lithosphere should be approached with caution.