Understanding subduction infancy to mature subduction in Southwest Japan via the self-consistent formation of a weak slab interface

The weak slab interface controls long-term subduction dynamics. A weak hydrous layer at the slab interface promotes mechanical decoupling between the forearc mantle and the subducting slab and converts a hot forearc mantle to a cold mantle. Often referred to as a cold nose, the cold forearc mantle, plays a key role in the transition from subduction infancy to mature subduction. This study was the first to numerically demonstrate the self-consistent formation of a weak hydrous layer with permeability anisotropy based on the Southwest Japan subduction zone case, where transition-related geological features were present. Our models showed that mechanical decoupling by spontaneous downdip growth of the weak hydrous layer created a cold nose by converting a hot forearc mantle to a cold mantle. The emergence of a cold nose explained the migration of the forearc-to-arc volcanic front, expressed as the formation of mid-Miocene forearc high-magnesium andesite and Quaternary arc adakite. Furthermore, the weak hydrous layer providing a pathway for free-water transport toward the mantle wedge tip elucidates slab/mantle-derived geochemical components in deep groundwater as well as large S-wave delay times and non-volcanic seismic tremors in the forearc.

Rheology of Naturally Deformed Antigorite Serpentinite: Strain and Strain-Rate Dependence at Mantle-Wedge Conditions

Antigorite serpentinite is expected to occur in parts of subduction plate boundaries, and may suppress earthquake slip, but the dominant deformation mechanisms and resultant rheology of antigorite are unclear. An exhumed plate boundary shear zone exposed near Nagasaki, Japan, contains antigorite deformed at 474°C ± 30°C. Observations indicate that a foliation defined by (001) crystal facets developed during plate-boundary shear. Microstructures indicating grain-scale dissolution at high-stress interfaces and precipitation in low-stress regions suggest that dissolution-precipitation creep contributed to foliation development. Analysis of crystal orientations indicate a small contribution from dislocation activity. We suggest a frictional-viscous rheology for antigorite, where dissolution-precipitation produces a foliation defined by (001) crystal facets and acts to resolve strain incompatibilities, allowing for efficient face-to-face sliding between facets. This rheology can not only explain aseismic behavior at ambient plate boundary conditions, but also some of the contrasting behaviors shown by previous field and laboratory studies.

In Situ Study on Dehydration and Phase Transformation of Antigorite

Antigorite is the main carrier of water in Earth’s subduction zones. The dehydration processes of antigorite were investigated by carrying out in situ phase transition experiments using a dynamic diamond anvil cell, with a time-resolved Raman scattering system, at 0.3–10 GPa and 396–1100 K. Three typical phase transformation reactions occurred within the P–T range of this study, corresponding to three reaction products. At low pressures (<0.7 GPa), antigorite transfers to talc and forsterite; as the temperature increases, the talc disappears and a combination of forsterite and clinoenstatite occurs. At moderate pressures (1.8–7.5 GPa), antigorite dehydrates into forsterite and clinoenstatite as temperatures increase; with the continuous increase in pressure, the dehydration products become clinoenstatite and phase A. At high pressures (>8.6 GPa), the products of the dehydration phase transition of antigorite are consistently clinoenstatite and phase A. Compared with the previous studies carried out by large-volume presses (such as a multi anvil press and a piston-cylinder press), the reaction to produce phase A occurs at higher P–T conditions, and the stable temperature region for talc as a dehydration product is narrower. Moreover, large quantities of pores with 5–10 μm in diameter formed in dehydration products, supporting the hypothesis that intermediate-depth earthquakes may result from dehydration embrittlement. The precise phase boundary determined by this in situ study provides a better understanding of the dehydration phase transition behavior and geological phenomena exhibited by antigorite under different pressure and temperature conditions.

Influence of Petrogenesis on the Engineering Properties of Ultramafic Aggregates and on Their Suitability in Concrete

