Acompanhamento e Projeção de óbitos pela Covid 19 no Brasil

 

Economia e Energia  

Nº 106, janeiro a junho de 2020
ISSN 1518-2932 Disponível em: http://ecen.com  
(números anteriores) http://ecen.com

ACOMPANHAMENTO E PROJEÇÃO DE ÓBITOS PELA COVID NO BRASIL 
Redação Preliminar em Discussão E&E 106

Baixar o arquivo em pdf

Palavra do Editor

Razões para a divulgação desse Acompanhamento

Um conhecimento só é realmente útil quando possibilita uma “previsão” do futuro.

Prever o futuro no curto e médio prazo é uma tarefa ingrata porque os pretensos profetas serão punidos também no curto e médio prazo.

Cientes desse risco, resolvemos, mesmo assim, ir divulgando os resultados parciais de nossas projeções e questões na esperança que eles possam ser úteis.

Na realidade, como não se sabe o comportamento que será adotado por governos e cidadãos quanto ao isolamento social, não podemos ter uma ideia precisa do que acontecerá. O que a metodologia realmente fornece é uma ideia para onde estamos indo. Para onde realmente vamos é uma questão onde todos somos os protagonistas, conscientes ou não disto.

Ou seja, os autores deste trabalho já temos a desculpa padrão para nossas prováveis falhas.

Acreditamos que o exame dos dados por diferentes ângulos, mesmo por pessoas treinadas no uso das estatísticas, mas sem conhecimento específico do problema médicos, pode ser útil aos especialistas e a população interessada para compreender as possibilidades e os riscos envolvidos. 

Estamos apresentando, nessa versão preliminar, os resultados comentados para o Rio de Janeiro e Brasil. Temos os mesmos dados disponíveis para as 5 Macrorregiões e para os estados.  Nosso objetivo é de irmos aprofundando e corrigindo dados e observações ao longo da pandemia. Estamos naturalmente abertos a colaborações externas e ao debate,

Carlos Feu Alvim

Texto e dados em discussão, redação preliminar de 12/06/2020

COVID 19

Abordagem logística da propagação do Novo Coronavirus

Carlos Feu Alvim, Olga Mafra e José Israel Vargas

O tratamento logístico foi aplicado por Cesare Marchetti[1] a inúmeros sistemas sociais e econômicos e tecnológicos. No Brasil, o professor José Israel Vargas o aplicou também a diversos sistemas, com ênfase em casos brasileiros.

Aqui mesmo, já aplicamos esse ajuste para vários casos e o assunto já foi tratado no número zero dessa revista,[2] em artigo de Omar Campos Ferreira, e uma visão completa sobre a metodologia está no exemplar de número 45 da E&E, em artigo de José Israel Vargas “Prospectiva Tecnológica, Previsão com um Simples Modelo Matemático”[3].

A metodologia, que nasceu na análise de Volterra e Lotka de sistemas de competição biológica, pode ser aplicada à Covid 19 que está se colocando como a maior pandemia desde a gripe espanhola, ocorrida há praticamente um século (1918 a 1920).

Fundamentalmente, essa metodologia se baseia em um modelo no qual a taxa de infecção em uma população depende do número de infectados N e do número dos que restam a infectar (N*-N) sendo N* o número final de atingidos, denominado de nicho. Matematicamente isso se traduz na equação na qual, ao longo do tempo o número de contaminados será proporcional ao produto de um termo crescente N e outro decrescente
N*-N.

T = dN/dt = a.N.(N*-N) ou

T = a.N*.N – a.N2       (1)

Onde dN/dt pode ser associado ao número de novos infectados (variação de N) por unidade de tempo (tomada aqui como um dia).  É fácil mostrar que o valor máximo da taxa de contaminações acontecerá quando N = N*/2. Na equação de T reconhecemos a equação de uma curva de segundo grau do tipo  
               [y = a1.x – a2.x2],     
mais precisamente, uma parábola. Um ajuste pelo método de mínimos quadrados nos permite estimar o número final de infectados N*; ao atingir metade desse nicho estaremos no máximo de contaminações por dia.

Supondo que o número de óbitos seja uma fração fixa do número de atingidos, o mesmo tipo de equação pode ser usado nos dois casos, tomando dados parciais disponíveis para estimar a quantidade de infeções ou de óbitos finais. Como o número de óbitos tende a ser mais bem avaliado, procuramos nos concentrar aqui nesses dados.

Um problema prático, no caso, é que as políticas de distanciamento social reduzem o tamanho deste nicho. Se essa política é mudada, o nicho também muda de tamanho. Outro problema, no caso do Brasil que tem dimensões continentais, é que a epidemia vai se difundindo aos poucos pelo País e alcança cada município e até cada bairro em diferentes datas. Outro problema prático é a necessidade de uma população relativamente grande para que o tratamento estatístico faça sentido. Nossa escolha foi tratar os dados em três níveis: país região e estados.

Da resolução da equação diferencial dN/dt = a.N.(N-N*) resulta a equação:

F = N/N* = 1/(1+ e(-at-b))                    (2)

Onde F é a fração do “nicho” ocupado e a função encontrada é da chamada curva em S, denominada de logística.

Tomando-se o logaritmo neperiano da equação encontra-se

Ln(F/(1-F)) = a.t + b              (3)

Essa equação facilita estimar os coeficientes a e b da função logística ajustando uma reta nos valores de Ln(F/(1-F))

Aplicação da metodologia aos números do Estado do Rio de Janeiro

Os dados, fornecidos pelo MS – Ministério da Saúde, para o Rio de Janeiro são aqui apresentados como exemplo de aplicação da metodologia. Eles se referem ao número de pessoas infectadas e no de óbitos registrados, em consequência da Covid 19. Como já mencionamos, preferimos concentrar a análise nos dados de óbitos, não obstante o retardo implícito das ocorrências e seu registro em relação à data de contaminação.

Ver e baixar arquivo E&E 106 parcial

Aplicação da metodologia aos números do Brasil

A aplicação da metodologia tem suas limitações em um país das dimensões do Brasil onde existem, vários sistemas em ambientes diferentes com uma grande variação de latitude, clima e ocupação do espaço físico.

Ver e baixar arquivo E&E 106 parcial

[1] Cesare Marchetti. Society as a learning system. Technological Forecasting and Social Change, 18:267-282, 1980

[2] https://ecen.com/eee0/eeezero.pdf

[3] http://w.ecen.com/eee45/eee45p/prospeccao_tecnologica.htm

[4] Média quadrática dos desvios relativos

[5] Naturalmente, os últimos três dados ainda estariam sujeitos às variações semanais porque média é centrada (o próprio dia e três anteriores e três posteriores); assim os três últimos dados é feita a correção da “semanalidade” .

[6] O Ministério da Saúde – MS tem divulgados os dados por 100 mil habitantes. No caso, preferimos adotar as mortes por milhão de habitantes que favorecem a comparação diária e não dificultam a acumulada.

[7] J.C. Fisher and R.H. Pry. A simple substitution model of technological change; Technological Forecasting and Social Change, 3:75-88, 1971. https://www.sciencedirect.com/science/article/pii/S0040162571800057

