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.     

Development of anisotropic structure in the Earth's lower mantle by solid-state convection

Seismological observations reveal highly anisotropic patches at the bottom of the Earth's lower mantle, whereas the bulk of the mantle has been observed to be largely isotropic. These patches have been interpreted to correspond to areas where subduction has taken place in the past or to areas where mantle plumes are upwelling, but the underlying cause for the anisotropy is unknown-both shape-preferred orientation of elastically heterogeneous materials and lattice-preferred orientation of a homogeneous material have been proposed. Both of these mechanisms imply that large-strain deformation occurs within the anisotropic regions, but the geodynamic implications of the mechanisms differ. Shape-preferred orientation would imply the presence of large elastic (and hence chemical) heterogeneity whereas lattice-preferred orientation requires deformation at high stresses. Here we show, on the basis of numerical modelling incorporating mineral physics of elasticity and development of lattice-preferred orientation, that slab deformation in the deep lower mantle can account for the presence of strong anisotropy in the circum-Pacific region. In this model-where development of the mineral fabric (the alignment of mineral grains) is caused solely by solid-state deformation of chemically homogeneous mantle material-anisotropy is caused by large-strain deformation at high stresses, due to the collision of subducted slabs with the core-mantle boundary.     

Texturing of the Earth's inner core by Maxwell stresses.

Elastic anisotropy in the Earth's inner core has been attributed to a preferred lattice orientation, which may be acquired during solidification of the inner core or developed subsequent to solidification as a result of plastic deformation. But solidification texturing alone cannot explain the observed depth dependence of anisotropy, and previous suggestions for possible deformation processes have all relied on radial flow, which is inhibited by thermal and chemical stratification. Here we investigate the development of anisotropy as the inner core deforms plastically under the influence of electromagnetic (Maxwell) shear stresses. We estimate the flow caused by a representative magnetic field using polycrystal plasticity simulations for epsilon-iron, where the imposed deformation is accommodated by basal and prismatic slip. We find that individual grains in an initially random polycrystal become preferentially oriented with their c axes parallel to the equatorial plane. This pattern is accentuated if deformation is accompanied by recrystallization. Using the single-crystal elastic properties of epsilon-iron at core pressure and temperature, we average over the simulated orientation distribution to obtain a pattern of elastic anisotropy which is similar to that observed seismologically.     

Coesite eclogite facies metamorphic tectonite: Evidence for ultrahigh pressure deformation

An ultrahigh pressure ductile shear mélange crops out on the beach of Yangkou Bay near Qingdao City. The mélange is composed of weakly deformed blocks in a highly ductilly flow mylonites. Ultrahigh pressure metamorphic (UHPM) tectonite includes strongly deformed eclogite and mylonitized eclogite. Coesite occurs in the tectonite as both interstitial mineral and inclusions in garnet and omphacite, indicating that the deformation took place in the stability field of coesite (800—850℃, >30 GPa) in the upper mantle. Coesite is rounded or short prismatic grains with undulatory extinction, and often fractured, suggesting brittle deformation. Garnet is also characterized by brittle fractures and sometimes necked and slightly elongated. Omphacite is elongated, with long axis preferred orientation. Undulatory extinction, subgrains and dynamically recrystallized grains suggest plastic flow of omphacite. Ultrahigh pressure metamorphic tectonite was probably formed in the ductile shear zone during the early stage of exhumation of the ultrahigh pressure metamorphic rocks. Its kinematic indicators point to the transport direction of the UHPM slab during the early stage of exhumation.     

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.     

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.     

Boundary-layer mantle flow under the Dead Sea transform fault inferred from seismic anisotropy

