Avaliação da influência de fatores hidroclimáticos sobre a estrutura populacional de Podocnemis expansa (Testudines: Podocnemididae) no Rio Formoso, sudeste da Amazônia brasileira

Evaluation of the influence of hydroclimatic factors on the population structure of Podocnemis expansa (Testudines: Podocnemididae) in Formoso River, southeastern Brazilian Amazon

RESUMO

O aumento do conhecimento sobre a ecologia populacional das espécies de quelônios amazônicos é essencial para a recuperação e manutenção dos estoques naturais desses répteis. Neste estudo foi avaliada a influência de fatores hidroclimáticos sobre parâmetros populacionais e a condição corporal de Podocnemis expansa (Schweigger, 1812) em um trecho do Rio Formoso, sudeste da Amazônia brasileira, caracterizado pelo intenso uso das águas superficiais na irrigação de culturas agrícolas. O trabalho de campo ocorreu em três campanhas realizadas entre o final de julho e início de agosto dos anos de 2016, 2017 e 2018, totalizando um esforço amostral de 205 horas. Os indivíduos foram capturados por meio dos métodos de arrasto de rede, pesca com anzol e rede de espera; e os dados hidroclimáticos obtidos no site da Agência Nacional das Águas - ANA. Um total de 466 espécimes foi capturado e marcado, incluindo 47 machos adultos, 76 fêmeas adultas e 343 juvenis- revelando-se uma população composta predominantemente por juvenis (73,60%) e com razão sexual funcional desviada para as fêmeas (0,62:1). A redução no número de indivíduos das classes de tamanhos maiores e melhores condições corporais estiveram associadas, respectivamente, a uma maior vazão média do rio na estação seca e ao aumento da vazão no período chuvoso. Estudos futuros de longo prazo devem buscar aprofundar os conhecimentos sobre como fatores hidroclimáticos e as atividades agrícolas desenvolvidas na região afetam a demografia da espécie.

PALAVRAS-CHAVE
Tartaruga-da-Amazônia; ecologia populacional; vazão; precipitação; agricultura irrigada

ABSTRACT

Evaluation of the influence of hydroclimatic factors on the population structure of Podocnemis expansa (Testudines: Podocnemididae) in Formoso River, southeastern Brazilian Amazon. It is essential to have more knowledge about the population ecology of Amazonian chelonian species in order to recover and maintain the natural stocks of these reptiles. This study assessed the influence of hydroclimatic factors on population parameters and the body condition of Podocnemis expansa (Schweigger, 1812) in a specific region of Formoso River, in the southeast of Brazilian Amazon rainforest, where there is an intense use of surface water for irrigation of agricultural crops. The fieldwork took place in three campaigns carried out between the end of July and the beginning of August in the years 2016, 2017 and 2018, with a total sample effort of 205 hours. The individuals were captured using the methods of net trawling, hook fishing and waiting net; and the hydroclimatic data obtained from the website of the Agência Nacional de Águas - ANA. A total of 466 specimens were captured and tagged, being 47 adult males, 76 adult females and 343 juveniles, demonstrating a population composed mostly of juveniles (73.60%) and with a deviated functional sex ratio for females (0.62:1). The reduction in the number of individuals of larger sizes and better body conditions were associated, respectively, with higher river flow average in the dry season and an increase in flow during the rainy period. Future long-term studies are necessary to deepen the knowledge about how hydroclimatic factors and agricultural activities developed in the region affect the demography of the species.

KEYWORDS
Giant Amazon River Turtle; population ecology; flow; precipitation; irrigated agriculture

A pressão antrópica sobre os quelônios amazônicos tem se intensificado e fatores como a poluição das águas, a expansão das atividades agropecuárias, queimadas, desmatamentos das várzeas e matas ciliares, barramentos dos corpos d’água e a caça ilegal estão contribuindo para o declínio das populações desses répteis na região amazônica (Rodrigues, 2005; Alho et al., 2015). As mudanças climáticas são um risco adicional para as populações de quelônios amazônicos, tendo em vista que para muitas espécies a temperatura de incubação determina o sexo dos filhotes (Mitchell & Janzen, 2010; Stanford et al., 2020), e que os padrões de comportamento e atividades desses répteis estão fortemente associados às condições climáticas (Souza, 2004). Como resultado desse cenário de múltiplas ameaças, seis (40,00%) das 15 espécies de tartarugas de água doce que ocorrem na Amazônia brasileira estão em alguma categoria de ameaça (IUCN, 2016).

A recuperação e manutenção dos índices populacionais de quelônios amazônicos a níveis desejáveis estão condicionadas à continuidade, ao aprimoramento e ampliação dos trabalhos de proteção, manejo e pesquisa com esses animais (Salera-Júnior et al., 2016). Informações sobre a abundância, razão sexual e estrutura etária das populações de quelônios são essenciais para avaliar práticas de manejo, programas de conservação, bem como os impactos de fatores como doenças, caça e conversão de hábitat (Litzgus & Mousseau, 2004; Balestra et al., 2016; Ferrara et al., 2017; Schneider et al., 2018). Por outro lado, o conhecimento sobre a ecologia populacional das espécies de quelônios amazônicos permanece escasso (Fachín-Terán et al., 2003; Fachín-Terán & Vogt, 2004; Bernhard &Vogt, 2012; Bernardes et al., 2014; Miorando et al., 2015). Características desses répteis como o crescimento lento, a maturidade sexual tardia e a longevidade, tornam o estudo de suas populações um desafio (Shine & Iverson, 1995; Stanford et al., 2020).

Práticas associadas às atividades agrícolas como o desmatamento, controle/desvio das águas dos rios para irrigação, uso extensivo de agrotóxicos e fertilizantes - na medida em que deterioram, fragmentam e eliminam habitats, são reconhecidas como ameaças às populações de tartarugas de água doce (Rodrigues, 2005; Chessman, 2011; Usuda et al., 2012). Para os quelônios amazônicos, os efeitos dessas atividades intensificam mudanças no regime hídrico da região, contribuindo com cheias ou secas prolongadas e assoreamento dos rios, que podem causar graves distúrbios nos processos reprodutivos, migratórios e alimentares desses répteis (Luz et al., 2019). Alterações nos ciclos de chuvas e inundações dos rios, por exemplo, podem afetar de maneira negativa a reprodução das tartarugas de água doce na Amazônia, uma vez que a nidificação e incubação dos ovos dessas espécies ocorrem nos bancos de areia expostos ao longo de corpos hídricos durante a estação seca, e a eclosão dos filhotes coincide com o aumento dos níveis da água (Alho & Pádua, 1982; Eisemberg et al., 2016).

Mudanças nos regimes de vazão dos rios também têm efeitos sobre a abundância, razão sexual, estrutura etária e condição corporal dos quelônios (Usuda et al., 2012; Eisemberg et al., 2016), tendo em vista que a variabilidade do fluxo das águas exerce um papel importante na determinação da composição biótica, estrutura, função e diversidade dos ecossistemas fluviais (Bunn & Arthington, 2002; Döll & Zhang, 2010). Ashton et al. (2015), por exemplo, verificaram que tartarugas de um rio com o fluxo regulado por barragem cresceram mais lentamente e foram menores em idades semelhantes, quando comparadas com indivíduos de um rio sem barramento, associando essas diferenças às alterações causadas pelo controle do fluxo do rio.

Reduções e interrupções no fluxo dos rios podem alterar padrões sazonais, estabelecer condições lênticas e restringir o acesso a habitats reprodutivos e alimentares, afetando espécies migrantes de longa distância como Podocnemis expansa (Schweigger, 1812), conhecida popularmente como Tartaruga-da-Amazônia (Alho, 2011; Castello & Macedo, 2016). Essa espécie possui ampla distribuição geográfica, ocorrendo em rios e lagos das bacias do Amazonas e Orinoco na Colômbia, Venezuela, Guianas, Equador, Peru, Bolívia e Brasil, onde pode atingir a bacia Araguaia-Tocantins na região central (Ferrara et al., 2017). Em consequência da sobre-exploração de suas populações naturais e da constante degradação de seu habitat, o estado de conservação de P. expansa no território brasileiro é de quase ameaçada de extinção (NT) (ICMBio, 2018). A espécie também é listada no Apêndice II da Convenção sobre Comércio Internacional das Espécies da Fauna e Flora Selvagens Ameaçadas de Extinção (CITES, 2017); classificada pela International Union for Conservation of Nature (IUCN, 2016) como de baixo risco de extinção (LC), mas dependente de ações de conservação; e considerada criticamente ameaçada (CR) pelo Tortoise and Freshwater Turtle Specialist Group (Rhodin et al., 2017).

Poucos são os trabalhos que avaliam as populações de Podocnemis expansa em locais diretamente afetados por práticas agrícolas e em áreas protegidas (Figueroa et al., 2013; Hernández & Espín, 2006; Mogollones et al., 2010; Portelinha et al., 2014). Mediante o exposto, este estudo objetiva descrever a abundância, razão sexual, estrutura etária e condição corporal de P. expansa em um trecho do Rio Formoso caracterizado pela intensa utilização de suas águas para irrigação, e avaliar se a curto prazo variáveis hidroclimáticas explicam as possíveis variações desses parâmetros populacionais. Espera-se que a abundância, quantidade de adultos, proporção de fêmeas e a condição corporal de P. expansa estejam positivamente correlacionadas com as variáveis de precipitação e vazão devido a um aumento de produtividade do ecossistema aquático. Adicionalmente, pretende-se verificar se indivíduos com lesões recentes ou cicatrizadas teriam condições corporais inferiores quando comparados a indivíduos sem lesões.

MATERIAL E MÉTODOS

Área de estudo. O estudo foi realizado em um trecho de aproximadamente 10 km do Rio Formoso (coordenadas 10°43ʼ01.1ˮ- 10°45ʼ28.9ˮS e 49°51ʼ26.2ˮ- 49°51ʼ50.8ˮW), localizado no município de Lagoa da Confusão, Estado do Tocantins, sudeste da Amazônia brasileira (Fig. 1). A bacia do Rio Formoso situa-se na Região Hidrográfica do Tocantins-Araguaia, apresenta uma área de drenagem de 20.508,79 km² e abrange parte do território de 21 municípios dos estados do Tocantins e Goiás (Marinho-Filho et al., 2013). O clima predominante é o tropical quente semi-úmido (Gusmão et al., 2013), marcado por um período chuvoso de outubro a abril, uma estação seca de maio a setembro e precipitação média anual de 1579,6 mm (Alves et al., 2016). A região é constituída por um complexo ecótono que apresenta elementos vegetacionais do Cerrado, Floresta Amazônica e Pantanal (Brito et al., 2006), caracterizado pela alta diversidade de fauna e flora, e incluído na lista de sítios Ramsar no Brasil (Lopes et al., 2017; Brasil, 2019). Também se localiza no entorno do Corredor Ecológico Araguaia-Bananal - instrumento de gestão territorial projetado para promover a conexão entre quatro unidades de conservação e seis reservas indígenas existentes na região (Galinkin et al., 2003).

