Paisagens epidemiológicas na Amazônia: explorando a associação entre uso da terra, indicadores ambientais e doenças vetoriais

Resumo:

Este estudo investiga, por meio do uso de indicadores ambientais, uso da terra e técnicas de análise de agrupamento, características de municípios da Amazônia Legal que podem estar associadas à incidência de doenças vetoriais. Identificamos e descrevemos seis agrupamentos de municípios amazônicos com características ambientais e agrárias similares. Com base em um conjunto de indicadores epidemiológicos, exploramos a incidência de doenças tropicais vetoriais negligenciadas ao longo desses agrupamentos. Os resultados evidenciam uma grande heterogeneidade ambiental dos municípios amazônicos, com perfis diferenciados e bem definidos nos agrupamentos obtidos. Além disso, demonstram a relação de atividades econômicas e degradação ambiental com a propagação de doenças. Essa abordagem permite uma visão abrangente dos desafios ambientais e de saúde enfrentados pela região, contribuindo para a elaboração de estratégias de conservação mais eficazes e adaptadas às necessidades específicas de cada perfil de município.

Palavras-chave:
Indicadores Ambientais; Ecossistema Amazônico; Vetores de Doenças

Abstract:

This study investigates the characteristics of municipalities in the Legal Amazon that may be associated with the incidence of vector-borne diseases using environmental and land use indicators and cluster analysis. We identified and described six groups of Amazonian municipalities with similar environmental and agrarian characteristics. Based on a set of epidemiological indicators, we explore the incidence of neglected vector-borne tropical diseases in these groups. Results show a great environmental heterogeneity in the Amazonian municipalities, with well-defined profiles in the obtained groups. Moreover, they show the relation of economic activities and environmental degradation with the spread of diseases. This approach can comprehensively show the environmental and health challenges of regions, contributing to the development of more effective conservation strategies to be adapted to the specific needs of each municipality profile.

Keywords:
Environmental Indicators; Amazonian Ecosystem; Diseases Vectors

Resumen:

A partir de indicadores ambientales, este estudio evalúa técnicas de uso del suelo y análisis de conglomerados, las características de los municipios de la Amazonía Legal que pueden estar asociadas con la incidencia de enfermedades vectoriales. Para ello, identifica y describe seis grupos de municipios amazónicos con características ambientales y agrarias similares. Con base en un conjunto de indicadores epidemiológicos, se evalúa la incidencia de enfermedades tropicales vectoriales desatendidas en estos grupos. Los resultados muestran una gran heterogeneidad ambiental de los municipios amazónicos, con perfiles diferenciados y bien definidos en los grupos evaluados. Además, demuestran la relación de las actividades económicas y la degradación ambiental con la propagación de enfermedades. Este enfoque permite una visión integral de los desafíos ambientales y sanitarios que enfrenta la región y contribuye al desarrollo de estrategias de conservación más efectivas y adaptadas a las necesidades específicas de cada perfil de municipio.

Palabras-clave:
Indicadores Ambientales; Ecosistema Amazónico; Vectores de Enfermedades

Introdução

A Amazônia é reconhecida mundialmente como uma das regiões mais biodiversas e ecologicamente importantes do planeta, desempenhando um papel crucial na regulação do clima regional e global ^1,2. A diversidade cultural, linguística e de assentamentos humanos na região amazônica brasileira é notável, englobando diversas etnias indígenas, comunidades quilombolas e ribeirinhas, além de apresentar uma variedade de núcleos urbanos, que vão desde pequenos núcleos urbanizados, pequenas e médias cidades até centros urbanos industrializados ligados ao Polo Industrial, como Manaus (Amazonas), ou às cadeias de matérias primas (minério e produtos agropecuários) exportadas, como Belém (Pará) ³. No entanto, a crescente pressão humana sobre o bioma tem levantado preocupações significativas sobre a sua sustentabilidade a longo prazo. A região está cada vez mais sujeita a impactos resultantes das formas de ocupação que suprimem a floresta, levando ao aumento da temperatura e eventos de secas extremas como o experimentado em 2023 ^4,5, e maior susceptibilidade a incêndios ^6,7,8. Estudos sugerem que partes da floresta Amazônica poderão sofrer mudanças estruturais significativas no futuro devido às mudanças climáticas e às intensas conversões de uso e cobertura da terra, acelerando as emissões de carbono e impactando a biodiversidade e a saúde humana, o que também afeta valores socioeconômicos e culturais da região ^9,10.

A atual situação ambiental da Região Amazônica é resultado de múltiplas e contínuas ações humanas que afetam a região historicamente e comprometem sua integridade ambiental ^11,12,13. Dentre os principais fatores que afetam a região e levam ao distúrbio e perda florestal, destacam-se a exploração madeireira, o uso do fogo no manejo de áreas com atividades agropecuárias, a grilagem de terras ilegais, a exploração mineral e a expansão da fronteira agropecuária e de projetos de infraestrutura sobre áreas de floresta ^14,15,16. A ocorrência e a intensidade desses fatores ameaçam a integridade do bioma e são modulados por uma complexa interação entre fatores sociais, fundiários e econômicos, associados à ocupação histórica dos municípios amazônicos.

Ao longo das últimas décadas, a ocupação da Amazônia ocorreu de maneira heterogênea tanto no tempo quanto no espaço, influenciada por políticas governamentais, interesses econômicos e pressões socioambientais. Enquanto algumas áreas foram colonizadas de forma planejada, por meio da implementação de projetos de assentamento agropecuários e programas de desenvolvimento regional criados pelo governo, outras foram ocupadas de modo espontâneo e desordenado, por meio de processos de migração das populações provenientes do Sul e do Nordeste, principalmente que vieram em busca de terra e de oportunidades econômicas ^17,18. As políticas de ocupação não deram conta de acomodar nos assentamentos rurais todo o contingente populacional que migrou para essa região na década de 1970, sejam eles do Instituto Nacional de Colonização e Reforma Agrária (INCRA) ou de projetos de colonização privada. Essa diversidade na ocupação da paisagem amazônica resultou em uma variedade de padrões de uso da terra, que incluem desde extensas áreas desmatadas para atividades agropecuárias até as áreas de ocupação tradicional, como as terras indígenas, quilombolas, ribeirinhas, além das de proteção integral e de manejo sustentável com uso extrativista.

