1A morfometria da rede de drenagem sempre chamou a atenção da comunidade científica brasileira, com muitos estudos nas mais variadas escalas de análise, sendo, até certo ponto, algo superado e sem grandes novidades desde a 1ª edição da obra “Geomorfologia” de Christofoletti (1980), amplamente divulgada e utilizada até os dias atuais. Christofoletti (1980) foi utilizado em demasia ao longo do tempo, sobretudo por sua capacidade de expor, em apenas uma obra, um compilado de metodologias elaboradas na década de 1940 e 1950, sobre vertentes, bacias hidrográficas e suas drenagens, relevo, entre outras temáticas. Obras estas, que devido à época, são de difícil acesso. Horton (1945), Strahler (1952) e tantos outros autores mencionados trouxeram inovações teórico-metodológicas, e que fizeram e ainda fazem parte do construto da ciência geomorfológica e, consequentemente, da Geomorfologia Fluvial.

2Desde então, muitos seguidores elaboraram trabalhos seja sobre a ótica tectônica, hidrológica, de movimentos de massa ou processos antrópicos, contudo, muitos são os trabalhos, mais pautados na aplicação da metodologia, do que propriamente na inovação dos conceitos, métodos e procedimentos metodológicos. Barros e Reis (2019) e Reis et al. (2023) destacam que no Brasil, os estudos para investigar padrões morfológicos e morfométricos do relevo e da rede de drenagem, quase sempre tem como unidades de estudo, as bacias hidrográficas. Sua análise partindo do princípio sistêmico é fundamental para atrelar tais estudos à perspectiva do planejamento e da gestão ambiental.

3Diante da exigência de uma obra que suprisse a necessidade de atualização das bases do Christofoletti (1980), Stevaux e Latrubesse (2017) objetivaram, ante a experiência adquirida por tais pesquisadores, compreender os processos fluviais e ampliaram a discussão sobre a análise morfométrica, que carecia de novos capítulos. Tal obra apresenta alguns parâmetros distintos e outras formas de análise. Mesmo com tal avanço, é importante destacar nesse contexto, que devido à pouca abrangência em território nacional, os sistemas cársticos têm sido mais estudados no contexto da espeleologia e geologia, e pouco se têm na literatura uma discussão sobre a análise morfométrica em sistemas cársticos. Aqui se tratando da bacia hidrográfica como a do rio Formoso, em que apresenta sistemas cárstico e não-cárstico (nessa pesquisa chamado de terrígeno), o que é ainda mais difícil de se deparar. Fatos que trazem consigo uma heterogeneidade de drenagens, de formas, de densidades e de características superficiais e subterrâneas.

4Antes de adentrar a essa questão, torna-se oportuno entender que a morfometria é um tipo de análise que abrange os componentes físicos associados ao relevo e à rede hidrográfica, de modo que influenciarão na velocidade e no tempo do escoamento das águas, além da variação na infiltração da água no solo (Leal e Tonello, 2016; Khanday e Javed, 2017; Stevaux e Latrubesse, 2017). Portanto, avaliar as características físicas de uma bacia hidrográfica, seja por meio de quantificações e/ou espacialização (mapeamento) de dados, fornecem subsídios para um melhor planejamento e gestão dos recursos hídricos.

5Em ambientes cársticos, os aspectos relacionados ao relevo são igualmente importantes, contudo, o substrato rochoso exerce influência direta sobre os demais componentes, seja pelo grau de coesão das rochas, resistência dos minerais e evolução geológica, o que determina uma influência sobre as águas superficiais e subterrâneas. Principalmente, pelos processos de dissolução química das rochas, influenciando na dinâmica, distribuição, geometria, forma, magnitude e frequência dos fluxos fluviais.

6Ao compreender Fritzsons e Mantovani (2010), torna-se possível afirmar que em sistemas cársticos como Bonito, no Mato Grosso do Sul, isto é, localizado em clima úmido e semiúmido, o substrato geológico torna-se fundamental, pois a precipitação exerce ação direta nos processos morfodinâmicos, moldando os canais de drenagem e interferindo na sua distribuição superficial e subterrânea. As rochas carbonáticas são permeáveis, porosas e solúveis, propensas à dissolução química, estando sujeitas a processos de abaulamento e com alta influências hidrogeológicas.

