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.     

Geophysical evidence from the MELT area for compositional controls on oceanic plates

Magnetotelluric and seismic data, collected during the MELT experiment at the southern East Pacific Rise, constrain the distribution of melt beneath this mid-ocean-ridge spreading centre and also the evolution of the oceanic lithosphere during its early cooling history. Here we focus on structures imaged at distances approximately 100 to 350 km east of the ridge crest, corresponding to seafloor ages of approximately 1.3 to 4.5 million years (Myr), where the seismic and electrical conductivity structure is nearly constant and independent of age. Beginning at a depth of about 60 km, we image a large increase in electrical conductivity and a change from isotropic to transversely anisotropic electrical structure, with higher conductivity in the direction of fast propagation for seismic waves. Conductive cooling models predict structure that increases in depth with age, extending to about 30 km at 4.5 Myr ago. We infer, however, that the structure of young oceanic plates is instead controlled by a decrease in water content above a depth of 60 km induced by the melting process beneath the spreading centre.     

Upper-mantle dynamics revealed by helium isotope variations along the southeast Indian ridge

Helium isotope variations in igneous rocks are important for relating isotopic heterogeneity to convective mixing in the Earth's mantle. High 3He/4He ratios at many ocean islands, along with lower and relatively uniform values in mid-ocean-ridge basalts (MORBs), are thought to result from a well mixed upper-mantle source for MORB and a distinct deeper-mantle source for ocean island basalts. At finer scales, 3He/4He variations along mid-ocean ridges have been related to underlying mantle heterogeneity, but relationships between the scales of geochemical segmentation and mantle convection remain enigmatic. Here we present helium isotope data for MORB glasses recovered along approximately 5,800 km of the southeast Indian ridge, and develop an approach to quantitatively relate spatial variations in geochemical and geophysical parameters at the Earth's surface. A point-to-point correlation analysis reveals structure in the helium isotope data at length scales of approximately 150 and approximately 400 km that appears to be related to secondary convection in the underlying mantle.     

Origin of a 'Southern Hemisphere' geochemical signature in the Arctic upper mantle

The Gakkel ridge, which extends under the Arctic ice cap for approximately 1,800 km, is the slowest spreading ocean ridge on Earth. Its spreading created the Eurasian basin, which is isolated from the rest of the oceanic mantle by North America, Eurasia and the Lomonosov ridge. The Gakkel ridge thus provides unique opportunities to investigate the composition of the sub-Arctic mantle and mantle heterogeneity and melting at the lower limits of seafloor spreading. The first results of the 2001 Arctic Mid-Ocean Ridge Expedition (ref. 1) divided the Gakkel ridge into three tectonic segments, composed of robust western and eastern volcanic zones separated by a 'sparsely magmatic zone'. On the basis of Sr-Nd-Pb isotope ratios and trace elements in basalts from the spreading axis, we show that the sparsely magmatic zone contains an abrupt mantle compositional boundary. Basalts to the west of the boundary display affinities to the Southern Hemisphere 'Dupal' isotopic province, whereas those to the east-closest to the Eurasian continent and where the spreading rate is slowest-display affinities to 'Northern Hemisphere' ridges. The western zone is the only known spreading ridge outside the Southern Hemisphere that samples a significant upper-mantle region with Dupal-like characteristics. Although the cause of Dupal mantle has been long debated, we show that the source of this signature beneath the western Gakkel ridge was subcontinental lithospheric mantle that delaminated and became integrated into the convecting Arctic asthenosphere. This occurred as North Atlantic mantle propagated north into the Arctic during the separation of Svalbard and Greenland.     

Whole-mantle convection and the transition-zone water filter

Because of their distinct chemical signatures, ocean-island and mid-ocean-ridge basalts are traditionally inferred to arise from separate, isolated reservoirs in the Earth's mantle. Such mantle reservoir models, however, typically satisfy geochemical constraints, but not geophysical observations. Here we propose an alternative hypothesis that, rather than being divided into isolated reservoirs, the mantle is filtered at the 410-km-deep discontinuity. We propose that, as the ascending ambient mantle (forced up by the downward flux of subducting slabs) rises out of the high-water-solubility transition zone (between the 660 km and 410 km discontinuities) into the low-solubility upper mantle above 410 km, it undergoes dehydration-induced partial melting that filters out incompatible elements. The filtered, dry and depleted solid phase continues to rise to become the source material for mid-ocean-ridge basalts. The wet, enriched melt residue may be denser than the surrounding solid and accordingly trapped at the 410 km boundary until slab entrainment returns it to the deeper mantle. The filter could be suppressed for both mantle plumes (which therefore generate wetter and more enriched ocean-island basalts) as well as the hotter Archaean mantle (thereby allowing for early production of enriched continental crust). We propose that the transition-zone water-filter model can explain many geochemical observations while avoiding the major pitfalls of invoking isolated mantle reservoirs.     

