📊 statistics

A Log-Gaussian Cox process framework for zone-free probabilistic seismic hazard assessment in Sumatra

Este artigo propõe uma estrutura de avaliação de perigo sísmico probabilística livre de zonas para Sumatra que integra um processo de Cox Log-Gaussiano com covariáveis tectônicas e um modelo de réplica ETAS, demonstrando que esta abordagem produz acelerações de projeto significativamente superiores aos códigos atuais da Indonésia ao eliminar fronteiras arbitrárias de zonas de fonte.

Autores originais: Irwan Endrayanto Aluicius, Nanang Susyanto, Wiwit Suryanto, Dwi Ertiningsih, Fajar Adi-Kusumo

Publicado 2026-09-03

📖 6 min de leitura🧠 Leitura aprofundada

CC BY 4.0

Autores originais: Irwan Endrayanto Aluicius, Nanang Susyanto, Wiwit Suryanto, Dwi Ertiningsih, Fajar Adi-Kusumo

Artigo original sob licença CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/). ✨ Esta é uma explicação gerada por IA do artigo abaixo. Não foi escrita nem endossada pelos autores. Para precisão técnica, consulte o artigo original. Ler aviso legal completo

Os terremotos estão entre as forças mais destrutivas do nosso planeta e, durante décadas, os engenheiros têm dependido de um método específico para projetar edifícios que possam sobreviver a eles. Este método, conhecido como análise probabilística de perigo sísmico, tenta prever quão forte será o tremor do solo em um local específico ao longo de um longo período. Tradicionalmente, os cientistas abordaram isso dividindo uma região em grandes zonas fixas, como desenhar caixas em um mapa. Dentro de cada caixa, eles assumem que o risco de terremoto é o mesmo em todos os lugares e que os terremotos ocorrem de forma aleatória e independente uns dos outros. No entanto, em locais onde placas tectônicas massivas colidem e se esmagam umas contra as outras, essas suposições frequentemente falham. A tensão que se acumula no subsolo não para abruptamente nas bordas de uma caixa desenhada; ela flui suavemente pelo cenário. Além disso, quando um terremoto gigante ocorre, ele não acontece isoladamente; ele desencadeia uma longa cadeia de réplicas menores e agrupadas que seguem um padrão diferente da simples aleatoriedade.

Em um novo estudo focado na ilha de Sumatra, na Indonésia, pesquisadores propuseram uma maneira de ir além dessas caixas rígidas. Em vez de forçar o risco de terremoto em zonas pré-definidas, eles usaram uma estrutura matemática flexível que permite que o mapa de risco emerja da realidade física do solo. Ao analisar décadas de registros de terremotos e combinar-nos com medições de quão firmemente as placas tectônicas estão travadas umas nas outras, a equipe criou um mapa de perigo sísmico contínuo e fluido. Esta abordagem revelou que o risco é muito mais concentrado em certas áreas do que se pensava anteriormente, particularmente em regiões que estiveram quietas por um longo tempo, mas que estão, na verdade, acumulando uma energia imensa. As descobertas sugerem que os códigos de construção atuais na região podem subestimar significativamente o tremor que poderia ocorrer, especialmente em áreas onde o solo está mais firmemente travado e esperando para deslizar.

Os pesquisadores aplicaram este novo método a Sumatra, uma ilha longa e estreita situada em uma das fronteiras de placas mais ativas e perigosas da Terra. Aqui, a enorme placa Indo-Australiana mergulha sob a placa de Sunda, criando uma zona de pressão intensa. A equipe começou com um catálogo de mais de 5.000 terremotos registrados entre 1998 e 2024. Para entender o verdadeiro risco de fundo, eles primeiro tiveram que separar os terremotos principais das réplicas que os seguem, de forma muito semelhante a distinguir o evento principal dos seus ecos. Uma vez isolados os choques principais, eles procuraram pelos fatores físicos que controlavam onde esses terremotos aconteciam. Eles descobriram dois impulsionadores principais: a distância até a Grande Falha de Sumatra, uma grande rachadura que percorre a ilha, e o "acoplamento intersísmico", uma medida de quão firmemente as duas placas tectônicas estão presas uma à outra em sua fronteira.