7Ao analisar a literatura brasileira, pouco se tem sobre recursos hídricos em ambiente cárstico, contudo, em estudos importantes pelo mundo, Burri, Castiglioni e Sauro (1999), Piló (2000), Ford e Williams (2007), Gutiérrez et al. (2014) e Parise (2012) discutem que em sistemas cársticos as drenagens sempre se vinculam ao subterrâneo, na chamada drenagem criptorrêica, tendo uma ausência marcante de drenagens superficiais. Apontando fragilidades e potencialidades disso, se por um lado as drenagens subterrâneas tem maior capacidade de abastecer as superficiais em épocas de estiagem, mantendo uma vazão média por longos períodos, por outro lado, uma intensidade de uso das terras, como a agricultura, poderá interferir de forma significativa no lençol freático e aquífero cárstico, diante da porosidade e capacidade de dissolução das rochas, o que favorece a infiltração e a contaminação hídrica.

8Essa dinâmica hidrogeológica e as informações que podem ser obtidas e demonstradas por meio de representações cartográficas, interpretam o relevo e as formas do modelado cárstico, e detalham sua evolução e sucessão espacial e temporal. Tais representações do relevo ou documentos morfométricos são importantes na compreensão dos processos que modificam as bacias hidrográficas.

9Atualmente, são inúmeros os trabalhos existentes sobre a análise morfométrica no estado de Mato Grosso do Sul, porém, a grande maioria está vinculado a uma bacia hidrográfica genuinamente terrígena (Berezuk et al., 2014; Miguel et al., 2014; Pirajá e Rezende Filho, 2019; Jesus, Berezuk e Brugnolli, 2020); poucos são os estudos em sistemas cársticos (Brugnolli, Pinto e Alves, 2020); mostrando a ausência de trabalhos em bacias hidrográficas com ambientes distintos, o que possibilita verificar tais distinções e acrescentar novos capítulos às discussões sobre esses sistemas, avaliando-os do ponto de vista hídrico, geológico e geomorfológico, além do interesse na temática do diagnóstico ambiental.

10De tal modo, para trabalhar com a Bacia Hidrográfica do Rio Formoso (BHRF) é preciso contextualizar que a mesma apresenta: sistemas cárstico e terrígeno. O que traz uma dinâmica distinta para a bacia e para os recursos hídricos em particular, muitos desses ligados ao turismo de natureza. Há, ainda, a existência de Unidades de Conservação (UC’s), como os Monumentos Naturais (MNE’s) do Rio Formoso e da Gruta do Lago Azul, as Reservas Particulares do Patrimônio Natural (RPPN’s) da Fazenda São Geraldo e de São Pedro da Barra, além do Parque Nacional (PARNA) da Serra da Bodoquena, todos estes auxiliam na preservação e manutenção da fauna e flora, pois impõem restrições e mantêm a integridade das UC’s; bem como um intenso uso antrópico das terras, especialmente a soja e pastagem, que vêm impactando nos recursos hídricos da BHRF.

11Diante dessas justificativas principais, essa pesquisa objetiva avaliar a morfometria da rede de drenagem da BHRF, abordando todas as sub-bacias existentes na BHRF, tanto em termos areais quanto lineares, de modo a oferecer subsídios para uma gestão dos recursos hídricos cársticos e terrígenos, buscando a manutenção do equilíbrio, quantidade e qualidade das águas superficiais e, por consequência, as subterrâneas.

12No contexto geográfico, a BHRF (Figura 1) está situada a sudoeste do estado de Mato Grosso do Sul, sua totalidade em Bonito, um município reconhecido pelo turismo de natureza ligado às águas superficiais e pelas várias UC’s com vegetações nativas (Mata Atlântica e Cerrado) relativamente preservadas. O rio Formoso tem suas nascentes na Serra da Bodoquena e apresenta vários afluentes de grande relevância, como os rios Sucuri e Mimoso e os córregos Formosinho, Anhuma, Taquaral e tantos outros.

Procedimentos Metodológicos

13A metodologia consistiu, em um primeiro momento, no mapeamento e breve caracterização geológica e geomorfológica da BHRF. Com relação à geologia, optou-se pelo uso de bases da Companhia de Pesquisa de Recursos Minerais (CPRM) em que foi realizada uma atualização e detalhamento diante de algumas incongruências observadas nas fontes de dados (CPRM) frente à realidade terrestre. Tais ajustes ocorreram mediante imagens de satélite (Sentinel 2A/MSI, órbita 135 e ponto 101, datada de 15 de março de 2017), cartas topográficas (folhas: Aldeia Lalima - SF.21-X-A-V; Jardim - SF.21-X-C-II; Rio Perdido - SF.21-X-C-I; Vila Campão - SF.21-X-A-IV; disponibilizadas pela Diretoria de Serviço Geográfico do Brasil) e saídas de campo (dezembro de 2016; março e setembro de 2017 e dezembro de 2018) por meio da visualização de afloramentos e manuseio de rochas ao longo de estradas, margens de mananciais e em meio a pastagens e vegetações em toda a BHRF.

