We study numerically the effects of fault roughness on the nucleation process during earthquake sequences. The faults are governed by a rate and state friction law. The roughness introduces local barriers that complicate the nucleation process and result in asymmetric expansion of the rupture, nonmonotonic increase in the slip rates on the fault, and the generation of multiple slip pulses. These complexities are reflected as irregular fluctuations in the moment rate. There is a large difference between first slip events in the sequences and later events. In the first events, for roughness amplitude br ≤ 0.002, there is a large increase in the nucleation length with increasing br. For larger values of br, slip is mostly aseismic. For the later events there is a trade‐off between the effects of the finite fault length and the fault roughness. For br ≤ 0.002, the finite length is a more dominant factor and the nucleation length barely changes with br. For larger values of br, the roughness plays a larger role and the nucleation length increases significantly with br. Using an energy balance approach, where the roughness is accounted for in the fault stiffness, we derive an approximate solution for the nucleation length on rough faults. The solution agrees well with the main trends observed in the simulations for the later events and provides an estimate of the frictional and roughness properties under which faults experience a transition between seismic and aseismic slip.

To observe and quantify the production of microfracturing from initial yield to failure, we deformed Carrara marble samples in uniaxial compression at 20, 105, and 180°C and continuously observed a region of about 1 mm2 on an exposed face with a long‐working distance microscope. Using image processing and microscale strain‐mapping techniques, we measured local strains over a length scale of tens of micrometers. By treating the images with various filters, we identified linear damage features, as well as the magnitude of localized strain and the mode of deformation, i.e., shear versus normal deformation. In general, shear deformation is more prevalent after initial yielding, while tensile deformation dominates closer to peak stress. Independent measurements of both stress and microcrack density at different stages of each experiment provide a unique opportunity to explicitly compare the data with damage models. The model of Ashby and Sammis (1990) significantly underestimated the damage that the rock could sustain before peak stress, perhaps owing to the influence of weak grain boundaries on the damage production. In these samples, microcracks tended to form near boundaries before yield stress. During strain hardening, the damage parameters increased rapidly as longer microcracks grew along the boundaries and finally transected grains as loading neared peak stress. The microcrack density can be empirically related to the reduction of Young's modulus; stiffness ratios decay exponentially with increasing microcrack density for T ≤ 105°C.

Large, destructive earthquakes often propagate along thrust faults including megathrusts. The asymmetric interaction of thrust earthquake ruptures with the free surface leads to sudden variations in fault-normal stress, which affect fault friction. Here, we present full-field experimental measurements of displacements, particle velocities, and stresses that characterize the rupture interaction with the free surface, including the large normal stress reductions. We take advantage of these measurements to investigate the dependence of dynamic friction on transient changes in normal stress, demonstrate that the shear frictional resistance exhibits a significant lag in response to such normal stress variations, and identify a predictive frictional formulation that captures this effect. Properly accounting for this delay is important for simulations of fault slip, ground motion, and associated tsunami excitation.

No abstract