Usando essas medições físicas, a equipe construiu um modelo que trata o risco de terremoto como uma superfície invisível e suave, em vez de uma coleção de zonas separadas. Neste modelo, o risco é maior onde as placas estão mais firmemente travadas e onde a falha está mais próxima. Esta abordagem permitiu que eles vissem padrões que os métodos tradicionais perderam. Por exemplo, o modelo identificou claramente uma "lacuna sísmica" na costa das Ilhas Mentawai. Esta área não sofre uma ruptura massiva há mais de dois séculos, mas as placas ali estão tão firmemente travadas que o modelo prevê um risco muito alto de um futuro terremoto gigante. Como o modelo não depende da localização de terremotos passados para definir as zonas, ele pode prever alto risco em áreas que estiveram quietas por um longo tempo, desde que as condições físicas — como o travamento firme das placas — sejam adequadas.

Quando os pesquisadores usaram este novo mapa de risco contínuo para calcular o tremor de solo esperado para um período de retorno de 475 anos — um marco padrão para o projeto de edifícios — eles encontraram resultados que são surpreendentemente diferentes do atual código nacional de construção. Na cidade de Padang, localizada diretamente acima da lacuna de Mentawai travada, o novo modelo previu acelerações de solo quase três vezes maiores do que o exigido pelo código atual. Em outras cidades como Banda Aceh, a diferença era menor, mas ainda significativa. Os pesquisadores rastrearam essa lacuna até duas escolhas específicas em sua modelagem: eles consideraram toda a gama de magnitudes possíveis sem cortar as maiores, e incluíram a variabilidade estatística total de como o tremor do solo se comporta, o que tende a produzir picos mais altos do que as estimativas médias.

O estudo também testou se este novo mapa poderia prever terremotos futuros. Eles treinaram o modelo com dados de 1998 a 2015 e depois pediram que ele previsse o padrão de terremotos de 2016 a 2024. O modelo previu com sucesso as localizações gerais onde os terremotos ocorreriam, confirmando que os fatores físicos que eles usaram — o travamento das placas e a distância da falha — são guias confiáveis de onde o solo provavelmente tremerá. No entanto, o modelo previu um número total de terremotos maior do que o que realmente aconteceu no período de teste. Os pesquisadores explicaram isso observando que o período de treinamento incluiu um tempo de atividade incomumente alta após um terremoto massivo de 2004, o que inflou temporariamente a taxa de fundo. À medida que a região retorna a um ritmo normal, espera-se que as taxas previstas diminuam, mas o padrão espacial de onde o risco é mais alto permanece estável.

As implicações destas descobertas são profundas para a forma como pensamos sobre segurança em regiões propensas a terremotos. O método tradicional de desenhar caixas em um mapa pode criar uma falsa segurança em áreas logo fora de uma zona de alto risco, ou perder o perigo em áreas quietas que estão fisicamente preparadas para uma ruptura. Ao usar uma estrutura que flui com a física da terra, os pesquisadores mostraram que o risco em Sumatra não é uniforme nem estático. O maior perigo reside na lacuna de Mentawai, um lugar que tem estado quieto por séculos, mas que é atualmente a parte mais firmemente travada da fronteira de placas. O estudo sugere que os códigos de construção precisam ser reavaliados para refletir esta realidade contín-nua baseada na física, garantindo que as estruturas sejam projetadas para suportar o verdadeiro potencial do solo sob elas. Este trabalho não oferece apenas um novo mapa; oferece uma nova maneira de ver a terra, uma que reconhece que o risco de um desastre é uma paisagem suave e mutável, em vez de um conjunto de fronteiras fixas.

Resumo Técnico: Uma Estrutura de Processo de Cox Log-Gaussiano para Avaliação Probabilística de Perigo Sísmico Livre de Zonas em Sumatra

