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.     

Laboratory models of the thermal evolution of the mantle during rollback subduction

The subduction of oceanic lithosphere plays a key role in plate tectonics, the thermal evolution of the mantle and recycling processes between Earth's interior and surface. Information on mantle flow, thermal conditions and chemical transport in subduction zones come from the geochemistry of arc volcanoes, seismic images and geodynamic models. The majority of this work considers subduction as a two-dimensional process, assuming limited variability in the direction parallel to the trench. In contrast, observationally based models increasingly appeal to three-dimensional flow associated with trench migration and the sinking of oceanic plates with a translational component of motion (rollback). Here we report results from laboratory experiments that reveal fundamental differences in three-dimensional mantle circulation and temperature structure in response to subduction with and without a rollback component. Without rollback motion, flow in the mantle wedge is sluggish, there is no mass flux around the plate and plate edges heat up faster than plate centres. In contrast, during rollback subduction flow is driven around and beneath the sinking plate, velocities increase within the mantle wedge and are focused towards the centre of the plate, and the surface of the plate heats more along the centreline.     

Crystallographic preferred orientation of akimotoite and seismic anisotropy of Tonga slab

The mineral akimotoite, ilmenite-structured MgSiO(3), exists at the bottom of the Earth's mantle transition zone and within the uppermost lower mantle, especially under low-temperature conditions. Akimotoite is thought to be a major constituent of the harzburgite layer of subducting slabs, and the most anisotropic mineral in the mantle transition zone. It has been predicted that if akimotoite crystals are preferentially oriented by plastic deformation, a cold subducted slab would be extremely anisotropic. However, there have been no studies of crystallographic preferred orientations and very few reports of plastic deformation experiments for MgSiO(3) ilmenite. Here we present plastic deformation experiments on polycrystalline akimotoite, which were conducted at confining pressures of 20-22 GPa and temperatures of 1,000-1,300 degrees C. We found a change in crystallographic preferred orientation pattern of akimotoite with temperature, where the c-axis maximum parallel to the compression direction develops at high temperature, whereas the c axes are preferentially oriented parallel to the shear direction or perpendicular to the compression direction at lower temperature. The previously reported difference in compressional-wave seismic anisotropy between the northern and southern segments of the Tonga slab at depths of the mantle transition zone can conceivably be attributed to the difference in the crystallographic preferred orientation pattern of akimotoite at varying temperature within the slab.     

Evidence for mantle metasomatism by hydrous silicic melts derived from subducted oceanic crust

The low concentrations of niobium, tantalum and titanium observed in island-arc basalts are thought to result from modification of the sub-arc mantle by a metasomatic agent, deficient in these elements, that originates from within the subducted oceanic crust. Whether this agent is an hydrous fluid or a silica-rich melt has been discussed using mainly a trace-element approach and related to variable thermal regimes of subduction zones. Melting of basalt in the absence of fluid both requires high temperatures and yields melt compositions unlike those found in most modern or Mesozoic island arcs. Thus, metasomatism by fluids has been thought to be the most common situation. Here, however, we show that the melting of basalt under both H2O-added and low-temperature conditions can yield extremely alkali-rich silicic liquids, the alkali content of which increases with pressure. These liquids are deficient in titanium and in the elements niobium and tantalum and are virtually identical to glasses preserved in mantle xenoliths found in subduction zones and to veins found in exhumed metamorphic terranes of fossil convergent zones. We also found that the interaction between such liquids and mantle olivine produces modal mineralogies that are identical to those observed in metasomatized Alpine-type peridotites. We therefore suggest that mantle metasomatism by slab-derived melt is a more common process than previously thought.     

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.     

Arc-parallel flow in the mantle wedge beneath Costa Rica and Nicaragua

Resolving flow geometry in the mantle wedge is central to understanding the thermal and chemical structure of subduction zones, subducting plate dehydration, and melting that leads to arc volcanism, which can threaten large populations and alter climate through gas and particle emission. Here we show that isotope geochemistry and seismic velocity anisotropy provide strong evidence for trench-parallel flow in the mantle wedge beneath Costa Rica and Nicaragua. This finding contradicts classical models, which predict trench-normal flow owing to the overlying wedge mantle being dragged downwards by the subducting plate. The isotopic signature of central Costa Rican volcanic rocks is not consistent with its derivation from the mantle wedge or eroded fore-arc complexes but instead from seamounts of the Galapagos hotspot track on the subducting Cocos plate. This isotopic signature decreases continuously from central Costa Rica to northwestern Nicaragua. As the age of the isotopic signature beneath Costa Rica can be constrained and its transport distance is known, minimum northwestward flow rates can be estimated (63-190 mm yr(-1)) and are comparable to the magnitude of subducting Cocos plate motion (approximately 85 mm yr(-1)). Trench-parallel flow needs to be taken into account in models evaluating thermal and chemical structure and melt generation in subduction zones.     