14Para a geomorfologia, optou-se pelo mapeamento de duas variáveis: hipsometria e declividade. Ambas tiveram como base o Modelo Digital de Terreno (MDT/SRTM) da National Aeronautics and Space Administration (NASA – Earth Explorer, 2000). Tal modelo passou por métodos de pré-processamento de modo a corrigir falhas, distorções e ruídos dessas imagens de radar, proporcionando confiabilidade no mapeamento. Com isso, foi possível identificar os índices hipsométricos (altitudes) e a declividade (inclinações do relevo). Esta última observando classificações do relevo, segundo o Sistema Brasileiro de Classificação de Solos (SIBCS) da Empresa Brasileira de Pesquisa Agropecuária (EMBRAPA, 2018).

15Em um segundo momento, avaliou-se a morfometria da rede de drenagem, alcançada mediante uma análise de cada sub-bacia dos principais recursos hídricos da BHRF. Tal compartimentação trouxe uma ideia clara a respeito dos mananciais e a forma com que eles se caracterizam em cada sub-bacia analisada, embasados na metodologia de Christofoletti (1980) e Stevaux e Latrubesse (2017). De tal forma, ainda foi adicionada na interpretação a abrangência dos ambientes terrígeno e cárstico em cada sub-bacia, pois esse fator traz uma dinâmica diferenciada no que tange aos recursos hídricos superficiais.

16As variáveis embasadas na metodologia foram: a) Hierarquia Fluvial: abrange a ordem dos canais, de acordo com Horton (1945); b) Análise Linear: envolve a relação de bifurcação, comprimento médio dos segmentos fluviais, relação entre o comprimento médio dos canais de cada ordem, relação entre o índice médio do comprimento médio dos canais e o índice de bifurcação e, por fim, a extensão do percurso superficial; c) Análise Areal: abarca a forma da bacia com seu índice de circularidade, densidade hidrográfica, densidade de drenagem e o coeficiente de manutenção; conforme pode ser visualizado no Quadro 1 e Quadro 2.

Quadro 1. Características da sub-bacia hidrográfica, hierarquia fluvial e aspectos lineares da rede de drenagem.

Fator de Análise		Conceito e Equação
Caracte-rísticas da Bacia Hidro-gráfica	Área de drenagem (A)	A corresponde a toda a área drenada pelo conjunto do sistema fluvial, projetada em plano horizontal do divisor de água (Christofoletti, 1980).
	Perímetro (P)	Corresponde à medida/comprimento dos limites da bacia hidrográfica.
	Porcentagem do Ambiente Cárstico e Terrígeno	Chrsitofoletti (1980) afirma que dentro da hidrologia cárstica é notável a ausência de cursos superficiais. Avaliar a porcentagem desses sistemas em cada sub-bacia ajudará a analisar se tal afirmação se aplica na BHRF.
Hierarquia Fluvial	Ordem dos Canais	Consiste no processo de classificação do curso d'água (ou da área drenada que lhe pertence) no conjunto total da bacia hidrográfica na qual se encontra (Christofoletti, 1980). Nessa pesquisa aplicou-se o método de Horton (1945).
Análise Linear	Relação de Bifurcação	Relação entre o número total de segmentos de determinada ordem e o número total dos segmentos da ordem imediatamente superior. Estes valores indicam o grau de dissecação da bacia hidrográfica, quanto maior o valor do índice de bifurcação maior será o grau de dissecação, valores geralmente abaixo de 2 indicam relevo colinoso (Castro e Carvalho, 2009). Agrandir Original (png, 880 octets) Nu é o número de seguimentos de determinada ordem; Nu+1 é o número de seguimentos da ordem imediatamente superior.
	Comprimento médio dos segmentos fluviais	É a média aritmética entre os comprimentos de todos os rios da bacia (Brugnolli, Pinto e Alves, 2020).
	Relação entre o comprimento médio dos canais de cada ordem	Analisa a média aritmética entre os comprimentos dos rios de cada ordem e a quantidade de canais dessa ordem. Agrandir Original (png, 644 octets) Lm é Relação entre o Comprimento Médio dos Canais de cada ordem (Christofoletti, 1980); Lu é a soma total dos comprimentos dos canais de cada ordem; Nu é o número de seguimentos encontrados na respectiva ordem
	Relação entre o índice médio do comprimento médio dos canais e o índice de bifurcação	Analisa o comprimento médio dos canais e sua relação com o índice de bifurcação, ou seja, relação entre a extensão total e a quantidade de canais (Brugnolli, Pinto e Alves, 2020).
	Extensão do percurso superficial	Representa a distância média percorrida pelas enxurradas entre o interflúvio e o canal permanente. Durante a evolução do sistema de drenagem, a extensão do percurso superficial está ajustada ao tamanho apropriado relacionado com as bacias de primeira ordem (Christofoletti, 1980). Agrandir Original (png, 954 octets) Eps representa a extensão do percurso superficial; Dd é o valor da densidade de drenagem.

