Multivariate stochastic modelling: towards forecasts of paroxysmal phases at Stromboli