Fig. 1.
Localização da área de estudo no Rio Formoso, município de Lagoa da Confusão, Tocantins, Brasil. Hidrografia representada em azul; limites do município da Lagoa da Confusão/TO em vermelho; APA, Área de Proteção Ambiental; UC, Unidade de Conservação.

Os quinze municípios que representam a maior abrangência territorial (aproximadamente 97% da área) da bacia do Rio Formoso têm em conjunto uma população de 175.691 habitantes e densidade demográfica de 3,34 habitantes/km² (IBGE, 2019). Essa região compõe uma imensa área de várzea, com mais de 500 mil hectares, onde a agropecuária é a principal atividade econômica e desde a década de 1970 desenvolve-se o Projeto Rio Formoso - considerado o maior projeto de irrigação de arroz em área contínua da América Latina (Morais et al., 2017). Os projetos hidroagrícolas nessa área, por meio de barramentos e das técnicas de inundação e/ou subirrigação, utilizam intensamente as águas do Rio Formoso para cultivo de arroz, melancia, feijão, soja e outras culturas (Lima et al., 2013).

Os quatro barramentos existentes no Rio Formoso, localizados próximos a área do presente estudo, tiveram autorização pelo Instituto Natureza do Tocantins - NATURATINS (órgão estadual de fiscalização ambiental) - de construção e operação, respectivamente, nos anos de 2014 e 2018 (Tocantins, 2014, 2018). Esses barramentos são autovertentes; possuem uma base de concreto que se estende pelo leito do rio por aproximadamente seis metros a montante da barragem e cinco metros a jusante; têm estruturas projetadas para serem montadas sobre um eixo principal e removidas após o período de retenção da água na calha principal do rio, previsto para acontecer a cada ano entre julho e setembro (três últimos meses da estação seca). Foram construídos pela Associação dos Produtores Rurais do Rio Formoso de Lagoa da Confusão, com a finalidade de aumentar a oferta hídrica para irrigação de lavouras no período de estiagem (Carneiro-Neto et al., 2017).

Procedimentos de campo. O trabalho de campo ocorreu em três campanhas, cada uma com duração de sete dias, realizadas entre o final de julho e início de agosto, nos anos de 2016, 2017 e 2018. As campanhas foram realizadas em meses da estação seca (julho e agosto) visando capturar o maior número possível de espécimes e minimizar os custos, uma vez que a captura de quelônios aquáticos no período chuvoso é dificultada pelo maior volume de água nos rios (Balestra et al., 2016). As capturas aconteceram durante o dia, entre 08:00 h e 17:00 h, mediante autorizações SISBIO (48309-1e 58557-3) concedidas pelo Instituto Chico Mendes de Conservação da Biodiversidade - ICMBio. Foram amostrados cinco locais ao longo do trecho estudado no rio, dois caracterizados como reentrâncias, ambientes lênticos próximos às praias, com profundidade aproximada de um a dois metros; e três como remansos, ambientes no leito principal, com águas pouco correntes e profundidade aproximada de um a três metros.

Os métodos de captura, conforme Balestra et al. (2016), consistiram em: rede de espera - utilizando-se rede malhadeira tipo “transmalha” com 30 m de comprimento por 3 m de altura e distância entre nós de 12 cm (malha interna) e 40 cm (malhas externas); arrasto de rede - rede de náilon com 35 m de comprimento, 3 metros de altura e malha entre nós de 18 cm; e pesca - com linha de náilon comprida e anzol sem fisga, iscado com melancia ou mandioca. Os métodos permitem a captura de animais de todas as classes de tamanhos, exceto filhotes.

Os animais coletados foram submetidos a procedimentos de identificação sexual, biometria, avaliação das condições físicas e marcação individualizada. Seguindo Mogollones et al. (2010), o sexo foi determinado para indivíduos com comprimento curvilíneo da carapaça superior a 175 mm; por meio da observação visual do escudo anal do plastrão, que possui abertura em forma de U nos machos e de V nas fêmeas; e do tamanho da cauda, proporcionalmente maior nos machos. Conforme Portelinha et al. (2014), o comprimento e largura da carapaça e do plastrão foram obtidos de maneira curvilínea com fita métrica de 1 mm de precisão, enquanto a massa corporal foi verificada com balança de 100 g de precisão. Quando existentes nos espécimes capturados, as lesões e cicatrizes foram classificadas como sendo provenientes de tentativa de predação natural ou de colisão (em hélices dos motores/estruturas de embarcações; ou nas estruturas dos barramentos do rio) (Saumure & Bider, 1998). A marcação foi realizada de acordo com Balestra et al. (2016), perfurando-se os escudos marginais e os ossos periféricos da carapaça com furadeira portátil e brocas de diferentes espessuras, dependendo do tamanho do espécime. Para evitar o sangramento mais intenso do animal, não foram perfurados os escudos e os ossos periféricos que fazem a união da carapaça com o plastrão (do 4° ao 7° e do 18° ao 21° escudos marginais que coincidem com os 3° ao 7° ossos periféricos em ambos os lados da carapaça) (Pritchard & Trebbau, 1984). Os demais escudos marginais receberam até dois furos, compondo um sistema numérico no sentido horário, com unidades, dezenas, centenas e milhares. Os escudos perfurados na marcação foram lavados com álcool iodado e os animais devolvidos ao mesmo local em que foram capturados.

Parâmetros populacionais. Neste estudo considerou-se a definição de população estabelecida por Futuyma (1993), qual seja, conjunto de indivíduos da mesma espécie que coexistem em determinada área e período de tempo, bem como exibem continuidade reprodutiva. A razão sexual foi determinada por meio da divisão do número de machos pelo número de fêmeas em uma amostragem. Segundo Gibbons (1990), a razão sexual funcional deve ser calculada com base nos indivíduos que alcançaram a maturidade. Nesse sentido, foram considerados adultos os indivíduos com comprimento da carapaça ≥ 40 cm para os machos (Y. S. L. Bataus, dados inéditos) e ≥ 50 cm para fêmeas (Alho & Pádua, 1982). Para estrutura de tamanhos foram estabelecidas classes de tamanho com intervalos de 5 cm comprimento da carapaça (Portelinha et al., 2014). No cálculo da abundância relativa foi utilizada a taxa de captura por unidade de esforço (CPUE), na qual se divide o número de indivíduos capturados pelo esforço de captura empregado em horas (Fachín-Terán et al., 2003). O esforço de captura foi o resultado da multiplicação do número de petrechos (aparatos de captura) pela quantidade de horas que estes permanecerem ativos (Balestra et al., 2016). Conforme Peig & Green (2009), a condição corporal foi definida como sendo o índice de massa escalada (Scaled Mass Index - SMI) e calculada como: $$\widehat{M_i}=M_i{\left(\frac{C_0}{C_i}\right)}^{b_{SMA}}$$ , onde M_i e C_i são, respectivamente, a massa corporal e o comprimento da carapaça do indivíduo i; C₀ é a média do comprimento da carapaça da população da espécie estudada; e b_SMA é um fator de escala. O b_SMA foi computado de uma regressão de eixo principal padronizado (Standardized Major Axis - SMA) do ln(M + 1) sobre ln(C). O SMI é frequentemente utilizado como um índice de condição corporal e, ao contrário da maioria de outros métodos mais convencionais, pode considerar com sucesso a relação variante entre massa e comprimento corporal durante a ontogenia (Peig & Green, 2010), além de ser um bom preditor das variações das reservas de gordura e proteína em diversos grupos de vertebrados (Peig & Green, 2009). Neste sentido, indivíduos com melhores condições corporais devem apresentar valores maiores para esse índice e vice-versa, evidenciando por exemplo as condições de saúde e fitness do animal.

Análise estatística. Para selecionar as variáveis hidroclimáticas que mais explicam a variação da abundância relativa, razão sexual, estrutura etária e condição corporal de Podocnemis expansa, foi realizada uma seleção manual de modelos com base no Critério de Informação de Akaike (AIC) (Burnham & Anderson, 2002), partindo-se de um modelo vazio e inserindo individualmente as variáveis (Crawley, 2013; Heinze et al., 2018; Zuur et al., 2007). Foram compilados os dados de vazão e precipitação para a área de estudo no site da Agência Nacional das Águas Hidroweb (Ana, 2019) - das estações meteorológicas de números 1249002, 1249003, 26720000 e 26730000, todas localizadas na Bacia do Rio Formoso. Foram calculadas as médias dos valores de precipitação e vazão média referentes ao ano hidrológico, estações seca (maio a setembro) e chuvosa (outubro a abril), bem como para o período reprodutivo de P. expansa que na região ocorre de julho a outubro (Malvasio et al., 2019). Conforme Vergara et al. (2013), considerou-se para a área estudada o ano hidrológico com início no mês de outubro e término em setembro. Sendo assim, foram usados nos modelos e testes estatísticos um total de oito variáveis hidroclimáticas: precipitação e vazão anual média (Precip_med e Vazão_med); precipitação e vazão na época seca (Precip_seca e Vazão_seca), precipitação e vazão na época chuvosa (Precip_chuva e Vazão_chuva) e; precipitação e vazão no período reprodutivo (Precip_reprod e Vazão_reprod)(Tab. I).

Tab. I.
Valores médios ($$\stackrel{-}{x}$$ ) mensais e desvio padrão ($$\sigma$$ ) das variáveis hidroclimáticas avaliadas no Rio Formoso, estado do Tocantins, Brasil.

Variável	2016		2017		2018
	$$\stackrel{-}{x}$$	$$\sigma$$	$$\stackrel{-}{x}$$	$$\sigma$$	$$\stackrel{-}{x}$$	$$\sigma$$
Precipmed (mm)	101,36	177,48	112,51	117,71	136,67	131,75
Vazãomed (m³/s)	48,06	87,28	28,29	42,11	63,46	91,83
Precipseca (mm)	7,16	16,01	0,16	0,22	5,81	8,50
Vazãoseca (m³/s)	2,31	2,49	5,14	8,71	9,07	14,18
Precipreprod (mm)	44,96	68,17	22,14	44,04	24,35	35,93
Vazãoreprod (m³/s)	1,02	0,11	0,61	0,76	1,19	0,66
Precipchuva (mm)	168,64	211,89	192,76	85,84	230,14	85,51
Vazãochuva (m³/s)	80,74	103,50	44,83	48,76	102,31	104,30

