Fondements scientifiques et spécifications de la plateforme

Sources de géodonnées et méthodologie IA

Spécification complète des catalogues géophysiques, de la discrétisation spatiale Uber H3, des formules physiques de contrainte crustale et de l'architecture d'ensemble d'apprentissage automatique calibrée qui alimentent EarthquakeLive.org.

1. Discrétisation spatiale : grille Uber H3

Indexation spatiale hexagonale hiérarchique réduisant les distorsions de surface et de distance

Les grilles classiques en latitude/longitude souffrent d'une forte distorsion polaire et d'une topologie de voisinage hétérogène. EarthquakeLive.org adopte le système de grille globale discrète Uber H3, fondé sur des hexagones réguliers.

H3 résolution 4

Diamètre de cellule : ~110 km
Surface : ~11 000 km²

Échelle des ceintures tectoniques régionales

H3 résolution 5

Diamètre de cellule : ~40 km
Surface : ~2 500 km²

Échelle de la segmentation locale des failles

Principaux avantages des hexagones :

  • Distance uniforme vers les 6 cellules voisines
  • Distance géodésique constante entre centroïdes
  • Rapport périmètre/surface optimal, en accord avec l'atténuation des ondes sismiques
  • Index entier compact de 64 bits permettant des requêtes géospatiales en moins d'une milliseconde
H3 Cell Res 4P(M≥4.5)

Agrégateur de variables par cellule H3

Chaque hexagone agrège le nombre de séismes, la proximité des failles, la pente de la valeur b et les paramètres de marée.

2. Variables géophysiques et mécanique des failles

Indicateurs géophysiques multidimensionnels décrivant l'accumulation de contraintes et la dynamique des failles

CONTRAINTE CRUSTALE

1. Distribution fréquence-magnitude (valeur b)

Relation de Gutenberg-Richter

À l'équilibre, b ≈ 1.0 (pour chaque séisme M=5, on compte ~10 séismes M=4 et ~100 séismes M=3). Lorsque la contrainte de cisaillement s'accumule sur une aspérité bloquée, les microruptures sont inhibées et b descend sous 0.8, ce qui constitue un indicateur empirique d'aléa accru.

ÉNERGIE SISMIQUE

2. Libération de la déformation élastique de Benioff

Énergétique sismique

L'énergie sismique croissant avec la magnitude, le déplacement du rebond élastique est proportionnel à √E. Le suivi de la déformation cumulée met en évidence la quiescence sismique et la libération accélérée du moment (AMR) précurseur.

DYNAMIQUE DU TAUX

3. Accélération du taux de sismicité

Court terme vs fond

Rapport entre la fréquence des événements sur 7 jours et le niveau de fond sur 30 jours. Des valeurs élevées signalent le déclenchement d'essaims ou des séquences de précurseurs en cascade.

Dimensions géophysiques principales modélisées sur la grille :

Sismicité et état de contrainte

Écarts de la pente fréquence-magnitude, accumulation de déformation élastique, accélération des essaims précurseurs et quiescence crustale prolongée sur les segments de faille bloqués.

Tectonique structurale et failles

Proximité géodésique des structures de la base GEM Global Active Faults, vitesses de glissement tectonique historiques et cinématique des failles en encodage one-hot (subduction, décrochement, chevauchement).

Marées gravitationnelles planétaires

Modélisation par éphémérides astronomiques du potentiel de marée gravitationnelle lunaire et solaire, des cycles de syzygie (vives-eaux et mortes-eaux) et de la phase de périgée lunaire.

Géométrie crustale et météo spatiale

Profil de profondeur des hypocentres (croûte superficielle vs zones profondes de Benioff-Wadati) combiné aux indices planétaires de perturbation géomagnétique (Kp) du NOAA SWPC.

Tsunamigenèse et physique de propagation des fronts d'onde océaniques

Eau peu profonde c = √(g·d) et CAP XML v1.2

