Fault-network geometry influences earthquake frictional behaviour.

Understanding the factors governing the stability of fault slip is a crucial problem in fault mechanics1-3. The importance of fault geometry and roughness on fault-slip behaviour has been highlighted in recent lab experiments4-7 and numerical models8-11, and emerging evidence suggests that large-scale complexities in fault networks have a vital role in the fault-rupture process12-18. Here we present a new perspective on fault creep by investigating the link between fault-network geometry and surface creep rates in California, USA. Our analysis reveals that fault groups exhibiting creeping behaviour show smaller misalignment in their fault-network geometry. The observation indicates that the surface fault traces of creeping regions tend to be simple, whereas locked regions tend to be more complex. We propose that the presence of complex fault-network geometries results in geometric locking that promotes stick-slip behaviour characterized by earthquakes, whereas simpler geometries facilitate smooth fault creep. Our findings challenge traditional hypotheses on the physical origins of fault creep explained primarily in terms of fault friction19-21 and demonstrate the potential for a new framework in which large-scale earthquake frictional behaviour is determined by a combination of geometric factors and rheological yielding properties.

Recovery and facets for deformation twins in minerals and metals.

Type II and IV twins with irrational twin boundaries are studied by high-resolution transmission electron microscopy in two plagioclase crystals. The twin boundaries in these and in NiTi are found to relax to form rational facets separated by disconnections. The topological model (TM), amending the classical model, is required for a precise theoretical prediction of the orientation of the Type II/IV twin plane. Theoretical predictions also are presented for types I, III, V, and VI twins. The relaxation process that forms a faceted structure entails a separate prediction from the TM. Hence, faceting provides a difficult test for the TM. Analysis of the faceting by the TM is in excellent agreement with the observations.

Type III and IV deformation twins in minerals and metals.

SignificanceHistorically, two types of twins (I and II) have been categorized for twinning in minerals and metals. When analyzed by the topological model, a crystallographic construction used to define the defect structure of interfaces, triclinic and some other low-symmetry crystals do not fall into either category and instead form two new twinning types, namely, III and IV. Aside from accurately describing twin structures, these concepts are important for understanding the deformation of minerals such as plagioclase and for deriving constitutive models for the deformation.

Disclinations and disconnections in minerals and metals.

A different type of defect, the coherency disclination, is added to disclination types. Disconnections that include disclination content are considered. A criterion is suggested to distinguish disconnections with dislocation content from those with disclination content. Electron microscopy reveals unit disconnections in a low albite grain boundary, defects important in grain boundary sliding. Disconnections of varying step heights are displayed and shown to define both deformed and recovered structures.

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 Mohorovicic 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.

High-resolution seismic constraints on flow dynamics in the oceanic asthenosphere.

Convective flow in the mantle and the motions of tectonic plates produce deformation of Earth's interior, and the rock fabric produced by this deformation can be discerned using the anisotropy of the seismic wave speed. This deformation is commonly inferred close to lithospheric boundaries beneath the ocean in the uppermost mantle, including near seafloor-spreading centres as new plates are formed via corner flow, and within a weak asthenosphere that lubricates large-scale plate-driven flow and accommodates smaller scale convection. Seismic models of oceanic upper mantle differ as to the relative importance of these deformation processes: seafloor spreading fabric is very strong just beneath the crust-mantle boundary (the Mohorovicic discontinuity, or Moho) at relatively local scales, but at the global and ocean-basin scales, oceanic lithosphere typically appears weakly anisotropic when compared to the asthenosphere. Here we use Rayleigh waves, recorded across an ocean-bottom seismograph array in the central Pacific Ocean (the NoMelt Experiment), to provide unique localized constraints on seismic anisotropy within the oceanic lithosphere-asthenosphere system in the middle of a plate. We find that azimuthal anisotropy is strongest within the high-seismic-velocity lid, with the fast direction coincident with seafloor spreading. A minimum in the magnitude of azimuthal anisotropy occurs within the middle of the seismic low-velocity zone, and then increases with depth below the weakest portion of the asthenosphere. At no depth does the fast direction correlate with the apparent plate motion. Our results suggest that the highest strain deformation in the shallow oceanic mantle occurs during corner flow at the ridge axis, and via pressure-driven or buoyancy-driven flow within the asthenosphere. Shear associated with motion of the plate over the underlying asthenosphere, if present, is weak compared to these other processes.