Para relacionar a abundância relativa com as variáveis hidroclimáticas, foram usados modelos lineares generalizados (GLM) com a distribuição de erros binomial negativa (Venables & Ripley, 2002). O esforço total de cada amostragem e o tipo de método utilizado (rede de arrasto, rede de espera ou pesca) foram incluídos como covariáveis nos modelos. Para a razão sexual, foram usados GLMs com a distribuição de erros binomiais (Venables & Ripley, 2002; McCullagh & Nelder, 1989). Para a estrutura etária, foram consideradas o comprimento da carapaça (usado como um proxy da estrutura etária), já que o comprimento de quelônios é diretamente proporcional à idade dos indivíduos (Zug, 1991).

Primeiramente, a estrutura etária foi comparada entre os anos e entre os sexos por meio de uma análise de variância (ANOVA) ranqueada (não-paramétrica; teste de Shapiro-Wilk: W = 0,81297, p <0,001) e posteriormente foi realizada a seleção das variáveis hidroclimáticas por meio de um modelo linear simples. Após a seleção das variáveis hidroclimáticas, foi realizada uma regressão quantílica (Koenker & Bassett-Júnior, 1978; Koenker & Hallock, 2001; veja a revisão do método em Koenker, 2017). A condição corporal dos indivíduos foi comparada por meio de um teste não-paramétrico de Wilcoxon (teste de Shapiro-Wilk: W = 0,89668, p <0,001) entre indivíduos com e sem marcas de lesões recentes ou cicatrizadas. Por meio de ANOVA ranqueada, foram testadas a relação da condição corporal com os anos do estudo, sexos e fases (adulto e juvenil) dos indivíduos. Em seguida, foi realizada a seleção manual de variáveis hidroclimáticas que melhor explicam a condição corporal por meio de modelos lineares simples, incluindo como covariável o comprimento da carapaça. Todas as análises estatísticas consideraram um nível de significância de 5% e foram realizados no software livre R v. 3.6.3 (R Core Team, 2020). Para realizar a maioria das análises estatísticas foram usadas as funções do pacote stats do R, padrão no programa. Para construir os modelos lineares generalizados com erros na distribuição binomial negativa, foi usado o pacote MASS (Venables & Ripley, 2002); e para construir as regressões quantílicas foi aplicado o pacote quantreg (Koenker, 2020). Outras informações sobre a análise estatística foram incluídas no ^{Material Suplementar}.

RESULTADOS

Capturas e abundância. Foram capturados e marcados 466 espécimes de Podocnemis expansa (Tab. II), tendo sido recapturados (uma única vez) sete indivíduos (três machos juvenis, três machos adultos e uma fêmea juvenil) - o que equivale a uma taxa de recaptura de aproximadamente 1.50%. A abundância relativa não variou significativamente entre os anos (χ²= 1,985; gl = 2; p = 0,371), sendo de 1.35 (indivíduos/hora) para 2016, 2.41 para 2017 e 3.23 para 2018 (Tab. II). Por outro lado, houve variação na abundância relativa entre os métodos (χ² = 75,529; gl = 4; p < 0,001), tendo sido capturados 0.09 (indivíduos/hora) por meio da rede de espera, 2.83 por pesca e 5.57 por arrasto de rede.

Tab. II.
Total de indivíduos capturados (n), abundância relativa, razão sexual e estrutura etária de Podocnemis expansa entre os anos de 2016 e 2018, no Rio Formoso, estado do Tocantins, Brasil.

Ano	Capturados			Abundância relativa	Razão sexual						Estrutura etária
					juvenis			adultos			juvenis	adultos
	♂	♀	n		♂	♀	♂:♀	♂	♀	♂:♀
2016	29	85	114	1,35	13	50	0,26:1	16	35	0,46:1	55,26% (n=63)	44,74% (n=51)
2017	35	76	111	2,41	23	54	0,43:1	12	22	0,55:1	69,37% (n=77)	30,63% (n=34)
2018	89	152	241	3,23	70	133	0,53:1	19	19	1:1	84,23% (n=203)	15,77% (n=38)
Total	153	313	466	2,27	106	237	0,45:1	47	76	0,62:1	73,61% (n=343)	26,39% (n=123)

Razão sexual e estrutura etária. Do total de exemplares capturados, verificou-se razão sexual funcional de 0.62:1 desviada para as fêmeas (χ² = 6,837; gl = 1; p = 0,0116) e de 0,45:1 entre os juvenis, também desviada para as fêmeas (χ² = 50,032; gl = 1; p < 0,001)(Tab. II). Foram 26,40% de indivíduos adultos (χ² = 103,863; gl = 1; p < 0,001), incluindo 47 machos e 76 fêmeas; e 73,60% de juvenis, compostos por 106 machos e 237 fêmeas (Tab. II). O comprimento médio da carapaça foi de 30,24 ± 5,14 cm (variação: 19,00 - 48,90 cm) para fêmeas juvenis e 68,42 ± 5,32 cm (variação: 52,50 - 78,50 cm) para adultas, enquanto entre os machos juvenis foi de 35,02 ± 2,70 cm (variação: 27,80 - 39,80 cm) e 43,81 ± 3,51 cm (variação: 40,00 - 53,30 cm) entre os adultos. Os métodos de captura não foram seletivos para o sexo (χ² = 0,11551; gl = 2; p = 0,9439), bem como não houve variação significativa na razão sexual entre os anos de amostragem (χ² = 2,5269; gl = 2; p = 0,2827). Apesar disso, destaca-se que em 2018 foi encontrada razão 1:1 em adultos (Tab. II).

A estrutura de tamanho apresentou um padrão unimodal para os machos, sendo a maior frequência encontrada entre os indivíduos com comprimento da carapaça correspondente à classe de tamanho de 35 a 40 cm (Fig. 2). Para as fêmeas a tendência foi de uma curva bimodal, com destaques nas classes de tamanhos de 25 a 30, 65 e 75 cm de comprimento da carapaça (Fig. 2). Apesar de terem sido utilizados três diferentes métodos de captura, recém-eclodidos e juvenis menores que 19 cm (classe de tamanho de 15 cm) não foram capturados.

Fig. 2.
Distribuição de frequências das classes de tamanhos das carapaças de machos e fêmeas de Podocnemis expansa (Schweigger, 1812) capturados entre os anos de 2016 e 2018, no Rio Formoso, Estado do Tocantins, Brasil.

Verificou-se que a estrutura de classes de tamanho variou entre os anos (F = 9,447; gl = 2; p < 0,001) (Fig. 3). Não houve interação significativa entre os anos e o sexo, ou seja, a estrutura etária variou igualmente entre anos em ambos os sexos (F = 0,062; gl = 2; p = 0,94). A vazão na estação seca (média ± desvio padrão em 2016 = 2,31 ± 2,49 m³/s; 2017 = 5,14 ± 8,71 m³/s; e 2018 = 9,07 ± 14,18 m³/s) (Tab. I) foi a variável que melhor explicou a variação na estrutura etária entre os anos (modelo linear simples: AIC = 2531,9; χ² = 6899,1; gl = 1; p < 0,001; regressão quantílica: F = 5,923; gl =1; p < 0,001). Quanto menor a vazão na estação seca maior a frequência de indivíduos de classes de tamanhos maiores (Fig. 3).

Fig. 3.
Relação entre a vazão na estação seca (linha azul) e o tamanho de Podocnemis expansa (Schweigger, 1812) entre os anos de 2016 e 2018, no Rio Formoso, Estado do Tocantins, Brasil.

Condição corporal. O percentual de exemplares com pelo menos uma lesão ou cicatriz foi de 34,55% (n = 161). As maiores frequências foram para lesões e/ou cicatrizes provenientes de tentativas de predação natural (82 injúrias no casco, 53 faltas de dedos e 29 faltas de membros). No entanto, foram encontrados cinco indivíduos com fratura por colisão na parte anterior do plastrão. Espécimes com lesões recentes e/ou cicatrizadas apresentaram maior condição corporal quando comparados aos que não tinham lesões (W = 23232; p= 0,04494; sem lesões - SMI = 5,14 ± 0,99; com lesões - SMI=5,31±0,98). A condição corporal dos indivíduos também variou entre as classes etárias, sendo maior para os adultos (F = 279,9; gl = 1; p < 0,001; juvenis - SMI = 4,81 ± 0,61; adultos - SMI = 6,27 ± 1,03).

A proporção de indivíduos com lesões e/ou cicatrizes não variou entre os sexos (χ² = 0,758; gl = 1; p = 0,4434), assim como não houve alteração para essa variável entre os anos (χ² = 2,0404; gl = 2; p = 0,3605). A condição corporal dos indivíduos independe do sexo (F = 0,166; gl = 1; p = 0,684; machos - SMI = 5,02 ± 0,48; fêmeas - SMI = 5,29 ± 1,15), mas variou entre os anos de estudo (F = 10,88, gl = 2, p < 0,001). O teste de Conovera posteriori mostrou que o ano de 2017 teve menor condição corporal (2016 - SMI = 5,53 ± 1,17; 2017 - SMI = 4,92 ± 0,99; 2018 - SMI = 5,16 ± 0,83; p < 0,001). A vazão no período chuvoso (média ± desvio padrão em 2016 = 80,74 ± 103,50 m³/s; 2017 = 44,83 ± 48,76 m³/s; e 2018 = 102,31 ± 104,30 m³/s) (Tab. I) foi a variável que melhor explicou a variação na condição corporal entre os anos (modelo linear simples: AIC=602,50; χ² =17,4684; gl=1; p<0,001; regressão quantílica: F = 3,122; gl = 1; p = 0,014), tendo sido verificado que quanto maior a vazão no período chuvoso maior a condição corporal (Fig. 4).

Fig. 4.
Relação entre a vazão na estação chuvosa (linha azul) e a condição corporal de Podocnemis expansa (Schweigger, 1812) entre os anos de 2016 e 2018, no Rio Formoso, Estado do Tocantins, Brasil.

DISCUSSÃO

