Vulnerabilidade
Territorial do Cerrado

Áreas protegidas sob pressão antrópica

Análise integrada da pressão antrópica sobre unidades
de conservação do bioma Cerrado.

Área de Cerrado

Divisas Estaduais

Unidade de Conservação

Uso do Solo (5km)

Uso do Solo (10km)

Contexto Regional (20km)

▲

🗺️ Explore o mapa

Clique em uma Unidade de Conservação para ver os detalhes

Vulnerabilidade Territorial do Cerrado

O Cerrado é o bioma com maior biodiversidade de savana do mundo e um dos mais ameaçados do planeta. Com uma área demarcada de 2.004.423,77 km², abriga nascentes das principais bacias hidrográficas brasileiras e uma fauna e flora de riqueza incomparável.

Apesar de sua importância, apenas 9,12% do bioma está sob proteção formal: as 615 Unidades de Conservação distribuídas pelas esferas federal, estadual e municipal cobrem 182.746,61 km² do território. Dentro dessas áreas protegidas, 71,95% ainda preserva vegetação nativa do bioma — evidência de que as UCs cumprem seu papel de contenção da conversão, mas protegem uma fração pequena demais do território.

Considerando o bioma inteiro, 50,85% da área total do Cerrado — 1.019.248,98 km² — ainda mantém cobertura vegetal nativa, sendo que 44,29% dessa vegetação — 887.758,84 km² — está completamente fora de qualquer proteção formal, configurando extensas lacunas de conservação em território altamente vulnerável ao avanço agropecuário e à expansão urbana.

Este projeto desenvolve um índice multivariado de pressão antrópica para as 615 UCs do Cerrado, integrando dados de uso e cobertura do solo, características geométricas das unidades e categorias de manejo. O objetivo é oferecer uma ferramenta analítica que permita identificar onde a vulnerabilidade é mais crítica, qual é a dinâmica de pressão em cada unidade e como essa pressão se distribui entre as esferas administrativas — informações essenciais para priorização de gestão, fiscalização e políticas de conservação.

O Índice de Pressão Antrópica

O índice é composto por sete variáveis que capturam dimensões complementares de vulnerabilidade. Cada variável foi normalizada entre 0 e 1 antes da combinação — exceto V7, que já está nessa escala por construção. O resultado final é um valor entre 0 e 1 por UC, onde valores próximos de 1 indicam maior pressão antrópica.

Classificação do índice:

Faixa	Intervalo
Muito Baixo	abaixo de 0,2
Baixo	entre 0,2 e 0,4
Moderado	entre 0,4 e 0,6
Alto	entre 0,6 e 0,8
Muito Alto	acima de 0,8

As variáveis de uso do solo utilizam uma ponderação que diferencia o tipo de pressão antrópica: pastagem recebe multiplicador 1,0 por ser o principal vetor histórico de conversão do Cerrado; agricultura recebe multiplicador 1,5 por seu impacto mais intenso sobre solo, água e biodiversidade; área urbana recebe multiplicador 2,0 por ser praticamente irreversível na escala de tempo relevante para conservação. Corpos d'água naturais e praias são feições do ecossistema e não entram no cálculo de pressão.

Variáveis do Índice

Variável 1 — Intensidade Antrópica no Entorno Imediato (peso: 25%)

Mede a pressão antrópica ponderada no buffer de 5km ao redor da UC em 2024. É a variável de maior peso no índice por capturar o estado atual da zona de maior impacto direto sobre a unidade — onde o efeito de borda é mais intenso, afetando temperatura, umidade, penetração de espécies invasoras e pressão sobre fauna e flora nativas. Valores altos indicam que o entorno imediato da UC está predominantemente ocupado por usos antrópicos, especialmente área urbana e agricultura intensiva.

Como interpretar: um valor de 0,65 no entorno imediato de 5km significa que a zona de maior vulnerabilidade ecológica da UC está fortemente comprometida, independente do que acontece mais longe.

