Originally published as:

Size: px
Start display at page:

Download "Originally published as:"

Transcription

1 Originally published as: Kuehn, N. M., Hainzl, S., Scherbaum, F. (2008): Non-Poissonian earthquake occurrence in coupled stress release models and its effect on seismic hazard. - Geophysical Journal International, 174, 2, pp DOI:

2 Geophys. J. Int. (2008) 174, doi: /j X x Non-Poissonian earthquake occurrence in coupled stress release models and its effect on seismic hazard N. M. Kuehn, 1 S. Hainzl 2 and F. Scherbaum 1 1 Institute of Geosciences, University of Potsdam, Potsdam, Germany. nico@geo.uni-potsdam.de 2 GeoForschungsZentrum Potsdam, Potsdam, Germany Accepted 2008 April 28. Received 2008 April 28; in original form 2007 September 7 SUMMARY Most seismic hazard estimations are based on the assumption of a Poisson process for earthquake occurrence, even though both observations and models indicate a departure of real seismic sequences from this simplistic assumption. Instrumental earthquake catalogues show earthquake clustering on regional scales while the elastic rebound theory predicts a periodic recurrence of characteristic earthquakes on longer timescales for individual events. Recent implementations of time-dependent hazard calculations in California and Japan are based on quasi-periodic recurrences of fault ruptures according to renewal models such as the Brownian Passage Time model. However, these renewal models neglect earthquake interactions and the dependence on the stressing history which might destroy any regularity of earthquake recurrences in reality. To explore this, we investigate the (coupled) stress release model, a stochastic version of the elastic rebound hypothesis. In particular, we are interested in the time-variability of the occurrence of large earthquakes and its sensitivity to the occurrence of Gutenberg Richter type earthquake activity and fault interactions. Our results show that in general large earthquakes occur quasi-periodically in the model: the occurrence probability of large earthquakes is strongly decreased shortly after a strong event and becomes constant on longer timescales. Although possible stress-interaction between adjacent fault zones does not affect the recurrence time distributions in each zone significantly, it leads to a temporal clustering of events on larger regional scales. The non-random characteristics, especially the quasi-periodic behaviour of large earthquakes, are even more pronounced if stress changes due to small earthquakes are less important. The recurrence-time distribution for the largest events is characterized by a coefficient of variation from 0.6 to 0.84 depending on the relative importance of small earthquakes. Key words: Seismic cycle; Earthquake interaction, forecasting, and prediction; Statistical seismology. GJI Seismology 1 INTRODUCTION A key aspect in probabilistic seismic hazard analysis is the accurate modelling of the spatiotemporal distribution of earthquakes and in particular their recurrence times. These calculations are usually done with the Poisson process for distributed seismicity, making the assumption that seismicity in a specific time interval is independent of previous seismicity. Advantages of the Poissonian model are its simplicity and the fact that it arises naturally in a superposition of event sequences from independent sources (Vere-Jones 1995). Because the Poisson model has an exponential recurrence time distribution and thus a constant hazard function, this assumption leads to time-independent seismic hazard estimates. However, observations show that earthquake occurrence usually does not follow a Poisson process. While short- and intermediateterm earthquake clustering such as aftershock activity is evident (Kagan & Jackson 1991a; Rhoades & Evison 2004; Lombardi & Marzocchi 2007), the recurrence of large earthquakes on individual faults is rather speculative because seismic catalogues are short compared to the underlying processes (Parsons 2005). Even in regions with excellent seismic monitoring, like the Parkfield segment in California, the number of large earthquakes is below 10. Additional information may be gained from palaeoseismology, but a reliable statistical evaluation based on observational data is not feasible at the moment (see e.g. the dispute about the seismic gap theory Kagan & Jackson 1991b, 1995; Wesnousky 1994, 1996; Kagan 1996; Wyss & Wiemer 1999). Despite direct observational evidence, quasi-periodic recurrence of characteristic events is expected from the classical elastic rebound theory (Reid 1910, 1911) and has been observed for the proxy-data of small repeating earthquakes (Nadeau & McEvilly 1997) as well as in numerical simulations of fault dynamics (e.g. Zöller et al. 2006). Therefore, recent implementations of time-dependent hazard assessment in California and Japan (WGCEP 2003; Stein et al. 2006) C 2008 The Authors 649