Diante desse contexto, observa-se que as transformações ambientais e da paisagem ocorridas na Amazônia estão intimamente ligadas à dinâmica da economia agrária regional que foi se estruturando ao longo dos anos ¹³. A heterogeneidade dos padrões de ocupação da terra, assim como a dinâmica da economia agrária, tem implicações significativas nos padrões de degradação ambiental na região. Além disso, a Região Amazônica já sofre os impactos da mudança climática global, evidenciados pelo aumento acentuado das temperaturas, pelas alterações nos padrões de chuva e pelas secas e incêndios intensos que têm afetado o bioma nos últimos anos ^10,19,20. Tais mudanças resultam na redução da resiliência da floresta, em mudanças na composição de espécies, mortalidade em massa das árvores e perda irreversível de biodiversidade ^21,22.

Por sua vez, os modos de produção associados aos diferentes sistemas agrários e as transformações na paisagem e na biodiversidade por eles promovidas podem ter impactos diretos na saúde humana ^13,23. Por exemplo, as áreas de manejo sustentável para uso extrativista, que mantém extensas áreas de floresta em pé, desempenham um papel crucial na conservação da biodiversidade e na manutenção dos serviços ecossistêmicos. Por outro lado, a conversão de florestas em terras agrícolas pode levar à perda de habitat de espécies animais, que irão procurar alimento e abrigo em áreas de ocupação humana, aumentando a exposição e o risco de doenças zoonóticas que podem ser transmitidas para os seres humanos. De modo geral, o desmatamento e a fragmentação da paisagem podem desencadear o aumento da abundância de reservatórios de doenças e vetores em contato com as comunidades humanas ^24,25,26. Estudos projetam que o risco de doenças transmitidas por vetores na Região Amazônica seja influenciado pelas mudanças climáticas e ambientais ^27,28. Assim, diferentes atividades econômicas podem impactar a paisagem, contribuindo para a emergência e reemergência de doenças ²⁹. Leishmanioses, malária, doença de Chagas e dengue são doenças tropicais negligenciadas prevalentes na Região Amazônica, causadas por agentes infecciosos ou parasitas que atingem predominantemente populações vulneráveis ³⁰. Elas são indicadores de vulnerabilidade social e ambiental ¹³. Alguns estudos apontam para a associação entre a dinâmica de doenças transmissíveis, como a malária, doença de Chagas, dengue e febre amarela, à perda de biodiversidade e às mudanças na paisagem ^28,31. De acordo com esses estudos, a intensificação da pressão sobre esse bioma gera condições socioambientais críticas que podem desencadear o surgimento de novas doenças infecciosas, provenientes de reservatórios silvestres. Portanto, a saúde ambiental e humana na Amazônia brasileira estão intrinsecamente ligadas, e a economia agrária da região desempenha um papel central nessa interconexão, afetando tanto a saúde das populações locais quanto a saúde do ecossistema amazônico como um todo. Entender a heterogeneidade dos fatores que podem levar a grandes transformações ambientais (dentre elas, a perda de floresta) na Região Amazônica é fundamental para a promoção de um desenvolvimento mais sustentável para toda a região.

Diante desse cenário, a compreensão detalhada dos fatores associados às mudanças ambientais que afetam os municípios da Amazônia brasileira é essencial para trazer informações que possam ajudar na orientação de políticas públicas de conservação e para apoiar tomadas de decisão voltadas à gestão sustentável dos recursos naturais e à saúde humana. Nesse contexto, o presente artigo visa responder a seguinte pergunta científica: municípios amazônicos que apresentam características ambientais e de uso da terra semelhantes podem ser associados a um conjunto específico de doenças vetoriais? A fim de elucidar a questão formulada, nós desenvolvemos uma tipologia com categorias de municípios da Amazônia brasileira que apresentam características ambientais e agrárias em comum, tendo como base um conjunto de indicadores ambientais e de uso da terra. Com essa abordagem, torna-se possível identificar, para cada categoria obtida, os principais fatores econômicos envolvidos nas transformações das paisagens, assim como suas condições ambientais, considerando principalmente os processos de perda florestal e suas implicações na transmissão de doenças vetoriais como malária, doença de Chagas, leishmanioses e dengue. Para alcançar esse objetivo, foram utilizados dados e informações disponíveis sobre uma variedade de indicadores ambientais, de uso da terra e epidemiológicos, incluindo desmatamento, degradação florestal, áreas queimadas por incêndios e mineração ³².

Método

Área de estudo

A área de estudo contemplou todos os municípios inteiramente dentro da região da Amazônia Legal brasileira. A Amazônia Legal brasileira é uma região político-administrativa que abrange cerca de 5 milhões de km² e corresponde a aproximadamente 58,9% do território brasileiro. Essa região compreende os estados do Acre, Amapá, Amazonas, Pará, Rondônia, Roraima, Mato Grosso, Tocantins e parte do Maranhão. Em relação aos biomas brasileiros, a Amazônia Legal abriga todo o bioma Amazônico (composto principalmente por florestas tropicais), 20% do bioma Cerrado (composto predominantemente por vegetação de savana) e cerca de 40% do bioma Pantanal (planície de inundação), no Estado de Mato Grosso. A floresta tropical amazônica desempenha papel central na manutenção do ciclo hidrológico regional, contribuindo por meio da reciclagem e transporte de umidade dentro e fora da região ³³. O desmatamento resultante de múltiplas transições de uso e cobertura da terra representa uma ameaça imediata às florestas, sendo o desmatamento para pastagens a transição mais frequente e de alto impacto na Amazônia brasileira ^34,35. Outras transições de alto impacto, que ocorrem com menor frequência ou com menor extensão de área, porém com grande importância, são as conversões para a agricultura ³⁴ e mineração ³⁶. De forma geral, a intensificação ou diminuição da degradação do bioma amazônico respondem às mudanças nas políticas ambientais, como a fiscalização do desmatamento e das atividades associadas ao mercado ilegal de terras, o incentivo para atividades econômicas, a construção de infraestruturas que possibilitam o escoamento da produção agropecuária (como estradas, portos e hidrelétricas) e a oscilação dos preços das commodities. Segundo o Projeto de Monitoramento do Desmatamento na Amazônica Brasileira por Satélite (PRODES), o valor estimado do desmatamento acumulado para o ano de 2023 foi de 9.001km², representando uma redução de 22,37% em relação à taxa de desmatamento consolidada pelo PRODES 2022, que foi de 11.594km² para os nove estados da Amazônia Legal, mas ainda se mantendo em um alto patamar. As grandes e rápidas transformações da paisagem florestal podem levar à desequilíbrios nesse ecossistema e, por consequência, provocar mudanças nos ciclos de transmissão de doenças vetoriais, devido à perda de habitat e à proximidade da população humana aos remanescentes florestais, vetores e reservatórios, levando à uma maior exposição à doenças vetoriais.

