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.     

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

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

Continental collision slowing due to viscous mantle lithosphere rather than topography

Because the inertia of tectonic plates is negligible, plate velocities result from the balance of forces acting at plate margins and along their base. Observations of past plate motion derived from marine magnetic anomalies provide evidence of how continental deformation may contribute to plate driving forces. A decrease in convergence rate at the inception of continental collision is expected because of the greater buoyancy of continental than oceanic lithosphere, but post-collisional rates are less well understood. Slowing of convergence has generally been attributed to the development of high topography that further resists convergent motion; however, the role of deforming continental mantle lithosphere on plate motions has not previously been considered. Here I show that the rate of India's penetration into Eurasia has decreased exponentially since their collision. The exponential decrease in convergence rate suggests that contractional strain across Tibet has been constant throughout the collision at a rate of 7.03?×?10(-16)?s(-1), which matches the current rate. A constant bulk strain rate of the orogen suggests that convergent motion is resisted by constant average stress (constant force) applied to a relatively uniform layer or interface at depth. This finding follows new evidence that the mantle lithosphere beneath Tibet is intact, which supports the interpretation that the long-term strain history of Tibet reflects deformation of the mantle lithosphere. Under conditions of constant stress and strength, the deforming continental lithosphere creates a type of viscous resistance that affects plate motion irrespective of how topography evolved.     

When and how did plate tectonics begin？ Theoretical and empirical considerations

Plate tectonics is the horizontal motion of Earth’s thermal boundary layer (lithosphere) over the convecting mantle (asthenosphere) and is mostly driven by lithosphere sinking in subduction zones. Plate tectonics is an outstanding example of a self organizing, far from equilibrium complex system (SOFFECS), driven by the negative buoyancy of the thermal boundary layer and controlled by dissipation in the bending lithosphere and viscous mantle. Plate tectonics is an unusual way for a silicate planet to lose heat, as it exists on only one of the large five silicate bodies in the inner solar system. It is not known when this mode of tectonic activity and heat loss began on Earth. All silicate planets probably experienced a short-lived magma ocean stage. After this solidified, stagnant lid behavior is the common mode of planetary heat loss, with interior heat being lost by delamination and “hot spot” volcanism and shallow intrusions. Decompression melting in the hotter early Earth generated a different lithosphere than today, with thicker oceanic crust and thinner mantle lithosphere; such lithosphere would take much longer than at present to become negatively buoyant, suggesting that plate tectonics on the early Earth occurred sporadically if at all. Plate tectonics became sustainable (the modern style) when Earth cooled sufficiently that decompression melting beneath spreading ridges made thin oceanic crust, allowing oceanic lithosphere to become negatively buoyant after a few tens of millions of years. Ultimately the question of when plate tectonics began must be answered by informa- tion retrieved from the geologic record. Criteria for the operation of plate tectonics includes ophiolites, blueschist and ultra-high pressure metamorphic belts, eclogites, passive margins, transform faults, paleomagnetic demonstration of different motions of different cratons, and the presence of diagnostic geochemical and isotopic indicators in igneous rocks. This record must be interpreted individually; I interpret the record to indicate a progression of tectonic styles from active Archean tectonics and magmatism to something similar to plate tectonics at ~1.9 Ga to sustained, modern style plate tectonics with deep subduction——and powerful slab pull——beginning in Neoproterozoic time.     

Subcontinental-scale crustal velocity changes along the Pacific-North America plate boundary

Transient tectonic deformation has long been noted within approximately 100 km of plate boundary fault zones and within active volcanic regions, but it is unknown whether transient motions also occur at larger scales within plates. Relatively localized transients are known to occur as both seismic and episodic aseismic events, and are generally ascribed to motions of magma bodies, aseismic creep on faults, or elastic or viscoelastic effects associated with earthquakes. However, triggering phenomena and systematic patterns of seismic strain release at subcontinental (approximately 1,000 km) scale along diffuse plate boundaries have long suggested that energy transfer occurs at larger scale. Such transfer appears to occur by the interaction of stresses induced by surface wave propagation and magma or groundwater in the crust, or from large-scale stress diffusion within the oceanic mantle in the decades following clusters of great earthquakes. Here we report geodetic evidence for a coherent, subcontinental-scale change in tectonic velocity along a diffuse approximately 1,000-km-wide deformation zone. Our observations are derived from continuous GPS (Global Positioning System) data collected over the past decade across the Basin and Range province, which absorbs approximately 25 per cent of Pacific-North America relative plate motion. The observed changes in site velocity define a sharp boundary near the centre of the province oriented roughly parallel to the north-northwest relative plate motion vector. We show that sites to the west of this boundary slowed relative to sites east of it by approximately 1 mm yr(-1) starting in late 1999.     

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.     

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.     

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.     

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.     