Lithospheric-scale transform faults play an important role in the dynamics of global plate motion. Near-surface deformation fields for such faults are relatively well documented by satellite geodesy, strain measurements and earthquake source studies, and deeper crustal structure has been imaged by seismic profiling. Relatively little is known, however, about deformation taking place in the subcrustal lithosphere--that is, the width and depth of the region associated with the deformation, the transition between deformed and undeformed lithosphere and the interaction between lithospheric and asthenospheric mantle flow at the plate boundary. Here we present evidence for a narrow, approximately 20-km-wide, subcrustal anisotropic zone of fault-parallel mineral alignment beneath the Dead Sea transform, obtained from an inversion of shear-wave splitting observations along a dense receiver profile. The geometry of this zone and the contrast between distinct anisotropic domains suggest subhorizontal mantle flow within a vertical boundary layer that extends through the entire lithosphere and accommodates the transform motion between the African and Arabian plates within this relatively narrow zone.     

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.     

The plastic deformation of iron at pressures of the Earth's inner core

Soon after the discovery of seismic anisotropy in the Earth's inner core, it was suggested that crystal alignment attained during deformation might be responsible. Since then, several other mechanisms have been proposed to account for the observed anisotropy, but the lack of deformation experiments performed at the extreme pressure conditions corresponding to the solid inner core has limited our ability to determine which deformation mechanism applies to this region of the Earth. Here we determine directly the elastic and plastic deformation mechanism of iron at pressures of the Earth's core, from synchrotron X-ray diffraction measurements of iron, under imposed axial stress, in diamond-anvil cells. The epsilon-iron (hexagonally close packed) crystals display strong preferred orientation, with c-axes parallel to the axis of the diamond-anvil cell. Polycrystal plasticity theory predicts an alignment of c-axes parallel to the compression direction as a result of basal slip, if basal slip is either the primary or a secondary slip system. The experiments provide direct observations of deformation mechanisms that occur in the Earth's inner core, and introduce a method for investigating, within the laboratory, the rheology of materials at extreme pressures.     

Hydroxyl induced eclogite fabric and deformation mechanism

Eclogites from orogens often show strong plastic deformation and high hydroxyl content. We have studied the correlation between crystallographic preferred orientations of garnet and omphacite from natural eclogites with their hydroxyl contents using the electron back-scattered diffraction technique. The results show: 1) Omphacite has typical L-type or SL-type crystrallographic preferred orientations, that is, [001] is distributed in a girdle in the foliation plane with a maximum parallel to lineation; (010) is distributed in a girdle normal to tbe lineation with a maximum parallel to the foliation plane, suggesting a shear dominant deformation regime. Omphacite fabrics do not vary significantly with hydroxyl content, although the hydrous component may cause lower flow strength. 2) Hydroxyl can influence significantly flow properties of garnet in eclogite. Garnets behave as rigid bodies under low temperature and dry conditions. Grain boundary processes will dominate the deformation and lower the flow strength of garnet under high water fugacity conditions. Garnets show no crystallographic preferred orientation in both cases.These results may have important implications for a better understanding of deformation mechanisms and associated fluid activities during deep subduction and exhumation processes.     

Deformation of the lowermost mantle from seismic anisotropy

The lowermost part of the Earth's mantle-known as D″-shows significant seismic anisotropy, the variation of seismic wave speed with direction. This is probably due to deformation-induced alignment of MgSiO(3)-post-perovskite (ppv), which is believed to be the main mineral phase present in the region. If this is the case, then previous measurements of D″ anisotropy, which are generally made in one direction only, are insufficient to distinguish candidate mechanisms of slip in ppv because the mineral is orthorhombic. Here we measure anisotropy in D″ beneath North and Central America, where material from subducting oceanic slabs impinges on the core-mantle boundary, using shallow as well as deep earthquakes to increase the azimuthal coverage in D″. We make more than 700 individual measurements of shear wave splitting in D″ in three regions from two different azimuths in each case. We show that the previously assumed case of vertical transverse isotropy (where wave speed shows no azimuthal variation) is not possible, and that more complicated mechanisms must be involved. We test the fit of different MgSiO(3)-ppv deformation mechanisms to our results and find that shear on (001) is most consistent with observations and the expected shear above the core-mantle boundary beneath subduction zones. With new models of mantle flow, or improved experimental determination of the dominant ppv slip systems, this method will allow us to map deformation at the core-mantle boundary and link processes in D″, such as plume initiation, to the rest of the mantle.     