Desflorestamento da Amazônia

Desmatamento da Amazônia ao longo do tempo

Matéria em discussão:                

NOVA ONDA DE DESFLORESTAMENTO NA AMAZÔNIA 
O QUE DIZEM OS DADOS OBJETIVOS

Carlos Feu Alvim e Olga Mafra
(autoria provisória dessa edição preliminar)

Resumo

Em 2012 e 2017, em artigos sob a coordenação de José Israel Vargas, apresentamos a análise e revisão da evolução dos dados de desflorestamento anuais na Amazônia, usando uma modelagem logística de Volterra. A conclusão da primeira análise na revisa E&E 86 é que o desflorestamento caiu, entre 2004 e 2012 , mais que o previsto pela expectativa histórica. Já na revisão  na E&E 95, constatou-se que, quase simultaneamente com a Conferência do Clima de Paris (COP 21 em 2015), iniciava-se uma aparente retomada do desflorestamento da Amazônia, levantando inquietudes sobre o comportamento futuro. Nas duas ocasiões anteriores, foi feita uma projeção da trajetória de desmatamento esperada dentro do modelo adotado. Nesta revisita, a extrapolação anterior é mantida e comparada com os novos dados disponíveis. Foi confirmada a retomada do desflorestamento da Amazônia em um nível que ainda pode ser considerado como uma oscilação na tendência de queda, mas que pode evoluir para um pulso de desmatamento como já aconteceu anteriormente. A intenção deste trabalho é contribuir com a racionalidade dos dados para uma discussão que está tomando rumo considerado destrutivo para as legítimas aspirações e interesses do Brasil e do próprio Planeta.

Palavras Chave:

Amazônia, desmatamento, desflorestamento, modelagem logística, análise quantitativa.

____________

1       Introdução

O desflorestamento[1]  da Amazônia no Brasil foi objeto de análise nesta revista, sob a coordenação de J. I. Vargas, na E&E 86 (Vargas , et al., 2012)  e revista na E&E 95 (Vargas, et al., 2017). Nelas foi feita a avaliação da evolução desse desflorestamento com auxílio de uma modelagem matemática simples, desenvolvida por Volterra para tratar a ocupação de um nicho ambiental por um organismo vivo.

Essa abordagem tem sido aplicada com êxito por C. Marchetti e J. I. Vargas em diversos fenômenos complexos envolvendo variáveis econômicas e sociais. Os detalhes sobre a metodologia estão amplamente descritos nos dois trabalhos anteriores.

Na primeira avaliação (setembro de 2012), os dados cumulativos representados em um gráfico mostravam uma curva ao longo do tempo em “S”, indicativa de tendência de estacionamento da área desmatada. Na segunda avaliação (junho de 2017), foi confirmada a tendência de contenção do desmatamento e constatada frenagem nesse processo mais rápida que a esperada. Essa queda no desmatamento, entre 2004 e 2012, teria resultado do reconhecimento da importância política do problema. Isso se refletiu nos compromissos internacionais assumidos pelo Brasil, reforçado na Conferência do Clima de Paris.

Na ocasião da segunda avaliação, os resultados dos anos anteriores já indicavam uma retomada do desmatamento e o artigo se encerrava com a questão: Uma nova onda de desflorestamento?

Agora, em 2019, uma grande polêmica sobre o assunto se instalou e nos estimulou a voltar a tratar do tema. A polêmica se estabeleceu baseada na divulgação de dados de alerta mensais que, isoladamente, são inadequados para avaliar tendências de desmatamento. É bom realçar que o gráfico de desmatamento anual já apresentara outros picos e vales no processo que causaram alarme na época, mas não reverteram a tendência de longo prazo que era a redução do desmatamento anual[2].

Para entender a base da polêmica criada, cabe esclarecer que os dados divulgados mensalmente pelo Instituto Nacional de Pesquisas Espaciais – INPE são de natureza diferente dos dados anuais. O dado mensal é um alerta de desmatamento baseado em imagens por satélite que assinalam áreas com indícios de desmatamento. O dado anual é de acréscimo da área desmatada e reflete dados extraídos de imagens de satélite, verificados estatisticamente no campo para evitar problemas de má interpretação.

 Deve-se realçar que a experiência acumulada no processo de análise e verificação no local vai tornando essa análise preliminar cada vez mais confiável. Existe outro aspecto a ser considerado nos dados mensais, que é a cobertura de nuvens que prejudica a coleta de dados em alguns meses. Isso pode fazer com que os resultados de um mês acumulem variações que não puderam ser computadas no mês ou meses anteriores, onde predominava a escassa visibilidade.

 Ao final do artigo anterior (E&E 95, abril-junho 2017) havia um chamado à atenção para os sinais de retomada do desmatamento amazônico. O destaque a seguir transcreve nossa advertência:

“Este comportamento deve ser monitorado para prevenir eventuais pressões que deslocariam o equilíbrio alcançado e poderiam desencadear um novo ciclo de desmatamento semelhante ao ocorrido nas últimas décadas que teve seu auge em torno da virada do século.

 Esta atenção é necessária já que, os resultados para os dois últimos anos (2015 e 2016) mostram aumentos de áreas desflorestadas nesses anos na Amazônia Brasileira. Recentemente, foram propostas medidas de afrouxamento das restrições ao desmatamento, como a da MP, aprovada pelo Senado em 25 de maio de 2017 e vetada pelo Presidente Temer, que alterava os limites da Floresta Nacional do Jamanxim. Isso pode sinalizar uma tendência de negligenciar restrições ao desmatamento e comprometer as metas voluntárias apresentadas pelo Brasil em Paris, dando início a um novo ciclo de desflorestamento na Amazônia.”

J. I. Vargas, R. Grandsire e C. Feu Alvim em E&E 95

 Esta nota mostra que esta tendência de aumento de desmatamento está confirmada, mas ainda pode ser considerada como uma oscilação dentro da tendência de longo prazo.

2        Os Polêmicos Resultados de 2019

Os últimos resultados apresentados pelo INPE de alertas de desmatamento mostravam um aumento de quase 100% na área atingida, como indicado na Figura 1. Esse aumento surge da comparação dos dados de junho de 2019 com junho de 2018.

Um comentário apressado do Presidente da República a uma pergunta de jornalista desencadeou uma série de discussões que atingiram escala internacional nas quais os dados históricos foram quase esquecidos.

Figura 1: Comparativo do alerta de desmatamento do mês de junho com o mesmo mês de anos anteriores na publicação do site do INPE[3]

Na mesma forma de apresentação da Figura 1, o INPE colocou em seu site gráficos similares, abrangendo mais meses, os quais também indicavam o crescimento do desmatamento, mas em uma intensidade muito menor. No ambiente de “discussão de bar” que tomou o País e o Mundo essas “tecnicalidades” foram também desprezadas.

O INPE fornece ainda “alertas de desmatamento + degradação” que inclui áreas onde há indicações de redução da densidade florestal. Essa área correspondente a “alertas de desmatamento + degradação” atingiu, no mês de junho, 2072 km2. Esse dado também aparece na discussão na mídia, confundido com o de desmatamento.