17Fonte: Elaborado pelos autores.

18Quadro 2. Aspectos areais da rede de drenagem e das sub-bacias hidrográficas.

Análise Areal	Forma da Bacia - Índice de Circularidade	Relaciona a forma da bacia com um retângulo, correlacionando a razão entre a largura média e o comprimento axial da bacia. Villela e Mattos (1975) citam que uma bacia com fator de forma baixo é menos suscetível a enchentes.
	Densidade Hidrográfica	Relação que se expressa entre o número de rios e a área da bacia. Importante para indicar a capacidade de gerar novos cursos de água na bacia hidrográfica em função das características pedológicas, geológicas e climáticas da área (Castro e Carvalho, 2009). Agrandir Original (png, 759 octets) Dh é a densidade de drenagem; N é o número total de rios; A é a área total da bacia hidrográfica.
	Densidade de Drenagem	É definida pela correlação entre o comprimento total dos canais e a área de drenagem. Quanto maior o índice menor é a capacidade de infiltrar água. (Castro e Carvalho, 2009). Agrandir Original (png, 746 octets) Dd é a densidade de drenagem; Lt é o comprimento total dos canais; A é a área total da bacia hidrográfica.
	Coeficiente de Manutenção	Relaciona a forma da bacia com um círculo, é um número adimensional que varia com a forma da bacia, independentemente de seu tamanho. De modo que, quanto maior a irregularidade da bacia maior será o coeficiente de compacidade (Santos et al., 2012). Agrandir Original (png, 1,1k) Cm é o coeficiente de manutenção; Dd é o valor da densidade de drenagem.

19Fonte: Elaborado pelos autores.

Resultados

Caracterização geológica e do relevo da BHRF

20Inserido na análise morfométrica da rede de drenagem, avaliar o contexto geológico e do relevo da BHRF transmite informações cruciais para entender e correlacionar parâmetros morfométricos à ocorrência e distribuição da rede de drenagem. Tais técnicas que permeiam os preceitos do geoprocessamento e mapeamento temático auxiliam no estudo dos ambientes cársticos e terrígenos de modo a compreender suas depressões (poljes), superfícies e pavimentos cársticos, morros e serras alongadas, além das planícies terrígenas e cársticas, que apresentam distinções entre suas características (Figura 2).

21A avaliação inicia por meio da apreensão do contexto geológico da região da Serra da Bodoquena, onde está localizada a BHRF, que se assenta sobre rochas carbonáticas do grupo Corumbá, relacionadas aos sistemas cársticos, que por sua vez, incidiram no desenvolvimento de depósitos sedimentares quaternários (Sallun Filho e Karmann, 2007).

22De modo específico para a BHRF, é possível entender que o seu alto e médio curso apresentam tais rochas carbonáticas, calcários calcíticos, oolíticos e dolomíticos, enquanto que o baixo curso exibe rochas do grupo Cuiabá, majoritariamente terrígenas, filitos, arenitos e xistos, bem como rochas das formações Aquidauana (arenito) e Rio Bacuri (xistos grafitosos).

23Essa diversidade traz à paisagem uma heterogeneidade, em que, por um lado estão feições cársticas, como dolinas, abismos, cavernas, poljes, sumidouros, (re)surgências, morros residuais e serras alongadas; por outro apresentam superfícies aplainadas a onduladas terrígenas. Até, por tais motivos, que a bacia como um todo, suas paisagens, estão diretamente relacionadas aos processos hidrogeoquímicos cársticos, que originaram uma formidável paisagem de contraste entre os patamares mais elevados, que podem alcançar até 720 metros na serra da Bodoquena, até patamares reduzidos, cerca de 200 metros na foz do rio Formoso.

24Os patamares mais elevados são cársticos e formados por rochas como o calcário, que devido à dissolução química, apresentam um contraste dos pavimentos cársticos como a serra (com vales em “v” que alcançam declives acima de 75%), até as áreas de poljes, com declives que não atingem mais que 20% (relevo considerado suavemente ondulado).

25No médio curso, em meio a grandes e extensas falhas, lineamentos estruturais e fraturas, ocorre a mudança dos sistemas cárstico para o terrígeno. É a partir dessa compreensão geológica que se começa a notar a influência do substrato rochoso na atividade morfodinâmica dos cursos fluviais superficiais. A mudança no relevo também é algo notável, as áreas de poljes dão lugar (sentido oeste-leste/alto-baixo curso) a superfícies aplainadas a suavemente onduladas originadas por rochas sedimentares do Grupo Cuiabá.

