Assessing the dependency between the magnitudes of earthquakes and the magnitudes of their aftershocks
|
|
- Philip Dalton
- 6 years ago
- Views:
Transcription
1 Research Article Received: 3 September 2013, Revised: 30 January 2014, Accepted: 31 January 2014, Published online in Wiley Online Library (wileyonlinelibrary.com) DOI: /env.2268 Assessing the dependency between the magnitudes of earthquakes and the magnitudes of their aftershocks Kevin Nichols a * and Frederic Paik Schoenberg b Current assumptions on the level of inter-dependency between the magnitudes of earthquakes and the magnitudes of their aftershocks vary significantly, ranging from assumed independence for models like epidemic type aftershock-sequence (ETAS) to assumed positive correlation, at least for the events leading up to the largest earthquakes, for models like accelerated moment release. Recently, a method was developed by Marsan and Lengline for model independent stochastic de-clustering (MISD) of a multidimensional Hawkes self-exciting process, or as is the case here, earthquake activity. After expanding the algorithm to allow for a spatially dependent background rate (henceforth SDMISD), repeated application of SDMISD to a global earthquake catalog allows for estimates with confidence bounds for the average magnitude of an aftershock conditioned upon the magnitude(s) of the earthquake(s) that caused the aftershock. 15 years of data on M 5.3 earthquakes from the global Centroid Moment Tensor earthquake catalog indicate that the magnitudes of aftershocks are dependent on the magnitudes of the earthquakes that cause them, and indeed, larger than average earthquakes are more likely to trigger larger than average aftershocks. Copyright 2014 John Wiley & Sons, Ltd. Keywords: earthquake magnitude; mainshock; aftershock; stochastic de-clustering; ETAS; MISD: SDMISD 1. INTRODUCTION The epidemic type aftershock-sequence (ETAS) model (Ogata 1988, 1998) parametrically estimates the rates of earthquake occurrences over a spatial temporal support. The ETAS model identifies two sources for earthquakes. They are either mainshocks, with rate determined by a spatially dependent background rate, or they are aftershocks, with rate determined by the spatial and temporal coordinates, as well as the magnitudes, of previously occurring earthquakes. Assumedly, both mainshocks and aftershocks are capable of triggering additional aftershocks. Ogata (1998) established the following extension of a Hawkes self-exciting point-process model (Hawkes, 1971) to model the spatial temporal conditional rate of earthquake occurrence:.t; x; y; H t / D.x; y/ C X g.t t i ;x x i ;y y i I M i / (1) iwt i <t where H t represents the history of earthquakes occurring prior to time.t/,.x; y/ represents the background rate of a mainshock earthquake occurring over the spatial temporal support, and P g.t t i ;x x i ;y y i I M i / represents the rate at which aftershock earthquakes occur. Parameterization of g.:/ is often difficult when considering that it is difficult to determine which earthquakes are aftershocks. Once de-clustered, nonparametric methods for estimating g.:/ will prove helpful in determining the level of dependency between earthquake magnitudes and the magnitudes of their aftershocks. One key assumption of most earthquake models is that the magnitudes of earthquakes roughly follow the exponential distribution (f ) (Gutenberg and Richter, 1944), meaning that Z 1 M 0 f.m/dm D e ˇM 0 (2) Parametrization of this density f varies from catalog to catalog but generally is estimated using expectation maximization techniques. * Correspondence to: Kevin Nichols, McCarthy Hall 154, Department of Mathematics, California State University, 800 N. State College Blvd., Fullerton, CA 92831, U.S.A. knichols@fullerton.edu a b McCarthy Hall 154, Department of Mathematics, California State University, 800 N. State College Blvd., Fullerton, CA 92831, U.S.A. Department of Statistics, University of California, Los Angeles, CA, U.S.A. (2014) Copyright 2014 John Wiley & Sons, Ltd.
2 K. NICHOLS AND F. P. SCHOENBERG A key feature of ETAS models, see Schoenberg (2003), is the assumption of separability between the magnitude of an earthquake and the magnitudes of any aftershocks produced by that earthquake. An extension of spatial ETAS models (1) to incorporate the magnitude of potential earthquakes is, 0 1.t; x; y; M; H t / D y/ C X g.t t i ;x x i ;y y i I M i / A (3) iwt i <t which inherently assumes that the magnitude of an aftershock earthquake is independent of the magnitude of the earthquake that produced it (Schoenberg, 2003). In contrast to the assumptions carried by ETAS models for earthquake occurrence, consider now the accelerated moment release model (AMR), for references see Bowman et al., (1998), Bowman and King (2001), Jaumé and Sykes (1999). AMR theory assumes that the potential energy, or more specifically the rate of seismic moment release in moderate earthquakes, increases prior to the occurrence of larger earthquakes (Jaumé and Sykes, 1999). Specifically, the cumulative Benioff strain.t/, which is a measure of seismic energy or seismic potential, can be fitted by a relation of the form: N.t/ X.t/ D E i.t/ 1 2 id1 (4) where N.t/ represents the number of earthquakes occurring prior to time t within a set spatial window, usually restricted to the range of earthquakes associated with the same fault, of the present location, and E i represents the energy of each previously occurring event (Bowman and King, 2001). Most studies supporting AMR theory focus on the earthquakes that precede the very largest earthquakes in historical record (Bowman et al., 1998; Bowman and King 2001; Jaume and Sykes 1999). The inherent assumption being that the magnitudes of the earthquakes preceding the largest earthquakes are in fact inseparable from each other, that is, there is some degree of dependence or positive correlation between the magnitudes of earthquakes that occur in clusters preceding the largest earthquakes. While this contradicts the assumption of ETAS models, the contradiction frames the debate in recent literature as to whether or not the magnitudes of earthquakes and the magnitudes of their aftershocks are correlated. The research performed by Lippiello et al., (2007), Lippiello et al., (2008), Lippiello et al., (2009) indicates there is a positive correlation between earthquakes and the earthquakes that occur after them. However, the research of (Davidsen and Green, 2011) concludes the hypothesis that earthquakes are correlated in magnitude with the earthquakes that follow them cannot be rejected, and any positive correlation that may have been noticed is mainly a result of short-term aftershock incompleteness (STAI). There may be a simple reason as to why the conclusions on the subject are so varied. To date, work on the subject has focused on correlations between one earthquake and the next earthquake(s) to follow it. This does not address the question of interest, namely are the magnitudes earthquakes correlated with the magnitudes of their aftershocks, but rather addresses the question of whether or not the magnitudes of earthquakes are correlated with the magnitudes of earthquakes that occur after them. Perhaps the reason why research has not focused on interdependency between earthquakes and their aftershocks is because it is difficult to determine which earthquakes are main shocks, which earthquakes are aftershocks, and which previously occurring earthquake triggered an aftershock. There have been several proposed methodologies for separating an earthquake catalog into background or primary events and secondary or aftershock events. Most of these methodologies essentially remove the earthquakes in a small spatial temporal window around larger events and assign remaining events as mainshocks; for examples, see Gardener and Knopoff (1990), and Utsu (1995). The obvious problems with this approach are that main shock earthquakes that occur near large earthquakes are misclassified as aftershocks and that smaller earthquakes are not considered self-exciting. In fact, these methods for stochastically de-clustering earthquake catalogs leads to the common misuse of the term aftershock as any earthquake occurring after a large earthquake and misuse of the term foreshock an any earthquake occurring before a large earthquake. Alternative to these window based methodologies are methods that consider the space-time distance between events and then parametrically determine if an event is a mainshock or an aftershock depending on its proximity in both space and time to preceding earthquakes (e.g. Reasenberg, 1985; Davis and Frohlich, 1991). These types of methods are still used today in determining spatial temporal windows for aftershock occurrence; for example, see how aftershocks are defined in (Lippiello et al., 2009). Unfortunately, these methods require parameterization for defining the shape of the aftershock windows or for determining the significance of the link between two events and ultimately underestimate the probability that smaller earthquakes trigger aftershocks. Marsan and Lengline (2008) developed a model independent stochastic de-clustering (MISD) data algorithm for nonparametrically estimating background rate and triggering function g.t t i ;x x i ;y y i I M i /. Consequently, MISD allows for probabilistic assignment of each event as a mainshock or an aftershock, with individual probabilities assigned to each previously occurring earthquake as the trigger for the aftershock. In this paper, we expand the capability of the algorithm to allow for a spatially dependent background rate.x; y/, hence forth spatially dependent model independent stochastic de-clustering (SDMISD). Model independent stochastic de-clustering is an expectation-maximization algorithm. It should be noted that some work has been carried out on improving the MISD algorithm, most notably through Maximized Penalized Likelihood Estimation to increase the accuracy of the algorithm (Lewis and Mohler, 2011). MISD can readily be applied to global earthquake catalogs. MISD does not allow for the error-free classification of events as either mainshocks or aftershocks. However, repeatedly assigning events as mainshocks or aftershocks, and in the case of aftershocks determining the earthquake that triggered the aftershock, will produce a set of pairings between earthquakes and their aftershocks in which the correct lineage will have higher probability and therefor greater representation than any incorrect pairing. This provides an opportunity to observe the dependency of earthquake magnitudes and their aftershock magnitudes. wileyonlinelibrary.com/journal/environmetrics Copyright 2014 John Wiley & Sons, Ltd. (2014)
3 MAGNITUDE DEPENDENCE 2. DATA Fifteen years of earthquake data, ranging from 1996 to 2010, was obtained from the Centroid Moment Tensor (CMT) global catalog (gcmt). The data is freely available online and can be found in NDK format at NDK format includes the following information of interest: latitudes, longitudes, dates, times, and scalar moments M sc, which can in turn easily be converted into moment magnitudes M using the US Geological Survery (USGS) earthquake magnitude policy conversion formula: M D log 10.M sc / 10:7 (5) 1:5 Figure 1 is a global plot of all earthquakes contained within the temporal bounds set for the catalog. There are several different measures of magnitude when it comes to earthquakes, including the Richter scale, measures of fault geometry and seismic moment, measures of energy released, and measures of wavelength such as surface wave length magnitudes and body wave length magnitudes. However, the reason that moment magnitudes are ideal when compared with other choices is that moment magnitudes do not saturate for and after larger events. The consequences of saturation are that typically other measures of earthquake magnitude tend to underestimate the magnitudes of the larger events. As such, for our measure of earthquake moment, we chose moment magnitudes though it is important to note that the forthcoming methods could easily accommodate a separate choice of measure for earthquake magnitude. Most earthquake catalogs are believed to be incomplete for smaller earthquakes, meaning that small earthquakes are often unnoticed or unrecorded in seismology data. Additionally, following the largest earthquakes, there is a period of STAI as seismic instruments for detecting earthquakes are temporarily saturated (Kagan 2004, Zhuang et al., 2004, Helmstetter et al., 2006, Lippiello et al., 2007, Lippiello et al., 2008, Davidsen and Green 2011). As a result, small earthquakes occurring during this period of saturation go undetected by seismic instruments. To address these issues with incompleteness of the catalog, the minimum magnitude of completeness M c is approximated by computing the low-magnitude breakpoint using median based analysis of the segment slope (MBASS) (Amorese, 2007). MBASS is a bootstrap method that essentially assesses at what point the slope of the log-survivor function is linear (as it should be for a complete realization of an Figure 1. Map with spatial occurrences of CMT global catalog earthquakes occurring between 1996 and 2010 Figure 2. The median based analysis of segment slope using a survivor function. Data appears approximately exponentially distributed for events > 5.3 (2014) Copyright 2014 John Wiley & Sons, Ltd. wileyonlinelibrary.com/journal/environmetrics
4 K. NICHOLS AND F. P. SCHOENBERG exponentially distributed random variable like earthquake magnitudes). The average (and median) low-magnitude breakpoint after 1,000 bootstrap iterations of MBASS for the catalog used in this analysis is M c D 5:3. It is fair to assume that the catalog is relatively complete for moment magnitudes (M > 5:3). Figure 2 is a plot of magnitude against the log frequency of earthquakes greater than that magnitude and helps to visually confirm completeness of the catalog for events with magnitude greater than 5.3 as M c D 5:3 is the lowest magnitude cut-off such that the log-survivor function plot of all data above M c is linear (Wiemer and Wyss, 2000). As incompleteness of events smaller than 5.3 in the catalog may contribute to spurious results any events with magnitude less than 5.3 have been deleted from the catalog in this analysis. Additionally, in consideration of STAI, the relationship between the minimum magnitude of completeness M c.t; m/ and the time t (in days) following an earthquake of magnitude m is M c.t; m/ D m 4:5 0:75log 10.t/ (Helmstetter et al., 2006). In this analysis, to consider earthquakes during periods in which M c is less than the overall minimum magnitude of completeness can lead to spurious results. As such, any earthquakes occurring during a time period in which M c.t; m/ < 5:3 were deleted from the analysis. It should be noted that the impact of removing earthquakes because of STAI was minimal as the minimum magnitude of completeness for the gcmt is relatively high compared to other catalogs. For more information on the gcmt, please see (Ekstrom et al., 2012). The result is a 15-year period in which 11,535.M > 5:3/ earthquakes are used for the forthcoming analyses. 3. METHODS The goal of the proposed methodology is to observe any trends of dependency between the magnitudes of aftershock earthquakes and the earthquakes that trigger them. This first requires stochastic de-clustering of the data. Once the catalog has been de-clustered, conditional averages for aftershock magnitude can be computed by grouping the magnitudes of aftershocks that come from earthquakes of comparable magnitude. While SDMISD allows us to achieve the desired classification of each earthquake as a mainshock or as an aftershock triggered by a previously occurring earthquake, admittedly there will be potential nonsystematic misclassifications. However, by repeatedly applying SDMISD to the same catalog, and assuming that catalog has a reasonably large number of events, on average the SDMISD will correctly compute the lineage between earthquakes and their aftershocks more often than any misclassification. By repeating the process, it is possible to compute the conditional average of aftershock magnitudes conditioned upon the magnitudes of the earthquakes that triggered the aftershock. Additionally, repeating the process allows for an examination on the distribution of conditional means. Using these distributions of conditional aftershock means, confidence intervals are computed for the average aftershock magnitude caused by earthquakes with similar magnitudes. Forthcoming, this section shall be split into three parts. The first is dedicated to an examination of MISD algorithm and how we can implement kernel smoothing to estimate a nonconstant background rate, thus developing SDMISD. In addition to developing SDMISD itself, we implement a methodology for examining results for SDMISD under two competing hypotheses, namely the null hypothesis in which we assume that earthquake magnitudes are independent of each other and the alternative hypothesis in which we assume earthquake magnitudes are dependent on each other. To simulate under the null hypothesis, we can simply re-create our catalog and impose independence by randomly shuffling the magnitudes of the earthquakes in the catalog. To simulate under various versions of the alternative hypothesis will require the creation of synthetic catalogs. The second part of this section is dedicated toward the creation of realistic multidimensional synthetic catalogs. The advantage of illustrating the method using synthetic catalogs before using the global CMT catalog is that we can obtain a descriptive sense of results that are typical for various levels of earthquake magnitude correlation. Finally, in the third part, a methodology is described for repeatedly using SDMISD to create estimates and confidence intervals for the conditional means of aftershock magnitudes depending on the magnitudes of the earthquakes that trigger them MISD and SDMISD The following is a summary of the methodology for probabilistically determining connectivity between earthquakes. For greater detail, please see (Marsan and Lengline, 2008). The methodology begins with a simplified version of (1):.x; y; t/ D 0 C X t i <t i.x;y;t;m/ (6) where.x; y; t/ represents the rate at which earthquakes occur in space time and is the sum of a constant background rate for mainshock occurrence, 0, and the rates for all previously occurring earthquakes to trigger an aftershock, i.x; y; t/. It is important to note the contribution of a previously occurring earthquake at location (x i,y i ), time (t i ), and with magnitude (M i ) to the rate of earthquake occurrence at location.x; y/ and time.t/ is a nonparametric function of distances j.x; y/.x i ;y i /j, t t i and magnitude M i. Model independent stochastic declustering for determining 0 and i.x; y; t/ begins with a set of distances between any pairing of events in the data as well as an aprioriestimate for.j.x; y/j; t;m/and for 0. The algorithm begins by estimating triggering weights w i;j = j.j.x j ;y j /.x i ;y i /j;t j t i ;M i / if t i <t j, w i;j D 0; otherwise, and background weights w 0;j D j 0. j is a normalizing coefficient such that P.w i;j / C w 0;j D 1 for events i occurring prior to event j. The rate contributions of previously occurring earthquakes are then updated as.j.x; y/j; t;m/d 1 N M ıt S.j.x; y/j;ır/ X i;j2a w i;j (7) wileyonlinelibrary.com/journal/environmetrics Copyright 2014 John Wiley & Sons, Ltd. (2014)
5 MAGNITUDE DEPENDENCE where A consists of all pairs such that j.x j ;y j /.x i ;y i /j2j.x; y/j ır, t j t i 2 t ıt and M i 2 M ım (ır;ıt;ım are discretization parameters), N M is the number of earthquakes such that M i 2 M ım,ands.j.x; y/j/ is the surface covered by the disk with radii jı.x; y/j ır. The updated background rate is 0 D 1 NX w 0;j T S j D1 (8) where T and S represent the length of the temporal support and the surface area of the spatial support, respectively. After an initialization of 0 and..x; y/; t; M /, the process of updating weights w 0;j and w i;j and then updating 0 and..x; y/; t; M / continues until convergence. Once MISD has converged the j th event is identified as a mainshock with probability w 0;j P wi;j C w 0;j (9) and as an aftershock triggered by the ith event with probability w i;j P wi;j C w 0;j (10) While MISD is desirable in the sense that it allows for nonparametric definition of the Poissonian (lambda) rate space triggering function g, the requirement for constant background rate is sub-optimal. Typically, to nonparametrically estimate the rate of earthquake occurrence, we would consider kernel smoothing a density function. Literature will show that often it is not the choice of the kernel that matters all that much for kernel smoothing but rather the choice of smoothing parameter. Similar to (Helmstetter et al., 2006), we use an isotropic kernel: K.d/ D C.d/ r 2 C d 2 1:5 (11) where C(d) is normalizing factor and d is the bandwidth determined on a per event basis using cross-validation; see (Helmstetter et al., 2006) for greater detail. Unfortunately, the desired result is a background rate for all earthquakes, when we desire a background rate for mainshock earthquakes. To account for this, we consider the following adaptation to (8). Let 0;j be the background rate at the j th catalog event, K j be the estimated kernel density at the location of the j th catalog event, and KN be the average kernel density across the locations of all events. We substitute for (8) the following: K j NX 0;j D w T S KN 0;j j D1 (12) and then update background rate w 0;j D j 0;j. Note that the mean of 0;j from (12) is equal to 0 from (8) Synthetic catalogs Multiple dimensional earthquake catalogs were simulated under various levels of inter-dependency between primary events and secondary events. Catalogs were simulated by first creating a spatially dependent background rate and randomly generating primary events over a uniform temporal domain. This was performed by kernel smoothing points tracing line segments that represent synthetic fault lines, then sampling coordinates from the kernel support. Magnitudes for the primary events were simulated in accordance with Gutenberg Richter laws by using inversion techniques with an exponential distribution in which the minimum magnitude of completeness set to 5.3, which is the same as our global catalog. The parameter of the exponential distribution used to sample primary event magnitudes was selected so that the expectation of earthquake magnitude in our synthetic catalog was equal to the average earthquake in our global catalog. After generating synthetic primary events, first generation secondary events were generated. The expected number of secondary events (N a ft.m/) over the temporal domain.t / for primary event with magnitude.m/, following a brief period of STAI.c/ was determined in accordance with generalized Omori law (see Helmstetter et al., 2005 for greater detail): N af t.m/ D K m ˇm d / T.1 p/ c.1 p/ 1 p (13) where p D 0:9, K 0 = , D ˇ D 1. Using this expectation as the parameter of a Poisson distribution, the number of secondary events associated with each primary event was generated. The temporal coordinates of secondary events were determined using generalized (2014) Copyright 2014 John Wiley & Sons, Ltd. wileyonlinelibrary.com/journal/environmetrics
6 K. NICHOLS AND F. P. SCHOENBERG Omori Law by sampling from (13) using the inversion method. The spatial coordinates of secondary events were generated by randomly sampling a uniform direction (), and a distance.d/ from an adapted version of (13): f.d/d K 010. m ˇm d /.d/ p (14) The magnitudes of secondary events triggered M af t.m i / triggered by event m i were generated using the following equation: M af t.m i / D.1 a/m r C am i (15) where.m r / is a new randomly generated primary event. As (a) increases, so does interdependency between primary events and secondary events. After first generation aftershocks are generated from mainshocks, second generation aftershocks are generated from first generation aftershocks and so forth until randomly no further aftershocks are generated Estimates and confidence intervals Repeated probabilistic assignment of earthquakes as either mainshocks or aftershocks creates a large set of earthquake pairings for us to examine. From these pairings, and for each probabilistic assignment of all the earthquakes in a catalog, we can easily compute the average magnitude of aftershocks caused by earthquakes of equal magnitude. After 1,000 such assignments, we can create frequency based confidence intervals for the average aftershock resulting from a mainshock of size.m /. After 1000 of such classifications, the 5th and 95th percentiles for the average aftershock caused by an earthquake of magnitude M serve as realistic 90% confidence bounds on the expected aftershock caused by an earthquake of magnitude M (Efron and Tibshirani, 1986). The conditional expected aftershock caused by an earthquake of size M can then be used inversely to estimate the parameter of a conditional exponential distribution f for earthquake aftershocks caused by an earthquake of size M. Average Aftershock (90% confidence) Average Aftershock (90% confidence) Magnitude of Parent Shock Magnitude of Parent Shock Average Aftershock (90% confidence) Average Aftershock (90% confidence) Magnitude of Parent Shock Magnitude of Parent Shock Figure 3. Simulation based 90% confidence intervals for the average aftershock magnitude (y-axis) corresponding to the magnitude the previously occurring earthquake that excited the aftershock (x-axis). Results were simulated for self-exciting earthquake and aftershock magnitudes are independent (Figure 3(a), top left, 0:0), weakly positively correlated (Figure 3(b), top right, 0:3), positively correlated (Figure 3(c), bottom left, 0:55), and strongly positively correlated (Figure 3(d), bottom right, 0:75) wileyonlinelibrary.com/journal/environmetrics Copyright 2014 John Wiley & Sons, Ltd. (2014)
7 MAGNITUDE DEPENDENCE 4. RESULTS Figure 3 is a plot of the conditional means with confidence intervals of aftershock magnitude as a function of the magnitudes of the earthquakes that come from for the four synthetic catalogs. Important features for Figure 3 include the increase in variability for the conditional averages as the magnitudes of the triggering earthquakes increase. Figure 4 is a display of the average magnitude of aftershocks with confidence intervals conditioned upon the magnitude of the earthquakes they come from. Figure 5 is the same display, only this time the described analyses were performed after randomly shuffling the magnitudes of the earthquakes. Contrasting Figure 4 both with the results of the synthetic catalogs and with the results from imposing independence on the magnitudes of earthquakes and their aftershocks through random shuffling yields the following observations. Table 1 contains the estimated 90% confidence bounds for various earthquake magnitude ranges. Table 1 indicated there is very strong evidence that the magnitudes of earthquakes are correlated with the magnitudes of their aftershocks. The correlation between the magnitudes of earthquakes and their aftershocks is moderate to strong (compare with results from synthetic catalogs in which the correlation was 0.55 or 0.75). It may appear that positive trend does not seem to hold for the largest earthquakes (M > 7). Note, however, that there were very few earthquakes with such magnitudes and with so little data SDMISD struggles to define the triggering function. The result is severe variation in the confidence interval estimates. These are significant results that indicate that the local history of earthquakes are indeed powerful components in modeling not only the rate of future earthquake occurrences but also the magnitudes of future earthquakes. Figure 4. Ninety percent confidence intervals for the average aftershock magnitude (y-axis) associated with varying foreshock magnitudes (M) such that M > 5:3 (x-axis). Results are based on 1000 applications of model independent stochastic declustering (MISD) to the CMT catalog. Horizontal line represents the average magnitude of all aftershocks Figure 5. Ninety percent confidence intervals for the average aftershock magnitude (y-axis) associated with varying foreshock magnitudes (M) such that M > 5:3 (x-axis). Results are based on 1000 applications of model independent stochastic declustering (MISD) to the CMT catalog with randomly shuffled magnitudes. Horizontal line represents the average magnitude of all aftershocks (2014) Copyright 2014 John Wiley & Sons, Ltd. wileyonlinelibrary.com/journal/environmetrics
8 K. NICHOLS AND F. P. SCHOENBERG Table 1. Ninety percent confidence interval bounds for the average aftershock magnitude associated with varying foreshock magnitudes (M) such that M > 5:3 M Lower bound Upper bound M Lower bound Upper bound DISCUSSION The methodology described in this article has yielded results that potentially could influence the ways in which popular models like ETAS model earthquake rates. An immediate impact could be the inclusion of a conditional term for the distribution of earthquake magnitude depending on whether or not the earthquake is a mainshock or an aftershock. Updating (3) to reflect the findings of this article could potentially yield a new form of the ETAS model in which.t; x; y; M; H t / D f 1.M /.x; y/ C X iwt i <t f 2.M jm i / g.t t i ;x x i ;y y i I M i / (16) where f 1.M / and f 2.M jm i / represent the distribution of mainshock earthquake magnitudes and the distribution of aftershock earthquake magnitudes conditioned upon their triggering earthquake magnitude, respectively. Note that both f 1 and f 2 could easily be nonparametrically estimated using the methodology prescribed in this article and a reasonably sized catalog. REFERENCES Amorese D Applying a change-point detection method on frequency-magnitude distributions. Bulletin of the Seismological Society of America 97: , DOI: / Bowman D, King G Accelarating seismicity and stress accumulation before large earthquakes. Geophysical Research Letters 28(21): Bowman D, Ouillon G, Sammis CG, Sornette D, Sornette A An observational test of the critical earthquake concept. Journal of Geophysics 103(24): Davidsen J, Green A Are earthquake magnitudes clustered? Physical Review Letters 106: , DOI: /PhysRevLett Davis SD, Frohlich C Single-link cluster analysis, synthetic earthquake catalogs, and aftershock identification. Geophysical Journal International 104: Efron B, Tibshirani R Bootstrap methods for standard errors, confidence intervals, and other measures of statistical accuracy. Statistical Science 1(1): Ekstrom G, Nettles M, Dziewonski A The global CMT project : centroid-moment tensors for 13,017 earthquakes. Physics of the Earth and Planetary Interiors 200:1 9, DOI: /j.pepi Gardner J, Knopoff L Is the sequence of earthquakes in southern California, with aftershock removed, Poissonian? Bulletin of the Seismological Society of America 64: Gutenberg B, Richter C Frequency of earthquakes in California. Bulletin of the Seismological Society of America 34: Hawkes AG Spectra of some self-exciting and mutually exciting point processes. Biometrika 58: Helmstetter A, Kagan Y, Jackson D Importance of small earthquakes for stress transfers and earthquake triggering. Journal of Geophysical Research 110(B05S08):1 13, DOI: /2004JB Helmstetter A, Kagan Y, Jackson D Comparison of short-term and time-independent earthquake forecast models for southern California. Bulletin of the Seismological Society of America 96: wileyonlinelibrary.com/journal/environmetrics Copyright 2014 John Wiley & Sons, Ltd. (2014)
9 MAGNITUDE DEPENDENCE Jaumé S, Sykes L Evolving towards a critical point: a review of accelerating moment/energy release prior to large and great earthquakes. Pure and Applied Geophysics 155: Kagan Y Short term properties of earthquake magnitudes and models of earthquake source. Bulletin of the Seismological Society of America 94(4): Lewis E, Mohler G A nonparametric EM algorithm for a multiscale Hawkes process. Journal of Nonparametric Statistics. inreview. Lippiello E, Arcangelis L, Godano C Influence of time and space correlations on earthquake magnitude. Physical Review Letters 100: Lippiello E, Arcangelis L, Godano C Role of static stress diffusion in the spatiotemporal organization of aftershocks. Physical Review Letters 103: Lippiello E, Godano C, Arcangelis L Dynamical scaling in branching models for seismicity. Physical Review Letters 98: Marsan D, Lengline O Extending earthquake s reach through cascading. Science 319: Ogata Y Statistical models for earthquake occurrences and residual analysis for point processes. JASA 83:9 27. Ogata Y Space-time point-process models for earthquake occurrences. Annals of the Institute of Statistical Mathematics 50: Reasenberg P Second-order moment of central California seismicity. Journal of Geophysical Research 90(B7): Schoenberg FP Multi-dimensional residual analysis of point process models for earthquake occurrences. JASA 98(464): Utsu T, Ogata Y, Matsu ura RS The centenary of the Omori formula for a decay law of aftershock activity. Journal of Physics of the Earth 43:1 33. Wiemer S, Wyss M Minimum magnitude of completeness in earthquake catalogs: examples from Alaska, the Western United States, and Japan. Bulletin of the Seismological Society of America 90(4): Zhuang J, Ogata Y, Vere-Jones D Analyzing earthquake clustering features by using stochastic reconstruction. Journal of Geophysical Research 109(B05301):1 17, DOI: /2003JB (2014) Copyright 2014 John Wiley & Sons, Ltd. wileyonlinelibrary.com/journal/environmetrics
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 informationImpact 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 informationStatistical Properties of Marsan-Lengliné Estimates of Triggering Functions for Space-time Marked Point Processes
Statistical Properties of Marsan-Lengliné Estimates of Triggering Functions for Space-time Marked Point Processes Eric W. Fox, Ph.D. Department of Statistics UCLA June 15, 2015 Hawkes-type Point Process
More informationarxiv: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 informationA 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 informationThe 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 informationComparison 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 informationQuantifying 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 informationComparison of short-term and long-term earthquake forecast models for southern California
Comparison of short-term and long-term earthquake forecast models for southern California A. Helmstetter, Yan Kagan, David Jackson To cite this version: A. Helmstetter, Yan Kagan, David Jackson. Comparison
More informationarxiv: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 informationMagnitude Of Earthquakes Controls The Size Distribution Of Their. Triggered Events
A c c e p t e d f o r p u b l i c a t i o n i n J G R : S o l i d E a r t h 1 Magnitude Of Earthquakes Controls The Size Distribution Of Their Triggered Events Shyam Nandan 1, Guy Ouillon 2, Didier Sornette
More informationShort-Term Properties of Earthquake Catalogs and Models of Earthquake Source
Bulletin of the Seismological Society of America, Vol. 94, No. 4, pp. 1207 1228, August 2004 Short-Term Properties of Earthquake Catalogs and Models of Earthquake Source by Yan Y. Kagan Abstract I review
More informationAre Declustered Earthquake Catalogs Poisson?
Are Declustered Earthquake Catalogs Poisson? Philip B. Stark Department of Statistics, UC Berkeley Brad Luen Department of Mathematics, Reed College 14 October 2010 Department of Statistics, Penn State
More informationTheory 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 informationMulti-dimensional residual analysis of point process models for earthquake. occurrences. Frederic Paik Schoenberg
Multi-dimensional residual analysis of point process models for earthquake occurrences. Frederic Paik Schoenberg Department of Statistics, University of California, Los Angeles, CA 90095 1554, USA. phone:
More informationAftershock 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 informationJournal 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 informationTesting for Poisson Behavior
Testing for Poisson Behavior Philip B. Stark Department of Statistics, UC Berkeley joint with Brad Luen 17 April 2012 Seismological Society of America Annual Meeting San Diego, CA Quake Physics versus
More informationMagnitude uncertainties impact seismic rate estimates, forecasts, and predictability experiments
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 113,, doi:10.1029/2007jb005427, 2008 Magnitude uncertainties impact seismic rate estimates, forecasts, and predictability experiments Maximilian J. Werner 1 and Didier
More informationResearch Article. J. Molyneux*, J. S. Gordon, F. P. Schoenberg
Assessing the predictive accuracy of earthquake strike angle estimates using non-parametric Hawkes processes Research Article J. Molyneux*, J. S. Gordon, F. P. Schoenberg Department of Statistics, University
More informationSelf-exciting point process modeling of crime
Self-exciting point process modeling of crime G. O. Mohler M. B. Short P. J. Brantingham F. P. Schoenberg G. E. Tita Abstract Highly clustered event sequences are observed in certain types of crime data,
More informationAppendix 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 informationDistribution 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 informationSecondary Aftershocks and Their Importance for Aftershock Forecasting
Bulletin of the Seismological Society of America, Vol. 93, No. 4, pp. 1433 1448, August 23 Secondary Aftershocks and Their Importance for Aftershock Forecasting by Karen R. Felzer, Rachel E. Abercrombie,
More informationSouthern California Earthquake Center Collaboratory for the Study of Earthquake Predictability (CSEP) Thomas H. Jordan
Southern California Earthquake Center Collaboratory for the Study of Earthquake Predictability (CSEP) Thomas H. Jordan SCEC Director & Professor, University of Southern California 5th Joint Meeting of
More informationProximity to Past Earthquakes as a Least-Astonishing Hypothesis for Forecasting Locations of Future Earthquakes
Bulletin of the Seismological Society of America, Vol. 101, No. 4, pp. 1618 1629, August 2011, doi: 10.1785/0120090164 Proximity to Past Earthquakes as a Least-Astonishing Hypothesis for Forecasting Locations
More informationEarthquake Clustering and Declustering
Earthquake Clustering and Declustering Philip B. Stark Department of Statistics, UC Berkeley joint with (separately) Peter Shearer, SIO/IGPP, UCSD Brad Luen 4 October 2011 Institut de Physique du Globe
More informationThe 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 informationPoint processes, spatial temporal
Point processes, spatial temporal A spatial temporal point process (also called space time or spatio-temporal point process) is a random collection of points, where each point represents the time and location
More informationAdaptive 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 informationUNIVERSITY 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 informationSelf-similar earthquake triggering, Båth s law, and foreshock/aftershock magnitudes: Simulations, theory, and results for southern California
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 117,, doi:10.1029/2011jb008957, 2012 Self-similar earthquake triggering, Båth s law, and foreshock/aftershock magnitudes: Simulations, theory, and results for southern
More informationEarthquake Likelihood Model Testing
Earthquake Likelihood Model Testing D. Schorlemmer, M. Gerstenberger 2, S. Wiemer, and D. Jackson 3 Swiss Seismological Service, ETH Zürich, Schafmattstr. 30, 8093 Zürich, Switzerland. 2 United States
More informationEarthquake Likelihood Model Testing
Earthquake Likelihood Model Testing D. Schorlemmer, M. Gerstenberger, S. Wiemer, and D. Jackson October 8, 2004 1 Abstract The Regional Earthquake Likelihood Models (RELM) project aims to produce and evaluate
More informationA hidden Markov model for earthquake declustering
Click Here for Full Article JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 115,, doi:10.1029/2008jb005997, 2010 A hidden Markov model for earthquake declustering Zhengxiao Wu 1 Received 8 August 2009; revised 18
More informationSupplementary Materials for. Time-advanced occurrence of moderate-size. earthquakes in a stable intraplate region after a
Supplementary Materials for Time-advanced occurrence of moderate-size earthquakes in a stale intraplate region after a megathrust earthquake and their seismic properties Tae-Kyung Hong a,, Junhyung Lee
More informationHigh-resolution Time-independent Grid-based Forecast for M >= 5 Earthquakes in California
High-resolution Time-independent Grid-based Forecast for M >= 5 Earthquakes in California A. Helmstetter, Yan Kagan, David Jackson To cite this version: A. Helmstetter, Yan Kagan, David Jackson. High-resolution
More informationOn Mainshock Focal Mechanisms and the Spatial Distribution of Aftershocks
On Mainshock Focal Mechanisms and the Spatial Distribution of Aftershocks Ka Wong 1 and Frederic Paik Schoenberg 1,2 1 UCLA Department of Statistics 8125 Math-Science Building, Los Angeles, CA 90095 1554,
More informationLimitations of Earthquake Triggering Models*
Limitations of Earthquake Triggering Models* Peter Shearer IGPP/SIO/U.C. San Diego September 16, 2009 Earthquake Research Institute * in Southern California Why do earthquakes cluster in time and space?
More informationALM: 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 informationProbabilistic 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 informationSelf-exciting point process modeling of crime
Self-exciting point process modeling of crime G. O. Mohler M. B. Short P. J. Brantingham F. P. Schoenberg G. E. Tita Abstract Highly clustered event sequences are observed in certain types of crime data,
More informationarxiv:physics/ v1 [physics.geo-ph] 19 Jan 2005
APS preprint Distribution of the Largest Aftershocks in Branching Models of arxiv:physics/0501102v1 [physics.geo-ph] 19 Jan 2005 Triggered Seismicity: Theory of Båth s law A. Saichev 1, 2 and D. Sornette
More informationSimulated 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 informationFIRST PAGE PROOFS. Point processes, spatial-temporal. Characterizations. vap020
Q1 Q2 Point processes, spatial-temporal A spatial temporal point process (also called space time or spatio-temporal point process) is a random collection of points, where each point represents the time
More informationGeophysical 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 informationAccuracy of modern global earthquake catalogs
Physics of the Earth and Planetary Interiors 135 (2003) 173 209 Accuracy of modern global earthquake catalogs Yan Y. Kagan Department of Earth and Space Sciences, University of California, Los Angeles,
More informationTime-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 informationEarthquake prediction: Simple methods for complex phenomena. Bradley Luen
Earthquake prediction: Simple methods for complex phenomena by Bradley Luen A dissertation submitted in partial satisfaction of the requirements for the degree of Doctor of Philosophy in Statistics in
More informationEdinburgh Research Explorer
Edinburgh Research Explorer Statistical Modeling of the 1997-1998 Earthquake Sequence Citation for published version: Touati, S, Naylor, M & Main, IG 2014, 'Statistical Modeling of the 1997-1998 Earthquake
More informationModeling Aftershocks as a Stretched Exponential Relaxation
Modeling Aftershocks as a Stretched Exponential Relaxation A. Mignan Institute of Geophysics, Swiss Federal Institute of Technology Zurich, Switzerland arnaud.mignan@sed.ethz.ch Abstract: The decay rate
More informationAN OVERVIEW AND GUIDELINES FOR PROBABILISTIC SEISMIC HAZARD MAPPING
CO 2 TRACCS INTERNATIONAL WORKSHOP Bucharest, 2 September, 2012 AN OVERVIEW AND GUIDELINES FOR PROBABILISTIC SEISMIC HAZARD MAPPING M. Semih YÜCEMEN Department of Civil Engineering and Earthquake Studies
More informationComment 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 informationAnalyzing earthquake clustering features by using stochastic reconstruction
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 109,, doi:10.1029/2003jb002879, 2004 Analyzing earthquake clustering features by using stochastic reconstruction Jiancang Zhuang and Yosihiko Ogata Institute of Statistical
More informationSpace-time clustering of seismicity in California and the distance dependence of earthquake triggering
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 117,, doi:10.1029/2012jb009471, 2012 Space-time clustering of seismicity in California and the distance dependence of earthquake triggering P. M. Shearer 1 Received
More informationSpace-time ETAS models and an improved extension
Space-time ETAS models and an improved extension Yosihiko Ogata and Jiancang Zhuang The Institute of Statistical Mathematics, Minami-Azabu 4-6-7, Minato-Ku, Tokyo, 106-8569 Japan, ogata@ism.ac.p Abstract.
More informationPerformance 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 informationEarthquake clusters in southern California I: Identification and stability
JOURNAL OF GEOPHYSICAL RESEARCH: SOLID EARTH, VOL. 118, 2847 2864, doi:10.1002/jgrb.50179, 2013 Earthquake clusters in southern California I: Identification and stability Ilya Zaliapin 1 and Yehuda Ben-Zion
More informationTheme V - Models and Techniques for Analyzing Seismicity
Theme V - Models and Techniques for Analyzing Seismicity Stochastic simulation of earthquake catalogs Jiancang Zhuang 1 Sarah Touati 2 1. Institute of Statistical Mathematics 2. School of GeoSciences,
More informationThe 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 informationlog (N) 2.9<M< <M< <M< <M<4.9 tot in bin [N] = Mid Point M log (N) =
Solution Set for Assignment Exercise : Gutenberg-Richter relationship: log() = a + b. M A) For a time period between January, 90 to December 3, 998 tot in bin [] = 450 6 57 22 7 5 Mid Point M 3.5 3.65
More informationDeclustering and Poisson Tests
Geophys. J. Int. (2012)???, 1 11 Declustering and Poisson Tests Brad Luen Department of Mathematics Reed College Portland, OR, USA Philip B. Stark Department of Statistics University of California Berkeley,
More informationFrom 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 informationEarthquake 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 informationBecause 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 informationSUPPLEMENTARY INFORMATION
SUPPLEMENTARY INFORMATION DOI: 1.138/NGEO177 The Long Precursory Phase of Most Large Interplate Earthquakes Supplementary Information Supplementary Methods 1. Data and Classification We select the earthquakes
More information(1) Istituto Nazionale di Geofisica e Vulcanologia, Roma, Italy. Abstract
Stacked analysis of earthquake sequences: statistical space-time definition of clustering and Omori law behavior Tosi P.(1), De Rubeis V.(1) and Sbarra P.(1) (1) Istituto Nazionale di Geofisica e Vulcanologia,
More informationScaling 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 informationTesting the stress shadow hypothesis
Testing the stress shadow hypothesis Karen Felzer Emily Brodsky March 14, 2005 1 Abstract A fundamental question in earthquake physics is whether aftershocks are predominantly triggered by static stress
More informationUCLA UCLA Electronic Theses and Dissertations
UCLA UCLA Electronic Theses and Dissertations Title Non Parametric Estimation of Inhibition for Point Process Data Permalink https://escholarship.org/uc/item/7p18q7d1 Author Beyor, Alexa Lake Publication
More informationTesting the stress shadow hypothesis
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 110,, doi:10.1029/2004jb003277, 2005 Testing the stress shadow hypothesis Karen R. Felzer and Emily E. Brodsky Department of Earth and Space Sciences, University of
More informationLimits of declustering methods for disentangling exogenous from endogenous events in time series with foreshocks, main shocks, and aftershocks
PHYSICAL REVIEW E 79, 061110 2009 Limits of declustering methods for disentangling exogenous from endogenous events in time series with foreshocks, main shocks, and aftershocks D. Sornette 1,2,3,4, * and
More informationChien-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 informationSmall-world structure of earthquake network
Small-world structure of earthquake network 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,
More informationStephen Hernandez, Emily E. Brodsky, and Nicholas J. van der Elst,
1 The Magnitude-Frequency Distribution of Triggered Earthquakes 2 Stephen Hernandez, Emily E. Brodsky, and Nicholas J. van der Elst, 3 Department of Earth and Planetary Sciences, University of California,
More informationSecond Annual Meeting
SHale gas Exploration l and Exploitation induced Risks Seismic hazard assessment considering fluid-induced seismicity: Inter-event time distribution of seismicity induced by hydraulic fracturing operations
More informationEarthquake Patterns in Diverse Tectonic Zones of the Globe
Pure Appl. Geophys. 167 (2010), 721 741 Ó 2010 The Author(s) This article is published with open access at Springerlink.com DOI 10.1007/s00024-010-0075-3 Pure and Applied Geophysics Earthquake Patterns
More informationRelationship between accelerating seismicity and quiescence, two precursors to large earthquakes
GEOPHYSICAL RESEARCH LETTERS, VOL. 35, L15306, doi:10.1029/2008gl035024, 2008 Relationship between accelerating seismicity and quiescence, two precursors to large earthquakes Arnaud Mignan 1 and Rita Di
More informationPostScript 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 informationSeismic 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 informationPostScript 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 informationMinimum 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 informationDISCLAIMER BIBLIOGRAPHIC REFERENCE
DISCLAIMER This report has been prepared by the Institute of Geological and Nuclear Sciences Limited (GNS Science) exclusively for and under contract to the Earthquake Commission. Unless otherwise agreed
More informationConnecting near and farfield earthquake triggering to dynamic strain. Nicholas J. van der Elst. Emily E. Brodsky
Connecting near and farfield earthquake triggering to dynamic strain Nicholas J. van der Elst Emily E. Brodsky Department of Earth and Planetary Science, University of California, Santa Cruz Abstract:
More informationPersistent link:
Analysis of aftershock and foreshock activity in stable continental regions: Implications for aftershock forecasting and the Hazard of strong earthquakes Author: J. E. Ebel Persistent link: http://hdl.handle.net/2345/bc-ir:107015
More informationEarthquake catalogues and preparation of input data for PSHA science or art?
Earthquake catalogues and preparation of input data for PSHA science or art? Marijan Herak Department of Geophysics, Faculty of Science University of Zagreb, Zagreb, Croatia e-mail: herak@irb.hr EARTHQUAKE
More informationQuantifying 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 informationSome tricks for estimating space-time point process models. Rick Paik Schoenberg, UCLA
Some tricks for estimating space-time point process models. Rick Paik Schoenberg, UCLA 1. Nonparametric estimation. 2. Parametric estimation of clustering. 3. Separable estimation for covariates. 1 1.
More informationA 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 informationOceanic Transform Fault Seismicity Earthquakes of a Different Kind
University of New Hampshire University of New Hampshire Scholars' Repository Center for Coastal and Ocean Mapping Center for Coastal and Ocean Mapping 4-2009 Oceanic Transform Fault Seismicity Earthquakes
More informationNon commercial use only
Time dependent seismicity along the western coast of Canada Evangelos V. Christou, George Karakaisis, Emmanuel Scordilis Department of Geophysics, School of Geology, Aristotle University of Thessaloniki,
More informationA Bayesian Change Point Model To Detect Changes In Event Occurrence Rates, With Application To Induced Seismicity
Vancouver, Canada, July -5, 5 A Bayesian Change Point Model To Detect Changes In Event Occurrence Rates, With Application To Induced Seismicity Abhineet Gupta Graduate Student, Dept. of Civil and Environmental
More informationETH 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 informationMonte Carlo simulations for analysis and prediction of nonstationary magnitude-frequency distributions in probabilistic seismic hazard analysis
Monte Carlo simulations for analysis and prediction of nonstationary magnitude-frequency distributions in probabilistic seismic hazard analysis Mauricio Reyes Canales and Mirko van der Baan Dept. of Physics,
More informationA 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 informationAccelerating 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 informationTesting long-term earthquake forecasts: likelihood methods and error diagrams
Testing long-term earthquake forecasts: likelihood methods and error diagrams Yan Y. Kagan 1 1 Department of Earth and Space Sciences University of California, Los Angeles, California, USA June 17, 2008
More informationSeismic gaps and earthquakes
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 108, NO. B10, 2471, doi:10.1029/2002jb002334, 2003 Seismic gaps and earthquakes Yufang Rong, 1 David D. Jackson, and Yan Y. Kagan Department of Earth and Space Sciences,
More informationOriginally published as:
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,
More informationOn the Role of Multiple Interactions in Remote Aftershock Triggering: The Landers and the Hector Mine Case Studies
Bulletin of the Seismological Society of America, Vol. 96, No. 1, pp. 80 89, February 2006, doi: 10.1785/0120050029 On the Role of Multiple Interactions in Remote Aftershock Triggering: The Landers and
More information