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

Technisches Resümee: Saisonale adaptive Analyse von Radon-Daten und Korrelation mit regionalen Erdbebenereignissen

Problemstellung
Die Studie befasst sich mit der Herausforderung, transiente geochemische Anomalien in atmosphärischen Radonemissionen zu identifizieren, die als Vorläufer regionaler seismischer Ereignisse dienen können. Eine primäre Schwierigkeit in diesem Bereich besteht darin, echte tektonische Signale von multiskaligen Umweltfluktuationen (einschließlich diurnaler meteorologischer Rauschfaktoren wie Temperatur und Luftfeuchtigkeit sowie mittelfristiger synoptischer Faktoren wie Luftdruck und Bodenfeuchtigkeit) zu unterscheiden. Die Arbeit postuliert, dass ohne robuste statistische Rahmenbedingungen zur Filterung dieser Artefakte die Korrelation zwischen krustaler Entgasung und seismischer Beschleunigung verborgen bleibt.

Methodik
Die Forschung implementiert ein Framework der „Saisonalen adaptiven Analyse“ unter Verwendung einer globalen Standardabweichung (SD)-Baseline.

• Datenerfassung: Die kontinuierliche Radonüberwachung erfolgte in einer Höhe von 540 m (Breitengrad 35,30°N, Längengrad 25,64°O) mittels eines kommerziellen Airthings Silizium-Photodioden-Detektors. Der Sensor, der in einem spezialisierten Gehäuse zum Schutz vor Feuchtigkeitseintritt und direkter Umgebungstrahlung untergebracht war, arbeitete mittels Alpha-Spektrometrie.
• Datenverarbeitung: Um das diurnale Rauschen zu unterdrücken, wurden die hochpräzisen Roh-Zeitreihendaten zu täglichen Mittelwerten aggregiert.
• Statistischer Rahmen:
  • Eine globale Standardabweichung ($\sigma_{global}$) wurde über den gesamten Zeitreihen-Hintergrund berechnet, um einen statistischen Benchmark zu etablieren, der mittelfristige meteorologische Modulationen herausfiltert.
  • Ein adaptiver Alarm-Schwellenwert wurde definiert als: $Threshold = \mu_{6-day} + 2,5 \times \sigma_{global}$.
  • Diese Konfiguration nutzt einen 6-Tage-Vergangenheitsmittelwert ($\mu_{6-day}$), um kurzfristige Trends zu verfolgen, während die globale SD sicherstellt, dass die Sensitivität für schnelle, stressbedingte krustale Entgasungsereignisse gewahrt bleibt und gleichzeitig Fehlalarme durch lokales Umgebungsrauschen minimiert werden.
• Korrelation: Die Studie integriert diese geochemische Zeitreihe mit regionalen seismischen Aktivitätsdaten (Magnitude $M_L$) über den Zeitraum vom 24. Juni bis zum 28. Juli 2026.

Kernergebnisse
Die Analyse des Datensatzes von Juni–Juli 2026 lieferte die folgenden Beobachtungen:

• Statistische Parameter: Über 35 Tage betrug der mittlere Radonwert 20,3 Bq/m³ mit einer Standardabweichung von 11,3, wobei die Werte von einem Minimum von 7 Bq/m³ bis zu einem Maximum von 56 Bq/m³ reichten.
• Validierte Anomalie (Präseismisch): Eine signifikante Anomalie wurde am 29. Juni 2026 detektiert, die eine Abweichung von etwa 2,9 SD erreichte. Dieses Ereignis trat etwa 72 Stunden vor einem schweren M 5,2 Erdbeben am 2. Juli 2026 auf. Die Studie interpretiert dies als kurzfristigen geochemischen Vorläufer, der durch kritische krustale Dehnungsbeschleunigung und Mikro-Frakturierung getrieben wurde.
• Postseismische Reaktion: Ein distinkter Sub-Schwellenwert-Peak am 3. Juli folgte dem M 5,2 Ereignis; dieser wird auf post-rupture Stressrelaxation und transiente hydrogeologische Störungen zurückgeführt, die tief liegendes Radon in die Atmosphäre forcieren.
• Sub-Schwellenwert-Vorläufer: Zwei intermediäre Peaks wurden am 9. Juli (~37 Bq/m³) und 11. Juli (~33 Bq/m³) identifiziert. Obwohl diese den strengen +2,5 SD-Schwellenwert nicht überschritten, deutet ihre zeitliche Nähe zu einem sekundären seismischen Cluster (M 3,8 am 16. Juli und M 4,4 am 18. Juli) auf eine stufenweise Akkumulation regionaler tektonischer Spannung hin.
• Aktueller Alarm: Eine kritische 2,9 SD Anomalie wurde am 28. Juli 2026 aufgezeichnet. Die Autoren klassifizieren dies als einen unvalidierten, aber hochprioritären Alarm und merken an, dass sich das statistische Verhalten des vorangegangenen Vorläufers vom 29. Juni spiegelt.

Bedeutung und Behauptungen
Die Arbeit behauptet, die Machbarkeit des Einsatzes kostengünstiger, kommerziell verfügbarer Silizium-Photodioden-Technologie (speziell Airthings-Sensoren) für die hochauflösende Überwachung von Erdbebenvorläufern nachzuweisen, sofern diese mit einer robusten statistischen Verarbeitung und Umgebungsschilderung kombiniert wird.

• Entkopplung des Rauschens: Die Integration eines mittels Random Forest abgeleiteten Wetter-Baselines mit einem unabhängigen 6-Tage-Vergangenheitsmittelwert entkoppelt das meteorologische Rauschen erfolgreich von echten tektonischen Signalen.
• Mechanismus-Validierung: Die zeitliche Übereinstimmung zwischen den Radon-Entgasungsspitzen und den regionalen Erdbebenmagnituden stützt die Prämisse, dass der lokale tektonische Spannungsaufbau die Mechanismen des unterirdischen Gastransports vor einer mechanischen Ruptur moduliert.
• Skalierbarkeit: Die Studie behauptet, dass dieser Ansatz das Potenzial für den Einsatz dichter, kosteneffizienter regionaler Warnnetzwerke validiert, die über Anomalie-Auflösungsfähigkeiten verfügen, die mit teurer wissenschaftlicher Instrumentierung vergleichbar sind.
• Unmittelbare Anwendung: Die Autoren kommen zu dem Schluss, dass die Anomalie vom 28. Juli die breitere Region Kreta in ein „hochkritisches Überwachungsfenster“ versetzt, was eine Echtzeitbeobachtung und eine dringende geophysikalische Evaluierung als potenzielles kurzfristiges Vorläufersignal erfordert.