Seasonal Adaptive Analysis of Radon Data and Correlation with Regional Earthquake Events

Résumé Technique : Analyse Adaptative Saisonnière des Données de Radon et Corrélation avec les Événements Sismiques Régionaux

Énoncé du Problème
L'étude traite du défi de l'identification des anomalies géochimiques transitoires des émissions de radon atmosphérique pouvant servir de précurseurs aux événements sismiques régionaux. Une difficulté majeure dans ce domaine est de distinguer les signaux tectoniques authentiques des fluctuations environnementales multi-échelles, notamment le bruit météorologique diurne (température, humidité) et les facteurs synoptiques à moyen terme (pression barométrique, humidité du sol). L'article postule que sans cadres statistiques robustes pour filtrer ces artefacts, la corrélation entre le dégazage crustal et l'accélération sismique reste obscurcie.

Méthodologie
La recherche met en œuvre un cadre d'« Analyse Adaptative Saisonnière » utilisant une Base de Référence d'Écart Type (SD) Globale.

• Acquisition des Données : Une surveillance continue du radon a été menée à une altitude de 540 m (Latitude 35,30°N, Longitude 25,64°E) à l'aide d'un détecteur à photodiodes de silicium commercial de type Airthings. Le capteur, protégé dans un boîtier spécialisé pour empêcher l'infiltration d'humidité et les interférences de rayonnement ambiant direct, fonctionnait par spectrométrie alpha.
• Traitement des Données : Pour supprimer le bruit diurne, les données brutes de séries temporelles à haute fidélité ont été agrégées en valeurs moyennes quotidiennes.
• Cadre Statistique :
  • Un Écart Type Global ($\sigma_{global}$) a été calculé sur l'ensemble de la série temporelle de fond pour établir un référentiel statistique qui filtre les modulations météorologiques à moyen terme.
  • Un Seuil d'Alerte Adaptatif a été défini comme suit : $Seuil = \mu_{6-jours} + 2,5 \times \sigma_{global}$.
  • Cette configuration utilise une moyenne des 6 jours passés ($\mu_{6-jours}$) pour suivre les tendances à court terme, tandis que l'écart type global assure la sensibilité aux événements de dégazage crustal rapides induits par le stress, tout en minimisant les faux positifs provenant du bruit ambiant localisé.
• Corrélation : L'étude intègre cette série temporelle géochimique avec les données d'activité sismique régionale (Magnitude $M_L$) sur la période du 24 juin au 28 juillet 2026.

Résultats Clés
L'analyse de l'ensemble de données de juin-juillet 2026 a produit les observations suivantes :

• Paramètres Statistiques : Sur 35 jours, la valeur moyenne du radon était de 20,3 Bq/m³ avec un écart type de 11,3, allant d'un minimum de 7 Bq/m³ à un maximum de 56 Bq/m³.
• Anomalie Validée (Pré-sismique) : Une anomalie significative a été détectée le 29 juin 2026, atteignant un écart d'environ 2,9 SD. Cet événement est survenu environ 72 heures avant un séisme majeur de M 5,2 le 2 juillet 2026. L'étude interprète cela comme un précurseur géochimique à court terme piloté par l'accélération critique de la déformation crustale et la micro-fracturation.
• Réponse Post-sismique : Un pic de sous-seuil distinct le 3 juillet a suivi l'événement M 5,2, attribué à la relaxation des contraintes post-rupture et aux perturbations hydrogéologiques transitoires forçant le radon profond vers l'atmosphère.
• Précurseurs de Sous-seuil : Deux pics intermédiaires ont été identifiés le 9 juillet (~37 Bq/m³) et le 11 juillet (~33 Bq/m³). Bien qu'ils n'aient pas franchi le seuil strict de +2,5 SD, leur proximité temporelle avec un second cluster sismique (M 3,8 le 16 juillet et M 4,4 le 18 juillet) suggère une accumulation par paliers du stress tectonique régional.
• Alerte Actuelle : Une anomalie critique de 2,9 SD a été enregistrée le 28 juillet 2026. Les auteurs la classent comme une alerte non validée mais de haute priorité, notant que son comportement statistique reflète l'anomalie validée du 29 juin.

Signification et Revendications
L'article prétend démontrer la faisabilité de l'utilisation de la technologie de photodiodes de silicium de faible coût et commercialement disponible (spécifiquement les capteurs Airthings) pour la surveillance des précurseurs de séismes à haute résolution, lorsqu'elle est couplée à un traitement statistique robuste et un blindage environnemental.

• Découplage du Bruit : L'intégration d'une base de référence météorologique dérivée de Forêts Aléatoires (Random Forest) avec une moyenne indépendante des 6 jours passés permet de découpler avec succès le bruit météorologique des signaux tectoniques authentiques.
• Validation du Mécanisme : L'alignement temporel entre les pics de dégazage de radon et les magnitudes sismiques régionales soutient la prémisse selon laquelle l'accumulation de stress tectonique localisé module les mécanismes de transport de gaz souterrains avant la rupture mécanique.
• Évolutivité : L'étude affirme que cette approche valide le potentiel de déploiement de réseaux de surveillance régionaux denses et rentables, offrant des capacités de résolution d'anomalies comparables à l'instrumentation scientifique coûteuse.
• Application Immédiate : Les auteurs concluent que l'anomalie du 28 juillet place la région de Crète élargie dans une « fenêtre de surveillance hautement critique », nécessitant une observation en temps réel et une évaluation géophysique urgente en tant que signal précurseur potentiel à court terme.