26Essa mudança de rochas e relevo associadas às suas características propiciaram, por meio das camadas intercaladas de calcários (frágil e poroso) e filitos e arenitos (resistentes e impermeáveis), o desenvolvimento de uma rica e diversificada rede de drenagem, exibindo um comportamento hidrológico único, que pode ser entendido e explicado por meio de uma série de fórmulas e valores que se vinculam à análise morfométrica.

Morfometria da rede de drenagem da BHRF

27A análise morfométrica procurou identificar as características dos recursos hídricos em cada sub-bacia hidrográfica mediante diversos índices (área de drenagem, perímetro, ordem dos canais, relação de bifurcação, comprimento médio dos segmentos fluviais, relação entre o comprimento médio dos canais de cada ordem, relação entre o índice médio do comprimento médio dos canais e o índice de bifurcação, extensão do percurso superficial, forma da bacia, índice de circularidade, densidade hidrográfica, densidade de drenagem e coeficiente de manutenção), (Tabela 1 e Figura 3).

28Todos esses índices foram expressos segundo as análises linear e areal que Christofoletti (1980) definiu como eficazes para a compreensão das particularidades da rede de drenagem das bacias hidrográficas, além de abranger alguns de seus aspectos físicos e como estes interferem nos mananciais, de modo a tornar a análise morfométrica mais aplicável ao planejamento e à gestão dos recursos hídricos e assim encontrar as respostas sobre a possível influência das rochas e relevo nas águas superficiais, não só na quantidade, mas em sua espacialização, densidade, etc.

29A análise morfométrica se iniciou pelo rio Sucuri (Figura 4), o primeiro afluente do rio Formoso e um dos mais reconhecidos pelo turismo de natureza, devido à translucidez de suas águas. A característica cárstica prevalece em sua totalidade na sub-bacia, evidenciando apenas um manancial em 6,69 km², com uma densidade de drenagem baixa (0,27 km/km²) se comparado com as demais bacias hidrográficas. Esse fato pode ser explicado por três motivos principais: o primeiro Christofoletti (1980) discute que é o relevo aplainado que reduz a ocorrência de drenagens; a segunda Christofoletti (1980) discute sobre a bacia hidrográfica estar situada em uma mesma zona climática; contudo, o terceiro e último fator que altera a densidade de drenagem são as rochas.

30No caso da BHRF, as rochas carbonatadas influenciam significativamente na densidade de drenagem e densidade hidrográfica, em praticamente todas as sub-bacias hidrográficas. Isso acontece pelo fato das rochas carbonatadas facilitarem a dissolução e percolação dos calcários por fendas e fissuras, portanto, seus mananciais correm, por muitas vezes, de forma subterrânea (criptorreica), com diversas ressurgências em outros locais. Isso ocorre, principalmente, em áreas como a serra da Bodoquena e nos morros residuais nas proximidades do médio curso do rio Mimoso.

31Christofoletti (1980) ressalta essa afirmação, discutindo que dentro da hidrologia cárstica é notável a ausência de cursos superficiais, podendo somente visualizar os traços deixados por antigos rios ou canais que desapareceram. A circulação da água no carste faz-se em pontos de absorção, que são pontos de fraqueza da massa rochosa, contudo, é muito comum que não possa se ver o caminho traçado, tampouco os locais de ressurgências, visto que o traçado subterrâneo é independente dos canais superficiais.

32Em Bonito, tais rochas são datadas de períodos ligados à idade Cambriana e Vendiana (630 milhões de anos), apresentam uma particularidade, que é a presença de fendas e fissuras, que resultam em drenagens criptorreicas, com a percolação das águas, formando bacias subterrâneas (Brugnolli e Salinas Chávez, 2021).

33De fato, tal afirmação pôde ser comprovada em todas as sub-bacias predominantemente cársticas. A sub-bacia do córrego Alegre não se difere, apresentando apenas o curso principal (1ª ordem) e como exibe a característica de rochas carbonatadas das formações serra da Bodoquena, Bocaina, Cerradinho e a subunidade Carbonática do grupo Cuiabá, além do relevo plano que reduz a quantidade de mananciais, há uma reduzida densidade de drenagem e densidade hidrográfica (0,18 km/km² e 0,01 cursos por km² respectivamente).