Variável 2 — Intensidade Antrópica no Entorno Próximo (peso: 20%)

Aplica a mesma lógica de V1 no buffer de 10km em 2024. Captura a pressão numa escala espacial mais ampla, onde dinâmicas de expansão agropecuária e urbana se manifestam antes de atingir o entorno imediato. A diferença entre V1 e V2 é analiticamente relevante — uma UC com V2 muito maior que V1 está num processo de cerco progressivo, com a pressão ainda distante mas avançando. Uma UC com V1 e V2 igualmente altos já está completamente inserida numa matriz antropizada.

Como interpretar: valores similares entre V1 e V2 indicam pressão uniforme em todo o entorno; valores muito maiores em V1 do que em V2 indicam pressão concentrada no entorno imediato com algum amortecimento mais longe.

Variável 3 — Taxa de Conversão Recente (peso: 20%)

Mede a perda percentual de vegetação nativa do Cerrado no buffer de 10km entre 2015 e 2024 — uma janela de nove anos que captura pressão ativa e recente. É a variável mais importante para identificar UCs sob ameaça dinâmica, diferenciando unidades com entorno em processo de degradação acelerada daquelas cuja degradação já estava consolidada antes do período de análise. Casos de recuperação de vegetação — valores negativos — foram tratados como zero.

Como interpretar: é a variável que mais importa para avaliar tendência futura. Uma UC com V3 alto está sendo pressionada agora — o entorno que existia em 2015 está sendo convertido ativamente. Uma UC com V3 próximo de zero pode estar preservada ou pode já estar tão degradada que não há mais vegetação para perder — por isso V3 deve sempre ser lido em conjunto com V1 e V2.

Variável 4 — Isolamento Regional (peso: 15%)

Aplica a área antropizada ponderada no buffer de 20km em 2024, medindo o contexto paisagístico amplo ao redor da UC. UCs inseridas numa paisagem regional relativamente preservada têm maior resiliência, maior conectividade com outras áreas naturais e maior capacidade de recuperação após perturbações. UCs isoladas numa matriz antropizada estão ecologicamente ilhadas — mesmo que internamente preservadas, sofrem pressão crescente por fragmentação e perda de conectividade.

Como interpretar: V4 alto com V1 e V2 baixos indica uma UC ainda bem preservada internamente mas progressivamente cercada numa escala regional — sinal de alerta para o médio prazo.

Variável 5 — Vulnerabilidade Intrínseca por Tamanho e Forma (peso: 10%)

Combina duas características geométricas da própria UC que determinam sua vulnerabilidade independente do entorno. O componente de tamanho — peso 60% — reflete que UCs menores têm maior proporção de borda em relação à área interna, amplificando o efeito de borda e reduzindo a área de interior efetivamente protegida. O componente de compacidade — peso 40% — quantifica a irregularidade da forma através da razão entre área e perímetro ao quadrado, onde valores próximos de 1 indicam forma mais regular e compacta e valores menores indicam formas alongadas ou muito recortadas.

Como interpretar: uma UC pequena e muito irregular combina os dois fatores de vulnerabilidade — máxima proporção de borda e máxima exposição ao entorno. Esta variável é intrínseca à UC e não muda com o tempo, servindo como linha de base de vulnerabilidade estrutural sobre a qual as pressões externas se somam.

Variável 6 — Isolamento Ecológico (peso: 7%)

Mede a proporção de vegetação nativa presente no buffer de 20km que está fora de qualquer UC — remanescentes naturais não protegidos que funcionam como corredores ecológicos, permitindo fluxo gênico, migração de fauna e recolonização após perturbações. Quanto maior essa proporção de vegetação não protegida no entorno, menor o isolamento efetivo da UC. A variável é invertida para que valores altos indiquem maior pressão — ou seja, menor conectividade com remanescentes externos.

Como interpretar: V6 alto significa que a UC está cercada predominantemente por matriz antropizada sem remanescentes naturais externos que funcionem como corredores — situação de isolamento ecológico mesmo que o buffer de 20km não seja inteiramente convertido.