Hydrous olivine unable to account for conductivity anomaly at the top of the asthenosphere

The oceanic asthenosphere is observed to have high electrical conductivity, which is highly anisotropic in some locations. In the directions parallel and normal to the plate motion, the conductivity is of the order of 10(-1) and 10(-2) S m(-1), respectively, which cannot be explained by the conductivity of anhydrous olivine. But because hydrogen can be incorporated in olivine at mantle pressures, this observation has been attributed to olivine hydration, which might cause anisotropically high conductivity by proton migration. To examine this hypothesis, here we report the effect of water on electrical conductivity and its anisotropy for hydrogen-doped and undoped olivine at 500-1,500 K and 3 GPa. The hydrous olivine has much higher conductivity and lower activation energy than anhydrous olivine in the investigated temperature range. Nevertheless, extrapolation of the experimental results suggests that conductivity of hydrous olivine at the top of the asthenosphere should be nearly isotropic and only of the order of 10(-2) S m(-1). Our data indicate that the hydration of olivine cannot account for the geophysical observations, which instead may be explained by the presence of partial melt elongated in the direction of plate motion.     

Evolution and diversity of subduction zones controlled by slab width

Subducting slabs provide the main driving force for plate motion and flow in the Earth's mantle, and geodynamic, seismic and geochemical studies offer insight into slab dynamics and subduction-induced flow. Most previous geodynamic studies treat subduction zones as either infinite in trench-parallel extent (that is, two-dimensional) or finite in width but fixed in space. Subduction zones and their associated slabs are, however, limited in lateral extent (250-7,400 km) and their three-dimensional geometry evolves over time. Here we show that slab width controls two first-order features of plate tectonics-the curvature of subduction zones and their tendency to retreat backwards with time. Using three-dimensional numerical simulations of free subduction, we show that trench migration rate is inversely related to slab width and depends on proximity to a lateral slab edge. These results are consistent with retreat velocities observed globally, with maximum velocities (6-16 cm yr(-1)) only observed close to slab edges (<1,200 km), whereas far from edges (>2,000 km) retreat velocities are always slow (<2.0 cm yr(-1)). Models with narrow slabs (< or =1,500 km) retreat fast and develop a curved geometry, concave towards the mantle wedge side. Models with slabs intermediate in width ( approximately 2,000-3,000 km) are sublinear and retreat more slowly. Models with wide slabs (> or =4,000 km) are nearly stationary in the centre and develop a convex geometry, whereas trench retreat increases towards concave-shaped edges. Additionally, we identify periods (5-10 Myr) of slow trench advance at the centre of wide slabs. Such wide-slab behaviour may explain mountain building in the central Andes, as being a consequence of its tectonic setting, far from slab edges.     

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.     

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.     

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.     

Seismic evidence for convection-driven motion of the North American plate

Since the discovery of plate tectonics, the relative importance of driving forces of plate motion has been debated. Resolution of this issue has been hindered by uncertainties in estimates of basal traction, which controls the coupling between lithospheric plates and underlying mantle convection. Hotspot tracks preserve records of past plate motion and provide markers with which the relative motion between a plate's surface and underlying mantle regions may be examined. Here we show that the 115-140-Myr surface expression of the Great Meteor hotspot track in eastern North America is misaligned with respect to its location at 200 km depth, as inferred from plate-reconstruction models and seismic tomographic studies. The misalignment increases with age and is consistent with westward displacement of the base of the plate relative to its surface, at an average rate of 3.8 +/- 1.8 mm yr(-1). Here age-constrained 'piercing points' have enabled direct estimation of relative motion between the surface and underside of a plate. The relative displacement of the base is approximately parallel to seismic fast axes and calculated mantle flow, suggesting that asthenospheric flow may be deforming the lithospheric keel and exerting a driving force on this part of the North American plate.     

Trench-parallel flow and seismic anisotropy in the Mariana and Andean subduction systems

Shear-wave splitting measurements above the mantle wedge of the Mariana and southern Andean subduction zones show trench-parallel seismically fast directions close to the trench and abrupt rotations to trench-perpendicular anisotropy in the back arc. These patterns of seismic anisotropy may be caused by three-dimensional flow associated with along-strike variations in slab geometry. The Mariana and Andean subduction systems are associated with the largest along-strike variations of slab geometry observed on Earth and are ideal for testing the link between slab geometry and solid-state creep processes in the mantle. Here we show, with fully three-dimensional non-newtonian subduction zone models, that the strong curvature of the Mariana slab and the transition to shallow slab dip in the Southern Andes give rise to strong trench-parallel stretching in the warm-arc and warm-back-arc mantle and to abrupt rotations in stretching directions that are accompanied by strong trench-parallel stretching. These models show that the patterns of shear-wave splitting observed in the Mariana and southern Andean systems may be caused by significant three-dimensional flow induced by along-strike variations in slab geometry.     

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.     

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.     