34Por outro lado, seu coeficiente de manutenção é o mais elevado da BHRF, ou seja, para se manter perene, a bacia necessita de 5.676,60 m² de área por metro de canal fluvial. Tal fato ocorre por ser inversamente proporcional à densidade de drenagem, pois, por apontar apenas um manancial e possuir uma área relativamente extensa, tais valores ficam elevados. Esse coeficiente é fundamental na análise morfométrica, tais dados constituem um parâmetro importante para se estimar o tempo necessário para que a água, precipitada nos limites da bacia, atinja o exutório ou a saída da mesma (Fritzsons e Mantovani, 2010).

35Mesmo estando relacionado às chuvas, percebe-se que o coeficiente de manutenção varia não conforme a precipitação na BHRF, mas sim conforme o carste. Pode-se notar que quanto maior a porcentagem de carste, maior é o coeficiente de manutenção, profundamente relacionada à baixa densidade de drenagem nas sub-bacias cársticas.

36A região de Bonito é controlada por massas tropicais e polares, ou seja, uma região úmida com inverno de temperaturas baixas, que ocorre uma atuação expressiva da precipitação no desenvolvimento do relevo cárstico, seja na esculturação da serra da Bodoquena ou até nas regiões de dissolução química do calcário, como as áreas de poljes e dolinas que se fazem presentes na BHRF (Zavatini, 1992). Entretanto, não apresenta grandes variações de precipitação dentro da própria bacia, o que inviabiliza avaliar a influência das chuvas no coeficiente de manutenção, ficando evidente apenas a influência do substrato rochoso.

37Partindo para a análise da sub-bacia do córrego Taquaral, amplamente discutida em Brugnolli, Pinto e Alves (2020), expõe 100% de suas terras sobre terrenos cársticos, além de ser um canal de 4ª ordem. Mesmo com uma ordem mais elevada, percebe-se que o córrego Taquaral é temporário, e que as águas variam de forma distinta ao longo dos anos analisados (2016 a 2020). Apesar dos períodos chuvosos serem na primavera e no verão, foi durante esse período que o córrego perdeu suas águas e encontrava-se seco. Por outro lado, era possível encontrar água em períodos secos, como outono e inverno (Figura 5). Fato que não ocorreu com a sub-bacia hidrográfica vizinha ao norte, a do córrego Anhuma. Essa irregularidade é muito bem explicada conforme Carvalho Júnior et al. (2008) que discute, de que as análises que relacionam a chuva com bacias hidrográficas cársticas são amplamente utilizadas, sobretudo, devido aos aquíferos cársticos que desempenham papel importante no regime das águas superficiais, visto que, a água é armazenada em períodos chuvosos, devido à porosidade das rochas calcárias, causando assim, uma irregularidade no regime dos rios superficiais.

38Com relação às sub-bacias da BHRF, o córrego Taquaral apresenta densidade de drenagem de 0,60 km/km², sobretudo, devido à extensão da sub-bacia, que possui 106,06 km² e exibe setores dotados de drenagens criptorreicas. Embora haja uma baixa densidade de drenagem, a sub-bacia do córrego Taquaral proporciona um dos menores coeficientes de manutenção (1.668,25 m²/m) e percurso superficial de todas as sub-bacias (1,67 km). O presente córrego é intermitente e transpassa relevos declivosos como em sua nascente e vales encaixados na serra da Bodoquena e, até mesmo, um relevo aplainado nas proximidades da sua foz no córrego Formosinho.

39Ainda que tenha uma abundância de sumidouros e ressurgências, seu manancial é de 4ª ordem, devido à quantidade de afluentes existentes na parte da serra da Bodoquena. Com isso, possui uma densidade hidrográfica de 0,15 segmentos fluviais por km², além de ser dotada de grandes extensões de culturas.

40O córrego Serradinho, que é afluente do córrego Anhuma, já possui um volume maior de água do que o anteriormente citado, mas foi classificado como 3ª ordem e ainda está 100% em terrenos cársticos, além disso, possui uma densidade de drenagem e densidade hidrográfica de 0,44 km/km² e de 0,09 cursos d’água por km² respectivamente. Este córrego não possui extensos tributários, assim, sua quantidade de bifurcações é baixa.

41O córrego Anhuma exibe uma característica diferenciada das demais supracitadas, pois existem rochas terrígenas em sua sub-bacia, contudo são áreas pequenas (11,74 km²), o que não traz impactos notáveis à morfometria de drenagem. Porém, já ocorre uma elevação para 0,16 cursos por km² de densidade hidrográfica e uma densidade de drenagem de 0,49 km/km². Há de se destacar uma grande concentração de mananciais de pouca extensão nas rochas terrígenas da formação Puga, que não são traduzidas na análise geral da sub-bacia devido a sua extensão. Esse córrego já é de 4ª ordem e possui uma das maiores relações de bifurcação de toda a BHRF (6,00) para com os seus canais de terceira e segunda ordem.