Carbon-dioxide-rich silicate melt in the Earth's upper mantle.

The onset of melting in the Earth's upper mantle influences the thermal evolution of the planet, fluxes of key volatiles to the exosphere, and geochemical and geophysical properties of the mantle. Although carbonatitic melt could be stable 250 km or less beneath mid-oceanic ridges, owing to the small fraction (∼0.03 wt%) its effects on the mantle properties are unclear. Geophysical measurements, however, suggest that melts of greater volume may be present at ∼200 km (refs 3-5) but large melt fractions are thought to be restricted to shallower depths. Here we present experiments on carbonated peridotites over 2-5 GPa that constrain the location and the slope of the onset of silicate melting in the mantle. We find that the pressure-temperature slope of carbonated silicate melting is steeper than the solidus of volatile-free peridotite and that silicate melting of dry peridotite + CO(2) beneath ridges commences at ∼180 km. Accounting for the effect of 50-200 p.p.m. H(2)O on freezing point depression, the onset of silicate melting for a sub-ridge mantle with ∼100 p.p.m. CO(2) becomes as deep as ∼220-300 km. We suggest that, on a global scale, carbonated silicate melt generation at a redox front ∼250-200 km deep, with destabilization of metal and majorite in the upwelling mantle, explains the oceanic low-velocity zone and the electrical conductivity structure of the mantle. In locally oxidized domains, deeper carbonated silicate melt may contribute to the seismic X-discontinuity. Furthermore, our results, along with the electrical conductivity of molten carbonated peridotite and that of the oceanic upper mantle, suggest that mantle at depth is CO(2)-rich but H(2)O-poor. Finally, carbonated silicate melts restrict the stability of carbonatite in the Earth's deep upper mantle, and the inventory of carbon, H(2)O and other highly incompatible elements at ridges becomes controlled by the flux of the former.

A pressure solution flow law for the seismogenic zone: Application to Cascadia.

We develop a linear viscous constitutive relationship for pressure solution constrained by models of deformed metasedimentary rocks and observations of exposed rocks from ancient subduction zones. We include pressure and temperature dependence on the solubility of silica in fluid by parameterizing a practical van't Hoff relationship. This general flow law is well suited for making predictions about interseismic behavior of subduction zones. We apply the flow law to Cascadia, where thermal structure, geometry, relative plate velocity, and Global Positioning System velocity field are well constrained. Results are consistent with the temperature conditions at which resolvable ductile strain is recorded in subducted mudstones (at depths near the updip limit of the seismogenic zone) and with relative plate motion accommodated completely by viscous deformation (at depths near the downdip limit of the seismogenic zone). The flow law also predicts the observed forearc tapering of slip rate deficit with depth.

A periodic shear-heating mechanism for intermediate-depth earthquakes in the mantle.

Intermediate-depth earthquakes, at depths of 50-300 km in subduction zones, occur below the brittle-ductile transition, where high pressures render frictional failure unlikely. Their location approximately coincides with 600 to 800 degrees C isotherms in thermal models, suggesting a thermally activated mechanism for their origin. Some earthquakes may occur by frictional failure owing to high pore pressure that might result from metamorphic dehydration. Because some intermediate-depth earthquakes occur approximately 30 to 50 km below the palaeo-sea floor, however, the hydrous minerals required for the dehydration mechanism may not be present. Here we present an alternative mechanism to explain such earthquakes, involving the onset of highly localized viscous creep in pre-existing, fine-grained shear zones. Our numerical model uses olivine flow laws for a fine-grained, viscous shear zone in a coarse-grained, elastic half space, with initial temperatures from 600-800 degrees C and background strain rates of 10(-12) to 10(-15) s(-1). When shear heating becomes important, strain rate and temperature increase rapidly to over 1 s(-1) and 1,400 degrees C. The stress then drops dramatically, followed by low strain rates and cooling. Continued far-field deformation produces a quasi-periodic series of such instabilities.