Seismic reflection imaging of two megathrust shear zones in the northern Cascadia subduction zone

At convergent continental margins, the relative motion between the subducting oceanic plate and the overriding continent is usually accommodated by movement along a single, thin interface known as a megathrust. Great thrust earthquakes occur on the shallow part of this interface where the two plates are locked together. Earthquakes of lower magnitude occur within the underlying oceanic plate, and have been linked to geochemical dehydration reactions caused by the plate's descent. Here I present deep seismic reflection data from the northern Cascadia subduction zone that show that the inter-plate boundary is up to 16 km thick and comprises two megathrust shear zones that bound a >5-km-thick, approximately 110-km-wide region of imbricated crustal rocks. Earthquakes within the subducting plate occur predominantly in two geographic bands where the dip of the plate is inferred to increase as it is forced around the edges of the imbricated inter-plate boundary zone. This implies that seismicity in the subducting slab is controlled primarily by deformation in the upper part of the plate. Slip on the shallower megathrust shear zone, which may occur by aseismic slow slip, will transport crustal rocks into the upper mantle above the subducting oceanic plate and may, in part, provide an explanation for the unusually low seismic wave speeds that are observed there.     

Rejuvenation of the lithosphere by the Hawaiian plume

The volcanism responsible for creating the chain of the Hawaiian islands and seamounts is believed to mark the passage of the oceanic lithosphere over a mantle plume. In this picture hot material rises from great depth within a fixed narrow conduit to the surface, penetrating the moving lithosphere. Although a number of models describe possible plume-lithosphere interactions, seismic imaging techniques have not had sufficient resolution to distinguish between them. Here we apply the S-wave 'receiver function' technique to data of three permanent seismic broadband stations on the Hawaiian islands, to map the thickness of the underlying lithosphere. We find that under Big Island the lithosphere is 100-110 km thick, as expected for an oceanic plate 90-100 million years old that is not modified by a plume. But the lithosphere thins gradually along the island chain to about 50-60 km below Kauai. The width of the thinning is about 300 km. In this zone, well within the larger-scale topographic swell, we infer that the rejuvenation model (where the plume thins the lithosphere) is operative; however, the larger-scale topographic swell is probably supported dynamically.     

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.     

Earthquake slip on oceanic transform faults

Oceanic transform faults are one of the main types of plate boundary, but the manner in which they slip remains poorly understood. Early studies suggested that relatively slow earthquake rupture might be common; moreover, it has been reported that very slow slip precedes some oceanic transform earthquakes, including the 1994 Romanche earthquake. The presence of such detectable precursors would have obvious implications for earthquake prediction. Here we model broadband seismograms of body waves to obtain well-resolved depths and rupture mechanisms for 14 earthquakes on the Romanche and Chain transform faults in the equatorial Atlantic Ocean. We found that earthquakes on the longer Romanche transform are systematically deeper than those on the neighbouring Chain transform. These depths indicate that the maximum depth of brittle failure is at a temperature of approximately 600 degrees C in oceanic lithosphere. We find that the body waves from the Romanche 1994 earthquake can be well modelled with relatively deep slip on a single fault, and we use the mechanism and depth of this earthquake to recalculate its source spectrum. The previously reported slow precursor can be explained as an artefact of uncertainties in the assumed model parameters.     

华南块体各向异性分区及上地幔动力学含义

利用SKS波分裂方法,对布设在华南地区的宽频带流动地震观测台阵数据进行分析。研究结果表明,华南块体(SCB)的SKS波分裂自西向东存在明显的变化。在华南块体西部克拉通岩石圈保留完整的四川盆地区域,SKS波分裂不明显,反映该区域岩石圈厚且缺乏大规模的一致性变形;在四川盆地东缘褶皱带区域,SKS波分裂快波方向主要为NNE向,分裂延迟时间可以达到1s左右,记录了该区域岩石圈地幔经历的一致的显著变形;华南块体东部的SKS波分裂快波方向主要为ENE向,分裂延迟时间为1s左右,该区域各向异性主要来自软流圈流动的贡献,与大地幔楔的流动方向契合;四川盆地东缘以东200km范围内,各向异性结果大部分显示为Null值,显示该区域在地质历史上可能发生过复杂的地球动力学过程。     

Mid-mantle deformation inferred from seismic anisotropy