42O córrego Formosinho tem como afluentes os córregos Taquaral e Anhuma e possui, portanto, a maior sub-bacia hidrográfica da BHRF, retratando assim, uma boa quantidade de mananciais, mesmo que suas características cársticas não favoreçam a presença destes. Com relação à rede de drenagem, nesta comparação entre essas áreas, a referida sub-bacia apresentou uma densidade de drenagem de 0,46 km/km² (ou 460 metros de canal hídrico para cada km²).

43Também foi digno de nota, que a bacia hidrográfica do córrego Formosinho possui canais de primeira ordem curtos e runoffs consideráveis que refletem um coeficiente de manutenção alto para uma bacia hidrográfica de 414 km², delineados por um perímetro altamente irregular, assim como a BHRF, em que toda a região oriental é delimitada pelos divisores da serra da Bodoquena. Esse perímetro irregular remete a um índice de circularidade baixo (0,30). Um índice circular baixo como esse também corrobora à alta influência topográfica da sub-bacia, sendo esta, influenciada pelo escoamento hídrico e menores tendências de inundação.

44O córrego Restinga, em sua sub-bacia, denota apenas um canal fluvial, sem a presença de afluentes e, portanto, tampouco de bifurcações. Sua densidade hidrográfica alcançou 0,06 cursos d’água por km² e 0,33 km/km² de densidade de drenagem. Esses valores ainda são baixos, mesmo que essa sub-bacia ainda esteja assentada sobre rochas terrígenas em 19,84% de sua área.

45O córrego Saladeiro, localizado no médio curso da BHRF, afluente do córrego Bonito, ostenta um equilíbrio das áreas com afloramentos de rochas carbonatadas e terrígenas, porém, ainda não é possível constatar a diferença entre esses sistemas no tocante à quantidade de mananciais. Porém, seus únicos três afluentes estão sobre rochas terrígenas, o que eleva a densidade de drenagem para 0,37 km/km², principalmente pela pouca extensão da sub-bacia.

46Já o córrego Bonito denota uma elevação na densidade de drenagem, que chega a 0,46 km/km². Toda sua região central exibe rochas terrígenas, contudo, em seu alto e baixo curso, as rochas carbonatadas são preponderantes. Esses valores contrastam com os encontrados no córrego São João, visto que este possui 100% de rochas terrígenas.

47Assim, já é possível identificar as alterações que as rochas acarretam na quantidade de mananciais no córrego São João. Apesar da pouca extensão, é de 4ª ordem, com 18 mananciais em sua sub-bacia. Isso se traduz em uma densidade de drenagem elevada para a área de estudo, que a faz alcançar a 0,96 km/km² e uma densidade hidrográfica de 0,39 cursos d’água por km².

48O coeficiente de manutenção é inversamente proporcional à densidade de drenagem, alcançando assim, 1.046,79 metros de coeficiente, o que significa, segundo Christofoletti (1980), que quanto menor é o índice, menor será a área mínima necessária para a manutenção de um metro de canal de escoamento. Ou seja, quanto maior a quantidade de mananciais por km², menor será a área mínima de manutenção. Essa relação entre as propriedades das rochas foi primordial para compreender a morfometria de drenagem e foi nesse viés que a análise é importante para essa pesquisa, pois auxiliou na compreensão das características de cada sub-bacia que também se manifestaram na qualidade das águas.

49Ao prosseguir com a descrição dos dados, notou-se que o córrego Espírito Santo exibe, em sua maioria, rochas carbonatadas no alto e médio curso, o que reduz consideravelmente os valores de densidade hidrográfica e de drenagem (0,30 cursos d’água por km² e 0,56 km/km², respectivamente).

50Outro importante afluente do rio Mimoso é o córrego Barranco, que denotam menor vazão de água que o córrego Espírito Santo, todavia, é mais extenso. Sua sub-bacia abrange 49,12 km², destes, 52,23% estão sobrepostos em rochas terrígenas. Esses fatores trazem impactos na análise morfométrica, já que ocorre uma elevação na quantidade de mananciais da sub-bacia (17 no total) e, consequentemente, a elevação da densidade de drenagem para 0,88 km/km². Exibe ainda, um relevo menos acidentado que a sub-bacia do córrego Espírito Santo, mas suas rochas terrígenas favorecem a formação de afluentes temporários de 1ª ordem.

