Em termos gerais, a metodologia empregada compreendeu três etapas específicas e complementares, a saber: 1ª) definição de anos-padrão secos e chuvosos, 2ª) mapeamento morfológico da quantidade e área de barras e ilhas fluviais, 3ª) tratamento estatístico dos dados e discussão dos resultados.

Considerando que em regiões semiáridas a precipitação é a variável mais importante a definir um padrão para um sistema fluvial, observou-se que os regimes pluvial e fluvial no rio Jaguaribe apresentam características semelhantes, ou seja, os anos mais chuvosos correspondem aos momentos de maior energia no canal. Em contrapartida, nos anos menos chuvosos, observam-se os menores fluxos de energia no rio. Desse modo, para o estudo da evolução morfológica em barras e ilhas no rio Jaguaribe, optou-se pela determinação de anos-padrão secos (baixa energia) e chuvosos (alta energia), com o intuito de compreender as repercussões das variações pluviométricas na morfologia de barras e ilhas fluviais.

Assim, o processo de eleição de anos padrão se desenvolveu levando em conta um tratamento estatístico em torno do método de Tavares (1976), denominado Desvio Percentual das Precipitações Sazonais (DPPS), que propõe por meio de dados mensais, sazonais e anuais uma análise temporal a cerca da realidade de distribuição pluviométrica durante uma série histórica. O método propõe uma classificação de três tipos de anos-padrão: o ano seco, o habitual e o chuvoso.

Incialmente foram corrigidas as falhas para a cada uma das estações consideradas no estudo, ao passo que, a etapa seguinte consistiu na organização dos dados em torno dos totais mensais, respeitando seus limites anuais. Em seguida, os dados foram agrupados em períodos trimestrais, correspondentes a cada uma das estações do ano (verão, outono, inverno e primavera), permitindo assim, o calculo da precipitação média de cada uma das estações. Na sequência, calculou-se o desvio médio para cada estação do ano, a partir do valor da média trimestral subtraído pelo valor da média geral. Pelo desvio médio, chegou-se ao desvio percentual, que é resultado da divisão do desvio médio pela média geral do ano, multiplicado por 100 Conforme Tavares (1976) os desvios próximos de 0 (zero) são considerados habituais; aqueles que apresentaram grandes desvios positivos são os chuvosos;  e os de grandes desvios negativos correspondem aos secos.

Estabelecidos os anos-padrão e considerando a disponibilidade de imagens de satélite para a área, optou-se por estudar os anos de 1958, 1984, 2001, 2009, 2014. A escolha dos anos citados permitiu avaliar a evolução morfológica de barras e ilhas numa perspectiva temporal de 56 anos, com ênfase nas variações do regime pluviométrico (anos secos e chuvosos), concomitantemente.

O ano de 1958 marca o início da disponibilidade de imagens para a área, e mesmo não tendo sido considerado na determinação dos anos-padrão (em razão da falha de dados na série histórica), estudos anteriores realizados Cavalcante (2012) o evidenciam como um ano bastante seco. Conforme a aplicação da metodologia dos anos-padrão, 1984 foi definido como um ano chuvoso (cuja precipitação média anual foi de 1144,9mm ano^-1 e a vazão média anual registrada foi 153,5m³ s^-1); o ano de 2001 foi considerado seco (cuja precipitação média anual foi de 437,8mm ano^-1 e a vazão média anual foi 8,62m³ s^-1); o ano de 2009 foi considerado chuvoso (cuja precipitação média anual foi elevada, em torno de 1305mm ano^-1, e a vazão média anual foi de 195, m³ s^-1); 2014 foi considerado um ano seco (cuja precipitação média anual foi de 448,2mm ano^-1 e a vazão média anual foi de 10,3m³ s^-1) e marca o mapeamento mais recente da área.

 Para o ano de 1958, foram utilizadas fotografias aéreas cedidas pelo Instituto do Desenvolvimento Agrário do Ceará-IDACE na escala de 1:25.000; para os anos de 1984, 2001, 2009  foram  utilizadas imagens de satélite do sensor Landsat 5/Thematic Mapper (TM),  disponibilizadas pelo Instituto Nacional de Pesquisas Espaciais-INPE; e para o ano de 2014 foram utilizadas imagens do Landsat 8, sensor Operational Land Imager (OLI),  obtidas junto ao Serviço Geológico Americano-USGS.

