ANDENNE (BE)

Méthode des filtres multiples (Dziewonski, Bloch & Landisman, 1969), appliquée dans votre navigateur au signal brut de la station (100 / 200 échantillons par seconde), le fichier même des sismogrammes ci-dessus : tendance retirée, extrémités adoucies sur 2 s, spectre divisé fréquence par fréquence par la réponse de l’instrument (sensibilité, pôles et zéros du capteur, et filtres numériques du numériseur, d’après l’inventaire d’ORFEUS), puis, pour chacune des 81 fréquences de 50 à 0,5 Hz, filtre gaussien exp(−α ((f − f₀)/f₀)²) avec α = 18 et enveloppe du signal analytique ; la puissance (carré de l’enveloppe) est lissée en temps par une gaussienne d’écart type égal à celui du filtre, √(2α) T₀ / 2π (opérateur de lissage de Torrence & Webster, 1999, repris par Grinsted et al., 2004), et l’amplitude affichée en est la racine, soit l’amplitude efficace à la résolution du filtre, sans les évanouissements de l’enveloppe instantanée ; un échantillon toutes les 0,04 s. Ce calcul équivaut, à la normalisation près, à une transformée en ondelettes continue avec l’ondelette de Morlet, ω₀ = √(2α) = 6, la valeur usuelle (Torrence & Compo, 1998). Le filtre est tenu pour nul là où le filtre anti-repliement du numériseur atténue le signal de plus de 80 % (ou, si l’inventaire ne le décrit pas, au-delà de 0,8 fois la fréquence de Nyquist) ; les lignes dont le filtre déborde de cette limite sont hachurées. Amplitudes en vitesse du sol, en décibels par rapport au maximum de la station ; sous −40 dB, la couleur du fond.

Échelle de couleur, au choix : « station », 0 dB au maximum de la station choisie (verticale et transversale ensemble) ; « séisme », au maximum de toutes les stations tracées, pour comparer les couleurs d’une station à l’autre ; « absolue », 0 dB à 1 µm/s, de −40 à +40 dB (de 10 nm/s à 100 µm/s), la même sur toutes les fiches, saturée au-delà. Les niveaux ne sont pas recalculés : seule leur correspondance avec les couleurs change.

Au-dessus de chaque spectrogramme, le signal de sa composante, traité comme les sismogrammes : tendance retirée, extrémités adoucies sur 2 s, passe-bande de Butterworth d’ordre 4 appliqué dans les deux sens entre 0,5 Hz et la limite de la station, ramené à son maximum dans la fenêtre affichée ; ce maximum, écrit à droite, est le signal divisé par la sensibilité de l’instrument, sans correction de sa réponse en fréquence. La fenêtre de temps est commune aux sismogrammes et au spectrogramme.

Transversale : −E·cos β + N·sin β, β l’azimut arrière, soit la direction de la station vers l’épicentre (convention d’ObsPy), calculé sur l’ellipsoïde depuis la position de l’instrument dans l’inventaire d’ORFEUS ; les deux horizontales sont ramenées en vitesse par leur sensibilité avant la rotation, puis corrigées de leur réponse commune. Elle n’est calculée que si l’incertitude sur la position de l’épicentre laisse cette direction connue à 10° près, atan(incertitude / distance) : au-delà, plus de 17 % de la composante radiale, où arrivent l’onde P et une partie de l’onde S, passerait dans la transversale. Les deux composantes partagent le même 0 dB.

Calcul contrôlé par une seconde implémentation indépendante (ObsPy et NumPy) sur des enregistrements du catalogue : réponse de l’instrument comparée à celle d’ObsPy (evalresp), enveloppes et niveaux de couleur à la NumPy. Réponse comparée à celle d’ObsPy (evalresp), filtres numériques compris, jusqu’à la limite de chaque station. Signaux : archive publique d’ORFEUS (EIDA), réseaux BE, NL.