Dynamic triggering of earthquakes is promoted by crustal heterogeneities and bimaterial faults

Remotely triggered earthquakes and aftershocks constitute a great challenge in assessing seismic risk. A growing body of observations indicates that significant earthquakes can be triggered by moderate to great earthquakes occurring at distances of up to thousands of kilometres. Currently we lack the knowledge to predict the location of triggered events. We present numerical simulations showing that dynamic interactions between material heterogeneities (e.g. compliant fault zones, sedimentary basins) and seismic waves focus and enhance stresses sufficiently to remotely trigger earthquakes. Numerical simulations indicate that even at great distances (>100 km), the amplified transient dynamic stress near heterogeneities is equivalent to stress levels near the source rupture tip

Analysis of slip-weakening frictional laws with static restrengthening and their implications on the scaling, asymmetry and mode of dynamic rupture on homogeneous and bi-material interfaces

Dynamic simulations of homogeneous and heterogeneous fault rupture using the finite element method are presented giving rise to both crack-like and pulse-like rupture. We employ various slip-weakening frictional laws to examine their effect on the resulting earthquake rupture speed, size and mode. More complex rupture characteristics were produced with more strongly slip-weakening frictional laws, and the degree of slip-weakening had to be finely tuned to reproduce realistic earthquake rupture characteristics. Rupture propagation on a fault is controlled by the constitutive properties of the fault. A dynamic elasto-plastic constitutive law for the interface friction at the fault is formulated based on the Coulomb failure criterion and applied in a way analogous to non-associated elasto-plasticity. We provide benchmark tests of our method against other reported solutions in the literature. We demonstrate the applicability of our elasto-plastic fault model for modeling dynamic rupture and wave propagation in fault systems, and the rich array of dynamic properties produced by our elasto-plastic finite element fault model. These are governed by a number of model parameters including: the spatial and material heterogeneity of the fault, the fault strength, and not least of all the frictional law employed. Asymmetric bilateral fault rupture was produced for the heterogeneous case, where the degree of heterogeneity influenced the rupture speed in the different propagation directions

esys User's guide: Solving partial differential equations with Escript and Finley. Release 3.0 (r2601)

esys.escript is a python-based environment for implementing mathematical models, in particular those based on coupled, non-linear, time-dependent partial differential equations. It consists of four major components: • esys.escript core library • finite element solver esys.finley (which uses fast vendor-supplied solvers or our paso linear solver library) • the meshing interface esys.pycad • a model library. The current version supports parallelization through both MPI for distributed memory and OpenMP for distributed shared memory. The esys.pyvisi module from previous releases has been deprecated. For more info on this and other changes from previous releases see Appendix A.2. If you use this software in your research, then we would appreciate (but do not require) a citation. Some relevant references can be found in Appendix A.3

Identification of supershear transition mechanisms due to material contrast at bimaterial faults

Numerical modelling of dynamic rupture is conducted along faults separating similar and dissimilar materials. Supershear transition is enhanced in the direction of slip of the stiffer material (the negative direction) due to the bimaterial effect whereby a decrease in normal stress in front of the crack tip supports yielding ahead of the rupture. In the direction of slip of the more compliant material (the positive direction), an increase in normal stress ahead of the rupture tip delays or prevents the supershear transition, whereas the impact of the bimaterial effect on subshear ruptures is to promote rupture in the positive direction due to the tensile stress perturbation behind the rupture tip in this direction. We demonstrate that the material contrast and the parameter S control whether the transition from sub- to supershear velocity (supershear transition) is smooth or follows the Burridge–Andrews mechanism. Supershear transition along interfaces separating dissimilar materials is possible for higher values of the parameter S than supershear transition along material interfaces separating similar materials. The difference between pulse-like and crack-like rupture is small with regard to the supershear transition type

Stress heterogeneities in earthquake rupture experiments with material contrasts

We investigate significant heterogeneous stresses along bimaterial interfaces in laboratory and numerical experiments. These stresses, partially induced by model or experimental configuration, affect the supershear transition length and rupture speed, mode and directivity in uniaxial compression tests and dynamic rupture experiments with bimaterial interfaces. Using numerical simulations we show that normal and tangential stresses at the fault are distorted by the different stress-strain relationships of the materials. This distortion leads to altered supershear transition lengths, higher rupture potencies and amplifies the preference for rupture in the direction of slip of the slower and more compliant material. We demonstrate how this stress-distortion can be decreased in laboratory experiments by using larger specimen samples and in numerical models by using periodic boundary conditions

Bridging the macro to mesoscale: evaluating the fourth-order anisotropic damage tensor parameters from ultrasonic measurements of an isotropic solid under triaxial stress loading

One of the most challenging problems which arises in continuum damage mechanics is the selection of variables to describe the internal damage. Many theories have been proposed and various types of damage variables ranging from scalar to vector to tensor quantities have been used. In this paper we consider anisotropic damage and the most general form for damage by using a fourth-order tensor for the damage variables. We demonstrate how experimentally measured quantities can be related to the internal tensorial damage variables. We apply this analysis to experiments of an initially isotropic solid becoming transverse isotropic under triaxial or uniaxial stress loading

A study of localization limiters and mesh dependency in earthquake rupture

No complete physically consistent model of earthquake rupture exists that can fully describe the rich hierarchy of scale dependencies and nonlinearities associated with earthquakes. We study mesh sensitivity in numerical models of earthquake rupture and demonstrate that this mesh sensitivity may provide hidden clues to the underlying physics generating the rich dynamics associated with earthquake rupture. We focus on unstable slip events that occur in earthquakes when rupture is associated with frictional weakening of the fault. Attempts to simulate these phenomena directly by introducing the relevant constitutive behaviour leads to mesh-dependent results, where the deformation localizes in one element, irrespective of size. Interestingly, earthquake models with oversized mesh elements that are ill-posed in the continuum limit display more complex and realistic physics. Until now, the mesh-dependency problem has been regarded as a red herring—but have we overlooked an important clue arising from the mesh sensitivity? We analyse spatial discretization errors introduced into models with oversized meshes to show how the governing equations may change because of these error terms and give rise to more interesting physics. © 2010 The Royal Society

Simulating damage evolution and fracture propagation in sandstone containing a preexisting 3-D surface flaw under uniaxial compression

Preexisting flaws and rock heterogeneity have important ramifications on the process of rock fracturing and on rock stability in many applications. Therefore, there is great interest in numerical modelling of rock fracture and the underlying mechanisms. We simulated damage evolution and fracture propagation in sandstone specimens containing a preexisting 3-D surface flaw under uniaxial compression. We applied the linear elastic damage model based on the unified strength theory following the rock failure process analysis code. However, in contrast to the rock failure process analysis code, we used the finite element method with tetrahedron elements on unstructured meshes. It provided higher geometrical flexibility and allowed for a more accurate representation of the disk-shaped flaw with various flaw depths, angles, and lengths through locally adapted meshes. The rock heterogeneity was modelled by sampling the initial local Young's modulus from a Weibull distribution over a cubic grid. The values were then interpolated to the computational finite element method mesh. This method introduced an additional length scale for the rock heterogeneity represented by the cell size in the sampling grid. The generation of three typical surface cracking patterns, called wing cracks, anti-wing cracks, and far-field cracks, were identified in the simulation results. These depend on the geometry of the preexisting surface flaw. The simulated fracture propagation, coalescence types, and failure modes for the specimens with preexisting surface flaw show good agreement with recent experimental studies