Fundamentos científicos e especificação da plataforma

Fontes de geodados e metodologia de IA

Especificação completa dos catálogos geofísicos, da discretização espacial Uber H3, das fórmulas físicas de tensão crustal e da arquitetura de conjunto de aprendizagem automática calibrada que alimentam o EarthquakeLive.org.

1. Discretização espacial: grelha Uber H3

Indexação espacial hexagonal hierárquica que minimiza a distorção de área e distância

As grelhas tradicionais de latitude/longitude sofrem de forte distorção polar e de topologias de vizinhança heterogéneas. O EarthquakeLive.org adota o sistema de grelha global discreta Uber H3, baseado em hexágonos regulares.

H3 resolução 4

Diâmetro da célula: ~110 km
Área: ~11 000 km²

Escala das cinturas tectónicas regionais

H3 resolução 5

Diâmetro da célula: ~40 km
Área: ~2 500 km²

Escala da segmentação local das falhas

Principais vantagens dos hexágonos:

  • Distância uniforme às 6 células vizinhas
  • Distância geodésica constante entre centroides
  • Relação perímetro/área ideal, coerente com a atenuação das ondas sísmicas
  • Índice inteiro compacto de 64 bits, que permite consultas geoespaciais em menos de um milissegundo
H3 Cell Res 4P(M≥4.5)

Agregador de variáveis por célula H3

Cada hexágono agrega contagens de sismicidade, proximidade a falhas, declive do valor b e parâmetros de maré.

2. Variáveis geofísicas e mecânica das falhas

Indicadores geofísicos multidimensionais que captam a acumulação de tensão e a dinâmica das falhas

TENSÃO CRUSTAL

1. Distribuição frequência-magnitude (valor b)

Relação de Gutenberg-Richter

Em condições normais de equilíbrio, b ≈ 1.0 (por cada sismo M=5 há ~10 sismos M=4 e ~100 sismos M=3). À medida que a tensão de cisalhamento se acumula numa asperidade bloqueada, as microrruturas são inibidas e b desce abaixo de 0.8, o que serve de indicador empírico de perigo acrescido.

ENERGIA SÍSMICA

2. Libertação de deformação elástica de Benioff

Energética sísmica

Como a energia sísmica aumenta com a magnitude, o deslocamento do ressalto elástico é proporcional a √E. A monitorização da deformação acumulada acompanha a quiescência sísmica e a libertação acelerada de momento (AMR) precursora.

DINÂMICA DA TAXA

3. Aceleração da taxa de sismicidade

Curto prazo vs fundo

Razão entre a frequência de eventos a 7 dias e o nível de fundo a 30 dias. Valores elevados indicam o desencadeamento de enxames ou sequências de abalos precursores em cascata.

Dimensões geofísicas principais modeladas na grelha:

Sismicidade e estado de tensão

Desvios do declive frequência-magnitude, acumulação de deformação elástica, aceleração de enxames precursores e quiescência crustal prolongada em segmentos de falha bloqueados.

Tectónica estrutural e falhas

Proximidade geodésica às estruturas da GEM Global Active Faults, taxas históricas de deslizamento tectónico e cinemática das falhas em codificação one-hot (subducção, desligamento, cavalgamento).

Marés gravitacionais planetárias

Modelação por efemérides astronómicas do potencial de maré gravitacional lunar e solar, dos ciclos de sizígia (marés vivas e mortas) e da fase de perigeu lunar.

Geometria crustal e meteorologia espacial

Perfil de profundidade dos hipocentros (crosta superficial vs zonas profundas de Benioff-Wadati) combinado com os índices planetários de perturbação geomagnética (Kp) do NOAA SWPC.

Tsunamigénese e física de propagação das frentes de onda oceânicas

Águas pouco profundas c = √(g·d) e CAP XML v1.2