This study focuses on the use of petrology as a useful tool in construction applications (i.e., concrete). More specifically, this study investigates how the petrogenetic characteristics of ultramafic rocks derived from ophiolite complexes (Veria–Naousa, Gerania) can act as a key tool for the prediction of the final behaviour of ultramafic aggregates as concrete aggregates. For this reason, their petrographic, chemical and mineralogical characteristics were examined and correlated with their engineering properties for evaluating their suitability as concrete aggregates. This study had come up, for the first time, that the genesis environment of the ultramafic rocks is the determinant factor for their physico-mechanical characteristics. Their suitability is relevant to the impact of their mineralogical and structural characteristics both from the two different ophiolite complexes (Veria–Naousa and Gerania). Except serpentinization, the basic alteration process-index of ultramafic rocks, there are also other chemical indices which can be used for ultramafic rocks that may determine their properties. In this context the term ‘fertility rate’ (FR) was introduced which may characterize ultramafic rocks as fertile or not. Furthermore, the Ultramafic Rock Health Index (U.R.H.I.) as well as the Normalized U.R.H.I. (U.R.H.I.N) was also introduced and correlated with the engineering properties of the investigated aggregate rocks. The last index aims to assess and quantify the overall health conditions, encompassing the two major modifying factors that include removal of primary mineral phases, as well as the extent of serpentinization. The main conclusion of this paper is that the genesis environment of the ultramafic rocks is the critical factor that determines their mineralogical, petrographic and chemical characteristics which consequently determines the basic engineering properties of rocks that determine their suitability or not as concrete aggregates.

Metamorphism-facilitated faulting in deforming orthopyroxene: Implications for global intermediate-depth seismicity

The exothermic metamorphic reaction in orthopyroxene (Opx), a major component of oceanic lithospheric mantle, is shown to trigger brittle failure in laboratory deformation experiments under conditions where garnet exsolution takes place. The reaction product is an extremely fine-grained material, forming narrow reaction zones that are mechanically weak, thereby facilitating macroscopic faulting. Oceanic subduction zones are characterized by two separate bands of seismicity, known as the double seismic zone. The upper band of seismicity, located in the oceanic crust, is well explained by dehydration-induced mechanical instability. Our newly discovered metamorphism-induced mechanical instability provides an alternative physical mechanism for earthquakes in the lower band of seismicity (located in the oceanic lithospheric mantle), with no requirement of hydration/dehydration processes.

Dynamic slab segmentation due to brittle-ductile damage in the outer rise

Subduction is the major plate driving force, and the strength of the subducting plate controls many aspects of the thermochemical evolution of Earth. Each subducting plate experiences intense normal faulting^1-9 during bending that accommodates the transition from horizontal to downwards motion at the outer rise at trenches. Here we investigate the consequences of this bending-induced plate damage using numerical subduction models in which both brittle and ductile deformation, including grain damage, are tracked and coupled self-consistently. Pervasive slab weakening and pronounced segmentation can occur at the outer-rise region owing to the strong feedback between brittle and ductile damage localization. This slab-damage phenomenon explains the subduction dichotomy of strong plates and weak slabs¹⁰, the development of large-offset normal faults^6,7 near trenches, the occurrence of segmented seismic velocity anomalies¹¹ and distinct interfaces imaged within subducted slabs^12,13, and the appearance of deep, localized intraplate areas of reduced effective viscosity¹⁴ observed at trenches. Furthermore, brittle-viscously damaged slabs show a tendency for detachment at elevated mantle temperatures. Given Earth's planetary cooling history¹⁵, this implies that intermittent subduction with frequent slab break-off episodes¹⁶ may have been characteristic for Earth until more recent times than previously suggested¹⁷.

Accurate post-mortem alignment for Focused Ion Beam and Scanning Electron Microscopy (FIB-SEM) tomography

Drifts in the three directions (X, Y, Z) during the FIB-SEM slice-and-view tomography is an important issue in 3D-FIB experiments which may induce significant inaccuracies in the subsequent volume reconstruction and further quantification of morphological volume parameters of the sample microstructure. Cross-correlation is frequently applied directly to the cross-section image series for aligning FIB sliced images. This solution is hazardous and can be flawed as it has been easily demonstrated by a dedicated test experiment. As a result, a novel aligning procedure based on the quantification of the topography of the sample surface has been developed. This new approach will be compared to the common cross-correlation methods, as well as another approach consisting in using an artificial reference marker fabricated during the FIB procedure. All these methods will then be discussed in terms of accuracy and liability.

What's down there? The structures, materials and environment of deep-seated slow slip and tremor