Indicadores ambientais e epidemiológicos

Nós utilizamos um conjunto de indicadores ambientais para caracterizar o perfil ambiental dos municípios amazônicos observando as seguintes dimensões: mudanças ambientais, uso e cobertura da terra (LULC − land use and land cover) e anomalias climáticas. Compondo a dimensão de mudança ambiental temos os indicadores de desmatamento, degradação florestal, incêndios, mineração e fragmentação florestal. A dimensão LULC é representada por indicadores de agricultura, pastagens, área urbanizada e vegetação secundária. Indicadores de anomalias de precipitação e temperatura representam a dimensão de anomalias climáticas. Os indicadores foram calculados para todos os municípios da Amazônia Legal e foram construídos a partir de bases de dados abertas e consolidadas. Nós também exploramos a associação dos agrupamentos de municípios com a predominância de doenças tropicais negligenciadas vetoriais. Para isso, utilizamos um conjunto de indicadores epidemiológicos de incidência de leishmaniose cutânea (tegumentar) e visceral, malária, doença de Chagas e dengue para o período de 2015-2019, a partir das bases de dados do DATASUS (Departamento de Informática do Sistema Único de Saúde). A taxa de incidência é uma medida da frequência da doença em populações dinâmicas, e foi estimada dividindo o número de casos pela população no meio do período e multiplicando por cinco (quantidade de anos de observação). Todos os indicadores ambientais e epidemiológicos utilizados neste estudo são descritos detalhadamente em Rorato et al. ³² e estão disponíveis para acesso em: https://zenodo.org/records/7098053. Os indicadores ambientais, de uso da terra e epidemiológicos estão resumidos no Quadro 1.

Quadro 1
Indicadores ambientais, de uso da terra e epidemiológicos.

DIMENSÃO	INDICADOR	DESCRIÇÃO	UNIDADE	PERÍODO	ID
Mudança ambiental	Desmatamento	a área total de desmatamento no município entre 2010 a 2016 é dividida pela sua área florestal original	km2/km2	2010-2016	desmat
Mudança ambiental	Degradação florestal	a área total de degradação florestal no município entre 2007 a 2017 é dividida pela sua área florestal original	km2/km2	2007-2017	degrad
Mudança ambiental	Fogo	proporção da área do município que sofreu pelo menos um incêndio entre 2012 e 2017	km2/km2	2012-2017	fogo
Mudança ambiental	Mineração	proporção da área do município utilizada para mineração (industrial e artesanal) em 2017	km2/km2	2017	mine
Mudança ambiental	Fragmentação	proporção da área do município classificada como núcleo de vegetação natural em 2017	km2/km2	2017	nucleo
Mudança ambiental	Fragmentação	razão entre a borda total (perímetro) das manchas de vegetação natural do município e a raiz quadrada da área total dessas manchas, em 2017	m/sqrt(m2)	2017	borda
LULC	Vegetação secundária	proporção da área do município coberta por vegetação secundária em 2017	km2/km2	2017	veg_sec
LULC	Pastagem	proporção da área do município classificada como pastagem em 2017	km2/km2	2017	pasto
LULC	Agricultura	proporção da área do município classificada como agricultura em 2017	km2/km2	2017	agric
LULC	Área urbana	proporção da área do município classificada como urbana em 2017	km2/km2	2017	urb
Anomalias climáticas	Precipitação	média da área de anomalia de precipitação positiva observada durante as estações secas de 2007-2017, dividida pela área total	m2/m2	2007-2017	precp+
Anomalias climáticas	Precipitação	média da área de anomalia de precipitação negativa observada durante as estações secas de 2007-2017, dividida pela área total	m2/m2	2007-2017	precp-
Anomalias climáticas	Temperatura	média da área de anomalia de temperatura positiva observada durante as estações frias de 2007-2017, dividida pela área total	m2/m2	2007-2017	temp+
Epidemiológico	Doenças transmitidas por vetores	incidência de doença de Chagas, leishmaniose cutânea, leishmaniose visceral, dengue, malária	Taxa de incidência	2015-2019	inc_

• ID: identificador; LULC: uso e cobertura da terra (land use and land cover, acrônimo em inglês).

Análise de agrupamentos

Com o propósito de identificar grupos de municípios na Amazônia que enfrentam mudanças ambientais similares, conduzimos uma análise de agrupamentos utilizando todas as variáveis ambientais acima descritas. Essa técnica não supervisionada visa dividir um conjunto de dados em grupos similares, comumente chamados de clusters, visando maximizar a homogeneidade intra-cluster e a heterogeneidade inter-cluster³⁷. Nós utilizamos o método de agrupamento k-means, uma técnica não hierárquica amplamente empregada devido à sua simplicidade e eficácia, em que o usuário determina a priori o número de clusters, k. Nesse método, são selecionados k objetos (municípios) como centros iniciais, e os demais são alocados com base na similaridade, garantindo a menor distância entre cada objeto e o centro do cluster. O centro de cada cluster, conhecido como centro de gravidade, é calculado como a média dos objetos do cluster. Neste estudo, a distância euclidiana quadrada foi utilizada como medida de similaridade. Esse método permitiu agrupar municípios com características ambientais semelhantes com base nos valores dos indicadores ambientais e agrários.

