Читать книгу Statistical Methods and Modeling of Seismogenesis - Eleftheria Papadimitriou - Страница 15

The kernel estimate of magnitude CDF and such an estimate based on the exponential model [1.20], both obtained for the actual seismicity case from the Northern Aegean Sea (Greece), are represented in Figure 1.4, taken from Lasocki and Papadimitriou (2006). We can see that the kernel estimate better fits the observed data, though the difference between these two estimates in the larger magnitude range is pretty tiny. However, due to the significant impact of the magnitude CDF on the parameters of hazard [1.17] and [1.18], the difference between the mean return period estimates is dramatic. For instance, for magnitude 6.5, the return period obtained from the exponential distribution model is close to 100 years. The kernel estimate of this return period is less than 10 years. Lasocki and Papadimitriou (2006) compared the “exponential” and “kernel” estimates of the mean return period for different magnitudes, with the estimates drawn from the actual observations done in the preceding 94 years (Figure 1.5). The differences between the “kernel” estimates and the assessments from actual observations were insignificant compared to the huge deviation of the “exponential” estimate. This led to the conclusion that the “kernel” estimate was much better than the “exponential” one.

Also, many other studies indicated big differences between the “exponential” and “kernel” estimates of hazard parameters. The mentioned Monte Carlo analyses by Kijko et al. (2001), and the actual data studies by Lasocki and Papadimitriou (2006), suggest that in the case when these estimates differ, the “kernel” estimate is more accurate. All of this favors the kernel estimation of magnitude distribution functions for the seismic hazard assessment.

The PSHA, which uses the kernel estimation of magnitude distribution as an alternative to the parametric model [1.19] and [1.20], has been implemented on the IS-EPOS Platform (tcs.ah-epos.eu, Orlecka-Sikora et al. 2020). The kernel estimation of magnitude distribution is also applied in the SHAPE software package for time-dependent seismic hazard analysis (Leptokaropoulos and Lasocki 2020). SHAPE is open-source, downloadable from https://git.plgrid.pl/projects/EA/repos/seraapplications/browse/SHAPE_Package.

Figure 1.4. Comparison of the estimates obtained from the kernel estimation method (solid lines) and from the exponential model of magnitude distribution (dashed lines) for the data from the Northern Aegean Sea. Left – the CDF estimates juxtaposed with the cumulative histogram of magnitude data. Right – the mean return period estimates. Reprinted from Lasocki and Papadimitriou (2006, Figures 5b and 6b)

Figure 1.5. “Kernel” (circles) and “exponential” (squares) estimates of the mean return periods juxtaposed with the mean return period estimates drawn from the actual 94 year-long observations of seismicity in the Northern Aegean Sea (diamonds). Reprinted from Lasocki and Papadimitriou (2006, Figure 7b)