Deep-seated slow slip and tremor (SST), including slow slip events, episodic tremor and slip, and low-frequency earthquakes, occur downdip of the seismogenic zone of numerous subduction megathrusts and plate boundary strike-slip faults. These events represent a fascinating and perplexing mode of fault failure that has greatly broadened our view of earthquake dynamics. In this contribution, we review constraints on SST deformation processes from both geophysical observations of active subduction zones and geological observations of exhumed field analogues. We first provide an overview of what has been learned about the environment, kinematics and dynamics of SST from geodetic and seismologic data. We then describe the materials, deformation mechanisms, and metamorphic and fluid pressure conditions that characterize exhumed rocks from SST source depths. Both the geophysical and geological records strongly suggest the importance of a fluid-rich and high fluid pressure habitat for the SST source region. Additionally, transient deformation features preserved in the rock record, involving combined frictional-viscous shear in regions of mixed lithology and near-lithostatic fluid pressures, may scale with the tremor component of SST. While several open questions remain, it is clear that improved constraints on the materials, environment, structure, and conditions of the plate interface from geophysical imaging and geologic observations will enhance model representations of the boundary conditions and geometry of the SST deformation process. This article is part of a discussion meeting issue 'Understanding earthquakes using the geological record'.

Extrusion of subducted crust explains the emplacement of far-travelled ophiolites

Continental subduction below oceanic plates and associated emplacement of ophiolite sheets remain enigmatic chapters in global plate tectonics. Numerous ophiolite belts on Earth exhibit a far-travelled ophiolite sheet that is separated from its oceanic root by tectonic windows exposing continental crust, which experienced subduction-related high pressure-low temperature metamorphism during obduction. However, the link between continental subduction-exhumation dynamics and far-travelled ophiolite emplacement remains poorly understood. Here we combine data collected from ophiolite belts worldwide with thermo-mechanical simulations of continental subduction dynamics to show the causal link between the extrusion of subducted continental crust and the emplacement of far-travelled ophiolites. Our results reveal that buoyancy-driven extrusion of subducted crust triggers necking and breaking of the overriding oceanic upper plate. The broken-off piece of oceanic lithosphere is then transported on top of the continent along a flat thrust segment and becomes a far-travelled ophiolite sheet separated from its root by the extruded continental crust. Our results indicate that the extrusion of the subducted continental crust and the emplacement of far-travelled ophiolite sheets are inseparable processes.

Hydrous oceanic crust hosts megathrust creep at low shear stresses

The rheology of the metamorphosed oceanic crust may be a critical control on megathrust strength and deformation style. However, little is known about the strength and deformation style of metamorphosed basalt. Exhumed megathrust shear zones exposed on Kyushu, SW Japan, contain hydrous metabasalts deformed at temperatures between ~300° and ~500°C, spanning the inferred temperature-controlled seismic-aseismic transition. Field and microstructural observations of these shear zones, combined with quartz grain-size piezometry, indicate that metabasalts creep at shear stresses <100 MPa at ~370°C and at shear stresses <30 MPa at ~500°C. These values are much lower than those suggested by viscous flow laws for basalt. The implication is that relatively weak, hydrous, metamorphosed oceanic crust can creep at low viscosities over a wide shear zone and have a critical influence on plate interface strength and deformation style around the seismic-aseismic transition.

Evidence for a serpentinized plate interface favouring continental subduction

The dynamics of continental subduction is largely controlled by the rheological properties of rocks involved along the subduction channel. Serpentinites have low viscosity at geological strain rates. However, compelling geophysical evidence of a serpentinite channel during continental subduction is still lacking. Here we show that anomalously low shear-wave seismic velocities are found beneath the Western Alps, along the plate interface between the European slab and the overlying Adriatic mantle. We propose that these seismic velocities indicate the stacked remnants of a weak fossilised serpentinite channel, which includes both slivers of abyssal serpentinite formed at the ocean floor and mantle-wedge serpentinite formed by fluid release from the subducting slab. Our results suggest that this serpentinized plate interface may have favoured the subduction of continental crust into the upper mantle and the formation/exhumation of ultra-high pressure metamorphic rocks, providing new constraints to develop the conceptual and quantitative understanding of continental-subduction dynamics.

The dynamics of continental subduction is largely controlled by the rheological properties of rocks involved along the subduction channel. Here, the authors reveal a prominent, yet previously undetected, low-velocity body beneath the Western Alps, along the plate interface between the European slab and the overlying Adriatic mantle, which they interpret as a serpentinite layer.