Antes da análise, aplicamos uma transformação logarítmica nos indicadores para melhorar a normalidade de distribuição das variáveis. Em seguida, padronizamos os indicadores logarítmicos usando o método z-score, visando melhorar a precisão do algoritmo k-means, o qual se baseou na distância euclidiana das observações ^38,39. O número ideal de clusters foi determinado utilizando o método de Elbow, que consiste em plotar a variação explicada em função do número de grupos. A escolha do número de grupos é feita observando o gráfico e escolhendo o número de grupos de acordo com o cotovelo da curva, ou seja, o número de grupos em que a curva plotada começa a se assemelhar a uma linha reta. Todas as análises e visualizações foram conduzidas utilizando os pacotes stats e factoextra no software R (http://www.r-project.org) e o aplicativo Quantum Gis 3.0 (https://qgis.org/en/site/).

Modelagem da distribuição de doenças transmitidas por vetores nos agrupamentos ambientais

Modelos lineares generalizados foram utilizados para identificar associações entre a incidência das doenças de transmissão vetorial e os agrupamentos de municípios gerados pela classificação ambiental. Foram testados modelos com distribuição de Poisson, binomial negativa e binomial negativa com inflação de zeros, no ambiente R, utilizando as bibliotecas MASS e pscl. Todos tiveram como desfecho o total de casos no município (período de 2015 a 2019), e como offset o logaritmo da população estimada em 2017, com a classificação do agrupamento ambiental sendo a única variável explicativa. Para a dengue e as leishmanioses, o modelo binomial negativo teve a melhor qualidade de ajuste, medido pelo critério de informação de Akaike (AIC). Para malária e Chagas, esse não convergiu, e foi necessário recorrer ao modelo inflacionado de zeros.

Resultados e discussão

Classificação ambiental e de uso da terra dos municípios

Com base no método de Elbow, nós encontramos 6 grupos (k = 6) em nossa análise. Cada grupo é composto por um conjunto de municípios com perfil ambiental similar, cuja distribuição espacial é apresentada na Figura 1. A Figura 2 apresenta a distribuição dos grupos de municípios no espaço multivariado, cujos eixos são as três primeiras componentes principais. Em relação a separação dos grupos, é possível observar que os grupos 1, 3 e 6 são bem separados em relação ao eixo 1, enquanto o cluster 6 se destaca em relação ao eixo 2 (Figura 2a). O eixo 2 também separa o grupo 5 de todos os demais. Já os grupos 2 e 4 são melhor evidenciados em relação ao eixo 3 (Figura 2b). Nossos resultados evidenciam uma grande heterogeneidade ambiental dos municípios Amazônicos (Figuras 3 e 4), com perfis bem definidos em cada agrupamento de maneira geral. As características ambientais de cada grupo são sumarizadas nas Figuras 3 e 4 e são discutidas a seguir.

Figura 1
Classificação ambiental dos municípios da Amazônia Legal brasileira em 2017, obtida pela análise de agrupamentos (k-means).

Figura 2
Distribuição dos municípios da Amazônia Legal brasileira no espaço multivariado gerado pelas variáveis ambientais do estudo.

Figura 3
Gráficos de radar indicando a distribuição das variáveis ambientais em cada agrupamento de municípios da Amazônia Legal brasileira.

Figura 4
Boxplots dos indicadores ambientais ao longo dos agrupamentos de municípios da Amazônia Legal brasileira.

O agrupamento 1 é composto por 169 municípios, sendo o segundo maior agrupamento. Os municípios desse agrupamento estão em sua maioria localizados nos Estados do Pará, Mato Grosso e Maranhão, assim como alguns municípios isolados no interior dos estados do Amapá, Roraima, Amazonas, Rondônia, Tocantins e Acre (Figura 1). Esse agrupamento se destaca por apresentar municípios com as maiores médias de desmatamento, degradação florestal e efeito de borda florestal entre todos os grupos. Fazem parte desse grupo municípios como Santarém (Oeste do Pará), Novo Progresso e Jacareacanga (Sudoeste do Pará) e Colniza (Noroeste do Mato Grosso), que apresentam grandes áreas florestais, muitas vezes concentradas em unidades de conservação ou Terras Indígenas. Entretanto, parte desses municípios, como Novo Progresso e Colniza, também figuram entre os 10 municípios que mais desmataram no período considerado até 2023 ⁴⁰. O grupo também apresenta a segunda maior média em áreas queimadas (Figuras 3 e 4). No caso do desmatamento, práticas como agricultura mecanizada em larga escala e pecuária extensiva são impulsionadas pelo mercado de terras e pela conversão de áreas florestais para atender à demanda por commodities.

O agrupamento 2 é formado por 201 municípios, sendo o maior agrupamento de todos. Os municípios desse agrupamento estão concentrados principalmente nos estados do Maranhão, Tocantins, Pará, Acre e Rondônia (Figura 1). São mais numerosos no Estado do Pará, concentrando-se na região nordeste (Zona da Bragantina) e fazendo vizinhança com os municípios do Estado do Maranhão, que também estão incluídos nesse grupo. São municípios em que grande parte dos sistemas agrícolas são realizados em sistemas de pousio. Nesses sistemas, a fertilidade do solo é recuperada por meio do pousio, no qual após a colheita há o abandono da terra e crescimento da vegetação secundária, em intervalos de tempos que podem variar de região para região ^5,41. Fazem parte desse agrupamento municípios ribeirinhos dos estados do Amazonas e do Pará, em que o sistema de pousio também é prática comum. Os municípios desse agrupamento possuem as maiores médias de área de vegetação secundária, urbanização e anomalia de temperatura positiva (Figuras 3 e 4). São municípios que apresentaram pouco desmatamento no período analisado, seja porque apresentam grandes áreas de floresta primária ou porque, ao contrário, apresentam pequenos remanescentes de floresta primária, mas uma proporção grande de vegetação secundária em sua substituição.

O agrupamento 3 é formado por 111 municípios localizados principalmente nos estados do Amazonas, Roraima, Acre, Pará e Amapá (Figura 1). Os municípios desse grupo se localizam na região ocidental da Amazônia, com maior proporção de área de floresta, destacando-se por apresentarem as maiores médias da porcentagem de área de núcleo florestal em comparação aos demais grupos. Em contrapartida, possuem os menores valores para fogo, pastagens, agricultura, urbanização e anomalia de precipitação negativa (Figuras 3 e 4), indicando que são municípios que apresentam menores distúrbios ambientais e floresta mais íntegra.

O agrupamento 4 é composto por 90 municípios localizados principalmente em áreas de uso agropecuário consolidado, concentrando-se principalmente nos estados do Mato Grosso e Rondônia (Figura 1). Esses municípios são caracterizados ambientalmente pelas maiores médias de áreas sob efeito de borda florestal (empatando com o grupo 1), agricultura e anomalias de precipitação positiva. No Estado do Mato Grosso, os municípios indicados podem apresentar grandes proporções de áreas de fitofisionomia não florestal ou de floresta estacional. Além disso, esse agrupamento é o que apresenta a segunda maior área de pastagens, degradação florestal e de anomalias de precipitação negativa (Figuras 3 e 4).

O agrupamento 5 é composto por 43 municípios, localizados principalmente nos estados do Pará, Rondônia, Mato Grosso e Amapá (Figura 1). Esse grupo se destaca por apresentar municípios com maior a média de porcentagem de área de mineração e a segunda maior média de porcentagem de desmatamento em relação aos demais grupos. Destacam-se também a predominância de áreas sob efeito de borda e de degradação florestal (Figuras 3 e 4). Nesse agrupamento, encontram-se municípios como Itaituba e Jacareacanga (oeste do Pará) e São Felix do Xingu e Ourilândia do Norte (leste do Pará), onde a atividade de mineração ilegal tem se intensificado muito nos últimos anos, inclusive no interior de Terras Indígenas como a Kayapó e Munduruku, com significativos impactos sobre as populações locais ⁴². Esses municípios apresentam grandes áreas florestais devido à presença de Unidades de Conservação e Terras Indígenas, e nas áreas desmatadas, a principal cobertura é a pastagem, ainda que o agrupamento não apresente a maior média de áreas de pastagem dos agrupamentos.

O agrupamento 6 é o terceiro maior agrupamento, composto por 159 municípios localizados principalmente nos estados do Tocantins, sul do Maranhão e do Mato Grosso (Figura 1). São municípios pequenos, que se encontram em região de transição com o bioma cerrado e que podem apresentar grandes manchas de vegetação com fitofisionomias do cerrado. Esse agrupamento se destaca por apresentar as maiores médias de anomalia de precipitação negativa, de áreas de pastagens e área queimada, e também por ter a segunda maior média de área de agricultura anual, vegetação secundária e anomalias de temperatura positiva (Figuras 3 e 4).

Associação das doenças de transmissão vetorial com a classificação ambiental

Nessa etapa, nós exploramos a suposição de que a classificação ambiental e de uso da terra produzida pelo estudo se relaciona com a distribuição da incidência de algumas doenças vetoriais na região. Para as cinco doenças estudadas, a inclusão da variável “grupo” resultou em melhor ajuste, medido por análise de Deviance. A Figura 5 e a Tabela 1 apresentam os odds ratios (OR) de ocorrência dos agravos por grupo ambiental, calculados tendo o grupo “3” como referência. Esse grupo é o conjunto de municípios com maior preservação do bioma.

Figura 5
Associação entre a incidência de doenças de transmissão vetorial e os agrupamentos gerados pela classificação ambiental dos municípios da Amazônia Legal brasileira.

Tabela 1
Resultado do ajuste dos modelos lineares generalizados.

Doença/Grupo	OR	IC95%
Dengue
1	1,50	1,12-2,01
2	0,69	0,52-0,91
3	1,00	1,00-1,00
4	1,72	1,23-2,42
5	1,29	0,85-2,02
6	1,50	1,12-2,02
Leishmaniose tegumentar
1	0,86	0,66-1,13
2	0,35	0,27-0,45
3	1,00	1,00-1,00
4	0,77	0,57-1,06
5	0,75	0,51-1,12
6	0,36	0,27-0,47
Leishmaniose visceral
1	3,08	1,88-4,98
2	3,82	2,37-6,08
3	1,00	1,00-1,00
4	0,98	0,55-1,79
5	2,44	1,25-5,03
6	3,19	1,93-5,21
Malária
1	0,07	0,04-0,13
2	0,02	0,01-0,03
3	1,00	1,00-1,00
4	0,02	0,01-0,05
5	0,20	0,09-0,42
6	0,00	0,00-0,01
Doença de Chagas
1	0,33	0,11-1,04
2	0,18	0,08-0,39
3	1,00	1,00-1,00
4	0,01	0,00-0,07
5	0,27	0,09-0,85
6	0	NO

• IC95%: intervalo de 95% de confiança; NO: não observado; OR: odds ratio.

A dengue e a leishmaniose tegumentar apresentaram alta incidência em todos os agrupamentos ambientais (Figura 6). Contudo, a dengue tem as maiores incidências nos grupos 1, 4 e 6 (Figura 5a), que se concentram a Sul e Sudeste da região amazônica, sendo áreas de grande transformação ambiental e que compartilham características como a alta porcentagem de área queimada no período, fragmentação das florestas, forte presença da agricultura, pastagem, maior proporção de vegetação secundária e maior ocorrência de anomalias de precipitação negativa dentre os grupos. Por outro lado, a dengue tem a menor incidência no grupo 2, seguida do grupo 3. Esses grupos compartilham menor área de queimadas, menor degradação e fragmentação florestal, além de menor extensão de área agrícola.

Figura 6
Boxplots dos indicadores epidemiológicos ao longo dos agrupamentos de municípios da Amazônia Legal brasileira.

A leishmaniose tegumentar teve incidência estatisticamente similar nos grupos 1, 3, 4 e 5, que representam 53% de todos os municípios. São municípios que apresentam grande proporção de remanescentes de florestas, com variados graus de desmatamento e fragmentação florestal. Apenas os grupos 2 e 6 apresentaram menor chance dessa doença (Figura 5b). Esses dois grupos compartilham algumas características: menores áreas de núcleos florestais, maiores áreas com vegetação secundária e maior ocorrência de anomalia de temperatura positiva. Historicamente, a leishmaniose cutânea é associada a municípios amazônicos com grande proporção de área florestal, especialmente em comunidades extrativistas. Porém, o avanço da ocupação humana e da urbanização na região levou a leishmaniose cutânea para outras paisagens ⁴³, o que explica a sua ampla distribuição pela região, mas ainda assim apresenta uma relação com a presença da floresta. Novos focos dessa doença são associados a áreas com desmatamento recente para implantação de agricultura ⁴⁴, o que é coerente com os resultados encontrados.

A leishmaniose visceral teve maior incidência nos municípios dos grupos 1, 2, 5 e 6, num padrão bem diferente da leishmaniose tegumentar (Figura 5c). Já os grupos 3 e 4 foram os que apresentaram as menores chances para essa doença. Interessante notar que esses dois grupos (3 e 4) são muito distintos em todas as variáveis (Figura 4), sugerindo que os determinantes associados à proteção a essa doença são distintos em cada paisagem. No grupo 3, tem-se a maior extensão de floresta contínua e menor fragmentação florestal; no grupo 4, tem-se área de transição com o cerrado e são áreas com poucos remanescentes florestais e que apresentam grandes áreas de produção de grãos, principalmente no sul dos estados do Mato Grosso e de Rondônia. A leishmaniose visceral é uma doença de histórico predominantemente rural, que nas últimas décadas apresentou expansão para áreas urbanas e para várias partes do país, incluindo a região oriental da Amazônia Legal brasileira ⁴⁵.

Ocorrência de malária e Chagas apresentou forte associação com o grupo 3, caracterizado pela predominância de floresta contínua (Figuras 5d e 5e). O grupo 5, caracterizado pela mineração, também apresentou risco para ambas, embora bem mais baixo que o grupo 3. A relação da malária com o grupo 3 é explicada pela presença de áreas extensas de floresta contínua, uma vez que o vetor é adaptado ao ambiente humano próximo da mata em áreas de inundação ⁴⁶. Quanto ao grupo 5, além das áreas de floresta contínua, a incidência de malária também pode estar associada às atividades de mineração e ao desmatamento, que são predominantes nesse grupo, configurando os municípios de fronteira de ocupação, como Itaituba e São Félix do Xingu ⁴⁷.

Municípios do grupo 1 não diferiram estatisticamente do grupo 3 quanto à incidência de Chagas. Essas áreas do grupo 1 de maior incidência de Chagas estão concentradas no Pará. Interessante notar que os grupos 3 e 1 apontam para paisagens distintas dessa doença, a primeira associada à floresta, com presença de palmeiras que são habitats do triatomíneo, e onde está o reservatório ambiental predominante do parasita; o segundo, com a paisagem de degradação ambiental associada à redução das fontes de alimentação silvestres do triatomíneo e à adaptação ao meio urbano ⁴⁸.

Nossos resultados estão alinhados com estudos que sugerem a influência de atividades econômicas na degradação dos ecossistemas locais e na propagação de doenças ^13,34,49. Isso evidencia que as ações de prevenção e controle das doenças vetoriais na Amazônia precisam estar alinhadas com o modo de viver das pessoas da região, observando sua relação com o ambiente, mediada pela economia e pelos modos de produzir ³⁰. Ao avaliar fatores ambientais e agrários e seu impacto na prevalência das doenças, as intervenções podem ser mais direcionadas e eficazes.

A proximidade geográfica observada nos agrupamentos reflete dinâmicas comuns de ocupação e uso da terra entre municípios vizinhos, indicando a influência de processos macrorregionais. Fatores como desmatamento, degradação florestal, fogo e mudanças no uso da terra podem apresentar forte dependência espacial, reforçando a conexão entre áreas próximas. Essa contiguidade territorial é relevante, não apenas para compreender os padrões históricos de ocupação, mas também para dar suporte ao planejamento territorial e, oportunamente, à uma estratégia de gestão ambiental mais integrada. Ao considerar estratégias conjuntas entre municípios, é possível enfrentar desafios com ações compartilhadas entre eles, o que acreditamos ser mais eficaz, incluindo ações de adaptação e mitigação às mudanças climáticas em escalas locais e regionais.

Cabe ressaltar que existem limitações nas análises apresentadas, já que os indicadores ambientais não incluem variáveis importantes como a presença de cães e outros hospedeiros relevantes para alguns desses agravos, e também não estão incluídos dados mais específicos de áreas urbanas. Além disso, reconhecemos que o uso dos municípios como unidades territoriais de análise constitui uma escolha metodológica que limita a análise dos resultados à escala municipal. Ao escolher os municípios da Amazônia como unidade de análise devido a questões operacionais relacionadas com as bases de dados, observa-se que a heterogeneidade intramunicipal não é capturada nos resultados, como era de se esperar. Em relação a heterogeneidade de tamanho entre os municípios e o possível viés relacionado ao tamanho dos mesmos, destacamos que todos os indicadores foram construídos levando em conta a proporção de área em relação ao município ou à área de floresta do município (Quadro 1), objetivando minimizar um possível efeito relativo à área. Mesmo com essas limitações, os resultados apontam para o potencial de produção de classificações ambientais de interesse para a análise epidemiológica da paisagem.

Conclusão

As associações observadas neste estudo entre indicadores ambientais e de uso da terra com a incidência de doenças realça a importância de se avaliar a saúde dos ecossistemas, resultante de processos de mudanças na paisagem, que afetam a provisão dos serviços ecossistêmicos e consequentemente o bem-estar e a saúde humana. Em relação à pergunta que norteou este trabalho, sobre a possibilidade dos municípios amazônicos que apresentam características ambientais e de uso da terra semelhantes poderem ser associados a um conjunto específico de doenças vetoriais, os indicadores explorados apontam que sim. Apontam também para a necessidade de aprofundamento das análises e para a ampliação do conjunto de indicadores, acrescentando outras dimensões como a social e a econômica. Em adição, destacamos a necessidade de colaboração interdisciplinar entre especialistas em saúde pública, ecologistas e formuladores de políticas para enfrentar os complexos desafios da Região Amazônica. Ao abordar os desafios interligados da epidemiologia, da biodiversidade e da ocupação da região, as intervenções futuras podem contribuir para melhores resultados de saúde, sustentabilidade ambiental e bem-estar da população.

Destacamos ainda que a conservação da floresta amazônica é de grande relevância devido à sua influência no clima local, regional e global, na manutenção da biodiversidade, provisão de alimentos e água e no bem-estar das comunidades locais e indígenas. Assim, no planejamento e no desenvolvimento de estratégias de intervenções em saúde pública, o estado de conservação, os modos de produção indicados pelos usos da terra e os processos de restauração dos ecossistemas devem ser considerados como componentes integrantes das ações de prevenção a doenças.

Agradecimentos

Este trabalho foi financiado pelo Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq, Centro de Síntese em Biodiversidade e Serviços Ecossistêmicos − SinBiose, processo: 442357/2019-2).

