[1] We develop a forecast model to reproduce the distribution of main shocks, aftershocks and surrounding seismicity observed during 1986-2003 in a 300 Â 310 km area centered on the 1992 M = 7.3 Landers earthquake. To parse the catalog into frames with equal numbers of aftershocks, we animate seismicity in log time increments that lengthen after each main shock; this reveals aftershock zone migration, expansion, and densification. We implement a rate/state algorithm that incorporates the static stress transferred by each M ! 6 shock and then evolves. Coulomb stress changes amplify the background seismicity, so small stress changes produce large changes in seismicity rate in areas of high background seismicity. Similarly, seismicity rate declines in the stress shadows are evident only in areas with previously high seismicity rates. Thus a key constituent of the model is the background seismicity rate, which we smooth from 1981 to 1986 seismicity. The mean correlation coefficient between observed and predicted M ! 1.4 shocks (the minimum magnitude of completeness) is 0.52 for 1986-2003 and 0.63 for 1992-2003; a control standard aftershock model yields 0.54 and 0.52 for the same periods. Four M ! 6.0 shocks struck during the test period; three are located at sites where the expected seismicity rate falls above the 92 percentile, and one is located above the 75 percentile. The model thus reproduces much, but certainly not all, of the observed spatial and temporal seismicity, from which we infer that the decaying effect of stress transferred by successive main shocks influences seismicity for decades. Finally, we offer a M ! 5 earthquake forecast for 2005-2015, assigning probabilities to 324 10 Â 10 km cells. Motivation and Goals[2] There is growing evidence that large earthquakes can inhibit or promote failure on nearby faults for decades to centuries, and that the transfer of stress plays a governing role in this interaction [Harris, 1998;Stein, 1999;King and Cocco, 2000;Freed, 2005]. Some progress has also been achieved in explaining the time dependence of observed seismicity by Coulomb stress transfer in an elastic medium coupled with the concepts of Dieterich's [1994] rate and state friction [Stein et al., 1997;Parsons et al., 1999;Toda and Stein, 2003]. Competing and complementary concepts of dynamic stress triggering [Kilb, 2003;Gomberg et al., 2003], pore fluid diffusion [Peltzer et al., 1996], and viscoelastic stress transfer [Pollitz et al., 2004] have also been explored to explain aftershock and postseismic geodetic deformation, but here we will restrict ourselves to static stress changes coupled to rate/state friction.[3] Here we seek to explore how successive earthquakes redistribute stress, causing large, sudden seismicity rate changes that decay and are subsequently perturbed by other large shocks. To make the process more visible and accessible, we simulate the expected seismicity and compare it with the observations in several animations. To test the ability of the model to capture seismicity obser...

Abstract. The Kobe earthquake struck at the edge of the densely populated Osaka-Kyoto corridor in southwest Japan. We investigate how the earthquake transferred stress to nearby faults, altering their proximity to failure and thus changing earthquake probabilities. We find that relative to the pre-Kobe seismicity, Kobe aftershocks were concentrated in regions of calculated Coulomb stress increase and less common in regions of stress decrease. We quantify this relationship by forming the spatial correlation between the seismicity rate change and the Coulomb stress change. The correlation is significant for stress changes greater than 0.2-1.0 bars (0.02-0.1 MPa), and the nonlinear dependence of seismicity rate change on stress change is compatible with a state-and rate-dependent formulation for earthquake occurrence. We extend this analysis to future mainshocks by resolving the stress changes on major faults within 100 km of Kobe and calculating the change in probability caused by these stress changes. Transient effects of the stress changes are incorporated by the state-dependent constitutive relation, which amplifies the permanent stress changes during the aftershock period. Earthquake probability framed in this manner is highly timedependent, much more so than is assumed in current practice. Because the probabilities depend on several poorly known parameters of the major faults, we estimate uncertainties of the probabilities by Monte Carlo simulation. This enables us to include uncertainties on the elapsed time since the last earthquake, the repeat time and its variability, and the period of aftershock decay. We estimate that a calculated 3-bar (0.3-MPa) stress increase on the eastern section of the Arima-Takatsuki Tectonic Line (ATTL) near Kyoto causes fivefold increase in the 30-year probability of a subsequent large earthquake near Kyoto; a 2-bar (0.2-MPa) stress decrease on the western section of the ATTL results in a reduction in probability by a factor of 140 to 2000. The probability of a M•,, = 6.9 earthquake within 50 km of Osaka during 1997-2007 is estimated to have risen from 5-6% before the Kobe earthquake to 7-11% afterward; during 1997-2027, it is estimated to have risen from 14-16% before Kobe to 16-22%.

