PaleoBAYES: novo modelo bayesiano calcula magnitudes de terremotos antigos no Apenino

Como analista de infraestrutura digital, observo que a reconstrução da história sísmica funciona como um diagnóstico dos alicerces geológicos que sustentam nossas cidades. Um grupo do INGV, em parceria com a SUSTech de Shenzhen, apresentou um avanço metodológico que transforma medições de campo em estimativas rigorosas: o software open‑source PaleoBAYES.

Publicado no Bulletin of the Seismological Society of America, o trabalho descreve um modelo probabilístico baseado em inferência bayesiana que, pela primeira vez, combina três medições geológicas fundamentais: o rigetto (o deslocamento observado na superfície), o comprimento da ruptura da crosta e a idade do evento. Essa integração é armazenada e testada no novo banco de dados PaleoEcA, que reúne 44 paleoterremotos documentados no Apenino central ao longo dos últimos 27.000 anos.

Do ponto de vista prático, a contribuição do trabalho é dupla. Primeiro, os autores harmonizaram décadas de levantamentos paleossismológicos, associando a cada medição sua incerteza — uma etapa crítica para que a estimativa final não pareça mais precisa do que realmente é. Como explica Deborah Di Naccio, pesquisadora do INGV e primeira autora do estudo, a transparência na atribuição de incertezas evita a ilusão de números definitivos em um campo onde o registro é fragmentário.

Em segundo lugar, PaleoBAYES modela dois fatores frequentemente negligenciados em estimativas anteriores: o efeito do tempo sobre a preservação das evidências e a distribuição estatística das magnitudes. Na prática, o algoritmo admite que traços mais antigos têm maior probabilidade de terem sido apagados pelos processos geológicos e que a ocorrência de terremotos obedece à lei de Gutenberg‑Richter, ou seja, eventos de maior magnitudo são progressivamente menos frequentes. Ao incluir esses elementos, o resultado não é um valor pontual, mas um intervalo de credibilidade: para os 44 eventos do banco, a incerteza ao nível de 95% costuma ficar entre ±0,3 e ±0,7 unidades de magnitudo.

O código foi validado em sismos históricamente documentados. Para o terremoto do Fucino de 1915, PaleoBAYES estimou M 6,83 ± 0,13 — coerente com as medições instrumentais e macrossísmicas conhecidas (M 6,7–7,0). Essa verificação demonstra que o modelo, ao tratar as medidas como camadas de informação com incerteza explícita, produz estimativas compatíveis com o registro instrumental quando disponível.

Uma implicação relevante do estudo é que os vestígios geológicos preservados tendem a amostrar preferencialmente eventos de maior magnitudo. Em linguagem de sistemas, as trincas e deslocamentos são como as vias principais do sistema nervoso das montanhas: deixam sinais mais persistentes que as perturbações menores.

PaleoBAYES e o banco PaleoEcA representam uma infraestrutura aberta para a paleossismologia — uma camada de inteligência que permite transformar dados fragmentários de campo em probabilidades úteis para modelagem de risco e planejamento. Para as cidades e redes críticas da Europa, essa ferramenta melhora a compreensão do passado sísmico e, portanto, a resiliência futura, sem recorrer a fantasias tecnológicas: trata‑se de integrar melhor os dados e explicitar as incertezas.

O código open‑source e o banco de dados estão prontos para serem usados por pesquisadores e gestores de risco que desejem incorporar estimativas de magnitudo históricas em modelos de vulnerabilidade e políticas públicas.