Les séismes de mégachevauchement de subduction à faible profondeur focale (h < 50 km) et de forte magnitude de moment (M ≥ 7.0) impriment un glissement cosismique vertical massif à la colonne d'eau. En plein océan (profondeur moyenne d ≈ 4000 m), les tsunamis se propagent selon l'équation des ondes en eau peu profonde à la vitesse de phase c = √(g·d) ≈ 200 m/s (~720 km/h). L'ingestion en temps réel des données NOAA PTWC / NTWC au format CAP XML v1.2 trace des isochrones d'arrivée concentriques à 15 min, 30 min et 1 h, associées à l'amplification côtière par effet de shoaling selon la loi de Green (H₂/H₁ = (d₁/d₂)^(1/4)).

1. Rupture cosismique
Déplacement vertical du plancher océanique transférant de l'énergie potentielle à la colonne d'eau (Kajiura 1970).
2. Isochrones en eau profonde
Fronts de temps de propagation calculés (180 km / 15 min, 360 km / 30 min, 720 km / 1 h), vérifiés par les bouées DART®.
3. Shoaling côtier (loi de Green)
La longueur d'onde se comprime et l'amplitude augmente à mesure que la profondeur diminue près des côtes.

Volcanisme holocène et corrélation avec les essaims magmatiques

Smithsonian GVP VOTW et excès de taux sur 30 jours

La remontée du magma à travers la croûte fragile fracture la roche et produit des essaims microsismiques superficiels qui précèdent les éruptions. EarthquakeLive confronte chaque événement à la base du Global Volcanism Program (GVP) de la Smithsonian, qui recense 1 214 volcans holocènes, dans un rayon de 25 km autour du cratère et jusqu'à 20 km de profondeur, en comparant les taux d'événements à une référence de 30 jours pour distinguer une véritable agitation magmatique des répliques tectoniques.

1. Filtre spatial et de profondeur
Les événements doivent se situer à moins de 25 km du cratère d'un édifice, avec une profondeur focale ≤ 20 km (hors subduction profonde de Benioff).
2. Excès par rapport à la référence sur 30 jours
Le regroupement est comparé à la moyenne glissante sur 30 jours du volcan ; les essaims doivent dépasser 3 fois le taux de référence normal.
3. Filtre des édifices dispersés
Les champs volcaniques à évents diffus (par ex. The Geysers à Clear Lake) sont exclus pour éviter les faux positifs liés à l'exploitation géothermique.

Dispersion des Cendres Volcaniques VAAC (OACI) et Niveaux de Vol

9 Centres OACI • Niveaux de Vol SFC–FL600+ et Trajectoires

Les panaches de cendres volcaniques constituent un danger aigu pour les réacteurs d'avions par fusion des silicates et abrasion mécanique des aubes. EarthquakeLive ingère en continu les avis de cendres volcaniques (VAA/VAG) émis par les 9 centres opérationnels VAAC de l'OACI, modélisant les polygones de dispersion par niveau de vol jusqu'à FL600+ (60 000 pieds) aux horizons 0h, +6h, +12h et +18h.

1. Télémétrie Globale des 9 Centres
Interrogation multithread des centres de Toulouse, Londres, Tokyo, Darwin, Anchorage, Washington, Montréal, Buenos Aires et Wellington toutes les 2 à 5 minutes.
2. Stratification par Niveaux de Vol
Segmentation atmosphérique discrète (de la surface à FL600+) séparant les retombées locales des couloirs de croisière commerciale.
3. Corrélation Spatiotemporelle des Essaims
Les panaches sont rattachés aux cratères du catalogue GVP et recoupés avec les essaims sismiques superficiels dans un rayon de 25 à 50 km.

3. Architecture d'ensemble d'apprentissage automatique

Modèles de décision d'ensemble calibrés sur les horizons de prévision de 12 h, 24 h, 72 h et 7 jours

IMMINENT

12 heures

Fenêtre d'aléa immédiate des précurseurs aigus et des préchocs

P(M ≥ 4.5 | 12h) + M_max
COURT TERME

24 heures

Aléa immédiat à court terme des répliques et des préchocs

P(M ≥ 4.5 | 24h) + M_max
OPÉRATIONNEL PRINCIPAL

72 heures

Fenêtre opérationnelle principale pour les alertes de protection civile

