Evolution of stress deficit and changing rates of seismicity in cellular automaton models of earthquake faults

We investigate the internal dynamics of two cellular automaton models with heterogeneous strength fields and differing nearest neighbour laws. One mod el is a crack-like automaton, transferring ail stress from a rupture zone t o the surroundings. The other automaton is a partial stress drop automaton, transferring only a fraction of the stress within a rupture zone to the su rroundings. To study evolution of stress, the mean spectral density. f(k(r) ) of a stress deficit held is: examined prior to, and immediately following ruptures in both models. Both models display a power-law relationship betw een f(k(r)) and spatial wavenumber (k(r)) of the form f(k(r)) similar tok(r )(-beta). In the crack model, the evolution of stress deficit is consistent with cyclic approach to, and retreat from a critical state in which large events occur. The approach to criticality is driven by tectonic loading. Sh ort-range stress transfer in the model does not affect the approach to crit icality of broad regions in the model. The evolution of stress deficit in t he partial stress drop model is consistent with small fluctuations about a mean state of high stress, behaviour indicative of a self-organised critica l system. Despite statistics similar to natural earthquakes these simplifie d models lack a physical basis. physically motivated models of earthquakes also display dynamical complexity similar to that of a critical point syste m. Studies of dynamical complexity in physical models of earthquakes may le ad to advancement towards a physical theory for earthquakes.