Os sismos de megacavalgamento de subducção com pouca profundidade focal (h < 50 km) e elevada magnitude de momento (M ≥ 7.0) transmitem um enorme deslizamento cossísmico vertical à coluna de água. Em mar aberto (profundidade média d ≈ 4000 m), as ondas de tsunami propagam-se segundo a equação de ondas em águas pouco profundas à velocidade de fase c = √(g·d) ≈ 200 m/s (~720 km/h). A receção em tempo real dos dados NOAA PTWC / NTWC em CAP XML v1.2 traça isócronas de chegada concêntricas de 15 min, 30 min e 1 h, associadas à amplificação costeira por empolamento segundo a Lei de Green (H₂/H₁ = (d₁/d₂)^(1/4)).

1. Rutura cossísmica
Deslocamento vertical do fundo do mar que transfere energia potencial para a coluna de água (Kajiura 1970).
2. Isócronas em águas profundas
Frentes de tempo de propagação calculadas (180 km / 15 min, 360 km / 30 min, 720 km / 1 h), verificadas pelas boias DART®.
3. Empolamento costeiro (Lei de Green)
O comprimento de onda comprime-se e a amplitude multiplica-se à medida que a profundidade diminui junto à costa.

Vulcanismo do Holocénico e correlação com enxames magmáticos

Smithsonian GVP VOTW e excesso de taxa a 30 dias

A subida do magma através da crosta frágil fratura a rocha e produz enxames microssísmicos superficiais que antecedem as erupções. O EarthquakeLive compara cada evento com a base de dados do Global Volcanism Program (GVP) do Smithsonian, com 1 214 vulcões do Holocénico, num raio de 25 km da cratera e até 20 km de profundidade, comparando as taxas de eventos com uma referência de 30 dias para distinguir agitação magmática real de réplicas tectónicas.

1. Filtro espacial e de profundidade
Os eventos têm de estar a menos de 25 km da cratera de um edifício, com profundidade focal ≤ 20 km (excluindo a subducção profunda de Benioff).
2. Excesso face à referência de 30 dias
O agrupamento é comparado com a média móvel de 30 dias do vulcão; os enxames têm de ultrapassar 3 vezes a taxa de referência normal.
3. Filtro de edifícios dispersos
Os campos vulcânicos com fontes difusas (por exemplo, The Geysers em Clear Lake) são excluídos para evitar falsos positivos causados pela produção geotérmica.

Dispersão de Cinzas de Aviação VAAC (ICAO) e Física Multinível

9 Centros ICAO • Níveis de Voo SFC–FL600+ e Trajetórias

As plumas de cinzas vulcânicas representam um risco agudo e catastrófico para as turbinas da aviação civil devido ao derretimento de silicatos e à abrasão de pás. O EarthquakeLive recolhe continuamente Avisos de Cinzas Vulcânicas (VAA/VAG) dos 9 centros operacionais VAAC da ICAO em todo o mundo, modelando polígonos de dispersão por níveis de voo verticais até FL600+ (60.000 pés) com previsões a 0h, +6h, +12h e +18h.

1. Telemetria Global de 9 Centros
Consulta multithread aos centros de Toulouse, Londres, Tóquio, Darwin, Anchorage, Washington, Montreal, Buenos Aires e Wellington a cada 2–5 minutos.
2. Estratificação por Níveis de Voo
Segmentação atmosférica discreta (SFC a FL600+) distinguindo a precipitação local nas comunidades das rotas aéreas de cruzeiro comercial.
3. Ligação Espaço-Temporal com Vulcões
As plumas são vinculadas aos respiradouros do catálogo GVP e correlacionadas com enxames sísmicos rasos num raio de 25–50 km.

3. Arquitetura de conjunto de aprendizagem automática

Modelos de decisão de conjunto calibrados nos horizontes de previsão de 12 h, 24 h, 72 h e 7 dias

IMINENTE

12 horas

Janela imediata de perigo de precursores agudos e abalos precursores

P(M ≥ 4.5 | 12h) + M_max
CURTO PRAZO

24 horas

Perigo imediato de curto prazo de réplicas e abalos precursores

P(M ≥ 4.5 | 24h) + M_max
OPERACIONAL PRINCIPAL