O passo seguinte consistiu no georreferenciamento das imagens para compatibilizar as fotografias aéreas com as imagens de satélite do Landsat 5 e do Landsat 8. A sobreposição destas imagens possibilitou uma análise da mobilidade das ilhas e barras fluviais, permitindo assim inferir algumas questões sobre a dinâmica morfológica do canal. Cabe salientar que as datas das imagens correspondem ao período pós-estação chuvosa na região, quando se verificou uma menor dinâmica de fluxos no rio, bem como uma menor cobertura de nuvens nas imagens. Ademais, o processamento das imagens de satélite e das fotografias aéreas, bem como a quantificação das feições morfológicas, foram realizados utilizando-se as ferramentas do software ArcGIS 10.1.

As diferentes tipologias de barras e ilhas fluviais identificadas no canal do rio Jaguaribe foram classificadas tomando como base as propostas de Charlton (2007), Stevaux et al. (1992), Stevaux (1993), Morais (2002, 2006) e Leli (2015), considerando, para tanto, as condições especificas de semiaridez do rio Jaguaribe.

3. Resultados e discussões

Inicialmente foram analisados os dados morfológicos referentes ao total de barras fluviais. Assim, no ano de 1958, considerado seco, foram contabilizadas 154 barras fluviais recobrindo uma área de aproximadamente 30,5km², ao passo que, no ano de 1984, considerado chuvoso, essas feições passaram a ocupar uma área de 37,1km², distribuída em 203 barras, evidenciando, portanto, um aumento de 31,8% na quantidade e 21,6% na área de barras fluviais.

No ano de 2001, considerado um ano seco e, portanto, de baixos fluxos no canal, as barras apresentaram uma redução em relação ao ano de 1984, de modo que se verificou a existência de 157 barras fluviais, ocupando uma área de cerca de 30 km², o que relevou uma redução percentual na quantidade e na área de 22,6% e 19,1%, respectivamente. Como se pode observar, o aumento da precipitação e da vazão entre os anos de 1958 e 1984 e, em seguida, a redução entre os anos de 1984 e 2001 mantiveram relação direta com a dinâmica de formação de barras fluviais no canal.

Em 2009, verificou-se um novo ciclo de agradação no canal, evidenciado pelo incremento quantitativo e areal das feições, de forma que foram identificadas 205 barras recobrindo uma área no canal de cerca de 36km², o que percentualmente revelou um aumento quantitativo e areal de 30,5% e 20%, respectivamente, em relação ao ano de 2001. Em anos chuvosos, como 2009, por exemplo, acredita-se que, mesmo após o período de picos de fluxos, os materiais depositados lateralmente e na parte superior (top bar) das barras fluviais possam ser expostos à ação de pequenas correntes erosivas que remodelam essas feições, levando-as a adquirir diferentes morfologias.  

No ano de 2014, ainda que as barras fluviais tenham apresentado uma redução quantitativa, totalizando 153 feições (52 barras fluviais a menos que no ano de 2009), registrou-se um sensível incremento areal de 1,2km², ocupando uma área de 37,2km². Assim, entre 2009 e 2014, as barras fluviais apresentaram um decréscimo percentual em quantidade de 25,3% e um acréscimo percentual em área de 3,3%. Nota-se que, apesar do intervalo temporal de análise, entre 2001 e 2014 ser reduzido (13 anos), o retrabalhamento das feições manteve-se praticamente igual àquele verificado no período de 1984 e 2001 (portanto, um intervalo temporal maior, de 20 anos), quando se verificou uma média de formação de aproximadamente 48 a 50 barras entre um período e outro. Essa dinâmica indica que no rio Jaguaribe o ciclo de vida das barras fluviais é curto, em função da rápida mobilidade dos sedimentos entre o período seco e chuvoso, fazendo com que essas feições não se estabilizem por longos períodos no canal, sobretudo em eventos de cheias. 

Sumarizando, nota-se que a dinâmica de formação de barras fluviais no baixo curso do rio Jaguaribe se dá em consonância com a dinâmica das precipitações e vazões, de modo que, nos anos considerados secos (2001 e 2014) observou-se uma redução das feições devido à ausência de fluxos significativos para promover o transporte de sedimentos, ao passo que, em 1984 e 2009 houve um incremento de barras fluviais, sobretudo pela entrada de novos fluxos no canal, que, pela ação cíclica dos processos erosivos e sedimentares, promoveram a redistribuição dos sedimentos no canal, culminando em feições morfológicas de diferentes tamanhos e formas geométricas, a saber, barras laterais, centrais, longitudinais, em pontal. 