É um erro relativamente comum confundir áreas para as quais são emitidos alertas de desmatamento com a área de desmatamento propriamente dita. Só a verificação local, realizada com critérios estatísticos adequados, permite avaliar a área de desmatamento. Esse é o dado que interessa do ponto de vista de danos à natureza. Esse procedimento que alia os dados de satélite a verificações no local, com auxílio de outras entidades como Embrapa e Ibama, é que torna os dados brasileiros particularmente confiáveis. Dados como o espectro de luz emitido ajudam a aperfeiçoar o processo.

A razão pela qual são divulgados os alertas é que eles constituem um instrumento muito útil na prevenção, mitigação e repressão do desmatamento. Eles são muitas vezes capazes de indicar áreas onde a ação predadora do homem deva ser reprimida e áreas onde o período de seca favorece o início da combustão e sua propagação.

A inclusão dos alertas de degradação aperfeiçoou a capacidade de detecção da extração prévia de madeira de lei ou de outras ações lesivas à floresta como, por exemplo, seria a aplicação de produtos químicos que facilitam o desmate. Ambos os alertas aparecem na imprensa confundidos com desmatamento. Também ambos estão sujeitos a incertezas, mas os de degradação são mais sutis e, em consequência, mais sujeitos a erros.

No site do INPE, são mostrados os gráficos, para quatro anos, dos alertas de desmatamento não só para o mês de junho (Figura 1), como para o trimestre (abril a junho) e os ocorridos após o início da estação de queimadas (mês de agosto do ano anterior até junho do ano corrente). Os dados estão mostrados na Tabela 1.

Tabela 1: Alertas de Desmatamento do INPE, por Períodos, em km2 , de 2016 até junho de 2019

 

Desmatamento

Alertas de Desmatamento

 

ANUAL

JUNHO[4].

ABR-JUN

AGO-JUN

2016

7393

969

1800

     4639

2017

6947

610

1099

     4182

2018

7536

480

1528

     3976

2019

 

920

1907

     4575

 
       

Fonte: site do INPE visitado em agosto de 2019

O aumento de 2018 para 2019 que foi de 92% (quase 100%) para junho, foi de 25% para o trimestre e apenas 15% para os 11 meses de agosto a junho. Ou seja, se os alertas indicam, para todos os períodos, um aumento de desmatamento em 2019, sua dispersão indica que a magnitude do crescimento da taxa de desmatamento está longe de ser estabelecida. Além disto, ainda resta a outra metade do ano para apurar.

A Tabela 1 mostra ainda que o dado de junho de 2018 foi excepcionalmente baixo para o mês e o valor deste ano apenas o segundo mais alto nos 4 anos. Ou seja, a alta taxa de crescimento de alertas para junho deste ano (relativa a junho do ano passado), parece corresponder muito mais ao baixo valor observado no ano passado (480 km2) que ao valor alto de 2019 (920 kn2), como pode ser visto na Tabela 1[5].  

O uso dos alertas para avaliar o desmatamento parte do princípio de que “onde tem fumaça tem fogo”. Mas, ainda usando essa analogia, só a presença dos bombeiros vai indicar se realmente houve o sinistro ou se o alarme, ou o alerta, era falso. A hipótese de que os alertas possam ser usados para avaliar o efetivo desmatamento pode ser testada usando os valores da Tabela 1 para os outros anos.

Isso foi feito com auxílio da Tabela 2 que usa a taxa de crescimento nos períodos para fazer “previsões” de variação anual de desmatamento. Esses valores podem ser comparados com a variação efetivamente estimada para o ano (primeira coluna). As previsões de desmatamento anual a partir dos dados intermediários de alerta são, nesses dois anos da tentativa, decepcionantes.

Tabela 2: Variação das Taxas de Desmatamento por Período

 

Taxa
Apurada

Taxas Estimadas a partir de dados por período

 

ANUAL

JUNHO

ABR-JUN

AGO-JUN

2017

-6%

-37%

-39%

-10%

2018

8%

-21%

39%

-5%

2019

 

92%

25%

15%

Se tomados os alertas de junho como se fossem a efetiva indicação do desmatamento anual, o valor para junho de 2017 (-37%) nos levaria a comemorar uma queda de 37% no desmatamento enquanto a queda efetiva no ano foi de apenas 6% (primeira coluna). Também em 2018, tomando os dados de junho, estaríamos comemorando uma queda de 21% quando, na verdade, a apuração anual revelaria um aumento de 8% no desmatamento.

A diferença para o ano de 2019, é que todas as projeções indicam uma maior taxa de desmatamento, ao contrário do que aconteceu nos anos anteriores, onde os sinais para os diferentes períodos eram contraditórios. Isso é um forte indício de que realmente deve haver um aumento de desmatamento no ano em curso, que refletirá no índice anual, a menos que medidas drásticas sejam adotadas no restante do ano.

O valor parcial que melhor reproduz a variação anual parece ser, como é natural, o do maior período. O problema é que ele inclui dados do ano anterior. Levando isto em conta, escolhemos usar, como prévia do desmatamento em 2019, o aumento de 25% da taxa de alertas de desmatamento do trimestre. A faixa de erro estaria em cerca de +/- 30% o que dá ideia da incerteza existente nestes dados.

O ensaio aqui realizado sugere a necessidade de um trabalho sistemático, sobre os dados mensais de alarme, ao longo do ano, para buscar estimar antecipadamente o desmatamento anual. Este indicador serviria para orientar as decisões políticas. Os dados mensais de alarme permitiriam determinar a projeção do valor do ano em curso. Na apresentação dos dados, deveria ser indicada a margem de erro esperada. Isso facilitaria as interpretações e evitaria, ou pelo menos atenuaria, conclusões apressadas.

Para os dados econômicos mais críticos o IBGE adota a política de processar os dados em sigilo e divulgar, com uma antecedência de algumas horas, os resultados para a autoridade econômica correspondente. Como o cronograma é anunciado com antecedência, os técnicos dos organismos econômicos responsáveis pela área têm a capacidade de, em algumas horas, oferecer aos jornalistas a interpretação do que é divulgado. No caso, o sigilo e o curto tempo de avaliação evitam interferências políticas que possam macular os resultados. O exemplo poderia ser adotado para dados ambientais que apresentam sensibilidade semelhante.

Ou seja, pode-se ter armado “uma tempestade em copo d’água” em torno de um resultado intermediário com pouca significação prática. Paradoxalmente, esse destaque exagerado que se deu ao episódio pode servir agora para chamar a atenção para uma retomada do desmatamento, menos espetacular que a suposta no recente episódio, mas que já se mantem por quatro anos. Isso será visto no próximo item.

3       Incorporação dos dados 2017 a 2019

O trabalho anterior (E&E 95) reviu as projeções para os próximos anos. Nesta nova análise, vamos manter essa projeção. Isto é, aliás, útil para testar a confiabilidade das projeções feitas. Naquele trabalho foram considerados os dados até 2016.

Nesta revisão, foi possível atualizar o dado de 2016 e incorporar os de desmatamento de 2017, 2018 e 2019 (valor preliminar para o último ano). Os dados sobre o desmatamento da Amazônia Legal, levantados pelo INPE e pela EMBRAPA[i], são mostrados na Figura 2.