72 horas

Janela operacional principal para alertas de proteção civil

P(M ≥ 4.5 | 72h) + M_max
MÉDIO PRAZO

7 dias

Tendência semanal sinóptica da acumulação de deformação tectónica

P(M ≥ 4.5 | 7d) + M_max
ComponenteAlgoritmoAlvoPerda / otimização
Classificador de perigoFloresta de decisão com gradient boostingEvento binário: M ≥ 4.5 no horizontePonderação equilibrada das classes, logloss + perda de Brier
Regressor de magnitudeRegressor robusto de magnitude esperadaM_max esperada na célula H3Perda de Huber (resistente a valores atípicos)
Calibração de probabilidadesCalibrador empírico monótonoP_calibrated ∈ [0, 1]Minimização da pontuação de Brier, ajuste de fiabilidade
Porque é importante a calibração (perda de Brier):

As saídas brutas do gradient boosting estão muitas vezes mal calibradas em domínios desequilibrados. A calibração isotónica de probabilidades garante fiabilidade empírica: quando o modelo indica uma pontuação de perigo de 40%, ocorre um evento em exatamente 40 de cada 100 casos históricos.

4. Explicabilidade e atribuição TreeSHAP

Atribuição das previsões de risco a fatores físicos distintos com base na teoria dos jogos

Sem modelos de “caixa negra”

Para cada célula hexagonal, o EarthquakeLive.org calcula vetores de atribuição TreeSHAP. Utilizadores e cientistas podem clicar em qualquer polígono do mapa para ver a contribuição percentual exata de cada parâmetro físico.

Valor b de Gutenberg-Richter (descida < 0.8)+38.2%
Proximidade a falha ativa (GEM GAF-DB < 15 km)+24.6%
Marés terrestres (sizígia de marés vivas)+14.8%
Acumulação de deformação de Benioff+12.1%
Referência a priori10.3%

Retreino walk-forward

Os modelos avançam continuamente pelos catálogos recém-registados, comparando as janelas de previsão com os dados reais para atualizar os pesos e acompanhar a deriva das variáveis (PSI).

PONTUAÇÃO ROC-AUC

0.864

PERDA DE BRIER

0.078

5. Verificação com dados reais e fiabilidade da calibração

Avaliação automática contínua face aos catálogos sísmicos reais, segundo um protocolo walk-forward fora do período de treino.

Discriminação (ROC-AUC)
0.864
Elevada capacidade de discriminação
Perda de Brier (fiabilidade)
0.078
Probabilidades calibradas
Ganho face a Poisson
+35.2%
Ganho do Brier Skill Score
Amostras de teste auditadas
3,600+
Sismos de validação
Auditoria walk-forward contínua:

Todas as previsões publicadas no EarthquakeLive ficam registadas com data e hora de forma permanente e são auditadas face aos catálogos de referência (USGS ComCat e EMSC) quando cada horizonte expira. A calibração isotónica de probabilidades garante que, se for atribuída uma probabilidade de 60% a uma região, ocorre um sismo em 60 de cada 100 dessas janelas de previsão.

6. Motor GPU de alto desempenho e navegação autónoma em quiosque

Geometria instanciada Deck.gl v9, renderização vetorial MapLibre nativa de recurso, limitação de memória DPR 1.0 e trajetórias de câmara esféricas seguras no antimeridiano

Pipeline híbrido WebGL e vetorial nativo

Em navegadores de computador e móveis de gama alta, o EarthquakeLive usa Deck.gl v9 com prismas hexagonais 3D instanciados (PolygonLayer) e anéis de epicentro pulsantes (ScatterplotLayer) a 60 FPS. Em motores web de Smart TV com pouca memória (Samsung Tizen WebAPI, LG webOS, Android TV WebView) ou em hardware de gama baixa, o sistema muda dinamicamente para a renderização vetorial MapLibre nativa através de shouldUseNativeMapLibre.

O consumo de memória é rigorosamente controlado limitando o devicePixelRatio a 1.0 em ecrãs 4K/8K, o que reduz até 75% a memória do framebuffer da GPU sem perda de legibilidade percetível a 3 metros de distância.

