A mechanism for preseismic steady rupture fronts observed in laboratory experiments

It has been shown that the onset of frictional instability is characterized by a transition from stable, quasi-static rupture growth to unstable, inertially-controlled high-speed rupture. In particular, slow rupture fronts propagating at a steady speed V_(slow) of the order of 5% of the S-wave speed have been observed prior to the onset of dynamic rupture in recent fault-friction laboratory experiments. However, the precise mechanism governing this V_(slow) stage is unknown. Here we reproduce this phenomenon in numerical simulations of earthquake sequences that incorporate laboratory-derived rate-and-state friction laws. Our simulations show that the V_(slow) stage originates from a stress concentration inherited from the coalescence of interseismic slow creep fronts. Its occurrence is limited to a narrow range of the parameter space but is found in simulations with two commonly-used state-variable evolution laws in the rate-and-state formulation. The sensitivity of the speed V_(slow) to the model parameters suggests that the propagation speed V_(slow) reported in laboratory experiments may also be sensitive to parameters of friction and stress conditions. Our results imply that time and space dimensions associated with the propagation of V_(slow) on natural faults can be as much as a few seconds and several hundred meters, respectively. Hence the detection of such preseismic signals may be possible with near-field high-resolution observations

Linear response for spiking neuronal networks with unbounded memory

We establish a general linear response relation for spiking neuronal networks, based on chains with unbounded memory. This relation allows us to predict the influence of a weak amplitude time-dependent external stimuli on spatio-temporal spike correlations, from the spontaneous statistics (without stimulus) in a general context where the memory in spike dynamics can extend arbitrarily far in the past. Using this approach, we show how linear response is explicitly related to neuronal dynamics with an example, the gIF model, introduced by M. Rudolph and A. Destexhe. This example illustrates the collective effect of the stimuli, intrinsic neuronal dynamics, and network connectivity on spike statistics. We illustrate our results with numerical simulations.Comment: 60 pages, 8 figure

Ambiguity of the Moment Tensor

An earthquake on a fault separating two dissimilar materials does not have a well-defined moment density tensor. We present a complete characterization of this bimaterial ambiguity in the general case of slip on a fault in an anisotropic medium. The ambiguity can be eliminated by utilizing a potency density rather than a moment density representation of a bimaterial source

Self-similar slip pulses during rate-and-state earthquake nucleation

For a wide range of conditions, earthquake nucleation zones on rate- and state-dependent faults that obey either of the popular state evolution laws expand as they accelerate. Under the “slip” evolution law, which experiments show to be the more relevant law for nucleation, this expansion takes the form of a unidirectional slip pulse. In numerical simulations these pulses often tend to approach, with varying degrees of robustness, one of a few styles of self-similar behavior. Here we obtain an approximate self-similar solution that accurately describes slip pulses growing into regions initially sliding at steady state. In this solution the length scale over which slip speeds are significant continually decreases, being inversely proportional to the logarithm of the maximum slip speed V_(max), while the total slip remains constant. This slip is close to D_c(1−a/b)^(−1), where D_c is the characteristic slip scale for state evolution and a and b are the parameters that determine the sensitivity of the frictional strength to changes in slip rate and state. The pulse has a “distance to instability” as well as a “time to instability,” with the remaining propagation distance being proportional to (1−a/b)^(−2) [ln(V_(max)Θ_(bg)/D_c)]^(−1), where Θ_(bg) is the background state into which the pulse propagates. This solution provides a reasonable estimate of the total slip for pulses growing into regions that depart modestly from steady state

Spectral-element simulations of long-term fault slip: Effect of low-rigidity layers on earthquake-cycle dynamics