With time, convective processes in the Earth's mantle will tend to align crystals, grains and inclusions. This mantle fabric is detectable seismologically, as it produces an anisotropy in material properties--in particular, a directional dependence in seismic-wave velocity. This alignment is enhanced at the boundaries of the mantle where there are rapid changes in the direction and magnitude of mantle flow, and therefore most observations of anisotropy are confined to the uppermost mantle or lithosphere and the lowermost-mantle analogue of the lithosphere, the D" region. Here we present evidence from shear-wave splitting measurements for mid-mantle anisotropy in the vicinity of the 660-km discontinuity, the boundary between the upper and lower mantle. Deep-focus earthquakes in the Tonga-Kermadec and New Hebrides subduction zones recorded at Australian seismograph stations record some of the largest values of shear-wave splitting hitherto reported. The results suggest that, at least locally, there may exist a mid-mantle boundary layer, which could indicate the impediment of flow between the upper and lower mantle in this region.     

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.     

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.     

Prediction of Emperor-Hawaii seamount locations from a revised model of global plate motion and mantle flow

The bend in the Hawaiian-Emperor seamount chain is a prominent feature usually attributed to a change in Pacific plate motion approximately 47 Myr ago. However, global plate motion reconstructions fail to predict the bend. Here we show how the geometry of the Hawaiian-Emperor chain and other hotspot tracks can be explained when we combine global plate motions with intraplate deformation and movement of hotspot plumes through distortion by global mantle flow. Global mantle flow models predict a southward motion of the Hawaiian hotspot. This, in combination with a plate motion reconstruction connecting Pacific and African plates through Antarctica, predicts the Hawaiian track correctly since the date of the bend, but predicts the chain to be too far west before it. But if a reconstruction through Australia and Lord Howe rise is used instead, the track is predicted correctly back to 65 Myr ago, including the bend. The difference between the two predictions indicates the effect of intraplate deformation not yet recognized or else not recorded on the ocean floor. The remaining misfit before 65 Myr ago can be attributed to additional intraplate deformation of similar magnitude.     

Magma-maintained rift segmentation at continental rupture in the 2005 Afar dyking episode

Seafloor spreading centres show a regular along-axis segmentation thought to be produced by a segmented magma supply in the passively upwelling mantle. On the other hand, continental rifts are segmented by large offset normal faults, and many lack magmatism. It is unclear how, when and where the ubiquitous segmented melt zones are emplaced during the continental rupture process. Between 14 September and 4 October 2005, 163 earthquakes (magnitudes greater than 3.9) and a volcanic eruption occurred within the approximately 60-km-long Dabbahu magmatic segment of the Afar rift, a nascent seafloor spreading centre in stretched continental lithosphere. Here we present a three-dimensional deformation field for the Dabbahu rifting episode derived from satellite radar data, which shows that the entire segment ruptured, making it the largest to have occurred on land in the era of satellite geodesy. Simple elastic modelling shows that the magmatic segment opened by up to 8 m, yet seismic rupture can account for only 8 per cent of the observed deformation. Magma was injected along a dyke between depths of 2 and 9 km, corresponding to a total intrusion volume of approximately 2.5 km3. Much of the magma appears to have originated from shallow chambers beneath Dabbahu and Gabho volcanoes at the northern end of the segment, where an explosive fissural eruption occurred on 26 September 2005. Although comparable in magnitude to the ten year (1975-84) Krafla events in Iceland, seismic data suggest that most of the Dabbahu dyke intrusion occurred in less than a week. Thus, magma intrusion via dyking, rather than segmented normal faulting, maintains and probably initiated the along-axis segmentation along this sector of the Nubia-Arabia plate boundary.     

Seismic evidence for catastrophic slab loss beneath Kamchatka

In the northwest Pacific Ocean, a sharp corner in the boundary between the Pacific plate and the North American plate joins a subduction zone running along the southern half of the Kamchatka peninsula with a region of transcurrent motion along the western Aleutian arc. Here we present images of the seismic structure beneath the Aleutian-Kamchatka junction and the surrounding region, indicating that: the subducting Pacific lithosphere terminates at the Aleutian-Kamchatka junction; no relict slab underlies the extinct northern Kamchatka volcanic arc; and the upper mantle beneath northern Kamchatka has unusually slow shear wavespeeds. From the tectonic and volcanic evolution of Kamchatka over the past 10 Myr (refs 3-5) we infer that at least two episodes of catastrophic slab loss have occurred. About 5 to 10 Myr ago, catastrophic slab loss shut down island-arc volcanic activity north of the Aleutian-Kamchatka junction. A later episode of slab loss, since about 2 Myr ago, seems to be related to the activity of the world's most productive island-arc volcano, Klyuchevskoy. Removal of lithospheric mantle is commonly discussed in the context of a continental collision, but our findings imply that episodes of slab detachment and loss are also important agents in the evolution of oceanic convergent margins.