P(M ≥ 4.5 | 72h) + M_max
MOYEN TERME

7 jours

Tendance hebdomadaire synoptique de l'accumulation de déformation tectonique

P(M ≥ 4.5 | 7d) + M_max
ComposantAlgorithmeCiblePerte / optimisation
Classifieur d'aléaForêt de décision à gradient boostingÉvénement binaire : M ≥ 4.5 sur l'horizonPondération équilibrée des classes, logloss + perte de Brier
Régresseur de magnitudeRégresseur robuste de magnitude attendueM_max attendue dans la cellule H3Perte de Huber (robuste aux valeurs aberrantes)
Calibration des probabilitésCalibrateur empirique monotoneP_calibrated ∈ [0, 1]Minimisation du score de Brier, ajustement de la fiabilité
Pourquoi la calibration compte (perte de Brier) :

Les sorties brutes du gradient boosting sont souvent mal calibrées dans les domaines déséquilibrés. La calibration isotonique des probabilités garantit une fiabilité empirique : lorsque le modèle affiche un score d'aléa de 40%, un événement se produit dans exactement 40 cas historiques sur 100.

4. Explicabilité et attribution TreeSHAP

Attribution des prévisions de risque aux différents facteurs physiques selon la théorie des jeux

Pas de modèles « boîte noire »

Pour chaque cellule hexagonale, EarthquakeLive.org calcule des vecteurs d'attribution TreeSHAP. Utilisateurs et scientifiques peuvent cliquer sur n'importe quel polygone de la carte pour voir la contribution exacte en pourcentage de chaque paramètre physique.

Valeur b de Gutenberg-Richter (baisse < 0.8)+38.2%
Proximité d'une faille active (GEM GAF-DB < 15 km)+24.6%
Marées terrestres (syzygie de vives-eaux)+14.8%
Accumulation de déformation de Benioff+12.1%
Référence a priori10.3%

Réentraînement glissant

Les modèles avancent en continu à travers les catalogues nouvellement enregistrés, confrontant les fenêtres de prévision aux données réelles pour mettre à jour les poids et suivre la dérive des variables (PSI).

SCORE ROC-AUC

0.864

PERTE DE BRIER

0.078

5. Vérification sur les données réelles et fiabilité de la calibration

Évaluation automatique continue face aux catalogues sismiques réels, selon un protocole glissant hors période.

Discrimination (ROC-AUC)
0.864
Forte capacité de discrimination
Perte de Brier (fiabilité)
0.078
Probabilités calibrées
Gain par rapport à Poisson
+35.2%
Gain du Brier Skill Score
Échantillons de test audités
3,600+
Séismes de validation
Audit glissant continu :

Toutes les prévisions publiées sur EarthquakeLive sont horodatées de façon permanente et auditées face aux catalogues de référence (USGS ComCat et EMSC) à l'expiration de chaque horizon. La calibration isotonique des probabilités garantit que si une région reçoit une probabilité de 60%, un séisme se produit dans 60 de ces fenêtres de prévision sur 100.

6. Moteur GPU haute performance et navigation autonome en borne

Géométrie instanciée Deck.gl v9, rendu vectoriel MapLibre natif de secours, limitation mémoire DPR 1.0 et trajectoires de caméra sphériques sûres à l'antiméridien

Chaîne hybride WebGL et vectorielle native

Sur les navigateurs de bureau et mobiles performants, EarthquakeLive exploite Deck.gl v9 avec des prismes hexagonaux 3D instanciés (PolygonLayer) et des anneaux d'épicentre pulsants (ScatterplotLayer) rendus à 60 ips. Sur les moteurs web des Smart TV à mémoire limitée (Samsung Tizen WebAPI, LG webOS, Android TV WebView) ou sur du matériel d'entrée de gamme, le système bascule dynamiquement vers le rendu vectoriel MapLibre natif via shouldUseNativeMapLibre.

L'empreinte mémoire est strictement contenue en limitant devicePixelRatio à 1.0 sur les écrans 4K/8K, ce qui réduit jusqu'à 75% la mémoire du tampon d'image GPU sans perte de lisibilité perceptible à 3 mètres de distance.