Definição do Problema
A PSHA (Avaliação Probabilística de Perigo Sísmico) convencional, seguindo a estrutura de Cornell (1968), baseia-se em duas premissas primárias que são difíceis de justificar em contextos de subducção como Sumatra: (1) a região sismogênica é particionada em zonas de fonte de área com taxas de fundo espacialmente uniformes, e (2) eventos sísmicos ocorrem independentemente no tempo. Na realidade, o carregamento tectônico varia suavemente no espaço, em vez de ocorrer em passos discretos, e grandes sismos principais geram sequências de réplicas fortemente agrupadas tanto no espaço quanto no tempo. Os autores argumentam que essas premissas obscurecem os drivers físicos da sismicidade e necessitam de uma alternativa "livre de zonas", onde o termo de fonte espacial é uma intensidade estocástica ajustada derivada de covariáveis físicas.

Metodologia
O estudo propõe uma estrutura que substitui as zonas de fonte predefinidas por um Processo de Cox Log-Gaussiano (LGCP) para a sismicidade de fundo, acoplado a um modelo de Sequência de Após-sismos do Tipo Epidêmico (ETAS) autoexcitável para o agrupamento desencadeado.

Dados e Pré-processamento:
- Catálogo: Foi utilizado um conjunto de dados de trabalho de 5.356 eventos ( $M \ge 4,5$ ) de 1998–2024 em Sumatra ( $95^\circ\text{E}–108^\circ\text{E}$ , $6^\circ\text{S}–6^\circ\text{N}$ ).
- Desagrupamento (Declustering): Para isolar a sismicidade de fundo, aplicou-se o método de janela espaço-temporal de Gardner–Knopoff (1974), reduzindo o catálogo para 693 sismos principais. Esta etapa foi necessária porque as sequências de réplicas do megathrust de 2004–2005 aumentariam a taxa aparente de fundo em uma ordem de magnitude.
Modelagem da Intensidade de Fundo (LGCP):
- A intensidade de fundo $\lambda(s)$ é modelada como um processo duplamente estocástico onde $\log \lambda(s) = \mu(s) + Z(s)$ .
- Covariáveis: A tendência média $\mu(s)$ $μ (s)$ é uma função log-linear de duas covariáveis tectônicas fisicamente interpretáveis:
  - Acoplamento intersísmico ( $\phi$ ): A razão de acoplamento na interface da subducção de Sunda (derivada de inversão geodética), que atinge o pico no gap sísmico de Mentawai, ainda travado.
  - Distância à falha ( $D_{fault}$ ): A distância perpendicular à Grande Falha de Sumatra (GSF).
- Campo Latente: $Z(s)$ é um Campo Aleatório Gaussiano (GRF) com uma estrutura de covariância exponencial para capturar a variabilidade residual de mesoescala não explicada pelas covariáveis.
- Estimação: Os parâmetros foram estimados via estimativa de máxima verossimilhança (MLE) log-linear de Poisson em uma grade de $60 \times 60$ .
Componente Desencadeado:
- O componente de agrupamento desencadeado é descrito separadamente por um modelo ETAS autoexcitável com um kernel temporal de Omori–Utsu. Este modelo contabiliza o gatilho evento-a-evento que a taxa de fundo impulsionada por covariáveis deixa sem explicação.
Integração do Perigo:
- A superfície de intensidade LGCP ajustada serve diretamente como o termo de fonte espacial na integral de perigo de Cornell, eliminando a necessidade de delineamento de zonas.
- A taxa anual de excedência da Aceleração de Pico do Solo (PGA) é calculada utilizando uma distribuição de magnitude de Gutenberg–Richter e a Equação de Predição de Movimento de Solo (GMPE) de Youngs et al. (1997) para a interface de subducção.
Validação e Comparação:
- Seleção de Modelo: O LGCP foi comparado contra um processo de cluster de Thomas de Neyman–Scott (uma linha de base estatística sem covariáveis físicas) usando o Critério de Informação de Akaike (AIC) e a função L não homogênea.
- Validação de Previsão: Um teste de retenção temporal (treino em 1998–2015, teste em 2016–2024) e a suíte de testes do Modelo de Probabilidade de Terremotos Regional (RELM) foram usados para avaliar o desempenho da previsão espacial.