A correlation between mid-ocean-ridge basalt chemistry and distance to continents

To fully understand the structure and dynamics of the Earth's convecting mantle, the origins of temperature variations within the mantle need to be resolved. Different hypotheses have been proposed to account for these temperature variations: for example, heat coming from the decay of radioactive elements or heat flowing out of the Earth's core. In addition, theoretical studies suggest that the thermal properties of continental masses can affect mantle convection, but quantitative data that could allow us to test these models are scarce. To address this latter problem, we have examined the chemistry of mid-ocean-ridge basalt--which reflects the temperature of the source mantle--as a function of the distance of the ridge from the closest continental margin. No correlation is observed for oceanic ridges close to subduction zones or hotspots; subduction zones probably inhibit thermal transfer between the mantle beneath continents and ocean, whereas hotspots influence the major-element chemistry of ridge basalts, which makes their interpretation with respect to mantle temperature more difficult. However, we do observe a significant correlation for mid-oceanic basalts from the Atlantic and Indian oceans. From this, we conclude that the location of continental masses relative to active ridges influences the large-scale thermal structure of the mantle and we estimate that the mantle cools by 0.05 to 0.1 degrees C per kilometre from the continental margins.     

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.     

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.     

Magmatic and amagmatic seafloor generation at the ultraslow-spreading Gakkel ridge,Arctic Ocean

A high-resolution mapping and sampling study of the Gakkel ridge was accomplished during an international ice-breaker expedition to the high Arctic and North Pole in summer 2001. For this slowest-spreading endmember of the global mid-ocean-ridge system, predictions were that magmatism should progressively diminish as the spreading rate decreases along the ridge, and that hydrothermal activity should be rare. Instead, it was found that magmatic variations are irregular, and that hydrothermal activity is abundant. A 300-kilometre-long central amagmatic zone, where mantle peridotites are emplaced directly in the ridge axis, lies between abundant, continuous volcanism in the west, and large, widely spaced volcanic centres in the east. These observations demonstrate that the extent of mantle melting is not a simple function of spreading rate: mantle temperatures at depth or mantle chemistry (or both) must vary significantly along-axis. Highly punctuated volcanism in the absence of ridge offsets suggests that first-order ridge segmentation is controlled by mantle processes of melting and melt segregation. The strong focusing of magmatic activity coupled with faulting may account for the unexpectedly high levels of hydrothermal activity observed.     

Melting in the Earth's deep upper mantle caused by carbon dioxide

The onset of partial melting beneath mid-ocean ridges governs the cycling of highly incompatible elements from the mantle to the crust, the flux of key volatiles (such as CO2, He and Ar) and the rheological properties of the upper mantle. Geophysical observations indicate that melting beneath ridges begins at depths approaching 300 km, but the cause of this melting has remained unclear. Here we determine the solidus of carbonated peridotite from 3 to 10 GPa and demonstrate that melting beneath ridges may occur at depths up to 330 km, producing 0.03-0.3% carbonatite liquid. We argue that these melts promote recrystallization and realignment of the mineral matrix, which may explain the geophysical observations. Extraction of incipient carbonatite melts from deep within the oceanic mantle produces an abundant source of metasomatic fluids and a vast mantle residue depleted in highly incompatible elements and fractionated in key parent-daughter elements. We infer that carbon, helium, argon and highly incompatible heat-producing elements (such as uranium, thorium and potassium) are efficiently scavenged from depths of approximately 200-330 km in the upper mantle.     

Osmium isotopic constraints on the nature of the DUPAL anomaly from Indian mid-ocean-ridge basalts