Em seguida, foram analisados dos dados morfológicos relacionados às barras laterais e centrais. Os resultados mostraram que, no ano de 1958, considerado seco, as barras laterais somavam 107 feições ocupando uma área de 19,1km², enquanto, em 1984, foram contabilizadas 147 feições recobrindo uma área de 29,6km², revelando um aumento percentual de 37,3% na quantidade e 55% na área de barras laterais.

 Em contrapartida, no ano de 2001, considerado seco, as barras laterais apresentaram um decréscimo quantitativo e areal, de modo que foram contabilizadas 113 feições laterais ocupando uma área de 27,2km². Em termos percentuais, os dados demonstraram uma diminuição de 23,1% na quantidade e 8,1% na área dessas feições em relação ao ano de 1984. Decerto, a redução de barras laterais entre 1984 e 2001 está relacionada ao fato de que os baixos fluxos propiciaram o desenvolvimento de vegetação na superfície das barras, favorecendo a anexação destas à planície de inundação.

 Contudo, no ano de 2009, as barras laterais voltaram a apresentar um aumento em função dos elevados fluxos registrados, de modo que foram identificadas 133 feições, abrangendo uma área de 27,7km². Esses valores equivalem a um aumento percentual de 17,7% na quantidade e 2,5% na área em relação ao ano de 2001.

No ano de 2014, considerado seco, as barras laterais apresentaram uma redução na quantidade, de modo que, foram mapeadas 122 feições, no entanto, com uma área de 30,7 km², representando em termos absolutos, uma redução quantitativa de 11 barras, e um aumento areal de 3km².

Outra investigação importante realizada nesse momento, diz respeito aos dados do coeficiente de correlação linear e de determinação estimados para a relação vazão, precipitação e quantidade de barras laterais. A relação entre a vazão média anual e a quantidade total de barras laterais revelou um coeficiente de correlação de 0,8052 e um coeficiente de determinação de 0,6484 (R²= 65%), ao passo que para a relação entre a precipitação média anual e a quantidade total de barras laterais, obteve-se um coeficiente de correlação de 0,8234 e um coeficiente de determinação de 0,678 (R²=68%).

Os valores observados para os anos de 1958 e 1984 revelaram que o retrabalhamento das barras centrais pela ação dos processos erosivos contribuiu para o incremento quantitativo, no entanto, as feições apresentaram menores dimensões em termos de área. A justificativa para isso pode estar associada ao fato de que as barras de tamanho médio e grande, desprovidas de vegetação, possam ter se fragmentado em barras menores, o que contribuiu para um incremento quantitativo de feições em 1984. De todo modo, essa dinâmica observada em 1984 revelou um canal bastante entrelaçado, cujo surgimento de novas barras centrais esteve diretamente associado à divisão da linha de fluxo.

Em 2001, cuja precipitação e vazão médias anuais foram reduzidas, as barras centrais apresentaram uma diminuição na quantidade e na área, de forma que foram contabilizadas 44 feições (12 feições a menos que em 1984) ocupando uma área de 2,9km², revelando, portanto, uma redução percentual de 21,4% na quantidade e 60,8% na área. Esse decréscimo na quantidade de barras centrais verificado em 2001 indica uma sensível tendência de retificação do canal, com a diminuição do número de fluxos secundários, e a incorporação de barras centrais às margens do rio.

No ano de 2009, considerado chuvoso, verificou-se uma nova fase de agradação no rio, evidenciada pelo acréscimo das feições, que passaram a ocupar uma área de 8,4km², correspondendo a 72 barras centrais, que equivalem a um aumento percentual de 63,6% na quantidade e de 189,6% na área, denotando, portanto, que o canal encontrava-se mais entrelaçado por barras e ilhas fluviais em relação ao ano de 2001. Porém, no ano de 2014, essas feições voltaram a apresentar uma redução, de maneira que foram mapeadas 31 feições centrais ocupando uma área de 6,4km², equivalente a uma redução percentual de 56,9% na quantidade e 23,8% na área de barras centrais.

