Dislocation creep in MgSiO3 perovskite at conditions of the Earth's uppermost lower mantle

Seismic anisotropy provides an important observational constraint on flow in the Earth's deep interior. The quantitative interpretation of anisotropy, however, requires knowledge of the slip geometry of the constitutive minerals that are responsible for producing rock fabrics. The Earth's lower mantle is mostly composed of (Mg, Fe)SiO3 perovskite, but as MgSiO3 perovskite is not stable at high temperature under ambient pressure, it has not been possible to investigate its mechanical behaviour with conventional laboratory deformation experiments. To overcome this limitation, several attempts were made to infer the mechanical properties of MgSiO3 perovskite on the basis of analogue materials. But perovskites do not constitute an analogue series for plastic deformation, and therefore the direct investigation of MgSiO3 perovskite is necessary. Here we have taken advantage of recent advances in experimental high-pressure rheology to perform deformation experiments on coarse-grained MgSiO3 polycrystals under pressure and temperature conditions of the uppermost lower mantle. We show that X-ray peak broadening measurements developed in metallurgy can be adapted to low-symmetry minerals to identify the elementary deformation mechanisms activated under these conditions. We conclude that, under uppermost lower-mantle conditions, MgSiO3 perovskite deforms by dislocation creep and may therefore contribute to producing seismic anisotropy in rocks at such depths.     

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.     

Seismic evidence of negligible water carried below 400-km depth in subducting lithosphere

Strong evidence exists that water is carried from the surface into the upper mantle by hydrous minerals in the uppermost 10-12?km of subducting lithosphere, and more water may be added as the lithosphere bends and goes downwards. Significant amounts of that water are released as the lithosphere heats up, triggering earthquakes and fluxing arc volcanism. In addition, there is experimental evidence for high solubility of water in olivine, the most abundant mineral in the upper mantle, for even higher solubility in olivine's high-pressure polymorphs, wadsleyite and ringwoodite, and for the existence of dense hydrous magnesium silicates that potentially could carry water well into the lower mantle (deeper than 1,000?km). Here we compare experimental and seismic evidence to test whether patterns of seismicity and the stabilities of these potentially relevant hydrous phases are consistent with a wet lithosphere. We show that there is nearly a one-to-one correlation between dehydration of minerals and seismicity at depths less than about 250?km, and conclude that the dehydration of minerals is the trigger of instability that leads to seismicity. At greater depths, however, we find no correlation between occurrences of earthquakes and depths where breakdown of hydrous phases is expected. Lastly, we note that there is compelling evidence for the existence of metastable olivine (which, if present, can explain the distribution of deep-focus earthquakes) west of and within the subducting Tonga slab and also in three other subduction zones, despite metastable olivine being incompatible with even extremely small amounts of water (of the order of 100?p.p.m. by weight). We conclude that subducting slabs are essentially dry at depths below 400?km and thus do not provide a pathway for significant amounts of water to enter the mantle transition zone or the lower mantle.     

The depth distribution of azimuthal anisotropy in the continental upper mantle

The most likely cause of seismic anisotropy in the Earth's upper mantle is the lattice preferred orientation of anisotropic minerals such as olivine. Its presence reflects dynamic processes related to formation of the lithosphere as well as to present-day tectonic motions. A powerful tool for detecting and characterizing upper-mantle anisotropy is the analysis of shear-wave splitting measurements. Because of the poor vertical resolution afforded by this type of data, however, it has remained controversial whether the splitting has a lithospheric origin that is 'frozen-in' at the time of formation of the craton, or whether the anisotropy originates primarily in the asthenosphere, and is induced by shear owing to present-day absolute plate motions. In addition, predictions from surface-wave-derived models are largely incompatible with shear-wave splitting observations. Here we show that this disagreement can be resolved by simultaneously inverting surface waveforms and shear-wave splitting data. We present evidence for the presence of two layers of anisotropy with different fast-axis orientations in the cratonic part of the North American upper mantle. At asthenospheric depths (200-400 km) the fast axis is sub-parallel to the absolute plate motion, confirming the presence of shear related to current tectonic processes, whereas in the lithosphere (80-200 km), the orientation is significantly more northerly. In the western, tectonically active, part of North America, the fast-axis direction is consistent with the absolute plate motion throughout the depth range considered, in agreement with a much thinner lithosphere.     

Natural examples of olivine lattice preferred orientation patterns with a flow-normal a-axis maximum