Variável 7 — Vulnerabilidade pela Categoria de Manejo (peso: 3%)

Variável qualitativa que reflete o nível de proteção formal estabelecido pelo Sistema Nacional de Unidades de Conservação — SNUC. Categorias de proteção integral recebem valores menores, indicando menor vulnerabilidade pelo regime de manejo. Categorias de uso sustentável recebem valores maiores, refletindo que permitem usos humanos diretos que reduzem a efetividade da proteção. A Área de Proteção Ambiental recebe o valor mais alto — 0,9 — por ser a categoria mais permissiva do SNUC, permitindo agricultura, pecuária e urbanização dentro de seus limites. Estações Ecológicas e Reservas Biológicas recebem o valor mais baixo — 0,1 — por terem o regime de proteção mais restritivo.

Como interpretar: esta variável não mede pressão externa, mas vulnerabilidade institucional — duas UCs com pressão externa idêntica têm vulnerabilidade real diferente se uma é Parque Nacional e a outra é APA municipal.

Resultados por Esfera Administrativa

A análise das 615 UCs revela uma hierarquia sistemática e consistente de pressão antrópica entre as esferas administrativas.

UCs Federais (229 unidades):

Apresentam a menor pressão média — 0,435 — com 64% das unidades nas faixas Muito Baixo e Baixo e nenhuma UC na faixa Muito Alto.

A UC federal com maior índice de pressão é a Reserva Particular do Patrimônio Natural Chácara Mangueiras — 0,742, seguida pela Reserva Extrativista da Mata Grande — 0,636.

O predomínio de RPPNs e Reservas Extrativistas entre as federais mais pressionadas é significativo — são categorias com menor capacidade de fiscalização direta pelo poder público federal, seja por dependerem de proprietários privados no caso das RPPNs, seja por estarem em regiões de fronteira agrícola ativa no caso das Resex.

UCs Estaduais (247 unidades):

Apresentam pressão média de 0,538 — 24% superior à federal — e o maior desvio padrão do conjunto — 0,170 —, refletindo a heterogeneidade de capacidade de gestão ambiental entre os estados do Cerrado.

A UC estadual com maior índice é a Área de Proteção Ambiental Morro de São Bento — 0,736, seguida pelo Parque Estadual Zé Bolo Flô — 0,715.

A presença de um Parque Estadual entre as mais pressionadas é preocupante — é uma categoria de proteção integral que deveria ter gestão mais rigorosa.

UCs Municipais (139 unidades):

Apresentam a maior pressão média — 0,627, com 86% das unidades nas faixas Moderado e Alto.

A UC municipal com maior índice é o Parque Natural Municipal das Capivaras — 0,794, seguido pela Área de Proteção Ambiental Jupiá — 0,782.

A concentração de APAs entre as municipais mais pressionadas reflete a tendência de municípios criarem UCs nas categorias mais permissivas do SNUC em contextos já urbanizados, onde a proteção efetiva é estruturalmente limitada.

A diferença de 44% entre a pressão média das municipais e das federais não é atribuível apenas à localização em contextos urbanos — reflete também diferenças reais de capacidade institucional, recursos de fiscalização e rigor na aplicação das normas de proteção. UCs municipais combinam localização estruturalmente vulnerável com menor capacidade de resposta institucional, configurando o grupo de maior vulnerabilidade sistêmica no Cerrado.

Tipologia de Pressão

Além do índice numérico, cada UC foi classificada numa tipologia qualitativa que contextualiza a natureza da pressão predominante.

Pressão Urbana Consolidada: caracteriza UCs em contexto urbano ou periurbano — proporção de área urbana superior a 20% no buffer de 10km — onde a pressão é estrutural e preexistente e o desafio de gestão é conter a expansão e manter conectividade ecológica em ambiente hostil.

Pressão Rural Crescente: caracteriza UCs com perda ativa de vegetação nativa no entorno entre 2015 e 2024 superior a 10% — são as unidades sob ameaça dinâmica mais preocupante em termos de tendência futura.