Referências bibliográficas

• ¹ Strand J, Soares-Filho B, Costa MH, Oliveira U, Ribeiro SC, Pires GF, et al. Spatially explicit valuation of the Brazilian Amazon Forest's ecosystem services. Nat Sustain 2018; 1:657-64.
• ² Nobre C, Encalada A, Anderson E, Neves EG. Science panel for the Amazon. Amazon Assessment Report 2021: executive summary. https://repositorio.usp.br/item/003095068 (accessed on 11/Apr/2024).
  » https://repositorio.usp.br/item/003095068
• ³ Becker BK. A urbe amazônida: a floresta e a cidade. Rio de Janeiro: Garamond Universitária; 2013.
• ⁴ Costa F, Marengo J. Statement on the 2023 Amazon drought and its unforeseen consequences. https://www.theamazonwewant.org/wp-content/uploads/2023/12/231207-AMAZON-DROUGHT-STATEMENT_ENGLISH.pdf (accessed on 11/Apr/2024).
  » https://www.theamazonwewant.org/wp-content/uploads/2023/12/231207-AMAZON-DROUGHT-STATEMENT_ENGLISH.pdf
• ⁵ Silva SS, Brown F, Sampaio AO, Silva ALC, Santos NCRS, Lima AC, et al. Amazon climate extremes: increasing droughts and floods in Brazil's state of Acre. Perspect Ecol Conserv 2023; 21:311-7.
• ⁶ Oliveira G, Mataveli G, Stark SC, Jones MW, Carmeta R, Brunsell NA, et al. Increasing wildfires threaten progress on halting deforestation in Brazilian Amazonia. Nat Ecol Evol 2023; 7:1945-6.
• ⁷ Barbosa MLF, Delgado RC, Andrade CF, Teodoro PE, Silva Junior CA, Wanderley HS, et al. Recent trends in the fire dynamics in Brazilian Legal Amazon: interaction between the ENSO phenomenon, climate and land use. Environ Dev 2021; 39:100648.
• ⁸ Dutra DJ, Anderson LO, Fearnside PM, Graça PMLA, Yanai AM, Dalagnol R, et al. Fire dynamics in an emerging deforestation frontier in southwestern Amazonia, Brazil. Fire 2023; 6:2.
• ⁹ Burton C, Kelley DI, Jones CD, Betts RA, Cardoso M, Anderson L. South American fires and their impacts on ecosystems increase with continued emissions. Clim Resil Sustain 2022; 1:e8.
• ¹⁰ Flores BM, Montoya E, Sakschewski B, Nascimento N, Staaç A, Betts RA, et al. Critical transitions in the Amazon forest system. Nature 2024; 626:555-64.
• ¹¹ Vieira ICG, Toledo PM, Silva JMC, Higuchi H. Deforestation and threats to the biodiversity of Amazonia. Braz J Biol 2008; 68:949-56.
• ¹² Nepstad DC, Stickler CM, Soares-Filho B, Merry F. Interactions among Amazon land use, forests and climate: prospects for a near-term forest tipping point. Philos Trans R Soc B Biol Sci 2008; 363:1737-46.
• ¹³ Codeço CT, Dal'Asta AP, Rorato AC, Lana RM, Neves TC, Andreazzi CS, et al. Epidemiology, biodiversity, and technological trajectories in the Brazilian Amazon: from malaria to COVID-19. Front Public Health 2021; 9:647754.
• ¹⁴ Pfaff ASP. What drives deforestation in the Brazilian Amazon?: evidence from satellite and socioeconomic data. J Environ Econ Manag 1999; 37:26-43.
• ¹⁵ Aguiar APD, Câmara G, Escada MIS. Spatial statistical analysis of land-use determinants in the Brazilian Amazonia: exploring intra-regional heterogeneity. Ecol Model 2007; 209:169-88.
• ¹⁶ Santos AM, Silva CFA, Almeida Junior PM, Rudke AP, Melo SN. Deforestation drivers in the Brazilian Amazon: assessing new spatial predictors. J Environ Manage 2021; 294:113020.
• ¹⁷ Machado L. A fronteira agrícola na amazônia brasileira. Revista Brasileira de Geografia 1992; 54:27-55.
• ¹⁸ Becker BK. Geopolítica da Amazônia. Estud Av 2005; 19:71-86.
• ¹⁹ Marengo JA, Jimenez JC, Espinoza JC, Cunha AP, Aragão LEO. Increased climate pressure on the agricultural frontier in the Eastern Amazonia - Cerrado transition zone. Sci Rep 2022; 12:457.
• ²⁰ Lapola DM, Pinho P, Barlow J, Aragão LE, Berenguer E, Carmenta R, et al. The drivers and impacts of Amazon forest degradation. Science 2023; 379:eabp8622.
• ²¹ Esquivel-Muelbert A, Baker TR, Dexter KG, Lewis SL, Brienen RJ, Feldpausch TR, et al. Compositional response of Amazon forests to climate change. Global Chang Biol 2019; 25:39-56.
• ²² Boulton CA, Lenton TM, Boers N. Pronounced loss of Amazon rainforest resilience since the early 2000s. Nat Clim Change 2022; 12:271-8.
• ²³ Rohr JR, Civitello DJ, Halliday FW, Hudson PJ, Lafferty KD, Wood CL, et al. Towards common ground in the biodiversity - disease debate. Nat Ecol Evol 2020; 4:24-33.
• ²⁴ Prist PR, Andrea PSD, Metzger JP. Landscape, climate and hantavirus cardiopulmonary syndrome outbreaks. EcoHealth 2017; 14:614-29.
• ²⁵ Bloomfield LSP, McIntosh TL, Lambin EF. Habitat fragmentation, livelihood behaviors, and contact between people and nonhuman primates in Africa. Landsc Ecol 2020; 35:985-1000.
• ²⁶ Wilkinson DA, Marshall JC, French NP, Hayman DTS. Habitat fragmentation, biodiversity loss and the risk of novel infectious disease emergence. J R Soc Interface 2018; 15:20180403.
• ²⁷ Carvalho BM, Perez LP, Oliveira BFA, Jacobson LS, Horta MA, Sobral A, et al. Vector-borne diseases in Brazil: climate change and future warming scenarios. Sustainability in Debate 2020; 11:361-404.
• ²⁸ Winck GR, Raimundo RLG, Fernandes-Ferreira H, Bueno MG, D'Andrea PS, Rocha FL, et al. Socioecological vulnerability and the risk of zoonotic disease emergence in Brazil. Sci Adv 2022; 8:eabo5774.
• ²⁹ Lowe R, Lee S, Lana RM, Codeço CT, Castro MC, Pascual M. Emerging arboviruses in the urbanized Amazon rainforest. BMJ 2020; 371:m4385.
• ³⁰ Codeço CT, Rorato AC, Dal'Asta AP, Escada MIS, Lana RM, Barbosa M,et al. O futuro sustentável da Amazônia: integrando saúde humana, economia e meio ambiente. https://www.gov.br/cnpq/pt-br/acesso-a-informacao/acoes-e-programas/programas/sinbiose-1/pdfs/brief_trajetorias_port_tela.pdf (accessed on 11/Apr/2024)
  » https://www.gov.br/cnpq/pt-br/acesso-a-informacao/acoes-e-programas/programas/sinbiose-1/pdfs/brief_trajetorias_port_tela.pdf