Tectonic plate motion is thought to cause solid-state plastic flow within the underlying upper mantle and accordingly lead to the development of a lattice preferred orientation of the constituent olivine crystals. The mechanical anisotropy that results from such preferred orientation typically produces a direction of maximum seismic wave velocity parallel to the plate motion direction. This has been explained by the existence of an olivine preferred orientation with an 'a-axis' maximum parallel to the induced mantle flow direction. In subduction zones, however, the olivine a axes have been inferred to be arranged roughly perpendicular to plate motion, which has usually been ascribed to localized complex mantle flow patterns. Recent experimental work suggests an alternative explanation: under conditions of high water activity, a 'B-type' olivine preferred orientation may form, with the a-axis maximum perpendicular to the flow direction. Natural examples of such B-type preferred orientation are, however, almost entirely unknown. Here we document widespread B-type olivine preferred orientation patterns from a subduction-type metamorphic belt in southwest Japan and show that these patterns developed in the presence of water. Our discovery implies that mantle flow above subduction zones may be much simpler than has generally been thought.     

Sound velocities of majorite garnet and the composition of the mantle transition region

The composition of the mantle transition region, characterized by anomalous seismic-wave velocity and density changes at depths of approximately 400 to 700 km, has remained controversial. Some have proposed that the mantle transition region has an olivine-rich 'pyrolite' composition, whereas others have inferred that it is characterized by pyroxene- and garnet-rich compositions ('piclogite'), because the sound velocities in pyrolite estimated from laboratory data are substantially higher than those seismologically observed. Although the velocities of the olivine polymorphs at these pressures (wadsleyite and ringwoodite) have been well documented, those of majorite (another significant high-pressure phase in the mantle transition region) with realistic mantle compositions have never been measured. Here we use combined in situ X-ray and ultrasonic measurements under the pressure and temperature conditions of the mantle transition region to show that majorite in a pyrolite composition has sound velocities substantially lower than those of earlier estimates, owing to strong nonlinear decreases at high temperature, particularly for shear-wave velocity. We found that pyrolite yields seismic velocities more consistent with typical seismological models than those of piclogite in the upper to middle parts of the region, except for the potentially larger velocity jumps in pyrolite relative to those observed at a depth of 410 km. In contrast, both of these compositions lead to significantly low shear-wave velocities in the lower part of the region, suggesting possible subadiabatic temperatures or the existence of a layer of harzburgite-rich material supplied by the subducted slabs stagnant at these depths.     

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.     

Redox freezing and melting in the Earth's deep mantle resulting from carbon-iron redox coupling

Very low seismic velocity anomalies in the Earth's mantle may reflect small amounts of melt present in the peridotite matrix, and the onset of melting in the Earth's upper mantle is likely to be triggered by the presence of small amounts of carbonate. Such carbonates stem from subducted oceanic lithosphere in part buried to depths below the 660-kilometre discontinuity and remixed into the mantle. Here we demonstrate that carbonate-induced melting may occur in deeply subducted lithosphere at near-adiabatic temperatures in the Earth's transition zone and lower mantle. We show experimentally that these carbonatite melts are unstable when infiltrating ambient mantle and are reduced to immobile diamond when recycled at depths greater than ～250?kilometres, where mantle redox conditions are determined by the presence of an (Fe,Ni) metal phase. This 'redox freezing' process leads to diamond-enriched mantle domains in which the Fe(0), resulting from Fe(2+) disproportionation in perovskites and garnet, is consumed but the Fe(3+) preserved. When such carbon-enriched mantle heterogeneities become part of the upwelling mantle, diamond will inevitably react with the Fe(3+) leading to true carbonatite redox melting at ～660 and ～250 kilometres depth to form deep-seated melts in the Earth's mantle.     

Resistance to mantle flow inferred from the electromagnetic strike of the Australian upper mantle

Seismic anisotropy is thought to result from the strain-induced lattice-preferred orientation of mantle minerals, especially olivine, owing to shear waves propagating faster along the a-axis of olivine crystals than along the other axes. This anisotropy results in birefringence, or 'shear-wave splitting', which has been investigated in numerous studies. Although olivine is also anisotropic with respect to electrical conductivity (with the a-axis being most conductive), few studies of the electrical anisotropy of the upper mantle have been undertaken, and these have been limited to relatively shallow depths in the lithospheric upper mantle. Theoretical models of mantle flow have been used to infer that, for progressive simple shear imparted by the motion of an overriding tectonic plate, the a-axes of olivine crystals should align themselves parallel to the direction of plate motion. Here, however, we show that a significant discrepancy exists between the electromagnetic strike of the mantle below Australia and the direction of present-day absolute plate motion. We infer from this discrepancy that the a-axes of olivine crystals are not aligned with the direction of the present-day plate motion of Australia, indicating resistance to deformation of the mantle by plate motion.     

