Search for: All records

Abstract Numerical simulations of Sequences of Earthquakes and Aseismic Slip (SEAS) have rapidly progressed to address fundamental problems in fault mechanics and provide self‐consistent, physics‐based frameworks to interpret and predict geophysical observations across spatial and temporal scales. To advance SEAS simulations with rigor and reproducibility, we pursue community efforts to verify numerical codes in an expanding suite of benchmarks. Here we present code comparison results from a new set of quasi‐dynamic benchmark problems BP6‐QD‐A/S/C that consider an aseismic slip transient induced by changes in pore fluid pressure consistent with fluid injection and diffusion in fault models with different treatments of fault friction. Ten modeling groups participated in problems BP6‐QD‐A and BP6‐QD‐S considering rate‐and‐state fault models using the aging (‐A) and slip (‐S) law formulations for frictional state evolution, respectively, allowing us to better understand how various computational factors across codes affect the simulated evolution of pore pressure and aseismic slip. Comparisons of problems using the aging versus slip law, and a constant friction coefficient (‐C), illustrate how aseismic slip models can differ in the timing and amount of slip achieved with different treatments of fault friction given the same perturbations in pore fluid pressure. We achieve excellent quantitative agreement across participating codes, with further agreement attained by ensuring sufficiently fine time‐stepping and consistent treatment of boundary conditions. Our benchmark efforts offer a community‐based example to reveal sensitivities of numerical modeling results, which is essential for advancing multi‐physics SEAS models to better understand and construct reliable predictive models of fault dynamics. 

Abstract We examine a simple mechanism for the spatiotemporal evolution of transient, slow slip. We consider the problem of slip on a fault that lies within an elastic continuum and whose strength is proportional to sliding rate. This rate dependence may correspond to a viscously deforming shear zone or the linearization of a nonlinear, rate‐dependent fault strength. We examine the response of such a fault to external forcing, such as local increases in shear stress or pore fluid pressure. We show that the slip and slip rate are governed by a type of diffusion equation, the solution of which is found using a Green's function approach. We derive the long‐time, self‐similar asymptotic expansion for slip or slip rate, which depend on both timetand a similarity coordinateη = x/t, wherexdenotes fault position. The similarity coordinate shows a departure from classical diffusion and is owed to the nonlocal nature of elastic interaction among points on an interface between elastic half‐spaces. We demonstrate the solution and asymptotic analysis of several example problems. Following sudden impositions of loading, we show that slip rate ultimately decays as 1/twhile spreading proportionally tot, implying both a logarithmic accumulation of displacement and a constant moment rate. We discuss the implication for models of postseismic slip as well as spontaneously emerging slow slip events.