Residual Analysis Methods for Space-time Point Processes with Applications to Earthquake Forecast Models in California

Robert A. Clements, Frederic P. Schoenberg, & Danijel Schorlemmer

Published 2011, SCEC Contribution #1490

Modern, powerful techniques for the residual analysis of spatial-temporal point process models are reviewed and compared. These methods are applied to California earthquake forecast models used in the Collaboratory for the Study of Earthquake Predictability (CSEP). Assessments of these earthquake forecasting models have previously been performed using simple, low-power means such as the L-test
and N-test. We instead propose residual methods based on rescaling, thinning, superposition, weighted K-functions, and deviance residuals. Rescaled residuals can be useful for assessing the overall fit of a model, but as with thinning and superposition, rescaling is generally impractical when λ is volatile. While residual thinning and superposition may be useful for identifying spatial locations where a model fits poorly, these methods have limited power when the modeled conditional intensity assumes extremely low or high values somewhere in the observation region, and this is commonly the case for earthquake forecasting models. A recently proposed hybrid method of thinning and superposition, called super-thinning, is a more powerful alternative. The weighted K-function is powerful for evaluating the degree of clustering or inhibition in a model. Competing models are
also compared using pixel-based approaches, such as Pearson residuals and deviance residuals. The different residual analysis techniques are demonstrated using the CSEP models and are used to highlight certain deficiencies in the models, such as the overprediction of seismicity in inter-fault zones for the model proposed by Helmstetter,
Kagan and Jackson (2007), the underprediction of the model proposed by Kagan, Jackson and Rong (2007) in forecasting seismicity around the Imperial, Laguna Salada, and Panamint clusters, and the underprediction of the model proposed by Shen, Jackson and Kagan (2007) in forecasting seismicity around the Laguna Salada, Baja, and Panamint clusters.

Citation
Clements, R. A., Schoenberg, F. P., & Schorlemmer, D. (2011). Residual Analysis Methods for Space-time Point Processes with Applications to Earthquake Forecast Models in California. Annals of Applied Statistics, 5(4), 2549-2571. doi: 10.1214/11-AOAS487.