Ductile deformation during carbonation of serpentinized peridotite.

Carbonated serpentinites (listvenites) in the Samail Ophiolite, Oman, record mineralization of 1-2 Gt of CO2, but the mechanisms providing permeability for continued reactive fluid flow are unclear. Based on samples of the Oman Drilling Project, here we show that listvenites with a penetrative foliation have abundant microstructures indicating that the carbonation reaction occurred during deformation. Folded magnesite veins mark the onset of carbonation, followed by deformation during carbonate growth. Undeformed magnesite and quartz overgrowths indicate that deformation stopped before the reaction was completed. We propose deformation by dilatant granular flow and dissolution-precipitation assisted the reaction, while deformation in turn was localized in the weak reacting mass. Lithostatic pore pressures promoted this process, creating dilatant porosity for CO2 transport and solid volume increase. This feedback mechanism may be common in serpentinite-bearing fault zones and the mantle wedge overlying subduction zones, allowing massive carbonation of mantle rocks.

Geophysical evidence from the MELT area for compositional controls on oceanic plates.

Magnetotelluric and seismic data, collected during the MELT experiment at the southern East Pacific Rise, constrain the distribution of melt beneath this mid-ocean-ridge spreading centre and also the evolution of the oceanic lithosphere during its early cooling history. Here we focus on structures imaged at distances approximately 100 to 350 km east of the ridge crest, corresponding to seafloor ages of approximately 1.3 to 4.5 million years (Myr), where the seismic and electrical conductivity structure is nearly constant and independent of age. Beginning at a depth of about 60 km, we image a large increase in electrical conductivity and a change from isotropic to transversely anisotropic electrical structure, with higher conductivity in the direction of fast propagation for seismic waves. Conductive cooling models predict structure that increases in depth with age, extending to about 30 km at 4.5 Myr ago. We infer, however, that the structure of young oceanic plates is instead controlled by a decrease in water content above a depth of 60 km induced by the melting process beneath the spreading centre.

Spreading-rate dependence of melt extraction at mid-ocean ridges from mantle seismic refraction data.

A variety of observations indicate that mid-ocean ridges produce less crust at spreading rates below 20 mm yr(-1) (refs 1-3), reflecting changes in fundamental ridge processes with decreasing spreading rate. The nature of these changes, however, remains uncertain, with end-member explanations being decreasing shallow melting or incomplete melt extraction, each due to the influence of a thicker thermal lid. Here we present results of a seismic refraction experiment designed to study mid-ocean ridge processes by imaging residual mantle structure. Our results reveal an abrupt lateral change in bulk mantle seismic properties associated with a change from slow to ultraslow palaeo-spreading rate. Changes in mantle velocity gradient, basement topography and crustal thickness all correlate with this spreading-rate change. These observations can be explained by variations in melt extraction at the ridge, with a gabbroic phase preferentially retained in the mantle at slower spreading rates. The estimated volume of retained melt balances the approximately 1.5-km difference in crustal thickness, suggesting that changes in spreading rate affect melt-extraction processes rather than total melting.

Trench-parallel anisotropy produced by foundering of arc lower crust.

Many volcanic arcs display fast seismic shear-wave velocities parallel to the strike of the trench. This pattern of anisotropy is inconsistent with simple models of corner flow in the mantle wedge. Although several models, including slab rollback, oblique subduction, and deformation of water-rich olivine, have been proposed to explain trench-parallel anisotropy, none of these mechanisms are consistent with all observations. Instead, small-scale convection driven by the foundering of dense arc lower crust provides an explanation for the trench-parallel anisotropy, even in settings with orthogonal convergence and no slab rollback.