Seismic Compression-Quiescence Anomalies as Earthquake Precursors: Prospective Prediction for the Kashmir Himalayan Fault Belt and Validation Across Four Tectonic Regimes

1. Definição do Problema

O arco de Caxemira-Himalaia representa uma das regiões sismicamente mais perigosas do mundo, impulsionada pela convergência norte da placa Indiana e da placa Euroasiática. Essa convergência armazena uma energia de deformação elástica significativa ao longo de sistemas de falhas paralelos (Falha Frontal Principal, Falha de Limite Principal, Falha de Panjal e Falha Central Principal). Apesar de catástrofes históricas como o terremoto de M 7,6 em Caxemira em 2005, os métodos atuais para a previsão de terremotos nesta região enfrentam três limitações específicas:

• Falta de um limiar quantitativo único válido para diferentes configurações tectônicas.
• Ausência de registros de previsão prospectivamente verificáveis especificamente para Caxemira.
• Incapacidade de atribuir sinais de carregamento a segmentos de falhas específicos, em vez de amplas regiões geográficas.

Métodos existentes, como a quiescência sísmica (estatística Z), o mapeamento do valor-b e a liberação acelerada de momento, não alcançaram simultaneamente esses objetivos.

2. Metodologia

O estudo introduz um framework baseado em dados centrado no Índice de Compressão de Magnitude ($C$) e no sinal de queda de C (C-drop), validado em quatro regimes tectônicos distintos (Caxemira, Japão, Turquia, Itália).

Métricas Principais:

• Índice de Compressão de Magnitude ($C$): Uma medida livre de parâmetros definida como $C(k) = \max(0, 1 - \sigma[W] / \sigma_0)$. Quantifica o estreitamento das distribuições de magnitude dos terremotos à medida que uma falha acumula tensão. Conforme o carregamento de tensão homogeneiza a sismicidade local, o desvio padrão ($\sigma$) das magnitudes diminui, fazendo com que $C$ suba em direção a 1.
• A Queda de C (C-drop): O principal sinal precursor. É definida como um declínio sistemático e sustentado pós-pico em $C$ (pelo menos 0,050 a partir de um pico $\ge$ 0,25). Esse declínio ocorre quando a falha começa a falhar, permitindo que terremotos maiores entrem na janela de observação e aumentem a variância da magnitude.
• Índice de Quiescência Sísmica ($Q$): Uma medida secundária de redução de taxa ($Q = \max(0, 1 - \lambda[W] / \lambda_0$). Um $Q$ elevado combinado com um $C$ elevado indica uma falha criticamente travada.

Processamento de Dados e Identificação de Zonas:

• Janela Móvel: Os cálculos utilizam uma janela móvel de $W=20$ terremotos mais recentes em uma zona, suavizada por uma média móvel de 12 eventos.
• Linha de Base: Valores de referência ($\sigma_0, \lambda_0$) são derivados do período 2000–2009.
• Zoneamento Baseado em Dados: Diferente de estudos anteriores que utilizam coordenadas arbitrárias, os centros das zonas são posicionados em picos naturais de densidade de sismicidade identificados via suavização Gaussiana do catálogo ISC (2000–2026). Os raios das zonas capturam 90% dos clusters de sismicidade local.
• Sistema de Pontuação de Seis Critérios: Para filtrar falsos positivos (ex: sequências de réplicas), uma zona é classificada como CRÍTICA apenas se atender a pelo menos 3 de 6 critérios simultaneamente:
  1. Compressão: $C \ge 0,25$.
  2. Quiescência: $Q \ge 0,20$.
  3. Queda de C (C-drop): Declínio confirmado a partir do pico.
  4. Valor-b: $b < 0,75$ (indicando concentração de tensão).
     \ 5. Baixa Atividade: $<15$ eventos nos últimos 90 dias.
  5. Frescor (Freshness): O pico ocorreu nos últimos 180 dias.

3. Principais Contribuições

• O Sinal de Queda de C (C-drop): Estabelece o declínio pós-pico na compressão de magnitude como um precursor primário e fisicamente fundamentado para terremotos significativos iminentes, distinto da simples quiescência baseada em taxa.
• Decomposição de Zonas Baseada em Dados: Demonstra que zonas sísmicas compostas podem ser decompostas para identificar segmentos de falhas ativos específicos. Por exemplo, o estudo isola o sinal de carregamento ativo na região de Anantnag para a Falha de Panjal (Ramban–Banihal), em vez da descarga do núcleo da Falha Central Principal.
• Verificação Prospectiva: O artigo apresenta uma previsão formal arquivada no Zenodo (18 de junho de 2026) antes da ocorrência de eventos ativadores, fornecendo um registro verificável da aplicação do método.
• Validação Trans-Tectônica: O método é testado retrospectivamente em quatro tipos de limites de placas (colisão, subducção, deslizamento lateral e arco compressivo), demonstrando a universalidade do mecanismo de C-drop.

4. Resultados

Validação Retrospectiva:
Em 53 ciclos confirmados de C-drop em quatro regimes tectônicos, um terremoto de magnitude alvo ocorreu dentro de 90 dias em 83,0% dos casos (44/53).

• Japão (Subducção): 88,9% de taxa de detecção.
• Turquia (Transformante): 91,7% de taxa de detecção.
• Itália (Compressão): 75,0% de taxa de detecção.
• Índia/Caxemira (Colisão): 73,3% de taxa de detecção.
• Subconjunto de Alta Confiança: Para ciclos com $C_{pico} \ge 0,30$ e queda $\ge 0,055$, a taxa de detecção atingiu 100% em todas as regiões.

Status da Zona de Caxemira (até 20 de junho de 2026):

• Zona Z21 (Anantnag): Pontuação 5/6. $C_{drop} = 0,129$. Dois eventos M 4,8 ocorreram em 11 de abril e 5 de junho de 2026, durante a fase ativa de C-drop. O ciclo permanece aberto com $C=0,260$ (67% do carregamento de pico retido).
• Zona Z15b (Caxemira Norte): Pontuação 5/6. $C_{drop} = 0,092$. Maior carregamento em 26 anos ($C_{pico}=0,585$). A taxa de detecção histórica para M $\ge$ 4,5 é de 96,3%. Previsto M 5,0–6,0 para julho–setembro de 2026.
• Zona Z15 (Muzaffarabad): Pontuação 4/6. $C=0,342$ com quiescência excepcionalmente profunda ($Q=0,795$). A confirmação do C-drop (limiar de 0,050) é esperada dentro de 2–4 semanas.
• Zona FZ-NPB (Gilgit): Pontuação 4/6. Maior valor de $C$ do estudo (0,530). O C-drop acabou de começar; uma previsão emergente foi sinalizada para agosto–outubro de 2026.

Eventos de Validação:
Três eventos validaram a janela de previsão prospectiva:

1. M 3,3 (19 de junho de 2026): Ocorreu 18 horas após o arquivo no Zenodo, quebrando um silêncio de 13 dias na Zona Z21.
2. M 4,8 (11 de abril de 2026) & M 4,8 (5 de junho de 2026): Ambos ocorreram durante a fase ativa de C-drop da Zona Z21, representando liberações parciais de tensão.

5. Significância e Alegações

O artigo afirma fornecer uma previsão de terremoto prospectiva e prospectivamente verificável para o cinturão de falhas de Caxemira-Himalaia, abordando uma lacuna crítica na prática sismológica atual.

• Mecanismo Físico: Os autores argumentam que o C-drop reflete um processo universal de mecânica de falhas onde a variância da magnitude diminui sob carga e aumenta na falha, independentemente do cenário tectônico.
• Precisão: Ao utilizar a decomposição de zonas baseada em dados, o estudo refina a atribuição espacial da tensão, identificando a Falha de Panjal como a fonte ativa para o sinal de Anantnag, em vez da anteriormente assumida Falha Central Principal.
• Avaliação de Risco: O estudo destaca que aproximadamente 10 milhões de pessoas vivem nas zonas afetadas entre Jammu e Caxemira, Caxemira administrada pelo Paquistão e Gilgit-Baltistan.
• Chamada para Ação: Os autores convidam a comunidade sismológica internacional (incluindo NCS e CSEP) para monitorar e avaliar independentemente a previsão, enfatizando que o método é autossuficiente para distinguir entre estados de falha carregados (predictíveis) e descarregados (não previsíveis).

O artigo conclui que o método C-drop, combinado com o zoneamento baseado em dados e um filtro de seis critérios, oferece um framework robusto para identificar estados críticos de carregamento de falhas e emitir previsões prospectivas com alta confiança estatística.