Upper-mantle water stratification inferred from observations of the 2012 Indian Ocean earthquake

Small effect of partial melt on electrical anomalies in the asthenosphere

High conductivity anomalies in the shallow mantle are frequently attributed to minor partial melt (basalt or carbonatite) in the olivine-dominated peridotites. Conductivity of a melt-mineral mixture depends on the configuration of melt that could be affected by grain size of the constitutive mineral(s), but this has rarely been explored. Here, we provide experimental evidence using a conductive carbonatite analog and olivine that the bulk conductivity decreases systematically with increasing olivine grain size. The required amount of melt for producing the geophysically resolved high conductivities in the asthenosphere is much greater than previously assumed. We suggest that the effect of partial melt on many conductive regions in the asthenosphere is small. Instead, the electrical anomalies (especially those away from mid-ocean ridges) originate more likely from subsolidus solid assemblages in the upper mantle. This reconciles well the geochemical and petrological constraints of the shallow mantle with its geophysically determined electrical properties.

An apparatus for measuring nonlinear viscoelasticity of minerals at high temperature

We describe a high-temperature, uniaxial creep apparatus designed to investigate nonlinear attenuation of materials over a wide range of temperatures (25-1300 °C) using forced oscillations combined with a bias stress. This apparatus is primarily designed for investigation of minerals and rocks with high melting temperatures. An oscillatory compressional stress is used to determine attenuation and Young's modulus at frequencies of 10^-1-10² Hz and high stress amplitudes (>0.1 MPa). Large bias stresses are applied in addition to the oscillatory stresses such that attenuation tests are conducted simultaneously with the ongoing creep. The complex compliance of the apparatus was characterized by conducting calibration tests on orientated crystals of sapphire. The real part of the apparatus compliance exhibits a dependence on sample length and frequency, whereas the imaginary part is only dependent on frequency. The complex compliance is not dependent on the oscillation amplitude or the bias stress. We assess the accuracy and precision of this calibration by comparing measurements of the attenuation and Young's modulus of aluminum and acrylic to previously published values. We outline a set of criteria defining the conditions over which this apparatus can precisely determine the attenuation and Young's modulus of a sample based on the sample length and expected values of attenuation and Young's modulus.

Dislocation interactions in olivine control postseismic creep of the upper mantle

Changes in stress applied to mantle rocks, such as those imposed by earthquakes, commonly induce a period of transient creep, which is often modelled based on stress transfer among slip systems due to grain interactions. However, recent experiments have demonstrated that the accumulation of stresses among dislocations is the dominant cause of strain hardening in olivine at temperatures ≤600 °C, raising the question of whether the same process contributes to transient creep at higher temperatures. Here, we demonstrate that olivine samples deformed at 25 °C or 1150–1250 °C both preserve stress heterogeneities of ~1 GPa that are imparted by dislocations and have correlation lengths of ~1 μm. The similar stress distributions formed at these different temperatures indicate that accumulation of stresses among dislocations also provides a contribution to transient creep at high temperatures. The results motivate a new generation of models that capture these intragranular processes and may refine predictions of evolving mantle viscosity over the earthquake cycle.

Models of the viscosity evolution of mantle rocks are central to analyses of postseismic deformation but constraints on underlying physical processes are lacking. Here, the authors present measurements of microscale stress heterogeneity in olivine suggesting that long-range dislocation interactions contribute to viscosity evolution.

Illuminating subduction zone rheological properties in the wake of a giant earthquake

Deformation associated with plate convergence at subduction zones is accommodated by a complex system involving fault slip and viscoelastic flow. These processes have proven difficult to disentangle. The 2010 M _w 8.8 Maule earthquake occurred close to the Chilean coast within a dense network of continuously recording Global Positioning System stations, which provide a comprehensive history of surface strain. We use these data to assemble a detailed picture of a structurally controlled megathrust fault frictional patchwork and the three-dimensional rheological and time-dependent viscosity structure of the lower crust and upper mantle, all of which control the relative importance of afterslip and viscoelastic relaxation during postseismic deformation. These results enhance our understanding of subduction dynamics including the interplay of localized and distributed deformation during the subduction zone earthquake cycle.

Coupled afterslip and transient mantle flow after the 2011 Tohoku earthquake

Modeling of postseismic deformation following great earthquakes has revealed the viscous structure of the mantle and the frictional properties of the fault interface. However, for giant megathrust events, viscoelastic flow and afterslip mechanically interplay with each other during the postseismic period. We explore the role of afterslip and viscoelastic relaxation and their interaction in the aftermath of the 2011 M _w (moment magnitude) 9.0 Tohoku earthquake based on a detailed model analysis of the postseismic deformation with laterally varying, experimentally constrained, rock rheology. Mechanical coupling between viscoelastic relaxation and afterslip notably modifies both the afterslip distribution and surface deformation. Thus, we highlight the importance of addressing mechanical coupling for long-term studies of postseismic relaxation, especially in the context of the geodynamics of the Japan trench across the seismic cycle.