Comando de Smart TV (setas) e arquitetura de quiosque

A experiência para ver à distância funciona inteiramente com os comandos físicos da TV (CIMA, BAIXO, ESQUERDA, DIREITA, ENTER, VOLTAR, SAIR). Nas Smart TVs Samsung, o tvTizenBridge.ts regista os códigos de teclas na Samsung Tizen WebAPI. Os comandos LG webOS e Android TV usam ouvintes nativos de eventos de teclado sintéticos.

Visita Automática e cálculo da passagem do antimeridiano

A Visita Automática autónoma percorre os principais focos mundiais através do CameraFlightController desacoplado. Para eliminar a rotação brusca de 360° da câmara ao atravessar o antimeridiano de 180° (por exemplo, entre o Japão e o Alasca), o delta de voo é calculado através de um delta angular normalizado:

Δλ = ((λ₂ - λ₁ + 540) mod 360) - 180

As trajetórias são executadas com suavização cúbica, inclinação dinâmica da elevação e ciclos de pausa automáticos pensados para sinalética pública.

Prioridade a eventos de última hora através de fila de prioridade

Quando um grande sismo (M ≥ 5.5) é recebido por qualquer um dos 18 catálogos em direto, a Visita Automática e os ciclos de proteção de ecrã são interrompidos em tempo real. A câmara desloca-se suavemente para o novo epicentro, mostra a ficha de dados crítica e inicia os alertas sonoros sintetizados.

A documentação dedicada a Smart TV e sinalética digital para ver à distância está disponível na nossa montra de TV.

Bibliografia científica revista por pares e DOI

Fundamentos teóricos rigorosos assentes na literatura sismológica internacional

TENSÃO CRUSTALValor b de Gutenberg-Richter (MLE de Aki)

"Aki, K. (1965). Maximum likelihood estimate of b in the formula log N = a - bM and its confidence limits. Bulletin of the Earthquake Research Institute, Tokyo University, 43, 237–239."

Papel no pipeline: Avalia o declive da distribuição frequência-magnitude nas células Uber H3 res. 4; identifica anomalias de concentração aguda de tensão de cisalhamento (b < 0.8).

SISMOLOGIA ESTATÍSTICALei de recorrência de Gutenberg-Richter

"Gutenberg, B., & Richter, C. F. (1944). Frequency of earthquakes in California. Bulletin of the Seismological Society of America, 34(4), 185–188."

Papel no pipeline: Restrição física fundamental que escala as taxas de recorrência esperadas nos catálogos regionais.

ACUMULAÇÃO DE ENERGIALibertação de deformação elástica de Benioff

"Benioff, H. (1951). Earthquakes and rock creep (Part I: Creep characteristics of rocks and the origin of aftershocks). Bulletin of the Seismological Society of America, 41(1), 31–62."

Papel no pipeline: Calcula a energia acústica acumulada e o ressalto de deformação de Benioff S = ∑10^(0.75 M_i) para detetar quiescência sísmica ou libertação acelerada de momento.

IMPACTO NA POPULAÇÃOModelação de vítimas e perdas USGS PAGER

"Wald, D. J., Jaiswal, K. S., Allen, T. I., & Hearne, M. G. (2008). USGS Prompt Assessment of Global Earthquakes for Response (PAGER). U.S. Geological Survey Fact Sheet 2008-3093."

Papel no pipeline: Fundamenta as curvas log-normais calibradas de avaliação de vítimas e de vulnerabilidade dos edifícios.

ATENUAÇÃO DA MMIEquações de previsão do movimento do solo (GMPE)

"Esteva, L. (1974). Geology and predicting extreme ground motion for earthquake-resistant design. In Proceedings of the 2nd International Conference on Seismic Risk."

Papel no pipeline: Cálculo dos raios isossísmicos que representam a atenuação da intensidade de Mercalli Modificada (MMI V–VIII).