We calculate the probability of strong shaking in Istanbul, an urban center of 10 million people, from the description of earthquakes on the North Anatolian fault system in the Marmara Sea during the past 500 years and test the resulting catalog against the frequency of damage in Istanbul during the preceding millennium. Departing from current practice, we include the time-dependent effect of stress transferred by the 1999 moment magnitude M = 7.4 Izmit earthquake to faults nearer to Istanbul. We find a 62 +/- 15% probability (one standard deviation) of strong shaking during the next 30 years and 32 +/- 12% during the next decade.

Magma intrusions and eruptions commonly produce abrupt changes in seismicity far from magma conduits that cannot be associated with the diffusion of pore fluids or heat. Such 'swarm' seismicity also migrates with time, and often exhibits a 'dog-bone'-shaped distribution. The largest earthquakes in swarms produce aftershocks that obey an Omori-type (exponential) temporal decay, but the duration of the aftershock sequences is drastically reduced, relative to normal earthquake activity. Here we use one of the most energetic swarms ever recorded to study the dependence of these properties on the stress imparted by a magma intrusion. A 1,000-fold increase in seismicity rate and a 1,000-fold decrease in aftershock duration occurred during the two-month-long dyke intrusion. We find that the seismicity rate is proportional to the calculated stressing rate, and that the duration of aftershock sequences is inversely proportional to the stressing rate. This behaviour is in accord with a laboratory-based rate/state constitutive law, suggesting an explanation for the occurrence of earthquake swarms. Any sustained increase in stressing rate--whether due to an intrusion, extrusion or creep event--should produce such seismological behaviour.

[1] Two M $ 6 well-recorded strike-slip earthquakes struck just 4 km and 48 days apart in Kagoshima prefecture, Japan, in 1997, providing an opportunity to study earthquake interaction. Aftershocks are abundant where the Coulomb stress is calculated to have been increased by the first event, and they abruptly stop where the stress is dropped by the second event. This ability of the main shocks to toggle seismicity on and off argues that static stress changes play a major role in exciting aftershocks, whereas the dynamic Coulomb stresses, which should only promote seismicity, appear to play a secondary role. If true, the net stress changes from a sequence of earthquakes might be expected to govern the subsequent seismicity distribution. However, adding the stress changes from the two Kagoshima events does not fully capture the ensuing seismicity, such as its rate change, temporal decay, or migration away from the ends of the ruptures. We therefore implement a stress transfer model that incorporates rate/state friction, in which seismicity is treated as a sequence of independent nucleation events that are dependent on the fault slip, slip rate, and elapsed time since the last event. The model reproduces the temporal response of seismicity to successive stress changes, including toggling, decay, and aftershock migration. Nevertheless, the match of observed to predicted seismicity is quite imperfect, due perhaps to inadequate knowledge of several model parameters. However, to demonstrate the potential of this approach, we build a probabilistic forecast of larger earthquakes on the expected rate of small aftershocks, taking advantage of the large statistical sample the small shocks afford. Not surprisingly, such probabilities are highly time-and location-dependent: During the first decade after the main shocks, the seismicity rate and the chance of successive large shocks are about an order of magnitude higher than the background rate and are concentrated exclusively in the stress triggering zones.INDEX TERMS: 7209 Seismology: Earthquake dynamics and mechanics; 7223 Seismology: Seismic hazard assessment and prediction; 7230 Seismology: Seismicity and seismotectonics; 7260 Seismology: Theory and modeling; KEYWORDS: Coulomb stress, stress triggering, seismicity Citation: Toda, S., and R. Stein, Toggling of seismicity by the 1997 Kagoshima earthquake couplet: A demonstration of timedependent stress transfer,