Anisotropy of Earth's D' layer and stacking faults in the MgSiO3 post-perovskite phase

The post-perovskite phase of (Mg,Fe)SiO3 is believed to be the main mineral phase of the Earth's lowermost mantle (the D' layer). Its properties explain numerous geophysical observations associated with this layer-for example, the D' discontinuity, its topography and seismic anisotropy within the layer. Here we use a novel simulation technique, first-principles metadynamics, to identify a family of low-energy polytypic stacking-fault structures intermediate between the perovskite and post-perovskite phases. Metadynamics trajectories identify plane sliding involving the formation of stacking faults as the most favourable pathway for the phase transition, and as a likely mechanism for plastic deformation of perovskite and post-perovskite. In particular, the predicted slip planes are {010} for perovskite (consistent with experiment) and {110} for post-perovskite (in contrast to the previously expected {010} slip planes). Dominant slip planes define the lattice preferred orientation and elastic anisotropy of the texture. The {110} slip planes in post-perovskite require a much smaller degree of lattice preferred orientation to explain geophysical observations of shear-wave anisotropy in the D' layer.     

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.     

The elastic constants of MgSiO3 perovskite at pressures and temperatures of the Earth's mantle.

The temperature anomalies in the Earth's mantle associated with thermal convection can be inferred from seismic tomography, provided that the elastic properties of mantle minerals are known as a function of temperature at mantle pressures. At present, however, such information is difficult to obtain directly through laboratory experiments. We have therefore taken advantage of recent advances in computer technology, and have performed finite-temperature ab initio molecular dynamics simulations of the elastic properties of MgSiO3 perovskite, the major mineral of the lower mantle, at relevant thermodynamic conditions. When combined with the results from tomographic images of the mantle, our results indicate that the lower mantle is either significantly anelastic or compositionally heterogeneous on large scales. We found the temperature contrast between the coldest and hottest regions of the mantle, at a given depth, to be about 800 K at 1,000 km, 1,500 K at 2,000 km, and possibly over 2,000 K at the core-mantle boundary.     

Efficacy of the post-perovskite phase as an explanation for lowermost-mantle seismic properties