3 650 N. M. Kuehn, S. Hainzl and F. Scherbaum assume quasi-periodic recurrence of characteristic events modelled by renewal processes such as the Brownian Passage Time model (Matthews et al. 2002; Zöller & Hainzl 2007). Although those renewal models account for some intrinsic variability, they neglect the dependence on the stressing history and the impact of earthquake interactions. However, stress transfer, especially from nearby faults and frequent smaller, Gutenberg Richter type earthquakes, can advance or delay the earthquake nucleation of large ruptures and it is an open question whether this will destroy any natural periodicity. In this work, we therefore, investigate numerical simulations of the stress release model which is a stochastic version of elastic rebound theory where the earthquake probability is a function of the stress level which is itself altered by each earthquake (Vere-Jones 1978). The stress release model is particularly interesting, since it is easily extendable into the spatial domain (the so-called coupled or linked stress release model) and can take stress interactions due to fault communication (Liu et al. 1999; Lu et al. 1999; Bebbington & Harte 2003). Several studies (for a review see Bebbington & Harte 2003) have applied the stress release model to reproduce the recurrence of large historical earthquakes in various seismically active regions showing that it is mostly superior to the Poissonian model (Zheng & Vere-Jones 1991, 1994; Zhuang & Ma 1998; Liu et al. 1999; Lu et al. 1999; Lu & Vere-Jones 2000; Lu 2005; Rotondi & Varini 2006, 2007; Imoto & Hurukawa 2006). Although this does not prove that the stress release model reflects the true evolution of seismicity on long timescales, it indicates that the stress release model can be regarded as a valid working hypothesis for studying the effect of stress interactions on recurrence time distributions. We investigate numerical simulations of the stress release model that consists of a large number of seismic cycles. Thereby, we explicitly take into account Gutenberg Richter distributed small and intermediate magnitude events because these might play an important role for stress redistribution and should not be neglected (Hanks 1992; Helmstetter et al. 2005; Marsan 2005). Additionally, we address the important question how coupling between different seismic zones affects the occurrence probability distributions on a local and regional scale. In the first sections, we will outline the model (Section 2) and the numerical setup (Section 3) in more detail. Subsequently, the resulting characteristics of the model are discussed (Section 4) and the main conclusions are summarized (Section 5). 2 THE MODEL The stress release model (Vere-Jones 1978) is a stochastic version of elastic rebound theory (Reid 1911), one of the classical models for earthquake occurrence. It assumes that the occurrence probability of an earthquake in a seismic active region depends on some quantity Z which may be interpreted as the mean stress level present in the region. The occurrence probability is assumed to increase with an increasing stress level and to drop when an earthquake occurs. Due to tectonic loading, the stress level is generally assumed to increase linearly between two earthquakes and to drop when an earthquake occurs, so its temporal evolution can be written as Z(t) = Z 0 + ρt t i <t Z i, (1) where Z 0 is the stress level at time t = 0 and ρ is the loading rate. t i <t Z i is the accumulated stress release of earthquakes that occurred prior to time t, with t i and Z i the occurrence time and stress release of the ith earthquake, respectively. It is important to note that the quantity Z, which is usually simply called stress, is a scalar and not the true stress tensor. Z is used to describe the underlying state of the system or region under study. It is assumed that the average stress release Z of an earthquake is related to some power of its energy release, Z E η. Combining the energy release with the magnitude using the relation m = 2 3 log 10 E + const (Kanamori & Anderson 1975), the following equation which connects the stress release with the magnitude is derived: Z = 10 h(m m 0), (2) where m 0 denotes the lower cut-off magnitude. For h = 0.75 and 1.5, the mean stress release is proportional to Benioff strain (Benioff 1951) and to seismic moment, respectively. Most applications use h = The occurrence probability of an earthquake is controlled by a stress-dependent conditional intensity function λ(z). For small t, λ(z) t represents the probability of an event occurring in the time interval (t, t + t). Obviously, λ has to be a non-decreasing function of Z to agree with the physical understanding that higher stress levels result in a higher earthquake risk. λ(z) = const would result in a random (Poisson process) model, whereas { 0, Z Zc λ(z) =, Z > Z c corresponds to a time-predictable model (Jaume & Bebbington 2004). We follow most of the articles in the literature and use λ(z) = exp [α + β Z]. (3) Our trials with a power law and a linear function for λ(z) indicate that the behaviour of the model does not depend strongly on its functional form. This is concordant with conclusions of Zheng & Vere-Jones (1994). Furthermore, it can be shown that with the exponential function, the stress release model is equivalent to the rate-state friction model of Dieterich (1994) as well as stress corrosion (Scholz 2002) in its instantaneous response to stress changes. Note that in eq. (3) the Poisson process (β = 0) is included. With eq. (1), the conditional intensity function becomes [ ( λ(t) = exp α + β Z 0 + ρt ) ] Z i. (4) t i <t Stress transfer and stress-triggering between distant faults and seismogenic zones can be included in a multivariate extension of the simple stress release model, leading to the so-called coupled (or linked ) stress release model (Liu et al. 1999; Lu et al. 1999; Bebbington & Harte 2003). Here, the local stress depends also on the earthquake activity in adjacent regions Z (i) (t) = Z (i) (0) + ρ (i) t j Z ( j) k θ ij S ( j) (t), (5) where S ( j) (t) = t k <t is the accumulated stress release of earthquakes in region j that occurred prior to time t. Superscript (i) denotes the ith region and θ ij measures the fraction of stress transferred from region j to region i. A positive or negative value of θ ij results in damping (stress shadowing) or excitation (stress triggering), respectively. The diagonal elements are θ ii = 1. The occurrence probability of the next event in each region is again modelled by a stress-dependent hazard function λ (i) (t) = exp [ α (i) + β (i) Z (i) (t) ]. (6)

4 Non-Poissonian earthquake occurrence 651 Note that in the simple model as described here, aftershocks that occur in the same region as the main shock are not taken into account. However, one might think of more elaborate models where one region is decomposed into a fault area describing the recurrence of main shocks which transfer stress to an aftershock zone (see, e.g. Borovkov & Bebbington 2003). Simulation of the stress release model requires specification of a probability distribution of the event sizes (i.e. magnitudes) to calculate the stress release for the events according to eq. (2). For that, we assume the well-known Gutenberg Richter distribution (Gutenberg & Richter 1944), empirically observed for regional seismicity. In particular, we use the double truncated Gutenberg Richter distribution where the maximum magnitude depends on the stress level. The probability density has the form ln 10 b 10 bm p(m) =, (7) 10 bm 0 10 bm max where m 0 and m max are the minimum and maximum magnitude, respectively. m max depends on the stress level in a way that the maximum possible earthquake just unloads the system. In the Appendix, it is shown that m max can be easily calculated by inverting eq. (2), m max = m 0 + log 10 Z. (8) h 3 EXPERIMENTAL DESIGN We investigate the (coupled) stress release model regarding its temporal properties and its consequences for seismic hazard. Therefore, we perform various simulations with different parameter sets. Parameters that can influence the temporal behaviour of the model are α, β and ρ in eq. (4) and the coupling strength θ in eq. (5). Furthermore, the parameters of the frequency magnitude distribution have to be specified. Since in seismic hazard analysis the frequency magnitude relation is often estimated independently, we assume that the a, b and m max value of the Gutenberg Richter relation are given. As a consequence, we can then fix the parameters α and ρ in eq. (4), as is shown in Appendix A. The temporal variability of the model is then only dependent on two free parameters, namely β and h. We perform simulations of both the uncoupled and the coupled model. Application of the coupled stress release model usually needs careful consideration concerning the number, shape and size of the regions taken into account (Bebbington & Harte 2003). However, to reduce complexity, we investigate the simple case of isotropically coupled regions. Stress is only transferred to the four neighbouring regions, and the coupling strength θ between different regions is always the same in a single simulation, while θ ii = 1. We require that the total stress that is transferred to neighbouring regions must be smaller than the stress released in the host region, which leads to the constraint θ To be as simple as possible, all regions of the model are assigned the same parameters. For the Gutenberg Richter relation, we use a maximum magnitude m max = 7.5 and its recurrence interval T r = a in all simulations. However, to study the parameter dependence of our results, we repeat the simulations with different values of θ, β and h. If values other than h = 0.75 are used, β needs to be adjusted in order to have comparable impacts of large earthquakes on the conditional intensity function. Therefore, we recalculate β by requiring that an event of magnitude m max has the same effect on the occurrence probability as before, namely β10 h(mmax m0) = β (mmax m0). (9) After neglecting the first events in each simulation to ensure that the system has evolved towards the equilibrium state, two catalogues are obtained. One contains events of any magnitude. Since we are particularly interested in the behaviour of the large magnitudes, the other catalogue consists of 1000 events per region with a magnitude m > 6.5. The simulations themselves are carried out as follows. (1) Using the conditional intensity function λ(t) (eq. 4), the next event time t i in each region (index i) is simulated using the transform method (Devroye 1986), namely by inverting u i = exp [ ti τ i e α+β Z i (s) ds ], (10) where u i is uniformly distributed within [0, 1] and τ i is the time of the last stress jump in this region (for the derivation, see Appendix B). (2) In the region with the smallest value, the next event occurs at time t i. (3) The magnitude of this event is simulated from distribution (7) and converted into a stress release Z by eq. (2). (4) The stress transfer is calculated by θ Z. In all affected regions, the stress level is updated and the subsequent event times are recomputed with eq. (10) according to the altered intensity functions. Then step 2 is again repeated and so on. 4 SIMULATION RESULTS We analyse the characteristics of large earthquakes occurring in this model. At first we investigate the uncoupled model describing an isolated seismic zone. Subsequently, we turn to the effects of interactions between different seismic regions. Our standard simulations are carried out with the parameters β = 0.1, h = 0.75 and b = 1in the Gutenberg Richter relation. 4.1 The uncoupled model In the case of the uncoupled model, a key result of our simulations is that the occurrence of small and intermediate earthquakes does not lead to a random recurrence of the largest events. Even if earthquakes occur with a fractal size distribution, related to a Gutenberg Richter distribution, the evolution of the stress level shows some regularity. A weak seismic cycle occurs which is illustrated for one example in Fig. 1. A systematic analysis of the interevent-time distribution is shown in Fig. 2 for a much longer simulation. While this distribution is monotonically decaying (with a power-law tail) if all events are considered, it is peaked if only larger magnitude earthquakes (m > 6.5) are taken into account. The peak corresponds to a preferred interevent-time, indicating quasi-periodicity. We have also fitted the exponential, gamma, log-normal, Weibull and Brownian passage time distributions to the interevent-time distribution, since these are widely used in seismic hazard calculations (Matthews et al. 2002; Yakovlev et al. 2006), using the maximum likelihood method. The highest likelihood value was achieved by the gamma distribution, with log-likelihood values of LL = 7961 for the gamma, LL = 7966 for the Weibull, LL = 7973 for the log-normal, LL = 7979 for the exponential and LL = 8008 for the BPT distribution. The right part of Fig. 2 shows the interevent time histogram of the larger magnitude events together with the best-fitting gamma distribution. A general characteristic of the interevent-times can be quantified by the coefficient of variation, which is defined by the ratio of

5 652 N. M. Kuehn, S. Hainzl and F. Scherbaum magnitude time [a] Figure 1. Example of evolving state of a 50-yr-long simulation with magnitudes (bars) and conditional intensity (dashed line) density density interevent time [a] interevent time [a] Figure 2. Left-hand side: Interevent-time distribution for all events with m > 3 (dashed) and fitted exponential distribution (solid); right-hand side: Intereventtime histogram for the large events (m > 6.5) and fitted gamma distribution. The parameters of the simulation are β = 0.1, h = 0.75, b = 1. standard deviation to mean, C V = SD/mean. For a Poisson (completely random) process the coefficient of variation has a value C V = 1, whereas values C V > 1 and C V < 1 correspond to a clustered and quasi-periodic process, respectively (Telesca et al. 2002). Fig. 3 illustrates how the gamma distribution the best fitting distribution for the interevent times of the larger magnitudes changes for the different C V -values of 0.5, 1 and 2, if the mean value is constant. For C V = 1, the gamma distribution corresponds to the exponential distribution. In Fig. 4, the resulting coefficient of variation is shown as a function of the lower magnitude cut-off for different values of h. Here we see that for small cut-off magnitudes the coefficient of variation takes values C V > 1, indicating that small events occur clustered. For larger cut-off magnitudes, the values of the coefficient of variation decrease, reaching a minimum at C V = 0.85 at m = 5.1 for h = 0.75 and C V = 0.61 at m = 7 for h = 1.5. Hence, earthquakes with a magnitude greater than these are characterized by quasi-periodic occurrence. The quasi-periodic behaviour stems from the fact that large earthquakes decrease the stress level by a significant amount, leaving it far from its equilibrium state, whereas small earthquakes cause only slight deviations in the stress level. However, the quasiperiodic behaviour applies only if all large events are considered together (i.e. m > 5.1 for h = 0.75 or m > 7 for h = 1.5), as they all produce a stress release high enough to cause significant density C V = 0.5 C V = 1 C V = Figure 3. Gamma distribution with parameter values corresponding to C V = 0.5, = 1 and = 2, where the mean value is fixed to 1. x

6 Non-Poissonian earthquake occurrence 653 CV cutoff magnitude CV cutoff magnitude Figure 4. Left-hand side: Coefficient of variation for a simulation with β = 0.1 and b = 1 as a function of the lower magnitude cut-off for h = 0.75 (stars), h = 1 (crosses) and h = 1.5 (rhombs); right-hand side: Coefficient of variation for seismicity of a single fault, taken from Zöller et al. (2006). stress release per magnitude h = 0.75 h = h = magnitude Figure 5. Fraction of stress released per magnitude for Benioff strain (h = 0.75), seismic moment (h = 1.5) and a value h = 1. perturbations in the stress level. The largest earthquakes, if considered separately, show an increasing trend of the coefficient of variation (right-hand part of Fig. 4) and hence a more random occurrence since slightly smaller events frequently occurring in between also significantly perturb the stress level. The result of a quasi-periodic behaviour of large events is in general agreement with the results for the fault model by Zöller et al. (2006) which simulates the seismicity on a single fault plane by spontaneous quasi-dynamic rupture propagation and stopping on a discretized fault embedded in the elastic crust. For this physically based fault model, the corresponding dependence of the C V -value on the lower magnitude cut-off is shown in Fig. 4. The general trend the decrease of C V for larger magnitudes cut-offs is very similar to the stochastic stress release model, although some differences occur for very large magnitude cut-offs: This discrepancy is probably related to the fact that the Gutenberg Richter frequency magnitude distribution breaks down in the physical model at a magnitude of approximately 6 due to the occurrence of so-called runaway ruptures. These are events which will rupture, after growing to a certain size, the whole fault because the stress concentration at the rupture tip becomes large enough to overcome all stress barriers. Thus independent magnitude 6 and 7 earthquakes do not occur in this physical model in contrast to our simulations of the stress release model. Despite this specific difference, the general similarity of the clustering properties between the stochastic stress release model and the computational expensive quasi-dynamic fault simulations indicates that the stress release model is able to mimic the fundamental fault dynamics in a first approximation. While it is not yet computationally feasible to study the effects of fault interactions in these more realistic models, this can be done with the linked stress release model (see Section 4.2). For the stress release model, the dependence of the recurrence time statistics on the h-value can be seen in Fig. 4. The dependence is caused by differences in the amount of stress that is released in different magnitude ranges. This is shown in Fig. 5. If Benioff strain is used (h = 0.75), the greater portion of stress is released by the smaller events, whereas for seismic moment (h = 1.5) the main stress release results from the large magnitude events. Hence, in the latter case, the magnitude level that can cause a large deviation from the normal stress level is enhanced. Together with the survivor function (Daley & Vere-Jones 2003), the interevent-time distribution can be used to calculate an empirical hazard function, which is defined as the event rate at time t given survival up to t. The hazard function is calculated as the probability density function divided by the survivor function. The Poisson process has a constant hazard function (Daley & Vere- Jones 2003). Fig. 6 shows the empirical hazard function (event rate) for all (left-hand side) and the large events (right-hand side).

7 654 N. M. Kuehn, S. Hainzl and F. Scherbaum hazard rate 3 2 hazard rate time [a] time [a] Figure 6. Empirical hazard function calculated from the interevent-time histogram with β = 0.1, h = 0.75 and b = 1. Left-hand side: All events with m > 3; right-hand side: All events with m > CV magnitude cutoff Figure 7. Coefficient of variation for a single region with (θ = 0.23, crosses) and without (θ = 0, stars) coupling (h = 0.75 and β = 0.1). The hazard function for all events decreases from high values at small interevent-times to an almost constant level. This corresponds to the clustering observed for the small magnitudes. By contrast, the hazard level for the large events (m > 6.5) is almost zero immediately after the last large earthquake and then starts to rise for increasing waiting times. This is due to the fact that the system needs to be reloaded before another large earthquake can happen again. On long timescales, the hazard function reaches a constant level. Changing the values of β does not alter the general characteristics of the model: events with a small magnitude occur clustered (C V > 1), whereas large ones occur quasi-periodically (C V < 1). However, for β 0, eq. (3) reduces to λ(t) = e α = const, corresponding to a Poisson process (C V = 1). Hence, the C V -values for low cutoff magnitudes decrease with decreasing β, but increase for higher values of β. However, the minimum of the coefficient of variation is quite stable: For 30 simulations (h = 0.75) with β-values varying between 0.01 and 0.3, the mean minimum takes place at C V = ± This indicates that the observed seismic cycle is mainly due to the seismic gap after large earthquakes which results from the time needed to reload the seismic zone such that it can produce the next large event. This reloading time depends on ρ but not on β. 4.2 The effect of coupling In the previous section, we examined the simple (uncoupled) stress release model concerning its temporal properties. Now we turn to the effects of stress interaction between different seismic regions. We find that the presence of coupling does not significantly affect the general temporal variability observed in each seismic zone, as can be seen in Fig. 7. Here, the coefficient of variation is shown for one specific region of the uncoupled (θ = 0, stars) and a strongly coupled (θ = 0.23, crosses) model, respectively. Parameters h = 0.75 and β = 0.1 are used in both cases. Since stress is transferred into the four neighbouring regions, a coupling strength θ = 0.23 means that 92 per cent of the stress released by an earthquake is redistributed in the system and 8 per cent is lost. As one can see in Fig. 7, the values of the coefficient of variation are very similar for both the uncoupled and coupled case. In the presence of coupling, the C V -values are just slightly increased, and the minimum of the curve is shifted a little to higher magnitudes. However, the general appearance of the graph is unchanged. Hence, equivalently to the uncoupled case, large earthquakes occur quasi-periodically in each individual seismic region. On the other hand, earthquakes in different regions start to correlate in the presence of stress interactions. This effect can be seen in

8 Non-Poissonian earthquake occurrence CV magnitude cutoff CV magnitude cutoff Figure 8. Dependence of the coefficient of variation on the coupling strength for the stacked activity, for θ = 0 (crosses), θ = 0.1 (stars) and θ = 0.23 (rhombs). Left-hand side: Benioff strain (h = 0.75); right-hand side: Seismic moment (h = 1.5) hazard rate normalized time Figure 9. Hazard function for large events (m > 6.5) in the uncoupled (solid) and strongly coupled (dashed) model. The long dashed line represents the hazard in the same region, the short dashed line in regions adjacent to that where the earthquake occurred. The times are normalized by the mean recurrence time. Fig. 8. Here, all events from all 100 regions are considered together as one catalogue for which the coefficient of variation is shown. This corresponds to a superposition of the individual stress release models for the different regions. The left-hand side of Fig. 8 shows simulations with h = 0.75, the right-hand part with h = 1.5. In the case of an uncoupled model (θ = 0), the different regions are independent, so their superposition leads to a Poisson process (Vere-Jones 1995) with C V = 1 (crosses in Fig. 8). However, if coupling is present, the processes in the different regions are not independent any more, so deviations from C V = 1 emerge. Coupling leads to a clustering of the events, independent of the magnitude, since each earthquake leads to an increase of the occurrence probability in adjacent zones. The clustering is also observed for small magnitudes because of triggered aftershock activity in neighbouring zones. Hence, on the global scale seismic cycling cannot be observed. One can see in Fig. 8 that the strongest clustering occurs for intermediate to large cut-off magnitudes. We explain the tendency that the largest events are again more random in the same way as before for the isolated region: The intermediate to large earthquakes can already perturb the stress field significantly and thus randomize the occurrence times of the largest events. Similar to the case of a single region, the clustering is enhanced for h = 1.5 when the stress is mainly released by the largest earthquakes (see Fig. 5). Fig. 9 compares the hazard functions in coupled and uncoupled seismic regions. On the one hand, the probability for the recurrence of large earthquakes in the same seismic zone is found to be almost identical in both cases. This indicates again that coupling does not significantly alter the seismic cycles in each region. On the other hand, the hazard level of a large earthquake is greatly enhanced in adjacent regions immediately after the event, if both regions are positively coupled. Thus correlated occurrence of major events in neighbouring regions will occur. However, the hazard function decays in adjacent regions quickly with time and, on long timescales, it reaches in all cases a constant level after approximately 20 per cent of the mean recurrence time of the large events. 5 SUMMARY AND CONCLUSIONS In probabilistic seismic hazard assessment (Cornell 1968), most studies assume a Poissonian occurrence of earthquakes, that is, a random occurrence in time and space. However, earthquake observations show that earthquakes cannot be treated as independent random events, because the occurrence of earthquakes is characterized by long- and short-term patterns (Scholz 2002). In seismically highly active regions where active faults have been identified and their characteristics are well known, often a quasi-periodic

9 656 N. M. Kuehn, S. Hainzl and F. Scherbaum recurrence of characteristic events combined with Poisson distributed background activity is deemed applicable (Matthews et al. 2002; WGCEP 2003). However, the validation of this assumption is difficult and it is thus controversially debated (see e.g. Kagan 1996). In particular, it is questionable whether or not the self-organized stress evolution due to fault communication destroy any regularity of the recurrence of large earthquakes. In particular, also small earthquakes seem to play an important role for the stress redistribution and should not be neglected (Hainzl et al. 2000; Helmstetter et al. 2005; Marsan 2005). In general, observational data could be used to directly estimate the hazard function, e.g. by non-parametric methods (Faenza et al. 2003). However, this usually cannot solve the problem because catalogues are mostly short compared to the underlying processes (Parsons 2005). In this work, we investigate the impact of stress evolution and interactions for seismic hazard estimations by means of the stress release model. In this model the earthquake activity depends on a scalar stress-variable denoting the mean stress level in the seismic region under consideration (Vere-Jones 1978). Tectonic loading increases this level steadily whereas it is decreased by earthquakes (dependent on the earthquake magnitude) based on the assumptions of classical models of earthquake mechanism (Reid 1911). The stress release model can be considered as a mean field model describing the net effect of earthquakes in a seismogenic volume. In previous studies, it has been shown that the stress release model achieves a fit to main shock occurrence that is usually superior to the Poisson process (e.g. Zheng & Vere-Jones 1994; Bebbington & Harte 2001; Imoto 2001; Imoto & Hurukawa 2006; Rotondi & Varini 2006). Hence, it seems worthwhile thinking about its application to seismic hazard analysis, although several issues pertaining the fitting of the model (e.g. the spatial zonation) have to be carefully considered (Bebbington & Harte 2003) and the model has not yet been rigorously tested (as opposed to other models, see e.g. Schorlemmer et al. (2007) and other articles in the same volume). Nevertheless, given its physical basis in elastic rebound theory and the number of good fits achieved, the stress release model can be regarded as a viable hypothesis for earthquake occurrence in seismically active regions. The aim of this work is to investigate the predictions of the stress release model in the general case of Gutenberg Richter type seismicity and fault interactions. Therefore, we perform numerical simulations of this model based on the well-known Gutenberg Richter relation for the frequency magnitude distribution. In practice, the parameters of the magnitude distribution can be determined independently, and thus most of the model parameters of the stress release model can be fixed (see Appendix). The remaining two parameters are describing the strength of the time dependence (β) and the scaling between stress release and magnitude (h). Both parameters are found to influence the strength but not the general characteristics of the observed non-poissonian behaviour. Our simulations show that the occurrence of earthquakes is generally correlated, in particular small earthquakes tend to cluster, while large ones occur quasiperiodically. Large earthquakes decrease the occurrence probability of the next event significantly within the first approximately 20 per cent of the mean recurrence time (see Fig. 9). On longer timescales, the hazard function reaches a constant level. The corresponding peaked recurrence-time distribution for the largest events is characterized by coefficients of variation from 0.61 (for h = 1.5) to 0.84 (for h = 0.75). This is similar to results from simulations of a physically based model of a single fault (Zöller et al. 2006). While the effect of fault interactions cannot be analysed appropriately by latter models due to computational limits, the linked stress release model can serve as a first approximation. We find that the above characteristics are robust if stress coupling between adjacent seismic zones is considered. Even very strong coupling does not significantly change the recurrence-time distributions of each zone, and major earthquakes occur quasi-periodically in each region. This is in agreement with conclusions drawn from the analysis of a number of real, but very short sequences of recurrent sequences on specific faults (Ellsworth et al. 1999). However, this does not imply that seismicity in adjacent seismic zones are independent in our simulations. Large events in different zones tend to cluster which is reflected in a strongly increased hazard level after a major failure in neighbouring regions (cf. Fig. 9). The role of small earthquakes can be studied in this model by varying the parameter h for a fixed Richter b-value. If Benioff strain is used (h = 0.75), a greater portion of stress is released by small magnitude events than by large earthquakes, whereas it is vice versa if seismic moment scaling (h = 1.5) is assumed (see Fig. 5). The non-poissonian behaviour mentioned above is found to be strongly enhanced for larger h-values when small earthquakes play only a minor role. Thus small earthquakes randomize the stress state evolution. However, even if they are dominating the stress release (in the case of h = 0.75), the seismicity shows a clear time-dependence. These results might have an important impact on seismic hazard assessment. Although direct application of the stress release model to specific seismic source zones might be problematic because of validation and fitting problems (cf. Bebbington & Harte 2003), our analysis shows that the occurrence of seismic cycles seems to be a robust feature which occur even in the presence of Gutenberg Richter distributed earthquake sizes and strong fault interactions. Thus, the numerically observed hazard functions for distributed seismicity could be a basis for future considerations of seismic cycles even in source regions where the active fault system is poorly known. ACKNOWLEDGMENTS This work was supported by the Deutsche Forschungsgemeinschaft (SCHE280/14.2). We are grateful to Warner Marzocchi and Jochen Woessner as well as to the editor Massimo Cocco for suggestions which helped to improve the manuscript. REFERENCES Beauval, C., Hainzl, S. & Scherbaum, F., Probabilistic seismic hazard estimation in low-seismicity regions considering non Poissonian seismic occurrence, Geophys. J. Int., 164, Bebbington, M.S. & Harte, D., On the statistics of the linked stress release model, J. Appli. Prob., 38A, Bebbington, M.S. & Harte, D., The linked stress release model for spatio-temporal seismicity: formulations, procedures and applications, Geophys. J. Int., 54, Benioff, H., Crustal strain characteristics derived from earthquake sequences, Trans. AGU, 32, Borovkov, K. & Bebbington, M.S., A stochastic two-node stress transfer model reproducing Omori s law, Pure appl. Geophys., 160, Cornell, C.A., Engineering seismic risk analysis, Bull. seism. Soc. Am., 58, Daley, D.J. & Vere-Jones, D., 2003, An Introduction to the Theory of Point Processes, 2nd edn, Vol. 1, Springer-Verlag, New York. Devroye, L., Non-Uniform Random Variate Generation, Springer- Verlag, New York. Dieterich, J., A constitutive law for rate of earthquake production and application to earthquake clustering, J. geophys. Res., 99, Ellsworth, W.L., Matthews, M.V., Nadeau, R.M., Nishenko, S.P., Reasenberg, P.A. & Simpson, R.W., 1999, A physically based earthquake recurrence

10 Non-Poissonian earthquake occurrence 657 model for estimation of long-term earthquake probabilities, U.S. Geol. Surv. Open-File Rept., Faenza, L., Marzocchi, W. & Boschi, E., A non-parametric hazard model to characterize the spatio-temporal occurrence of large earthquakes: an application to the Italian catalogue, Geophys. J. Int., 155, Gutenberg, B. & Richter, C.F., Frequency of earthquakes in California, Bull. seism. Soc. Am., 34, Hainzl, S., Zöller, G., Kurths, J. & Zschau, J., Seismic quiescence as an indicator for large earthquakes in a system of self-organized criticality, Geophys. Res. Lett., 27, Hanks, T.C., Small earthquakes, tectonic forces, Science, 256, Helmstetter A., Kagan, Y.Y. & Jackson D.D., Importance of small earthquakes for stress transfers and earthquake triggering, J. geophys. Res., 110, B05S08. Imoto, M., Application of the stress release model to the Nankai earthquake sequence, southwest Japan, Tectonophysics, 338, Imoto, M. & Hurukawa, N., Assessing potential seismic activity in Vrancea, Romania, using a stress release model, Earth Planets Space, 58, Imoto, M., Maeda, K. & Yoshida, A., Use of statistical models to analyze periodic seismicity observed for clusters in the Kanto region, Central Japan, Pure appl. Geophys., 155, Jaume, S.C. & Bebbington, M.S., Accelerating seismic release from a self-correcting stochastic model, J. geophys. Res., 109, Kagan, Y.Y., Comment on The Gutenberg-Richter or characteristic earthquake distribution, which is it? by Steven G. Wesnousky, Bull. seism. Soc. Am., 86, Kagan, Y.Y. & Jackson, D.D., 1991a. Long-term earthquake clustering, Geophys. J. Int., 104, Kagan, Y.Y. & Jackson, D.D., 1991b. Seismic gap hypothesis, ten years after, J. geophys. Res., 96, Kagan, Y.Y. & Jackson D.D., New seismic gap hypothesis, five years after, J. geophys. Res., 100, Kanamori, H. & Anderson, D.L., Theoretical basis of some empirical relations in seismology, Bull. seism. Soc. Am., 65, Lombardi, A.M. & Marzocchi, W., Evidence of clustering and nonstationarity in the time distribution of large worldwide earthquakes, J. geophys. Res., 112, B Liu, J., Chen, Y., Shi, Y. & Vere-Jones, D., Coupled stress release model for time-dependent seismicity, Pure appl. Geophys., 155, Lu, C., The degree of predictability of earthquakes in several regions of China: statistical analysis of historical data, J. Asian Earth Sci., 25, Lu, C. & Vere-Jones, D., Application of linked stress release model to historical earthquake data: comparison between two kinds of tectonic seismicity, Pure appl. Geophys., 157, Lu, C., Harte, D. & Bebbington, M.C., A linked stress release model for historical Japanese earthquakes: coupling among major seismic regions, Earth Planets Space, 51, Marsan, D., The role of small earthquakes in redistributing crustal elastic stress, Geophys. J. Int., 163, McCann, W.R., Nishenko, S.P., Sykes L.R. & Krause, J., Seismic gaps and plate tectonics: seismic potential for major boundaries, Pure appl. Geophys., 217, Matthews, M.V., Ellsworth, W.L. & Reasenberg, P.A., A Brownian model for recurrent earthquakes, Bull. seism. Soc. Am., 92, Nadeau, R.M. & McEvilly, T.V., Seismological studies at Parkfield V: characteristic microearthquake sequences as fault-zone drilling targets, Bull. seism. Soc. Am., 87, Ogata, Y., Statistical models of point occurrences and residual analysis for point processes, J. Am. Stat, Assoc., 83, Parsons, T., Significance of stress-transfer in time-dependent earthquake probability calculations, J. geophys. Res., 110, B05S02, doi: /2004jb Reid, H.F., The mechanics of the earthquake, The California Earthquake of April 18, 1906, Report of the State Investigation Commission, Vol. 2, Carnegie Institution of Washington, Washington, DC. Reid, H.F., The elastic rebound theory of earthquakes, Bull. Dep. Geol. Univ. Calif., 6, Rhoades, D.A. & Evison, F.F., Long-range earthquake forecasting with every earthquake a precursor according to scale, Pure appl. Geophys., 161, Rotondi, R. & Varini, E., Bayesian analysis of marked stress release models for time-dependent hazard assesment in the western Gulf of Corinth, Tectonophysics, 423, Rotondi, R. & Varini, E., Bayesian inference of stress release models applied to some Italian seismogenic zones, Geophys. J. Int., 169, Scholz, C.H., The Mechanics of Earthquakes and Faulting, Cambridge University Press, New York. Schorlemmer, M., Gerstenberger, C., Wiemer, S., Jackson, D.D. & Rhoades, D.A., Earthquake likelihood model testing, SRL 78, Stein, R.S., Toda, S., Parsons, T. & Grunewald, E., A new probabilistic seismic hazard assessment for greater Tokyo, Philos. Trans. Royal SOC. A Math. Phys. Eng. Sci., 364, Telesca, L., Cuomo, V., Lapenna, V. & Macchiato, M., On the methods to identify clustering properties in sequences of seismic time-occurrences, J. Seism., 6, Utsu, T., Ogata, Y. & Matsu ura, R.S., 1995, The centenary of the Omori formula for a decay law of aftershock activity, J. Phys. Earth, 43, Vere-Jones, D., Earthquake prediction a statistician s view, J. Phys. Earth, 26, Vere-Jones, D., Forecasting earthquakes and earthquake risk, Int. J. Forecast., 11, Vere-Jones, D., Robinson, R. & Wang, W., Remarks on the accelerated moment release model: problems of model formulation, simulation and estimation, Geophys. J. Int., 144, Yakovlev, G., Turcotte, D.L., Rundle, J.B. & Rundle, P.B., Simulation based distributions or earthquake recurrence times on the San Andreas fault system, Bull. seism. Soc. Am., 96, Wesnousky, S.G., The Gutenberg-Richter or characteristic earthquake distribution, which is it?, Bull. seism. Soc. Am., 84, Wesnousky, S.G., Reply to Yan Kagan s comment on The Gutenberg- Richter or characteristic earthquake distribution, which is it?, Bull. seism. Soc. Am., 86, Working Group on California Earthquake Probabilities, Earthquake probabilities in the San Francisco Bay region, US Geol. Survey Open File Report , US Geological Survey. Wyss, M. & Wiemer, S., How can one test the seismic gap hypothesis? The case of repeated ruptures in the Aleutians, Pure appl. Geophys., 157, 155, Zheng, X. & Vere-Jones, D., Applications of stress release models to earthquakes from North China, Pure appl. Geophys., 157, Zheng, X. & Vere-Jones, D., Further applications of the stochastic stress release model to historical earthquake data, Tectonophysics, 229, Zhuang, J. & Ma, L., The stress release model and results from modelling features of some seismic regions in China, Acta Seismol. Sin., 11, Zöller, G. & Hainzl, S., Recurrence time distributions of large earthquakes in a stochastic model for coupled fault systems: the role of fault interaction, Bull. seism. Soc. Am., 97, , doi: / Zöller, G., Hainzl, S., Ben-Zion, Y. & Holschneider, M., Earthquake activity related seismic cycles in a model for a heterogeneous strike-slip fault, Tectonophysics, 423, APPENDIX A Here we derive how the parameters α and ρ in the conditional intensity function (eq. 5) can be fixed if a Gutenberg Richter relation for the frequency magnitude relationship is assumed. In the stress release model, two processes oppose each other which both affect the stress level: stress gain due to tectonic loading and

11 658 N. M. Kuehn, S. Hainzl and F. Scherbaum stress release due to earthquakes. The system reaches an equilibrium stress state Z eq when these two processes equal each other. Tectonical loading is expressed by the loading rate ρ, which is the time-derivative of the stress function (eq. 1). The average stress release due to earthquakes is the product of the average stress loss of an earthquake Z times the conditional intensity function λ(t). Since the stress level trends to the equilibrium state Z eq, the value of the conditional intensity function at Z eq can serve as an approximation of λ(t). Hence, we get as a condition for the stable state the equation ρ = λ(z eq ) Z. (A1) However, the equilibrium state is also a good approximation for the maximum stress level, Z eq = Z max, since larger values of Z are unstable. The maximum stress level can be calculated by Z max = 10 h(mmax m0). On the other hand, the average stress loss Z can be calculated by summing up the average stress loss per magnitude for all magnitudes. Thus, Z is given by Z = mmax m 0 p(m) Z(m)dm, (A2) ln 10 b 10 bm where p(m) = is the probability density function 10 bm 0 10 bmmax of magnitudes and Z = 10 h(m m0) is the stress loss due to an earthquake of magnitude m. The solution of the integral gives b Z = (h b) 10(h b)(mmax m0) 1. (A3) 1 10 b (mmax m 0) Hence, the condition for the stable state (eq. A1) can be written as α+β b Zmax ρ = e (h b) 10(h b)(mmax m0) 1. (A4) 1 10 b (mmax m 0) Assuming that the stress level is near its stable state, the earthquake rate for all events with a magnitude m > m 0 can be constrained by λ(z max ) = 10 a bm 0. (A5) Thus, using λ(z) = e α+β Z as the conditional intensity function, this becomes e α+β Zmax = 10 a bm 0. leading to α = ln 10(a bm 0 ) β Z max. (A6) (A7) Thus, assuming a maximum magnitude as well as the a- and b- value of the Gutenberg Richter relation for a seismic active region, eqs (A7) and (A4) can be used to fix the parameters α and ρ, leaving β free. APPENDIX B In this section, the inverse transform method for calculating the next event time of the stress release model is explained. Therefore, first the conditional survivor and the conditional intensity function of a point process are defined. The conditional survivor function S n (t H τn 1 ) represents the probability that no event occurs in the interval [τ n 1, t], given the history H τn 1 of the process. It is related to the conditional probability density p n (t H τn 1 ), representing the probability of an event occurring at time t given H τn 1,via t S n (t H τn 1 ) = 1 p n (s H τn 1 )ds. (B1) τ n 1 For t (τ n 1, τ n ] the conditional intensity function can be expressed by λ(t) = p n(t H τn 1 ) S n (t H τn 1 ) = d dt log [ S n (t H τn 1 ) ], (B2) where λ(t) is also called the risk or hazard function. The inverse transform method employs that the survivor function S n (t H τn 1 ) of a point process is uniformly distributed in [0, 1]. Hence, if u is a sample from a uniform distribution on [0, 1], t = S 1 n (u) is a simulated point of the process, where S 1 n denotes the inverse of the survivor function. The hazard function for the stress release model is [ ] λ(t) = exp α + β(z 0 + ρt t i <t Z i ). (B3) Hence, using eq. (B2), the survivor function for t >τcan be calculated: [ t ] S(t H τ ) = exp e α+β Z(s) ds = exp { A[e βρ(t τ) 1] }, τ (B4) with A = 1 βρ eα+β Z(τ). So, if u is a sample from a uniform distribution on [0, 1], solving u = S(t H τ ) (B5) gives for the next event time t = τ + 1 [ βρ ln 1 ln u ]. (B6) A

arxiv:physics/ v2 [physics.geo-ph] 18 Aug 2003

arxiv:physics/ v2 [physics.geo-ph] 18 Aug 2003 Is Earthquake Triggering Driven by Small Earthquakes? arxiv:physics/0210056v2 [physics.geo-ph] 18 Aug 2003 Agnès Helmstetter Laboratoire de Géophysique Interne et Tectonophysique, Observatoire de Grenoble,

More information

Recurrence Times for Parkfield Earthquakes: Actual and Simulated. Paul B. Rundle, Donald L. Turcotte, John B. Rundle, and Gleb Yakovlev

Recurrence Times for Parkfield Earthquakes: Actual and Simulated. Paul B. Rundle, Donald L. Turcotte, John B. Rundle, and Gleb Yakovlev Recurrence Times for Parkfield Earthquakes: Actual and Simulated Paul B. Rundle, Donald L. Turcotte, John B. Rundle, and Gleb Yakovlev 1 Abstract In this paper we compare the sequence of recurrence times

More information

Impact of earthquake rupture extensions on parameter estimations of point-process models

Impact of earthquake rupture extensions on parameter estimations of point-process models 1 INTRODUCTION 1 Impact of earthquake rupture extensions on parameter estimations of point-process models S. Hainzl 1, A. Christophersen 2, and B. Enescu 1 1 GeoForschungsZentrum Potsdam, Germany; 2 Swiss

More information

The Role of Asperities in Aftershocks

The Role of Asperities in Aftershocks The Role of Asperities in Aftershocks James B. Silva Boston University April 7, 2016 Collaborators: William Klein, Harvey Gould Kang Liu, Nick Lubbers, Rashi Verma, Tyler Xuan Gu OUTLINE Introduction The

More information

A GLOBAL MODEL FOR AFTERSHOCK BEHAVIOUR

A GLOBAL MODEL FOR AFTERSHOCK BEHAVIOUR A GLOBAL MODEL FOR AFTERSHOCK BEHAVIOUR Annemarie CHRISTOPHERSEN 1 And Euan G C SMITH 2 SUMMARY This paper considers the distribution of aftershocks in space, abundance, magnitude and time. Investigations

More information

3D MODELING OF EARTHQUAKE CYCLES OF THE XIANSHUIHE FAULT, SOUTHWESTERN CHINA

3D MODELING OF EARTHQUAKE CYCLES OF THE XIANSHUIHE FAULT, SOUTHWESTERN CHINA 3D MODELING OF EARTHQUAKE CYCLES OF THE XIANSHUIHE FAULT, SOUTHWESTERN CHINA Li Xiaofan MEE09177 Supervisor: Bunichiro Shibazaki ABSTRACT We perform 3D modeling of earthquake generation of the Xianshuihe

More information

The largest aftershock: How strong, how far away, how delayed?

The largest aftershock: How strong, how far away, how delayed? GEOPHYSICAL RESEARCH LETTERS, VOL. 39,, doi:10.1029/2011gl050604, 2012 The largest aftershock: How strong, how far away, how delayed? M. Tahir, 1 J.-R. Grasso, 1 and D. Amorèse 2 Received 15 December 2011;

More information

Comment on Systematic survey of high-resolution b-value imaging along Californian faults: inference on asperities.

Comment on Systematic survey of high-resolution b-value imaging along Californian faults: inference on asperities. Comment on Systematic survey of high-resolution b-value imaging along Californian faults: inference on asperities Yavor Kamer 1, 2 1 Swiss Seismological Service, ETH Zürich, Switzerland 2 Chair of Entrepreneurial

More information

Coulomb stress changes due to Queensland earthquakes and the implications for seismic risk assessment

Coulomb stress changes due to Queensland earthquakes and the implications for seismic risk assessment Coulomb stress changes due to Queensland earthquakes and the implications for seismic risk assessment Abstract D. Weatherley University of Queensland Coulomb stress change analysis has been applied in

More information

arxiv:physics/ v1 6 Aug 2006

arxiv:physics/ v1 6 Aug 2006 The application of the modified form of Båth s law to the North Anatolian Fault Zone arxiv:physics/0608064 v1 6 Aug 2006 1. INTRODUCTION S E Yalcin, M L Kurnaz Department of Physics, Bogazici University,

More information

Distribution of volcanic earthquake recurrence intervals

Distribution of volcanic earthquake recurrence intervals JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 114,, doi:10.1029/2008jb005942, 2009 Distribution of volcanic earthquake recurrence intervals M. Bottiglieri, 1 C. Godano, 1 and L. D Auria 2 Received 21 July 2008;

More information

Simulated and Observed Scaling in Earthquakes Kasey Schultz Physics 219B Final Project December 6, 2013

Simulated and Observed Scaling in Earthquakes Kasey Schultz Physics 219B Final Project December 6, 2013 Simulated and Observed Scaling in Earthquakes Kasey Schultz Physics 219B Final Project December 6, 2013 Abstract Earthquakes do not fit into the class of models we discussed in Physics 219B. Earthquakes

More information

Probabilistic seismic hazard estimation in low-seismicity regions considering non-poissonian seismic occurrence

Probabilistic seismic hazard estimation in low-seismicity regions considering non-poissonian seismic occurrence Geophys. J. Int. (26) 164, 543 55 doi: 1.1111/j.1365-246X.26.2863.x Probabilistic seismic hazard estimation in low-seismicity regions considering non-poissonian seismic occurrence Céline Beauval, Sebastian

More information

Earthquake predictability measurement: information score and error diagram

Earthquake predictability measurement: information score and error diagram Earthquake predictability measurement: information score and error diagram Yan Y. Kagan Department of Earth and Space Sciences University of California, Los Angeles, California, USA August, 00 Abstract

More information

COULOMB STRESS CHANGES DUE TO RECENT ACEH EARTHQUAKES

COULOMB STRESS CHANGES DUE TO RECENT ACEH EARTHQUAKES COULOMB STRESS CHANGES DUE TO RECENT ACEH EARTHQUAKES Madlazim Physics Department, Faculty Mathematics and Sciences of Surabaya State University (UNESA) Jl. Ketintang, Surabaya 60231, Indonesia. e-mail:

More information

On the nucleation of creep and the interaction between creep and seismic slip on rate- and state-dependent faults

On the nucleation of creep and the interaction between creep and seismic slip on rate- and state-dependent faults Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 34, L15303, doi:10.1029/2007gl030337, 2007 On the nucleation of creep and the interaction between creep and seismic slip on rate- and state-dependent

More information

Performance of national scale smoothed seismicity estimates of earthquake activity rates. Abstract

Performance of national scale smoothed seismicity estimates of earthquake activity rates. Abstract Performance of national scale smoothed seismicity estimates of earthquake activity rates Jonathan Griffin 1, Graeme Weatherill 2 and Trevor Allen 3 1. Corresponding Author, Geoscience Australia, Symonston,

More information

Does Aftershock Duration Scale With Mainshock Size?

Does Aftershock Duration Scale With Mainshock Size? GEOPHYSICAL RESEARCH LETTERS, VOL.???, NO., PAGES 1 16, Does Aftershock Duration Scale With Mainshock Size? A. Ziv A. Ziv, Ben-Gurion University of the Negev, Beer-Sheva, 84105, Israel. (e-mail: zival@bgu.ac.il)

More information

Comment on the paper. Layered seismogenic source model and probabilistic seismic-hazard analyses in. Central Italy

Comment on the paper. Layered seismogenic source model and probabilistic seismic-hazard analyses in. Central Italy Comment on the paper Layered seismogenic source model and probabilistic seismic-hazard analyses in Central Italy by Pace B., L. Peruzza, G. Lavecchia, P. Boncio. BSSA, vol. 96, No. 1, pp. 107-132, February

More information

Time-varying and long-term mean aftershock hazard in Wellington

Time-varying and long-term mean aftershock hazard in Wellington Time-varying and long-term mean aftershock hazard in Wellington A. Christophersen, D.A. Rhoades, R.J. Van Dissen, C. Müller, M.W. Stirling, G.H. McVerry & M.C. Gerstenberger GNS Science, Lower Hutt, New

More information

Pattern Dynamics and Forecast Methods in Seismically Active Regions

Pattern Dynamics and Forecast Methods in Seismically Active Regions Pattern Dynamics and Forecast Methods in Seismically Active Regions K.F. Tiampo (1), J.B. Rundle (1), S. McGinnis (1) and W. Klein (2) (1) CIRES, University of Colorado, Boulder, CO USA (e-mail: kristy@fractal.colorado.edu;

More information

Minimum preshock magnitude in critical regions of accelerating seismic crustal deformation

Minimum preshock magnitude in critical regions of accelerating seismic crustal deformation Bollettino di Geofisica Teorica ed Applicata Vol. 44, n. 2, pp. 103-113; June 2003 Minimum preshock magnitude in critical regions of accelerating seismic crustal deformation C.B. Papazachos University

More information

PostScript le created: August 6, 2006 time 839 minutes

PostScript le created: August 6, 2006 time 839 minutes GEOPHYSICAL RESEARCH LETTERS, VOL.???, XXXX, DOI:10.1029/, PostScript le created: August 6, 2006 time 839 minutes Earthquake predictability measurement: information score and error diagram Yan Y. Kagan

More information

Aftershock From Wikipedia, the free encyclopedia

Aftershock From Wikipedia, the free encyclopedia Page 1 of 5 Aftershock From Wikipedia, the free encyclopedia An aftershock is a smaller earthquake that occurs after a previous large earthquake, in the same area of the main shock. If an aftershock is

More information

Multifractal Analysis of Seismicity of Kutch Region (Gujarat)

Multifractal Analysis of Seismicity of Kutch Region (Gujarat) P-382 Multifractal Analysis of Seismicity of Kutch Region (Gujarat) Priyanka Midha *, IIT Roorkee Summary The geographical distribution of past earthquakes is not uniform over the globe. Also, in one region

More information

Accelerating energy release prior to large events in simulated earthquake cycles: implications for earthquake forecasting

Accelerating energy release prior to large events in simulated earthquake cycles: implications for earthquake forecasting Accelerating energy release prior to large events in simulated earthquake cycles: implications for earthquake forecasting Peter Mora and David Place QUAKES, Department of Earth Sciences, The University

More information

Renewal models of seismic recurrence applied to paleoseismological. data. Abstract. I. Mosca 1, R. Console 2,3 and G. D Addezio 2

Renewal models of seismic recurrence applied to paleoseismological. data. Abstract. I. Mosca 1, R. Console 2,3 and G. D Addezio 2 * Manuscript Click here to download Manuscript: manuscript.doc Renewal models of seismic recurrence applied to paleoseismological data I. Mosca 1, R. Console 2,3 and G. D Addezio 2 1 Utrecht University,

More information

Quantifying early aftershock activity of the 2004 mid-niigata Prefecture earthquake (M w 6.6)

Quantifying early aftershock activity of the 2004 mid-niigata Prefecture earthquake (M w 6.6) Click Here for Full Article JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 112,, doi:10.1029/2006jb004629, 2007 Quantifying early aftershock activity of the 2004 mid-niigata Prefecture earthquake (M w 6.6) Bogdan

More information

Testing aftershock models on a time-scale of decades

Testing aftershock models on a time-scale of decades Testing aftershock models on a time-scale of decades A. Christophersen D.A. Rhoades S. Hainzl D.S. Harte GNS Science Consultancy Report 6/6 March 6 DISCLAIMER This report has been prepared by the Institute

More information

Adaptive Kernel Estimation and Continuous Probability Representation of Historical Earthquake Catalogs

Adaptive Kernel Estimation and Continuous Probability Representation of Historical Earthquake Catalogs Bulletin of the Seismological Society of America, Vol. 92, No. 3, pp. 904 912, April 2002 Adaptive Kernel Estimation and Continuous Probability Representation of Historical Earthquake Catalogs by Christian

More information

Journal of Asian Earth Sciences

Journal of Asian Earth Sciences Journal of Asian Earth Sciences 7 (0) 09 8 Contents lists available at SciVerse ScienceDirect Journal of Asian Earth Sciences journal homepage: www.elsevier.com/locate/jseaes Maximum magnitudes in aftershock

More information

Hitoshi Hirose (1), and Kazuro Hirahara (2) Abstract. Introduction

Hitoshi Hirose (1), and Kazuro Hirahara (2) Abstract. Introduction Three dimensional simulation for the earthquake cycle at a subduction zone based on a rate- and state-dependent friction law: Insight into a finiteness and a variety of dip-slip earthquakes Hitoshi Hirose

More information

Originally published as:

Originally published as: Originally published as: Roth, F., Dahm, T., Hainzl, S. (2017): Testing stress shadowing effects at the South American subduction zone. Geophysical Journal International, 211, 2, pp. 1272 1283. DOI: http://doi.org/10.1093/gji/ggx362

More information

PostScript file created: June 11, 2011; time 985 minutes RANDOM STRESS AND OMORI S LAW. Yan Y. Kagan

PostScript file created: June 11, 2011; time 985 minutes RANDOM STRESS AND OMORI S LAW. Yan Y. Kagan PostScript file created: June 11, 2011; time 985 minutes RANDOM STRESS AND OMORI S LAW Yan Y. Kagan Department of Earth and Space Sciences, University of California, Los Angeles, California 90095-1567,

More information

Seismic Hazard Assessment for Specified Area

Seismic Hazard Assessment for Specified Area ESTIATION OF AXIU REGIONAL AGNITUDE m At present there is no generally accepted method for estimating the value of the imum regional magnitude m. The methods for evaluating m fall into two main categories:

More information

Mechanical origin of aftershocks: Supplementary Information

Mechanical origin of aftershocks: Supplementary Information Mechanical origin of aftershocks: Supplementary Information E. Lippiello Department of Mathematics and Physics, Second University of Naples, Via Vivaldi 43, 81100 Caserta, Italy & Kavli Institute for Theoretical

More information

6 Source Characterization

6 Source Characterization 6 Source Characterization Source characterization describes the rate at which earthquakes of a given magnitude, and dimensions (length and width) occur at a given location. For each seismic source, the

More information

Source parameters II. Stress drop determination Energy balance Seismic energy and seismic efficiency The heat flow paradox Apparent stress drop

Source parameters II. Stress drop determination Energy balance Seismic energy and seismic efficiency The heat flow paradox Apparent stress drop Source parameters II Stress drop determination Energy balance Seismic energy and seismic efficiency The heat flow paradox Apparent stress drop Source parameters II: use of empirical Green function for

More information

Forecasting Earthquakes

Forecasting Earthquakes Forecasting Earthquakes Lecture 9 Earthquake Recurrence ) Long-term prediction - Elastic Rebound Theory Stress Stress & strain accumulation 0 0 4 F slow accumulation of stress & strain that deforms rock

More information

Chien-chih Chen Institute of Geophysics, National Central University, Chungli, Taiwan

Chien-chih Chen Institute of Geophysics, National Central University, Chungli, Taiwan Geophys. J. Int. (2003) 155, F1 F5 FAST TRACK PAPER Accelerating seismicity of moderate-size earthquakes before the 1999 Chi-Chi, Taiwan, earthquake: Testing time-prediction of the self-organizing spinodal

More information

Quantifying the effect of declustering on probabilistic seismic hazard

Quantifying the effect of declustering on probabilistic seismic hazard Proceedings of the Ninth Pacific Conference on Earthquake Engineering Building an Earthquake-Resilient Society 14-16 April, 2011, Auckland, New Zealand Quantifying the effect of declustering on probabilistic

More information

Synthetic Seismicity Models of Multiple Interacting Faults

Synthetic Seismicity Models of Multiple Interacting Faults Synthetic Seismicity Models of Multiple Interacting Faults Russell Robinson and Rafael Benites Institute of Geological & Nuclear Sciences, Box 30368, Lower Hutt, New Zealand (email: r.robinson@gns.cri.nz).

More information

Y. Y. Kagan and L. Knopoff Institute of Geophysics and Planetary Physics, University of California, Los Angeles, California 90024, USA

Y. Y. Kagan and L. Knopoff Institute of Geophysics and Planetary Physics, University of California, Los Angeles, California 90024, USA Geophys. J. R. astr. SOC. (1987) 88, 723-731 Random stress and earthquake statistics: time dependence Y. Y. Kagan and L. Knopoff Institute of Geophysics and Planetary Physics, University of California,

More information

A rate-state model for aftershocks triggered by dislocation on a rectangular fault: a review and new insights

A rate-state model for aftershocks triggered by dislocation on a rectangular fault: a review and new insights A rate-state model for aftershocks triggered by dislocation on a rectangular fault: a review and new insights Rodolfo Console and Flaminia Catalli Istituto Nazionale di Geofisica e Vulcanologia, Rome,

More information

Effect of an outer-rise earthquake on seismic cycle of large interplate earthquakes estimated from an instability model based on friction mechanics

Effect of an outer-rise earthquake on seismic cycle of large interplate earthquakes estimated from an instability model based on friction mechanics Effect of an outer-rise earthquake on seismic cycle of large interplate earthquakes estimated from an instability model based on friction mechanics Naoyuki Kato (1) and Tomowo Hirasawa (2) (1) Geological

More information

Appendix O: Gridded Seismicity Sources

Appendix O: Gridded Seismicity Sources Appendix O: Gridded Seismicity Sources Peter M. Powers U.S. Geological Survey Introduction The Uniform California Earthquake Rupture Forecast, Version 3 (UCERF3) is a forecast of earthquakes that fall

More information

M 7.0 earthquake recurrence on the San Andreas fault from a stress renewal model

M 7.0 earthquake recurrence on the San Andreas fault from a stress renewal model Click Here for Full Article JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 111,, doi:10.1029/2006jb004415, 2006 M 7.0 earthquake recurrence on the San Andreas fault from a stress renewal model Tom Parsons 1 Received

More information

THE DOUBLE BRANCHING MODEL FOR EARTHQUAKE FORECAST APPLIED TO THE JAPANESE SEISMICITY

THE DOUBLE BRANCHING MODEL FOR EARTHQUAKE FORECAST APPLIED TO THE JAPANESE SEISMICITY 1 2 3 THE DOUBLE BRANCHING MODEL FOR EARTHQUAKE FORECAST APPLIED TO THE JAPANESE SEISMICITY 4 Anna Maria Lombardi Istituto Nazionale di Geofisica e Vulcanologia,Via di Vigna Murata 605, 00143 Roma, Italy,

More information

ALM: An Asperity-based Likelihood Model for California

ALM: An Asperity-based Likelihood Model for California ALM: An Asperity-based Likelihood Model for California Stefan Wiemer and Danijel Schorlemmer Stefan Wiemer and Danijel Schorlemmer 1 ETH Zürich, Switzerland INTRODUCTION In most earthquake hazard models,

More information

Generic Earthquake Simulator

Generic Earthquake Simulator Generic Earthquake Simulator by Terry E. Tullis, Keith Richards-Dinger, Michael Barall, James H. Dieterich, Edward H. Field, Eric Heien, Louise H. Kellogg, Fred Pollitz, John Rundle, Michael Sachs, Donald

More information

GLOBAL SOURCE PARAMETERS OF FINITE FAULT MODEL FOR STRONG GROUND MOTION SIMULATIONS OR PREDICTIONS

GLOBAL SOURCE PARAMETERS OF FINITE FAULT MODEL FOR STRONG GROUND MOTION SIMULATIONS OR PREDICTIONS 13 th orld Conference on Earthquake Engineering Vancouver, B.C., Canada August 1-6, 2004 Paper No. 2743 GLOBAL SOURCE PARAMETERS OF FINITE FAULT MODEL FOR STRONG GROUND MOTION SIMULATIONS OR PREDICTIONS

More information

Comparison of Short-Term and Time-Independent Earthquake Forecast Models for Southern California

Comparison of Short-Term and Time-Independent Earthquake Forecast Models for Southern California Bulletin of the Seismological Society of America, Vol. 96, No. 1, pp. 90 106, February 2006, doi: 10.1785/0120050067 Comparison of Short-Term and Time-Independent Earthquake Forecast Models for Southern

More information

Mid-Continent Earthquakes As A Complex System

Mid-Continent Earthquakes As A Complex System SRL complex earthquakes 5/22/09 1 Mid-Continent Earthquakes As A Complex System Niels Bohr once observed How wonderful that we have met with a paradox. Now we have some hope of making progress. This situation

More information

UNIVERSITY OF CALGARY. Nontrivial Decay of Aftershock Density With Distance in Southern California. Javad Moradpour Taleshi A THESIS

UNIVERSITY OF CALGARY. Nontrivial Decay of Aftershock Density With Distance in Southern California. Javad Moradpour Taleshi A THESIS UNIVERSITY OF CALGARY Nontrivial Decay of Aftershock Density With Distance in Southern California by Javad Moradpour Taleshi A THESIS SUBMITTED TO THE FACULTY OF GRADUATE STUDIES IN PARTIAL FULFILLMENT

More information

ETH Swiss Federal Institute of Technology Zürich

ETH Swiss Federal Institute of Technology Zürich Swiss Federal Institute of Technology Zürich Earthquake Statistics using ZMAP Recent Results Danijel Schorlemmer, Stefan Wiemer Zürich, Swiss Seismological Service, Switzerland Contributions by: Matt Gerstenberger

More information

Friction. Why friction? Because slip on faults is resisted by frictional forces.

Friction. Why friction? Because slip on faults is resisted by frictional forces. Friction Why friction? Because slip on faults is resisted by frictional forces. We first describe the results of laboratory friction experiments, and then discuss the implications of the friction constitutive

More information

A TESTABLE FIVE-YEAR FORECAST OF MODERATE AND LARGE EARTHQUAKES. Yan Y. Kagan 1,David D. Jackson 1, and Yufang Rong 2

A TESTABLE FIVE-YEAR FORECAST OF MODERATE AND LARGE EARTHQUAKES. Yan Y. Kagan 1,David D. Jackson 1, and Yufang Rong 2 Printed: September 1, 2005 A TESTABLE FIVE-YEAR FORECAST OF MODERATE AND LARGE EARTHQUAKES IN SOUTHERN CALIFORNIA BASED ON SMOOTHED SEISMICITY Yan Y. Kagan 1,David D. Jackson 1, and Yufang Rong 2 1 Department

More information

An Empirical Model for Earthquake Probabilities in the San Francisco Bay Region, California,

An Empirical Model for Earthquake Probabilities in the San Francisco Bay Region, California, Bulletin of the Seismological Society of America, Vol. 93, No. 1, pp. 1 13, February 2003 An Empirical Model for Earthquake Probabilities in the San Francisco Bay Region, California, 2002 2031 by Paul

More information

Interactions between earthquakes and volcano activity

Interactions between earthquakes and volcano activity GEOPHYSICAL RESEARCH LETTERS, VOL. 34, L24303, doi:10.1029/2007gl031438, 2007 Interactions between earthquakes and volcano activity Nicolas Lemarchand 1 and Jean-Robert Grasso 1 Received 29 July 2007;

More information

PostScript file created: April 24, 2012; time 860 minutes WHOLE EARTH HIGH-RESOLUTION EARTHQUAKE FORECASTS. Yan Y. Kagan and David D.

PostScript file created: April 24, 2012; time 860 minutes WHOLE EARTH HIGH-RESOLUTION EARTHQUAKE FORECASTS. Yan Y. Kagan and David D. PostScript file created: April 24, 2012; time 860 minutes WHOLE EARTH HIGH-RESOLUTION EARTHQUAKE FORECASTS Yan Y. Kagan and David D. Jackson Department of Earth and Space Sciences, University of California,

More information

Di#erences in Earthquake Source and Ground Motion Characteristics between Surface and Buried Crustal Earthquakes

Di#erences in Earthquake Source and Ground Motion Characteristics between Surface and Buried Crustal Earthquakes Bull. Earthq. Res. Inst. Univ. Tokyo Vol. 2+,**0 pp.,/3,00 Di#erences in Earthquake Source and Ground Motion Characteristics between Surface and Buried Crustal Earthquakes Paul Somerville* and Arben Pitarka

More information

Rate and State-Dependent Friction in Earthquake Simulation

Rate and State-Dependent Friction in Earthquake Simulation Rate and State-Dependent Friction in Earthquake Simulation Zac Meadows UC Davis - Department of Physics Summer 2012 REU September 3, 2012 Abstract To better understand the spatial and temporal complexity

More information

Theory of earthquake recurrence times

Theory of earthquake recurrence times Theory of earthquake recurrence times Alex SAICHEV1,2 and Didier SORNETTE1,3 1ETH Zurich, Switzerland 2Mathematical Department, Nizhny Novgorod State University, Russia. 3Institute of Geophysics and Planetary

More information

Ground displacement in a fault zone in the presence of asperities

Ground displacement in a fault zone in the presence of asperities BOLLETTINO DI GEOFISICA TEORICA ED APPLICATA VOL. 40, N. 2, pp. 95-110; JUNE 2000 Ground displacement in a fault zone in the presence of asperities S. SANTINI (1),A.PIOMBO (2) and M. DRAGONI (2) (1) Istituto

More information

Heterogeneous Coulomb stress perturbation during earthquake cycles in a 3D rate-and-state fault model

Heterogeneous Coulomb stress perturbation during earthquake cycles in a 3D rate-and-state fault model Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 35, L21306, doi:10.1029/2008gl035614, 2008 Heterogeneous Coulomb stress perturbation during earthquake cycles in a 3D rate-and-state fault

More information

On the earthquake predictability of fault interaction models

On the earthquake predictability of fault interaction models On the earthquake predictability of fault interaction models Warner Marzocchi, and Daniele Melini Istituto Nazionale di Geofisica e Vulcanologia, Rome, Italy StatSei 9 meeting, Potsdam, June 15-17, 2015

More information

San Francisco Bay Area Earthquake Simulations: A step toward a Standard Physical Earthquake Model

San Francisco Bay Area Earthquake Simulations: A step toward a Standard Physical Earthquake Model San Francisco Bay Area Earthquake Simulations: A step toward a Standard Physical Earthquake Model Steven N. Ward Institute of Geophysics and Planetary Physics, University of California, Santa Cruz, CA,

More information

A mathematical formulation of accelerating moment release based on the stress accumulation model

A mathematical formulation of accelerating moment release based on the stress accumulation model Click Here for Full Article JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 112,, doi:10.1029/2006jb004671, 2007 A mathematical formulation of accelerating moment release based on the stress accumulation model A.

More information

Velocity Changes of Seismic Waves and Monitoring Stress in the Crust

Velocity Changes of Seismic Waves and Monitoring Stress in the Crust Bull. Earthq. Res. Inst. Univ. Tokyo Vol. 12,**- pp.,+3,,0 + * +, +, Velocity Changes of Seismic Waves and Monitoring Stress in the Crust Muneyoshi Furumoto + *, Yoshihiro Hiramatsu + and Takashi Satoh,

More information

EFFECTS OF NON-LINEAR WEAKENING ON EARTHQUAKE SOURCE SCALINGS

EFFECTS OF NON-LINEAR WEAKENING ON EARTHQUAKE SOURCE SCALINGS Extended abstract for the 11th International Conference on Fracture 2005 1 EFFECTS OF NON-LINEAR WEAKENING ON EARTHQUAKE SOURCE SCALINGS J.-P. Ampuero Geosciences Department, Princeton University, USA

More information

UCERF3 Task R2- Evaluate Magnitude-Scaling Relationships and Depth of Rupture: Proposed Solutions

UCERF3 Task R2- Evaluate Magnitude-Scaling Relationships and Depth of Rupture: Proposed Solutions UCERF3 Task R2- Evaluate Magnitude-Scaling Relationships and Depth of Rupture: Proposed Solutions Bruce E. Shaw Lamont Doherty Earth Observatory, Columbia University Statement of the Problem In UCERF2

More information

Observational analysis of correlations between aftershock productivities and regional conditions in the context of a damage rheology model

Observational analysis of correlations between aftershock productivities and regional conditions in the context of a damage rheology model Geophys. J. Int. (2009) 177, 481 490 doi: 10.1111/j.1365-246X.2009.04145.x FAST TRACK PAPER Observational analysis of correlations between aftershock productivities and regional conditions in the context

More information

Dependence of Possible Characteristic Earthquakes on Spatial Samph'ng: Illustration for the Wasatch Seismic Zone, Utah

Dependence of Possible Characteristic Earthquakes on Spatial Samph'ng: Illustration for the Wasatch Seismic Zone, Utah Dependence of Possible Characteristic Earthquakes on Spatial Samph'ng: Illustration for the Wasatch Seismic Zone, Utah Seth Stein Northwestern University Anke Friedrich Universitiit Potsdam* Andrew Newman

More information

The magnitude distribution of earthquakes near Southern California faults

The magnitude distribution of earthquakes near Southern California faults JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 116,, doi:10.1029/2010jb007933, 2011 The magnitude distribution of earthquakes near Southern California faults Morgan T. Page, 1 David Alderson, 2 and John Doyle 3

More information

Module 7 SEISMIC HAZARD ANALYSIS (Lectures 33 to 36)

Module 7 SEISMIC HAZARD ANALYSIS (Lectures 33 to 36) Lecture 35 Topics Module 7 SEISMIC HAZARD ANALYSIS (Lectures 33 to 36) 7.4.4 Predictive Relationships 7.4.5 Temporal Uncertainty 7.4.6 Poisson Model 7.4.7 Other Models 7.4.8 Model Applicability 7.4.9 Probability

More information

Scaling of apparent stress from broadband radiated energy catalogue and seismic moment catalogue and its focal mechanism dependence

Scaling of apparent stress from broadband radiated energy catalogue and seismic moment catalogue and its focal mechanism dependence Earth Planets Space, 53, 943 948, 2001 Scaling of apparent stress from broadband radiated energy catalogue and seismic moment catalogue and its focal mechanism dependence Z. L. Wu Institute of Geophysics,

More information

Geophysical Journal International

Geophysical Journal International Geophysical Journal International Geophys. J. Int. (2011) 184, 759 776 doi: 10.1111/j.1365-246X.2010.04857.x Global earthquake forecasts Yan Y. Kagan and David D. Jackson Department of Earth and Space

More information

Numerical study on multi-scaling earthquake rupture

Numerical study on multi-scaling earthquake rupture GEOPHYSICAL RESEARCH LETTERS, VOL. 31, L266, doi:1.129/23gl1878, 24 Numerical study on multi-scaling earthquake rupture Hideo Aochi Institut de Radioprotection et de Sûreté Nucléaire, France Satoshi Ide

More information

From the Testing Center of Regional Earthquake Likelihood Models. to the Collaboratory for the Study of Earthquake Predictability

From the Testing Center of Regional Earthquake Likelihood Models. to the Collaboratory for the Study of Earthquake Predictability From the Testing Center of Regional Earthquake Likelihood Models (RELM) to the Collaboratory for the Study of Earthquake Predictability (CSEP) Danijel Schorlemmer, Matt Gerstenberger, Tom Jordan, Dave

More information

On the validity of time-predictable model for earthquake generation in north-east India

On the validity of time-predictable model for earthquake generation in north-east India Proc. Indian Acad. Sci. (Earth Planet. Sci.), Vol. 101, No. 4, December 1992, pp. 361-368. 9 Printed in India. On the validity of time-predictable model for earthquake generation in north-east India V

More information

Numerical simulation of seismic cycles at a subduction zone with a laboratory-derived friction law

Numerical simulation of seismic cycles at a subduction zone with a laboratory-derived friction law Numerical simulation of seismic cycles at a subduction zone with a laboratory-derived friction law Naoyuki Kato (1), Kazuro Hirahara (2) and Mikio Iizuka (3) (1) Earthquake Research Institute, University

More information

A forward test of the Decelerating-Accelerating Seismic Strain model in the Mediterranean

A forward test of the Decelerating-Accelerating Seismic Strain model in the Mediterranean Bollettino di Geofisica Teorica ed Applicata Vol. 50, n. 3, pp. 235-254; September 2009 A forward test of the Decelerating-Accelerating Seismic Strain model in the Mediterranean B.C. PAPAZACHOS, G.F. KARAKAISIS,

More information

Aspects of risk assessment in power-law distributed natural hazards

Aspects of risk assessment in power-law distributed natural hazards Natural Hazards and Earth System Sciences (2004) 4: 309 313 SRef-ID: 1684-9981/nhess/2004-4-309 European Geosciences Union 2004 Natural Hazards and Earth System Sciences Aspects of risk assessment in power-law

More information

Sensitivity study of forecasted aftershock seismicity based on Coulomb stress calculation and rate and state dependent frictional response

Sensitivity study of forecasted aftershock seismicity based on Coulomb stress calculation and rate and state dependent frictional response Click Here for Full Article JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 115,, doi:10.1029/2009jb006838, 2010 Sensitivity study of forecasted aftershock seismicity based on Coulomb stress calculation and rate

More information

Because of its reputation of validity over a wide range of

Because of its reputation of validity over a wide range of The magnitude distribution of declustered earthquakes in Southern California Leon Knopoff* Department of Physics and Astronomy and Institute of Geophysics and Planetary Physics, University of California,

More information

Preparation of a Comprehensive Earthquake Catalog for Northeast India and its completeness analysis

Preparation of a Comprehensive Earthquake Catalog for Northeast India and its completeness analysis IOSR Journal of Applied Geology and Geophysics (IOSR-JAGG) e-issn: 2321 0990, p-issn: 2321 0982.Volume 2, Issue 6 Ver. I (Nov-Dec. 2014), PP 22-26 Preparation of a Comprehensive Earthquake Catalog for

More information

A double branching model for earthquake occurrence

A double branching model for earthquake occurrence Click Here for Full Article JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 113,, doi:10.1029/2007jb005472, 2008 A double branching model for earthquake occurrence Warner Marzocchi 1 and Anna Maria Lombardi 1 Received

More information

Clock Advance in Earthquake Triggering of Volcanic Eruptions

Clock Advance in Earthquake Triggering of Volcanic Eruptions Clock Advance in Earthquake Triggering of Volcanic Eruptions Mark Bebbington 1 and Warner Marzocchi 2 1 Volcanic Risk Solutions, Massey University, Palmerston North, New Zealand 2 Istituto Nazionale di

More information

The Centenary of the Omori Formula for a Decay Law of Aftershock Activity

The Centenary of the Omori Formula for a Decay Law of Aftershock Activity The Centenary of the Omori Formula for a Decay Law of Aftershock Activity Author; Tokuji Utsu, Yosihiko Ogata, and Ritsuko S. Matsu'ura Presentater; Okuda Takashi 8. p Values from Superposed Sequences

More information

Depth variation of coseismic stress drop explains bimodal earthquake magnitude-frequency distribution

Depth variation of coseismic stress drop explains bimodal earthquake magnitude-frequency distribution Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 35, L24301, doi:10.1029/2008gl036249, 2008 Depth variation of coseismic stress drop explains bimodal earthquake magnitude-frequency distribution

More information

Supporting Information for Break of slope in earthquake-size distribution reveals creep rate along the San Andreas fault system

Supporting Information for Break of slope in earthquake-size distribution reveals creep rate along the San Andreas fault system GEOPHYSICAL RESEARCH LETTERS Supporting Information for Break of slope in earthquake-size distribution reveals creep rate along the San Andreas fault system Inessa Vorobieva, 1,2 Peter Shebalin, 1,2 Clément

More information

ON NEAR-FIELD GROUND MOTIONS OF NORMAL AND REVERSE FAULTS FROM VIEWPOINT OF DYNAMIC RUPTURE MODEL

ON NEAR-FIELD GROUND MOTIONS OF NORMAL AND REVERSE FAULTS FROM VIEWPOINT OF DYNAMIC RUPTURE MODEL 1 Best Practices in Physics-based Fault Rupture Models for Seismic Hazard Assessment of Nuclear ON NEAR-FIELD GROUND MOTIONS OF NORMAL AND REVERSE FAULTS FROM VIEWPOINT OF DYNAMIC RUPTURE MODEL Hideo AOCHI

More information

Seismic Quiescence before the 1999 Chi-Chi, Taiwan, M w 7.6 Earthquake

Seismic Quiescence before the 1999 Chi-Chi, Taiwan, M w 7.6 Earthquake Bulletin of the Seismological Society of America, Vol. 96, No. 1, pp. 321 327, February 2006, doi: 10.1785/0120050069 Seismic Quiescence before the 1999 Chi-Chi, Taiwan, M w 7.6 Earthquake by Yih-Min Wu

More information

Accelerating Seismic Crustal Deformation in the Southern Aegean Area

Accelerating Seismic Crustal Deformation in the Southern Aegean Area Bulletin of the Seismological Society of America, Vol. 92, No. 2, pp. 570 580, March 2002 Accelerating Seismic Crustal Deformation in the Southern Aegean Area by C. B. Papazachos, G. F. Karakaisis, A.

More information

Damage Localization, Sensitivity of Energy Release and Catastrophe Transition

Damage Localization, Sensitivity of Energy Release and Catastrophe Transition Damage Localization, Sensitivity of Energy Release and Catastrophe Transition Y.L.Bai (1), H.L.Li (1), F.J.Ke (1,2) and M.F.Xia (1,3) (1) LNM, Institute of Mechanics, Chinese Academy of Sciences, Beijing

More information

Probabilities for Jumping Fault Segment Stepovers

Probabilities for Jumping Fault Segment Stepovers GEOPHYSICAL RESEARCH LETTERS, VOL.???, XXXX, DOI:1.129/, Probabilities for Jumping Fault Segment Stepovers Bruce E. Shaw Lamont-Doherty Earth Observatory, Columbia University, New York James H. Dieterich

More information

A. Akinci 1, M. Murru 1, R. Console 2, G. Falcone 1, S. Pucci 1 1. Istituto Nazionale di Geofisica e Vulcanologia, Rome, Italy 2

A. Akinci 1, M. Murru 1, R. Console 2, G. Falcone 1, S. Pucci 1 1. Istituto Nazionale di Geofisica e Vulcanologia, Rome, Italy 2 Earthquake Forecasting Models and their Impact on the Ground Motion Hazard in the Marmara Region, Turkey A. Akinci 1, M. Murru 1, R. Console 2, G. Falcone 1, S. Pucci 1 1 Istituto Nazionale di Geofisica

More information

Mechanics of Earthquakes and Faulting

Mechanics of Earthquakes and Faulting Mechanics of Earthquakes and Faulting Lecture 20, 30 Nov. 2017 www.geosc.psu.edu/courses/geosc508 Seismic Spectra & Earthquake Scaling laws. Seismic Spectra & Earthquake Scaling laws. Aki, Scaling law

More information

Scale-free network of earthquakes

Scale-free network of earthquakes Scale-free network of earthquakes Sumiyoshi Abe 1 and Norikazu Suzuki 2 1 Institute of Physics, University of Tsukuba, Ibaraki 305-8571, Japan 2 College of Science and Technology, Nihon University, Chiba

More information

Stress triggering and earthquake probability estimates

Stress triggering and earthquake probability estimates JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 109,, doi:10.1029/2003jb002437, 2004 Stress triggering and earthquake probability estimates Jeanne L. Hardebeck 1 Institute for Geophysics and Planetary Physics, Scripps

More information