A comunidade científica nacional e internacional, têm uma atenção especial para a evolução do desflorestamento na região amazônica. Nas avaliações feitas pelo Governo, o desmatamento, chegou a ser a principal responsável pela contribuição brasileira aos gases de efeito estufa. Além de sua contribuição para o aquecimento global, estima-se que a destruição da floresta tenha influência marcante no clima do continente como um todo e, particularmente, no regime pluvial de nosso país.

No gráfico da Figura 2, constam os dados a partir de 1988. O tratamento anual do desmatamento na Amazônia vinha sendo predominantemente qualitativo, até 1989, quando se iniciou a fotointerpretação de imagens obtidas pelo satélite Landsat 5, pelo INPE. Esse trabalho contribuiu para a elaboração da Comunicação Nacional para a Convenção Quadro das Nações Unidas sobre Mudança do Clima da Comunicação do Brasil, em 1994[ii]. Segundo este relatório, a contribuição do desmatamento, ocorrido até então (470 mil km2), responderia por cerca de 50% das emissões de gases causadores do efeito estufa pelo Brasil no período.

Figura 2: Evolução do desmatamento anual da Amazônia Legal

Com os dados anuais e a estimativa do desmatamento até então, foi possível avaliar, no trabalho anterior, a evolução do desmatamento histórico acumulado na Amazônia Legal.

A Figura 3 mostra a evolução do desmatamento acumulado e o ajuste realizado para a descrição da trajetória provável ao longo do tempo. O desmatamento evoluiria para atingir uma área total de 891 mil km2 da Amazônia Legal. Esse valor máximo (nicho) é um dos parâmetros principais a ser ajustado. O outro parâmetro (delta) relaciona-se com o tempo entre atingir 10% de desmatamento e 90%, estimado em 47 anos. Ambos os fatores são indicados na Figura 3. O eixo da esquerda representa a fração F do nicho ocupado e no eixo à direita está o valor do desmatamento acumulado.

A metodologia de estabelecer a curva de tendência é simplesmente aplicar aos dados disponíveis o melhor ajuste para o tipo de função matemática utilizada (logística).  

Figura 3: Desmatamento acumulado na Amazônia Legal e ajuste baseado em uma curva logística, mostrando a tendência de se limitar a expansão do desmatamento.

Pelo ajuste, o total desmatado da região amazônica até 2018 é de 796 mil km2, correspondentes a 19,4% da área total da Amazônia (estimada em 4,109 milhões de km2) [iii]. A previsão da área final desmatada corresponde a 22% da área da Amazônia.

Dentro dessa perspectiva de médio prazo os dados sobre o desflorestamento mostram uma evolução relativamente positiva já que haveria uma preservação de quase 80% de uma floresta nativa o que é excepcional em termos mundiais. Já a ação predatória sobre a biodiversidade ainda permanece largamente ignorada. Isto se deve não só à enorme variedade e complexidade desse bioma, mas também aos ainda limitados esforços científicos realizados para elucidá-la.

Pode-se observar que os dados reais apresentam uma oscilação em torno do ajuste e estaríamos passando de um período em que a evolução real estava abaixo da tendência para uma onde haveria uma fase um pouco acima da previsão.

O fenômeno dessas oscilações em torno da média já foi muito bem explorado por C. Marchetti e J. I. Vargas conforme comentado no artigo anterior.

Essas oscilações ficam mais visíveis quando comparadas com na representação de Fisher Pry mostrada na Figura 4. No caso, o ajuste foi automaticamente refeito quando foram introduzidos os dados dos últimos anos. O ajuste obtido é, no entanto, muito próximo do anterior. A Figura 4 ilustra também a metodologia usada para o ajuste e projeção.

Figura 4: Representação Fisher-Pry da função logística do desmatamento da Floresta Amazônica (1977 – 2018).

A curva mostrada nas Figuras 3 e 4 corresponde aos dados acumulados de desmatamento, ou seja, à integral dos dados de acréscimos anuais de desmatamento mostrados na Figura 2, acrescida de uma estimativa do valor até a data inicial. A curva das variações, dita diferencial, é mostrada na Figura 5, juntamente com o ajuste usado no artigo anterior (dados até 2016).

A tendência observada na Figura 3 parece coerente com uma estabilização na área coberta pela Floresta Amazônica dentro do prometido pelas autoridades brasileiras na Conferência de Paris. Apenas o tempo de estabilização é um pouco maior. Uma melhor ideia do que significa esse ajuste, em termos de desflorestamento, surge da Figura 5 que representa os valores anuais e o ajuste[6].

Figura 5: Comparação do desflorestamento anual verificado com o ajuste logístico

A Figura 5 mostra o comportamento da taxa de desmatamento, permitindo distinguir as variações ocorridas, muito sujeitas a circunstâncias econômicas e políticas, sobrepostas a uma tendência de longo prazo que passou por um máximo em 1992 e se aproximaria do zero menos drasticamente do que se poderia desejar, mas mais coerente com o processo normal de evolução de fenômenos sociais complexos.

Vale a pena analisar os picos e vales do processo e suas circunstâncias: O ano de 1991 é o segundo do Governo Collor que apresentou forte recessão no primeiro ano e foi precedido por grave crise econômica e financeira. Os dados da Figura 5, não mostram o possível pico no desmatamento. O ano de 1995, pico de desmatamento, foi o primeiro ano do Governo FHC e o ano 2004 o segundo do Governo Lula tendo sido precedido por profunda crise financeira na transição de governos. Estamos assistindo, não há mais dúvida agora, a uma retomada do desmatamento que já vem do governo anterior e estamos ainda imersos em uma profunda crise econômica; o início do desmatamento coincide praticamente com o da crise. Tentativas de associar atividade econômica e desmatamento da Amazônia também foram feitas como parte de nossas análises anteriores, com êxito limitado.

4       Enfrentando a Crise

Mesmo tendo sido provavelmente provocada pelo uso equivocado do indicador, houve uma escalada de declarações internas e externas que já configuram a Crise.

A atual retomada do desmatamento iniciou-se em 2015 e tem aparente seguimento em 2019. Até agora, os dados existentes deixam dúvida se assistimos uma simples retomada da trajetória natural indicada neste trabalho pela curva “de ajuste” ou se tratar de um novo ciclo de desmatamento aproveitando o vácuo de poder que ocorre frequentemente em mudanças de governo ou mesmo a expectativa de afrouxamento da política de contenção do desflorestamento.

O alarme atual em torno do tema não se deve a uma indicação clara fornecida pelos índices, mas a uma interpretação precipitada da comparação de um mês de junho com índice anormalmente baixo de 2018 com um possível índice alto do mesmo mês em 2019. As amostragens de outros períodos não confirmam a possibilidade de ser real um crescimento de quase 100% no desmatamento anual como aventado.

A situação criada tem motivações políticas no Brasil e no exterior, agravadas por erros de comunicação do Governo. Com efeito, existe, do lado do Governo Federal, um discurso que pode ser interpretado como de perda de importância da atitude ambiental considerada “politicamente correta”. Do outro lado, existe a legítima preocupação, interna e externa, com a preservação da natureza e uma somatória de interesses políticos da oposição de contrariar o governo, além de interesses externos de disseminar a incapacidade do Brasil de controlar a proteção da Amazônia.