Global anisotropy and the thickness of continents

For decades there has been a vigorous debate about the depth extent of continental roots. The analysis of heat-flow, mantle-xenolith and electrical-conductivity data all indicate that the coherent, conductive part of continental roots (the 'tectosphere') is at most 200-250 km thick. Some global seismic tomographic models agree with this estimate, but others suggest that a much thicker zone of high velocities lies beneath continental shields, reaching a depth of at least 400 km. Here we show that this disagreement can be reconciled by taking into account seismic anisotropy. We show that significant radial anisotropy, with horizontally polarized shear waves travelling faster than those that are vertically polarized, is present under most cratons in the depth range 250-400 km--similar to that found under ocean basins at shallower depths of 80-250 km. We propose that, in both cases, the anisotropy is related to shear in a low-viscosity asthenospheric channel, located at different depths under continents and oceans. The seismically defined 'tectosphere' is then at most 200-250 km thick under old continents. The 'Lehmann discontinuity', observed mostly under continents at about 200-250 km, and the 'Gutenberg discontinuity', observed under oceans at depths of about 60-80 km, may both be associated with the bottom of the lithosphere, marking a transition to flow-induced asthenospheric anisotropy.     

山东地区地壳及上地幔结构研究

利用布设在山东省境内的宽频带流动地震观测台阵和国家地震局固定地震观测台站记录的地震数据,应用接收函数和SKS波分裂方法,研究山东地区的地壳与上地幔结构,得到该区域的地壳厚度、地壳平均P波与S波的波速比以及SKS波分裂延迟的分布情况。结果表明,山东地区地壳厚度范围为28～39km;胶南隆起的北段和南段以及鲁西隆起北侧济阳凹陷的地壳厚度小于32km,鲁西隆起下方的地壳比较厚。研究区P波与S波的波速比主要分布在1.67～1.94之间,鲁西隆起西南部和胶南隆起北段该比值小于1.75,可能是由中上地壳增厚以及下地壳减薄和拆沉造成。鲁西隆起南北P波与S波的波速比差异反映地壳活动的差异。地幔物质的各向异性显示,山东地区西部的地壳减薄和拆沉可能仍在进行。     

Striped iron zoning of olivine induced by dislocation creep in deformed peridotites.

Deformation of solid materials affects not only their microstructures, but also their microchemistries. Although chemical unmixing of initially homogeneous multicomponent solids is known to occur during deformation by diffusion creep, there has been no report on their chemical zoning due to deformation by dislocation creep, in either natural samples or laboratory experiments. Here we report striped iron zoning of olivine ((Mg,Fe)2SiO4) in deformed peridotites, where the iron concentration increases at subgrain boundaries composed of edge dislocations. We infer that this zoning is probably formed by alignment of edge dislocations dragging a so-called Cottrell 'atmosphere' of solute atoms (iron in this case) into subgrain boundaries during deformation of the olivine by dislocation creep. We have found that the iron zoning does not develop in laboratory experiments of high strain rates where dislocations move too fast to drag the Cottrell atmosphere. This phenomenon might have important implications for the generation of deep-focus earthquakes, as transformation of olivine to high-pressure phases preferentially occurs in high-iron regions, and therefore along subgrain boundaries which would be preferentially aligned in plastically deformed mantle peridotites.     

Inferring nonlinear mantle rheology from the shape of the Hawaiian swell

The convective circulation generated within the Earth's mantle by buoyancy forces of thermal and compositional origin is intimately controlled by the rheology of the rocks that compose it. These can deform either by the diffusion of point defects (diffusion creep, with a linear relationship between strain rate and stress) or by the movement of intracrystalline dislocations (nonlinear dislocation creep). However, there is still no reliable map showing where in the mantle each of these mechanisms is dominant, and so it is important to identify regions where the operative mechanism can be inferred directly from surface geophysical observations. Here we identify a new observable quantity--the rate of downstream decay of the anomalous seafloor topography (swell) produced by a mantle plume--which depends only on the value of the exponent in the strain rate versus stress relationship that defines the difference between diffusion and dislocation creep. Comparison of the Hawaiian swell topography with the predictions of a simple fluid mechanical model shows that the swell shape is poorly explained by diffusion creep, and requires a dislocation creep rheology. The rheology predicted by the model is reasonably consistent with laboratory deformation data for both olivine and clinopyroxene, suggesting that the source of Hawaiian lavas could contain either or both of these components.     