We develop a spectral element method for the simulation of long-term histories of spontaneous seismic and aseismic slip on faults subjected to tectonic loading. Our approach reproduces all stages of earthquake cycles: nucleation and propagation of earthquake rupture, postseismic slip and interseismic creep. We apply the developed methodology to study the effects of low-rigidity layers on the dynamics of the earthquake cycle in 2-D. We consider two cases: small (M ~ 1) earthquakes on a fault surrounded by a damaged fault zone and large (M ~ 7) earthquakes on a vertical strike-slip fault that cuts through shallow low-rigidity layers. Our results indicate how the source properties of repeating earthquakes are affected by the presence of a damaged fault zone with low rigidity. Compared to faults in homogeneous media, we find (1) reduction in the earthquake nucleation size, (2) amplification of slip rates during dynamic rupture propagation, (3) larger recurrence interval, and (4) smaller amount of aseismic slip. Based on linear stability analysis, we derive a theoretical estimate of the nucleation size as a function of the width and rigidity reduction of the fault zone layer, which is in good agreement with simulated nucleation sizes. We further examine the effects of vertically-stratified layers (e.g., sedimentary basins) on the nature of shallow coseismic slip deficit. Our results suggest that low-rigidity shallow layers alone do not lead to coseismic slip deficit. While the low-rigidity layers result in lower interseismic stress accumulation, they also cause dynamic amplification of slip rates, with the net effect on slip being nearly zero

A window into the complexity of the dynamic rupture of the 2011 Mw 9 Tohoku-Oki earthquake

The 2011 Mw 9 Tohoku-Oki earthquake, recorded by over 1000 near-field stations and multiple large-aperture arrays, is by far the best recorded earthquake in the history of seismology and provides unique opportunities to address fundamental issues in earthquake source dynamics. Here we conduct a high resolution array analysis based on recordings from the USarray and the European network. The mutually consistent results from both arrays reveal rupture complexity with unprecedented resolution, involving phases of diverse rupture speed and intermittent high frequency bursts within slow speed phases, which suggests spatially heterogeneous material properties. The earthquake initially propagates down-dip, with a slow initiation phase followed by sustained propagation at speeds of 3 km/s. The rupture then slows down to 1.5 km/s for 60 seconds. A rich sequence of bursts is generated along the down-dip rim of this slow and roughly circular rupture front. Before the end of the slow phase an extremely fast rupture front detaches at about 5 km/s towards the North. Finally a rupture front propagates towards the south running at about 2.5 km/s for over 100 km. Key features of the rupture process are confirmed by the strong motion data recorded by K-net and KIK-net. The energetic high frequency radiation episodes within a slow rupture phase suggests a patchy image of the brittle-ductile transition zone, composed of discrete brittle asperities within a ductile matrix. The high frequency is generated mainly at the down-dip edge of the principal slip regions constrained by geodesy, suggesting a variation along dip of the mechanical properties of the mega thrust fault or their spatial heterogeneity that affects rise time

Spectral-element modeling of spontaneous earthquake rupture on rate and state faults: Effect of velocity-strengthening friction at shallow depths

We develop a spectral-element methodology (SEM) for simulating dynamic rupture on rate and state faults and use it to study how the rupture is affected by a shallow fault region of steady-state velocity-strengthening friction. Our comparison of the developed SEM and a spectral boundary-integral method (BIM) for an anti-plane (two-dimensional) test problem shows that the two methods produce virtually identical solutions for the finest resolution we use and that the convergence with grid reduction of the developed SEM methodology is comparable to that of BIM. We also use the test problem to compare numerical resolution required for different state evolution laws and for linear slip-weakening friction. Using our three-dimensional implementation of the methodology, we find that a shallow velocity-strengthening fault region can significantly alter dynamic rupture and ground motion. The velocity-strengthening region suppresses supershear propagation at the free surface occurring in the absence of such region, which could explain the lack of universally observed supershear rupture near the free surface. In addition, the velocity-strengthening region promotes faster fall-off of slip velocity behind the rupture front and decreases final slip throughout the entire fault, causing a smaller average stress drop. The slip decrease is largest in the shallow parts of the fault, resulting in the depth profile of slip qualitatively consistent with observations of shallow co-seismic slip deficit. The shallow velocity-strengthening region also reduces the amplification of strong ground motion due to a low-velocity bulk structure

