O modelo, em aberto

Como funciona

Para cada incêndio activo, o browser corre uma simulação Rothermel 1972 ao longo de 12 direcções a partir do ponto de ignição. Cada direcção avança em passos de 5 minutos, lendo o terreno, o combustível e o tempo a cada passo.

• 01Lê declive e exposição (aspecto) do terreno Cesium ~30 m no ponto de ignição.
• 02Procura a classe de combustível (mato/pinhal/eucalipto/floresta mista/agrícola/urbano/água/rocha) num mapa nacional construído a partir do OpenStreetMap.
• 03Interpola vento, temperatura e humidade pelas 5 estações IPMA mais próximas (IDW vectorial — preserva direcção do vento, não só intensidade).
• 04Aplica RCM do concelho (IPMA 1-5) como amplificador global do ROS.
• 05Quando disponível, aplica FWI DMC/DC (Van Wagner 1987) para captar a secura cumulativa que a humidade snapshot não vê.
• 06Re-amostra terreno e combustível a cada 4 passos (≈ 20 min de fogo) — captura mudanças de paisagem ao longo do caminho.
• 07Os pontos a 30 min / 1 h / 3 h / 6 h em todas as 12 direcções formam os perímetros previstos.

Componentes activos

Precisão histórica

Para cada incêndio do corpus, comparámos o perímetro previsto a 24 h com o perímetro real (ICNF Áreas Ardidas). O índice Jaccard mede a sobreposição — 0 = sem sobreposição, 1 = perímetro perfeito.

Buckets agrupados por região × tipo de combustível dominante × intervalo de RCM. Mínimo 3 amostras por bucket.

Comparação com outros sistemas

Honestamente: em precisão pura, sistemas calibrados comercialmente (Tecnosylva) ou operados por analistas (FARSITE da BMAP da ANEPC) batem-nos. O que nos diferencia é acessibilidade — somos o único sistema público em PT que mostra previsão de propagação ao cidadão em tempo real.

Limitações conhecidas

O que ainda não fazemos bem, e por porquê:

• ⚠Sem spotting (transporte de embers a longa distância). Em incêndios extremos, pavesas saltam centenas de metros à frente da frente principal e iniciam novos focos. Modelo de Albini (1979) está na roadmap.
• ⚠Sem feedback loop em tempo real. Quando os bombeiros reportam perímetro real, não corrigimos o modelo automaticamente. Sistemas profissionais fazem isto.
• ⚠Mapa de combustíveis OSM, não ICNF. A Cartografia da Vegetação Florestal do ICNF é mais precisa, mas tem barreiras de licenciamento. Em PT, eucalipto e pinhal estão particularmente subestimados no OSM.
• ⚠Sem vento variável no tempo. Usamos uma snapshot da estação no momento da previsão. Em incêndios de 6 h o vento pode mudar 90°.
• ⚠Calibração só contra incêndios ≥ 1000 ha. Para incêndios pequenos (< 100 ha) o modelo não foi validado contra ground truth.
• ⚠Sem validação cega contra BMAP/ANEPC. Os analistas profissionais correm modelos no terreno; comparar com os deles seria o teste definitivo.

Próximos passos

Trabalho explícito para subir o Jaccard. Em ordem de impacto/esforço:

1. 01Re-calibração com fuel map + FWI activo (Jaccard esperado 0,55-0,60).
2. 02Modelo de spotting Albini 1979 (+5-10 % em incêndios extremos).
3. 03Vento variável no tempo via IPMA forecast (não apenas observação).
4. 04Mapa de combustíveis ICNF (Cartografia da Vegetação Florestal).
5. 05Validação cega independente: pedir à ANEPC ou BMAP para correr o modelo num incêndio real e comparar.
6. 06Feedback loop: incorporar perímetros Copernicus EMS quando activados.