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.     

Contrasting origins of the upper mantle revealed by hafnium and lead isotopes from the Southeast Indian Ridge

The origin of the isotopic signature of Indian mid-ocean ridge basalts has remained enigmatic, because the geochemical composition of these basalts is consistent either with pollution from recycled, ancient altered oceanic crust and sediments, or with ancient continental crust or lithosphere. The radiogenic isotopic signature may therefore be the result of contamination of the upper mantle by plumes containing recycled altered ancient oceanic crust and sediments, detachment and dispersal of continental material into the shallow mantle during rifting and breakup of Gondwana, or contamination of the upper mantle by ancient subduction processes. The identification of a process operating on a scale large enough to affect major portions of the Indian mid-ocean ridge basalt source region has been a long-standing problem. Here we present hafnium and lead isotope data from across the Indian-Pacific mantle boundary at the Australian-Antarctic discordance region of the Southeast Indian Ridge, which demonstrate that the Pacific and Indian upper mantle basalt source domains were each affected by different mechanisms. We infer that the Indian upper-mantle isotope signature in this region is affected mainly by lower continental crust entrained during Gondwana rifting, whereas the isotope signature of the Pacific upper mantle is influenced predominantly by ocean floor subduction-related processes.     

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.     

Re-Os isotopic evidence for long-lived heterogeneity and equilibration processes in the Earth's upper mantle

The geochemical composition of the Earth's upper mantle is thought to reflect 4.5 billion years of melt extraction, as well as the recycling of crustal materials. The fractionation of rhenium and osmium during partial melting in the upper mantle makes the Re-Os isotopic system well suited for tracing the extraction of melt and recycling of the resulting mid-ocean-ridge basalt. Here we report osmium isotope compositions of more than 700 osmium-rich platinum-group element alloys derived from the upper mantle. The osmium isotopic data form a wide, essentially gaussian distribution, demonstrating that, with respect to Re-Os isotope systematics, the upper mantle is extremely heterogeneous. As depleted and enriched domains can apparently remain unequilibrated on a timescale of billions of years, effective equilibration seems to require high degrees of partial melting, such as occur under mid-ocean ridges or in back-arc settings, where percolating melts enhance the mobility of both osmium and rhenium. We infer that the gaussian shape of the osmium isotope distribution is the signature of a random mixing process between depleted and enriched domains, resulting from a 'plum pudding' distribution in the upper mantle, rather than from individual melt depletion events.     

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.     

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.     

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.     

Geochemical tracing of Pacific-to-Atlantic upper-mantle flow through the Drake passage

The Earth's convecting upper mantle can be viewed as comprising three main reservoirs, beneath the Pacific, Atlantic and Indian oceans. Because of the uneven global distribution and migration of ridges and subduction zones, the surface area of the Pacific reservoir is at present contracting at about 0.6 km2 x y(r-1), while the Atlantic and Indian reservoirs are growing at about 0.45 km2 x yr(-1) and 0.15 km2 x yr(-1), respectively. Garfunkel and others have argued that there must accordingly be net mantle flow from the Pacific to the Atlantic and Indian reservoirs (in order to maintain mass balance), and Alvarez further predicted that this flow should be restricted to the few parts of the Pacific rim (here termed 'gateways') where there are no continental roots or subduction zones that might act as barriers to shallow mantle flow. The main Pacific gateways are, according to Alvarez, the southeast Indian Ocean, the Caribbean Sea and the Drake passage. Here we report geochemical data which confirm that there has been some outflow of Pacific mantle into the Drake passage--but probably in response to regional tectonic constraints, rather than global mass-balance requirements. We also show that a mantle domain boundary, equivalent to the Australian-Antarctic discordance, must lie between the Drake passage and the east Scotia Sea.     

North China sub-craton lithospheric structure elucidated through coal mine blasting

We present a super-range seismic observation along the >1300-km-long profile passing through the Yinchuan basin and the Ordos block from the blasting point towards the southeast triggered by a large-(dynamite) scale coal blast in the Ningxia Hui Autonomous Region’s Helan Mountain (yielded a profile). The seismic wave information from the uppermost mantle reflecting different depths was obtained by the China continental seismic survey. Pn refracted waves from the uppermost mantle were effectively traced up t...     

Evidence of lower-mantle slab penetration phases in plate motions

It is well accepted that subduction of the cold lithosphere is a crucial component of the Earth's plate tectonic style of mantle convection. But whether and how subducting plates penetrate into the lower mantle is the subject of continuing debate, which has substantial implications for the chemical and thermal evolution of the mantle. Here we identify lower-mantle slab penetration events by comparing Cenozoic plate motions at the Earth's main subduction zones with motions predicted by fully dynamic models of the upper-mantle phase of subduction, driven solely by downgoing plate density. Whereas subduction of older, intrinsically denser, lithosphere occurs at rates consistent with the model, younger lithosphere (of ages less than about 60 Myr) often subducts up to two times faster, while trench motions are very low. We conclude that the most likely explanation is that older lithosphere, subducting under significant trench retreat, tends to lie down flat above the transition to the high-viscosity lower mantle, whereas younger lithosphere, which is less able to drive trench retreat and deforms more readily, buckles and thickens. Slab thickening enhances buoyancy (volume times density) and thereby Stokes sinking velocity, thus facilitating fast lower-mantle penetration. Such an interpretation is consistent with seismic images of the distribution of subducted material in upper and lower mantle. Thus we identify a direct expression of time-dependent flow between the upper and lower mantle.     

