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

Resumen Técnico: Análisis Adaptativo Estacional de Datos de Radón y Correlación con Eventos Sísmicos Regionales

Planteamiento del Problema
El estudio aborda el desafío de identificar anomalías geoquímicas transitorias en las emisiones atmosféricas de radón que puedan servir como precursores de eventos sísmicos regionales. Una dificultad primordial en este dominio es distinguir las señales tectónicas genuinas de las fluctuaciones ambientales multiescala, incluyendo el ruido meteorológico diarío (temperatura, humedad) y los factores sinópticos de mediano plazo (presión barométrica, humedad del suelo). El artículo postula que, sin marcos estadísticos robustos para filtrar estos artefactos, la correlación entre la desgasificación cortical y la aceleración sísmica permanece oscurecida.

Metodología
La investigación implementa un marco de "Análisis Adaptativo Estacional" utilizando una Línea Base de Desviación Estándar (DE) Global.

• Adquisición de Datos: El monitoreo continuo de radón se llevó a cabo a una altitud de 540 m (Latitud 35.30°N, Longitud 25.64°E) utilizando un detector de fotodiodo de silicio comercial Airthings. El sensor, protegido dentro de un recinto especializado para evitar la entrada de humedad e interferencias de radiación ambiental directa, operó mediante espectrometría alfa.
• Procesamiento de Datos: Para suprimir el ruido diarío, los datos de series temporales de alta fidelidad se agregaron en valores medios diarios.
• Marco Estadístico:
  • Se computó una Desviación Estándar Global ($\sigma_{global}$) a través de toda la serie temporal de fondo para establecer un punto de referencia estadístico que filtre las modulaciones meteorológicas de mediano plazo.
  • Se definió un Umbral de Alerta Adaptativo como: $Umbral = \mu_{6-días} + 2.5 \times \sigma_{global}$.
  • Esta configuración utiliza una media de los 6 días pasados ($\mu_{6-días}$) para rastrear tendencias a corto plazo, mientras que la DE global asegura la sensibilidad ante eventos de desgasificación cortical rápidos inducidos por estrés, minimizando al mismo tiempo los falsos positivos provenientes del ruido ambiental localizado.
• Correlación: El estudio integra esta serie temporal geoquímica con datos de actividad sísmica regional (Magnitud $M_L$) durante el periodo del 24 de junio al 28 de julio de 2026.

Resultados Clave
El análisis del conjunto de datos de junio–julio de 2026 arrojó las siguientes observaciones:

• Parámetros Estadísticos: Durante 35 días, el valor medio de radón fue de 20.3 Bq/m³ con una desviación estándar de 11.3, variando desde un mínimo de 7 Bq/m³ hasta un máximo de 56 Bq/m³.
• Anomalía Validada (Pre-sísmica): Se detectó una anomalía significativa el 29 de junio de 2026, alcanzando una desviación de aproximadamente 2.9 DE. Este evento ocurrió aproximadamente 72 horas antes de un terremoto mayor de M 5.2 el 2 de julio de 2026. El estudio interpreta esto como un precursor geoquímico de corto plazo impulsado por la aceleración de la deformación cortical crítica y la microfracturación.
• Respuesta Post-sísmica: Un pico distintivo por debajo del umbral el 3 de julio siguió al evento M 5.2, atribuido a la relajación del estrés post-ruptura y a perturbaciones hidrogeológicas transitorias que fuerzan el radón de origen profundo hacia la atmósfera.
• Precursores Sub-umbral: Se identificaron dos picos intermedios el 9 de julio (~37 Bq/m³) y el 11 de julio (~33 Bq/m³). Aunque no superaron el estricto umbral de +2.5 DE, su proximidad temporal a un clúster sísmico secundario (M 3.8 el 16 de julio y M 4.4 el 18 de julio) sugiere una acumulación escalonada de estrés tectónico regional.
• Alerta Actual: Se registró una anomalía crítica de 2.9 DE el 28 de julio de 2026. Los autores clasifican esto como una alerta de alta prioridad pero no validada, señalando que su comportamiento estadístico refleja el precursor validado del 29 de junio.

Significancia y Reivindicaciones
El artículo afirma demostrar la viabilidad de utilizar tecnología de fotodiodo de silicio de bajo costo y comercialmente disponible (específicamente sensores Airthings) para el monitoreo de alta resolución de precursores de terremotos cuando se combina con un procesamiento estadístico robusto y blindaje ambiental.

• Desacoplamiento de Ruido: La integración de una línea base meteorológica derivada de Random Forest con una media independiente de 6 días pasados desacopla con éxito el ruido meteorológico de las señales tectónicas genuinas.
• Validación del Mecanismo: La alineación temporal entre los picos de desgasificación de radón y las magnitudes de los terremotos regionales respalda la premisa de que la acumulación de estrés tectónico localizado modula los mecanismos de transporte de gases subsuperficiales antes de la ruptura mecánica.
• Escalabilidad: El estudio afirma que este enfoque valida el potencial de desplegar redes de alerta regional densas y rentables, ofreciendo capacidades de resolución de anomalías comparables a la instrumentación científica costosa.
• Aplicación Inmediata: Los autores concluyen que la anomalía del 28 de julio sitúa a la región de Creta en una "ventana de monitoreo altamente crítica", lo que requiere observación en tiempo real y una evaluación geofísica urgente como una posible señal precursora a corto plazo.