Principais Resultados

Desempenho do Modelo: O LGCP impulsionado por covariáveis superou decisivamente a linha de base de cluster de Thomas, com uma melhoria de razão de verossimilhança de $2\Delta\ell = 645,7$ ( $\chi^2_2, p < 10^{-4}$ ) e um $\Delta$ AIC de 643,7. O LGCP reproduziu com sucesso o gradiente de perigo observado, incluindo o alto risco do gap de Mentawai, sem fronteiras predefinidas.
Significância das Covariáveis: Ambas as covariáveis foram altamente significativas (Cox–Snell pseudo- $R^2 \approx 0,61$ ). O coeficiente de acoplamento foi positivo (maior acoplamento $\rightarrow$ maior taxa), e o coeficiente de distância à falha foi negativo (mais próximo da GSF $\rightarrow$ maior taxa).
Estimativas de Perigo: A estrutura gerou acelerações de projeto de 475 anos significativamente superiores ao código nacional indonésio (SNI 1726:2019). As razões de PGA derivadas do LGCP em relação ao PGA do código foram aproximadamente:
- Padang: 2,7x (1,79 g)
- Pekanbaru: 2,5x
- Medan: 2,1x
- Banda Aceh: 1,5x
Padrão Espacial: O maior perigo concentrou-se no noroeste e ao longo do corredor de Padang, correspondendo diretamente ao gap de alto acoplamento de Mentawai. Em contraste, Banda Aceh mostrou menor perigo devido ao baixo acoplamento residual após a ruptura de 2004.
Validação: A retenção temporal confirmou que o padrão espacial do LGCP é preservado em dados não vistos (quantile do teste L $\gamma = 0,21$ ). No entanto, a taxa absoluta foi superestimada por um fator de 3,3 no período de teste, atribuído ao fato de a janela de treinamento incluir uma atividade elevada de pós-2004.

Significância e Alegações
O artigo afirma demonstrar um acoplamento matemático direto entre uma intensidade de Cox ajustada e impulsionada por covariáveis e a integral de perigo de Cornell, uma combinação não demonstrada anteriormente em um contexto de subducção. A principal significância reside na capacidade de derivar movimentos de solo de projeto a partir de um modelo de fonte estocástica contínua, em vez de zonas predefinidas.

Os autores enfatizam que a estrutura:

Reproduz gradientes físicos: Captura o gap sísmico de Mentawai e o gradiente de perigo através de Sumatra usando apenas covariáveis físicas (acoplamento e distância à falha), sem inserir histórico de ruptura ou limites de zona.
Oferece uma orientação conservadora: Os valores de perigo mais elevados em comparação ao código nacional são atribuídos a duas escolhas de modelagem específicas e defensáveis: (a) o uso de um piso de completude de magnitude ( $M_c=4,5$ ) consistente com o catálogo, que pondera eventos maiores de forma mais intensa, e (b) a aplicação de um termo de variabilidade aleatória não truncado no GMPE.
Distingue escalas de tempo: Os autores esclarecem que, embora o LGCP seja apropriado para o design de engenharia de longo prazo (décadas/séculos), onde o carregamento tectônico governa a taxa, o modelo ETAS autoexcitável é melhor adequado para previsões operacionais de curto prazo (dias/meses), onde o agrupamento desencadeado domina.

O artigo conclui que a estrutura LGCP livre de zonas fornece uma alternativa geologicamente sólida e reproduzível à tradicional PSHA de fonte de área, particularmente em regiões onde o carregamento tectônico varia suavemente e os dados de ruptura histórica são escassos (ex: o gap travado de Mentawai).

Resumo Técnico: Uma Estrutura de Processo de Cox Log-Gaussiano para Avaliação Probabilística de Perigo Sísmico Livre de Zonas em Sumatra

Mais como este