Seismic Radiation From Simple Models of Earthquakes

We review some basic features of shear wave generation and energy balance for a 2D anti plane rupture. We first study the energy balance for a flat fault, and for a fault that contains a single localized kink. We determine an exact expression for the partition between strain energy flow released from the elastic medium surrounding the fault, radiated energy flow and energy release rate. This balance depends only on the rupture speed and the residual stress intensity factor. When the fault contains a kink, the energy available for fracture is reduced so that the rupture speed is reduced. When rupture speed changes abruptly, the radiated energy flow also changes abruptly. As rupture propagates across the kink, a shear wave is emitted that has a displacement spectral content that decreases like ω^(-2) at high frequencies. We then use spectral elements to model the propagation of an antiplane crack with a slip-weakening friction law. Since the rupture front in this case has a finite length scale, the wave emitted by the kink is smoothed at very high frequencies but its general behavior is similar to that predicted by the simple sharp crack model. A model of a crack that has several kinks and wanders around a mean rupture directions, shows that kinks reduce the rupture speed along the average rupture direction of the fault. Contrary to flat fault models, a fault with kinks produces high frequency waves that are emitted every time the rupture front turns at a kink. Finally, we discuss the applicability of the present results to a 3D rupture model

Deep ductile shear localization facilitates near-orthogonal strike-slip faulting in a thin brittle lithosphere

Some active fault systems comprise near-orthogonal conjugate strike-slip faults, as highlighted by the 2019 Ridgecrest and the 2012 Indian Ocean earthquake sequences. In conventional failure theory, orthogonal faulting requires a pressure-insensitive rock strength, which is unlikely in the brittle lithosphere. Here, we conduct 3D numerical simulations to test the hypothesis that near-orthogonal faults can form by inheriting the geometry of deep ductile shear bands. Shear bands nucleated in the deep ductile layer, a pressure-insensitive material, form at 45 degree from the maximum principal stress. As they grow upwards into the brittle layer, they progressively rotate towards the preferred brittle faulting angle, ~30 degree, forming helical shaped faults. If the brittle layer is sufficiently thin, the rotation is incomplete and the near-orthogonal geometry is preserved at the surface. The preservation is further facilitated by a lower confining pressure in the shallow portion of the brittle layer. For this inheritance to be effective, a thick ductile fault root beneath the brittle layer is necessary. The model offers a possible explanation for orthogonal faulting in Ridgecrest, Salton Trough, and Wharton basin. Conversely, faults nucleated within the brittle layer form at the optimal angle for brittle faulting and can cut deep into the ductile layer before rotating to 45 degree. Our results thus reveal the significant interactions between the structure of faults in the brittle upper lithosphere and their deep ductile roots

Properties of Dynamic Earthquake Ruptures With Heterogeneous Stress Drop

Earthquake rupture is a notoriously complex process, at all observable scales. We introduce a simplified semi-dynamic crack model to investigate the connection between the statistical properties of stress and those of macroscopic source parameters such as rupture size, seismic moment, apparent stress drop and radiated energy. Rupture initiation is treated consistently with nucleation on a linear slip-weakening fault, whereas rupture propagation and arrest are treated according to the Griffith criterion. The available stress drop is prescribed as a spatially correlated random field and is shown to potentially sustain a broad range of magnitudes. By decreasing the amplitude of the stress heterogeneities or increasing their correlation length the distribution of earthquake sizes presents a transition from Gutenberg- Richter to characteristic earthquake behavior. This transition is studied through a mean-field analysis. The bifurcation to characteristic earthquake behavior is sharp, reminiscent of a first-order phase transition. A lower roll-off magnitude observed in the Gutenberg-Richter regime is shown to depend on the correlation length of the available stress drop, rather than being a direct signature of the nucleation process. More generally, we highlight the possible role of the stress correlation length scale on deviations from earthquake source self-similarity. The present reduced model is a building block towards understanding the effect of structural and dynamic fault heterogeneities on the scaling of source parameters and on basic properties of seismicity