Matrice des liens Tous les phénomènes mensuels croisés deux à deux, avec un contrôle des fausses découvertes

Données. 24 séries mensuelles de l’observatoire, depuis janvier 1979 (début de l’ère satellitaire) ou leur premier mois ; moyennes mensuelles des séries quotidiennes et horaires ; nombre mensuel de séismes M ≥ 6 du catalogue USGS ; mois incomplets écartés.

Prétraitement. Trous de trois mois au plus interpolés ; CO₂ et CH₄ pris comme croissance sur 12 mois ; anomalie par rapport à la moyenne de chaque mois de l’année (cycle saisonnier retiré) ; moyenne glissante centrée de 121 mois retranchée (tendances et variations de plus de 10 ans retirées) ; séries centrées réduites.

Statistique. Pour chaque paire, corrélation de Pearson pour chaque décalage de −12 à +12 mois, sur la plus longue période commune sans trou (au moins 180 mois) ; on retient la plus forte en valeur absolue.

Significativité. 2 999 substituts de la première série à phases aléatoires (Ebisuzaki, 1997) : ils gardent son spectre, donc son autocorrélation, et l’on refait la même recherche du meilleur décalage ; p = (1 + nombre de substituts au moins aussi corrélés) / 3 000. Cela tient compte à la fois de la persistance des séries et du choix du décalage.

Tests multiples. 276 paires testées : à 5 %, une douzaine sortiraient par hasard. Le taux de fausses découvertes de Benjamini et Hochberg (1995), recommandé pour ces criblages (Wilks, 2016), limite à 5 % en moyenne la part de faux liens parmi ceux retenus (q ≤ 0,05).

Sans l’ENSO. Chaque série est régressée par moindres carrés sur Niño 3.4 décalé de 0 à 6 mois, puis la matrice est recalculée sur les résidus. Méthode linéaire simple ; Compo et Sardeshmukh (2010) proposent une version dynamique plus complète.

Calcul par le pipeline Python de l’observatoire (numpy, data-pipeline/py/obspipe/links.py), graine fixe : le résultat est reproductible.