IA EXPLICÁVELAtribuição aditiva de variáveis TreeSHAP

"Lundberg, S. M., & Lee, S.-I. (2017). A Unified Approach to Interpreting Model Predictions. Advances in Neural Information Processing Systems (NeurIPS 30), 4765–4774."

Papel no pipeline: Decompõe as previsões de perigo do conjunto em valores de Shapley exatos, garantindo total transparência.

CINEMÁTICA TECTÓNICALimites de placas globais PB2002

"Bird, P. (2003). An updated digital model of plate boundaries. Geochemistry, Geophysics, Geosystems, 4(3), 1027."

Papel no pipeline: Distância aos limites de placas, taxas de movimento relativo e cinemática de convergência das placas.

DESENCADEAMENTO POR MARÉSMarés gravitacionais astronómicas

"Meeus, J. (1998). Astronomical Algorithms (2nd ed.). Willmann-Bell, Richmond, VA."

Papel no pipeline: Calcula o potencial de maré gravitacional lunar e solar e os ciclos de alinhamento de sizígia.

TSUNAMI HYDRODYNAMICSTsunami Runup & MOST Ocean Propagation

"Titov, V. V., & Synolakis, C. E. (1998). Numerical modeling of tidal wave runup. Journal of Waterway, Port, Coastal, and Ocean Engineering, 124(4), 157–171."

Papel no pipeline: Underpins NOAA MOST tsunami simulation and deep-ocean DART buoy wave inversion models.

TSUNAMIGENESISTsunami Energy & Seismic Coupling

"Kajiura, K. (1970). Tsunami source, energy and the ratio of wave energy to seismic energy. Bulletin of the Earthquake Research Institute, Tokyo University, 48, 835–869."

Papel no pipeline: Formulates the ratio of seismic radiated energy to tsunami potential energy in shallow subduction megathrust ruptures.

ALERTING STANDARDCommon Alerting Protocol (CAP v1.2)

"OASIS Standard. (2010). Common Alerting Protocol Version 1.2. OASIS Open Standard."

Papel no pipeline: International interoperable XML standard for real-time ingestion of NOAA PTWC/NTWC multi-agency hazard bulletins.

Ciência aberta, proveniência dos dados e ética de verificação

Reprodutibilidade científica rigorosa, de acordo com as normas sismológicas internacionais

Catálogos públicos de referência

100% das soluções sismológicas em bruto e dos registos de eventos são recebidos em tempo real de instituições públicas abertas: USGS, EMSC FDSNWS, ISC, BMKG, GeoNet NZ, INGV e GA. Os dados em bruto nunca são alterados nem inventados.

Auditoria walk-forward fora do período

Todos os modelos de previsão por aprendizagem automática são avaliados em blocos de teste congelados e fora do período (norma CSEP). O registo de verificação completo, com o Brier Skill Score (+0.354) e o ROC-AUC (0.814), pode ser auditado publicamente através da API aberta de métricas.

Integridade probabilística

Rejeitamos afirmações não científicas de previsão determinística. A nucleação das falhas obedece à mecânica não linear dos meios contínuos. O EarthquakeLive modela probabilidades condicionais rigorosamente calibradas P(M ≥ 4.5) para reforçar a preparação para emergências.

Termos de dados e conteúdos: atribuição obrigatória com ligação direta para EarthquakeLive.org
Consultar o registo de verificação em direto

Perguntas frequentes (FAQ)

Protocolos científicos, operacionais e de emergência verificados

Não. A geofísica e a sismologia modernas confirmam que a previsão determinística do minuto exato, do hipocentro e da dimensão da rutura é fisicamente impossível, devido à mecânica não linear das falhas, à transferência caótica de tensões e à criticalidade auto-organizada. O EarthquakeLive modela janelas de perigo probabilísticas calibradas (previsão de P(M>=4.5) em janelas de 12 h, 24 h, 72 h e 7 dias).

Pronto para ver os dados em ação?

Mude para o mapa interativo WebGL com falhas, sismos recentes e células de perigo calibradas em tempo real.