Although models of homogeneous faults de-velop seismicity that has a Gutenberg-Richter distribution, this is only a transient state that is followed by events that are strongly influenced by the nature of the boundaries. Models with geometrical inhomogeneities of fracture thresholds can limit the sizes of earthquakes but now favor the characteristic earthquake model for large earthquakes. Before discussing recent progress toward our goal of studying simulated seismicity on a network of faults (3), we first summarize some of the relevant results in the modeling of seismicity on a single fault as a preliminary step toward the more difficult problem. The focus on simulations of the power-law part of the distribution has led a number of authors to suggest that a scale-independent physics regulates the self-organization. The simplest scale-independent fault model is a homogeneous one, in which it is assumed that the structure of a given seismic region does not play an important part in the earthquake process at any scale, and that the seismicity is dominated by the mechanics of the self-organization due to the stress redistribution on a homogeneous landscape of structure (4-11). The total stress on the fault in scalar one-and two-dimensional (antiplane) elastostatic fractures does not decrease with time; thus, ultimately the stress in the system must exceed the fracture strength everywhere and a fracture must occur that is larger than any given size. If we require that the largest earthquakes be of finite length and that their growth stops by the same mechanism as the smaller ones, then these quasistatic models must ultimately develop a fracture whose length is greater than the largest that is geophysically possible.Thus, an event must develop that is equal to the size of any finite computational lattice; at this point, the fracture interacts with the boundaries, and the system is no longer homogeneous: the evolutionary development of fractures is subsequently influenced by the interaction with the boundaries, by the Abbreviations: G-R, Gutenberg-Richter; SAF, San Andreas fault; B-K, Burridge-Knopoff; 1D, one-dimensional.Proc. Natl. Acad. Sci. USA 93 (1996) There are several scales of nonplanar geometric features: (i) topographic irregularities on fault surfaces and the presence of fault gouge that are the cause of friction on faults; (ii) larger-scale geometrical fluctuations of faults such as bends, stepovers, bifurcations, etc., of faults; (iii) the dendritic nature of the network of secondary faults associated with plate boundaries and a possibly tegular character of the space between elements of the network; and (iv) the influence of the curvature of the earth and the distances between triple junctions of the tectonic plates. We have no contribution to make to this discussion on the influence of this last item on the organization of earthquakes worldwide.While there is much to be said concerning the dynamics of friction at the smallest scales, limitations of space do not allow for a full discussion of ...