Capturas e abundância. Assim como relatado por Portelinha et al. (2014) (recaptura de 0,70%), no presente estudo se verificou um baixo índice de recaptura (aproximadamente 1,50%) para P. expansa. Baixos índices de recaptura também foram registrados para P. erythrocephala (Bernhard & Vogt, 2012) (< 10,00%), P. sextuberculata (Fachín-Terán et al.,2003) (3,5%) e P. unifilis (Fachín-Terán & Vogt, 2004) (5,00%). De acordo com Bernhard & Vogt (2012), um baixo número de recapturas pode indicar que apenas uma pequena fração de indivíduos foi capturada em relação ao tamanho total da população. Isso pode ocorrer devido a fatores como baixa eficiência dos métodos de captura, amostragem de apenas uma pequena parte da área ocupada pela população estudada ou ainda em razão de altas taxas de imigração/migração (Balestra et al., 2016).

A abundância relativa variou de acordo com o método. Considerando o método mais eficiente e o período total amostrado, o resultado registrado no presente estudo (5,57 espécimes capturados por hora de amostragem com o arrasto de rede) foi inferior ao relatado por Portelinha et al. (2014) (10,9 tartarugas capturadas por hora, também por meio do arrasto de rede), o que pode ser explicado pelo fato do trabalho desses autores ter sido realizado em uma área protegida. Diversos autores relatam que locais protegidos, quando comparados a áreas sem proteção, apresentam maiores abundâncias relativas de tartarugas (Fachín-Terán et al., 2003; Tomillo et al., 2008; Alcantara et al., 2013; Sung et al.2013). Hernández & Espín (2006) por sua vez, após três anos de amostragem em estudo realizado na Venezuela, relataram diminuição na abundância relativa de P. expansa (variação de 3,7 para 1,2 indivíduos capturados/100 m de rede de arrasto, respectivamente, no primeiro e terceiro ano). Eles explicaram que essa redução na taxa de captura poderia ser resultado de diferenças nos métodos de amostragem e nas atividades de caça. Usuda et al. (2012) relataram acentuada diminuição da densidade de Chinemys reevesii (Gray, 1831) (Testudines, Geoemydidae) em consequência de práticas artificiais de alteração dos rios para controle de inundação, na província de Chiba, Japão. De modo parecido, Chessman (2011) apontou a diminuição do fluxo dos rios (devido a um longo período de seca e o uso da água na agricultura irrigada) como provável causa do declínio da abundância das espécies de Chelidae Chelodina longicollis (Shaw, 1794) e Emydura macquarii (Gray, 1830), em estudo realizado na Austrália.

Não foi possível verificar se as variáveis hidroclimáticas tinham relação com mudanças na abundância relativa de Podocnemis expansa, pois não houve variação significativa nesse parâmetro populacional ao longo do período avaliado. Ao avaliarem a influência de fatores ambientais na reprodução de P. expansa no Rio Javaés, Estado do Tocantins (TO), Simoncini et al. (2019) constataram que quanto maior a precipitação na nascente do rio, maior a quantidade de fêmeas desovando nas praias do local. Em uma lagoa no Parque Nacional de Brasília/DF, Fraxe-Neto et al. (2011) associaram maiores taxas de captura de Acanthochelys spixii (Duméril & Bibron, 1835) (Testudines, Chelidae) ao aumento da abundância e riqueza de insetos que fazem parte da dieta dessa espécie, sugerindo que isso estaria relacionado ao aumento dos níveis de água no local. Desse modo, estudos futuros devem avaliar a longo prazo os efeitos da vazão e precipitação sobre abundância de P. expansa no Rio Formoso, especialmente porque as águas superficiais da região são cada vez mais utilizadas para irrigação de culturas agrícolas (Morais et al., 2017).

Razão sexual e estrutura etária. A razão sexual de P. expansa desviada para as fêmeas (0,45♂:1♀ entre os juvenis; e 0,62♂:1♀ entre os adultos) verificada no presente estudo também foi observada por Hernández & Espín (2006) (0,03♂:1♀), Mogollones et al. (2010) (0,09♂:1♀) e Peñaloza et al. (2013) (0,16♂:1♀) - todas com capturas de P. expansa por arrasto de rede no médio Orinoco, Venezuela. Por outro lado, Portelinha et al. (2014) em estudo realizado no Rio Javaés, sul da Amazônia brasileira, utilizando os métodos de mergulho e rede de arrasto, encontraram uma razão sexual em P. expansa desviada para os machos entre os adultos (1,4♂:1♀) e para as fêmeas entre os juvenis (0,4♂:1♀). Existem registros de proporções sexuais desviadas para os machos em estudos populacionais de P. sextubercularta (Fachín-Terán et al., 2003) e P. unifilis (Fachín-Terán & Vogt, 2004; Miorando et al., 2015), bem como de desvio para fêmeas em P. erythrocephala (Bernhard & Vogt, 2012). Desvios na razão sexual em populações de quelônios podem ser explicados por diferenças na proporção sexual dos indivíduos recém-eclodidos, devido à determinação do sexo dependente da temperatura (Valenzuela, 2001; Bonach et al., 2011). Proporções sexuais enviesadas também podem ser resultado de diferenças na idade/tamanho de maturação sexual, nas taxas de movimentação/migração, mortalidade e sobrevivência de machos e fêmeas, e/ou consequência dos métodos de captura, período do ano e locais de amostragem utilizados (Gibbons, 1990; Steen et al., 2006).

Os métodos de captura utilizados no presente estudo não foram seletivos para o sexo. Nesse sentido, entre os fatores citados anteriormente, o mais provável é que o maior número de fêmeas em relação aos machos verificado na população estudada seja resultado de uma razão sexual diferencial dos indivíduos recém-eclodidos (também em favor das fêmeas) ao longo do tempo. Razões sexuais primárias de Podocnemis expansa desviadas para fêmeas foram verificadas por Malvasio et al. (2002) (n = 494; 0,7♂:1♀) e Malvasio et al. (2012) (n = 187; 0,2♂:1♀) em estudos realizados no Rio Javáes/TO, do qual o Rio Formoso é afluente; e por Valenzuela et al. (1997) (n = 173; 0,1♂:1♀) na Colômbia. Do ponto de vista da conservação, a razão sexual desviada para as fêmeas pode ser considerada positiva, uma vez que um único macho é capaz de fertilizar várias fêmeas e assim manter ou incrementar uma população (Vogt, 1994). No entanto, uma alta diminuição no número de machos devido ao aumento da temperatura de incubação dos ninhos pode dificultar o encontro intersexual para acasalamento, possivelmente diminuindo a persistência das populações a longo prazo (Janzen, 1994). Esse fenômeno também é preocupante entre as espécies de tartarugas marinhas, sobretudo em razão do risco de feminização das populações dessas espécies associado ao aumento das temperaturas devido às mudanças climáticas (Glen & Mrosovsky, 2004; Hawkes et al., 2007; Fuller et al., 2013; Laloë et al., 2016; Jensen et al., 2018; Stanford et al., 2020).

Não foi possível avaliar se as variáveis hidroclimáticas explicavam variações na razão sexual da população estudada, tendo em vista que não houve mudança significativa na proporção dos sexos entre os anos de amostragem. Nesse sentido, estudos futuros devem testar quais os efeitos de variáveis hidroclimáticas (além das que testamos nesse estudo, por exemplo aquelas relacionadas à temperatura) na proporção dos sexos da população de Podocnemis expansa, sobretudo entre os recém-eclodidos e a longo prazo.

Alcantara et al. (2013) afirmam que as variações climáticas e sua influência nas proporções sexuais dos filhotes devem ser monitoradas continuamente, devendo-se tomar cuidado para interpretar os resultados em estudos de curto prazo. Eisemberg et al. (2016) também apontam que mudanças nos ciclos de chuvas e inundações têm impactos diretos nas tartarugas de água doce da Amazônia e citam, por exemplo, a redução no número de dias da estação reprodutiva e alterações nos níveis das águas dos rios como responsáveis por afetar as características térmicas dos ninhos e, em consequência, distorcer a proporção dos sexos dos filhotes, uma vez que muitas espécies possuem determinação sexual dependente da temperatura. Nesse sentido, a manutenção de níveis elevados das águas dos rios por um maior período de tempo poderia diminuir a temperatura nos locais de nidificação de espécies como P. expansa, desviando a razão sexual dos filhotes para os machos, uma vez que nessa espécie temperaturas de incubação mais altas favorecem a geração de fêmeas e temperaturas mais baixas a produção de machos (Valenzuela, 2001; Bonach et al., 2011).

A população de P. expansa estudada é composta em sua maioria por juvenis (73,60%), assim como o que foi verificado para a mesma espécie em estudos realizados na Venezuela (Mogollones et al., 2010; Peñaloza et al., 2013) e Colômbia (Figueroa et al., 2013). Nos trabalhos conduzidos na Venezuela, os autores argumentaram que ações de proteção das praias de nidificação e de reintrodução de filhotes criados em cativeiro, realizadas na área estudada (médio Orinoco), eram as prováveis causas da maior proporção de tartarugas juvenis na população; enquanto na Colômbia, o desvio para os juvenis teria ocorrido por viés na seleção do local e período de amostragem. Portelinha et al. (2014) encontraram predominância de adultos (55,47%) ao avaliarem períodos em que indivíduos de P. expansa foram capturados apenas na estação reprodutiva e argumentaram que isso poderia ser explicado pelas diferenças no comportamento migratório dos grupos etários. Em geral, a abundância de juvenis é uma das características da estrutura populacional de répteis (Cagle, 1954). Gibbs & Amato (2000) apontaram que a classe juvenil é a mais abundante nas populações de muitas espécies de tartarugas. No presente estudo, não foi demonstrado que populações de tartarugas afetadas por distúrbios tiveram seu recrutamento reduzido, conforme mencionado em outras pesquisas (Marchand & Litvaitis, 2004; Chessman, 2011). Vale ressaltar que porcentagens mais altas de jovens podem não resultar necessariamente em populações estáveis (Crouse & Frazer, 1995), especialmente para a população de P. expansa avaliada na região em estudo, onde as taxas de sobrevivência para os diferentes grupos etários são desconhecidas.

