Emulation of simulated earthquake catalogues

• Autorzy: Gallucci M. , Guillas S. , Vallianatos F.
• Języki publikacji: EN

Abstrakty

EN

In earthquake occurrence studies, the so-called q value can be considered both as one of the parameters describing the distribution of inter-event times and as an index of non-extensivity. Using simulated datasets, we compare four kinds of estimators, based on principle of maximum entropy (POME), method of moments (MOM), maximum likelihood (MLE), and probability weighted moments (PWM) of the parameters (q and τ0) of the distribution of inter-events times, assumed to be a generalized Pareto distribution (GPD), as defined by Tsallis (1988) in the frame of non-extensive statistical physics. We then propose to use the unbiased version of PWM estimators to compute the q value for the distribution of inter-event times in a realistic earthquake catalogue simulated according to the epidemic type aftershock sequence (ETAS) model. Finally, we use these findings to build a statistical emulator of the q values of ETAS model. We employ treed Gaussian processes to obtain partitions of the parameter space so that the resulting model respects sharp changes in physical behaviour. The emulator is used to understand the joint effects of input parameters on the q value, exploring the relationship between ETAS model formulation and distribution of inter-event times.

Słowa kluczowe

EN

probability weighted moments; q-statistics; inter-event times; ETAS model; emulator

• Wydawca: Instytut Geofizyki PAN ; Springer
• Czasopismo: Acta Geophysica
• Rocznik: 2012
• Tom: Vol. 60, no. 3
• Strony: 699--718
• Opis fizyczny: Bibliogr. 35 poz.

Twórcy

autor

Gallucci M.

autor

Guillas S.

autor

Vallianatos F.

• Dipartimento di Statistica "G. Parenti", Universita di Firenze, Firenze, Italy,

Bibliografia