Os dados ambientais necessários à defesa da posição brasileira existem, mas foram colocados em dúvida pelas discussões dentro do próprio Governo. Esse desentrosamento também dificulta a tomada de medidas efetivas de repressão ao desmatamento, já complexas pela falta de recursos para a verificação nos locais indicados pelo sistema e a própria contenção dos incêndios.

Os dados históricos disponíveis a respeito do desmatamento ao longo de três décadas mostram um comportamento que é típico de situações em que o equilíbrio se estabelece por um jogo de forças e interesses. A tendência de longo prazo ao longo de 30 anos é a de êxito na contenção do desmatamento amazônico como assinala o ajuste mostrado na Figura 5 e anteriores,

Sempre existe o perigo real de que um descuido institucional na defesa das florestas possa vir a causar uma grave perda da cobertura florestal como nos episódios de 1995 e 2004, ambos em início de governo.

Por outro lado, os compromissos brasileiros na Conferência de Paris que prometem, para o final de 2030, zerar a contribuição florestal à emissão de CO2 podem não ser realizáveis.

Pode-se perceber na Figura 5 que a tendência projetada não prevê o desmatamento zero para a Amazônia. Medidas de recomposição de florestas na própria Amazônia e em outros biomas poderiam compensar as emissões e possibilitar o alcance das metas. Isso quase certamente não acontecerá em um quadro de recessão econômica e em uma atividade onde a reposição, mesmo na presença de recursos é lenta. Não interessa também a ninguém que as metas globais de emissão no Brasil se cumpram a custa da recessão prolongada que compense as prováveis emissões florestais.

Na mesma revista em que foi feita a análise do desflorestamento (E&E 95), também foram apontadas as previsíveis dificuldades do cumprimento das metas estabelecidas, sobretudo se forem criadas condições de uma nova expansão da economia que é absolutamente vital para o bem-estar e a paz social da Nação. Vale repetir as conclusões do artigo “As Metas brasileiras de emissões de Gases de Efeito Estufa e a Contribuição Nacionalmente Determinada – CND do Brasil”:

“A intensidade de emissões (relativas ao PIB) no Brasil já se encontra em um patamar muito baixo em razão da presença dos renováveis. Manter os atuais coeficientes já é um desafio para muitos setores. Não parece racional a passividade nos setores produtivos na aceitação da redução de emissões adicionais em alguns destes itens. Talvez muitos acreditam que recursos externos ou do próprio governo resolverão os problemas: isto é certamente uma ilusão.

As metas que se transformaram em compromisso, são ambiciosas como solicitadas aos países (pela ONU). As metas Emissões/PIB setoriais podem se revelar incompatíveis com o crescimento. A recente tentativa de modificação da legislação sobre o uso da terra na Amazônia é talvez a primeira reação organizada de setores econômicos contra medidas associadas às emissões de GEE. Seria melhor para o conceito do País que as metas merecessem uma discussão mais profunda antes de serem assumidas.”

Carlos Feu Alvim e Olga Mafra na E&E 95

Infelizmente, isso não ocorreu. A atitude das forças produtivas brasileiras foi a “alegre aceitação” sem contestação dos compromissos assumidos. Deve-se reconhecer que as metas foram fixadas sem o necessário diálogo com os setores interessados. Isso aconteceu não obstante as tentativas do Ministério de Meio Ambiente de propiciar os foros de discussão que esbarraram na incapacidade de diálogo fundamentado, devido ao escasso tempo disponível e da falta da preparação técnica que possibilitasse uma discussão fundamentada com os setores produtivos.  

Para atenuar as possíveis preocupações, alegou-se que estavam sendo fixadas apenas intenções de contribuição voluntária e não compromissos. A realidade está demonstrando que mesmo compromissos voluntários assumidos pelo País acabam se tornando obrigações passíveis muitas vezes de justificar retaliações.

Outro fator a considerar é que mesmo sendo mantidos os compromissos, são esperadas oscilações em torno da trajetória para alcançar as metas estabelecidas. Em um sistema de equilíbrio dinâmico de forças socioeconômicas, nem sempre é possível atingir as metas. Isto deve ser encarado de forma natural e a apuração constante dos progressos alcançados e as justificativas do que está acontecendo auxiliam a alcançar este equilíbrio. Outros países deverão rever suas metas e isto pode ser feito em um clima de concertação.

5       Conclusão

A situação presente encerra ameaças ao conceito do País e pode acarretar danos à Soberania Nacional e deverá ser enfrentada com a seriedade implícita na ameaça existente e preferentemente com recursos próprios. A entrada desordenada de equipes de “salvamento” da Amazônia pode produzir danos importantes à visão de nossa capacidade de lidar com o problema. A eventual ajuda externa é bem vinda, desde que sob o controle nacional.

Deve-se reconhecer que o Governo parece ter compreendido a gravidade da situação, embora ela tenha sido ativada com o que nos parece um falso gatilho. Isso motivou um pronunciamento presidencial e reuniões de coordenação no mais alto nível governamental.

É muito importante que seja mantida a credibilidade de nossas instituições encarregadas da apuração e verificação dos compromissos assumidos e dos mecanismos, técnicos ou de força que tem o Estado para corrigir os eventuais desvios. Do contrário, seria propiciar aos críticos externos novos motivos de alimentar suspeições e fortalecer sua cobiça. Nossos órgãos ambientais competentes e independentes fazem parte da primeira linha de nossa defesa. Essa consciência deve existir tanto nesses órgãos como nas demais autoridades do Governo e, sobretudo, na sociedade nacional.

6      Bibliografia

Vargas , J. I. e Gorgozinho, P. M. 2012. Modelagem Matemática Simples do Desmatamento da Amazônia. Economia e Energia E&E. [Online] Setembro de 2012. http://eee.org.br/wp-content/uploads/2017/02/eee86p.pdf.

Vargas, J. I., Grandsire, R. e Alvim, C Feu. 2017. Acompanhamento da evolução do eesflorestamento da Amazônia usando modelagem matemática simples. Economia Energia – E&E 95. [Online] setembro de 2017. http://ecen.com.br/wp-content/uploads/2017/06/eee95p.pdf.

[i] PRODES/IMPE-EMBRAPA Taxas anuais do desmatamento – 1988 até 2016: Taxa de desmatamento anual (km2/ano) em http://www.obt.inpe.br/prodes/prodes_1988_2016n.htm

[ii] Governo Brasileiro: Comunicação Nacional para a Convenção Quadro das Nações Unidas sobre Mudança do Clima, Brasilia, 1994. em

http://sirene.mcti.gov.br/documents/1686653/1706391/205854.pdf/5eadb8ca-f316-49ec-9dd1-7ba80754b20d

[iii] Marta Salomon e Tânia Monteiro, Maior parte de área desmatada da Amazônia virou pasto, Estado de São Paulo, 03/09/2011 citando relatório Instituto Nacional de Pesquisas Espaciais (Inpe) e pela Empresa Brasileira de Pesquisa Agropecuária (Embrapa) apresentado ontem ao Palácio do Planalto no dia anterior.