O comprimento médio da carapaça verificado para as fêmeas adultas de P. expansa capturadas no presente estudo (68,42 ± 5,32 cm) foi semelhante aos apresentados por Alho & Pádua (1982) (66,00 ± 3,06 cm) e Portelinha et al. (2014) (64,70 ± 7,10 cm). Para os machos o resultado do nosso estudo (comprimento médio da carapaça igual a 43,81 ± 3,51 cm) também foi similar ao relatado por Portelinha et al. (2014) (41,90 ± 5,70 cm). Embora os tamanhos médios de P. expansa verificados no presente estudo tenham variado pouco em relação aos apresentados pelos autores citados, destaca-se que aspectos da história de vida, como o tamanho corporal, tendem a variar na distribuição de uma espécie (Greaves & Litzgus, 2009). Para muitas espécies de tartarugas o tamanho corporal pode estar relacionado à latitude e a fatores ambientais (temperatura, produtividade e sazonalidade) (Ashton & Feldman, 2003; Werner et al., 2016), sendo esperado que existam tartarugas maiores em latitudes mais baixas e menores em latitudes mais altas (Angielczyk et al., 2015). De acordo com Sung et al. (2013), as diferenças no tamanho médio ao longo da área de distribuição de uma espécie também podem ser devidas à pressão de caça, sendo que locais com histórico de caça ilegal de tartarugas apresentam indivíduos de tamanhos corporais menores do que os de locais protegidos. Na região do Corredor Ecológico Araguaia-Bananal, onde se realizou o presente estudo, Faria & Malvasio (2018) verificaram que os quelônios sofrem intensa pressão de caça, especialmente P. expansa.

O número reduzido de indivíduos de P. expansa (n = 1) com comprimento da carapaça menor ou igual à classe de tamanho de 15 cm pode ser devido ao uso de habitats diferentes do restante da população, conforme observado por Portelinha et al. (2014) para a mesma espécie, e por Fachín-Terán & Vogt (2004) para Podocnemis unifilis (Troschel, 1848). Miorando et al. (2015) explicaram que os métodos de captura muitas vezes não são adequados na amostragem dos microhabitats utilizados pelos juvenis de tamanho pequeno, motivo pelo qual estimar a quantidade desses indivíduos nas populações é um grande desafio para os estudos demográficos de Podocnemis. As fêmeas das classes de tamanhos 45, 50 e 55 cm também tiveram poucas representantes capturadas (apenas uma para cada classe), o que pode ser um indício de maior pressão de caça sobre esse sexo e classes de tamanho. Ataídes et al. (2010) relataram que os espécimes de tamanho médio são preferencialmente consumidos pela maioria das pessoas que fazem uso da carne dos quelônios no Rio Javaés/TO. A sobre-exploração de ovos, captura de adultos e comercialização de quelônios amazônicos têm afetado de maneira negativa a distribuição das classes de tamanho nas populações (Fachín-Terán & Vogt, 2004; Hernández & Espín, 2006; Miorando et al., 2015).

A vazão média do rio no período seco (maio a setembro) foi a variável que melhor explicou a variação na estrutura de classes de tamanhos entre os anos (Fig. 3). Era esperado que o aumento da disponibilidade de alimentos e a facilidade de deslocamento, proporcionados pela maior vazão no período seco, contribuíssem de maneira diretamente proporcional com as variações na estrutura de classes de tamanhos de P. expansa. No entanto, os resultados demonstraram que a redução no número de indivíduos das classes de tamanhos maiores está relacionada de maneira inversa com a vazão média do rio na estação seca.Os quatro barramentos construídos no rio Formoso, a partir do ano de 2014, podem ter impedido que parte dos indivíduos de P. expansa completasse o processo de migração dos habitats alimentares para a área do presente estudo, tendo em vista que as estruturas de represamento são operadas de julho a setembro (meses que fazem parte do período reprodutivo). As práticas de alteração de rios, incluindo canalização, barramentos e retirada de água para irrigação, afetam negativamente as características demográficas das populações de tartarugas (Alho, 2011; Chessman, 2011; Usuda et al., 2012). A pressão de caça diferencial sobre os indivíduos de tamanhos maiores também pode ser umas das causas da redução no número de espécimes maiores. Indícios dessa pressão diferencial foram verificados nos dois primeiros anos do presente estudo, quando em um dos locais de amostragem foram encontrados cascos de 30 indivíduos de P. expansa (Fig. 5), dos quais identificou-se que oito eram machos e 20 fêmeas, respectivamente, com 43,20 cm e 68,21 cm de tamanho médio da carapaça, todos predados por humanos.

Fig. 5.
Carapaça e plastrão de treze espécimes adultos de Podocnemis expansa (Schweigger, 1812) predados por humanos no ano 2017, em um dos pontos de amostragem no Rio Formoso, Estado do Tocantins, Brasil.

Condição corporal. A condição corporal de um animal é definida como uma medida do seu estado energético, especialmente associada ao tamanho relativo das reservas de gordura e proteína (Schulte-Hostedde et al., 2001, 2005). Muitos aspectos comportamentais e processos fisiológicos que influenciam direta ou indiretamente o sucesso reprodutivo, como a seleção de habitat, forrageamento, crescimento e migração, podem ser fortemente afetados pela condição corporal (Peig & Green, 2009; Thomson et al., 2009). Variações nas condições corporais de quelônios têm sido correlacionadas a componentes climáticos (Henen, 1997; McCoy et al., 2011), diferenças nos níveis de conservação do habitat (Polo-Cavia et al., 2010) e a diferenças entre os grupos etários nas atividades de migração e reprodução (Jessop et al., 2004).

A variação entre os anos na condição corporal dos espécimes de P. expansa capturados foi associada à vazão no período chuvoso (outubro a abril) (Fig. 4). Nos anos em que houve aumento da vazão do rio no período chuvoso, os indivíduos apresentaram melhores índices de condição corporal. Howard et al. (2017) atribuíram o aumento da condição corporal de Chelodina longicollis (Shaw, 1794) (Testudines, Chelidae) ao restabelecimento da vazão em um rio anteriormente submetido à seca severa, no sudeste da Austrália. Assim como relatado por esses autores, o aumento na condição corporal verificada no presente estudo pode ser resultado de uma maior produção de alimentos e/ou de uma redução na competição intraespecífica por alimentos, possibilitados pelo aumento da vazão no período chuvoso. O presente estudo não avaliou de maneira direta os efeitos dos barramentos do Rio Formoso sobre a condição corporal dos indivíduos de P. expansa, mas os resultados obtidos podem servir de linha de base para estudos futuros. Ashton et al. (2015), por exemplo, verificaram que alterações do fluxo e regimes térmicos causadas por represamento tiveram efeitos negativos significativos e de longo prazo na condição corporal e reprodução de tartarugas Actinemys marmorata (Baird & Girard, 1852) (Testudines, Emydidae), no rio Trinity, noroeste da Califórnia.

Indivíduos de P. expansa adultos e com lesões apresentaram melhores condições corporais. Outros estudos também apontaram que quelônios adultos exibiram índice de condição corporal mais elevado do que os juvenis e explicaram que isso se deve ao comportamento de acúmulo de energia por parte dos adultos para investimento na reprodução (Henen, 1997; Jessop et al., 2004). Cecala et al. (2009) verificaram que espécimes de Malaclemys terrapin (Schoepff, 1793) (Testudines, Emydidae) com lesões no casco não sofreram redução na condição corporal em comparação a indivíduos sem lesões. De modo de semelhante, Meek (2007) e Biaggini & Corti (2018) observaram que indivíduos de Testudo hermanni (Gmelin, 1789) (Testudines, Testudinidae) não tiveram a condição corporal afetada por lesões no casco ou perda de membros. No entanto, vale destacar que indivíduos com lesões podem apresentar desvantagens nos conflitos intraespecíficos, tempo de forrageamento e capacidade de fuga, afetando negativamente a aptidão reprodutiva ou longevidade (Smith, 1995; Downes & Shine 2001; Taylor & Jackson, 2003).

A maior probabilidade de indivíduos adultos de P. expansa apresentarem lesões recentes e/ou cicatrizadas pode ser explicada pelo fato desse grupo etário possuir taxa de sobrevivência maior (Mogollones et al., 2010), especialmente porque a maioria das injúrias/cicatrizes verificadas neste estudo era proveniente de tentativas de predação natural. Por outro lado, as fraturas na parte anterior do plastrão, verificadas em cinco indivíduos adultos de P. expansa e que não eram oriundas de tentativas de predação natural, podem ter sido causadas pelas estruturas dos barramentos existentes no rio. Os espécimes podem ter fraturado o plastrão ao transporem as barragens (nos períodos em que estavam vertendo água), uma vez que as estruturas de represamento são autovertentes, têm até dois metros de altura de lâmina d’água a montante e, imediatamente a jusante, possuem uma base de concreto e nível da água reduzido. Em estudo realizado no rio Mary, Austrália, Limpus (2008) verificou que 43% dos indivíduos capturados (n=456) de Emydura macquariikrefftii (Gray, 1871) (Testudines, Chelidae) apresentavam fraturas na carapaça ou plastrão resultantes de colisão com superfícies duras (concreto, metal) e que ocorriam quando esses animais passavam por cima da infraestrutura de represamento do rio. É importante que estudos futuros avaliem a longo prazo as condições corporais e injúrias para P. expansa no Rio Formoso, tendo em vista que os efeitos de lesões e/ou más condições corporais podem aumentar a probabilidade de declínio das populações (Litzgus et al., 2008; Bennett & Litzgus, 2014).

Este estudo foi o primeiro a avaliar parâmetros da população de uma espécie de quelônio no Rio Formoso. Os resultados evidenciam que a população de P. expansa na área estudada é composta predominantemente por juvenis, possui razão sexual desviada para fêmeas e abundância inferior à verificada em locais protegidos. Além disso, indivíduos com lesões recentes e/ou cicatrizadas apresentaram melhores condições corporais, indicando que as lesões não prejudicaram a aptidão de P. expansa. Os resultados demonstram ainda que a redução no número de indivíduos das classes de tamanhos maiores e melhores condições corporais estiveram associadas, respectivamente, a uma maior vazão média do rio na estação seca e ao aumento da vazão no período chuvoso.

Embora os resultados revelem que a população de P. expansa no Rio Formoso mantém seu recrutamento e que a razão sexual desviada para as fêmeas pode ser positiva para a conservação, os indícios de pressão de caça diferencial sobre os indivíduos adultos na área estudada e a provável diminuição do número de machos em razão do aumento na temperatura de incubação dos ovos (Jansen, 1994; Mitchell & Jansen, 2010) podem contribuir para o declínio da população a longo prazo. O intenso uso das águas superficiais para irrigação de culturas agrícolas na área de estudo (Lima et al., 2013; Morais et al., 2017) e alterações no ciclo hidrológico devido às mudanças climáticas (Guimberteau et al., 2013) também podem comprometer a manutenção das populações de P. expansa, na medida em que afetam de maneira negativa a vazão do rio. Desse modo, diagnósticos para a implementação e/ou aperfeiçoamento de planos de manejo e conservação também precisam ter em vista os efeitos das variáveis hidroclimáticas sobre a dinâmica das populações da referida espécie. Ademais, considerando que a bacia do Rio Formoso é constituída por ambientes heterogêneos (rios, lagos, praias, igapós e várzeas) importantes para o ciclo de vida de P. expansa e outras espécies de quelônios, estudos futuros de longo prazo devem buscar entender melhor a influência das variáveis hidroclimáticas sobre os parâmetros populacionais, conhecer os padrões de migração, movimentação e uso do espaço, bem como a demografia da espécie, sobretudo para compreender e propor medidas que mitiguem os impactos deletérios do avanço das atividades agrícolas na região.