Oceanic crust recycling controlled by weakening at slab edges

Retreating subduction zones such as the Lesser Antilles, Gibraltar and Scotia have been migrating towards the Atlantic Ocean by cutting their way through the oceanic crust. This spontaneously retreating subduction is enabled by the development of faults at the edges of the slab, but the physical mechanisms controlling fault propagation and direction remain unknown. Here, using 3D numerical subduction models we show that oceanic lithosphere recycling is mainly controlled by the intensity of strain-induced weakening of fractures forming at the edges of the slab. Intense strain-induced weakening causes predominantly brittle fault propagation and slab narrowing until detachment. Without weakening, preponderantly ductile slab edge propagation occurs, which causes slab widening. This rheological control is not affected by the proximity of non-weakened passive continental margins. Natural examples suggest that slab edges follow convergent paths that could be controlled by fractures weakening due to deep water penetration into the oceanic lithosphere.

In Situ Quantitative Tensile Testing of Antigorite in a Transmission Electron Microscope

The determination of the mechanical properties of serpentinites is essential toward the understanding of the mechanics of faulting and subduction. Here we present the first in situ tensile tests on antigorite in a transmission electron microscope. A push-to-pull deformation device is used to perform quantitative tensile tests, during which force and displacement are measured, while the evolving microstructure is imaged with the microscope. The experiments have been performed at room temperature on 2 × 1 × 0.2 μm3 beams prepared by focused ion beam. The specimens are not single crystals despite their small sizes. Orientation mapping indicated that several grains were well oriented for plastic slip. However, no dislocation activity has been observed even though the engineering tensile stress went up to 700 MPa. We show also that antigorite does not exhibit a purely elastic-brittle behavior since, despite the presence of defects, the specimens accumulate permanent deformation and did not fail within the elastic regime. Instead, we observe that strain localizes at grain boundaries. All observations concur to show that under these experimental conditions, grain boundary sliding is the dominant deformation mechanism. This study sheds a new light on the mechanical properties of antigorite and calls for further studies on the structure and properties of grain boundaries in antigorite and more generally in phyllosilicates.

Insight into the microphysics of antigorite deformation from spherical nanoindentation

The mechanical behaviour of antigorite strongly influences the strength and deformation of the subduction interface. Although there is microstructural evidence elucidating the nature of brittle deformation at low pressures, there is often conflicting evidence regarding the potential for plastic deformation in the ductile regime at higher pressures. Here, we present a series of spherical nanoindentation experiments on aggregates of natural antigorite. These experiments effectively investigate the single-crystal mechanical behaviour because the volume of deformed material is significantly smaller than the grain size. Individual indents reveal elastic loading followed by yield and strain hardening. The magnitude of the yield stress is a function of crystal orientation, with lower values associated with indents parallel to the basal plane. Unloading paths reveal more strain recovery than expected for purely elastic unloading. The magnitude of inelastic strain recovery is highest for indents parallel to the basal plane. We also imposed indents with cyclical loading paths, and observed strain energy dissipation during unloading-loading cycles conducted up to a fixed maximum indentation load and depth. The magnitude of this dissipated strain energy was highest for indents parallel to the basal plane. Subsequent scanning electron microscopy revealed surface impressions accommodated by shear cracks and a general lack of dislocation-induced lattice misorientation. Based on these observations, we suggest that antigorite deformation at high pressures is dominated by sliding on shear cracks. We develop a microphysical model that is able to quantitatively explain Young's modulus and dissipated strain energy data during cyclic loading experiments, based on either frictional or cohesive sliding of an array of cracks contained in the basal plane. This article is part of a discussion meeting issue 'Serpentinite in the earth system'.

Impact of interseismic deformation on phase transformations and rock properties in subduction zones

Phase transformations greatly affect physical properties of rocks and impose a first-order control on geodynamic processes. Under high deformation rates, rheological heterogeneities cause large spatial variations of stress in materials. Until now, the impact of higher deformation rates, rock heterogeneity and stress build up on phase transformations and material properties is not well understood. Here we show, that phase transitions are controlled by the stress build-up during fast deformation. In a deformation experiment (600 °C, 1.47 GPa), rock heterogeneity was simulated by a strong elliptical alumina inclusion in a weak calcite matrix. Under deformation rates comparable to slow earthquakes, calcite transformed locally to aragonite matching the distribution of maximum principal stresses and pressure (mean stress) from mechanical models. This first systematic investigation documents that phase transformations occur in a dynamic system during deformation. The ability of rocks to react during fast deformation rates may have serious consequences on rock rheology and thus provide unique information on the processes leading to giant ruptures in subduction zones.