• ³¹ Keesing F, Ostfeld RS. Impacts of biodiversity and biodiversity loss on zoonotic diseases. Proc Natl Acad Sci U S A 2021; 118:e2023540118.
• ³² Rorato AC, Dal'Asta AP, Lana RM, Santos RBN, Escada MIS, Vogt CM, et al. Trajetorias: a dataset of environmental, epidemiological, and economic indicators for the Brazilian Amazon. Sci Data 2023; 10:65.
• ³³ Marengo JA, Souza CMJ, Thonicke K, Burton C, Halladay K, Betts RA, et al. Changes in climate and land use over the Amazon Region: current and future variability and trends. Front Earth Sci 2018; 6:228.
• ³⁴ Nunes CA, Berenguer E, França F, Ferreira J, Lees AC, Louzada J, et al. Linking land-use and land-cover transitions to their ecological impact in the Amazon. Proc Natl Acad Sci 2022; 119:e2202310119.
• ³⁵ Almeida CA, Coutinho AC, Esquerdo JCDM, Adami M, Venturieri A, Diniz CG, et al. High spatial resolution land use and land cover mapping of the Brazilian Legal Amazon in 2008 using Landsat-5/TM and MODIS data. Acta Amaz 2016; 46:291-302.
• ³⁶ Sonter LJ, Herrera D, Barrett DJ, Galford GL, Moran CJ, Soares-Filho BS. Mining drives extensive deforestation in the Brazilian Amazon. Nat Commun 2017; 8:1013.
• ³⁷ Fernandez P, Mourato S, Moreira M, Pereira L. A new approach for computing a flood vulnerability index using cluster analysis. Physics and Chemistry of the Earth 2016; 94:47-55.
• ³⁸ Rorato AC, Picoli MCA, Verstegen JA, Camara G, Silva Bezerra FG, Escada MIS. Environmental threats over Amazonian Indigenous lands. Land 2021; 10:267.
• ³⁹ Marôco J. Análise estatística com utilização do SPSS. Lisboa: Edições Sílabo; 2007.
• ⁴⁰ Sobre TerraBrasillis. https://terrabrasilis.dpi.inpe.br/app/map/deforestation?hl=pt-br (accessed on 11/Apr/2024).
  » https://terrabrasilis.dpi.inpe.br/app/map/deforestation?hl=pt-br