The isotopic compositions of mid-ocean-ridge basalts (MORB) from the Indian Ocean have led to the identification of a large-scale isotopic anomaly relative to Pacific and Atlantic ocean MORB. Constraining the origin of this so-called DUPAL anomaly may lead to a better understanding of the genesis of upper-mantle heterogeneity. Previous isotopic studies have proposed recycling of ancient subcontinental lithospheric mantle or sediments with oceanic crust to be responsible for the DUPAL signature. Here we report Os, Pb, Sr and Nd isotopic compositions of Indian MORB from the Central Indian ridge, the Rodriguez triple junction and the South West Indian ridge. All measured samples have higher (187)Os/(188)Os ratios than the depleted upper-mantle value and Pb, Sr and Nd isotopic compositions that imply the involvement of at least two distinct enriched components in the Indian upper-mantle. Using isotopic and geodynamical arguments, we reject both subcontinental lithospheric mantle and recycled sediments with oceanic crust as the cause of the DUPAL anomaly. Instead, we argue that delamination of lower continental crust may explain the DUPAL isotopic signature of Indian MORB.     

Spreading rate dependence of gravity anomalies along oceanic transform faults

Mid-ocean ridge morphology and crustal accretion are known to depend on the spreading rate of the ridge. Slow-spreading mid-ocean-ridge segments exhibit significant crustal thinning towards transform and non-transform offsets, which is thought to arise from a three-dimensional process of buoyant mantle upwelling and melt migration focused beneath the centres of ridge segments. In contrast, fast-spreading mid-ocean ridges are characterized by smaller, segment-scale variations in crustal thickness, which reflect more uniform mantle upwelling beneath the ridge axis. Here we present a systematic study of the residual mantle Bouguer gravity anomaly of 19 oceanic transform faults that reveals a strong correlation between gravity signature and spreading rate. Previous studies have shown that slow-slipping transform faults are marked by more positive gravity anomalies than their adjacent ridge segments, but our analysis reveals that intermediate and fast-slipping transform faults exhibit more negative gravity anomalies than their adjacent ridge segments. This finding indicates that there is a mass deficit at intermediate- and fast-slipping transform faults, which could reflect increased rock porosity, serpentinization of mantle peridotite, and/or crustal thickening. The most negative anomalies correspond to topographic highs flanking the transform faults, rather than to transform troughs (where deformation is probably focused and porosity and alteration are expected to be greatest), indicating that crustal thickening could be an important contributor to the negative gravity anomalies observed. This finding in turn suggests that three-dimensional magma accretion may occur near intermediate- and fast-slipping transform faults.     

Contrasting crustal production and rapid mantle transitions beneath back-arc ridges

The opening of back-arc basins behind subduction zones progresses from initial rifting near the volcanic arc to seafloor spreading. During this process, the spreading ridge and the volcanic arc separate and lavas erupted at the ridge are predicted to evolve away from being heavily subduction influenced (with high volatile contents derived from the subducting plate). Current models predict gradational, rather than abrupt, changes in the crust formed along the ridge as the inferred broad melting region beneath it migrates away from heavily subduction-influenced mantle. In contrast, here we show that across-strike and along-strike changes in crustal properties at the Eastern Lau spreading centre are large and abrupt, implying correspondingly large discontinuities in the nature of the mantle supplying melt to the ridge axes. With incremental separation of the ridge axis from the volcanic front of as little as 5?km, seafloor morphology changes from shallower complex volcanic landforms to deeper flat sea floor dominated by linear abyssal hills, upper crustal seismic velocities abruptly increase by over 20%, and gravity anomalies and isostasy indicate crustal thinning of more than 1.9?km. We infer that the abrupt changes in crustal properties reflect rapid evolution of the mantle entrained by the ridge, such that stable, broad triangular upwelling regions, as inferred for mid-ocean ridges, cannot form near the mantle wedge corner. Instead, the observations imply a dynamic process in which the ridge upwelling zone preferentially captures water-rich low-viscosity mantle when it is near the arc. As the ridge moves away from the arc, a tipping point is reached at which that material is rapidly released from the upwelling zone, resulting in rapid changes in the character of the crust formed at the ridge.     

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.     

Multiple seismic discontinuities near the base of the transition zone in the Earth's mantle

The seismologically defined boundary between the transition zone in the Earth's mantle (410-660 km depth) and the underlying lower mantle is generally interpreted to result from the breakdown of the gamma-spinel phase of olivine to magnesium-perovskite and magnesiowustite. Laboratory measurements of these transformations of olivine have determined that the phase boundary has a negative Clapeyron slope and does indeed occur near pressures corresponding to the base of the transition zone. But a computational study has indicated that, because of the presence of garnet minerals, multiple seismic discontinuities might exist near a depth of 660 km (ref. 4), which would alter the simple negative correlation of changes in temperature with changes in the depth of the phase boundary. In particular, garnet minerals undergo exothermic transformations near this depth, acting to complicate the phase relations and possibly effecting mantle convection processes in some regions. Here we present seismic evidence that supports the existence of such multiple transitions near a depth of 660 km beneath southern California. The observations are consistent with having been generated by garnet transformations coupling with the dissociation of the gamma-spinel phase of olivine. Temperature anomalies calculated from the imaged discontinuity depths--using Clapeyron slopes determined for the various transformations--generally match those predicted from an independent P-wave velocity model of the region.     