Considerando que o índice de entrelaçamento do canal por barras fluviais aumenta com a carga de sedimentos e diminui com a presença de vegetação (COULTHARD, 2005), é possível presumir que em anos secos, como 1958, 2001 e 2014, por exemplo, os baixos fluxos associados ao desenvolvimento da vegetação sobre as barras centrais tenderam a estabilizar o padrão de canal, diminuindo a erosão e a disponibilidade de sedimentos (inibindo a formação de novas barras centrais), e, portanto reduzindo o índice de entrelaçamento. Por outro lado, em anos chuvosos como 1984 e 2009, o entrelaçamento do canal é maior em virtude do maior aporte de água e sedimentos.

A relação entre a vazão média anual e quantidade total de barras centrais revelou um coeficiente de correlação de 0,9262 e um coeficiente de determinação de 0,8579 (R²= 85%), ao passo que para a relação precipitação média anual e a quantidade total de barras laterais obteve-se um coeficiente de correlação de 0,9178 e um coeficiente de determinação de 0,8425 (R²=84%). Tanto para as barras laterais, quanto para as barras centrais, os valores da correlação linear de Pearson foram superiores a 0,5, indicando uma associação linear forte positiva entre as variáveis.

O mapeamento mostrou a predominância de barras laterais e centrais ao longo de toda a escala temporal, especialmente nos segmentos que drenam os municípios de São João do Jaguaribe, Tabuleiro do Norte (Distrito de Peixe Gordo) e Limoeiro do Norte. Em adição, observou-se que as barras laterais e centrais de menor porte apresentaram uma maior mobilidade ao longo da escala temporal estudada, ao passo que as de médio e grande porte, especialmente as que se encontravam vegetadas, demonstraram um retrabalhamento preferencial nas bordas, exceto em períodos de grandes fluxos no canal. Do ponto de vista comparativo, acredita-se que, mesmo nos períodos de cheias, as barras de grande porte em rios de áreas semiáridas sofram de modo mais acentuado a erosão nas bordas quando comparadas às barras de grande porte em rios de áreas úmidas/subúmidas, haja vista a intensidade dos fluxos fluviais durante o período chuvoso. Esse processo contribui, inclusive, para que ocorram alterações na morfologia das barras centrais e laterais, mudando completamente a geometria entre os períodos seco e chuvoso.

 Essa dinâmica observada no rio Jaguaribe corrobora com a ideia de Stevaux (1992) de que, dependendo da dimensão areal das barras, o comportamento deposicional parece ser mais significativo para a evolução de barras laterais, enquanto, para as barras centrais, os processos erosivos são mais importantes, contribuindo para um retrabalhamento contínuo.

As barras em pontal foram mapeadas em trechos onde o canal apresenta um padrão meândrico, permitindo avaliar as condições hidrodinâmicas que se desenvolvem nas margens côncavas e convexas das curvas do rio. No caso do rio Jaguaribe, a dinâmica de sedimentação observada em ambas as margens (côncava e convexa) contradiz o esquema teórico tradicional acerca do comportamento morfológico dos meandros, demonstrando que a acreção lateral não se produz exclusivamente nas margens convexas.

Quanto aos dados, verificou-se que no ano de 1958 foram contabilizadas 19 barras em pontal, ocupando uma área de 4,2km², ao passo que, em 1984, ocorreu um incremento de 26 feições ocupando uma área de aproximadamente 5,7km², distribuídas ao longo dos trechos sinuosos do canal. Em termos percentuais, verificou-se um aumento de 35,7% na área e 36,4% no número de feições.

Contudo, no ano de 2009, as barras em pontal voltaram a apresentar significativo acréscimo, de maneira que foram mapeadas 36 feições, abrangendo uma área de 9,4km², revelando, portanto, um aumento percentual em quantidade e área da ordem de 125% e 100%, respectivamente, em relação ao ano de 2001. No ano de 2014, essas feições apresentaram um novo acréscimo, sendo contabilizadas 39 feições, ocupando uma área de 9,6km², o que indicou um aumento percentual de 8,3% na quantidade e 2,1% na área dessas feições. Acredita-se que uma parcela dos sedimentos que adentraram o canal para a formação dessas barras em pontal sejam provenientes do próprio aporte do rio (sobretudo em 2009, por se tratar de um ano chuvoso, cuja dinâmica de fluxos e sedimentos foi intensa), bem como em decorrência de processos erosivos nas margens do rio em trechos desprovidos de vegetação.