Constraining the chemical, rheological and electromagnetic properties of the lowermost mantle (D') is important to understand the formation and dynamics of the Earth's mantle and core. To explain the origin of the variety of characteristics of this layer observed with seismology, a number of theories have been proposed, including core-mantle interaction, the presence of remnants of subducted material and that D' is the site of a mineral phase transformation. This final possibility has been rejuvenated by recent evidence for a phase change in MgSiO3 perovskite (thought to be the most prevalent phase in the lower mantle) at near core-mantle boundary temperature and pressure conditions. Here we explore the efficacy of this 'post-perovskite' phase to explain the seismic properties of the lowermost mantle through coupled ab initio and seismic modelling of perovskite and post-perovskite polymorphs of MgSiO3, performed at lowermost-mantle temperatures and pressures. We show that a post-perovskite model can explain the topography and location of the D' discontinuity, apparent differences in compressional- and shear-wave models and the observation of a deeper, weaker discontinuity. Furthermore, our calculations show that the regional variations in lower-mantle shear-wave anisotropy are consistent with the proposed phase change in MgSiO3 perovskite.     

High-pressure polymorphs of olivine and the 660-km seismic discontinuity

It had long been accepted that the 400-km seismic discontinuity in the Earth's mantle results from the phase transition of (Mg,Fe)2-SiO4-olivine to its high-pressure polymorph beta-spinel (wadsleyite), and that the 660-km discontinuity results from the breakdown of the higher-pressure polymorph gamma-spinel (ringwoodite) to MgSiO3-perovskite and (Mg,Fe)O-magnesiowüstite. An in situ multi-anvil-press X-ray study indicated, however, that the phase boundary of the latter transition occurs at pressures 2 GPa lower than had been found in earlier studies using multi-anvil recovery experiments and laser-heated diamond-anvil cells. Such a lower-pressure phase boundary would be irreconcilable with the accuracy of seismic measurements of the 660-km discontinuity, and would thus require a mineral composition of the mantle that is significantly different from what is currently thought. Here, however, we present measurements made with a laser-heated diamond-anvil cell which indicate that gamma-Mg2SiO4 is stable up to pressure and temperature conditions equivalent to 660-km depth in the Earth's mantle (24 GPa and 1,900 K) and then breaks down into MgSiO3-perovskite and MgO (periclase). We paid special attention to pressure accuracy and thermal pressure in our experiments, and to ensuring that our experiments were performed under nearly hydrostatic, inert pressure conditions using a variety of heating methods. We infer that these factors are responsible for the different results obtained in our experiments compared to the in situ multi-anvil-press study.     

各向异性油藏中五点井网的变形作用

针对各向异性油藏中水驱开发的研究主要集中于各向异性方向垂直或平行于井网单元中任意两口井的连线方向等特殊的情形,而对于其他布井方向的研究相对缺乏的这种情况,采用坐标变换方法将各向异性油藏转化为各向同性油藏,然后利用流线方法计算正五点井网的开发指标,分析了各向异性程度和方向对水驱开发效果的影响。结果表明,在各向异性油藏中,各向异性程度越大,井网的变形程度也越大;当井网单元中任意两口井的连线与渗透率主轴方向既不平行又不垂直时,井网开发效果对各向异性程度的变化反应十分敏感,井网单元可能会发生的破坏和重组现象,这种破坏和重组现象是各向异性方向和大小两者共同作用的结果,从而造成了井网开发效果的不确定性。为了有效避免可能发生的井网的破坏和重组现象,布井时最好选择注水井生产井方向或者注水井注水井方向平行或者垂直于渗透率主轴方向。     

Atom-by-atom substitution of Mn in GaAs and visualization of their hole-mediated interactions

The discovery of ferromagnetism in Mn-doped GaAs has ignited interest in the development of semiconductor technologies based on electron spin and has led to several proof-of-concept spintronic devices. A major hurdle for realistic applications of Ga(1-x)Mn(x)As, or other dilute magnetic semiconductors, remains that their ferromagnetic transition temperature is below room temperature. Enhancing ferromagnetism in semiconductors requires us to understand the mechanisms for interaction between magnetic dopants, such as Mn, and identify the circumstances in which ferromagnetic interactions are maximized. Here we describe an atom-by-atom substitution technique using a scanning tunnelling microscope (STM) and apply it to perform a controlled study at the atomic scale of the interactions between isolated Mn acceptors, which are mediated by holes in GaAs. High-resolution STM measurements are used to visualize the GaAs electronic states that participate in the Mn-Mn interaction and to quantify the interaction strengths as a function of relative position and orientation. Our experimental findings, which can be explained using tight-binding model calculations, reveal a strong dependence of ferromagnetic interaction on crystallographic orientation. This anisotropic interaction can potentially be exploited by growing oriented Ga(1-x)Mn(x)As structures to enhance the ferromagnetic transition temperature beyond that achieved in randomly doped samples.     

Structure characterization of melt drawn polyethylene ultrathin films

Highly oriented ultrathin polyethylene (PE) films were prepared by a melt-draw technique. Transmission electron microscopy study on the obtained ultrathin films indicates that the melt drawn PE thin films consist of highly oriented edge-on lamellae aligned perpendicular to the drawing direction. Electron diffraction confirms that the PE chains in crystal-line phase are highly oriented in film plane along the drawing direction, while only a random orientation of the crystallographic a-and b-axes can be described through electron diffraction. The IR results on the melt-drawn ultrathin PE films demonstrate that the PE molecular chains in both crystalline and amorphous phases of the melt drawn thin films are well oriented along the drawing direction. Moreover, the IR results indicate that in the crystalline phase, the crystallographic b-axis tend to lie in the thin film plane, while in the amorphous phase, the skeleton plane of some local chains prefers to parallel the film plane. The alignment of b-axis of PE crystals in ultrathin films originates from the fact that the b-axis is the fastest growth direction of the PE crystals.