• Abe, S., and N. Suzuki (2003), Law for the distance between successive earthquakes, J. Geophys. Res. 108, B2, 2113-2116, DOI: 10.1029/2002JB 002220.
• Abe, S., and N. Suzuki (2005), Scale-free statistics of time interval between successive earthquakes, Physica A 350, 2-4, 588-596, DOI: 10.1016/j.physa.2004.10.040.
• Gramacy, R.B. (2007), tgp: An R package for Bayesian nonstationary, semiparametric nonlinear regression and design by treed Gaussian process models, J. Stat. Softw. 19, 9, 1-46.
• Gramacy, R.B., and H.K.H. Lee (2008), Bayesian treed Gaussian process models with an application to computer modeling, J. Am. Stat. Assoc. 103, 483, 1119-1130, DOI: 10.1198/016214508000000689.
• Gramacy, R.B., and M. Taddy (2010), Categorical inputs, sensitivity analysis, optimization and importance tempering with tgp version 2, an R package for treed Gaussian process models, J. Stat. Softw. 33, 6, 1-48.
• Greenwood, J.A., J.M. Landwehr, N.C. Matalas, and J.R. Wallis (1979), Probability weighted moments: Definition and relation to parameters of several distributions expressable in inverse form, Water Resour. Res. 15, 5, 1049-1054, DOI: 10.1029/WR015i005p01049.
• Harte, D. (2010), PtProcess: An R package for modelling marked point processes indexed by time, J. Stat. Softw. 35, 8, 1-32.
• Hosking, J.R.M., and J.R. Wallis (1987), Parameter and quantile estimation for the generalized Pareto distribution, Technometrics 29, 3, 339-349, DOI: 10.2307/1269343.
• Kagan, Y.Y., and F. Shoenberg (2001), Estimation of the upper cutoff parameter for the tapered Pareto distribution, J. Appl. Probab. 38A, 158-175, DOI: 10.1239/jap/1085496599.
• McKay, M.D., R.J. Beckman, and W.J. Conover (1979), A comparison of three methods for selecting values of input variables in the analysis of output from a computer code, Technometrics 21, 2, 239-245, DOI: 10.2307/1268522.
• Naylor, M., I.G. Main, and S. Touati (2009), Quantifying uncertainty in mean earthquake interevent times for a finite sample, J. Geophys. Res. 114, B01316, DOI: 10.1029/2008JB005870.
• O’Hagan, A. (2006), Bayesian analysis of computer code outputs: A tutorial, Reliab. Eng. Syst. Safe. 91, 10-11, 1290-1300, DOI: 10.1016/j.ress.2005.11.025.
• Ogata, Y. (1988), Statistical models for earthquake occurrences and residual analysis for point processes, J. Am. Stat. Assoc. 83, 401, 9-27, DOI: 10.2307/2288914.
• Ogata, Y. (1998), Space-time point-process models for earthquake occurrences, Ann. Inst. Statist. Math. 50, 2, 379-402, DOI: 10.1023/A:1003403601725.
• Ogata, Y., and J. Zhuang (2006), Space-time ETAS models and an improved extension, Tectonophysics 413, 1-2, 13-23, DOI: 10.1016/j.tecto.2005.10.016.
• Pisarenko, V.F., D. Sornette, and M.V. Rodkin (2010), Distribution of maximum earthquake magnitudes in future time intervals: application to the seismicity of Japan (1923-2007), Earth Planets Space 62, 7, 567-578, DOI: 10.5047/eps.2010.06.003.
• Rougier, J., S. Guillas, A. Maute, and A.D. Richmond (2009), Expert knowledge and multivariate emulation: The thermosphere-ionosphere electrodynamics general circulation model (TIE-GCM), Technometrics 51, 4, 414-424, DOI: 10.1198/TECH.2009.07123.
• Sacks, J., W.J. Welch, T.J. Mitchell, and H.P. Wynn (1989), Design and analysis of computer experiments, Stat. Sci. 4, 4, 409-423, DOI: 10.1214/ss/1177012413.
• Shalizi, C.R. (2007), Maximum likelihood estimation for q-exponential (Tsallis) distributions, arXiv:math/0701854v2 [math.ST].
• Singh, V.P., and H. Guo (1997), Parameter estimation for 2-parameter generalized pareto distribution by POME, Stoch. Hydrol. Hydraul. 11, 3, 211-227, DOI: 10.1007/BF02427916.
• Telesca, L. (2010a), Nonextensive analysis of seismic sequences, Physica A 389, 9, 1911-1914, DOI: 10.1016/j.physa.2010.01.012.
• Telesca, L. (2010b), A non-extensive approach in investigating the seismicity of L’Aquila area (central Italy), struck by the 6 April 2009 earthquake (ML = 5.8), Terra Nova 22, 2, 87-93, DOI: 10.1111/j.1365-3121.2009. 00920.x.
• Touati, S., M. Naylor, and I.G. Main (2009), Origin and nonuniversality of the earthquake interevent time distribution, Phys. Rev. Lett. 102, 16, 168501, DOI: 10.1103/PhysRevLett.102.168501.
• Touati, S., M. Naylor, I.G. Main, and M. Christie (2011), Masking of earthquake triggering behavior by a high background rate and implications for epidemic-type aftershock sequence inversions, J. Geophys. Res. 116, B03304, DOI: 10.1029/2010JB007544.
• Tsallis, C. (1988), Possible generalization of Boltzmann–Gibbs statistics, J. Stat. Phys. 52, 1-2, 479-487, DOI: 10.1007/BF01016429.
• Tsallis, C. (1999), Nonextensive statistics: theoretical, experimental and computanional evidences and connections, Braz. J. Phys. 29, 1, DOI: 10.1590/S0103-97331999000100002.
• Tsallis, C. (2009), Introduction to Nonextensive Statistical Mechanics: Approaching a Complex World, Springer, New York, DOI: 10.1007/978-0-387-85359-8.
• Vallianatos, F. (2009), A non-extensive approach to risk assessment, Nat. Hazards Earth Syst. Sci. 9, 1, 211-216, DOI: 10.5194/nhess-9-211-2009.
• Vallianatos, F. (2011), A non extensive statistical physics approach to the polarity reversals of the geomagnetic field, Physica A 390, 10, 1773-1778, DOI: 10.1016/j.physa.2010.12.040.
• Vallianatos, F., and P. Sammonds (2010), Is plate tectonics a case of non-extensive thermodynamics?, Physica A 389, 21, 4989-4993, DOI:10.1016/j.physa.2010.06.056.
• Vallianatos, F., and P. Sammonds (2011), A non-extensive statistics of the faultpopulation at the Valles Marineris extensional province, Mars, Tectonophysics 509, 1-2, 50-54, DOI: 10.1016/j.tecto.2011.06.001.
• Vallianatos, F., E. Kokinou, and P. Sammonds (2011a), Non-extensive statistical physics approach to fault population distribution. A case study from the Southern Hellenic Arc (Central Crete), Acta Geophys. 59, 4, 770-784, DOI: 10.2478/s11600-011-0015-3.
• Vallianatos, F., D. Triantis, and P. Sammonds (2011b), Non-extensivity of the isothermal depolarization relaxation currents in uniaxial compressed rocks, Europhys. Lett. 94, 6, 68008, DOI: 10.1209/0295-5075/94/68008.
• Wang, Q., D.D. Jackson, and J. Zhuang (2010a), Missing links in earthquake clustering models, Geophys. Res. Lett. 37, L21307, DOI: 10.1029/2010GL044858.
• Wang, Q, F.P. Schoenberg, and D.D. Jackson (2010b), Standard errors of parameter estimates in the ETAS model, Bull. Seismol. Soc. Am. 100, 5A, 1989-2001, DOI: 10.1785/0120100001.