A Bayesian Change Point Model To Detect Changes In Event Occurrence Rates, With Application To Induced Seismicity
|
|
- Elvin Hart
- 6 years ago
- Views:
Transcription
1 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 Engineering, Stanford University, Stanford, CA 9435, USA Jack W. Baker Associate Professor, Dept. of Civil and Environmental Engineering, Stanford University, Stanford, CA 9435, USA ABSTRACT: A significant increase in earthquake occurrence rates has been observed in recent years in parts of Central and Eastern US. There is a possibility that this increased seismicity is anthropogenic and is referred to as induced seismicity. In this paper, a Bayesian change point model is implemented to evaluate whether temporal features of observed earthquakes support the hypothesis that a change in seismicity rates has occurred. This model is then used to estimate when the change is likely to have occurred. The magnitude of change is also quantified by estimating the distributions of seismicity rates before and after the change. These calculations are validated using a simulated data set with a known change point and event occurrence rates; and then applied to earthquake occurrence data for a site in Oklahoma.. INTRODUCTION The level of seismicity in the Central and Eastern US (CEUS) has increased markedly since approximately 9 (Ellsworth, 3). For Oklahoma, the cumulative number of earthquakes with magnitude 3 since 974 is shown in Fig.. The figure shows a marked increase in seismicity rate in Oklahoma after 8. The magnitude 3 threshold was chosen since Coppersmith et al. () described that there is catalog completeness in CEUS at this magnitude level. (All earthquake data from Oklahoma used in this paper has been obtained from catalogs developed and maintained by Oklahoma Geological Survey s Leonard Geophysical Observatory.) There is a possibility that this increased seismicity is a result of underground wastewater injection (e.g., Ellsworth 3; Keranen et al. 3, 4). Seismicity generated as a result of human activities is referred to as induced or triggered seismicity. Induced seismicity in a region can alter its seismic hazard and risk. One of the important components in calculating seismic hazard through Probabilistic Seismic Hazard Analysis (PSHA) is the activity Number of earthquakes Date Figure : Cumulative number of earthquakes in Oklahoma with magnitude 3 from 974 to Sept 4 rate at a seismic source. For a given seismic source, we thus need to establish whether earthquake occurrence data indicates that a change in rates has occurred, or whether the earthquake activity is consistent with normal features of a process having a constant rate.
2 Vancouver, Canada, July -5, 5 In this study, we describe a Bayesian changepoint model that uses event occurrence data to indicate whether a change in event rates occurred. We assume that the event occurrences belong to a Poisson distribution. If a change in rates is detected, we also obtain the probability of change occurring at any given time, and the probability distributions of rates before and after the change. Along with the description of the model, an algorithmic implementation is provided. The model is validated through its application to a simulated data set with known properties. The Bayesian model is then implemented on a region in Oklahoma to evaluate whether change in rates has occurred in this region. Finally, the distributions for seismicity rates before and after the change are calculated. This detection could be used to inform seismic hazard analysis and can serve as a decision support tool for operations that may be linked to induced seismicity.. BAYESIAN MODEL FOR CHANGE POINT DE- TECTION Change point models are used to detect changes in occurrence rates of events. A Bayesian model for change point detection is implemented here to quantify changes in seismicity rates. The unknowns in our problem are the date of change, event rate before the change and rate after the change. Unlike parameter estimation techniques where a single value of a parameter is obtained, using a Bayesian model yields a probability distribution for the parameters. The probability distribution for rates before and after the change is helpful in seismic hazard calculation as it allows a more rigorous accounting of uncertainties. For our model, it is assumed that occurrence of earthquakes is a Poisson process. Declustering an earthquake catalog (i.e. removing dependent aftershocks and foreshocks and only preserving independent mainshocks) accomplishes the requirements of this assumption as described by Gardner and Knopoff (974) and van Stiphout et al. (). Furthermore, it is assumed that the seismicity rates before and after the change, as well as the date of change are mutually independent. Using Bayesian analysis to detect change point in a Poisson process has been described by Raftery and Akman (986)... Marginal posterior distributions of time of change and event rates We define data in an observation period [, T] as a vector of inter-event times (i.e. the time between successive events) t. The first event occurs at time and the n + th event occurs at time T. t = {t,t,...,t n } s.t. t i = T () i The variables τ, λ, and λ define the date of change, event occurrence rate before the change, and occurrence rate after the change, respectively. Since the events belong to a Poisson process, the inter-event times are exponentially distributed. such that f T λ(s) (t) = λ(s)e λ(s)t () λ(s) = { λ, λ, s τ τ < s T (3) A Gamma distribution with parameters k j and θ j is used to define the conjugate prior of λ j. π(λ j ) λ k j j e λ j/θ j (4) The prior distribution for the time of change, τ is assumed to be uniformly distributed over the observation period. This implies that change is equally likely to occur at any time during the observation period. π(τ) = T, τ T (5) The likelihood function for the unknown parameters can be written as L (τ,λ,λ t) t τ t n = λ e λ t λ e λ t t=t t=t τ+ = λ N(τ) e λ τ. N(T ) N(τ) λ e λ (T τ) (6) where N(t) represents the number of events between [, t]. Using the fact that all parameters are
3 Vancouver, Canada, July -5, 5 mutually independent, the posterior density is calculated as π(τ,λ,λ t) L (τ,λ,λ t)π(λ )π(λ )π(τ) (7) Then the marginal distributions for each of τ, λ and λ can be obtained by integrating the posterior density over the remaining two variables. The marginal posterior distribution of τ is calculated as π(τ t) where = π(τ).. π(λ,λ,τ t)dλ dλ λ N(τ)+k e λ ( λ N(T ) N(τ)+k (.e λ τ+ θ ) dλ T τ+ θ ) dλ = T. Γ(r (τ))γ(r (τ)) S (τ) r (τ) S (τ) r (τ) (8) r (τ) = N(τ) + k S (τ) = τ + θ (9) r (τ) = N(T ) N(τ) + k S (τ) = T τ + θ The marginal posterior distribution of λ is calculated as shown below. A closed form solution for this double integration does not exist. Hence, to evaluate the probability distribution, the time range is discretized on a per day basis, and summed up to approximate the marginal distribution. π(λ t) = T T π(λ,λ,τ t)dλ dτ λ r (τ) e λ S (τ) dλ.π(τ)λ r (τ) e λ S (τ) dτ T [ τ= T.λ r (τ) e λ S (τ). Γ(r (τ))s (τ) r (τ) ] () Similarly, the marginal posterior distribution of λ can be calculated as [ π(λ t) T.λ r (τ) e λ S (τ) T τ=. Γ(r (τ))s (τ) r (τ) ] () Another quantity of interest is the ratio of prechange event rate to post-change event rate (or viceversa), defined as β = λ /λ. Lindley (965) describes the following function of β conditional on τ to follow the F-distribution with d.o.f. p and p, where p k = r k (τ). S (τ)r (τ) S (τ)r (τ) β F p,p () The above equation can be used to compute the probability distribution of β conditional on τ. ( ) r (τ) r (τ) p(β τ, t) = B(r (τ),r (τ)) r (τ) ( ) S (τ) + S (τ)β (r (τ)+r (τ)). S (τ) ( ) S (τ)r (τ) r (τ). S (τ)r (τ) β. S (τ)r (τ) S (τ)r (τ) (3) The marginal distribution of β is then calculated as p(β t) = T p(β τ, t)π(τ t)dτ T π(τ)β r (τ).(s (τ) + S (τ)β) (r (τ)+r (τ)) dτ T τ= T.β r (τ).(s (τ) + S (τ)β) (r (τ)+r (τ)) (4) It is noted after the evaluation of eq. 4 that there is some difference in the equation terms compared to those obtained by Raftery and Akman (986). We are able to verify and validate our equation with the same data as Raftery and Akman (986), and believe this formulation to be correct. 3
4 Vancouver, Canada, July -5, 5.. Bayes factor The change-point model described in the previous section does not indicate whether the data supports the presence of a change in the observed time range. When the model is applied to observed data, it assumes that there is a change and calculates the probability of change on any given date. To indicate whether data favors a change-point model, the Bayes factor is used. The Bayes factor is a Bayesian alternative to hypothesis testing. In this case, it is defined as the ratio of the likelihood function for a constant rate model (H ) to that for a change model (H ). Hence, it is used to compare which model better describes the data. B (t) = L (H t) (5) L (H t) The constant rate model is characterized by a single unknown parameter - the rate of occurrence, λ. A gamma distribution with parameters k and θ, similar to eq. 4, is used as its prior distribution. Then L (H t) = L (λ t)π(λ )dλ = Γ(N(T ) + k ( ) ) k Γ(k ) θ (. + T θ and L (H t) = T ) (N(T )+k ) L (τ,λ,λ t) (6).π(λ )π(λ )π(τ)dλ dλ dτ ( ) k ( ) k θ θ Γ(k )Γ(k ) T [ ] π(τ)γ(r (τ))γ(r (τ)). S (τ) r(τ) S (τ) r (7) (τ) τ= If the value of parameters for gamma conjugate priors are k j =.5 and θ j f or j =,,, then it is shown by Raftery and Akman (986) that the equation for Bayes factor can be simplified to B (t) = 4 π T n Γ(n + /) [ T. Γ(r (τ))γ(r (τ)) τ=. S (τ) r (τ) S (τ) r (τ) ] (8) Smaller values of Bayes factor (less than ) imply that the change model is more strongly supported by the data. Bayes factor values of less than. (or greater than ) are typically used to indicate decisive preference for one model or the other (Kass and Raftery, 995). In this study, a Bayes factor of smaller than 3 is used to indicate that data favors the change model compared to the constant rate model (i.e if the Bayes factor for an observed period is calculated to be less than 3, it indicates that a change has occurred within this observed period). 3. ALGORITHMIC IMPLEMENTATION OF THE MODEL Since all the unknown variables in the model described above (τ, λ, and λ ) are continuous and closed form solutions are not possible for marginal distributions of λ, and λ, some approximations are required to implement the algorithm for the change-point model. Additionally, overflow conditions are encountered in the algorithm, for instance when computing Γ(x) function for big values. The approximations and inputs required for the implementation of our algorithm are described below. Later, the verification and validation processes for the algorithm are described. 3.. Approximations and prior parameter inputs for the algorithm Under the first approximation, the continuous variables are discretized. Since it is not computationally possible to calculate the marginal distribution of τ over a continuous range, the algorithm is implemented on a per-day basis. Hence, probability of change happening at time τ is calculated for every day of the observation period. The proportionality in eq. 8 is converted to a probability by dividing the values for each day by the sum of values for all 4
5 Vancouver, Canada, July -5, 5 days in the observation period, such that the sum of probabilities equals one. To compute the marginal distributions for λ, λ, and β, the integral over τ is discretized to a summation by again dividing the observation period on a per-day basis. Then the marginal distributions are approximated as shown in the last step of eqs., and 4. The results are calculated at discrete values of λ, λ, and β over a certain range [X low,x high ], and the proportionality is converted to probability again by dividing each result by the sum of results over the complete range. It is ensured that the selected range [X low,x high ] for each of the variables is wide enough such that probabilities at the extreme points are essentially zero. Secondly, to address the overflow problem, the results are calculated in the log domain and then converted back to the original domain. One of the reasons for the overflow problem is computation of Γ(x) which grows very rapidly with increasing values of x. In the log domain, algorithms are available to calculate log[γ(x)] directly without first needing to calculate Γ(x). Another overflow problem is encountered when the results obtained in log domain at each discrete value are converted back to original space and summed together. Their sum could be very large and outside the floating point range of the computer. Hence, to calculate probability at each discrete value, the log results are scaled such that the sum is within the range of computation. As an example of the implementation of approximations and modifications described above, the algorithm for calculating the probability of change happening on any given day is written as: (a) Discretize the observation period on a per-day basis, τ = {τ i } T i=. (b) At each τ i, calculate prob i = log(t ) + log[γ(r (τ i ))] + log[γ(r (τ i ))] r (τ i )log[s (τ i )] (c) Find a scale such that r (τ i )log[s (τ i )]. e probi scale realmax i while preserving as many smallest values as possible. Here realmax is the largest finite floating-point number in IEEE double precision. Update prob i = prob i scale. (d) Normalize to obtain probability, prob i = eprob i i e prob i. The final step in the implementation of the algorithm is selecting the values for the prior parameters (k j and θ j ) for the gamma distribution. The values are selected as the same ones that are used in the development of the formulation for Bayes factor in eq. 8 (i.e. k j =.5 and θ j f or j =,,). Brief analysis was performed to assess the sensitivity of the prior parameter values on results, and it was found that marginal distributions of τ, λ, λ, and β did not vary significantly with different prior parameters. 3.. Verification and validation The algorithm developed above for the changepoint model is verified by comparing the results with those presented in Raftery and Akman (986). (The data used for comparison is the coal-mining disasters data described by Jarrett (979) and extended per Raftery and Akman (986)). The results are in very good agreement. The algorithms are further validated through application on simulated data. To obtain the simulated data, inter-event times are randomly generated from an exponential distribution with different rates and a known change point. The change-point model is then implemented on this data and the results for τ, λ, and λ are compared with the inputs. Here, we describe one simulation case for validation. For this validation case, inter-event duration are generated at a rate of.5 events/day and another 5 are generated at a rate of.5 events/day. A plot of cumulative number of events is shown in Fig.. The first event is assumed to occur on --. Since the event rate reduces starting from the st event, a change occurs between the th and the st event. Without any additional 5
6 Vancouver, Canada, July -5, 5 information, we cannot determine a specific date of change. Hence, it is assumed that change occurred the day before the st event which is known from the data as Number of events 5 5 Known date of change 5 5 Date Figure : Cumulative number of events for simulated data with known change on 5-6- The change-point model yields a Bayes factor of for this data, implying that data strongly supports the change-point model. The probability of change happening on a given date is plotted on Fig. 3. A 95% credible interval for dates of change is determined between and Hence, the actual date of change is within the 95% credible interval. Additionally, the date of change with highest associated probability is found to be This date exactly matches our input date of change. Thus, the Bayesian change point model correctly estimates the date of change on simulated data. The rates of event occurrence are also estimated using the model. The results of probability distributions of rates before and after the change are shown in Fig. 4. The modes of posterior distributions for rates are.54 and.57 respectively, which are in good agreement with the input rates. 4. CHANGE POINT MODEL IMPLEMENTATION IN OKLAHOMA The change-point model described in previous sections, is used to quantify the change in seismicity rates in a local region of Oklahoma. The local region that is used in this paper is a 5 km radius Number of events 9 Known date of change 95% credible interval Date Number of events x Probability of change Figure 3: Probability of change on simulated data represented by the thick solid line, with its 95% credible interval between and Probability of rate Pre change rate Post change rate Input rates λ (events/day) Figure 4: Probability of pre-change and post-change rates area centered around the location of "Well " described in Keranen et al. (3). The center of this region is located km north-west of Prague, OK at (35.56 N, W). The cumulative number of events within the chosen region from the declustered catalog are shown with the dashed line in Fig. 5. No events are observed between 974 and 9 in this region. A Bayes factor of is calculated for this data, implying that the data almost certainly supports the change-point model over a constant rate model. The probability of change occurring on any given day is then calculated and is shown in Fig. 5. 6
7 Vancouver, Canada, July -5, 5 The highest probability for change is observed on A 95% credible interval is calculated to be between 8-- and --4. This implies that there is a high chance that change in seismicity rates in this local region occurs between late 8 and early. The ratio of pre-change to post-change rate, β is shown in Fig. 7. The highest probability is calculated at a ratio of This implies that the post-change rate is about 3 times greater than the pre-change rate. Number of earthquakes % credible interval Cumulative number of earthquakes. 5 Date Figure 5: Probability of change on any given day for the local region, represented by the solid line Probability of change on date Probability β = λ /λ Figure 7: Probability of ratio of pre-change to postchange rates in the local region The probabilities of seismicity rates before and after the change are also calculated and are shown in Fig. 6. Since the pre-change rate is governed by the period of no seismicity between 974 and 9, its distribution is left-tailed. The modes for prechange and post-change rates are and.8, respectively. Probability of rate Pre change rate Post change rate Seismic Rate Figure 6: Probability of activity rates for earthquakes in the local region 5. CONCLUSIONS An algorithm for change-point detection in event rates using Bayesian statistics was developed in this paper. The Bayes factor was defined in the model to determine whether event occurrence data supported that a change had occurred in event rates. If the data supported a change model, the probability of change happening on a given day could be computed. The model could also be used to estimate the probability distributions of event rates before and after the change. The model was validated through its application on a simulated data set with known properties. After validating the model on simulated data, it was implemented on a region in Oklahoma. The model detected that a change in seismicity rates occurred sometime between late 8 and mid. This period agreed well with the time of change expected through a visual inspection of the data. The post-change rate was estimated to be approximately 3 times the pre-change rate. This high increase in activity rates can substantially increase the seismic hazard. Assuming that this local region 7
8 Vancouver, Canada, July -5, 5 is a seismic source and that the magnitude distribution is held constant from pre-change to postchange, this would imply an increase of about 3 times in rate of exceeding any ground motion at this site from this source. However, more research needs to be carried out before calculating hazard from this increased seismicity. The magnitude distribution of induced earthquakes could be different than tectonic earthquakes due to difference in b- values (in a Gutenberg-Richter relation) or due to an upper bound on earthquake magnitudes. Additionally, there could be a difference in ground motions from induced earthquakes compared to natural earthquakes (Hough, 4). The change point model does not make any association with the causes of rate change. Hence, after determining that a change has occurred, it should be linked with a physical phenomenon that might have caused this change. In the case of induced seismicity in Oklahoma, some of these physical phenomena could be change in number of injection wells, cumulative injection volume or basement pore pressure. The information about dates of change provided by the change-point model can be used to identify the causes of induced seismicity. This identification can assist in decision-making for operations potentially linked with induced seismicity and can thus be used as a tool for risk mitigation. Although the change-point model described here was applied for the case of induced seismicity, this is a versatile model with other potential applications. Some of the other applications of this model could be to detect change in storm occurrence rates to inform decisions about climate change, or to detect change in population migration rates to inform decisions about developing urban infrastructure. Thus the change-point model described in this paper can serve as a decision-support tool for a variety of applications involving occurrences of potentially non-stationary events. for nuclear facilities project. Report no., Electric Power Research Institute (EPRI). Ellsworth, W. L. (3). Injection-induced earthquakes. Science, 34(64), 594. Gardner, J. K. and Knopoff, L. (974). Is the sequence of earthquakes in southern California, with aftershocks removed, Poissonian?. Bulletin of the Seismological Society of America, 64(5), Hough, S. E. (4). Shaking from injection-induced earthquakes in the Central and Eastern United States. Bulletin of the Seismological Society of America, 4(5), Jarrett, R. G. (979). A note on the intervals between coal-mining disasters. Biometrika, 66(), Kass, R. E. and Raftery, A. E. (995). Bayes factors. Journal of the American Statistical Association, 9(43), Keranen, K. M., Savage, H. M., Abers, G. A., and Cochran, E. S. (3). Potentially induced earthquakes in Oklahoma, USA: Links between wastewater injection and the Mw 5.7 earthquake sequence. Geology, 4(6), Keranen, K. M., Weingarten, M., Abers, G. A., Bekins, B. A., and Ge, S. (4). Sharp increase in central Oklahoma seismicity since 8 induced by massive wastewater injection. Science, 345(695), Lindley, D. V. (965). Introduction to probability and statistics from Bayesian viewpoint. part inference, Vol.. Cambridge University Press Archive. Raftery, A. and Akman, V. (986). Bayesian analysis of a Poisson process with a change-point.. Biometrika, 73(), van Stiphout, T., Zhuang, J., and Marsan, D. (). Seismicity declustering. Community Online Resource for Statistical Seismicity Analysis,. 6. REFERENCES Coppersmith, K. J., Salomone, L. A., Fuller, C. W., Glaser, L. L., Hanson, K. L., Hartleb, R. D., Lettis, W. R., Lindvall, S. C., McDuffie, S. M., McGuire, R. K., and others (). Central and Eastern United States (CEUS) seismic source characterization (SSC) 8
RISK ANALYSIS AND RISK MANAGEMENT TOOLS FOR INDUCED SEISMICITY
Eleventh U.S. National Conference on Earthquake Engineering Integrating Science, Engineering & Policy June 25-29, 2018 Los Angeles, California RISK ANALYSIS AND RISK MANAGEMENT TOOLS FOR INDUCED SEISMICITY
More informationBayesian Treatment of Induced Seismicity in Probabilistic Seismic Hazard Analysis
Baker, J. W., and Gupta, A. (2016). Bayesian Treatment of Induced Seismicity in Probabilistic Seismic Hazard Analysis. Bulletin of the Seismological Society of America, (in press). Bayesian Treatment of
More informationBayesian Treatment of Induced Seismicity in Probabilistic Seismic-Hazard Analysis
Bulletin of the Seismological Society of America, Vol. 106, No. 3, pp. 860 870, June 2016, doi: 10.1785/0120150258 Bayesian Treatment of Induced Seismicity in Probabilistic Seismic-Hazard Analysis by Jack
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 informationHow will induced seismicity in Oklahoma respond to decreased saltwater injection rates?
SEISMOLOGY How will induced seismicity in Oklahoma respond to decreased saltwater injection rates? Cornelius Langenbruch* and Mark D. Zoback In response to the marked number of injection-induced earthquakes
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 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 informationRisk Treatment. Todd Shipman PhD, Alberta Geological Survey/Alberta Energy Regulator November 17 th,2017 Induced Seismicity Workshop, Yellowknife NWT
Risk Treatment Todd Shipman PhD, Alberta Geological Survey/Alberta Energy Regulator November 17 th,2017 Induced Seismicity Workshop, Yellowknife NWT Risk Management Approach to Induced Seismicity Establishing
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 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 ECAT SOFTWARE PACKAGE TO ANALYZE EARTHQUAKE CATALOGUES
THE ECAT SOFTWARE PACKAGE TO ANALYZE EARTHQUAKE CATALOGUES Tuba Eroğlu Azak Akdeniz University, Department of Civil Engineering, Antalya Turkey tubaeroglu@akdeniz.edu.tr Abstract: Earthquakes are one of
More informationA comparison of seismicity rates and fluid-injection operations in Oklahoma and California: Implications for crustal stresses
A comparison of seismicity rates and fluid-injection operations in Oklahoma and California: Implications for crustal stresses Thomas Göbel 1 Downloaded 06/03/15 to 130.118.44.150. Redistribution subject
More informationOccurrence of negative epsilon in seismic hazard analysis deaggregation, and its impact on target spectra computation
Occurrence of negative epsilon in seismic hazard analysis deaggregation, and its impact on target spectra computation Lynne S. Burks 1 and Jack W. Baker Department of Civil and Environmental Engineering,
More informationUnraveling the Mysteries of Seismicity in Oklahoma. A Story Map
Unraveling the Mysteries of Seismicity in Oklahoma A Story Map Background Why create a story map? Visualization of various types of data (graphics, text, symbols, etc.) within a geographical context to
More informationSEISMIC HAZARD CHARACTERIZATION AND RISK EVALUATION USING GUMBEL S METHOD OF EXTREMES (G1 AND G3) AND G-R FORMULA FOR IRAQ
13 th World Conference on Earthquake Engineering Vancouver, B.C., Canada August 1-6, 2004 Paper No. 2898 SEISMIC HAZARD CHARACTERIZATION AND RISK EVALUATION USING GUMBEL S METHOD OF EXTREMES (G1 AND G3)
More informationPSHA results for the BSHAP region
NATO Science for Peace and Security Programme CLOSING CONFERENCE OF THE NATO SfP 983054 (BSHAP) PROJECT Harmonization of Seismic Hazard Maps for the Western Balkan Countries October 23, 2011 Ankara, Turkey
More informationModule 7 SEISMIC HAZARD ANALYSIS (Lectures 33 to 36)
Lecture 34 Topics Module 7 SEISMIC HAZARD ANALYSIS (Lectures 33 to 36) 7.3 DETERMINISTIC SEISMIC HAZARD ANALYSIS 7.4 PROBABILISTIC SEISMIC HAZARD ANALYSIS 7.4.1 Earthquake Source Characterization 7.4.2
More informationA comparison of seismicity rates and fluid injection operations in Oklahoma and California: Implications for crustal stresses
A comparison of seismicity rates and fluid injection operations in Oklahoma and California: Implications for crustal stresses Thomas Goebel California Institute of Technology Seismological Laboratory,
More informationSupplemental Information for. Geospatial analysis of Oklahoma earthquakes ( ): Quantifying the limits of regional-scale earthquake mitigation
GSA Data Repository 8 Supplemental Information for Geospatial analysis of Oklahoma earthquakes ( - 6): Quantifying the limits of regional-scale earthquake mitigation measures Ryan M. Pollyea, Neda Mohammadi,
More informationCHARACTERIZING SPATIAL CROSS-CORRELATION BETWEEN GROUND- MOTION SPECTRAL ACCELERATIONS AT MULTIPLE PERIODS. Nirmal Jayaram 1 and Jack W.
Proceedings of the 9th U.S. National and 10th Canadian Conference on Earthquake Engineering Compte Rendu de la 9ième Conférence Nationale Américaine et 10ième Conférence Canadienne de Génie Parasismique
More informationIGC. 50 th INDIAN GEOTECHNICAL CONFERENCE PROBABILISTIC SEISMIC HAZARD ANALYSIS FOR WARANGAL CONSIDERING SINGLE SEISMOGENIC ZONING
50 th IGC 50 th INDIAN GEOTECHNICAL CONFERENCE 17 th 19 th DECEMBER 2015, Pune, Maharashtra, India Venue: College of Engineering (Estd. 1854), Pune, India PROBABILISTIC SEISMIC HAZARD ANALYSIS FOR WARANGAL
More informationAspects 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 informationCEUS Seismic Source Characterization (SSC) Project
CEUS Seismic Source Characterization (SSC) Project 2012 EERI Annual Meeting and National Earthquake Conference April 12, 2012 L. A. Salomone CEUS SSC Project Manager DRIVERS: CEUS SSC Project and NGA-East
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 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 informationPreliminary Statistics Lecture 2: Probability Theory (Outline) prelimsoas.webs.com
1 School of Oriental and African Studies September 2015 Department of Economics Preliminary Statistics Lecture 2: Probability Theory (Outline) prelimsoas.webs.com Gujarati D. Basic Econometrics, Appendix
More informationPotential Injection-Induced Seismicity Associated With Oil & Gas Development
Potential Injection-Induced Seismicity Associated With Oil & Gas Development A primer on technical & regulatory considerations informing risk management & mitigation Second Edition September, 2017 www.statesfirstinitiative.org
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 informationConsideration of prior information in the inference for the upper bound earthquake magnitude - submitted major revision
Consideration of prior information in the inference for the upper bound earthquake magnitude Mathias Raschke, Freelancer, Stolze-Schrey-Str., 6595 Wiesbaden, Germany, E-Mail: mathiasraschke@t-online.de
More informationWhat We Know (and don t know)
What We Know (and don t know) about the M5.1 La Habra Earthquake Dr. Robert Graves U.S. Geological Survey In collaboration with: Dr. Egill Hauksson and Dr. Thomas Göbel Caltech Dr. Elizabeth Cochran U.S.
More informationStandardized Tests of RELM ERF s Against Observed Seismicity
SSS Standardized Tests of RELM ERF s Against Observed Seismicity Danijel Schorlemmer 1, Dave D. Jackson 2, Matt Gerstenberger 3, Stefan Wiemer 1 1 Zürich, Swiss Seismological Service, Switzerland 2 Department
More informationSIGNAL PROCESSING AND PROBABILISTIC SEISMIC HAZARD ANALYSIS TOOLS FOR CHARACTERIZING THE IMPACT OF NEAR-FAULT DIRECTIVITY
JOINT CONFERENCE PROCEEDINGS 7th International Conference on Urban Earthquake Engineering (7CUEE) & 5th International Conference on Earthquake Engineering (5ICEE) March 3-5, 2010, Tokyo Institute of Technology,
More informationUncertainties in estimates of the occurrence rate of rare space weather events
Uncertainties in estimates of the occurrence rate of rare space weather events Jeffrey J. Love Geomagnetism Program USGS Natural Hazards Mission jlove@usgs.gov Love, J. J., 2012. Credible occurrence probabilities
More informationNew USGS Maps Identify Potential Ground- Shaking Hazards in 2017
New USGS Maps Identify Potential Ground- Shaking Hazards in 2017 New USGS Maps Identify Potential Ground-Shaking Hazards in 2017 The central U.S. faces continued hazards from human-induced earthquakes
More informationProbabilistic Seismic Hazard Analysis of Nepal considering Uniform Density Model
Proceedings of IOE Graduate Conference, 2016 pp. 115 122 Probabilistic Seismic Hazard Analysis of Nepal considering Uniform Density Model Sunita Ghimire 1, Hari Ram Parajuli 2 1 Department of Civil Engineering,
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, SMALL EARTHQUAKES
EARTHQUAKE CLUSTERS, SMALL EARTHQUAKES AND THEIR TREATMENT FOR HAZARD ESTIMATION Gary Gibson and Amy Brown RMIT University, Melbourne Seismology Research Centre, Bundoora AUTHORS Gary Gibson wrote his
More informationRecent Earthquakes in Oklahoma. and Potential for Induced Seismicity Austin Holland Oklahoma State Seismologist
Recent Earthquakes in Oklahoma and the Mid-Continent: Significance and Potential for Induced Seismicity Austin Holland Oklahoma State Seismologist Outline Seismicity Rates in Mid-continent and Oklahoma
More informationModule 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 information5. Probabilistic Seismic Hazard Analysis
Probabilistic Seismic Hazard Analysis (PSHA) proposed by C.A. Cornell (1968) used to determine the design earthquake for all locations in USA. PSHA gives a relative quantification i of the design earthquake,
More informationENGINEERING-SEISMOLOGICAL ASPECTS OF EARTHQUAKE SCENARIO DEVELOPMENT ON THE EXAMPLE OF TASHKENT, UZBEKISTAN
International Journal of Geology, Earth & Environmental Sciences ISSN: 2277-281 (Online) 218 Vol. 8 (2) May-August, pp. 3-35/Alixanovich ENGINEERING-SEISMOLOGICAL ASPECTS OF EARTHQUAKE SCENARIO DEVELOPMENT
More informationPreliminary test of the EEPAS long term earthquake forecast model in Australia
Preliminary test of the EEPAS long term earthquake forecast model in Australia Paul Somerville 1, Jeff Fisher 1, David Rhoades 2 and Mark Leonard 3 Abstract 1 Risk Frontiers 2 GNS Science 3 Geoscience
More informationSENSITIVITY OF INDUCED SEISMICITY RISK TO SOURCE CHARACTERIZATION, GROUND MOTION PREDICTION, AND EXPOSURE
Paper N 3025 Registration Code: S-J1463158968 SENSITIVITY OF INDUCED SEISMICITY RISK TO SOURCE CHARACTERIZATION, GROUND MOTION PREDICTION, AND EXPOSURE A. Gupta (1), J. W. Baker (2) (1) PhD Candidate,
More informationEstimation of reliability parameters from Experimental data (Parte 2) Prof. Enrico Zio
Estimation of reliability parameters from Experimental data (Parte 2) This lecture Life test (t 1,t 2,...,t n ) Estimate θ of f T t θ For example: λ of f T (t)= λe - λt Classical approach (frequentist
More informationA hybrid probabilistic seismic hazard analysis of a low and moderate seismic region: Sri Lanka a case study
Proceedings of the Tenth Pacific Conference on Earthquake Engineering Building an Earthquake-Resilient Pacific 6-8 November 2015, Sydney, Australia A hybrid probabilistic seismic hazard analysis of a low
More informationHow to predict the probability of a major nuclear accident after Fukushima Da
How to predict the probability of a major nuclear accident after Fukushima Dai-ichi? CERNA Mines ParisTech March 14, 2012 1 2 Frequentist approach Bayesian approach Issues 3 Allowing safety progress Dealing
More informationBayesian Inference for Clustered Extremes
Newcastle University, Newcastle-upon-Tyne, U.K. lee.fawcett@ncl.ac.uk 20th TIES Conference: Bologna, Italy, July 2009 Structure of this talk 1. Motivation and background 2. Review of existing methods Limitations/difficulties
More informationIllustrating a Bayesian Approach to Seismic Collapse Risk Assessment
2th International Conference on Applications of Statistics and Probability in Civil Engineering, ICASP2 Vancouver, Canada, July 2-5, 25 Illustrating a Bayesian Approach to Seismic Collapse Risk Assessment
More informationDisaggregation of seismic drift hazard
Disaggregation of seismic drift hazard J.W. Baker, C.A. Cornell & P. Tothong Dept. of Civil and Environmental Engineering, Stanford University, Stanford, CA, USA Keywords: Seismic Demand, Probabilistic
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 informationOriginally published as:
Originally published as: Zöller, G., Holschneider, M., Hainzl, S. (2013): The maximum earthquake magnitude in a time horizon: Theory and case studies. - Bulletin of the Seismological Society of America,
More informationChallenges and Strategies for Monitoring Induced Seismic Activity
Challenges and Strategies for Monitoring Induced Seismic Activity Designing and operating induced seismic monitoring networks to meet regulations Dario Baturan Geophysical Society of Tulsa 2015 Introduction
More informationSeismic Analysis of Structures Prof. T.K. Datta Department of Civil Engineering Indian Institute of Technology, Delhi. Lecture 03 Seismology (Contd.
Seismic Analysis of Structures Prof. T.K. Datta Department of Civil Engineering Indian Institute of Technology, Delhi Lecture 03 Seismology (Contd.) In the previous lecture, we discussed about the earth
More informationUniform Hazard Spectrum(UHS) for performance based seismic design
Uniform Hazard Spectrum(UHS) for performance based seismic design *Jun-Kyoung Kim 1), Soung-Hoon Wee 2) and Seong-Hwa Yoo 2) 1) Department of Fire Protection and Disaster Prevention, Semyoung University,
More informationDependence 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 informationThe Length to Which an Earthquake will go to Rupture. University of Nevada, Reno 89557
The Length to Which an Earthquake will go to Rupture Steven G. Wesnousky 1 and Glenn P. Biasi 2 1 Center of Neotectonic Studies and 2 Nevada Seismological Laboratory University of Nevada, Reno 89557 Abstract
More informationRegional Workshop on Essential Knowledge of Site Evaluation Report for Nuclear Power Plants.
Regional Workshop on Essential Knowledge of Site Evaluation Report for Nuclear Power Plants. Development of seismotectonic models Ramon Secanell Kuala Lumpur, 26-30 August 2013 Overview of Presentation
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 information10.1 A summary of the Virtual Seismologist (VS) method for seismic early warning
316 Chapter 10 Conclusions This final Chapter is made up of the following: a summary of the Virtual Seismologist method for seismic early warning, comments on implementation issues, conclusions and other
More informationAssessment and Mitigation of Ground Motion Hazards from Induced Seismicity. Gail M. Atkinson
Assessment and Mitigation of Ground Motion Hazards from Induced Seismicity Gail M. Atkinson (with acknowledgement to many collaborators, especially Ghofrani and Assatourians) NSERC/TransAlta/Nanometrics
More informationOPTIMIZATION OF RESPONSE SIMULATION FOR LOSS ESTIMATION USING PEER S METHODOLOGY
13 th World Conference on Earthquake Engineering Vancouver, B.C., Canada August 1-6, 2004 Paper No. 1066 OPTIMIZATION OF RESPONSE SIMULATION FOR LOSS ESTIMATION USING PEER S METHODOLOGY Hesameddin ASLANI
More informationNew Madrid and Central U.S. Region Earthquake Hazard
New Madrid and Central U.S. Region Earthquake Hazard Rob Williams U.S. Geological Survey, Golden, CO March 16, 2017 Reelfoot Lake, TN U.S. Department of the Interior U.S. Geological Survey 1811-12 New
More informationSEISMIC HAZARD ANALYSIS. Instructional Material Complementing FEMA 451, Design Examples Seismic Hazard Analysis 5a - 1
SEISMIC HAZARD ANALYSIS Instructional Material Complementing FEMA 451, Design Examples Seismic Hazard Analysis 5a - 1 Seismic Hazard Analysis Deterministic procedures Probabilistic procedures USGS hazard
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 informationA METHODOLOGY FOR ASSESSING EARTHQUAKE-INDUCED LANDSLIDE RISK. Agency for the Environmental Protection, ITALY (
A METHODOLOGY FOR ASSESSING EARTHQUAKE-INDUCED LANDSLIDE RISK Roberto W. Romeo 1, Randall W. Jibson 2 & Antonio Pugliese 3 1 University of Urbino, ITALY (e-mail: rwromeo@uniurb.it) 2 U.S. Geological Survey
More informationRisk Evaluation. Todd Shipman PhD, Alberta Geological Survey/Alberta Energy Regulator November 17 th,2017 Induced Seismicity Workshop, Yellowknife NWT
Risk Evaluation Todd Shipman PhD, Alberta Geological Survey/Alberta Energy Regulator November 17 th,2017 Induced Seismicity Workshop, Yellowknife NWT Risk Management Approach to Induced Seismicity Establishing
More informationEnhanced remote earthquake triggering at fluid injection sites in the Midwestern US
Enhanced remote earthquake triggering at fluid injection sites in the Midwestern US Nicholas van der Elst 1 Heather Savage 1 Katie Keranen 2 Geoff Abers 1 1 Lamont-Doherty Earth Observatory 2 Cornell University
More informationStochastic Renewal Processes in Structural Reliability Analysis:
Stochastic Renewal Processes in Structural Reliability Analysis: An Overview of Models and Applications Professor and Industrial Research Chair Department of Civil and Environmental Engineering University
More informationBayesian Inference for Normal Mean
Al Nosedal. University of Toronto. November 18, 2015 Likelihood of Single Observation The conditional observation distribution of y µ is Normal with mean µ and variance σ 2, which is known. Its density
More information(1) Introduction to Bayesian statistics
Spring, 2018 A motivating example Student 1 will write down a number and then flip a coin If the flip is heads, they will honestly tell student 2 if the number is even or odd If the flip is tails, they
More informationSEISMIC-HAZARD ASSESSMENT: Conditional Probability
Conditional Probability SEISMIC-HAZARD ASSESSMENT: Conditional Probability Supplies Needed calculator PURPOSE Previous exercises in this book have outlined methods for inferring the patterns and history
More informationBC HYDRO SSHAC LEVEL 3 PSHA STUDY METHODOLOGY
BC HYDRO SSHAC LEVEL 3 PSHA STUDY METHODOLOGY M. W. McCann, Jr. 1, K. Addo 2 and M. Lawrence 3 ABSTRACT BC Hydro recently completed a comprehensive Probabilistic Seismic Hazard Analysis (PSHA) to evaluate
More informationI N T R O D U C T I O N T O P R O B A B I L I S T I C S E I S M I C H A Z A R D A N A LY S I S
I N T R O D U C T I O N T O P R O B A B I L I S T I C S E I S M I C H A Z A R D A N A LY S I S J A C K W. B A K E R Copyright 2015 Jack W. Baker Preferred citation for this document: Baker, Jack W. (2015)
More informationProbabilistic Earthquake Risk Assessment of Newcastle and Lake Macquarie Part 1 Seismic Hazard.
Probabilistic Earthquake Risk Assessment of Newcastle and Lake Macquarie Part 1 Seismic Hazard. T. Dhu, D. Robinson, C. Sinadinovski, T. Jones, A. Jones & J. Schneider Geoscience Australia, Canberra, Australia.
More informationEARTHQUAKE HAZARD ASSESSMENT IN KAZAKHSTAN
EARTHQUAKE HAZARD ASSESSMENT IN KAZAKHSTAN Dr Ilaria Mosca 1 and Dr Natalya Silacheva 2 1 British Geological Survey, Edinburgh (UK) imosca@nerc.ac.uk 2 Institute of Seismology, Almaty (Kazakhstan) silacheva_nat@mail.ru
More informationAn Introduction to Probabilistic Seismic Hazard Analysis (PSHA) Jack W. Baker
An Introduction to Probabilistic Seismic Hazard Analysis (PSHA) Jack W. Baker Version 1.3 October 1 st, 2008 Acknowledgements This document is based in part on an excerpt from a report originally written
More informationOperational Earthquake Forecasting in Italy: perspectives and the role of CSEP activities
Operational Earthquake Forecasting in Italy: perspectives and the role of CSEP activities Warner Marzocchi, Istituto Nazionale di Geofisica e Vulcanologia, Italy The research was developed partially within
More informationBayesian network for infrastructure seismic risk assessment and decision support
Bayesian network for infrastructure seismic risk assessment and decision support Michelle T. Bensi Armen Der Kiureghian University of California, Berkeley Daniel Straub Technical University of Munich Research
More informationStep-Stress Models and Associated Inference
Department of Mathematics & Statistics Indian Institute of Technology Kanpur August 19, 2014 Outline Accelerated Life Test 1 Accelerated Life Test 2 3 4 5 6 7 Outline Accelerated Life Test 1 Accelerated
More informationI N T R O D U C T I O N T O P R O B A B I L I S T I C S E I S M I C H A Z A R D A N A LY S I S
I N T R O D U C T I O N T O P R O B A B I L I S T I C S E I S M I C H A Z A R D A N A LY S I S J A C K W. B A K E R Copyright 2013 Jack W. Baker Preferred citation for this document: Baker, Jack W. (2013)
More informationprocesses Sai Xiao, Athanasios Kottas and Bruno Sansó Abstract
Nonparametric Bayesian modeling and inference for renewal processes Sai Xiao, Athanasios Kottas and Bruno Sansó Abstract We propose a flexible approach to modeling and inference for renewal processes.
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 informationUGRC 144 Science and Technology in Our Lives/Geohazards
UGRC 144 Science and Technology in Our Lives/Geohazards Session 3 Understanding Earthquakes and Earthquake Hazards Lecturer: Dr. Patrick Asamoah Sakyi Department of Earth Science, UG Contact Information:
More informationEstimation of Regional Seismic Hazard in the Korean Peninsula Using Historical Earthquake Data between A.D. 2 and 1995
Bulletin of the Seismological Society of America, Vol. 94, No. 1, pp. 269 284, February 2004 Estimation of Regional Seismic Hazard in the Korean Peninsula Using Historical Earthquake Data between A.D.
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 informationCoulomb 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 informationIntroduction to Bayesian Statistics
Bayesian Parameter Estimation Introduction to Bayesian Statistics Harvey Thornburg Center for Computer Research in Music and Acoustics (CCRMA) Department of Music, Stanford University Stanford, California
More informationACCOUNTING FOR SITE EFFECTS IN PROBABILISTIC SEISMIC HAZARD ANALYSIS: OVERVIEW OF THE SCEC PHASE III REPORT
ACCOUNTING FOR SITE EFFECTS IN PROBABILISTIC SEISMIC HAZARD ANALYSIS: OVERVIEW OF THE SCEC PHASE III REPORT Edward H FIELD 1 And SCEC PHASE III WORKING GROUP 2 SUMMARY Probabilistic seismic hazard analysis
More informationA BAYESIAN MATHEMATICAL STATISTICS PRIMER. José M. Bernardo Universitat de València, Spain
A BAYESIAN MATHEMATICAL STATISTICS PRIMER José M. Bernardo Universitat de València, Spain jose.m.bernardo@uv.es Bayesian Statistics is typically taught, if at all, after a prior exposure to frequentist
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 informationDevelopment of Probabilistic Seismic Hazard Analysis for International Sites, Challenges and Guidelines
Development of Probabilistic Seismic Hazard Analysis for International Sites, Challenges and Guidelines ABSTRACT Dr. Antonio Fernandez Ares Paul C. Rizzo Associates, Inc. 500 Penn Center Boulevard, Suite
More informationSEISMIC INPUT FOR CHENNAI USING ADAPTIVE KERNEL DENSITY ESTIMATION TECHNIQUE
SEISMIC INPUT FOR CHENNAI USING ADAPTIVE KERNEL DENSITY ESTIMATION TECHNIQUE G. R. Dodagoudar Associate Professor, Indian Institute of Technology Madras, Chennai - 600036, goudar@iitm.ac.in P. Ragunathan
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 informationGEM's community tools for probabilistic seismic hazard modelling and calculation
GEM's community tools for probabilistic seismic hazard modelling and calculation Marco Pagani, GEM Secretariat, Pavia, IT Damiano Monelli, GEM Model Facility, SED-ETH, Zürich, CH Graeme Weatherill, GEM
More information6 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 informationBayesian Statistics as an Alternative for Analyzing Data and Testing Hypotheses Benjamin Scheibehenne
Bayesian Statistics as an Alternative for Analyzing Data and Testing Hypotheses Benjamin Scheibehenne http://scheibehenne.de Can Social Norms Increase Towel Reuse? Standard Environmental Message Descriptive
More informationSeismic Source Characterization in Siting New Nuclear Power Plants in the Central and Eastern United States
Seismic Source Characterization in Siting New Nuclear Power Plants in the Central and Eastern United States ABSTRACT : Yong Li 1 and Nilesh Chokshi 2 1 Senior Geophysicist, 2 Deputy Director of DSER Nuclear
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 informationThomas Bayes versus the wedge model: An example inference using a geostatistical prior function
Thomas Bayes versus the wedge model: An example inference using a geostatistical prior function Jason M. McCrank, Gary F. Margrave, and Don C. Lawton ABSTRACT The Bayesian inference is used to estimate
More informationCharacterizing and Responding to Seismic Risk Associated with Earthquakes Potentially Triggered by Saltwater Disposal and Hydraulic Fracturing
Characterizing and Responding to Seismic Risk Associated with Earthquakes Potentially Triggered by Saltwater Disposal and Hydraulic Fracturing R. J. Walters 1, M. D. Zoback 1, J. W. Baker 2 and G. C. Beroza
More information