A discontinuity in mantle composition beneath the southwest Indian ridge

The composition of mid-ocean-ridge basalt is known to correlate with attributes such as ridge topography and seismic velocity in the underlying mantle, and these correlations have been interpreted to reflect variations in the average extent and mean pressures of melting during mantle upwelling. In this respect, the eastern extremity of the southwest Indian ridge is of special interest, as its mean depth of 4.7 km (ref. 4), high upper-mantle seismic wave velocities and thin oceanic crust of 4-5 km (ref. 6) suggest the presence of unusually cold mantle beneath the region. Here we show that basaltic glasses dredged in this zone, when compared to other sections of the global mid-ocean-ridge system, have higher Na(8.0), Sr and Al2O3 compositions, very low CaO/Al2O3 ratios relative to TiO2 and depleted heavy rare-earth element distributions. This signature cannot simply be ascribed to low-degree melting of a typical mid-ocean-ridge source mantle, as different geochemical indicators of the extent of melting are mutually inconsistent. Instead, we propose that the mantle beneath approximately 1,000 km of the southwest Indian ridge axis has a complex history involving extensive earlier melting events and interaction with partial melts of a more fertile source.     

Lithospheric structure of the Rio Grande rift

A high-resolution, regional passive seismic experiment in the Rio Grande rift region of the southwestern United States has produced new images of upper-mantle velocity structure and crust-mantle topography. Synthesizing these results with geochemical and other geophysical evidence reveals highly symmetric lower-crustal and upper-mantle lithosphere extensional deformation, suggesting a pure-shear rifting mechanism for the Rio Grande rift. Extension in the lower crust is distributed over a region four times the width of the rift's surface expression. Here we propose that the laterally distributed, pure shear extension is a combined effect of low strain rate and a regionally elevated geotherm, possibly abetted by pre-existing lithospheric structures, at the time of rift initiation. Distributed extension in the lower crust and mantle has induced less concentrated vertical mantle upwelling and less vigorous small-scale convection than would have arisen from more localized deformation. This lack of highly focused mantle upwelling may explain a deficit of rift-related volcanics in the Rio Grande rift compared to other major rift systems such as the Kenya rift.     

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.     

Global systematics of arc volcano position

Global systematics in the location of volcanic arcs above subduction zones are widely considered to be a clue to the melting processes that occur at depth, and the locations of the arcs have often been explained in terms of the release of hydrous fluids near the top of the subducting slab (see, for example, refs 3-6). Grove et al. conclude that arc volcano location is controlled by melting in the mantle at temperatures above the water-saturated upper-mantle solidus and below the upper limit of stability of the mineral chlorite and in particular, that the arc fronts lie directly above the shallowest point of such melt regions in the mantle. Here we show that this conclusion is incorrect because the calculated arc locations of Grove et al. are in error owing to the inadequate spatial resolution of their numerical models, and because the agreement that they find between predicted and observed systematics arises from a spurious correlation between calculated arc location and slab dip. A more informative conclusion to draw from their experiments is that the limits of chlorite stability (figure 1b of ref. 7) cannot explain the global systematics in the depth to the slab beneath the sharply localized arc fronts.     

Anisotropy of thermal diffusivity in the upper mantle.

Heat transfer in the mantle is a key process controlling the Earth's dynamics. Upper-mantle mineral phases, especially olivine, have been shown to display highly anisotropic thermal diffusivity at ambient conditions, and seismic anisotropy data show that preferred orientations of olivine induced by deformation are coherent at large scales (>50 km) in the upper mantle. Thus heat transport in the upper mantle should be anisotropic. But the thermal anisotropy of mantle minerals at high temperature and its relationship with deformation have not been well constrained. Here we present petrophysical modelling and laboratory measurements of thermal diffusivity in deformed mantle rocks between temperatures of 290 and 1,250 K that demonstrate that deformation may induce a significant anisotropy of thermal diffusivity in the uppermost mantle. We found that heat transport parallel to the flow direction is up to 30 per cent faster than that normal to the flow plane. Such a strain-induced thermal anisotropy implies that the upper-mantle temperature distribution, rheology and, consequently, its dynamics, will depend on deformation history. In oceans, resistive drag flow would result in lower vertical diffusivities in both the lithosphere and asthenosphere and hence in less effective heat transfer from the convective mantle. In continents, olivine orientations frozen in the lithosphere may induce anisotropic heating above mantle plumes, favouring the reactivation of pre-existing structures.     

Mantle wedge control on back-arc crustal accretion