http://www.estadao.com.br/noticias/geral,maior-parte-de-area-desmatada-da-amazonia-virou-pasto-mostra-estudo-imp-,767823

_________________

Notas

[1] O termo desflorestamento é mais apropriado para o que se quer  medir na Amazônia, por força de uso, o termo desmatamento, correntemente usado no mesmo sentido, será muitas vezes usado neste trabalho em substituição ao desflorestamento.

[2] Os dados recentes aqui analisados parecem, no entanto, não se enquadrar inteiramente como oscilações de curta duração, como as anteriormente observadas, e devem ser acompanhados com atenção.

[3] http://www.obt.inpe.br/OBT/noticias/alertas-do-deter-na-amazonia-em-junho-somam-2-072-03-km2

[4] Dados referentes ao mês de junho estimados a partir da Figura 1.

[5] Essa dificuldade reforça o problema de se tomar um mês isolado para avaliar a variação de desmatamento. O INPE costuma advertir em todos os seus boletins de desmatamento para um problema adicional que é a forte influência das condições de visibilidade sobre o valor da avaliação mensal.

[6] Corresponde á Figura 2 com a curva de ajuste.

Desflorestamento da Amazônia


Economia e Energia Nº 95, Abril a Junho de 2017 – Ano XXI
ISSN 1518-2932

ACOMPANHAMENTO DA EVOLUÇÃO DO DESFLORESTAMENTO DA AMAZÔNIA USANDO MODELAGEM MATEMÁTICA SIMPLES

José Israel Vargas, Rafael Grandsire e Carlos Feu Alvim

Resumo

Repete-se aqui a análise realizada em 2012 dos dados de desflorestamento  na Amazônia usando uma modelagem logística de Volterra e aplicada extensamente por Cesare Marchetti e José Israel Vargas. A queda no ritmo do desflorestamento pós Copenhagen superou as expectativas da análise anterior. Nos últimos anos, entretanto, houve uma retomada do desflorestamento levantando uma inquietude sobre o comportamento futuro.

Palavras Chave:

Amazônia, desmatamento, desflorestamento, modelagem logística, análise quantitativa.

1       Introdução

O desflorestamento de modo geral, mas em particular o da Amazônia, segue sendo o foco principal das preocupações ambientais, no Brasil. Sua evolução foi abordada no N° 86[i] dessa revista com auxílio de uma modelagem matemática simples como um fenômeno que estaria em fase de saturação, assumindo a forma de uma curva em S e, tendendo a um valor um pouco inferior ao previsto inicialmente quando não havia compromisso explícito brasileiro com a contribuição do País em conter o efeito estufa.

O Brasil se comprometeu na iNDC (Intended Nationally Determined Contribution), cujo título completo em português é “Pretendida Contribuição Nacionalmente Determinada para Consecução do Objetivo da Convenção Quadro das Nações Unidas sobre Mudanças do Clima” [ii] a:

“- fortalecer políticas e medidas com vistas a alcançar, na Amazônia brasileira, o desmatamento ilegal zero até 2030 e a compensação das emissões de gases de efeito de estufa provenientes da supressão legal da vegetação até 2030”;

Este comportamento deve ser monitorado para prevenir eventuais pressões que deslocariam o equilíbrio alcançado e poderiam desencadear um novo ciclo de desmatamento semelhante ao ocorrido nas últimas décadas que teve seu auge em torno da virada do século. Esta atenção é necessária já que, os resultados para os dois últimos anos (2015 e 2016) mostram aumentos de áreas desflorestadas nesses anos na Amazônia Brasileira. Recentemente, foram propostas medidas de afrouxamento das restrições ao desmatamento, como a da MP, aprovada pelo Senado em 25 de maio de 2017 e vetada pelo Presidente Temer, que alterava os limites da Floresta Nacional do Jamanxim. Isso pode sinalizar uma tendência de negligenciar restrições ao desmatamento e comprometer as metas voluntárias apresentadas pelo Brasil em Paris, dando início a um novo ciclo de desflorestamento na Amazônia. Essa região tem importante papel, tanto no que diz respeito às mudanças climáticas, a nível global, quanto na presumida extinção da imensa biodiversidade regional que, estima-se, seja a mais rica do planeta.

2       Incorporação dos dados 2008 a 2015

Os dados do trabalho anterior alcançavam o ano de 2007 e nesta revisão, foi possível incorporar dados de desflorestamento até 2016 (valor preliminar para o último ano). Os dados sobre o desmatamento da Amazônia Legal, levantados pelo INPE e pela EMBRAPA[iii], são mostrados na Figura № 1 e mostram uma sensível redução nos últimos sete anos. O dado para 1988 representa a média avaliada de 1977 a 1988 e o último ano mostrado (2016) é uma avaliação preliminar.

Figura 1: Evolução do desflorestamento anual da Amazônia Legal

Existe, na comunidade científica nacional e internacional, uma atenção especial com a evolução do desflorestamento na região amazônica que, segundo avaliações feitas pelo governo, chegou a ser a principal responsável pela contribuição brasileira aos gases de efeito estufa. Além de sua contribuição no aquecimento global, estima-se que a destruição da floresta tenha influência marcante no clima do continente como um todo e, particularmente, no regime pluvial de nosso país.

O tratamento destes efeitos vinha sendo predominantemente qualitativo, até 1989, quando se iniciou a fotointerpretação de imagens obtidas pelo satélite Landsat 5, pelo INPE – Instituto Nacional de Pesquisas Espaciais. Esse trabalho contribuiu para a elaboração da Comunicação Nacional para a Convenção Quadro das Nações Unidas sobre Mudança do Clima da Comunicação do Brasil, em 1994[iv]. Segundo este relatório, a contribuição do desmatamento, ocorrido até então (470 mil km2), responderia por cerca de 50% das emissões de gases causadores do efeito estufa pelo Brasil.

P. Aguiar et al.[v]., em trabalho conjunto do Instituto Nacional de Pesquisas Espaciais – INPE e do Museu Paraense Emílio Goeldi – MPEG, estimaram as emissões de CO2 por desmatamento na Amazônia Brasileira utilizando um modelo matemático que combina mapas anuais de novas áreas desmatadas com informações espacialmente explícitas sobre a distribuição da quantidade de biomassa nos diferentes tipos de vegetação do Bioma Amazônia. O modelo leva também em conta as diferenças inter-regionais em termos de biomassa, com a localização dos desmatamentos, e considera suas diferentes etapas, a saber: o corte, a queima no solo, a decomposição das raízes e as queimadas nos anos sucessivos, subtraindo-se o que foi acumulado com o crescimento e corte da vegetação secundária.

Segundo os dados anuais do Governo Federal, disponíveis anualmente a partir de 1988, e considerando-se a estimativa do total acumulado até 1994, acima mencionada, o total desmatado da região amazônica até 2015 seria de 777 mil km2, correspondentes a 18,9% da área total da Amazônia (estimada em 4,109 milhões de km2) [vi]