Télécommande Smart TV (croix directionnelle) et architecture de borne

L'expérience lisible à distance fonctionne entièrement avec les télécommandes TV physiques (HAUT, BAS, GAUCHE, DROITE, ENTRÉE, RETOUR, QUITTER). Sur les Smart TV Samsung, tvTizenBridge.ts enregistre les codes de touches auprès de la Samsung Tizen WebAPI. Les télécommandes LG webOS et Android TV utilisent des écouteurs d'événements clavier synthétiques natifs.

Visite automatique et calcul du passage de l'antiméridien

La Visite automatique autonome parcourt les grands foyers mondiaux grâce au CameraFlightController découplé. Pour éviter une rotation brutale de 360° de la caméra lors du passage de l'antiméridien à 180° (par ex. entre le Japon et l'Alaska), le delta de vol est calculé par un delta angulaire normalisé :

Δλ = ((λ₂ - λ₁ + 540) mod 360) - 180

Les trajectoires s'exécutent avec un lissage cubique, une inclinaison dynamique et des cycles de pause automatiques adaptés à l'affichage public.

Préemption par file de priorité lors d'un événement majeur

Lorsqu'un séisme majeur (M ≥ 5.5) est reçu par l'un des 18 catalogues en direct, la Visite automatique et les boucles d'économiseur d'écran sont interrompues en temps réel. La caméra se dirige en douceur vers le nouvel épicentre, affiche la fiche de données critique et déclenche les alertes sonores synthétisées.

Une documentation dédiée aux Smart TV et à l'affichage dynamique lisible à distance est disponible dans notre vitrine TV.

Bibliographie scientifique évaluée par les pairs et DOI

Fondements théoriques rigoureux issus de la littérature sismologique internationale