At mid-ocean ridges, plate separation leads to upward advection and pressure-release partial melting of fertile mantle material; the melt is then extracted to the spreading centre and the residual depleted mantle flows horizontally away. In back-arc basins, the subducting slab is an important control on the pattern of mantle advection and melt extraction, as well as on compositional and fluid gradients. Modelling studies predict significant mantle wedge effects on back-arc spreading processes. Here we show that various spreading centres in the Lau back-arc basin exhibit enhanced, diminished or normal magma supply, which correlates with distance from the arc volcanic front but not with spreading rate. To explain this correlation we propose that depleted upper-mantle material, generated by melt extraction in the mantle wedge, is overturned and re-introduced beneath the back-arc basin by subduction-induced corner flow. The spreading centres experience enhanced melt delivery near the volcanic front, diminished melting within the overturned depleted mantle farther from the corner and normal melting conditions in undepleted mantle farther away. Our model explains fundamental differences in crustal accretion variables between back-arc and mid-ocean settings.     

High 3He/4He ratios in picritic basalts from Baffin Island and the role of a mixed reservoir in mantle plumes

The high 3He/4He ratio of volcanic rocks thought to be derived from mantle plumes is taken as evidence for the existence of a mantle reservoir that has remained largely undegassed since the Earth's accretion. The helium isotope composition of this reservoir places constraints on the origin of volatiles within the Earth and on the evolution and structure of the Earth's mantle. Here we show that olivine phenocrysts in picritic basalts presumably derived from the proto-Iceland plume at Baffin Island, Canada, have the highest magmatic 3He/4He ratios yet recorded. A strong correlation between 3He/4He and 87Sr/86Sr, 143Nd/144Nd and trace element ratios demonstrate that the 3He-rich end-member is present in basalts that are derived from large-volume melts of depleted upper-mantle rocks. This reservoir is consistent with the recharging of depleted upper-mantle rocks by small volumes of primordial volatile-rich lower-mantle material at a thermal boundary layer between convectively isolated reservoirs. The highest 3He/4He basalts from Hawaii and Iceland plot on the observed mixing trend. This indicates that a 3He-recharged depleted mantle (HRDM) reservoir may be the principal source of high 3He/4He in mantle plumes, and may explain why the helium concentration of the 'plume' component in ocean island basalts is lower than that predicted for a two-layer, steady-state model of mantle structure.     

Evolution of helium isotopes in the Earth's mantle

Degassing of the Earth's mantle through magmatism results in the irreversible loss of helium to space, and high (3)He/(4)He ratios observed in oceanic basalts have been considered the main evidence for a 'primordial' undegassed deep mantle reservoir. Here we present a new global data compilation of ocean island basalts, representing upwelling 'plumes' from the deep mantle, and show that island groups with the highest primordial signal (high (3)He/(4)He ratios) have striking chemical and isotopic similarities to mid-ocean-ridge basalts. We interpret this as indicating a common history of mantle trace element depletion through magmatism. The high (3)He/(4)He in plumes may thus reflect incomplete degassing of the deep Earth during continent and ocean crust formation. We infer that differences between plumes and the upper-mantle source of ocean-ridge basalts reflect isolation of plume sources from the convecting mantle for approximately 1-2 Gyr. An undegassed, primordial reservoir in the mantle would therefore not be required, thus reconciling a long-standing contradiction in mantle dynamics.     

Pressure sensitivity of olivine slip systems and seismic anisotropy of Earth's upper mantle

The mineral olivine dominates the composition of the Earth's upper mantle and hence controls its mechanical behaviour and seismic anisotropy. Experiments at high temperature and moderate pressure, and extensive data on naturally deformed mantle rocks, have led to the conclusion that olivine at upper-mantle conditions deforms essentially by dislocation creep with dominant [100] slip. The resulting crystal preferred orientation has been used extensively to explain the strong seismic anisotropy observed down to 250 km depth. The rapid decrease of anisotropy below this depth has been interpreted as marking the transition from dislocation to diffusion creep in the upper mantle. But new high-pressure experiments suggest that dislocation creep also dominates in the lower part of the upper mantle, but with a different slip direction. Here we show that this high-pressure dislocation creep produces crystal preferred orientations resulting in extremely low seismic anisotropy, consistent with seismological observations below 250 km depth. These results raise new questions about the mechanical state of the lower part of the upper mantle and its coupling with layers both above and below.     

Neon isotopes constrain convection and volatile origin in the Earth's mantle

Identifying the origin of primordial volatiles in the Earth's mantle provides a critical test between models that advocate magma-ocean equilibration with an early massive solar-nebula atmosphere and those that require subduction of volatiles implanted in late accreting material. Here we show that neon isotopes in the convecting mantle, resolved in magmatic CO2 well gases, are consistent with a volatile source related to solar corpuscular irradiation of accreting material. This contrasts with recent results that indicated a solar-nebula origin for neon in mantle plume material, which is thought to be sampling the deep mantle. Neon isotope heterogeneity in different mantle sources suggests that models in which the plume source supplies the convecting mantle with its volatile inventory require revision. Although higher than accepted noble gas concentrations in the convecting mantle may reduce the need for a deep mantle volatile flux, any such flux must be dominated by the neon (and helium) isotopic signature of late accreting material.