Abstract
Most point process models for earthquakes in the literature assume that the magnitude is independent and identically distributed. This potentially hinders the ability of the model to describe the main features of data sets containing multiple earthquake mainshock aftershock sequences in succession. This study presents a novel multivariate fractional Hawkes process model designed to capture magnitude dependent triggering behaviour by incorporating history dependence into the magnitude distribution. This is done by discretising the magnitude range into disjoint intervals and modelling events with magnitude in these ranges as the subprocesses of a mutually exciting Hawkes process using the Mittag-Leffler density as the kernel function so that the point process has a history dependent mark distribution. We apply this model to two data sets, Japan and the Middle America Trench, both containing multiple mainshock aftershock sequences and compare it to the existing ETAS model by using information criteria, residual diagnostics and retrospective prediction performance. We find that for both data sets all metrics indicate that the multivariate fractional Hawkes process performs favourably against the ETAS model due to its history dependent magnitude distribution. Furthermore, we are able to infer characteristics of the data sets that cannot be inferred from the ETAS model.