Agradecimentos

Agradecemos à Universidade Federal do Tocantins - UFT pelas oportunidades promovidas por meio do Programa de Pós-Graduação em Ciências do Ambiente. Ao Programa Nacional de Cooperação Acadêmica - PROCAD/CAPES, pelo apoio logístico no projeto “Ecologia da Transição Cerrado-Floresta Amazônica” - processo nº 88881.068430-2014-01, edital nº 071/2013, coordenado pelo Dr. Guarino Rinaldi Colli, da Universidade de Brasília - UNB. Ao Instituto Chico Mendes de Conservação da Biodiversidade - ICMBio pelas autorizações no Sistema de Autorização e Informação em Biodiversidade - SISBIO de nº 48309-1e nº 58557-3. Aos proprietários da Fazenda e Pousada Praia Alta, Eloi Bernardon e Eduardo Bernardon, pelo apoio logístico. À Prof. Natani Cruz pela revisão da língua inglesa do Abstract. Aos amigos do grupo de pesquisa Crocodilianos e Quelônios da Região Norte (CROQUE-UFT) pela indispensável ajuda em campo.

REFERÊNCIAS

• Alcantara, A. S.; Silva, D. F. & Pezzuti, J. C. B. 2013. Effects of the hydrological cycle and human settlements on the population status of Podocnemis unifilis (Testudines: Podocnemididae) in the Xingu River, Brazil. Chelonian Conservation and Biology 12(1):134-142.
• Alho, C. J. & Pádua, L. F. 1982. Reproductive parameters and nesting behavior of the Amazon turtle Podocnemis expansa (Testudinata: Pelomedusidae) in Brazil. Canadian Journal of Zoology 60(1):97-103.
• Alho, C. J.; Reis, R. E. & Aquino, P. P. 2015. Amazonian freshwater habitats experiencing environmental and socioeconomic threats affecting subsistence fisheries. Ambio 44(5):412-425.
• Alho, C. J. R. 2011. Environmental effects of hydropower reservoirs on wild mammals and freshwater turtles in Amazonia: a review. Oecologia Australis 15(3):593-604.
• Alves, K. C. C. L. F.; Viola, M. R.; Mello, C. R.; Giongo, M. & Santos, A. F. 2016. Distribuição da precipitação mensal, anual e máxima diária anual na bacia hidrográfica do rio Formoso, Tocantins. Ambiência 12(1):49-70.
• Ana - Agência Nacional das Águas. 2019. HIDROWEB v3.1.1. Disponível em <Disponível em http://www.snirh.gov.br/hidroweb/apresentacao >. Acesso em 03.12.2019.
  » http://www.snirh.gov.br/hidroweb/apresentacao
• Angielczyk, K. D.; Burroughs, R. W. & Feldman, C. R. 2015. Do turtles follow the rules? Latitudinal gradients in species richness, body size, and geographic range area of the worldʼs turtles. Journal of Experimental Zoology, Part B: Molecular and Developmental Evolution 324(3):270-294.
• Ashton, K. G. & Feldman, C. R. 2003. Bergmann’s rule in non avian reptiles: turtles follow it, lizards and snakes reverse it. Evolution 57(5):1151-1163.
• Ashton, D. T.; Bettaso, J. B. & Welsh Jr, H. H. 2015. 1across a decade in size, growth, and body condition of western pond turtle (Actinemys marmorata) populations on free-flowing and regulated forks of the Trinity River in northwest California. Copeia 103(3):621-633.
• Ataídes, A. G.; Malvasio, A. & Parente, T. G. 2010. Percepções sobre o consumo de quelônios no entorno do Parque Nacional do Araguaia, Tocantins: conhecimentos para conservação. Gaia Scientia 4(1):7-20.
• Balestra, R. A. M.; Valadão, R. M.; Vogt, R. C.; Bernhard, R.; Ferrara, C. R.; Brito, E. S.; Arias, R. B.; Malvasio, A.; Lustosa, A. P. G.; Souza, F. L.; Drummond, G. M.; Bassetti, L. A. B.; Coutinho, M. E.; Ferreira-Júnior, P. D.; Campos, Z. M. S.; Mendonça, S. H. S. T.; Rocha, J. M. N. & Luz, V. L. F. 2016. Roteiro para Inventários e Monitoramentos de Quelônios Continentais. Biodiversidade Brasileira 6(1):114-152.
• Bennett, A. M. & Litzgus, J. D. 2014. Injury rates of freshwater turtles on a recreational waterway in Ontario, Canada. Journal of Herpetology 48(2):262-266.
• Bernardes, V. C. D.; Ferrara, C. R.; Vogt, R. C. & Schneider, L. 2014. Abundance and population structure of Podocnemis erythrocephala (Testudines, Podocnemididae) in the Unini River, Amazonas. Chelonian Conservation and Biology 13(1):89-95.
• Bernhard, R. & Vogt, R. C. 2012. Population structure of the turtle Podocnemis erythrocephala in the Rio Negro Basin, Brazil. Herpetologica 68(4):491-504.
• Biaggini, M. & Corti, C. 2018. Facing habitat reduction in your own shell: patterns of non-lethal injuries in the endangered tortoise Testudo Hermanni in Italy. Herpetological Conservation and Biology 1(3):539-550.
• Bonach, K.; Malvasio, A.; Matushima, E. R. & Verdade, L. M. 2011. Temperature-sex determination in Podocnemis expansa (Testudines, Podocnemididae). Iheringia, Série Zoologia 101(3):151-155.
• Brasil. Ministério do Meio Ambiente. 2019. Sítios Ramsar do Brasil. Disponível em <Disponível em https://www.mma.gov.br/areas-protegidas/instrumentos-de-gestao/s%C3%ADtios-ramsar.html >. Acesso em 03.12.2019.
  » https://www.mma.gov.br/areas-protegidas/instrumentos-de-gestao/s%C3%ADtios-ramsar.html
• Brito, E. R.; Martins, S. V.; Oliveira-Filho, A. T. D.; Silva, E. & Silva, A. F. D. 2006. Estrutura fitossociológia de um fragmento natural de floresta inundável em área de orizicultura irrigada, município de Lagoa da Confusão, Tocantins. Revista Árvore 30(5):829-836.
• Burnham, K. P. & Anderson, D. R. 2002. Model selection and multimodel inference: a practical information-theoretic approach. 2ed. New York, Springer-Verlag. 488p.
• Bunn, S. E. & Arthington, A. H. 2002. Basic principles and ecological consequences of altered flow regimes for aquatic biodiversity. Environmental Management 30(4):492-507.
• Cagle, F. R. 1954. Observation on the life cycles of painted turtles (Genus Chrysemys). The American Midland Naturalist 52:225-235.
• Carneiro-Neto, M. B.; Mendonça, K. L. & Mendonça, J. C. 2017. Avaliação da Implantação de uma Barragem Articulada no Rio Formoso na Cidade de Lagoa da Confusão - TO. Revista Integralização Universitária 11(17):181-200.
• Castello, L. & Macedo, M. N. 2016. Large-scale degradation of Amazonian freshwater ecosystems. Global Change Biology 22(3):990-1007.
• Cecala, K. K.; Gibbons, J. W. & Dorcas, M. E. 2009. Ecological effects of major injuries in diamondback terrapins: implications for conservation and management. Aquatic Conservation: Marine and Freshwater Ecosystems 19(4):421-427.
• Chessman, B. C. 2011. Declines of freshwater turtles associated with climatic drying in Australia’s Murray-Darling Basin. Wildlife Research 38(8):664-671.
• CITES - Convention on International Trade in Endangered Species of Wild Fauna and Flora. 2017. Appendices I, II and III. Disponível em <Disponível em https://cites.org/eng >. Acesso em 01.01.2017.
  » https://cites.org/eng