Usando dados até 2015, a estimativa da área desflorestada ao final do processo é de cerca de 0,88 milhões de km2. Este valor resulta do melhor ajuste com os últimos dados disponíveis a partir da metodologia mostrada em edição anterior que usou dados até 2007 (E&E № 86). A estimativa deste trabalho, mostrada a seguir, é inferior à apresentada naquela edição que era de 1,11 milhão de km2. Nesta reavaliação, foram também corrigidos dados do desmatamento acumulado.

Se os dados sobre o desflorestamento mostram uma evolução positiva, a ação predatória sobre a biodiversidade ainda permanece largamente ignorada. Isto se deve não só à enorme variedade e complexidade desse bioma, mas também aos ainda limitados esforços científicos realizados para elucidá-la.

3       Discussão

Parece conveniente recordar que a metodologia descrita no trabalho anterior deriva do tratamento matemático dado à Teoria da Evolução e logo, à competição entre espécies vivas, por Vito Volterra, em seu famoso trabalho intitulado “Leçons sur la théorie mathématique de la lute pour la vie”[vii]  e extensamente explorado por Cesare Marchetti no International Institute for Applied System Analysis[viii]. Aqui usaremos somente uma solução particular da equação diferencial não linear de Volterra, a chamada versão Malthusiana. Ela descreve a dinâmica da ocupação, por predadores, de uma área limitada fisicamente, em recursos de todos os tipos (ou em extensão), definindo assim o chamado nicho. Detalhes dessa abordagem encontram-se nas referências[ix].

A equação Malthusiana (1), apresentada abaixo:

dN/dt = aN.(N* – N)                    (1)

mostra que a intensidade da ocupação (destruição) da floresta, por unidade de tempo dt (anual), é proporcional ao que já foi anteriormente desmatado, N, no tempo t e ao que resta a ser desmatado, N*- N, onde N* é o nicho total desflorestável no início do processo.

O a é uma constante de proporcionalidade, que mede a taxa de desflorestamento. A solução da equação (1) é conhecida como equação logística (ou epidemiológica, por descrever corretamente a evolução temporal de epidemias)[x] e é expressa a seguir:

N(t) = N*/[1+exp(at+b)]             (2)

Chamando de F a fração de ocupação do nicho N*, ou seja, F = N/N*, ao correr do tempo t a equação logística assume a forma:

F/(1-F) = exp(at+b)                    (3)

sendo b um localizador temporal, como se pode ver em sua representação dita de Fisher-Pry, obtida extraindo-se o logaritmo natural de (3) e consequentemente linearizando a equação:

ln[F/(1-F)] = at + b                    (4)

As representações gráficas das equações estão mostradas mais adiante, para o caso em estudo. A equação (1) tem forma de sino e varia desde valores negativos próximos de zero (de fato, desde menos infinito), passa pelo valor máximo 0,5 do valor total, (considerado este igual a 1 ou 100%) e assume valores positivos próximos de 1 (de fato, valores positivos muito grandes), no seu limite superior. Ela é chamada de “curva da vida”, por descrever quaisquer processos evolutivos, que consistam de ações de nascimento, crescimento até um ápice, seguido por decréscimo e morte, com a cessação do processo; ou, alternativamente, descreve uma mutação inicial, a difusão competitiva (Darwinista) da nova espécie e seu crescimento, até à saturação do nicho, mediante sua ocupação completa por senectude final do referido organismo, ou do processo considerado, por falta de “alimento” no nicho; ou por pressão de forças socioeconômicas emergentes contrárias. No caso presente, supõe-se justamente que essas forças limitariam o nível de desmatamento tolerável que não seria o total da área original. Assim, o nicho a ser ocupado aqui considerado não é o total da área (cerca de 4 milhões de km2), mas estaria limitado a cerca de 1 milhão de km2.

No presente caso, para adequar os dados à representação das equações acima mencionadas, estimou-se inicialmente o desmatamento acumulado, considerando-se os dados mensais disponíveis e o valor acumulado para 1991. Em seguida, a partir de um valor tentativo do nicho N* a ser ocupado estimou-se a fração F=N/N* alcançada a cada ano, sendo N o total desmatado acumulado até o ano.

A equação (1) permite determinar avaliar o “nicho” a partir dos dados do desflorestamento anual (representação de dN/dt) e dos valores acumulados N de desflorestamento.

Figura 2: Determinação do Nicho N* usando-se ajuste para dados de dN/dt e N que correspondem ao desflorestamento anual e acumulado, respectivamente.

Os dados mostrados na Figura 2 permitem estimar o nicho em 882 mil km2. Uma metodologia alternativa (usada no trabalho anterior) é procurar o valor máximo ajustando um polinômio de segundo grau e identificar esse máximo com N*/2.

Coeficientes muitos semelhantes resultam do ajuste na representação Fisher Pry mostrada na Figura 3 que corresponde a representação da equação (4). O Melhor ajuste para a reta nesta representação determina o nicho de 891 mil km2 para o desflorestamento limite e uma constante de tempo (de 10% para 90% de ocupação do nicho) Δ=48 anos. A representação abaixo usa o valor do nicho de 882 mil km2, obtido usando a metodologia anteriormente mostrada e possui uma constante de tempo de 47 anos.

Figura 3: Representação Fisher-Pry da função logística do desmatamento da Floresta Amazônica (1977 – 2015).

A Figura 4 mostra os dados acumulados do desflorestamento em fração do nicho e em valores em km2.

Figura 4 Desflorestamento acumulado e ajuste mostrando a tendência à limitação do valor desflorestado da Amazônia Legal

A tendência observada parece coerente com uma estabilização na área coberta pela Floresta Amazônica dentro do prometido pelas autoridades brasileiras na Conferência de Paris. Uma melhor ideia do que significa esse ajuste em termos de desflorestamento anual é o mostrado na Figura 5 que representa os valores anuais.

Figura 5: Comparação do desflorestamento mensal verificado com o ajuste logístico

A Figura 4 fornece uma melhor ideia da tendência do desflorestamento. A Figura 5 mostra as variações ocorridas que estão naturalmente muito sujeitas a circunstâncias econômicas e sociais. O ano de 1991 é o segundo do Governo Collor que apresentou forte recessão no primeiro ano. 1995 foi o primeiro ano do Governo FHC e 2004 o segundo do Governo Lula. 1994 e 2004 foram anos de crescimento econômico acima de 5% ao ano. A análise dos dados não permite, no entanto, obter uma clara correlação entre o crescimento econômico. Sabe-se que as mudanças de governo, principalmente quando envolvem mudança de tendência política, costumam propiciar uma inércia nos mecanismos do Estado o que poderia reforçar os picos de 1995 e 2004, também relacionados aos da atividade econômica. Na circunstância vigente nestes dois anos, se confirmada a expansão do desmatamento, certamente deverá ser atribuída a controles mais frouxos uma vez que não estamos em período de expansão da atividade econômica.

O compromisso brasileiro da Conferência de Paris não fixa uma trajetória anual para o desflorestamento. A trajetória mostrada na Figura 5 pode ser considerada a “natural” para fenômenos complexos como o desflorestamento. Mesmo sendo mantido o compromisso, esperam-se oscilações em torno dessa trajetória. Em um sistema de equilíbrio dinâmico de forças socioeconômicas, isto pode ser encarado de forma natural e a apuração constante e a pronta publicação do que está acontecendo auxilia alcançar este equilíbrio.

