Robust estimation of magnitude-frequency relationship parameters

The precise estimation of the a and b parameters of Richter's magnitude-frequency relationship is of primary importance, since the evaluation of seismicity and assessment of seismic hazard depend on these two parameters. In the literature two popular methods of estimation are available for the estimation of these parameters, namely: least squares and maximum likelihood. However, in implementing these statistical methods, engineers very seldom check the validity of the underlying assumptions with respect to the available data and this may lead to serious problems. Under non-normality least squares estimators (LSEs) are neither efficient nor robust and maximum likelihood estimators (MLEs) are elusive due to numerous complexities. A robust estimation procedure, the modified maximum likelihood method (MML), can be utilized to estimate the unknown parameters a and b in such situations. The resulting estimators are explicit functions of sample observations and are shown to be considerably more efficient than the commonly used least squares estimators. In addition, we demonstrate that the MML estimators are more appropriate to estimate the parameters of Richter's magnitude-frequency relationship based on the comparison of their performance with those of the least squares estimators by using the seismic database on earthquakes recorded in Turkey.


Fault-based probabilistic seismic hazard assessment of the eastern Makran subduction and the Chaman transform fault, Pakistan: Emphasis on the source characterization of megathrust
Shah, Syed Tanvir; Özacar, Atilla Arda; Gülerce, Zeynep (Elsevier BV, 2021-01-01)
Seismic source characterization (SSC) for probabilistic seismic hazard assessment (PSHA) in regions characterized by subduction megathrust involves a considerable ambiguity. Lack of detailed geologic, seismic, and geodetic data increases the uncertainties. The enigma is enhanced in regions where thin-skinned accretionary prism faults are part of active deformation. In this study, a planar SSC model for seismically active eastern Makran subduction zone, its associated accretionary prism faults and Chaman tra...
Quasi-fission and fusion-fission reactions in Ca-48+Pb-208 collisions at E-c.m.=190 MeV
Yilmaz, O.; Turan, Gürsevil; Yilmaz, B. (Springer Science and Business Media LLC, 2020-02-05)
The yields of quasifission fragments in the Ca-48+Pb-208 are investigated in the framework of the stochastic mean field (SMF) approach based on the Skyrme-TDHF calculations. A method which combines the microscopic time-dependent Hartree-Fock (TDHF) theory with a statistical de-excitation model, GEMINI++, is used for calculating the yields of fusion-fission fragments after de-excitation processes of primary reaction products which is assumed to be dominated by neutron emission in the present system. Results ...
Effective Stress Based Large Strain Modulus Degradation Response of Fully Saturated Clean Sands
Bilge, Habib Tolga; Çetin, Kemal Önder (2015-9-11)
Current state of practice to assess ground response under seismic loading is mostly based on well-known and widely-used total stress-based equivalent linear assessment methodology, as part of which, mostly, strain compatible shear modulus and damping relationships need to be used. As implied by its name, "total stress", accumulation of excess pore water pressure which is a significant component of cyclic response, is not considered. Neglection of this component will produce unconservatively biased deformati...
Numerical calculation of rough turbulent boundary layer for incompressible fluids
Atay, Görkem; Bayer, Özgür; Çıray, Cahit; Department of Mechanical Engineering (2018)
In this thesis, a numerical code for calculating turbulent boundary layer parameters is developed. Two-dimensional turbulent flow over a rough flat plate with zero pressure gradient is numerically solved by using the integral method. Fluid mechanics formulations, such as momentum integral equation are coupled with the results of Nikuradse’s experiment. Throughout the calculation, turbulent boundary layer parameters, such as boundary layer thickness, displacement thickness, momentum thickness, local skin-fri...
Nonnormal regression. I. Skew distributions
İslam, Muhammed Qamarul; Yildirim, F (2001-01-01)
In a linear regression model of the type y = thetaX + e, it is often assumed that the random error e is normally distributed. In numerous situations, e.g., when y measures life times or reaction times, e typically has a skew distribution. We consider two important families of skew distributions, (a) Weibull with support IR: (0, infinity) on the real line, and (b) generalised logistic with support IR: (-infinity, infinity). Since the maximum likelihood estimators are intractable in these situations, we deriv...
Citation Formats
M. S. Yücemen and A. Akkaya, “Robust estimation of magnitude-frequency relationship parameters,” STRUCTURAL SAFETY, pp. 32–39, 2012, Accessed: 00, 2020. [Online]. Available: