Giorni orsono Guido Guidi, nel post dedicato all’ottobrata ed al nuovo indice circolatorio OPI, ha segnalato un sito della Rutgers University dedicato alle coperture nevose. In tale sito, alla pagina, è presente una interessante tabella che riporta i dati mensili di fonte NOAA delle superfici innevate nell’emisfero boreale per il periodo 1966 – 2013 su cui ho mi è parso interessante condurre alcune analisi del tutto preliminari di cui qui sotto riporto i metodi impiegati ed i risultati ottenuti.
Anzitutto per ragioni di qualità dei dati (i dati dei primi anni solo parzialmente presenti e con alcuni spikes), ho preferito analizzare i soli dati dal gennaio 1973 al dicembre 2012. Da tali dati mensili ho ricavato gli innevamenti medi annuali espressi in milioni di kmq. Le medie, ad una prima analisi visuale, mostrano la presenza di un evidente discontinuità (alias change point, alias breakpoint), con l’innevamento che cala drasticamente a decorrere da fine anni ’80. Per sostituire all’analisi visiva un criterio più oggettivo ho applicato il test di Bai e Perron presente nel software Strucchange di R Cran e che è deputato all’individuazione di discontinuità singole o multiple. Il test, i cui risultati sono illustrati in figura 1, individua un’unica discontinuità che con una confidenza del 99% ricade fra l’aprile 1986 ed il settembre 1988 (linea orizzontale rossa) e che ha come momento più probabile di accadimento il marzo 1987. Il livello di confidenza molto alto (99%) ed il fatto che la banda di confidenza sia così ristretta rendono il test assai probante. Si noti che in virtù del “gradino” del 1987 la nevosità media annua passa dai 25.6 milioni di kmq del periodo gennaio 1973- marzo 1987 ed i 24.7 milioni di kmq del periodo aprile 1987-dicembre 2012.