Earthquake faults occur in networks that have dynamical modes not displayed by single isolated faults. Using simulations of the network of strike-slip faults in southern California, we find that the physics depends critically on both the interactions among the faults, which are determined by the geometry of the fault network, as well as on the stress dissipation properties of the nonlinear frictional physics, similar to the dynamics of integrate-and-fire neural networks.

Earthquakes on a specified fault (or fault segment) with magnitudes greater than a specified value have a statistical distribution of recurrence times. The mean recurrence time can be related to the rate of strain accumulation and the strength of the fault. Very few faults have a recorded history of earthquakes that define a distribution well. For hazard assessment, in general, a statistical distribution of recurrence times is assumed along with parameter values. Assumed distributions include the Weibull (stretched exponential) distribution, the lognormal distribution, and the Brownian passage-time (inverse Gaussian) distribution. The distribution of earthquake waiting times is the conditional probability that an earthquake will occur at a time in the future if it has not occurred for a specified time in the past. The distribution of waiting times is very sensitive to the distribution of recurrence times. An exponential distribution of recurrence times is Poissonian, so there is no memory of the last event. The distribution of recurrence times must be thinner than the exponential if the mean waiting time is to decrease as the time since the last earthquake increases. Neither the lognormal or the Brownian passage time distribution satisfies this requirement. We use the "Virtual California" model for earthquake occurrence on the San Andreas fault system to produce a synthetic distribution of earthquake recurrence times on various faults in the San Andreas system. We find that the synthetic data are well represented by Weibull distributions. We also show that the Weibull distribution follows from both damage mechanics and statistical physics.

We discuss the problem of earthquake forecasting in the context of new models for the dynamics based on statistical physics. Here we focus on new, topologically realistic system-level approaches to the modeling of earthquake faults. We show that the frictional failure physics of earthquakes in these complex, topologically realistic models leads to self-organization of the statistical dynamics, and produces statistical distributions characterizing the activity, notably the Gutenberg-Richter magnitude frequency distribution, that are similar to those observed in nature. In particular, we show that a parameterization of friction that includes a simple representation of a dynamic stress intensity factor is needed to organize the dynamics. We also show that the slip distributions for synthetic events obtained in the model are also similar to those observed in nature

No abstract

Virtual California is a topologically realistic simulation of the interacting earthquake faults in California. Inputs to the model arise from field data, and typically include realistic fault system topologies, realistic long-term slip rates, and realistic frictional parameters. Outputs from the simulations include synthetic earthquake sequences and space-time patterns together with associated surface deformation and strain patterns that are similar to those seen in nature. Here we describe details of the data assimilation procedure we use to construct the fault model and to assign frictional properties. In addition, by analyzing the statistical physics of the simulations, we can show that that the frictional failure physics, which includes a simple representation of a dynamic stress intensity factor, leads to selforganization of the statistical dynamics, and produces empirical statistical distributions (probability density functions: PDFs) that characterize the activity. One type of distribution that can be constructed from empirical measurements of simulation data are PDFs for recurrence intervals on selected faults. Inputs to simulation dynamics are based on the use of time-averaged event-frequency data, and outputs include PDFs representing measurements of dynamical variability arising from fault interactions and spacetime correlations. As a first step for productively using model-based methods for earthquake forecasting, we propose that simulations be used to generate the PDFs for recurrence intervals instead of the usual practice of basing the PDFs on standard forms (Gaussian, Log-Normal, Pareto, Brownian Passage Time, and so forth). Subsequent development of simulation-based methods should include model enhancement, data assimilation and data mining methods, and analysis techniques based on statistical physics.