Metal saturation in the upper mantle

The oxygen fugacity f(O2)of the Earth's mantle is one of the fundamental variables in mantle petrology. Through ferric-ferrous iron and carbon-hydrogen-oxygen equilibria, f(O2) influences the pressure-temperature positions of mantle solidi and compositions of small-degree mantle melts. Among other parameters, f(O2) affects the water storage capacity and rheology of the mantle. The uppermost mantle, as represented by samples and partial melts, is sufficiently oxidized to sustain volatiles, such as H2O and CO2, as well as carbonatitic melts, but it is not known whether the shallow mantle is representative of the entire upper mantle. Using high-pressure experiments, we show here that large parts of the asthenosphere are likely to be metal-saturated. We found that pyroxene and garnet synthesized at >7 GPa in equilibrium with metallic Fe can incorporate sufficient ferric iron that the mantle at >250 km depth is so reduced that an (Fe,Ni)-metal phase may be stable. Our results indicate that the oxidized nature of the upper mantle can no longer be regarded as being representative for the Earth's upper mantle as a whole and instead that oxidation is a shallow phenomenon restricted to an upper veneer only about 250 km in thickness.     

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.     

Magma-assisted rifting in Ethiopia

The rifting of continents and evolution of ocean basins is a fundamental component of plate tectonics, yet the process of continental break-up remains controversial. Plate driving forces have been estimated to be as much as an order of magnitude smaller than those required to rupture thick continental lithosphere. However, Buck has proposed that lithospheric heating by mantle upwelling and related magma production could promote lithospheric rupture at much lower stresses. Such models of mechanical versus magma-assisted extension can be tested, because they predict different temporal and spatial patterns of crustal and upper-mantle structure. Changes in plate deformation produce strain-enhanced crystal alignment and increased melt production within the upper mantle, both of which can cause seismic anisotropy. The Northern Ethiopian Rift is an ideal place to test break-up models because it formed in cratonic lithosphere with minor far-field plate stresses. Here we present evidence of seismic anisotropy in the upper mantle of this rift zone using observations of shear-wave splitting. Our observations, together with recent geological data, indicate a strong component of melt-induced anisotropy with only minor crustal stretching, supporting the magma-assisted rifting model in this area of initially cold, thick continental lithosphere.     

Lithospheric thermal-rheological structures of the continental margin in the northern South China Sea

Thermal structures of three deep seismic profiles in the continental margin in the northern South China Sea are calculated, their "thermal" lithospheric thicknesses are evaluated based on the basalt dry solidus, and their rheological structures are evaluated with linear frictional failure criterion and power-law creep equation. "Thermal" lithosphere is about 90 km in thickness in shelf area, and thins toward the slope, lowers to 60-65 km in the lower slope, ocean crust and Xisha Trough. In the mid-west of the studied area, the lithospheric rheological structure in shelf area and Xisha Islands is of four layers: brittle, ductile, brittle and ductile. Because of uprising of heat mantle and thinning of crust and lithosphere in Xisha Trough, the bottom of the upper brittle layer is only buried at 16 km. In the eastern area, the bottom of the upper brittle layer in the north is buried at 20 km or so, while in lower slope and ocean crust, the rheological structure is of two layers of brittle and ductile, and crust and uppermost mantle form one whole brittle layer whose bottom is buried at 30-32 km. Analyses show that the characteristics of rheological structure accord with the seismic result observed. The character of rheological stratification implies that before the extension of the continent margin, there likely was a ductile layer in mid-lower crust. The influence of the existence of ductile layer to the evolution of the continent margin and the different extensions of ductile layer and brittle layer should not be overlooked. Its thickness, depth and extent in influencing continent margin's extension and evolution should be well evaluated in building a dynamic model for the area.     

中国东北上地幔各向异性研究

利用中国数字地震台网记录的地震波形,观测其穿过地球外核和地幔的SKS波分裂特征,确定了快波S波的振动方向和快慢波时间延迟。这种经过地幔传播的SKS震相的横波分裂主要来源于上地幔的各向异性。可以用应变引起的上地幔中矿物质的结晶优势排列来解释。结果表明,对于所分析的10个观测台站均发现了明显的横波分裂现象。时间延迟在0.4～1.8 s。在东部地区,快波振动方向与中国东北西太平洋俯冲带的板块相对运动方向基本一致,说明俯冲的太平洋板块地幔流对各向异性有重要意义。在西部地区,快波振动方向和时间延迟都具有较大差异,说明其总体形变的起因是复杂的,各向异性受到地幔流和历史构造运动的双重影响。     

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.