Mechanism of subsidence of the Northeast Japan forearc during the late period of a gigantic earthquake cycle

The forearc in Northeast Japan subsided (3–4 mm/year) in the interseismic ~100 years before the 2011 Tohoku earthquake (M_W9.1) just like it did during this event. This study attempts to understand the mechanism of the vertical displacement of the forearc during gigantic earthquake cycles via numerical modeling. The results suggest that the interseismic subsidence rate in the forearc increases with the duration of the locking of the asperity of the gigantic earthquake over several hundred years, due to the increasing slip deficit rate on the deeper parts of the plate interface. The increasing slip deficit rate is caused by both the decreasing the shear stress in the shear zone owing to the continuous locking of the asperity and the increasing the mobility of the continental lithosphere owing to the viscoelastic relaxation in the mantle wedge. The deep slip deficit rate extending to ~100 km depth of the plate interface is necessary to explain the observed interseismic forearc subsidence rate. The results also suggest hundreds of years of continuous locking of the asperities of a gigantic earthquake in the western Kuril subduction zone, where fast forearc subsidence has been observed as well.

The equation of motion for supershear frictional rupture fronts

The rupture fronts that mediate the onset of frictional sliding may propagate at speeds below the Rayleigh wave speed or may surpass the shear wave speed and approach the longitudinal wave speed. While the conditions for the transition from sub-Rayleigh to supershear propagation have been studied extensively, little is known about what dictates supershear rupture speeds and how the interplay between the stresses that drive propagation and interface properties that resist motion affects them. By combining laboratory experiments and numerical simulations that reflect natural earthquakes, we find that supershear rupture propagation speeds can be predicted and described by a fracture mechanics-based equation of motion. This equation of motion quantitatively predicts rupture speeds, with the velocity selection dictated by the interface properties and stress. Our results reveal a critical rupture length, analogous to Griffith's length for sub-Rayleigh cracks, below which supershear propagation is impossible. Above this critical length, supershear ruptures can exist, once excited, even for extremely low preexisting stress levels. These results significantly improve our fundamental understanding of what governs the speed of supershear earthquakes, with direct and important implications for interpreting their unique supershear seismic radiation patterns.

Slow-slip events in semi-brittle serpentinite fault zones

Slow-slip events are earthquake-like events only with much lower slip rates. While peak coseismic velocities can reach tens of meters per second, slow-slip is on the order of 10^−7±2 m/s and may last for days to weeks. Under the rate-and-state model of fault friction, slow-slip is produced only when the asperity size is commensurate with the critical nucleation size, a function of frictional properties. However, it is unlikely that all subduction zones embody the same frictional properties. In addition to friction, plastic flow of antigorite-rich serpentinite may significantly influence the dynamics of fault slip near the mantle wedge corner. Here, we show that the range of frictional parameters that generate slow slip is widened in the presence of a serpentinized layer along the subduction plate interface. We observe increased stability and damping of fast ruptures in a semi-brittle fault zone governed by both brittle and viscoelastic constitutive response. The rate of viscous serpentinite flow, governed by dislocation creep, is enhanced by high ambient temperatures. When effective viscosity is taken to be dynamic, long-term slow slip events spontaneously emerge. Integration of rheology, thermal effects, and other microphysical processes with rate-and-state friction may yield further insight into the phenomenology of slow slip.

High-pressure studies with x-rays using diamond anvil cells