Nota-se que, no ano de 2014, considerado seco, o comportamento das barras em pontal não seguiu o ritmo dos demais anos secos, como, por exemplo, 1958 e 2001, onde foi registrada uma redução das barras em pontal, indicando que a formação dessas tipologias de barras está associada ainda a outros fatores, sobretudo de caráter geomorfológico. Os dados de correlação linear e do coeficiente de determinação corroboraram para o entendimento desta questão. Assim, para a relação entre a vazão e a quantidade de barras em pontal verificou-se um coeficiente de correlação de 0,2659 e um coeficiente de determinação de 0,0708, ao passo que para a relação entre a precipitação e a quantidade de barras em pontal, o coeficiente de correlação e o coeficiente de determinação obtidos foram da ordem de 0,2560 e 0,0656. Para ambos os casos, os dados revelaram uma correlação positiva fraca (< 0,5), ou seja, vazão e precipitação não são variáveis determinantes na formação de barras em pontal.

Considerando que os dados de correlação linear são insuficientes para justificar o aumento na quantidade de barras em pontal no baixo curso do rio Jaguaribe, efetuou-se a estimativa do índice de sinuosidade para cada um dos anos-padrão em estudo. Assim, o índice de sinuosidade foi de 1,32 para o ano de 1958, 1,33 para 1984, 1,34 para o ano de 2001, 1,38 para 2009 e 1,40 para o ano de 2014. Com base nos dados obtidos, pode-se afirmar que o aumento na quantidade de barras em pontal entre 1958 e 2014 se justifica em função do aumento da sinuosidade do rio ao longo da escala temporal, revelando um canal mais curvilíneo, portanto de fluxos mais lentos, favorecendo, assim, a sedimentação de barras do tipo em pontal. 

A análise da evolução morfológica das ilhas fluviais no baixo curso do rio Jaguaribe ao longo da escala temporal considerada evidenciou mudanças significativas, como o aparecimento de novas ilhas, a presença de ilhas residuais que aumentaram e reduziram suas áreas, e o desaparecimento de ilhas por incorporação à planície de inundação.

Em 1984, foram identificadas 13 ilhas fluviais, recobrindo uma área de 2,3km², ao passo que, em 2001, apesar da redução do número de ilhas, observou-se um aumento areal de 0,4 km² (40 ha), uma vez que foram identificadas 12 feições, abrangendo uma área de 2,7km². Assim, entre os anos de 1984 e 2001, notou-se um decréscimo percentual na quantidade de 7,6% e um acréscimo na área de 17,3%. Possivelmente, o desaparecimento de uma única ilha fluvial entre os anos de 1984 e 2001 está relacionado a dois processos pontuais: 1) aos eventos de alta energia verificados nesse intervalo temporal, especificamente nos anos de 1985 e 1989, cujas precipitações médias anuais registradas foram de 1907,7mm ano^-1 e 1046mm ano^-1, respectivamente, e que podem ter ocasionado a remoção da cobertura vegetal da ilha e, consequentemente, a remobilização do material sedimentar e 2) que a ilha fluvial tenha se anexado à planície de inundação nesse intervalo temporal.

No ano de 2009, foram mapeadas 18 feições, recobrindo uma área 2,8km², representando, portanto, um aumento quantitativo de 6 ilhas fluviais e um acréscimo areal de 0,1km²  em relação ao ano de 2001. Entre 2001 e 2009, houve um acréscimo percentual na quantidade e na área de 50% e 3,6%, respectivamente. No ano de 2014, observou-se um aumento considerável de ilhas fluviais, uma vez que foram contabilizadas 28 feições recobrindo uma área de 3,2km², representando, em relação a 2009, um aumento percentual quantitativo de 55,6% e 14,3% na área. De modo direto, nota-se que a evolução de ilhas fluviais nos anos-padrão não corresponde às variações de precipitação e vazão, o que releva a existência de outros fatores na formação dessas feições.

