Modeling Anisotropy in the Earthquake Cycle: a Numerical Multiscale Model

We have developed a methodology for incorporating and studying the effects of anisotropy when simulating the full earthquake cycle. The method is developed for a vertical strike-slip fault in two-dimensions, with antiplane motion. Inertial terms are dropped from the elastic anisotropic wave equation to obtain a steady state problem. This resulting equilibrium equation is discretized with a finite difference method. A nonlinear rate-and-state friction law is enforced at the fault. Time stepping is adaptive to capture highly varying time scales, and as such is able to produce self-consistent initial conditions

Community Code Verification Exercise for Simulating Sequences of Earthquakes and Aseismic Slip (SEAS)

Numerical simulations of sequences of earthquakes and aseismic slip (SEAS) have made great progress over past decades to address important questions in earthquake physics. However, significant challenges in SEAS modeling remain in resolving multiscale interactions between earthquake nucleation, dynamic rupture, and aseismic slip, and understanding physical factors controlling observables such as seismicity and ground deformation. The increasing complexity of SEAS modeling calls for extensive efforts to verify codes and advance these simulations with rigor, reproducibility, and broadened impact. In 2018, we initiated a community code‐verification exercise for SEAS simulations, supported by the Southern California Earthquake Center. Here, we report the findings from our first two benchmark problems (BP1 and BP2), designed to verify different computational methods in solving a mathematically well‐defined, basic faulting problem. We consider a 2D antiplane problem, with a 1D planar vertical strike‐slip fault obeying rate‐and‐state friction, embedded in a 2D homogeneous, linear elastic half‐space. Sequences of quasi‐dynamic earthquakes with periodic occurrences (BP1) or bimodal sizes (BP2) and their interactions with aseismic slip are simulated. The comparison of results from 11 groups using different numerical methods show excellent agreements in long‐term and coseismic fault behavior. In BP1, we found that truncated domain boundaries influence interseismic stressing, earthquake recurrence, and coseismic rupture, and that model agreement is only achieved with sufficiently large domain sizes. In BP2, we found that complexity of fault behavior depends on how well physical length scales related to spontaneous nucleation and rupture propagation are resolved. Poor numerical resolution can result in artificial complexity, impacting simulation results that are of potential interest for characterizing seismic hazard such as earthquake size distributions, moment release, and recurrence times. These results inform the development of more advanced SEAS models, contributing to our further understanding of earthquake system dynamics

Community Code Verification Exercise for Simulating Sequences of Earthquakes and Aseismic Slip (SEAS)

Numerical simulations of sequences of earthquakes and aseismic slip (SEAS) have made great progress over past decades to address important questions in earthquake physics. However, significant challenges in SEAS modeling remain in resolving multiscale interactions between earthquake nucleation, dynamic rupture, and aseismic slip, and understanding physical factors controlling observables such as seismicity and ground deformation. The increasing complexity of SEAS modeling calls for extensive efforts to verify codes and advance these simulations with rigor, reproducibility, and broadened impact. In 2018, we initiated a community code‐verification exercise for SEAS simulations, supported by the Southern California Earthquake Center. Here, we report the findings from our first two benchmark problems (BP1 and BP2), designed to verify different computational methods in solving a mathematically well‐defined, basic faulting problem. We consider a 2D antiplane problem, with a 1D planar vertical strike‐slip fault obeying rate‐and‐state friction, embedded in a 2D homogeneous, linear elastic half‐space. Sequences of quasi‐dynamic earthquakes with periodic occurrences (BP1) or bimodal sizes (BP2) and their interactions with aseismic slip are simulated. The comparison of results from 11 groups using different numerical methods show excellent agreements in long‐term and coseismic fault behavior. In BP1, we found that truncated domain boundaries influence interseismic stressing, earthquake recurrence, and coseismic rupture, and that model agreement is only achieved with sufficiently large domain sizes. In BP2, we found that complexity of fault behavior depends on how well physical length scales related to spontaneous nucleation and rupture propagation are resolved. Poor numerical resolution can result in artificial complexity, impacting simulation results that are of potential interest for characterizing seismic hazard such as earthquake size distributions, moment release, and recurrence times. These results inform the development of more advanced SEAS models, contributing to our further understanding of earthquake system dynamics

A Computational Method for Earthquake Cycles within Anisotropic Media

We present a numerical method for the simulation of earthquake cycles on a 1-D fault interface embedded in a 2-D homogeneous, anisotropic elastic solid. The fault is governed by an experimentally motivated friction law known as rate-and-state friction which furnishes a set of ordinary differential equations which couple the interface to the surrounding volume. Time enters the problem through the evolution of the ordinary differential equations along the fault and provides boundary conditions for the volume, which is governed by quasi-static elasticity. We develop a time-stepping method which accounts for the interface/volume coupling and requires solving an elliptic partial differential equation for the volume response at each time step. The 2-D volume is discretized with a second-order accurate finite difference method satisfying the summation-by-parts property, with boundary and fault interface conditions enforced weakly. This framework leads to a provably stable semi-discretization. To mimic slow tectonic loading, the remote side-boundaries are displaced at a slow rate, which eventually leads to earthquake nucleation at the fault. Time stepping is based on an adaptive, fourth-order Runge–Kutta method and captures the highly varying timescales present. The method is verified with convergence tests for both the orthotropic and fully anisotropic cases. An initial parameter study reveals regions of parameter space where the systems experience a bifurcation from period one to period two behaviour. Additionally, we find that anisotropy influences the recurrence interval between earthquakes, as well as the emergence of aseismic transients and the nucleation zone size and depth of earthquakes

The Community Code Verification Exercise for Simulating Sequences of Earthquakes and Aseismic Slip (SEAS)

International audienceAbstract Numerical simulations of sequences of earthquakes and aseismic slip (SEAS) have made great progress over past decades to address important questions in earthquake physics. However, significant challenges in SEAS modeling remain in resolving multiscale interactions between earthquake nucleation, dynamic rupture, and aseismic slip, and understanding physical factors controlling observables such as seismicity and ground deformation. The increasing complexity of SEAS modeling calls for extensive efforts to verify codes and advance these simulations with rigor, reproducibility, and broadened impact. In 2018, we initiated a community code-verification exercise for SEAS simulations, supported by the Southern California Earthquake Center. Here, we report the findings from our first two benchmark problems (BP1 and BP2), designed to verify different computational methods in solving a mathematically well-defined, basic faulting problem. We consider a 2D antiplane problem, with a 1D planar vertical strike-slip fault obeying rate-and-state friction, embedded in a 2D homogeneous, linear elastic half-space. Sequences of quasi-dynamic earthquakes with periodic occurrences (BP1) or bimodal sizes (BP2) and their interactions with aseismic slip are simulated. The comparison of results from 11 groups using different numerical methods show excellent agreements in long-term and coseismic fault behavior. In BP1, we found that truncated domain boundaries influence interseismic stressing, earthquake recurrence, and coseismic rupture, and that model agreement is only achieved with sufficiently large domain sizes. In BP2, we found that complexity of fault behavior depends on how well physical length scales related to spontaneous nucleation and rupture propagation are resolved. Poor numerical resolution can result in artificial complexity, impacting simulation results that are of potential interest for characterizing seismic hazard such as earthquake size distributions, moment release, and recurrence times. These results inform the development of more advanced SEAS models, contributing to our further understanding of earthquake system dynamics