Water-rich basalts at mid-ocean-ridge cold spots

Although water is only present in trace amounts in the suboceanic upper mantle, it is thought to play a significant role in affecting mantle viscosity, melting and the generation of crust at mid-ocean ridges. The concentration of water in oceanic basalts has been observed to stay below 0.2 wt%, except for water-rich basalts sampled near hotspots and generated by 'wet' mantle plumes. Here, however, we report unusually high water content in basaltic glasses from a cold region of the mid-ocean-ridge system in the equatorial Atlantic Ocean. These basalts are sodium-rich, having been generated by low degrees of melting of the mantle, and contain unusually high ratios of light versus heavy rare-earth elements, implying the presence of garnet in the melting region. We infer that water-rich basalts from such regions of thermal minima derive from low degrees of 'wet' melting greater than 60 km deep in the mantle, with minor dilution by melts produced by shallower 'dry' melting--a view supported by numerical modelling. We therefore conclude that oceanic basalts are water-rich not only near hotspots, but also at 'cold spots'.     

Mapping the Hawaiian plume conduit with converted seismic waves

The volcanic edifice of the Hawaiian islands and seamounts, as well as the surrounding area of shallow sea floor known as the Hawaiian swell, are believed to result from the passage of the oceanic lithosphere over a mantle hotspot. Although geochemical and gravity observations indicate the existence of a mantle thermal plume beneath Hawaii, no direct seismic evidence for such a plume in the upper mantle has yet been found. Here we present an analysis of compressional-to-shear (P-to-S) converted seismic phases, recorded on seismograph stations on the Hawaiian islands, that indicate a zone of very low shear-wave velocity (< 4 km s(-1)) starting at 130-140 km depth beneath the central part of the island of Hawaii and extending deeper into the upper mantle. We also find that the upper-mantle transition zone (410-660 km depth) appears to be thinned by up to 40-50 km to the south-southwest of the island of Hawaii. We interpret these observations as localized effects of the Hawaiian plume conduit in the asthenosphere and mantle transition zone with excess temperature of approximately 300 degrees C. Large variations in the transition-zone thickness suggest a lower-mantle origin of the Hawaiian plume similar to the Iceland plume, but our results indicate a 100 degrees C higher temperature for the Hawaiian plume.     

Upper mantle SH velocity structure beneath Qiangtang Terrane by modeling triplicated phases

We constrain SH wave velocity structure for the upper mantle beneath western Qiangtang Terrane by comparing regional distance seismic triplicated waveforms with synthetic seismograms, based on an intermediate event （-220 km） recorded by the INDEPTH-Ⅲ seismic array. The ATIP model reveals a low-velocity anomaly with up to -4% variation at the depth of 190-270 km and a relatively small velocity gradient above the depth of 410 km in the upper mantle, which is in agreement with previous results. In combination with other geological studies, we suggest that the depth of top asthenosphere is 190 km and no large-scale lithosphere thinning occurs in western Qiangtang Terrane, besides, Qiangtang Terrane has the same kind of upper mantle structure as the stable Eurasia.     

Helium solubility in olivine and implications for high 3He/4He in ocean island basalts

High 3He/4He ratios found in ocean island basalts are the main evidence for the existence of an undegassed mantle reservoir. However, models of helium isotope evolution depend critically on the chemical behaviour of helium during mantle melting. It is generally assumed that helium is strongly enriched in mantle melts relative to uranium and thorium, yet estimates of helium partitioning in mantle minerals have produced conflicting results. Here we present experimental measurements of helium solubility in olivine at atmospheric pressure. Natural and synthetic olivines were equilibrated with a 50% helium atmosphere and analysed by crushing in vacuo followed by melting, and yield a minimum olivine-melt partition coefficient of 0.0025 +/- 0.0005 (s.d.) and a maximum of 0.0060 +/- 0.0007 (s.d.). The results indicate that helium might be more compatible than uranium and thorium during mantle melting and that high 3He/4He ratios can be preserved in depleted residues of melting. A depleted source for high 3He/4He ocean island basalts would resolve the apparent discrepancy in the relative helium concentrations of ocean island and mid-ocean-ridge basalts.     

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.