Esse mesmo processo foi verificado por Ham e a Church (2002) estudando o rio Fraser, no Canadá, onde se constatou que a redução do regime de cheias e a ausência de grandes fluxos entre 1977 e 2002, portanto 25 anos, possibilitaram o crescimento da vegetação na superfície das barras, propiciando a formação de novas ilhas fluviais. Outros estudos desenvolvidos no Rio Grande, no Novo México (MEYER e HEPLER, 2007) evidenciaram esse mesmo processo, ou seja, os autores constataram que os fluxos elevados entre o período de 1973 a 1985 foram capazes de perturbar a estabilidade das ilhas fluviais, ocasionando a remoção da vegetação e a redistribuição dos sedimentos no rio. Contudo, os baixos fluxos registrados a partir de 1986 possibilitaram a fixação da vegetação, ocasionando um aumento de ilhas fluviais de 50 (em 1985) para 203 em 1992. Os baixos fluxos registrados em anos subsequentes, especialmente entre o período de 1990 a 2000, favoreceram ainda mais a estabilização das ilhas fluviais.

Em suma, considerando os fatores de formação de ilhas fluviais citados por Osterkamp (1998), pode-se afirmar que, no caso do rio Jaguaribe, em seu baixo curso, três destes fatores parecem ser mais relevantes, a saber: 1) a partir da evolução de barras centrais; 2) processos de avulsão no canal e 3) a construção de reservatórios, a exemplo da barragem do Castanhão.

Por último, verificou-se que as ilhas fluviais de pequenas dimensões ocorrem predominantemente no inicio do baixo curso, enquanto as de maiores dimensões ocorrem a partir do município de Itaiçaba em direção à foz do rio. Assim, quando comparadas às grandes ilhas, as pequenas ilhas apresentam menor estabilidade no sistema fluvial, sofrendo um retrabalhamento mais intenso pela ação dos processos erosivos. Em contrapartida, as grandes ilhas, como a do Pinto, do Caldereiro, dos Veados e a ilha Grande, são praticamente estáveis, de modo que o retrabalhamento se dá frequentemente nas bordas dessas feições. Assim, não apenas a vegetação presente na ilha como também a existência de troncos, galhos e folhas auxiliam na estabilização dos sedimentos mais finos, fazendo com que a área de deposição da ilha aumente significativamente, levando-a a adquirir maior estabilidade no sistema fluvial.

4. Conclusões

Os resultados obtidos neste estudo indicaram que a dinâmica quantitativa de barras fluviais no baixo curso do rio Jaguaribe se dá em consonância com a dinâmica de vazões, de modo que, nos anos considerados secos, observou-se uma redução das feições devido à ausência de fluxos significativos para promover o transporte de sedimentos, ao passo que, em anos chuvosos, houve um incremento na quantidade de barras fluviais (laterais e centrais) em decorrência de novos ciclos de agradação no canal fluvial. Tal afirmação pôde ser ratificada com base nos dados de correlação linear de Pearson, que evidenciaram uma correlação linear positiva em boa parte dos cenários analisados.

A exceção a essa regra foi observada nas barras em pontal, cujo comportamento quantitativo não correspondeu às variações de vazão nos diversos anos estudados, de modo que os dados de correlação linear revelaram uma associação positiva fraca, indicando a influência de outros fatores de caráter geomorfológico na elaboração dessas feições. Logo, tal cenário pode ser justificado, em parte, pelo aumento da sinuosidade do rio ao longo da escala temporal, revelando um canal mais curvilíneo, portanto, de fluxos mais lentos, o que favoreceu a sedimentação de barras do tipo em pontal. 

No caso das ilhas fluviais, constatou-se um incremento dessas feições nos últimos anos no baixo curso do rio Jaguaribe, fato que pode estar relacionado à regularização das vazões exercida pela barragem do Castanhão e pelas diversas passagens molhadas presentes ao longo do rio, associada, sobretudo, à ausência de grandes cheias entre esse período. Assim, a ausência de grandes fluxos nos últimos anos contribuiu para que a tensão de cisalhamento fosse mínima, sendo insuficiente para superar os efeitos de ligação da vegetação e, consequentemente, a redistribuição dos sedimentos na ilha. Assim, a redução dos picos de cheias e as descargas de baixos fluxos liberadas pela barragem forneceram a água necessária para o crescimento da vegetação e a colonização e estabilização das barras fluviais, favorecendo, assim, a formação de novas ilhas fluviais.

5. Referências bibliográficas