PORTULAN/parlamento-pt huggingface dataset:
Figura 3.13. Linhas isócronas de propagação do complexo de incêndios de Sandomil-Seia (1a ocorrência, sector centro-leste), Esculca-Côja-Arganil (2a ocorrência, com duas ignições junto ao limite sul) e Vide-Seia (3a ocorrência, sudeste). O fogo de Sandomil eclodiu no centro leste da mancha queimada, numa paisagem com ocupação florestal mais fragmentada comparativamente às ignições da serra do Açor. Como tal, o seu desenvolvimento foi mais lento, e também mais difícil de caracterizar. A sul, a ocorrência designada como Esculca-Côja-Arganil procede de dois reacendimentos distintos separados por 4 km, respetivamente na imediação da mata da Margaraça, junto a Relva Velha, e mais a oeste próximo de Monte Redondo. Estas duas ignições tiveram desenvolvimentos distintos. A primeira progrediu principalmente para NNE e deslocou-se mais rapidamente, atingindo cerca de 4,5 km/h das 15 às 17 horas. A segunda avançou em direção a Côja (NNW) à taxa de 1-2 km/h, tendo, todavia, aumentado para quase 4 km/h no período 18-19 horas quando atingiu a periferia leste de Oliveira do Hospital. Os impactos mais notórios em Oliveira do Hospital, nomeadamente na zona industrial, e nas aldeias a sul, resultaram, portanto, da ignição de Relva Velha. 17 DE ABRIL DE 2018_____________________________________________________________________________________________________________________67
Figura 3.14. Linhas isócronas de propagação do complexo de incêndios de Ponte das Portelinhas-Figueiredo-Sertã (1a ignição) e Nespereira-Pedrogão Pequeno-Sertã (canto sudoeste). Complexo de Pataias Os incêndios do litoral são consideravelmente mais simples na sua evolução comparativamente às ocorrências do interior. Para tal concorreram a vegetação relativamente uniforme, terreno plano e condução pelo vento. Dois reacendimentos na freguesia de Pataias em Praia da Légua (14:01 horas) e Burinhosa (14:33 horas) deram origem ao incêndio que devastou a quase totalidade das matas nacionais do Pedrogão e Leiria e aproximadamente metade da mata nacional do Urso. O incêndio foi predominantemente de superfície (mas com crestamento total das copas) nos pinhais de estatura superior a 10-12 m, e de copas nos restantes. 17 DE ABRIL DE 2018_____________________________________________________________________________________________________________________69
A direção de propagação do incêndio manteve-se paralela à linha de costa, portanto NNE, com alguma rotação para leste após Vieira de Leiria, que aliás funcionou como um obstáculo parcial, estreitando o incêndio e provavelmente reduzindo a sua velocidade de deslocação. Nas duas primeiras horas os fogos propagaram-se a cerca de 1 km/h, valor este que quadruplicou à 3.a hora e se manteve durante a 4.a hora quando as frentes se juntaram. No entanto, a velocidade de propagação viria ainda a aumentar para 6 km/h entre as 18 e as 21 horas, atingindo um máximo de 6,9 km/h na mancha da mata de Leiria que havia ardido em 2003; valores tão elevados só se viriam a registar entre as 22 e as 23 horas, 5,8 km/h, de novo numa área essencialmente desarborizada. Globalmente, o maior crescimento do incêndio deu-se entre as 19 e as 21 horas, com uma média de 2249 ha/h. O crescimento do incêndio abrandou substancialmente a partir da 2.a hora do dia 16 de outubro diminuindo progressivamente durante o resto da noite. Figura 3.15. Linhas isócronas de propagação do complexo de incêndios de Légua-Pataias (1.a ignição, extremo sul) e Burinhosa-Pataias (2.a ignição). II SÉRIE-E — NÚMERO 16_____________________________________________________________________________________________________________________70
Quiaios-Figueira da Foz O incêndio de Quiaios teve causa intencional, com início declarado às 14:36 horas. Os seus padrões de crescimento e orientação são similares aos do incêndio de Pataias. Assinale-se que durante quatro horas manteve uma velocidade de propagação consistentemente elevada (média 4,8 km/h das 16 às 20 horas, com um máximo de 5,4 km/h), apesar de não alcançar os picos registados na mata nacional de Leiria. A área de fogo de copas foi, no entanto mais expressiva em Quiaios, dada a representatividade de pinhal baixo (~8 m) com muito mato e acácias sob coberto, um legado do grande incêndio de 1993; foi possível observar no terreno a mudança para fogo de superfície aquando da entrada do fogo em pinhal adulto com descontinuidade vertical. Figura 3.16. Linhas isócronas de propagação do incêndio de Quiaios-Figueira da Foz. 3.5.3. Influência dos incêndios do passado O legado do regime de fogo e dos fogos antecedentes, através da sua influência na paisagem – padrões espaciais de acumulação de combustível e de tipo de vegetação – condiciona as 17 DE ABRIL DE 2018_____________________________________________________________________________________________________________________71
(piro-cúmulo-nimbo). Para que este fenómeno possa suceder há normalmente dois pré-requisitos: ● Uma base da atmosfera alta e extremadamente seca (V invertido entre a adiabática seca e a húmida) que favorece os fenómenos de microburst, alimentados por precipitações escassas que criam no processo de evaporação fenómenos de outflowou downdraft; ● Presença de humidade na atmosfera, contribuída pela superfície ou por frentes frias próximas. 3.5.4.1. Piro-convectividade no dia 15 de outubro de 2017 A estrutura atmosférica a 15 de outubro claramente indica evolução para as condições necessárias, primeiro na estrutura de V invertido e depois com o aumento de humidade ao fim do dia, modificando-se de 20-25% para 75-85%. A estrutura em V invertido identifica-se também pelos valores elevados (10-11) do índice C-HAINES. A radiosondagem de Lisboa descreve como a estrutura vertical da atmosfera se alterou entre 14 e 16 de outubro. Mostra também como no dia 15 a atmosfera evolui para condições mais facilitadoras da transição de PyroCu para PyroCb, devido ao perfil de V invertido, mas com subida substancial da parte mais baixa em altitude. Assim, nas horas finais do dia, o aumento da humidade relativa do ar conjugado com a energia do incêndio permite que o LCL se situe até 10-11 km, muito próximo da zona de plena convecção (Figura 3.17). Figura 3.17. Radiosondagem sobre Lisboa, 14 de outubro de 2017. Como é visível nas radiosondagens, a estrutura do dia 14 exige muita energia para ascender do LCL para o Nível de Condensação por Convecção - CCL (Convective Condensation Level), o que o torna improvável. A do dia 15 (Figura 3.18), e ainda mais a do dia 16 (Figura 3.19), já se acerca de uma condição na qual o LCL tende a subir devido à energia dos incêndios, aproximando-se perigosamente do CCL e consequentemente da transição de PyroCu para PyroCb a partir da qual seriam possíveis fenómenos piro-convectivos de downdraft com alteração do crescimento dos incêndios. Destaque-se o LCL subiu de 1,2 km para quase 3 km num dia e na zona de Viseu para 3,4 km (IPMA 2018). Nesta situação e na presença de incêndio, 17 DE ABRIL DE 2018_____________________________________________________________________________________________________________________73
que sobe o LCL quase 1 km, a transição a CCL é mais que provável (CLEMENTS e LAREAU 2016). É claro que as condições descritas serão mais favoráveis à medida que o dia progride e a chegada de mais humidade facilite ainda mais a capacidade da coluna de aceder à zona de convecção livre. Estas condições são confirmadas pela análise radar das plumas dos incêndios efetuada pelo IPMA, sendo identificados oito casos de piro-convecção durante o episódio do dia 15 ao dia 16 de outubro (Figura 3.19). Como é visível, a coluna da maioria acerca-se dos 10 km identificados como sendo a altura de convecção livre estimada por radiosondagem. Figura 3.18. Radiosondagem sobre Lisboa, 15 de outubro de 2017. II SÉRIE-E — NÚMERO 16_____________________________________________________________________________________________________________________74
Figura 3.19. Quadro de PyroCb e mapa de PyroCu e PyroCb detetados pelo IPMA com radar. Reproduzido de IPMA (2018). 3.5.4.2. Momentos de convectividade Os dados de IPMA identificam dois momentos básicos de piro-convecção, a tarde de 15 e a noite de 15 para 16 de outubro. O dia 15 de outubro está claramente sob a influência do vento de sul impelido pelo Ophelia o que faz com que estes incêndios sejam conduzidos pelo vento. A piro-convecção dá-se em incêndios sem vento significativo e com a atmosfera instável (ROTHERMEL 1991), mas também com vento, quando o incêndio consegue romper a estratificação ou dispersão que aquele impõe e cria uma coluna convectiva, alterando a meteorologia local. Operacionalmente, esta capacidade é observável através dos perfis de vento de BYRAM (1954) e dos jatos troposféricos de baixa altitude (low jets). A análise dos perfis atmosféricos, tanto daqueles previstos pelo modelo GFS (Global Forecast System) como do que foi extraído da radiosondagem efetuada sobre Lisboa, confirmam a presença de um low jet que vai enfraquecendo ao longo do tempo. Portanto, inicialmente e nas horas posteriores ao meio dia, os incêndios teriam essa capacidade piro-convectiva no âmbito do fluxo geral do vento, a qual subsequentemente desapareceria com a eliminação do low jet. Tal condiz com a primeira onda de piro-convecção detetada pelo IPMA e que coincide com o aumento do low jet na previsão das 15:00 e o seu posterior desvanecimento em superfície. O II SÉRIE-E — NÚMERO 16_____________________________________________________________________________________________________________________76
low jet da Figura 3.20 condiz com o perfil de piro-convecção de Byram tipo 1a às 12 e às 15 horas, mas já não às 18 horas. Este tipo de piro-convectividade gera incêndios estreitos e alongados, como aqueles observados na primeira fase de crescimento da área ardida no dia 15 de outubro, nos quais a pluma de convecção gera indrafts potentes que estendem mais rapidamente a cabeça do incêndio. Observando as radiosondagens de Lisboa (Figura 3.21) é visível a diminuição da velocidade do vento acima do CCL na zona de convecção livre. Figura 3.20. Perfil de vento para cinco posições entre Coimbra e Vigo para três horas do dia 15 de outubro de 2017 segundo o modelo GFS. Figura 3.21. Perfil real extraído da radiosondagem de Lisboa do dia 15 de outubro de 2017, na qual se observa um low jet condizente com o perfil de piro-convecção de Byram tipo 1a. A descrição anterior é bem visível nas imagens do satélite Meteosat (Figura 3.22). A leitura da altitude das nuvens permite diferenciar as 17 UTC, com topos de nuvens até 9 km de altitude nas duas regiões de incêndios, das 18:30 horas, que reflete a dispersão ocorrida com a perda 17 DE ABRIL DE 2018_____________________________________________________________________________________________________________________77
de intensidade do low jet e à altura de 5 km, já no LCL normal e esperado onde o piro-cúmulo se dispersa. Figura 3.22. Imagem Meteosat do topo de nuvens em 15 de outubro de 2017, presente às 17 horas UTC até 9 km de altitude (cima) e ausente às 18:30 horas UTC (baixo). Os dados meteorológicos da enorme maioria das estações do IPMA correspondem a máxima intensidade de frente de chamas entre as 12 e as 15 horas UTC (Figura 3.23), período durante o qual é detetado o maior número de PyroCb no quadro de um incêndio dominado pelo vento. A segunda onda de piro-convecção no Centro de Portugal ocorre entre as 10 horas do dia 15 de outubro e as 3 horas do dia 16 de outubro e coincide com o máximo crescimento do conjunto de II SÉRIE-E — NÚMERO 16_____________________________________________________________________________________________________________________78
incêndios em área e em velocidade. Mas este 2.o momento piro-convectivo não é explicável pela meteorologia. O fenómeno é agora menos influenciado pelo vento e mais pela grande superfície que arde com frentes de chamas intensas e extensas, ainda que tenha passado a fase de máxima intensidade, tal como determinada pelas condições meteorológicas. Este tipo de piro-convecção é mais perigoso, já que se associa a downdrafts da pluma de convecção e ao crescimento do incêndio bastante mais em largura e em várias direções, em contraste com o fenómeno anterior (ROTHERMEL 1991). Figura 3.23. Intensidade potencial da frente de chamas (kW/m) com base na meteorologia das diferentes estações do IPMA e para as condições dominantes de combustível. A linha vermelha a tracejado assinala o limite da capacidade de extinção. É grande a disparidade entre a intensidade do incêndio modelada a partir dos dados meteorológicos observados (Figura 3.23) e a intensidade reconstruída com base na análise de crescimento e velocidade de propagação e efeitos observados no terreno. A intensidade máxima da frente de chamas tal como determinada pela meteorologia dá-se a meio do dia (1.o momento de piro-convecção), porém o crescimento máximo ocorre de noite (2.o momento de piro-convecção). Esta diferença entre os picos meteorológicos e real de crescimento do fogo está de acordo com o que foi observado em Pedrogão Grande em junho de 2017. Para os dois incêndios de Burinhosa-Pataias (Figura 3.24) e Figueiredo-Sertã (Figura 3.25) vê-se claramente que a intensidade máxima se verificou após o pico teórico ditado pela meteorologia. Cabe destacar que todos os momentos máximos de crescimento do incêndio e de intensidade máxima real da frente de chamas ocorreram quando a humidade começou a aumentar. Esta é uma das condições básicas para a formação e/ou transição para PyroCb e está amplamente documentada na literatura, e.g. ROTHERMEL (1991). Há, no entanto, que diferenciar dois aspetos críticos: ● Nos incêndios litorais nas matas públicas a piro-convecção com vento dirige o vento e dissipa-se quando o low jet desaparece, já que a estrutura atmosférica sem low jetdissipa a piro-convecção, o que se observa bem na imagem de satélite; ● No incêndio de Figueiredo-Sertã, e ainda que haja piro-convecção por efeito do vento inicial, a expansão mais importante dá-se já de noite e longe dos picos de vento e intensidade máxima modelada, sendo uma expressão da piro-convecção com downdraft e transição da pluma de PyroCu para PyroCb. 17 DE ABRIL DE 2018_____________________________________________________________________________________________________________________79
Figura 3.24. Intensidade real e modelada do incêndio de Burinhosa-Pataias das 15 horas do dia 15 de outubro às 3 horas do dia 16 de outubro. Repare-se na evolução da humidade do ar (HR). Figura 3.25. Intensidade real e modelada do incêndio de Figueiredo-Sertã das 15 horas do dia 15 às 4 horas do dia 16 de outubro. Repare-se na evolução da humidade do ar (HR). A calma súbita que antecede o “colapso” da pluma constitui uma indicação clara de downdraftgeneralizado nos incêndios convectivos (ROTHERMEL 1991), quando o indraft cessa antes da aparição do downdraft. Essa diminuição momentânea da velocidade do vento antes do seu pico máximo é claramente visível nos dados da estação de Nelas (Figura 3.26) e do parque eólico de Fajão, Pampilhosa da Serra (Figura 3.27), situado na região onde obtivemos testemunhos pessoais e observámos evidências de downdraft. No caso de Nelas, a intensidade mínima do vento coincide com o início da subida da humidade relativa do ar, condição necessária para o processo na análise de radiosondagem. 010203040506070010.00020.00030.00040.00050.00060.00015 16 17 18 19 20 21 22 23 24 1 2 3kw/m meteokw/m realHR0102030405060010000200003000040000500006000070000800009000013 14 15 16 17 18 19 20 21 22 23 24 1 2 3 4kw/m meteokw/m realHRII SÉRIE-E — NÚMERO 16_____________________________________________________________________________________________________________________80
Figura 3.26. Velocidade do vento máxima e humidade relativa do ar (intervalos de 10 minutos) na estação meteorológica de Nelas (IPMA) durante os dias 15 e 16 de outubro de 2017. A elipse assinala a coincidência entre a velocidade mínima do vento, que antecede os valores máximos, e o início da subida da humidade do ar. Figura 3.27. Velocidade do vento (médias de 10 minutos) ao nível das turbinas no parque eólico de Fajão, Pampilhosa da Serra, durante o dia 15 de outubro e até às 8 horas do dia 16 de outubro de 2017. A linha resulta de um spline cúbico ajustado ao conjunto de observações individuais. Veja-se a redução que antecede os valores máximos registados durante o período. Conclui-se que se verificaram duas tipologias de PyroCb no dia 15 de outubro, no litoral dominada pelo vento (potência do vento Pw superior à potência do fogo Pf, de acordo com a conceptualização de Byram) e no interior dominada pela energia do incêndio ao final da tarde (Pw● Ventos fortes determinaram o crescimento rápido do incêndio; ● A posterior redução do vento permitiu maior verticalidade da pluma do incêndio; ● Os incêndios litorais foram dominados por um low jet; ● A chegada de uma massa de ar mais húmida permitiu a transição de PyroCu a PyroCb, primeiro no litoral e depois no interior. 17 DE ABRIL DE 2018_____________________________________________________________________________________________________________________81
3.5.4.3. Zona de piro-convecção A piro-convecção situou-se no centro do país, sobre os incêndios de Vilarinho-Lousã, Figueiredo-Sertã e Esculca-Côja-Arganil. As três colunas exibem o fenómeno por volta das 17 horas, sendo evidente a existência de interações entre elas e a criação de uma zona de tempestade de fogo. Com a diminuição do vento as três colunas dispersam-se, e quando chega a noite a grande área em combustão mantém o LCL elevado e próximo dos 4 km. O início da recuperação de humidade favorece agora o crescimento vertical da convecção livre, com o posterior downdraft ou “colapso” de coluna. O mapa do índice de aerossóis do satélite Suomi NPP mostra que a concentração e espessura do fumo se dava em altitude (Figura 3.28). Figura 3.28. Índice de aerossol medido pelo satélite Suomi NPP, com valores elevados no Norte de Portugal e sobre a zona de ventos fracos que são indicativos das plumas de fumo dos incêndios. A concentração elevada de fumo dos incêndios, porém sem dispersão, confirma a piro-convecção.  https://ozoneaq.gsfc.nasa.gov/data/omps/ 3.5.4.4. Sequência espácio-temporal da piro-convecção A análise do momento piro-convectivo da tarde de 15 e noite de 16 de outubro permite descrever mais detalhadamente a sequência de acontecimentos. Tal como descrito e agora explicitamente apresentado na Figura 3.29 nem todos os incêndios coincidem no seu crescimento, havendo que procurar a explicação nas suas localizações no terreno e interações. Para isso utilizámos a análise do fator dominante do incêndio, usando o conceito de Byram e Rothermel para classificar a força dominante, isto é, a energia do vento (Pw) ou a energia do incêndio (Pf). A Figura 3.30mostra esta evolução para os fogos da Sertã, Arganil, Pataias e Quiaios. Esta análise confirma os momentos distintos de PyroCb no litoral e no interior, ao mesmo tempo que se reforça a ideia das diferentes interações intuídas na reconstrução da piro-convecção, assim como na simulação de cenários. Na Figura 3.29 observam-se quatro grandes momentos no desenvolvimento dos eventos, e que o vento do Ophelia está muito longe de ser o único criador do megaevento de incêndios. O Ophelia iniciou a expansão dos incêndios, mas a posição dos mesmos, o estado do combustível e as interações entre colunas convectivas acabaram por gerar uma sequência complexa que o Quadro 3.4. descreve. II SÉRIE-E — NÚMERO 16_____________________________________________________________________________________________________________________82
Figura 3.29. Crescimento horário dos mega-incêndios estudados. No eixo das abcissas a hora acumula desde as 9 horas do dia 15 de outubro. Figura 3.30. Evolução da relação entre a potência do vento (Pw) e a potência do fogo (Pf) para quatro dos incêndios. No caso do complexo de Arganil-Seia apresenta-se também a relação a 1500 m de altitude; note-se o domínio de Pf cerca das 20 horas, momento de PyroCb e downdraft nesse incêndio. 17 DE ABRIL DE 2018_____________________________________________________________________________________________________________________83
Quadro 3.4. Fases e controlos da expansão dos mega-incêndios em análise e interação entre eles, diferenciando o litoral do interior. Hora Fase Litoral Interior 11:00 a 14:00 Condução pelo vento O vento domina a propagação e gera incêndios compridos.  14:00 a 15:00 O vento domina a propagação e gera incêndios compridos.  No incêndio de Pataias Pw>Pf, porém em Quiaios Pw e Pf já são equivalentes.  Arganil e Lousã desenvolvem os primeiros PyroCb “à sombra” da coluna de convecção da Sertã. O indraft assim gerado alonga o fogo da Sertã em sua direção, o que mantém o incêndio com Pw próximo de Pf e com severidade muito elevada. 16:00 a 19:00 1a vaga de PyroCbcom vento a baixar O vento perde força à superfície e mantém valores equivalentes de Pw e Pf. A coluna ascende e há verticalização da piro-convecção, mas a situação quebra aos 5000 m de altitude e origina muitas projeções, com incêndios alongados (Figura 3.30). É o mais baixo PyroCu detetado pelo IPMA (6500 m). Situação idêntica ao litoral. O vento enfraquece e as colunas verticalizam-se até que quebram em altitude. Os incêndios alongam-se por focos secundários. Forte convectividade sob a coluna e 1.o episódio generalizado de PyroCu transitando para PyroCb, pequenos e separados. 19:00 a 23:00 2a vaga de PyroCb com a chegada de humidade As colunas geram PyroCu, mas não PyroCb e dispersam-se com a desaparição do low jet. A humidade relativa aumenta e a situação desagrava-se pontualmente. O fogo é conduzido pelo vento, com projeções a longa distância. A coluna de Pataias adquire piro-convectividade com o acréscimo de humidade, mas não a coluna de Quiaios, à “sombra” da 1.a. A entrada de humidade permite que os incêndios acedam à altitude de convecção livre na atmosfera com vento já reduzido. Criam-se grandes PyroCb sendo Pf o dobro de Pw. Os incêndios continuam a gerar indrafts poderosos e interagem entre si. Sertã contínua debaixo da influência de Lousã e Arganil, e o último exibe o seu máximo crescimento. Lousã, em posição mais litoral, continua a alongar-se com Pw ligeiramente acima de Pf. Das 19 às 20 horas os três incêndios têm atividade piro-convectiva. É uma situação baseada em outflows repetidos e sustentados. 00:00 a 04:00 Colapso final Fogo totalmente conduzido pelo vento sob humidade alta.  Ligeira rotação para NE. Lousã colapsa sobre todo o seu flanco leste. Este incêndio foi o 1.o a surgir e manteve piro-convecção elevada com Pw similar a Pf. É o último PyroCb a norte e desenvolve o seu máximo crescimento para NE. O seu downdraft faz com que Arganil passe a ser conduzido pelo vento (Pw>>Pf). O máximo Pw de Arganil dá-se nesta fase, quando o vento geral era mínimo – é o vento do downdraft.  II SÉRIE-E — NÚMERO 16_____________________________________________________________________________________________________________________84
Foto de Hélio Madeiras Figura 3.31. Coluna de convecção geral do incêndio de Pataias, inclinada com o vento, mas erguida. O fumo ascende verticalmente sob a coluna em situação de piro-convecção que se desenvolve protegida do vento. Esta situação ocorre quando as potências do vento e do fogo se (quase) equivalem.  3.5.5. Simulação do complexo de incêndios A simulação dos incêndios permite validar as sequências de factos identificados, e projetar diferentes cenários para o futuro, assim como comprovar a capacidade de previsão destes eventos de que se dispõe com a tecnologia atual. Utilizou-se para tal o software Wildfire Analyst (RAMIREZ et al. 2015), aplicado aos incêndios de Arganil, Lousã, Sertã, Quiaios e Pataias. A simulação teve que ser ajustada para os diferentes modelos de combustível, especialmente para aqueles de floresta, com fatores de 2,98 a 5,63. Este intervalo é coerente com a literatura referente a fogo de copas, que indica subestima por parte dos modelos do USDA Forest Service (ALEXANDER e CRUZ 2006), mas o grande número de colunas convectivas e as suas interações explica também a magnitude dos valores de ajustamento, que se deve portanto à grande capacidade de propagação por copas nestes incêndios sob alta instabilidade. A análise de sensibilidade dos dados permitiu observar que: ● Pataias. Simula-se melhor com os dados de Cabo Carvoeiro do que com a mais próxima estação, Leiria, que está “à sombra” do incêndio e cuja medição do vento é, portanto, afetada; ● Quiaios. Com os dados da estação mais próxima, Figueira da Foz, o incêndio não responde, exigindo fatores de ajustamento superiores a 8. Usando os dados de simulação de Pataias é possível alcançar 91% da realidade. Isto indica-nos que este incendio se situa sob a influência do 1.o e os dois alteram a atmosfera criando condições similares. De facto, a estação de Figueira da Foz é afetada pelo incêndio de Pataias como se observa na Figura 3.32. 17 DE ABRIL DE 2018_____________________________________________________________________________________________________________________85
Figura 3.32. Campo de ventos do dia 15 de outubro às 18:00 horas. A estação da Figueira da Foz é a única discordante e aponta ao quadrante sul. ● Sertã. Simula-se melhor com os dados da Lousã do que com os da zona. Basicamente tal indica que rapidamente ficou sob a influência do indraft do potente incêndio da Lousã. De facto, como se vê na Figura 3.33, a simulação com os dados teóricos da zona do incêndio fá-lo propagar-se para NE em vez de “procurar” o conjunto da Lousã e Arganil-Seia. II SÉRIE-E — NÚMERO 16_____________________________________________________________________________________________________________________86
Figura 3.32. Simulação sem interação do conjunto Lousã e Arganil-Seia (norte) e Sertã (sul), imagem superior. Observe-se como o fogo da Sertã se dirige para NE, o que não aconteceu na realidade. A simulação na imagem inferior é corrigida com os valores meteorológicos da Lousã assumindo que devido ao indraft estavam a afetar o incêndio a sul. ● Lousã. Inicialmente responde bem aos dados locais, mas a partir das 18 horas é a estação de Nelas aquela que melhor permite simular o incêndio. Nelas reflete o efeito dos downdrafts dos diferentes momentos do PyroCb, enquanto a Lousã permanece dentro do indraft. ● Arganil-Seia. Responde perfeitamente aos dados da Lousã primeiro e aos de Nelas mais tarde. É exatamente a mesma dinâmica e por tanto a mesma interação. Cabe somente destacar que a ignição de Vide responde ao outflow do poderoso downdrafta norte, determinante para o episódio que se gera no flanco esquerdo do incêndio da Lousã (registado por Nelas). Para simular esta parte haveria que usar a meteorologia das Penhas Douradas. Os mega-incêndios do dia 15 de outubro interagiram fortemente entre si devido à instabilidade atmosférica e todos exibiram momentos de piro-convectividade dominada pela pluma, superando a intensidade do vento. Esta interação extrema é confirmável através das modificações necessárias para ajustar a simulação á realidade, alterações que consistiram em atribuir a cada incêndio as condições de certos incêndios vizinhos. 17 DE ABRIL DE 2018_____________________________________________________________________________________________________________________87
Figura 3.33. Simulação ajustada dos 5 incêndios estudados mostrando o seu crescimento máximo possível durante 26 horas. Da simulação dos incêndios conclui-se que: 1. Os incêndios do litoral funcionaram como uma só coluna, movendo-se dominados pela primeira (Pataias) no fluxo de vento e criando “sombra” que distorceu os dados locais, os quais não explicam a propagação dos fogos. 2. Os incêndios do interior foram dominados pelo desenrolar do incêndio do primeiro incêndio, Lousã. A sua coluna convectiva dominou por interação o complexo de Arganil-Seia, exceto o incêndio de Vide, e “atraiu” o intenso incêndio da Sertã para si, o qual de outra forma se teria desviado para nordeste. 3. Estes dois complexos de incêndios formaram uma sombra sobre o Norte de Portugal, provavelmente reforçando a menor velocidade do vento na região, facto que ajuda a explicar a não formação de maiores incêndios nessa região. De facto, o conjunto de estações na projeção do PyroCb de Lousã-Arganil registaram à hora crítica uma inversão da direção do vento no sentido SE, onde estavam os incêndios (Figura 3.34). 4. Seriam expectáveis incêndios maiores no Norte de Portugal e especialmente na Galiza na ausência da extensa nuvem de convecção gerada pelo complexo interior de incêndios de Lousã, Sertã e Arganil. II SÉRIE-E — NÚMERO 16_____________________________________________________________________________________________________________________88
Figura 3.34. Campo de ventos do dia 16 de outubro de 2017 às 03:00 horas. Observe-se como as estações mais a norte indicam orientação para sudeste ou este, quando as restantes têm componente noroeste a nordeste. A fim de ter uma melhor sensação do significado deste episódio de incêndios efetuámos simulações adicionais para um cenário de ausência de precipitação e usando os dados de cada estação meteorológica. Os resultados e conclusões são: 1. Sem chuva no pós-Ophelia, 18 horas de simulação: área ardida de 504 137 ha. Sem a ocorrência de precipitação a área ardida duplicaria. 2. Sem chuva no pós-Ophelia, 48 horas de simulação: 532 010 ha. Para lá de 24 horas, e perdidas as condições de piro-convecção, seriam muito limitadas as possibilidades de crescimento dos incêndios (em valor relativo). 3. Sem interação: 364 048 ha. É a situação de área ardida mais mitigada, o que mostra claramente a importância das interações para explicar o crescimento extraordinário dos incêndios no Centro do país na noite de 15 para 16 de outubro. A análise do incêndio da Sertã leva a conclusões inquietantes ao comparar o crescimento ajustado à realidade com aquele que teria sido possível transcorridas 48 horas, e com a expansão máxima sob as condições criadas durante a interação e na ausência da limitação imposta pelos incêndios anteriores (de 2017). Neste último caso a superfície ardida ascenderia a 504 mil hectares. Desta simulação ressalta a importância das interações para originar episódios de tempestade de fogo e a gestão da paisagem (aqui como em incêndios anteriores) como a única ferramenta de limitação de mega-incêndios desta ordem de magnitude no futuro. 17 DE ABRIL DE 2018_____________________________________________________________________________________________________________________89
Figura 3.35. Simulação do incêndio da Sertã, ajustada à realidade (esquerda) e projetada para 48 horas (centro), e simulação livre (condições do dia 15 de outubro durante 48 horas) se o incêndio tivesse ocorrido sem que os incêndios de junho e agosto tivessem sucedido. 3.5.6. Severidade dos mega-incêndios A severidade do fogo avalia o grau de modificação na vegetação em consequência da passagem do fogo. Indica, portanto, o impacto agregado do fogo nos vários estratos de vegetação existentes, das árvores à manta morta sobre o solo; indiretamente indicará também o aquecimento potencial do solo atingido durante a combustão e a sua subsequente vulnerabilidade a fenómenos erosivos. Aqui consideramos a severidade do fogo dos sete mega-incêndios que ocorreram a 15 de outubro, sabendo desde já que apenas um não reflete um complexo de duas ou mais ignições, Vilarinho-Lousã. Globalmente, é manifesto o enorme impacto destes incêndios, uma vez que as classes de severidade Muito Elevada e Extrema perfazem ligeiramente mais de metade da área ardida total, sendo que a classe Extrema representa quase um terço do total (Quadro 3.5). As classes Baixa e Moderada constituem ainda assim um quarto da área ardida. A distribuição da severidade do fogo variou entre os incêndios. Embora a amostra seja insuficiente para detetar padrões consistentes são evidentes algumas associações: ● Albitelhe-Vouzela apresenta a maior % do somatório de severidades Baixa e Moderada (1/3 do total), para o que pode ter contribuído a % comparativamente elevada de combustíveis “jovens” e elevada recorrência do fogo; ● Sabugueiro-Seia é distinto na grande prevalência de severidade Extrema, sendo-o também no que respeita à ocupação por matos; ● Figueiredo-Sertã corresponde à segunda mais elevada fração de área nas classes Muito Elevada ou Extrema, conjugando % intermédias de ocupação por matos e por combustível “velho”, bem como grandes áreas de pinhal jovem e homogéneo; ● A severidade em Vilarinho-Lousã é apenas ligeiramente inferior à de Figueiredo-Sertã, sendo também o incêndio com maior representatividade de espaços florestais com combustível mais “velho”. Quadro 3.5. Distribuição (%) da área ardida por classe de severidade nos mega-incêndios alvos de estudo, designados pelo nome do município em que se deu a ignição. Classe de severidade Alcobaça Figueira da Foz Lousã Sertã Arganil Seia Vouzela Total Baixa 7.98 11.73 9.60 7.62 10.36 7.02 14.70 9,71 Moderada 16.83 17.28 13.86 13.14 22.25 7.35 18.81 16,22 Elevada 26.78 24.27 19.88 21.87 27.55 11.96 21.10 22,58 Muito Elevada 17.62 20.70 17.05 27.07 18.74 12.13 15.70 19,03 II SÉRIE-E — NÚMERO 16_____________________________________________________________________________________________________________________90
Extrema 30.80 26.02 39.61 30.30 21.09 61.54 29.69 32,47 Figura 3.36. Relação entre a área ardida com severidade Extrema (% do total) e a ocupação do solo por eucalipto (R2=0,76, p=0,0109) e por matos (R2=0,76, p=0,0027) para os mega-incêndios individuais. Procurámos analisar a relação entre o nível de severidade do fogo e os vários descritores de uso do solo e combustibilidade. As únicas relações significativas constam da Figura 3.36 e mostram diminuição da severidade do fogo para maiores taxas de ocupação por eucalipto, e o seu aumento quando os matos estão mais representados. No primeiro caso, a explicação deverá residir no efeito combinado destes fatores: ● Gestão ativa mais provável, implicando operações episódicas ou frequentes de redução do combustível; ● Povoamentos tendencialmente mais jovens, portanto com menor carga de combustível; ● Estrutura do povoamento menos conducente a fogo de copas, isto é, com maior descontinuidade vertical e horizontal. No caso dos matos há essencialmente um estrato de combustível, cuja remoção pelo fogo é proporcionalmente mais completa do que em formações florestais por a fração de biomassa fina ser largamente predominante. Acresce que a folhada sob matos é frequentemente incipiente, portanto totalmente consumida pelo fogo e expondo mais solo mineral, e que estas formações frequentemente revestem áreas rochosas. 17 DE ABRIL DE 2018_____________________________________________________________________________________________________________________91
Figura 3.37. Classificação da severidade dos incêndios nas matas nacionais costeiras: Buronhosa-Pataias (esquerda) e o conjunto formado por Quiaios-Figueira da Foz e as ocorrências de Vagos (direita). Os mapas de severidade do fogo (Figuras 3.37-3.40) permitem apreciar a respetiva variação espacial em cada um dos casos. Em incêndios desta natureza é expectável que na sua grande maioria as manchas de severidade mais baixa estejam confinadas a pequenas áreas ditadas pelas circunstâncias de desenvolvimento do fogo. Nos incêndios de 15 de outubro observa-se que as manchas menos severamente afetadas correspondem tipicamente: ● Às áreas de desenvolvimento inicial do incêndio; ● Às orlas dos incêndios, especialmente nos flancos e retaguarda (onde a propagação se dá contra o vento); ● A efeitos topográficos locais de (des)alinhamento da força e direção do vento com a morfologia do terreno que quebram a propagação do fogo nas encostas mais inclinadas ou pelo contrário o “canalizam” através dos vales e poupam o terreno elevado. É expectável também uma elevada associação entre as fases mais ativas de expansão dos incêndios e os níveis mais elevados de severidade, uma vez que o contributo do estrato arbóreo para a progressão do fogo (portanto de copas) aumenta durante esses períodos. Esta ligação é bastante evidente nos vários mega-incêndios, exceto na serra da Estrela, dada a já referida predominância de matos. Da mesma forma, as áreas que arderam de forma mais dispersa e difusa por eclosão de focos secundários tendem a apresentar menor severidade, o que resulta também da fragmentação do espaço florestal. II SÉRIE-E — NÚMERO 16_____________________________________________________________________________________________________________________92
Figura 3.38. Classificação da severidade dos mega-incêndios de Vilarinho-Lousã (esquerda) e Figueiredo-Sertã (direita). Destacamos alguns aspetos relevantes e particulares relativamente aos vários mega-incêndios em análise: 1. Burinhosa-Pataias. A área ardida em 2003 (~2500 ha) ardeu uniformemente com severidade Extrema e Muito Elevada formando uma mancha distinta pela continuidade e extensão, presumivelmente devido à estrutura de combustível, contínua e com elevada carga. São visíveis linhas de severidade mais baixa que não correspondem a atributos observáveis no terreno e que muito provavelmente refletem os flancos de propagação das duas ignições das quais resulta a área queimada e, eventualmente, projeções que se tenham desenvolvido significativamente. 2. Quiaios-Figueira da Foz. O sinal é pouco evidente comparativamente a Burinhosa-Pataias, mas a área correspondente ao grande incêndio de 1993 (3000 ha) ardeu agora com mais severidade, dada a baixa estatura do pinhal e a elevada quantidade de matos e acácias no sub-bosque. Também são observáveis algumas linhas de menor severidade resultantes dos padrões de progressão e alguns efeitos sombra de redução da severidade do fogo causados por ilhas não ardidas (por descontinuidade do combustível) individualizando-se claramente uma área queimada em 2013 que o incêndio flanqueou. 17 DE ABRIL DE 2018_____________________________________________________________________________________________________________________93
Figura 3.39. Classificação da severidade dos mega-incêndios de Esculca-Côja-Arganil e Sandomil-Seia (esquerda) e Sabugueiro-Seia e Casal Boavista-Folgosinho-Gouveia (direita). 5. Vilarinho-Lousã. Sendo o maior mega-incêndio e um dos que globalmente teve maior impacto na vegetação é também aquele onde se observa a maior alternância entre classes de severidade a escalas espaciais mais pequenas. Tal é particularmente notório no ramo esquerdo do incêndio a norte da barragem da Aguieira e explica-se pela heterogeneidade da ocupação do solo e pela grande importância da propagação por projeções. Os espaços florestais são razoavelmente fragmentados por áreas agrícolas e sociais, e sendo o eucalipto predominante é também de esperar heterogeneidade e um mosaico nas condições de acumulação e estrutura do combustível. 6. Figueiredo-Sertã. Grandes áreas contínuas de severidade Muito Elevada ou Extrema caracterizam este incêndio, particularmente a norte do rio Zêzere, mas também nas cumeadas que o antecedem a sul, revestidas continuamente por pinhal jovem, eucaliptal e, mais a norte, matos. No canto sudoeste, queimado pela ignição de Figueiredo-Sertã, a severidade foi mitigada pela presença de áreas agrícolas e aldeias. Após a cumeada do Fajão, no sul da serra do Açor, observa-se uma área longitudinal em que a severidade se reduziu substancialmente, correspondendo à expansão do fogo no sentido descendente das encostas e com o seu abrandamento, havendo bastante maior heterogeneidade espacial da vegetação e presença significativa de folhosas caducifólias; esta situação repete-se em duas áreas junto ao término norte do incêndio. 7. Esculca-Côja-Arganil. O padrão de severidade deste complexo de incêndios mostra a sul uma vasta área de severidade continuamente Muito Elevada ou Extrema, interrompida localmente nos espaços agrícolas em torno das aldeias e correspondente à serra do Açor; uma mancha distinta de 146 ha junto a Vale de Maceira e à qual foi atribuída severidade reduzida corresponde a um incêndio de 13 II SÉRIE-E — NÚMERO 16_____________________________________________________________________________________________________________________94
de agosto de 2017 que muito provavelmente não voltou a arder. Severidades mais moderadas têm alguma expressão no sector sudoeste, a sul de Côja, e coincidem grosso modo com ocupação por floresta caducifólia. De Oliveira do Hospital para norte o padrão de severidade tem alguma semelhança com aquele descrito para a generalidade de Vilarinho-Lousã. Entre Tábua e Oliveira do Hospital é visível uma extensão alongada de menor severidade, aparentemente resguardada por aldeias e campos agrícolas. 8. Sabugueiro-Seia e Casal Boavista–Folgosinho-Gouveia. Tal como já referido este conjunto distingue-se pela grande prevalência da classe de severidade Extrema. Severidades nas classes reduzida e moderada ocorreram localmente em afloramentos rochosos com ou sem pinhal disperso, mas também na imediação de aldeias e em pequenas manchas de caducifólias. A exceção mais relevante à regra geral é, contudo, uma mancha relativamente grande perto de Salgueiro de Baixo que ardeu com severidade apenas. Elevada e que havia sido percorrida por fogos em 2011, 2012 e 2013. 9. Albitelhe-Vouzela. A forma da área queimada e a direção conhecida dos ventos sugere fortemente que se trata de dois fogos, no entanto apenas há registos de ocorrências no sector ocidental. Neste a severidade é geralmente bastante inferior relativamente ao sector oriental, onde é similar à dos incêndios do ponto anterior, predominando largamente a classe Extrema. A menor severidade na parte ocidental é explicável pela grande fragmentação do espaço florestal, especialmente nas imediações de Oliveira de Frades. Tal como nos incêndios de Vilarinho-Lousã e Esculca-Côja-Arganil, nomeadamente nas suas secções norte, trata-se de uma tipologia de território nunca antes percorrido pelo fogo nesta escala. Pelo contrário, o nível de severidade atingido na parte leste é surpreendente, já que uma parte substancial tinha ardido há menos de cinco anos. Figura 3.40. Classificação da severidade do mega-incêndio de Albitelhe-Vouzela. 17 DE ABRIL DE 2018_____________________________________________________________________________________________________________________95
RELATÓRIO  |  COMISSÃO TÉCNICA INDEPENDENTE  visíveis faixas de gestão de combustível, mosaicos de parcelas orientadas linearmente e parcelas individuais (Figura 3.41). A escala da gestão reflete o constrangimento imposto pela pequena dimensão da propriedade, incluindo a dificuldade de aceitação de intervenções de silvicultura preventiva que acarretem perda de rendimento. Nas ZIF da CAULE o esforço de intervenção é algo inferior à recomendação de tratamento anual de 5% do território para efeitos de DFCI, sendo que apenas uma parte das intervenções são de gestão do combustível. Ainda assim foi conseguido localmente algum grau de sobreposição das faixas de gestão de combustível, essencial para perturbar a propagação de um grande incêndio e oferecer oportunidades consecutivas de combate. O potencial de crescimento de um grande incêndio era muito elevado na serra do Açor, dado que as últimas manchas queimadas significativas datavam de 1991, 1995, 1996 e 2005, dependendo da localização. Assumindo um esforço anual de tratamento de 3%, efetivo durante quatro anos, para as áreas ocupadas por formações arbustivas imediatamente a norte dos locais de reacendimento do incêndio referenciado como Esculca-Côja-Arganil, a sua velocidade de expansão seria reduzida entre <5% e 30% à escala da paisagem, dependendo do posicionamento mais aleatório ou mais estratégico das intervenções. A taxa de intervenção era, portanto, insuficiente para perturbar (passivamente) o potencial de propagação, ainda que localmente (nas áreas tratadas) o fogo pudesse abrandar, aliás como foi reportado nas imediações de Moura da Serra pelas equipas da CAULE, embora não seja suportado pelas isócronas de desenvolvimento do fogo (Figura 3.13). Nessa fase a gestão de combustíveis terá criado oportunidades efetivas de combate, mas o comportamento do fogo exigia meios de combate suficientes e adequados no local; refira-se que nas circunstâncias do dia 15 de outubro o ATI teria que ser bastante reforçado para ser eficaz. Logo, o efeito das áreas intervencionadas no incêndio adviria fundamentalmente do seu aproveitamento pelas forças de combate, o qual não aconteceu. Os próprios recursos de apoio ao combate a incêndios da CAULE, depois de numa fase inicial estarem dedicadas ao combate ao incêndio florestal, tiveram de redirecionar a ação para a sua autodefesa e defesa das populações, na situação que se viveu de simultaneidade de situações de socorro e de ausência de apoio de outras entidades de proteção civil. Este caso ilustra claramente a complementaridade e sinergia entre prevenção e combate. Figura 3.41. Sobreposição entre intervenções silvícolas efetuadas no quadrante sudeste das ZIF da CAULE (2011-2016) e classes de severidade do fogo. Informação de base fornecida pela CAULE. 17 DE ABRIL DE 2018_____________________________________________________________________________________________________________________97