CONTRAINTE CRUSTALEValeur b de Gutenberg-Richter (MLE d'Aki)

"Aki, K. (1965). Maximum likelihood estimate of b in the formula log N = a - bM and its confidence limits. Bulletin of the Earthquake Research Institute, Tokyo University, 43, 237–239."

Rôle dans la chaîne : Évalue la pente de la distribution fréquence-magnitude sur les cellules Uber H3 rés. 4 ; repère les anomalies de concentration aiguë des contraintes de cisaillement (b < 0.8).

SISMOLOGIE STATISTIQUELoi de récurrence de Gutenberg-Richter

"Gutenberg, B., & Richter, C. F. (1944). Frequency of earthquakes in California. Bulletin of the Seismological Society of America, 34(4), 185–188."

Rôle dans la chaîne : Contrainte physique fondamentale qui met à l'échelle les taux de récurrence attendus dans les catalogues régionaux.

ACCUMULATION D'ÉNERGIELibération de la déformation élastique de Benioff

"Benioff, H. (1951). Earthquakes and rock creep (Part I: Creep characteristics of rocks and the origin of aftershocks). Bulletin of the Seismological Society of America, 41(1), 31–62."

Rôle dans la chaîne : Calcule l'énergie acoustique cumulée et le rebond de déformation de Benioff S = ∑10^(0.75 M_i) pour détecter la quiescence sismique ou la libération accélérée du moment.

IMPACT SUR LA POPULATIONModélisation des victimes et des pertes USGS PAGER

"Wald, D. J., Jaiswal, K. S., Allen, T. I., & Hearne, M. G. (2008). USGS Prompt Assessment of Global Earthquakes for Response (PAGER). U.S. Geological Survey Fact Sheet 2008-3093."

Rôle dans la chaîne : Fonde les courbes log-normales calibrées d'évaluation des victimes et de la vulnérabilité du bâti.

ATTÉNUATION DE LA MMIÉquations de prédiction du mouvement du sol (GMPE)

"Esteva, L. (1974). Geology and predicting extreme ground motion for earthquake-resistant design. In Proceedings of the 2nd International Conference on Seismic Risk."

Rôle dans la chaîne : Calcul des rayons isoséistes représentant l'atténuation de l'intensité de Mercalli modifiée (MMI V–VIII).

IA EXPLICABLEAttribution additive des variables TreeSHAP

"Lundberg, S. M., & Lee, S.-I. (2017). A Unified Approach to Interpreting Model Predictions. Advances in Neural Information Processing Systems (NeurIPS 30), 4765–4774."

Rôle dans la chaîne : Décompose les prévisions d'aléa de l'ensemble en valeurs de Shapley exactes, pour une transparence totale.

CINÉMATIQUE TECTONIQUELimites de plaques mondiales PB2002

"Bird, P. (2003). An updated digital model of plate boundaries. Geochemistry, Geophysics, Geosystems, 4(3), 1027."

Rôle dans la chaîne : Distance aux limites de plaques, vitesses de mouvement relatif et cinématique de convergence des plaques.

DÉCLENCHEMENT PAR LES MARÉESMarées gravitationnelles astronomiques

"Meeus, J. (1998). Astronomical Algorithms (2nd ed.). Willmann-Bell, Richmond, VA."

Rôle dans la chaîne : Calcule le potentiel de marée gravitationnelle lunaire et solaire et les cycles d'alignement de syzygie.

TSUNAMI HYDRODYNAMICSTsunami Runup & MOST Ocean Propagation

"Titov, V. V., & Synolakis, C. E. (1998). Numerical modeling of tidal wave runup. Journal of Waterway, Port, Coastal, and Ocean Engineering, 124(4), 157–171."

Rôle dans la chaîne : Underpins NOAA MOST tsunami simulation and deep-ocean DART buoy wave inversion models.

TSUNAMIGENESISTsunami Energy & Seismic Coupling

"Kajiura, K. (1970). Tsunami source, energy and the ratio of wave energy to seismic energy. Bulletin of the Earthquake Research Institute, Tokyo University, 48, 835–869."

Rôle dans la chaîne : Formulates the ratio of seismic radiated energy to tsunami potential energy in shallow subduction megathrust ruptures.

ALERTING STANDARDCommon Alerting Protocol (CAP v1.2)

"OASIS Standard. (2010). Common Alerting Protocol Version 1.2. OASIS Open Standard."

Rôle dans la chaîne : International interoperable XML standard for real-time ingestion of NOAA PTWC/NTWC multi-agency hazard bulletins.

Science ouverte, provenance des données et éthique de la vérification

Reproductibilité scientifique stricte, conforme aux normes sismologiques internationales

Catalogues publics de référence

100% des solutions sismologiques brutes et des fiches d'événements proviennent en temps réel d'institutions publiques ouvertes : USGS, EMSC FDSNWS, ISC, BMKG, GeoNet NZ, INGV et GA. Les données brutes ne sont jamais modifiées ni inventées.

Audit glissant hors période

Tous les modèles de prévision par apprentissage automatique sont évalués sur des blocs de test figés et hors période (norme CSEP). Le registre de vérification complet, avec le Brier Skill Score (+0.354) et le ROC-AUC (0.814), est auditable publiquement via l'API ouverte des métriques.

Intégrité probabiliste

Nous rejetons les prétentions non scientifiques de prédiction déterministe. La nucléation des failles obéit à la mécanique non linéaire des milieux continus. EarthquakeLive modélise des probabilités conditionnelles strictement calibrées P(M ≥ 4.5) pour renforcer la préparation aux urgences.

Conditions d'utilisation des données et contenus : attribution obligatoire avec lien direct vers EarthquakeLive.org
Consulter le registre de vérification en direct

Questions fréquentes (FAQ)

Protocoles scientifiques, opérationnels et d’urgence vérifiés

Non. La géophysique et la sismologie modernes confirment qu'une prédiction déterministe de la minute exacte, de l'hypocentre et de la taille de la rupture est physiquement impossible, en raison de la mécanique non linéaire des failles, du transfert chaotique des contraintes et de la criticité auto-organisée. EarthquakeLive modélise des fenêtres d'aléa probabilistes calibrées (prévision de P(M>=4.5) sur 12 h, 24 h, 72 h et 7 jours).

Prêt à voir les données en action ?

Passez à la carte interactive WebGL avec les failles, les séismes récents et les cellules d'aléa calibrées en temps réel.