SOCIEDADE BRASILEIRA PARA
O ESTUDO DE QUIRÓPTEROS
Brazilian Bat Research Society
Modelos de distribuição potencial de espécies de morcegos nos biomas brasileiros
A distribuição geográfica de uma espécie é a manifestação de complexas interações entre as características intrínsecas de um organismo (e.g. fisiologia, história de vida, aspectos demográficos e capacidade de dispersão) com o ambiente, que por sua vez varia no espaço e tempo, limitando a distribuição e abundância das espécies (Brown et al. 1996). Portanto, a distribuição geográfica e sua variação ao longo do tempo são características ecológicas e evolutivas fundamentais de uma espécie, e delimitá-las no período presente é uma das etapas fundamentais para o seu completo conhecimento (Brown et al. 1996).
No entanto, a distribuição geográfica é muitas vezes desconhecida para um número imenso de taxa, principalmente no caso de espécies menos conhecidas e/ou crípticas, ou ainda que ocupam áreas de difícil acesso. Nos últimos 20 anos, a partir da crescente disponibilidade de dados ambientais georreferenciados em boa qualidade e resolução espacial, o fácil acesso a bancos de dados públicos da biodiversidade, e o aumento de poder e velocidade dos computadores, bem como o desenvolvimento de diversos algoritmos, a modelagem de distribuição de espécies ganhou imenso destaque em diversas áreas da biologia, e vem sendo amplamente utilizada para prever a distribuição geográfica das espécies (Elith e Leathwick 2009b; Araújo et al. 2019). O uso destas técnicas é particularmente valioso em países ou áreas mal amostradas, extensas e biologicamente diversas, com recursos limitados (Cooper-Bohannon et al., 2016; Hernandez et al., 2008). Este é o caso do Brasil, com uma dimensão continental, com alta riqueza de espécies, mas com extensas áreas nunca amostradas (Bernard et al. 2011).
A modelagem de distribuição de espécies inclui diversas técnicas, que de forma genérica correlacionam pontos de presença em localidades conhecidas das espécies com informações ambientais e espaciais dessas localidades (Elith e Leathwick 2009a, 2009b). Comumente esses dados de distribuição são encontrados na forma de pontos de ocorrência e/ou "Extensão de Ocorrência – EOO" (IUCN, 2014). Os pontos de ocorrência geralmente subestimam a distribuição das espécies (os chamados erros de omissão), pois as localidades estudadas são geralmente aquelas facilmente acessadas pelos pesquisadores (Brito et al., 2009; Fernández & Nakamura, 2015; Rondinini et al., 2006).
Por outro lado, os mapas de EOO por serem polígonos com o intervalo plausível da espécie, muitas vezes construídos por pontos de interpolação ou desenhando um polígono que inclui os pontos mais externos (técnica de polígono convexo mínimo – MCP), geralmente falham ao representar a distribuição das espécies incluindo áreas inadequadas (i.e., áreas onde a espécie está ausente – Elith & Leathwick, 2009; Rodrigues, 2011). As distribuições estimadas usando esta metodologia geralmente resultam em altos níveis de erros de falsa presença (os chamados erros de comissão) e, assim, as espécies são assumidas como protegidas em locais onde na verdade não ocorrem (Fourcade, 2016; Gaston & Rodrigues, 2003; Rondinini et al., 2005).
Estimar a distribuição de uma espécie da maneira mais acurada possível ou determinar como as espécies irão responder às mudanças climáticas são questões fundamentais no manejo e conservação da biodiversidade. Nesse contexto, avanços têm sido feitos na habilidade de coleta, estocagem e disponibilização de dados temporal e espacialmente referenciados em plataformas públicas. Das 181 espécies de morcegos que ocorrem no Brasil, modelos de distribuição de espécies já haviam sido gerados para 135 delas (Delgado-Jaramillo et al. 2020). Estes modelos foram posteriormente atualizados com mais de 20 mil registros de ocorrência obtidos em todo Brasil (Aguiar et al. 2021). Com o intuito de facilitar a divulgação de dados para além da comunidade acadêmica, a SBEQ, a partir da colaboração com diferentes pesquisadores e instituições, optou por atualizar os modelos de distribuição para as espécies brasileiras de morcegos e disponibilizá-los gratuitamente de forma unificada em sua plataforma online. Aqui você poderá visualizar os modelos gerados e baixar os dados usados na modelagem.
Entretanto, é extremamente importante considerar que modelos de distribuição de espécies são passíveis de erros e de melhorias. Estes modelos podem e devem passar por constante atualização, agregando novos registros e corrigindo registros equivocados. Por isso, entre em contato com a SBEQ caso você tenha algum registro que não tenha sido considerado nesta modelagem, ou caso faça um novo registro que posa corrigir/aprimorar/refinar os modelos aqui disponibilizados. Sua colaboração será muito bem-vinda.
Metodologia
A partir dos estudos de Delgado et al. (2020) e Aguiar et al. (2020), 27.318 registros de morcegos de 172 espécies reconhecidas para o território brasileiro foram obtidos. Neste conjunto inicial de registros, filtros foram aplicados de modo a manter nas etapas subsequentes de análise apenas um registro de ocorrência por quadrante de ~10 x 10 km. Um conjunto de variáveis ambientais também com resolução espacial de ~10 x 10 km e com extensão de todo território brasileiro (porção continental) foi criado a partir de um conjunto inicial de 12 variáveis não-correlacionadas (Tabela 1). O algoritmo Maxent (Phillips et al. 2006) implementado no pacote maxnet (Phillips et al. 2017) foi utilizado para gerar modelos de distribuição para as espécies de quirópteros para cada bioma brasileiro.
É importante destacar que preferiu-se restringir a área de estudo dos modelos por bioma, uma vez que a grande extensão do território brasileiro impõe uma série de dificuldades que afetam negativamente o poder de predição e acurácia dos modelos (Elith et al. 2010, Merow et al. 2013). Os modelos foram gerados com o mínimo de 10 ocorrências por bioma para cada espécie, e um script customizado foi desenvolvido para selecionar o melhor conjunto de variáveis ambientais para cada bioma onde a espécie é encontrada (Zurell et al. 2020).
Uma vez que cada bioma possui características ambientais distintas, a partir da metodologia utilizada espera-se que o espaço ambiental ocupado pelas espécies em cada bioma encontre-se bem representado nos modelos, gerando maior confiabilidade nos resultados. Os parâmetros necessários para gerar cada modelo foram definidos de forma espécie-específica a partir de testes estatísticos. A performance dos modelos e predições foram estimadas utilizando-se a técnica de validação cruzada, a partir de um script com funções adaptadas de Zurell et al. (2020). A performance dos modelos foi avaliada a partir da habilidade estatística verdadeira (TSS) e a área sob a curva característica do operador receptor (AUC), sendo os melhores modelos aqueles com valores de AUC próximos de 1 e de TSS próximos de 0,5.
Os modelos foram gerados em ambiente R (R Core Team 2021), e para mais detalhes sobre como os modelos foram gerados e pacotes R utilizados consulte Ávila et al. (2022). A partir dos modelos gerados, mapas de distribuição potencial de cada espécie para cada bioma foram criados. Os mapas foram gerados a partir da distribuição potencial contínua resultante dos modelos criados, e em formato binário, onde áreas de presença são determinadas a partir de um limiar de corte. Uma vez que o trabalho desenvolvido no escopo deste projeto visa o auxílio na conservação das espécies de morcegos no Brasil, optou-se por incluir mapas binários obtidos a partir de dois limiares de corte: um limiar mais restritivo (Máxima sensibilidade e especificidade), e um limiar menos restritivo (Probabilidade média). Para mais informações sobre limiares, consulte Liu et al. (2013).