Pressure profoundly alters all states of matter. The symbiotic development of ultrahigh-pressure diamond anvil cells, to compress samples to sustainable multi-megabar pressures; and synchrotron x-ray techniques, to probe materials' properties in situ, has enabled the exploration of rich high-pressure (HP) science. In this article, we first introduce the essential concept of diamond anvil cell technology, together with recent developments and its integration with other extreme environments. We then provide an overview of the latest developments in HP synchrotron techniques, their applications, and current problems, followed by a discussion of HP scientific studies using x-rays in the key multidisciplinary fields. These HP studies include: HP x-ray emission spectroscopy, which provides information on the filled electronic states of HP samples; HP x-ray Raman spectroscopy, which probes the HP chemical bonding changes of light elements; HP electronic inelastic x-ray scattering spectroscopy, which accesses high energy electronic phenomena, including electronic band structure, Fermi surface, excitons, plasmons, and their dispersions; HP resonant inelastic x-ray scattering spectroscopy, which probes shallow core excitations, multiplet structures, and spin-resolved electronic structure; HP nuclear resonant x-ray spectroscopy, which provides phonon densities of state and time-resolved Mössbauer information; HP x-ray imaging, which provides information on hierarchical structures, dynamic processes, and internal strains; HP x-ray diffraction, which determines the fundamental structures and densities of single-crystal, polycrystalline, nanocrystalline, and non-crystalline materials; and HP radial x-ray diffraction, which yields deviatoric, elastic and rheological information. Integrating these tools with hydrostatic or uniaxial pressure media, laser and resistive heating, and cryogenic cooling, has enabled investigations of the structural, vibrational, electronic, and magnetic properties of materials over a wide range of pressure-temperature conditions.

Seismic evidence for flow in the hydrated mantle wedge of the Ryukyu subduction zone

It is widely accepted that water-rich serpentinite domains are commonly present in the mantle above shallow subducting slabs and play key roles in controlling the geochemical cycling and physical properties of subduction zones. Thermal and petrological models show the dominant serpentine mineral is antigorite. However, there is no good consensus on the amount, distribution and alignment of this mineral. Seismic velocities are commonly used to identify antigorite-rich domains, but antigorite is highly-anisotropic and depending on the seismic ray path, its properties can be very difficult to distinguish from non-hydrated olivine-rich mantle. Here, we utilize this anisotropy and show how an analysis of seismic anisotropy that incorporates measured ray path geometries in the Ryukyu arc can constrain the distribution, orientation and amount of antigorite. We find more than 54% of the wedge must consist of antigorite and the alignment must change from vertically aligned to parallel to the slab. This orientation change suggests convective flow in the hydrated forearc mantle. Shear wave splitting analysis in other subduction zones indicates large-scale serpentinization and forearc mantle convection are likely to be more widespread than generally recognized. The view that the forearc mantle of cold subduction zones is dry needs to be reassessed.

Dehydration of lawsonite could directly trigger earthquakes in subducting oceanic crust

Intermediate-depth earthquakes in cold subduction zones are observed within the subducting oceanic crust, as well as the mantle. In contrast, intermediate-depth earthquakes in hot subduction zones predominantly occur just below the Mohorovičić discontinuity. These observations have stimulated interest in relationships between blueschist-facies metamorphism and seismicity, particularly through dehydration reactions involving the mineral lawsonite. Here we conducted deformation experiments on lawsonite, while monitoring acoustic emissions, in a Griggs-type deformation apparatus. The temperature was increased above the thermal stability of lawsonite, while the sample was deforming, to test whether the lawsonite dehydration reaction induces unstable fault slip. In contrast to similar tests on antigorite, unstable fault slip (that is, stick-slip) occurred during dehydration reactions in the lawsonite and acoustic emission signals were continuously observed. Microstructural observations indicate that strain is highly localized along the fault (R1 and B shears), and that the fault surface develops slickensides (very smooth fault surfaces polished by frictional sliding). The unloading slope during the unstable slip follows the stiffness of the apparatus at all experimental conditions, regardless of the strain rate and temperature ramping rate. A thermomechanical scaling factor for the experiments is within the range estimated for natural subduction zones, indicating the potential for unstable frictional sliding within natural lawsonite layers.

Time-lapse nanoscopy of friction in the non-Amontons and non-Coulomb regime