Rapid mantle flow with power-law creep explains deformation after the 2011 Tohoku mega-quake

The deformation transient following large subduction zone earthquakes is thought to originate from the interaction of viscoelastic flow in the asthenospheric mantle and slip on the megathrust that are both accelerated by the sudden coseismic stress change. Here, we show that combining insight from laboratory solid-state creep and friction experiments can successfully explain the spatial distribution of surface deformation in the first few years after the 2011 M_w 9.0 Tohoku-Oki earthquake. The transient reduction of effective viscosity resulting from dislocation creep in the asthenosphere explains the peculiar retrograde displacement revealed by seafloor geodesy, while the slip acceleration on the megathrust accounts for surface displacements on land and offshore outside the rupture area. Our results suggest that a rapid mantle flow takes place in the asthenosphere with temporarily decreased viscosity in response to large coseismic stress, presumably due to the activation of power-law creep during the post-earthquake period.

Lower-crustal rheology and thermal gradient in the Taiwan orogenic belt illuminated by the 1999 Chi-Chi earthquake

The strength of the lithosphere controls tectonic evolution and seismic cycles, but how rocks deform under stress in their natural settings is usually unclear. We constrain the rheological properties beneath the Taiwan orogenic belt using the stress perturbation following the 1999 Chi-Chi earthquake and fourteen-year postseismic geodetic observations. The evolution of stress and strain rate in the lower crust is best explained by a power-law Burgers rheology with rapid increases in effective viscosities from ~10¹⁷ to ~10¹⁹ Pa s within a year. The short-term modulation of the lower-crustal strength during the seismic cycle may alter the energy budget of mountain building. Incorporating the laboratory data and associated uncertainties, inferred thermal gradients suggest an eastward increase from 19.5±2.5°C/km in the Coastal Plain to 32±3°C/km in the Central Range. Geodetic observations may bridge the gap between laboratory and lithospheric scales to investigate crustal rheology and tectonic evolution.

Dynamics of fault motion and the origin of contrasting tectonic style between Earth and Venus

Plate tectonics is one mode of mantle convection that occurs when the surface layer (the lithosphere) is relatively weak. When plate tectonics operates on a terrestrial planet, substantial exchange of materials occurs between planetary interior and its surface. This is likely a key in maintaining the habitable environment on a planet. Therefore it is essential to understand under which conditions plate tectonics operates on a terrestrial planet. One of the puzzling observations in this context is the fact that plate tectonics occurs on Earth but not on Venus despite their similar size and composition. Factors such as the difference in water content or in grain-size have been invoked, but these models cannot easily explain the contrasting tectonic styles between Earth and Venus. We propose that strong dynamic weakening in friction is a key factor. Fast unstable fault motion is found in cool Earth, while slow and stable fault motion characterizes hot Venus, leading to substantial dynamic weakening on Earth but not on Venus. Consequently, the tectonic plates are weak on Earth allowing for their subduction, while the strong plates on Venus promote the stagnant lid regime of mantle convection.

Transient rheology of the Sumatran mantle wedge revealed by a decade of great earthquakes

Understanding the rheological properties of the upper mantle is essential to develop a consistent model of mantle dynamics and plate tectonics. However, the spatial distribution and temporal evolution of these properties remain unclear. Here, we infer the rheological properties of the asthenosphere across multiple great megathrust earthquakes between 2004 and 2014 along the Sumatran subduction zone, taking advantage of decade-long continuous GPS and tide-gauge measurements. We observe transient mantle wedge flow following these earthquakes, and infer the temporal evolution of the effective viscosity. We show that the evolution of stress and strain rate following these earthquakes is better matched by a bi-viscous than by a power-law rheology model, and we estimate laterally heterogeneous transient and background viscosities on the order of ~10¹⁷ and ~10¹⁹ Pa s, respectively. Our results constitute a preliminary rheological model to explain stress evolution within earthquake cycles and the development of seismic hazard in the region.

The rheology of the upper mantle is key to understanding how plate tectonics may evolve. Here, using GPS and tide-gauge measurements along the Sumatran subduction zone, the authors’ show that a bi-viscous rheology model is needed to explain the stress and strain evolution of the upper mantle following earthquakes.

Imaging the distribution of transient viscosity after the 2016 Mw 7.1 Kumamoto earthquake

The deformation of mantle and crustal rocks in response to stress plays a crucial role in the distribution of seismic and volcanic hazards, controlling tectonic processes ranging from continental drift to earthquake triggering. However, the spatial variation of these dynamic properties is poorly understood as they are difficult to measure. We exploited the large stress perturbation incurred by the 2016 earthquake sequence in Kumamoto, Japan, to directly image localized and distributed deformation. The earthquakes illuminated distinct regions of low effective viscosity in the lower crust, notably beneath the Mount Aso and Mount Kuju volcanoes, surrounded by larger-scale variations of viscosity across the back-arc. This study demonstrates a new potential for geodesy to directly probe rock rheology in situ across many spatial and temporal scales.