• ⁴¹ Do Vale I, Miranda IS, Mitja D, Grimaldi M, Nelson BW, Desjardins T, et al. Tree regeneration under different land-use mosaics in the Brazilian Amazon's "Arc of Deforestation". Environ Manage 2015; 56:342-54.
• ⁴² Rorato AC, Camara G, Escada MIS, Picoli MCA, Moreira T, Verstegen JA. Brazilian amazon indigenous peoples threatened by mining bill. Environ Res Lett 2020; 15:1040a3.
• ⁴³ Lainson R, Shaw JJ, Ryan L, Ribeiro RS, Silveira FT. Leishmaniasis in Brazil. XXI. Visceral leishmaniasis in the Amazon Region and further observations on the role of Lutzomyia longipalpis (Lutz & Neiva, 1912) as the vector. Trans R Soc Trop Med Hyg 1985; 79:223-6.
• ⁴⁴ Martins LM, Rebêlo JMM, Santos MCFV, Costa JML, Silva AR, Ferreira LA. Ecoepidemiologia da leishmaniose tegumentar no Município de Buriticupu, Amazônia do Maranhão, Brasil, 1996 a 1998. Cad Saúde Pública 2004; 20:735-43.
• ⁴⁵ Conti RV, Moura Lane VF, Montebello L, Pinto Junior VL. Visceral leishmaniasis epidemiologic evolution in timeframes, based on demographic changes and scientific achievements in Brazil. J Vector Borne Dis 2016; 53:99-104.
• ⁴⁶ Reis IC, Honório NA, Barros FSM, Barcellos C, Kitron U, Camara DCP, et al. Epidemic and endemic malaria transmission related to fish farming ponds in the Amazon frontier. PLoS One 2015; 10:e0137521.
• ⁴⁷ Caldas RJC, Nogueira LMV, Rodrigues ILA, Andrade EGR, Costa CML, Trindade LNM. Incidence of malaria among indigenous people associated with the presence of artisanal mining. Rev Gaúcha Enferm 2023; 44:e20220098.
• ⁴⁸ Coura JR, Junqueira AC. Risks of endemicity, morbidity and perspectives regarding the control of Chagas disease in the Amazon Region. Mem Inst Oswaldo Cruz 2012; 107:145-54.
• ⁴⁹ Murta FLG, Marques LLG, Santos APC, Batista TSB, Mendes MO, Silva ED, et al. Perceptions about malaria among Brazilian gold miners in an Amazonian border area: perspectives for malaria elimination strategies. Malar J 2021; 20:286.

Datas de Publicação

• Publicação nesta coleção
  09 Jun 2025
• Data do Fascículo
  2025

Histórico

• Recebido
  12 Abr 2024
• Revisado
  12 Maio 2024
• Aceito
  14 Jan 2025