51O rio Mimoso, principal afluente do rio Formoso, retrata uma quantidade considerável de afluentes (74 ao todo), porém, em sua extensão de 250,93 km², não há grandes alterações na densidade de drenagem, que chega a 0,80 km/km². A sub-bacia aponta 57,59% de seus terrenos em sistemas cársticos, porém, o relevo exerce atuação expressiva, pois seus terrenos acidentados favorecem o aumento da quantidade de cursos fluviais, além de ser apontado como o menor índice de circularidade da BHRF, o que, segundo Christofoletti (1980), significa identificar e interpretar o processo de alargamento ou alongamento da bacia hidrográfica.

52Já o córrego Retiro (Figura 6) possui características distintas, com 100% de rochas terrígenas, além de uma morfometria igualmente díspar, em que a grande quantidade de recursos hídricos traz uma densidade de drenagem elevada (1,38 km/km²), a maior dentre todas as sub-bacias. Sua densidade hidrográfica também ficou elevada para a área de estudo, com 0,83 cursos por km², e seu coeficiente de manutenção chega a menores valores, ou seja, 724,62 metros são necessários para manutenção da perenidade dos cursos fluviais da bacia hidrográfica. Estes valores reduzidos traduzem a grande quantidade de mananciais existentes, que chega a 51 segmentos fluviais, dos quais o córrego Retiro (principal) é classificado como de 5ª ordem, segundo Horton (1945 apud Christofoletti, 1980).

53Segundo os dados provenientes da Tabela 1 e Figura 3, a BHRF possui 278 segmentos fluviais, destes, 73% são segmentos de primeira ordem com a ampla maioria de córrego temporários. Já o rio Formoso é um curso fluvial de 6ª ordem com uma relação de bifurcação maior entre os segmentos de 1ª e 2ª ordem, com uma densidade de drenagem de 6,64 km/km² e uma densidade hidrográfica de 0,21 cursos d’água por km². Esses baixos valores de densidade foram potencializados pelo alto e médio curso da BHRF que, por pertencerem a um sistema cárstico, ostentam menos recursos hídricos superficiais.

Conclusões

54Os dados morfométricos encontrados auxiliaram, portanto, na identificação das características e sua disponibilização ao longo das bacias hidrográficas, sobretudo na íntima relação entre as rochas terrígenas e carbonatadas com a quantidade de cursos fluviais. Por mais que na literatura se tenham afirmações dessa relação, são escassas as pesquisas que buscam analisar uma área com os dois sistemas e compará-los, o que levaria – ou não – a refutar tais afirmações.

55Atualmente, muito se discute sobre a morfometria da rede de drenagem, inversamente proporcional, pouco se têm sobre morfometrias em regiões cársticas, ainda menos buscando correlacionar essa análise com as rochas e o modelado cárstico. Na literatura, é amplamente discutido algumas características principais, como a predominância de drenagens criptorreicas em carstes bem desenvolvidos. Tais drenagens correm pelas fendas e fissuras das rochas, favorecendo uma rica e diversificada rede de drenagem superficial e subterrânea.

56Foi possível constatar que, mesmo o carste de Bonito ainda ser considerado recente, já é possível notar uma alteração significativa das drenagens superficiais nas bacias hidrográficas cársticas e terrígenas. Constatando que, no carste, as drenagens são escassas e espaçadas, exibindo grandes áreas sem qualquer tipo de drenagem superficial (estas correm em predominância, no subterrâneo), vide as áreas de poljes das regiões do médio curso da BHRF. Por outro lado, no baixo curso (sub-bacias predominantemente terrígenas) a quantidade de mananciais superficiais é notável, isso ficou evidente pela densidade de drenagem e hidrográfica de tais sub-bacias hidrográficas.

57A morfometria da rede de drenagem possibilita trazer para a bacia hidrográfica analisada, uma melhor gestão dos recursos hídricos. Compreendendo que no médio curso ocorrem poucos mananciais, são nessas áreas que ocorrem plantações, muitas delas atingindo as áreas de preservação, o que reduziria ainda mais os mananciais superficiais, interferindo no volume e na vazão dos mananciais de maior porte. O regime irregular das águas superficiais ficaria mais comprometido, o que impactaria não apenas nas questões ambientais, mas também sociais e econômicos, visto que Bonito é, reconhecidamente, um dos maiores polos turísticos brasileiros, sobretudo pela ligação com as águas superficiais cênicas.

58Diante disso, essa pesquisa tem como propósito, além de analisar e trazer informações sobre os recursos hídricos superficiais da BHRF, auxiliar na instituição de um banco de dados que poderá ser utilizado por órgãos responsáveis em gerir e planejar o município de Bonito, sendo um auxílio para possíveis tomadas de decisões, salvaguardando os recursos hídricos, sobretudo aqueles poucos situados em ambiente cárstico, que abastecem e mantêm o equilíbrio da BHRF.