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

Resumo Técnico: Análise Adaptativa Sazonal de Dados de Radônio e Correlação com Eventos Sísmicos Regionais

Definição do Problema
O estudo aborda o desafio de identificar anomalias geoquímicas transitórias nas emissões atmosféricas de radônio que possam servir como precursores de eventos sísmicos regionais. Uma dificuldade primária neste domínio é distinguir sinais tectônicos genuínos de flutuações ambientais multiescala, incluindo ruído meteorológico diurno (temperatura, umidade) e fatores sinóticos de médio prazo (pressão barométrica, umidade do solo). O artigo postula que, sem estruturas estatísticas robustas para filtrar esses artefatos, a correlação entre a desgaseificação cortical e a aceleração sísmica permanece obscurecida.

Metodologia
A pesquisa implementa uma estrutura de "Análise Adaptativa Sazonal" utilizando uma Linha de Base de Desvio Padrão (DP) Global.

• Aquisição de Dados: O monitoramento contínuo de radônio foi conduzido a uma altitude de 540m (Latitude 35,30°N, Longitude 25,64°E) utilizando um detector de fotodiodo de silício comercial da Airthings. O sensor, protegido dentro de um invólucro especializado para evitar a entrada de umidade e interferência direta de radiação ambiente, operou via espectrometria alfa.
• Processamento de Dados: Para suprimir o ruído diurno, os dados de séries temporais brutos de alta fidelidade foram agregados em valores médios diários.
• Estrutura Estatística:
  • Um Desvio Padrão Global ($\sigma_{global}$) foi computado ao longo de toda a série temporal de fundo para estabelecer um parâmetro estatístico que filtre as modulações meteorológicas de médio prazo.
  • Um Limiar de Alerta Adaptativo foi definido como: $Threshold = \mu_{6-day} + 2,5 \times \sigma_{global}$.
  • Esta configuração utiliza uma média dos últimos 6 dias ($\mu_{6-day}$) para rastrear tendências de curto prazo, enquanto o DP global garante a sensibilidade a eventos de desgaseificação cortical rápidos e induzidos por estresse, minimizando falsos positivos provenientes de ruído ambiente localizado.
• Correlação: O estudo integra esta série temporal geoquímica com dados de atividade sísmica regional (Magnitude $M_L$) ao longo do período de 24 de junho a 28 de julho de 2026.

Principais Resultados
A análise do conjunto de dados de junho–julho de 2026 gerou as seguintes observações:

• Parâmetros Estatísticos: Ao longo de 35 dias, o valor médio de radônio foi de 20,3 Bq/m³ com um desvio padrão de 11,3, variando de um mínimo de 7 Bq/m³ a um máximo de 56 Bq/m³.
• Anomalia Validada (Pré-sísmica): Uma anomalia significativa foi detectada em 29 de junho de 2026, atingindo um desvio de aproximadamente 2,9 DP. Este evento ocorreu cerca de 72 horas antes de um grande terremoto de M 5,2 em 2 de julho de 2026. O estudo interpreta este evento como um precursor geoquímico de curto prazo impulsionado pela aceleração da deformação cortical crítica e microfraturamento.
• Resposta Pós-sísmica: Um pico distinto abaixo do limiar foi identificado em 3 de julho após o evento de M 5,2, atribuído ao relaxamento de estresse pós-ruptura e distúrbios hidrogeológicos transitórios que forçam o radônio profundo para a atmosfera.
• Precursores Sub-limiar: Dois picos intermediários foram identificados em 9 de julho (~37 Bq/m³) e 11 de julho (~33 Bq/m³). Embora não tenham ultrapassado o rigoroso limiar de +2,5 DP, sua proximidade temporal com um segundo cluster sísmico (M 3,8 em 16 de julho e M 4,4 em 18 de julho) sugere uma acumulação em degraus do estresse tectônico regional.
• Alerta Atual: Uma anomalia crítica de 2,9 DP foi registrada em 28 de julho de 2026. Os autores classificam esta como um alerta não validado, porém de alta prioridade, observando que seu comportamento estatístico espelha o precursor validado de 29 de junho.

Significância e Alegações
O artigo afirma demonstrar a viabilidade do uso de tecnologia de fotodiodo de silício de baixo custo e comercialmente disponível (especificamente sensores Airthings) para o monitoramento de alta resolução de precursores de terremotos quando pareada com processamento estatístico robusto e blindagem ambiental.

• Desacoplamento de Ruído: A integração de uma linha de base meteorológica derivada de Random Forest com uma média independente dos últimos 6 dias desacopla com sucesso o ruído meteorológico de sinais tectônicos genuínos.
• Validação do Mecanismo: O alinhamento temporal entre os picos de desgaseificação de radônio e as magnitudes dos terremotos regionais apoia a premissa de que o acúmulo de estresse tectônico localizado modula os mecanismos de transporte de gases subsuperficiais antes da ruptura mecânica.
• Escalabilidade: O estudo afirma que esta abordagem valida o potencial de implantação de redes de alerta regionais densas e de baixo custo, oferecendo capacidades de resolução de anomalias comparáveis à instrumentação científica dispendiosa.
• Aplicação Imediata: Os autores concluem que a anomalia de 28 de julho coloca a região de Creta em uma "janela de monitoramento altamente crítica", necessitando de observação em tempo real e avaliação geofísica urgente como um potencial sinal precursor de curto prazo.