4       Uma Nova Onda de Desflorestamento?

Pouco se comentou no Brasil sobre a retomada do desflorestamento verificada em 2015 e 2016 mostrada na Figura 5 previamente à viagem do Presidente Temer à Noruega. Só em 21/06/2017 a Folha de São Paulo [xi] comentou a tendência preocupante da política ambiental brasileira em virtude dos resultados para os dois últimos anos. A entrevista do ministro norueguês, na realidade, faz referências muito elogiosas aos resultados dos anos imediatamente procedentes. O episódio teve ainda desdobramentos em pronunciamentos oficiais durante a visita presidencial.

Um desdobramento da metodologia aqui apresentada já foi usado na análise de outros fenômenos de longo prazo, onde, como no aqui tratado, intervêm muitos fatores e ocorrem oscilações em torno do ajuste da logística.

As Figuras 3, 4 e 5 mostram que o desflorestamento real vem oscilando em torno de uma tendência de longo prazo de ocupação de um nicho que corresponderia a cerca de um quarto da floresta original.

A mesma metodologia aqui aplicada para aos dados globais (de longo prazo) pode ser aplicada a períodos mais restritos. Isso é feito dividindo-se o intervalo de tempo e ajustando uma curva logística a cada um dos ciclos observados.

A aplicação da metodologia ao ciclo observado entre 1992 a 2012 resulta em um ajuste muito bom como o mostrado na Figura 6.

Figura 6: Ajuste dos dados entre 1997 e 2012

O resultado do ajuste é que o ciclo estudado já estaria encerrado e o valor acumulado em 2012 já representava 97% do nicho e o desmatamento anual esperado no ciclo tenderia rapidamente para zero. A Figura 7 mostra o ajuste em relação as emissões no período 1997/2012 e as que seguiram nos anos 2012/2016 que parecem indicar o início de um novo ciclo ou oscilação.

Figura 7: Dados Observados do desflorestamento anual e ajuste para o período 1997 e 2012

Como pode ser observado, o ajuste logístico é capaz de oferecer uma boa descrição para o ciclo observado. No caso, não foi possível identificar um padrão nos ciclos anteriores que pudesse ser útil para a previsão de um próximo.

5       Conclusão

O desflorestamento da Amazônia Brasileira vinha sendo paulatinamente reduzido e constitui um fenômeno que pode ser, em sua tendência global, reproduzido no longo prazo por uma modelagem matemática simples. O processo está sujeito a variações anuais que precisam ser monitoradas. Uma reversão de tendência é possível se for modificado o atual equilíbrio socioeconômico que levou à Política Ambiental que, em termos de contenção do efeito estufa, foi resumida na “Pretendida Contribuição Nacional Determinada” mais conhecida pela sigla em inglês iNDC que o Brasil apresentou na Conferência do Clima da ONU em Paris em 2016 e que foi recentemente ratificada pelo Brasil.

O objetivo deste trabalho foi apenas modelar o que vem ocorrendo. A identificação de uma tendência de longo prazo e as oscilações identificadas podem subsidiar uma melhor compreensão do que vem ocorrendo. Uma política para a Amazônia, deveria explicitar os objetivos a serem alcançados de maneira que o País pudesse dispor de um planejamento onde os recursos naturais da região fossem usados para alcançar um equilíbrio realmente sustentável dos pontos de vista ecológico, social e econômico. Desse mapeamento é que deveria surgir o conjunto de metas. No caso do desflorestamento, as rígidas metas estabelecidas parecem satisfazer apenas um dos pontos de vista, o que, normalmente, conduz a não sustentabilidade.

 

Referências

[i] http://ecen.com.br/wp-content/uploads/2017/02/eee86p.pdf

[ii] GOVERNO BRASILEIRO. Pretendida Contribuição Nacional para Consecução do Objetivo da Convenção Quadro das Nações Unidas sobre Mudanças do Clima. Brasília, 2016 em

http://www.mma.gov.br/images/arquivo/80108/BRASIL%20iNDC%20portugues%20FINAL.pdf .

[iii] PRODES/IMPE-EMBRAPA Taxas anuais do desmatamento – 1988 até 2016: Taxa de desmatamento anual (km2/ano) em http://www.obt.inpe.br/prodes/prodes_1988_2016n.htm

[iv] Governo Brasileiro: Comunicação Nacional para a Convenção Quadro das Nações Unidas sobre Mudança do Clima, Brasilia, 1994. em

http://sirene.mcti.gov.br/documents/1686653/1706391/205854.pdf/5eadb8ca-f316-49ec-9dd1-7ba80754b20d.

[v] Aguiar, A. P.; Ometto, J. O.; Nobre, C.; Câmara, G.; Longo, K.; Alvalá, R., Araujo, R. J.; Soares, V.; Vareliano, D.; Almeida, C.; Vieira, I. e Almeida, A. Estimativa das Emissões de CO2 por desmatamento na Amazônia Brasileira; Relatório Técnico Sintético INPE, 2009.

[vi] Marta Salomon e Tânia Monteiro, Maior parte de área desmatada da Amazônia virou pasto, Estado de São Paulo, 03/09/2011 citando relatório Instituto Nacional de Pesquisas Espaciais (Inpe) e pela Empresa Brasileira de Pesquisa Agropecuária (Embrapa) apresentado ontem ao Palácio do Planaltono dia anterior. Em

http://www.estadao.com.br/noticias/geral,maior-parte-de-area-desmatada-da-amazonia-virou-pasto-mostra-estudo-imp-,767823

[vii] Volterra, V.; Leçons sur la théorie mathématique de la lute pour la vie. Paris, Gauthier-Villars, 1931.

[viii] Marchetti, C.; A Forecasting Model for Research and Innovation Activities in Selected Areas: A Support for Strategic Choices, 1991. International Institute of Applied Systems Analysis, Laxenburg, Austria.

[ix] Goel, N. S.; Maitra, S. C.; Montroll, E. W.; On the Volterra and Other Nonlinear Models of Interacting Populations. Reviews of Modern Physics, 1971, 43(2), 231-276.

Vargas, J. I.; A Prospectiva tecnológica: a previsão com simples modelo matemático. Economia e Energia, nº 45 (julho/agosto de 2004) e nº 46 (setembro/outubro de 2004) em http://ecen.com/eee45/eee45p/ecen_45p.htm

[x] Vargas, J. I.; Ferreira, O. C.; Alvim, C. F.; Corgozinho, P. M.; Tratamento Logístico das Ocorrências Anuais de Dengue no Rio de Janeiro (1985 – 2008); Economia e Energia, nº 71 (dezembro 2008/janeiro 2009). Em http://ecen.com/eee71/eee71p/dengue_no_rj.htm

[xi] Vemos tendências preocupantes, diz ministro norueguês do ambiente. Folha de São Paulo 21/02/2017 em http://www1.folha.uol.com.br/ambiente/2017/06/1894546-noruega-critica-politica-ambiental-de-temer-as-vesperas-de-visita-oficial.shtml