EUPEN (BE)

Séismes tectoniques et induits du catalogue de l’Observatoire, présumés compris ; les tirs de carrière, explosions et autres événements ne sont pas comptés. Distances épicentrales sur l’ellipsoïde WGS84.

Commune : limites communales de 2025 (Statbel), appliquées à tous les séismes, y compris ceux antérieurs à la fusion des communes. Pour un épicentre situé hors de Belgique, seule la zone circulaire est donnée. Le nom de chaque séisme est celui du catalogue : il peut désigner une localité voisine de l’épicentre.

Magnitude : celle que donne la fiche de chaque séisme (ML de préférence, sinon MW ou MS) ; des magnitudes de types différents sont comparées telles quelles.

Les épicentres des séismes historiques ne sont connus qu’à quelques kilomètres près, parfois davantage : leur appartenance à la zone est indicative. Les petits séismes n’ont été enregistrés qu’à mesure que le réseau de stations s’est étoffé : les séismes anciens de la zone sont surtout les plus forts, et l’activité récente paraît plus fréquente.

Méthode des filtres multiples (Dziewonski, Bloch & Landisman, 1969), appliquée dans votre navigateur au signal brut de la station (100 é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.