Pressão Rural Estabilizada: caracteriza UCs com entorno degradado mas sem conversão recente significativa — a pressão histórica foi alta mas aparentemente arrefeceu.

Pressão Moderada: caracteriza UCs com entorno relativamente preservado e sem dinâmica recente de conversão acelerada — geralmente as federais mais remotas do Cerrado.

Lacunas de Conservação

Os 44,29% de vegetação nativa remanescente fora de qualquer UC — 887.758,84 km² — representam o dado mais crítico do projeto em termos de política pública. São áreas que ainda preservam o bioma mas estão completamente expostas à conversão agropecuária e urbana sem qualquer proteção formal.

Quando somados à vegetação dentro das UCs — 50,85% da área protegida —, o Cerrado ainda mantém uma cobertura vegetal nativa expressiva, mas distribuída de forma altamente vulnerável: mais da metade dessa vegetação está fora de qualquer proteção.

A identificação espacial dessas lacunas — onde estão, em quais estados, próximas a quais UCs — é uma das contribuições analíticas centrais deste projeto para subsidiar decisões de criação de novas unidades de conservação ou corredores ecológicos.

1. Origem do Projeto

Este projeto nasceu de uma conversa com o assistente de inteligência artificial Claude, da Anthropic, durante uma busca por projetos de geoprocessamento com potencial de aplicação prática no mercado ambiental brasileiro. A partir da constatação de que o Cerrado — o bioma mais ameaçado do Brasil — conta com menos de 10% do território sob proteção formal, surgiu a ideia de desenvolver uma análise que fosse além da simples visualização das Unidades de Conservação existentes, quantificando a pressão antrópica que essas unidades enfrentam e identificando padrões de vulnerabilidade por esfera administrativa. A concepção e estruturação da metodologia foram desenvolvidas em diálogo com o Claude ao longo de várias sessões de trabalho. Todas as decisões analíticas, a implementação técnica, o tratamento dos dados e a validação dos resultados foram realizados pelo autor.

O processo não foi linear — envolveu descobertas, revisões de abordagem e ajustes metodológicos sucessivos, alguns dos quais estão descritos na aba sobre o projeto.

2. Obtenção dos dados

Os dados utilizados no projeto foram obtidos integralmente de fontes públicas e oficiais brasileiras.

Os polígonos das Unidades de Conservação foram obtidos do portal de dados abertos do Ministério do Meio Ambiente e Mudança do Clima.
O limite oficial do bioma Cerrado foi obtido do IBGE.
Os rasters de uso e cobertura do solo dos anos de 2000, 2015 e 2024 foram obtidos do MapBiomas Coleção 10 — o sistema de monitoramento anual da cobertura e uso da terra no Brasil.
A malha estadual utilizada para contextualização geográfica também foi obtida do IBGE.

Uma particularidade relevante foi identificada nessa etapa: o limite vetorial do bioma Cerrado disponibilizado pelo IBGE é ligeiramente maior do que o limite implícito nos rasters do MapBiomas, o que exigiu atenção nos recortes para garantir consistência espacial entre as fontes.

3. Preparação e tratamento dos dados

Após a obtenção dos dados, todos os vetores foram carregados e reprojetados para o EPSG 5880 — SIRGAS 2000 Brasil Policônico — projeção escolhida por cobrir todo o território brasileiro com distorção mínima de área, necessária porque as UCs do Cerrado se distribuem por múltiplas zonas UTM. As UCs foram separadas por esfera administrativa — federal, estadual e municipal — e recortadas pelo limite do bioma Cerrado, garantindo que apenas as porções efetivamente inseridas no bioma entrassem nos cálculos.

UCs parcialmente inseridas no bioma foram mantidas com sua porção recortada, decisão metodológica que preserva a representação completa da realidade de proteção do bioma.