Originally discovered by Leonard da Vinci in the 15th century, the force of friction is directly proportional to the applied load (known as Amontons' first law of friction). Furthermore, kinetic friction is independent of the sliding speed (known as Coulomb's law of friction). These empirical laws break down at high normal pressure (due to plastic deformation) and low sliding speed (in the transition regime between static friction and kinetic friction). An important example of this phenomenon is friction between the asperities of tectonic plates on the Earth. Despite its significance, little is known about the detailed mechanism of friction in this regime due to the lack of experimental methods. Here we demonstrate in situ time-lapse nanoscopy of friction between asperities sliding at ultralow speed (∼0.01 nm/s) under high normal pressure (∼GPa). This is made possible by compressing and rubbing a pair of nanometer-scale crystalline silicon anvils with electrostatic microactuators and monitoring its dynamical evolution with a transmission electron microscope. Our analysis of the time-lapse movie indicates that superplastic behavior is induced by decrystallization, plastic deformation, and atomic diffusion at the asperity-asperity interface. The results hold great promise for a better understanding of quasi-static friction under high pressure for geoscience, materials science, and nanotechnology.

Subduction of fracture zones controls mantle melting and geochemical signature above slabs

For some volcanic arcs, the geochemistry of volcanic rocks erupting above subducted oceanic fracture zones is consistent with higher than normal fluid inputs to arc magma sources. Here we use enrichment of boron (B/Zr) in volcanic arc lavas as a proxy to evaluate relative along-strike inputs of slab-derived fluids in the Aleutian, Andean, Cascades and Trans-Mexican arcs. Significant B/Zr spikes coincide with subduction of prominent fracture zones in the relatively cool Aleutian and Andean subduction zones where fracture zone subduction locally enhances fluid introduction beneath volcanic arcs. Geodynamic models of subduction have not previously considered how fracture zones may influence the melt and fluid distribution above slabs. Using high-resolution three-dimensional coupled petrological-thermomechanical numerical simulations of subduction, we show that enhanced production of slab-derived fluids and mantle wedge melts concentrate in areas where fracture zones are subducted, resulting in significant along-arc variability in magma source compositions and processes.

Slow earthquakes associated with fault healing on a serpentinized plate interface

Slow earthquakes that occur at subduction zones are distinct from regular earthquakes in terms of their slip behavior. We consider this difference to relate to localized hydration reactions at the plate interface that influence the frictional properties. The results of laboratory friction experiments indicate that simulated serpentine faults are characterized by a low healing rate and large slip-weakening distance compared with unaltered dry fault patches. These results are consistent with the slip mechanism of slow earthquakes, indicating that a locally serpentinized plate interface could trigger slow earthquakes, assisted by pore pressure build-up, whereas unaltered dry patches that remain strongly coupled are potential sites of regular earthquakes.

High-pressure experimental geosciences: state of the art and prospects

This paper aims at reviewing the current advancements of high pressure experimental geosciences. The angle chosen is that of in situ measurements at the high pressure (P) and high temperature (T) conditions relevant of the deep Earth and planets, measurements that are often carried out at large facilities (X-ray synchrotrons and neutron sources). Rather than giving an exhaustive catalogue, four main active areas of research are chosen: the latest advancements on deep Earth mineralogy, how to probe the properties of melts, how to probe Earth dynamics, and chemical reactivity induced by increased P-T conditions. For each area, techniques are briefly presented and selected examples illustrate their potentials, and what that tell us about the structure and dynamics of the planet.

Electron back-scattering diffraction (EBSD) measurements of antigorite lattice-preferred orientations (LPO)

Lattice preferred orientations of serpentines induce a strong anisotropy of various properties in serpentine bearing-rocks. Lattice preferred orientations had so far been obtained only by X-ray diffraction techniques. We have applied electron back-scattering diffraction to the measurement of the lattice preferred orientations of antigorite in a naturally deformed high-pressure serpentinite. This technique is very sensitive to sample preparation that can lead to surface amorphization in the case of serpentine. A polishing procedure is described that avoids amorphization and allows accurate electron back-scattering diffraction measurements with optimized experimental conditions in a variable pressure scanning electron microscope. Results indicate that deformation leads to lattice preferred orientations characterized by extremely strong c-axis clustering perpendicular to the foliation, as expected for a layered silicate. In the foliation plane, a significant clustering of the a-axis is observed and tentatively attributed to intracrystalline deformation mechanisms. These data suggest that antigorite deforms mostly by gliding along the basal plane of the layered phyllosilicate structure, but that gliding may occur along directions favouring a-axis alignment. Electron back-scattering diffraction appears to be a reliable method for determining phyllosilicate lattice preferred orientations in deformed rocks, with potential applications for determining anisotropy of properties like seismic velocities or thermal and electrical conductivities.