• Crawley, M. J. 2013.The R Book. 2ed. Chichester, John Wiley & Sons. 975p.
• Crouse, D. & Frazer, N. 1995. Population models and structure. In: Bjorndal, K. ed. Biology and conservation of sea turtles. Washington D.C., Smithsonian Institution Press, p. 601-603.
• Döll, P. & Zhang, J. 2010. Impact of climate change on freshwater ecosystems: a global-scale analysis of ecologically relevant river flow alterations. Hydrology & Earth System Sciences Discussions 7(1):1305-1342.
• Downes, S. & Shine, R. 2001. Why does tail loss increase a lizardʼs later vulnerability to snake predators? Ecology 82(5):1293-1303.
• Eisemberg, C. C.; Balestra, R. A. M.; Famelli, S.; Pereira, F. F.; Bernardes, V. C. D. & Vogt, R. C. 2016. Vulnerability of Giant South American Turtle (Podocnemis expansa) nesting habitat to climate-change-induced alterations to fluvial cycles. Tropical Conservation Science 9(4):1-12.
• Fachín-Terán, A. & Vogt, R. C. 2004. Estrutura populacional, tamanho e razão sexual de Podocnemis unifilis (Testudines, Podocnemididae) no rio Guaporé (RO), norte do Brasil. Phyllomedusa 3(1):29-42.
• Fachín-Terán, A.; Vogt, R. C. & Thorbjarnarson, J. B. 2003. Estrutura populacional, razão sexual e abundância de Podocnemis sextuberculata (Testudines, Podocnemididae) na reserva de desenvolvimento sustentável Mamirauá, Amazonas, Brasil. Phyllomedusa 2(1):43-63.
• Faria, V. A. & Malvasio, A. 2018. Aspectos sobre a caça, comercialização e consumo de quelônios na região do Corredor Ecológico Araguaia Bananal no Estado do Tocantins. Revista Ouricuri 8(2):26-48.
• Ferrara, C. R.; Fagundes, C. K.; Morcatty, T. Q. & Vogt, R. C. 2017. Quelônios Amazônicos Guia de identificação e distribuição. Manaus, Wildlife Conservation Society Brasil.182p.
• Figueroa, I. C.; Fachín-Terán, A. & Duque, S. R. 2013. Habitat, population structure and community consume of Podocnemis unifilis and P. expansa (Testudines: Podocnemididae) turtles in Curare-Los Ingleses reserve, La Pedrera, Amazonas, Colombia. Mundo Amazónico 4:153-173.
• Fraxe-Neto, H. J.; Brasil, M. A.; Horta, G. F.; Barros, T. O.; Falcon, G. B. & Colli, G. R. 2011. Demografia de Acanthochelys spixii (Testudines, Chelidae) no Cerrado brasileiro. Demography of Acanthochelys spixii (Testudines, Chelidae) in the Brazilian Cerrado. Chelonian Conservation and Biology 10(1):82-90.
• Fuller, W. J.; Godley, B. J.; Hodgson, D. J.; Reece, S. E.; Witt, M. J. & Broderick, A. C. 2013. Importance of spatio-temporal data for predicting the effects of climate change on marine turtle sex ratios. Marine Ecology Progress Series 488:267-274.
• Futuyma, D. J. 1993. Biologia Evolutiva. 2ed. Ribeirão Preto, SGB. 620p.
• Galinkin, M.; Dias, A.; Latrubesse, E. M.; Scardua, F. P.; Mendonça, A. F. & Arruda, M. B. 2003. Projeto corredor ecológico Araguaia - Bananal. In: Arruda, M. B. & Sá, L. F. S. N. orgs. Corredores ecológicos: uma abordagem integradora de ecossistemas no Brasil. Brasília, IBAMA, p. 81-132.
• Gibbons, J. W. 1990. Sex ratios and their significance among turtle populations. In: Gibbons, J. W. ed. Life History and Ecology of the Slider Turtle. Washington D.C., Smithsonian Institution Press, p. 171-182.
• Gibbs, J. & Amato, G. 2000. Genetics and demography in turtle conservation. In: Klemens, M. W. ed. 2000. Turtle Conservation. Washington D.C., Smithsonian Institution Press, p. 207-217.
• Glen, F. & Mrosovsky, N. 2004. Antigua revisited: the impact of climate change on sand and nest temperatures at a hawksbill turtle (Eretmochelys imbricata) nesting beach. Global Change Biology 10(12):2036-2045.
• Greaves, W. F. & Litzgus, J. D. 2009. Variation in life history characteristics among populations of North American wood turtles: a view from the north. Journal of Zoology 279(3):298-309.
• Guimberteau, M.; Ronchail, J.; Espinoza, J. C.; Lengaigne, M.; Sultan, B.; Polcher, J.; Drapeau, G.; Guyot, J. L.; Ducharne, A. & Ciais, P. 2013. Future changes in precipitation and impacts on extreme stream flow over Amazonian sub-basins. Environmental Research Letters 8(1):014035.
• Gusmão, A. C. V. L.; Silva, B. D.; Montenegro, S. M. G. L.; Galvíncio, J. D. & Oliveira, L. D. 2013. Índice de vegetação e temperatura da superfície no ecótono Ilha do Bananal por sensoriamento remoto. Revista de Geografia 30(3):209-225.
• wkes, L. A.; Broderick, A. C.; Godfrey, M. H. & Godley, B. J. 2007. Investigating the potential impacts of climate change on a marine turtle population. Global Change Biology 13(5):923-932.
• Heinze, G.; Wallisch, C. & Dunkler, D. 2018. Variable selection - A review and recommendations for the practicing statistician. Biometrical Journal 60(3):431-449.
• Henen, B. T. 1997. Seasonal and annual energy budgets of female desert tortoises (Gopherus agassizii). Ecology 78(1):283-296.
• Hernández, O. & Espín, R. 2006. Efectos del reforzamiento sobre la población de tortuga Arrau (Podocnemis expansa) enel Orinoco medio, Venezuela. Interciencia 31(6):424-430.
• Howard, K.; Beesley, L.; Ward, K. & Stokeld, D. 2017. Preliminary evidence suggests freshwater turtles respond positively to an environmental water delivery during drought. Australian Journal of Zoology 64(5):370-373.
• IBGE - Instituto Brasileiro de Geografia e Estatística. 2019. Cidades. Disponível em <Disponível em https://cidades.ibge.gov.br/brasil/panorama >. Acesso em 12.01.2019.
  » https://cidades.ibge.gov.br/brasil/panorama
• ICMBio - Instituto Chico Mendes de Conservação da Biodiversidade. 2018. Livro Vermelho da Fauna Brasileira Ameaçada de Extinção. Brasília, ICMBio. 492p.
• IUCN- International Union for Conservation of Nature. 2016. The IUCN Red List of Threatened Species. Version 2016-3. Disponível em <Disponível em http://www.iucnredlist.org >. Acesso em 01.01.2017.
  » http://www.iucnredlist.org
• Janzen, F. J. 1994. Climate change and temperature-dependent sex determination in reptiles. Proceedings of the National Academy of Sciences 91:7487-7490.
• Jensen, M. P.; Allen, C. D.; Eguchi, T.; Bell, I. P.; LaCasella, E. L.; Hilton, W. A.; Hof, C. A. M. & Dutton, P. H. 2018. Environmental warming and feminization of one of the largest sea turtle populations in the world. Current Biology 28(1):154-159.
• Jessop, T. S.; Hamann, M. & Limpus, C. J. 2004. Body condition and physiological changes in male green turtles during breeding. Marine Ecology Progress Series 276:281-288.
• Koenker, R. 2017. Quantile regression: 40 years on. Annual Review of Economics 9:155-176.
• Koenker, R. 2020. quantreg: Quantile Regression. R package version 5.61. Disponível em <https://CRAN.R-project.org/package=quantreg>.
  » https://CRAN.R-project.org/package=quantreg