Os buffers de 5km, 10km e 20km foram gerados ao redor de cada UC, com a área interna da própria UC subtraída — garantindo que os cálculos de uso do solo nos buffers reflitam exclusivamente o entorno, não o interior das unidades. As áreas dos buffers foram calculadas usando geometria elipsoidal sobre o datum SIRGAS 2000 para máxima precisão. Todo o processamento tabular e vetorial foi realizado em ambiente Jupyter Notebook com Python, GeoPandas, Pandas e Rasterio.

4. Processamento dos rasters

Os rasters do MapBiomas de 2000, 2015 e 2024 foram recortados no QGIS pelo limite vetorial do Cerrado do IBGE. Em seguida foram reclassificados em seis categorias: cerrado — vegetação nativa —, pastagem, agricultura, área urbana, rios e corpos d'água naturais, e praias e áreas sem vegetação. Essa última categoria foi incluída para garantir a correta representação visual de feições naturais específicas sem distorcer os cálculos de pressão antrópica.

Para cada UC e cada buffer foram extraídas as áreas — em ha e em proporção do buffer total — de cada classe de uso do solo nos anos de 2015 e 2024. As proporções foram a base de cálculo de todas as variáveis do índice, garantindo comparabilidade entre UCs de tamanhos muito diferentes. Para identificar o percentual de cobertura vegetal nativa dentro e fora das UCs em relação ao bioma total, todos os polígonos das UCs foram dissolvidos numa geometria única — eliminando sobreposições — e cruzados com os rasters de uso do solo.

5. Construção do índice de pressão antrópica

O índice foi construído de forma iterativa, com sucessivos ajustes nas variáveis e nos pesos a partir da análise dos resultados intermediários. A versão inicial utilizava apenas variáveis de variação histórica de longo prazo, mas os primeiros resultados revelaram um problema metodológico relevante — UCs já muito degradadas no início do período de análise apareciam com baixa pressão por não terem vegetação remanescente para perder. Esse diagnóstico levou à reformulação das variáveis, incorporando medidas de estado atual do entorno — V1 e V2 — ao lado da medida de dinâmica recente — V3.

A versão final do índice combina sete variáveis com pesos diferenciados: intensidade antrópica ponderada no entorno imediato de 5km com peso de 25%, intensidade antrópica ponderada no entorno próximo de 10km com peso de 20%, taxa de conversão recente entre 2015 e 2024 no buffer de 10km com peso de 20%, isolamento regional no buffer de 20km com peso de 15%, vulnerabilidade intrínseca por tamanho e forma da UC com peso de 10%, isolamento ecológico por conectividade com remanescentes externos com peso de 7%, e vulnerabilidade pela categoria de manejo segundo o SNUC com peso de 3%.

As variáveis de uso do solo utilizam uma ponderação interna que diferencia o tipo de pressão: pastagem com multiplicador 1,0, agricultura com multiplicador 1,5 e área urbana com multiplicador 2,0 — refletindo o grau crescente de irreversibilidade de cada uso. Todas as variáveis foram normalizadas entre 0 e 1 antes da combinação, exceto V7 — categoria de manejo — que já está nessa escala por construção.

6. Desenvolvimento da aplicação

Os mapas interativos foram desenvolvidos com a biblioteca Folium em Python, com sucessivos ciclos de teste e ajuste visual iniciados ainda na fase de análise de dados no Jupyter Notebook. A aplicação web foi construída com o framework Flask, integrando os mapas, os indicadores estatísticos, os gráficos e a interface de navegação por esfera administrativa e categoria de pressão.

Para a definição do layout visual da aplicação, imagens do protótipo funcional foram submetidas aos assistentes de IA ChatGPT e Qwen, cujas sugestões de modernização foram avaliadas e seletivamente incorporadas pelo autor. O resultado visual final foi alcançado após múltiplos ciclos de refinamento, buscando equilíbrio entre clareza analítica, riqueza de informação e experiência de uso.

O sistema integra ferramentas de análise de dados, geoprocessamento e desenvolvimento web para transformar dados demográficos em visualizações interativas.