Extreme Value Methods for Modelling Historical Series of Large Volcanic Magnitudes

Volcanic eruptions are among the most extreme events on Earth and it seems natural to make use of the theory of extreme values to improve understanding of volcanic processes. The dataset we use is a catalogue of large eruptions over the last two millennia, in which the date of occurrence and magnitude are recorded. The dataset is affected by a recording bias, mostly for eruptions of lower magnitude, though this under-recording process largely disappears in the most recent 400 years. Coles and Sparks (2006) modelled these data, via maximum likelihood, using a Poisson process motivated by extreme value theory, with an intensity function that takes into account the recording bias. Nevertheless, the fitted model did not seem entirely consistent with the observed data, since this intensity function does not represent effectively the temporal evolution of the censoring effect in the recording process.