• Koenker, R. & Bassett-Júnior, G. 1978. Regression quantiles. Econometrica: Journal of the Econometric Society 46(1):33-50.
• Koenker, R. & Hallock, K. F. 2001. Quantile regression. Journal of Economic perspectives 15(4):143-156.
• Laloë, J. O.; Esteban, N.; Berkel, J. & Hays, G. C. 2016. Sand temperatures for nesting sea turtles in the Caribbean: Implications for hatchling sex ratios in the face of climate change. Journal of Experimental Marine Biology and Ecology 474:92-99.
• Lima, R. A. M.; Parente, T. G. & Malvasio, A. 2013. Concepções dos trabalhadores sobre o uso de agrotóxicos em projetos de irrigação, Lagoa da Confusão - TO, Brasil. Gaia Scientia 7(1):31-41.
• Limpus, C. 2008. Freshwater Turtles in the Mary River: Review of biological data for turtles in the Mary River, with emphasis on Elusor macrurus and Elseya albagula. Brisbane, Queensland Government. 84p.
• Litzgus, J. D. & Mousseau, T. A. 2004. Demography of a southern population of the spotted turtle (Clemmys guttata). Southeastern Naturalist 3(3):391-400.
• Litzgus, J. D.; Bolton, F. & Schulte-Hostedde, A. I. 2008. Reproductive output depends on body condition in spotted turtles (Clemmys guttata). Copeia (1):86-92.
• Lopes, M. H.; Franco, J. L. A. & Costa, K. S. 2017. Expressões da Natureza no Parque Nacional do Araguaia: Processos geoecológicos e diversidade da vida. Historia Ambiental Latinoamericana y Caribeña (HALAC) revista de la Solcha 7(2):65-100.
• Luz, V. L. F.; Malvasio, A.; Balestra, R. A. M.; Salera-Júnior, G.; Souza, V. L.; Portelinha, T. C. G.; Uhlig, V. M. & Portal, R. R. 2019. Conservação dos quelônios amazônicos no Brasil. In: Lacava, R. V. & Balestra, R. A. M. orgs. Plano de Ação Nacional para a Conservação dos Quelônios Amazônicos. Brasília, IBAMA, p.11-18.
• Malvasio, A.; Nascimento-Rocha, J. M.; Santos, H. D.; Ataídes, A. G. & Portelinha, T. C. G. 2012. Morfometria e histologia das gônadas de machos e fêmeas recém-eclodidos de Podocnemis expansa e Podocnemis unifilis (Testudines, Podocnemididae. Acta Scientiarum, Biological Sciences 34(1):105-112.
• Malvasio, A.; Souza, A. M.; Ferreira-Júnior, P. D.; Reis, E. S. & Sampaio, F. A. A. 2002. Temperatura de incubação dos ovos e granulometria dos sedimentos das covas relacionadas à determinação sexual em Podocnemis expansa (Schweigger, 1812) e P. unifilis (Troschel, 1848) (Testudines, Pelomedusidae). Publicações Avulsas do Instituto Pau Brasil de História Natural 5(1):11-25.
• Malvasio, A.; Sousa-Segundo, J. P. B.; Salera-Júnior, G.; Ataídes, A. G.; Montelo, K. M.; Karajá, A. I. D. A.; Lopes, T. K. M.; Prado, T. R. L. & Portelinha, T. C. G. 2019. Biologia Populacional e Reprodutiva de Podocnemis expansa e Podocnemis unifilis no Parque Nacional do Araguaia, Tocantins: Um Histórico Entre os Anos de 2004 e 2012. In: Pinheiro, R. T. org. Biodiversidade na região da ilha do Bananal/Cantão. Palmas, EDUFT, p. 77-102.
• Marchand, M. N. & Litvaitis, J. A. 2004. Effects of habitat features and landscape composition on the population structure of a common aquatic turtle in a region undergoing rapid development. Conservation Biology 18(3):758-767.
• Marinho-Filho, G. M.; Maciel, G. F.; Dias, R. R.; Magalhães-Filho, L. L.; Rezende, C. D. S. A.; Figueroa, F. E. V. & Oliveira, L. M. 2013. Avaliação de Características Morfométricas da Bacia Hidrográfica do Rio Formoso - TO. Revista de Ciências Ambientais 7(1):37-48.
• McCoy, E. D.; Moore, R. D.; Mushinsky, H. R. & Popa, S. C. 2011. Effects of rainfall and the potential influence of climate change on two congeneric tortoise species. Chelonian Conservation and Biology 10(1):34-41.
• McCullagh, P. & Nelder, J. A. 1989. Generalized Linear Models. London, Chapman and Hall. 511p.
• Meek, R. 2007. Non-lethal injury in Hermannʼs tortoise, Testudo hermanni, in Croatia and Montenegro. Herpetological Bulletin 100:23-29.
• Miorando, P. S.; Giarrizzo, T. & Pezzuti, J. C. 2015. Population structure and allometry of Podocnemis unifilis (Testudines, Podocnemididae) in a protected area upstream Belo Monte dam in Xingu River, Brazil. Anais da Academia Brasileira de Ciências 87(4):2067-2079.
• Mitchell, N. J. & Janzen, F. J. 2010. Temperature-dependent sex determination and contemporary climate change. Sexual Development 4(1-2):129-140.
• Mogollones, S. C.; Rodríguez, D. J.; Hernández, O. & Barreto, G. R. 2010. A demographic study of the arrau turtle (Podocnemis expansa) in the Middle Orinoco River, Venezuela. Chelonian Conservation and Biology 9(1):79-89.
• Morais, P. B.; Nolêto-Júnior, S. & Martins, I. C. M. 2017. Análise de Sustentabilidade do Projeto Hidroagrícola Javaés/Lagoa, no Estado do Tocantins. Cadernos de Ciência &Tecnologia 34(1):83-111.
• Peig, J. & Green, A. J. 2009. New perspectives for estimating body condition from mass/length data: the scaled mass index as an alternative method. Oikos 118(12):1883-1891.
• Peig, J. & Green, A. J. 2010. The paradigm of body condition: a critical reappraisal of current methods based on mass and length. Functional Ecology 24(6):1323-1332.
• Peñaloza, C. L.; Hernández, O.; Espín, R.; Crowder, L. B. & Barreto, G. R. 2013. Harvest of endangered side neck river turtles (Podocnemis spp.) in the middle Orinoco, Venezuela. Copeia 2013(1):111-120.
• Polo-Cavia, N.; Engstrom, T.; López, P. & Martín, J. 2010. Body condition does not predict immunocompetence of western pond turtles in altered versus natural habitats. Animal Conservation13(3):256-264.
• Portelinha, T. C.; Malvasio, A.; Piña, C. I. & Bertoluci, J. 2014. Population Structure of Podocnemis expansa (Testudines: Podocnemididae) in Southern Brazilian Amazon. Copeia4:707-715.
• Pritchard, P. C. H. & Trebbau, P. 1984. The turtles of Venezuela. Ohio, Society for the Study of Amphibians and Reptiles. 414p.
• R Core Team. 2020. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria.
• Rhodin, A. G. J.; Iverson, J. B.; Bour, R.: Fritz, U.; Georges, A.; Shaffer, H. B. & van Dijk, P. P. 2017. Turtles of the World: Annotated Checklist and Atlas of Taxonomy, Synonymy, Distribution, and Conservation Status. 8ed. In: Rhodin, A. G. J.; Iverson, J. B.; van Dijk, P. P.; Saumure, R. A.; Buhlmann, K. A.; Pritchard, P. C. H. & Mittermeier, R. A. eds. Conservation Biology of Freshwater Turtles and Tortoises: A Compilation Project of the IUCN/SSC Tortoise and Freshwater Turtle Specialist Group. Chelonian Research Monographs 7:1-292.
• Rodrigues, M. T. 2005. The Conservation of Brazilian Reptiles: Challenges for a Megadiverse Country. Conservation Biology 19(3):659-664.
• Salera-Júnior, G.; Balestra, R. A. M. & Luz, V. L. F. 2016. Breve histórico da conservação dos quelônios amazônicos no Brasil. In: Balestra, R. A. M. org. Manejo conservacionista e monitoramento populacional de quelônios amazônicos. Brasília, IBAMA , p.11-14.
• Saumure, R. A. & Bider, J. R. 1998. Impact of agricultural development on a population of wood turtles (Clemmys insculpta) in southern Quebec, Canada. Chelonian Conservation and Biology 3:37-45.
• Schneider, A. C.; Arnold, T. W.; Huber, P. W. & Lewis, T. L. 2018. An 18-year mark-recapture study of wood turtles (Glyptemys insculpta) in Michigan. Journal of Herpetology 52(2):193-200.
• Schulte-Hostedde, A. I.; Millar, J. S. & Hickling, G. J. 2001. Evaluating body condition in small mammals. Canadian Journal of Zoology 79:1021-1029.
• Schulte-Hostedde, A. I.; Zinner, B.; Millar, J. S. & Hickling, G. J. 2005. Restitution of mass-size residuals: validating body condition indices. Ecology 86:155-163.
• Shine, R. & Iverson, J. B. 1995. Patterns of survival, growth and maturation in turtles. Oikos 72:343-348.
• Simoncini, M. S.; Portelinha, T. C. G.; Montelo, K. M.; Falcon, G. B.; Collicchio, E.; Balestra, R. A. M.; Luz, V. L. F. & Malvasio, A. 2019. Avaliação da influência dos fatores ambientais sazonais na reprodução da tartaruga-da-amazônia (Podocnemis expansa): um estudo de caso no Tocantins. In: Lacava, R. V. & Balestra, R. A. M. orgs. 2019. Plano de Ação Nacional para a Conservação dos Quelônios Amazônicos. Brasília, IBAMA , p. 89-99.
• Smith, L. D. 1995. Effects of limb autotomy and tethering on juvenile blue crab survival from cannibalism. Marine Ecology Progress Series 116(1):65-74.
• Souza, F. L. 2004. Uma revisão sobre padrões de atividade, reprodução e alimentação de cágados brasileiros (Testudines, Chelidae). Phyllomedusa 3:15-27.
• Stanford, C. B.; Iverson, J. B.; Rhodin, A. G.; van Dijk, P. P.; Mittermeier, R. A.; Kuchling, G.; Berry, K. H.; Bertolero, A.; Bjorndal, K. A.; Blanck, T. E. G.; Buhlmann, K. A.; Burke, R. L.; Congdon, J. D.; Diagne, T.; Edwards, T.; Eisemberg, C. C.; Ennen, J. R.; Forero-Medina, G.; Frankel, M.; Fritz, U.; Gallego-García, N.; Georges, A.; Gibbons, J. W.; Gong, S.; Goode, E. V.; Shi, H. T.; Hoang, H.; Hofmeyr, M. D.; Horne, B. D.; Hudson, R.; Juvik, J. O.; Kiester, R. A.; Koval, P.; Le, M.; Lindeman, P. V.; Lovich, J. E.; Luiselli, L.; McCormack, T. E.M.; Meyer, G. A.; Páez, V. P.; Platt, K.; Platt, S. G.; Pritchard, P. C.H.; Quinn, H. R.; Roosenburg, W. M.; Seminoff, J. A.; Shaffer, H. B.; Spencer, R.; Van Dyke, J. U.; Vogt, R. C. & Walde, A. D. 2020. Turtles and tortoises are in trouble. Current Biology 30(12):R721-R735.
• Steen, D. A.; Aresco, M. J.; Beilke, S. G.; Compton, B. W.; Condon, E. P.; Dodd Jr, C. K.; Forrester, H.; Gibbons, J. W.; Greene, J. L.; Johnson, G.; Langen, T. A.; Oldham, M. J.; Oxier, D. N.; Saumure, R. A.; Schueler, F. W.; Sleeman, J. M.; Smith, L. L.; Tucker, J. K.; Gibbs, J. P. 2006. Relative vulnerability of female turtles to road mortality. Animal Conservation 9(3):269-273.
• Sung, Y. H.; Karraker, N. E. & Hau, B. C. 2013. Demographic evidence of illegal harvesting of an endangered Asian turtle. Conservation Biology 27(6):1421-1428.
• Taylor, P. W. & Jackson, R. R. 2003. Interacting effects of size and prior injury in jumping spider conflicts. Animal Behaviour 65(4):787-794.
• Thomson, J. A.; Burkholder, D.; Heithaus, M. R. & Dill, L. M. 2009. Validation of a rapid visual-assessment technique for categorizing the body condition of green turtles (Chelonia mydas) in the field. Copeia 2009(2):251-255.
• Tocantins - Governo do Estado do Tocantins. 2014. Portarias números 318, 319, 320 e 321, de 23 de julho de 2014. Diário Oficial do Estado do Tocantins, n. 4.176, p. 24-25. Disponível em <Disponível em https://diariooficial.to.gov.br/busca/?por=edicao&edicao=4176 >. Acesso em 02.01.2020.
  » https://diariooficial.to.gov.br/busca/?por=edicao&edicao=4176
• Tocantins - Governo do Estado do Tocantins. 2018. Extrato da Licença de Operação número 2855-2018. Diário Oficial do Estado do Tocantins, n. 5.165, p. 56. Disponível em <Disponível em https://diariooficial.to.gov.br/busca/?por=edicao&edicao=5165 >. Acesso em 02.01.2020.
  » https://diariooficial.to.gov.br/busca/?por=edicao&edicao=5165
• Tomillo, P. S.; Saba, V. S.; Piedra, R.; Paladino, F. V. & Spotila, J. R. 2008. Effects of illegal harvest of eggs on the population decline of leatherback turtles in Las Baulas Marine National Park, Costa Rica. Conservation Biology 22(5):1216-1224.
• Usuda, H.; Morita, T. & Hasegawa, M. 2012. Impacts of river alteration for flood control on freshwater turtle populations. Landscape and Ecological Engineering 8(1):9-16.
• Valenzuela, N. 2001. Constant, shift, and natural temperature effects on sex determination in Podocnemis expansa turtles. Ecology 82(11):3010-3024.
• Valenzuela, N.; Botero, R. & Martínez, E. 1997. Field study of sex determination in Podocnemis expansa from Colombian Amazonia. Herpetologica 53(3):390-398.
• Venables, W. N. & Ripley, B. D. 2002. Modern Applied Statistics with S. 4ed. New York, Springer. 498p.
• Vergara, F. E.; Reis, F. C. & Magalhães, L. N. L. 2013. Proposta de vazão de referência Q 90 para o Rio Formoso na Bacia do Araguaia. Engenharia Ambiental 10(1):84-102.
• Vogt, R. C. 1994. Temperature Controlled Sex Determination as a Tool for Turtle Conservation. Chelonian Conservation and Biology 1(2):159-162.
• Werner, Y. L.; Korolker, N.; Sion, G. & Göçmen, B. 2016. Bergmannʼs and Renschʼs rules and the spur-thighed tortoise (Testudo graeca). Biological Journal of the Linnean Society 117(4):796-811.
• Zug, G. R. 1991. Age determination in turtles. Ohio, Society for the Study of Amphibians and Reptiles . 28p.
• Zuur, A. F.; Ieno, E. N.; Smith, G. M. 2007. Analysing ecological data. New York, Springer Science + Business Media. 672p.

Material suplementar.

O seguinte material online está disponível para este artigo:

Análise estatística; Condição corporal; Número de indivíduos capturados e abundância.

Datas de Publicação

• Publicação nesta coleção
  18 Out 2021
• Data do Fascículo
  2021

Histórico

• Recebido
  27 Abr 2020
• Aceito
  04 Maio 2021