Effects of contact guidance and gravity on L929 cell orientation

Contact guidance and external force field are two important factors that influence cell orientation. Microgroove-like patterns on PDMS substrates were fabricated using soft-lithography, and the effects of gravitational field direction relative to the pattern, on L929 cell orientation were investigated. The majority of cells were aligned along the grooves when gravity was parallel to the groove direction, because of contact guidance. When gravity was perpendicular to the groove direction, cell alignment decreased noticeably. Although contact guidance induced by the pattern played a dominant role in determining cell orientation, gravity also had an important influence when the field direction was perpendicular to the groove direction. No synergistic influence on cell alignment was observed when the field direction was parallel to the groove direction.     

基于横波分裂方法的海南地幔柱研究

利用横波分裂方法,对北京大学、南方科技大学及桂林理工大学在广东和广西布设的宽频带流动台阵以及国家地震台网在广东、广西和海南的固定台站记录的XKS震相进行计算,得到海南及邻近地区的上地幔各向异性参数,分裂延迟时间在0.35～1.90s之间.综合研究区域内前人发表的和本文的横波分裂结果,认为自新生代以来,伴随东部太平洋板块...     

西太平洋俯冲带岩石圈变形研究

全球海洋岩石圈的最大弯曲与大地震发生在俯冲带。当弯曲应力超过岩石承受范围，就会产生正断层和地震，海水沿着正断层进入上地幔并发生蛇纹石化，引发浅源地震并可能造成灾难性海啸。选取西太平洋最具代表性的日本、伊豆-小笠原、马里亚纳和雅浦俯冲带以及汤加-克马德克俯冲带，归纳近些年的地球物理观测及地球动力学模拟的结果，对比分析了不同俯冲带挠曲正断层的分布特征，并探讨了俯冲板块变形与地震之间的相关性,以揭示俯冲板片弯曲变形及相应的正断层与潜在板块水化特征。     

Destruction geodynamics of the North China craton and its Paleoproterozoic plate tectonics

Much attention has been paid in the last two decades to the physical and chemical processes as well as temporal-spatial variations of the lithospheric mantle beneath the North China Craton. In order to provide insights into the geodynamics of this variation, it is necessary to thoroughly study the state and structure of the lithospheric crust and mantle of the North China Craton and its adjacent regions as an integrated unit. Based on the velocity structure of the crust and upper mantle constrained from seismological studies, this paper presents various available geophysical results regarding the lithosphere thickness, the nature of crust-mantle boundary, the upper mantle structure and deformation characteristics as well as their tectonic features and evolution systematics. Combined with the obtained data from petrology and geochemistry, a mantle flow model is proposed for the tectonic evolution of the North China Craton during the Mesozoic-Cenozoic. We suggest that subduction of the Pacific plate made the mantle underneath the eastern Asian continent unstable and able to flow faster. Such a regional mantle flow system would cause an elevation of melt/fluid content in the upper mantle of the North China Craton and the lithospheric softening, which, subsequently resulted in destruction of the North China Craton in different ways of delamination and thermal erosion in Yanshan, Taihang Mountains and the Tan-Lu Fault zone. Multiple lines of evidence recorded in the crust of the North China Craton, such as the amalgamation of the Archean eastern and western blocks, the subduction of Paleo-oceanic crust and Paleo-continental residue, indicate that the Earth in the Paleoproterozoic had already evolved into the plate tectonic system similar to the present plate tectonics.     

Measuring the onset of locking in the Peru-Chile trench with GPS and acoustic measurements

The subduction zone off the west coast of South America marks the convergence of the oceanic Nazca plate and the continental South America plate. Nazca-South America convergence over the past 23 million years has created the 6-km-deep Peru-Chile trench, 150 km offshore. High pressure between the plates creates a locked zone, leading to deformation of the overriding plate. The surface area of this locked zone is thought to control the magnitude of co-seismic release and is limited by pressure, temperature, sediment type and fluid content. Here we present seafloor deformation data from the submerged South America plate obtained from a combination of Global Positioning System (GPS) receivers and acoustic transponders. We estimate that the measured horizontal surface motion perpendicular to the trench is consistent with a model having no slip along the thrust fault between 2 and 40 km depth. A tsunami in 1996, 200 km north of our site, was interpreted as being the result of an anomalously shallow interplate earthquake. Seismic coupling at shallow depths, such as we observe, may explain why co-seismic events in the Peruvian subduction zone create large tsunamis.